• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.10
• /
• pp.3293-3303
• /
• 1996

Numerical analyses have been performed to estimate the absorption heat and mass transfer coefficients in absorption process of the LiBr aqueous solution and the total heat and mass transfer rates in a vertical tube absorber which is coolING ed by air. Axisymmetric cylindrical coordinate system was adopted to model the circular tube and the transport equations were solved by the finite volume method. Absorption behaviors of heat and mass transfer were analyzed through falling film of the LiBr aqueous solution contacted by water vapor in tube. Effects of film Reynolds number on heat and mass transfer coefficients have been also investigated. Especially, effects of tube diameter have been considered to observe the total heat and mass transfer rates through falling film along the tube. Based on the analysis it has been found that the total mass transfer rate increases rapidly in a region with low film Reynolds number(10 ~ 40) as the film Reynolds number increases, while decreases beyond that region. The total heat and mass transfer rates increase with increasing the tube diameter.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.155-158
• /
• 2002

In general, TVC(Thermal Vapor Compressor) is used to boost/compress a low pressure vapor to a higher pressure for further utilization. The one-dimensional method is simple and reasonably accurate, but cannot realize the detail as like the back flow and recirculation in the mixing chamber, viscous shear effect, and etc. In this study, the axisymmetric How simulations have been performed to reveal the detailed flow characteristics for the various ejector shapes. The Navier-Stokes and energy equations are solved together with the continuity equation In the compressible flow fields. The standard $k-{\epsilon}$ model is selected for the turbulence modeling. The commercial computational fluid dynamic code FLUENT software is used for the simulation. The results contain the entrainment ratio under the various motive, suction and discharge pressure conditions. The numerical results are compared with the experimental data, and the comparison shows the good agreement. The three different flow regimes (double chocking, single chocking and back flow) have been clearly distinguished according to each boundary pressure values. Also the effects of the various shape variables (nozzle position, nozzle outlet diameter, mixing tube diameter, mixing tube converging angle, and etc.) are quantitatively discussed.

• Journal of Mechanical Science and Technology
• /
• v.18 no.2
• /
• pp.313-324
• /
• 2004

A rich phenomenon in the dynamics of azimuthal vortices in a circular cylinder caused by the inertial oscillation is investigated numerically at high Reynolds numbers and moderate Rossby numbers. In the actual spin-up flow where both the Ekman circulation and the bottom friction effects are included, the first appearance of a seed vortex is generated by the Ekman boundary-layer on the bottom wall and the subsequent roll-up near the corner bounded by the side wall. The existence of the small vortex then rapidly propagates toward the inviscid region and induces a complicated pattern in the distribution of azimuthal vorticity, i.e. inertial oscillation. The inertial oscillation however does not deteriorate the classical Ekman-pumping model in the time scale larger than that of the oscillatory motion. Motions of single vortex and a pair of vortices are further investigated under a slip boundary-condition on the solid walls. For the case of single vortex, repeated change of the vorticity sign is observed together with typical propagation of inertial waves. For the case of a pair of vortices with a two-step profile in the initial azimuthal velocity, the vortices' movement toward the outer region is resisted by the crescent-shape vortices surrounding the pair. After touching the border between the core and outer regions, the pair vortices weaken very fast.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.7 no.4
• /
• pp.46-52
• /
• 2003

This paper presents the numerical thermal analysis for regeneratively cooled rocket thrust chambers. An integrated numerical model incorporates computational fluid dynamics for the hot-gas thermal environment, and thermal analysis for the liner and coolant channels. The flow and temperature fields in rocket thrust chambers is assumed to be axisymmetric steady state which is presumed to the combustion liner. The heat flux computed from nozzle flow is used to predict the temperature distribution of the combustion liner As a result, we present the wall temperature of combustion liner and the temperature change of coolant.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.14 no.3
• /
• pp.9-15
• /
• 2010

Numerical investigation was carried out on axisymmetric over-expanded rocket nozzle to predict flow fields of transitional shock separation patterns. The unsteady, compressible N-S equations with k-$\omega$ SST for turbulence model closure were solved using a fully implicit finite volume scheme. Computed results were in good agreement with previous experimental works. It was found that strong side-loads were generated during the transition of RSS to FSS due to the development of a vortex ring in the inviscid jet core region. Hysteresis phenomenon exhibited by the shock-separation patterns was also found during the start-up and shut-down processes.

• Journal of ILASS-Korea
• /
• v.5 no.1
• /
• pp.13-22
• /
• 2000

The present study has numerically analyzed the vaporization characteristics of fuel droplets in the high temperature convective flow field. The axisymmetric governing equations for mass, momentum, energy, and species are solved by an iterative and implicite unstructured finite-volume method. The moving boundary due to vaporization is handled by the deformable unstructured grid technique. The pressure-velocity coupling in the density-variable flows is treated by the SIMPLEC algorithm. In terms of the matrix solver, Bi-CGSTAB is employed for the numerically efficient and stable convergence. The n-decane is used as a liquid fuel and the initial droplet temperature is 300K. Computations are performed for the nonevaporating and evaporating droplets with the relative interphase velocity(25m/s). The unsteady vaporization process has been simulated up to the nondimensional time, 25. Numerical results indicate that the mathematical model developed in this study succesfully simulates the main features of the droplet vaporization process in the convective environment.

• Journal of Ship and Ocean Technology
• /
• v.1 no.2
• /
• pp.19-30
• /
• 1997

The accurate prediction of the flow and the pressure distribution near the tip of the blade is crucial in determining the tip vortex cavitation inception which usually occurs on the blade tip or inside the core of the tip vortex just downstream of the blade tip. An improved boundary element method is applied to the prediction of the flow around propeller blades, with emphasis at the tip region. In the method, the Blow adapted grid and a higher order panel method, which combines a hyperboloidal panel geometry with a hi-quadratic dipole distribution, are used in order to accurately model the trailing wake geometry and the highly rolled-up regions in the wake. The method is applied to several propeller geometries and the results have been found to agree well to the existing experimental data. Inviscid flow methods are able to predict the pressures at the tip as well as the shape of the trailing wake. On the other hand, they are unable to determine the flow inside the viscous core of the tip vortex, where cavitation inception often occurs. Thus, a method is presented that treats the flow inside the viscous core. The inner flow is treated with a 2-D Clavier-stokes solution without making any assumptions for axisymmetric flow and conicity of the flow along the tip trajectory. The method can thus allow the treatment of general propeller blade configurations. The velocity and pressure distributions inside the core are shown and compared to those from other numerical methods.

• Journal of the Korean Society of Safety
• /
• v.21 no.6 s.78
• /
• pp.7-13
• /
• 2006

A transient heat transfer and thermal stress analysis for a brake drum of commercial vehicles have been performed by ANSYS code in the cases of single braking and the repeated braking condition. The temperature and thermal stress distributions in the brake drum under various braking conditions were obtained using a two-dimensional axisymmetric model. In case of deceleration of 0.3G with an initial vehicle speed of 60km/h, the maximum temperature in the drum was $87.6^{\circ}C$ after braking application. The maximum stress of 78.7MPa in the drum occurred at the intersection between the flange and hat under a condition in which repeated 15 cycles braking with an initial vehicle speed of 60km/h and a deceleration of 0.3G is applied to according to KS R1129. The maximum stress value is much lower than the yield strength of drum material(FC250).

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.10
• /
• pp.2712-2723
• /
• 1994

Numerical simulation of an axisymmetric ethylene-air jet diffusion flame has been carried out in order to investigate flame dynamics and soot formation. The model solves the time-dependent Navier-Stokes equations and includes models for soot formation, chemical reaction, molecular diffusion, thermal conduction, and radiation. Numerically FCT(Flux Corrected Transport) and DOM(Discrete Ordinate Method) methos are used for convection and radiation trasport respectively. Simulation was conducted for a 5 cm/sec fuel jet flowing into a coflowing air stream. The maximum flame temperature was found to be approximately 2100 K, and was located at an axial position of approximately 5 cm from the base of the flame. The maximum soot volume fraction was about $7{\times}10^{-7}$, and was located within the high temperature region where the fuel mole fraction ranges from 0.01 to 0.1. The buoyancy-driven low-frequency(12~13 Hz) structures convected along the outer region of the flame were captured. In case without radiation trasport, the maximum temperature was higher by 150 K than in case with radiation. Also the maximum soot volume fraction reached about $8{\times}10^{-6}$. As the the hydrocarbon fuel forms many soot particles, the radiation transport becomes to play a more important role.

• Journal of computational fluids engineering
• /
• v.18 no.4
• /
• pp.69-73
• /
• 2013

Mass manufacturing of optical fiber includes the process of very thin glass fiber drawing by heating and softening the high purity silica preform and applying the draw tension on the softened tip of preform neck-down profile in a draw furnace. In this computational study, this process is numerically modeled with simplified geometry of the draw furnace which is comprised of essential parts such as concentric graphite heater, muffle tube, and insulation surrounding the heater. The iterative computational scheme is employed between one-dimensional model of neck-down profile prediction and two-dimensional axisymmetric thermo-fluid CFD computation of radiative heating and working gas convection. The computational results show the experimentally observed neck-down profile in heated section of preform, while yielding the reasonable values of draw tension and heater wattage. Also, this study analyzes and discusses the effects of heating conditions such as heater length and temperature on several important aspects of glass fiber drawing process.