• 한국전산유체공학회지
• /
• 제21권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.

• 한국전산유체공학회지
• /
• 제19권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.

• 대한수학회지
• /
• 제59권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.

• 대한기계학회논문집B
• /
• 제24권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.

• 대한조선학회논문집
• /
• 제51권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.

• 한국해안·해양공학회논문집
• /
• 제24권4호
• /
• pp.277-286
• /
• 2012

기체와 액체의 유동 및 고체의 변형을 동시에 고려할 수 있는 혼상류모델을 이용하여 파동장을 해석하는 경우, 서로 다른 비혼합의 유체 경계면의 시간변형을 고정도로 추적하는 것이 대단히 중요하다. 본 연구에서는 경계면의 추적에 있어서 VOF(Volume of Fluid)법으로 대표되는 경계면 형상의 재구축이 필요한 Geometrical-type의 경계추적법의 대신에 Algebraic-type의 경계추적법인 THINC(Tangent of Hyperbola for INterface-Capturing)법을 적용하였다. THINC법은 경계면에 대한 형상의 구축이 필요하지 않으므로 VOF법에 비해 비교적 간단한 알고리즘을 가지며, 기존의 Navier-Stokes solver에로 적용성이 용이한 장점을 갖는다. 본 연구에서는 THINC법의 기본적인 이류특성을 고찰하고, 혼상류수치모델인 TWOPM(one-field Model for immiscible TWO-Phase flows)과 결합한 수치모델을 파동장에 적용하여 비혼합 혼상류에서 경계면의 추적능을 검토하였다. 그 결과, 혼상류의 경계면 추적에 있어서 상대적으로 간단한 알고리즘의 THINC법이 기존의 VOF법과 유사한 정도를 갖는 해석법이라는 것을 확인할 수 있었고, 따라서 향후 기포의 연행을 동반하는 쇄파 및 쇄파력의 해석 등에 그의 적용성이 기대된다.

• 한국전산구조공학회논문집
• /
• 제28권1호
• /
• pp.71-78
• /
• 2015

LNG 선박에서 발생하는 슬로싱 충격하중은 다상유동 및 기체의 압축효과에 따라 CCS에서 발생하는 압력과 구조응답에 큰 영향을 미칠 수 있다. 본 연구에서는 슬로싱 운동 시 LNG의 유동에 의해 발생하는 슬로싱 충격을 시뮬레이션하기 위해서 다상유동을 적용한 수치해석 모델을 제시하였으며, 그 결과를 실험과 비교하여 타당성을 검토하였다. 또한 효율적인 구조 응답 계산을 위해 분사모델을 이용한 유체구조 연성해석 방법에 대해서 검토하고 멤브레인형 Mark III 화물창의 강도평가에 적용하여 LNG 화물창의 강도평가를 위한 가능성을 검토하였다.

• Structural Engineering and Mechanics
• /
• 제64권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.

• 한국전산유체공학회지
• /
• 제12권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.

• 한국유체기계학회 논문집
• /
• 제7권6호
• /
• pp.55-61
• /
• 2004

A numerical procedure for optimizing the shape of three-dimensional sedimentation tank is presented to maximize its sedimentation efficiency. The response surface based optimization is used as an optimization technique with Reynolds-averaged Navier-Stokes analysis for multi-phase flow. Standard $k-{\epsilon}$ model is used as a turbulence closure. Three design variables such as, tank height to center feed wall diameter ratio, blockage ratio of center feed wall and angle of distributor are chosen as design variables. Sedimentation efficiency is defined as an objective function. Full-factorial method is used to determine the training points as a means of design of experiment. Sensitivity of each design variable on the objective function has been evaluated. And, optimal values of the design variables have been obtained.