• 대한기계학회논문집B
• /
• 제21권12호
• /
• pp.1624-1634
• /
• 1997

A numerical study on two-dimensional unsteady transonic cascade flow has been performed by adopting dual time stepping and the k-.omega. turbulence model. An explicit 4 stage Runge-Kutta scheme for the compressible Navier-Stokes equations and an implicit Gauss-Seidel iteration scheme for the k-.omega. turbulence model are proposed for fictitious time stepping. This mixed time stepping scheme ensures the stability of numerical computation and exhibits a good convergence property with less computation time. Typical steady-state convergence accelerating schemes such as local time stepping, residual smoothing and multigrid combined with dual time stepping shows good convergence properties. Numerical results are presented for unsteady laminar flow past a cylinder and turbulent shock buffeting problem for bicircular arc cascade flow is discussed.

• 대한기계학회논문집B
• /
• 제23권10호
• /
• pp.1229-1239
• /
• 1999

A computational study on unsteady compressible flows has been performed by adopting a low Reynolds number $k-{\omega}$ turbulence model in conjunction with dual time stepping scheme. An explicit four-stage Runge-Kutta scheme for the Navier-Stokes equations and an approximate factorization scheme for the $k-{\omega}$ turbulence model equations are used. Computational results obtained for blade surface pressure distributions in the process of rotor-stator interaction in a turbine stage are in good agreement with extant experimental data. The effects of the wake from the stator on the boundary-layer transition over the rotor blade surface are discussed by showing that high intensity turbulence of the stator wake induces an early transition.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.59-62
• /
• 2002

In the present study, solution algorithms for the computation of unsteady flows on an unstructured mesh are presented. Dual time stepping is incorporated to achieve the 2-nd order temporal accuracy while reducing the linearization and the factorization errors associated with a linear solver. Hence, any time step can be used by only considering physical phenomena. Gauss-Seidel scheme is used to solve linear system of equations. Rigid motion and spring analogy method fur moving mesh are all considered and compared. Special treatments of spring analogy for high aspect ratio cells are presented. Finally, numerical results for oscillating wing are compared with experimental data.

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

In the present study, solution algorithms for the connotation of unsteady flows on an unstructured mesh me presented Dual time stepping is incorporated to achieve the 2-nd order temporal accuracy while reducing the linearization and the factorization errors associated with a linear solver. Hence, any time step can be used by only considering physical phenomena. Gauss-Seidel scheme is used to solve linear system of equations. Rigid motion and suing analogy method for moving mesh are all considered and compared. Special treatments of suing analogy for high aspect ratio cells are presented. Finally, numerical results for oscillating ing are compared with experimental data.

• 한국전산유체공학회지
• /
• 제9권2호
• /
• pp.52-61
• /
• 2004

In the present study, a solution algorithm for the computation of unsteady flows on unstructured meshes is presented. Dual time stepping is incorporated to achieve the second-order temporal accuracy while reducing errors associated with linearization and factorization. This allows any time step size, which is suitable for considering physical phenomena of interest. The Gauss-Seidel scheme is used to solve the linear system of equations. A special treatment based on spring analogy is made to handle meshes with high aspect-ratio cells. The present method was validated by comparing the results with experimental data and those obtained from rigid motion.

• 한국항공우주학회지
• /
• 제33권9호
• /
• pp.34-40
• /
• 2005

본 연구에서는 비정상 유동해석을 위한 CFD 코드의 개발을 위해 대각화 ADI 기법을 적용한 정상 해석기법과 내재적 이중시간 전진기법을 도입하였다. 정상상태 Navier-Stokes 방정식의 Jacobian 행렬은 비점성항에 대해서만 적용하였고 여기에 내재적 인공점성 연산자를 첨가하여 블록 5대각 행렬을 유도하였다. 시간단축을 위해 스칼라 5대각 행렬로 대체하였다. 가상시간에 대한 정상상태기법에 실시간에 대한 미분항이 포함된 새로운 잔류항을 정의하였다. 가상시간에 대해 수렴된 해로부터 실시간 해를 구하고 시간에 대해 적분을 수행하는 내재적 이중시간 전진기법을 이용한 비정상 Navier-Stokes 코드를 개발하였다. 이에 대한 검증으로 정지한 유체속에 진동하는 평판문제, 원기둥 후방의 주기적인 Karman 와류생성, 이중원호 익형주위의 충격파 진동문제등을 수치해석하여 이론치, 실험치, 타연구자의 계산결과와 비교, 분석하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.181-182
• /
• 2003

A preconditioned numerical method for gas-liquid to-phase flow is applied to solve cavitating flow. The present method employs a density based finite-difference method of dual time-stepping integration procedure and Roe's flux difference splitting approximation with MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. The method permits simple treatment of the whole gas-liquid two-phase flow field including wave propagation, large density changes and incompressible flow characteristics at low Mach number. By this method, two-dimensional internal flows through a venturi tuve and decelerating cascades are computed and discussed.

• 한국추진공학회지
• /
• 제6권2호
• /
• pp.75-84
• /
• 2002

고체 추진 로켓 내부 연소실의 비정상 유동을 수치적으로 해석하였다. 완전 보존식을 이용하여 2 차원 축-대칭 연소실 안의 연소 불안정을 해석하기 위한 수치 기법을 구성하였는데 비정상 유동을 해석하기 위한 수정된 $\kappa$-$\varepsilon$ 난류 모델이 사용되었다. 이산화한 지배 방정식은 연관된 경계 조건을 포함하여 dual time-stepping 방법으로 시적분 하였다. 정상 상태의 계산을 기반으로 연소실 내의 천이 압력파의 비정상 상태를 수치적으로 모사하기 위하여 압력 펄스 및 압력 변동을 연소실 상단에 부과하였다. 로켓 모터 연소실 내의 다양한 정상 상태 및 비정상 상태의 특성을 계산 및 해석하였다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• pp.1903-1908
• /
• 2004

A preconditioned numerical method for gas-liquid two-phase flows is applied to solve cavitating flow. The present method employs a finite-difference dual time-stepping integration procedure and the MUSCLTVD scheme. A homogeneous equilibrium cavitation model is used. The present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation, large density changes and incompressible flow characteristics at low Mach number. Some internal flows such as convergent-divergent nozzles are computed using this method. Comparisons of predicted and experimental results are provided and discussed.

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

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.