Developing a hybrid flux function suitable for hypersonic flow simulation with high‐order methods

In this paper, we develop a new hybrid Euler flux function based on Roe's flux difference scheme, which is free from shock instability and still preserves the accuracy and efficiency of Roe's flux scheme. For computational cost, only 5% extra CPU time is required compared with Roe's FDS. In hypersonic flow simulation with high‐order methods, the hybrid flux function would automatically switch to the Rusanov flux function near shock waves to improve the robustness, and in smooth regions, Roe's FDS would be recovered so that the advantages of high‐order methods can be maintained. Multidimensional dissipation is introduced to eliminate the adverse effects caused by flux function switching and further enhance the robustness of shock‐capturing, especially when the shock waves are not aligned with grids. A series of tests shows that this new hybrid flux function with a high‐order weighted compact nonlinear scheme is not only robust for shock‐capturing but also accurate for hypersonic heat transfer prediction. Copyright © 2015 John Wiley & Sons, Ltd.     

An unstructured finite‐volume method for coupled models of suspended sediment and bed load transport in shallow‐water flows

The aim of this work is to develop a well‐balanced finite‐volume method for the accurate numerical solution of the equations governing suspended sediment and bed load transport in two‐dimensional shallow‐water flows. The modelling system consists of three coupled model components: (i) the shallow‐water equations for the hydrodynamical model; (ii) a transport equation for the dispersion of suspended sediments; and (iii) an Exner equation for the morphodynamics. These coupled models form a hyperbolic system of conservation laws with source terms. The proposed finite‐volume method consists of a predictor stage for the discretization of gradient terms and a corrector stage for the treatment of source terms. The gradient fluxes are discretized using a modified Roe's scheme using the sign of the Jacobian matrix in the coupled system. A well‐balanced discretization is used for the treatment of source terms. In this paper, we also employ an adaptive procedure in the finite‐volume method by monitoring the concentration of suspended sediments in the computational domain during its transport process. The method uses unstructured meshes and incorporates upwinded numerical fluxes and slope limiters to provide sharp resolution of steep sediment concentrations and bed load gradients that may form in the approximate solutions. Details are given on the implementation of the method, and numerical results are presented for two idealized test cases, which demonstrate the accuracy and robustness of the method and its applicability in predicting dam‐break flows over erodible sediment beds. The method is also applied to a sediment transport problem in the Nador lagoon.Copyright © 2013 John Wiley & Sons, Ltd.     

A robust high‐resolution finite volume scheme for the simulation of long waves over complex domains

The propagation, runup and rundown of long surface waves are numerically investigated, initially in one dimension, using a well‐balanced high‐resolution finite volume scheme. A conservative form of the nonlinear shallow water equations with source terms is solved numerically using a high‐resolution Godunov‐type explicit scheme coupled with Roe's approximate Riemann solver. The scheme is also extended to handle two‐dimensional complex domains. The numerical difficulties related to the presence of the topography source terms in the model equations along with the appearance of the wet/dry fronts are properly treated and extended. The resulting numerical model accurately describes breaking waves as bores or hydraulic jumps and conserves volume across flow discontinuities. Numerical results show very good agreement with previously presented analytical or asymptotic solutions as well as with experimental benchmark data. Copyright © 2007 John Wiley & Sons, Ltd.     

Solution of the hyperbolic mild‐slope equation using the finite volume method

A finite volume solver for the 2D depth‐integrated harmonic hyperbolic formulation of the mild‐slope equation for wave propagation is presented and discussed. The solver is implemented on unstructured triangular meshes and the solution methodology is based upon a Godunov‐type second‐order finite volume scheme, whereby the numerical fluxes are computed using Roe's flux function. The eigensystem of the mild‐slope equations is derived and used for the construction of Roe's matrix. A formulation that updates the unknown variables in time implicitly is presented, which produces a more accurate and reliable scheme than hitherto available. Boundary conditions for different types of boundaries are also derived. The agreement of the computed results with analytical results for a range of wave propagation/transformation problems is very good, and the model is found to be virtually paraxiality‐free. Copyright © 2003 John Wiley & Sons, Ltd.     

