• 한국해양공학회지
• /
• 제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.

• 전기화학회지
• /
• 제14권4호
• /
• pp.231-238
• /
• 2011

수치해석을 이용하여 아연공기전지의 전기화학적 성능을 예측하였다. KOH 수계전해질 내부의 이동현상을 예측하기 위하여 Nernst-Planck식을 사용하였고, 전극 표면의 활성화손실을 모사하기 위해 아연극(음극)에는 Butler-Volmer식을, 공기극(양극)에는 Tafel식을 적용하였다. 정상상태해석을 통하여 아연/공기전지의 I-V곡선을 도출하였으며, 실험결과와의 I-V곡선 비교/분석을 통하여 수치 해석 모델의 타당성을 검증하였다. 전지반응 진행에 따른 전해질 내부의 이온 이동 및 분포 특성을 조사하기 위해 과도상태해석을 수행하였으며, 전극주변에서의 ${Zn(OH)_4}^{-2}$, $OH^-$, $K^+$ 이온들의 농도변화를 확인할 수 있었다. 또한, 다양한 전지전압조건 하에서 반응시간에 따라 아연극 표면에서의 ${ZnOH_4}^{2-}$의 농도 변화를 해석한 결과, 반응진행시 아연극 표면에서 ${ZnOH_4}^{2-}$의 농도가 최고성능을 나타내는 운전전압 0.63 V에서 약 1초 만에 포화농도에 도달하였으며, 일반적인 운전조건인 1.04 V에서는 약 13초 만에 포화농도에 접근하는 것으로 나타났다.

• 한국전산유체공학회지
• /
• 제15권2호
• /
• pp.7-13
• /
• 2010

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 5.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}10^6$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental datum in the low aspect ratio cavity (L/D = ~4.5). In the high aspect ratio cavity, however, there are other low dominant frequencies of the leading edge shear layer with the dominant frequencies of the feedback mechanism.

• 대한기계학회논문집B
• /
• 제37권3호
• /
• pp.229-235
• /
• 2013

유연한 혈관벽을 가진 경동맥 분지관을 흐르는 혈액의 유동을 해석하기 위하여 비정상상태, 비압축성, 뉴턴 유체를 가정한 3차원 유한요소해석을 수행하였다. 유체영역은 P2P1 유한요소를 사용하였으며, 격자의 움직임을 모사하기 위하여 arbitrary Lagrangian-Eulerian 기법을 적용하였다. Newmark 관계식을 이용하여 고체영역의 선형탄성 방정식의 변수들을 속도에 관한 방정식으로 간략화하였으며, 유체와 고체의 운동에 관하여 완전 결합된 공식을 얻었다. 맥동의 한 주기 동안에 혈관벽의 유연성이 유동장에 큰 영향을 미치며, 경동맥 분지각이 커짐에 따라 경동맥 공동에서 유동장의 정체영역이 더 넓게 분포한다는 연구결과를 얻었다.

• Wind and Structures
• /
• 제23권5호
• /
• pp.465-483
• /
• 2016

This study aims to enhance the understanding of the surface pressure distribution around rectangular bodies, by considering aspects such as the suction pressure at the leading edge on the top and side faces when the body aspect ratio and wind direction are changed. We carried out wind tunnel measurements and numerical simulations of flow around a series of rectangular bodies (a cube and two rectangular bodies) that were placed in a deep turbulent boundary layer. Based on a modern numerical platform, the Navier-Stokes equations with the typical two-equation model (i.e., the standard $k-{\varepsilon}$ model) were solved, and the results were compared with the wind tunnel measurement data. Regarding the turbulence model, the results of the $k-{\varepsilon}$ model are in overall agreement with the experimental results, including the existing data. However, because of the blockage effects in the computational domain, the pressure recovery region is underpredicted compared to the experimental data. In addition, the $k-{\varepsilon}$ model sometimes will fail to capture the exact flow features. The primary emphasis in this study is on the flow characteristics around rectangular bodies with various aspect ratios and approaching wind directions. The aspect ratio and wind direction influence the type of wake that is generated and ultimately the structural loading and pressure, and in particular, the structural excitation. The results show that the surface pressure variation is highly dependent upon the approaching wind direction, especially on the top and side faces of the cube. In addition, the transverse width has a substantial effect on the variations in surface pressure around the bodies, while the longitudinal length has less influence compared to the transverse width.

• 대한기계학회논문집B
• /
• 제24권2호
• /
• pp.204-213
• /
• 2000

