Semi‐Lagrangian semi‐implicit time‐splitting scheme for the shallow water equations

Time‐splitting technique applied in the context of the semi‐Lagrangian semi‐implicit method allows the use of extended time steps mainly based on physical considerations and reduces the number of numerical operations at each time step such that it is approximately proportional to the number of the points of spatial grid. To control time growth of the additional truncation errors, the standard stabilizing correction method is modified with no penalty for accuracy and efficiency of the algorithm. A linear analysis shows that constructed scheme is stable for time steps up to 2h. Numerical integrations with actual atmospheric fields of pressure and wind confirm computational efficiency, extended stability and accuracy of the proposed scheme. Copyright © 2006 John Wiley & Sons, Ltd.     

An operator splitting scheme for the stream‐function formulation of unsteady Navier–Stokes equations

A fictitious time is introduced into the unsteady equation of the stream function rendering it into a higher‐order ultra‐parabolic equation. The convergence with respect to the fictitious time (we call the latter ‘internal iterations’) allows one to obtain fully implicit nonlinear scheme in full time steps for the physical‐time variable. For particular choice of the artificial time increment, the scheme in full time steps is of second‐order of approximation. For the solution of the internal iteration, a fractional‐step scheme is proposed based on the splitting of the combination of the Laplace, bi‐harmonic and advection operators. A judicious choice for the time staggering of the different parts of the nonlinear advective terms allows us to prove that the internal iterations are unconditionally stable and convergent. We assess the number of operations needed per time step and show computational effectiveness of the proposed scheme. We prove that when the internal iterations converge, the scheme is second‐order in physical time and space, nonlinear, implicit and absolutely stable. The performance of the scheme is demonstrated for the flow created by oscillatory motion of the lid of a square cavity. All theoretical findings are demonstrated practically. Copyright © 2006 John Wiley & Sons, Ltd.     

An operator splitting algorithm for the three-dimensional advection–diffusion equation

Operator splitting algorithms are frequently used for solving the advection–diffusion equation, especially to deal with advection dominated transport problems. In this paper an operator splitting algorithm for the three-dimensional advection–diffusion equation is presented. The algorithm represents a second-order-accurate adaptation of the Holly and Preissmann scheme for three-dimensional problems. The governing equation is split into an advection equation and a diffusion equation, and they are solved by a backward method of characteristics and a finite element method, respectively. The Hermite interpolation function is used for interpolation of concentration in the advection step. The spatial gradients of concentration in the Hermite interpolation are obtained by solving equations for concentration gradients in the advection step. To make the composite algorithm efficient, only three equations for first-order concentration derivatives are solved in the diffusion step of computation. The higher-order spatial concentration gradients, necessary to advance the solution in a computational cycle, are obtained by numerical differentiations based on the available information. The simulation characteristics and accuracy of the proposed algorithm are demonstrated by several advection dominated transport problems. © 1998 John Wiley & Sons, Ltd.     

An upstream flux‐splitting finite‐volume scheme for 2D shallow water equations

An upstream flux‐splitting finite‐volume (UFF) scheme is proposed for the solutions of the 2D shallow water equations. In the framework of the finite‐volume method, the artificially upstream flux vector splitting method is employed to establish the numerical flux function for the local Riemann problem. Based on this algorithm, an UFF scheme without Jacobian matrix operation is developed. The proposed scheme satisfying entropy condition is extended to be second‐order‐accurate using the MUSCL approach. The proposed UFF scheme and its second‐order extension are verified through the simulations of four shallow water problems, including the 1D idealized dam breaking, the oblique hydraulic jump, the circular dam breaking, and the dam‐break experiment with 45° bend channel. Meanwhile, the numerical performance of the UFF scheme is compared with those of three well‐known upwind schemes, namely the Osher, Roe, and HLL schemes. It is demonstrated that the proposed scheme performs remarkably well for shallow water flows. The simulated results also show that the UFF scheme has superior overall numerical performances among the schemes tested. Copyright © 2005 John Wiley & Sons, Ltd.     

An accurate and efficient semi‐implicit method for section‐averaged free‐surface flow modelling

