Higher‐order upwind leapfrog methods for multi‐dimensional acoustic equations

A non‐dissipative and very accurate one‐dimensional upwind leapfrog method was successfully extended to higher‐order and multi‐dimensional acoustic equations. The governing equations in characteristic form and staggered grid were utilized to preserve the accuracy. Fourier analysis was performed to find the accurate scheme for acoustics and the resultant two‐dimensional methods were successfully applied to several classical test cases. Copyright © 2004 John Wiley & Sons, Ltd.     

A robust multi‐grid pressure‐based algorithm for multi‐fluid flow at all speeds

This paper reports on the implementation and testing, within a full non‐linear multi‐grid environment, of a new pressure‐based algorithm for the prediction of multi‐fluid flow at all speeds. The algorithm is part of the mass conservation‐based algorithms (MCBA) group in which the pressure correction equation is derived from overall mass conservation. The performance of the new method is assessed by solving a series of two‐dimensional two‐fluid flow test problems varying from turbulent low Mach number to supersonic flows, and from very low to high fluid density ratios. Solutions are generated for several grid sizes using the single grid (SG), the prolongation grid (PG), and the full non‐linear multi‐grid (FMG) methods. The main outcomes of this study are: (i) a clear demonstration of the ability of the FMG method to tackle the added non‐linearity of multi‐fluid flows, which is manifested through the performance jump observed when using the non‐linear multi‐grid approach as compared to the SG and PG methods; (ii) the extension of the FMG method to predict turbulent multi‐fluid flows at all speeds. The convergence history plots and CPU‐times presented indicate that the FMG method is far more efficient than the PG method and accelerates the convergence rate over the SG method, for the problems solved and the grids used, by a factor reaching a value as high as 15. Copyright © 2003 John Wiley & Sons, Ltd.     

A high‐order upwind control‐volume method based on integrated RBFs for fluid‐flow problems

This paper is concerned with the development of a high‐order upwind conservative discretization method for the simulation of flows of a Newtonian fluid in two dimensions. The fluid‐flow domain is discretized using a Cartesian grid from which non‐overlapping rectangular control volumes are formed. Line integrals arising from the integration of the diffusion and convection terms over control volumes are evaluated using the middle‐point rule. One‐dimensional integrated radial basis function schemes using the multiquadric basis function are employed to represent the variations of the field variables along the grid lines. The convection term is effectively treated using an upwind scheme with the deferred‐correction strategy. Several highly non‐linear test problems governed by the Burgers and the Navier–Stokes equations are simulated, which show that the proposed technique is stable, accurate and converges well. Copyright © 2010 John Wiley & Sons, Ltd.     

A fast nested multi‐grid viscous flow solver for adaptive Cartesian/Quad grids

A nested multi‐grid solution algorithm has been developed for an adaptive Cartesian/Quad grid viscous flow solver. Body‐fitted adaptive Quad (quadrilateral) grids are generated around solid bodies through ‘surface extrusion’. The Quad grids are then overlapped with an adaptive Cartesian grid. Quadtree data structures are employed to record both the Quad and Cartesian grids. The Cartesian grid is generated through recursive sub‐division of a single root, whereas the Quad grids start from multiple roots—a forest of Quadtrees, representing the coarsest possible Quad grids. Cell‐cutting is performed at the Cartesian/Quad grid interface to merge the Cartesian and Quad grids into a single unstructured grid with arbitrary cell topologies (i.e., arbitrary polygons). Because of the hierarchical nature of the data structure, many levels of coarse grids have already been built in. The coarsening of the unstructured grid is based on the Quadtree data structure through reverse tree traversal. Issues arising from grid coarsening are discussed and solutions are developed. The flow solver is based on a cell‐centered finite volume discretization, Roe's flux splitting, a least‐squares linear reconstruction, and a differentiable limiter developed by Venkatakrishnan in a modified form. A local time stepping scheme is used to handle very small cut cells produced in cell‐cutting. Several cycling strategies, such as the saw‐tooth, W‐ and V‐cycles, have been studies. The V‐cycle has been found to be the most efficient. In general, the multi‐grid solution algorithm has been shown to greatly speed up convergence to steady state—by one to two orders. Copyright © 2000 John Wiley & Sons, Ltd.     

A Riemann solver and upwind methods for a two‐phase flow model in non‐conservative form

We present a theoretical solution for the Riemann problem for the five‐equation two‐phase non‐conservative model of Saurel and Abgrall. This solution is then utilized in the construction of upwind non‐conservative methods to solve the general initial‐boundary value problem for the two‐phase flow model in non‐conservative form. The basic upwind scheme constructed is the non‐conservative analogue of the Godunov first‐order upwind method. Second‐order methods in space and time are then constructed via the MUSCL and ADER approaches. The methods are systematically assessed via a series of test problems with theoretical solutions. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

