An unstructured finite volume technique for the 3D Poisson equation on arbitrary geometry using a σ‐coordinate system

We present a solver for a three‐dimensional Poisson equation issued from the Navier–Stokes equations applied to model rivers, estuaries, and coastal flows. The three‐dimensional physical domain is composed of an arbitrary domain in the horizontal direction and is bounded by an irregular free surface and bottom in the vertical direction. The equations are transformed vertically to the σ‐coordinate system to obtain an accurate representation of top and bottom topographies. The method is based on a second‐order finite volume technique on prisms consisting of triangular grids in the horizontal direction. The algorithm is accompanied by an analysis of different linear system solvers in order to achieve fast solutions. Numerical experiments are conducted to test the numerical accuracy and the computational efficiency of the proposed method. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical modelling of 3D stratified free surface flows: a case study of sediment dumping

A three‐dimensional numerical model has been developed to simulate stratified flows with free surfaces. The model is based on the Reynolds‐averaged Navier–Stokes (RANS) equations with variable fluid density. The equations are solved in a transformed σ‐coordinate system with the use of operator‐splitting method (Int. J. Numer. Meth. Fluids 2002; 38 :1045–1068). The numerical model is validated against the one‐dimensional diffusion problem and the two‐dimensional density‐gradient flow. Excellent agreements are obtained between numerical results and analytical solutions. The model is then used to study transport phenomena of dumped sediments into a water body, which has been modelled as a strongly stratified flow. For the two‐dimensional problem, the numerical results compare well with experimental data in terms of mean particle falling velocity and spreading rate of the sediment cloud for both coarse and medium‐size sediments. The model is also employed to study the dumping of sediments in a three‐dimensional environment with the presence of free surface. It is found that during the descending process an annulus‐like cloud is formed for fine sediments whereas a plate‐like cloud for medium‐size sediments. The model is proven to be a good tool to simulate strongly stratified free surface flows. Copyright © 2005 John Wiley & Sons, Ltd.     

Calculation of three‐dimensional turbulent flow with a finite volume multigrid method

A numerical method for the efficient calculation of three‐dimensional incompressible turbulent flow in curvilinear co‐ordinates is presented. The mathematical model consists of the Reynolds averaged Navier–Stokes equations and the k–ε turbulence model. The numerical method is based on the SIMPLE pressure‐correction algorithm with finite volume discretization in curvilinear co‐ordinates. To accelerate the convergence of the solution method a full approximation scheme‐full multigrid (FAS‐FMG) method is utilized. The solution of the k–ε transport equations is embedded in the multigrid iteration. The improved convergence characteristic of the multigrid method is demonstrated by means of several calculations of three‐dimensional flow cases. Copyright © 1999 John Wiley & Sons, Ltd.     

Dilation‐based shock capturing for high‐order methods

In this paper, we introduce a shock‐capturing artificial viscosity technique for high‐order unstructured mesh methods. This artificial viscosity model is based on a non‐dimensional form of the divergence of the velocity. The technique is an extension and improvement of the dilation‐based artificial viscosity methods introduced in Premasuthan et al., 15 and further extended in Nguyen and Peraire 27 . The approach presented has a number attractive properties including non‐dimensional analytical form, sub‐cell resolution, and robustness for complex shock flows on anisotropic meshes. We present extensive numerical results to demonstrate the performance of the proposed approach. Copyright © 2016 John Wiley & Sons, Ltd.     

Central finite volume methods for ideal and shallow water magnetohydrodynamics

We propose two‐dimensional central finite volume methods based on our multidimensional extensions of Nessyahu and Tadmor's one‐dimensional non‐oscillatory central scheme and a constrained transport‐type method to solve ideal magnetohydrodynamic problems (MHD) and shallow water magnetohydrodynamic problems (SMHD). The main numerical scheme is second‐order accurate both in space and time and uses an original Cartesian grid coupled to a Cartesian‐ or diamond‐staggered dual grid to by‐pass the resolution of the Riemann problems at the cell interfaces. To treat the non‐vanishing magnetic field/flux divergence we have constructed an adaptation of Evans and Hawley's constrained transport method specifically designed for central schemes. Our numerical results show the efficiency and the potential of the scheme. Copyright © 2010 John Wiley & Sons, Ltd.     

