• Journal of ILASS-Korea
• /
• v.11 no.2
• /
• pp.89-97
• /
• 2006

An inviscid axisymmetric model capable of predicting droplet bouncing and the detailed pre-impact motion, influenced by the ambient pressure, has been developed using boundary element method (BEM). Because most droplet impact simulations of previous studies assumed that a droplet was already in contact with the impacting substrate at the simulation start, the previous simulations could not accurately describe the effect of the gas compressed between a failing droplet and the impacting substrate. To properly account for the surrounding gas effect, an effect is made to release a droplet from a certain height. High gas pressures are computationally observed in the region between the droplet and the impact surface at instances just prior to impact. The current simulation shows that the droplet retains its spherical shape when the surface tension energy is dominant over the dissipative energy. When increasing the Weber number, the droplet surface structure is highly deformed due to the appearance of the capillary waves and, consequently, a pyramidal surface structure is formed; this phenomenon was verified with our experiment. Parametric studies using our model include the pre-impact behavior which varies as a function of the Weber number and the surrounding gas pressure.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.640-645
• /
• 2001

The present study is to investigate the characteristics of the impulse waves discharged from the exit of the convergent and divergent pipes. An experiment is carried out using a shock tube with an open end and is compared to the computation of the axisymmetric, compressible, unsteady Euler equations, which are solved by the second-order total variation diminishing(TVD) scheme. For the computational work, some initial compression waves are assumed inside the pipe so that those are identical to the experimental ones of the shock tube. The results show that the peak pressures of the impulse waves discharged from the exit of convergent and divergent pipes decrease with an increase in the wavelength of the initial compression wave. All of the impulse waves have a strong directivity toward the pipe axis, regardless of the exit type of the pipe employed. The impulse waves discharged from the divergent pipe are stronger than those from the straight pipe, while the impulse waves of the convergent pipe are weaker than those from the straight pipe. It is believed that the convergent pipe can playa role of a passive control to reduce the peak pressure of the impulse wave. The present computations represent the experimented impulse waves with a good accuracy.

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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.7
• /
• pp.1039-1046
• /
• 2002

Co-flow axisymmetric laminar premixed flame of methane was used to study the influence of air temperature and $N_2$ addition on the flame structure, temperature field and emission characteristics. OH 2-D images and temperatures along the centerline were measured experimentally by PLIF and CARS techniques respectively to observe the influences of dilution and thermal effects of $N_2$ in the gas mixture. Also, the concentration of NOx was measured at each condition by gas analyser to see the suppression effect of N2 addition on NOx emissions. It was found that OH concentrations distribute widely as air temperature goes higher, while the effect of $N_2$ addition is not significant. But $N_2$ addition highly contributes to the flame front and NOx emissions which was argued to be due to the reduction of flame temperature. In accordance with experimental study, numerical simulation using CHEMKlN code was carried out to compare the temperature results with those acquired by CARS measurement, and we could find that there is good agreement between those results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.6
• /
• pp.778-787
• /
• 1999

In this study, the autoignition process of liquid fuel, injected into hot and stagnant air in a 2-D axisymmetric confined cylindrical combustor, has been investigated. Eulerian-Lagrangian scheme was adopted to analyze the two-phase flow and combustion. The unsteady conservation equations were used to solve the transition of the gas field. Interactions between two phases were accounted by using the particle source in cell (PSI-Cell) model, which was used for detailed consideration of the finite rates of transports between phases. And infinite conduction model was adopted for the vaporization of droplets. The results have shown that the process of the autoignition consists of heating up of droplets, vaporization, mixing and ignition. The ignition criteria could be determined by the temporal variations of temperature, reaction rate and species mass fraction. And the effects of various parameters on ignition phenomena are examined. These have shown that the increasing the reaction rate and/or the vaporization rate can reduce the ignition delay time.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.8
• /
• pp.1087-1096
• /
• 2001

