INT: reacting*EulerFoam: correct handling of K for energy equation (#2033)

• Integration from OpenFOAM.org: (link).
• It was reported in #2033 (closed) that reactingTwoPhaseEulerFoam and reactingMultiphaseEulerFoamcould not reproduce the closed-form solutions available for Sod's shock tube test case.
• To this end, two test cases were setup (the link will be broken in the future):
  • reactingTwoPhaseEulerFoam
  • reactingMultiphaseEulerFoam
• Effects of the bug fix on these test cases are illustrated below (left: previous state, right: new state):

reactingTwoPhaseEulerFoam

reactingMultiphaseEulerFoam

Tutorials

• The bug fix changes the output of the existing tutorials:

new-logs.gz old-logs.gz

new-testLoopReport.gz old-testLoopReport.gz

Environment (development)

    base0    = base
    base1    = develop
    api      = 2102
    patch    = 210210
    HEAD     = 8abac1dc4d
    version  = com
    compiler = Clang (system)
             = clang version 9.0.1 
    mpi      = SYSTEMOPENMPI
             = mpirun (Open MPI) 1.10.7.0.5e373bf1fd
    OS       = Description:     openSUSE Leap 15.1
    opts     = linux64ClangDPInt32Opt