Two‐dimensional dam break flow simulation

Numerical modelling of shallow water flow in two dimensions is presented in this work with the results obtained in dam break tests. Free surface flow in channels can be described mathematically by the shallow‐water system of equations. These equations have been discretized using an approach based on unstructured Delaunay triangles and applied to the simulation of two‐dimensional dam break flows. A cell centred finite volume method based on Roe's approximate Riemann solver across the edges of the cells is presented and the results are compared for first‐ and second‐order accuracy. Special treatment of the friction term has been adopted and will be described. The scheme is capable of handling complex flow domains as shown in the simulation corresponding to the test cases proposed, i.e. that of a dam break wave propagating into a 45° bend channel (UCL) and in a channel with a constriction (LNEC‐IST). Comparisons of experimental and numerical results are shown. Copyright © 2000 John Wiley & Sons, Ltd.     

Numerical methods for conjunctive two‐dimensional surface and three‐dimensional sub‐surface flows

Sophisticated catchment runoff problems necessitate conjunctive modeling of overland flow and sub‐surface flow. In this paper, finite difference numerical methods are studied for simulation of catchment runoff of two‐dimensional surface flow interacting with three‐dimensional unsaturated and saturated sub‐surface flows. The equations representing the flows are mathematically classified as a type of heat diffusion equation. Therefore, two‐ and three‐dimensional numerical methods for heat diffusion equations were investigated for applications to the surface and sub‐surface flow sub‐models in terms of accuracy, stability, and calculation time. The methods are the purely explicit method, Saul'yev's methods, the alternating direction explicit (ADE) methods, and the alternating direction implicit (ADI) methods. The methods are first examined on surface and sub‐surface flows separately; subsequently, 12 selected combinations of methods were investigated for modeling the conjunctive flows. Saul'yev's downstream (S‐d) method was found to be the preferred method for two‐dimensional surface flow modeling, whereas the ADE method of Barakat and Clark is a less accurate, stable alternative. For the three‐dimensional sub‐surface flow model, the ADE method of Larkin (ADE‐L) and Brian's ADI method are unconditionally stable and more accurate than the other methods. The calculations of the conjunctive models utilizing the S‐d surface flow sub‐model give excellent results and confirm the expectation that the errors of the surface and sub‐surface sub‐models interact. The surface sub‐model dominates the accuracy and stability of the conjunctive model, whereas the sub‐surface sub‐model dominates the calculation time, suggesting the desirability of using a smaller time increment for the surface sub‐model. Copyright © 2000 John Wiley & Sons, Ltd.     

Simulation of non‐hydrostatic,density‐stratified flow in irregular domains

A numerical model has been developed for simulating density‐stratified flow in domains with irregular but simple topography. The model was designed for simulating strong interactions between internal gravity waves and topography, e.g. exchange flows in contracting channels, tidally or convectively driven flow over two‐dimensional sills or waves propagating onto a shoaling bed. The model is based on the non‐hydrostatic, Boussinesq equations of motion for a continuously stratified fluid in a rotating frame, subject to user‐configurable boundary conditions. An orthogonal boundary fitting co‐ordinate system is used for the numerical computations, which rely on a fourth‐order compact differentiation scheme, a third‐order explicit time stepping and a multi‐grid based pressure projection algorithm. The numerical techniques are described and a suite of validation studies are presented. The validation studies include a pointwise comparison of numerical simulations with both analytical solutions and laboratory measurements of non‐linear solitary wave propagation. Simulation results for flows lacking analytical or laboratory data are analysed a posteriori to demonstrate satisfaction of the potential energy balance. Computational results are compared with two‐layer hydraulic predictions in the case of exchange flow through a contracting channel. Finally, a simulation of circulation driven by spatially non‐uniform surface buoyancy flux in an irregular basin is discussed. Copyright © 2000 John Wiley & Sons, Ltd.     

