scholarly journals A Computational Method for Free Surface Hydrodynamics

1981 ◽  
Vol 103 (2) ◽  
pp. 136-141 ◽  
Author(s):  
C. W. Hirt ◽  
B. D. Nichols

There are numerous flow phenomena in pressure vessel and piping systems that involve the dynamics of free fluid surfaces. For example, fluid interfaces must be considered during the draining or filling of tanks, in the formation and collapse of vapor bubbles, and in seismically shaken vessels that are partially filled. To aid in the analysis of these types of flow phenomena, a new technique has been developed for the computation of complicated free-surface motions. This technique is based on the concept of a local average volume of fluid (VOF) and is embodied in a computer program for two-dimensional, transient fluid flow called SOLA-VOF. The basic approach used in the VOF technique is briefly described, and compared to other free-surface methods. Specific capabilities of the SOLA-VOF program are illustrated by generic examples of bubble growth and collapse, flows of immiscible fluid mixtures, and the confinement of spilled liquids.

Author(s):  
Jeremy Rice ◽  
Amir Faghri

A new technique, called the surface velocity correction technique (SVC), is developed to track a free surface such as a liquid-vapor interface. SVC is a computationally inexpensive, and accurate method to capture interfacial fluid phenomena. This method uses a finite volume technique to discretize the governing equations, and a semi-Legrangian mesh to locate the interface between two fluids. The effectiveness of this technique is demonstrated through several classical examples and the results are also compared to both analytical and VOF solutions. The examples include: the shape of a meniscus in a capillary tube in mechanical equilibrium, the rise of a meniscus in a capillary tube, and the instability growth of a free flowing cylindrical column of fluid.


Ocean Science ◽  
2011 ◽  
Vol 7 (6) ◽  
pp. 733-743 ◽  
Author(s):  
J. Stark ◽  
J. Gorman ◽  
M. Hennessey ◽  
F. Reseghetti ◽  
J. Willis ◽  
...  

Abstract. A new technique for determining the depth of expendable bathythermographs (XBTs) is developed. This new method uses a forward-stepping calculation which incorporates all of the forces on the XBT devices during their descent. Of particular note are drag forces which are calculated using a new drag coefficient expression. That expression, obtained entirely from computational fluid dynamic modeling, accounts for local variations in the ocean environment. Consequently, the method allows for accurate determination of depths for any local temperature environment. The results, which are entirely based on numerical simulation, are compared with the experiments of LM Sippican T-5 XBT probes. It is found that the calculated depths differ by less than 3% from depth estimates using the standard fall-rate equation (FRE). Furthermore, the differences decrease with depth. The computational model allows an investigation of the fluid flow patterns along the outer surface of the probe as well as in the interior channel. The simulations take account of complex flow phenomena such as laminar-turbulent transition and flow separation.


2011 ◽  
Vol 8 (4) ◽  
pp. 1777-1802
Author(s):  
J. Stark ◽  
J. Gorman ◽  
M. Hennessey ◽  
F. Reseghetti ◽  
J. Willis ◽  
...  

Abstract. A new technique for determining the depth of expendable bathythermographs (XBTs) is developed. This new method combines a forward-stepping calculation which incorporates all of the forces on the XBT devices during their descent. Of particular note are drag forces which are calculated using a new drag coefficient expression. That expression, obtained entirely from computational fluid dynamic modeling, accounts for local variations in the ocean environment. Consequently, the method allows for accurate determination of depths for any local temperature environment. The results, which are entirely based on numerical simulation, are compared with an experimental descent of an LM-Sippican T-5 XBT. It is found that the calculated depths differ by less than 3 % from depth estimates using the industry standard FRE. Furthermore, the differences decrease with depth. The computational model allows an investigation of the fluid patterns along the outer surface of the probe as well as in the interior channel. The simulations take account of complex flow phenomena such as laminar-turbulent transition and flow separation.


Author(s):  
Francisco de Melo Viríssimo ◽  
Paul A. Milewski

