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

• Atmosphere
• /
• v.27 no.4
• /
• pp.455-466
• /
• 2017

In Korea, property and human damages occur annually due to heavy precipitation during the summer. On August 8, 2015, heavy rainfall occurred in the Seoul metropolitan area due to an outflow boundary, and $77mmhr^{-1}$ rainfall was recorded in Gwangju, Gyeonggi Province. In this study, the simulation of the WRF numerical model is performed to understand the cause and characteristics of heavy rainfall using the Conditional Instability of the Second Kind (CISK), potential vorticity (PV), frontogenesis function, and convective available potential energy (CAPE) analyses, etc. Convective cells initiated over the Shandong Peninsula and located on the downwind side of an upper level trough. Large amounts of water vapor were supplied to the Shandong Peninsula along the southwestern edge of a high pressure system, and from the remnants of typhoon Soudelor. The mesoscale convective system (MCS) developed through CISK process and moved over to the Yellow Sea. The outflow boundary from the MCS progressed east and pushed cold pool eastward. The warm and humid air over the Korean Peninsula further enhanced convective development. As a result, a new MCS developed rapidly over land. Because of the latent heat release due to convection and precipitation, strong potential vorticity was generated in the lower atmosphere. The rapid development of MCS and the heavy rainfall occurred in an area where the CAPE value was greater than $1300Jkg^{-1}$ and the fronto-genesis function value of 1.5 or greater coincided. The analysis result shows that the MCS driven by an outflow boundary can be identified using CISK process.

• Atmosphere
• /
• v.26 no.1
• /
• pp.141-158
• /
• 2016

The objective of this study is to examine the impact of urban canopy and the horizontal resolution on simulated meteorological variables such as 10-m wind speed, 2-m temperature and precipitation using WRF model for a local, convective rainfall case. We performed four sensitivity tests by varying the use of urban canopy model (UCM) and the horizontal resolution, then compared the model results with observations of AWS network. The focus of our study is over the Seoul metropolitan area for a convective rainfall that occurred on 16 August 16 2015. The analysis shows that mean diurnal variation of temperature is better simulated by the model runs with UCM before the convective rainfall. However, after rainfall, model shows significant difference in air temperature among sensitivity tests depending on the simulated rainfall amount. The rainfall amount is significantly underestimated in 0.5 km resolution model run compared to 1.5 km resolution, particularly over the urban areas. This is due to earlier occurrence of light rainfall in 0.5 km resolution model. Earlier light rainfall in the afternoon eliminates convective instability significantly, which prevents occurrence of rainfall later in the evening. The use of UCM results in a higher maximum rainfall in the domain, which is due to higher temperature in model runs with urban canopy. Earlier occurrence of rainfall in 0.5 km resolution model is related to rapid growth of PBL. Enhanced mixing and higher temperature result in rapid growth of PBL, which provides more favorable conditions for convection in the 0.5 km resolution run with urban canopy. All sensitivity tests show dry bias, which also contributes to the occurrence of light precipitation throughout the simulation period.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.17 no.3
• /
• pp.256-261
• /
• 2005

The objective of the present study is to study the Soret effect of both nanoparticles and solute on the convective instabilities in binary nanofluids. A new stability criterion is obtained based on the linear stability theory. The results show that the Soret effect of solute(${\psi}_{bf}$) makes the binary nanofluids unstable significantly and the convective motion in a binary nanofluid sets in easily as the ratio of Soret coefficient of nanofluid to that of binary basefluid ${\delta}_4$ increases for ${\delta}_4$ > -1. It is also found that as an increase of the volume fraction of nanoparticles, nanofluid becomes stable but at a separation ratio of ${\psi}=-0.3$ the state of fluid changes from stable to unstable.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.285-287
• /
• 2001

