Velocity-equilibrium model
Velocity-equilibrium and mechanical-equilibrium flows are solved in ECOGEN using the velocity-equilibrium model [SCazeP+21]. For N phases involved and without any extra physics (surface tension, viscosity…), this model reads:
where \(\alpha_k\), \(\rho_k\), \(p_k\) and \(e_k\) are the volume fraction, density, pressure and internal energy of each phase, respectively, and for which \(k\) indicates the phase index. The mixture density and pressure are
while the mixture total energy is
where \(e\) is the mixture specific internal energy
\(e_k \left( \rho_k , p_k \right)\) is defined via an equation of state (EOS) and \(Y_k\) are the mass fractions
The relaxation of pressures between the phases is
where \(j\) are phases different from \(k\) and \(\mu_{k,j}\) are the pressure-relaxation coefficients related to the \(k\)–\(j\) interactions. Herein, the pressure-relaxation coefficient \(\mu\) is considered the same for each phase combination. The interfacial pressure is defined as
where \(z_k = \rho_k c_k\) and \(c_k\) are the acoustic impedance and speed of sound of the phase \(k\), respectively.
Since pressures are in disequilibrium here, the total energy equation of the mixture is replaced by the internal-energy equation for each phase. Nevertheless, conservation of the mixture total energy can be written in its usual form
We note that this equation is redundant when the internal energy equations are also computed. However, in practice, we include it in our computations to ensure that the total energy is numerically conserved, and thus preserve a correct treatment of shock waves.
This model is solved thanks to the numerical method presented in [SCazeP+21] where infinite as well as finite pressure-relaxation rates are possible.
Tests cases
Tests are provided with ECOGEN package and may be described in details later.
<testCase>libTests/referenceTestCases/UEq/1D/transports/interfaceWaterAir/</testCase>
<testCase>libTests/referenceTestCases/UEq/1D/shockTubes/interfaceAirHelium/</testCase>
<testCase>libTests/referenceTestCases/UEq/1D/shockTubes/interfaceWaterAir/</testCase>
<testCase>libTests/referenceTestCases/UEq/1D/shockTubes/interfaceWaterAirNASG/</testCase>
<testCase>libTests/referenceTestCases/UEq/1D/shockTubes/epoxySpinel/</testCase>
<testCase>libTests/referenceTestCases/UEq/1D/shockTubes/mixtures/</testCase>
<testCase>libTests/referenceTestCases/UEq/1D/mixture/waterAir/</testCase>
<testCase>libTests/referenceTestCases/UEq/1D/shockOnInterface/sharpInterfaceWaterAir/</testCase>
<testCase>libTests/referenceTestCases/UEq/1D/shockOnInterface/diffusedInterfaceWaterAir/</testCase>
<testCase>libTests/referenceTestCases/UEq/1D/cavitation/</testCase>
<testCase>libTests/referenceTestCases/UEq/1D/sphericalCollapse/Pratio1427/</testCase>
<testCase>libTests/referenceTestCases/UEq/2D/nonSphericalCollapseNearWall/</testCase>
<testCase>libTests/referenceTestCases/UEq/2D/squareToCircleSymmetry/</testCase>