A numerical scheme for strong blast wave driven by explosion

After the detonation of a solid high explosive, the material has extremely high pressure keeping the solid density and expands rapidly driving strong shock wave. In order to simulate this blast wave, a stable and accurate numerical scheme is required due to large density and pressure changes in time and space. The compressible fluid equations are solved by a fractional step procedure which consists of the advection phase and non‐advection phase. The former employs the Rational function CIP scheme in order to preserve monotone signals, and the latter is solved by interpolated differential operator scheme for achieving the accurate calculation. The procedure is categorized into the fractionally stepped semi‐Lagrangian. The accuracy of our scheme is confirmed by checking the one‐dimensional plane shock tube problem with 10³ times initial density and pressure jump in comparison with the analytic solution. The Sedov–Taylor blast wave problem is also examined in the two‐dimensional cylindrical coordinate in order to check the spherical symmetry and the convergence rates. Two‐ and three‐dimensional simulations for the blast waves from the explosion in the underground magazine are carried out. It is found that the numerical results show quantitatively good agreement with the experimental data. Copyright © 2006 John Wiley & Sons, Ltd.     

A robust low‐dissipation AUSM‐family scheme for numerical shock stability on unstructured grids

The simple low‐dissipation advection upwind splitting method (SLAU) scheme is a parameter‐free, low‐dissipation upwind scheme that has been applied in a wide range of aerodynamic numerical simulations. In spite of its successful applications, the SLAU scheme could be showing shock instabilities on unstructured grids, as many other contact resolved upwind schemes. Therefore, a hybrid upwind flux scheme is devised for improving the shock stability of SLAU scheme, without compromising on accuracy and low Mach number performance. Numerical flux function of the hybrid scheme is written in a general form, in which only the scalar dissipation term is different from that of the SLAU scheme. The hybrid dissipation term is defined by using a differentiable multidimensional‐shock‐detection pressure weight function, and the dissipation term of SLAU scheme is combined with that of the Van Leer scheme. Furthermore, the hybrid dissipation term is only applied for the solution of momentum fluxes in numerical flux function. Based on the numerical test results, the hybrid scheme is deemed to be a successful improvement on the shock stability of SLAU scheme, without compromising on the efficiency and accuracy. Copyright © 2016 John Wiley & Sons, Ltd.     

High‐order shock capturing schemes for turbulence calculations

This work investigates high‐order central compact methods for simulating turbulent supersonic flows that include shock waves. Several different types of previously proposed characteristic filters, including total variation diminishing, monotone upstream‐centered scheme for conservation laws, and weighted essentially non‐oscillatory filters, are investigated in this study. Similar to the traditional shock capturing schemes, these filters can eliminate the numerical instability caused by large gradients in flow fields, but they also improve efficiency compared with classical shock‐capturing schemes. Adding the nonlinear dissipation part of a classical shock‐capturing scheme to a central scheme makes the method suitable for incorporation into any existing central‐based high‐order subsonic code. The amount of numerical dissipation to add is sensed by means of the artificial compression method switch. In order to improve the performance of the characteristic filters, we propose a hybrid approach to minimize the dissipation added by the characteristic filter. Through several numerical experiments (including a shock/density wave interaction, a shock/vortex interaction, and a shock/mixing layer interaction) we show that our hybrid approach works better than the original method, and can be used for future turbulent flow simulations that include shocks. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical assessment of a shock‐detecting sensor for low dissipative high‐order simulation of shock–vortex interactions

Considering the importance of high‐order schemes implementation for the simulation of shock‐containing turbulent flows, the present work involves the assessment of a shock‐detecting sensor for filtering of high‐order compact finite‐difference schemes for simulation of this type of flows. To accomplish this, a sensor that controls the amount of numerical dissipation is applied to a sixth‐order compact scheme as well as a fourth‐order two‐register Runge–Kutta method for numerical simulation of various cases including inviscid and viscous shock–vortex and shock–mixing‐layer interactions. Detailed study is performed to investigate the performance of the sensor, that is, the effect of control parameters employed in the sensor are investigated in the long‐time integration. In addition, the effects of nonlinear weighting factors controlling the value of the second‐order and high‐order filters in fine and coarse non‐uniform grids are investigated. The results indicate the accuracy of the nonlinear filter along with the promising performance of the shock‐detecting sensor, which would pave the way for future simulations of turbulent flows containing shocks. Copyright © 2014 John Wiley & Sons, Ltd.     

