• Journal of Korean Society for Atmospheric Environment
• /
• v.16 no.4
• /
• pp.339-350
• /
• 2000

The three-dimensional new corrdinate system over a single hill double hills and complex terrain with a single hill and a rectangular obstacle was generated using a body-fitted coordinate system. Control of the coordinate line distribution in the field was executed by generalizing the elliptic generating system to Poisson equation. ▽2ξ=P. The new coordinate system was well fitted to the surface boundary of single hill and double hills. But in the case of complex terrain with hill and rectangular obstacle there was smoothing tendency around the rectangular obstacle. In order to show the validity of the body-fitted coordinate system the heat diffusion equation was transformed and the temperature distribution was calculated over the various terrain. The results showed the temperature distribution was very symmetrical and stable around hills and obstacle. As a result the couple of a body-fitted coordinate system and the heat diffusion equation were executed successfully. Wind field over complex terrain with hill and rectangular obstacle which represent urban area was simulated stably in body-fitted coordinate system. The qualitative result show the enhancement of wind speed at the upwind direction of a hill and a rectangular obstacle and the recirculation zone at the downwind direction.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.3
• /
• pp.1015-1027
• /
• 1996

In this paper, a numerical procedure for the calculation of the incompressible Navier-Stokes equations in a generalized nonorthogonal body fitted coordinate system is proposed and is validated through three test problems. Present numerical procedure derives the pressure equation by using the pressure substitution method on the regular grid system, and discretized momentum equations are based on the covariant velocity components. Cavity flow, backward facing step flow, and two dimensional channel flow with a sinusoidal wavy wall are chosen as three test problems. Numerical solutions obtained by present procedure shows a good agreement with previous numerical and/or experimental results. Convergence rate is also satisfactory.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.64-70
• /
• 1996

Three dimensional turbulent flow fields around ships are simulated by a numerical method. Reynolds Averaged Navier-Stokes equations are used where Reynolds stresses are approximated by Baldwin-Lomax and Sub-Grid Scale(SGS) turbulence models. Body-fitted coordinate system is introduced to conform three dimensional ship geometries. The governing equations are discretized by a finite volume method. Temporal derivatives are approximated by the forward differencing and the convection terms are approximated by the QUICK or Kawamura scheme. The 2nd-order centered differencing is used for other spatial derivatives. Pressure and velocity fields are simultaneously iterated by the Highly Simplified Marker-And-Cell method. To verify the numerical method and turbulence models, flow fields around ships are simulated and compared to the experiments.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.2
• /
• pp.39-47
• /
• 1997

In the present paper, viscous flow fields around a wigley hull with appendages are analysed to study interactions between the hull and appendages. Navier-Stokes and continuity equations are solved by a finite volume method in a body-fitted coordinate system which conforms three dimensional ship geometries with appendages. A Sub-Grid Scale(SGS) turbulent model is used for a calculation of high Reynolds number flow. Numerical computations has been done for a Wigley hull form at $Rn=1.0{\times}10^6$. The results show that the present approach can predict, at least in qualitative sense, the influence of the appendages upon the flow field around a ship.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.3
• /
• pp.253-258
• /
• 2012

The Cartesian grid system has generally been used in casting simulations, even though it does not represent sloped and curved surfaces very well. These distorted boundaries cause several problems, and special treatment is necessary to resolve them. A cut cell method on a Cartesian grid has been developed for the simulation of threedimensional mold filling. Cut cells at a cast/mold interface are generated on Cartesian grids, and the governing equations are computed using the volume and areas of the cast at the cut cells. In this paper, we propose a new method based on the partial cell treatment (PCT) that can consider the cutting cells which are cut by the cast and the mold. This method provides a better representation of the surface geometry, and will be used in the computation of velocities that are defined on the cell boundaries in the Cartesian gird system. Various test examples for several casting process are computed and validated.

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

• Transactions of the Korean Society of Automotive Engineers
• /
• v.1 no.1
• /
• pp.32-40
• /
• 1993

Flow fields in model engine cooling passages are studied numerically by using TURBO-3D program, a finite volume based 3-D turbulent flow program adopting a general body fitted coordinate system. The effects of exit position on mass flow rate at each gasket hole are examined for a model cooling passage in order to understand the flow distribution inside the water jacket. The results of the present study can be applied to the design of high performance, high reliability engine.

• Nuclear Engineering and Technology
• /
• v.32 no.2
• /
• pp.169-179
• /
• 2000

A detailed numerical analysis of the evolution of thermal stratification in a curved piping system in a nuclear power plant is performed. A finite volume based thermal-hydraulic computer code has been developed employing a body-fitted, non-orthogonal curvilinear coordinate for this purpose. The cell-centered, non-staggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is prevented by the application of momentum interpolation method. The SIMPLE algorithm is employed for the pressure and velocity coupling, and the convection terms are approximated by a higher-order bounded scheme. The thermal-hydraulic computer code developed in the present study has been applied to the analysis of thermal stratification in a curved duct and some of the predicted results are compared with the available experimental data. It is shown that the predicted results agree fairly well with the experimental measurements and the transient formation of thermal stratification in a curved duct is also well predicted.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.2
• /
• pp.313-324
• /
• 1992

Two-dimensional Navier-Stokes code has been developed for analysis of turbomachinery blade rows and other internal flows. The Navier-Stokes equations are written in a Cartesian coordinate system, then mapped into a generalized body-fitted coordinate system. All direction of viscous terms are incorporated and turbulent effects are modeled using the Baldwin-Lomax algebraic model. Equation are discretized using finite difference method on the C-type grids and solved using implicit LU-ADI decomposition scheme. Calculations are made at a VKI turbine cascade flow in a transonic wind-tunnel and compared to experimental data. Present numerical scheme is shown to be in good agreement with the previous experimental results and simulates the two-dimensional viscous flow phenomena.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.2
• /
• pp.575-583
• /
• 1991

The weld pool convection problem that occurs during the stationary GTA welding has been studied, considering the four driving forces for weld pool convection, i.e., the electromagnetic force, the buoyancy force, the aerodynamic drag force, and the surface tension force at the weld pool surface. In the numerical simulation, the difficulties associated with the irregular moving liquid-solid interface have been successfully overcome by adopting a Boundary-Fitted Coordinate system. In the experiments to show the validity of the numerical analysis, a deep periphery and shallow centerpentrated weld pool shape was observed from the etched specimen. It could be revealed that this type of weld pool shape could be simulated, only when some of aerodynamic drag force distributions are considered. Although slight disagreement arose, the calculated and the observed weld pool shapes were in a reasonable agreement.