A low Mach correction able to deal with low Mach acoustics

International audienceThis article deals with acoustic computations in low Mach number flows with density based solvers. For ensuring a good resolution of the low Mach number base flow, a scheme able to deal with stationary low Mach number flows is necessary. Previously proposed low Mach number fixes are tested with acoustic computations. Numerical results prove that they are not accurate for acoustic computations. The issues raised with acoustic computations with low Mach number fixes are discussed, and a new scheme is developed, in order to be accurate not only for steady low Mach number flows, but also for acoustic computations. Numerical tests show the improvement of the proposed scheme with respect to the state of the art

Rescaling of the Roe Scheme in Low Mach-Number Flow Regions II: Artificial Speed of Sound and Low Mach Number Fix

We look at two simple modifications of the Roe scheme in the incompressible limit, based on different ideas: the Rossow's artificial speed of sound and the Rieper's low Mach number fix. Both schemes modify the eigenspaces of the dissipation matrix. The analysis emphasizes the properties of the dissipation matrix for the Von Neumann stability, the asymptotic behaviour and the solution accuracy in the incompressible limit. Numerical results in the very low-speed limit are discussed for robustness, consistency and accuracy issues of the numerical procedure. Possible occurrence of checkerboard pressure modes, when using a collocated arrangement for velocity components and pressure in the finite-volume scheme, and spurious acoustic modes, is also illustrated for both schemes

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

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal