• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.9 s.240
• /
• pp.988-996
• /
• 2005

Flame structure and extinction mechanism of counterflow methane/air non-premixed flame diluted with nitrogen are studied by NASA 2.2 s drop tower experiments and two-dimensional numerical simulations with finite rate chemistry and transport properties. Extinction mechanism at low strain rate is examined through the comparison among results of microgravity experiment, 1D and 2D simulations with a finite burner diameter. A two-dimensional simulation in counterflow flame especially with a finite burner diameter is shown to be very important in explaining the importance of multidimensional effects and lateral heat loss in flame extinction, effects that cannot be understood using a one-dimensional flamelet model. Extinction mechanism at low strain rate is quite different from that at high strain rate. Low strain rate flame is extinguished initially at the outer flame edge, the flame shrinks inward, and finally is extinguished at the center. It is clarified from the overall fractional contribution by each term in energy equation to heat release rate that the contribution of radiation fraction with 1D and 2D simulations does not change so much and the overall fractional contribution is decisively attributed to radial conduction ('lateral heat loss'). The experiments by Maruta et at. can be only completely understood if multi-dimensional heat loss effects are considered. It is, as a result, verified that the turning point, which is caused only by pure radiation heat loss, has to be shifted towards much lower global strain rate in microgravity flame.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.6
• /
• pp.688-696
• /
• 2004

Nonlinear dynamic behavior of diffusive-thermal instability in diluted CH$_4$/O$_2$ diffusion flames is numerically investigated by adopting detailed chemistry and transport. Counterflow diffusion flame is adopted as a model flamelet. Particular attention is focused on the pulsating-instability regime, which arises for Lewis numbers greater than unity, and the instability occurs at high strain rate near extinction condition in this flame configuration. Once a steady flame structure is obtained for a prescribed value of initial strain rate, transient solution of the flame is calculated after a finite amount of strain-rate perturbation is imposed on the steady flame. Transient evolution of the flame depends on the initial strain rate and the amount of perturbed strain rate. Basically, the dynamic behaviors can be classified into two types, namely non-oscillatory decaying solution and diverging solution leading to extinction. The peculiar oscillatory solution, which has been found in the previous study adopting one-step chemistry and constant Lewis numbers, is net observed in this study, which is attributed to both convective flow and preferential diffusion effects.

• 한국연소학회:학술대회논문집
• /
• 2004.06a
• /
• pp.95-101
• /
• 2004

Dynamic behavior of diffusive-thermal instability in diluted $CH_4/O_2$ diffusion flames is numerically investigated by adopting detailed chemistry and transport. Counterflow diffusion flame is adopted as a model flamelet. Particular attention is focused on the pulsating-instability regime, which arises for Lewis numbers greater than unity, and the instability occurs at high strain rate near extinction condition in this flame configuration. Once a steady flame structure is obtained for a prescribed value of initial strain rate. transient solution of the flame is calculated after a finite amount of strain-rate perturbation is imposed Oil the steady flame. Transient evolution of the flame depends on the initial strain rate and the amount of perturbed strain rate. Basically, the dynamic behaviors can be classified into two types, namely non-oscillatory decaying solution and diverging solution leading to extinction. The peculiar oscillatory solution. which has been found in the previous study adopting one-step chemistry and constant Lewis numbers, is not observed in this study, which is attributed to both convective flow and preferential diffusion effects.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.1
• /
• pp.29-35
• /
• 2001

Radiation-induced oscillatory instability in CH$_4$/Air diffusion flames is numerically investigated by adopting detailed chemistry. Counterflow diffusion flame is employed as a model flamelet and optically thin gas-phase radiation is assumed. Attention is focused on the extinction regime induced by radiative heat loss, which occurs at low strain rate. Once a steady flame structure is obtained for a prescribed value of initial strain rate, transient solution of the flame is calculated after a finite amount of strain-rate perturbation is imposed on the steady flame. Depending on the initial strain rate and the amount of perturbed strain rate, transient evolution of the flame exhibits various types of flame-evolution behaviors. Basically, the dynamic behaviors can be classified into two types, namely oscillatory decaying solution and diverging solution leading to extinction.

• Korea-Australia Rheology Journal
• /
• v.17 no.3
• /
• pp.99-110
• /
• 2005

This work presents results of finite element analysis of isothermal incompressible creeping viscoelastic flows with the tensor-logarithmic formulation of the Leonov model especially for the planar geometry with singular comers in the domain. In the case of 4:1 contraction flow, for all 5 meshes we have obtained solutions over the Deborah number of 100, even though there exists slight decrease of convergence limit as the mesh becomes finer. From this analysis, singular behavior of the comer vortex has been clearly seen and proper interpolation of variables in terms of the logarithmic transformation is demonstrated. Solutions of 4:1:4 contraction/expansion flow are also presented, where there exists 2 singular comers. 5 different types spatial resolutions are also employed, in which convergent solutions are obtained over the Deborah number of 10. Although the convergence limit is rather low in comparison with the result of the contraction flow, the results presented herein seem to be the only numerical outcome available for this flow type. As the flow rate increases, the upstream vortex increases, but the downstream vortex decreases in their size. In addition, peculiar deflection of the streamlines near the exit comer has been found. When the spatial resolution is fine enough and the Deborah number is high, small lip vortex just before the exit comer has been observed. It seems to occur due to abrupt expansion of the elastic liquid through the constriction exit that accompanies sudden relaxation of elastic deformation.

