A conservative semi‐implicit method for coupled surface–subsurface flows in regional scale

A semi‐implicit method for coupled surface–subsurface flows in regional scale is proposed and analyzed. The flow domain is assumed to have a small vertical scale as compared with the horizontal extents. Thus, after hydrostatic approximation, the simplified governing equations are derived from the Reynolds averaged Navier–Stokes equations for the surface flow and from the Darcy's law for the subsurface flow. A conservative free‐surface equation is derived from a vertical integral of the incompressibility condition and extends to the whole water column including both, the surface and the subsurface, wet domains. Numerically, the horizontal domain is covered by an unstructured orthogonal grid that may include subgrid specifications. Along the vertical direction a simple z‐layer discretization is adopted. Semi‐implicit finite difference equations for velocities and a finite volume approximation for the free‐surface equation are derived in such a fashion that, after simple manipulation, the resulting discrete free‐surface equation yields a single, well‐posed, mildly nonlinear system. This system is efficiently solved by a nested Newton‐type iterative method that yields simultaneously the pressure and a non‐negative fluid volume throughout the computational grid. The time‐step size is not restricted by stability conditions dictated by friction or surface wave speed. The resulting algorithm is simple, extremely efficient, and very accurate. Exact mass conservation is assured also in presence of wetting and drying dynamics, in pressurized flow conditions, and during free‐surface transition through the interface. A few examples illustrate the model applicability and demonstrate the effectiveness of the proposed algorithm. Copyright © 2015 John Wiley & Sons, Ltd.     

A semi‐implicit numerical model for urban drainage systems

In this paper, a semi‐implicit numerical model for one‐dimensional urban drainage networks is formulated in such a fashion as to intrinsically account for arbitrary cross sections, for the occurrence of dry areas, for free surface, and for pressurized flows. The governing differential equations are discretized with a consistent mass conservative scheme that naturally applies to all flow regimes. The resulting mildly nonlinear system, at every time step, is efficiently solved with a converging, properly devised, nested Newton‐type algorithm. It will be shown that with the proposed semi‐implicit model, high accuracy can be achieved at a moderate computational cost. Copyright © 2013 John Wiley & Sons, Ltd.     

A two‐dimensional numerical scheme of dry/wet fronts for the Saint‐Venant system of shallow water equations

We propose a new two‐dimensional numerical scheme to solve the Saint‐Venant system of shallow water equations in the presence of partially flooded cells. Our method is well balanced, positivity preserving, and handles dry states. The latter is ensured by using the draining time step technique in the time integration process, which guarantees non‐negative water depths. Unlike previous schemes, our technique does not generate high velocities at the dry/wet boundaries, which are responsible for small time step sizes and slow simulation runs. We prove that the new scheme preserves ‘lake at rest’ steady states and guarantees the positivity of the computed fluid depth in the partially flooded cells. We test the new scheme, along with another recent scheme from the literature, against the analytical solution for a parabolic basin and show the improved simulation performance of the new scheme for two real‐world scenarios. Copyright © 2014 John Wiley & Sons, Ltd.     

A robust and well‐balanced scheme for the 2D Saint‐Venant system on unstructured meshes with friction source term

In the following lines, we propose a numerical scheme for the shallow‐water system supplemented by topography and friction source terms, in a 2D unstructured context. This work proposes an improved version of the well‐balanced and robust numerical model recently introduced by Duran et al. (J. Comp. Phys., 235 , 565–586, 2013) for the pre‐balanced shallow‐water equations, accounting for varying topography. The present work aims at relaxing the robustness condition and includes a friction term. To this purpose, the scheme is modified using a recent method, entirely based on a modified Riemann solver. This approach preserves the robustness and well‐balanced properties of the original scheme and prevents unstable computations in the presence of low water depths. A series of numerical experiments are devoted to highlighting the performances of the resulting scheme. Simulations involving dry areas, complex geometry and topography are proposed to validate the stability of the numerical model in the neighbourhood of wet/dry transitions. Copyright © 2015 John Wiley & Sons, Ltd.     

Development of temporal and spatial high‐order schemes for two‐fluid seven‐equation two‐pressure model and its applications in two‐phase flow benchmark problems

