In this work, numerical parametric studies on spray combustion have been conducted. In simulation of turbulence, RNG ${\kappa}-{\varepsilon}model$ is adopted. Initial spray distribution is specified by Rosin-Rammler distribution function. Eddy break-up model is adopted as a combustion model. The parameters considered are inlet air temperature, swirl number, and SMD. With higher inlet air temperature, the axial velocities are increased and penetration of primary jet is stronger than that of lower inlet air temperature and temperature at the exit of combustor is more uniform. Combustion efficiency is improved with high inlet air temperature. The effect of swirl number on flow field is not significant. It affect only recirculation zone. So temperature at upstream of combustor is influenced. Combustion efficiency deteriorate as SMD of fuel spray increase.

Experimental studies on effects of the interaction of duplex swirl injector and the liquid temperature on the spray characteristics were conducted. Water and fuel were used as a test fluid for the experiments. The drop size distributions of the liquid spray were measured with Malvern particle sizer. The liquid temperature and distance of injectors were adopted as the operating parameter. The results show that SMD decrease as the increases of fuel temperature and pressure. The spray angle increase as fluid temperature increases. For fuel spray, SMD of impacting surface increase as the distance of injectors is lengthened.

Electrostatic pressure-swirl nozzle for practical oil burner application has been designed. The charge injection method has been used in this design, where the nozzle consists of a sharp pointed tungsten wire as a charge injector and the nozzle body grounded. The spray characteristics of the nozzles have been investigated by using an insulating liquid, i.e. kerosene without active surface agent. Breakup length of liquid decreased with an increase in applied voltage and injection pressure, while the spray angle increased with an increased in both applied voltage and injection pressure. An empirical equations have been suggested to predict the breakup length for electrostatic pressure-swirl atomizer. The experimental result was within the range of the predicted equations. The SMD decreased between the ranges of 2.8 ${\sim}$ 33% when the conventional nozzle was compared to the electrostatic with -10 kV applied to the electrode at a radial distance from 5 to 20 mm.

The autoignition and subsequent flame propagation of initially nonpremixed turbulent system have been numerically analyzed. The unsteady flamelet modeling based on the RIF (Representative Interactive Flamelet) concept has been employed to account for the influences of turbulence on these essentially transient combustion processes. In this RIF approach, the partially premixed burning, diffusive combustion and formation of pollutants(NOx, soot) can be consistently modeled by utilizing the comprehensive chemical mechanism. To treat the spatially distributed inhomogeneity of scalar dissipation rate, the multiple RIFs are employed in the framework of EPFM(Eulerian Particle Flamelet Model) approach. Computations are made for the various initial conditions of pressure, temperature, and fuel composition. The present turbulent combustion model reasonably well predicts the essential features of autoignition process in the transient gaseous fuel jets injected into high pressure and temperature environment.

An optical single cylinder diesel engine equipped with a common-rail injection system has been built to investigate diesel combustion and emission characteristics. Three optical widows (piston crown quartz for bottom view of the cylinder, upper piston quartz for allowing laser sheet and liner quartz for side view) have been placed in the optical engine to visualize spray characteristics and combustion process inside the cylinder. Before doing further research using various optical diagnostics with the optical engine, fundamental combustion experiments and flame visualization incorporating a high speed motion analyser have been carried out with a wide range of engine operating conditions.

