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.

[1]  J. A. Sethian,et al.  Fast Marching Methods , 1999, SIAM Rev..

[2]  Frederick Stern,et al.  Unsteady RANS simulation of the ship forward speed diffraction problem , 2006 .

[3]  Emilio F. Campana,et al.  A level set technique applied to unsteady free surface flows , 2001 .

[4]  J A Sethian,et al.  A fast marching level set method for monotonically advancing fronts. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[5]  J. López,et al.  An improved PLIC-VOF method for tracking thin fluid structures in incompressible two-phase flows , 2005 .

[6]  Frederick Stern,et al.  Towing tank PIV measurement system, data and uncertainty assessment for DTMB Model 5512 , 2001 .

[7]  M. Sussman.  A second order coupled level set and volume-of-fluid method for computing growth and collapse of vapor bubbles , 2003 .

[8]  D. Drew,et al.  Theory of Multicomponent Fluids , 1998 .

[9]  Ian M. Mitchell,et al.  A hybrid particle level set method for improved interface capturing , 2002 .

[10]  D. Juric,et al.  A front-tracking method for the computations of multiphase flow , 2001 .

[11]  H. Udaykumar,et al.  Sharp interface Cartesian grid method II: A technique for simulating droplet interactions with surfaces of arbitrary shape , 2005 .

[12]  Anthony F. Molland,et al.  Quantifying the Airflow Distortion over Merchant Ships. Part II: Application of the Model Results , 2006 .

[13]  Frederick Stern,et al.  An unsteady single‐phase level set method for viscous free surface flows , 2007 .

[14]  R. Fedkiw,et al.  A Boundary Condition Capturing Method for Poisson's Equation on Irregular Domains , 2000 .

[15]  F. Menter.  Two-equation eddy-viscosity turbulence models for engineering applications , 1994 .

[16]  J. Sethian,et al.  LEVEL SET METHODS FOR FLUID INTERFACES , 2003 .

[17]  Wei-Cheng Wang.  A Jump Condition Capturing Finite Difference Scheme for Elliptic Interface Problems , 2004, SIAM J. Sci. Comput..

[18]  Anthony F. Molland,et al.  Quantifying the Airflow Distortion over Merchant Ships. Part I: Validation of a , 2006 .

[19]  Phillip Colella,et al.  An efficient second-order projection method for viscous incompressible flow , 1991 .

[20]  Mark Sussman,et al.  An Efficient, Interface-Preserving Level Set Redistancing Algorithm and Its Application to Interfacial Incompressible Fluid Flow , 1999, SIAM J. Sci. Comput..

[21]  Sangmook Shin,et al.  Internal wave computations using the ghost fluid method on unstructured grids , 2005 .

[22]  S. Osher,et al.  Level set methods: an overview and some recent results , 2001 .

[23]  Ronald Fedkiw,et al.  A Boundary Condition Capturing Method for Multiphase Incompressible Flow , 2000, J. Sci. Comput..

[24]  Ronald Fedkiw,et al.  Regular Article: The Ghost Fluid Method for Deflagration and Detonation Discontinuities , 1999 .

[25]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[26]  Matthias Klemt.  RANSE Simulation of Ship Seakeeping using Overlapping Grids , 2005 .

[27]  Shan Zhao,et al.  High order matched interface and boundary method for elliptic equations with discontinuous coefficients and singular sources , 2006, J. Comput. Phys..

[28]  D. Jacqmin.  Regular Article: Calculation of Two-Phase Navier–Stokes Flows Using Phase-Field Modeling , 1999 .

[29]  G. Syms.  Numerical Simulation of Frigate Airwakes , 2004 .

[30]  Frederick Stern,et al.  Uncertainty Assessment for Towing Tank Tests With Example for Surface Combatant DTMB Model 5415 , 2005 .

[31]  Joe F. Thompson,et al.  Numerical grid generation , 1985 .

[32]  Volker Bertram,et al.  Benchmarking of computational fluid dynamics for ship flows: the Gothenburg 2000 workshop , 2003 .

[33]  Heinz Pitsch,et al.  A ghost-fluid method for large-eddy simulations of premixed combustion in complex geometries , 2007, J. Comput. Phys..

[34]  Riccardo Broglia,et al.  On the application of the single-phase level set method to naval hydrodynamic flows , 2007 .

[35]  Frederick Stern,et al.  Unsteady RANS method for ship motions with application to roll for a surface combatant , 2006 .

[36]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[37]  K. R. Reddy,et al.  Numerical simulation of ship airwake , 2000 .

[38]  P. Wesseling,et al.  A mass‐conserving Level‐Set method for modelling of multi‐phase flows , 2005 .

[39]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[40]  J. Sethian.  Numerical algorithms for propagating interfaces: Hamilton-Jacobi equations and conservation laws , 1990 .

[41]  R Nangia,et al.  Novel Vortex Flow Devices - Columnar Vortex Generators Studies for Airwakes , 2004 .

[42]  Cheng Jian-long,et al.  Natural vegetation recovery on waste dump in opencast coalmine area , 2005, Journal of Forestry Research.

[43]  D. Fletcher,et al.  A New Volume of Fluid Advection Algorithm , 2000 .

[44]  James A. Sethian,et al.  The Fast Construction of Extension Velocities in Level Set Methods , 1999 .

[45]  Woo Il Lee,et al.  A new VOF‐based numerical scheme for the simulation of fluid flow with free surface. Part I: New free surface‐tracking algorithm and its verification , 2003 .