• Journal of computational fluids engineering
• /
• v.20 no.4
• /
• pp.93-101
• /
• 2015

This paper provides numerical results of the simulation for the flow around the hull and the propeller of KCS model ship advancing in shallow water conditions. A finite volume method is used to solve the unsteady Reynolds averaged Navier-Stokes(RANS) equations, where the wave-making problem is solved by using a volume-of-fluid(VOF) method. The wave formed near the hull surface in shallow water conditions shows a deep trough dominant pattern that causes the loss of buoyancy followed by hull squat. The flow past the hull increases as the depth of water decreases. However, the axial flow velocity around the stern shows a reduction in magnitude by the effect of shallow water accompanied by the hull-propeller interaction. As a results, the thrust and torque coefficient increase about 8.3% and 6.2%, respectively for a depth of h/T=3.0 corresponding to a depth Froude number of $F_h=0.693$. The resistance coefficient increases about 11.6% at this Froude number condition.

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

• Journal of computational fluids engineering
• /
• v.19 no.4
• /
• pp.52-60
• /
• 2014

Recently, coastal structures have been built to protect coastal areas. However, if a tidal wave caused by an earthquake hits the coast, it would cause catastrophic damages. It is important to analyze the basics and the characteristics of a tsunami to reduce damages caused by natural disasters. In this study, a tsunami passing over different topographical changes is simulated with VOF method based on FVM(Finite Volume Method). The reduction of both scale and velocity is accomplished by similarity analysis, and an initial energy is generated by increasing the water level as needed to create a tsunami as if it is caused by a crustal movement. It is found that the present method is appropriate to simulate the tsunami with its mechanism.

• Journal of the Korean Mathematical Society
• /
• v.59 no.3
• /
• pp.519-548
• /
• 2022

In this paper, a stabilized-penalized collocated finite volume (SPCFV) scheme is developed and studied for the stationary generalized Navier-Stokes equations with mixed Dirichlet-traction boundary conditions modelling an incompressible biological fluid flow. This method is based on the lowest order approximation (piecewise constants) for both velocity and pressure unknowns. The stabilization-penalization is performed by adding discrete pressure terms to the approximate formulation. These simultaneously involve discrete jump pressures through the interior volume-boundaries and discrete pressures of volumes on the domain boundary. Stability, existence and uniqueness of discrete solutions are established. Moreover, a convergence analysis of the nonlinear solver is also provided. Numerical results from model tests are performed to demonstrate the stability, optimal convergence in the usual L² and discrete H¹ norms as well as robustness of the proposed scheme with respect to the choice of the given traction vector.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.9
• /
• pp.1166-1174
• /
• 2000

This work presents a two-dimensional quasi-steady state model to study the fluid flow and heat transfer in high-power density welding process of thin AISI-304 stainless steel plates. The enthalpy method and the finite volume method were used for a numerical analysis of the mushy region phase change as well as the heat flow at the weld pool and the heat-affected zone. The results show that the mushy region distributed around the weld pool becomes wider downstream and the surface heat losses by convection and radiation can be significant factors in welding process especially when a welding speed is relatively low.

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

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.24 no.4
• /
• pp.277-286
• /
• 2012

In the numerical simulation of wave fields using a multi-phase flow model that considers simultaneous flows of materials with different states such as gas, liquid and solid, there is need of an accurate representation of the interface separating the fluids. We adopted an algebraic interface capturing method called tangent of hyperbola for interface-capturing(THINC) method for the capture of the free-surface in computations of multi-phase flow simulations instead of geometrical-type methods such a volume of fluid(VOF) method. The THINC method uses a hyperbolic tangent functions to represent the surface, and compute the numerical flux for the fluid fraction functions. One of the remarkable advantages of THINC method is its easy applicability to incorporate various numerical codes based on Navier-Stokes solver because it does not require the extra geometric reconstruction needed in most of VOF-type methods. Several tests were carried out in order to investigate the advection of interfaces and to verify the applicability of the THINC method to wave fields based on the one-field model for immiscible two-phase flows (TWOPM). The numerical results revealed that the THINC method is able to track the interface between air and water separating the fluids although its algorithm is fairly simple.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.1
• /
• pp.71-78
• /
• 2015

The reliable sloshing assessment methods for LNG CCS(cargo containment system) are important to satisfy the structural strength of the systems. Multiphase fluid flow of LNG and Gas Compressibility may have a large effect on excited pressures and structural response. Impinging jet model has been introduced to simulate the impact of the LNG sloshing and analyze structural response of LNG CCS as a practical FSI(fluid structure interaction) method. The practical method based on fluid structure interaction analysis is employed in order to evaluate the structural strength in actual scale for Mark III CCS. The numerical model is based on an Euler model that employs the CVFEM(control volume based finite element method). It includes the particle motion of gas to simulate not only the interphase interaction between LNG liquid and gas and the impact load on the LNG insulation box. The analysis results by proposed method are evaluated and discussed for an effectiveness of FSI analysis method.

• Structural Engineering and Mechanics
• /
• v.64 no.4
• /
• pp.515-526
• /
• 2017

Dynamic response of functionally graded Carbon nanotubes (FG-CNT) reinforced pipes conveying viscous fluid under accelerated moving load is presented. The mixture rule is used for obtaining the material properties of nano-composite pipe. The radial force induced by viscous fluid is calculated by Navier-Stokes equation. The material properties of pipe are considered temperature-dependent. The structure is simulated by Reddy higher-order shear deformation shell theory and the corresponding motion equations are derived by Hamilton's principal. Differential quadrature (DQ) method and the Integral Quadrature (IQ) are applied for analogizing the motion equations and then the Newmark time integration scheme is used for obtaining the dynamic response of structure. The effects of different parameters such as boundary conditions, geometrical parameters, velocity and acceleration of moving load, CNT volume percent and distribution type are shown on the dynamic response of pipe. Results indicate that increasing CNTs leads to decrease in transient deflection of structure. In accelerated motion of the moving load, the maximum displacement is occurred later with respect to decelerated motion of moving load.

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