Comments on ‘Adiabatic shock capturing in perfect gas hypersonic flows’

Some comments are provided on the shock‐capturing techniques and stabilization parameters used in a recent paper (B.S. Kirk, Int. J. Numer. Meth. Fluids 2009; DOI: 10.1002/fld.2195 ) in conjunction with the SUPG formulation of compressible flows. Copyright © 2010 John Wiley & Sons, Ltd.     

An enthalpy‐preserving shock‐capturing term for residual distribution schemes

In this contribution, we investigate strategies to perform shock‐capturing computation of steady hypersonic flow fields by means of residual distribution schemes. The ultimate objective is the computation of flow solutions for which the correct upstream enthalpy value is recovered in the postshock region. To this end, the parallelism existing between the classical Bx scheme and the stabilized finite element techniques is exploited. The simple Lax‐Friedrichs dissipation term is leveraged to build two new residual distribution schemes. Upon testing on both inviscid and viscous steady problems, solutions obtained with one of the two schemes are shown to recover the correct upstream total enthalpy level in the postshock region. This last scheme provides also improved wall pressure and skin friction predictions; heat transfer predictions are, unfortunately, similar to those offered by the Bx scheme. A conjecture for explaining this behavior is exposed.     

Development and validation of a SUPG finite element scheme for the compressible Navier–Stokes equations using a modified inviscid flux discretization

This paper considers the streamline‐upwind Petrov–Galerkin (SUPG) method applied to the unsteady compressible Navier–Stokes equations in conservation‐variable form. The spatial discretization, including a modified approach for interpolating the inviscid flux terms in the SUPG finite element formulation, and the second‐order accurate time discretization are presented. The numerical method is discussed in detail. The performance of the algorithm is then investigated by considering inviscid flow past a circular cylinder. Validation of the finite element formulation via comparisons with experimental data for high‐Mach number perfect gas laminar flows is presented, with a specific focus on comparisons with experimentally measured skin friction and convective heat transfer on a 15° compression ramp. Copyright © 2007 John Wiley & Sons, Ltd.     

Rheology of rarefied gas flow in hypersonic shock and boundary layers

Using the Maxwell method, transfer equations describing molecular gas flows in viscous shock and hypersonic boundary layers are obtained. It is shown that, in contrast to the Navier-Stokes approximation, the kinetic model proposed makes it possible correctly to describe hypersonic flow around bodies under conditions of strong nonequilibrium of the internal and translational degrees of freedom of the gas particles.     

A high‐resolution scheme for compressible multicomponent flows with shock waves

A simple methodology for a high‐resolution scheme to be applied to compressible multicomponent flows with shock waves is investigated. The method is intended for use with direct numerical simulation or large eddy simulation of compressible multicomponent flows. The method dynamically adds non‐linear artificial diffusivity locally in space to capture different types of discontinuities such as a shock wave, contact surface or material interface while a high‐order compact differencing scheme resolves a broad range of scales in flows. The method is successfully applied to several one‐dimensional and two‐dimensional compressible multicomponent flow problems with shock waves. The results are in good agreement with experiments and earlier computations qualitatively and quantitatively. The method captures unsteady shock and material discontinuities without significant spurious oscillations if initial start‐up errors are properly avoided. Comparisons between the present numerical scheme and high‐order weighted essentially non‐oscillatory (WENO) schemes illustrate the advantage of the present method for resolving a broad range of scales of turbulence while capturing shock waves and material interfaces. Also the present method is expected to require less computational cost than popular high‐order upwind‐biased schemes such as WENO schemes. The mass conservation for each species is satisfied due to the strong conservation form of governing equations employed in the method. Copyright © 2010 John Wiley & Sons, Ltd.     

A note on compressive limiting for two‐material flows

In this short note we describe a simple extension to the multi‐material shock‐capturing algorithm presented in (J. Comput. Phys. 2007; 223 :262–297) that can be used to maintain sharp material interfaces. The method takes the form of an artificial compression which is designed so that the material indicator jumps across only a few cells but which does not excite physical instabilities in the flow. The advantages of the approach include its simplicity and flexibility in that it provides a parameter that effectively determines the captured interface thickness. Copyright © 2009 John Wiley & Sons, Ltd.     