The zone-conditioned CMC equations are derived by taking an unconditional average of the generic conservation equations multiplied by delta and Heaviside functions in terms of mixture fraction and reaction progress variable. The resulting equations are essentially in the same form as the single zone CMC equations except for separate flow fields for burned and unburned gas. The zone-conditioned two-fluid equations are applied to a stagnating turbulent premixed flame brush of Cheng and Shepherd[5l. It is shown that the flame stretch factor is of crucial importance to accurately reproduce the measured mean reaction progress variable and conditional velocities. Further work is in progress for the relationship between surface and volume averages and extension to partially premixed combustion on the basis of a triple flame structure, e. g. in a lifted turbulent diffusion flame.

The measurement of velocity and stain rate field has been conducted in opposed impinging jet combustion. When a smaller diameter (5mm) orifice of pre-chamber was used, previous studies had reported that the combustion phase showed a shift from weak turbulent combustion to moderate turbulent combustion in the modified Borghi Diagram. In the case with smaller orifice diameter (5mm), NOx emission was substantially reduced by a factor 1/2 while the combustion pressure remains at the same as that in the conventional combustion. Hence, in this study, the experiment setup using PIV technique was designed to identify the relation of the strain rate distribution and NOx reduction associated with moderate turbulent combustion.

A two-dimensional direct numerical simulation is performed to investigate the flame structure of $CH_4/N_2-Air$ counterflow nonpremixed flame interacting with a single vortex. The detailed transport properties and a modified 16-step augmented reduced mechanism based on Miller and Bowman's detailed reaction mechanism are adopted in this calculation. To quantify the strain on flame induced by a vortex, a scalar dissipation rate (SDR) is introduced. Results show that the fuel and air-side vortex cause an unsteady extinction. In this case, the flame interacting with a vortex is extinguished at much larger SDR than steady flame. It is also found that air-side vortex extinguishes a flame more rapidly than fuel-side vortex. The unsteady effect induced by flame-vortex interaction does not lead to a transient OH overshoot of the maximum steady concentration observed in experiment, while $HO_2$ radical increases more than the maximum steady concentration with increasing SDR. In addition, it is seen that NO and $NO_2$ are not sensitive to the unsteady variation of SDR.

Characteristics of lifted flames in axisymmetric laminar coflow jets have been investigated experimentally. Approximate equations for velocity and concentration with virtual origins have been proposed to analyze the behavior of flames in coflow jets. Measuring Rayleigh intensity to investigate the concentration field. proposed approximate equations were confirmed. By using the results of OH PLIF, direct photography and Rayleigh scattering measurement, it is shown that the locations of maximum intensity in direct photography coincide with the tribrachial points in axisymmetric jets and the tribrachial points travel on the stoichiometric contour. For coflow jets, the experimental results of liftoff height have been successfully correlated with nozzle exit velocity using predicted behavior from proposed approximated equations. These results substantiate the stabilization mechanism in coflow jet is based on the balance between flame propagation speed and axial flow velocity, same as for the free jets.

The dynamic behaviors of the single vortex and flame-vortex interaction in a $CH_4/Air$ diffusion flame with addition of $CO_2$ were investigated numerically. The numerical method was based on a predictor-corrector for low Mach number flow. A two-step global reaction mechanism was adopted as a combustion model. Through comparison of results by effect of $CO_2$ added either on the fuel or oxidizer side, it was found that the growth of single vortex and entrainment of surrounding fluid by $CO_2$addition on the fuel side are larger than those by $CO_2$ addition on oxidizer side. Also, when $CO_2$ is added on fuel side, flame-vortex interaction becomes more significant than on air side.

The catalytic heat exchanger was designed which employs the regenerative preheating system of combustion air. The characteristics of the catalytic heat exchanger have been experimentally studied at the various operating parameters. The results showed that the mixture velocity did not affect significantly the performance of catalytic combustor whereas the preheating temperature of combustion air affected significantly the conversion rate. The complete conversion was achieved in the catalyzed honeycomb at a preheating temperature of $370-390^{\circ}C$, a mixture velocity of 0.53 $^{\sim}$ 0.75 m/s and an equivalence ratio of 0.19 $^{\sim}$ 0.27. The heat exchange efficiency of the catalytic heat exchanger appeared to be about 75 % when the air of room temperature was used as a working fluid. The results showed that both the heat balance of the system and the mixture conditions determine its stable catalytic combustion.

Computational flow dynamics(CFD) has been frequently applied to the waste incinerators to understand the flow performance for various design and operating parameters. Though it needs many simplifications and complicated flow models, the reasonability of its results is not fully evaluated. For example, the inlet condition is calculated from an arbitrarily assumed properties of combustion gas release from the waste bed, since the combustion in the bed is difficult to be predicted. In this study, the computational modeling and calculation procedures of CFD for the grate type waste incinerator were evaluated using comparative simulations. Though the assumption method on the generation of the combustion gas directly affected the temperature and gas species concentrations, the overall flow pattern was dominated by the secondary air jets. The gaseous reaction could be included by assuming the release of the products of incomplete combusion from the bed. However, the reaction effficiency cannot not be directly evaluated from the species concentration, since it is not possible to simulate the actual co-existence of fuel rich or oxygen rich puffs over the bed. In predicting the turbulence, the higher order model, such as Reynolds stress model, gave difference shape of local recirculation zones, but similar results was acquired from the standard $k-{\varepsilon}$ model. Introducing radiation model was required for accurate temperature prediction, but it also caused heat imbalance due to the fixed temperature of the inlet, i.e. the waste bed. Thus, the computational modeling procedures on incinerators and the analysis of the predicted results should be progressed carefully. Though not validated experimentally, current simulation method is capable of comparative evaluation on the flow-related parameters such as the furnace shape and secondary air injection using identical inlet conditions. Quantitative analysis using measures of the residence time and mixing is essential to compare the flow performance efficiently.

A numerical investigation on the flame propagation and extinction in a micro combustor is described. Previous measurements of $H_2-air$ flame propagation in a submilimeter scale combustor exhibited significance of wall effects on burning velocity and extinction. The heat transfer to wall becomes important not only in the cooling of burnt gases but also during the flame ropagation, which has be by and large ignored in macro scale combustor calculations. In order to take the heat loss into account the combustion calculation, we developed a numerical code with a heat transfer model that was determined empirically from measured data. PISO algorithm was used for differencing of conservation equations. $H_2-air$ reaction was modeled with 10 species - 16 steps. Comparison with measured data showed good agreement in flame propagation speed. Also the pressure decrease after flame extinction was accurately predicted by the model. A further study is desirable for a better quenching model that can predict the quenching location.

The concentration gradient effects on triple flame have been studied experimentally using a slot burner in order to stabilize stationary triple flame in coflowing stream. By means of contraction we generate the coflowing stream with uniform velocity and linear concentration gradient at the outlet of the slot. In this paper we investigated the response of the triple flame. to the concentration gradient, like the stability, the liftoff height, and the structure of triple flame. Flow velocity is measured with Laser Doppler Velocimetry. As the concentration gradient increases. the flame propagation velocity in immediately upstream triple point increases until the liftoff height of triple flame becomes minimum, and then decreases.

Trichlorethane (TCE) is known as one of major carcinogens. TCE is difficult to be incinerated environmentally friendly, so pyrolysis is suggested for TCE treatment. In this study, we examined effects of cylinder-type pyrolysis reactor design parameters like existence of baffle inside reactor and reactor operating condition like heating reactor wall temperature and residence time numerically using CFX 4.3, a commercial computational fluid dynamic program.

The purpose of this study is to verify that the modified Lagrangian model can predict temperature, flow and scalar fields in the high temperature recirculation region of swirling confined diffusion flame. In the meantime numerical results from EBU and Equilibrium PDF models as well as experimental results are compared with those from the modified Lagrangian model. Adaption of three different turbulent models were accompanied with this procedure. Look-up table of the ignition characteristic time scale which is one of important factors of the Lagrangian model was referred to the 11-step reduced mechanism. Eventually, results with the Lagrangian model show a good accordance with experimental results, which shows the validity of this model. Results from Chen's model differ from those of the others. Numerical results of ${\widetilde{k}$ show significant deviation from experimental results for three models.

A numerical study was conducted to investigate the combustion phenomena of normal start and unstart processes based on ISL's RAMAC 30 experiments with different diluent amounts and fill pressures in a ram accelerator. The initial projectile launching speed was 1.8 km/s which corresponded to the superdetonative speed of the stoichiometric $H_2/O_2$ mixture diluted with 5 $CO_2$ or 4 $CO_2$. Experiments with same condition except for projectile surface material demonstrated that ignition was successful with an aluminum projectile, but no combustion was observed in case of a steel projectile. In this study, it was found that neither shock nor viscous heating was sufficient to ignite the mixture at a low speed of 1.8 km/s, as was found in the experiments using a steel projectile. However, we could succeed in igniting the mixtures by imposing a minimal amount of additional heat to the combustor section and simulate the normal start and unstart processes found in the experiments with an aluminum projectile. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations coupled with a Baldwin-Lomax turbulence model and detailed chemistry reaction equations of $H_2/O_2/CO_2$ suitable for high-pressure gaseous combustion were considered. The governing equations were discretized by a high order accurate upwind scheme and solved in a fully coupled manner with a fully implicit, time accurate integration method. The numerical results matched almost exactly to the experimental results. As a result, it was found that the normal start and unstart processes depended on the strength of gas mixture, development of shock-induced combustion wave stabilized by the first separation bubble, and its size and location.

The lean premixed technique has been proven very efficient in reducing NOx emissions from gas turbine combustors. However combustion instability is susceptible to occur in lean premixed combustor. So laboratory-scale dump combustor was used to understanding the underlying mechanisms causing combustion instabilities. In this study, tests were conducted at atmospheric pressure and inlet air was up to $360^{\circ}C$ with natural gas. The observed instability was a longitudinal mode with a frequency of ${\sim}341.8Hz$. At selected unstable conditions, phase-resolved OH chemiluminescence images were captured to investigate flame structure with various equivalence ratio. Combustion instability was observed to occur at higher value of equivalence ratio(>0.69). This study was performed to investigate the effects of equivalence ratio and fuel split measuring NOx and acoustic wave. The results reveal the effect of fuel-air unmixedness on lean premixed combustor.

The tulip-inversion of flames in half-open tubes was investigated experimentally. Experiments was carried out in tubes with various shapes. The image of a flame propagation were pictured by HICCD(High speed intensified CCD) and the dynamic pressure of tubes was measured by a piezo pressure sensor. By analyzing the images of the flame propagation, we found the time and the distance for the occurrence of tulip-inversion. Regardless of the shapes of tubes, time of tulip-inversion are similar and inversely proportional to the burning velocity. But distances have different tendency. In a straight tube, the distance of tulip-inversion increases when the burning velocity increases. But in a converging tube, the distance of tulip-inversion decreases when a burning velocity increases. And the distance of tulip-inversion in a converging tube is much smaller than the distance of tulip-inversion in a straight tube. These results are caused by the deceleration of a flame when the diameter of a hole in open-side of a tube is small. The deceleration causes little effect on the time of tulip-inversion.

Experimental investigations are performed on the stability and the structure of bluff-body stabilized hydrogen flames. The velocities of coflow air are varied from subsonic to supersonic velocity of Mach 1.8 and OH PLIF images and Schilieren images are used for analysis. Three characteristic flame modes are classified into three regimes with the variation of fuel-air velocity ratio; a jet like flame, a central-jet dominated flame and a recirculation zone flame. Stability curves are drawn to find the blowout regimes and to show that flame stability is improved by increasing the lip thickness of fuel nozzle that works as bluff-body. $Damk{\ddot{o}hler$ number is adopted in order to scale the blowout curves of each flame obtained at different sizes of the bluff-body and all blowout curves are scaled successfully regardless of its bluff-body size.

The characteristics of DeNOx conversion process by plasma/post-heating system with the simulated gas containing ethene is investigated experimentally. Without plasma treatment, $NO-NO_2$ conversion doesn't occur by $400^{\circ}C$ in a mixture of $N_2/O_2$ with a trace gas of ethene. But $NO-NO_2$ conversion occurs as temperature increases above $400^{\circ}C$. The NO can, however, be converted to $NO_2$ at lower temperatures by treating the gas mixture with non-thermal plasma. The $NO-NO_2$ conversion enhances further by passing the plasma treated gas through the post-heating furnace. Results show that 20%${\sim}50%$ more conversion of NO to $NO_2$ is observed when the temperatures of the post-heating furnace are maintained at $300^{\circ}C$ or $400^{\circ}C$. The additional $NO-NO_2$ conversion by post-heating is due to the reaction of ethene with the byproducts or radicals generated from the plasma reaction.

The main objective of this study is to investigate the correlation of chemical structure and cetane number of reformed diesel fuels by ultrasonic irradiation. In order to analyze the effect of the chemical structure and the cetane number of reformed diesel fuels by ultrasonic irradiation, $^1H-NMR$ was used. From the study, following conclusive remarks can be made. 1) BI(=Branch Index), aromatics percentages, and $H_{\alpha}(={\alpha}-methyl$ functional group) of the reformed diesel fuels by ultrasonic irradiation decreased more than those of the conventional diesel fuel. 2) All the cetane numbers which were calculated from carbon type structure and hydrogen type distribution of the reformed diesel fuels increased more than those of the conventional diesel fuel. 3) Using predicated equation of cetane number caculated from carbon type structure is more reasonable than that caculated from hydrogen type distribution 4) BI, aromatics percentages, and $H_{\alpha}$ on both of conventional fuel and reformed diesel fuels by ultrasonic irradiation are inversely proportional to cetane number on these fuels.

