Improvements to compressible Euler methods for low‐Mach number flows

In the present study improvements to numerical algorithms for the solution of the compressible Euler equations at low Mach numbers are investigated. To solve flow problems for a wide range of Mach numbers, from the incompressible limit to supersonic speeds, preconditioning techniques are frequently employed. On the other hand, one can achieve the same aim by using a suitably modified acoustic damping method. The solution algorithm presently under consideration is based on Roe's approximate Riemann solver [Roe PL. Approximate Riemann solvers, parameter vectors and difference schemes. Journal of Computational Physics 1981; 43 : 357–372] for non‐structured meshes. The numerical flux functions are modified by using Turkel's preconditioning technique proposed by Viozat [Implicit upwind schemes for low Mach number compressible flows. INRIA, Rapport de Recherche No. 3084, January 1997] for compressible Euler equations and by using a modified acoustic damping of the stabilization term proposed in the present study. These methods allow the compressible Euler equations at low‐Mach number flows to be solved, and they are consistent in time. The efficiency and accuracy of the proposed modifications have been assessed by comparison with experimental data and other numerical results in the literature. Copyright © 2000 John Wiley & Sons, Ltd.     

A time‐accurate scheme for the calculations of unsteady reactive flows at low Mach number

The artificial compressibility method is extended to the case of unsteady turbulent reacting flows at low Mach number. The resulting scheme is applied to the calculation of a propagating one‐dimensional (1D) planar turbulent flame with a realistic heat release parameter. An eddy break‐up‐like approach, with a conventional gradient expression for the turbulent fluxes, is retained to model this reacting turbulent flow. A quenched form of the mean reaction rate is used to ensure the existence of a steady regime of propagation, for which the present results are compared with those obtained by a steady analysis of the mean flame brush structure, with excellent agreement. A sensitivity analysis of the convergence rate to the values of the artificial compressibility factor and the pseudo‐time is carried out. It is shown that a reduced artificial compressibility factor of 5–10, combined with a pseudo‐Courant number of ≈1000, represents a good compromise to optimize the convergence rate. Copyright © 1999 John Wiley & Sons, Ltd.     

SOLUTION ALGORITHM FOR LOW-PECLET-NUMBER REACTING FLOW

An enhanced solution strategy based on the SIMPLER algorithm is presented for low-Peclet-number mass transport calculations with applications in low-pressure material processing. The accurate solution of highly diffusive flows requires boundary conditions that preserve specified chemical species mass fluxes. The implementation of such boundary conditions in the standard SIMPLER solution procedure leads to degraded convergence that scales with the Peclet number. Modifications to both the non-linear and linear parts of the solution algorithm remove the slow convergence problem. In particular, the linearized species transport equations must be implicitly coupled to the boundary condition equations and the combined system must be solved exactly at each non-linear iteration. The pressure correction boundary conditions are reformulated to ensure that continuity is preserved in each finite volume at each iteration. The boundary condition scaling problem is demonstrated with a simple linear model problem. The enhanced solution strategy is implemented in a baseline computer code that is used to solve the multicomponent Navier–Stokes equations on a generalized, multiple-block grid system. Accelerated convergence rates are demonstrated for several material-processing example problems. © 1997 John Wiley & Sons, Ltd.     

Transition between regular and Mach reflection of shock waves: new numerical and experimental results

New numerical and experimental results on the transition between regular and Mach reflections of steady shock waves are presented. The influence of flow three-dimensionality on transition between steady regular and Mach reflection has been studied in detail both numerically and experimentally. Characteristic features of 3D shock wave configuration, such as peripheral Mach reflection, non-monotonous Mach stem variation in transverse direction, the existence of combined Mach-regular-peripheral Mach shock wave configuration, have been found in the numerical simulations. The application of laser sheet imaging technique in streamwise direction allowed us to confirm all the details of shock wave configuration in the experiments. Close agreement of the numerical and experimental data on Mach stem heights is shown. Received 23 November 2000 / Accepted 25 April 2001     

马赫-曾德尔干涉仪消光比和单通道光开关