The capturing of free surfaces in incompressible multi‐fluid flows

By treating it as a contact discontinuity in the density field, a free surface between two immiscible fluids can be automatically ‘captured’ by the enforcement of conservation laws. A surface‐capturing method of this kind requires no special tracking or fitting treatment for the free surface, thereby offering the advantage of algorithm simplicity over the surface‐tracking or the surface‐fitting method. A surface‐capturing method based on a new multi‐fluid incompressible Navier–Stokes formulation is developed. It is applied to a variety of free‐surface flows, including the Rayleigh–Taylor instability problem, the ship waves around a Wigley hull and a model bubble‐rising problem to demonstrate the validity and versatility of the present method. Copyright © 2000 John Wiley & Sons, Ltd.     

A high‐order finite difference method for numerical simulations of supersonic turbulent flows

This paper describes a new class of three‐dimensional finite difference schemes for high‐speed turbulent flows in complex geometries based on the high‐order monotonicity‐preserving (MP) method. Simulations conducted for various 1D, 2D, and 3D problems indicate that the new high‐order MP schemes can preserve sharp changes in the flow variables without spurious oscillations and are able to capture the turbulence at the smallest computed scales. Our results also indicate that the MP method has less numerical dissipation and faster grid convergence than the weighted essentially non‐oscillatory method. However, both of these methods are computationally more demanding than the COMP method and are only used for the inviscid fluxes. To reduce the computational cost for reacting flows, the scalar equations are solved by the COMP method, which is shown to yield similar results to those obtained by the MP in supersonic turbulent flows with strong shock waves. Copyright © 2011 John Wiley & Sons, Ltd.     

A lattice boltzmann model for compressible perfect gas

A new lattice Boltzmann model for compressible perfect gas is proposed. The numerical example shows that it can be used to simulate shock wave and contact discontinuity. The results are comparable with those obtained by traditional methods. The ratio of specific heats γ may be chosen according to the requirement of problems. The project supported by the National Natural Science Foundation of China     

A study of time integration schemes for the numerical modelling of free surface flows

A semi‐discrete finite element methodology for the modelling of transient free surface flows in the context of Eulerian interface capturing is proposed. The focus of this study is put on the choice of an appropriate time integration strategy for the accurate modelling of the dynamics of free surfaces and of interfacial physics. It is composed of an adaptive time integration scheme for the Navier–Stokes equations, and of the implicit midpoint rule for the transport equation of the Eulerian marker variable. The adaptive scheme allows the automatic determination of a time‐step size that follows the physics of the problem under study, which facilitates the accurate modelling of stiff free surface flows. It is shown that the implicit midpoint rule reduces mass loss for each fluid. Various free surface flow problems are studied to verify and validate the proposed time integration strategy. Copyright © 2005 John Wiley & Sons, Ltd.     

The robustness issue on multigrid schemes applied to the Navier–Stokes equations for laminar and turbulent,incompressible and compressible flows

The paper's leitmotiv is condensed in one word: robustness. This is a real hindrance for the successful implementation of any multigrid scheme for solving the Navier–Stokes set of equations. In this paper, many hints are given to improve this issue. Instead of looking for the best possible speed‐up rate for a particular set of problems, at a given regime and in a given condition, the authors propose some ideas pursuing reasonable speed‐up rates in any situation. In a previous paper, the authors presented a multigrid method for solving the incompressible turbulent RANS equations, with particular care in the robustness and flexibility of the solution scheme. Here, these concepts are further developed and extended to compressible laminar and turbulent flows. This goal is achieved by introducing a non‐linear multigrid scheme for compressible laminar (NS equations) and turbulent flow (RANS equations), taking benefit of a convenient master–slave implementation strategy. Copyright © 2004 John Wiley & Sons, Ltd.     

Controlling the disturbances in a hypersonic flat-plate shock layer by unsteady action from the surface

