Coupled ghost fluid/two‐phase level set method for curvilinear body‐fitted grids

A coupled ghost fluid/two‐phase level set method to simulate air/water turbulent flow for complex geometries using curvilinear body‐fitted grids is presented. The proposed method is intended to treat ship hydrodynamics problems. The original level set method for moving interface flows was based on Heaviside functions to smooth all fluid properties across the interface. We call this the Heaviside function method (HFM). The HFM requires fine grids across the interface. The ghost fluid method (GFM) has been designed to explicitly enforce the interfacial jump conditions, but the implementation of the jump conditions in curvilinear grids is intricate. To overcome these difficulties a coupled GFM/HFM method was developed in which approximate jump conditions are derived for piezometric pressure and velocity and pressure gradients based on exact continuous velocity and stress and jump in momentum conditions with the jump in density maintained but continuity of the molecular and turbulent viscosities imposed. The implementation of the ghost points is such that no duplication of memory storage is necessary. The level set method is adopted to locate the air/water interface, and a fast marching method was implemented in curvilinear grids to reinitialize the level set function. Validations are performed for three tests: super‐ and sub‐critical flow without wave breaking and an impulsive plunging wave breaking over 2D submerged bumps, and the flow around surface combatant model DTMB 5512. Comparisons are made against experimental data, HFM and single‐phase level set computations. The proposed method performed very well and shows great potential to treat complicated turbulent flows related to ship flows. Copyright © 2007 John Wiley & Sons, Ltd.     

Semi‐coupled air/water immersed boundary approach for curvilinear dynamic overset grids with application to ship hydrodynamics

For many problems in ship hydrodynamics, the effects of air flow on the water flow are negligible (the frequently called free surface conditions), but the air flow around the ship is still of interest. A method is presented where the water flow is decoupled from the air solution, but the air flow uses the unsteady water flow as a boundary condition. The authors call this a semi‐coupled air/water flow approach. The method can be divided into two steps. At each time step the free surface water flow is computed first with a single‐phase method assuming constant pressure and zero stress on the interface. The second step is to compute the air flow assuming the free surface as a moving immersed boundary (IB). The IB method developed for Cartesian grids (Annu. Rev. Fluid Mech. 2005; 37 :239–261) is extended to curvilinear grids, where no‐slip and continuity conditions are used to enforce velocity and pressure boundary conditions for the air flow. The forcing points close to the IB can be computed and corrected under a sharp interface condition, which makes the computation very stable. The overset implementation is similar to that of the single‐phase solver (Comput. Fluids 2007; 36 :1415–1433), with the difference that points in water are set as IB points even if they are fringe points. Pressure–velocity coupling through pressure implicit with splitting of operators or projection methods is used for water computations, and a projection method is used for the air. The method on each fluid is a single‐phase method, thus avoiding ill‐conditioned numerical systems caused by large differences of fluid properties between air and water. The computation is only slightly slower than the single‐phase version, with complete absence of spurious velocity oscillations near the free surface, frequently present in fully coupled approaches. Validations are performed for laminar Couette flow over a wavy boundary by comparing with the analytical solution, and for the surface combatant model David Taylor Model Basin (DTMB) 5512 by comparing with Experimental Fluid Dynamics (EFD) and the results of two‐phase level set computations. Complex flow computations are demonstrated for the ONR Tumblehome DTMB 5613 with superstructure subject to waves and wind, including 6DOF motions and broaching in SS7 irregular waves and wind. Copyright © 2008 John Wiley & Sons, Ltd.     

Finite‐element/level‐set/operator‐splitting (FELSOS) approach for computing two‐fluid unsteady flows with free moving interfaces

The present work is devoted to the study on unsteady flows of two immiscible viscous fluids separated by free moving interface. Our goal is to elaborate a unified strategy for numerical modelling of two‐fluid interfacial flows, having in mind possible interface topology changes (like merger or break‐up) and realistically wide ranges for physical parameters of the problem. The proposed computational approach essentially relies on three basic components: the finite element method for spatial approximation, the operator‐splitting for temporal discretization and the level‐set method for interface representation. We show that the finite element implementation of the level‐set approach brings some additional benefits as compared to the standard, finite difference level‐set realizations. In particular, the use of finite elements permits to localize the interface precisely, without introducing any artificial parameters like the interface thickness; it also allows to maintain the second‐order accuracy of the interface normal, curvature and mass conservation. The operator‐splitting makes it possible to separate all major difficulties of the problem and enables us to implement the equal‐order interpolation for the velocity and pressure. Diverse numerical examples including simulations of bubble dynamics, bifurcating jet flow and Rayleigh–Taylor instability are presented to validate the computational method. Copyright © 2004 John Wiley & Sons, Ltd.     

