58 int main(
int argc,
char *argv[])
74 Info<<
"\nStarting time loop\n" <<
endl;
139 amaxSf =
max(
mag(aphiv_pos),
mag(aphiv_neg));
147 Info<<
"Time = " << runTime.timeName() <<
nl <<
endl;
152 (aphiv_pos*rhoU_pos + aphiv_neg*rhoU_neg)
153 + (a_pos*p_pos + a_neg*p_neg)*
mesh.
Sf();
156 aphiv_pos*(rho_pos*(e_pos + 0.5*
magSqr(U_pos)) + p_pos)
157 + aphiv_neg*(rho_neg*(e_neg + 0.5*
magSqr(U_neg)) + p_neg)
158 + aSf*p_pos - aSf*p_neg;
168 U.dimensionedInternalField() =
169 rhoU.dimensionedInternalField()
170 /
rho.dimensionedInternalField();
171 U.correctBoundaryConditions();
172 rhoU.boundaryField() =
rho.boundaryField()*
U.boundaryField();
194 & (a_pos*U_pos + a_neg*U_neg)
207 rhoE.boundaryField() =
227 p.dimensionedInternalField() =
228 rho.dimensionedInternalField()
230 p.correctBoundaryConditions();
235 Info<<
"ExecutionTime = " << runTime.elapsedCpuTime() <<
" s"
236 <<
" ClockTime = " << runTime.elapsedClockTime() <<
" s"