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.     

Simple treatment of non‐aligned boundaries in a Cartesian grid shallow flow model

A simple method is proposed for treating curved or irregular boundaries in Cartesian grid shallow flow models. It directly evaluates fictional values in ‘ghost’ cells adjacent to boundary cells and requires no interpolation or generation of cut cells. The boundary treatment is implemented in a dynamically adaptive quadtree grid‐based solver of the hyperbolic shallow water equations and validated against several test cases with analytical or alternative numerical solutions. The method is easy to code, accurate, and demonstrably effective in dealing with irregular computational domains in shallow flow simulations. Results are presented for still water in a basin of complicated geometry, steady hydraulic jump in an open channel with a converging sidewall, wind‐induced circulation in a circular shallow lake, and shock wave diffraction in a channel containing a contraction and expansion. Copyright © 2007 John Wiley & Sons, Ltd.     

A structured but non‐uniform Cartesian grid‐based model for the shallow water equations

In large‐scale shallow flow simulations, local high‐resolution predictions are often required in order to reduce the computational cost without losing the accuracy of the solution. This is normally achieved by solving the governing equations on grids refined only to those areas of interest. Grids with varying resolution can be generated by different approaches, e.g. nesting methods, patching algorithms and adaptive unstructured or quadtree gridding techniques. This work presents a new structured but non‐uniform Cartesian grid system as an alternative to the existing approaches to provide local high‐resolution mesh. On generating a structured but non‐uniform Cartesian grid, the whole computational domain is first discretized using a coarse background grid. Local refinement is then achieved by directly allocating a specific subdivision level to each background grid cell. The neighbour information is specified by simple mathematical relationships and no explicit storage is needed. Hence, the structured property of the uniform grid is maintained. After employing some simple interpolation formulae, the governing shallow water equations are solved using a second‐order finite volume Godunov‐type scheme in a similar way as that on a uniform grid. Copyright © 2010 John Wiley & Sons, Ltd.     

Simulation of viscous water column collapse using adapting hierarchical grids

An adaptive hierarchical grid‐based method for predicting complex free surface flows is used to simulate collapse of a water column. Adapting quadtree grids are combined with a high‐resolution interface‐capturing approach and pressure‐based coupling of the Navier–Stokes equations. The Navier–Stokes flow solution scheme is verified for simulation of flow in a lid‐driven cavity at Re=1000. Two approaches to the coupling of the Navier–Stokes equations are investigated as are alternative face velocity and hanging node interpolations. Collapse of a water column as well as collapse of a water column and its subsequent interaction with an obstacle are simulated. The calculations are made on uniform and adapting quadtree grids, and the accuracy of the quadtree calculations is shown to be the same as those made on the equivalent uniform grids. Results are in excellent agreement with experimental and other numerical data. A sharp interface is maintained at the free surface. The new adapting quadtree‐based method achieves a considerable saving in the size of the computational grid and CPU time in comparison with calculations made on equivalent uniform grids. Copyright © 2005 John Wiley & Sons, Ltd.     

多尺度嵌入式离散裂缝模型模拟方法

天然裂缝性油藏和人工压裂油藏内裂缝形态多样,分布复杂,传统的离散裂缝模型将裂缝作为基岩网格的边界,采用非结构化网格进行网格划分,其划分过程复杂,计算量大。嵌入式离散裂缝模型划分网格时不需要考虑油藏内的裂缝形态,只需对基岩系统进行简单的网格剖分,可以大大降低网格划分的复杂度,从而提高计算效率。然而,在油藏级别的数值模拟和人工压裂裂缝下的产能分析中,仍然存在计算量巨大、模拟时间过长的问题。本文提出嵌入式离散裂缝模型的多尺度数值计算格式,使用多尺度模拟有限差分法研究嵌入式离散裂缝模型渗流问题。通过在粗网格上求解局部流动问题计算多尺度基函数,多尺度基函数可以捕捉裂缝与基岩间的相互关系,反映单元内的非均质性,因此该方法既有传统尺度升级法的计算效率,又可以保证计算精度,数值结果表明这是一种有效的裂缝性油藏数值模拟方法。     

An unstructured grid,three‐dimensional model based on the shallow water equations

