• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 1996년도 제7회 학술강연회논문집
• /
• pp.107-116
• /
• 1996

Numerical simulations have been conducted to characterize unsteady flow structures in an axisymmetric supersonic inlet diffuser with sinusoidal pressure oscillations at the diffuser exit. The formulation is based on the unsteady Navier-Stokes equations and turbulence closure is achieved using a two-layer model with a too-Reynolds-number scheme for the near-wall treatment. The governing equations are formulated in an integral form, and are discretized by the four-stage Runge-Kutta scheme for temporal terms and the Harten-Yee upwind TVD scheme for convective terms. Results indicated that the inlet shock characteristics are significantly modified by acoustic oscillations originating from the combustor. The characteristics of shock/boundarv-layer interactions (such as the size of separation bubble, terminal shock shape, and vorticity intensity) are also greatly iufluenced by the shock oscillation due to acoustic waves.

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

A numerical study is carried out to investigate the transient process of combustion phenomena associated with hypersonic propulsion devices. Reynolds averaged Navier-Stokes equations for reactive flows are used as governing equations with a detailed chemistry mechanism of hydrogen-air mixture and two-equation SST turbulence modeling. The governing equations are discretized by a high order accurate upwind scheme and solved in a fully coupled manner with a fully implicit time accurate method. At first, oscillating shock-induced combustion is analyzed and the comparison with experimental result gives the validity of present computational modeling. Secondly, the model ram accelerator experiment was simulated and the results show the detailed transient combustion mechanisms. Thirdly, the evolution of oblique detonation wave is simulated and the result shows transient and final steady state behavior at off-stability condition. Finally, shock wave/boundary layer interaction in combustible mixture is studied and the criterion of boundary layer flame and oblique detonation wave is identified.

• 대한기계학회논문집B
• /
• 제26권3호
• /
• pp.375-383
• /
• 2002

The flowfield of transverse jet in a supersonic air stream subjected to shock wave turbulent boundary layer interactions is simulated numerically by Generalized Taylor Galerkin(GTG) finite element methods. Effects of turbulence are taken into account with a two-equation (k-$\varepsilon$) model with a compressibility correction. Injection pressures and slot widths are varied in the present study. Pressure, separation extents, and penetration heights are compared with experimental data. Favorable comparisons with experimental measurements are demonstrated.

• Advances in aircraft and spacecraft science
• /
• 제9권5호
• /
• pp.415-431
• /
• 2022

Secondary flows have a huge impact on losses generation in modern low pressure gas turbines (LPTs). At design point, the interaction of the blade profile with the end-wall boundary layer is responsible for up to 40% of total losses. Therefore, predicting accurately the end-wall flow field in a LPT is extremely important in the industrial design phase. Since the inlet boundary layer profile is one of the factors which most affects the evolution of secondary flows, the first main objective of the present work is to investigate the impact of two different inlet conditions on the end-wall flow field of the T106A, a well known LPT cascade. The first condition, labeled in the paper as C1, is represented by uniform conditions at the inlet plane and the second, C2, by a flow characterized by a defined inlet boundary layer profile. The code used for the simulations is based on the Discontinuous Galerkin (DG) formulation and solves the Reynolds-averaged Navier-Stokes (RANS) equations coupled with the Spalart Allmaras turbulence model. Secondly, this work aims at estimating the influence of viscosity and turbulence on the T106A end-wall flow field. In order to do so, RANS results are compared with those obtained from an inviscid simulation with a prescribed inlet total pressure profile, which mimics a boundary layer. A comparison between C1 and C2 results highlights an influence of secondary flows on the flow field up to a significant distance from the end-wall. In particular, the C2 end-wall flow field appears to be characterized by greater over turning and under turning angles and higher total pressure losses. Furthermore, the C2 simulated flow field shows good agreement with experimental and numerical data available in literature. The C2 and inviscid Euler computed flow fields, although globally comparable, present evident differences. The cascade passage simulated with inviscid flow is mainly dominated by a single large and homogeneous vortex structure, less stretched in the spanwise direction and closer to the end-wall than vortical structures computed by compressible flow simulation. It is reasonable, then, asserting that for the chosen test case a great part of the secondary flows details is strongly dependent on viscous phenomena and turbulence.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 1999년도 춘계학술대회논문집
• /
• pp.253-256
• /
• 1999

