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 Godunov‐type fractional semi‐implicit method based on staggered grid for dam‐break modeling

A semi‐implicit finite volume model based upon staggered grid is presented for solving shallow water equation. The model employs a time‐splitting scheme that uses a predictor–corrector method for the advection term. The fluxes are calculated based on a Riemann solver in the prediction step and a downwind scheme in the correction step. A simple TVD scheme is employed for shock capturing purposes in which the Minmond limiter is used for flux functions. As a consequence of using staggered grid, an ADI method is adopted for solving the discretized equations for 2‐D problems. Several 1‐D and 2‐D flows have been modeled with satisfactory results when compared with analytical and experimental test cases. The model is also capable of simulating supercritical as well as subcritical flow. Copyright © 2010 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.     

A Godunov method for the computation of erosional shallow water transients

A Godunov method is proposed for the computation of open‐channel flows in conditions of rapid bed erosion and intense sediment transport. Generalized shallow water equations govern the evolution of three distinct interfaces: the water free‐surface, the boundary between pure water and a sediment transport layer, and the morphodynamic bottom profile. Based on the HLL scheme of Harten, Lax and Van Leer (1983), a finite volume numerical solver is constructed, then extended to second‐order accuracy using Strang splitting and MUSCL extrapolation. Lateralisation of the momentum flux is adopted to handle the non‐conservative product associated with bottom slope. Computational results for erosional dam‐break waves are compared with experimental measurements and semi‐analytical Riemann solutions. Copyright © 2003 John Wiley & Sons, Ltd.     

A semi‐implicit finite difference method for non‐hydrostatic,free‐surface flows

In this paper a semi‐implicit finite difference model for non‐hydrostatic, free‐surface flows is analyzed and discussed. It is shown that the present algorithm is generally more accurate than recently developed models for quasi‐hydrostatic flows. The governing equations are the free‐surface Navier–Stokes equations defined on a general, irregular domain of arbitrary scale. The momentum equations, the incompressibility condition and the equation for the free‐surface are integrated by a semi‐implicit algorithm in such a fashion that the resulting numerical solution is mass conservative and unconditionally stable with respect to the gravity wave speed, wind stress, vertical viscosity and bottom friction. Copyright © 1999 John Wiley & Sons, Ltd.     

Further experiences with computing non‐hydrostatic free‐surface flows involving water waves

A semi‐implicit, staggered finite volume technique for non‐hydrostatic, free‐surface flow governed by the incompressible Euler equations is presented that has a proper balance between accuracy, robustness and computing time. The procedure is intended to be used for predicting wave propagation in coastal areas. The splitting of the pressure into hydrostatic and non‐hydrostatic components is utilized. To ease the task of discretization and to enhance the accuracy of the scheme, a vertical boundary‐fitted co‐ordinate system is employed, permitting more resolution near the bottom as well as near the free surface. The issue of the implementation of boundary conditions is addressed. As recently proposed by the present authors, the Keller‐box scheme for accurate approximation of frequency wave dispersion requiring a limited vertical resolution is incorporated. The both locally and globally mass conserved solution is achieved with the aid of a projection method in the discrete sense. An efficient preconditioned Krylov subspace technique to solve the discretized Poisson equation for pressure correction with an unsymmetric matrix is treated. Some numerical experiments to show the accuracy, robustness and efficiency of the proposed method are presented. Copyright © 2004 John Wiley & Sons, Ltd.     

A high‐order fully explicit flux‐form semi‐Lagrangian shallow‐water model

A new approach is proposed for constructing a fully explicit third‐order mass‐conservative semi‐Lagrangian scheme for simulating the shallow‐water equations on an equiangular cubed‐sphere grid. State variables are staggered with velocity components stored pointwise at nodal points and mass variables stored as element averages. In order to advance the state variables in time, we first apply an explicit multi‐step time‐stepping scheme to update the velocity components and then use a semi‐Lagrangian advection scheme to update the height field and tracer variables. This procedure is chosen to ensure consistency between dry air mass and tracers, which is particularly important in many atmospheric chemistry applications. The resulting scheme is shown to be competitive with many existing numerical methods on a suite of standard test cases and demonstrates slightly improved performance over other high‐order finite‐volume models. Copyright © 2014 John Wiley & Sons, Ltd.     

