scholarly journals An FEM-Level-set Numerical Model for Potential Flow with Free Surface

2015 ◽  
Vol 126 ◽  
pp. 237-241
Author(s):  
Longfei Cong ◽  
Bin Teng
2002 ◽  
Vol 456 ◽  
pp. 113-135 ◽  
Author(s):  
H. KAGEMOTO ◽  
M. MURAI ◽  
M. SAITO ◽  
B. MOLIN ◽  
š. MALENICA

A row of fifty identical, truncated vertical cylinders is submitted to regular head waves, with wave periods in a narrow range around the period of the so-called Neumann trapped mode. The free-surface elevation is measured at 14 locations along the array. Response amplitude operators of the free-surface motion are compared with numerical predictions from a potential flow model. Resonance effects, at wave periods equal to or larger than the critical one, are found to be much less than given by the numerical model. It is advocated that these discrepancies are due to dissipative effects taking place in the boundary layers at the cylinder walls. An artificial means is devised to incorporate dissipation in the potential flow model, whereby the cylinder walls are made slightly porous; the inward normal velocity of the flow is related to the dynamic pressure. The coefficient of proportionality is based on existing knowledge for circular cylinders in oscillatory flows. With this modification in the numerical code, excellent agreement is obtained with the experiments. The numerical model is further used for the case of a very long array composed of 1000 cylinders; it is found that with dissipation at the cylinder walls, the wave action steadily decreases along the array, even for wave periods substantially larger than the critical one. On the other hand, at wave periods less than the critical one, dissipation plays a negligible role; the observed decay is solely due to diffraction effects. Implications of these results for very large structures such as column-supported floating airports are discussed. In particular, it is concluded that scale effects may be an important issue in the experimental analysis of such multi-column structures.


Author(s):  
Erlend Liavåg Grotle ◽  
Hans Bihs ◽  
Eilif Pedersen ◽  
Vilmar Æsøy

In this paper, numerical simulations of non-linear sloshing in rectangular tanks are presented. Model implementations in the open source software REEF3D are tested and results compared with experimental data. Three different conditions are compared with experiments in 2D. First, the free surface time-evolution is compared for both linear and non-linear sloshing. In the last case, video images from the SPHERIC project are compared with simulations images of the free surface. A condition with lateral wave impacts in sloshing, with a frequency closer to the natural frequency of the first mode, can be found in this case. The non-linear sloshing, case 2, is also simulated in 3D. The numerical model is solving the RANS equations with the k-ω turbulence model. The level set method is used to capture the interface. Higher order discretization schemes are implemented to handle time-evolution and convective fluxes. A ghost cell method is used to account for solid boundaries and multiple grids for parallel computations. It is found that the limiting factor for the eddy-viscosity has significant influence in case 2 and 3. As the sloshing becomes more violent, the increased strain at the gas-liquid interface overproduces turbulence energy with unrealistically high damping of the motion. 3D simulations are only performed in case 2, which shows slightly better comparison than with 2D. Due to non-linearities and small damping, the time to reach steady-state may take several cycles, but no information is given in the paper [1]. The last case shows promising results for the global motion. As expected, the break up of the liquid surface makes it difficult to resolve each phase. But overall, the numerical model predicts the sloshing motion reasonably well.


2016 ◽  
Author(s):  
Arun Kamath ◽  
Hans Bihs ◽  
Øivind A. Arnsten