A numerical model for run‐up of breaking waves

In the present paper, the numerical method for the three‐dimensional run‐up, given in Johnsgard and Pedersen [‘A numerical model for three‐dimensional run‐up’, Int. J. Numer. Methods Fluids, 24 , 913–931 (1997)], is extended to include wave breaking. In the fundamental problem of run‐up of a uniform bore, the present model is compared with analytical solutions from the literature. The numerical solutions converge, but very slowly. This is not due to the numerical model, but rather to the structure of the solutions themselves. Numerical results for two realistic but simplified tsunami cases are also presented. In the first case, two‐dimensional simulations are performed concerning the run‐up of a tsunami in Portugal, in the second case, the three dimensional wave pattern generated after a slide in Tafjord, Norway in 1931, is studied. A discussion of different aspects of the model is summarized at the end of the paper. Copyright © 1999 John Wiley & Sons, Ltd.     

A finite element method for the two‐dimensional extended Boussinesq equations

A new numerical method for Nwogu's (ASCE Journal of Waterway, Port, Coastal and Ocean Engineering 1993; 119 :618)two‐dimensional extended Boussinesq equations is presented using a linear triangular finite element spatial discretization coupled with a sophisticated adaptive time integration package. The authors have previously presented a finite element method for the one‐dimensional form of these equations (M. Walkley and M. Berzins (International Journal for Numerical Methods in Fluids 1999; 29 (2):143)) and this paper describes the extension of these ideas to the two‐dimensional equations and the application of the method to complex geometries using unstructured triangular grids. Computational results are presented for two standard test problems and a realistic harbour model. Copyright © 2002 John Wiley & Sons, Ltd.     

A non‐iterative implicit algorithm for the solution of advection–diffusion equation on a sphere

A numerical algorithm for the solution of advection–diffusion equation on the surface of a sphere is suggested. The velocity field on a sphere is assumed to be known and non‐divergent. The discretization of advection–diffusion equation in space is carried out with the help of the finite volume method, and the Gauss theorem is applied to each grid cell. For the discretization in time, the symmetrized double‐cycle componentwise splitting method and the Crank–Nicolson scheme are used. The numerical scheme is of second order approximation in space and time, correctly describes the balance of mass of substance in the forced and dissipative discrete system and is unconditionally stable. In the absence of external forcing and dissipation, the total mass and L₂‐norm of solution of discrete system is conserved in time. The one‐dimensional periodic problems arising at splitting in the longitudinal direction are solved with Sherman–Morrison's formula and Thomas's algorithm. The one‐dimensional problems arising at splitting in the latitudinal direction are solved by the bordering method that requires a prior determination of the solution at the poles. The resulting linear systems have tridiagonal matrices and are solved by Thomas's algorithm. The suggested method is direct (without iterations) and rapid in realization. It can also be applied to linear and nonlinear diffusion problems, some elliptic problems and adjoint advection–diffusion problems on a sphere. Copyright © 2015 John Wiley & Sons, Ltd.     

Applications of the artificial compressibility method for turbulent open channel flows

A three‐dimensional (3‐D) numerical method for solving the Navier–Stokes equations with a standard k–ε turbulence model is presented. In order to couple pressure with velocity directly, the pressure is divided into hydrostatic and hydrodynamic parts and the artificial compressibility method (ACM) is employed for the hydrodynamic pressure. By introducing a pseudo‐time derivative of the hydrodynamic pressure into the continuity equation, the incompressible Navier–Stokes equations are changed from elliptic‐parabolic to hyperbolic‐parabolic equations. In this paper, a third‐order monotone upstream‐centred scheme for conservation laws (MUSCL) method is used for the hyperbolic equations. A system of discrete equations is solved implicitly using the lower–upper symmetric Gauss–Seidel (LU‐SGS) method. This newly developed numerical method is validated against experimental data with good agreement. Copyright © 2005 John Wiley & Sons, Ltd.     

