• 천문학회보
• /
• 제38권2호
• /
• pp.101.2-101.2
• /
• 2013

We present a new computational fluid dynamic (CFD) simulation code. The code employs the moving and polyhedral unstructured mesh scheme, which is known as a superior approach to the conventional SPH (smoothed particle hydrodynamics) and AMR (adaptive mesh refinement) schemes. The code first generates unstructured meshes by the Voronoi tessellation at every time step, and then solves the Riemann problem for surfaces of every Voronoi cell to update the hydrodynamic states as well as to move former generated meshes. For the second-order accuracy, the MUSCL-Hancock scheme is implemented. To increase efficiency for generating Voronoi tessellation we also develop the incremental expanding method, by which the CPU time is turned out to be just proportional to the number of particles, i.e., O(N). We will discuss the applications of our code in the context of cosmological simulations as well as numerical experiments for galaxy formation.

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

Several high resolution schemes such as OSHER, MUSCL, SMART, GAMMA, WACEB and CUBISTA are comparatively studied with respect to the accurate capturing of fluid interfaces throughout the application to two typical test cases of a translation test and a collapsing water column problem with a return wave. 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 and an interface capturing method, which is based on the finite-volume technique and fully conservative. The calculated results show that SMART scheme gives the best performance with respect to accuracy and robustness.

• International Journal of Fluid Machinery and Systems
• /
• 제2권3호
• /
• pp.248-253
• /
• 2009

This paper presents a numerical study for unsteady flows in a high-pressure steam turbine with a partial admission stage. Compressible Navier-Stokes equations are solved by the high-order high-resolution finite-difference method based on the fourth-order compact MUSCL TVD scheme, Roe's approximate Riemann solver, and the LU-SGS scheme. The SST-model is also solved for evaluating the eddy-viscosity. The unsteady two-dimensional flows through whole nozzle-rotor cascade channels considering a partial admission are numerically investigated. 108 nozzle passages with two blockages and 60 rotor passages are simultaneously calculated. The influence of the flange in the nozzle box to the lift of rotors is predicted. Also the efficiency of the partial admission stage changing the number of blockages and the number of nozzles is parametrically predicted.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 추계 학술대회논문집
• /
• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2005년도 춘계 학술대회논문집
• /
• 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년도 춘계 학술대회논문집
• /
• 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.

• 한국추진공학회지
• /
• 제13권3호
• /
• pp.27-33
• /
• 2009

본 논문에서는 초음속 공동유동에서 발생하는 압력진동에 미치는 수직 다공벽의 영향을 조사하기 위하여 수치해석적 연구를 수행하였다. 본 연구에서는 2차원 사각공동내부에 수직다공벽을 설치하여, 기류 마하수를 1.50, 1.83, 2.50로 변화시켰으며, 다공벽의 기공율을 변화시켰다. 수치계산에서는 2차원 비정상 압축성 Navier-Stokes 방정식을 수치적으로 풀기 위하여, TVD 유한 차분 MUSCL법을 사용하였다. 본 수치계산 결과에 의하면, 본 연구에서 적용된 수직 다공벽은 공동 상류단에 발생하는 비정상 전단층의 특성을 상당히 변화시켰으며, 공동내부에서 발생하는 압력진동을 크게 줄이는 것으로 알려졌다. 이와 같은 수직다공벽을 이용한 피동제어 효과는 기류마하수 그리고 다공벽의 기공율에 의존하는 것으로 나타났다.

• 한국추진공학회지
• /
• 제13권1호
• /
• pp.27-33
• /
• 2009

충격파관 유동의 측정결과는 격막의 유한 파막 시간에 의하여 영향을 받게 된다. 그러나 지금까지 이에 관한 구체적인 연구사례는 많지 않다. 본 연구에서는 저압력비에서 작동하는 충격파관의 유동을 조사하기 위하여, 실험 및 수치해석을 수행하였다. 수치해석에서는 2차원 비정상 압축성 Navier-Stokes 방정식에 TVD MUSCL 유한차분법을 적용하였다. 본 수치해석 결과는 충격파관의 실험결과를 잘 예측하였으며, 충격파관의 파막 과정을 시간의 함수로 적절히 모사할 수 있었다. 본 연구에서는 유한의 파막 시간으로 인하여 발생하는 Non-centered 팽창파의 특성을 정량화하기 위하여 가상중심의 개념을 적용하였다. 본 연구로부터 충격파관의 압력비가 증가할수록 파막 시간은 감소하였으며, 충격파관 유동에 미치는 파막 시간의 영향은 저압력비에서 현저하게 나타나게 됨을 알았다.

• 한국자동차공학회논문집
• /
• 제6권2호
• /
• pp.12-20
• /
• 1998

Exter ballistics of a typical high-speed projectile is studied through a flow-visualization experiment and an unstructured grid Navier-Srokes computation. Experiment produced a schlieren photograph that adequately shows the characteristic features of this complex flow, namely two kinds of oblique cone shocks and turbulent wake developing into the downstream. A hybrid scheme of finite volume-element method is used to simulate the compressible Reynolds-Averaged Navier-Stok- es solution on unstructured grids. Osher's approximate Riemann solver is used to discretize the cinvection term. Higher-order spatial accuracy is obtained by MUSCL extension and van Albada ty- pe flux limiter is used to stabilize the numerical oscillation near the solution discontinuity. Accurate Gakerkin method is used to discretize the viscous term. Explict fourth-order Runge-Kutta method is used for the time-stepping, which simplifies the application of MUSCL extension. A two-layer k-$\varepsilon$ turbulence model is used to simulate the turbulent wakes accurately. Axisymmetric folw and two-dimensional flow with an angle of attack have been computed. Grid-dependency is also checked by carrying out the computation with doubled meshes. 2-D calculation shows that effect of angle of attack on the flow field is negligible. Axi-symmetric results of the computation agrees well with the flow visualization. Primary oblique shock is represented within 2-3 meshes in numerical results, and the varicose mode of the vortex shedding is clearly captured in the turbulent wake region.

• 한국항공우주학회지
• /
• 제30권8호
• /
• pp.21-29
• /
• 2002

본 연구는 다공벽과 공동을 사용한 피동제어법을 천음속 습공기 유동에서 발생하는 충격파와 경계층 간섭에 적응하였다. 지배방정식은 액적성장 방정식과 완전히 결합된 2차원, 비정상, 압축성 Navier-Stokes 방정식이며, 3차 오더 MUSCL 타입의 TVD 기법을 사용하였다. 또 난류모델로는 Baldwin-Lomax 모델을 적용하였다. 본 연구에서 적용한 제어법의 유용성을 조사하기 위해 유동의 전압손실과 충격파 변위의 시간의존성 거동을 해석하였다. 수치계산 결과로부터 본 연구의 피동제어기법을 통해 천음속 습공기 유동에서 발생하는 충격파/경계층 간섭으로 인한 전압손실이 상당히 감소하였고, 익에서 충격파 운동을 억제하는 것으로 나타났다. 또 다공영역의 위치가 본 연구의 제어법의 효과에 상당한 영향을 준다는 것이 발견하였다.