On the efficiency of semi‐implicit and semi‐Lagrangian spectral methods for the calculation of incompressible flows

Classical semi‐implicit backward Euler/Adams–Bashforth time discretizations of the Navier–Stokes equations induce, for high‐Reynolds number flows, severe restrictions on the time step. Such restrictions can be relaxed by using semi‐Lagrangian schemes essentially based on splitting the full problem into an explicit transport step and an implicit diffusion step. In comparison with the standard characteristics method, the semi‐Lagrangian method has the advantage of being much less CPU time consuming where spectral methods are concerned. This paper is devoted to the comparison of the ‘semi‐implicit’ and ‘semi‐Lagrangian’ approaches, in terms of stability, accuracy and computational efficiency. Numerical results on the advection equation, Burger's equation and finally two‐ and three‐dimensional Navier–Stokes equations, using spectral elements or a collocation method, are provided. Copyright © 2001 John Wiley & Sons, Ltd.     

A semi‐implicit scheme for 3D free surface flows with high‐order velocity reconstruction on unstructured Voronoi meshes

In this paper, we present a computationally efficient semi‐implicit scheme for the simulation of three‐dimensional hydrostatic free surface flow problems on staggered unstructured Voronoi meshes. For each polygonal control volume, the pressure is defined in the cell center, whereas the discrete velocity field is given by the normal velocity component at the cell faces. A piecewise high‐order polynomial vector velocity field is then reconstructed from the scalar normal velocities at the cell faces by using a new high‐order constrained least‐squares reconstruction operator. The reconstructed high‐order piecewise polynomial velocity field is used for trajectory integration in a semi‐Lagrangian approach to discretize the nonlinear convective terms in the governing PDE. For that purpose, a high‐order Taylor method is used as ODE integrator. The resulting semi‐implicit algorithm is extensively validated on a large set of different academic test problems with exact analytical solution and is finally applied to a real‐world engineering problem consisting of a curved channel upstream of two micro‐turbines of a hydroelectric power plant. For this realistic case, some experimental reference data are available from field measurements. Copyright © 2012 John Wiley & Sons, Ltd.     

High‐order 𝒞1 finite‐element interpolating schemes—Part I: Semi‐Lagrangian linear advection

This paper is devoted to the development of accurate high‐order interpolating schemes for semi‐Lagrangian advection. The characteristic‐Galerkin formulation is obtained by using a semi‐Lagrangian temporal discretization of the total derivative. The semi‐Lagrangian method requires high‐order interpolators for accuracy. A class of ??¹ finite‐element interpolating schemes is developed and two semi‐Lagrangian methods are considered by tracking the feet of the characteristic lines either from the interpolation or from the integration nodes. Numerical stability and analytical results quantifying the amount of artificial viscosity induced by the two methods are presented in the case of the one‐dimensional linear advection equation, based on the modified equation approach. Results of test problems to simulate the linear advection of a cosine hill illustrate the performance of the proposed approach. Copyright © 2007 John Wiley & Sons, Ltd.     

An entropy fixed cell‐centered Lagrangian scheme

On the basis of the work [P.‐H. Maire, R. Abgrall, J. Breil, J. Ovadia, SIAM J. Sci. Comput. 29 (2007), 1781–1824], we present an entropy fixed cell‐centered Lagrangian scheme for solving the Euler equations of compressible gas dynamics. The scheme uses the fully Lagrangian form of the gas dynamics equations, in which the primary variables are cell‐centered. And using the nodal solver, we obtain the nodal viscous‐velocity, viscous‐pressures, antidissipation velocity, and antidissipation pressures of each node. The final nodal velocity is computed as a weighted sum of viscous‐velocity and antidissipation velocity, so do nodal pressures, whereas these weights are calculated through the total entropy conservation for isentropic flows. Consequently, the constructed scheme is conservative in mass, momentum, and energy; preserves entropy for isentropic flows, and satisfies a local entropy inequality for nonisentropic flows. One‐ and two‐dimensional numerical examples are presented to demonstrate theoretical analysis and performance of the scheme in terms of accuracy and robustness.Copyright © 2013 John Wiley & Sons, Ltd.     

