• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.299-300
• /
• 2003

• 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.

• 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.

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

• Journal of Mechanical Science and Technology
• /
• v.20 no.11
• /
• pp.1961-1971
• /
• 2006

A soft recovery system (SRS) is a device that stops a high speed projectile without damaging the projectile. The SRS is necessary to verify the shock resistant requirements of microelectronics and electro-optic sensors in smart munitions, where the projectiles experience over 20,000 g acceleration inside the barrel. In this study, a computer code for the performance evaluation of a SRS based on ballistic compression decelerator concept has been developed. It consists of a time accurate compressible one-dimensional Euler code with use of deforming grid and a projectile motion analysis code. The Euler code employs Roe's approximate Riemann solver with a total variation diminishing (TVD) method. A fully implicit dual time stepping method is used to advance the solution in time. In addition, the geometric conservation law (GCL) is applied to predict the solutions accurately on the deforming mesh. The equation of motion for the projectile is solved with the four-stage Runge-Kutta time integration method. A small scale SRS to catch a 20 mm bullet fired at 500 m/s within 1,600 g-limit has been designed with the proposed method.

• 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.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• 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.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.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.27-41
• /
• 2008

An unstructured hybrid mesh flow solver has been developed for the simulation of three-dimensional steady and unsteady incompressible flow fields. The incompressible Navier-Stokes equations with an artificial compressibility method were discretized by using a node-based finite-volume method. For the unsteady time-accurate computation, a dual-time stepping method was adopted to satisfy a divergence-free flow field at each physical time step. An implicit time integration method with local time stepping was implemented to accelerate the convergence in the pseudo-time sub-iteration procedure. The one-equation Spalart-Allmaras turbulence model has been adopted to solve high-Reynolds number flow fields. The flow solver was parallelized to minimize the CPU time and to overcome the computational overhead. This method has been applied to calculate steady and unsteady flow fields around submarine configurations and a 3-D infinite cylinder. Validations were made by comparing the predicted results with those of experiments or other numerical results. It was demonstrated that the present method is efficient and robust for the prediction of steady and unsteady incompressible flow fields.