• Journal of Aerospace System Engineering
• /
• v.1 no.2
• /
• pp.37-42
• /
• 2007

This paper proposes a numerical modeling approach to simulate the hybrid combustion phenomena. From the physical understandings of hybrid combustion, the computational domain was separated into three regions: the solid fuel, gas phase reactive flow, and the interface between solid and fluid. Moreover, for the accurate calculation, computational grids for these regions was generated at every time step considering the instantaneous moving interface which are governed by the balance equations using thermal pyrolysis. In the domain of reactive flow, by virtue of diffusion flame structure, turbulent combustion modeling was introduced using either mixture fraction approach or mean reaction rate approach.

• Bulletin of the Korean Chemical Society
• /
• v.13 no.4
• /
• pp.398-404
• /
• 1992

We have investigated the kinetics of diffusion-influenced bimolecular reactions in which one reactant has an internal mode, called the gating mode, that activates or deactivates its reactivity intermittently. The rate law and an expression for the time-dependent rate coefficient have been obtained from the general formalism based on the hierarchy of kinetic equations involving reactant distribution functions. The analytic expression obtained for the steady-state reaction rate constant coincides with the one obtained by Szabo et al., who derived the expression by employing the conventional concentration-gradient approach. For the time-dependent reaction rate coefficient, we obtained for the first time an exact analytic expression in the Laplace domain which was then inverted numerically to give the time-domain results.

• Journal of Industrial Technology
• /
• v.17
• /
• pp.241-247
• /
• 1997

A numerical model has been proposed for a diffusion flame reactor to manufacture ultrafine $TiO_2$ powders. The model equations such as mass balance equation, the 0th, 1st, and 2nd moment equations of aerosols were considered. The phenomena such as $TiCl_4$ reaction rate, $TiO_2$ nucleation rate and the coagulation of $TiO_2$ powders were included in the aerosol dynamic equation. It is found that the $TiO_2$ particle concentration becomes higher, as the inlet $TiCl_4$ concentration and the total gas flow rate increase, and also as the flame temperature decreases. The $TiO_2$ particle size increases, as the flame temperature and the inlet $TiCl_4$ concentration increase and the total gas flow rate decreases.

• Journal of applied mathematics & informatics
• /
• v.28 no.5_6
• /
• pp.1035-1054
• /
• 2010

We consider a mixed type singularly perturbed one dimensional elliptic problem with discontinuous source term. The domain under consideration is partitioned into two subdomains. A convection-diffusion and a reaction-diffusion type equations are posed on the first and second subdomains respectively. Two hybrid difference schemes on Shishkin mesh are constructed and we prove that the schemes are almost second order convergence in the maximum norm independent of the diffusion parameter. Error bounds for the numerical solution and its numerical derivative are established. Numerical results are presented which support the theoretical results.

• 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 the Korean Society of Combustion
• /
• v.21 no.4
• /
• pp.23-29
• /
• 2016

Reactivity of gasification defined by bouardard reaction is critical parameter in efficiency of the gasifier. In this study, char reactivity of the gasification was derived from the experiments using the intrinsic reaction kinetics model. Pressurized wire mesh heating reactor (PWMR) can produce high temperature and high pressure conditions up to 50 atm and 1750 K, respectively and PWMR was designed to evaluate the intrinsic reaction kinetics of $CO_2$ gasification. In this study, Kideco and KCH (sub-bituminous Indonesian coal) were pulverized and converted into char. Experiments used the PWMR were conducted and the conditions of the temperature and pressure were 1373~1673 K, 1~40 atm. To distinguish the pressure effect from high pressurized condition, internal and external effectiveness factors were considered. Finally, the intrinsic kinetics of the Kideco and KCH coal char were derived from $n^{th}$ order reaction rate equations.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.26 no.1
• /
• pp.20-27
• /
• 2022

Present modeling study of nanoenergetics focuses on the numerical simulation of reaction wave propagation in normal direction across nanoscale multilayers of aluminum and nickel combination. The governing equations for atomic and thermal diffusion are employed in one-dimensional semi-infinitely alternating Al/Ni multilayered structures and the numerical results show the established patterns of quasi-steady intermetallic reaction waves. Also, the reaction wave speed is confirmed to be highly independent of reaction wave directions in such nanoenergetic structures.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.16 no.3
• /
• pp.133-140
• /
• 2007

A 2-D model of fluid flow, mass transport and electrochemistry is analysed to examine the effect of current density at the current collector depending on active layer thickness of catlyst in polymer elecrolyte fuel cells. The finite element method is used to solve the continuity, potential and Maxwell-Stefan equations in the flow channel and gas diffusion electrode regions. For the material behavior of electrode reactions in the active catalyst layers, the agglomerate model is implemented to solve the diffusion-reaction problem. The calculated model results are described and compared with the different thickness of active catalyst layers. The significance of the results is discussed in the viewpoint of the current collecting capabilities as well as mass transportation phenomena, which is inferred that the mass transport of reactants dictates the efficiency of the electrode in the present analysis.

• Journal of Korean Society for Atmospheric Environment
• /
• v.12 no.5
• /
• pp.577-585
• /
• 1996

A screening study was performed in order to resolve the flow, combustion and emission characteristics of the gas furmace with co-axial diffusion flane burner. A control-valume based finite-difference method with the power-law scheme was employed for discretization. Numerical procedure for the differential equation was used by SIMPLEST to enclosute rapid converge. A k-.varepsilon. model was incorporated for the closure of turbulence. The mass fraction and mixture fraction were calculated by cinserved scalar method. An equilibrium analysis was employed to determine the concentration of radicals in the product stream and conserbation equations were them solved for N amd NO by Zelovich reaction scheme. The method was exercised in a simple one-dimensional case first, to determine the effects of air ratio, temperature and residence time on NO formation and applied to a furnace with co-axial diffusion flame burner.

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

This paper investigates the linear stability of reacting mixing layers with special emphasis on the existence of multiple unstable modes. The governing equations for laminar flows are from two-dimensional compressible boundary-layer equations. The chemistry is a finite rate single step irreversible reaction with Arrhenius kinetics. For the incompressible reacintg mixing layer with variable density. A necessary condition for instability has been derived. The condition requires that the angular momentum, not the vorticity, to have a maximum in the flow domain. New inflectional modes of instability are found to exist in the outer part of the mixing layer. For the compressible reacting mixing layer, supersonic unstable modes may exist in the abscence of a generalized inflection point. The outer modes at high Mach numbers in the reacting mixing layer are continuations of the inflectional modes of low Mach number flows. However, the generalized inflection point is less important at supersonic flows.