A stabilized Nitsche‐type extended embedding mesh approach for 3D low‐ and high‐Reynolds‐number flows

This paper presents a stabilized extended finite element method (XFEM) based fluid formulation to embed arbitrary fluid patches into a fixed background fluid mesh. The new approach is highly beneficial when it comes to computational grid generation for complex domains, as it allows locally increased resolutions independent from size and structure of the background mesh. Motivating applications for such a domain decomposition technique are complex fluid‐structure interaction problems, where an additional boundary layer mesh is used to accurately capture the flow around the structure. The objective of this work is to provide an accurate and robust XFEM‐based coupling for low‐ as well as high‐Reynolds‐number flows. Our formulation is built from the following essential ingredients: Coupling conditions on the embedded interface are imposed weakly using Nitsche's method supported by extra terms to guarantee mass conservation and to control the convective mass transport across the interface for transient viscous‐dominated and convection‐dominated flows. Residual‐based fluid stabilizations in the interior of the fluid subdomains and accompanying face‐oriented fluid and ghost‐penalty stabilizations in the interface zone stabilize the formulation in the entire fluid domain. A detailed numerical study of our stabilized embedded fluid formulation, including an investigation of variants of Nitsche's method for viscous flows, shows optimal error convergence for viscous‐dominated and convection‐dominated flow problems independent of the interface position. Challenging two‐dimensional and three‐dimensional numerical examples highlight the robustness of our approach in all flow regimes: benchmark computations for laminar flow around a cylinder, a turbulent driven cavity flow at Re = 10000 and the flow interacting with a three‐dimensional flexible wall. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of preconditioned iterative solvers for incompressible flow problems

Solving efficiently the incompressible Navier–Stokes equations is a major challenge, especially in the three‐dimensional case. The approach investigated by Elman et al. (Finite Elements and Fast Iterative Solvers. Oxford University Press: Oxford, 2005) consists in applying a preconditioned GMRES method to the linearized problem at each iteration of a nonlinear scheme. The preconditioner is built as an approximation of an ideal block‐preconditioner that guarantees convergence in 2 or 3 iterations. In this paper, we investigate the numerical behavior for the three‐dimensional lid‐driven cavity problem with wedge elements; the ultimate motivation of this analysis is indeed the development of a preconditioned Krylov solver for stratified oceanic flows which can be efficiently tackled using such meshes. Numerical results for steady‐state solutions of both the Stokes and the Navier–Stokes problems are presented. Theoretical bounds on the spectrum and the rate of convergence appear to be in agreement with the numerical experiments. Sensitivity analysis on different aspects of the structure of the preconditioner and the block decomposition strategies are also discussed. Copyright © 2011 John Wiley & Sons, Ltd.     

Resolution to a three‐dimensional physical oceanographic problem using the non‐linear Galerkin method

This paper presents the numerical results concerning the adaptation of the non‐linear Galerkin method to three‐dimensional geophysical fluid equations. This method was developed by Marion and Temam to solve the Navier–Stokes two‐dimensional equations. It allows a substantial decrease in calculation costs due to the application of an appropriate treatment to each mode based on its position in the spectrum. The large scales involved in the study of geophysical flow require that the earth's rotational effects and the existence of a high degree of stratification be taken into account. These phenomena play an important role in the distribution of the energy spectrum. It is shown here that the non‐linear Galerkin method is very well‐suited to the treatment of these phenomena. First, the method for the particular situation of a rigid‐lid with a flat bottom is validated, for which the functional basis used is particularly well‐adapted. Then the more general case of a domain exhibiting variable bathymetry is presented, which necessitates the use of the transformation σ, thus providing a study domain with a cylindrical configuration. Copyright © 1999 John Wiley & Sons, Ltd.     