Current existing main nuclear thermal‐hydraulics (T‐H) system analysis codes, such as RALAP5, TRACE, and CATHARE, play a crucial role in the nuclear engineering field for the design and safety analysis of nuclear reactor systems. However, two‐fluid model used in these T‐H system analysis codes is ill posed, easily leading to numerical oscillations, and the classical first‐order methods for temporal and special discretization are widely employed for numerical simulations, yielding excessive numerical diffusion. Two‐fluid seven‐equation two‐pressure model is of particular interest due to the inherent well‐posed advantage. Moreover, high‐order accuracy schemes have also attracted great attention to overcome the challenge of serious numerical diffusion induced by low‐order time and space schemes for accurately simulating nuclear T‐H problems. In this paper, the semi‐implicit solution algorithm with high‐order accuracy in space and time is developed for this well‐posed two‐fluid model and the robustness and accuracy are verified and assessed against several important two‐phase flow benchmark tests in the nuclear engineering T‐H field, which include two linear advection problems, the oscillation problem of the liquid column, the Ransom water faucet problem, the reversed water faucet problem, and the two‐phase shock tube problem. The following conclusions are achieved. (1) The proposed semi‐implicit solution algorithm is robust in solving two‐phase flows, even for fast transients and discontinuous solutions. (2) High‐order schemes in both time and space could prevent excessive numerical diffusion effectively and the numerical simulation results are more accurate than those of first‐order time and space schemes, which demonstrates the advantage of using high‐order schemes.     

An implicit two‐dimensional non‐hydrostatic model for free‐surface flows

An implicit finite volume model in sigma coordinate system is developed to simulate two‐dimensional (2D) vertical free surface flows, deploying a non‐hydrostatic pressure distribution. The algorithm is based on a projection method which solves the complete 2D Navier–Stokes equations in two steps. First the pressure term in the momentum equations is excluded and the resultant advection–diffusion equations are solved. In the second step the continuity and the momentum equation with only the pressure terms are solved to give a block tri‐diagonal system of equation with pressure as the unknown. This system can be solved by a direct matrix solver without iteration. A new implicit treatment of non‐hydrostatic pressure, similar to the lower layers is applied to the top layer which makes the model free of any hydrostatic pressure assumption all through the water column. This treatment enables the model to evaluate both free surface elevation and wave celerity more accurately. A series of numerical tests including free‐surface flows with significant vertical accelerations and nonlinear behaviour in shoaling zone are performed. Comparison between numerical results, analytical solutions and experimental data demonstrates a satisfactory performance. Copyright © 2006 John Wiley & Sons, Ltd.     

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

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

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.     

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.     

Finite‐volume multi‐stage schemes for shallow‐water flow simulations

A finite‐volume multi‐stage (FMUSTA) scheme is proposed for simulating the free‐surface shallow‐water flows with the hydraulic shocks. On the basis of the multi‐stage (MUSTA) method, the original Riemann problem is transformed to an independent MUSTA mesh. The local Lax–Friedrichs scheme is then adopted for solving the solution of the Riemann problem at the cell interface on the MUSTA mesh. The resulting first‐order monotonic FMUSTA scheme, which does not require the use of the eigenstructure and the special treatment of entropy fixes, has the generality as well as simplicity. In order to achieve the high‐resolution property, the monotonic upstream schemes for conservation laws (MUSCL) method are used. For modeling shallow‐water flows with source terms, the surface gradient method (SGM) is adopted. The proposed schemes are verified using the simulations of six shallow‐water problems, including the 1D idealized dam breaking, the steady transcritical flow over a hump, the 2D oblique hydraulic jump, the circular dam breaking and two dam‐break experiments. The simulated results by the proposed schemes are in satisfactory agreement with the exact solutions and experimental data. It is demonstrated that the proposed FMUSTA schemes have superior overall numerical accuracy among the schemes tested such as the commonly adopted Roe and HLL schemes. Copyright © 2007 John Wiley & Sons, Ltd.     

A coupled surface‐subsurface model for hydrostatic flows under saturated and variably saturated conditions

