A computationally efficient 3D finite‐volume scheme for violent liquid–gas sloshing

We describe a semi‐implicit volume‐of‐fluid free‐surface‐modelling methodology for flow problems involving violent free‐surface motion. For efficient computation, a hybrid‐unstructured edge‐based vertex‐centred finite volume discretisation is employed, while the solution methodology is entirely matrix free. Pressures are solved using a matrix‐free preconditioned generalised minimum residual algorithm and explicit time‐stepping is employed for the momentum and interface‐tracking equations. The high resolution artificial compressive (HiRAC) volume‐of‐fluid method is used for accurate capturing of the free surface in violent flow regimes while allowing natural applicability to hybrid‐unstructured meshes. The code is parallelised for solution on distributed‐memory architectures and evaluated against 2D and 3D benchmark problems. Good parallel scaling is demonstrated, with almost linear speed‐up down to 6000 cells per core. Finally, the code is applied to an industrial‐type problem involving resonant excitation of a fuel tank, and a comparison with experimental results is made in this violent sloshing regime. Copyright © 2015 John Wiley & Sons, Ltd.     

A weakly compressible MPS method for modeling of open‐boundary free‐surface flow

A mesh‐free particle method, based on the moving particle semi‐implicit (MPS) interaction model, has been developed for the simulation of two‐dimensional open‐boundary free‐surface flows. The incompressibility model in the original MPS has been replaced with a weakly incompressible model. The effect of this replacement on the efficiency and accuracy of the model has been investigated. The new inflow–outflow boundary conditions along with the particle recycling strategy proposed in this study extend the application of the model to open‐boundary problems. The final model is able to simulate open‐boundary free surface flow in cases of large deformation and fragmentation of free surface. The models and proposed algorithms have been validated and applied to sample problems. The results confirm the model's efficiency and accuracy. Copyright © 2009 John Wiley & Sons, Ltd.     

CFD simulation of two‐ and three‐dimensional free‐surface flow

The paper describes the implementation of moving‐mesh and free‐surface capabilities within a 3‐d finite‐volume Reynolds‐averaged‐Navier–Stokes solver, using surface‐conforming multi‐block structured meshes. The free‐surface kinematic condition can be applied in two ways: enforcing zero net mass flux or solving the kinematic equation by a finite‐difference method. The free surface is best defined by intermediate control points rather than the mesh vertices. Application of the dynamic boundary condition to the piezometric pressure at these points provides a hydrostatic restoring force which helps to eliminate any unnatural free‐surface undulations. The implementation of time‐marching methods on moving grids are described in some detail and it is shown that a second‐order scheme must be applied in both scalar‐transport and free‐surface equations if flows driven by free‐surface height variations are to be computed without significant wave attenuation using a modest number of time steps. Computations of five flows of theoretical and practical interest—forced motion in a pump, linear waves in a tank, quasi‐1d flow over a ramp, solitary wave interaction with a submerged obstacle and 3‐d flow about a surface‐penetrating cylinder—are described to illustrate the capabilities of our code and methods. Copyright © 2003 John Wiley & Sons, Ltd.     

A massively parallel GPU‐accelerated model for analysis of fully nonlinear free surface waves

We implement and evaluate a massively parallel and scalable algorithm based on a multigrid preconditioned Defect Correction method for the simulation of fully nonlinear free surface flows. The simulations are based on a potential model that describes wave propagation over uneven bottoms in three space dimensions and is useful for fast analysis and prediction purposes in coastal and offshore engineering. A dedicated numerical model based on the proposed algorithm is executed in parallel by utilizing affordable modern special purpose graphics processing unit (GPU). The model is based on a low‐storage flexible‐order accurate finite difference method that is known to be efficient and scalable on a CPU core (single thread). To achieve parallel performance of the relatively complex numerical model, we investigate a new trend in high‐performance computing where many‐core GPUs are utilized as high‐throughput co‐processors to the CPU. We describe and demonstrate how this approach makes it possible to do fast desktop computations for large nonlinear wave problems in numerical wave tanks (NWTs) with close to 50/100 million total grid points in double/single precision with 4 GB global device memory available. A new code base has been developed in C++ and compute unified device architecture C and is found to improve the runtime more than an order in magnitude in double precision arithmetic for the same accuracy over an existing CPU (single thread) Fortran 90 code when executed on a single modern GPU. These significant improvements are achieved by carefully implementing the algorithm to minimize data‐transfer and take advantage of the massive multi‐threading capability of the GPU device. Copyright © 2011 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 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.     

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.     

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.     

Depth‐averaged two‐dimensional curvilinear explicit finite analytic model for open‐channel flows

