Contributions to the Development of Entropy-Stable Schemes for Compressible Flows

Entropy-Stable (ES) schemes have gathered considerable attention over the last decade, especially in the context of under-resolved simulations of compressible turbulent flows, where achieving both high-order accuracy and robustness is difficult. ES schemes provide stability in a nonlinear and integral sense: the total entropy of the discrete solution can be made non-decreasing, in agreement with the second principle of thermodynamics. Additionally, the amount of entropy produced by the scheme is known and can be modified, making room for analysis and improvements. This thesis delves into some of the challenges currently limiting their use in practice. The current state of the art solves the compressible Navier-Stokes equations for a single-component perfect gas in chemical and thermal equilibrium. This model is inappropriate in aerospace engineering applications such as hypersonics and combustion, which typically involve chemically reacting gas mixtures far from equilibrium. As a first step towards enabling their use for these applications, we formulated ES schemes for the multicomponent compressible Euler equations. Special care had to be taken as we found out that the theoretical foundations of ES schemes begin to crumble in the limit of vanishing partial densities. The realization that ES schemes can only go as far as their theory led us to review some of it. A fundamental result supporting the development of limiting strategies for high-order methods is the minimum entropy principle for the compressible Euler equations. It states that the specific entropy of the physically relevant weak solution does not decrease. We proved that the same result holds for the specific entropy of the gas mixture in the multicomponent case. While entropy-stability is a valuable property, it does not imply a well-behaved solution. One must recall that the second principle is a prescription on the correct behavior of a system at the global level only. To better understand how ES schemes may or may not improve the quality of the numerical solution, we revisited two classical problems encountered in the development of shock-capturing techniques. First, we studied the receding flow problem, which is a simple setup used to study the anomalous temperature rise, termed "overheating", typically observed in shock reflection and shock interaction calculations. Previous studies showed that the anomaly can be cured if conservation of entropy is enforced, but at the considerable price of total energy conservation. Entropy-Conservative (EC) schemes, a particular instance of ES schemes, can achieve both simultaneously and therefore appeared as a potential solution. We showed that while the overheating is correlated to entropy production, entropy conservation does not necessarily prevent it. Second, we studied the behavior of ES schemes in the low Mach number regime, where shock-capturing schemes are known to suffer from severe accuracy degradation issues. A classic remedy to this problem is the flux-preconditioning technique, which consists in modifying artificial dissipation terms to enforce consistent low Mach behavior. We showed that ES schemes suffer from the same issues and that the flux-preconditioning technique can improve their behavior without interfering with entropy-stability. Furthermore, we demonstrated analytically that these issues stem from an acoustic entropy production field which scales improperly with the Mach number, generating spatial fluctuations that are inconsistent with the equations. An important outgrowth of this effort is the discovery that skew-symmetric dissipation operators can alter the way entropy is produced or conserved locally.PHDAerospace EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/155304/1/gouasmia_1.pd

超臨界炭化水素におけるマルチフィジックス熱流動の数値解法

Tohoku University博士（情報科学）thesi

Direct Numerical Simulation of Interfacial Flows: Implicit Sharp-Interface Method (I-SIM)

In recent work (Nourgaliev, Liou, Theofanous, JCP in press) we demonstrated that numerical simulations of interfacial flows in the presence of strong shear must be cast in dynamically sharp terms (sharp interface treatment or SIM), and that moreover they must meet stringent resolution requirements (i.e., resolving the critical layer). The present work is an outgrowth of that work aiming to overcome consequent limitations on the temporal treatment, which become still more severe in the presence of phase change. The key is to avoid operator splitting between interface motion, fluid convection, viscous/heat diffusion and reactions; instead treating all these non-linear operators fully-coupled within a Newton iteration scheme. To this end, the SIM’s cut-cell meshing is combined with the high-orderaccurate implicit Runge-Kutta and the “recovery” Discontinuous Galerkin methods along with a Jacobian-free, Krylov subspace iteration algorithm and its physics-based preconditioning. In particular, the interfacial geometry (i.e., marker’s positions and volumes of cut cells) is a part of the Newton-Krylov solution vector, so that the interface dynamics and fluid motions are fully-(non-linearly)-coupled. We show that our method is: (a) robust (L-stable) and efficient, allowing to step over stability time steps at will while maintaining high-(up to the 5th)-order temporal accuracy; (b) fully conservative, even near multimaterial contacts, without any adverse consequences (pressure/velocity oscillations); and (c) highorder-accurate in spatial discretization (demonstrated here up to the 12th-order for smoothin-the-bulk-fluid flows), capturing interfacial jumps sharply, within one cell. Performance is illustrated with a variety of test problems, including low-Mach-number “manufactured” solutions, shock dynamics/tracking with slow dynamic time scales, and multi-fluid, highspeed shock-tube problems. We briefly discuss preconditioning, and we introduce two physics-based preconditioners – “Block-Diagonal” and “Internal energy-Pressure-Velocity Partially Decoupled”, demonstrating the ability to efficiently solve all-speed flows with strong effects from viscous dissipation and heat conduction