69 int main(
int argc,
char *argv[])
73 "Solver for two incompressible, isothermal, immiscible fluids with" 75 "using VOF (volume of fluid) phase-fraction based interface capturing," 76 " with optional dynamic mesh motion (including overset)\n" 77 "and mesh topology changes including adaptive re-meshing." 87 #include "initContinuityErrs.H" 88 #include "createFields.H" 97 IOobject::READ_IF_PRESENT,
105 #include "CourantNo.H" 106 #include "setInitialDeltaT.H" 115 Info<<
"\nStarting time loop\n" <<
endl;
126 #include "CourantNo.H" 127 #include "setDeltaT.H" 138 scalar timeBeforeMeshUpdate =
runTime.elapsedCpuTime();
144 Info<<
"Execution time for mesh.update() = " 145 <<
runTime.elapsedCpuTime() - timeBeforeMeshUpdate
168 #include "alphaControls.H" 184 #include "alphaEqnSubCycle.H"
Sets blocked cells mask field.
constexpr char nl
The newline '\n' character (0x0a)
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Ostream & endl(Ostream &os)
Add newline and flush stream.
interfaceProperties interface(alpha1, U, thermo->transportPropertiesDict())
tmp< volScalarField > rAU
Sets blocked cells mask field.
const surfaceScalarField & ghf
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Creates and initialises the velocity velocity field Uf.
const uniformDimensionedVectorField & g
Info<< "Reading field U\"<< 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\"<< endl;autoPtr< compressible::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
Correct phi on new faces C-I faces.
Calculates and outputs the mean and maximum Courant Numbers.
tmp< surfaceScalarField > absolute(const tmp< surfaceScalarField > &tphi, const volVectorField &U)
Return the given relative flux in absolute form.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const volScalarField & gh
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
messageStream Info
Information stream (stdout output on master, null elsewhere)
const dimensionSet dimMass(1, 0, 0, 0, 0, 0, 0)
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Execute application functionObjects to post-process existing results.
void makeRelative(surfaceScalarField &phi, const volVectorField &U)
Make the given flux relative.
CMULES: Multidimensional universal limiter for explicit corrected implicit solution.
Adjust the balance of fluxes on the faces between interpolated and calculated to obey continuity...
static constexpr const zero Zero
Global zero (0)