• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 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.

• 한국전산유체공학회지
• /
• 제13권3호
• /
• pp.35-43
• /
• 2008

A methodology for the simulation of compressible high Reynolds number flow over rigid and moving bodies on a structured Cartesian grid is described in this paper. The approach is based on a modified version of the Brinkman Penalization method. To avoid oscillations in the vicinity of the body and to simulate shcok-containing flows, a Weighted Essentially Non-Oscillatory scheme is used to discretize the spatial flux derivatives. For high Reynolds number viscous flow, two turbulence models of the two-equation Menter's SST URANS model and a two-equation Detached Eddy Simulation are implemented. Some simple flow examples are given to assess the accuracy of the technique. Finally, a moving grid capability is demonstrated.

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

A methodology for the simulation of compressible high Reynolds number flow over rigid and moving bodies on a structured Cartesian grid is described in this paper. The approach is based on a modified version of the Brinkman Penalization method. To avoid oscillations in the vicinity of the body and to simulate shcok-containing flows, a Weighted Essentially Non-Oscillatory scheme is used to discretize the spatial flux derivatives. For high Reynolds number viscous flow, two turbulence models of the two-equation Menter's SST URANS model and a two-equation Detached Eddy Simulation are implemented. Some simple flow examples are given to assess the accuracy of the technique. Finally, a moving grid capability is demonstrated.

• 한국전산유체공학회지
• /
• 제11권4호
• /
• pp.39-47
• /
• 2006

In this paper, the authors developed a multidimensional interpolation method inside a finite volume cell in the computation of high-order accurate numerical flux such as the fifth order WEND (weighted essentially non-oscillatory) scheme. This numerical method starts from a simple Taylor series expansion in a proper spatial order of accuracy, and the WEND filter is used for the reconstruction of sharp nonlinear waves like shocks in the compressible flow. Two kinds of interpolations are developed: one is for the cell-averaged values of conservative variables divided in one mother cell (Type 1), and the other is for the vertex values in the individual cells (Type 2). The result of the present study can be directly used to the cell refinement as well as the convective flux between finer and coarser cells in the Cartesian adaptive grid system (Type 1) and to the post-processing as well as the viscous flux in the Navier-Stokes equations on any types of structured and unstructured grids (Type 2).

• 대한설비공학회:학술대회논문집
• /
• 대한설비공학회 2006년도 하계학술발표대회 논문집
• /
• pp.116-121
• /
• 2006

Numerical simulations for dendritic growth of crystals are conducted in this study by the level set method. The effect of order of difference is tested for reinitialization error in simple problems and authors founded in case of 1st order of difference that very fine grids have to be used to minimize the error and higher order of difference is desirable to minimize the reinitialization error The 2nd and 4th order Runge-Kutta scheme in time and 3rd and 5th order of WENO schemes with Godunov scheme are applied for space discretization. Numerical results are compared with the analytical theory, phase-field method and other researcher's level set method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제27권4호
• /
• pp.207-231
• /
• 2023

This paper aims to improve the alternative formulation of the fifth- and sixth-order accurate weighted essentially non-oscillatory (AWENO) finite difference schemes. The first is to derive the AWENO scheme with sixth-order accuracy in the smooth region of the solution. Second, a new weighted polynomial functions combining the perturbed forms with conserved variable to the AWENO is constructed; the new form of tunable functions are invented to maintain non-oscillatory property. Detailed numerical experiments are presented to illustrate the behavior of the new perturbational AWENO schemes. The performance of the present scheme is evaluated in terms of accuracy and resolution of discontinuities using a variety of one and two-dimensional test cases. We show that the resulted perturbational AWENO schemes can achieve fifth- and sixth-order accuracy in smooth regions while reducing numerical dissipation significantly near singularities.

• 한국추진공학회지
• /
• 제22권2호
• /
• pp.66-73
• /
• 2018

밀폐용기 내 Zirconium/Potassium Perchlorate의 연소를 수치적 모델링을 통해 전산해석을 수행하였다. 5차 WENO 공간차분법과 improved delayed detached eddy (IDDES) 난류모델을 사용하여 충격파가 동반되는 내부 유동구조를 모사하였고, 라그랑지안 연소모델을 통해 화약 입자를 계산하였다. 옆면 중앙에 센서가 설치된 원통형 밀폐용기 내부 유동분석을 통해 압력 진동이 발생하는 원인을 규명하였다. 또한 센서 다이어프램 깊이 변화에 따라 측정되는 압력 데이터를 실험값과 비교분석 하였다. 그 결과 센서 탭의 깊이가 약 2.36 mm 이상으로 커지면 유동속도가 아음속으로 감쇠하고 복잡한 eddy가 발생하여 측정값에 큰 불규칙성을 야기하는 현상을 관측하였다.

• 천문학회보
• /
• 제46권1호
• /
• pp.36.3-36.3
• /
• 2021

We study ultra-relativistic jets on several tens kpc scales through three-dimensional relativistic hydrodynamic (RHD) simulations using a new RHD code based on the weighted essentially non-oscillatory (WENO) scheme. Utilizing the high-resolution and high-accuracy capabilities of the new code, we especially explore the structures and energetics of nonlinear flows, such as shocks, turbulence, velocity shear in different parts of jets. We find that the mildly relativistic shocks which form in the jet backflow are most effective for the shock dissipation of the jet energy, while the turbulent dissipation is largest either in the backflow or in the shocked ICM, depending on the jet parameter. The velocity shear is strongest across the jet flow to the cocoon boundary. Our results should have important implications for the studies of high-energy cosmic-ray production in radio galaxies.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• pp.426-432
• /
• 2004

Large-Eddy Simulation (LES) is applied for the simulation of compressible flat plate boundary with Reynolds number up to 5 X 10$^{5}$ . Numerical examples include shock/boundary layer interaction and boundary layer transition, aiming future application to the analysis of transonic fan/compressor cascades. The present LES code uses hybrid com-pact/WENO scheme for the spatial discretization and compact diagonalized implicit scheme for the time integration. The present code successfully predicted the bypass transition of subsonic boundary layer. As for supersonic turbulent boundary layer, mean and fluctuation velocity of the attached boundary, as well as the evolution of the friction coefficient and the displacement thickness both upstream and downstream of the separation region are all in good agreement with experiment. The separation point also agreed with the experiment. In the simulation of the shock/laminar boundary layer interaction, the dependence of the transition upon the shock strength is reproduced qualitatively, but the extent of the separation region is overpredicted. These numerical examples show that LES can predict the behavior of boundary layer including transition and shock interaction, which are hardly managed by the conventional Reynolds-averaged Navier-Stokes approach, although there needs to be more effort before achieving quantitative agreement.

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

Considerable efforts are required to develop a monotone, robust and stable high-order numerical scheme for solving the hyperbolic system. The discontinuous Galerkin(DG) method is a natural choice, but elimination of the spurious oscillations from the high-order solutions demands a new development of proper limiters for the DG method. There are several available limiters for controlling or removing unphysical oscillations from the high-order approximate solution; however, very few studies were directed to analyze the exact role of the limiters in the hyperbolic systems. In this study, the performance of the several well-known limiters is examined by comparing the high-order($p^1$, $p^2$, and $p^3$) approximate solutions with the exact solutions. It is shown that the accuracy of the limiter is in general problem-dependent, although the Hermite WENO limiter and maximum principle limiter perform better than the TVD and generalized moment limiters for most of the test cases. It is also shown that application of the troubled cell indicators may improve the accuracy of the limiters under some specific conditions.