A domain decomposition approach to compute wave breaking (wave‐breaking flows)

A heterogeneous domain decomposition approach is followed to simulate the unsteady wavy flow generated by a body moving beneath a free surface. Attention being focused on complex free surface configurations, including wave‐breaking phenomena, a two‐fluid viscous flow model is used in the free surface region to capture the air–water interface (via a level‐set technique), while a potential flow approximation is adopted to describe the flow far from the interface. Two coupling strategies are investigated, differing in the transmission conditions. Both the adopted approaches make use of the inviscid velocity field as boundary condition in the Navier–Stokes solution. For validation purposes, two different two‐dimensional non‐breaking flows are simulated. Domain decomposition results are compared with both fully viscous and fully inviscid results, obtained by solving the corresponding equations in the whole fluid domain, and with available experimental data. Finally, the unsteady evolution of a steep breaking wave is followed and some of the physical phenomena, experimentally observed, are reproduced. Copyright © 2003 John Wiley & Sons, Ltd.     

A level set technique applied to unsteady free surface flows

An unsteady Navier–Stokes solver for incompressible fluid is coupled with a level set approach to describe free surface motions. The two‐phase flow of air and water is approximated by the flow of a single fluid whose properties, such as density and viscosity, change across the interface. The free surface location is captured as the zero level of a distance function convected by the flow field. To validate the numerical procedure, two classical two‐dimensional free surface problems in hydrodynamics, namely the oscillating flow in a tank and the waves generated by the flow over a bottom bump, are studied in non‐breaking conditions, and the results are compared with those obtained with other numerical approaches. To check the capability of the method in dealing with complex free surface configurations, the breaking regime produced by the flow over a high bump is analyzed. The analysis covers the successive stages of the breaking phenomenon: the steep wave evolution, the falling jet, the splash‐up and the air entrainment. In all phases, numerical results qualitatively agree with the experimental observations. Finally, to investigate a flow in which viscous effects are relevant, the numerical scheme is applied to study the wavy flow past a submerged hydrofoil. Copyright © 2001 John Wiley & Sons, Ltd.     

Unified one‐fluid formulation for incompressible flexible solids and multiphase flows: Application to hydrodynamics using the immersed structural potential method (ISPM)

In this paper, we present a two‐dimensional computational framework for the simulation of fluid‐structure interaction problems involving incompressible flexible solids and multiphase flows, further extending the application range of classical immersed computational approaches to the context of hydrodynamics. The proposed method aims to overcome shortcomings such as the restriction of having to deal with similar density ratios among different phases or the restriction to solve single‐phase flows. First, a variation of classical immersed techniques, pioneered with the immersed boundary method (IBM), is presented by rearranging the governing equations, which define the behaviour of the multiple physics involved. The formulation is compatible with the “one‐fluid” formulation for two‐phase flows and can deal with large density ratios with the help of an anisotropic Poisson solver. Second, immersed deformable structures and fluid phases are modelled in an identical manner except for the computation of the deviatoric stresses. The numerical technique followed in this paper builds upon the immersed structural potential method developed by the authors, by adding a level set–based method for the capturing of the fluid‐fluid interfaces and an interface Lagrangian‐based meshless technique for the tracking of the fluid‐structure interface. The spatial discretisation is based on the standard marker‐and‐cell method used in conjunction with a fractional step approach for the pressure/velocity decoupling, a second‐order time integrator, and a fixed‐point iterative scheme. The paper presents a wide d range of two‐dimensional applications involving multiphase flows interacting with immersed deformable solids, including benchmarking against both experimental and alternative numerical schemes.     

Non-conservative pressure-based compressible formulation for multiphase flows with heat and mass transfer