The problem of two layers of immiscible fluid, bordered above by an unbounded layer of passive fluid and below by a flat bed, is formulated and discussed. The resulting equations are given by a first-order, four-dimensional system of PDEs of mixed-type. The relevant physical parameters in the problem are presented and used to write the equations in a non-dimensional form. The conservation laws for the problem, which are known to be only six, are explicitly written and discussed in both non-Boussinesq and Boussinesq cases. Both dynamics and nonlinear stability of the Cauchy problem are discussed, with focus on the case where the upper unbounded passive layer has zero density, also called the free surface case. We prove that the stability of a solution depends only on two ‘baroclinic’ parameters (the shear and the difference of layer thickness, the former being the most important one) and give a precise criterion for the system to be well-posed. It is also numerically shown that the system is nonlinearly unstable, as hyperbolic initial data evolves into the elliptic region before the formation of shocks. We also discuss the use of simple waves as a tool to bound solutions and preventing a hyperbolic initial data to become elliptic and use this idea to give a mathematical proof for the nonlinear instability.


2002 ◽  
Vol 46 (03) ◽  
pp. 186-200 ◽  
Author(s):  
Pierre C. Sames ◽  
Delphine Marcouly ◽  
Thomas E. Schellin

To validate an existing finite volume computational method, featuring a novel scheme to capture the temporal evolution of the free surface, fluid motions in partially filled tanks were simulated. The purpose was to compare computational and experimental results for test cases where measurements were available. Investigations comprised sloshing in a rectangular tank with a baffle at 60% filling level and in a cylindrical tank at 50% filling level. The numerical study started with examining effects of systematic grid refinement and concluded with examining effects of three-dimensionality and effects of variation of excitation period and amplitude. Predicted time traces of pressures and forces compared favorably with measurements.


2000 ◽  
Vol 402 ◽  
pp. 109-136 ◽  
Author(s):  
AMY WARNCKE LANG ◽  
MORTEZA GHARIB

This experimental investigation into the nature of free-surface flows was to study the effects of surfactants on the wake of a surface-piercing cylinder. A better understanding of the process of vorticity generation and conversion at a free surface due to the absence or presence of surfactants has been gained. Surfactants, or surface contaminants, have the tendency to reduce the surface tension proportionally to the respective concentration at the free surface. Thus when surfactant concentration varies across a free surface, surface tension gradients occur and this results in shear stresses, thus altering the boundary condition at the free surface. A low Reynolds number wake behind a surface-piercing cylinder was chosen as the field of study, using digital particle image velocimetry (DPIV) to map the velocity and vorticity field for three orthogonal cross-sections of the flow. Reynolds numbers ranged from 350 to 460 and the Froude number was kept below 1.0. In addition, a new technique was used to simultaneously map the free surface deformation. Shadowgraph imaging of the free surface was also used to gain a better understanding of the flow. It was found that, depending on the surface condition, the connection of the shedding vortex filaments in the wake of the cylinder was greatly altered with the propensity for surface tension gradients to redirect the vorticity near the free surface to that of the surface-parallel component. This result has an impact on the understanding of turbulent flows in the vicinity of a free surface with varying surface conditions.


Author(s):  
G. D. Gkikas ◽  
F. van Walree

A computational method for the seakeeping behavior of a cruise ship at zero speed and under severe wind and oblique wave loads is presented. The proposed methodology is a time-domain panel method where the transient Green functions used for the estimation and implementation of the free surface effects on the vessel’s motions are estimated assuming constant low lateral speed, instead of the common practice zero speed influence functions. For the evaluation of the overall hydrodynamic forces, the so called “blended approach” is followed in the sense that the induced hydrodynamic pressures due to the scattering and radiation phenomena are calculated over the linearized position of the body, ignoring any displacements with respect to its mean position, while the hydrostatic and non-linear Froude-Krylov forces are considered at the actual body location and taking into account the free surface elevation at each time step. For the validation of the proposed methodology, heave and roll motions, the drift velocity as well as lateral accelerations of the vessel were investigated for two cases of severe beam seas combined with a constant strong wind load and the results were compared against experimental model tests. The model tests were performed to investigate the vessel’s behavior under extreme weather conditions. The low lateral speed Green functions were estimated for a speed similar to the one that the vessel was expected to drift, an estimation based on the model tests, as well as for the case where the input speed corresponded to the half of the expected speed. Good agreement was presented for both cases, showing that accurate and computationally efficient numerical simulations of the vessel’s motions under severe wind and wave excitations can be obtained by using low lateral speed transient Green functions.


Author(s):  
E.Y. Kenig ◽  
A.A. Ganguli ◽  
T. Atmakidis ◽  
P. Chasanis

2008 ◽  
Vol 141 (1-3) ◽  
pp. 204-221 ◽  
Author(s):  
Vinay R. Gopala ◽  
Berend G.M. van Wachem

Sign in / Sign up

Export Citation Format

Share Document