Numerical simulations of the condition in the iron ore sintering bed are performed for various parameters. The sintering bed is modelled as an unsteady one-dimensional progress of solid material, containing cokes and iron ore. Bed temperature, solid mass and gas species distributions are predicted for various parameters of moisture contents, cokes contents and air suction rates, along with the various particle diameters of the solid for sensitivity analysis. Calculation results show that influences of these parameters on the bed condition should be carefully evaluated for achievement of the self-sustaining combustion without the high temperature section, which can cause the excessive melting in the bed. It suggests that the model should be extended to consider the bed structural change and multiple solid phase, which can treat the inerts and fuel particles separately.

Interaction between flames and non-thermal plasmas of DBD type has been experimentally investigated. Vigorous streamers were observed under flame conditions because of the increase of reduced field (electrical) at high temperature as well as the seeding of free electrons and ions generated inside the flame. Flame lengths were significantly shortened as the applied voltage increased on account of intense mixing by ionic winds and soot-induced flows. Flame luminosities severely decreased under plasma conditions, which means the reduction of soot, since the residence time was reduced because of the flame shortening. Temperature and major species concentrations measured by FTIR were not changed despite the plasma generation. which shows overall chemistries were not affected by non-thermal plasmas.

The influence of N2 addition on soot formation, flame temperature and NOx emissions is investigated experimentally with methane fuel co-flow diffusion flames. The motivation of the present investigation is the differences in NOx reduction reported between fuel-side and oxidizer-side introduction of N2. To determine the influence of dilution alone, fuel was diluted with nitrogen while keeping the adiabatic flame temperature fixed by changing the temperature of the reactants. And to see the thermal effect only, air was supplied at different temperature without N2 addition. N2 addition into fuel side suppressed the soot formation than the case of oxidizer-side, while flame temperature enhanced the soot formation almost linearly. These results reveals the relative influences of the thermal, concentration effects of N2 additives on soot formation In accordance with experimental study, numerical simulation using CHEMKIN 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. Emission test revealed that NOx emissions were affected by not only flame temperature but also N2 addition.