Healing of shock instability for Roe's flux‐difference splitting scheme on triangular meshes

Healing of nonphysical flow solutions and shock instability from the use of Roe's flux‐difference splitting scheme is presented. The proposed method heals nonphysical flow solutions such as the carbuncle phenomenon, the shock instability from the odd–even decoupling problem, and the expansion shock generated from the violated entropy condition. The performance and efficiency of the proposed method are evaluated by solving several benchmark and complex high‐speed compressible flow problems. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

Euler solutions using flux‐based wave decomposition

The Euler equations are solved for non‐hydrostatic atmospheric flow problems in two dimensions using a high‐resolution Godunov‐type scheme. The Riemann problem is solved using a flux‐based wave decomposition suggested by LeVeque. This paper describes in detail, the design and implementation of the Riemann solver used for computing the Godunov fluxes. The methodology is then validated against benchmark cases for non‐hydrostatic atmospheric flows. Comparisons are made with solutions obtained from the National Center for Atmospheric Research's state‐of‐the‐art numerical model. The method shows promise in simulating non‐hydrostatic flows, which are characterized by steep gradients on the meso‐, micro‐ and urban‐scales. Copyright © 2006 John Wiley & Sons, Ltd.     

Godunov‐type adaptive grid model of wave–current interaction at cuspate beaches

This paper presents a second‐order accurate Godunov‐type numerical scheme for depth‐ and period‐averaged wave–current interaction. A flux Jacobian is derived for the wave conservation equations and its eigensystem determined, enabling Roe's approximate Riemann solver to be used to evaluate convective fluxes. Dynamically adaptive quadtree grids are used to focus on local hydrodynamic features, where sharp gradients occur in the flow variables. Adaptation criteria based on depth‐averaged vorticity, wave‐height gradient, wave steepness and the magnitude of velocity gradients are found to produce accurate solutions for nearshore circulation at a half‐sinusoidal beach. However, the simultaneous combination of two or more separate criteria produces numerical instability and interference unless all criteria are satisfied for mesh depletion. Simulations of wave–current interaction at a multi‐cusped beach match laboratory data from the United Kingdom Coastal Research Facility (UKCRF). A parameter study demonstrates the sensitivity of nearshore flow patterns to changes in relative cusp height, angle of wave incidence, bed roughness, offshore wave height and assumed turbulent eddy viscosity. Only a small deviation from normal wave incidence is required to initiate a meandering longshore current. Nearshore circulation patterns are highly dependent on the offshore wave height. Reduction of the assumed eddy viscosity parameter causes the primary circulation cells for normally incident waves to increase in strength whilst producing rip‐like currents cutting diagonally across the surf zone. Copyright © 2004 John Wiley & Sons, Ltd.     

An unstructured finite volume model for dam‐break floods with wet/dry fronts over complex topography

A robust, well‐balanced, unstructured, Godunov‐type finite volume model has been developed in order to simulate two‐dimensional dam‐break floods over complex topography with wetting and drying. The model is based on the nonlinear shallow water equations in hyperbolic conservation form. The inviscid fluxes are calculated using the HLLC approximate Riemann solver and a second‐order spatial accuracy is achieved by implementing the MUSCL reconstruction technique. To prevent numerical oscillations near shocks, slope‐limiting techniques are used for controlling the total variation of the reconstructed field. The model utilizes an explicit two‐stage Runge–Kutta method for time stepping, whereas implicit treatments for friction source terms. The novelties of the model include the flux correction terms and the water depth reconstruction method both for partially and fully submerged cells, and the wet/dry front treatments. The proposed flux correction terms combined with the water depth reconstruction method are necessary to balance the bed slope terms and flux gradient in the hydrostatical steady flow condition. Especially, this well‐balanced property is also preserved in partially submerged cells. It is found that the developed wet/dry front treatments and implicit scheme for friction source terms are stable. The model is tested against benchmark problems, laboratory experimental data, and realistic application related to dam‐break flood wave propagation over arbitrary topography. Numerical results show that the model performs satisfactorily with respect to its effectiveness and robustness and thus has bright application prospects. Copyright © 2010 John Wiley & Sons, Ltd.     

