66int main(
int argc,
char *argv[])
70 "Solver for two compressible, non-isothermal immiscible fluids"
71 " using VOF phase-fraction based interface capturing.\n"
72 "With optional mesh motion and mesh topology changes including"
73 " adaptive re-meshing."
81 #include "initContinuityErrs.H"
82 #include "createDyMControls.H"
83 #include "createFields.H"
85 #include "CourantNo.H"
86 #include "setInitialDeltaT.H"
94 Info<<
"\nStarting time loop\n" <<
endl;
105 #include "CourantNo.H"
106 #include "alphaCourantNo.H"
107 #include "setDeltaT.H"
119 scalar timeBeforeMeshUpdate =
runTime.elapsedCpuTime();
121 if (isA<dynamicRefineFvMesh>(
mesh))
133 if (isA<dynamicRefineFvMesh>(
mesh))
137 alpha2.correctBoundaryConditions();
139 rho.correctBoundaryConditions();
146 Info<<
"Execution time for mesh.update() = "
147 <<
runTime.elapsedCpuTime() - timeBeforeMeshUpdate
157 #include "correctPhi.H"
160 fvc::makeRelative(
phi,
U);
171 #include "alphaControls.H"
172 #include "compressibleAlphaEqnSubCycle.H"
194 runTime.printExecutionTime(Info);
CMULES: Multidimensional universal limiter for explicit corrected implicit solution.
const uniformDimensionedVectorField & g
const surfaceScalarField & ghf
const volScalarField & gh
const volScalarField & alpha1
const volScalarField & psi2
const volScalarField & alpha2
const volScalarField & psi1
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< scalar, fvPatchField, volMesh > volScalarField
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
constexpr char nl
The newline '\n' character (0x0a)
Execute application functionObjects to post-process existing results.
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
isoAdvection advector(alpha1, phi, U)