Shallow flow simulation on dynamically adaptive cut cell quadtree grids

A computationally efficient, high‐resolution numerical model of shallow flow hydrodynamics is described, based on dynamically adaptive quadtree grids. The numerical model solves the two‐dimensional non‐linear shallow water equations by means of an explicit second‐order MUSCL‐Hancock Godunov‐type finite volume scheme. Interface fluxes are evaluated using an HLLC approximate Riemann solver. Cartesian cut cells are used to improve the fit to curved boundaries. A ghost‐cell immersed boundary method is used to update flow information in the smallest cut cells and overcome the time step restriction that would otherwise apply. The numerical model is validated through simulations of reflection of a surge wave at a wall, a low Froude number potential flow past a circular cylinder, and the shock‐like interaction between a bore and a circular cylinder. The computational efficiency is shown to be greatly improved compared with solutions on a uniform structured grid implemented with cut cells. Copyright © 2006 John Wiley & Sons, Ltd.     

Sparse element for shallow water wave equations on triangle meshes

Within the mixed FEM, the mini‐element that uses a bubble shape function for the solution of the shallow water wave equations on triangle meshes is simplified to a sparse element formulation. The new formulation has linear shape functions for water levels and constant shape functions for velocities inside each element. The suppression of decoupled spurious solutions is excellent with the new scheme. The linear dispersion relation of the new element has similar advantages as that of the wave equation scheme (generalised wave continuity scheme) proposed by Lynch and Gray. It is shown that the relation is monotonic over all wave numbers. In this paper, the time stepping scheme is included in the dispersion analysis. In case of a combined space–time staggering, the dispersion relation can be improved for the shortest waves. The sparse element is applied in the flow model Bubble that conserves mass exactly. At the same time, because of the limited number of degrees of freedom, the computational efficiency is high. The scheme is not restricted to orthogonal triangular meshes. Three test cases demonstrate the very good accuracy of the proposed scheme. The examples are the classical quarter annulus test case for the linearised shallow water equations, the hydraulic jump and the tide in the Elbe river mouth. Copyright © 2013 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.     

Improved multislope MUSCL reconstruction on unstructured grids for shallow water equations

In shallow water flow and transport modeling, the monotonic upstream‐centered scheme for conservation laws (MUSCL) is widely used to extend the original Godunov scheme to second‐order accuracy. The most important step in MUSCL‐type schemes is MUSCL reconstruction, which calculate‐extrapolates the values of independent variables from the cell center to the edge. The monotonicity of the scheme is preserved with the help of slope limiters that prevent the occurrence of new extrema during reconstruction. On structured grids, the calculation of the slope is straightforward and usually based on a 2‐point stencil that uses the cell centers of the neighbor cell and the so‐called far‐neighbor cell of the edge under consideration. On unstructured grids, the correct choice for the upwind slope becomes nontrivial. In this work, 2 novel total variation diminishing schemes are developed based on different techniques for calculating the upwind slope and the downwind slope. An additional treatment that stabilizes the scheme is discussed. The proposed techniques are compared to 2 existing MUSCL reconstruction techniques, and a detailed discussion of the results is given. It is shown that the proposed MUSCL reconstruction schemes obtain more accurate results with less numerical diffusion and higher efficiency.     

Exact pressure integrations on submerged bodies in waves using a quadtree adaptive mesh algorithm

The development of an adaptive free surface, mesh cutting, methodology, in order to analytically integrate pressures on varying wet parts of partially submerged surfaces in the presence of waves, is presented. Given a function of free‐surface elevation, the algorithm checks for the intersection of the body with the free surface and, based on user‐defined parameters, modifies the initial mesh, by subdividing the elements where necessary and eliminating others, via a quadtree approach. Redundant sub‐divisions, generated in the quad‐division process, are partially eliminated, but the quadrilateral nature of the elements is always kept. The free‐surface function must be single‐valued and its definition domain simply connected. Hydrostatic and Froude–Krylov forces are computed exactly on each panel by means of analytical formulations, which are derived and presented, based on the theory of linear gravity waves and from applying Green's theorem. Copyright © 2014 John Wiley & Sons, Ltd.     

The Lagrange–Galerkin method for the two‐dimensional shallow water equations on adaptive grids