A fractional two‐step implicit algorithm using curvilinear collected grids

An efficient fractional two‐step implicit algorithm is reported to simulate incompressible fluid flows in a boundary‐fitted curvilinear collocated grid system. Using the finite volume method, the convection terms are discretized by the high‐accuracy Roe's scheme to minimize numerical diffusion. An implicitness coefficient Π is introduced to accelerate the rate of convergence. It is demonstrated that the proposed algorithm links the fractional step method to the pressure correction procedure, and the SIMPLEC method could be considered as a special case of the fractional two‐step implicit algorithm (when Π=1). The proposed algorithm is applicable to unsteady flows and steady flows. Three benchmark two‐dimensional laminar flows are tested to evaluate the performance of the proposed algorithm. Performance is measured by sensitivity analyses of the efficiency, accuracy, grid density, grid skewness and Reynolds number on the solutions. Results show that the model is efficient and robust. Copyright © 2006 John Wiley & Sons, Ltd.     

A numerical study of an unsteady laminar flow in a doubly constricted 3D vessel

Unsteady flow dynamics in doubly constricted 3D vessels have been investigated under pulsatile flow conditions for a full cycle of period T. The coupled non‐linear partial differential equations governing the mass and momentum of a viscous incompressible fluid has been numerically analyzed by a time accurate Finite Volume Scheme in an implicit Euler time marching setting. Roe's flux difference splitting of non‐linear terms and the pseudo‐compressibility technique employed in the current numerical scheme makes it robust both in space and time. Computational experiments are carried out to assess the influence of Reynolds' number and the spacing between two mild constrictions on the pressure drop across the constrictions. The study reveals that the pressure drop across a series of mild constrictions can get physiologically critical and is also found to be sensitive both to the spacing between the constrictions and the oscillatory nature of the inflow profile. The flow separation zone on the downstream constriction is seen to detach from the diverging wall of the constriction leading to vortex shedding with 3D features earlier than that on the wall in the spacing between the two constrictions. Copyright © 2002 John Wiley & Sons, Ltd.     

Weighted essential non‐oscillatory schemes for tidal bore on unstructured meshes

In this paper, the third‐order weighted essential non‐oscillatory (WENO) schemes are used to simulate the two‐dimensional shallow water equations with the source terms on unstructured meshes. The balance of the flux and the source terms makes the shallow water equations fit to non‐flat bottom questions. The simulation of a tidal bore on an estuary with trumpet shape and Qiantang river is performed; the results show that the schemes can be used to simulate the current flow accurately and catch the stronger discontinuous in water wave, such as dam break and tidal bore effectively. Copyright © 2008 John Wiley & Sons, Ltd.     

A Cartesian grid technique based on one‐dimensional integrated radial basis function networks for natural convection in concentric annuli

This paper reports a radial basis function (RBF)‐based Cartesian grid technique for the simulation of two‐dimensional buoyancy‐driven flow in concentric annuli. The continuity and momentum equations are represented in the equivalent stream function formulation that reduces the number of equations from three to one, but involves higher‐order derivatives. The present technique uses a Cartesian grid to discretize the problem domain. Along a grid line, one‐dimensional integrated RBF networks (1D‐IRBFNs) are employed to represent the field variables. The capability of 1D‐IRBFNs to handle unstructured points with accuracy is exploited to describe non‐rectangular boundaries in a Cartesian grid, while the method's ability to avoid the reduction of convergence rate caused by differentiation is instrumental in improving the quality of the approximation of higher‐order derivatives. The method is applied to simulate thermally driven flows in annuli between two circular cylinders and between an outer square cylinder and an inner circular cylinder. High Rayleigh number solutions are achieved and they are in good agreement with previously published numerical data. Copyright © 2007 John Wiley & Sons, Ltd.     

