• Journal of Ocean Engineering and Technology
• /
• v.29 no.6
• /
• pp.427-435
• /
• 2015

This paper provides the numerical results of a simulation of the flow around a propeller working beneath the free surface. A finite volume method is used to solve the unsteady Reynolds averaged Navier-Stokes (URANS) equations, where the wave-making problem is solved using a volume-of-fluid (VOF) method. The numerical analysis focuses on the propeller wake structure affected by the free surface, where we consider another free surface boundary condition that treats the free surface as a rigid wall surface. The propeller wake under the effect of these two free surface conditions shows a reduction in the magnitude of the longitudinal and vertical flow velocities, and its vortical structures strongly interact with the free surface. The thrust and torque coefficient under the free surface effect decrease about 3.7% and 3.1%, respectively. Finally, the present numerical results show a reasonable agreement with the available experimental data.

• Journal of computational fluids engineering
• /
• v.21 no.1
• /
• pp.110-116
• /
• 2016

The present study numerically investigated the deformation of the interface of two-phase fluids flow in a blast furnace. To simulate three-dimensional(3D) incompressible viscous two-phase flow in the furnace filled with the air and molten iron, the volume of fluid(VOF) method based on the finite volume method has been utilized. In addition, the porous medium with the porosity has been considered as the bed of the particles such as cokes and char etc. For the comparison, the single phase flow and the two-phase flow without the porosity have been simulated. The two-phase flow without porosity condition revealed the smooth parabolic profile of the free surface near the outlet. However, the free surface under the porosity condition formed the viscous finger when the free surface was close to the outlet. This viscous finger accelerated the velocity of the free surface falling and the outflow velocity of the fluids near the outlet.

• Nuclear Engineering and Technology
• /
• v.48 no.2
• /
• pp.360-367
• /
• 2016

A windowless spallation target can provide a neutron source and maintain neutron chain reaction for a subcritical reactor, and is a key component of China's nuclear waste transmutation of coupling accelerator and subcritical reactor. The main issue of the windowless target design is to form a stable and controllable free surface that can ensure that energy spectrum distribution is acquired for the neutron physical design when the high energy proton beam beats the lead-bismuth eutectic in the spallation target area. In this study, morphology and flow characteristics of the free surface of the windowless target were analyzed through the volume of fluid model using computational fluid dynamics simulation, and the results show that the outlet cross section size of the target is the key to form a stable and controllable free surface, as well as the outlet with an arc transition. The optimization parameter of the target design, in which the radius of outlet cross section is $60{\pm}1mm$, is verified to form a stable and controllable free surface and to reduce the formation of air bubbles. This work can function as a reference for carrying out engineering design of windowless target and for verification experiments.

• The KSFM Journal of Fluid Machinery
• /
• v.15 no.5
• /
• pp.60-66
• /
• 2012

One of commonly physical phenomena encountered in pump sump systems in which its significant influence to the hydraulic performance of pump system plays an important role in the field of fluid engineering, is the appearance of free surface and submerged vortices. In this paper, a study of the vortices behavior and their formative mechanism of asymmetry is considered in this paper by using numerical approach. The Reynolds-Averaged Navier-Stokes (RANS) equations and k-omega Shear Stress Transport turbulence model used to describe the properties of turbulent flows, in company with VOF multiphase model, are implemented by Fluent code with multi-block structured grid system. In the numerical simulation, the calculated elevation of air-water interface and vortex core contours are used to classify visually surface vortices as well as submerged vortices. It is shown that the free surface vortex is identified by the concavity of liquid region from the free surface and swirling flow at that own plane. To investigate the distinctive behavior of these vortices corresponding to each given flow rate at the same water level, some numerical testing of them are considered here in such a manner that the flow pattern of surface vortex are obtained similarly to the obtained results from experiment. Furthermore, the influence due to the change of grid refinement and the variation of depth of the concavity are also considered in this paper. From that, these influential factors will be implemented to design a good pump sump with higher performance in the future.

• Journal of the korean society of oceanography
• /
• v.31 no.4
• /
• pp.173-182
• /
• 1996

A median-density fluid was injected into the upper layer of a two-layered fluid in a rotating cylindrical container. Several sets of the top and bottom boundary configurations were employed and the flow pattern of each layer including the injected fluid was observed to determine the factors that affect the path of the injected intermediate fluid. The axisymmetric path of the intermediate fluid when the upper layer had a free surface, changed into the asymmetric path with bulged-shape radial spreading whenever either the upper layer or the lower layer had ${\beta}$-effect. The internal Fronds number that controls the shape of the interface turned out to be the most important parameter that determines the radial spreading in terms of location and strength. When the upper and lower layer had the ${\beta}$-effect, convective overturning produced anticyclonic vortices at the frontal edge of the intermediate fluid, and that could enhance the vertical mixing of different density fluids. The intermediate fluid did not produce any topographic effect on the upper-layer motion during its spreading over the interface, since its thickness was very small. However, its anticyclonic motion within the bulged-shape produced a cyclonic motion in the lower layer just beneath the bulge.

