Simulation of free surface flows using volume of fluid method and genetic algorithm

2014 ◽  
Vol 16 (5) ◽  
pp. 1110-1124 ◽  
Author(s):  
N. Soleiman Beygi ◽  
H. Hakimzadeh ◽  
M. R. Chenaglou

In this paper, details of a numerical model development for simulation of fluid flows with moving free surface are presented. The unsteady incompressible Navier–Stokes equations on a fixed grid system are used to obtain velocity and pressure values in the computational domain and volume of fluid (VOF) method is used to determine free surface location. In order to reduce numerical smearing and increase accuracy of numerical modeling of fluid flow with moving free surface, a new free surface-tracking method is proposed. The proposed method is a combination of genetic algorithm and free surface tracking method based on donor and acceptor scheme. The specification of the new combinational method can be summarized in determining orientation vector and plane constant to represent the free surface orientation in each cell. The proposed algorithm can be easily used in any unstructured grids. In this method, the fluid flow equations are explicitly discretized with the finite volume method and the projection method is used to determine the velocity and pressure magnitude in computational domain. Validity of the new solution algorithm is demonstrated through its application to the dam break and the bore motion examples.

2012 ◽  
Vol 134 (9) ◽  
Author(s):  
Sajad Pooyan ◽  
Mohammad Passandideh-Fard

In this paper, a numerical model is developed that can simulate the unsteady axisymmetric free-surface flow of a perfectly conductive liquid under an electrostatic field. The effect of the electrostatic field is modeled by a force distributed on the liquid free surface. Assuming the liquid as a perfect conductor makes it possible to reduce the general electromagnetic equations to electrostatic equations. The Navier–Stokes equations are solved to find the velocity and pressure fields. The free surface advection and reconstruction are performed based on the volume-of-fluid method using Youngs’ algorithm. To evaluate the effect of the electric field on the free surface, the electrostatic potential is first solved for the entire computational domain. Next, the electric field intensity and the surface density of the electric charge are calculated on the free surface after which the electric force can be determined. The computational method for treating this force is similar to that of the surface tension using the continuum surface force method. The developed model is validated by a comparison between the calculated results with those of the analytics as well as experiments for an electrowetting scenario.


Author(s):  
Hans Bihs ◽  
Mayilvahanan Alagan Chella ◽  
Arun Kamath ◽  
Øivind Asgeir Arntsen

For the stability of offshore structures, such as offshore wind foundations, extreme wave conditions need to be taken into account. Waves from extreme events are critical from the design perspective. In a numerical wave tank, extreme waves can be modeled using focused waves. Here, linear waves are generated from a wave spectrum. The wave crests of the generated waves coincide at a preselected location and time. Focused wave generation is implemented in the numerical wave tank module of REEF3D, which has been extensively and successfully tested for various wave hydrodynamics and wave–structure interaction problems in particular and for free surface flows in general. The open-source computational fluid dynamics (CFD) code REEF3D solves the three-dimensional Navier–Stokes equations on a staggered Cartesian grid. Higher order numerical schemes are used for time and spatial discretization. For the interface capturing, the level set method is selected. In order to test the generated waves, the time series of the free surface elevation are compared with experimental benchmark cases. The numerically simulated free surface elevation shows good agreement with experimental data. In further computations, the impact of the focused waves on a vertical circular cylinder is investigated. A breaking focused wave is simulated and the associated kinematics is investigated. Free surface flow features during the interaction of nonbreaking focused waves with a cylinder and during the breaking process of a focused wave are also investigated along with the numerically captured free surface.


1995 ◽  
Vol 117 (4) ◽  
pp. 683-690 ◽  
Author(s):  
Peter E. Raad ◽  
Shea Chen ◽  
David B. Johnson

A new method of calculating the pressure field in the simulation of two-dimensional, unsteady, incompressible, free surface fluid flow by use of a marker and cell method is presented. A critical feature of the new method is the introduction of a finer mesh of cells in addition to the regular mesh of finite volume cells. The smaller (micro) cells are used only near the free surface, while the regular (macro) cells are used throughout the computational domain. The movement of the free surface is accomplished by the use of massless surface markers, while the discrete representation of the free surface for the purpose of the application of pressure boundary conditions is accomplished by the use of micro cells. In order to exploit the advantages offered by micro cells, a new general equation governing the pressure field is derived. Micro cells also enable the identification and treatment of multiple points on the free surface in a single surface macro cell as well as of points on the free surface that are located in a macro cell that has no empty neighbors. Both of these situations are likely to occur repeatedly in a free surface fluid flow simulation, but neither situation has been explicitly taken into account in previous marker and cell methods. Numerical simulation results obtained both with and without the use of micro cells are compared with each other and with theoretical solutions to demonstrate the capabilities and validity of the new method.


