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

• KSTLE International Journal
• /
• v.5 no.2
• /
• pp.52-56
• /
• 2004

The hard disk drive performance depends strongly on air bearing characterisitcs of the head slider. The objective of the slider design is to provide accurate positioning of the magnetic read/write element at the very small height above the disk. Application of the numerical methods is required due to complexity of the air bearing surface of the slider. The Boundary-Fitted Coordinate System Divergence Formulation method can be used for calculation of pressure distribution in the case of steep film thickness gradients. In the present work, the interpolating functions used in the expression for the Couette flow are modified in order to improve the solution characteristics in the extremely high compressibility number region. The advantages of the modified method are demonstrated on example of the flat skewed slider. Finally, the modi.ed method is applied to analysis of the static characteristics of the femto-slider. The analysis results indicate the effect of the silder's air bearing surface crown on the flying height and the pitching angle in steady state position.

• Tribology and Lubricants
• /
• v.16 no.6
• /
• pp.432-439
• /
• 2000

The present work is an attempt to calculate the steady state pressure and perturbed pressure of herringbone grooved journal bearings. A generalized coordinate system is introduced to handle the complex bearing geometry. The coordinates are fitted to the groove boundary and the Reynold's equation is transformed to be fitted to this coordinate system using the Gauss divergence theorem. This method makes it possible to deal with an arbitrary configuration of a lubricated surface. The caharacteristics of finite herringbone groove journal bearing are well calculated using this method.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2002.10b
• /
• pp.105-106
• /
• 2002

An expanded scheme of direct numerical solution method for solving the Reynolds' equation in the boundary fitted coordinate systems for the gas lubrication with ultra low clearance is presented. Skewed slider is calculated by this scheme and results are compared to the original direct numerical solution. The modified scheme has advantages in stability in high compressibility number region. At the lower A region the difference in results of original and modified method is several percents.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.3
• /
• pp.125-131
• /
• 1998

The paper deals with the application of Boundary-Fitted Coordinate System (BFCS) to the two wave models of parabolic and hyperbolic types developed on a rectangular grid system. Since the BFCS conforms the boundaries of the region in such wary that boundary conditions or calculation process can be accurately represented, improvement in predicting the wave fields can be achieved. The numerical results show a good agreement with the analytical results for either waves propagating or reflecting along a circular channel of constant depth. Simulation of reflecting waves in a parabolic wave model is accomplished by the backward calculation as if waves approached at the cross wall take a turn in the opposite direction and propagate against a channel.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 1999.06a
• /
• pp.317-324
• /
• 1999

The present work is an attempt to calculate the steady state pressure and perturbed pressure of herringbone grooved journal bearings. A generalized coordinate system is introduced to handle the complex bearing geometry. The coordinates are fitted to the groove boundary and the Reynold's equation is transformed to be fitted to this coordinates system using the Gauss divergence theorem. This method makes it possible to deal with an arbitrary configuration of a lubricated surface. The characteristics of finite herringbone grooved journal are well calculated using this method.

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

Turbulent flows around tube banks and in the diffuser were studied using a non-orthogonal boundary fitted coordinate system and the modified K-.epsilon. turbulence model. In these cases, many problems emerge which stem from the geometrical complexity of the flow domain and the physical complexity of turbulent flow itself. To treat the complex geometry, governing equations were reformulated in a non-orthogonal coordinate system with Cartesian velocity components and discretised by the finite volume method with a non-staggered variable arrangement. The modified K-.epsilon. model of Hanjalic and Launer was applied to solve above two cases under the condition of strong and mild pressure gradient. The results using the modified K-.epsilon. model results in both test cases.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1996.11a
• /
• pp.27-41
• /
• 1996

In this study the heat transfer and fluid flow of the molten pool in stationary gas tungsten arc welding using argon shielding gas were investigated. Transporting phenomena from the welding arc to the base material surface, such as current density, heat flux, arc pressure and shear stress acting on the weld pool surface, were taken from the simulation results of the corresponding welding arc. Various driving forces for the weld pool convection were considered, self-induced electromagnetic, surface tension, buoyancy, and impinging plasma arc forces. Furthermore, the effect of surface depression due to the arc pressure acting on the molten pool surface was considered. Because fusion boundary has a curved and unknown shape during welding, a boundary-fitted coordinate system was adopted to precisely describe the boundary for the momentum equation. The numerical model was applied to AISI 304 stainless steel and compared with the experimental results.

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

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