• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.9
• /
• pp.1185-1191
• /
• 1999

Radiation-induced oscillatory instability in diffusion flames is numerically investigated with nonlinear dynamics considered. As the simplest flame model, a diffusion flame established in the stagnant mixing layer is employed with optically thin gas-phase radiation and unity Lewis numbers for all species. Attention is focused on the radiation-induced extinction regime, which occurs at large $Damk\ddot{o}hler$ number. Once the steady flame structure is obtained for a prescribed value of the initial $Damk\ddot{o}hler$ number, transient solution of the flame is calculated after a finite amount of the $Damk\ddot{o}hler$-number perturbation is imposed on the steady flame. Transient evolution of the flame exhibits three types of flame-evolution behaviors, namely decaying oscillatory solution, diverging solution to extinction and stable limit-cycle solution. A dynamic extinction boundary is identified for laminar flamelet library.

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

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.2
• /
• pp.129-138
• /
• 2013

A spatio-temporal models as systems of ODE which describe two-species Beddington - DeAngelis type predator-prey system living in a habitat of two identical patches linked by migration is investigated. It is assumed in the model that the per capita migration rate of each species is influenced not only by its own but also by the other one's density, i.e. there is cross diffusion present. We show that a standard (self-diffusion) system may be either stable or unstable, a cross-diffusion response can stabilize an unstable standard system and destabilize a stable standard system. For the diffusively stable model, numerical studies show that at a critical value of the bifurcation parameter the system undergoes a Turing bifurcation and the cross migration response is an important factor that should not be ignored when pattern emerges.

• Journal of the Korean Society of Combustion
• /
• v.17 no.3
• /
• pp.9-16
• /
• 2012

Nonlinear dynamics of pulsating instability in radiating counterflow diffusion flames is numerically investigated by imposing Damk$\ddot{o}$hler number perturbation. Stable limit-cycle solutions occur in small ranges of Damk$\ddot{o}$hler numbers past bifurcation point of instability. Period doubling cascade and chaotic behaviors appear just before dynamic extinction occurs. Nonlinear dynamics is also studied when large disturbances are imposed to flames. For weak steady flames, the dynamic extinction range shrinks as the magnitudes of disturbances are increased. However, strong steady flames can overcome relatively large disturbances, thereby the dynamic extinction range extending. Stable limit-cycle behaviors reappears prior to dynamic extinction when the steady flames are strong enough.

• Journal of The Korean Astronomical Society
• /
• v.37 no.5
• /
• pp.601-603
• /
• 2004

We present the results of the linear analysis for the Parker-Jeans instability in the magnetized gas disks including the effect of cosmic-ray diffusion along the magnetic field lines. We adopted an uni-formly rotating two temperature layered disk with a horizontal magnetic fields and solved the perturbed equations numerically. Fragmentation of gases takes place and filamentary structures are formed by the growth of the instability. Nagai et al. (1998) showed that the direction of filaments being formed by the Parker-Jeans instability depends on the strength of pressure outside the unperturbed gas disk. We found that at some range of external pressures the direction of filaments is also governed by the value of the diffusion coefficient of CR along the magnetic field lines k.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.8
• /
• pp.857-864
• /
• 2012

A linear stability analysis of a diffusion flame with radiation heat loss is performed to identify linearly unstable conditions for the Damk$\ddot{o}$hler number and radiation intensity. We adopt a counterflow diffusion flame with unity Lewis number as a model. Near the kinetic limit extinction regime, the growth rates of disturbances always have real eigenvalues, and a neutral stability condition perfectly falls into the quasi-steady extinction. However, near the radiative limit extinction regime, the eigenvalues are complex, which implies pulsating instability. A stable limit cycle occurs when the temperatures of the pulsating flame exceed the maximum temperature of the steady-state flame with real positive eigenvalues. If the instantaneous temperature of the pulsating flame is below the maximum temperature, the flame cannot recover and goes to extinction. The neutral stability curve of the radiation-induced instability is plotted over a broad range of radiation intensities.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.26 no.7
• /
• pp.915-921
• /
• 2020

Thermoacoustic instability caused by air conditioning in a combustion chamber has emerged as a problem that must be solved to establish a stable combustion system. Thermoacoustic instability is largely divided into primary and secondary acoustic instability. In this study, an experimental study of the effects of heat losses was conducted to investigate the mechanism of secondary acoustic instability. To generate the secondary acoustic instability, a quarter-wavelength resonator with one open end and one closed end was used, and the inside of the resonator was filled with premixed gases. Subsequently, secondary acoustic instability with downward-propagating flames could be realized via thermal expansion on the burnt side. To control heat losses qualitatively, an additional co-axial tube was installed in the resonator with air or nitrogen supply. Therefore, additional diffusion flames can be formed at the top of the resonator depending on the injection of the oxidizer into the co-axial tube when rich premixed flames are used. Consequently, secondary acoustic instability could not be achieved by increasing heat losses to the ambient when the additional diffusion flame was not formed, and the opposite result was obtained with the additional diffusion flame.

• The Pure and Applied Mathematics
• /
• v.15 no.3
• /
• pp.329-342
• /
• 2008

There have been many experimental and observational evidences which indicate the predator response to prey density needs not always monotone increasing as in the classical predator-prey models in population dynamics. Holling type functional response depicts situations in which sufficiently large number of the prey species increases their ability to defend or disguise themselves from the predator. In this paper we investigated the stability and instability property for a Holling type predator-prey system of a generalized form. Hopf type bifurcation properties of the non-diffusive system and the diffusion effects on instability and bifurcation values are studied.

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

Computational experiments on fundamental un stretched laminar burning velocities and flame response to stretch (represented by the Markstein number) of hydrogen-air flames at high temperatures and pressures were conducted in order to understand the dynamics of the flames including hydrogen as an attractive energy carrier in conditions encountered in practical applications such as internal combustion engines. Outwardly propagating spherical premixed flames were considered for a fuel-equivalence ratio of 0.6, pressures of 5 to 50 atm, and temperatures of 298 to 1000 K. For these conditions, ratios of unstretched-to-stretched laminar burning velocities varied linearly with flame stretch (represented by the Karlovitz number), similar to the flames at normal temperature and normal to moderately elevated pressures, implying that the "local conditions" hypothesis can be extended to the practical conditions. Increasing temperatures tended to reduce tendencies toward preferential-diffusion instability behavior (increasing the Markstein number) whereas increasing pressures tended to increase tendencies toward preferential-diffusion instability behavior (decreasing the Markstein number).

• Bulletin of the Society of Naval Architects of Korea
• /
• v.19 no.4
• /
• pp.19-29
• /
• 1982

Galerkin's finite element method is applied to a two-dimensional heat convection-diffusion problem arising in the hydrodynamic lubrication of thrust bearings used in naval vessels. A parabolized thermal energy equation for the lubricant, and thermal diffusion equations for both bearing pad and the collar are treated together, with proper juncture conditions on the interface boundaries. it has been known that a numerical instability arises when the classical Galerkin's method, which is equivalent to a centered difference approximation, is applied to a parabolic-type partial differential equation. Probably the simplest remedy for this instability is to use a one-sided finite difference formula for the first derivative term in the finite difference method. However, in the present coupled heat convection-diffusion problem in which the governing equation is parabolized in a subdomain(Lubricant), uniformly stable numerical solutions for a wide range of the Peclet number are obtained in the numerical test based on Galerkin's classical finite element method. In the present numerical convergence errors in several error norms are presented in the first model problem. Additional numerical results for a more realistic bearing lubrication problem are presented for a second numerical model.