• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.1
• /
• pp.60-66
• /
• 2000

In order to improve cavitation characteristics for a high-speed propeller, leading edge shape of a 2-D hydrofoil is investigated numerically and experimentally. For flowfield analysis around the leading edge, the incompressible Reynolds Averaged Navier-Stokes(RANS) equation is solved using the standard $k-\varepsilon$ turbulence model and a finite volume method(FVM). The cavitation thickness, which is occurred on hydrofoil surface, is predicted using the panel code. It is shown that the calculation codes predict the experimental trend fairly well. From these results, new hydrofoils are designed

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.3 no.2
• /
• pp.161-171
• /
• 1999

Aerodynamic sensitivity analysis codes are developed via the hand-differentiation using a direct differentiation method and an adjoint method respectively from discrete two-dimensional compressible Navier-Stokes equations. Unlike previous other researches, Baldwin-Lomax algebraic turbulence model is also differentiated by hand to obtain design sensitivities with respect to design variables of interest in turbulent flows. Discrete direct sensitivity equations and adjoint equations are efficiently solved by the same time integration scheme adopted in the flow solver routine. The required memory for the adjoint sensitivity code is greatly reduced at the cost of the computational time by allowing the large banded flux jacobian matrix unassembled. Direct sensitivity code results are found to be exactly coincident with sensitivity derivatives obtained by the finite difference. Adjoint code results of a turbulent flow case show slight deviations from the exact results due to the limitation of the algebraic turbulence model in implementing the adjoint formulation. However, current adjoint sensitivity code yields much more accurate sensitivity derivatives than the adjoint code with the turbulence eddy viscosity being kept constant, which is a usual assumption for the prior researches.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.15 no.3
• /
• pp.233-251
• /
• 2011

The Gauge-Uzawa method [GUM] in [9] which is a projection type algorithm to solve evolution Navier-Stokes equations has many advantages and superior performance. But this method has been studied for backward Euler time discrete scheme which is the first order technique, because the classical second order GUM requests rather strong stability condition. Recently, the second order time discrete GUM was modified to be unconditionally stable and estimated errors in [12]. In this paper, we contemplate several GUMs which can be derived by the same manner within [12], and we dig out properties of them for both stability and accuracy. In addition, we evaluate an stability condition for the classical GUM to construct an adaptive GUM for time to make free from strong stability condition of the classical GUM.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2005.11a
• /
• pp.450-454
• /
• 2005

A numerical analysis is performed to analyze the reacting flow-field of an axisymetric coaxial ramjet combustor. Two dimensional Navier-Stokes equation with low Reynolds number $k-\varepsilon$ turbulence model is utilized and finite-rate chemistry model is adopted. Eddy dissipation model is applied for a modeling of turbulent combustion. Two different types of combustors (combustor with a suddenly expanded dump and combustor with wedge-shaped flame holders) are compared in a view point of flame stabilizing.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.127-132
• /
• 1998

The present study is aimed to investigate flow characteristics of Two dimensional backward-facing step by numerical approach. A convection conservative difference scheme based upon SOLA algorithm is used for the solution of the two-dimensional incompressible Navier-Stokes equations to simulate the laminar, transitional and turbulent flow conditions at which the experimental data can be available for the backward-facing step. The twenty kinds of Reynolds number are used for the calculations. In an effort to demonstrate that the reported solutions are dependent on the mesh refinement, computations are performed on seven different meshes of uniformly increasing refinement. Also to investigate the result of inflow dependence, two kinds of the inflow profile are chosen for the laminar flow. As criterion of benchmarking the result of numerical simulation, reattachment length is used for the selected Reynolds numbers.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.133-140
• /
• 1998

A Navier-Stokes code based on a unstructured finite volume method is used to simulate the MIT flapping foil experiment. A low Reynolds number $k-{\varepsilon}$ turbulence model is used to close the Reynolds averaged Navier-Stokes equations. Computations are carried out for a domain involving two flapping foils and a downstream hydrofoil. The computational domain is meshed with unstructured quadrilateral elements, partly structured. Numerical solutions show good agreement with experiment. Unsteadiness inside boundary layer is entrained when a unsteady vortex impinge on the blade surface. It shoves that local peak value inside the boundary layer and also local minimum near the edge of boundary layer as it developes along the blade surface. The unsteadiness inside the boundary layer is almost isolated from the free stream unsteadiness and being convected at local boundary layer speed, less than the free stream value.

• Journal of computational fluids engineering
• /
• v.7 no.4
• /
• pp.19-27
• /
• 2002

An efficient implicit multigrid method is presented for the Navier-Stokes and k-ω turbulence equations. Freezing and limiting strategies are applied to improve the robustness and convergence of the multigrid method. The eddy viscosity and strongly nonlinear production terms of turbulence are frozen in the coarser grids by passing down the values without update of them. The turbulence equations together with the Navier-Stokes equations, however, are consecutively solved on the coarser grids in a loosely coupled fashion. A simple limit for k is also introduced to circumvent slow-down of convergence. Numerical results for the unseparated and separated transonic airfoil flows show that all computations converge well without any robustness problem and the computing time is reduced to a factor of about 3 by the present multigrid method.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.5
• /
• pp.53-65
• /
• 1993

In this work the dynamic response of 3-D arbitrarily shaped rigid massless foundation is numerically obtained using boundary element under non-relaxed boundary condition. The problem is formulated in time domain by the boundary element method. The fundamental solutions used in this work are the Stokes solutions of the three dimensional elastodynamics. This method has advantages over frequency domain techniques in that it provides in a natural and direct way the time history of the response and forms the basis for elct:ension to nonlinear problems. This work is verified and can be used for practical purpose.

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