Some considerations for high‐order ‘incremental remap’‐based transport schemes: edges,reconstructions, and area integration

The problem of two‐dimensional tracer advection on the sphere is extremely important in modeling of geophysical fluids and has been tackled using a variety of approaches. A class of popular approaches for tracer advection include ‘incremental remap’ or cell‐integrated semi‐Lagrangian‐type schemes. These schemes achieve high‐order accuracy without the need for multistage integration in time, are capable of large time steps, and tend to be more efficient than other high‐order transport schemes when applied to a large number of tracers over a single velocity field. In this paper, the simplified flux‐form implementation of the Conservative Semi‐LAgrangian Multi‐tracer scheme (CSLAM) is reformulated using quadratic curves to approximate the upstream flux volumes and Gaussian quadrature for integrating the edge flux. The high‐order treatment of edge fluxes is motivated because of poor accuracy of the CSLAM scheme in the presence of strong nonlinear shear, such as one might observe in the midlatitudes near an atmospheric jet. Without the quadratic treatment of upstream edges, we observe at most second‐order accuracy under convergence of grid resolution, which is returned to third‐order accuracy under the improved treatment. A shallow‐water barotropic instability also reveals clear evidence of grid imprinting without the quadratic correction. Consequently, these tests reveal a problem that might arise in tracer transport near nonlinearly sheared regions of the real atmosphere, particularly near cubed‐sphere panel edges. Although CSLAM is used as the foundation for this analysis, the conclusions of this paper are applicable to the general class of incremental remap schemes. Copyright © 2012 John Wiley & Sons, Ltd.     

Particle trajectory calculations with a two‐step three‐time level semi‐Lagrangian scheme well suited for curved flows

This study proposes a new two‐step three‐time level semi‐Lagrangian scheme for calculation of particle trajectories. The scheme is intended to yield accurate determination of the particle departure position, particularly in the presence of significant flow curvature. Experiments were performed both for linear and non‐linear idealized advection problems, with different flow curvatures. Results for simulations with the proposed scheme, and with three other semi‐Lagrangian schemes, and with an Eulerian method are presented. In the linear advection problem the two‐step three‐time level scheme produced smaller root mean square errors and more accurate replication of the angular displacement of a Gaussian hill than the other schemes. In the non‐linear advection experiments the proposed scheme produced, in general, equal or better conservation of domain‐averaged quantities than the other semi‐Lagrangian schemes, especially at large Courant numbers. In idealized frontogenesis simulations the scheme performed equally or better than the other schemes in the representation of sharp gradients in a scalar field. The two‐step three‐time level scheme has some computational overhead as compared with the other three semi‐Lagrangian schemes. Nevertheless, the additional computational effort was shown to be worthwhile, due to the accuracy obtained by the scheme in the experiments with large time steps. The most remarkable feature of the scheme is its robustness, since it performs well both for small and large Courant numbers, in the presence of weak as well strong flow curvatures. Copyright © 2009 John Wiley & Sons, Ltd.     

A multi‐moment transport model on cubed‐sphere grid

A conservative, single‐cell‐based semi‐Lagrangian transport model is proposed in this paper. Using multi‐moment concept, an additional moment, i.e. volume‐integrated average (VIA), is treated as the model variable besides the point value (PV) updated in the traditional semi‐Lagrangian schemes. A quadratic interpolation function is constructed based on local degrees of freedom defined within each single cell. The PV moment is advanced by the semi‐Lagrangian formulation, whereas the VIA moment is updated by a finite volume formulation to rigorously ensure the numerical conservation. The numerical fluxes are computed from the PV moments defined along the boundary edges of the control volume. The scheme is extended to the spherical geometry through the application of the cubed‐sphere grid that eliminates the polar singularity in the conventional longitude/latitude coordinates by using the quasi‐uniform grid spacing covering the whole sphere. The single‐cell‐based scheme is well suited for the treatment of the connections between different patches. A simple quasi‐monotone limiter to the PV moment is applied to suppress non‐physical oscillations. The proposed scheme has been validated via representative benchmark tests and the performance is competitive to other existing transport schemes. Copyright © 2010 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.     

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.     