Model reduction using L1‐norm minimization as an application to nonlinear hyperbolic problems

We are interested in the model reduction techniques for hyperbolic problems, particularly in fluids. This paper, which is a continuation of an earlier paper of Abgrall et al, proposes a dictionary approach coupled with an L¹ minimization approach. We develop the method and analyze it in simplified 1‐dimensional cases. We show in this case that error bounds with the full model can be obtained provided that a suitable minimization approach is chosen. The capability of the algorithm is then shown on nonlinear scalar problems, 1‐dimensional unsteady fluid problems, and 2‐dimensional steady compressible problems. A short discussion on the cost of the method is also included in this paper.     

A discontinuous Galerkin scheme for the numerical solution of flow problems with discontinuities

In this paper, a new discontinuous Galerkin finite element method for the numerical solution of flow problems with discontinuities is presented. The method is based on the limitation in every cell of the difference between the extrema values and the mean value of the numerical solution. The algorithm and technical details for the implementation of the method are presented in one‐and two‐dimensional problems. Numerical experiments for classical test problems are solved on unstructured triangulations to demonstrate the performance of the proposed method. Copyright © 2011 John Wiley & Sons, Ltd.     

A convergence accelerator of a linear system of equations based upon the power method

This paper considers the convergence rate of an iterative numerical scheme as a method for accelerating at the post‐processor stage. The methodology adapted here is: (1) residual eigenmodes included in the origin of the convex hull are eliminated; (2) remaining residual terms are smoothed away by the main convergence algorithm. For this purpose, the polynomial matrix approach is employed for deriving the characteristic equation by two different methods. The first method is based on vector scaling and the second is based on the normal equations approach. The input for both methods is the solution difference between two consecutive iteration/cycle levels obtained from the main program. The singular value decomposition was employed for both methods due to the ill‐conditioned structure of the matrices. The use of the explicit form of the Richardson extrapolation in the present work overrules the need to employ the Richardson iteration with a Leja ordering. The performance of these methods was compared with the GMRES algorithm for three representative problems: two‐dimensional boundary value problem using the Laplace equation, three‐dimensional multi‐grid, potential solution over a sphere and the one‐dimensional steady state Burger equation. In all three examples both methods have the same rate of convergence, or better, as that of the GMRES method in terms of computer operational count. However, in terms of storage requirements, the method based upon vector scaling has a significant advantage over the normal equations approach as well as the GMRES method, in which only one vector of the N grid‐points is required. Copyright © 2001 John Wiley & Sons, Ltd.     

Divergence preserving reconstruction of the nodal components of a vector field from its normal components to edges

We have developed a new divergence preserving method for the reconstruction of the Cartesian components of a vector field from the orthogonal projection of a vector field to the normals to edges in two dimensional. In this method, discrete divergences computed from the nodal components and from the normal ones are exactly the same. Our new method consists of two stages. At the first stage, we use an extended version of the local procedure described in [J. Comput. Phys., 139 :406–409, 1998] to obtain a ‘reference’ nodal vector. This local procedure is exact for linear vector fields; however, the discrete divergence is not preserved. Then, we formulate a constrained optimization problem, in which this reference vector plays the role of a target, and the divergence constraints are enforced by using Lagrange multipliers. It leads to the solution of ‘elliptic’ like discrete equations for the cell‐centered Lagrange multipliers. The new global divergence preserving method is exact for linear vector fields. We describe all details of our new method and present numerical results, which confirm our theory. Copyright © 2016 John Wiley & Sons, Ltd.     

An efficient and accurate algebraic interface capturing method for unstructured grids in 2 and 3 dimensions: The THINC method with quadratic surface representation