2018 ◽  
Vol 2018 ◽  
pp. 1-9 ◽  
Author(s):  
Deepak Garg ◽  
Antonella Longo ◽  
Paolo Papale

This work aims to develop a numerical wave tank for viscous and inviscid flows. The Navier-Stokes equations are solved by time-discontinuous stabilized space-time finite element method. The numerical scheme tracks the free surface location using fluid velocity. A segregated algorithm is proposed to iteratively couple the fluid flow and mesh deformation problems. The numerical scheme and the developed computer code are validated over three free surface problems: solitary wave propagation, the collision between two counter moving waves, and wave damping in a viscous fluid. The benchmark tests demonstrate that the numerical approach is effective and an attractive tool for simulating viscous and inviscid free surface flows.


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.


Author(s):  
Gerasimos A. Kolokythas ◽  
Athanassios A. Dimas

In the present study, numerical simulations of the free-surface flow, developing by the propagation of nonlinear water waves over a rippled bottom, are performed assuming that the corresponding flow is two-dimensional, incompressible and viscous. The simulations are based on the numerical solution of the Navier-Stokes equations subject to the fully-nonlinear free-surface boundary conditions and the suitable bottom, inflow and outflow boundary conditions. The equations are properly transformed so that the computational domain becomes time-independent. For the spatial discretization, a hybrid scheme with finite-differences and Chebyshev polynomials is applied, while a fractional time-step scheme is used for the temporal discretization. A wave absorption zone is placed at the outflow region in order to efficiently minimize reflection of waves by the outflow boundary. The numerical model is validated by comparison to the analytical solution for the laminar, oscillatory, current flow which develops a uniform boundary layer over a horizontal bottom. For the propagation of finite-amplitude waves over a rigid rippled bed, the case with wavelength to water depth ratio λ/d0 = 6 and wave height to wavelength ratio H0/λ = 0.05 is considered. The ripples have parabolic shape, while their dimensions — length and height — are chosen accordingly to fit laboratory and field data. Results indicate that the wall shear stress over the ripples and the form drag forces on the ripples increase with increasing ripple height, while the corresponding friction force is insensitive to this increase. Therefore, the percentage of friction in the total drag force decreases with increasing ripple height.


Author(s):  
Y. G. Chen ◽  
W. G. Price ◽  
P. Temarel

This investigation continues the development of an anti-diffusive volume of fluid method [1] by improving accuracy through the addition of an artificial diffusion term, with a negative diffusion coefficient, to the original advection equation describing the evolution of the fluid volume fraction. The advection and diffusion processes are split into a set of two partial differential equations (PDEs). The improved anti-diffusive Volume of Fluid (VOF) method is coupled with a two-fluid flow solver to predict free surface flows and illustrated by examples given in two-dimensional flows. The first numerical example is a solitary wave travelling in a tank. The second example is a plunging wave generated by flow over a submerged obstacle of prescribed shape on a horizontal floor. The computational results are validated against available experimental data.


2011 ◽  
Vol 233-235 ◽  
pp. 701-706
Author(s):  
Bing Tao Zhao ◽  
Yi Xin Zhang ◽  
Kai Bin Xiong

The numerical simulation of the fluid flow is presented by CFD technique to characterize the flow pattern of cyclone spray scrubber. In this process, the Reynolds-averaged Navier-Stokes equations with the Reynolds stress turbulence model (RSM) for fluid flow are solved by use of the finite volume method based on the SIMPLE pressure correction algorithm in the fluid computational domain. According to the computational results, the tangential velocity, axial velocity and turbulence intensity of the gas flow are addressed in the different flowrate. The results indicate that the CFD method can effectively reveal the mechanism of gas flow in the cyclone spray scrubber.


Sign in / Sign up

Export Citation Format

Share Document