overRhoSimpleFoam.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | www.openfoam.com
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8  Copyright (C) 2017 OpenCFD Ltd
9 -------------------------------------------------------------------------------
10 License
11  This file is part of OpenFOAM.
12 
13  OpenFOAM is free software: you can redistribute it and/or modify it
14  under the terms of the GNU General Public License as published by
15  the Free Software Foundation, either version 3 of the License, or
16  (at your option) any later version.
17 
18  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
19  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
20  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21  for more details.
22 
23  You should have received a copy of the GNU General Public License
24  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
25 
26 Application
27  overRhoSimpleFoam
28 
29 Group
30  grpCompressibleSolvers
31 
32 Description
33  Overset steady-state solver for compressible turbulent flow.
34 
35 \*---------------------------------------------------------------------------*/
36 
37 #include "fvCFD.H"
38 #include "dynamicFvMesh.H"
39 #include "fluidThermo.H"
41 #include "simpleControl.H"
42 #include "pressureControl.H"
43 #include "fvOptions.H"
44 #include "cellCellStencilObject.H"
45 #include "localMin.H"
46 #include "oversetAdjustPhi.H"
47 
48 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
49 
50 int main(int argc, char *argv[])
51 {
52  argList::addNote
53  (
54  "Overset steady-state solver for compressible turbulent flow"
55  );
56 
57  #define CREATE_MESH createUpdatedDynamicFvMesh.H
58  #include "postProcess.H"
59 
60  #include "setRootCaseLists.H"
61  #include "createTime.H"
62  #include "createUpdatedDynamicFvMesh.H"
63  #include "createControl.H"
64  #include "createFields.H"
65  #include "createFieldRefs.H"
66  #include "createFvOptions.H"
67  #include "initContinuityErrs.H"
68 
69  turbulence->validate();
70 
71  // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
72 
73  Info<< "\nStarting time loop\n" << endl;
74 
75  while (simple.loop())
76  {
77  Info<< "Time = " << runTime.timeName() << nl << endl;
78 
79  // Pressure-velocity SIMPLE corrector
80  #include "UEqn.H"
81  #include "EEqn.H"
82  #include "pEqn.H"
83 
84  turbulence->correct();
85 
86  runTime.write();
87 
88  runTime.printExecutionTime(Info);
89  }
90 
91  Info<< "End\n" << endl;
92 
93  return 0;
94 }
95 
96 
97 // ************************************************************************* //
runTime
engineTime & runTime
Definition: createEngineTime.H:13
simple
const dictionary & simple
Definition: readFluidMultiRegionSIMPLEControls.H:1
fvOptions.H
turbulence
Info<< "Reading field U\n"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar("pos", dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar("neg", dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\n"<< endl;autoPtr< compressible::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
Definition: createFields.H:94
fluidThermo.H
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:369
createFvOptions.H
setRootCaseLists.H
Foam::Info
messageStream Info
Information stream (stdout output on master, null elsewhere)
simpleControl.H
postProcess.H
Execute application functionObjects to post-process existing results.
pressureControl.H
localMin.H
Foam::nl
constexpr char nl
Definition: Ostream.H:404
createTime.H
dynamicFvMesh.H
fvCFD.H
cellCellStencilObject.H
oversetAdjustPhi.H
Adjust the balance of fluxes on the faces between interpolated and calculated to obey continuity.
turbulentFluidThermoModel.H