A semi‐implicit finite difference model based on the three‐dimensional shallow water equations is modified to use unstructured grids. There are obvious advantages in using unstructured grids in problems with a complicated geometry. In this development, the concept of unstructured orthogonal grids is introduced and applied to this model. The governing differential equations are discretized by means of a semi‐implicit algorithm that is robust, stable and very efficient. The resulting model is relatively simple, conserves mass, can fit complicated boundaries and yet is sufficiently flexible to permit local mesh refinements in areas of interest. Moreover, the simulation of the flooding and drying is included in a natural and straightforward manner. These features are illustrated by a test case for studies of convergence rates and by examples of flooding on a river plain and flow in a shallow estuary. Copyright © 2000 John Wiley & Sons, Ltd.     

A Comparison Study Between an Adaptive Quadtree Grid and Uniform Grid Upscaling for Reservoir Simulation

An adaptive quadtree grid generation algorithm is developed and applied for tracer and multiphase flow in channelized heterogeneous porous media. Adaptivity was guided using two different approaches. In the first approach, wavelet transformation was used to generate a refinement field based on permeability variations. The second approach uses flow information based on the solution of an initial-time fine-scale problem. The resulting grids were compared with uniform grid upscaling. For uniform upscaling, two commonly applied methods were used: renormalization upscaling and local-global upscaling. The velocities obtained by adaptive grid and uniformly upscaled grids, were downscaled. This procedure allows us to separate the upscaling errors, on adaptive and uniform grids, from the numerical dispersion errors resulting from solving the saturation equation on a coarse grid. The simulation results obtained by solving on flow-based adaptive quadtree grids for the case of a single phase flow show reasonable agreement with more computationally demanding fine-scale models and local-global upscaled models. For the multiphase case, the agreement is less evident, especially in piston-like displacement cases with sharp frontal movement. In such cases a non-iterative transmissibility upscaling procedure for adaptive grid is shown to significantly reduce the errors and make the adaptive grid comparable to iterative local-global upscaling. Furthermore, existence of barriers in a porous medium complicates both upscaling and grid adaptivity. This issue is addressed by adapting the grid using a combination of flow information and a permeability based heuristic criterion.     

Adaptive Q‐tree Godunov‐type scheme for shallow water equations

This paper presents details of a second‐order accurate, Godunov‐type numerical model of the two‐dimensional shallow water equations (SWEs) written in matrix form and discretized using finite volumes. Roe's flux function is used for the convection terms and a non‐linear limiter is applied to prevent unwanted spurious oscillations. A new mathematical formulation is presented, which inherently balances flux gradient and source terms. It is, therefore, suitable for cases where the bathymetry is non‐uniform, unlike other formulations given in the literature based on Roe's approximate Riemann solver. The model is based on hierarchical quadtree (Q‐tree) grids, which adapt to inherent flow parameters, such as magnitude of the free surface gradient and depth‐averaged vorticity. Validation tests include wind‐induced circulation in a dish‐shaped basin, two‐dimensional frictionless rectangular and circular dam‐breaks, an oblique hydraulic jump, and jet‐forced flow in a circular reservoir. Copyright © 2001 John Wiley & Sons, Ltd.     

A mesh patching method for finite volume modelling of shallow water flow

A new mesh‐patching model is presented for shallow water flow described by the 2D non‐linear shallow water (NLSW) equations. The mesh‐patching model is based on AMAZON, a high‐resolution NLSW engine with an improved HLLC approximate Riemann solver. A new patching algorithm has been developed, which not only provides improved spatial resolution of flow features in particular parts of the mesh, but also simplifies and speeds up the (structured) grid generation process for an area with complicated geometry. The new patching technique is also compatible with increasingly popular parallel computing and adaptive grid techniques. The patching algorithm has been tested with moving bores, and results of test problems are presented and compared to previous work. Copyright © 2005 John Wiley & Sons, Ltd.     

Simulation of interface and free surface flows in a viscous fluid using adapting quadtree grids

