A stable second‐order mass‐weighted upwind scheme for unstructured meshes

In this paper, an original second‐order upwind scheme for convection terms is described and implemented in the context of a Control‐Volume Finite‐Element Method (CVFEM). The proposed scheme is a second‐order extension of the first‐order MAss‐Weighted upwind (MAW) scheme proposed by Saabas and Baliga (Numer. Heat Transfer 1994; 26B :381–407). The proposed second‐order scheme inherits the well‐known stability characteristics of the MAW scheme, but exhibits less artificial viscosity and ensures much higher accuracy. Consequently, and in contrast with nearly all second‐order upwind schemes available in the literature, the proposed second‐order MAW scheme does not need limiters. Some test cases including two pure convection problems, the driven cavity and steady and unsteady flows over a circular cylinder, have been undertaken successfully to validate the new scheme. The verification tests show that the proposed scheme exhibits a low level of artificial viscosity in the pure convection problems; exhibits second‐order accuracy for the driven cavity; gives accurate reattachment lengths for low‐Reynolds steady flow over a circular cylinder; and gives constant‐amplitude vortex shedding for the case of high‐Reynolds unsteady flow over a circular cylinder. Copyright © 2005 John Wiley & Sons, Ltd.     

Convergence issues in using high‐resolution schemes and lower–upper symmetric Gauss–Seidel method for steady shock‐induced combustion problems

This paper reports numerical convergence study for simulations of steady shock‐induced combustion problems with high‐resolution shock‐capturing schemes. Five typical schemes are used: the Roe flux‐based monotone upstream‐centered scheme for conservation laws (MUSCL) and weighted essentially non‐oscillatory (WENO) schemes, the Lax–Friedrichs splitting‐based non‐oscillatory no‐free parameter dissipative (NND) and WENO schemes, and the Harten–Yee upwind total variation diminishing (TVD) scheme. These schemes are implemented with the finite volume discretization on structured quadrilateral meshes in dimension‐by‐dimension way and the lower–upper symmetric Gauss–Seidel (LU–SGS) relaxation method for solving the axisymmetric multispecies reactive Navier–Stokes equations. Comparison of iterative convergence between different schemes has been made using supersonic combustion flows around a spherical projectile with Mach numbers M = 3.55 and 6.46 and a ram accelerator with M = 6.7. These test cases were regarded as steady combustion problems in literature. Calculations on gradually refined meshes show that the second‐order NND, MUSCL, and TVD schemes can converge well to steady states from coarse through fine meshes for M = 3.55 case in which shock and combustion fronts are separate, whereas the (nominally) fifth‐order WENO schemes can only converge to some residual level. More interestingly, the numerical results show that all the schemes do not converge to steady‐state solutions for M = 6.46 in the spherical projectile and M = 6.7 in the ram accelerator cases on fine meshes although they all converge on coarser meshes or on fine meshes without chemical reactions. The result is based on the particular preconditioner of LU–SGS scheme. Possible reasons for the nonconvergence in reactive flow simulation are discussed.Copyright © 2012 John Wiley & Sons, Ltd.     

Applications of a higher‐order bounded numerical scheme to turbulent flows

This paper applies the higher‐order bounded numerical scheme Weighted Average Coefficients Ensuring Boundedness (WACEB) to simulate two‐ and three‐dimensional turbulent flows. In the scheme, a weighted average formulation is used for interpolating the variables at cell faces and the weighted average coefficients are determined from a normalized variable formulation and total variation diminishing (TVD) constraints to ensure the boundedness of the solution. The scheme is applied to two turbulent flow problems: (1) two‐dimensional turbulent flow around a blunt plate; and (2) three‐dimensional turbulent flow inside a mildly curved U‐bend. In the present study, turbulence is evaluated by using a low‐Reynolds number version of the k–ω model. For the flow simulation, the QUICK scheme is applied to the momentum equations while either the WACEB scheme (Method 1) or the UPWIND scheme (Method 2) is used for the turbulence equations. The present study shows that the WACEB scheme has at least second‐order accuracy while ensuring boundedness of the solutions. The present numerical study for a pure convection problem shows that the ‘TVD’ slope ranges from 2 to 4. For the turbulent recirculating flow, two different mixed procedures (Method 1 and Method 2) produce a substantial difference for the mean velocities as well as for the turbulence kinetic energy. Method 1 predicts better results than Method 2 does, comparing the analytical solution and the experimental data. For the turbulent flow inside the mildly curved U‐bend, although the predictions of velocity distributions with two procedures are very close, a noticeable difference of turbulence kinetic energy is exhibited. It is noticed that the discrepancy exists between numerical results and the experimental data. The reason is the limit of the two‐equation turbulence model to such complex turbulent flows with extra strain‐rates. Copyright © 2001 John Wiley & Sons, Ltd.     

A Roe scheme for a compressible six‐equation two‐fluid model

We derive a partially analytical Roe scheme with wave limiters for the compressible six‐equation two‐fluid model. Specifically, we derive the Roe averages for the relevant variables. First, the fluxes are split into convective and pressure parts. Then, independent Roe conditions are stated for these two parts. These conditions are successively reduced while defining acceptable Roe averages. For the convective part, all the averages are analytical. For the pressure part, most of the averages are analytical, whereas the remaining averages are dependent on the thermodynamic equation of state. This gives a large flexibility to the scheme with respect to the choice of equation of state. Furthermore, this model contains nonconservative terms. They are a challenge to handle right, and it is not the object of this paper to discuss this issue. However, the Roe averages presented in this paper are fully independent from how those terms are handled, which makes this framework compatible with any treatment of nonconservative terms. Finally, we point out that the eigenspace of this model may collapse, making the Roe scheme inapplicable. This is called resonance. We propose a fix to handle this particular case. Numerical tests show that the scheme performs well. Copyright © 2012 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.     

A novel weighting switch function for uniformly high‐order hybrid shock‐capturing schemes

Hybrid schemes are very efficient for complex compressible flow simulation. However, for most existing hybrid schemes in literature, empirical problem‐dependent parameters are always needed to detect shock waves and hence greatly decrease the robustness and accuracy of the hybrid scheme. In this paper, based on the nonlinear weights of the weighted essentially non‐oscillatory (WENO) scheme, a novel weighting switch function is proposed. This function approaches 1 with high‐order accuracy in smooth regions and 0 near discontinuities. Then, with the new weighting switch function, a seventh‐order hybrid compact‐reconstruction WENO scheme (HCCS) is developed. The new hybrid scheme uses the same stencil as the fifth‐order WENO scheme, and it has seventh‐order accuracy in smooth regions even at critical points. Numerical tests are presented to demonstrate the accuracy and robustness of both the switch function and HCCS. Comparisons also reveal that HCCS has lower dissipation and less computational cost than the seventh‐order WENO scheme. Copyright © 2016 John Wiley & Sons, Ltd.     

A ‘well‐balanced’ finite volume scheme for blood flow simulation

We are interested in simulating blood flow in arteries with a one‐dimensional model. Thanks to recent developments in the analysis of hyperbolic system of conservation laws (in the Saint‐Venant shallow water equations context) we will perform a simple finite volume scheme. We focus on conservation properties of this scheme, which were not previously considered. To emphasize the necessity of this scheme, we present how a too simple numerical scheme may induce spurious flows when the basic static shape of the radius changes. On the contrary, the proposed scheme is ‘well‐balanced’: it preserves equilibria of Q = 0. Then examples of analytical or linearized solutions with and without viscous damping are presented to validate the calculations. The influence of abrupt change of basic radius is emphasized in the case of an aneurism. Copyright © 2012 John Wiley & Sons, Ltd.     

Improved third‐order weighted essentially nonoscillatory schemes with new smoothness indicators

In this article, we present two improved third‐order weighted essentially nonoscillatory (WENO) schemes for recovering their design‐order near first‐order critical points. The schemes are constructed in the framework of third‐order WENO‐Z scheme. Two new global smoothness indicators, τ_L3 and τ_L4, are devised by a nonlinear combination of local smoothness indicators (IS_k) and reference values (IS_G) based on Lagrangian interpolation polynomial. The performances of the proposed schemes are evaluated on several numerical tests governed by one‐dimensional linear advection equation or one‐ and two‐dimensional Euler equations. Numerical results indicate that the presented schemes provide less dissipation and higher resolution than the original WENO3‐JS and subsequent WENO3‐N scheme.     

A fourth‐order compact finite difference scheme for the steady stream function–vorticity formulation of the Navier–Stokes/Boussinesq equations

A fourth‐order compact finite difference scheme on the nine‐point 2D stencil is formulated for solving the steady‐state Navier–Stokes/Boussinesq equations for two‐dimensional, incompressible fluid flow and heat transfer using the stream function–vorticity formulation. The main feature of the new fourth‐order compact scheme is that it allows point‐successive overrelaxation (SOR) or point‐successive underrelaxation iteration for all Rayleigh numbers Ra of physical interest and all Prandtl numbers Pr attempted. Numerical solutions are obtained for the model problem of natural convection in a square cavity with benchmark solutions and compared with some of the accurate results available in the literature. Copyright © 2003 John Wiley & Sons, Ltd.     

The multi‐stage centred‐scheme approach applied to a drift‐flux two‐phase flow model

For two‐phase flow models, upwind schemes are most often difficult do derive, and expensive to use. Centred schemes, on the other hand, are simple, but more dissipative. The recently proposed multi‐stage (MUSTA ) method is aimed at coming close to the accuracy of upwind schemes while retaining the simplicity of centred schemes. So far, the MUSTA approach has been shown to work well for the Euler equations of inviscid, compressible single‐phase flow. In this work, we explore the MUSTA scheme for a more complex system of equations: the drift‐flux model, which describes one‐dimensional two‐phase flow where the motions of the phases are strongly coupled. As the number of stages is increased, the results of the MUSTA scheme approach those of the Roe method. The good results of the MUSTA scheme are dependent on the use of a large‐enough local grid. Hence, the main benefit of the MUSTA scheme is its simplicity, rather than CPU ‐time savings. Copyright © 2006 John Wiley & Sons, Ltd.     

Depth‐integrated,non‐hydrostatic model with grid nesting for tsunami generation,propagation, and run‐up

Tsunamis generated by earthquakes involve physical processes of different temporal and spatial scales that extend across the ocean to the shore. This paper presents a shock‐capturing dispersive wave model in the spherical coordinate system for basin‐wide evolution and coastal run‐up of tsunamis and discusses the implementation of a two‐way grid‐nesting scheme to describe the wave dynamics at resolution compatible to the physical processes. The depth‐integrated model describes dispersive waves through the non‐hydrostatic pressure and vertical velocity, which also account for tsunami generation from dynamic seafloor deformation. The semi‐implicit, finite difference model captures flow discontinuities associated with bores or hydraulic jumps through the momentum‐conserved advection scheme with an upwind flux approximation. The two‐way grid‐nesting scheme utilizes the Dirichlet condition of the non‐hydrostatic pressure and both the horizontal velocity and surface elevation at the inter‐grid boundary to ensure propagation of dispersive waves and discontinuities across computational grids of different resolution. The inter‐grid boundary can adapt to bathymetric features to model nearshore wave transformation processes at optimal resolution and computational efficiency. A coordinate transformation enables application of the model to small geographic regions or laboratory experiments with a Cartesian grid. A depth‐dependent Gaussian function smoothes localized bottom features in relation to the water depth while retaining the bathymetry important for modeling of tsunami transformation and run‐up. Numerical experiments of solitary wave propagation and N‐wave run‐up verify the implementation of the grid‐nesting scheme. The 2009 Samoa Tsunami provides a case study to confirm the validity and effectiveness of the modeling approach for tsunami research and impact assessment. Copyright © 2010 John Wiley & Sons, Ltd.     

On a higher‐order bounded discretization scheme

