58int main(
int argc,
char *argv[])
62 "Transient solver for compressible turbulent flow.\n"
63 "With optional mesh motion and mesh topology changes."
69 #include "createDyMControls.H"
71 #include "initContinuityErrs.H"
72 #include "createFields.H"
73 #include "createMRF.H"
76 #include "createControls.H"
82 #include "compressibleCourantNo.H"
83 #include "setInitialDeltaT.H"
88 Info<<
"\nStarting time loop\n" <<
endl;
92 #include "readControls.H"
99 autoPtr<volScalarField> divrhoU;
107 fvc::div(fvc::absolute(
phi,
rho,
U))
114 #include "setRDeltaT.H"
118 #include "compressibleCourantNo.H"
119 #include "setDeltaT.H"
147 rhoUf() *= faceMaskOld;
152 rhoUf() += (1-faceMaskOld)*rhoUfint;
157 rhoUf().boundaryFieldRef()[patchI] =
158 rhoUfint.boundaryField()[patchI];
166 #include "correctPhi.H"
179 fvc::makeRelative(
phi,
rho,
U);
213 runTime.printExecutionTime(Info);
Bound the given scalar field if it has gone unbounded.
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
surfaceScalarField faceMask(localMin< scalar >(mesh).interpolate(cellMask))
Creates and initialises the velocity field rhoUf if required.
autoPtr< surfaceVectorField > rhoUf
compressible::turbulenceModel & turbulence
Provides functions smooth spread and sweep which use the FaceCellWave algorithm to smooth and redistr...
Calculates and outputs the mean and maximum Courant Numbers.
GeometricField< vector, fvsPatchField, surfaceMesh > surfaceVectorField
messageStream Info
Information stream (stdout output on master, null elsewhere)
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Ostream & endl(Ostream &os)
Add newline and flush stream.
constexpr char nl
The newline '\n' character (0x0a)
Sets blocked cells mask field.
#define forAll(list, i)
Loop across all elements in list.