In this paper, the effects of water contents on combustion characteristics of paper sludge and capacity of fluidized bed combustor(FBC) were investigated using 0.26m diameter, 1.75m height pilot-plant scale combustor. Combustion tests of paper sludge containing water contents between 40wt% and 50wt% were performed. The temperature and emission variation, the pressure inside combustor were measured to monitor the fluidization quality. The experimental results showed that 30kg/hr feeding rate of sludge containing water up to 45wt% was preferable for this system. Sludge loading rate, heat release rate were calculated from experimental data as major parameters showing FBC capacity. Comparsion with sludge loading rate from other source was also performed

Combustion experiments are carried out to obtain the fundamental data for application of a internally circulating fluidized bed combustor to the combustion of paper sludge wastes. Experimental parameters are identified as secondary air ratio, sludge weight and water contents. The secondary air ratio was varied from 0 to 20% and water content was 14%, 25%, 35%, 45%, 62% and sludge weight was 30g, 60g, 90g. As a result, carbon conversion ratio was higher than injection of primary air. Emission of CO and NOx reduced with an increase of secondary air injection.

Waste fuels, which originate from different sources, have unique combustion characteristics. The characteristics should be considered in applying FBC(fluidized bed combustor) technology to those fuels. The effects of fuel properties and operating conditions on FBC reactivity were investigated by means of carbon based parameter called mean carbon conversion time, rate of carbon conversion, fraction of carbon conversion and carbon recovery. And the basic physical and chemical mechanisms taking place in a fluidized bed were summarized. Major parameters in designing and operating FBC were evaluated in terms of the fuel properties and the combustion environment.

