• 대한기계학회논문집B
• /
• 제39권7호
• /
• pp.557-566
• /
• 2015

본 논문에서는 $600{\times}600{\mu}m$ 사각 마이크로 채널에서의 기-액체 테일러(슬러그) 흐름의 특성을 실험을 통해 살펴보았다. 실험 유체로는 질소와 물을 사용하였으며, 액체 및 기체 겉보기 속도는 각각 0.01 ~ 3 m/s, 0.1 ~ 3 m/s 의 범위에서 테일러 흐름이 나타나는 구간에서만 데이터를 얻었다. 기포 길이, 액체 슬러그 길이, 기포 속도 그리고 기포 생성 주파수를 고속 카메라를 사용하여 이미지 분석을 통해 측정하였다. 제시된 측정값(기포 길이, 액체 슬러그 길이, 기포 속도)과 이전 문헌에서 제안된 경험적 모델의 비교결과 대부분 오차가 50% 이상으로 나타났다. 따라서 기포와 액체 슬러그 길이 그리고 기포 속도에 대한 개선된 모델을 제시하였고, ${\pm}20%$ 이내의 비교적 우수한 결과를 볼 수 있었다. 또한 기포 생성 주파수는 기포 길이, 액체슬러그 길이 그리고 기포 속도의 관계를 이용하여 ${\pm}20%$ 이내에서 예측가능함을 알 수 있었다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2008년도 추계학술대회B
• /
• pp.2558-2561
• /
• 2008

In the present study, a parallel Taylor-Galerkin/level set based two-phase flow code was developed using finite element discretization and domain decomposition method based on MPI (Message Passing Interface). The proposed method can be utilized for the analysis of a large scale free surface problem in a complex geometry due to the feature of FEM and domain decomposition method. Four-step fractional step method was used for the solution of the incompressible Navier-Stokes equations and Taylor-Galerkin method was adopted for the discretization of hyperbolic type redistancing and advection equations. A Parallel ILU(0) type preconditioner was chosen to accelerate the convergence of a conjugate gradient type iterative solvers. From the present parallel numerical experiments, it has been shown that the proposed method is applicable to the simulation of large scale free surface flows.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제16권4호
• /
• pp.225-234
• /
• 2012

We present two high-order potential flow models for the evolution of the interface in the Rayleigh-Taylor instability in two dimensions. One is the source-flow model and the other is the Layzer-type model which is based on an analytic potential. The late-time asymptotic solution of the source-flow model for arbitrary density jump is obtained. The asymptotic bubble curvature is found to be independent to the density jump of the fluids. We also give the time-evolution solutions of the two models by integrating equations numerically. We show that the two high-order models give more accurate solutions for the bubble evolution than their low-order models, but the solution of the source-flow model agrees much better with the numerical solution than the Layzer model.

• 한국가시화정보학회:학술대회논문집
• /
• 한국가시화정보학회 2002년도 추계학술대회 논문집
• /
• pp.85-88
• /
• 2002

In order to visually analyze the flow characteristics, gas was injected into the liquid bath through nozzle installed at the center of bottom of the bath. When gas was injected into the liquid bath, several flow patterns were observed bubble-liquid plumb, the spout flow that occurred at the free surface, liquid circulation flow by bubble's behavior, etc. Various bubbles, from small bubbles to Taylor bubbles, consisted of the bubble-liquid plumb. In the pure liquid region, the large and small several vortices were formed and irregularly circulated. These irregular repetition and circulation play a important role of mixing in the bath. The vortices were developed in the upper and the side wall regions and the movement of flow in the low region was very small. It is known as 'dead zone'.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.223-227
• /
• 2009

In the present study, a least square weighted residual method and Taylor-Galerkin method were formulated and tested for the discretization of the two hyperbolic type equations of level set method; advection and reinitialization equations. The two approaches were compared by solving a time reversed vortex flow and three-dimensional broken dam flow by employing a four-step splitting finite element method for the solution of the incompressible Navier-Stokes equations. From the numerical experiments, it was shown that the least square method is more accurate and conservative than Taylor-Galerkin method and both methods are approximately first order accurate when both advection and reinitialization phase are involved in the evolution of free surface.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.389-392
• /
• 2002

In many industrial applications, multiphase flow analysis is the norm rather than an exception as compared to more-conventional single-phase investigation. This paper describes the implementation of the multiphase flow simulation capability in the general purpose CFD software AVL FIRE/SWIFT. The governing equations are discretized based on a finite volume method (FVM) suitable fur very complex geometry, The pressure field is obtained using the SIMPLE algorithm. Depending on the characteristics of the multiphase flow to be examined, the user can choose either the two-fluid model or an explicit interface-tracking model based on the Volume-of-Fluid approach. For truly 'multi'-phase flow problems, it is also possible to apply a hybrid model where certain phases are explicitly tracked while the other phases are handled by the two fluid model. In order to demonstrate the capability of the method, applications to the Taylor bubble flow simulations are presented.

• 한국해양공학회지
• /
• 제35권1호
• /
• pp.38-49
• /
• 2021