对马赫-曾德尔(M—Z)干涉仪的滤光消光比特性进行了研究，提出了一种在其一端加上光纤全反镜的方法，大大改善了马赫一曾德尔干涉仪的输出谱特性，由此提出并设计了一种单通道光开关，实验得到了约37．4dB的开关动态范围。     

光纤光栅水听器技术实验研究

采用光纤布喇格光栅作为基本传感器件, 设计制作了光纤光栅水听器. 利用匹配光纤光栅解调技术, 实现了高灵敏度、高速度的动态信号测量, 并进行了相关的水下声音测量实验. 实验系统可测量水声信号频率范围为10 Hz~3 kHz, 测量动态范围为60 dB, 其结果具有很高的线性度.     

A hybrid pressure–density‐based algorithm for the Euler equations at all Mach number regimes

In the present work, we propose a reformulation of the fluxes and interpolation calculations in the PISO method, a well‐known pressure‐correction solver. This new reformulation introduces the AUSM⁺ ? up flux definition as a replacement for the standard Rhie and Chow method of obtaining fluxes and central interpolation of pressure at the control volume faces. This algorithm tries to compatibilize the good efficiency of a pressure based method for low Mach number applications with the advantages of AUSM⁺ ? up at high Mach number flows. The algorithm is carefully validated using exact solutions. Results for subsonic, transonic and supersonic axisymmetric flows in a nozzle are presented and compared with exact analytical solutions. Further, we also present and discuss subsonic, transonic and supersonic results for the well known bump test‐case. The code is also benchmarked against a very tough test‐case for the supersonic and hypersonic flow over a cylinder. Copyright © 2011 John Wiley & Sons, Ltd.     

An L2‐stable approximation of the Navier–Stokes convection operator for low‐order non‐conforming finite elements

We develop in this paper a discretization for the convection term in variable density unstationary Navier–Stokes equations, which applies to low‐order non‐conforming finite element approximations (the so‐called Crouzeix–Raviart or Rannacher–Turek elements). This discretization is built by a finite volume technique based on a dual mesh. It is shown to enjoy an L² stability property, which may be seen as a discrete counterpart of the kinetic energy conservation identity. In addition, numerical experiments confirm the robustness and the accuracy of this approximation; in particular, in L² norm, second‐order space convergence for the velocity and first‐order space convergence for the pressure are observed. Copyright © 2010 John Wiley & Sons, Ltd.     

Mach configuration in pseudo-stationary compressible flow

This paper is devoted to studying the local structure of Mach reflection, which occurs in the problem of the shock front hitting a ramp. The compressible flow is described by the full unsteady Euler system of gas dynamics. Because of the special geometry, the motion of the fluid can be described by self-similar coordinates, so that the unsteady flow becomes a pseudo-stationary flow in this coordinate system. When the slope of the ramp is less than a critical value, the Mach reflection occurs. The wave configuration in Mach reflection is composed of three shock fronts and a slip line bearing contact discontinuity. The local existence of a flow field with such a configuration under some assumptions is proved in this paper. Our result confirms the reasonableness of the corresponding physical observations and numerical computations in Mach reflection.

In order to prove the result, we formulate the problem to a free boundary value problem of a pseudo-stationary Euler system. In this problem two unknown shock fronts are the free boundary, and the slip line is also an unknown curve inside the flow field. The proof contains some crucial ingredients. The slip line will be transformed to a fixed straight line by a generalized Lagrange transformation. The whole free boundary value problem will be decomposed to a fixed boundary value problem of the Euler system and a problem to updating the location of the shock front. The Euler system in the subsonic region is an elliptic-hyperbolic composite system, which will be decoupled to the elliptic part and the hyperbolic part at the level of principal parts. Then some sophisticated estimates and a suitable iterative scheme are established. The proof leads to the existence and stability of the local structure of Mach reflection.

     

On the RR→MR transition of asymmetric shock waves in steady flows

In this paper, RR→MR transition of asymmetric shock waves has been theoretically studied. The transition can occur between the sonic-point and maximum-deflection criteria due to the the effects of expansion fans which are inherent flow structures. Comparison shows a better agreement among experiments and the analytical results. Some discrepancies reported in previous studies among experiments and theory have also been explained based on the threshold for RR→MR transition.