Many offshore constructions and operations involve water impact problems such as water slamming onto a structure or free fall of objects with subsequent water entry and emergence. Wave slamming on semi-submersibles, vertical members of jacket structures, crane operation of a diving bell and dropping of free fall lifeboats are some notable examples. The slamming and water entry problems lead to large instantaneous impact pressures on the structure, accompanied with complex free surface deformations. These need to be studied in detail in order to obtain a better understanding of the fluid physics involved and develop safe and economical design. In the special case of free-fall lifeboats, model testing can be expensive and time consuming. Here, numerical modelling can make useful contributions to the design process. The slamming of a free falling body into water involves several complex hydrodynamic features after its free-fall such as water entry, submergence into water and resurfacing. The water entry and submergence lead to formation of water jets and air cavities in the water resulting in large impact forces on the object. In order to evaluate the forces and hydrodynamics involved, the numerical model should be able to account for the complex free surface features, the instantaneous pressure changes around the lifeboat and accurately evaluate the loads on the lifeboat. As a step towards simulating free-fall lifeboats, water entry of a free-falling wedge into water is studied in this paper using a CFD model. The vertical velocity of the wedge during the process of free fall and water impact are calculated for different cases and the free surface deformations are captured in detail. Numerical results are compared with experimental data and a good agreement is seen. The open-source CFD model REEF3D is used in this study. The model solves the Reynolds-Averaged Navier-Stokes equations to evaluate the fluid flow. The convective terms are discretized using a 5th-order conservative finite difference WENO scheme. Time discretization is carried out using a 3rd-order Runge-Kutta scheme. Pressure discretization is carried out using Chorins projection method. The Poisson pressure equation is solved using a pre-conditioned BiCGStab algorithm. A sharp representation of the free surface is obtained using the level set method. The falling wedge is represented using the level set paradigm as well, avoiding the need for re-meshing during the simulation. Turbulence modeling is carried out using the k-ω model. Computational performance of the numerical model is improved by parallel processing using the MPI library.


2021 ◽  
Vol 36 (3) ◽  
pp. 165-176
Author(s):  
Kirill Nikitin ◽  
Yuri Vassilevski ◽  
Ruslan Yanbarisov

Abstract This work presents a new approach to modelling of free surface non-Newtonian (viscoplastic or viscoelastic) fluid flows on dynamically adapted octree grids. The numerical model is based on the implicit formulation and the staggered location of governing variables. We verify our model by comparing simulations with experimental and numerical results known from the literature.


2003 ◽  
Vol 41 (1) ◽  
pp. 110-112
Author(s):  
ZhixiaN. Cao ◽  
Rodney Day ◽  
Sarah Liriano

Author(s):  
C-E Janson

A potential-flow panel method is used to compute the waves and the lift force from surface-piercing and submerged bodies. In particular the interaction between the waves and the lift produced close to the free surface is studied. Both linear and non-linear free-surface boundary conditions are considered. The potential-flow method is of Rankine-source type using raised source panels on the free surface and a four-point upwind operator to compute the velocity derivatives and to enforce the radiation condition. The lift force is introduced as a dipole distribution on the lifting surfaces and on the trailing wake, together with a flow tangency condition at the trailing edge of the lifting surface. Different approximations for the spanwise circulation distribution at the free surface were tested for a surface-piercing wing and it was concluded that a double-model approximation should be used for low speeds while a single-model, which allows for a vortex at the free surface, was preferred at higher speeds. The lift force and waves from three surface-piercing wings, a hydrofoil and a sailing yacht were computed and compared with measurements and good agreement was obtained.


2017 ◽  
Vol 14 (06) ◽  
pp. 1750063 ◽  
Author(s):  
A. M. Hegab ◽  
S. A. Gutub ◽  
A. Balabel

This paper presents the development of an accurate and robust numerical modeling of instability of an interface separating two-phase system, such as liquid–gas and/or solid–gas systems. The instability of the interface can be refereed to the buoyancy and capillary effects in liquid–gas system. The governing unsteady Navier–Stokes along with the stress balance and kinematic conditions at the interface are solved separately in each fluid using the finite-volume approach for the liquid–gas system and the Hamilton–Jacobi equation for the solid–gas phase. The developed numerical model represents the surface and the body forces as boundary value conditions on the interface. The adapted approaches enable accurate modeling of fluid flows driven by either body or surface forces. The moving interface is tracked and captured using the level set function that initially defined for both fluids in the computational domain. To asses the developed numerical model and its versatility, a selection of different unsteady test cases including oscillation of a capillary wave, sloshing in a rectangular tank, the broken-dam problem involving different density fluids, simulation of air/water flow, and finally the moving interface between the solid and gas phases of solid rocket propellant combustion were examined. The latter case model allowed for the complete coupling between the gas-phase physics, the condensed-phase physics, and the unsteady nonuniform regression of either liquid or the propellant solid surfaces. The propagation of the unsteady nonplanar regression surface is described, using the Essentially-Non-Oscillatory (ENO) scheme with the aid of the level set strategy. The computational results demonstrate a remarkable capability of the developed numerical model to predict the dynamical characteristics of the liquid–gas and solid–gas flows, which is of great importance in many civilian and military industrial and engineering applications.


Sign in / Sign up

Export Citation Format

Share Document