가압유동층연소로(bed dia. 0.17m, freeboard dia. 0.25, total height 5m)에서 국내무연탄을 연소시켜 이에 대한 탈황특성을 고찰하였다. 실험은 압력($2{\sim}6atm$), 요동층온도(($850{\sim}950^{\circ}C$), 과잉공기(10, 20, 30,%)등의 조건과, 탈황을 위한 Ca/S몰비(($0.8{\sim}4.8$)가 탈황특성에 미치는 영향을 고찰하였다. 결과적으로, 본 연구의 실험범위에서 연소효율은 $80{\sim}99%$를 보였고, 연소온도, 압력 그리고 과잉공기가 증가 할수록 증가하였다. 배가스중의 $SO_2$배출농도는 압력, Ca/S몰비가 증가함에 따라 감소하였다. 탈황율은 상압에서 층(bed)온도의 증가에 따라 감소하였다. 운전압력이 증가함에 따라 탈황율의 감소폭이 둔화되었다. 과잉공기가 증가함에 따라 탈황율은 증가하는 경향을 나타냈다. 각각의 운전압력에서 과잉공기의 증가에 따라 약 10%의 증가폭을 보였다. 국내무연탄을 연소하는 경우 운전압력 4atm일 때 Ca/S몰비는 4이상 주입하여야 하고, 6atm일 때 Ca/S 몰비가 2.5이상주입하여야 150ppm이하를 보여 배출규제치를 만족하는 것을 보였다.