Extension of the Entropy Viscosity Method to the Multi-D Euler Equations and the Seven-Equation Two-Phase Model
The work presented in this dissertation focuses on the application of the entropy viscosity method to low-Mach single- and two-phase flow equations discretized using a continuous Galerkin finite element method with implicit time integration. The technique has been implemented and tested using the multiphysics simulation environment MOOSE (D Gaston, C Newsman, G Hansen and D Lebrun-Grandie. A parallel computational framework for coupled systems of nonlinear equations. Journal of Nucl. Eng. Design, 239, 1768-1778, 2009).
First, the entropy viscosity method, developed by Guermond et al. (J-L Guermond, R Pasquetti and B Popov. Entropy viscosity method for nonlinear conservation laws. Journal of Comput. Phys., 230, 4248-4267, 2011), is extended to the multi-dimensional Euler equations for both subsonic (very low Mach numbers) and supersonic flows. We show that the current definition of the viscosity coefficients is not adapted to low-Mach flows and we provide a robust alternate definition valid for any Mach number value. The new definitions are derived from a low-Mach asymptotic study, is valid for a wide range of Mach numbers and no longer requires an analytical expression of the entropy function. In addition, the entropy minimum principle is used to derive the viscous regularization terms for Euler equations with variable area for nozzle flow problems and was proved valid for any equation of state with a concave entropy. The new definition of the entropy viscosity method is tested on various 1-D and 2-D numerical benchmarks employing the ideal and the stiffened gas equation of states: flow in a converging-diverging nozzle, Leblanc shock tube, slow moving shock, strong shock for liquid phase, subsonic flows around a 2-D cylinder and over a circular hump, and supersonic flow in a compression corner. Convergence studies are performed using analytical solutions in 1-D and proved the entropy viscosity method to be second-order accurate for smooth solutions.
In a second part, the entropy viscosity method is applied to the seven-equation two-phase flow model. After deriving the dissipative terms using the same procedure as for the multi-D Euler equations, a low-Mach asymptotic study is performed in order to obtain a definition for the viscosity coefficients. Because the seven-equation model is derived by assuming that each phase obeys the Euler equations, the dissipative terms and the definition of the viscosity coefficients are analogous to the ones obtained for the single-phase system of equations. Then, 1-D numerical tests were performed to demonstrate that the entropy viscosity method properly stabilizes the flow simulations based on the seven-equation model.
Another focus of this work was to investigate the impact of source terms (gravity, friction, etc) onto the entropy viscosity method. The theoretical approach adopted here consists of deriving the entropy residual when accounting for the source terms and investigate the sign of the new terms in order to adapt the definition of the viscosity coefficients. Numerical 1-D tests are performed to validate this approach for both single- and two-phase flow models.
In the last part of this dissertation, the entropy viscosity method is applied to the 1-D grey radiation-hydrodynamic equations where the 1-D Euler equations are coupled to a radiation diffusion equation through relaxation terms. The method of manufactured solutions was used to prove second-order accuracy of the numerical stabilization method and also show that the entropy viscosity method yields the correct asymptotic diffusion limit. 1-D tests for inlet Mach number ranging from 1.2 to 50 are presented and show good agreement with semi-analytical solutions.