A pressure-based compressible multiphase flow solver has been developed based on non-conservative discretization of the mixture continuity equation. The formulation is an extension of the single phase incompressible pressure-correction approach, such that it can be applied to both two-phase flows using interface resolving methods and general n-phase ensemble-averaged mixture flows. The formulation is currently presented with the single pressure and single temperature assumption, but extension to multiple temperatures is straightforward. A robust treatment of phase change allows the method to model conditions with rapid phase change such as expansion through nozzles and valves. The method has been validated thoroughly using canonical single phase problems such as the shock tube, tank filling and sudden valve closure problems. Multiphase flow validation has been carried out for sound propagation in mixtures using the ensemble-averaged model and pressure wave transmission and reflection across an air-water interface, using the level set interface tracking method. The method has been used to study sound propagation in saturated steam-water systems under thermodynamic non-equilibrium, where the expected drastic reduction in the speed of sound is reproduced. Finally the method is applied to the problem of critical (choked) flow in a nozzle for a saturated steam-water system.     

A geometry‐based level set method for curvilinear overset grids with application to ship hydrodynamics

In a previous work (Int. J. Numer. Meth. Fluids 2007; 55 :867–897), we presented a two‐phase level set method to simulate air/water turbulent flows using curvilinear body‐fitted grids for ship hydrodynamics problems. This two‐phase level set method explicitly enforces jump conditions across the interface, thus resulting in a fully coupled representation of the air/water flow. Though the method works well with multiblock curvilinear grids, severe robustness problems were found when attempting to use it with overset grids. The problem was tracked to small unphysical level set discontinuities across the overset grids with large differences in curvature. Though negligible for single‐phase approaches, the problem magnifies with large density differences between the phases, causing computation failures. In this paper, we present a geometry‐based level set method for curvilinear overset grids that overcomes these difficulties. The level set transport and reinitialization equations are not discretized along grid coordinates, but along the upwind streamline and level set gradient directions, respectively. The method is essentially an unstructured approach that is transparent to the differences between overset grids, but still the discretization is under the framework of a finite differences approach. As a result, significant improvements in robustness and to a less extent in accuracy are achieved for the level set function interpolation between overset grids, especially with big differences in grid curvature. Example tests are shown for the case of bow breaking waves around the surface combatant model David Taylor Model Basin (DTMB) 5415 and for the steady‐state ONR Tumblehome DTMB 5613 with superstructure. In the first case, the results are compared against experimental data available and in the second against results of a semi‐coupled method. Copyright © 2011 John Wiley & Sons, Ltd.     

A discontinuous Galerkin‐based sharp‐interface method to simulate three‐dimensional compressible two‐phase flow

A numerical method for the simulation of compressible two‐phase flows is presented in this paper. The sharp‐interface approach consists of several components: a discontinuous Galerkin solver for compressible fluid flow, a level‐set tracking algorithm to follow the movement of the interface and a coupling of both by a ghost‐fluid approach with use of a local Riemann solver at the interface. There are several novel techniques used: the discontinuous Galerkin scheme allows locally a subcell resolution to enhance the interface resolution and an interior finite volume Total Variation Diminishing (TVD) approximation at the interface. The level‐set equation is solved by the same discontinuous Galerkin scheme. To obtain a very good approximation of the interface curvature, the accuracy of the level‐set field is improved and smoothed by an additional P_NP_M‐reconstruction. The capabilities of the method for the simulation of compressible two‐phase flow are demonstrated for a droplet at equilibrium, an oscillating ellipsoidal droplet, and a shock‐droplet interaction problem at Mach 3. Copyright © 2015 John Wiley & Sons, Ltd.     

Convected level set method for the numerical simulation of fluid buckling

‘Fluid buckling’ is a phenomenon characterized mainly by the existence of fluid toroidal oscillations during flow. It appears when a high viscosity fluid flows vertically against a flat surface and may occur in industrial applications, as in injection molding of a propergol in complex‐shaped cavities. These coiling or folding oscillations appear during the mold filling stage, leading to air entrapment. To understand and to model this free surface flow problem, a convected level set method is proposed. First, a sinus filter is applied to the distance function to get a smooth truncation far from the interface. Second, the reinitialization is embedded in the transport equation model, avoiding it as a separate step during calculation. In order to validate the method, numerical results are presented on classical interface capturing benchmarks. Finally, results are shown on two‐dimensional and three‐dimensional viscous jet buckling problems. Copyright © 2010 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 SPH multi‐fluid model based on quasi buoyancy for interface stabilization up to high density ratios and realistic wave speed ratios

