Bug: reactingTwoPhaseEulerFoam does not reproduce the analytical solution for the Sod Shock tube problem

Description:SodShockTube-issue-v2012.zip

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).

Solution:

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:

https://bugs.openfoam.org/view.php?id=3634

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