A simple high‐resolution advection scheme

A simple, robust, mass‐conserving numerical scheme for solving the linear advection equation is described. The scheme can estimate peak solution values accurately even in regions where spatial gradients are high. Such situations present a severe challenge to classical numerical algorithms. Attention is restricted to the case of pure advection in one and two dimensions since this is where past numerical problems have arisen. The authors' scheme is of the Godunov type and is second‐order in space and time. The required cell interface fluxes are obtained by MUSCL interpolation and the exact solution of a degenerate Riemann problem. Second‐order accuracy in time is achieved via a Runge–Kutta predictor–corrector sequence. The scheme is explicit and expressed in finite volume form for ease of implementation on a boundary‐conforming grid. Benchmark test problems in one and two dimensions are used to illustrate the high‐spatial accuracy of the method and its applicability to non‐uniform grids. Copyright © 2007 John Wiley & Sons, Ltd.     

Evaluation of an energy relaxation method for the simulation of unsteady,viscous, real gas flows

This study investigates a new energy relaxation method designed to capture the dynamics of unsteady, viscous, real gas flows governed by the compressible Navier–Stokes equations. We focus on real gas models accounting for inelastic molecular collisions and yielding temperature‐dependent heat capacities. The relaxed Navier–Stokes equations are discretized using a mixed finite volume/finite element method and a high‐order time integration scheme. The accuracy of the energy relaxation method is investigated on three test problems of increasing complexity: the advection of a periodic set of vortices, the interaction of a temperature spot with a weak shock, and finally, the interaction of a reflected shock with its trailing boundary layer in a shock tube. In all cases, the method is validated against benchmark solutions and the numerical errors resulting from both discretization and energy relaxation are assessed independently. Copyright © 2004 John Wiley & Sons, Ltd.     

Entropy dissipation scheme and minimums‐increase‐and‐maximums‐decrease slope limiter

In this paper, we continue to study the entropy dissipation scheme developed in former. We start with a numerical study of the scheme without the entropy dissipation term on the linear advection equation, which shows that the scheme is stable and numerical dissipation and numerical dispersion free for smooth solutions. However, the numerical results for discontinuous solutions show nonlinear instabilities near jump discontinuities. This is because the scheme enforces two related conservation properties in the computation. With this study, we design a so‐called ‘minimums‐increase‐and‐maximums‐decrease’ slope limiter in the reconstruction step of the scheme and delete the entropy dissipation in the linear fields and reduce the entropy dissipation terms in the nonlinear fields. Numerical experiments show improvements of the designed scheme compared with the results presented in former. However, the minimums‐increase‐and‐maximums‐decrease limiter is still not perfect yet, and better slope limiters are still sought. Copyright © 2011 John Wiley & Sons, Ltd.     

A high‐order PIC method for advection‐dominated flow with application to shallow water waves

A hybrid Eulerian‐Lagrangian particle‐in‐cell–type numerical method is developed for the solution of advection‐dominated flow problems. Particular attention is given over to the high‐order transfer of flow properties from the particles to the grid. For smooth flows, the method presented is of formal high‐order accuracy in space. The method is applied to solve the nonlinear shallow water equations resulting in a new, and novel, shock capturing shallow water solver. The approach is able to simulate complex shallow water flows, which can contain an arbitrary number of discontinuities. Both trivial and nontrivial bottom topography is considered, and it is shown that the new scheme is inherently well balanced, exactly satisfying the ‐property. The scheme is verified against several one‐dimensional benchmark shallow water problems. These include cases that involve transcritical flow regimes, shock waves, and nontrivial bathymetry. In all the test cases presented, very good results are obtained.     