We introduce a smoothed particle hydrodynamics (SPH) concept for the stabilization of the interface between 2 fluids. It is demonstrated that the change in the pressure gradient across the interface leads to a force imbalance. This force imbalance is attributed to the particle approximation implicit to SPH. To stabilize the interface, a pressure gradient correction is proposed. In this approach, the multi‐fluid pressure gradients are related to the (gravitational and fluid) accelerations. This leads to a quasi‐buoyancy correction for hydrostatic (stratified) flows, which is extended to nonhydrostatic flows. The result is a simple density correction that involves no parameters or coefficients. This correction is included as an extra term in the SPH momentum equation. The new concept for the stabilization of the interface is explored in 5 case studies and compared with other multi‐fluid models. The first case is the stagnant flow in a tank: The interface remains stable up to density ratios of 1:1000 (typical for water and air), in combination with artificial wave speed ratios up to 1:4. The second and third cases are the Rayleigh‐Taylor instability and the rising bubble, where a reasonable agreement between SPH and level‐set models is achieved. The fourth case is an air flow across a water surface up to density ratios of 1:100, artificial wave speed ratios of 1:4, and high air velocities. The fifth case is about the propagation of internal gravity waves up to density ratios of 1:100 and artificial wave speed ratios of 1:4. It is demonstrated that the quasi‐buoyancy model may be used to stabilize the interface between 2 fluids up to high density ratios, with real (low) viscosities and more realistic wave speed ratios than achieved by other weakly compressible SPH multi‐fluid models. Real wave speed ratios can be achieved as long as the fluid velocities are not very high. Although the wave speeds may be artificial in many cases, correct and realistic wave speed ratios are essential in the modelling of heat transfer between 2 fluids (eg, in engineering applications such as gas turbines).     

Optimization of unsteady incompressible Navier–Stokes flows using variational level set method

This paper presents the optimization of unsteady Navier–Stokes flows using the variational level set method. The solid–liquid interface is expressed by the level set function implicitly, and the fluid velocity is constrained to be zero in the solid domain. An optimization problem, which is constrained by the Navier–Stokes equations and a fluid volume constraint, is analyzed by the Lagrangian multiplier based adjoint approach. The corresponding continuous adjoint equations and the shape sensitivity are derived. The level set function is evolved by solving the Hamilton–Jacobian equation with the upwind finite difference method. The optimization method can be used to design channels for flows with or without body forces. The numerical examples demonstrate the feasibility and robustness of this optimization method for unsteady Navier–Stokes flows.Copyright © 2012 John Wiley & Sons, Ltd.     

An adaptive numerical scheme for solving incompressible 2‐phase and free‐surface flows

In this paper, we present a numerical scheme for solving 2‐phase or free‐surface flows. Here, the interface/free surface is modeled using the level‐set formulation, and the underlying mesh is adapted at each iteration of the flow solver. This adaptation allows us to obtain a precise approximation for the interface/free‐surface location. In addition, it enables us to solve the time‐discretized fluid equation only in the fluid domain in the case of free‐surface problems. Fluids here are considered incompressible. Therefore, their motion is described by the incompressible Navier‐Stokes equation, which is temporally discretized using the method of characteristics and is solved at each time iteration by a first‐order Lagrange‐Galerkin method. The level‐set function representing the interface/free surface satisfies an advection equation that is also solved using the method of characteristics. The algorithm is completed by some intermediate steps like the construction of a convenient initial level‐set function (redistancing) as well as the construction of a convenient flow for the level‐set advection equation. Numerical results are presented for both bifluid and free‐surface problems.     

An explicit level set approach for generalized shape optimization of fluids with the lattice Boltzmann method

This study is concerned with a generalized shape optimization approach for finding the geometry of fluidic devices and obstacles immersed in flows. Our approach is based on a level set representation of the fluid–solid interface and a hydrodynamic lattice Boltzmann method to predict the flow field. We present an explicit level set method that does not involve the solution of the Hamilton–Jacobi equation and allows using standard nonlinear programming methods. In contrast to previous works, the boundary conditions along the fluid–structure interface are enforced by second‐order accurate interpolation schemes, overcoming shortcomings of flow penalization methods and Brinkman formulations frequently used in topology optimization. To ensure smooth boundaries and mesh‐independent results, we introduce a simple, computationally inexpensive filtering method to regularize the level set field. Furthermore, we define box constraints for the design variables that guarantee a continuous evolution of the boundaries. The features of the proposed method are studied by two numeric examples of two‐dimensional steady‐state flow problems. Copyright © 2009 John Wiley & Sons, Ltd.     

A direct reinitialization approach of level‐set/splitting finite element method for simulating incompressible two‐phase flows

