• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.7-14
• /
• 2004

It is well known that preconditioning methods are efficient for convergence acceleration at compressible low Mach number flows. In this study, the original Euler equations and three preconditioners nondimensionalized differently are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as flux discretization and time integration respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning one produces Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.51-56
• /
• 2006

Several high resolution schemes such as OSHER, MUSCL, SMART, GAMMA, WACEB and CUBISTA are applied to two typical test cases of a translation test and a collapsing water column problem for the accurate capturing of fluid interfaces. It is accomplished by implementing the high resolution schemes in the in-house CFD code(PowerCFD) for computing 3-D flow with an unstructured cell-centered method, which is based on the finite-volume technique and fully conservative. The calculated results are found to show that SMART scheme gives the best performance with respect to accuracy and robustness.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.1
• /
• pp.10-17
• /
• 2008

In this paper an accurate and stable gridless method that can be applied to multi-dimensional convection problems is developed on a flow directional local grid. A two dimensional pure convection problem is calculated and more accurate and stable solution is obtained compared with other schemes in grid method. The tested numerical schemes include 1st-order upwind scheme, 2nd-order Leith scheme, 3rd-order MUSCL, and QUICK scheme. It is seen that more accurate results are expected when the schemes combined with a MMT control limiter.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.128-134
• /
• 2005

In the CFD area, the numerical analysis of high Mach number flow over a blunt-body poses many issues. Various numerical schemes have been developed to cover the issues, but the traditional schemes are still used widely due to the complexities of new schemes and intricacy of modifying the established codes. In the present study, the well-known Roe's FDS based on TVD-MUSCL scheme is used for the solution of very high Mach number three-dimensional flows posing carbuncle and non-physical phenomena in numerical analysis. A parametric study was carried out to account for the effects of the entropy fixing, grid configurations and initial condition. The carbuncle phenomena could be easily overcome by the entropy fixing, and the non-physical solution could be eliminated by the use of the modified initial condition regardless of entropy fixing and grid configurations.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.177-181
• /
• 2007

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma}$ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.586-592
• /
• 2000

터보머신 태부에 존재하는 정익 - 동익의 상호작용 유동현상을 수치모사 하는 코드를 병렬화 하였다 정익 - 동익의 상호작용을 해석하는 데에 편리하도륵 Multi-Block Grid System을 도입하여 계산영역을 형성하였고, 동익의 움직임으로 인해 발생하는 Sliding Interface부분은 Patched 알고리즘을 적용하여 해석하였다. 정익과 동익의 수를 1대 1로 단순화시켜 수치모사한 결과와 정익과 동익의 수를 실제 조건과 더 비슷하게 설정한 3대 4의 비율로 맞추어 수치모사한 결과를 비교하였다. 또한, 병렬컴퓨팅으로 인해 단축된 계산시간을 다른 연구에서의 계산시간들과 서로 비교하였다. 2차원 비정상 압축성 Navier-Stokes 방정식이 이용되었고, 난류모델링에는 K-w SST 모델링이 적응되었다. Roe의 FDS 기법을 사용하여 플럭스를 계산하였고, MUSCL 기법을 적용하여 3차의 공간정확도를 갖도록 하였다. 시간적분에는 이보성의 DP-SGS를 사용하였다. 해석결과의 분석에는 Time-averaged pressure distribution과 Pressure amplitude distribution 데이터를 사용했다.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.514-519
• /
• 2003

A computational study is performed to clarify the characteristics of supersonic moist air jet issuing from a simple sonic nozzle. The effects of the initial supersaturation on the Mach disk diameter and location, the barrel shock wave and jet boundary structures are investigated in details. The axisymmetric, compressible, Navier-Stokes equations, coupled with droplet growth equation, are solved using a third-order MUSCL type TVD finite-difference scheme. It is found that the Mach disk diameter increases with an increase in relative humidity of moist air. while its location is not significantly dependent on the relative humidity. As the relative humidity increases, the barrel shock wave and jet boundary are more expanded due to the local static pressure rise of nonequilibrium condensation.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.562-567
• /
• 2003

A computational study is performed to clarify the characteristics of supersonic moist air jet issuing from a simple sonic nozzle. The effects of the initial supersaturation on the Mach disk diameter and location, the barrel shock wave and jet boundary structures are investigated in details. The axisymmetric, compressible, Navier-Stokes equations, coupled with droplet growth equation, are solved using a third-order MUSCL type TVD finite-difference scheme. It is found that the Mach disk diameter increases with an increase in relative humidity of moist air. while its location is not significantly dependent on the relative humidity. As the relative humidity increases, the barrel shock wave and jet boundary are more expanded due to the local static pressure rise of nonequilibrium condensation.

• 한국연소학회:학술대회논문집
• /
• 2006.04a
• /
• pp.205-212
• /
• 2006

Three-dimensional numerical study was carried out for the investigation of the detonation wave structures propagating in tubes. Fluid dynamics equations and conservation equation of reaction progress variable were analyzed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Chemical reaction was modeled by using a simplified one-step irreversible kinetics model. The variable gas properties between unburned and burned states were considered by using variable specific heat ratio formulation. The unsteady computational results in three-dimension show the detailed mechanisms of rectangular and diagonal mode of detonation wave instabilities resulting same cell length but different cell width in smoked-foil record. The results for the small reaction constant shows the spinning mode of three-dimensional detonation wave dynamics, which was rarely observed in the previous numerical simulation of the detonation waves.

• Journal of Energy Engineering
• /
• v.24 no.2
• /
• pp.110-114
• /
• 2015

A numerical investigation has been carried out on recirculation zone formed downstream of a compressible flow over a backward facing step. The study has been performed by solving Two-Dimensional Navier-Stokes equations. The system of governing equations has been solved, using an explicit Harten-Yee Non- MUSCL Modified flux type TVD scheme and a zero-equation algebraic turbulence model to calculate the eddy viscosity coefficient. The recirculation region dimensions are characterized over a range of Mach numbers of fixed step height 5mmThe detail of recirculation zone such as pressure, temperature, recirculation length, strength etc are reported. The variations of these characteristics due to change of Mach number are also presented.