We have studied the nonlinear evolution of a magnetized disk of isothermal gas, which is sustained by its self-gravity. Our objective is to investigate how the Jeans, Parker, and convective instabilities compete with each other in structuring/de-structuring large scale condensations in such disk. The Poisson equation for the self-gravity has been solved with a fourth-order accurate Fourier method along with the Green function, and the MHD part has been handled by an isothermal TVD code. When large wavelength perturbations are applied, the combined action of the Jeans and Parker instabilities suppresses the development of the convection and forms a dense core of prolate shape in the mid-plane. Peripheral structures around it are filamentary. The low density filaments connect the dense core to the diffuse upper region. On the other hand, when small wavelength perturbations are applied, the disk develops into an equilibrium state which is reminiscent of the Mouschovias's 2-D non-linear equilibrium of the classical Parker instability under an externally given gravity.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.4
• /
• pp.429-438
• /
• 2017

Cranke-Nicolson scheme in native OpenFOAM source libraries was not able to provide 2nd order temporal accuracy of velocity and pressure since the volume flux of convective nonlinear terms was 1st accurate in time. In the present study the simplest way of getting the volume flux with 2nd order accuracy was proposed by using old fluxes. A possible numerical instability originated from an explicit estimation of volume fluxes could be handled by introducing a weighting factor which was determined by observing the ratio of the finally corrected volume flux to the intermediate volume flux at the previous step. The new calculation of volume fluxes was able to provide temporally accurate velocity and pressure with 2nd order. The improvement of temporal accuracy was validated by performing numerical simulations of 2D Taylor-Green vortex of which an exact solution was known and 2D vortex shedding from a circular cylinder.

• Korean Chemical Engineering Research
• /
• v.57 no.1
• /
• pp.142-147
• /
• 2019

The onset of buoyancy-driven convection in an initially isothermal and quiescent horizontal fluid layer was analyzed theoretically. It is well-known that at the critical Rayleigh number $Ra_c=669$ convective motion sets in with a constant-heat-flux cooling through the upper boundary. Here, based on the momentary instability concept, the dimensionless critical time ${\tau}_m$ to mark the onset of convective motion for Ra > 669 was analyzed theoretically. The energy method under the momentary stability concept was used to find the critical conditions as a function of the Rayleigh number Ra and the Prandtl number Pr. The predicted critical conditions were compared with the previous theoretical and experimental results. The momentary stability criterion gives more reasonable wavenumber than the conventional energy method.

• Journal of Mechanical Science and Technology
• /
• v.17 no.9
• /
• pp.1388-1396
• /
• 2003

The hydrodynamic instability of the three-dimensional boundary layer on a rotating disk introduces a periodic modulation of the mean flow in the form of stationary cross flow vortices. The instability labeled Type II by Faller occurs first at lower Reynolds number than that of well known Type I instability. Detailed numerical values of the amplification rates, neutral curves and other characteristics of the two instabilities have been calculated over a wide range of parameters. Presented are the neutral stability results concerning the two instability modes by solving the appropriate linear stability equations reformulated not only by considering whole convective terms but also by correcting some errors in the previous stability equations. The present stability results agree with the previously known ones within reasonable limit. Consequently, the flow is found to be always stable for a disturbance whose dimensionless wave number is greater than 0.75. Some spatial amplification contours have been computed for the stationary disturbance wave, whose azimuth angle $\varepsilon$= 11.29$^{\circ}$ to 15$^{\circ}$ and for the moving disturbance wave, whose azimuth angle $\varepsilon$ = 12.5$^{\circ}$ to 15$^{\circ}$. Also, some temporal amplification contours have been computed for the stationary disturbance wave, whose azimuth angle $\varepsilon$= 11.29$^{\circ}$ to 15$^{\circ}$ and for the moving disturbance wave, whose azimuth angle $\varepsilon$= 12$^{\circ}$ to 15$^{\circ}$. The flow instability was observed by using a white titanium tetrachloride gas over rotating disk system. When the numerical results are compared to the present experimental data, the numerical results agree quantitatively, indicating the existence of the selective frequency mechanism.

• Proceedings of the SAREK Conference
• /
• 2004.06a
• /
• pp.175-175
• /
• 2004