A depth‐averaged two‐dimensional model has been developed in the curvilinear co‐ordinate system for free‐surface flow problems. The non‐linear convective terms of the momentum equations are discretized based on the explicit–finite–analytic method with second‐order accuracy in space and first‐order accuracy in time. The other terms of the momentum equations, as well as the mass conservation equation, are discretized by the finite difference method. The discretized governing equations are solved in turn, and iteration in each time step is adopted to guarantee the numerical convergence. The new model has been applied to various flow situations, even for the cases with the presence of sub‐critical and supercritical flows simultaneously or sequentially. Comparisons between the numerical results and the experimental data show that the proposed model is robust with satisfactory accuracy. Copyright © 2000 John Wiley & Sons, Ltd.     

An efficient edge‐based level set finite element method for free surface flow problems

We present an efficient technique for the solution of free surface flow problems using level set and a parallel edge‐based finite element method. An unstructured semi‐explicit solution scheme is proposed. A custom data structure, obtained by blending node‐based and edge‐based approaches is presented so to allow a good parallel performance. In addition to standard velocity extrapolation (for the convection of the level set function), an explicit extrapolation of the pressure field is performed in order to impose both the pressure boundary condition and the volume conservation. The latter is also improved with a modification of the divergence free constrain. The method is shown to allow an efficient solution of both simple benchmark cases and complex industrial examples. Copyright © 2012 John Wiley & Sons, Ltd.     

Semi‐implicit subgrid modelling of three‐dimensional free‐surface flows

In this paper, a semi‐implicit numerical model for two‐ and three‐dimensional free‐surface flows will be formulated in such a fashion as to intrinsically account for subgrid bathymetric details. It will be shown that with the proposed subgrid approach the model accuracy can be substantially improved without increasing the corresponding computational effort. Copyright © 2010 John Wiley & Sons, Ltd.     

Portable data‐parallel surface reconstruction on a uniform rectilinear grid

With the increasing heterogeneity and on‐node parallelism of high‐performance computing hardware, a major challenge is to develop portable and efficient algorithms and software. In this work, we present our implementation of a portable code to perform surface reconstruction using NVIDIA's Thrust library. Surface reconstruction is a technique commonly used in volume tracking methods for simulations of multimaterial flow with interfaces. We have designed a 3D mesh data structure that is easily mapped to the 1D vectors used by Thrust and at the same time is simple to use and uses familiar data structure terminology (such as cells, faces, vertices, and edges). With this new data structure in place, we have implemented a piecewise linear interface reconstruction algorithm in 3 dimensions that effectively exploits the symmetry present in a uniform rectilinear computational cell. Finally, we report performance results, which show that a single implementation of these algorithms can be compiled to multiple backends (specifically, multi‐core CPUs, NVIDIA GPUs, and Intel Xeon Phi processors), making efficient use of the available parallelism on each. We also compare performance of our implementation to a legacy FORTRAN implementation in Message Passing Interface (MPI) and show performance parity on single and multi‐core CPU and achieved good parallel speed‐ups on GPU. Our research demonstrates the advantage of performance portability of the underlying data‐parallel programming model.     

GENSMAC3D: a numerical method for solving unsteady three‐dimensional free surface flows

A numerical method for solving three‐dimensional free surface flows is presented. The technique is an extension of the GENSMAC code for calculating free surface flows in two dimensions. As in GENSMAC, the full Navier–Stokes equations are solved by a finite difference method; the fluid surface is represented by a piecewise linear surface composed of quadrilaterals and triangles containing marker particles on their vertices; the stress conditions on the free surface are accurately imposed; the conjugate gradient method is employed for solving the discrete Poisson equation arising from a velocity update; and an automatic time step routine is used for calculating the time step at every cycle. A program implementing these features has been interfaced with a solid modelling routine defining the flow domain. A user‐friendly input data file is employed to allow almost any arbitrary three‐dimensional shape to be described. The visualization of the results is performed using computer graphic structures such as phong shade, flat and parallel surfaces. Results demonstrating the applicability of this new technique for solving complex free surface flows, such as cavity filling and jet buckling, are presented. Copyright © 2001 John Wiley & Sons, Ltd.     

An unconditionally stable semi‐Lagrangian method for the spherical atmospherical shallow water equations

A semi‐implicit, semi‐Lagrangian, mixed finite difference–finite volume model for the shallow water equations on a rotating sphere is introduced and discussed. Its main features are the vectorial treatment of the momentum equation and the finite volume approach for the continuity equation. Pressure and Coriolis terms in the momentum equation and velocity in the continuity equation are treated semi‐implicitly. Moreover, a splitting technique is introduced to preserve symmetry of the numerical scheme. An alternative asymmetric scheme (without splitting) is also introduced and the efficiency of both is discussed. The model is shown to be conservative in geopotential height and unconditionally stable for 0.5≤θ≤1. Numerical experiments on two standard test problems confirm the performance of the model. Copyright © 2000 John Wiley & Sons, Ltd.     