Improvements to compressible Euler methods for low‐Mach number flows

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

Simulation of free surface flows using the flux‐difference splitting scheme on the hybrid Cartesian/immersed boundary method

In this study, a method is developed to simulate the interaction between free surface flows and moving or deforming boundaries using the flux‐difference splitting scheme on the hybrid Cartesian/immersed boundary method. At each physical time step, the boundary is defined by an unstructured triangular surface grid. Immersed boundary (IB) nodes are distributed inside an instantaneous fluid domain based on edges crossing the boundary. At an IB node, dependent variables are reconstructed along the local normal line to the boundary. Inviscid fluxes are computed using Roe's flux‐difference splitting scheme for immiscible and incompressible fluids. The free surface is considered as a contact discontinuity in the density field. The motion of free surface is captured without any additional treatment along the fluid interface. The developed code is validated by comparisons with other experimental and computational results for a piston‐type wave maker, impulsive motion of a submerged circular cylinder, flow around a submerged hydrofoil, and Rayleigh–Taylor instability. The developed code is applied to simulate wave generation due to a continuously deforming bed beneath the free surface. The violent motion of a free surface caused by sloshing in a spherical tank is simulated. In this case, the free surface undergoes breakup and reconnection. Copyright © 2011 John Wiley & Sons, Ltd.     

A stabilized Nitsche‐type extended embedding mesh approach for 3D low‐ and high‐Reynolds‐number flows

This paper presents a stabilized extended finite element method (XFEM) based fluid formulation to embed arbitrary fluid patches into a fixed background fluid mesh. The new approach is highly beneficial when it comes to computational grid generation for complex domains, as it allows locally increased resolutions independent from size and structure of the background mesh. Motivating applications for such a domain decomposition technique are complex fluid‐structure interaction problems, where an additional boundary layer mesh is used to accurately capture the flow around the structure. The objective of this work is to provide an accurate and robust XFEM‐based coupling for low‐ as well as high‐Reynolds‐number flows. Our formulation is built from the following essential ingredients: Coupling conditions on the embedded interface are imposed weakly using Nitsche's method supported by extra terms to guarantee mass conservation and to control the convective mass transport across the interface for transient viscous‐dominated and convection‐dominated flows. Residual‐based fluid stabilizations in the interior of the fluid subdomains and accompanying face‐oriented fluid and ghost‐penalty stabilizations in the interface zone stabilize the formulation in the entire fluid domain. A detailed numerical study of our stabilized embedded fluid formulation, including an investigation of variants of Nitsche's method for viscous flows, shows optimal error convergence for viscous‐dominated and convection‐dominated flow problems independent of the interface position. Challenging two‐dimensional and three‐dimensional numerical examples highlight the robustness of our approach in all flow regimes: benchmark computations for laminar flow around a cylinder, a turbulent driven cavity flow at Re = 10000 and the flow interacting with a three‐dimensional flexible wall. Copyright © 2016 John Wiley & Sons, Ltd.     

A γ‐DGBGK scheme for compressible multi‐fluids

The paper presents a Discontinuous Galerkin γ‐BGK (γ‐DGBGK) method for compressible multicomponent flow simulations by coupling the discontinuous Galerkin method with a γ‐BGK scheme based on WENO limiters. In this γ‐DGBGK method, the construction of the flux in the DG method is based on the kinetic scheme which not only couples the convective and dissipative terms together, but also includes both discontinuous and continuous terms in the flux formulation at cell interfaces. WENO limiters are used to obtain uniform high‐order accuracy and sharp non‐oscillatory shock transition, and time accuracy obtained by integration for the flux function at the cell interface. Numerical examples in one and two space dimensions are presented to illustrate the robust and accuracy of the present scheme. Copyright © 2010 John Wiley & Sons, Ltd.