We present in this paper an efficient and accurate volume of fluid (VOF) type scheme to compute moving interfaces on unstructured grids with arbitrary quadrilateral mesh elements in 2D and hexahedral elements in 3D. Being an extension of the multi‐dimensional tangent of hyperbola interface capturing (THINC) reconstruction proposed by the authors in Cartesian grid, an algebraic VOF scheme is devised for arbitrary quadrilateral and hexahedral elements. The interface is cell‐wisely approximated by a quadratic surface, which substantially improves the numerical accuracy. The same as the other THINC type schemes, the present method does not require the explicit geometric representation of the interface when computing numerical fluxes and thus is very computationally efficient and straightforward in implementation. The proposed scheme has been verified by benchmark tests, which reveal that this scheme is able to produce high‐quality numerical solutions of moving interfaces in unstructured grids and thus a practical method for interfacial multi‐phase flow simulations. Copyright © 2014 John Wiley & Sons, Ltd.     

A σ‐coordinate three‐dimensional numerical model for surface wave propagation

A three‐dimensional numerical model based on the full Navier–Stokes equations (NSE) in σ‐coordinate is developed in this study. The σ‐coordinate transformation is first introduced to map the irregular physical domain with the wavy free surface and uneven bottom to the regular computational domain with the shape of a rectangular prism. Using the chain rule of partial differentiation, a new set of governing equations is derived in the σ‐coordinate from the original NSE defined in the Cartesian coordinate. The operator splitting method (Li and Yu, Int. J. Num. Meth. Fluids 1996; 23 : 485–501), which splits the solution procedure into the advection, diffusion, and propagation steps, is used to solve the modified NSE. The model is first tested for mass and energy conservation as well as mesh convergence by using an example of water sloshing in a confined tank. Excellent agreements between numerical results and analytical solutions are obtained. The model is then used to simulate two‐ and three‐dimensional solitary waves propagating in constant depth. Very good agreements between numerical results and analytical solutions are obtained for both free surface displacements and velocities. Finally, a more realistic case of periodic wave train passing through a submerged breakwater is simulated. Comparisons between numerical results and experimental data are promising. The model is proven to be an accurate tool for consequent studies of wave‐structure interaction. Copyright © 2002 John Wiley & Sons, Ltd.     

Simulation of non‐linear free surface motions in a cylindrical domain using a Chebyshev–Fourier spectral collocation method

When a liquid is perturbed, its free surface may experience highly non‐linear motions in response. This paper presents a numerical model of the three‐dimensional hydrodynamics of an inviscid liquid with a free surface. The mathematical model is based on potential theory in cylindrical co‐ordinates with a σ‐transformation applied between the bed and free surface in the vertical direction. Chebyshev spectral elements discretize space in the vertical and radial directions; Fourier spectral elements are used in the angular direction. Higher derivatives are approximated using a collocation (or pseudo‐spectral) matrix method. The numerical scheme is validated for non‐linear transient sloshing waves in a cylindrical tank containing a circular surface‐piercing cylinder at its centre. Excellent agreement is obtained with Ma and Wu's [Second order transient waves around a vertical cylinder in a tank. Journal of Hydrodynamics 1995; Ser. B4 : 72–81] second‐order potential theory. Further evidence for the capability of the scheme to predict complicated three‐dimensional, and highly non‐linear, free surface motions is given by the evolution of an impulse wave in a cylindrical tank and in an open domain. Copyright © 2001 John Wiley & Sons, Ltd.     

Incompressible multiphase flow and encapsulation simulations using the moment‐of‐fluid method

A moment‐of‐fluid method is presented for computing solutions to incompressible multiphase flows in which the number of materials can be greater than two. In this work, the multimaterial moment‐of‐fluid interface representation technique is applied to simulating surface tension effects at points where three materials meet. The advection terms are solved using a directionally split cell integrated semi‐Lagrangian algorithm, and the projection method is used to evaluate the pressure gradient force term. The underlying computational grid is a dynamic block‐structured adaptive grid. The new method is applied to multiphase problems illustrating contact‐line dynamics, triple junctions, and encapsulation in order to demonstrate its capabilities. Examples are given in two‐dimensional, three‐dimensional axisymmetric (R–Z), and three‐dimensional (X–Y–Z) coordinate systems. Copyright © 2015 John Wiley & Sons, Ltd.