• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 추계 학술대회논문집
• /
• pp.1-9
• /
• 2001

In this paper, some numerical methods recently developed for gas-liquid two-phase flows are reviewed. And then, a preconditioning method to solve cavitating flow by the author is introduced. This method employs a finite-difference Runge-Kutta method combined with MUSCL TVD scheme, and a homogeneous equilibrium cavitation model. So that it permits to treat simply the whole gas-liquid two-phase flow field including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Finally, numerical results such as detailed observations of the unsteady cavity flows, a sheet cavitation break-off phenomena and some data related to performance characteristics of hydrofoils are shown.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2000년도 춘계 학술대회논문집
• /
• pp.84-90
• /
• 2000

Ground effect of three-dimensional wing is studied. LU-factored Implicit upwind TVD scheme and Baldwin-Lomax turbulence model are used for this calculation. To investigate ground effect, NACA 4415 wing at M=0.5 calculated. Two different angles of attack and three cases of flight height are calculated. As increasing angle of attack, the ground effect becomes strong. In case of NACA 4415 wing in ground effect, strength of wing tip vortex becomes stronger than that of free flight.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.181-182
• /
• 2003

A preconditioned numerical method for gas-liquid to-phase flow is applied to solve cavitating flow. The present method employs a density based finite-difference method of dual time-stepping integration procedure and Roe's flux difference splitting approximation with MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. The method permits simple treatment of the whole gas-liquid two-phase flow field including wave propagation, large density changes and incompressible flow characteristics at low Mach number. By this method, two-dimensional internal flows through a venturi tuve and decelerating cascades are computed and discussed.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.215-220
• /
• 2008

A numerical method for gas-liquid two-phase flow is applied to solve shock-bubble interaction problems. The present method employs a finite-difference Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. By this method, a Riemann problem for shock tube was computed for validation. Then, shock-bubble interaction problems between cylindrical bubbles located in the liquid and incident liquid shock wave are computed.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년 추계학술대회논문집
• /
• pp.215-220
• /
• 2008

A numerical method for gas-liquid two-phase flow is applied to solve shock-bubble interaction problems. The present method employs a finite-difference Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. By this method, a Riemann problem for shock tube was computed for validation. Then, shock-bubble interaction problems between cylindrical bubbles located in the liquid and incident liquid shock wave are computed.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.49-54
• /
• 1998

3-dimensional Euler solver is parallelized. The spatial discretization method is the 2nd order TVD scheme and DADI method with multigrid is used as a time integration. In order to parallelize this solver, the domain decomposition method with overlapped grid and message passing techniques are used. The informations on the each inter-processor bound-aries are communicated with MPI library. Finally, the parallel performance repsented by calculating the ONERA M6 wing at transonic flow condition using CRAY T3E and C90.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.29-34
• /
• 1998

An efficient 3-dimensional compressible solver is developed using the second-order upwind TVD scheme and the multigrid diagonalized ADI method. The multigrid method is improved so that the present DADI algorithm obtains better convergence rates. Results are computed on Cray C90 computer for transonic unsaperated flows past ONERA-M6 wing to demonstrate the accuracy and efficiency. The results show good agreement with experimetal data. A reduction of four orders of residual for 3-dimensional transonic flow is obtained about 99 seconds.

• 천문학회보
• /
• 제36권2호
• /
• pp.56.2-56.2
• /
• 2011

Building a relativistic magnetohydrodynamic (RMHD) codes based on upwind schemes is a challenging project, because the characteristic wave structures for RMHDs has not yet been analytically given. We obtained an analytic expression of eigenvalues and eigenvectors of the flux Jacobian matrix of RMHDs for one-dimensional, isothermal flows with two velocity and magnetic field components (that is, x and y components only), which can be used to build numerical codes. The degeneracies were taken into account. Here, we present preliminary test results with an RMHD code based on the total variation diminishing (TVD) scheme.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년도 제30회 춘계학술대회논문집
• /
• pp.285-288
• /
• 2008

본 연구는 연료(수소)층과 산화제(공기)층의 사이에 불활성기체(질소)또는 연료(수소)를 평행분사하는 수치해석을 다루고 있다. 수치해석을 위해서 완전 보존되는 비정상 2차 시간정확도법과 2차 TVD방법이 유한 체적법과 사용되었다. 결과는 3가지 종류로 구성되어있다. 첫째는 연료와 산화제의 단일 혼합층이고, 둘째는 연료와 산화제의 사이에 불활성기체를 분사하는 방식이며, 세 번째는 연료와 산화제의 사이에 연료를 분사하는 방식이다. 전체 유동층의 수직두께는 4cm이며 삽입된 중간층의 두께는 1,2,4mm의 세가지 경우에 대하여 계산하였다.

• 천문학회보
• /
• 제44권2호
• /
• pp.47.2-47.2
• /
• 2019

Observations indicate that turbulence in molecular clouds of the interstellar medium (ISM) is highly supersonic (M >> 1) and strongly magnetized (β ≈ 0.1), while in the intracluster medium (ICM) it is subsonic (M <~1) and weakly magnetized (β ≈ 100). Here, M is the turbulent Mach number and β is the ratio of the gas to magnetic pressures. Although magnetohydrodynamic (MHD) turbulence in such environments has been previously studied through numerical simulations, some of its properties as well as its consequences are not yet fully described. In this talk, we report a study of MHD turbulence in molecular clouds and the ICM using a newly developed code based the high-order accurate, WENO (Weighted Essentially Non-Oscillatory) scheme. The simulation results using the WENO code are generally in agreement with those presented in the previous studies with, for instance, a TVD code (Porter et al. 2015 &, Park & Ryu 2019), but reveal more detailed structures on small scales. We here present and compare the properties of simulated turbulences with WENO and TVD codes, such as the spatial distribution of density, the density probability distribution functions, and the power spectra of kinetic and magnetic energies. We also describe the populations of MHD shocks and the energy dissipation at the shocks. Finally, we discuss the implications of this study on star formation processes in the ISM and shock dissipation in the ICM.