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

A high resolution numerical method aimed at solving gas-liquid two-phase flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

• 한국전산유체공학회지
• /
• 제12권2호
• /
• pp.53-62
• /
• 2007

A high resolution numerical method aimed at solving cavitating flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media at isothermal condition and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

• 한국해양공학회지
• /
• 제20권4호
• /
• pp.31-36
• /
• 2006

Numerical schemes for spacing and time are tested to accurately simulate the wave propagation. The tested numerical schemesinclude 2nd-order central differencing, l-order upwind scheme, 2nd-order Leith scheme, 3rd-order MUSCLE, QUICK and QUICKEST schemes in spacing and the Euler and 4th-order Runge-Kutta(R-K) schemes in time. It is seen that more accurate results are expected when the higher-order schemes, especially the schemes combined with a TVD control limiter, are used for solving the wave equation. The 3rd-order upwind scheme with limiter and the 4th-order R-K scheme in tim￡ are finally applied to the wave-making simulation in a digital wave tank.

• 천문학회지
• /
• 제29권2호
• /
• pp.195-205
• /
• 1996

The evolution of the stellar debris after tidal disruption due to the super massive black hole's tidal force is difficult to solve numerically because of the large dynamical range of the problem. We developed an SPH (Smoothed Particle Hydrodynamics) - TVD (Total Variation Diminishing) hybrid code in which the SPH is used to cover a widely spread debris and the TVD is used to compute the stream collision more accurately. While the code in the present form is not sufficient to obtain desired resoultion, it could provide a useful tool in studying the aftermath of the stellar disruption by a massive black hole.

• 한국항공우주학회지
• /
• 제31권2호
• /
• pp.10-16
• /
• 2003

초음속 유동장에서 비정상 직열배치공동 유동 연구를 Navier-Stokes 방정식과 k-$\varepsilon$난류모델을 이용하여 수치 모사하였다. 공동주위의 비정상유동은 자유 전단층과 공동 내부 유동의 상호 작용에 의한 주기성으로 특정지어진다. 수치적 방법은 van Leer의 유량한계계수를 이용한 유량벡터분리법에 기반을 둔 TVD방법을 사용하였다. 먼저, 단일 공동에 대한 압력비 그래프를 통한 주진동 주파수를 비교하였다. 직열배치공동유동에서는 전후방 공동을 몇가지 세장비로 조합하여, 전방 공동의 영향에 의한 후방 공동 유동장의 주진동 주파수 특성등을 분석하였다. 전방공동의 세장비가 작은 경우에는 강한 전방공동의 전단층이 후방공동의 전단츨과 합하여 하나의 전단층 특성을 줌으로 전후방 공동의 주진동특성이 비슷하게 나타나고, 전방공동의 세장비가 큰경우에는 전방공동의 전단층이 확산되면서 내부에 후방공동의 전단층이 발생함으로 주진동수의 특성이 다르게 나타나는 것을 볼 수 있었다.

• 한국항공우주학회지
• /
• 제32권9호
• /
• pp.12-26
• /
• 2004

수치점성을 효과적으로 줄이기 위한 새로운 AUSM-type 수치기법을 개발하였다. TVD 제한자 분석을 통하여 제어면에서의 불성치를 보다 정확하게 예측할 수 있는 기준을 찾아내고 이를 이용하여 제어면의 물성치를 적절히 정의하는 M-AUSMPW+를 개발하였다. M-AUSMPW+의 장점은 다차원 유동을 해석하는데 있어 분명히 나타난다. M-AUSMPW+ 는 유동과 격자계가 일치하지 않는 격자계에서 수치점성을 효과적으로 제거할 수 있기 때문에 다차원 유동을 계산하는데 있어 기존의 어떠한 수치기법보다 정확하고 효과적인 계산이 가능하다. 접촉불연속면, 와류 유동, 충격파 경계층 상호작용, 정성 충격파관 문제를 통해 이를 확인하였다.

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

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

• 천문학회보
• /
• 제21권
• /
• pp.6.1-6.1
• /
• 1996

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

The unsteady supersonic flow over two- and three-Dimensional cavities has been analyzed by the integration of unsteady Reynolds-Averaged Navier-Stokes(RANS) with the k - w turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in cavities. Numerical method is upwind TVD scheme based on the flux vector split with the Van Leer limiters, and time accuracy is used explicit 4th stage Runge-Kutta scheme. Cavity flows are Comparison of two- and three-dimensional. The cavity has a L/D ratio of 3 for two-dimensional case. and same L/D and W/D ratio is 1 for three-dimensional case. The Mach and Reynolds numbers are held constant at 1.5 and 450000 respectively. For the three-dimensional case, the flow field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follow Rossiter's formula. On the other hand, the self-sustained oscillating flow transitions to a 'wake mode' for the two-dimensional simulation, with more violent fluctuations inside the cavity.

• 한국전산유체공학회지
• /
• 제10권4호통권31호
• /
• pp.51-58
• /
• 2005

The unsteady supersonic flow over two- and three-dimensional cavities has been analyzed by the integration of unsteady Reynolds-Averaged Navier-Stokes(RANS) with the k-$\omega$ turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in the cavity. An explicit 4th order Runge-Kutta scheme and an upwind TVD scheme based on the flux vector split with the van Leer limiters are used for time and space discritizations, respectively. The cavity has a L/D ratio of 3 for two-dimensional case, and same L/D and W/D ratio of I for three-dimensional case. The Mach and Reynolds numbers are 1.5 and 450000 respectively. In the three-dimensional flow, the field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follows Rossiter's formula. In the two-dimensional simulation, the self-sustained oscillating flow has more violent fluctuation inside the cavity. The primary fluctuating frequencies of two- and three- dimensional flow agree very well with the 2nd mode of Rossiter's frequency. In the three-dimensional flow, the 1st mode of frequency could be seen.