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

An Euler/Navier Stokes solver has been developed for the analysis of steady and unsteady flows. The $q-{\omega}$ turbulent model has been incorporated into the solver in strongly coupled manner for stability and robustness. A new Chimera hole cutting algorithm, Cut-paste algorithm, has been devised for automatic Chimera hole cutting. Number of viscous/inviscid numerical computations demonstrate the accuracy and the versatility of the solver.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.63-65
• /
• 2003

General classes of boundary-pressure-driven flows of incompressible Newtonian fluids in three­dimensional (3D) channels and in 3D pipes with known steady laminar realizations are investigated respectively. The characteristic physical and geometrical quantities of the flows are subsumed in the kinetic Reynolds number Re and a parameter $\psi$, which involves the energetic ratio and the directions of the boundary-driven part and the pressure-driven part of the laminar flow. The solution of non-stationary dimension-free Navier-Stokes equations is sought in the form $\underline{u}=u_{L}+U,\;where\;u_{L}$ is the scaled laminar velocity and periodical conditions are prescribed for U in the unbounded directions. The objects of our numerical investigations are autonomous systems (S) of ordinary differential equations for the time-dependent coefficients of the spatial Stokes eigenfunction, where these systems (S) were received by application of the Galerkin-method to the dimension-free Navier-Stokes equations for u.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.112-122
• /
• 1996

A finite element scheme using the concept of PISO method has been developed to solve the Navier-Stokes viscous flows in all speed range. This scheme includes development of new pressure equation that retains both the hyperbolic term related with the density variation and the elliptic term reflecting the incompressibility constraint. The present method is applied to the incompressible two-dimensional driven cavity flow problems(Re=100, 400 and 1,000). For compressible flows, the Carter plate problem(M=3 and Re=1,000) is computed. Finally, we have simulated the shock-boundary layer interaction(M=2 and Re=2.96×10/sup 5/), a more difficult problem, and compared its results with the experiment to demonstrate the shock capturing capability of the present solution algorithm.

• 유체기계공업학회:학술대회논문집
• /
• 2000.12a
• /
• pp.169-175
• /
• 2000

A design method for transonic turbine blades is developed based on Navier-Stokes equations. The present computing process is done on the four separate steps, 1.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. Applications are made to the VKI transonic rotor blades. Computed results are analyzed with respect to the aerodynamic performance and are compared with the experimental data.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.9
• /
• pp.34-40
• /
• 2005

In present study, a two dimensional unsteady Navier-Stokes solver has been developed using the Diagonalized ADI (DADI) method and implicit dual time stepping method. The jacobian matrices in steady state Navier-Stokes equations are introduced from inviscid flux terms. The implicit treatment of artificial dissipation terms results in a block penta-diagonal matrix system and it becomes a scalar penta-diagonal matrix by diagonalization. In steady state equations about fictitious time, a new residual including a real time derivative term is introduced. From a converged solution about fictitious time, a real time unsteady solution can be obtained, which is called 'implicit dual time stepping method'. For code validation, an oscillating flat plate, a regular Karman vortices past a circular cylinder and shock buffeting around a bicircular airfoil problems are numerically solved. And they are compared with a theoretical solution, experiments and other researcher's computations.

• Journal of Navigation and Port Research
• /
• v.28 no.8
• /
• pp.753-757
• /
• 2004

The free-surface flow is simulated to make clear the viscous interaction of stem waves and the sub-breaking phenomena around a high speed vehicle. The Navier-Stokes equation is solved by a finite difference method where the body-fitted coordinate system, the wall function and the triple-grid system are invoked They are applied to study precisely on the stem flow of S-103 as to which extensive experimental data are available. Computations are extended to the submerged revolutional body. The numerical result shows that the gradient of M/Us is greatly influenced by the submerged depth And the stem wave is influenced by the separation due to the bow wave.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.5
• /
• pp.755-761
• /
• 1987

Two-Dimensional incompressible laminar boundary layer with the reversed flow region is computed using the parially parabolized Navier-Stokes equations in primitive variables. The velocities and the pressure are explicity coupled in the difference equation and the resulting penta-diagonal matrix equations are solved by a streamwise marching technique. The test calculations for the trailing edge region of a finite flat plate and Howarth's linearly retarding flows demonstrate that the method is accurate, efficient and capable of predicting the reversed flow region.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.123-141
• /
• 1996

In this paper the upwind flux difference splitting Navier-Stokes method has been applied to study the perfect gas and the equilibrium chemically reacting hypersonic flow over an axisymmetric sphere-cone(5°) geometry. The effective gamma(γ), enthalpy to internal energy ratio was used to couple chemistry with the fluid mechanics for equilibrium chemically reacting air. The test case condition was at altitude(30km) and Mach number(15). The equilibrium shock thickness over the blunt body region was much thinner than that of perfect gas shock. The pressure difference between perfect gas and equilibrium gas was about 3 ∼ 5 percent. The heat transfer coefficient were also calculated. The results were compared with VSL results in order to validate the current numerical analysis. The results from current method were compared well VSL results ; however, not well at near nose. The proper boundary condition and grid system will be studied to improve the solution quality.

• Journal of computational fluids engineering
• /
• v.4 no.2
• /
• pp.67-73
• /
• 1999

In the present paper, three types of the flow solvers were used to investigate the influence on the airfoil shape optimization. The adopted equations, i.e., Euler, thin layer Navier-Stokes and full Navier-Stokes ones. are solved using implicit LU-ADI decomposition scheme. The gradient projection method with the sinusoidal function was used as an optimization algorithm. The present numerical method was applied to the drag minimization problems under the initial shape of NACA0012 airfoils.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.13-16
• /
• 2008

After the collapse of the Tacoma bay bridge at Tacoma Washington, the accurate prediction of aerodynamics became crucial to the sound design of bridges. CFD(Computational Fluid Dynamics) becomes important tool for the prediction on wind effects on the bridge due to the recent development of CFD. The usage of CFD is further prompted by the advantages in using CFD, such as low-cost and fast feed-back of design. In this paper, an unsteady compressible Reynolds averaged Navier-Stokes code is used for the computation of the flow over bridges. Coakley's ��q-${\omega}$ �� two-equation turbulence model is used for the turbulent eddy viscosity. For accurate and stable computations, the local preconditioning method is adapted to the code. Aerodynamic characteristics of a couple bridges are presented to show the validity and the accuracy of the method.