Consistent hybrid finite volume/element formulations: model and complex viscoelastic flows

The accuracy and consistency of a new cell‐vertex hybrid finite element/volume scheme are investigated for viscoelastic flows. Finite element (FE) discretization is employed for the momentum and continuity equation, with finite volume (FV) applied to the constitutive law for stress. Here, the interest is to explore the consequences of utilizing conventional cell‐vertex methodology for an Oldroyd‐B model and to demonstrate resulting drawbacks in the presence of complex source terms on structured and unstructured grids. Alternative strategies worthy of consideration are presented. It is demonstrated how high‐order accuracy may be achieved in steady state by respecting consistency in the formulation. Both FE and FV spatial discretizations are embedded in the scheme, with FV triangular sub‐cells referenced within parent triangular finite elements. Both model and complex flow problems are selected to quantify and assess accuracy, appealing to analysis and experimental validation. The test problem is that of steady sink flow, a pure extensional flow, which reflects some of the numerical difficulties involved in solving more generalized viscoelastic flows, where both source and flux terms may contribute equally to stress propagation. In addition, a complex transient filament‐stretching flow is chosen to compute the evolution of stress fields within liquid bridges. Shortcomings of the various stress upwinding schemes are discussed in this context, whilst dealing with such free‐surface type problems. Here, stress fluctuation distribution alone is advocated, and a Lax‐scheme is found to deliver accuracy and stability to the computational results, comparing well with the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical modelling of gas–solid fluidized beds using the two‐fluid approach

This paper details an approach to modelling gas–solid fluidized beds using the two‐fluid granular temperature model. Details concerning the difficulties associated with the boundary conditions, particularly for curved boundaries, are described along with a novel means of obtaining the internal stress of the solid‐phase, in part, by solving an implicit equation. This results in a scheme that is stable even when the solid volume fraction is close to maximum packing. A transient, mixed finite element discretization is used to solve the multi‐phase equations with a discontinuous finite element representation of the granular temperature and continuity equations. A new solution method is proposed to solve the coupled momentum and continuity equations based on Arnoldi iteration. Two fluidized beds are modelled, one in the bubbling regime and the other in the slugging regime. These simulations are compared with experiments. Copyright © 2001 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.     

An efficient finite difference scheme for free‐surface flows in narrow rivers and estuaries

This paper presents a free‐surface correction (FSC) method for solving laterally averaged, 2‐D momentum and continuity equations. The FSC method is a predictor–corrector scheme, in which an intermediate free surface elevation is first calculated from the vertically integrated continuity equation after an intermediate, longitudinal velocity distribution is determined from the momentum equation. In the finite difference equation for the intermediate velocity, the vertical eddy viscosity term and the bottom‐ and sidewall friction terms are discretized implicitly, while the pressure gradient term, convection terms, and the horizontal eddy viscosity term are discretized explicitly. The intermediate free surface elevation is then adjusted by solving a FSC equation before the intermediate velocity field is corrected. The finite difference scheme is simple and can be easily implemented in existing laterally averaged 2‐D models. It is unconditionally stable with respect to gravitational waves, shear stresses on the bottom and side walls, and the vertical eddy viscosity term. It has been tested and validated with analytical solutions and field data measured in a narrow, riverine estuary in southwest Florida. Model simulations show that this numerical scheme is very efficient and normally can be run with a Courant number larger than 10. It can be used for rivers where the upstream bed elevation is higher than the downstream water surface elevation without any problem. Copyright © 2003 John Wiley & Sons, Ltd.     

A semi‐implicit numerical method for the free‐surface Navier–Stokes equations