An accurate, efficient and robust numerical method for the solution of the section‐averaged De St. Venant equations of open channel flow is presented and discussed. The method consists in a semi‐implicit, finite‐volume discretization of the continuity equation capable to deal with arbitrary cross‐section geometry and in a semi‐implicit, finite‐difference discretization of the momentum equation. By using a proper semi‐Lagrangian discretization of the momentum equation, a highly efficient scheme that is particularly suitable for subcritical regimes is derived. Accurate solutions are obtained in all regimes, except in presence of strong unsteady shocks as in dam‐break cases. By using a suitable upwind, Eulerian discretization of the same equation, instead, a scheme capable of describing accurately also unsteady shocks can be obtained, although this scheme requires to comply with a more restrictive stability condition. The formulation of the two approaches allows a unified implementation and an easy switch between the two. The code is verified in a wide range of idealized test cases, highlighting its accuracy and efficiency characteristics, especially for long time range simulations of subcritical river flow. Finally, a model validation on field data is presented, concerning simulations of a flooding event of the Adige river. Copyright © 2009 John Wiley & Sons, Ltd.     

An operator splitting numerical scheme for thermal/isothermal incompressible viscous flows

A numerical scheme for time‐dependent incompressible viscous fluid flow, thermally coupled under the Boussinesq approximation is presented. The scheme combines an operator splitting in the time discretization and linear finite elements in the space discretization, and is an extension of one previously applied for isothermal incompressible viscous flow governed by the Navier–Stokes equations. To show the efficiency of the scheme, numerical results are presented for mixed convection, and natural convection at high Rayleigh numbers. Restricting the scheme to the isothermal case, some numerical results at high Reynolds numbers are included, i.e. the scheme is tested for a small viscosity and a large force term, which are not trivial tasks to deal with. Copyright © 1999 John Wiley & Sons, Ltd.     

Solution of the Navier–Stokes equations in velocity–vorticity form using a Eulerian–Lagrangian boundary element method

This paper describes the Eulerian–Lagrangian boundary element model for the solution of incompressible viscous flow problems using velocity–vorticity variables. A Eulerian–Lagrangian boundary element method (ELBEM) is proposed by the combination of the Eulerian–Lagrangian method and the boundary element method (BEM). ELBEM overcomes the limitation of the traditional BEM, which is incapable of dealing with the arbitrary velocity field in advection‐dominated flow problems. The present ELBEM model involves the solution of the vorticity transport equation for vorticity whose solenoidal vorticity components are obtained iteratively by solving velocity Poisson equations involving the velocity and vorticity components. The velocity Poisson equations are solved using a boundary integral scheme and the vorticity transport equation is solved using the ELBEM. Here the results of two‐dimensional Navier–Stokes problems with low–medium Reynolds numbers in a typical cavity flow are presented and compared with a series solution and other numerical models. The ELBEM model has been found to be feasible and satisfactory. Copyright © 2000 John Wiley & Sons, Ltd.     

Interfacial flow computations using adaptive Eulerian–Lagrangian method for spacecraft applications

Understanding the interfacial dynamics and fluid physics associated with the operation of spacecraft is important for scientific as well as engineering purposes. To help address the issues associated with moving boundaries, interfacial dynamics, and spatial‐temporal variations in time and length scales, a 3‐D adaptive Eulerian–Lagrangian method is extended and further developed. The stationary (Eulerian) Cartesian grid is adopted to resolve the fluid flow, and the marker‐based triangulated moving (Lagrangian) surface meshes are utilized to treat the phase boundary. The key concepts and numerical procedures, along with the selected interfacial flow problems are presented. Specifically, the liquid fuel draining dynamics in different flow regimes, and the liquid surface stability under vertically oscillating gravitational acceleration are investigated. Direct assessment of experimental measurement and scaling analysis is made to highlight the computational performance of the present approach as well as the key fluid physics influenced by the given flow parameters. Copyright © 2010 John Wiley & Sons, Ltd.     

A method for computing compressible,highly stratified flows in astrophysics based on operator splitting

A new numerical approach based on consistent operator splitting is presented for computing compressible, highly stratified flows in astrophysics. The algorithm is particularly designed to search for steady or almost steady solutions for the time-dependent Navier–Stokes equations, describing viscous flow under the influence of a strong gravitational field. The algorithm proposed is multidimensional and works in Cartesian, cylindrical or spherical co-ordinates. It uses a second-order finite volume scheme with third-order upwinding and a second-order time discretization. An adaptive time step control and monotonic multilevel grid distribution has been incorporated to speed up convergence. This method has been incorporated into a hydrodynamical code by which, for the first time, for two-dimensional models the dynamics of the boundary layer in the accretion disk around a compact star could be computed over the whole viscous time scale. © 1998 John Wiley & Sons, Ltd.     

A well‐balanced upstream flux‐splitting finite‐volume scheme for shallow‐water flow simulations with irregular bed topography