In this paper, the governing differential equations for hydrostatic surface‐subsurface flows are derived from the Richards and from the Navier‐Stokes equations. A vertically integrated continuity equation is formulated to account for both surface and subsurface flows under saturated and variable saturated conditions. Numerically, the horizontal domain is covered by an unstructured orthogonal grid that may include subgrid specifications. Along the vertical direction, a simple z‐layer discretization is adopted. Semi‐implicit finite difference equations for velocities, and a finite volume approximation for the vertically integrated continuity equation, are derived in such a fashion that, after simple manipulation, the resulting discrete pressure equation can be assembled into a single, two‐dimensional, mildly nonlinear system. This system is solved by a nested Newton‐type method, which yields simultaneously the (hydrostatic) pressure and a nonnegative fluid volume throughout the computational grid. The resulting algorithm is relatively simple, extremely efficient, and very accurate. Stability, convergence, and exact mass conservation are assured throughout also in presence of wetting and drying, in variable saturated conditions, and during flow transition through the soil interface. A few examples illustrate the model applicability and demonstrate the effectiveness of the proposed algorithm.     

The solution of two‐dimensional free‐surface problems using automatic mesh generation

A new method is described for the iterative solution of two‐dimensional free‐surface problems, with arbitrary initial geometries, in which the interior of the domain is represented by an unstructured, triangular Eulerian mesh and the free surface is represented directly by the piecewise‐quadratic edges of the isoparametric quadratic‐velocity, linear‐pressure Taylor–Hood elements. At each time step, the motion of the free surface is computed explicitly using the current velocity field and, once the new free‐surface location has been found, the interior nodes of the mesh are repositioned using a continuous deformation model that preserves the original connectivity. In the event that the interior of the domain must be completely remeshed, a standard Delaunay triangulation algorithm is used, which leaves the initial boundary discretisation unchanged. The algorithm is validated via the benchmark viscous flow problem of the coalescence of two infinite cylinders of equal radius, in which the motion is due entirely to the action of capillary forces on the free surface. This problem has been selected for a variety of reasons: the initial and final (steady state) geometries differ considerably; in the passage from the former to the latter, large free‐surface curvatures—requiring accurate modelling—are encountered; an analytical solution is known for the location of the free surface; there exists a large body of literature on alternative numerical simulations. A novel feature of the present work is its geometric generality and robustness; it does not require a priori knowledge of either the evolving domain geometry or the solution contained therein. Copyright © 1999 John Wiley & Sons, Ltd.     

Simulation of dam‐ and dyke‐break hydrodynamics on dynamically adaptive quadtree grids

Flooding due to the failure of a dam or dyke has potentially disastrous consequences. This paper presents a Godunov‐type finite volume solver of the shallow water equations based on dynamically adaptive quadtree grids. The Harten, Lax and van Leer approximate Riemann solver with the Contact wave restored (HLLC) scheme is used to evaluate interface fluxes in both wet‐ and dry‐bed applications. The numerical model is validated against results from alternative numerical models for idealized circular and rectangular dam breaks. Close agreement is achieved with experimental measurements from the CADAM dam break test and data from a laboratory dyke break undertaken at Delft University of Technology. Copyright © 2004 John Wiley Sons, Ltd.     

Well‐balanced positivity preserving central‐upwind scheme for the shallow water system with friction terms

Shallow water models are widely used to describe and study free‐surface water flow. While in some practical applications the bottom friction does not have much influence on the solutions, there are still many applications, where the bottom friction is important. In particular, the friction terms will play a significant role when the depth of the water is very small. In this paper, we study shallow water equations with friction terms and develop a semi‐discrete second‐order central‐upwind scheme that is capable of exactly preserving physically relevant steady states and maintaining the positivity of the water depth. The presence of the friction terms increases the level of complexity in numerical simulations as the underlying semi‐discrete system becomes stiff when the water depth is small. We therefore implement an efficient semi‐implicit Runge‐Kutta time integration method that sustains the well‐balanced and sign preserving properties of the semi‐discrete scheme. We test the designed method on a number of one‐dimensional and two‐dimensional examples that demonstrate robustness and high resolution of the proposed numerical approach. The data in the last numerical example correspond to the laboratory experiments reported in [L. Cea, M. Garrido, and J. Puertas, Journal of Hydrology, 382 (2010), pp. 88–102], designed to mimic the rain water drainage in urban areas containing houses. Since the rain water depth is typically several orders of magnitude smaller than the height of the houses, we develop a special technique, which helps to achieve a remarkable agreement between the numerical and experimental results. Copyright © 2015 John Wiley & Sons, Ltd.     

