43 int main(
int argc,
char *argv[])
51 #include "createFields.H" 52 #include "createFieldRefs.H" 66 scalar meanCoNum = 0.0;
68 Info<<
"\nStarting time loop\n" <<
endl;
111 max(
max(phiv_pos + cSf_pos, phiv_neg + cSf_neg), v_zero)
116 min(
min(phiv_pos - cSf_pos, phiv_neg - cSf_neg), v_zero)
141 amaxSf =
max(
mag(aphiv_pos),
mag(aphiv_neg));
148 #include "setRDeltaT.H" 152 #include "setDeltaT.H" 157 Info<<
"Time = " << runTime.timeName() <<
nl <<
endl;
159 phi = aphiv_pos*rho_pos + aphiv_neg*rho_neg;
163 (aphiv_pos*rhoU_pos + aphiv_neg*rhoU_neg)
164 + (a_pos*p_pos + a_neg*p_neg)*
mesh.Sf()
170 aphiv_pos*(rho_pos*(e_pos + 0.5*
magSqr(U_pos)) + p_pos)
171 + aphiv_neg*(rho_neg*(e_neg + 0.5*
magSqr(U_neg)) + p_neg)
172 + aSf*p_pos - aSf*p_neg
187 U.correctBoundaryConditions();
188 rhoU.boundaryFieldRef() ==
rho.boundaryField()*
U.boundaryField();
209 & (a_pos*U_pos + a_neg*U_neg)
220 e.correctBoundaryConditions();
222 rhoE.boundaryFieldRef() ==
225 e.boundaryField() + 0.5*
magSqr(
U.boundaryField())
242 p.correctBoundaryConditions();
243 rho.boundaryFieldRef() ==
psi.boundaryField()*
p.boundaryField();
249 Info<<
"ExecutionTime = " << runTime.elapsedCpuTime() <<
" s" 250 <<
" ClockTime = " << runTime.elapsedClockTime() <<
" s" autoPtr< compressible::turbulenceModel > turbulence
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
GeometricField< tensor, fvPatchField, volMesh > volTensorField
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
static tmp< GeometricField< typename innerProduct< vector, Type >::type, fvsPatchField, surfaceMesh > > dotInterpolate(const surfaceVectorField &Sf, const GeometricField< Type, fvPatchField, volMesh > &tvf)
Interpolate field onto faces.
GeometricField< vector, fvsPatchField, surfaceMesh > surfaceVectorField
dimensionedScalar sqrt(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Read the control parameters used by setDeltaT.
tmp< surfaceScalarField > interpolate(const RhoType &rho)
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
dimensionedScalar neg(const dimensionedScalar &ds)
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
psiReactionThermo & thermo
dimensionedScalar pos(const dimensionedScalar &ds)
const dimensionSet dimVolume(pow3(dimLength))
const dimensionedScalar e
Elementary charge.
dimensioned< scalar > magSqr(const dimensioned< Type > &)
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
dimensionedSymmTensor dev2(const dimensionedSymmTensor &dt)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const dimensionedScalar c
Speed of light in a vacuum.
Calculates the mean and maximum wave speed based Courant Numbers.
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
dimensioned< scalar > mag(const dimensioned< Type > &)
const volScalarField & psi
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Execute application functionObjects to post-process existing results.
word fluxScheme("Kurganov")
Provides functions smooth spread and sweep which use the FaceCellWave algorithm to smooth and redistr...
Read the control parameters used by setDeltaT.
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)