The dynamic structures of unsteady CH$_4$/Air jet diffusion flame with a flame-vortex interaction were numerically investigated. A timed-dependent, axisymmetric computational model and a low mach number approximation were employed in the present calculation. A two-step global reaction mechanism which considers 6 species, was used to calculate the reaction rates. The predicted results including the gravitational effect show that the large outer vortices and the small inner vortices can be well simulated without any additional disturbances near nozzle tip. It was found that the temperature and species concentrations have deviated values even for the same mixture fraction in the flame-vortex interaction region. It was also shown that the flame surface is not deformed by the inner vortex in upstream region, while in downstream region, the flame surface is compressed or stretched by the outer vortex roll-up. The present unsteady jet flame configuration accompanying a flame-vortex interaction is expected to give good implications for the unsteady structures of turbulent flames.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.3
• /
• pp.386-398
• /
• 1998

The present experiment is conducted to investigate heat transfer characteristics on the impinging surface with secondary flows around circular nozzle jets. The changed vortex pattern around jet affects significantly the flow characteristics and heat transfer coefficients on the impinging surface. The effects of the jet vortex control are also considered with jet nozzle-to-plate distances and main jet velocities. The vortex pattern around a jet is changed from a convective instability to an absolute instability with a velocity suction ratio of the main jet and the secondary counterflow. With the absolute instability condition, the jet potential core length increases and the heat transfer on the impinging surface is increased by small scale eddies. The region of high heat transfer coefficients is enlarged with the high Reynolds number due to increasing secondary peak values. The effect of suction flows is influenced largely with collars attached the exit of the jet nozzle because the attached collar guides well the counterflow around the main jet.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.2
• /
• pp.697-709
• /
• 1996

Supersonic axisymmetric jets issuing from various kinds of nozzles with a throat diameter of a few millimeters were experimentally investigated. The exit Mach number and Reynolds number based on the throat diameter of nozzle were in the range of 1.0 ~ 5.9 and 8.4$\times$ $10^4$ ~ 2.9$\times$$10^6$, respectively. The nozzle pressure ratio was varied from 5 to 85. Present paper aims to offer fundamental information of the supersonic free-jets, with an emphasis to give data with which the shape of the free-jets can be depicted under a specified condition. Experimental data are summarized to enable an estimation of the shape of the supersonic free-jets. The result shows that the shape of free-jets is dependent on only the nozzle pressure ratio.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.490-495
• /
• 2003

In steel-making process of iron and steel industry, the purity and quality of steel can be dependent on the amount of CO contained in the molten metal. Recently, the supersonic oxygen jet is being applied to the molten metal in the electric furnace and thus reduces the CO amount through the chemical reactions between the oxygen jet and molten metal, leading to a better quality of steel. In this application, the supersonic oxygen jet is limited in the distance over which the supersonic velocity is maintained. In order to get longer supersonic jet propagation into the molten metal, a supersonic coherent jet is suggested as one of the alternatives which are applicable to the electric furnace system. It has a flame around the conventional supersonic jet and thus the entrainment effect of the surrounding gas into the supersonic jet is reduced, leading to a longer propagation of the supersonic jet. In this regard, gasdynamics mechanism about why the combustion phenomenon surrounding the supersonic jet causes the jet core length to be longer is not yet clarified. The present study investigates the major characteristics of the supersonic coherent jet, compared with the conventional supersonic jet. A computational study is carried out to solve the compressible, axisymmetric Navier-Stokes equations. The computational results of the supersonic coherent jet are compared with the conventional supersonic jets.

• Korea-Australia Rheology Journal
• /
• v.13 no.3
• /
• pp.141-148
• /
• 2001

In this research, experimental studies leave been performed on the hydrodynamic interaction between a spherical particle and a plane wall by measuring the force between the particle and wall. To approach the system as a resistance problem, a servo-driving system was set-up by assembling a microstepping motor, a ball screw and a linear motion guide for the particle motion. Glycerin and dilute solution of polyacrylamide in glycerin were used as Newtonian and non-Newtonian fluids, respectively. The polymer solution behaves like a Boger fluid when the concentration is 1,000 ppm or less. The experimental results were compared with the asymptotic solution of Stokes equation. The result shows that fluid inertia plays all important role in the particle-wall interaction in Newtonian fluid. This implies that the motion of two particles in suspension is not reversible even in Newtonian fluid. In non-Newtonian fluid, normal stress difference and viscoelasticity play important roles as expected. In the dilute solution weak shear thinning and the migration of polymer molecules in the inhomogeneous flow field also affect the physic of the problem.