68int main(
int argc,
char *argv[])
72 "Solver for two compressible, non-isothermal, immiscible fluids"
73 " using VOF phase-fraction based interface capturing approach.\n"
74 "Supports dynamic mesh motions including overset cases."
86 #include "createDyMControls.H"
87 #include "createFields.H"
94 #include "correctPhi.H"
99 #include "CourantNo.H"
100 #include "setInitialDeltaT.H"
108 Info<<
"\nStarting time loop\n" <<
endl;
112 #include "readControls.H"
116 #include "setRDeltaT.H"
120 #include "CourantNo.H"
121 #include "alphaCourantNo.H"
122 #include "setDeltaT.H"
134 scalar timeBeforeMeshUpdate =
runTime.elapsedCpuTime();
140 Info<<
"Execution time for mesh.update() = "
141 <<
runTime.elapsedCpuTime() - timeBeforeMeshUpdate
146 if (
mesh.topoChanging())
159 localMin<scalar>(
mesh).interpolate(cellMask.oldTime());
171 Uf.boundaryFieldRef()[patchI] =
172 Uint.boundaryField()[patchI];
180 #include "correctPhi.H"
186 faceMask = localMin<scalar>(
mesh).interpolate(cellMask);
192 fvc::makeRelative(
phi,
U);
202 #include "alphaControls.H"
203 #include "compressibleAlphaEqnSubCycle.H"
231 runTime.printExecutionTime(Info);
CMULES: Multidimensional universal limiter for explicit corrected implicit solution.
const uniformDimensionedVectorField & g
const surfaceScalarField & ghf
const volScalarField & gh
const volScalarField & psi2
const volScalarField & psi1
surfaceScalarField faceMask(localMin< scalar >(mesh).interpolate(cellMask))
tmp< surfaceScalarField > talphaPhi1Corr0
Read the control parameters used by setDeltaT.
autoPtr< surfaceVectorField > Uf
Creates and initialises the velocity velocity field Uf.
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
GeometricField< scalar, fvPatchField, volMesh > volScalarField
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)
Adjust the balance of fluxes on the faces between interpolated and calculated to obey continuity.
Execute application functionObjects to post-process existing results.
Sets blocked cells mask field.
Sets blocked cells mask field.
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
#define forAll(list, i)
Loop across all elements in list.