• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• pp.723-728
• /
• 2008

We performed numerical analysis of base drag phenomenon, when a projectile with backward step flies into atmosphere at supersonic speed. We compared with other researchers. From our previous studies that were 2-dimensional simulation, we found out from sophisticated simulations that need dense mesh points to compare base pressure and velocity profile after from base with experimental data. Therefore, we focus on high order spatial disceretization over 3rd order with TVD such as MUSCL TVD 3rd, 5th, and WENO 5th order, and Limiters such as minmod, Triad. Moreover, we enforce to flux averaging schemes such as Roe, RoeM, HLLE, AUSMDV. In present, one dimensional result of Euler tests, there are Sod, Lax, Shu-Osher and interacting blast wave problems. AUSMDV as a flux averaging scheme with MUSCL TVD 5th order as spatial resolution is good agreement with exact solutions than other combinations. We are carrying out the same approaches into 3-dimensional base flow only candidate flux schemes that are Roe, AUSMDV. Additionally, turbulence models are used in 3-dimensional flow, one is Menter s SST DES model and another is Sparlat-Allmaras DES/DDES model in Navier-Stokes equations.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제15권4호
• /
• pp.291-306
• /
• 2011

We consider the second order Strang splitting method to approximate the solution to the KdV equation. The model equation is split into three sets of initial value problems containing convection and dispersal terms separately. TVD MUSCL or MUSCL scheme is applied to approximate the convection term and the second order centered difference method to approximate the dispersal term. In time stepping, explicit third order Runge-Kutta method is used to the equation containing convection term and implicit Crank-Nicolson method to the equation containing dispersal term to reduce the CFL restriction. Several numerical examples of weakly and strongly dispersive problems, which produce solitons or dispersive shock waves, or may show instabilities of the solution, are presented.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2005년도 제24회 춘계학술대회논문집
• /
• pp.343-346
• /
• 2005

전산 해석과 병렬처리를 이용하여 정사각 관 내부를 전파하는 데토네이션파의 삼차원 구조를 살펴보았다. 가변 비열비 공식과 간략화 된 일단계 Arrhenius 반응 모델과 연계된 비점성 유체 방정식을 MUSCL 기반 TVD 해법과 4단계 Runge-Kutta 적분 방법을 이용하여 해석하였다. 삼차원에서의 비정상 해석 결과로부터 그을음 막 기록(smoked-foil record)에서 같은 길이와 다른 폭을 가지는 수평 및 대각 방향 불안정에 의한 상세한 파면 구조를 파악할 수 있었다.

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

A high resolution scheme for solving gas-liquid two-phase flows with cavitation is described. This scheme uses the curvilinear coordinate grid and solves the density based momentum equations for mixture of gas-liquid medium with a preconditioning method to treat both compressible and incompressible flow characteristics. The present preconditioned method is based on the Runge-Kutta explicit finite-difference scheme, and is improved by using the diagonalization, the flux difference splitting and the MUSCL-TVD schemes to save computational effort and to increase stability and resolvability, especially at gas-liquid contact surfaces. A homogeneous equilibrium cavitation model is used to treat the gas-liquid two-phase medium in cavitating flow as a locally homogeneous pseudo-single-phase medium. Therefore, it is easy to solve cavitating flow, including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Some numerical results obtained by the present scheme are shown.

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

본 연구의 목적은 수소 및 탄화수소를 이용한 bluff-body 화염 안정화 및 hyper-mixer와 같은 초음속 연소기 연료 분사기 주변의 유동해석 이다. 대규모 박리영역을 갖는 후방단 형상의 초음속 주변 유동을 해석하기 위한 검증 단계로 다양한 수치기법을 평가해 보았다. 고차 공간 정확도로 확장된 5차 TVD-MUSCL기법, Roe FDS, S-A DES/DDES 조합이 기저유동과 기저유출 유동에 좋은 결과를 보여 주었다.

• Journal of Advanced Marine Engineering and Technology
• /
• 제22권5호
• /
• pp.670-679
• /
• 1998

The solution of hyperbolic equation with wave nature has sharp discontinuties in the medium at the wave front. Difficulties encounted in the numrtical solution of such problem in clude among oth-ers numerical oscillation and the representation of sharp discontinuities with good resolution at the wave front. In this work inviscid Burgers equation and modified heat conduction equation is intro-duced as hyperboic equation. These equations are caculated by numerical methods(explicit method MacCormack method Total Variation Diminishing(TVD) method) along various Courant numbers and numerical solutions are compared with the exact analytic solution. For inviscid Burgers equa-tion TVD method remains stable and produces high resolution at sharp wave front but for modified heat Conduction equation MacCormack method is recommmanded as numerical technique.

• Water Engineering Research
• /
• 제3권1호
• /
• pp.23-30
• /
• 2002

A numerical model for the solution of two-dimensional dam break problems using fractional step method is developed on unstructured grid. The model is based on second-order Weighted Averaged Flux(WAF) scheme with HLLC approximate Riemann solver. To control the nonphysical oscillations associated with second-order accuracy, TVD scheme with SUPERBEE limiter is used. The developed model is verified by comparing the computational solutions with analytic solutions in idealized test cases. Very good agreements have been achieved in the verifications.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2003년도 제20회 춘계학술대회 논문집
• /
• pp.41-44
• /
• 2003

현재까지 초음속 제트에 관한 많은 연구들이 수행되어져 왔지만, 초음속 제트의 3차원 구조에 관해서는 상세히 알려져 있지 않다. 본 연구에서는 수치해석을 통하여 초음속 제트의 3차원 구조를 상세히 연구한다. 수치계산에서는 TVD법을 비정상, 3차원, 압축성 오일러 방정식에 적용하였으며, 얻어진 수치해석 결과를 이용하여, 초음속 제트의 유동장을 가시화 하였다. 본 수치해석 결과로부터 초음속 제트의 3차원 구조는 2차원적 구조와 매우 다르며, 압력비에 크게 의존한다는 것을 알았다.

• 한국방재학회:학술대회논문집
• /
• 한국방재학회 2008년도 정기총회 및 학술발표대회
• /
• pp.233-236
• /
• 2008

To analyze the saline intrusion in the place, such as an estuary, the three-dimensional numerical model is developed. In this study, the advection terms of the governing equations are discretized by upwind scheme. By using an explicit scheme for the longitudinal direction and an implicit scheme for the vertical direction, the numerical model is free from the restriction of temporal step size caused by a relatively small grid ratio. The equation of state is used to consider the density, and the scalar transport equation for salinity is employed the third order TVD to scheme to prevent unphysical oscillation near discontinuity. In order to verify saline intrusion, the numerical model is conducted to compare the previous model in the lock exchange. The present model generally show a good agreement with the previous one.

• 천문학회지
• /
• 제29권spc1호
• /
• pp.459-460
• /
• 1996