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.     

Accuracy and efficiency of two numerical methods of solving the potential flow problem for highly nonlinear and dispersive water waves

The accuracy and efficiency of two methods of resolving the exact potential flow problem for nonlinear waves are compared using three different one horizontal dimension (1DH) test cases. The two model approaches use high‐order finite difference schemes in the horizontal dimension and differ in the resolution of the vertical dimension. The first model uses high‐order finite difference schemes also in the vertical, while the second model applies a spectral approach. The convergence, accuracy, and efficiency of the two models are demonstrated as a function of the temporal, horizontal, and vertical resolutions for the following: (1) the propagation of regular nonlinear waves in a periodic domain; (2) the motion of nonlinear standing waves in a domain with fully reflective boundaries; and (3) the propagation and shoaling of a train of waves on a slope. The spectral model approach converges more rapidly as a function of the vertical resolution. In addition, with equivalent vertical resolution, the spectral model approach shows enhanced accuracy and efficiency in the parameter range used for practical model applications. Copyright © 2015 John Wiley & Sons, Ltd.     

A finite volume unstructured multigrid method for efficient computation of unsteady incompressible viscous flows

An unstructured non‐nested multigrid method is presented for efficient simulation of unsteady incompressible Navier–Stokes flows. The Navier–Stokes solver is based on the artificial compressibility approach and a higher‐order characteristics‐based finite‐volume scheme on unstructured grids. Unsteady flow is calculated with an implicit dual time stepping scheme. For efficient computation of unsteady viscous flows over complex geometries, an unstructured multigrid method is developed to speed up the convergence rate of the dual time stepping calculation. The multigrid method is used to simulate the steady and unsteady incompressible viscous flows over a circular cylinder for validation and performance evaluation purposes. It is found that the multigrid method with three levels of grids results in a 75% reduction in CPU time for the steady flow calculation and 55% reduction for the unsteady flow calculation, compared with its single grid counterparts. The results obtained are compared with numerical solutions obtained by other researchers as well as experimental measurements wherever available and good agreements are obtained. Copyright © 2004 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.     

Implicit preconditioned high‐order compact scheme for the simulation of the three‐dimensional incompressible Navier–Stokes equations with pseudo‐compressibility method

This paper combines the pseudo‐compressibility procedure, the preconditioning technique for accelerating the time marching for stiff hyperbolic equations, and high‐order accurate central compact scheme to establish the code for efficiently and accurately solving incompressible flows numerically based on the finite difference discretization. The spatial scheme consists of the sixth‐order compact scheme and 10th‐order numerical filter operator for guaranteeing computational stability. The preconditioned pseudo‐compressible Navier–Stokes equations are marched temporally using the implicit lower–upper symmetric Gauss–Seidel time integration method, and the time accuracy is improved by the dual‐time step method for the unsteady problems. The efficiency and reliability of the present procedure are demonstrated by applications to Taylor decaying vortices phenomena, double periodic shear layer rolling‐up problem, laminar flow over a flat plate, low Reynolds number unsteady flow around a circular cylinder at Re = 200, high Reynolds number turbulence flow past the S809 airfoil, and the three‐dimensional flows through two 90°curved ducts of square and circular cross sections, respectively. It is found that the numerical results of the present algorithm are in good agreement with theoretical solutions or experimental data. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical simulation of free‐surface flow using the level‐set method with global mass correction

A new numerical method that couples the incompressible Navier–Stokes equations with the global mass correction level‐set method for simulating fluid problems with free surfaces and interfaces is presented in this paper. The finite volume method is used to discretize Navier–Stokes equations with the two‐step projection method on a staggered Cartesian grid. The free‐surface flow problem is solved on a fixed grid in which the free surface is captured by the zero level set. Mass conservation is improved significantly by applying a global mass correction scheme, in a novel combination with third‐order essentially non‐oscillatory schemes and a five stage Runge–Kutta method, to accomplish advection and re‐distancing of the level‐set function. The coupled solver is applied to simulate interface change and flow field in four benchmark test cases: (1) shear flow; (2) dam break; (3) travelling and reflection of solitary wave and (4) solitary wave over a submerged object. The computational results are in excellent agreement with theoretical predictions, experimental data and previous numerical simulations using a RANS‐VOF method. The simulations reveal some interesting free‐surface phenomena such as the free‐surface vortices, air entrapment and wave deformation over a submerged object. Copyright © 2009 John Wiley & Sons, Ltd.     