The characteristic‐based split (CBS) meshfree method for free surface flow problems in ALE formulation

A semi‐implicit characteristic‐based split (CBS) meshfree algorithm in the arbitrary Lagrangian Eulerian (ALE) framework is proposed for the numerical solution of incompressible free surface flow problem in the paper. The algorithm is the extension of general CBS method which was initially introduced in finite element framework, this is due to the fact that CBS method not only can enhance the stability, but also avoid LBB condition when equal order basis function is used to approximate velocity and pressure variables. Meanwhile, a simple way for node update and node speed calculation is developed which is used to capture the free surface exactly. The numerical solutions are compared with available analytical and numerical solutions, which shows that the proposed method has better ability to simulate the free surface incompressible flow problem. Copyright © 2009 John Wiley & Sons, Ltd.     

A semi‐implicit scheme for large Eddy simulation of piston engine flow and combustion

A semi‐implicit scheme is presented for large eddy simulation of turbulent reactive flow and combustion in reciprocating piston engines. First, the governing equations in a deforming coordinate system are formulated to accommodate the moving piston. The numerical scheme is made up of a fourth‐order central difference for the diffusion terms in the transport equations and a fifth‐order weighted essentially nonoscillatory (WENO) scheme for the convective terms. A second‐ order Adams–Bashforth scheme is used for time integration. For higher density ratios, it is combined with a predictor–corrector scheme. The numerical scheme is explicit for time integration of the transport equations, except for the continuity equation which is used together with the momentum equation to determine the pressure field and velocity field by using a Poisson equation for the pressure correction field. The scheme is aimed at the simulation of low Mach number flows typically found in piston engines. An efficient multigrid method that can handle high grid aspect ratio is presented for solving the pressure correction equation. The numerical scheme is evaluated on two test engines, a laboratory four‐stroke engine with rectangular‐shaped engine geometry where detailed velocity measurements are available, and a modified truck engine with practical cylinder geometry where lean ethanol/air mixture is combusted under a homogeneous charge compression ignition (HCCI) condition. Copyright © 2012 John Wiley & Sons, Ltd.     

An unconditionally stable semi‐Lagrangian method for the spherical atmospherical shallow water equations

A semi‐implicit, semi‐Lagrangian, mixed finite difference–finite volume model for the shallow water equations on a rotating sphere is introduced and discussed. Its main features are the vectorial treatment of the momentum equation and the finite volume approach for the continuity equation. Pressure and Coriolis terms in the momentum equation and velocity in the continuity equation are treated semi‐implicitly. Moreover, a splitting technique is introduced to preserve symmetry of the numerical scheme. An alternative asymmetric scheme (without splitting) is also introduced and the efficiency of both is discussed. The model is shown to be conservative in geopotential height and unconditionally stable for 0.5≤θ≤1. Numerical experiments on two standard test problems confirm the performance of the model. Copyright © 2000 John Wiley & Sons, Ltd.     

A conservative local interface sharpening scheme for the constrained interpolation profile method

A conservative local interface sharpening scheme has been developed for the constrained interpolation profile method with the conservative semi‐Lagrangian scheme, because the conservative semi‐Lagrangian scheme does not feature a mechanism to control the interface thickness, thus causing an increase of numerical error with the advance of the time step. The proposed sharpening scheme is based on the conservative level set method proposed by Olsson and Kreiss. However, because their method can cause excessive deformation of the free‐surface in certain circumstances, we propose an improvement of the method by developing a local sharpening technique. Several advection tests are presented to assess the correctness of the advection and the improved interface sharpening scheme. This is followed by the validations of dam‐breaking flow and the rising bubble flows. The mass of the fluid is exactly conserved and the computed terminal velocity of the rising bubble agrees well with the experiments compared with other numerical methods such as the volume of fluid method (VOF), the front tracking method, and the level set method. Copyright © 2011 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.     

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.     

