The effect of using non-orthogonal boundary-fitted grids for solving the shallow water equations

The purpose of this paper is to examine the effects of using non-orthogonal boundary-fitted grids for the numerical solution of the shallow water equations. Two geometries with well known analytical solutions are introduced in order to investigate the accuracy of the numerical solutions. The results verify that a reasonable departure from orthogonality can be allowed when the rate of change of cell areas is kept sufficiently small (i.e. it is not necessary to create a strictly orthogonal grid when the grid is sufficiently smooth).     

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.     

Implementation of some higher-order convection schemes on non-uniform grids

A generalized formulation is applied to implement the quadratic upstream interpolation (QUICK) scheme, the second-order upwind (SOU) scheme and the second-order hybrid scheme (SHYBRID) on non-uniform grids. The implementation method is simple. The accuracy and efficiency of these higher-order schemes on non-uniform grids are assessed. Three well-known bench mark convection-diffusion problems and a fluid flow problem are revisited using non-uniform grids. These are: (1) transport of a scalar tracer by a uniform velocity field; (2) heat transport in a recirculating flow; (3) two-dimensional non-linear Burgers equations; and (4) a two-dimensional incompressible Navier-Stokes flow which is similar to the classical lid-driven cavity flow. The known exact solutions of the last three problems make it possible to thoroughly evaluate accuracies of various uniform and non-uniform grids. Higher accuracy is obtained for fewer grid points on non-uniform grids. The order of accuracy of the examined schemes is maintained for some tested problems if the distribution of non-uniform grid points is properly chosen.     

A comparison of spatial grids for numerical modelling of flows in near-coastal seas

Three-dimensional algorithms for the numerical computation of flows caused by tides or meteorological forcing are developed for four of Arakawa's spatial grid types using a spectral method in the vertical dimension. Three of the grids, in which the velocity components are computed at the same grid points, offer potential advantages over the commonly used C-grid. The computed results from the four grids are compared for three test problems based on the linearized hydrodynamical equations. It is concluded that the B-grid provides a viable alternative to the C-grid, with significant advantages when a spectral method is used.     

Computational properties of vertical grids for a nonhydrostatic anelastic model

Computational dispersion properties of all vertically staggered grids, which are presently available, are analysed in terms of frequency and group velocity components using the second-order centre difference scheme for a nonhydrostatic anelastic approximation system with a general method. The inertial-gravitational waves with a horizontal scale of a hundred-, ten- and one-kilometres are considered. The comparison analysis shows that the Charny-Phillips (CP) and Lorenz grids are suitable for waves at all abovementioned horizontal scales, while the Lorenz time staggered and Charny-Phillips time staggered grids are applicable only to waves with a horizontal scale less than 10 km. The unstaggered (N) grid is not suitable for simulating waves at any horizontal scale. In an idealised flow numerical test, the result on the CP grid has much less error than that on the N grid.     

An Adaptive Grid Eulerian Method for the Computation of Free Surface Flows

A numerical scheme for the prediction of free surface flows is presented and investigated. The method is based on an adaptive grid Eulerian finite-volume method, where non-orthogonal boundary-fitted moving grids are employed to follow the free surface. The underlying flow solver consists in a pressure-correction scheme of SIMPLE type with multigrid acceleration, which is iteratively combined with the moving grid technique. Several numerical examples are considered to illustrate the capabilities of the approach.     

动态混合网格生成及隐式非定常计算方法

建立了一种基于动态混合网格的非定常数值计算方法. 混合网格由贴体的四边形网格、外场 的多层次矩形网格和中间的三角形网格构成. 当物体运动时，贴体四边形网格随物体运动而 运动，而外场的矩形网格保持静止，中间的三角形网格随之变形；当物体运动位移较大，导 致三角形网格的质量降低，甚至导致网格相交时，在局部重新生成网格. 新网格上的物理量 由旧网格上的物理量插值而得. 为了提高计算效率，采用了双时间步和子迭代相结合的隐式 有限体积格式计算非定常Navier-Stokes方程. 子迭代采用高效的块LU-SGS方法. 利用该 方法数值模拟了NACA0012振荡翼型的无黏和黏性绕流，得到了与实验和他人计算相当一致 的结果.     

OVERLAPPING GRIDS FOR THE SIMULATION OF PARTICLE INTERACTION AT THE MICRO SCALE

This paper describes the use of overlapping grids for the calculation of flow around single and multipleparticle configurations at the micro scale. The basic equations for calculation are those for conservation of mass and momentum which are solved using a common Finite-Volume formulation. The hydrodynamic particle-particle and particle-wall interaction can be calculated by using an overlapping or Chimera grid scheme. With the grid structuring procedure it is possible to use simple and structured grids around the particles and the overall main grid geometry. The particle grids are lapped over the main grid such that they can move independently after each time step without remeshing the whole geometry. The paper gives results for the validation of the code developed for general test cases, for a rotating ellipsoid in simple shear flow, the flow around particles attached to a wall, the motion of a particle in the vicinity of a wall and some results for the flow through a packed bed configuration.     

