• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.99-108
• /
• 2010

Fractional Step Methods(FSM) are popular in simulation of unsteady incompressible flow. In this study, we demonstrate that FSM, combined with a Volume-Of-Fluid method, can be further applied to simulation of multiphase flow. The interface between the fluids is constructed by the effective least squares volume-of-fluid interface reconstruction algorithm and advected by the velocity using the operator split advection algorithm. To verify our numerical methodology, our results are compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. The present FSM sheds light on accurate simulation of turbulent multiphase flow which is found in many engineering applications.

• Journal of computational fluids engineering
• /
• v.20 no.2
• /
• pp.103-108
• /
• 2015

Fluid-body interaction analysis of floating body with six degree-of-freedom motion is presented. In this study, three-dimensional incompressible Navier-Stokes equations are employed as a governing equation. The numerical method is based on a finite-volume approach on a cartesian grid together with a fractional-step method. To represent the body motion, the immersed boundary method for direct forcing is employed. In order to simulate the coupled six degree-of-freedom motion, Euler's equations based on rigid body dynamics are utilized. To represent the complex body shape, level-set based algorithm is utilized. In order to describe the free surface motion, the volume of fluid method utilizing the tangent of hyperbola for interface capturing scheme is employed. This study showed three different continuums(air, water and body) are simultaneously simulated by newly developed code. To demonstrate the applicability of the current approach, two different problems(dam-breaking with stationary obstacle and water entry) are simulated and all results are validated.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.9-10
• /
• 2006

A numerical method for simulating tree surface flows including the surface tension is presented. Numerical scheme is based an a fractional-step method with a finite volume formulation and the interface between liquid and gas is tracked by Volume of Fluid (VOF) method. Piecewise Linear Interface Calculation (PLIC) method is used to reconstruct the interface and the surface tension is considered using a Continuum Surface Force (CSF) model. Several free surface flow phenomena were simulated to show its effectiveness to find such phenomena.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.910-915
• /
• 2003

A Navier-Stokes equation solver for incompressible viscous flows with free surface is developed and tested. This is based upon a fractional time step method and a non-staggered finite volume formulation for unstructured meshes. For time advancement scheme, Adams -Bashforth method for convective term and Crank-Nicolson method for diffusive term are applied. The interface between two fluids with different fluid properties is tracked with Piecewise Linear Interface Calculation(PLIC) Volume-of-Fluid(VOF) methods. Computational results are presented for some test problems: the broken dam, the sloshing in a rectangular tank, the filling of a cylindrical tank.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.12
• /
• pp.1555-1569
• /
• 2000

Numerical simulation of fluid flow with moving free surface has been carried out. For the free surface flow, a VOF(Volume of Fluid)-based algorithm utilizing a fixed grid system has been investigated. In order to reduce numerical smearing at the free surface represented on a fixed grid system, a new free surface tracking algorithm based on the donor-acceptor scheme has been presented. Novel features of the proposed algorithm are characterized as two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed algorithm can be easily implemented in any irregular non-uniform grid systems that are usual in finite element method (FEM). Moreover, the proposed algorithm can be extended and applied to the 3-D free surface flow problem without additional efforts. For computation of unsteady incompressible flow, a finite element approximation based on the explicit fractional step method has been adopted. In addition, the SUPG(streamline upwind/Petrov-Galerkin) method has been implemented to deal with convection dominated flows. Combination of the proposed free surface tracking scheme and explicit fractional step formulation resulted in an efficient solution algorithm. Validity of the present solution algorithm was demonstrated from its application to the broken dam and the solitary wave propagation problems.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.3
• /
• pp.248-257
• /
• 2011

In the present study, we proposed a new numerical wave tank model to analyze the vertical tension-leg circular floating bodies, using a 2-D Navier-Stokes solver. An IBM(Immersed Boundary Method) capable of handling interactions between waves and moving structures with complex geometry on a standard regular Cartesian grid system is coupled to the VOF(Volume of Fluid) method for tracking the free surface. Present numerical results for the motions of the floating body were compared with existing experimental data as well as numerical results based on FAVOR(Fractional Area Volume Obstacle Representation) algorithm. For detailed examinations of the present model, the additional hydraulic experiments for floating motions and free surface transformations were conducted. Further, the versatility of the proposed numerical model was verified via the numerical and physical experiments for the general rectangular floating bodies. Numerical results were compared with experiments and good agreement was archived.