The weak Lagrange–Galerkin finite element method for the two‐dimensional shallow water equations on adaptive unstructured grids is presented. The equations are written in conservation form and the domains are discretized using triangular elements. Lagrangian methods integrate the governing equations along the characteristic curves, thus being well suited for resolving the non‐linearities introduced by the advection operator of the fluid dynamics equations. An additional fortuitous consequence of using Lagrangian methods is that the resulting spatial operator is self‐adjoint, thereby justifying the use of a Galerkin formulation; this formulation has been proven to be optimal for such differential operators. The weak Lagrange–Galerkin method automatically takes into account the dilation of the control volume, thereby resulting in a conservative scheme. The use of linear triangular elements permits the construction of accurate (by virtue of the second‐order spatial and temporal accuracies of the scheme) and efficient (by virtue of the less stringent Courant–Friedrich–Lewy (CFL) condition of Lagrangian methods) schemes on adaptive unstructured triangular grids. Lagrangian methods are natural candidates for use with adaptive unstructured grids because the resolution of the grid can be increased without having to decrease the time step in order to satisfy stability. An advancing front adaptive unstructured triangular mesh generator is presented. The highlight of this algorithm is that the weak Lagrange–Galerkin method is used to project the conservation variables from the old mesh onto the newly adapted mesh. In addition, two new schemes for computing the characteristic curves are presented: a composite mid‐point rule and a general family of Runge–Kutta schemes. Results for the two‐dimensional advection equation with and without time‐dependent velocity fields are illustrated to confirm the accuracy of the particle trajectories. Results for the two‐dimensional shallow water equations on a non‐linear soliton wave are presented to illustrate the power and flexibility of this strategy. Copyright © 2000 John Wiley & Sons, Ltd.     

Godunov-type solution of the shallow water equations on adaptive unstructured triangular grids

A Godunov-type upwind finite volume solver of the non-linear shallow water equations is described. The shallow water equations are expressed in a hyperbolic conservation law formulation for application to cases where the bed topography is spatially variable. Inviscid fluxes at cell interfaces are computed using Roe's approximate Riemann solver. Second-order accurate spatial calculations of the fluxes are achieved by enhancing the polynomial approximation of the gradients of conserved variables within each cell. Numerical oscillations are curbed by means of a non-linear slope limiter. Time integration is second-order accurate and implicit. The numerical model is based on dynamically adaptive unstructured triangular grids. Test cases include an oblique hydraulic jump, jet-forced flow in a flat-bottomed circular reservoir, wind-induced circulation in a circular basin of non-uniform bed topography and the collapse of a circular dam. The model is found to give accurate results in comparison with published analytical and alternative numerical solutions. Dynamic grid adaptation and the use of a second-order implicit time integration scheme are found to enhance the computational efficiency of the model.     

一个求解浅水波方程的时域分段展开算法

发展了一种时域分段展开自适应方法求解一维非线性浅水波方程。通过时域分段展开,将一个非线性的时空耦合初边值问题转化为一系列的线性空间边值问题,并采用有限元方法递推求解;通过展开阶数的递进,实现了分段时域的自适应计算,当不同步长时可保持稳定的计算精度。研究结果表明,当步长较大而Heun’s法、四阶Runge-Kutta法不能得到合理结果时,本文算法仍能保证足够的计算精度。     

Finite‐volume‐type VOF method on dynamically adaptive quadtree grids

This paper describes a finite‐volume volume‐of‐fluid (VOF) method for simulating viscous free surface flows on dynamically adaptive quadtree grids. The scheme is computationally efficient in that it provides relatively fine grid resolution at the gas–liquid interface and coarse grid density in regions where flow variable gradients are small. Special interpolations are used to ensure volume flux conservation where differently sized neighbour cells occur. The numerical model is validated for advection of dyed fluid in unidirectional and rotating flows, and for two‐dimensional viscous sloshing in a rectangular tank. Copyright © 2004 John Wiley & Sons, Ltd.     

Contradiction between the C‐property and mass conservation in adaptive grid based shallow flow models: cause and solution

When performing shallow flow simulations on adaptive grids, the C‐property (i.e. conservation property) and the mass conservation may not be simultaneously preserved, that is, either C‐property or mass conservation is likely to be violated following grid refining or coarsening. The cause of such a contradiction is analyzed in detail in this work, which essentially links to the reconstruction of bed and flow information in those newly created cells during grid adaptation. An effective approach is subsequently proposed to resolve the contradiction by locally modifying the bed elevation in certain problematic cells when reconstructing flow information using linear interpolation as part of the grid adaptation procedure. Copyright © 2015 John Wiley & Sons, Ltd.     