This study extends the upstream flux‐splitting finite‐volume (UFF) scheme to shallow water equations with source terms. Coupling the hydrostatic reconstruction method (HRM) with the UFF scheme achieves a resultant numerical scheme that adequately balances flux gradients and source terms. The proposed scheme is validated in three benchmark problems and applied to flood flows in the natural/irregular river with bridge pier obstructions. The results of the simulations are in satisfactory agreement with the available analytical solutions, experimental data and field measurements. Comparisons of the present results with those obtained by the surface gradient method (SGM) demonstrate the superior stability and higher accuracy of the HRM. The stability test results also show that the HRM requires less CPU time (up to 60%) than the SGM. The proposed well‐balanced UFF scheme is accurate, stable and efficient to solve flow problems involving irregular bed topography. Copyright © 2009 John Wiley & Sons, Ltd.     

A multigrid method for steady incompressible navier-stokes equations based on flux difference splitting

The steady Navier–Stokes equations in primitive variables are discretized in conservative form by a vertex-centred finite volume method Flux difference splitting is applied to the convective part to obtain an upwind discretization. The diffusive part is discretized in the central way. In its first-order formulation, flux difference splitting leads to a discretization of so-called vector positive type. This allows the use of classical relaxation methods in collective form. An alternating line Gauss–Seidel relaxation method is chosen here. This relaxation method is used as a smoother in a multigrid method. The components of this multigrid method are: full approximation scheme with F-cycles, bilinear prolongation, full weighting for residual restriction and injection of grid functions. Higher-order accuracy is achieved by the flux extrapolation method. In this approach the first-order convective fluxes are modified by adding second-order corrections involving flux limiting. Here the simple MinMod limiter is chosen. In the multigrid formulation the second-order discrete system is solved by defect correction. Computational results are shown for the well known GAMM backward-facing step problem and for a channel with a half-circular obstruction.     

A well‐balanced explicit/semi‐implicit finite element scheme for shallow water equations in drying–wetting areas

Numerical solutions of the shallow water equations can be used to reproduce flow hydrodynamics occurring in a wide range of regions. In hydraulic engineering, the objectives include the prediction of dam break wave propagation, fluvial floods and other catastrophic flooding phenomena, the modeling of estuarine and coastal circulations, and the design and optimization of hydraulic structures. In this paper, a well‐balanced explicit and semi‐implicit finite element scheme for shallow water equations over complex domains involving wetting and drying is proposed. The governing equations are discretized by a fractional finite element method using a two‐step Taylor–Galerkin scheme. First, the intermediate increment of conserved variable is obtained explicitly neglecting the pressure gradient term. This is then corrected for the effects of pressure once the pressure increment has been obtained from the Poisson equation. In order to maintain the ‘well‐balanced’ property, the pressure gradient term and bed slope terms are incorporated into the Poisson equation. Moreover, a local bed slope modification technique is employed in drying–wetting interface treatments. The proposed model is well validated against several theoretical benchmark tests. Copyright © 2014 John Wiley & Sons, Ltd.     

A multigrid method for steady incompressible Navier–Stokes equations based on partial flux splitting

Flux splitting is applied to the convective part of the steady Navier–Stokes equations for incompressible flow. Partial upwind differences are introduced in the split first-order part, while central differences are used in the second-order part. The discrete set of equations obtained is positive, so that it can be solved by collective variants of relaxation methods. The partial upwinding is optimized in the same way as for a scalar convection–diffusion equation, but involving several Peclet numbers. It is shown that with the optimum partial upwinding accurate results can be obtained. A full multigrid method in W-cycle form, using red–black successive under-relaxation, injection and bilinear interpolation, is described. The efficiency of this method is demonstrated.     

Multirate time stepping for accelerating explicit discontinuous Galerkin computations with application to geophysical flows

This paper presents multirate explicit time‐stepping schemes for solving partial differential equations with discontinuous Galerkin elements in the framework of Large‐scale marine flows. It addresses the variability of the local stable time steps by gathering the mesh elements in appropriate groups. The real challenge is to develop methods exhibiting mass conservation and consistency. Two multirate approaches, based on standard explicit Runge–Kutta methods, are analyzed. They are well suited and optimized for the discontinuous Galerkin framework. The significant speedups observed for the hydrodynamic application of the Great Barrier Reef confirm the theoretical expectations. Copyright © 2012 John Wiley & Sons, Ltd.     

A flux‐splitting solver for shallow water equations with source terms

This paper introduces a stable flux‐splitting solver for one‐dimensional (1D) shallow water equations. This solver is specifically designed to satisfy a strengthened consistency condition for stationary solutions that ensures the stability and accuracy of the scheme. It applies to channels with variable depth and width, including terms modelling friction at bottom and vertical walls. Some numerical tests by comparison to both analytical solutions and experimental measurements show the good performances of the scheme. Copyright © 2003 John Wiley & Sons, Ltd.     

