• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.8
• /
• pp.1122-1133
• /
• 2003

Numerical characteristics of implicit upwind schemes, such as upwind ADI, line Gauss-Seidel (LGS) and point Gauss-Seidel (LU) algorithms, for Navier-Stokes equations have been investigated. Time-derivative preconditioning method was applied for efficient convergence at low Mach/Reynolds number regime as well as at large grid aspect ratios. All the algorithms were expressed in approximate factorization form and von Neumann stability analysis was performed to identify stability characteristics of the above algorithms in the presence of high grid aspect ratios. Stability analysis showed that for high aspect ratio computations, the ADI and LGS algorithms showed efficient damping effect up to moderate aspect ratio if we adopt viscous preconditioning based on min-CFL/max-VNN time-step definition. The LU algorithm, on the other hand, showed serious deterioration in stability characteristics as the grid aspect ratio increases. Computations for several practical applications also verified these results.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.81-89
• /
• 2006

Two dimensional turbulent flow simulations on the low Mach number - high Reynolds number flow about the NACA 4412 airfoil are carried out as the airfoil approaches a ground. It has turned out that angle of attack between 2 and 8 degrees is recommended for the airfoil to utilize the benefit of ground effect. For the large angle of attack, the increment of lift due to the ground effect is faded away and negative aerodynamic effect such as destabilizing aspect in static longitudinal stability occurs and the separation point moves to forward as the airfoil approaches a ground.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.199-205
• /
• 2006

Two dimensional turbulent flow simulations on the low Mach number - high Reynolds number flow about the NACA 4412 airfoil are carried out as the airfoil approaches a ground. It has been turned out that angle of attack between 2 and 8 is recommended for the airfoil to utilize the benefit of ground effect. For the large angle of attack, the increment of lift due to the ground effect is eliminated and negative aerodynamic effect such as destabilizing aspect in static longitudinal stability are occurred as the airfoil approaches a ground.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.106-112
• /
• 1999

In this article, combination of the FAS-FMG multi-grid method and the Krylov subspace method was presented in solving two dimensional driven-cavity flows. Three algorithms of the Krylov subspace method, CG, CGSTAB(Bi-CG Stabilized) and GMRES method were tested with MILU preconditioner. As a smoother of the pressure correction equation, the MILU-CG is recommended rather than MILU-GMRES(k) or MILU-CGSTAB, since the MILU-GMRES(k) preconditioner has too much computation on the coarse grid compared to the MILU-CG one. As for the momentum equation, relatively cheap smoother like SIP solver may be sufficient.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.178-183
• /
• 2000

An efficient Newton-GMRES algorithm is presented for computing two-dimensional steady compressible viscous flows on unstructured hybrid meshes. The scheme is designed on cell-centered finite volume method which accepts general polygonal meshes. Steady-state solution is obtained with pseudo-transient continuation strategy. The preconditioned, restarted general minimum residual(GMRES) method is employed in matrix-free form to solve the linear system arising at each Newton iteration. The incomplete LU fartorization is employed for the preconditioning of linear system. The Spalart-Allmars one equation turbulence model is fully coupled with the flow equations to simulate turbulence effect. The accuracy, efficiency and robustness of the presently developed method are demonstrated on various test problems including laminar and turbulent flows over flat plate and airfoils.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.2
• /
• pp.69-78
• /
• 2009

Cavitating flow simulation is of practical importance for many hydraulic engineering systems, such as pump, turbine, nozzle, injector, etc. In the present work, a solver for cavitating flow has been developed and applied to simulate the flows past axisymmetric cylinders. Governing equations are the two-phase Navier-Stokes equations, comprised of continuity equation of liquid and vapor phase. The momentum equation is in the mixture phase. The solver employed an implicit, dual time, preconditioned algorithm in curvilinear coordinates. Computations were carried out for three axisymmetric cylinders: hemispherical, ogive, and caliber-0 forebody shape. Then, the present calculations were compared with experiments and other numerical results to validate the present solver. Also, the code has shown its capability to accurately simulate the re-entrant jet phenomena and ventilated cavitation. Hence, it has been found that the present numerical code has successfully accounted for cavitating flows past axisymmetric cylinders.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.8
• /
• pp.661-669
• /
• 2007

Study on the feasibility assessment for applying the overlap grid method to numerical calculations on the ground effect has been performed. The objective of the present study is to settle the problem in the grid generation process. A low Mach number preconditioned turbulent flow solver using the overlap grid and the multi-block grid methods has been developed and applied to the ground effect simulation around the RAE 101 airfoil. It has been verified that the overlap grid method not only can provide sufficiently accurate solutions but also work out the grid generation problem in the ground effect simulations.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.193-198
• /
• 2006

It takes a lot of time and effort to generate grids for numerical analysis of problems with ground effect because the relative attitude and height of airfoil should be maintained to the ground as well as the inflow. A low Mach number preconditioned turbulent flow solver using the overlap grid technique has been developed and applied to the ground effect simulation. It has been validated that the present method using the multi-block grid gives us highly accurate solutions comparing with the experimental data of the RAE 101 airfoil in an unbounded condition. Present numerical method has been extended to simulate ground effect problems by using the overlapped grid system to avoid tedious work in generating multi-block grid system. An extended method using the overlapped grid has been verified and validated by comparing with results of multi-block method and experimental data as well. Consequently, the overlapped grid method can provide not only sufficiently accurate solutions but also the efficiency to simulate ground effect problems. It is shown that the pressure and aerodynamic centers move backward by the ground effect as the airfoil approaches to the ground.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.251-255
• /
• 2008

Flow field around helicopter involves incompressible flow near the blade root and compressible flow at the blade tip. A problem occurs for low Mach number flow due to the stiffness of the governing equations. Time-derivative preconditioning techniques have been incorporated to reduce the stiffness that occurs at low speed region. The preconditioned form of the compressible Navier-Stokes and Euler equations is used. Computations are performed for the Caradonna-Tung's hovering and non-lifting forward flight case. Computational results are in good agreement with the experimental data.