Semi‐implicit methods are known for being the basis of simple, efficient, accurate, and stable numerical algorithms for simulating a large variety of geophysical free‐surface flows. Geophysical flows are typically characterized by having a small vertical scale as compared with their horizontal extents. Hence, the hydrostatic approximation often applies, and the free surface can be conveniently represented by a single‐valued function of the horizontal coordinates. In the present investigation, semi‐implicit methods are extended to complex free‐surface flows that are governed by the full incompressible Navier–Stokes equations and are delimited by solid boundaries and arbitrarily shaped free‐surfaces. The primary dependent variables are the velocity components and the pressure. Finite difference equations for momentum, and a finite volume discretization for continuity, are derived in such a fashion that, after simple manipulation, the resulting pressure equation yields a well‐posed piecewise linear system from which both the pressure and the fluid volume within each computational cell are naturally derived. This system is efficiently solved by a nested Newton type iterative scheme, and the resulting fluid volumes are assured to be nonnegative and bounded from above by the available cell volumes. The time step size is not restricted by stability conditions dictated by surface wave speed, but can be freely chosen just to achieve the desired accuracy. Several examples illustrate the model applicability to a large range of complex free‐surface flows and demonstrate the effectiveness of the proposed algorithm. Copyright © 2013 John Wiley & Sons, Ltd.     

Global volume conservation in unsteady free surface flows with energy absorbing far‐end boundaries

A wave absorption filter for the far‐end boundary of semi‐infinite large reservoirs is developed for numerical simulation of unsteady free surface flows. Mathematical model is based on finite volume solution of the Navier–Stokes equations and depth‐integrated continuity equation to track the free surface. The Sommerfeld boundary condition is applied at the far‐end of the truncated computational domain. A dissipation zone is formed by applying artificial pressure on water surface to dissipate the kinetic energy of the outgoing waves. The computational scheme is tested to verify the conservation of total fluid volume in the domain for long simulation durations. Combination of the Sommerfeld boundary and dissipation zone can effectively minimize reflections and prevent cumulative changes in total fluid volume in the domain. Solitary wave, nonlinear periodic waves and irregular waves are simulated to illustrate the numerical developments. Earthquake excited surface waves and nonlinear hydrodynamic pressures in a dam–reservoir are computed. Copyright © 2009 John Wiley & Sons, Ltd.     

Parallel implementation of a non‐hydrostatic model for free surface flows with semi‐Lagrangian advection treatment

The parallel implementation of an unstructured‐grid, three‐dimensional, semi‐implicit finite difference and finite volume model for the free surface Navier–Stokes equations (UnTRIM ) is presented and discussed. The new developments are aimed to make the code available for high‐performance computing in order to address larger, complex problems in environmental free surface flows. The parallelization is based on the mesh partitioning method and message passing and has been achieved without negatively affecting any of the advantageous properties of the serial code, such as its robustness, accuracy and efficiency. The key issue is a new, autonomous parallel streamline backtracking algorithm, which allows using semi‐Lagrangian methods in decomposed meshes without compromising the scalability of the code. The implementation has been carefully verified not only with simple, abstract test cases illustrating the application domain of the code but also with advanced, high‐resolution models presently applied for research and engineering projects. The scheme performance and accuracy aspects are researched and discussed. Copyright © 2008 John Wiley & Sons, Ltd.     

Concurrent use of finite element and finite volume methods for shallow water flows in locally 1‐D channel networks

A numerical model is developed for shallow water equation in locally 1‐D channel networks. The model concurrently uses the standard Galerkin finite element method for the continuity equation and the finite volume method with an upwind scheme for the momentum equation. The surface gradient method is consistently employed. A minimum treatment is given for channel junctions so that application to multiply connected channels do not require any special consideration The model is capable of computing different types of transcritical flows, wet and dry flows, and flows with complex source terms. Standardized test problems and laboratory experimental data are used for verifying the model. Applicability of the models is validated in a multiply connected channel network draining hydromorphic farmlands located in a West African savanna, and Manning's roughness coefficient is identified, so that the steady solution is consistent with field observations. Unsteady simulation demonstrates that the model is capable of stably reproducing shifts of hydraulic jumps in flows of sub‐millimeter water depths. Copyright © 2011 John Wiley & Sons, Ltd.     