Extension of an explicit finite volume method to large time steps (CFL>1): application to shallow water flows

In this work, the explicit first order upwind scheme is presented under a formalism that enables the extension of the methodology to large time steps. The number of cells in the stencil of the numerical scheme is related to the allowable size of the CFL number for numerical stability. It is shown how to increase both at the same time. The basic idea is proposed for a 1D scalar equation and extended to 1D and 2D non‐linear systems with source terms. The importance of the kind of grid used is highlighted and the method is outlined for irregular grids. The good quality of the results is illustrated by means of several examples including shallow water flow test cases. The bed slope source terms are involved in the method through an upwind discretization. Copyright © 2005 John Wiley & Sons, Ltd.     

复杂无粘流场数值模拟的矩形/三角形混合网格技术

建立了一套模拟复杂无粘流场的矩形／三角形混合网格技术，其中三角形仅限于物面附近，发挥非结构网格的几何灵活性，以少量的网格模拟复杂外型；同时在以外的区域采用矩形结构网格，发挥矩形网格计算简单快速的优势，有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足．混合网格系统由修正的四分树方法生成．将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算，消除了全矩形网格模拟曲边界的台阶效应，同时保证了网格间的通量守恒．数值实验表明本方法在模拟复杂无粘流场方面的灵活性和高效性．     

A staggered conservative scheme for every Froude number in rapidly varied shallow water flows

This paper proposes a numerical technique that in essence is based upon the classical staggered grids and implicit numerical integration schemes, but that can be applied to problems that include rapidly varied flows as well. Rapidly varied flows occur, for instance, in hydraulic jumps and bores. Inundation of dry land implies sudden flow transitions due to obstacles such as road banks. Near such transitions the grid resolution is often low compared to the gradients of the bathymetry. In combination with the local invalidity of the hydrostatic pressure assumption, conservation properties become crucial. The scheme described here, combines the efficiency of staggered grids with conservation properties so as to ensure accurate results for rapidly varied flows, as well as in expansions as in contractions. In flow expansions, a numerical approximation is applied that is consistent with the momentum principle. In flow contractions, a numerical approximation is applied that is consistent with the Bernoulli equation. Both approximations are consistent with the shallow water equations, so under sufficiently smooth conditions they converge to the same solution. The resulting method is very efficient for the simulation of large‐scale inundations. Copyright © 2003 John Wiley & Sons, Ltd.     

Improved total variation diminishing schemes for advection simulation on arbitrary grids

The requirements for flux limiter functions preserving total variation diminishing (TVD) are derived based on a 1D nonuniform grid, and a new TVD region is determined to fit arbitrary 1D grids. Some second‐order TVD schemes called improved TVD schemes are developed, such as modified Van Leer scheme, modified Van Albada scheme, and modified SUPERBEE scheme. Then, they are extended to 2D grids. Because the element sizes and face positions are taken into account, good behaviors are observed in the implementations in both 1D and 2D cases for pure advection simulation. That is, good conservation, better monotonicity, and higher accuracy are maintained by the improved TVD schemes compared with the present ones deduced from uniform grids, and they keep superiorities even when implemented on poor grids. Among all the improved TVD schemes, the modified SUPERBEE is only recommended for poor 1D grids, but the modified Van Leer scheme can suit both poor 1D and 2D grids. Copyright © 2011 John Wiley & Sons, Ltd.     

A PARALLEL BLOCK-STRUCTURED MULTIGRID METHOD FOR THE PREDICTION OF INCOMPRESSIBLE FLOWS

In this paper a parallel multigrid finite volume solver for the prediction of steady and unsteady flows in complex geometries is presented. For the handling of the complexity of the geometry and for the parallelization a unified approach connected with the concept of block-structured grids is employed. The parallel implementation is based on grid partitioning with automatic load balancing and follows the message-passing concept, ensuring a high degree of portability. A high numerical efficiency is obtained by a non-linear multigrid method with a pressure correction scheme as smoother. By a number of numerical experiments on various parallel computers the method is investigated with respect to its numerical and parallel efficiency. The results illustrate that the high performance of the underlying sequential multigrid algorithm can largely be retained in the parallel implementation and that the proposed method is well suited for solving complex flow problems on parallel computers with high efficiency.     

A parallel adaptive mesh approach for flowfields with shock waves

A simple and efficient implementation of Adaptive Mesh Refinement (AMR) with distributed-memory approach is presented. Introducing a lookup table including grid connectivity information and simplified algorithms for AMR, the procedures for reconstructing adaptive grids are carried out in parallel, with local data to a large extent. A simple static load-balancing scheme is adopted, and the grids are not repartitioned and no data redistribution is performed. A numerical example on two different parallel computers shows that the proposed implementation of AMR is effective to reduce the computational time for unsteady flows with shock waves. Received 23 October 2000 / Accepted 30 March 2001 Published online 11 June 2002     