The Streamline Upwind Petrov-Galerkin(SUPG) finite element method is used to solve the two-dimensional laminar and turbulent flow. The flow is simulated by averaged Navier-Stokes equations with a penalty function approach and the lograithmic(k-$\varepsilon$) turbulent model is employed to take into account its turbulent behavior. The near-wall viscous sub-layer model is employed to approach the dominant viscous effects in the near wall zones. To find a good-enough initial guess of the Newton-Raphson iteration solving Nonlinear Matrix the Incremental method is considered for momentum and the Incomplete logarithmic turbu-lent equations for Turbulence. The validation of our method is investigated in comparision with published experimental data.

• 한국전산유체공학회지
• /
• 제2권2호
• /
• pp.26-34
• /
• 1997

An implicit viscous turbulent flow solver is developed for two-dimensional geon unstructured triangular meshes. The flux terms are discretized based on a cell-centered formulation with the Roe's flux-difference splitting. The solution is advanced in time us backward-Euler time-stepping scheme. At each time step, the linear system of equation approximately solved wi th the Gauss-Seidel relaxation scheme. The effect of turbulence is with a standard k-ε two-equation model which is solved separately from the mean flow equation the same backward-Euler time integration scheme. The triangular meshes are generated advancing-front/layer technique. Validations are made for flows over the NACA 0012 airfoil. Douglas 3-element airfoil. Good agreements are obtained between the numerical result experiment.

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

An implicit viscous turbulent flow solver is developed for two-dimensional geometries on unstructured triangular meshes. The flux terms are discretized based on a cell-centered finite-volume formulation with the Roe's flux-difference splitting. The solution is advanced in time using an implicit backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with the Gauss-Seidel relaxation scheme. The effect of turbulence effects is approximated with a standard $k-{\varepsilon}$ two-equation model which is solved separately from the mean flow equations using the same backward-Euler time integration scheme. The triangular meshes are generated using an advancing-front/layer technique. Validations are made for flows over the NACA0012 airfoil and the Douglas 3-element airfoil. Good agreements are obtained between the numerical results and the experiment.

• 대한기계학회논문집
• /
• 제13권5호
• /
• pp.999-1011
• /
• 1989

본 연구의 목적은 첫째, 재발달 경계층에서의 난류에서지 및 난류전단응력에 대한 전달방정식들의 각 항의 측정자료들을 보다 정확히 제시하고 항들간의 균형을 비교 평가함으로써 비평형 유동으로부터 평형유동으로 회복되는 과정을 검토하고, 둘째, 대표적인 난류 모델들로써 표존 k-.epsilon.모델 및 레이놀즈 응력 모델을 사용한 수치계산을 수행함으로써 이와같은 모델들이 비평형 유동을 서술함에 있어 발생될 수 있는 문제점들을 고찰하는데 있다.

• 한국화재소방학회논문지
• /
• 제13권4호
• /
• pp.20-29
• /
• 1999

본 연구는 두 가지 유형의 아트리움 공간에 대해서 Zone모델과 Field모델을 비교하였으며 특히 천장에 열 유속을 갖는 아트리움 화재에 대해서 SMEP화재 모델을 적용하여 연기거동을 수치해석 하였다. Zone 모델로는 NIST에서 개발된 CFAST 그리고 CSIRO에서 개발된 NBTC 1-room 모델을 사용하여 SMEP Field모델을 검증하였다. PISO 알고리즘과 부력항을 포함한 수정 k-e epsilon 난류모델을 사용한 SMEP은 연속, 운동, 에너지 그리고 농도 방정식을 풀었으며, 상용 Zone 모델들과의 비교는 서로 유사한 결과를 보였다. 천장이 유리로 만들어진 아트리움의 경우, 태양열에 의한 천장 열 유속을 고려함으로써 좀더 실제적인 화재현상을 규명할 수 있다. 수치해석결과 천장 열 유속 조건은 연층의 온도 분포에는 영향을 미치지만 연기의 하강과 거동에는 커다란 영향을 미치지 않고 있음을 확인 할 수 있었다. 따라서 화재 감지기나 배연 설비 시스템의 설치시 이러한 점들이 고려되어야만 한다.

• 한국자동차공학회논문집
• /
• 제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.