Depth‐integrated nonhydrostatic free‐surface flow modeling using weighted‐averaged equations

In this study, a depth‐integrated nonhydrostatic flow model is developed using the method of weighted residuals. Using a unit weighting function, depth‐integrated Reynolds‐averaged Navier‐Stokes equations are obtained. Prescribing polynomial variations for the field variables in the vertical direction, a set of perturbation parameters remains undetermined. The model is closed generating a set of weighted‐averaged equations using a suitable weighting function. The resulting depth‐integrated nonhydrostatic model is solved with a semi‐implicit finite‐volume finite‐difference scheme. The explicit part of the model is a Godunov‐type finite‐volume scheme that uses the Harten‐Lax‐van Leer‐contact wave approximate Riemann solver to determine the nonhydrostatic depth‐averaged velocity field. The implicit part of the model is solved using a Newton‐Raphson algorithm to incorporate the effects of the pressure field in the solution. The model is applied with good results to a set of problems of coastal and river engineering, including steady flow over fixed bedforms, solitary wave propagation, solitary wave run‐up, linear frequency dispersion, propagation of sinusoidal waves over a submerged bar, and dam‐break flood waves.     

An efficient edge‐based level set finite element method for free surface flow problems

We present an efficient technique for the solution of free surface flow problems using level set and a parallel edge‐based finite element method. An unstructured semi‐explicit solution scheme is proposed. A custom data structure, obtained by blending node‐based and edge‐based approaches is presented so to allow a good parallel performance. In addition to standard velocity extrapolation (for the convection of the level set function), an explicit extrapolation of the pressure field is performed in order to impose both the pressure boundary condition and the volume conservation. The latter is also improved with a modification of the divergence free constrain. The method is shown to allow an efficient solution of both simple benchmark cases and complex industrial examples. Copyright © 2012 John Wiley & Sons, Ltd.     

Third‐order‐accurate semi‐implicit Runge–Kutta scheme for incompressible Navier–Stokes equations

A semi‐implicit three‐step Runge–Kutta scheme for the unsteady incompressible Navier–Stokes equations with third‐order accuracy in time is presented. The higher order of accuracy as compared to the existing semi‐implicit Runge–Kutta schemes is achieved due to one additional inversion of the implicit operator I‐τγL, which requires inversion of tridiagonal matrices when using approximate factorization method. No additional solution of the pressure‐Poisson equation or evaluation of Navier–Stokes operator is needed. The scheme is supplied with a local error estimation and time‐step control algorithm. The temporal third‐order accuracy of the scheme is proved analytically and ascertained by analysing both local and global errors in a numerical example. Copyright © 2005 John Wiley & Sons, Ltd.     

A dual‐weighted trust‐region adaptive POD 4‐D Var applied to a finite‐volume shallow water equations model on the sphere

In this paper we study solutions of an inverse problem for a global shallow water model controlling its initial conditions specified from the 40‐yr ECMWF Re‐analysis (ERA‐40) data sets, in the presence of full or incomplete observations being assimilated in a time interval (window of assimilation) with or without background error covariance terms. As an extension of the work by Chen et al. (Int. J. Numer. Meth. Fluids 2009), we attempt to obtain a reduced order model of the above inverse problem, based on proper orthogonal decomposition (POD), referred to as POD 4D‐Var for a finite volume global shallow water equation model based on the Lin–Rood flux‐form semi‐Lagrangian semi‐implicit time integration scheme. Different approaches of POD implementation for the reduced inverse problem are compared, including a dual‐weighted method for snapshot selection coupled with a trust‐region POD adaptivity approach. Numerical results with various observational densities and background error covariance operator are also presented. The POD 4‐D Var model results combined with the trust‐region adaptivity exhibit similarity in terms of various error metrics to the full 4D Var results, but are obtained using a significantly lesser number of minimization iterations and require lesser CPU time. Based on our previous and current work, we conclude that POD 4‐D Var certainly warrants further studies, with promising potential of its extension to operational 3‐D numerical weather prediction models. Copyright © 2011 John Wiley & Sons, Ltd.