• Journal of Biomedical Engineering Research
• /
• v.37 no.2
• /
• pp.84-89
• /
• 2016

In this study, we proposed ameliorated method for region of interest (ROI) study to improve its accuracy using partial volume effect (PVE). PVE which arose in volumetric images when more than one tissue type occur in a voxel, could be used to reduce an amount of gray matter and cerebrospinal fluid within ROI of diffusion tensor image (DTI). In order to define ROIs, individual b0 image was spatially aligned to the JHU DTI-based atlas using linear and non-linear registration (http://cmrm.med.jhmi.edu/). Fractional anisotropy (FA) and mean diffusivity (MD) maps were estimated by fitting diffusion tensor model to each image voxel, and their mean values were computed within each ROI with PVE threshold. Participants of this study consisted of 20 healthy controls, 27 Alzheimer's disease and 27 normal-pressure hydrocephalus patients. The result showed that the mean FA and MD of each ROI were increased and decreased respectively, but standard deviation was significantly decreased when PVE was applied. In conclusion, the proposed method suggested that PVE was indispensable to improve an accuracy of DTI ROI study.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.1-10
• /
• 1998

서로 섞이지 않는 두 비압축성 유체의 유동을 해석하기 위하여 VOF 방법에 기초한 수치 기법을 개발하였다. 유체간의 계면형상의 거동은 유동장내의 유체의 점유체적비의 변화에 의해 묘사되는데 이를 지배하는 이동방정식을 풀기 위한 새로운 대류항 계산법을 고안하였다. 대류항은 유체계면의 방향에 따라 풍상법과 역풍상법의 적절한 조합을 취하여 계산하는데 여기에 대각방향의 상류효과를 포함시켜 시간에 대한 2차 정확도를 갖도록 하였다. 또한 이 방법을 유량보정수송(FCT)법과 결합시켜 해의 단조성을 보장하였다. 몇 가지 단순 문제에 대한 시험 결과 이 기법이 수치오차에 의한 계면형상의 변형과 파손을 감소시킴을 확인하였다.

• Wind and Structures
• /
• v.38 no.4
• /
• pp.261-275
• /
• 2024

Before performing an experimental study on the downburst-generated wave, it is necessary to examine the scale effects and corresponding corrections or compensations. Analysis of similarity is conducted to conclude the non-dimensional force ratios that account for the dynamic similarity in the interaction of downburst with wave between the prototype and the scale model, along with the corresponding scale factors. The fractional volume of fluid (VOF) method in association with the impinging jet model is employed to explore the characteristics of the downburst-generated wave numerically, and the validity of the proposed scaling method is verified. The study shows that the location of the maximum radial wind velocity in a downburst-wave field is a little higher than that identified in a downburst over the land, which might be attributed to the presence of the wave which changes the roughness of the underlying surface of the downburst. The impinging airflow would generate a concavity in the free surface of the water around the stagnation point of the downburst, with a diameter of about two times the jet diameter (D_jet). The maximum wave height appears at the location of 1.5D_jet from the stagnation point. Reynolds number has an insignificant influence on the scale effects, in accordance with the numerical investigation of the 30 scale models with the Reynolds number varying from 3.85 × 10⁴ to 7.30 × 10⁹. The ratio of the inertial force of air to the gravitational force of water, which is denoted by G, is found to be the most significant factor that would affect the interaction of downburst with wave. For the correction or compensation of the scale effects, fitting curves for the measures of the downburst-wave field (e.g., wind profile, significant wave height), along with the corresponding equations, are presented as a function of the parameter G.

• Journal of Mechanical Science and Technology
• /
• v.16 no.5
• /
• pp.720-731
• /
• 2002

Flow with moving free surfaces is analyzed with an the Eulerian coordinate system. This study proposes a semi-implicit filling algorithm using VOF in which the PLIC (Piecewise Linear Interface Calculation) -type interface reconstruction method and the donor-acceptor-type front advancing scheme are adopted. Also, a new scheme using extrapolation of the stream function is proposed to find the velocity of the node that newly enters the computational domain. The effect of wall boundary conditions on the flow field and temperature field is examined by numerically solving a two-dimensional casting process.