A Stable Semi-Implicit Method for Free Surface Flows

2005 ◽  
Vol 73 (6) ◽  
pp. 940-947 ◽  
Author(s):  
Cassio M. Oishi ◽  
José A. Cuminato ◽  
Valdemir G. Ferreira ◽  
Murilo F. Tomé ◽  
Antonio Castelo ◽  
...  

The present work is concerned with a semi-implicit modification of the GENSMAC method for solving the two-dimensional time-dependent incompressible Navier-Stokes equations in primitive variables formulation with a free surface. A projection method is employed to uncouple the velocity components and pressure, thus allowing the solution of each variable separately (a segregated approach). The viscous terms are treated by the implicit backward method in time and a centered second order method in space, and the nonlinear convection terms are explicitly approximated by the high order upwind variable-order nonoscillatory scheme method in space. The boundary conditions at the free surface couple the otherwise segregated velocity and pressure fields. The present work proposes a method that allows the segregated solution of free surface flow problems to be computed by semi-implicit schemes that preserve the stability conditions of the related coupled semi-implicit scheme. The numerical method is applied to both the simulation of free surface and to confined flows. The numerical results demonstrate that the present technique eliminates the parabolic stability restriction required by the original explicit GENSMAC method, and also found in segregated semi-implicit methods with time-lagged boundary conditions. For low Reynolds number flows, the method is robust and very efficient when compared to the original GENSMAC method.

Author(s):  
A. Ganguly ◽  
V. Shigunov ◽  
O. Turan

A finite volume method with a multiphase type free surface description is employed to calculate the flow around ships in shallow and restricted channels. The flows at critical and supercritical depth Froude numbers (Fnd = 1.0 and Fnd = 1.18) are calculated for Series–60 monohull and a medium speed catamaran. A steady state solution for Reynolds-averaged Navier-Stokes equations with a k-ε turbulence model is obtained by time marching. Computed wave profiles are in good agreement with model tests in the near field of the ship. The computed and measured resistance agree fairly well.


Author(s):  
Vedanth Srinivasan ◽  
De Ming Wang

This paper presents a numerical method that couples the incompressible Navier-Stokes equations with the Volume of Fluid method in a Cartesian co-ordinate system for tracking immiscible interfaces in multiple dimensions. The governing equations are discretized based on a finite volume method on a non-staggered fixed grid. The free surface flow problem is solved as a single phase flow system in which the free surface is captured using a Switching Technique for Advection and Capturing of Surfaces (STACS) scheme. The effects of surface tension at the interfaces are treated using a Continuum Surface Force (CSF) model. The pressure velocity coupling is achieved using a SIMPLE strategy. The coupled system, implemented in the commercial CFD software, AVL FIRE/SWIFT, is applied to a two dimensional dam breaking problem. The simulation results reveal a multitude of phenomena such as, free surface vortex generation, air entrapment and splashing of the liquid surge front. The computational results are in good agreement with experimental data, wherever available. The effects of time and grid resolution on the solution behavior are elaborated in detail. Different convection schemes are tested and the current method is compared to another existing interface capturing methodology.


2005 ◽  
Vol 4 (2) ◽  
Author(s):  
C. M. Oishi ◽  
J. A. Cuminato ◽  
V. G. Ferreira ◽  
M. F. Tomé ◽  
A. Castelo ◽  
...  

The present work is concerned with a numerical method for solving the two-dimensional time-dependent incompressible Navier-Stokes equations in the primitive variables formulation. The diffusive terms are treated by Implicit Backward and Crank-Nicolson methods, and the non-linear convection terms are, explicitly, approximated by the high order upwind VONOS (Variable-Order Non-Oscillatory Scheme) scheme. The boundary conditions for the pressure field at the free surface are treated implicitly, and for the velocity field explicitly. The numerical method is then applied to the simulation of free surface and confined flows. The numerical results show that the present technique eliminates the stability restriction in the original explicit method. For low Reynolds number flow dynamics, the method is robust and produces numerical results that compare very well with the analytical solutions.


2005 ◽  
Vol 4 (2) ◽  
pp. 106
Author(s):  
C. M. Oishi ◽  
J. A. Cuminato ◽  
V. G. Ferreira ◽  
M. F. Tomé ◽  
A. Castelo ◽  
...  