This paper presents a new higher‐order bounded scheme, weighted‐average coefficient ensuring boundedness (WACEB), for approximating the convective fluxes in solving transport equations with the finite volume difference method (FVDM). The weighted‐average formulation is used for interpolating the variables at cell faces and the weighted‐average coefficient is determined from normalized variable formulation and total variation diminishing (TVD) constraints to ensure the boundedness of solution. The new scheme is tested by solving three problems: (1) a pure convection of a box‐shaped step profile in an oblique velocity field, (2) a sudden expansion of an oblique velocity field in a cavity, and (3) a laminar flow over a fence. The results obtained by the present WACEB are compared with the UPWIND and the QUICK schemes and it is shown that this scheme has at least second‐order accuracy, while ensuring boundedness of solutions. Moreover, it is demonstrated that this scheme produces results that better agree with the experimental data in comparison with other schemes. Copyright © 2000 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.     

Depth‐integrated free‐surface flow with a two‐layer non‐hydrostatic formulation

This paper presents the derivation of a depth‐integrated wave propagation and runup model from a system of governing equations for two‐layer non‐hydrostatic flows. The governing equations are transformed into an equivalent, depth‐integrated system, which separately describes the flux‐dominated and dispersion‐dominated processes. The depth‐integrated system reproduces the linear dispersion relation within a 5 error for water depth parameter up to kd = 11, while allowing direct implementation of a momentum conservation scheme to model wave breaking and a moving‐waterline technique for runup calculation. A staggered finite‐difference scheme discretizes the governing equations in the horizontal dimension and the Keller box scheme reconstructs the non‐hydrostatic terms in the vertical direction. An semi‐implicit scheme integrates the depth‐integrated flow in time with the non‐hydrostatic pressure determined from a Poisson‐type equation. The model is verified with solitary wave propagation in a channel of uniform depth and validated with previous laboratory experiments for wave transformation over a submerged bar, a plane beach, and fringing reefs. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical solution of the two‐phase expansion of a metastable flashing liquid jet using the dispersion‐controlled dissipative scheme

This paper presents a study of the stationary phenomenon of superheated or metastable liquid jets, flashing into a two‐dimensional axisymmetric domain, while in the two‐phase region. In general, the phenomenon starts off when a high‐pressure, high‐temperature liquid jet emerges from a small nozzle or orifice expanding into a low‐pressure chamber, below its saturation pressure taken at the injection temperature. As the process evolves, crossing the saturation curve, one observes that the fluid remains in the liquid phase reaching a superheated condition. Then, the liquid undergoes an abrupt phase change by means of an oblique evaporation wave. Across this phase change the superheated liquid becomes a two‐phase high‐speed mixture in various directions, expanding to supersonic velocities. In order to reach the downstream pressure, the supersonic fluid continues to expand, crossing a complex bow shock wave. The balance equations that govern the phenomenon are mass conservation, momentum conservation, and energy conservation, plus an equation‐of‐state for the substance. A false‐transient model is implemented using the shock capturing scheme: dispersion‐controlled dissipative (DCD), which was used to calculate the flow conditions as the steady‐state condition is reached. Numerical results with computational code DCD‐2D v1 have been analyzed. Copyright © 2009 John Wiley & Sons, Ltd.     

Higher‐order bounded differencing schemes for compressible and incompressible flows

In recent years, three higher‐order (HO) bounded differencing schemes, namely AVLSMART, CUBISTA and HOAB that were derived by adopting the normalized variable formulation (NVF), have been proposed. In this paper, a comparative study is performed on these schemes to assess their numerical accuracy, computational cost as well as iterative convergence property. All the schemes are formulated on the basis of a new dual‐formulation in order to facilitate their implementations on unstructured meshes. Based on the proposed dual‐formulation, the net effective blending factor (NEBF) of a high‐resolution (HR) scheme can now be measured and its relevance on the accuracy and computational cost of a HR scheme is revealed on three test problems: (1) advection of a scalar step‐profile; (2) 2D transonic flow past a circular arc bump; and (3) 3D lid‐driven incompressible cavity flow. Both density‐based and pressure‐based methods are used for the computations of compressible and incompressible flow, respectively. Computed results show that all the schemes produce solutions which are nearly as accurate as the third‐order QUICK scheme; however, without the unphysical oscillations which are commonly inherited from the HO linear differencing scheme. Generally, it is shown that at higher value of NEBF, a HR scheme can attain better accuracy at the expense of computational cost. Copyright © 2006 John Wiley & Sons, Ltd.     