A hybrid FVM–LBM method for single and multi‐fluid compressible flow problems

The lattice Boltzmann method (LBM) has established itself as an alternative approach to solve the fluid flow equations. In this work we combine LBM with the conventional finite volume method (FVM), and propose a non‐iterative hybrid method for the simulation of compressible flows. LBM is used to calculate the inter‐cell face fluxes and FVM is used to calculate the node parameters. The hybrid method is benchmarked for several one‐dimensional and two‐dimensional test cases. The results obtained by the hybrid method show a steeper and more accurate shock profile as compared with the results obtained by the widely used Godunov scheme or by a representative flux vector splitting scheme. Additional features of the proposed scheme are that it can be implemented on a non‐uniform grid, study of multi‐fluid problems is possible, and it is easily extendable to multi‐dimensions. These features have been demonstrated in this work. The proposed method is therefore robust and can possibly be applied to a variety of compressible flow situations. Copyright © 2009 John Wiley & Sons, Ltd.     

A fully conservative high‐order upwind multi‐moment method using moments in both upwind and downwind cells

We propose a fully conservative high‐order upwind multi‐moment method for the conservation equation. The proposed method is based on a third‐order polynomial interpolation function and semi‐Lagrangian formulation and is a variant of the constrained interpolation profile conservative semi‐Lagrangian scheme with third‐order polynomial function method. The third‐order interpolation function is constructed based on three constraints in the upwind cell (two boundary values and a cell average) and a constraint in the downwind cell (a cell center value). The proposed method shows fourth‐order accuracy in a benchmark problem (sine wave propagation). We also propose a less oscillatory formulation of the proposed method. The less oscillatory formulation can minimize numerical oscillations. These methods were validated through scalar transport problems, and compressible flow problems (shock tube and 2D explosion problems). Copyright © 2016 John Wiley & Sons, Ltd.     

Finite difference scheme for the solution of fluid flow problems on non‐staggered grids

A new finite difference methodology is developed for the solution of computational fluid dynamics problems that do not require the use of staggered grid systems. Previous successful and robust non‐staggered methods, which used primitive variables and mass conservation in order to solve the pressure field, either interpolate cell‐face velocities or interpolate the pressure gradients in a special way, usually with an upwind‐bias to avoid the problem of odd–even coupling between the velocity and pressure fields. The new methodology presented does not detail a ‘special interpolation procedure for a primitive variable’, however, it manages to avoid the problem of odd–even coupling. The odd–even coupling is avoided by applying fourth‐order dissipation to the pressure field. It is shown that this approach can be regarded as a modified Rhie and Chow scheme. The method is implemented using a SIMPLE‐type algorithm and is applied to two test problems: laminar flow over a backward‐facing step and laminar flow in a square cavity with a driven lid. Good agreement is obtained between the numerical solutions and the corresponding benchmark solutions. The pressure dissipation term was found to successfully suppress wiggles in the pressure field. Copyright © 2000 John Wiley & Sons, Ltd.     

A pressure‐based unstructured grid method for all‐speed flows

An all‐speed algorithm based on the SIMPLE pressure‐correction scheme and the ‘retarded‐density’ approach has been formulated and implemented within an unstructured grid, finite volume (FV) scheme for both incompressible and compressible flows, the latter involving interaction of shock waves. The collocated storage arrangement for all variables is adopted, and the checkerboard oscillations are eliminated by using a pressure‐weighted interpolation method, similar to that of Rhie and Chow [Numerical study of the turbulent flow past an airfoil with trailing edge separation. AIAA Journal 1983; 21 : 1525]. The solution accuracy is greatly enhanced when a higher‐order convection scheme combined with adaptive mesh refinement (AMR) are used. Copyright © 2000 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.     

Global and local central non‐upwind finite volume schemes for hyperbolic conservation laws in porous media

The focus of this paper is on the development of convective flow approximation schemes that offer the benefits of an upwind formalism without actually upwinding. Development of robust schemes that remove or reduce upwind dependence would represent a significant step forward leading to a fundamental simplification of current methods for flow in porous media and reservoir simulation. Copyright © 2009 John Wiley & Sons, Ltd.     

An implicit numerical algorithm for solving non‐hydrostatic free‐surface flow problems

Details are given of the development of a two‐dimensional vertical numerical model for simulating unsteady free‐surface flows, using a non‐hydrostatic pressure distribution. In this model, the Reynolds equations and the kinematic free‐surface boundary condition are solved simultaneously, so that the water surface elevation can be integrated into the solution and solved for, together with the velocity and pressure fields. An efficient numerical algorithm has been developed, deploying implicit parameters similar to those used in the Crank–Nicholson method, and generating a block tri‐diagonal algebraic system of equations. The model has been applied to simulate a range of unsteady flow problems involving relatively strong vertical accelerations. The results show that the numerical algorithm described is able to produce accurate predictions and is also easy to apply. Copyright © 2001 John Wiley & Sons, Ltd.     

