• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.392-394
• /
• 2011

This paper is concerned about the hydrodynamic coefficients of hydrofoil. We discretized the incompressible Navier-Stokes equation with second order Runge-kutta for the time in the second order compact scheme for the spatial. The three-dimensional CFD code based on hybrid mesh on the finite volume method is used to simulate flow around NACA series foils. Lift and drag coefficient is calculated for several NACA series foils using different mesh types. Our aim is to obtain the lift and drag coefficient to evaluate the robustness of the solver and to shaw the advantage of using trim tab at the trailing edge. It concludes with a discussion of results and recommendations for future work.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.181-189
• /
• 1998

A numerical study based on the three-dimensional Reynolds averaged Navier-Stokes equations is presented to analyze the transonic flowfield through two-stage axial compressor. Explicit four-step Runge-Kutta scheme is used for solution algorithm, and local time step and implicit residual averaging are introduced for enhancing the convergency. Artificial dissipation model is adopted to assure the stability of solution. The solver is coupled with Baldwin-Lomax model to describe turbulence. To avoid calculating the unsteady flow, a mixing process is modeled at a station between rotating and stationary blade rows. Results show a variety of important physical phenomena. Comparison of the flowfields with and without tip clearance shows that the effect is considerable in this flowfield. Comparisons with experimental data carried out to validate the calculational results show reasonable agreements. Some remedies are also suggested to improve the revealed problems.

• Kyungpook Mathematical Journal
• /
• v.52 no.2
• /
• pp.167-177
• /
• 2012

In this article, we propose exponentially fitted error correction methods(EECM) which originate from the error correction methods recently developed by the authors (see [10, 11] for examples) for solving nonlinear stiff initial value problems. We reduce the computational cost of the error correction method by making a local approximation of exponential type. This exponential local approximation yields an EECM that is exponentially fitted, A-stable and L-stable, independent of the approximation scheme for the error correction. In particular, the classical explicit Runge-Kutta method for the error correction not only saves the computational cost that the error correction method requires but also gives the same convergence order as the error correction method does. Numerical evidence is provided to support the theoretical results.

• 한국연소학회:학술대회논문집
• /
• 2006.04a
• /
• pp.205-212
• /
• 2006

Three-dimensional numerical study was carried out for the investigation of the detonation wave structures propagating in tubes. Fluid dynamics equations and conservation equation of reaction progress variable were analyzed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Chemical reaction was modeled by using a simplified one-step irreversible kinetics model. The variable gas properties between unburned and burned states were considered by using variable specific heat ratio formulation. The unsteady computational results in three-dimension show the detailed mechanisms of rectangular and diagonal mode of detonation wave instabilities resulting same cell length but different cell width in smoked-foil record. The results for the small reaction constant shows the spinning mode of three-dimensional detonation wave dynamics, which was rarely observed in the previous numerical simulation of the detonation waves.

• Journal of Power System Engineering
• /
• v.9 no.1
• /
• pp.23-29
• /
• 2005

Numerical analysis of three-dimensional viscous flow-fields in the turbine rotor passages was carried out to investigate flow physics including the interaction between secondary vortices, tip leakage vortex, and the rotor wake. The blade tip geometry was accurately modeled adopting the embedded H grid system. An explicit four-stage Runge-Kutta scheme was used for the time integration of both the mean flow and turbulence equations. The computational results for the entire turbine rotor flows, particularly the tip clearance flow and the secondary flows, were interpreted and compared with the experimental data from the Penn State turbine stage. The predictions for major features of the flow field have been found to be in good agreement with the experimental data.

• Kyungpook Mathematical Journal
• /
• v.60 no.3
• /
• pp.629-645
• /
• 2020

A uniformly convergent numerical method is developed for solving singularly perturbed 1-D parabolic convection-diffusion problems. The developed method applies a non-standard finite difference method for the spatial derivative discretization and uses the implicit Runge-Kutta method for the semi-discrete scheme. The convergence of the method is analyzed, and it is shown to be first order convergent. To validate the applicability of the proposed method two model examples are considered and solved for different perturbation parameters and mesh sizes. The numerical and experimental results agree well with the theoretical findings.

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

A numerical simulation method has been under development for solving turbulent flows around a ship model in maneuvering motion using the Reynolds Averaged Navier-Stokes equations. The method used second-order finite differences, collocated grids, pressure-Poisson equation and four-stage Runge-Kutta scheme as key components of the solution method. A modified Baldwin-Lomax model is used for the turbulence closure. This paper presents a preliminary result of the computational study on turbulent flows past a ship model in drift motion. Calculations are carried out for a Series 60 $C_B=0.6$ ship model, for which detailed experimental data are available. The results of the present calculations are compared with the experimental data for hydrodynamic forces acting on the model as well as velocity distributions at longitudinal sections. Only fair agreements has been achieved. The computational results show the complex asymmetrical shear flow patterns including three-dimensional separations followed by formation of bilge vortices both in bow and stern regions.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.171-175
• /
• 2003

Numerical analysis of three-dimensional viscous flow-fields in the turbine rotor passages is carried out to investigate flow physics including the interaction between secondary vortices, tip leakage vortex, and the rotor wake. The blade tip geometry is accurately modeled adopting the embedded H grid topology. An explicit four-stage Runge-Kutta scheme is used for the time integration of both the mean flow and turbulence equations. The computational results for the entire turbine rotor flows, particularly the tip clearance flow and the secondary flows, are interpreted and compared with the experimental data from the Penn State turbine stage. Good agreement between the experimental data and the numerical prediction was achieved in the sense of the major features of the flow fields.

• Journal of Hydrospace Technology
• /
• v.2 no.2
• /
• pp.1-13
• /
• 1996

This paper presents the result of a computational study on the wake characteristics of two tanker models, i.e. HSVA and Mystery hull forms. The focus of the study is on the distributions of axial, radial and tangential velocities of the two hull forms in way of the propeller, especially over the propeller disk. The effect of bilge vortices on the velocity distribution is also concerned. For the computation of stern and wake flows of the two hull farms, the incompressible Reynolds-Averaged Wavier-Stokes(RANS) equations are numerically solved by the second order finite difference method, which employs a four stage Runge-Kutta scheme with a residual averaging technique and the Baldwin-Lomax model. The calculated pressure distributions on the hull surface and the axial, radial and tangential velocity distributions over the propeller disk are presented for the two hull forms. Finally, the result of wake analysis for the computed wake distribution over the propeller disk is given in comparison with those for the experimental wake distribution fur the both hull forms.

• 유체기계공업학회:학술대회논문집
• /
• 1999.12a
• /
• pp.27-32
• /
• 1999

A numerical study based on the three-dimensional thin-layer Navier-Stokes solver is carried out to analyze the flowfield through a single stage transonic compressor. Explicit four-step Runge-Kutta scheme with spatially variable time step and implicit residual smoothing is used. The governing equations are discretized with exploit finite difference method. Mixed-out average method is used at the interface between rotor and stator. And, an artificial dissipation model is used to assure the stability of solution. The results with k-$\omega$ turbulence model were compared to the results with Baldwin-Lomax model, and physical phenomena of transonic compressor are presented. The two turbulence models give the results that show reasonably good agreements with experimental data.