An experimental and numerical investigation of container filling with viscous liquids

This work is concerned with a study of container filling, with particular reference to the food industry. A computer code was developed and an experimental rig was built, the main purpose being to validate the software. The computational fluid dynamic (CFD) code, called GENSMAC, was specifically designed for relatively slow viscous flow and was capable of capturing multiple free surfaces. This paper focuses on the design of the experimental rig and how it functions. The visual output of the code is then compared with high‐speed photographic shots of glucose syrup being jetted into a tub for a selected number of flow regimes. Copyright © 1999 John Wiley & Sons, Ltd.     

A semi‐implicit scheme for 3D free surface flows with high‐order velocity reconstruction on unstructured Voronoi meshes

In this paper, we present a computationally efficient semi‐implicit scheme for the simulation of three‐dimensional hydrostatic free surface flow problems on staggered unstructured Voronoi meshes. For each polygonal control volume, the pressure is defined in the cell center, whereas the discrete velocity field is given by the normal velocity component at the cell faces. A piecewise high‐order polynomial vector velocity field is then reconstructed from the scalar normal velocities at the cell faces by using a new high‐order constrained least‐squares reconstruction operator. The reconstructed high‐order piecewise polynomial velocity field is used for trajectory integration in a semi‐Lagrangian approach to discretize the nonlinear convective terms in the governing PDE. For that purpose, a high‐order Taylor method is used as ODE integrator. The resulting semi‐implicit algorithm is extensively validated on a large set of different academic test problems with exact analytical solution and is finally applied to a real‐world engineering problem consisting of a curved channel upstream of two micro‐turbines of a hydroelectric power plant. For this realistic case, some experimental reference data are available from field measurements. Copyright © 2012 John Wiley & Sons, Ltd.     

Analysis of hybrid algorithms for the Navier–Stokes equations with respect to hydrodynamic stability theory

Hybrid three‐dimensional algorithms for the numerical integration of the incompressible Navier–Stokes equations are analyzed with respect to hydrodynamic stability in both linear and nonlinear fields. The computational schemes are mixed—spectral and finite differences—and are applied to the case of the channel flow driven by constant pressure gradient; time marching is handled with the fractional step method. Different formulations—fully explicit convective term, partially and fully implicit viscous term combined with uniform, stretched, staggered and non‐staggered meshes, x‐velocity splitted and non‐splitted in average and perturbation component – are analyzed by monitoring the evolution in time of both small and finite amplitude perturbations of the mean flow. The results in the linear field are compared with correspondent solutions of the Orr–Sommerfeld equation; in the nonlinear field, the comparison is made with results obtained by other authors. Copyright © 2002 John Wiley & Sons, Ltd.     

A higher resolution edge‐based finite volume method for the simulation of the oil–water displacement in heterogeneous and anisotropic porous media using a modified IMPES method

In this article, we present a higher‐order finite volume method with a ‘Modified Implicit Pressure Explicit Saturation’ (MIMPES) formulation to model the 2D incompressible and immiscible two‐phase flow of oil and water in heterogeneous and anisotropic porous media. We used a median‐dual vertex‐centered finite volume method with an edge‐based data structure to discretize both, the elliptic pressure and the hyperbolic saturation equations. In the classical IMPES approach, first, the pressure equation is solved implicitly from an initial saturation distribution; then, the velocity field is computed explicitly from the pressure field, and finally, the saturation equation is solved explicitly. This saturation field is then used to re‐compute the pressure field, and the process follows until the end of the simulation is reached. Because of the explicit solution of the saturation equation, severe time restrictions are imposed on the simulation. In order to circumvent this problem, an edge‐based implementation of the MIMPES method of Hurtado and co‐workers was developed. In the MIMPES approach, the pressure equation is solved, and the velocity field is computed less frequently than the saturation field, using the fact that, usually, the velocity field varies slowly throughout the simulation. The solution of the pressure equation is performed using a modification of Crumpton's two‐step approach, which was designed to handle material discontinuity properly. The saturation equation is solved explicitly using an edge‐based implementation of a modified second‐order monotonic upstream scheme for conservation laws type method. Some examples are presented in order to validate the proposed formulation. Our results match quite well with others found in literature. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical simulation of standing wave with 3D predictor-corrector finite difference method for potential flow equations

A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.