This study deals with jet impingement, which is extensively used in the process industries to achieve intense heating, cooling or drying rates and also widely employed as a test flow for turbulent models due to its complex flow configuration, on a flat plate by numerical methods. In this calculation, the finite volume method was employed to solve the Navier-stokes equation based on the non-orthogonal coordinate with non-staggered variable arrangement. To get a better understanding for the fluid flow and heat transfer characteristics of the turbulent jet impingements, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model was adapted and compared with the experimental data and the result of standard $k-{\varepsilon}$ turbulent model. Numerical calculations were carried out with various flow rates, nozzle to plate distances. In the case of the axisymmetric jet impingement on a flat plate, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model showed better agreement with the experimental data than the standard $k-{\varepsilon}$ turbulent model in the prediction of the mean velocity profiles, the turbulent velocity profiles. the turbulent shear stress and the heat transfer rate. The highest heat transfer rate can be obtained when the impingement occurs within the potential core..

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

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought about by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 2.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}106$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental data in the low aspect ratio cavity (L/D = ~ 4.5). In the large aspect ratio cavity, however, there are other low dominant frequencies due to the leading edge shear layer with the dominant frequencies of the feedback mechanism. The characteristics of the acoustic wave propagation are analyzed using the Correlation of Pressure Distribution (CPD).

• 한국전산유체공학회지
• /
• 제21권2호
• /
• pp.112-119
• /
• 2016

NEXTSat-1 is the next-generation small-size artificial satellite system planed by the Satellite Technology Research Center(SatTReC) in Korea Advanced Institute of Science and Technology(KAIST). For the control of attitude and transition of the orbit, the system has adopted a RHM(Resisto-jet Head Module), which has a very simple geometry with a reasonable efficiency. An axisymmetric model is devised with two coil-resistance heaters using xenon(Xe) gas, and the minimum required specific impulse is 60 seconds under the thrust more than 30 milli-Newton. To design the module, seven basic parameters should be decided: the nozzle shape, the power distribution of heater, the pressure drop of filter, the diameter of nozzle throat, the slant length and the angle of nozzle, and the size of reservoir, etc. After quasi one-dimensional analysis, a theoretical value of specific impulse is calculated, and the optima of parameters are found out from the baseline with a series of multi-physical numerical simulations based on the compressible Navier-Stokes equations for gas and the heat conduction energy equation for solid. A commercial code, COMSOL Multiphysics is used for the computation with a FEM (finite element method) based numerical scheme. The final values of design parameters indicate 5.8% better performance than those of baseline design after the verification with all the tuned parameters. The present method should be effective to reduce the time cost of trial and error in the development of RHM, the thruster of NEXTSat-1.

• Journal of Advanced Marine Engineering and Technology
• /
• 제36권8호
• /
• pp.1091-1096
• /
• 2012

유한체적법을 기반으로 나비에 스톡스 방정식을 비구조격자로 풀어 실린더 주위의 공력특성을 규명하였다. 보텍스, 속도, 압력, 잔차, 항력계수 등의 데이터를 가지고 분석하였고 레이놀즈 수는 50, 100이다. 유동특성은 Re>50에서 주기적으로 진동하는 소용돌이를 후류에 형성하며 이 현상은 이동하는 실린더에서도 유사한 현상을 보여 주었다. 지면효과는 실린더 위쪽에서 형성된 소용돌이가 벽면에 근접할수록 실린더의 후방으로 길게 늘어나는 형상을 보이고, 실린더와 근접벽면 사이의 유속이 정체되어 실린더와 벽면 사이의 간격이 0.6 에서는 근접 평판과 실린더 사이의 유동이 거의 끊겨짐을 알 수 있었다. 본 수치계산의 검증을 위하여 항력계수를 타 연구결과와 비교하였다.

• Wind and Structures
• /
• 제23권6호
• /
• pp.505-521
• /
• 2016

The objective of this study is to investigate numerically the effect of building roof shaps on wind flow and pollutant dispersion in a street canyon with one row of trees of pore volume, $P_{vol}=96%$. A three-dimensional computational fluid dynamics (CFD) model is used to evaluate air flow and pollutant dispersion within an urban street canyon using Reynolds-averaged Navier-Stokes (RANS) equations and the Explicit Algebraic Reynolds Stress Models (EARSM) based on k-${\varepsilon}$ turbulence model to close the equation system. The numerical model is performed with ANSYS-CFX code. Vehicle emissions were simulated as double line sources along the street. The numerical model was validated by the wind tunnel experiment results. Having established this, the wind flow and pollutant dispersion in urban street canyons (with six roof shapes buildings) are simulated. The numerical simulation results agree reasonably with the wind tunnel data. The results obtained in this work, indicate that the flow in 3D domain is more complicated; this complexity is increased with the presence of trees and variability of the roof shapes. The results also indicated that the largest pollutant concentration level for two walls (leeward and windward wall) is observed with the upwind wedge-shaped roof. But the smallest pollutant concentration level is observed with the dome roof-shaped.