• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.35-38
• /
• 2002

The convergence acceleration methods for the compressible Wavier-Stokes equations are studied ,which are multigrid method and implicit preconditioned multistage time stepping method. In this paper, the performance of implicit preconditioning methods are studied for the full-coarsening multigrid methods on the high Reynolds number compressible flow computations. The effect of numerical flux on the convergence are investigated for the inviscid and viscous calculations.

• Journal of computational fluids engineering
• /
• v.8 no.4
• /
• pp.16-25
• /
• 2003

The design method for transonic turbine blades has been developed based on Wavier-Stokes equations. The present computing process is done on the four separate steps, i.e., determination of the blade profile, generation of the computational grids, cascade flow simulation and analysis of the computed results in the sense of the aerodynamic performance. The blade shapes are designed using the cubic polynomials under the control of the design parameters. Numerical methods for the flow equations are based on Van-Leer's FVS with an upwind TVD scheme on the finite volume. In the present study, numerical simulation has been done to investigate the effects of the design parameters on the aerodynamic peformance of the axial-flow turbine blades. Applications are made to the VKI transonic rotor blades. Computed results are analyzed with respect to four parameters and compared with the experimental data.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.2
• /
• pp.202-212
• /
• 1997

In this study, performance of the multi-stage explicit methods for numerical computation of the unsteady Navier-Stokes equations is investigated. Three methods under consideration are 1 st-, 2 nd-, and 4 th-order Runge-Kutta (R-K) methods. Compared in this estimation is stability, accuracy, and CPU time of each method. The computational codes developed are applied to the two-dimensional flow in a square cavity driven by an oscillating lid. It turned out that at Reynolds number 400, the 1 st-order R-K method is the best, while at 3200 the 2 nd-order R-K is recommended. At higher Reynolds numbers, it is conjectured that the 4 th-order R-K method will be the best algorithm among three due to its highest stability.

• Bulletin of the Korean Mathematical Society
• /
• v.46 no.4
• /
• pp.627-644
• /
• 2009

The gauge-Uzawa method which has been constructed in [11] is a projection type method to solve the evolution Navier-Stokes equations. The method overcomes many shortcomings of projection methods and displays superior numerical performance [11, 12, 15, 16]. However, we have obtained only suboptimal accuracy via the energy estimate in [11]. In this paper, we study semi-discrete gauge-Uzawa method to prove optimal accuracy via energy estimate. The main key in this proof is to construct the intermediate equation which is formed to gauge-Uzawa algorithm. We will estimate velocity errors via comparing with the intermediate equation and then evaluate pressure errors via subtracting gauge-Uzawa algorithm from Navier-Stokes equations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.1
• /
• pp.62-70
• /
• 2001

Three-dimensional dynamic stall over an oscillating wing has been analyzed by using a compressible Navier-Stokes code. The code solved the thin-layer Navirer-Stokes equations with a second-order time accuracy for a semispan wing with 0.3048m chord, a NACA 0015 airfoil section, and zero twist Computations were made for a freestream Mach number of 0.29, a chord Reynolds number of 1.95$\times$10(sup)6 and a reduced frequency equal to 0.1. Numerical results were compared with experimental data which include the hysteresis of lift, drag and moment at various wing span. The comparison reveals the quantitative as well as qualitative nature of the three-dimensional dynamic stall.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.6
• /
• pp.525-531
• /
• 2008

• Journal of Korea Water Resources Association
• /
• v.43 no.9
• /
• pp.801-811
• /
• 2010

A three-dimensional numerical model called NEWTANK is employed to investigate solitary wave run-up with an internal wave-maker on a steep slope. The numerical model solves the spatially averaged Navier-Stokes equations for two-phase flows. The LES (large-eddy-simulation) approach is adopted to model the turbulence effect by using the Smagorinsky SGS (sub-grid scale) closure model. A two-step projection method is adopted in numerical solutions, aided by the Bi-CGSTAB (Bi-Conjugate Gradient Stabilized) method to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate VOF (volume-of-fluid) method is used to track the distorted and broken free surface. A solitary wave is first internally generated and propagated over a constant water depth in the three-dimensional domain. Numerically predicted results are compared with analytical solutions and numerical errors are analyzed in detail. The model is then applied to study solitary wave run-up on a steep slope and the obtained results are compared with available laboratory measurements.

• Journal of the Korean Mathematical Society
• /
• v.50 no.5
• /
• pp.1129-1163
• /
• 2013

Based on finite element discretization, two linearization approaches to the defect-correction method for the steady incompressible Navier-Stokes equations are discussed and investigated. By applying $m$ times of Newton and Picard iterations to solve an artificial viscosity stabilized nonlinear Navier-Stokes problem, respectively, and then correcting the solution by solving a linear problem, two linearized defect-correction algorithms are proposed and analyzed. Error estimates with respect to the mesh size $h$, the kinematic viscosity ${\nu}$, the stability factor ${\alpha}$ and the number of nonlinear iterations $m$ for the discrete solution are derived for the linearized one-step defect-correction algorithms. Efficient stopping criteria for the nonlinear iterations are derived. The influence of the linearizations on the accuracy of the approximate solutions are also investigated. Finally, numerical experiments on a problem with known analytical solution, the lid-driven cavity flow, and the flow over a backward-facing step are performed to verify the theoretical results and demonstrate the effectiveness of the proposed defect-correction algorithms.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.120-126
• /
• 1991

An expanded scheme of a direct numerical solution method for solving the Navier-Stokes equation considering the modified boundary conditions for gas lubrication with ultra low clearance at high .LAMBDA. region is presented. Many examples are calculated by this scheme and their results are compared to the previous solutions using P$^{2}$H$^{[-992]}$ . This scheme has the advantages of fast calculation time and stable convergence in high .LAMBDA. region, and gives very good results in the case of fluid film thickness discontinuity.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.392-404
• /
• 1997

The flow induced by the oscillatory motion of a solid body is important in a number of practical problems. As the solid boundary oscillates harmonically, there is steady streaming motion invoked by the Reynolds stresses, which could cause extensive migration of the fluid during a period of fluid motion. We here analyzed the flow in a square cavity with an oscillating top wall for the parameters which make the time derivatives and the convective terms equally important in the entire cavity flow. The full Navier-Stokes equations are solved by the second-order time accurate Momentum Coupling Method which is devised by the authors. The particular numerical scheme does not need subiteration at each time step which is usually a required process to calculate the incompressible Navier-Stokes equations. The effect of two parameters, the Reynolds number and the frequency parameter, on the oscillatory flow has been investigated.