• Coupled systems mechanics
• /
• v.2 no.4
• /
• pp.389-410
• /
• 2013

Three-dimensional Lagrangian fluid finite element is applied to seismic response analysis of an oil storage tank with a floating roof. The fluid element utilized in the present analysis is formulated based on the displacement finite element method considering only volumetric elasticity and its element stiffness matrix is derived by using one-point integration method in order to avoid volumetric locking. The method usually adds a rotational penalty stiffness to satisfy the irrotational condition for fluid motion and modifies element mass matrices through the projected mass method to suppress spurious hourglass-mode appeared in compensation for one-point integration. In the fluid element utilized in the present paper, a small hourglass stiffness is employed. The fluid and structure domains for the objective oil storage tank are modeled by eight-node solid elements and four-node shell elements, respectively, and the transient response of the floating roof structure or the free surface are evaluated by implicit direct time integration method. The results of seismic response analyses are compared with those by other method and the validation of the present analysis using three-dimensional Lagrangian fluid finite elements is shown.

• Special Issue of the Society of Naval Architects of Korea
• /
• 2015.09a
• /
• pp.56-62
• /
• 2015

In ship structure, many parts are in contact with inner or outer fluid as stern, ballast and oil tanks. Fatigue damages are sometimes observed in these tanks which seem to be caused by resonance with exciting force of engine and propeller. Vibration characteristics of these tanks in contact with fluid are significantly affected by fluid coupling effect. Therefore it is important to exactly predict vibration characteristics of tank structure. In order to estimate the vibration characteristics, the fluid-structure interaction(FSI) problem should be solved precisely. But it is difficult to estimate exactly the magnitude of the fluid coupling effect because it has some problems such as a fluid-structure interaction, influence by the free surface, vibration modes of structural panels and depth of water. In this paper, with fluid coupling effect, the effect of structural constraint between panels on the vibration characteristics are investigated numerically and discussed.

• Ocean Systems Engineering
• /
• v.13 no.1
• /
• pp.57-77
• /
• 2023

In the present study, the total hydrodynamic pressure exerted by the fluid on walls of rectangular tanks due to horizontal excitations of different frequencies, is investigated by pressure based finite element method. Fluid within the tanks is invisid, compressible and its motion is considered to be irrotational and it is simulated by two dimensional eight-node isoparametric. The walls of the tanks are assumed to be rigid. The total hydrodynamic pressure increases with the increase of exciting frequency and has maximum value when the exciting frequency is equal to the fundamental frequency. However, the hydrodynamic pressure has decreasing trend for the frequency greater than the fundamental frequency. Hydrodynamic pressure at the free surface is independent to the height of fluid. However, the pressure at base and mid height of vertical wall depends on height of fluid. At these two locations, the hydrodynamic pressure decreases with the increase of fluid depth. The depth of undisturbed fluid near the base increases with the increase of depth of fluid when it is excited with fundamental frequency of fluid. The sloshing of fluid with in the tank increases with the increase of exciting frequency and has maximum value when the exciting frequency is equal to the fundamental frequency of liquid. However, this vertical displacement is quite less when the exciting frequency is greater than the fundamental frequency.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.9 no.2 s.42
• /
• pp.71-80
• /
• 2005

The solution to a liquid sloshing problem is challenge to the field of engineering. This is not only because the dynamic boundary condition at the free surface is nonlinear, but also because the position of the free surface varies with time in a manner not known a priori. Therefore, this nonlinear phenomenon, which is characterized by the oscillation of the unrestrained free surface of the fluid, is a difficult mathematical problem to solve numerically and analytically. In this study, three-dimensional boundary element method(BEM), which is based on the so-called an arbitrary Lagrangian-Eulerian(ALE) approach for the fluid flow problems with a free surface, was formulated to solve the behavior of the nonlinear free surface motion. An ALE-BEM has the advantage to track the free surface along any prescribed paths by using only one displacement variable, even for a three-dimensional problem. Also, some numerical examples were presented to demonstrate the validity and the applicability of the developed procedure.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2000.04a
• /
• pp.427-434
• /
• 2000

This paper presents a computer program for a 2-D fluid-structure-soil interaction analysis. With this computer program the fluid can be modeled by a spurious free 4-node displacement-based fluid element which uses rotational penalty and mass projection technique in conjunction with the one point reduced integration scheme to remove the spurious zero energy modes. The structure and near field soil are discretized by the standard finite elements while the unbounded far field soil are discretized by the standard finite elements while the unbounded far field soil is represented by the frequency dependent dynamic infinite elements. Sine this method models directly the fluid-structure-soil system it can be applied to the dynamci analysis of 2-D liquid storage structure with complex geometry. For the purpose of verification dynamic analyses for tanks on a rigid foundation and on compliant embankment are carried out. Comparison of the present results with those by ANSYS program shows good agreement.