• Journal of Electrochemical Science and Technology
• /
• v.6 no.4
• /
• pp.146-151
• /
• 2015

Here I report an electrochemical simulation work that compares voltammetric current and resistance of a complex electrochemical reaction over a potential scan. For this work, the finite element method is employed which are frequently used for voltammetry but rarely for impedance spectroscopy. Specifically, this method is used for simulation of a complex reaction where a heterogeneous faradaic reaction is followed by a homogeneous chemical reaction. By tracing the current and its polarization resistance, I learn that their relationship can be explained in terms of rate constants of charge transfer and chemical change. An unexpected observation is that even though the resistance is increased by the rate of the following chemical reaction, the current can be increased due to the potential shift of the resistance made by the proceeding faradaic reaction. This report envisions a possibility of the FEM-based resistance simulation to be applied to understand a complex electrochemical reaction. Until now, resistance simulations are mostly based on equivalent circuits or complete mathematical equations and have limitations to find proper models. However, this method is based on the first-principles, and is expected to be complementary to the other simulation methods.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.320-325
• /
• 2003

Acoustic-Pressure Response of diluted hydrogen-air diffusion flames is investigated numerically by adopting a fully unsteady analysis of flame structures. In the low-pressure regime, the amplification index remains low and constant at low frequencies. As acoustic frequency increases, finite-rate chemistry is enhanced through a nonlinear accumulation of heat release rate, leading to a high amplification index. Finally, the flame responses decrease at high frequency due to the response lag of the transport zone. For a medium-pressure operation and low-frequency excitation, the amplification index is low and constant. It then decreases at moderate frequencies. As frequency increases further, the amplification index increases appreciably due to an intense accumulation effect.

• Journal of Mechanical Science and Technology
• /
• v.15 no.1
• /
• pp.88-97
• /
• 2001

The present study numerically investigates the fuel-air mixing characteristics, flame structure, and pollutant emission inside a double-swirler combustor. A PSR(Perfectly Stirred Reactor) based microstructural model is employed to account for the effects of finite rate chemistry on the flame structure and NO formation. The turbulent combustion model is extended to nonadiabatic flame condition with radiation by introducing an enthalpy variable, and the radiative heat loss is calculated by a local, geometry-independent model. The effects of turbulent fluctuation are taken into account by the joint assumed PDFs. Numerical model is based on the non-orthogonal body-fitted coordinate system and the pressure/velocity coupling is handled by PISO algorithm in context with the finite volume formulation. The present PSR-based turbulent combustion model has been applied to analyze the highly intense turbulent nonpremixed flame field in the double swirler combustor. The detailed discussions were made for the flow structure, combustion effects on flow structure, flame structure, and emission characteristics in the highly intense turbulent swirling flame of the double swirler burner.

• Applied Chemistry for Engineering
• /
• v.8 no.4
• /
• pp.667-672
• /
• 1997

Thermal behavior of high capacity Nickel/Metal hybride battery in analyzed using the NISA software which is based on the three dimensional finite element method. Differential energy balance equation is used for the conduction heat transfer of the battery, while convective heat transfer equation is used for the interface between the battery and air. Heat generation rate and convective heat transfer coefficient are tested as variables to investigate thermal behavior, and the generalized equation for maximum temperature inside the battery is developed. The abrupt rise of the battery temperature due to the quick charge or discharge can be prevented from the use of metallic cooling fin. In addition, temperature augmentation of the battery is negligible when the low thermal conductive and thin insulating material is used outside of the battery case.

• Journal of Mechanical Science and Technology
• /
• v.15 no.4
• /
• pp.510-521
• /
• 2001

A nonlinear acoustic instability of subcritical liquid-oxygen droplet flames burning in gaseous hydrogen environment are investigated numerically. Emphases are focused on the effects of finite-rate kinetics by employing a detailed hydrogen-oxygen chemistry and of the phase change of liquid oxygen. Results show that if nonlinear harmonic pressure oscillations are imposed, larger flame responses occur during the period that the pressure passes its temporal minimum, at which point flames are closer to extinction condition. Consequently, the flame response function, normalized during one cycle of pressure oscillation, increases nonlinearly with the amplitude of pressure perturbation. This nonlinear response behavior can be explained as a possible mechanism to produce the threshold phenomena for acoustic instability, often observed during rocket-engine tests.