• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.1
• /
• pp.193-201
• /
• 1995

The role of thermal contact resistance between a casting and a metal mold as well as natural convection in the melt during solidification of a pure metal is numerically studied. Numerical simulation is performed for a rectangular cavity using the coordinate transformation by boundary-fitted coordinate and pure aluminum is used as the phase- change material. The influences of thermal contact resistance on the interface shape and position, solidified volume fraction, temperature field and local heat transfer are investigated.

• Water for future
• /
• v.27 no.3
• /
• pp.95-105
• /
• 1994

The boundary layer flow under a sluice gate is numerically solved by the random vortex sheet method combined with the vortex-in-cell method in a boundary-fitted coordinate system. The numerical solution shows that the boundary layer developed along the vertical sluice gate wall is the primary cause for the discrepancy in the contraction ratio between the laboratory experiments and inviscid theory; the bottom boundary layer plays much a smaller role in the discrepancy. By dimensional analysis it is concluded that the discrepancy is inversely proportional to the 3/4th power of the gate opening, as analyzed by Benjamin(1956). The results of the numerical simulation and dimensional analysis show a good agreement with experimental results obtained by Benjamin(1956).

• Journal of Korean Society for Atmospheric Environment
• /
• v.14 no.3
• /
• pp.219-228
• /
• 1998

In the present investigation, a numerical model developed for the prediction of the wind flow over complex terrain is validated by comparing with the field experiments. For the solution of the Reynolds - Averaged Clavier- stokes equations which are the governing equations of the microscale atmospheric flow, the model is constructed based on the finite-volume formulation and the SIMPLEC pressure-correction algorithm for the hydrodynamic computation. The boundary- fitted coordinate system is employed for the detailed depiction of topography. The boundary conditions and the modified turbulence constants suitable for an atmospheric boundary- layer are applied together with the k- s turbulence model. The full- scale experiments of Cooper's Ridge, Kettles Hill and Askervein Hill are chosen as the validation cases . Comparisons of the mean flow field between the field measurements and the predicted results show good agreement. In the simulation of the wind flow over Askervein Hill , the numerical model predicts the three dimensional flow separation in the downslope of the hill including the blockage effect due to neighboring hills . Such a flow behavior has not been simulated by the theoretical predictions. Therefore, the present model may offer the most accurate prediction of flow behavior in the leeside of the hill among the existing theoretical and numerical predictions.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.10
• /
• pp.2418-2425
• /
• 1993

A two-dimensional transient heat transfer model for reactive gas assisted laser cutting process with a moving Gaussian heat source is developed using a numerical finite difference technique. The kerf width, melting front shape and temperature distribution were calculated by using the boundary-fitted coordinate system to handle the ejection of workpiece material and heat input from reaction and evaporation. An analytical solution for cutting front movement was adopted and numerical simulation was performed to calculate the temperature distribution and melting front thickness. To calculate the moving velocity of cutting front, the normal distribution of the cutting gas velocity was used. The kerf width was revealed to be dependent on the cutting velocity, laser power and cutting gas velocity.

• International Union of Geodesy and Geophysics Korean Journal of Geophysical Research
• /
• v.26 no.1
• /
• pp.1-14
• /
• 1998

Various numerical methods for the two dimensional shallow water equations have been applied to the problems of flood routing, tidal circulation, storm surges, and atmospheric circulation. These methods are often based on the Alternating Direction Implicity(ADI) method. However, the ADI method results in inaccuracies for large time steps when dealing with a complex geometry or bathymetry. Since this method reduces the performance considerably, a fully implicit method developed by Wilders et al. (1998) is used to improve the accuracy for a large time step. Finite Difference Methods are defined on a rectangular grid. Two drawbacks of this type of grid are that grid refinement is not possibile locally and that the physical boundary is sometimes poorly represented by the numerical model boundary. Because of the second deficiency several purely numerical boundary effects can be involved. A boundary fitted curvilinear coordinate transformation is used to reduce these difficulties. It the curvilinear coordinate transformation is used to reduce these difficulties. If the coordinate transformation is orthogonal then the transformed shallow water equations are similar to the original equations. Therefore, an orthogonal coorinate transformation is used for defining coordinate system. A multigrid (MG) method is widely used to accelerate the convergence in the numerical methods. In this study, a technique using a MG method is proposed to reduce the computing time and to improve the accuracy for the orthogonal to reduce the computing time and to improve the accuracy for the orthogonal grid generation and the solutions of the shallow water equations.

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

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.37-44
• /
• 2000

Numerical investigation is made to study the effects of natural convection on the formation of macrosegregation of a Pb-Sn alloy solidification process in a 2-D confined rectangle mold. The governing equations are calculated using previous continuum models with SIMPLE algorithm doring the solidification process. In addition. to track the solid-liquid interface with time variations. the moving boundary condition Is adopted and irregular interface shapes are treated with Boundary-Fitted Coordinate system. As the temperature reduce from the liquidus to the solidus, the liquid concentration of Sn. the lighter constituent, increases. Then the buoyancy-driven flow due to temperature and liquid composition gradients, called thermosolutal convection or double diffusion, occurs in the mushy region and forms the complicated macrosegregation maps. Related to this phnomena, effects on the macrosegregation formation depending on the cooling condition and gravity values are described.

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

Complex flow field around the Darius turbine rotating stationally are simulated by solving the three dimensional incompressible Navier-Stokes equation numerically. The rotating coordinate system is employed so that the boundary conditions on the blades of the rotor become simple. In order to impose the boundary condition on the blades precisely, the boundary fitted coordinate system is employed. Fractional step method is used to solve the basic equations. The complex flow fields due to the three dimensionality of the geometry of the turbine and the rotation of the turbine are obtained and they are visualized effectively by using the technique of the computer graphics.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.3
• /
• pp.368-375
• /
• 1987

Effective thermal conductivity of the porous nuclear fuel has been investigated numerically. Difficulties associated with irregular shape of pore have been overcome by using the Body Fitted Coordinate Systems. A computer code has been developed to solve the governing equation with appropriate boundary conditions by transforming from the Cartesian coordinates to the nonorthogonal curvilinear coordinates. The effects of the porosity have been investigated. For a convenient use of the result, a correlation equation was suggested under the assumption of circular pore. The computation results by the assumption of randomly oriented elliptic pore has been agreed more closely to existing experimental result than that by the assumption of circular pore.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 1997.10a
• /
• pp.31-38
• /
• 1997

The interest in numerically-generated, Boundary-Fitted Coordinate Systems (BFCS) arises from the need for conforming the boundaries of the region in such way that boundary conditions can be accurately represented. The parabolic approximation method in solving wave phenomena is known to have a great merit as time-saving method. However, the method shows a disagreement for the wide angle and behind the structure (omitted)