Typo in twoPhaseEulerFoam

Summary

A typo seems to be placed in twoPhaseEulerFoam, EEqns.H file, on the energy equation for the second phase.

What is the current bug behaviour?

In openfoam/applications/solvers/multiphase/twoPhaseEulerFoam/EEqns.H

Lines 50 to 55

      + (
            he2.name() == thermo2.phasePropertyName("e")
          ? fvc::div(fvc::absolute(alphaPhi2, alpha2, U2), p)
          + p*fvc::ddt(**alpha1**)
          : -alpha2*dpdt
        )

What is the expected correct behavior?

The following correction should be made to match the equation for the first energy field.

Lines 50 to 55

      + (
            he2.name() == thermo2.phasePropertyName("e")
          ? fvc::div(fvc::absolute(alphaPhi2, alpha2, U2), p)
          + p*fvc::ddt(**alpha2**)
          : -alpha2*dpdt
        )

Relevant logs and/or images

None

Environment information

Not relevant

Possible fixes

Quite straightforward

assigned to @andy

closed via commit 7904a273

There will be slight differences due to the change. On multiphase/twoPhaseEulerFoam/laminar/fluidisedBed tutorial:

Original code:

Fixed code:

added community label

mentioned in commit 7904a273