A new adaptive quadtree method for simulating laminar viscous fluid problems with free surfaces and interfaces is presented in this paper. The Navier–Stokes equations are solved with a SIMPLE‐type scheme coupled with the Compressive Interface Capturing Scheme for Arbitrary Meshes (CICSAM) (Numerical prediction of two fluid systems with sharp interfaces, Ph.D. Thesis, Imperial College of Science, Technology and Medicine, London, 1997) volume of fluid (VoF) method and PLIC reconstruction of the volume fraction field during refinement and derefinement processes. The method is demonstrated for interface advection cases in translating and shearing flow fields and found to provide high interface resolution at low computational cost. The new method is also applied to simulation of the collapse of a water column and the results are in excellent agreement with other published data. The quadtree grids adapt to follow the movement of the free surface, whilst maintaining a band of the smallest cells surrounding the surface. The calculation is made on uniform and adapting quadtree grids and the accuracy of the quadtree calculation is shown to be the same as that made on the equivalent uniform grid. Copyright © 2004 John Wiley & Sons, Ltd.     

Three‐dimensional modeling of non‐hydrostatic free‐surface flows on unstructured grids

A three‐dimensional numerical model is presented for the simulation of unsteady non‐hydrostatic shallow water flows on unstructured grids using the finite volume method. The free surface variations are modeled by a characteristics‐based scheme, which simulates sub‐critical and super‐critical flows. Three‐dimensional velocity components are considered in a collocated arrangement with a σ‐coordinate system. A special treatment of the pressure term is developed to avoid the water surface oscillations. Convective and diffusive terms are approximated explicitly, and an implicit discretization is used for the pressure term to ensure exact mass conservation. The unstructured grid in the horizontal direction and the σ coordinate in the vertical direction facilitate the use of the model in complicated geometries. Solution of the non‐hydrostatic equations enables the model to simulate short‐period waves and vertically circulating flows. Copyright © 2015 John Wiley & Sons, Ltd.     

A robust algorithm for computing fluid flows on highly non‐smooth staggered grids

A new method for computing the fluid flow in complex geometries using highly non‐smooth and non‐orthogonal staggered grid is presented. In a context of the SIMPLE algorithm, pressure and physical tangential velocity components are used as dependent variables in momentum equations. To reduce the sensitivity of the curvature terms in response to coordinate line orientation change, these terms are exclusively computed using Cartesian velocity components in momentum equations. The method is then used to solve some fairly complicated 2‐D and 3‐D flow field using highly non‐smooth grids. The accuracy of results on rough grids (with sharp grid line orientation change and non‐uniformity) was found to be high and the agreement with previous experimental and numerical results was quite good. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

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.     

Large eddy simulation of free surface shallow‐water flow

Shallow‐water flow with free surface frequently occurs in ambient water bodies, in which the horizontal scale of motion is generally two orders of magnitude greater than the water depth. To accurately predict this flow phenomenon in more detail, a three‐dimensional numerical model incorporating the method of large eddy simulation (LES) has been developed and assessed. The governing equations are split into three parts in the finite difference solution: advection, dispersion and propagation. The advection part is solved by the QUICKEST scheme. The dispersion part is solved by the central difference method and the propagation part is solved implicitly using the Gauss–Seidel iteration method. The model has been applied to free surface channel flow for which ample experimental data are available for verification. The inflow boundary condition for turbulence is generated by a spectral line processor. The computed results compare favourably with the experimental data and those results obtained by using a periodic boundary condition. The performance of the model is also assessed for the case in which anisotropic grids and filters with horizontal grid size of the order of the water depth are used for computational efficiency. The coarse horizontal grid was found to cause a significant reduction in the large‐scale turbulent motion generated by the bottom turbulence, and the turbulent motion is predominately described by the sub‐grid scale (SGS) terms. The use of the Smagorinsky model for SGS turbulence in this situation is found inappropriate. A parabolic mixing length model, which accounts for the filtered turbulence, is then proposed. The new model can reproduce more accurately the flow quantities. Copyright © 2000 John Wiley & Sons, Ltd.     

A two‐grid fictitious domain method for direct simulation of flows involving non‐interacting particles of a very small size

The full resolution of flows involving particles whose scale is hundreds or thousands of times smaller than the size of the flow domain is a challenging problem. A naive approach would require a tremendous number of degrees of freedom in order to bridge the gap between the two spatial scales involved. The approach used in the present study employs two grids whose grid size fits the two different scales involved, one of them (the micro‐scale grid) being embedded into the other (the macro‐scale grid). Then resolving first the larger scale on the macro‐scale grid, we transfer the so obtained data to the boundary of the micro‐scale grid and solve the smaller size problem. Since the particle is moving throughout the macro‐scale domain, the micro‐scale grid is fixed at the centroid of the moving particle and therefore moves with it. In this study we combine such an approach with a fictitious domain formulation of the problem resulting in a very efficient algorithm that is also easy to implement in an existing CFD code. We validate the method against existing experimental data for a sedimenting sphere, as well as analytical results for motion of an inertia‐less ellipsoid in a shear flow. Finally, we apply the method to the flow of a high aspect ratio ellipsoid in a model of a human lung airway bifurcation. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical analysis of turbulent structure in compound meandering open channel by algebraic Reynolds stress model

