• 한국추진공학회지
• /
• 제8권1호
• /
• pp.27-37
• /
• 2004

예조건화된 오일러 방정식의 수렴특성에 대한 연구를 수행하였다. 지배방정식의 거동을 이해하기 위하여 섭동 해석을 수행하였다. 중앙부에 10％ 원호를 가진 2차원 관을 통과하는 다양한 마하수의 비점성 유동장에 대해 수치 계산을 수행하였다. 공간차분은 Roe의 FDS를 사용하고 시간적분은 LU-SGS 기법을 사용하였다. 압력 및 속도의 수렴특성은 마하수와 상관없이 일정하게 유지되었으나, 온도의 수렴성은 마하수가 작아질수록 악화되는 것으로 나타났다. 섭동 해석을 통해 이러한 지배방정식의 수렴특성을 설명할 수 있었으며, 수렴특성이 예조건화 행렬의 거동 특성에 의해 결정된다는 사실을 알 수 있었다.

• 한국항공우주학회지
• /
• 제35권12호
• /
• pp.1067-1074
• /
• 2007

예조건화 오일러 방정식의 수렴특성을 개선하기 위하여 밀도에 대한 온도의 미분항의 크기를 조절하는 온도예조건화 방법을 제시하고자 한다. 다양한 마하수의 원형 턱이 있는 이차원 관내 유동을 계산하였다. 온도예조건화에 맞는 새로운 Roe의 FDS 기법의 수치 소산항을 유도하여 사용하였다. 온도예조건화는 예조건화 오일러 방정식의 수렴 속도를 증가시키는 것으로 나타났다

• 한국항공우주학회지
• /
• 제36권2호
• /
• pp.115-122
• /
• 2008

예조건화 오일러 방정식의 수렴 특성에 미치는 특성 조건수의 영향을 조사하였다. 그리고 Choi와 Merkle 예조건화를 적용한 경우와 온도 예조건화를 적용한 경우의 수렴 특성을 분석하였다. 공간차분을 위해 예조건화 Roe의 FDS 기법을 적용하였고, 시간적분을 위해 예조건화 LU-SGS 기법을 적용하였다. 지배방정식의 수렴 특성은 특성 조건수에 크게 영향을 받는 것으로 나타났으며, 최적의 특성 조건수가 존재하는 것으로 나타났다. 그리고 최적의 특성 조건수는 Choi 와 Merkle 예조건화를 적용한 경우와 온도 예조건화를 적용한 경우가 서로 다른 것으로 나타났다.

• 한국전산유체공학회지
• /
• 제10권1호
• /
• pp.99-106
• /
• 2005

Performance of first, second and third order accurate methods for calculation of in viscid fluxes in fluid flow governing equations are investigated here. For the purpose, an upwind method based on Roe's scheme is used to solve 2-dimensional Euler equations. To increase the accuracy of the method two different schemes are applied. The first one is a second and third order upwind-based algorithm with the MUSCL extrapolation Van Leer (1979), based on primitive variables. The other one is an upwind-based algorithm with the Chakravarthy extrapolation to the fluxes of mass, momentum and energy. The results show that the thickness of shock layer in the third order accuracy is less than its value in second order. Moreover, applying limiter eliminates the oscillations near the shock while increases the thickness of shock layer especially in MUSCL method using Van Albada limiter.

• 한국전산유체공학회지
• /
• 제12권1호
• /
• pp.43-52
• /
• 2007

A comparative study on flux functions for the 2-dimensional Euler equations has been conducted. Explicit 4-stage Runge-Kutta method is used to integrate the equations. Flux functions used in the study are Steger-Warming's, van Leer's, Godunov's, Osher's(physical order and natural order), Roe's, HLLE, AUSM, AUSM+, AUSMPW+ and M-AUSMPW+. The performance of MUSCL limiters and MLP limiters in conjunction with flux functions are compared extensively for steady and unsteady problems.

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

A comparative study on flux functions for the 2-dimensional Euler equations has been conducted. Explicit 4-stage Runge-Kutta method is used to integrate the equations. Flux functions used in the study are Steger-Warming's, van Leer's. Godunov's, Osher's(physical order and natural order), Roe's, HILE, AUSM, AUSM+ and AUSMPW+. The performance of MUSCL limiters and MLP limiters in conjunction with flux functions are compared extensively for steady and unsteady problems.

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

An Euler solver is developed to predict accurate aerodynamic data such as lift coefficient, drag coefficient, and moment coefficient. The conservation law form of the compressible Euler equations are used in the generalized curvilinear coordinates system. The Euler solver uses a finite volume method and the second order Roe's flux difference splitting scheme with min-mod flux limiter to calculate the fluxes accurately. An implicit scheme which includes the boundary conditions is implemented to accelerate the convergence rate. The multi-block grid is integrated into the flow solver for complex geometry. The flowfields are analyzed around NACA 0012 airfoil in the cases of $M_{\infty}=0.75,\;\alpha=2.0\;and\;M_{\infty}=0.80,\;\alpha=1.25$. The numerical results are compared with other numerical results from the literature. The final goal of this research is to prepare a robust and an efficient Navier-Stokes solver eventually.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 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 Mechanical Science and Technology
• /
• 제20권11호
• /
• pp.1961-1971
• /
• 2006

A soft recovery system (SRS) is a device that stops a high speed projectile without damaging the projectile. The SRS is necessary to verify the shock resistant requirements of microelectronics and electro-optic sensors in smart munitions, where the projectiles experience over 20,000 g acceleration inside the barrel. In this study, a computer code for the performance evaluation of a SRS based on ballistic compression decelerator concept has been developed. It consists of a time accurate compressible one-dimensional Euler code with use of deforming grid and a projectile motion analysis code. The Euler code employs Roe's approximate Riemann solver with a total variation diminishing (TVD) method. A fully implicit dual time stepping method is used to advance the solution in time. In addition, the geometric conservation law (GCL) is applied to predict the solutions accurately on the deforming mesh. The equation of motion for the projectile is solved with the four-stage Runge-Kutta time integration method. A small scale SRS to catch a 20 mm bullet fired at 500 m/s within 1,600 g-limit has been designed with the proposed method.

• International Journal of Aeronautical and Space Sciences
• /
• 제3권1호
• /
• pp.74-85
• /
• 2002

The numerical experiment has been conducted to investigate the unsteady shock wave reflecting phenomena. The cell-vertex finite-volume, Roe's upwind flux difference splitting method with unstructured grid is implemented to solve unsteady Euler equations. The $4^{th}$-order Runge-Kutta method is applied for time integration. A linear reconstruction of the flux vector using the least-square method is applied to obtain the $2^{nd}$-order accuracy for the spatial derivatives. For a better resolution of the shock wave and slipline, the dynamic grid adaptation technique is adopted. The new concept of grid adaptation technique, which is much simpler than that of conventional techniques, is introduced for the current study. Three error indicators (divergence and curl of velocity, and gradient of density) are used for the grid adaptation procedure. Considering the quality of the solution and the numerical efficiency, the grid adaptation procedure was updated up to $2^{nd}$ level at every 20 time steps. For the convenience of comparison with other experimental and analytical results, the case of interaction between the straight incoming shock wave and a sharp wedge is simulated for various flow conditions. The numerical results show good agreement with other experimental and analytical results, in the shock wave reflecting structure, slipline, and the trajectory of the triple points. Some critical cases show disagreement with the analytical results, but these cases also have been proven to show hysteresis phenomena.