A consistent splitting scheme for unsteady incompressible viscous flows I. Dirichlet boundary condition and applications

A well‐recognized approach for handling the incompressibility constraint by operating directly on the discretized Navier–Stokes equations is used to obtain the decoupling of the pressure from the velocity field. By following the current developments by Guermond and Shen, the possibilities of obtaining accurate pressure and reducing boundary‐layer effect for the pressure are analysed. The present study mainly reports the numerical solutions of an unsteady Navier–Stokes problem based on the so‐called consistent splitting scheme (J. Comput. Phys. 2003; 192 :262–276). At the same time the Dirichlet boundary value conditions are considered. The accuracy of the method is carefully examined against the exact solution for an unsteady flow physics problem in a simply connected domain. The effectiveness is illustrated viz. several computations of 2D double lid‐driven cavity problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Use of two‐step split‐operator approach for 2D shallow water flow computation

The objective of this paper is to present a methodology of using a two‐step split‐operator approach for solving the shallow water flow equations in terms of an orthogonal curvilinear co‐ordinate system. This approach is in fact one kind of the so‐called fractional step method that has been popularly used for computations of dynamic flow. By following that the momentum equations are decomposed into two portions, the computation procedure involves two steps. The first step (dispersion step) is to compute the provisional velocity in the momentum equation without the pressure gradient. The second step (propagation step) is to correct the provisional velocity by considering a divergence‐free velocity field, including the effect of the pressure gradient. This newly proposed method, other than the conventional split‐operator methods, such as the projection method, considers the effects of pressure gradient and bed friction in the second step. The advantage of this treatment is that it increases flexibility, efficiency and applicability of numerical simulation for various hydraulic problems. Four cases, including back‐water flow, reverse flow, circular basin flow and unsteady flow, have been demonstrated to show the accuracy and practical application of the method. Copyright © 1999 John Wiley & Sons, Ltd.     

Edge‐based finite element method for shallow water equations

This paper describes an edge‐based implementation of the generalized residual minimum (GMRES) solver for the fully coupled solution of non‐linear systems arising from finite element discretization of shallow water equations (SWEs). The gain in terms of memory, floating point operations and indirect addressing is quantified for semi‐discrete and space–time analyses. Stabilized formulations, including Petrov–Galerkin models and discontinuity‐capturing operators, are also discussed for both types of discretization. Results illustrating the quality of the stabilized solutions and the advantages of using the edge‐based approach are presented at the end of the paper. Copyright © 2001 John Wiley & Sons, Ltd.     

A modular,operator‐splitting scheme for fluid–structure interaction problems with thick structures

We present an operator‐splitting scheme for fluid–structure interaction (FSI) problems in hemodynamics, where the thickness of the structural wall is comparable to the radius of the cylindrical fluid domain. The equations of linear elasticity are used to model the structure, while the Navier–Stokes equations for an incompressible viscous fluid are used to model the fluid. The operator‐splitting scheme, based on the Lie splitting, separates the elastodynamics structure problem from a fluid problem in which structure inertia is included to achieve unconditional stability. We prove energy estimates associated with unconditional stability of this modular scheme for the full nonlinear FSI problem defined on a moving domain, without requiring any sub‐iterations within time steps. Two numerical examples are presented, showing excellent agreement with the results of monolithic schemes. First‐order convergence in time is shown numerically. Modularity, unconditional stability without temporal sub‐iterations, and simple implementation are the features that make this operator‐splitting scheme particularly appealing for multi‐physics problems involving FSI. Copyright © 2013 John Wiley & Sons, Ltd.     

An eigenvector‐based linear reconstruction approach for time stepping in discontinuous Galerkin scheme used to solve shallow water equations

Discontinuous Galerkin (DG) methods have shown promising results for solving the two‐dimensional shallow water equations. In this paper, the classical Runge–Kutta (RK) time discretisation is replaced by the eigenvector‐based reconstruction (EVR) that allows the second‐order time accuracy to be achieved within a single time‐stepping procedure. Moreover, the EVRDG approach yields stable solutions near drying and wetting fronts, whereas the classical RKDG approach yields instabilities. The proposed EVRDG technique is compared with the original RKDG approach on various test cases with analytical solutions. The EVRDG solutions are shown to be as accurate as those obtained with the RKDG scheme. Besides, the EVRDG scheme is 1.6 times faster than the RKDG method. Simulating dambreaks involving dry beds confirms that EVRDG scheme gives correct solutions, whereas the RKDG method yields instabilities. Copyright © 2012 John Wiley & Sons, Ltd.