The results of a wind-tunnel experiment on the joint action of periodic acoustic fast-mode disturbances of the outer flow and disturbances generated at the leading edge of a plate on the hypersonic (M_∞ = 21) viscous shock layer on the plate are presented. The possibility of positively controlling the intensity of density fluctuations in the plate shock layer by means of disturbances introduced from the leading edge is shown. Direct numerical simulation of the suppression (enhancement) of disturbances under the simultaneous action on the shock layer of the two-dimensional fast-mode acoustic waves in the outer flow and the source of two-dimensional suction/injection disturbances near the leading edge of the plate is performed under the experimental conditions. The experimental and calculated results are shown to be in good agreement.     

Evaluation of some high‐order shock capturing schemes for direct numerical simulation of unsteady two‐dimensional free flows

The present study addresses the capability of a large set of shock‐capturing schemes to recover the basic interactions between acoustic, vorticity and entropy in a direct numerical simulation (DNS) framework. The basic dispersive and dissipative errors are first evaluated by considering the advection of a Taylor vortex in a uniform flow. Two transonic cases are also considered. The first one consists of the interaction between a temperature spot and a weak shock. This test emphasizes the capability of the schemes to recover the production of vorticity through the baroclinic process. The second one consists of the interaction of a Taylor vortex with a weak shock, corresponding to the framework of the linear theory of Ribner. The main process in play here is the production of an acoustic wave. The results obtained by using essentially non‐oscillatory (ENO), total variation diminishing (TVD), compact‐TVD and MUSCL schemes are compared with those obtained by means of a sixth‐order accurate Hermitian scheme, considered as reference. The results are as follows; the ENO schemes agree pretty well with the reference scheme. The second‐order accurate Upwind‐TVD scheme exhibits a strong numerical diffusion, while the MUSCL scheme behavior is very sensitive to the value on the parameter β in the limiter function minmod. The compact‐TVD schemes do not yield improvement over the standard TVD schemes. Copyright © 2000 John Wiley & Sons, Ltd.     

A stabilized finite element formulation for high‐speed inviscid compressible flows using finite calculus

This paper aims at the development of a new stabilization formulation based on the finite calculus (FIC) scheme for solving the Euler equations using the Galerkin FEM on unstructured triangular grids. The FIC method is based on expressing the balance of fluxes in a space–time domain of finite size. It is used to prevent the creation of instabilities typically present in numerical solutions due to the high convective terms and sharp gradients. Two stabilization terms, respectively called streamline term and transverse term, are added via the FIC formulation to the original conservative equations in the space–time domain. An explicit fourth‐order Runge–Kutta scheme is implemented to advance the solution in time. The presented numerical test examples for inviscid flows prove the ability of the proposed stabilization technique for providing appropriate solutions especially near shock waves. Although the derived methodology delivers precise results with a nearly coarse mesh, a mesh refinement technique is coupled to the solution process for obtaining a suitable mesh particularly in the high‐gradient zones. Copyright © 2014 John Wiley & Sons, Ltd.     

A finite element model for free surface flows on fixed meshes

In this paper, we present a finite element model for free surface flows on fixed meshes. The main novelty of the approach, compared with typical fixed mesh finite element models for such flows, is that we take advantage of the particularities of free surface flow, instead of considering it a particular case of two‐phase flow. The fact that a given free surface implies a known boundary condition on the interface, allows us to solve the Navier–Stokes equations on the fluid domain uncoupled from the solution on the rest of the finite element mesh. This, together with the use of enhanced integration allows us to model low Froude number flows accurately, something that is not possible with typical two‐phase flow models applied to free surface flow. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical study of the oscillations induced by shock/shock interaction in hypersonic double-wedge flows