The present work is concerned with a numerical method for solving the two-dimensional time-dependent incompressible Navier-Stokes equations in the primitive variables formulation. The diffusive terms are treated by Implicit Backward and Crank-Nicolson methods, and the non-linear convection terms are, explicitly, approximated by the high order upwind VONOS (Variable-Order Non-Oscillatory Scheme) scheme. The boundary conditions for the pressure field at the free surface are treated implicitly, and for the velocity field explicitly. The numerical method is then applied to the simulation of free surface and confined flows. The numerical results show that the present technique eliminates the stability restriction in the original explicit method. For low Reynolds number flow dynamics, the method is robust and produces numerical results that compare very well with the analytical solutions.


2012 ◽  
Vol 1 (33) ◽  
pp. 71
Author(s):  
Takahide Honda ◽  
Peter Wellens ◽  
Marcel Van Gent

COMFLOW is a general 3D free-surface flow solver. The numerical method is based on the Navier-Stokes equations in a porous medium, with additional force terms to represent the (turbulent) interaction of the flow with the medium. The free surface is displaced by means of the Volume-Of-Fluid method. The main objective in this paper is to validate the permeable flow model in 3D. Tailor-made physical model tests were performed for this purpose. In the experiment surface elevations are measured inside and around a permeable structure with 18 wave gauges in total. The measurements are represented well by the simulation results.


Fluids ◽  
2021 ◽  
Vol 6 (10) ◽  
pp. 350
Author(s):  
Georgios A. Leftheriotis ◽  
Iason A. Chalmoukis ◽  
Guillermo Oyarzun ◽  
Athanassios A. Dimas

An advanced numerical model is presented for the simulation of wave-induced free-surface flow, utilizing an efficient hybrid parallel implementation. The model is based on the solution of the Navier–Stokes equations using large-eddy simulation of large-scale coastal free-surface flows. The three-dimensional immersed boundary method was used for the enforcement of the no-slip boundary condition on the bed surface. The water-air interface was tracked using the level-set method. The numerical model was effectively validated against laboratory measurements involving wave propagation over a flatbed with an elliptical shoal, whose presence induces combined wave refraction and diffraction phenomena. The parallel implementation of the model enabled the efficient simulation of depth-resolved, wave-induced, three-dimensional, free-surface flow; the model parallel efficiency and strong scaling are quantitatively demonstrated.


2007 ◽  
Vol 51 (01) ◽  
pp. 47-64
Author(s):  
James C. Huan ◽  
Thomas T. Huang

A fast turnaround and an accurate computational fluid dynamics (CFD) approach for ship total resistance prediction is developed. The approach consists of a nonlinear free surface potential flow solver (PShip code) with a wet-or-dry transom stern model, and a Reynolds-averaged Navier-Stokes (RANS) equation solver that solves viscous free surface flow with a prescribed free surface given from the PShip. The prescribed free surface RANS predicts a viscous correction to the pressure resistance (viscous form) and viscous flow field around the hull. The viscous free surface flow solved this way avoids the time-consuming RANS iterations to resolve the free surface profile. The method, however, requires employing a flow characteristic-based nonreflecting boundary condition at the free surface. The approach can predict the components of ship resistance, the associated wave profile around the hull, and the sinkage and trim of the ship. Validation of the approach is presented with Wigley, Series 60 (CB = 0.6), and NSWCCD Model 5415 hulls. An overall accuracy of ±2% for ship total resistance prediction is achieved. The approach is applied to evaluating the effects of a stern flap on a DD 968 model on ship performance. An empirical viscous form resistance formula is also devised for a quick ship total resistance estimate.


Author(s):  
I. L. Collings

AbstractSolutions are found to two cusp-like free-surface flow problems involving the steady motion of an ideal fluid under the infinite-Froude-number approximation. The flow in each case is due to a submerged line source or sink, in the presence of a solid horizontal base.


Author(s):  
Iraj Saeedpanah ◽  
M. Shayanfar ◽  
E. Jabbari ◽  
Mohammad Haji Mohammadi

Free surface flows are frequently encountered in hydraulic engineering problems including water jets, weirs and around gates. An iterative solution to the incompressible two-dimensional vertical steady Navier-Stokes equations, comprising momentum and continuity equations, is used to solve for the priori unknown free surface, the velocity and the pressure fields. The entire water body is covered by a unstructured finite element grid which is locally refined. The dynamic boundary condition is imposed for the free surface where the pressure vanishes. This procedure is done continuously until the normal velocities components vanish. To overcome numerical errors and oscillations encountering in convection terms, the SUPG (streamline upwinding Petrov-Galerkin) method is applied. The solution method is tested for different discharges onto a standard spillway geometries. The results shows good agreement with available experimental data.


Sign in / Sign up

Export Citation Format

Share Document