Parallel adaptive refinement for unsteady flow calculations on 3D unstructured grids

A parallel adaptive refinement algorithm for three‐dimensional unstructured grids is presented. The algorithm is based on an hierarchical h‐refinement/derefinement scheme for tetrahedral elements.The algorithm has been fully parallelized for shared‐memory platforms via a domain decomposition of the mesh at the algebraic level. The effectiveness of the procedure is demonstrated with applications which involve unsteady compressible fluid flow. A parallel speedup study of the algorithm also is included. Published in 2004 by John Wiley & Sons, Ltd.     

Solution of shallow water equations using fully adaptive multiscale schemes

The concept of fully adaptive multiscale finite volume methods has been developed to increase spatial resolution and to reduce computational costs of numerical simulations. Here grid adaptation is performed by means of a multiscale analysis based on biorthogonal wavelets. In order to update the solution in time we use a local time stepping strategy that has been recently developed for hyperbolic conservation laws. The adaptive multiresolution scheme is now applied to two‐dimensional shallow water equations with source terms. The efficiency of the scheme is demonstrated on several problems with a general geometry, including circular damp breaks, oblique hydraulic jump, supercritical channel flows encountering sudden change in cross‐section, and, finally, the bore wave and its interactions. Copyright © 2005 John Wiley & Sons, Ltd.     

Simulation of shallow flows over variable topographies using unstructured grids

Simulation of shallow flows over variable topographies is a challenging case for most available shock‐capturing schemes. This problem arises because the source terms and flux gradients are not balanced in the numerical computations. Treatments for this problem generally work well on structured grids, but they are usually too expensive, and most of them are not directly applicable to unstructured grids. In this paper we propose two efficient methods to treat the source terms without upwinding and to satisfy the compatibility condition on unstructured grids. In the first method, the calculation of the bed slope source term is performed by employing a compatible approximation of water depth at the cell interfaces. In the second one, different components of the bed slope term are considered separately and a compatible discretization of the components is proposed. The present treatments are applicable for most schemes including the Roe's method without changing the performance of the original scheme for smooth topographies. Copyright © 2006 John Wiley & Sons, Ltd.     

A simplified adaptive Cartesian grid system for solving the 2D shallow water equations

This paper presents a new simplified grid system that provides local refinement and dynamic adaptation for solving the 2D shallow water equations (SWEs). Local refinement is realized by simply specifying different subdivision levels to the cells on a background uniform coarse grid that covers the computational domain. On such a non‐uniform grid, the structured property of a regular Cartesian mesh is maintained and neighbor information is determined by simple algebraic relationships, i.e. data structure becomes unnecessary. Dynamic grid adaptation is achieved by changing the subdivision level of a background cell. Therefore, grid generation and adaptation is greatly simplified and straightforward to implement. The new adaptive grid‐based SWE solver is tested by applying it to simulate three idealized test cases and promising results are obtained. The new grid system offers a simplified alternative to the existing approaches for providing adaptive mesh refinement in computational fluid dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

A method for analysing nesting techniques for the linearized shallow water equations

A method for analysing different nesting techniques for the linearized shallow water equations is presented. The problem is formulated as an eigenvector–eigenvalue problem. A necessary condition for stability is that the spectral radius of the propagation matrix is less than or equal to one. Two test cases are presented. The first test case is analysed, and effects of enforcing volume conservation and nudging in time are studied. A nesting technique is found that causes no growth of any eigenvectors for reasonable time steps. This nesting technique is then used on both test cases, and results are compared to an everywhere refined model and a coarse grid model. Copyright © 2002 John Wiley & Sons, Ltd.     

A fast nested multi‐grid viscous flow solver for adaptive Cartesian/Quad grids