Slug flow is the most common multi-phase flow encountered in oil and gas industry. In this study, the hydrodynamic features of flow in pipes investigated numerically using computational fluid dynamic (CFD) simulations for the effect of slug flow on the vertical and bent pipeline. The compressible Reynold averaged Navier-Stokes (RANS) equation was used as the governing equation, with the volume of fluid (VOF) method to capture the outline of the bubble in a pipeline. The simulations were tested for the grid and time step convergence, and validated with the experimental and theoretical results for the main hydrodynamic characteristics of the Taylor bubble, i.e., bubble shape, terminal velocity of bubble, and the liquid film velocity. The slug flow was simulated with various air and water injection velocities in the pipeline. The simulations revealed the effect of slug flow as the pressure occurring in the wall of the pipeline. The peak pressure and pressure oscillations were observed, and those magnitudes and trends were compared with the change in air and water injection velocities. The mechanism of the peak pressures was studied in relation with the change in bubble length, and the maximum peak pressures were investigated for the different positions and velocities of the air and water in the pipeline. The pressure oscillations were investigated in comparison with the bubble length in the pipe and the oscillation was provided with the application of damping. The pressures were compared with the case of a bent pipe, and a 1.5 times higher pressures was observed due to the compression of the bubbles at the corner of the bent. These findings can be used as a basic data for further studies and designs on pipeline systems with multi-phase flow.

• 한국지하수토양환경학회:학술대회논문집
• /
• 한국지하수토양환경학회 2001년도 총회 및 춘계학술발표회
• /
• pp.89-92
• /
• 2001

This study aims at investigating the relationship between dispersion coefficient ratio to molecular diffusion coefficient (D$_{l}$ /D$_{m}$) and Peclet number (Pe) for multi-solute system in non-Darcian flow regime. Existing understanding on solute dispersion is primarily derived from one-solute system in Darcian flow regime. We found that solute dispersion in rock fractures can be characterized by the mechanism of both macrodispersion and Taylor dispersion, even for non-Darcian f]ow domain. For the Darcian flow regime even different solutes lead to the same D$_{l}$ /D$_{m}$ at same Pe. However, as the flow becomes non-Darcian, solute with a higher molecular diffusion coefficient result in higher D$_{l}$ /D$_{m}$ at tile same Pe than that with a lower diffusion coefficient.cient.

• 대한기계학회논문집B
• /
• 제28권7호
• /
• pp.789-799
• /
• 2004

An experimental study is conducted to investigate the effects of duct corrugation angle on heat/mass transfer characteristics in wavy ducts of a primary surface heat exchanger application. Local heat/mass transfer coefficients on the wavy duct sidewalls are determined by using a naphthalene sublimation technique. The corrugation angles(${\alpha}$) of the wavy ducts are 145$^{\circ}$, 130$^{\circ}$, 115$^{\circ}$ and 100$^{\circ}$. And the Reynolds numbers based on the duct hydraulic diameter vary from 300 to 3,000. The results show that at the low Re(Re $\leq$1000), the secondary vortices called Taylor-Gortler vortices perpendicular to the main flow direction are generated due to effect of duct curvature. By these secondary vortices, high heat/mass transfer regions are formed on both pressure-side and suction-side walls. At the high Re(Re $\geq$ 1000), these secondary flows are vanished with helping flow transition to turbulent flow and the regions which show high heat/mass coefficients by flow reattachment are formed on suction side. As corrugation angle decreases, the local peak Sh induced by Taylor-Gortler vortices increase at Re $\leq$1000. At high Re(Re $\geq$ 1000), by the existence of different kind of secondary flows called Dean vortices, non-uniform Sh distribution appears along spanwise direction at the narrow corrugation angle (${\alpha}$=100$^{\circ}$). Average Sh also increase by the enhanced effect of secondary vortices and flow reattachment. More pumping power (pressure loss) is required with the smaller corrugation angle due to the enhancement of flow instability.

• 대한기계학회논문집
• /
• 제19권7호
• /
• pp.1731-1740
• /
• 1995

Experimental investigations have been made to study the double-diffusive nature of convection of an initially stratified salt-water solution due to heating from below in a rotating cylindrical cavity. The objective is to examine the flow phenomena and the heat transfer characteristics according to the changes in temperature gradient, concentration gradient and rotating velocity of cavity. Thermal and solutal boundary conditions at side wall are adiabatic and impermeable, respectively. The top and bottom plate are maintained each at constant temperature and concentration. The cavity is put into a state of solid body rotation. Like the stationary case, the types of initially-formed flow pattern are classified into three regimes depending on the effective Rayleigh number and Taylor number; stagnant flow regime, single mixed-layer flow regime and successively formed multi-mixed layer flow regime. At the same effective Rayleigh number, the number of initially-formed mixed layer and its growth rate decrease as the effect of rotation increases. The temperature and concentration profiles are both uniform in each layer due to convective mixing in the layered-flow regime, but look both liner in stagnant flow regime and single mixed-layer flow regime. At the interface between adjacent layers, the temperature changes smoothly but the concentration changes rapidly.