A mass conservative well‐balanced reconstruction at wet/dry interfaces for the Godunov‐type shallow water model

This paper presents a novel mass conservative, positivity preserving wetting and drying treatment for Godunov‐type shallow water models with second‐order bed elevation discretization. The novel method allows to compute water depths equal to machine accuracy without any restrictions on the time step or any threshold that defines whether the finite volume cell is considered to be wet or dry. The resulting scheme is second‐order accurate in space and keeps the C‐property condition at fully flooded area and also at the wet/dry interface. For the time integration, a second‐order accurate Runge–Kutta method is used. The method is tested in two well‐known computational benchmarks for which an analytical solution can be derived, a C‐property benchmark and in an additional example where the experimental results are reproduced. Overall, the presented scheme shows very good agreement with the reference solutions. The method can also be used in the discontinuous Galerkin method. Copyright © 2016 John Wiley & Sons, Ltd.     

A kinetic interpretation of the section‐averaged Saint‐Venant system for natural river hydraulics

The classical Saint‐Venant system is well suited for the modeling of dam breaks, hydraulic jumps, reservoir emptying, flooding etc. For many applications, the extension of the Saint‐Venant system to the case of non‐rectangular channels is necessary and this section‐averaged Saint‐Venant system exhibits additional source terms. The main difficulty of these equations consists of the discretization of these source terms. In this paper we propose a kinetic interpretation for the section averaged Saint‐Venant system and derive an associated numerical scheme. The numerical scheme—2nd order in space and time—preserves the positivity of the water height, and is well‐balanced. Numerical results including comparisons with analytic and experimental test problems illustrate the accuracy and the robustness of the numerical algorithm. Copyright © 2010 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.     

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.     

Three‐dimensional numerical modelling of free surface flows with non‐hydrostatic pressure

A three‐dimensional numerical model is developed for incompressible free surface flows. The model is based on the unsteady Reynolds‐averaged Navier–Stokes equations with a non‐hydrostatic pressure distribution being incorporated in the model. The governing equations are solved in the conventional sigma co‐ordinate system, with a semi‐implicit time discretization. A fractional step method is used to enable the pressure to be decomposed into its hydrostatic and hydrodynamic components. At every time step one five‐diagonal system of equations is solved to compute the water elevations and then the hydrodynamic pressure is determined from a pressure Poisson equation. The model is applied to three examples to simulate unsteady free surface flows where non‐hydrostatic pressures have a considerable effect on the velocity field. Emphasis is focused on applying the model to wave problems. Two of the examples are about modelling small amplitude waves where the hydrostatic approximation and long wave theory are not valid. The other example is the wind‐induced circulation in a closed basin. The numerical solutions are compared with the available analytical solutions for small amplitude wave theory and very good agreement is obtained. Copyright © 2002 John Wiley & Sons, Ltd.     

Method of fundamental solutions for partial‐slip fibrous filtration flows

In this study a Stokeslet‐based method of fundamental solutions (MFS) for two‐dimensional low Reynolds number partial‐slip flows has been developed. First, the flow past an infinitely long cylinder is selected as a benchmark. The numerical accuracy is investigated in terms of the location and the number of the Stokeslets. The benchmark study shows that the numerical accuracy increases when the Stokeslets are submerged deeper beneath the cylinder surface, as long as the formed linear system remains numerically solvable. The maximum submergence depth increases with the decrease in the number of Stokeslets. As a result, the numerical accuracy does not deteriorate with the dramatic decrease in the number of Stokeslets. A relatively small number of Stokeslets with a substantial submergence depth is thus chosen for modeling fibrous filtration flows. The developed methodology is further examined by application to Taylor–Couette flows. A good agreement between the numerical and analytical results is observed for no‐slip and partial‐slip boundary conditions. Next, the flow about a representative set of infinitely long cylindrical fibers confined between two planar walls is considered to represent the fibrous filter flow. The obtained flowfield and pressure drop agree very well with the experimental data for this setup of fibers. The developed MFS with submerged Stokeslets is then applied to partial‐slip flows about fibers to investigate the slip effect at fiber–fluid interface on the pressure drop. The numerical results compare qualitatively with the analytical solution available for the limit case of infinite number of fibers. Copyright © 2008 John Wiley & Sons, Ltd.