The turbulent flow in a compound meandering channel with a rectangular cross section is one of the most complicated turbulent flows, because the flow behaviour is influenced by several kinds of forces, including centrifugal forces, pressure‐driven forces and shear stresses generated by momentum transfer between the main channel and the flood plain. Numerical analysis has been performed for the fully developed turbulent flow in a compound meandering open‐channel flow using an algebraic Reynolds stress model. The boundary‐fitted coordinate system is introduced as a method for coordinate transformation in order to set the boundary conditions along the complicated shape of the meandering open channel. The turbulence model consists of transport equations for turbulent energy and dissipation, in conjunction with an algebraic stress model based on the Reynolds stress transport equations. With reference to the pressure–strain term, we have made use of a modified pressure–strain term. The boundary condition of the fluctuating vertical velocity is set to zero not only for the free surface, but also for computational grid points next to the free surface, because experimental results have shown that the fluctuating vertical velocity approaches zero near the free surface. In order to examine the validity of the present numerical method and the turbulent model, the calculated results are compared with experimental data measured by laser Doppler anemometer. In addition, the compound meandering open channel is clarified somewhat based on the calculated results. As a result of the analysis, the present algebraic Reynolds stress model is shown to be able to reasonably predict the turbulent flow in a compound meandering open channel. Copyright © 2005 John Wiley & Sons, Ltd.     

Flux and source term discretization in two‐dimensional shallow water models with porosity on unstructured grids

Two‐dimensional shallow water models with porosity appear as an interesting path for the large‐scale modelling of floodplains with urbanized areas. The porosity accounts for the reduction in storage and in the exchange sections due to the presence of buildings and other structures in the floodplain. The introduction of a porosity into the two‐dimensional shallow water equations leads to modified expressions for the fluxes and source terms. An extra source term appears in the momentum equation. This paper presents a discretization of the modified fluxes using a modified HLL Riemann solver on unstructured grids. The source term arising from the gradients in the topography and in the porosity is treated in an upwind fashion so as to enhance the stability of the solution. The Riemann solver is tested against new analytical solutions with variable porosity. A new formulation is proposed for the macroscopic head loss in urban areas. An application example is presented, where the large scale model with porosity is compared to a refined flow model containing obstacles that represent a schematic urban area. The quality of the results illustrates the potential usefulness of porosity‐based shallow water models for large scale floodplain simulations. Copyright © 2005 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.     

Dispersion and stability analyses of the linearized two‐dimensional shallow water equations in boundary‐fitted co‐ordinates

In the present investigation, a Fourier analysis is used to study the phase and group speeds of a linearized, two‐dimensional shallow water equations, in a non‐orthogonal boundary‐fitted co‐ordinate system. The phase and group speeds for the spatially discretized equations, using the second‐order scheme in an Arakawa C grid, are calculated for grids with varying degrees of non‐orthogonality and compared with those obtained from the continuous case. The spatially discrete system is seen to be slightly dispersive, with the degree of dispersivity increasing with an decrease in the grid non‐orthogonality angle or decrease in grid resolution and this is in agreement with the conclusions reached by Sankaranarayanan and Spaulding (J. Comput. Phys., 2003; 184 : 299–320). The stability condition for the non‐orthogonal case is satisfied even when the grid non‐orthogonality angle, is as low as 30° for the Crank Nicolson and three‐time level schemes. A two‐dimensional wave deformation analysis, based on complex propagation factor developed by Leendertse (Report RM‐5294‐PR, The Rand Corp., Santa Monica, CA, 1967), is used to estimate the amplitude and phase errors of the two‐time level Crank–Nicolson scheme. There is no dissipation in the amplitude of the solution. However, the phase error is found to increase, as the grid angle decreases for a constant Courant number, and increases as Courant number increases. Copyright © 2003 John Wiley & Sons, Ltd.