THE PROGRESS OF THE OVERLAPPING GRID TECHNIQUES

从网格装配和插值计算两个主要方面对现有的重叠网格方法进行了综述。首先,从挖洞方法和建立嵌入网格关系环节的寻点技术出发归纳和介绍了网格装配方法;其次,介绍了数值迭代过程中的插值计算方法,并特别讨论了插值守恒性以及插值计算精度等问题;另外,对重叠网格方法的并行计算和应用成果也作了介绍;最后,通过总结认为重叠网格方法在改进网格装配方法、改善插值和并行计算效率等方面仍需进一步研究。     

Grid‐independent depth‐averaged simulations with a collocated unstructured finite volume scheme

In this article, the depth‐averaged transport equations are written in a new way so that it is possible to solve the transport equations for very small water depths. Variables are interpolated into the cell face with two different schemes and, the schemes are compared in terms of computational cost and accuracy. The bed source terms are computed using two different assumptions. The effect of these assumptions on numerical simulations is then investigated. Solutions of transport equations on different types of unstructured triangular grids are compared and, an appropriate choice of grid is suggested. Copyright © 2011 John Wiley & Sons, Ltd.     

An immersed boundary method for the incompressible Navier–Stokes equations in complex geometry

A nodally exact convection–diffusion–reaction scheme developed in Cartesian grids is applied to solve the flow equations in irregular domains within the framework of immersed boundary (IB) method. The artificial momentum forcing term applied at certain points in the flow and inside the body of any shape allows the imposition of no‐slip velocity condition to account for the body of complex boundary. Development of an interpolation scheme that can accurately lead to no‐slip velocity condition along the IB is essential since Cartesian grid lines generally do not coincide with the IB. The results simulated from the proposed IB method agree well with other numerical and experimental results for several chosen benchmark problems. The accuracy and fidelity of the IB flow solver to predict flows with irregular IBs are therefore demonstrated. Copyright © 2007 John Wiley & Sons, Ltd.     

Stream function–potential function coordinates for aeroacoustics and unsteady aerodynamics

‘Stream function as a coordinate approach’ (SFC) combined with compact high-order finite difference schemes has been developed and applied to aeroacoustics and unsteady aerodynamics problems. Straightforward implementation of SFC creates coarse grids at the vicinity of stagnation points that smears high-order numerical computations. Grid clustering is employed to resolve coarse grid near stagnations points. The agreement between numerical results and particle image velocimetry (PIV) measurements for flapping airfoil shows the robustness of the current approach for performing high-order computations.     

Cell-vertex Based Parallel and Adaptive Explicit 3D Flow Solution on Unstructured Grids

A parallel adaptive Euler flow solution algorithm is developed for 3D applications on distributed memory computers. Significant contribution of this research is the development and implementation of a parallel grid adaptation scheme together with an explicit cell vertex-based finite volume 3D flow solver on unstructured tetrahedral grids. Parallel adaptation of grids is based on grid-regeneration philosophy by using an existing serial grid generation program. Then, a general partitioner repartitions the grid. An adaptive sensor value, which is a measure to refine or coarsen grids, is calculated considering the pressure gradients in all partitioned blocks of grids. The parallel performance of the present study was tested. Parallel computations were performed on Unix workstations and a Linux cluster using MPI communication library. The present results show that overall adaptation scheme developed in this study is applicable to any pair of a flow solver and grid generator with affordable cost. It is also proved that parallel adaptation is necessary for accurate and efficient flow solutions.     

A finite volume Navier-Stokes algorithm for adaptive grids

A compact, finite volume, time-marching scheme for the two-dimensional Navier-Stokes equations of viscous fluid flow is presented. The scheme is designed for unstructured (locally refined) quadrilateral meshes. An earlier inviscid equation (Euler) scheme is employed for the convective terms and the emphasis is on treatment of the viscous terms. An essential feature of the algorithm is that all necessary operations are restricted to within each cell, which is very important when dealing with unstructured grids. Numerical issues which have to be addressed when developing a Navier-Stokes scheme are investigated. These issues are not limited to the particular Navier-Stokes scheme developed in the present work but are general problems. Specifically, the extent of the numerical molecule, which is related to the compactness of the scheme and to its suitability for unstructured grids, is examined. An approach which considers suppression of odd-even mode decoupling of the solution when designing a scheme is presented. In addition, accuracy issues related to grid stretching as well as boundary layer solution contamination due to artificial dissipation are addressed. Although the above issues are investigated with respect to the specific scheme presented, the conclusions are valid for an entire class of finite volume algorithms. The Navier-Stokes solver is validated through test cases which involve comparisons with analytical, numerical and experimental results. The solver is coupled to an adaptive algorithm for high-Reynolds-number aerofoil flow computations.