Optimized sixth‐order monotonicity‐preserving scheme by nonlinear spectral analysis

In this paper, sixth‐order monotonicity‐preserving optimized scheme (OMP6) for the numerical solution of conservation laws is developed on the basis of the dispersion and dissipation optimization and monotonicity‐preserving technique. The nonlinear spectral analysis method is developed and is used for the purpose of minimizing the dispersion errors and controlling the dissipation errors. The new scheme (OMP6) is simple in expression and is easy for use in CFD codes. The suitability and accuracy of this new scheme have been tested through a set of one‐dimensional, two‐dimensional, and three‐dimensional tests, including the one‐dimensional Shu–Osher problem, the two‐dimensional double Mach reflection, and the Rayleigh–Taylor instability problem, and the three‐dimensional direct numerical simulation of decaying compressible isotropic turbulence. All numerical tests show that the new scheme has robust shock capturing capability and high resolution for the small‐scale waves due to fewer numerical dispersion and dissipation errors. Moreover, the new scheme has higher computational efficiency than the well‐used WENO schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

Simulation of multiple shock–shock interference using implicit anti‐diffusive WENO schemes

Accurate computations of two‐dimensional turbulent hypersonic shock–shock interactions that arise when single and dual shocks impinge on the bow shock in front of a cylinder are presented. The simulation methods used are a class of lower–upper symmetric‐Gauss–Seidel implicit anti‐diffusive weighted essentially non‐oscillatory (WENO) schemes for solving the compressible Navier–Stokes equations with Spalart–Allmaras one‐equation turbulence model. A numerical flux of WENO scheme with anti‐diffusive flux correction is adopted, which consists of first‐order and high‐order fluxes and allows for a more flexible choice of first‐order dissipative methods. Experimental flow fields of type IV shock–shock interactions with single and dual incident shocks by Wieting are computed. By using the WENO scheme with anti‐diffusive flux corrections, the present solution indicates that good accuracy is maintained and contact discontinuities are sharpened markedly as compared with the original WENO schemes on the same meshes. Computed surface pressure distribution and heat transfer rate are also compared with experimental data and other computational results and good agreement is found. Copyright © 2009 John Wiley & Sons, Ltd.     

A horizontally curvilinear non‐hydrostatic model for simulating nonlinear wave motion in curved boundaries

A horizontally curvilinear non‐hydrostatic free surface model that embeds the second‐order projection method, the so‐called θ scheme, in fractional time stepping is developed to simulate nonlinear wave motion in curved boundaries. The model solves the unsteady, Navier–Stokes equations in a three‐dimensional curvilinear domain by incorporating the kinematic free surface boundary condition with a top‐layer boundary condition, which has been developed to improve the numerical accuracy and efficiency of the non‐hydrostatic model in the standard staggered grid layout. The second‐order Adams–Bashforth scheme with the third‐order spatial upwind method is implemented in discretizing advection terms. Numerical accuracy in terms of nonlinear phase speed and amplitude is verified against the nonlinear Stokes wave theory over varying wave steepness in a two‐dimensional numerical wave tank. The model is then applied to investigate the nonlinear wave characteristics in the presence of dispersion caused by reflection and diffraction in a semicircular channel. The model results agree quantitatively with superimposed analytical solutions. Finally, the model is applied to simulate nonlinear wave run‐ups caused by wave‐body interaction around a bottom‐mounted cylinder. The numerical results exhibit good agreement with experimental data and the second‐order diffraction theory. Overall, it is shown that the developed model, with only three vertical layers, is capable of accurately simulating nonlinear waves interacting within curved boundaries. Copyright © 2011 John Wiley & Sons, Ltd.     

Solution to Euler equations by high‐resolution upwind compact scheme based on flux splitting

A high‐resolution upwind compact method based on flux splitting is developed for solving the compressive Euler equations. The convective flux terms are discretized by using the modified advection upstream splitting method (AUSM). The developed scheme is used to compute the one‐dimensional Burgers equation and four different example problems of supersonic compressible flows, respectively. The results show that the high‐resolution upwind compact scheme based on modified AUSM+ flux splitting can capture shock wave and other discontinuities, obtain higher resolution and restrain numerical oscillation. Copyright © 2007 John Wiley & Sons, Ltd.     

