• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.9
• /
• pp.34-40
• /
• 2005

In present study, a two dimensional unsteady Navier-Stokes solver has been developed using the Diagonalized ADI (DADI) method and implicit dual time stepping method. The jacobian matrices in steady state Navier-Stokes equations are introduced from inviscid flux terms. The implicit treatment of artificial dissipation terms results in a block penta-diagonal matrix system and it becomes a scalar penta-diagonal matrix by diagonalization. In steady state equations about fictitious time, a new residual including a real time derivative term is introduced. From a converged solution about fictitious time, a real time unsteady solution can be obtained, which is called 'implicit dual time stepping method'. For code validation, an oscillating flat plate, a regular Karman vortices past a circular cylinder and shock buffeting around a bicircular airfoil problems are numerically solved. And they are compared with a theoretical solution, experiments and other researcher's computations.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.1000-1002
• /
• 2017

Here, the unsteady Navier-Stokes solver has been developed using implicit dual time stepping method. The implicit dual time stepping method introduced the pseudo time step for solving the new residual including the steady state residual and real time derivative. For the validation of code, Stokes 2nd problem, the laminar flow on the oscillating flat plate was selected and compare the calculating results with analytic solutions. The calculating velocity profile and skin friction has a good agreement with analytic solutions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.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.

• Journal of computational fluids engineering
• /
• v.9 no.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.

• Journal of computational fluids engineering
• /
• v.4 no.3
• /
• pp.44-51
• /
• 1999

A modified version of SST turbulence model is suggested to simulate unsteady separated flows over oscillating airfoils. The original SST model, which shows good performance in predicting various steady flows, often results in oscillatory behavior of aerodynamic loads in large separated flow regions. It is shown that this oscillatory behavior is due to the adoption of the absolute value of vorticity in generalizing the original model. As a remedy, a modification is made such that the vorticity in the original SST model is replaced by strain rate. The present model is verified for a mild separated airfoil flow at fixed angle of incidence and for unsteady flowfields about oscillating airfoils. The results are compared with BSL model and original SST model. It is illustrated that the present model gives a better agreement with the experimental results than other two models.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.1
• /
• pp.24-31
• /
• 2006

In the present study, a strongly coupled analysis code is developed for transonic flutter analysis. For aerodynamic analysis, two dimensional Reynolds-Averaged Navier-Stokes equation was used for governing equation, and ε-SST for turbulence model, DP-SGS(Data Parallel Symmetric Gauss Seidel) Algorithm for parallelization algorithm. 2 degree-of-freedom pitch and plunge model was used for structural analysis. To obtain flutter response in the time domain, dual time stepping method was applied to both flow and structure solver. Strongly coupled method was implemented by successive iteration of fluid-structure interaction in pseudo time step. Computed results show flutter speed boundaries and limit cycle oscillation phenomena in addition to typical flutter responses - damped, divergent and neutral responses. It is also found that the accuracy of transonic flutter analysis is strongly dependent on the methodology of fluid-structure interaction as well as on the choice of turbulence model.

• Journal of Satellite, Information and Communications
• /
• v.10 no.2
• /
• pp.34-39
• /
• 2015

In this paper, a design of an offset Cassegrain antenna is proposed for Ku and Ka dual-band without increasing the antenna size. For Efficiency of computation and implementation, the frequency selectivity surface (FSS) of reflecting the Ka-band signal and passing the Ku-band is provided for the sub-reflector instead of the main reflector. The proposed FSS hyperboloid sub-reflector is the periodic structure of a unit cell comprising octagon metal rings embedded in the multiple layers. The proposed design is verified for 19 GHz and 45 GHz bands by the use of precise electromagneitc-field simulations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.6
• /
• pp.788-802
• /
• 1998

A numerical method using dual time stepping and preconditioning procedure for efficient computations of unsteady low speed compressible flow problems is developed. The time-derivative preconditioning method which is valid at low speed flow conditions cannot maintain temporal accuracy because of the modification of the time-derivative term in Navier-Stokes equations. The dual time stepping procedure is incorporated to enable the time accurate computations and this procedure introduces a pseudo-time derivative in addition to the physical time derivative. At a given physical time, an inner iteration can be carried out until a steady state in pseudo-time is achieved. This will effectively yield a time accurate solution. Computational capabilities of the above algorithm are demonstrated through computation of a variety of practical fluid flows and it is shown that the algorithms is efficient in the essentially incompressible flows and low Mach number compressible flows with heat source.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• 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.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2002.04a
• /
• pp.59-59
• /
• 2002

램제트 추진기관은 압축과정을 별도의 부품 없이 형상에 의해서 감속하여 연소 압력비를 얻는다. 따라서 구동 마하수와 형상에 의해 흡입과정의 압축 효율이 결정된다. 설계점은 충분한 유량을 확보 할 수 있는 유량과 충격파 각을 조절하여 전압력 손실을 줄이도록 고려되어야 한다. 또한 연소가 일어나면 연소실 압력이 배압으로 작용하고 비행시에 받음각은 변하므로 이에 따른 성능 분석도 고려 되어야 할 사항이다. 본 연구는 국내에서 실험한 형상에 대해 수치계산을 수행하여 코드의 검증과 아울러 램제트 유동장의 수치적 시뮬레이션도 설계단계에서 하나의 도구로 이용할 수 있음을 보여준다. 실험에서는 배압 조건을 얻기 위해 유동 블록키지를 유로 내에 두어 상응하는 배압을 얻었지만 본 계산에서는 압력 경계조건을 직접 사용하였다. 유동이 비정상 특성을 가지므로 시간 정확도를 이차로 가지도록 이중시간 전진법을 사용하였다. 사용한 압력비는 충격파가 카울 끝에 닿는 임계상태에 가까운 12, 13, 14에 대해 계산을 수행하였고 부스터모드로 흡입구 끝이 막혀 있다가 램제트 모드로 바뀌어 연소실 압력이 위의 압력비라고 가정할 때의 비정상 천이 과정을 계산해 보았다. 본 계산은 흡입구 부분만을 떼어놓고 적절한 가정 하에 수행되었지만 연소실 내부도 비정상 특성을 가지므로 흡입구와 연소실을 동시에 같이 계산해야한다. 추후에 전체적인 계산을 진행하기 위한 전 단계로 흡입구 계산만을 수행하여 실험과 잘 일치하는 계산 결과를 얻었고 전체 계산을 위한 연구는 진행 중에 있다.