Boundary‐fitted non‐linear dispersive wave model for regions of arbitrary geometry

A vertically integrated non‐linear dispersive wave model is expressed in non‐orthogonal curvilinear co‐ordinate system for simulating shallow or deep water wave motions in regions of arbitrary geometry. Both dependent and independent variables are transformed so that an irregular physical domain is converted into a rectangular computational domain with contravariant velocities. Thus, the wall condition for enclosures surrounding a typical physical domain, such as a channel, port or harbor, is satisfied accurately and easily. The numerical scheme is based on staggered grid finite‐difference approximations, which result in implicit formulations for the momentum equations and semi‐explicit formulation for the continuity equation. Test cases of linear wave propagation in converging, diverging and circular channels are performed to check the reliability of model simulations against the analytical solutions. Cnoidal waves of different steepness values in a circular channel are also considered as examples to non‐linear wave propagation within curved walls. In closing, remarks concerning versatility and practical uses of the numerical model are made. Copyright © 2004 John Wiley & Sons, Ltd.     

An improved technique for solving two‐dimensional shallow water problems

A numerical model for solving the 2D shallow water equations is proposed herewith. This model is based on a finite volume technique in a generalized co‐ordinate system, coupled with a semi‐implicit splitting algorithm in which a Helmholtz equation is used for the surface elevation. Several benchmark problems have proven the good accuracy of this method in complex geometries. Nevertheless, several numerical perturbations were noted in the surface elevation. After finding the origin, a new numerical technique is suggested, to avoid these perturbations. Several severe tests are proposed to validate this technique. Copyright © 1999 John Wiley & Sons, Ltd.     

A semi‐Lagrangian level set method for incompressible Navier–Stokes equations with free surface

In this paper, we formulate a level set method in the framework of finite elements‐semi‐Lagrangian methods to compute the solution of the incompressible Navier–Stokes equations with free surface. In our formulation, we use a quasi‐monotone semi‐Lagrangian scheme, which is both unconditionally stable and essentially non oscillatory, to compute the advective terms in the Navier–Stokes equations, the transport equation and the equation of the reinitialization stage for the level set function. The method we propose is quite robust and flexible with regard to the mesh and the geometry of the domain, as well as the magnitude of the Reynolds number. We illustrate the performance of the method in several examples, which range from a benchmark problem to test the volume conservation property of the method to the flow past a NACA0012 foil at high Reynolds number. Copyright © 2005 John Wiley & Sons, Ltd.     

An immersed‐boundary finite‐volume method for direct simulation of flows with suspended paramagnetic particles

In the present study, we have proposed an immersed‐boundary finite‐volume method for the direct numerical simulation of flows with inertialess paramagnetic particles suspended in a nonmagnetic fluid under an external magnetic field without the need for any model such as the dipole–dipole interaction. In the proposed method, the magnetic field (or force) is described by the numerical solution of the Maxwell equation without current, where the smoothed representation technique is employed to tackle the discontinuity of magnetic permeability across the particle–fluid interface. The flow field, on the other hand, is described by the solution of the continuity and momentum equations, where the discrete‐forcing‐based immersed‐boundary method is employed to satisfy the no‐slip condition at the interface. To validate the method, we performed numerical simulations on the two‐dimensional motion of two and three paramagnetic particles in a nonmagnetic fluid subjected to an external uniform magnetic field and then compared the results with the existing finite‐element and semi‐analytical solutions. Comparison shows that the proposed method is robust in the direct simulation of such magnetic particulate flows. This method can be extended to more general flows without difficulty: three‐dimensional particulate flows, flows with a great number of particles, or flows under an arbitrary external magnetic field. Copyright © 2010 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.