A nested multi‐grid solution algorithm has been developed for an adaptive Cartesian/Quad grid viscous flow solver. Body‐fitted adaptive Quad (quadrilateral) grids are generated around solid bodies through ‘surface extrusion’. The Quad grids are then overlapped with an adaptive Cartesian grid. Quadtree data structures are employed to record both the Quad and Cartesian grids. The Cartesian grid is generated through recursive sub‐division of a single root, whereas the Quad grids start from multiple roots—a forest of Quadtrees, representing the coarsest possible Quad grids. Cell‐cutting is performed at the Cartesian/Quad grid interface to merge the Cartesian and Quad grids into a single unstructured grid with arbitrary cell topologies (i.e., arbitrary polygons). Because of the hierarchical nature of the data structure, many levels of coarse grids have already been built in. The coarsening of the unstructured grid is based on the Quadtree data structure through reverse tree traversal. Issues arising from grid coarsening are discussed and solutions are developed. The flow solver is based on a cell‐centered finite volume discretization, Roe's flux splitting, a least‐squares linear reconstruction, and a differentiable limiter developed by Venkatakrishnan in a modified form. A local time stepping scheme is used to handle very small cut cells produced in cell‐cutting. Several cycling strategies, such as the saw‐tooth, W‐ and V‐cycles, have been studies. The V‐cycle has been found to be the most efficient. In general, the multi‐grid solution algorithm has been shown to greatly speed up convergence to steady state—by one to two orders. Copyright © 2000 John Wiley & Sons, Ltd.     

A mass‐conservative staggered immersed boundary model for solving the shallow water equations on complex geometries

In this work, an approach is proposed for solving the 3D shallow water equations with embedded boundaries that are not aligned with the underlying horizontal Cartesian grid. A hybrid cut‐cell/ghost‐cell method is used together with a direction‐splitting implicit solver: Ghost cells are used for the momentum equations in order to prescribe the correct boundary condition at the immersed boundary, while cut cells are used in the continuity equation in order to conserve mass. The resulting scheme is robust, does not suffer any time step limitation for small cut cells, and conserves fluid mass up to machine precision. Moreover, the solver displays a second‐order spatial accuracy, both globally and locally. Comparisons with analytical solutions and reference numerical solutions on curvilinear grids confirm the quality of the method. Copyright © 2015 John Wiley & Sons, Ltd.     

Development and testing of a simple 2D finite volume model of sub‐critical shallow water flow

Many environmental applications of shallow water flow modelling can be characterized as only slowly varying and everywhere sub‐critical. A simplified finite volume model is therefore developed that is capable of describing pertinent shallow water flow processes more efficiently than the usual Godunov/ Riemann characteristics approaches. The model is tested against a number of analytical and numerical solutions to the governing equations. The model reproduces accurately flow round a circular bend, flow over topography, flow up an initially dry beach and floodwave propagation down a meandering river reach, with mass conservative solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Coupling between shallow water and solute flow equations: analysis and management of source terms in 2D

A two‐dimensional model for the simulation of solute transport by convection and diffusion into shallow water flow over variable bottom is presented. It is based on a finite volume method over triangular unstructured grids. A first order upwind technique is applied to solve the flux terms in both the flow and solute equations and the bed slope source terms and a centred discretization is applied to the diffusion and friction terms. The convenience of considering the fully coupled system of equations is indicated and the methodology is well explained. Three options are suggested and compared in order to deal with the diffusion terms. Some comparisons are carried out in order to show the performance in terms of accuracy and computational effort of the different options. Copyright © 2005 John Wiley & Sons, Ltd.     

Boiling flow simulations on adaptive octree grids

Boiling flow simulations are conducted on adaptive octree grids. A phase change model consistent with the mixture formulation, in conjunction with the Volume-of-Fluid (VOF) model, is used to track the liquid–vapor interface. Test cases including Rayleigh Taylor instability and bubble growth in a uniform superheat are conducted to validate the phase change model on adaptive grids. The validated model is then used to conduct film boiling simulations on both two-dimensional and three-dimensional adaptive grids. The average wall Nusselt number agrees well with the widely accepted correlations of Berenson (1961) and Klimenko (1981) and Klimenko and Shelepen (1982) for film boiling on a horizontal surface. For the test cases presented, the efficiency of the adaptive technique, as measured by the adaptive mesh refinement (AMR) efficiency, is mostly in the range of 50–80%. Although this efficiency is a function of the nature and dimensionality of the problem, this range of efficiency is comparable to those obtained in the simulations of primary jet atomization conducted by Fuster et al. (2009). This work opens the prospect of conducting more realistic (three-dimensional) multi-modal boiling flow simulations, and problems of similar complexity, in an efficient manner.