In this paper, the shock pattern oscillations induced by shock/shock interactions over double-wedge geometries in hypersonic flows were studied numerically by solving 2D inviscid Euler equations for a multi-species system. Laminar viscous effects were considered in some cases. Temperature-dependent thermodynamic properties were employed in the state and energy equations for consideration of the distinct change of the thermodynamic state. It was shown that the oscillation results in high-frequency fluctuations of heating and pressure loads over wedge surfaces. In a case with a relatively lower free-stream Mach number, the shock/shock interaction structure maintains a seven-shock configuration during the entire oscillation process. On the other hand, the oscillation is accompanied by a transition between a six-shock configuration (regular interaction) and a seven-shock configuration (Mach interaction) in a case with a higher free-stream Mach number. Numerical results also indicate that the critical wedge angle for the transition from a steady to an oscillation solution is higher compared to the corresponding value in earlier numerical research in which the perfect diatomic gas model was used.      

A fully implicit scheme for simulating ionized gas flows using the gas dynamics electrodynamics coupled system

The flow of ionized gases under the influence of electromagnetic fields is governed by the coupled system of the compressible flow equations and the Maxwell equations. In this system, coupling of the flow with the electromagnetic field is obtained through nonlinear and stiff source terms, which may cause difficulties with the numerical solution of the coupled system. The discontinuous Galerkin finite element method is used for the numerical solution of this system. For the magnetic field vector, discontinuous Galerkin discretization is performed using a divergence‐free vector base for the magnetic field to preserve zero divergence in the element and retain the implicit constraint of a divergence‐free magnetic field vector down to very low level both globally and locally. To circumvent difficulties resulting from the presence of the stiff source terms, implicit time marching is used for the fully coupled system to avoid wrong wave shapes and propagation speeds that are obtained when the coupling source terms are lagged in time or by using splitting iterative schemes. Numerical solutions for benchmark problems computed on collocated meshes for the flow and electromagnetic field variables with this fully coupled monolithic approach showed good agreement with other numerical solutions and exact results. Copyright © 2014 John Wiley & Sons, Ltd.     

Accurate particle splitting for smoothed particle hydrodynamics in shallow water with shock capturing

The solution for the shallow water equations using smoothed particle hydrodynamics is attractive, being a mesh‐free, automatically adaptive method without special treatment for wet–dry interfaces. However, the relatively new method is limited by the variable kernel size or smoothing length being inversely proportional to water depth causing poor resolution at small depths. Boundary conditions at solid walls have also not been well resolved. To solve the resolution problem in small depths, a particle splitting procedure was developed (conveniently into seven particles), which conserves mass and momentum by varying the smoothing length, velocity and acceleration of each refined particle. This improves predictions in the shallowest depths where the error associated with splitting is reduced by one order of magnitude in comparison to other published works. To provide good shock capturing behaviour, particle interactions are treated as a Riemann problem with Monotone Upstream‐centred Scheme for Conservation Laws (MUSCL) reconstruction providing stability. For solid boundaries, the recent modified virtual boundary particle method was developed further to enable the zeroth moment to be accurately conserved where the smoothing length of particles is changing rapidly during particle splitting. The resulting method is applied to the one‐dimensional and the two‐dimensional axisymmetric wet‐bed dam break problems showing close agreement with analytical solutions, demonstrating the need for particle splitting. To demonstrate wetting and drying in a more complex case, the scheme is applied to oscillating water in a two‐dimensional parabolic basin and produces good agreement with the analytical solution. The method is finally applied to the European Concerted Action on DAm break Modelling dam‐break test case representative of realistic conditions and good predictions are made of experimental measurements with a 40% reduction in the computational time when particle splitting is employed. The overall method has thus become quite sophisticated but its generality and versatility will be attractive for various shallow water problems. Copyright © 2011 John Wiley & Sons, Ltd.     

Similarity solutions for the flow behind an exponential shock in a non-ideal gas

Similarity solutions for the flow of a non-ideal gas behind a strong exponential shock driven out by a piston (cylindrical or spherical) moving with time according to an exponential law are obtained. Similarity solutions exist only when the surrounding medium is of constant density. Solutions are obtained, in both the cases, when the flow between the shock and the piston is isothermal or adiabatic. It is found that the assumption of zero temperature gradient brings a profound change in the density distribution as compare to that of the adiabatic case. Effects of the non-idealness of the gas on the flow-field between the shock and the piston are investigated. The variations of density-ratio across the shock and the location of the piston with the parameter of non-idealness of the gas are also obtained.