Computation of unsteady viscous incompressible flows in generalized non‐inertial co‐ordinate system using Godunov‐projection method and overlapping meshes

Time‐dependent incompressible Navier–Stokes equations are formulated in generalized non‐inertial co‐ordinate system and numerically solved by using a modified second‐order Godunov‐projection method on a system of overlapped body‐fitted structured grids. The projection method uses a second‐order fractional step scheme in which the momentum equation is solved to obtain the intermediate velocity field which is then projected on to the space of divergence‐free vector fields. The second‐order Godunov method is applied for numerically approximating the non‐linear convection terms in order to provide a robust discretization for simulating flows at high Reynolds number. In order to obtain the pressure field, the pressure Poisson equation is solved. Overlapping grids are used to discretize the flow domain so that the moving‐boundary problem can be solved economically. Numerical results are then presented to demonstrate the performance of this projection method for a variety of unsteady two‐ and three‐dimensional flow problems formulated in the non‐inertial co‐ordinate systems. Copyright © 2002 John Wiley & Sons, Ltd.     

High‐order compact numerical schemes for non‐hydrostatic free surface flows

The development of a numerical scheme for non‐hydrostatic free surface flows is described with the objective of improving the resolution characteristics of existing solution methods. The model uses a high‐order compact finite difference method for spatial discretization on a collocated grid and the standard, explicit, single step, four‐stage, fourth‐order Runge–Kutta method for temporal discretization. The Cartesian coordinate system was used. The model requires the solution of two Poisson equations at each time‐step and tridiagonal matrices for each derivative at each of the four stages in a time‐step. Third‐ and fourth‐order accurate boundaries for the flow variables have been developed including the top non‐hydrostatic pressure boundary. The results demonstrate that numerical dissipation which has been a problem with many similar models that are second‐order accurate is practically eliminated. A high accuracy is obtained for the flow variables including the non‐hydrostatic pressure. The accuracy of the model has been tested in numerical experiments. In all cases where analytical solutions are available, both phase errors and amplitude errors are very small. Copyright © 2006 John Wiley & Sons, Ltd.     

A further work on multi‐phase two‐fluid approach for compressible multi‐phase flows

This paper is to continue our previous work Niu (Int. J. Numer. Meth. Fluids 2001; 36 :351–371) on solving a two‐fluid model for compressible liquid–gas flows using the AUSMDV scheme. We first propose a pressure–velocity‐based diffusion term originally derived from AUSMDV scheme Wada and Liou (SIAM J. Sci. Comput. 1997; 18 (3):633—657) to enhance its robustness. The scheme can be applied to gas and liquid fluids universally. We then employ the stratified flow model Chang and Liou (J. Comput. Physics 2007; 225 :240–873) for spatial discretization. By defining the fluids in different regions and introducing inter‐phasic force on cell boundary, the stratified flow model allows the conservation laws to be applied on each phase, and therefore, it is able to capture fluid discontinuities, such as the fluid interfaces and shock waves, accurately. Several benchmark tests are studied, including the Ransom's Faucet problem, 1D air–water shock tube problems, 2D shock‐water column and 2D shock‐bubble interaction problems. The results indicate that the incorporation of the new dissipation into AUSM⁺‐up scheme and the stratified flow model is simple, accurate and robust enough for the compressible multi‐phase flows. Copyright © 2008 John Wiley & Sons, Ltd.     

A new weighted upwind finite volume element method based on non‐standard covolume for time‐dependent convection–diffusion problems

In modern numerical simulation of problems in energy resources and environmental science, it is important to develop efficient numerical methods for time‐dependent convection–diffusion problems. On the basis of nonstandard covolume grids, we propose a new kind of high‐order upwind finite volume element method for the problems. We first prove the stability and mass conservation in the discrete forms of the scheme. Optimal second‐order error estimate in L₂‐norm in spatial step is then proved strictly. The scheme is effective for avoiding numerical diffusion and nonphysical oscillations and has second‐order accuracy. Numerical experiments are given to verify the performance of the scheme. Copyright © 2013 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.     

Computing non‐Newtonian fluid flow with radial basis function networks

This paper is concerned with the application of radial basis function networks (RBFNs) for solving non‐Newtonian fluid flow problems. Indirect RBFNs, which are based on an integration process, are employed to represent the solution variables; the governing differential equations are discretized by means of point collocation. To enhance numerical stability, stress‐splitting techniques are utilized. The proposed method is verified through the computation of the rectilinear and non‐rectilinear flows in a straight duct and the axisymmetric flow in an undulating tube using Newtonian, power‐law, Criminale–Ericksen–Filbey (CEF) and Oldroyd‐B models. The obtained results are in good agreement with the analytic and benchmark solutions. Copyright © 2005 John Wiley & Sons, Ltd.