A pressure correction‐volume of fluid method for simulation of two‐phase flows

A pressure correction method coupled with the volume of fluid (VOF) method is developed to simulate two‐phase flows. A volume fraction function is introduced in the VOF method and is governed by an advection equation. A modified monotone upwind scheme for a conservation law (modified MUSCL) is used to solve the solution of the advection equation. To keep the initial sharpness of an interface, a slope modification scheme is introduced. The continuum surface tension (CST) model is used to calculate the surface tension force. Three schemes, central‐upwind, Parker–Youngs, and mixed schemes, are introduced to compute the interface normal vector and the gradient of the volume fraction function. Moreover, a height function technique is applied to compute the local curvature of the interface. Several basic test problems are performed to check the order of accuracy of the present numerical schemes for computing the interface normal vector and the gradient of the volume fraction function. Three physical problems, two‐dimensional broken dam problem, static drop, and spurious currents, and three‐dimensional rising bubble, are performed to demonstrate the efficiency and accuracy of the pressure correction method. Copyright © 2010 John Wiley & Sons, Ltd.     

A new weighting method for improving the WENO‐Z scheme

The calculation of the weight of each substencil is very important for a weighted essentially nonoscillatory (WENO) scheme to obtain high‐order accuracy in smooth regions and keep the essentially nonoscillatory property near discontinuities. The weighting function introduced in the WENO‐Z scheme provides a straightforward method to analyze the accuracy order in smooth regions. In this paper, we construct a new sixth‐order global smoothness indicator (GSI‐6) and a function about GSI‐6 and the local smoothness indicators (IS_k) to calculate the weights. The analysis and numerical results show that, with the new weights, the scheme satisfies the sufficient condition for the fifth‐order convergence in smooth regions even at critical points. Meanwhile, it can also maintain low dissipation for discontinuous solutions due to relative large weights assigned to discontinuous substencils.     

A depth‐integrated non‐hydrostatic finite element model for wave propagation

This paper presents a two‐dimensional finite element model for simulating dynamic propagation of weakly dispersive waves. Shallow water equations including extra non‐hydrostatic pressure terms and a depth‐integrated vertical momentum equation are solved with linear distributions assumed in the vertical direction for the non‐hydrostatic pressure and the vertical velocity. The model is developed based on the platform of a finite element model, CCHE2D. A physically bounded upwind scheme for the advection term discretization is developed, and the quasi second‐order differential operators of this scheme result in no oscillation and little numerical diffusion. The depth‐integrated non‐hydrostatic wave model is solved semi‐implicitly: the provisional flow velocity is first implicitly solved using the shallow water equations; the non‐hydrostatic pressure, which is implicitly obtained by ensuring a divergence‐free velocity field, is used to correct the provisional velocity, and finally the depth‐integrated continuity equation is explicitly solved to satisfy global mass conservation. The developed wave model is verified by an analytical solution and validated by laboratory experiments, and the computed results show that the wave model can properly handle linear and nonlinear dispersive waves, wave shoaling, diffraction, refraction and focusing. Copyright © 2013 John Wiley & Sons, Ltd.     

Three‐dimensional modeling of non‐hydrostatic free‐surface flows on unstructured grids

A three‐dimensional numerical model is presented for the simulation of unsteady non‐hydrostatic shallow water flows on unstructured grids using the finite volume method. The free surface variations are modeled by a characteristics‐based scheme, which simulates sub‐critical and super‐critical flows. Three‐dimensional velocity components are considered in a collocated arrangement with a σ‐coordinate system. A special treatment of the pressure term is developed to avoid the water surface oscillations. Convective and diffusive terms are approximated explicitly, and an implicit discretization is used for the pressure term to ensure exact mass conservation. The unstructured grid in the horizontal direction and the σ coordinate in the vertical direction facilitate the use of the model in complicated geometries. Solution of the non‐hydrostatic equations enables the model to simulate short‐period waves and vertically circulating flows. Copyright © 2015 John Wiley & Sons, Ltd.