A dissipation‐free numerical method to solve one‐dimensional hyperbolic flow equations

In this paper, a numerical method to capture the shock wave propagation in 1‐dimensional fluid flow problems with 0 numerical dissipation is presented. Instead of using a traditional discrete grid, the new numerical method is built on a range‐discrete grid, which is obtained by a direct subdivision of values around the shock area. The range discrete grid consists of 2 types: continuous points and shock points. Numerical solution is achieved by tracking characteristics and shocks for the movements of continuous and shock points, respectively. Shocks can be generated or eliminated when triggering entropy conditions in a marking step. The method is conservative and total variation diminishing. We apply this new method to several examples, including solving Burgers equation for aerodynamics, Buckley‐Leverett equation for fractional flow in porous media, and the classical traffic flow. The solutions were verified against analytical solutions under simple conditions. Comparisons with several other traditional methods showed that the new method achieves a higher accuracy in capturing the shock while using much less grid number. The new method can serve as a fast tool to assess the shock wave propagation in various flow problems with good accuracy.     

High resolution flux‐difference‐splitting scheme on adaptive grid for open‐channel flows

The effectiveness and usefulness of further enhancing the shock resolution of a second‐order accurate scheme for open‐channel flows by using an adaptive grid is investigated. The flux‐difference‐splitting (FDS) scheme based on the Lax–Wendroff numerical flux is implemented on a fixed as well as on a self‐adjusting grid for this purpose. The grid‐adjusting procedure, developed by Harten and Hyman, adjusts the grid by averaging the local characteristic velocities with respect to the signal amplitude in such a way that a shock always lies on a mesh point. This enables a scheme capable of perfectly resolving a stationary shock to capture a shock that moves from mesh point to mesh point. The Roe's approximate Jacobian is used for conservation and consistency, while theoretically sound treatment for satisfying entropy inequality conditions ensures physically realistic solutions. Details about inclusion of source terms, often left out of analyses for the homogeneous part of governing equations, are also explained. The numerical results for some exacting problems are compared with analytical as well as experimental results for examining improvements in resolution of discontinuities by the adaptive grid. Copyright © 2001 John Wiley & Sons, Ltd.     

An adaptive multiresolution semi‐intrusive scheme for UQ in compressible fluid problems

This paper deals with the introduction of a multiresolution strategy into the semi‐intrusive scheme, recently introduced by the authors, aiming to propagate uncertainties in unsteady compressible fluid applications. The mathematical framework of the multiresolution setting is presented for the cell‐average case, and the coupling with the semi‐intrusive scheme is described from both the theoretical and algorithmic point‐of‐view. Some reference test cases are performed to demonstrate the convergence properties and the efficiency of the overall scheme: the linear advection problem for both smooth and discontinuous initial conditions, the inviscid Burgers equation, and an uncertain shock tube problem obtained by modifying the well‐known Sod shock problem. For all the cases, the convergence curves are computed with respect to semi‐analytical (exact) solutions. In the case of the shock tube problem, an original technique to obtain a reference highly‐accurate numerical stochastic solution has also been developed. Copyright © 2015 John Wiley & Sons, Ltd.     

Time‐linearized time‐harmonic 3‐D Navier–Stokes shock‐capturing schemes

In the present paper, a numerical method for the computation of time‐harmonic flows, using the time‐linearized compressible Reynolds‐averaged Navier–Stokes equations is developed and validated. The method is based on the linearization of the discretized nonlinear equations. The convective fluxes are discretized using an O(Δx) MUSCL scheme with van Leer flux‐vector‐splitting. Unsteady perturbations of the turbulent stresses are linearized using a frozen‐turbulence‐Reynolds‐number hypothesis, to approximate eddy‐viscosity perturbations. The resulting linear system is solved using a pseudo‐time‐marching implicit ADI‐AF (alternating‐directions‐implicit approximate‐factorization) procedure with local pseudo‐time‐steps, corresponding to a matrix‐successive‐underrelaxation procedure. The stability issues associated with the pseudo‐time‐marching solution of the time‐linearized Navier–Stokes equations are discussed. Comparison of computations with measurements and with time‐nonlinear computations for 3‐D shock‐wave oscillation in a square duct, for various back‐pressure fluctuation frequencies (180, 80, 20 and 10 Hz), assesses the shock‐capturing capability of the time‐linearized scheme. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

