Bug: reactingTwoPhaseEulerFoam does not reproduce the analytical solution for the Sod Shock tube problem
reactingTwoPhaseEulerFoam (OpenFOAM v2012) is not reproducing the exact solution for the single phase Sod shock tube benchmark. In particular, the temperature field is significantly different (see figure, dashed line representing the exact solution). reactingMultiPhaseEulerFoam is also affected and twoPhaseEulerFoam should not (although not tested for this repository)This issue is not observed in old OpenFOAM versions (OF<5).
This Bug is now solved in multiphaseEulerFoam (OpenFOAM 8) and I strongly suggest to implement this fix also for reactingTwo/MultiPhaseEulerFoam (OF-v2012). Here all the details:
Steps to Reproduce:
1- Enter Sod shock tube test case for reactingTwoPhaseEulerFoam (OF-v2012,v2006). The analytical solution is also included. 2- ./Allclean; ./Allrun 3- paraview -> load state file -> ParaviewState/vsAnalytical.pvsm (paraview version 5.0-5.6) 4- select “thermo:rho.gas” and click “apply” to check gas density