• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.21-25
• /
• 2007

A general purpose program NUFLEX has been extended for two-phase flows with topologically complex interface and cavitation flows with liquid-vapor phase change caused by large pressure drop. In analysis of two-phase flow, the phase interfaces are tracked by employing a LS(Level Set) method. Compared with the VOF(Volume-of-Fluid) method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. Also, the cavitation process is computed by including the effects of evaporation and condensation for bubble formation and collapse as well as turbulence in flows. The volume-faction and continuity equations are adapted for cavitation models with phase change. The LS and cavitation formulation are implemented into a general purpose program for 3-D flows and verified through several test problems.

• Journal of Ocean Engineering and Technology
• /
• v.4 no.1
• /
• pp.55-61
• /
• 1990

For the tracking of three dimensional free surface motions, a method referred to as the Volume of Fluid(VOF) algorithm is extended. In order to calculate the slope of three dimensional free surface which is the most important for the advection algorithm that decides the amount of fluid from cell to cell and for the application of free surface boundary condition, a simple method utilizing two dimensional slope informations is introduced. The extended algroithm is tested by demonstrating the simulation of a propagating sinusoidal wave through the channel whose width changes abruptly.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.1
• /
• pp.1-13
• /
• 2008

In the present work, wave transformation by vertical columns and its wave forces acting on them are discussed using a direct 3-D numerical model based on the VOF (Volume Of Fluid) method. The numerical results for wave transformations and wave forces are critically compared to an advanced experimental data, and provide the verification of the numerical model used in the present study. Overall model-data comparisons are good. After verification of the numerical model, it is used to simulate wave fields around dual vertical columns with arbitrary cross section, and the characteristics of nonlinear wave forces and wave transformations according to the variations of different cross section types of vertical columns, an interval of vertical columns and incident wave angle are discussed.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2005.10a
• /
• pp.269-272
• /
• 2005

For the three decades, the mold-filling of injection molding process was modeled as Hele-Shaw model. However, this model can not consider the 3D effect. In this paper, numerical simulations of three dimensional mold-filling during the filling phase were performed. The governing equations were discretized by segregated finite element method, which used equal order interpolation for pressure and velocity fields. The iterative linear equation solver (JCG, SOR) was employed for the solution of the momentum and pressure equations. Volume of Fluid (VOF) was employed for the melt front advancement. To check the validity of the numerical results, the results were compared with the experimental ones. The agreements between the experiment and the numerical results were found to be satisfactory.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.409-418
• /
• 2014

Simulation of flow past a complex marine structure requires a fine resolution in the vicinity of the structure, whereas a coarse resolution is enough far away from it. Therefore, a lot of grid cells may be wasted, when a simple Cartesian grid system is used for an Immersed Boundary Method (IBM). To alleviate this problems while maintaining the Cartesian frame work, we adopted an Adaptive Mesh Refinement (AMR) scheme where the grid system dynamically and locally refines as needed. In this study, We implemented a moving IBM and an AMR technique in our basic 3D incompressible Navier-Stokes solver. A Volume Of Fluid (VOF) method was used to effectively treat the free surface, and a recently developed Lagrangian Dynamic Subgrid-scale Model (LDSM) was incorporated in the code for accurate turbulence modeling. To capture vortex induced vibration accurately, the equation for the structure movement and the governing equations for fluid flow were solved at the same time implicitly. Also, We have developed an interface by using AutoLISP, which can properly distribute marker particles for IBM, compute the geometrical information of the object, and transfer it to the solver for the main simulation. To verify our numerical methodology, our results were compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. Using the verified code, we investigated the following cases. (1) simulating flow around a floating sphere. (2) simulating flow past a marine structure.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.219-235
• /
• 2014

With the development of the technology of underwater moving bodies, the need for developing the knowledge of surface effect interaction of free surface and underwater moving bodies is increased. Hence, the two-phase flow is a subject which is interesting for many researchers all around the world. In this paper, the non-linear free surface deformations which occur during the water-exit of a circular cylinder due to its buoyancy are solved using finite volume discretization based code, and using Volume of Fluid (VOF) scheme for solving two phase flow. Dynamic mesh model is used to simulate dynamic motion of the cylinder. In addition, the effect of cylinder mass in presence of an external force is studied. Moreover, the oblique exit and entry of a circular cylinder with two exit angles is simulated. At last, water-exit of a circular cylinder in six degrees of freedom is simulated in 3D using parallel processing. The simulation errors of present work (using VOF method) for maximum velocity and height of a circular cylinder are less than the corresponding errors of level set method reported by previous researchers. Oblique exit shows interesting results; formation of waves caused by exit of the cylinder, wave motion in horizontal direction and the air trapped between the waves are observable. In 3D simulation the visualization of water motion on the top surface of the cylinder and the free surface breaking on the front and back faces of the 3D cylinder at the exit phase are observed which cannot be seen in 2D simulation. Comparing the results, 3D simulation shows better agreement with experimental data, specially in the maximum height position of the cylinder.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.2
• /
• pp.149-159
• /
• 2016

In this study, the solitary wave-induced tsunami force acting on an onshore bridges in coastal area was numerically modelled by means of TWOPM-3D based on Navier-Stokes solver and VOF method which can track free surface effectively. The validity of numerical analysis was verified by comparing the experimental tsunami bore force acting on vertical wall and column structure. In particular, the characteristics of tsunami force with the changing tsunami intensity were surveyed through numerical experiments. The availability of 3-dimensional numerical analysis was reviewed through the comparison between the existing numerical results and design criteria for each drag force coefficient by applying Morison equation considering only drag force. As reasonable and high-precision estimation method of tsunami force, it was suggested to apply the estimation method taking drag and inertial force into consideration at the same time.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.5
• /
• pp.405-414
• /
• 2015

Fluid transport is a key issue in the development of microfluidic systems. Recently, Myong (2014) has proposed a new concept for droplet transport without external power sources and numerically validated the results for a hypothetical 2D, initially having a hemicylindrical droplet. In this paper, the movement of an actual water droplet, initially having a 3D hemispherical shape, on horizontal hydrophilic/hydrophobic surfaces is simulated using a commercial computational fluid dynamics (CFD) package, Fluent, with VOF (volume of fluid) method. The results are compared with the 2D analysis of Myong (2014), and the transport mechanism for the actual water droplet is examined based on the numerical results of the time evolution of the droplet shape, as well as the total kinetic, gravitational, surface free and pressure energies inside the droplet.

• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.9
• /
• pp.571-577
• /
• 2017

In this numerical study, the separation efficiency of three-phase separator in an oil-sand plant was studied with various inlet head configurations. The free water knockout (FWKO) vessel was designed with a flow rate of $15.89m^3/day$ (100 bbl/day) and the SOR(stream-to-oil ratio)=3.5 was derived using Stokes' law. For modularization, optimization of the design of the inlet head configuration was performed with parallel-connected dual FWKO vessels. The feed condition of bitumen emulsion was API=17, $T_{in}=150^{\circ}C$ and $P_{in}=50bar$. A mean residence time was determined the time when 95% of the oil and water in FWKO vessel was separated. The combination between the volume of fluid (VOF) and the discrete phase model (DPM) was used to simulate the phase separation phenomenon in a multi-phase separator. Furthermore, in order to calculate multi-phase flow the pseudo-transient method was adopted.