• Journal of the Society of Naval Architects of Korea
• /
• v.44 no.1 s.151
• /
• pp.16-25
• /
• 2007

Studies of unsteady-airfoil flows have been motivated mostly by efforts to avoid. or reduce such undesirable effects as flutter, noise and vibrations, dynamic stall. In this paper, we carry out a computational study of viscous flows around a two-dimensional oscillating airfoil to investigate unsteady effects in these important and challenging flows. A fully implicit incompressible RANS solver has been used for calculating unsteady viscous flows around an airfoil. The cell-centered End order finite volume method is utilized to discretize governing equations. in order to ease the flow computation for fluid region changing in time, improve the qualify of solution and simplify the grid generation for an oscillating airfoil flow, the computational method adopts a moving and deforming grid generation technique based on the multi-block grid topology. The numerical method is applied for calculating viscous flows of an oscillating NACA 0012 in uniform flow. The computational results are compared with available experimental data. Computed results are compared with experimental data and flow characteristics of the experiment are reproduced well In the computed results.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.900-907
• /
• 2000

The three-dimensional turbulent flow in curved pipes susceptible to flow-accelerated corrosion has been analyzed numerically to predict the pressure and shear stress distributions on the inner surface of the pipes. The analysis employs the body-fitted non-orthogonal curvilinear coordinate system and a standard $ {\kappa}-{\varepsilon}$ turbulence model with wall function method. The finite volume method is used to discretize the governing equations. The convection term is approximated by a high-resolution and bounded discretization scheme. The cell-centered, non-staggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is prevented by the application of a modified version of momentum interpolation scheme. The SIMPLE algorithm is employed for the pressure and velocity coupling. The numerical calculations have been performed for two curved pipes with different bend angles and curvature radii, and discussions have been made on the distributions of the primary and secondary flow velocities, pressure and shear stress on the inner surface of the pipe to examine applicability of the present analysis method. As the result it is seen that the method is effective to predict the susceptible systems or their local areas where the fluid velocity or local turbulence is so high that the structural integrity can be threatened by wall thinning degradation due to flow-accelerated corrosion.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.86-95
• /
• 2008

For analyses of multi-phase flows in a water-cooled nuclear power plant, a three-dimensional SIMPLE-algorithm based hydrodynamic solver CUPID-S has been developed. As governing equations, it adopts a two-fluid three-field model for the two-phase flows. The three fields represent a continuous liquid, a dispersed droplets, and a vapour field. The governing equations are discretized by a finite volume method on an unstructured grid to handle the geometrical complexity of the nuclear reactors. The phasic momentum equations are coupled and solved with a sparse block Gauss-Seidel matrix solver to increase a numerical stability. The pressure correction equation derived by summing the phasic volume fraction equations is applied on the unstructured mesh in the context of a cell-centered co-located scheme. This paper presents the numerical method and the preliminary results of the calculations.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.201-210
• /
• 2009

A volume capturing method using unstructured grid system for numerical analysis of multiphase flows is introduced in the present paper. This method uses an interface capturing method (CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The novelty of CICSAM lies in the adaptive combination of high resolution discretization scheme which ensures the preservation of the sharpness and shape of the interface while retaining boundedness of the field, and no explicit interface reconstruction which is perceived to be difficult to implement on unstructured grid system. Several typical test cases for multiphase flows are presented, which are simulated by an in-house solution code(PowerCFD). This code employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with CICSAM. It is found that the present method simulates efficiently and accurately complex free surface flows such as multiphase flows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.6
• /
• pp.443-453
• /
• 2009

The Rayleigh-Taylor instability, the bubble rising in both partially and fully filled containers and the droplet splash are simulated by an in-house solution code(PowerCFD), which are typical benchmark problems among multiphase flows with material interface due to density difference. The present method(code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with interface capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The present results are compared with other numerical solutions found in the literature. It is found that the present method simulates efficiently and accurately complex free surface flows such as multiphase flows with material interface due to both density difference and instability.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.7
• /
• pp.2347-2355
• /
• 1996

Purpose of the present study is to develop a computational design program for shape optimization, combining the numerical optimization technique with the flow analysis code. The present methodology is then validated in three cases of aerodynamic shape optimization. In the numerical optimization, a feasible direction optimization algorithm and shape functions are considered. In the flow analysis, the Navier-Stokes equations are discretized by a cell-centered finite volume method, and Roe's flux difference splitting TVD scheme and ADI method are used. The developed design code is applied to a transonic channel flow over a bump, and an external flow over a NACA0012 airfoil to minimize the wave drag induced by shock waves. Also a separated subsonic flow over a NACA0024 airfoil is considered to determine a maximum allowable thickness of the airfoil without separation.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.13 no.4
• /
• pp.267-280
• /
• 2009

In this paper we investigate a simple two-level additive Schwarz preconditioner for the P1 symmetric interior penalty Galerkin method of the Poisson equation on rectangular meshes. The construction is based on the decomposition of the global space of piecewise linear polynomials into the sum of local subspaces, each of which corresponds to an element of the underlying mesh, and the global coarse subspace consisting of piecewise constants. This preconditioner is a direct combination of the block Jacobi iteration and the cell-centered finite difference method, and thus very easy to implement. Explicit upper and lower bounds for the maximum and minimum eigenvalues of the preconditioned matrix system are derived and confirmed by some numerical experiments.

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

An adaptive Cartesian grid method having the best elements of structured, unstructured, and Cartesian grids is developed to solve the steady two-dimensional Euler equations. The solver is based on a cell-centered finite-volume method with Roe's flux-difference splitting and implicit point Gauss-seidel time integration method. Calculations of several compressible flows are carried out to show the efficiency of the developed computer code. The results were generally in good agreements with existing data in the literature and the developed code has the good ability to capture important feature of the flows.

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

• Journal of computational fluids engineering
• /
• v.16 no.1
• /
• pp.14-21
• /
• 2011

Super-cavitating flows around under-water bodies are being studied for drag reduction and dramatic speed increase. In this paper, high speed super-cavitating flow around a two-dimensional symmetric wedge-shaped body were studied using an unsteady Reynolds-averaged Navier-Stokes equations solver based on a cell-centered finite volume method. To verify the computational method, flow over a hemispherical head-form body was simulated and validated against existing experimental data. Various computational conditions, such as different wedge angles and caviation numbers, were considered for the super-cavitating flow around the wedge-shaped body. Super-cavity begins to form in the low pressure region and propagates along the wedge body. The computed cavity lengths and velocities on the cavity boundary with varying cavitation number were validated by comparing with analytic solution.