View Revisions: Issue #3634

Summary 0003634: multiphaseEulerFoam does not reproduce the analytical solution for the Sod Shock tube problem
Revision 2021-02-24 10:56 by henry
Description multiphaseEulerFoam (OpenFOAM 8) is not reproducing the exact solution for the single phase Sod shock tube benchmark. In particular, the temperature field is significantly different (see figure). This issue is not observed in older OpenFOAM versions. TwoPhaseEulerFoam is not presenting this problem and reactingEulerFoam is affected since OpenFOAM 6.

Please find attached the results for the test performed with different OpenFOAM versions ( dashed line representing the exact solution). The “first bad commit” is also reported in the “Additional information” section.

Summary for reactingTwoPhaseEulerFoam:
OF-4 - ok.
OF-5 - convergence issue (easily solved just re-adding fluid.correct() at the beginning of the loop in pEqn.H). Then the solution is ok.
OF-6 - not good (higher temperature)
OF-7 - not good (higher temperature)
OF-8 - multiphasePhaseEulerFoam: not good (higher temperature)
Revision 2021-02-24 09:54 by federico
Description multiphaseEulerFoam (OpenFOAM 8) is not reproducing the exact solution for the single phase Sod shock tube benchmark. In particular, the temperature field is significantly different (see figure). This issue is not observed in older OpenFOAM versions. TwoPhaseEulerFoam is not presenting this problem and reactingEulerFoam is affected since OpenFOAM 6. This issue is also observed with OpenFOAM v2006 (same behaviour).

Please find attached the results for the test performed with different OpenFOAM versions ( dashed line representing the exact solution). The “first bad commit” is also reported in the “Additional information” section.

Summary for reactingTwoPhaseEulerFoam:
OF-4 - ok.
OF-5 - convergence issue (easily solved just re-adding fluid.correct() at the beginning of the loop in pEqn.H). Then the solution is ok.
OF-6 - not good (higher temperature)
OF-7 - not good (higher temperature)
OF-8 - multiphasePhaseEulerFoam: not good (higher temperature)