60 int main(
int argc,
char *argv[])
76 Info<<
"\nStarting time loop\n" <<
endl;
121 max(
max(phiv_pos + cSf_pos, phiv_neg + cSf_neg), v_zero);
123 min(
min(phiv_pos - cSf_pos, phiv_neg - cSf_neg), v_zero);
147 amaxSf =
max(
mag(aphiv_pos),
mag(aphiv_neg));
155 Info<<
"Time = " << runTime.timeName() <<
nl <<
endl;
158 phiv_pos = U_pos &
mesh.
Sf();
159 phiv_neg = U_neg &
mesh.
Sf();
165 aphiv_pos = phiv_pos - aSf;
166 aphiv_neg = phiv_neg + aSf;
168 phi = aphiv_pos*rho_pos + aphiv_neg*rho_neg;
171 (aphiv_pos*rhoU_pos + aphiv_neg*rhoU_neg)
172 + (a_pos*p_pos + a_neg*p_neg)*
mesh.
Sf();
175 aphiv_pos*(rho_pos*(e_pos + 0.5*
magSqr(U_pos)) + p_pos)
176 + aphiv_neg*(rho_neg*(e_neg + 0.5*
magSqr(U_neg)) + p_neg)
177 + aSf*p_pos - aSf*p_neg;
188 U.dimensionedInternalField() =
189 rhoU.dimensionedInternalField()
190 /
rho.dimensionedInternalField();
191 U.correctBoundaryConditions();
192 rhoU.boundaryField() =
rho.boundaryField()*
U.boundaryField();
212 & (a_pos*U_pos + a_neg*U_neg)
225 rhoE.boundaryField() =
245 p.dimensionedInternalField() =
246 rho.dimensionedInternalField()
248 p.correctBoundaryConditions();
255 Info<<
"ExecutionTime = " << runTime.elapsedCpuTime() <<
" s"
256 <<
" ClockTime = " << runTime.elapsedClockTime() <<
" s"