Some considerations for high‐order ‘incremental remap’‐based transport schemes: edges,reconstructions, and area integration

The problem of two‐dimensional tracer advection on the sphere is extremely important in modeling of geophysical fluids and has been tackled using a variety of approaches. A class of popular approaches for tracer advection include ‘incremental remap’ or cell‐integrated semi‐Lagrangian‐type schemes. These schemes achieve high‐order accuracy without the need for multistage integration in time, are capable of large time steps, and tend to be more efficient than other high‐order transport schemes when applied to a large number of tracers over a single velocity field. In this paper, the simplified flux‐form implementation of the Conservative Semi‐LAgrangian Multi‐tracer scheme (CSLAM) is reformulated using quadratic curves to approximate the upstream flux volumes and Gaussian quadrature for integrating the edge flux. The high‐order treatment of edge fluxes is motivated because of poor accuracy of the CSLAM scheme in the presence of strong nonlinear shear, such as one might observe in the midlatitudes near an atmospheric jet. Without the quadratic treatment of upstream edges, we observe at most second‐order accuracy under convergence of grid resolution, which is returned to third‐order accuracy under the improved treatment. A shallow‐water barotropic instability also reveals clear evidence of grid imprinting without the quadratic correction. Consequently, these tests reveal a problem that might arise in tracer transport near nonlinearly sheared regions of the real atmosphere, particularly near cubed‐sphere panel edges. Although CSLAM is used as the foundation for this analysis, the conclusions of this paper are applicable to the general class of incremental remap schemes. Copyright © 2012 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.     

A two‐dimensional vertical non‐hydrostatic σ model with an implicit method for free‐surface flows

An implicit finite difference model in the σ co‐ordinate system is developed for non‐hydrostatic, two‐dimensional vertical plane free‐surface flows. To accurately simulate interaction of free‐surface flows with uneven bottoms, the unsteady Navier–Stokes equations and the free‐surface boundary condition are solved simultaneously in a regular transformed σ domain using a fully implicit method in two steps. First, the vertical velocity and pressure are expressed as functions of horizontal velocity. Second, substituting these relationship into the horizontal momentum equation provides a block tri‐diagonal matrix system with the unknown of horizontal velocity, which can be solved by a direct matrix solver without iteration. A new treatment of non‐hydrostatic pressure condition at the top‐layer cell is developed and found to be important for resolving the phase of wave propagation. Additional terms introduced by the σ co‐ordinate transformation are discretized appropriately in order to obtain accurate and stable numerical results. The developed model has been validated by several tests involving free‐surface flows with strong vertical accelerations and non‐linear waves interacting with uneven bottoms. Comparisons among numerical results, analytical solutions and experimental data show the capability of the model to simulate free‐surface flow problems. Copyright © 2004 John Wiley & Sons, Ltd.     

On the use of characteristic‐based split meshfree method for solving flow problems

This study presents characteristic‐based split (CBS) algorithm in the meshfree context. This algorithm is the extension of general CBS method which was initially introduced in finite element framework. In this work, the general equations of flow have been represented in the meshfree context. A new finite element and MFree code is developed for solving flow problems. This computational code is capable of solving both time‐dependent and steady‐state flow problems. Numerical simulation of some known benchmark flow problems has been studied. Computational results of MFree method have been compared to those of finite element method. The results obtained have been verified by known numerical, analytical and experimental data in the literature. A number of shape functions are used for field variable interpolation. The performance of each interpolation method is discussed. It is concluded that the MFree method is more accurate than FEM if the same numbers of nodes are used for each solver. Meshfree CBS algorithm is completely stable even at high Reynolds numbers. Copyright © 2007 John Wiley & Sons, Ltd.     

Solving viscoelastic free surface flows of a second‐order fluid using a marker‐and‐cell approach

This work is concerned with the numerical simulation of two‐dimensional viscoelastic free surface flows of a second‐order fluid. The governing equations are solved by a finite difference technique based on the marker‐and‐cell philosophy. A staggered grid is employed and marker particles are used to represent the fluid free surface. Full details for the approximation of the free surface stress conditions are given. The resultant code is validated and convergence is demonstrated. Numerical simulations of the extrudate swell and flow through a planar 4:1 contraction for various values of the Deborah number are presented. Copyright © 2006 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.