Computation of a moving interface by the level‐set (LS) method typically requires reinitialization of LS function. An inaccurate execution of reinitialization results in incorrect free surface capturing and thus errors such as mass gain/loss so that an accurate and robust reinitialization process in the LS method is essential for the simulation of free surface flows. In the present study, we pursue further development of the reinitialization process, which directly corrects the LS function after advection is carried out by using the normal vector to the interface instead of solving the reinitialization equation of hyperbolic type. The Taylor–Galerkin method is adopted to discretize the advection equation of the LS function and the P1P1 splitting finite element method is applied to solve the Navier–Stokes equation. The proposed algorithm is validated with the well‐known benchmark problems, i.e. stretching of a circular fluid element, time‐reversed single‐vortex, solitary wave propagation, broken dam flow and filling of a container. The simulation results of these flows are in good agreement with previously existing experimental and numerical results. Copyright © 2010 John Wiley & Sons, Ltd.     

A hybrid Lagrangian–Eulerian particle‐level set method for numerical simulations of two‐fluid turbulent flows

A coupled Lagrangian interface‐tracking and Eulerian level set (LS) method is developed and implemented for numerical simulations of two‐fluid flows. In this method, the interface is identified based on the locations of notional particles and the geometrical information concerning the interface and fluid properties, such as density and viscosity, are obtained from the LS function. The LS function maintains a signed distance function without an auxiliary equation via the particle‐based Lagrangian re‐initialization technique. To assess the new hybrid method, numerical simulations of several ‘standard interface‐moving’ problems and two‐fluid laminar and turbulent flows are conducted. The numerical results are evaluated by monitoring the mass conservation, the turbulence energy spectral density function and the consistency between Eulerian and Lagrangian components. The results of our analysis indicate that the hybrid particle‐level set method can handle interfaces with complex shape change, and can accurately predict the interface values without any significant (unphysical) mass loss or gain, even in a turbulent flow. The results obtained for isotropic turbulence by the new particle‐level set method are validated by comparison with those obtained by the ‘zero Mach number’, variable‐density method. For the cases with small thermal/mass diffusivity, both methods are found to generate similar results. Analysis of the vorticity and energy equations indicates that the destabilization effect of turbulence and the stability effect of surface tension on the interface motion are strongly dependent on the density and viscosity ratios of the fluids. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical simulation of cavitating flow in 2D and 3D inducer geometries

A computational method is proposed to simulate 3D unsteady cavitating flows in spatial turbopump inducers. It is based on the code FineTurbo, adapted to take into account two‐phase flow phenomena. The initial model is a time‐marching algorithm devoted to compressible flow, associated with a low‐speed preconditioner to treat low Mach number flows. The presented work covers the 3D implementation of a physical model developed in LEGI for several years to simulate 2D unsteady cavitating flows. It is based on a barotropic state law that relates the fluid density to the pressure variations. A modification of the preconditioner is proposed to treat efficiently as well highly compressible two‐phase flow areas as weakly compressible single‐phase flow conditions. The numerical model is applied to time‐accurate simulations of cavitating flow in spatial turbopump inducers. The first geometry is a 2D Venturi type section designed to simulate an inducer blade suction side. Results obtained with this simple test case, including the study of its general cavitating behaviour, numerical tests, and precise comparisons with previous experimental measurements inside the cavity, lead to a satisfactory validation of the model. A complete three‐dimensional rotating inducer geometry is then considered, and its quasi‐static behaviour in cavitating conditions is investigated. Numerical results are compared to experimental measurements and visualizations, and a promising agreement is obtained. Copyright © 2004 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.     

Numerical prediction of droplet dynamics in turbulent flow,using the level set method

In the present article, the droplet dynamics in turbulent flow is numerically predicted. The modelling is based on an interfacial marker-level set (IMLS) method, coupled with the Reynolds-averaged Navier–Stokes (RANS) equations to predict the dynamics of turbulent two-phase flow. The governing equations for time-dependent, two-dimensional and incompressible two-phase flow are described in both phases and solved separately using a control volume approach on structured cell-centred collocated grids. The topological changes of the interface are predicted by applying the level set approach. The kinematic and dynamic conditions on the interface separating the two phases are satisfied. The numerical method proposed is validated against a well-known computational fluid dynamics problem. Further, the deformation and breakup of a single droplet either suddenly moved in air or exposed to turbulent stream are numerically investigated. In general, the developed numerical method demonstrates remarkable capability in predicting the characteristics of complex turbulent two-phase flows.