Augmented upwind numerical schemes for the groundwater transport advection‐dispersion equation with local operators

When solute transport is advection‐dominated, the advection‐dispersion equation approximates to a hyperbolic‐type partial differential equation, and finite difference and finite element numerical approximation methods become prone to artificial oscillations. The upwind scheme serves to correct these responses to produce a more realistic solution. The upwind scheme is reviewed and then applied to the advection‐dispersion equation with local operators for the first‐order upwinding numerical approximation scheme. The traditional explicit and implicit schemes, as well as the Crank‐Nicolson scheme, are developed and analyzed for numerical stability to form a comparison base. Two new numerical approximation schemes are then proposed, namely, upwind–Crank‐Nicolson scheme, where only for the advection term is applied, and weighted upwind‐downwind scheme. These newly developed schemes are analyzed for numerical stability and compared to the traditional schemes. It was found that an upwind–Crank‐Nicolson scheme is appropriate if the Crank‐Nicolson scheme is only applied to the advection term of the advection‐dispersion equation. Furthermore, the proposed explicit weighted upwind‐downwind finite difference numerical scheme is an improvement on the traditional explicit first‐order upwind scheme, whereas the implicit weighted first‐order upwind‐downwind finite difference numerical scheme is stable under all assumptions when the appropriate weighting factor (θ) is assigned.     

Robustness of the hybrid DRP‐WENO scheme for shock flow computations

This paper describes a new variant of hybrid scheme that is constructed by a wave‐capturing scheme and a nonoscillatory scheme for flow computations in the presence of shocks. The improved fifth‐order upwind weighted essentially nonoscillatory scheme is chosen to be conjugated with the seven‐point dispersion‐relation‐preserving scheme by means of an adaptive switch function of grid‐point type. The new hybrid scheme can achieve a better resolution than the hybrid scheme which is based on the classical weighted essentially scheme. Ami Harten's multiresolution analysis algorithm is applied to density field for detecting discontinuities and setting point values of the switch function adaptively. Moreover, the tenth‐order central filter is applied in smooth part of the flow field for damping dispersion errors. This scheme can promote overall computational efficiency and yield oscillation‐free results in shock flows. The resolution properties and robustness of the new hybrid scheme are tested in both 1D and 2D linear and nonlinear cases. It performs well for computing flow problems with rich structures of weak/strong shocks and large/small vortices, such as the shock‐boundary layer interaction problem in a shock tube, which illustrates that it is very robust and accurate for direct numerical simulation of gas‐dynamics flows. Copyright © 2011 John Wiley & Sons, Ltd.     

Mass,momentum and total excess energy transported by a weak planeN wave

A weakly nonlinear plane acoustic wave is emitted into an ideal gas of semi-infinite extent from an infinite plate by its sinusoidal motion of single period. The wave develops into anN wave in the far field, as long as the energy dissipation is negligible everywhere except for discontinuous shock fronts. The third-order effects at shock fronts are evaluated, i.e., the generation of reflected acoustic wave as a result of the interaction of shock and expansion wave and the production of entropy by the energy dissipation at shock fronts. Consideration of these effects enables one to estimate the whole mass, momentum and total excess energy (sum of the kinetic energy and excess of internal energy over an initial undisturbed value) transported by theN wave to the accuracy of third order of wave amplitude. It is shown that the mass and total excess energy transported by theN wave increase and the momentum decreases to asymptotic limits as the wave propagates. The result shows good agreement with a numerical result obtained by solving the Euler equations with a high-resolution TVD upwind scheme.