To solve the problems of energy and environment conservation issued recently, mainly in mechanical engineering point of view, R&D's on the thermal, fluid and environmental engineering technology have been carried out by two R&D departments in the Korea Institute of Machinery & Materials (KIMM). Now there are 65 researchers in the two. The representative projects in the field of thermal and fluid engineering are development of an inactive gas generator and development of a cryogenic cooler for electronic sensors. Pyrolysis and melting of wastes, gas treatment using nonthermal plasma, and desalination are important technology to be developed in environmental R&D areas. To reduce the emission from the existing diesel engines for buses, an LPG direct injection type of bus engine is being developed supported by LPG supply companies. These several R&D projects which have been carried out in KIMM are introduced briefly.

The liquid fuel film on the cylinder liner is believed to be a major source of engine-out hydrocarbon emissions in SI engines, especially during cold start and warm-up period. Quantifying the liquid fuel film on the cylinder liner is essential to understand the engine-out hydrocarbon emissions formation in SI engines. In this research, two-dimensional visualization was carried out to quantify liquid fuel film on the quartz liner in an SI engine test rig. The visualization was based on laser-induced fluorescence and total reflection. Using a quartz liner and a special lens, only the liquid fuel on the liner was visualized. The calibration technique was developed to quantify the fluorescence signal with the thickness gage and the calibration device. The fluorescence intensity increases linearly with increase in the fuel film thickness on the quartz liner. Using this technique, the distribution of the fuel film thickness on the cylinder liner was measured quantitatively for different valve lifts and injected fuel mass in the test rig.

As a part of micro engine development feasibility estimation was done through fabrication and test of down scaled combustor and MEMS fabricated spark electrode. In an experimental observation of the down scaled combustion phenomena where flame propagation was observed by optical method and pressure change in combustor which gives the information about the reaction generated thermal energy was recorded and analyzed. Optimal combustor scale was derived to be about 2mm considering increased heat loss effect and thermal energy generation capability. Through the fabrication and discharge test of MEMS electrode effects of electrode width and gap was investigated. Electrode was fabricated by thick PR mold and electroplating. From the result discharge voltage characteristic in sub millimeter scale electrode having thickness of $40{\mu}m$ was obtained. From the result base technology for design and fabrication of micro engine was obtained.

This paper computes the bluff-body stabilized jet and flame. This study numerically investigates the nonpremixed $C_{2}H_{4}-air$ jet for the nonreacting case and the nonpremixed $CH_{3}OH-air$ turbulent flames for the reacting case using the laminar flamelet model on modified KIVA2 code. And this study predicts $NO_{x}$ formation characteristics using Eulerian Particle Flamelet Model. In the present study, the turbulent combustion model is applied to analyze both nonreacting and reacting case. And both standard $k-{\varepsilon}$ model and modified $k-{\varepsilon}$ model are used in nonreacting case. Calculations are compared with experimental data in terms of velocity, mixture fraction, mixture fraction Root Mean Square and Temperature. The present model correctly predicts the essential features of flame structures and $NO_{x}$ formation characteristics in the bluff-body stabilized flames.

Combustor design technique is established by reverse engineering of existing combustor and applying heat & mass balance equations for the combustion process. The ratio of entrained air for each air slot is found to be almost proportional to the area ratio from the result of numerical simulation. The shape of the combustor is modified by the numerical analysis to get circumferentially uniform flow inside the combustion chamber required for the flame stability.

A conventional flame type gas combustion major portion of heat is transferred to the body by convection due to small radiant ability of the gas flame. Increasing the radiation component of heat flux in the combustion zone allows to augment the efficiency of gas utilization. Such effect can be reached by using radiative gas burner applied to metal mesh combustion. Basically the gas radiant burner consists of metallic mesh of high heat resisting steels. In terms of this regards, we have made the burner consisted of metal mesh and measured the radiative flame stability of natural gas/air mixture on the metal mesh burner. The pressure loss through the metal mesh is defined by pressure-velocity slope. The more increased the pressure-velocity slope of the metal mesh is, the wider the stable zone of radiave flame on the metal mesh burner is. And the augmentation of mixture flowrate through the metal mesh make narrow the permissible range of equivalence ratio.

OH radical concentration have been measured in a methane-air partially premixed flames using PLIF. Excitation lines were selected $Q_{1}(6)$ branch, (1,0) band. The system is consisted of Nd:YAG laser, dye laser and frequency doubler to make pump beam for OH radical. On the direct photographs, flame height increases as fuel flow rate and equivalence ratio increase. And on the PLIF images, OH radical is distributed from premixed flame front to nonpremixed flame front through the flame structure with all equivalence ratio. OH overall concentrations increase with equivalence ratio. At the stoichiometric equivalence ratio, the peak of OH radical concentration exists strongly near the inner cone. As equivalence ratio is changed to richer, OH radical distribution goes thinly and the peak is increased as longitudinal direction. As the flow goes to the downstream, OH radical concentration decreases and broadens, because OH radical reacts with another species after OH formation at the initial oxidization. This phenomenon resembles radial distribution. At the l00cc fuel flowrate, the radial peak of OH radical exists from x/R=l.0 to 1.5.

Bomb calorimeter was developed for measuring the calorific value of combustible matter such as wastes. The calorimeter consist of bomb, stirred-water type bucket, thermometer and ignition circuit. Operation and performance of the calorimeter have been tested experimentally. In the present study, calorific values of light oil, lamp oil and bunker C oil is measured using the bomb calorimeter. Mass of the sample is fixed at lg, and oxygen pressure in the bomb is used as an experimental parameter. Sample in the oxygen bomb is burned with electrically heated Ni-Cr wire of 100mm in length, and temperature of water in the bucket become increased by $5^{\circ}C$ during about 30min. Calorific value of the sample is calculated with the temperature difference of water. Combustion tests, such as the record of temperature history and the inspection of remnants, are performed at 4, 6, 8 and 10 atm of the oxygen pressure. From the test results, oxygen pressure in the bomb must be over 10atm for complete combustion.

The combustion characteristics have been investigated to develop the low $NO_{x}$ gas turbine combustor. The lean premixed combustion technology was applied to reduce the $NO_{x}$ emission. Also, the conventional combustor was designed and tested for the baseline of low $NO_{x}$ combustor performance. The test was conducted at the condition of high temperature and ambient pressure. The combustion air which has the temperature of 500K were supplied to the combustor through the air preheater. The temperature and emissions of $NO_{x}$ and CO were measured at the exit of combustor. The premixing chamber can be operated very lean condition of equivalence ratio around 0.35. The $NO_{x}$ was decreased with decreasing the equivalence ration. The CO was decreased with decreasing the equivalence ratio, but the CO was increased with decreasing the equivalence ratio below 0.45. But, at the very lean condition of equivalence ratio below 0.35 both NOx and CO were increased because of the flame unstability. The $NO_{x}$ was decreased slightly and CO was increased with increasing inlet air flowrate. This results can be used to determine the size of combustor. The low $NO_{x}$ combustor has lower values of $NO_{x}$ and CO compared with conventional one. Consequently the performance of combustor shows the possibility of the application to the gas turbine system.

It was numerically studied that characteristics of fluid flow and heat transfer in a tube with disk and annular baffle for heat exchanger of condensing boiler. Using a finite volume technique and CFD code, STAR-CD, the governing equations were solved and the temperature and flow fields were investigated. The interval between tube and annular baffle, height and diameter of baffle were selected as important design parameters, and the effects of these parameters on heat transfer and fluid flow were studied. As a result, in the case of with interval, the pressure was decreased but heat transfer was increased. Also heat transfer was slowly increased as the size of disk and annular baffle were increased and the distance between baffles were decreased.

The effect of spray characteristics on desulfurization yield was evaluated by performing experiments with pilot spray drying sorber (SDS). Among the variables of operating conditions, the slurry-spraying conditions were chosen as major parameters; Stoichiometric ratio and Sauter mean diameter of slurry droplet were varied for the different gas temperatures and $SO_{2}$ concentrations in the inlet gas flow. From the experimental results, we proposed semi-empirical models of desulfurization yield for both Stoichiometric ratio and Sauter mean diameter of droplets. The optimal condition of spray can be determined based on these results, which might be applied to the design or scale-up of SDS systems.

This work presents experimental results of filter media test by using particles from $0.02\;to\;0.6{\mu}m$ in diameter and by applying different charging states. In order to investigate the electret filter performance, the collection efficiency and the pressure drop of filter were measured. The face velocities of test filters varied from 2.4 to 20.4 cm/s. Another experiment setup for the cartridge cabin air filter was used to get an collection efficiency in submicron region. After charging level of electret filter severely decreased, the change of collection efficiency was verified. Experimental results show that the reliability of electret filters can be poor in some conditions.

Ultrafine particles have been widely used in many high technology industrial areas. The spherical nonagglomerated and uniform nanometer-size $SiO_{2}$ particles are synthesized by the injection of TEOS vapor, ions and reaction gas in furnace. Ions were generated by Corona discharge electrode and these ions charge $SiO_{2}$ particles. As a result, spherical, nonagglomerated and ultrafine particles are generated in various conditions, it's morphology, charging portion and size distribution are examined by using TEM, ESP and SMPS. As the applied voltage of electrode changes from 0 kV to 5.0 kV, it is observed that the mean diameter of $SiO_{2}$ particle decreases from 94 nm to 42 nm.

The effects of the electrohydrodynamic (EHD) flow and turbulent diffusion on the collection efficiency of a model ESP composed of the plates with a cavity were studied through numerical computation. The electric field and ion space charge density were calculated by the Poisson equation of the electrical potential and the current continuity equation. The EHD flow field was solved by the continuity and momentum equations of the gas phase including the electrical body force induced by the movement of ions under the electric field. The RNG $k-{\varepsilon}$ model was used to analyze the turbulent flow. The particle concentration distribution was calculated from the convective diffusion equation of the particle phase. As the ion space charge increased, the particulate collection efficiency increased because the electrical potential increased over the entire domain in the ESP. The collection efficiency decreased and then increased, i.e. had a minimum value, as the EHD circulating flow became stronger when the electrical migration velocity of the charged particle was low. However, the collection efficiency decreased with the stronger EHD flow when the electrical migration of the particle was higher relatively. The collection efficiency of the model ESP increased as the turbulent diffusivity of the particle increased when the electrical migration velocity of the particle was low. However, the collection efficiency decreased for increasing the turbulent diffusivity when the electrical migration of the particle was higher relatively.

The time to measure the size distribution using Condensation Nuclei Counter(CNC) and Differential Mobility Analyzer(DMA) can be shortened by classifying particles ramping the DMA voltage exponentially and continuously. In measurement, particles sampled at different time are mixed together going through sampling tube and CNC. Because the size distribution is inversed by using detector responses to sampling time intervals in this accelerated method, the mixing effects give inversion errors to the size distribution. The mixing effects can be considered by appling the transfer function with mixing effects to the data inversion. The inversion considering this effects gives birth to the size distribution shifted to the opposite direction of the size scanning.

A numerical analysis has been perfonned to estimate the effect of turbulent penetration and thermal stratified flow in the branch lines piping. This phenomenon of thermal stratification are usually observed in the piping lines of the safety related systems and may be identified as the source of fatigue in the piping system due to the thermal stress loading which are associated with plant operating modes. The turbulent penetration length reaches to $1^{st}$ valve in safety injection piping from reactor coolant system (RCS) at normal operation for nuclear power plant when a coolant does not leak out through valve. At the time, therefore, the thermal stratification does not appear in the piping between RCS piping and $1^{st}$ valve of safety injection piping. When a coolant leak out through the $1^{st}$ valve by any damage, however, the thermal stratification can occur in the safety injection piping. At that time, the maximum temperature difference of fluid between top and bottom in the piping is estimated about $50^{\circ}C$.

The present study examined the possibility of $NO_{x}$ reduction in the high temperature industrial furnaces. duct burner of gas turbine cogeneration and two-stage gas turbine combustor. The experimental study was carried out for the diffusion flame of second stage combustor with the variations of oxygen concentration and supplying rate of hot exhaust gas from first stage combustor. It also examined the flammability range and $NO_{x}$ formation of the second stage combustor in which the fuel is supplying into the mixture of oxygen hot exhaust gas from first stage combustor. The results show that the enrichment of oxygen and increase of exhaust gas lead to increase the $NO_{x}$ up to 50 ppm with 23% $O_{2}$ condition.

In this study, reduced-scale experiments were conducted to understand smoke movements in tunnel fires with the natural ventilation. The 1/20 scale experiments were conducted under the Froude scaling since the smoke movement in tunnels is governed by buoyancy force. Six cases of experiments(pool diameter is 6.5cm, 7.3cm, 8.3cm, 10cm, l2.5cm and l5.4cm), in which vertical vents positioned 1m from the fire source symmetrically, were conducted in order to evaluate the effect of the vent on smoke movement. In case of heat release rate under 2MW, smoke front reached to the tunnel exit about 20 see delayed with ventilation and the smoke velocity was proportional to the power of the heat release rate. Temperature after the vent was lower than without vent. In case of l5.4cm pool, the temperature difference was about $50^{\circ}C$. It was confirmed that the thickness of smoke layer was maintained uniformly under the 35% height of tunnel through the visualized smoke flow by a laser sheet and the digital camcoder.

An experimental study on penetration integrity of the reactor vessel has been performed under external vessel cooling during a core melt accident. In this study a series of experiments are performed for the verification of the effects of coolant in the annulus between the ICI(In-Core Instrumentation) nozzle and the thimble tube and also the effects of external vessel cooling on the integrity of the penetration using the test section including only one penetration and $Al_{2}O_{3}$ melt as a corium simulant. The experimental results have shown that penetration is more damaged in the case of no external vessel cooling compared with the case of external vessel cooling. It is preliminarily concluded that the external vessel cooling is very effective measure for the improvement of the penetration integrity. Also it is confirmed from the experimental results that the coolant in the annulus reduces the melt penetration distance through the annulus and enhance the integrity of the reactor vessel penetration in the end.

In this paper, we have investigated the drying characteristics of waste sludge in a rotary dryer equipped with disintegration device. The sludge samples($moisture\;contents\;:\;70{\sim}85%\;W.B.$) used in the experiment were mainly a soybean and an alcohol fermentation sludges in food industry, or sewage sludge. It was studied by measuring the change of moisture contents that the shaft speed($operating\;range\;:\;100{\sim}250rpm$) affected on the drying rate and the particle size distribution of the dried solid. The dried granular products of the diameter of 1 to 5mm are obtained from sludge materials. In the results, as the increased of shaft speed, the particle size and the moisture contents of dried solid was decreased, and the retention time was increased.

The purpose of present study is to find design parameters and operating conditions of the HVAC system in a subway platform. The simulation was carried out for the flow, heat and mass transfer for heating, ventilating and air-conditioning(HVAC) environments in the subway platform. The steady-state. incompressible flow assumption and standard $k-{\varepsilon}$ turbulence model are adopted. The location of HVAC air inlet above platform and the volume flow rate of curtain air released from inlet B are chosen as main parameters in this study. The results of present study are following: In the case of existence of train, the heat and contaminant released under the train have no effect on the average temperature and mass fraction of contaminant in the platform, but heat released on the train has influence on the average temperature in the platform. Train acts as an obstacle to exhaust the contaminant in the platform, but has good effect on the average temperature in the platform.

Several design parameters for a 100 kW molten carbonate fuel cell stack was described. Approximately 170 cells are required to generate 100 kW at a current density of $125\;mA/cm^{2}$ with $6000\;cm^{2}$ cells. An overall heat balance was calculated to predict exit temperature. In order to limit the stack temperature in the range of $600-700^{\circ}C$, current load cannot exceed $75\;mA/cm^{2}$ at atmospheric operation. The 100 kW power is expected only under pressurization. Recycle of cathode gas by more than 50% is recommended to run the stack at $125\;mA/cm^{2}$ and 3 atm. Manifolds should be designed based on gas flow rates for the suggested operating condition.

UTC(Unglazed Transpired Collector) system has recently emerged as a new solar air heating technology. It is relatively inexpensive because it has not a glazed material. And it demonstrates efficient particularly for the applications in which larger wall area facilities with a high outdoor air requirement. Mathematical algorithm for UTC thermal modeling has been understood for further improvement of the system. EES and TRNSYS model of actual solar wall panel could be developed for computer simulations under other conditions. Computer models could be validated with the measured data from fixed outdoor test cell in KIER(Korea Institute of Energy Research). Major design parameters could be identified such as panel configuration and absorptivity and emissivity values for UTC design.

Performance of the water direct contact air conditioning system, in which heat and mass are transferred directly between air and water droplet, is simulated by semi-empirical method. This system improves transport efficiency compared to conventional indirect contact system and cooling, heating, dehumidification and humidification are attained with one unit. In this study, temperature and flowrate for air and water are measured in the various cooling and heating conditions, and correlations for $h_{c}A/c_{pm}$ are derived from these data. Cooling and heating characteristics of the water direct contact air conditioning system are investigated using correlations.

There is a chilled ceiling panel which carries out the air conditioning by radiation and convection between the room and cold ceiling panel surface. In order to verify heat transfer characteristics between them in cooling system with radiant chilled ceiling panel, analytical and experimental studies were performed for various design and operating parameters such as tube space and diameter, inlet water temperature, mass flow rate, cooling load, and so on. In this study, we found that the tube space and inlet water temperature were more important elements than the tube diameter and water flow rate for the performance of radiant chilled ceiling panel. The cooling capacity of the radiant chilled ceiling panel had the maximum value of $65W/m^{2}$ because the highest cooling capacity was limited by the condensation on the panel surface. The results of comparison between numerical analysis and experiment showed a resonable agreement qualitatively, especially for low cooling capacity.

In this paper, CFD technique has been used at design stage to predict space air distribution in a cycle stadium with air circulation system. An air circulation flow of 0.67 rev./min was observed at computed results in the stadium space with and without air circulation system. Comparing the thermal comfort of the two models with or without air circulation system showed that the thermal environment in the former was superior in the latter. Energy savings could be achieved for the model with air circulation due to its lower air inflow temperature.

In this study, the minimum heat flux conditions are experimentally investigated for the spray cooling of hot plate. The hot plates are cooled down from the initial temperature of about $900^{\circ}C$, and the local heat flux and surface temperatures are calculated from the measured temperature-time history. The results show that the minimum heat flux point temperatures increase linearly resulting from the propagation of wetting front with the increase of the distance from the stagnation point of spray flow. However, in the wall region, the minimum heat flux point temperature becomes independent of the distance. Also, the experimental results show that the velocity of wetting front increases with the increase of the droplet flow rate.

The present research is an experimental study of subcooled flow boiling behavior using flat, microporousenhanced square heater surfaces in pure FC-72. Two $1-cm^{2}$ copper surfaces, one highly polished (plain) and one microporous coated, were flush-mounted into a 12.7 mm square, horizontal flow channel. Testing was performed for fluid velocities ranging from 0.5 to 4 m/s (Reynolds numbers from 18,700 to 174,500) and pure subcooling levels from 4 to 20 K. Results showed both surfaces' nucleate flow boiling curves collapsed to one line showing insensitivity to fluid velocity and subcooling. The log-log slope of the microporous surface nucleate boiling curves was lower than the plain surface due to the conductive thermal resistance of the microporous coating layer. Both, increased fluid velocity and subcooling, increase the CHF values for both surfaces, however, the already enhanced boiling characteristics of the microporous coating appear dominant and require higher fluid velocities to provide additional enhancement of CHF to the microporous surface.

A simultaneous visualization of the behavior of bubbles and dry spots has been carried out for pool boiling of R-1l3 on a horizontal heater. From video imaging and image processing analysis, the formation of bubbles and dry spots occurs simultaneously, and therefore they should be considered as a synchronized concept rather than independent identities. The dry spot density is equivalent to the active site density in the region before CHF. At CHF point, large dry areas due to the coalescence of neighboring dry areas cover the heater surface.

An experimental study on critical heat flux (CHF) has been performed for water flow in a uniformly heated vertical 3 by 3 rod bundle under low flow and a wide range of pressure conditions. The objective of this study is to investigate the parametric trends of CHF with 3 by 3 rod bundle test section where three unheated rods exist. The general trends of the CHF are coincident with previous understandings. At low flow and system pressure above 3 MPa, some critical qualities are larger than 1.0 due to counter-current flow in test sections. Since there is a supply of water to the heated section from unheated section, the maximum CHFs at system pressure between 2 and 4 MPa are not shown.

We modeled the liquid film dryout(LFD) process for both tube and annulus which have uniformly heated vertical channels. We set phenomenological initial conditions in the model. The initial void fraction on the onset of the annular flow location is derived from the physical chum-to-annular flow criterion with the help of the drift-flux-model. The initial thermodynamic-equilibrium-quality is calculated by iteration with the flow quality to find the onset of the annular-flow location. Present model tends to predict very well at the lower exit quality but under-estimates at the higher exit quality. We found that the prediction error of the present model is gradually bigger as the inlet subcooling approaches near the saturation. We obtained excellent results for both tube and annulus channels as the mean of 0.97 and root-mean-square error of 11% for the number of 3883 experimental data on tubes, and of 0.96 and of 12% for 593 on annuli. The present model extended the applicable range to the relatively low exit quality region than previous LFD models.

The minimum heat flux conditions are experimentally investigated for the subcooled liquid film flow on the horizontal plate. The experimental results show that the minimum heat flux point temperature becomes higher with the increase of the velocity and the subcooling of the liquid film flow. However, the effect of distance from the leading edge of the heat transfer plate on the minimum heat flux is almost negligible. Also, the experimental results show that the propagation velocity of wetting front increase with increasing the velocity and the subcooling of the liquid film flow.

In this study, numerical investigation has been performed on the spreading and solidification of a droplet impacting onto a solid substrate in the thermal spray process. The finite difference method with volume-of-fluid approach is used to analyze the free surface flow and the source-based enthalpy method is employed to model the latent heat release during the solidification. In this work, the numerical model is validated through the comparison of the present numerical result with experimental data available for the flat substrate.

The plate heat exchanger is characterized. by low pressure drop and high heat transfer coefficient. The experimental study has been performed on the condensation heat transfer and pressure drop characteristics of the plate heat exchangers in this study. In the present study, a brazed type plate heat exchanger was investigated at a chevron angle of $45^{\circ},\;55^{\circ},\;and\;70^{\circ}$ with R410A. Condensation temperatures were varied from $20^{\circ}C\;and\;30^{\circ}C$, and mass flux was ranged from $13{\sim}34\;kg/m^{2}s$ with constant heat flux ($=5\;kw/m^{2}$). The heat transfer coefficient and pressure drop increased with the chevron angle. Average condensation heat transfer coefficients and pressure drops are decreased with increasing condensation tempeature.

The purpose of this paper is to propose the experimental method of thermal properties of Phase Change Materials (PCMs) by using T-history method. As far, in order to measure the heat of fusion and specific heat of PCMs, conventional thermal analysis methods such as DSC and DTA have been used. Because these methods test very small samples, thermal properties of samples are usually different from those of materials consisting of several components. For these reasons, T-history method, the simple measurement method of the heat of fusion and specific heat of PCMs have been performed. In this paper, we investigated the thermal properties of low temperature PCMs(below $0^{\circ}C$) under the charging process by using T-history method. The results are compared to those of DSC method. The T-history method will be useful for selection of the best PCM from lots of candidates and development of new PCMs.

The present study is an experimental investigation of nucleate boiling heat transfer mechanism in pool boiling from wire heaters immersed in saturated FC-72 coolant and water. The vapor volume flow rate departing from a wire during nucleate boiling was determined by measuring the volume of bubbles, varying $25{\mu}m,\;75{\mu}m,\;and\;390{\mu}m$, from a wire utilizing the consecutive-photo method. The effects of the wire size on heat transfer mechanism during a nucleate boiling were investigated by measuring vapor volume flow rate and the frequency of bubbles departing from a wire immersed in saturated FC-72. One wire diameter of $390{\mu}m$ was selected and tested in saturated water to investigate the fluid effect on the nucleate boiling heat transfer mechanism. Results of the study showed that an increase in nucleate boiling heat transfer coefficients with reductions in wire diameter was related to the decreased latent heat contribution. The latent heat contribution of boiling heat transfer for the water test was found to be higher than that of FC-72. The frequency of departing bubbles was correlated as a function of bubble diameters.

Phenomena of direct contact condensation (DCC) heat transfer between steam and water are characterized by the transport of heat and mass through a moving steam/water interface. Application of the phenomena of DCC heat transfer to the engineering industries provides some advantageous features in the viewpoint of enhanced heat transfer. This study proposes a simple condensation model on the steam jets discharging into subcooled water from a single horizontal pipe for the prediction of the steam jet shapes. The analysis model was derived from the mass, momentum and energy equations as well as a thermal balance equation with condensing characteristics at the steam/water interface for the axi-symmetric coordinates. The extremely large heat transfer rate at the steam/water interface was reflected in the effective thermal conductivity estimated from the previous experimental results. The analysis results were compared with the experimental ones. The analysis model predicted that the steam jet shape (i. e. radius and length) was increasing as the steam mass flux and the pool temperature were increasing, which was similar in trend to that observed in the experiment.

Solid-liquid phase change (i.e. melting or solidification) occurs in a number of situations of practical interest. Some common examples include the melting of edible oil, metallurgical process such as casting and welding, and materials science applications such as crystal growth. Therefore, due to the practical importance of the subject, there have been a large number of experimental and numerical studies of problems involving phase change during the past few decades. Also, this study presented the effective way to enhance phase change heat transfer.

It is derived for the temperature profile in a cylindrical and planar shape capillary pumped loop evaporators subject to a uniform heat flux prior to the initiation of boiling using the finite difference method. The results of the analysis allow for the determination of applied power levels for which nucleation is likely to occur only within the vapor grooves of the evaporator while maintaining subcooling in the liquid core, thereby increasing the likelihood of a successful startup. Also, limits are found for which additional increases in the applied heat flux do not increase the temperature difference between the vapor grooves and the wick-liquid core interface. Several advantages of larger diameter evaporators observed experimentally in startup are explained and quantified by the model. This analysis is appropriate for standard capillary pumped loop evaporators during a fully-flooded startup as well as starter pump designs and loop heat pipes.

This paper describes total flow control of an ice slurry system for pump energy saving. Similar turbo machinery has a characteristic that input power ratio is proportional to the three time of revolution speed ratio. To reduce the energy cost of brine pump in ice slurry storage systems, inverter is adapted instead of 3-way valve to control the speed of brine pump motor. One type of cooling load profile was used as driving load of the system, generated by a boiler and warm water storage tank. As results of the laboratory test, energy consumption and cost of the pump were reduced by 11.4%.

In this study, the performance of the simultaneous heating & cooling water making system using R134a was investigated by simulation. The most important effect upon heating COP was intermediate pressure depending on input water temperature. With the input water temperature of $10^{\circ}C\;and\;20^{\circ}C$, optimum intermediate pressure were 923 and 1040kPa, respectively. At that optimum intermediate pressure, the maximum heating COP of the system operated between $0^{\circ}C$ evaporating temperature and $70^{\circ}C$ condensing temperature were 4.15 and 3.83. With installation of the subcoolers in high or low pressure section, the system COP was increased by reducing the refrigerant mass flow rate. Under the optimum pressure and $10^{\circ}C$ input water temperature, it was found that heating COP was maximized when the low-subcooler and high-subcooler capacity rate were taken by 14% and 13%, respectively.

The objective of this paper is the investigation of the relationships between the surface roughness and film evaporative characteristics of the surface. For example, when the droplet of liquid is in contact with the solid surface, its behavior strongly depends on the surface characteristics. The material properties and geometry - profile shape, waviness, roughness - of the surfaces have strongly influenced on the wettability of the droplet. To investigate the effect of the surface roughness on the film evaporation, firstly, the characteristics of wettabilities were studied according to contact angle and surface tree energy of specimens with various roughness heights. Secondly, the experimental test were carried out on capacities of the tubes diversly roughened by using different kinds of emery papers. Finally, the relationships between the film evaporation characteristics and surface roughness were explained by means of the correlation of contact angle and surface free energy with surface roughness and the influences of surface tree energy on the heat transfer performance.

This paper presents a modelling of thermal discharge performance for a static ice-on-coil ice-storage tank. Through the present study, discharging characteristics were examined with the existing results of theoretical and numerical heat transfer analyses. Also, an experiment was conducted to obtain a real set of discharge performance. The thermal effectiveness, the ratio of the actual heat transfer rate to the maximum possible heat transfer rate, decreased when the stored energy decreased during discharging period. And the effectiveness increased as the coolant flow rate through the storage increased, of which increasing rate decreased abruptly near the maximum and the minimum stored energy. An empirical correlation was obtained from the experimental and the numerical analysis data.

This paper describes an optimal scheduling for ice slurry systems for energy cost saving. The optimization technique applied in the study is the dynamic programming method, for which the state variable is the storage in the ice storage tank and the control variable is the state of chiller's on-off switching. Though the costs during charge period is included in optimization by taking the average cost of ice per hour for slurry making, the time horizon for the simulation is limited building cooling period because accurate charge rate from the ice maker into the ice storage tank cannot be estimated during the charge period. In the operating simulation after optimizing procedure, energy consumption and operating cost for the optimal control are calculated and compared with them for a conventional control with one case of cooling load profile.

This paper discusses the effect of varying LiBr solution circuits flow rates for a direct fired double effect commercial absorption chiller in the parallel flow configuration. The effects of solution flow rates have been investigated for generator, condenser, solution heat exchanger, absorber and evaporator. According to the result of this work, it was found that sensible heat rate of generator increases and refrigerant vapor generated in that decreases when inlet solution flow rate of that increases. As solution flow rate of absorber increases, the degree of superheat increases because of decreasing solution heat exchanger efficiency. The flashing vapor at the top of absorber increases in proportion to the degree of superheat whileas decreases cooling capacity inversely.

A fundamental study on the under water harvest-type ice storage system and its temperature characteristics in ice storage system was performed experimentally. The experiments were conducted by changing the inlet refrigerant temperature of an evaporator to analyzing the thermal fluid motion inside the ice storage tank. From the experimental results, the cold storage characteristics were investigated by measuring the axial and radial temperature variations inside the ice storage tank with respect to the inlet and outlet refrigerant temperatures of an evaporator. In case of the under water harvest-type ice storage system, thermal fluid. motion inside the ice storage tank was shown differently in comparison with that of other ice storage systems. During the cooling storage process, there was no supercooling phenomenon in the ice storage tank. These results show the characteristic of this system and the possibility of application to other fields.

A numerical model which simulates the heat and mass transfer processes within a counter-current plate type generator for ammonia/water absorption refrigerators was developed. Ammonia/water solution flows downward under gravity and ammonia/water vapor generated by flow boiling flows upward. The flow pattern within the generator was assumed to be a bubbly flow, and the liquid and vapor phase were assumed to be saturated. It was shown that the boiling of ammonia occurred mainly in the upper part of the generator. The effects of the generator length, the wall temperature and the mass flow rate of ammonia/water solution into the generator on the generation of ammonia/water vapor were investigated.

A three-dimensional numerical analysis of the flow and heat transfer characteristic of wood-flour-filled polypropylene melt in an extrusion die was carried out Used for this analysis were Finite Concept Method based on FVM, unstructured grid and non-Newtonian fluid viscosity model. Temperature and flow fields are closely coupled through temperature dependent viscosity and viscous dissipation. With large Peclet, Nahme, Brinkman numbers, viscous heating caused high temperature belt near die housing, Changing taper plate thickness and examining some predefined parameters at die exit investigated the effect of taper plate on velocity and temperature uniformities. In the presence of taper plate, uniformity at die exit could be improved and there existed an optimum thickness to maximize it.

A finite element model for the process of squeeze casting for metal matrix composites (MMCs) in cylindrical mold is developed. The fluid flow and the heat transfer are the fundamental phenomena in the squeeze casing process. To describe heat transfer with solidification of molten aluminum, the energy equation in terms of temperature and enthalpy are applied to two dimensional axisymmetric model which is similar to the experimental system. And one dimensional flow model is employed to simulate the transient metal flow. The direct iteration technique was used to solve the resulting nonlinear algebraic equations. A computer program is developed to calculate the enthalpy, temperature and fluid velocity. Cooling curves and temperature distribution during infiltration and solidification are calculated for pure aluminum. The temperature is measured and recorded experimentally. At two points of the perform inside and one point of the mold outside, thermocouple wire are installed. The time-temperature data are compared with the calculated cooling curves. The experimental results show that the finite element model can estimate the solidification time and predict the cooling process.

One of the principal components of the KSTAR (Korea Superconducting Tokamak Advanced Research) tokamak structure is the vacuum vessel, which acts as the high vacuum boundary for the plasma and also provides the structural support for internal components. Hyundai Heavy Industries Inc. has performed the engineering design of the vacuum vessel. Here the overall configuration of the KSTAR vacuum vessel was briefly described and then the design methodology and the analysis results were presented. The vacuum vessel consists of double walls, several ports, leaf spring style supports. Double walls are separated by reinforcing ribs and filled with baking/shielding water. The overall external dimensions of the main body are 3.39 m high, 1.11 m inner radius, 2.99 m outer radius, and made of SA240-316LN. The vacuum vessel was designed to be capable of achieving the base pressure of $1\times10^{-8}$ Torr, and also to be structurally capable of sustaining the vacuum pressure, the electromagnetic and thermal loads during plasma disruption and bakeout, respectively. The vacuum vessel will be baked out maximum $150^{\circ}C$ by hot pressurized water through the channels formed between double walls and the reinforcing ribs. A 3-D temperature distribution and the resulting thermal loads in the vessel were calculated during bakeout. It was found that the vacuum vessel and its supports were structurally rigid based on the thermal stress analysis. The maximum electromagnetic loads on the vacuum vessel induced by eddy and halo currents resulting from the engineering plasma radial and vertical disruption scenarios have been estimated. The stress analyses have been performed based on these electromagnetic loads and the resulting stresses at he critical locations of the vacuum vessel were within the allowable stresses.

Theoretical approach was taken to the whole CNG refueling process. In particular, this study was focused on the prediction of flow rate at any given piping configuration of CNG system, in order that a simulation program for the CNG refueling system should be developed. The simulation result of refueling process was compared with experimental result obtained from various kinds of fueling configuration. The simulation results showed a satisfactory agreement within 10% errors in fueling time, fueling amount, and residual pressure. The developed program would be used a good engineering tools for estimating fueling performance for a any given CNG station.

A numerical study of three-dimensional fluid flow and heat transfer in the metering section of a single screw extruder has been performed. The mathematical model for the screw channel is simplified by unwound channel and fixing the coordinate system to the screw. The pressure boundary and the prescribed mass flow rate conditions are imposed on the inlet and outlet, respectively. The commercial code STAR-CD based on the finite volume method is used to obtain the results of the present work. The computation of the reverse flow, which cannot be computed by the marching-type 3-D model, is performed in the present study.

Advanced Pressurized Reactor 1400(APR-1400), which is a standard evolutionary advanced light water reactor(ALWR), has been developed from 1992 as one of long-term Government Project(G-7). The APR-1400 is designed to operate at the rated output of 4000MWt to produce an electric power output of around 1450MWe. The balance of plant (BOP) for the secondary system consists of main steam, feedwater, condensate, turbine generator and auxiliary system. In this paper, we describe the major design features of secondary component, balance of plant configuration, and then the turbine cycle thermal performance evaluation using PEPSE code.

Experimental schemes that enable characterization of phase-change phenomena in the micro scale regime is essential for understanding the phase-change kinetics. Particularly, monitoring rapid vaporization on a submicron length scale is an important yet challenging task in a variety of laser-processing applications, including steam laser cleaning and liquid-assisted material ablation. This paper introduces a novel technique based on Michelson interferometry for probing the liquid-vaporization process on a solid surface heated by a KrF excimer laser pulse (${\lambda}=248nm,\;FWHM=24\;ns$) in water. The effective thickness of a microbubble layer has been measured with nanosecond time resolution. The maximum bubble size and growth rate are estimated to be of the order of $0.1{\mu}m\;and\;1\;m/s$, respectively. The results show that the acoustic enhancement in the laser induced vaporization process is caused by bubble expansion in the initial growth stage, not by bubble collapse. This work demonstrates that the interference method is effective for detecting bubble nucleation and microscale vaporization kinetics.

As the temperature of liquid under negative pressure approaches the absolute zero, the nucleation process due to thermal fluctuations hardly occurs. Instead of this mechanism, quantum fluctuations may lead the formation of nucleus for new phase in metastable state. In this study, the thermal as well as quantum nucleation bubble in liquid helium under negative pressure was investigated theoretically. The energy barrier against nucleation was estimated by molecular interaction due to the Londom dispersion force. It is shown that the phase transition from liquid to vapor in is possible due to the quantum tunneling below 0.2 K for Helium-4 and 0.1 K for Helium-3, at negative pressures close to the ideal tensile strength at which every liquid molecules become bubbles simultaneously.

Numerical Aanlysis is made on the thermal performance of micro heat pipe in a axial flat grooved channel. The flow of liquid and vapor is investigated in trapezoidal grooves and the effect of variable shear stress along the interface of the liquid and vapor considered. The results from this study are obtained in the axial variation of pressure difference between vapor and liquid, contact angle, velocity of liquid and vapor and so forth. In addition, maximum heat transport capacity of micro-heat pipe is provided by varying the operation temperature, and compared with that from Schneider and Devos's model in which the interfacial shear stress is neglected.

There are various wick types for heat pipe. In the present study, the manufacturing technology of a sintered wick among various wick types is discussed. The sintering technology using metal has been applied broadly in the field of electronic-telecommunication as well as heat pipes. A study of manufacturing procedure and characteristic of sintered wick for heat pipe have been performed. Copper powder was used as wick material and stainless steel as a mandrel. A manufacturing technology of the mandrel for arranging vapor core in heat pipe, a sintering technology by first or second times and operating temperature for sintering, the measurements of a porosity, pore size, and pore distribution of sintered wick were considered. In the meantime, a heat pipe with sintered wick has been manufactured and a performance test of the heat pipe has been performed in order to review cooling performance. The performance test results for the 4mm diameter heat pipe with the sintered wick shows the stability since the temperature difference between a evaporator and a condenser of the heat pipe is less than $4.4^{\circ}C$, and thermal resistance is less than $0.7^{\circ}C/W$.

Radiative transfer by nongray gas mixtures with nonuniform concentration and temperature profiles were studied by using the statistical narrow-band model and ray-tracing method with the sufficiently accurate $T_{60}$ quadrature set. Transmittances through the nonhomogeneous gas mixtures were calculated by using the Curtis-Godson approximation. Three different cases with different temperature and concentration profiles were considered to obtain benchmark solutions for nongray gas mixtures with nonuniform concentration and temperature profiles. The solutions obtained from this study were verified and found to be very well matched with the previous solutions for uniform gas mixtures. The results presented in this paper can be used in developing various solution methods for radiative transfer by nongray gas mixtures.

Film cooling characteristics has been investigated numerically with the aid of FLUENT software for the sunk or the lifted upstream wall from the slot injection exit. In this study, with the fixed blowing ratio of 1 and the fixed coolant injection angle of $30^{\circ}$, the downstream flow field and the downstream temperature field were examined in terms of velocity vector, turbulent kinetic energy, temperature contours, and downstream wall temperature. Upstream wall was sunk or lifted from 1d to 5d(d=slot width). The result shows that the up-Id upstream wall has the best film cooling performance. This is due to the fact that the up-1d upstream wall configuration reduces velocity gradient just enough to minimize the turbulent mixing between the mainstream and the coolant just off the slot exit.

In the present study, local heat transfer rates for co-rotating disks with two modified hubs having ventilation holes are investigated for Rossby number of 0.04, 0.1 and 0.35 to evaluate the influence of incoming flows through hub holes. A naphthalene sublimation technique is employed to determine the detailed local heat/mass transfer coefficients on the rotating disks using the heat and mass transfer analogy. Flow field measurements are conducted using Laser Doppler Anemometry (LDA) and numerical calculations are performed simultaneously to analyze the flow patterns induced by the disk rotation. The basic flow structure in a cavity between co-rotating disks consists of three regions; the solid-body rotating inner region, the outer region with turbulence vortices and the shroud boundary layer region. The heat/mass transfer. rates on the co-rotating disks are very low near the hub due to the solid-body rotation and those increase rapidly in the outer region due to turbulence mixing. The modified hubs with ventilation holes enhances significantly the heat/mass transfer rates on the region near the hub. The results also show that the heat transfer of Hub-2 is superior to that of Hub-1, but Hub-1 is more profitable for destructing the solid-body rotating inner region.

Cooling characteristics using convection and conduction heat transfer in a parallel channel with extruding heat sources are studied numerically. A two-dimensional model has been developed for numerical prediction of transient, compressible, viscous, laminar flow, and conjugate heat transfer between parallel plates with uniform block heat sources. The finite volume method is used to solve this problem. The considered assembly consists of two channels formed by two covers and one PCB which has three uniform heat source blocks. Five different cooling methods are considered to find efficient cooling method in a given geometry and heat source. The velocity and temperature fields, local temperature distribution along surface of blocks, and the maximum temperature in each block are obtained.

Heat transfer analysis of a iron core type linear motor for surface mounting device applications was considered in this study. In order to avoid the complex conjugate problem a fluid flow regime and a solid regime were considered separately. First, film coefficients of the moving parts were evaluated from computational fluid dynamic analysis and those of the stationary parts from the existing empirical or analytic correlations. And then, by applying them, internal and external temperatures of the linear motor pal1s were computed through finite element analysis. Both computation and measurement were carried out with respect to motor driving power. The measurement did not exhibit a linear temperature variation trend with respect to motor power while the computation revealed a linear correlation. Nonetheless, the computations agreed with the measurements within an error range of 20%. It indicates that an adequate heat transfer model for the reciprocative coil assembly may help more exact prediction.

This paper reports a result of numerical heat transfer analysis for heat spreader with various ribs. Four different ribs are compared in this study. In general, the heat transfer on a vertical plate is enhanced when a rib is attached as the surface area increases, and the growth of the boundary layer is interrupted. However, for a low flow less than 0.1m/s, it is observed that the heat transfer is sensitive to the height of a rib: it decreases as the height increases. Among the four ribs, the H-shaped rib showed better performance than other ribs.

Numerical analysis on the three-dimensional laminar flows (Re=1000) and heat transfer in a rectangular channel with punched longitudinal vortex generator have been conducted to explore the heat transfer enhancement and the combined effect of the angle of attack ${\alpha}$ and the lovour angle ${\beta}$. Rectangular winglets have been used as vortex generators. Velocity and temperature fields and spanwise averaged Nu and friction factor were presented. Enhancement of heat transfer and flow loss penalty are evidenced. The results show performance characteristics allowing a reduction in heat transfer surface area of 62% for fixed heat duty and for fixed pumping power compared with that of channel flow without vortex generator. However, adding lovour angle to the vortex generator shows no positive effect on the heat transfer enhancement.

Flow and thermal characteristics of cooling system for the robot controller were numerically as well as experimentally investigated. To obtain the overall flows within controller, the system level solutions were analysed at first and then the board level solutions were pursued to understand the detailed flow and temperature fields near the main board which have a significant influence on the cooling of electronic components. The evaluation for a performance of the heat exchanger was conducted on the basis of the obtained flow and temperature patterns. The results showed that the heat exchanger made a small contribution to the cooling of controller and caused an increase of the temperature in CPU.

This paper deals with the development of a new method that can obtain heat transfer coefficient and reference tree stream temperature simultaneously. The method is based on transient heat transfer experiments using two narrow-band TLCs. The method is validated through error analysis in terms of the random uncertainties in the measured temperatures. It is shown how the uncertainties in heat transfer coefficient and tree stream temperature can be reduced. The general method described in this paper is applicable to many heat transfer models with unknown free stream temperature.

Single-phase heat transfer coefficients and pressure drops of R-22 were measured in smooth, horizontal copper tubes with inner diameters of 3.36 mm, 5.35 mm, 6.54 mm and 8.12 mm, respectively. The experiments were conducted in the closed loop, which was driven by a magnetic gear pump. Data are presented for the following range of variables: Reynolds from 1000 to 20000. Single-phase heat transfer coefficients increased by $10{\sim}30%$ as the inner diameter of tube was reduced and it was found that a well-known previous correlation, Gnielinski's correlation, was not suitable for the small diameter tubes. But the pressure drop in the small diameter tubes have been shown slightly deviations with Blauius' correlation. Based on an analogy between heat and mass transfer, the new heat transfer correlation is proposed to predict the experimental data successfully.

An experimental study was carried out on the characteristics of the mixed-convection heat transfer from a protruding heat source module which had uniform heat flux and was located on a flat plate in the inclined channel. The effects of the inclined channel(${\varphi}=0{\sim}90^{\circ}$) was studied for the input power($Q=3,\;7W$) and inlet air velocities($V_{i}=0.1{\sim}0.9m/s$). Experimental results indicate that the input power was most effective parameter on the temperature differences between inlet air and module. The effects of the inclined angle was negligible when the inlet velocities were above 0.5m/s and 0.9m/s at Q = 3W, 7W respectively. As the inclined angle of the channel increases, the temperatures of the module are decreased. So we obtained the best condition on the adiabatic board at the vertical channel.

The pressure drop and heat transfer of the spirally indented tubes were measured and analyzed. Eight sample tubes of indentation depth 0.4, 0.7mm and indentation pitch 10, 14, 20, 26mm were used in this experimental tests. And all the tubes have same outer diameter of 16mm, and same indentation start number of I. Air was used as the internal fluid from 10000 to 50000 for Reynolds Number. The friction factors and heat transfer coefficients have increased when indentation depths increase and indentation pitches decrease. Finally, the correlations were made between the effect of the tube geometry and characteristics of tubes for the pressure drop and heat transfer.

Quenching phenomena of hemispherical downward facing convex surfaces with narrow gaps have been investigated experimentally. Experiments employed test sections having 1 and 2 mm in gap thickness and 1 atm in system pressure. From interpretations of the temperature and the heat flux history, it was found that the flooding inside the gap was restricted by CCFL phenomena and quenching process was propagated from lower to upper region of the internal copper shell. The ratio of the maximum heat fluxes at 1 mm to 2mm in gap thickness was the almost same that obtained by steady state experiments. The quenching scenario of the hemispherical downward facing surface with narrow gap has been suggested.

Analysis of radiation heating system for producing 60" size PDP panels was carried out using $RADCAD^{TM}$ software. Optimum arrangement of infrared heating elements was found to obtain uniform temperature distribution in PDP panel during heating. Heating capacity of each heater was determined to obtain an appropriate maximum panel temperature. Parametric study to find the effect of design parameters such as the thermophysical and optical properties of glass and cooling system was carried out. As a reference system, about 35 kW heating capacity was chosen to obtain about 800 K maximum panel temperature after 30 minute heating. The maximum temperature difference in panel was below 20 K. The maximum/minimum and its difference in the panel were very sensitive to the variation of the emissivity of glass and cooling block.

A mathematical model has been developed to describe the turbulent and reversed flow with convective heat transfer in a cylindrical combustion chamber. By using the mathematical model for high temperature flow enables the trends in overall heat transfer rates to be predicted. The model was applied to the design of the combustion chamber. The influences of the size of air inlet and inlet velocity were investigated for process optimization. Through modelling work it is found that the heat transfer rate to the chamber wall may be enhanced by adjusting the air flow and heat transfer pattern through selecting the air inlet condition. Internal plate has less influence to the heat transfer characteristics.

The entrainment of air into spray jets has been considered. Entrainment is defined as the quantity of ambient gas that is drawn into a spray. Numerical study is performed to investigate an air entrainment into spray jets and compared with results of experiment of air entrainment. Experimental measurements were performed with PDA and PIV system. Experimental and numerical results show that the air entrainment was affected droplet size and velocity.

The characteristics of conjugated heat transfer in the inclined channel was experimentally investigated. The simulated module is attached to the bottom of the inclined channel and is heated with constant heat flux. The experimental parameters of this study are input power (Q = 3, 7W), inlet air velocity ($V_{i}=0.1{\sim}0.9m/s$) and inclined channel angle (${\varphi}=0{\sim}90^{\circ}$). The results show that input power was most effective parameter on the temperature differences between module and air. As the inclined channel angle increases, the temperatures of the module are increased. And we obtained the best condition on the conductive board when ${\varphi}=0^{\circ}$.

An experimental study on the performance evaluation of a brazed plate heat-exchanger with 10RT of normal cooling capacity has been carried out. In the present study, a brazed type plate heat exchanger was tested at a chevron angle $25^{\circ}$ with refrigerant R-22. Mass flux was ranged from $23\;to\;58kg/m^{2}s$ in condensation, and from $22\;to\;53kg/m^{2}s$ in evaporation. The heat transfer coefficient and pressure drop increased with the mass flux increases. The water side pressure drop increased with the cooling water flow rate and chilled water flow rate increases, while mass flux has little affect. It is also shown that the system performance can be improved by enlarging condensation heat transfer area.

In the present study a numerical simulation is performed on a natural convection inside a square cavity with a vibrating wall. The study has been conducted varying the heat transfer rate, wall excitation frequency and also the orientation of the cavity. The temperature and velocities inside the cavity was observed and also, the heat transfer coefficients on the heating wall was seen. From the results, it can be seen that the temperature inside the cavity decreases when excited with the proper frequency and the heat transfer coefficient increased with cavity inclination angle, ${\theta}$. It is also found from the results that flow resonance is occurred near the inclination angle ${\theta}=90^{\circ}$.

Two perforated plates are used to investigate local heat/mass transfer characteristics in an impingement/effusion cooling system. A naphthalene sublimation method is conducted to determine the local heat/mass transfer coefficients on the upward facing surface of the effusion plate. The two plates are placed in parallel position with gap distances of 1, 2, 4 and 6 times of effusion hole diameter. The effects of hole arrangements of the plates are studied for staggered, square, and hexagonal arrays. The experiments are conducted at Reynolds number of 10,000 based on the effusion hole diameter. The results show that the smaller hole size in the staggered array has the higher transfer coefficients on the stagnation region due to the formation of higher momentum flows through the impingement holes. In the square array, heat/mass transfer on the target plate is more uniform as the number of impingement holes increases. High and uniform heat/mass transfer coefficients are obtained in the hexagonal array.

Three-dimensional natural convection from a discrete flush-mounted circular heat source on the bottom of a cubic enclosure was studied by using a holographic interferometric tomography. The heat source was located at the off-center of the bottom plate so that three-dimensional temperature field can be achieved. A set of multidirectional holographic interferogram was recorded by employing a double-reference beam, double-exposure holographic technique in order to eventually reconstruct the temperature fields. The recorded interferometric data appear good enough to be further processed to extract optical pathlength data from them and finally reconstruct the temperature fields. A complete analysis of the temperature fields including the field reconstructions and comparison with thermocouple measurements is underway and will be reported shortly.

Experiments on the condensation and evaporation heat transfer characteristics inside plate heat exchanger with R134A are performed in this study. The test plate heat exchangers in 45o, 55o and 70o shevron angle are used. Varying the mass flux of the refrigerant and the saturation temperatures, the average heat transfer coefficients are investigated. It is shown that the heat transfer is increased with increasing shevron angle. Experiments results show that average condensation heat transfer coefficients are decreased with increasing condensation temperature but those of evaporation are increased with increasing evaporation temperature.

In a square cavity, the flow phenomena in the surrounding of the bubble attached at the upper cooled solid wall were studied by using a thermo-sensitive liquid-crystal tracer and image processing techniques. This method offers the advantage of measuring the entire flow field in a selected plane within the fluid at a given instant of time in contrast to point by point method like T/C. Quantitative data of the temperature were obtained by applying a colour-image-processing to the. visualized image. As the growing of a bubble, In a bubble size appears the flow phenomena which the direction of flow is reversed in the entire temperature and flow field. The observed phenomena are described with regard to thermocapillary convection.

Present study deals with combustion characteristics and performance of U type radiation tube burner with fin which combustion capacity is 30,000kcal/hr and the maximum capacity of supply fuel is $30Nm^{3}/hr$. Temperature difference of radiation tube is about $173^{\circ}C$ at 25% capacity and this show relatively small temperature difference for convenient type. Thermal efficiency is satisfactory as $72{\sim}81%$. Also, radiative efficiency of radiation tube is $52{\sim}73%$. The efficiency of heat exchanger is $27{\sim}37%$. Therefore, radiative efficiency is improved to $1{\sim}10%$ after installing fin.

Comprehensive numerical computations are made of side-heated squire cavity which is exposed to zero mean g-jitter. Numerical solutions are acquires to the governing two-dimensional Navier-Stokes equations for a Boussinesq fluid. When the system is exposed to pure sinusoidal g-jitter inclined to the vertical axis, in spite of zero mean gravity there exist non zero net flow fields [8]. The resonance phenomenon are observed in moderate Rayleigh number. And, by comprehensive numerical work, unlike[5], it is found that they are related with the overshoot phenomenon of the sudden gravity up problem.

An experimental study is carried out to investigate the effects of air and water mass flow rates on cooling characteristics of mist impinging jet on a flat plate. Experiments are conducted with air mass flow rates from 0.0 to 3.0 g/s, and water mass flow rates from 5.0 to 20.0 g/s. An air-atomizing nozzle is used for the purpose of controlling air and water mass flow rates. In this study, a new test section is designed to obtain local heat transfer coefficient distributions. Heat transfer characteristics of the mist impinging jet are explained with the aid of flow visualization. Surface temperature and heat transfer coefficient distributions become more uniform as air mass flow rate increases, and that the increases in water flow rate mainly enhance cooling performance. Air mass flow rate weakly influences averaged heat transfer coefficient when water mass flow rate is low, but averaged heat transfer coefficient increases remarkably as air mass flow rate in case of high water mass flow rate.

An experimental investigation of heat transfer characteristics on confined jet impinging plate using multiple slot jets has been performed. The effects of jet Reynolds numbers(Re=2000, 3950, 5900, 7900), dimensionlesss slot-to-plate distances(H/B=2, 4, 6, 8) and jet-to-jet distances(S=16B, 20B, 24B, 30B) on the local and average heat transfer coefficients have been examined. To clarify local heat transfer characteristics, naphthalene sublimation technique were used. From the experimental results, it was found that the local and average heat transfer rates increase with increasing jet Reynolds number. Measurements of local heat transfer coefficients produced by multiple of slot jets have given an indication of the nature of the interaction between jets and of the uniformity of heat transfer obtainable with various arrangements. At S/B=20, Re=7900 and H/B=6, maximum average Nusselt number is obtained.

An experimental study was conducted to compare the heat transfer characteristics of an impinging slot jet and three kinds of impinging circular jets. Thermochromic liquid crystal with an image processing system was employed to measure the temperature of impinging wall where constant heat flux condition was applied. The distribution of convective heat transfer coefficients were then evaluated for eight nozzle-to-surface distance settings for each jet cases. The cooling effect was linearly proportional to the number of nozzles for circular jet cases at the same nozzle exit speed. However, the heat transfer under constant volume flow rate was the most at single circular jet. It was concluded that the overall convective heat transfer was better at the circular jets than the slot jet.

The heat transfer and flow measurements on a pedestal encountered in chip cooling. A uniform wall temperature boundary condition at the plate surface and on a pedestal was created using shroud method. Liquid crystal was used to measure the plate surface temperature. The jet Reynolds number (Re) ranges from 11,000 to 50,000, the dimensionless nozzle-to-surface distance (L/d) from 2 to 10, and the dimensionless pedestal diameter-to-height (H/D) from 0 to 1.0. The results show that the Nusselt number distributions at the near the pedestal exhibit secondary maxima at $r/d{\cong}1.0\;and\;1.5$. The formation of the secondary maxima is attributed to an create in the vortex by the pedestal.

To find the effective combinations of blowing ratio and injection angle for a straight slot film cooling, film cooling characteristics was investigated using both flow visualization experiment and numerical simulation. Injection angles from $15^{\circ}\;to\;50^{\circ}$ and blowing ratios from 0.2 to 3.0 were selected for the simulation. Comparison between experimental and numerical results shows a good agreement, for the case of the injection angle of $30^{\circ}$ and blowing ratio ranging from 0.55 to 2.0. Film cooling effectiveness was found to be an increasing function of blowing ratio. The effects of injection angle became prominent as the blowing ratio increases. An interesting phenomenon was found for the injection angle of $15^{\circ}$ : the lowest film cooling effectiveness for the blowing ratio smaller than 1.0, but the highest film cooling effectiveness for the blowing ratio greater than 2.0 within wide range of downstream region. There exist optimum injection angles corresponding to maximum film cooling effectiveness : injection angle of $25^{\circ}$ for the blowing ratio from 0.2 to 2.0, and injection angle of $15^{\circ}$ for the blowing ratio of 3.0. Present study provides a design combination among film cooling effectiveness, blowing ratio, and injection angle.

Advanced gas turbine engines employ turbine entry temperatures so high that cooling of the turbine blades is essential. The coolant flow introduces losses which need to be minimized, and therefore it is important that the minimum amount of coolant is used. This work presents the result of the one-dimensional analysis and the effect of the boundary conditions on coolant flow rate in gas turbine blades.

The objective of this study was to investigate the characteristics of air flow and heat transfer caused by trapezoid rods array in impinging air jet system. In this study, trapezoid rods have been set up on front of flat plate to act as a turbulence promoter. Local Nusselt numbers were determined as a function of three parameters: (a) the space from rods to heating surface(C=1, 2, 4mm), (b) the pitch between each rods(P=30, 40, 50mm), (c) the distance from nozzle exit to flat plate(H/B=2, 6, 10). And this research compared the above with the experiment without trapezoid rods. As a result, heat transfer performance was best under the condition of C=1mm and as the pitch is 30mm. In this case, maximum rate of heat transfer augmentation is about 1.9 times greater compared to that without trapezoid rods.

This study is to research the heat transfer characteristics in copper-water heat pipes with screen wick, #100. Recently, the semiconductor capacity of an electronic unit has been larger, on the contrary, its size is smaller than before. As a result, a high-performance cooling system is needed. Experimental variables are inclination angle and temperature of cooling water. The distilled water was used for the working fluid. At a inclination angle ${-6}^{\circ}$, #100 2layer screen mesh is shown the best heat transfer performance.

Agglomerated and nonagglomerated SiO2 particles are synthesized in furnace by the electrohydrodynamic spraying method and the vapor feeding method for the test particle generator this study. These polydispersed particles are classified with DMA to extract equal mobility particles. Then these particles are introduced into CNC (Condensation Nuclei Counter) to see the pulse height using Multi-channel Analyzer. The response characteristics of these two kinds of particles in CNCs (TSI CNC 3022 and 3025A) have been studied as a function of particle size using mono disperse particles classified by DMA. The results show that the higher drag resistance particles, so called agglomerated particles have generated a lower CNC pulse height than the spherical particles for these two different CNCs, which means the nonagglomerated particles may start to grow larger than the agglomerated particles.

Particle contamination on the data storage disks has been a serious problem for magnetic hard disk drive manufacturers. For high storage optical disks, such as DVD-ROM/RAM or NFR (near field recording) system, particle-induced damages can be also detected because only a few micrometer particles can prevent read/write signal from optical lens. The increasing areal density and smaller bit size accelerates particle induced damages on the optical disk. One of the methods to prevent particle contamination on the optical disk surface is to handle the disk enclosed in a cartridge like a modern DVD-RAM disk. However, even for a perfectly sealed disk drive, particles are found inside the drive. The other method is to improve disk surface characteristics. Particle contamination on the surface can be reduced by proper selection of disk coating materials. [n this paper, particle detachment ratios for CD (compact disk), DVD (digital versatile disk), HD (magnetic hard disk), HD with Jut lubricant, and aluminosilicate substrate HD were investigated. Surface roughness and surface energy of the test disks were compared with the particle detachment ratios. Proper substrate and lubricant characteristics to reduce particle contamination on the disk surface were found.

The present study conducted heat transfer analysis in multi-layer of CD-R. It is necessary to analyze heat transfer during the recording process to find optimum power and write strategy in CD-R. This study investigated effects of several parameters such as recording speed, laser power, layer thickness and thermal property. The calculated results presented temperature distribution in the multi-layer and detailed information of recording characteristics. Optimum laser power was estimated, comparing an optimum mark length with the calculated mark lengths. The results showed that the optimum laser power was influenced significantly by the layer thickness and the thermal properties of the dye.

Recently, as the environmental pollution becomes an important issue, aerosols which are the main components of the atmospheric pollution become the subject of a lot of researches. An impactor is a kind of sampling and measuring equipment for aerosols since 1970s. The impactor uses coating materials to obtain high performance. However, there are a lot of situations when coating materials are not available and in this case the classification efficiency of the impactor decreases significantly. In this study, the impaction plate of the impactor is cooled, therefore the impactor performs efficiently when coating materials are not available.

The relation between the aerodynamic diameter and some morphological parameters was studied for flame-generated aggregates. $SiO_{2}$ aggregates were generated from $SiCl_{4}$ in premixed methane/air flames. These aggregates were sampled and classified according to their aerodynamic diameter by a cascade impactor; moreover, computer program was developed and tested to find the equivalent area diameter and the fractal dimension of the aggregates. We calculated the parameters from the digitized images of the aggregate TEM micrographs. The aerodynamic diameters of the sampled aggregates were larger than $0.4{\mu}m$ in this experiment. In most cases, fractal dimension of their projection image was very close to 2.0 for these large aggregates. It was found that the equivalent area diameter of these aggregates was approximately three times larger than the Stokes' diameter of them.

A numerical technique for simulating the aggregation of charged particles was presented with a Brownian dynamic simulation in the free molecular regime. The Langevin equation was used for tracking each particle making up an aggregate. A periodic boundary condition was used for calculation of the aggregation process in each cell with 500 primary particles of 16 nm in diameter. We considered the thermal force and the electrostatic force for the calculation of the particle motion. The morphological shape of aggregates was described in terms of the fractal dimension. The fractal dimension for the uncharged aggregate was $D_{f}=1.761$. The fractal dimension changed slightly for the various amounts of bipolar charge. However, in case of unipolar charge, the fractal dimension decreased from 1.641 to 1.537 with the increase of the average number of charges on the particles from 0.2 to 0.3 in initial states.

This paper reports the cooling performance of a PDP(plasma display panel) with a heat spreader by means of numerical analysis. Due to the simplifications and assumptions inherent in the analysis, computed results are found to differ those of the experiment by 13%. Calculation shows a maximum temperature of $65^{\circ}C$ for the plasma glass, as opposed to the allowable temperature of $90^{\circ}C$, producing a temperature difference of $25^{\circ}C$ between the upper and lower regions. This is enough to cause cracks in the plasma glass. In order to avoid this, more ventholes are added at the upper center region of the back cover, thereby causing a $3^{\circ}C$ drop in the maximum temperature, which reduces the temperature difference to $12^{\circ}C$. The new design gives more uniform temperature distribution across the plasma glass.

The transport of charged particles in electrostatic precipitator is investigated by Eulerian numerical analysis. Collection efficiencies are calculated using various combinations of the assumptions about flow field, turbulent diffusivity and boundary condition at collecting electrode. The characteristics of calculated collection efficiencies are compared with the trends of published experimental results. It is found that the collection efficiency for the case using nonuniform turbulent flow field, nonuniform turbulent diffusivity and zero concentration boundary condition at collecting electrode is the most suitable for the prediction of collection efficiency of electrostatic precipitator.

An experimental study was performed to investigate adiabatic wall temperature and heat transfer coefficient around on a module with longitudinal fin heat sink cooled by forced air flow. In the first method, inlet air flow(1-7m/s) and input power(3-5W) was varied after a heated module were placed on an adiabatic floor($320{\times}550{\times}1mm^{3}$). An adiabatic wall temperature was determinated to use liquid crystal film(LCF). In the second method to determinate heat transfer coefficient, inlet air flow(1-7m/s) and the heat flux of rubber heater($0.031-0.062\;W/cm^{2}$) was varied after an adiabatic module was placed on rubber heater covering up an adiabatic floor. In addition, surface oil-film visualization were performed to characterize the macroscopic flow-field around a module.

The stem thrust required to unwedging a gate valve is influenced by the pressure and temperature when the valve is closed and by the changes in these conditions between closure and opening. "Pressure Locking" and "Thennal Binding" refer to situations where pressure and temperature effects cause the unwedging load to be much higher than normal. A model of these phenomena has been developed. The effects of pressure and temperature are analyzed to determine the change in this disk-to-seat "interference". Flexibilities or Stiffness of the disk and body strongly influence the unwedging thrust. Calculation and limited comparison to data have been performed for the RHR motor operated valve designs and scenario. Pressure changes can increase the unwedging thrust when bonnet pressure exceeds the pressure in the adjacent piping and temperature changes can increase the unwedging thrust when a temperature change after closure produces an increase in the disk-to-seat interference.

In this paper, we constructed four valve type pulse tube refrigerator and found the characteristic of orifice (needle valve) opening for using phase shifter and optimum cycle time - The valve timing was controlled by the digital timers. The experimental results showed the optimum frequence and cycle time at operating conditions. The results showed that the optimum process time existed and the rate was same at each operating frequence. The no- load temperature of the refrigerator was 85K.

This study represents numerical analysis on the thermal and fluid flow of the air layer in a solar collector. The boundary conditions was assumed that the top and bottom wall of the air layer have a heating and cooling surface. respectively. and this calculation model have a solid body with a cooling temperature of $20^{\circ}C$. As the results of simulations. the magnitudes of the velocity vectors and isotherms are increased proportionally to the tilt angles. As the tilt angle is increased. the mean Nusselt numbers are increased and the maximum value of the mean Nusselt number was appeared at tilt angle ${\theta}=75^{\circ}$.

This study represents the numerical analysis of natural convection of a microenvironments with a air permeability in the clothing air-layer. The clothing air layer of shoulder and arm was used for numerical analysis model. As a numerical analysis method, we adopted a finite volume method for two-dimensional laminar flow, and analyzed the flow and thermal characteristics of velocity, temperature and concentration in the air layer between body and clothing. As a temperature boundary conditions, we considered that a body skin has a high temperature with $34^{\circ}C$ the environmental temperatures are $5,\;15\;and\;25^{\circ}C$ for various permeability coefficients. The distributions of concentration, temperature and velocity were showed that two large cells were. formed at horizontal and vertical air layer, respectively. As the temperature difference between body skin and environment decrease, the heat transfer was decreased rapidly.

A new non-contact method of the frost thickness measurement has been developed. The method is based on the digital image processing technique to identify the reflection edge of the image captured by a CCD camera under laser sheet light illumination. To insure the accuracy of frost layer thickness, an in-situ calibration procedure is carried out with a calibration target with 0.5mm holes. Using the mapping function obtained by the calibration procedure, the contour of frost surface can be estimated with sub-pixel resolutions. The developed method is applied to study the effect of cooling plate temperature on the frost thickness in a small low speed wind tunnel.

Ringbom Stirling engine which is a kind of low temperature difference model Stirling engine is manufactured and its characteristics are measured at some temperature differences. Pressure, displacer position and rotation speed are measured. Displacer position and rotation speed are detected by photo-sensor. The hot side of Ringbom Stirling engine is warmed by electric heater. The cold side of Ringbom Stirling engine is cooled by the air. This result may be useful for further design and manufacture of Ringbom Stirling engine. Also, it would be used as an educational material for mechanical engineering students.

The purpose of this study is to develop a new type dryer, which is to proceed mixing and drying of wet-materials at the same time and drying process is carried out in a closed system. In this drying system, thermal contact occurs, when a fluidized zone is created by paddle mechanism. Accordingly, wet-materials is dried in a short time without any damage. Also wet-materials could be dried uniformly to low moisture contents. It is suitable to dry a small quantity of multi-kind materials. And drying process is carried out in a closed system, so as for environmental pollution dust not to be emitted into the atmosphere. Accordingly, it could be used to dry not only food and chemical materials, but also environmental pollution materials.

This study are development results of solar hot water system for anchovy proceeding. The heat amounts of boiling vessel are calculated 292.66W at forward and backward direction, and surface direction are calculated 373.14W. The polyenoic rate of anchovy are measured lower as high temperature, but monoenoic and polyacid are higher. Then the others. The maximum radiation of anchovy fishing grounds are shown $350kcal/m^2h$ at pm. 13:00, Chungdo in CHINA. Distributions of Velocity and temperature in boiling vessel are calculation. Solar collector and boiling vessel for anchovy proceeding ship are developed to automatic control system.

Catalytic combustion is the environmental-friendly technology, which has been applied to a variety of areas for industrial and domestic use in recent years. Accordingly, this study performed the development of the catalytic manufacturing technology for the high temperature and of the catalytic combustor in priority, which were aimed to be applied to a commercialized boiler. Paliadium(Pd) of a noble metal was used as a catalyst for the high temperature and supported on alumina($Al_[2}O_{3}$) and zirconia($ZrO_{2}$) in constant weight ratio. Activity of Pd catalysts is compared and analyzed in the catalytic combustion of natural gas. The ratio of $Pd/Al_{2}O_{3}=4$ was found to be better than any other weight ratios in activity and durability. The performance examination of catalysts and of combustion through the plate-type combustor made it possible to be developed the cylindrical-type combustor which has increased combustion area. Catalytic combustion boiler of 25,000 kcal/hr class was also developed, which had the optimum combustion condition at the nozzle of 5.95mm and the orifice of 21mm. This condition was determined through the performance experiments of catalytic combustion boiler to which the cylindrical-type catalytic combustor was applied.

Cryogenic LNG(Liquefied Natural Gas) which is stored in the cylindrical storage tanks of $100,000m^{3}$ has very complex flow phenomena and the changes of thermal properties with exterior conditions and operation mdoes. These complex thermofluid behaviors are affected by the storage, exterior conditions of LNG, design specifications and heat transfer characteristics of tanks, Also, those have influence on the stable storage and supply of LNG in the storage tanks. Thus this study performed the analysis on the 2-D heat transfer of the tank with exterior conditions, on the Cool Down Process in order to cool down the LNG Storage Tank at the initial normal state, and on the Filling Process considered for incoming and rising of LNG. The analysis on the Mixing LNG Storage was studied too. At last, the visualized program on the complex thermofluidodynamic analysis was developed on the basis of the above analyses. The development of this program becomes to be used to the basic design of the commercial tanks as well as to assure technical skill of the analysis on the thermal stability of the stored LNG in the LNG Storage Tank.

Exergetic and thermoeconomic analysis were performed for a 200kW Phosphoric Acid Fuel Cell(PAFC) plant which offers many advantage for cogeneration in the aspect of high electrical efficiency and low emission. This analytical study was based on the data obtained by in-field measurement of PC25 fuel cell plant to find whether this system is viable economically. For 100% load condition, the electrical efficiency and the unit cost of electricity are about 45% and 0.032 $/kWh respectively, which turn out to be much better than those for the 1000kW gas turbine cogeneration plant. Further, at lower loads, the unit costs of electricity and hot water increase slightly and consequently more economic operation is possible at any loads.

Numerical simulations of 1,2 dichloroethane(EDC) pyrolysis are conducted to understand the process on the production of the vinyl chloride monomer(VCM) and by-products. A chemical kinetic mechanism is developed, the adopted scheme involving 44 gas-phase species and 260 elementary forward and backward reactions. Detailed sensitivity analyses and the rates of production analysis are performed on each of the reactions and the various species, respectively. The concentrations of EDC, VCM, and HCI predicted by this mechanism are in good agreement with those deduced from experiments of commercial and laboratory scale. The mechanism is found to accurately predict the VCM yield and the production of by-products by varying the ranges of pyrolysis temperature, residence time, and pressure which impact on the pyrolysis of 1,2 dichloroethane. The influence of reactions related to H atom on the relative sensitivity of EDC becomes important as the residence time increases. The pyrolysis of EDC mainly occurs through $C_{2}H_{4}Cl_{2}+Cl=CH_{2}ClCHCl$.

A performance analysis of a SOFC/MGT hybrid system has been carried out for concept design. Thermo-dynamic models for each component being able to describe electrochemical characteristics and heat and mate-rial balance are proposed. Estimated is the power capacity of a SOFC suitable for the hybrid operation with a 5kW class MGT. Effects of current density and operating pressure are also investigated. Electric efficiency showed weak dependence on operating pressure and current density. It is desirable that the SOFC operates at high current density in manufacturing cost's point of view though operating with high current density slightly decreases the electric efficiency find specific power.

Although LP gas has lots of advantages, there has been limitation in application for automotive engine due to the several disadvantages, such as power decrease, complex fuel supply unit, and back fire etc. However LP gas direct injection engine has possibility to solve the problems above mentioned. LEM(Light Extinction Method) was employed for analysis of spacial and temporal distribution of LP gas which is directly injected into combustion chamber under various pressure and temperature conditions. The results from CVC(Constant Volume Chamber) were compared to those of RICEM(Rapid Induction, Compression and Expansion Machine) which simulate early- and late injection of direct injection engine. LPG fuel spray is affected by temperature and pressure in evaporation characteristics but it is more benefit to direct injection engine in every way such as, fuel distribution, evaporating speed and well wetting reduction.

The main reason for applying positive pressure variable clearance packing in fossil power plant is high efficiency and energy saving movement in the government. This study intends to analyze the turbine efficiency through the shaft packing improvement in thermal power plant and makes its comparison to that of the used packing.

The main reason for applying pulverizer classifier-rotary type in fossil power plant is boiler high efficiency and energy saving movement in the government. This study intends to analyze the boiler efficiency through the pulverizer classifier-rotary type improvement in thermal power plant and makes its comparison to that of the used fixed type.

A thermal design software is developed for the heat recovery steam generator(HRSG) of combined cogeneration systems. The heat transfer is calculated by using the element method to account for the varying thermal properties across the heat transfer elements. The circulation balance is computed for the evaporator to accurately estimate the steam generation rate and to check the proper circulation of the boiler water through the tubes. The software developed can be used to simulate HRSG systems with various combinations of auxiliary burner, wall superheater, superheater, reheater, evaporator, and economizer. Systems with several different combinations of the system components are successfully tested. And it is concluded that the developed software can be used for the design of heat recovery steam generators with various combinations of heat transfer components.

The mathematical method was developed and numerical analyses were carried out with various parameters to provide substantial data for optimal design and operation of urban utility systems. The composition of systems and their specifications, such as co-generation system, heat pump system, incineration system and other heating and cooling system could be obtained through these analyses for various resource and energy requirements in urban area. As results the system constituents and operating characteristics, and their economic performances such as the value of objective function, initial and an operating costs were discussed for various load patterns. The effective system design method and the excepted effects of the several unused energy recovery systems were also briefly discussed with the variation of the buildings and facilities species and their capacities.

A computer program, capable of performing thermal design analysis of the triple pressure bottoming system of combined cycle power plants, was developed. The program is based on thermal analysis of the heat recovery steam generator and estimation of its size and steam turbine power. The program is applicable to various parametric analyses including optimized design calculation. This paper presents examples of analysis results for the effects of arrangement of heat exchanger units, steam pressures and deaerating sources on design performance indices such as steam turbine power and the size of heat recovery steam generator.

In now a days, the concern to environment and energy saving problem is increased worldly. So many countries are developing the wind power system as clean energy system. In our country, Cheju local government has the plan of the Cheju Island wind farm and 600kW class 2 wind turbines, 660kW class 2 turbines, 225kW class 1 turbine and 750kW class 2 turbines has been operated at Hangwon. In this paper the field operation data of the wind turbines was analyzed and was compared with the characteristics & performance of each turbines. As the results, we would find the possibility of wind turbine in domestic and suggest the direction of developing technology.

Exergy concept is applied to the analysis of part-load performance of gas turbine engine. Exergy is a useful tool to find the source of irreversibility in thermal system. In this study, details of the performance characteristics of a heavy-duty gas turbine, l50MW-class GE 7FA model, are described by theoretical investigations with exergy analysis. Result shows that exergy destruction rate of gas turbine increases with decreased load, which means increase of irreversibility. Also, it is found that variations of IGV angle and amount of cooling air for turbine blades are closely related to the inefficiencies of compressor and turbine, respectively.

Heat (mass) transfer characteristics have been investigated on the endwall of a large-scale linear turbine cascade passage under a combustor-level high free-stream turbulence with a large length scale. Local heat (mass) transfer coefficients are measured by using the naphthalene sublimation technique. The result shows that local heat (mass) transfer on the endwall is greatly enhanced in the central region of the turbine passage, but there is no noticeable change in the local heat (mass) transfer in the region suffering severe heat load. Under the high free-stream turbulence, the local heat (mass) transfer coefficient shows more uniform distribution and its average value across the whole endwall region is increased by 26% of that at low turbulence condition. The heat (mass) transfer data on the endwall strongly supports that well-organized vortices near the endwall tends to suffer an suppression by the high free-stream turbulence.

The objective of the present study is to investigate the effects of the various inlet boundary layer thickness on convective heat transfer distribution in a turbine cascade endwall and blade suction surface. In addition, the proper height of the boundary layer fences for various inlet boundary layer thickness were applied to turbine cascade endwall in order to reduce the secondary flow, and to verify its influence on the heat transfer process within the turbine cascade. Convective heat transfer distributions on the experimental regions were measured by the image processing system. The results show that heat transfer coefficients on the blade suction surface were increased with an augmentation of inlet boundary layer thickness. However, in a turbine cascade endwall, magnitude of heat transfer coefficients did not change with variation of inlet boundary layer thickness. The results also present that the boundary layer fence is effective in reducing heat transfer on the suction surface. On the other hand, in the endwall region, boundary layer fence brought about the subsidiary heat transfer increment.

By matching a well established fast throughflow code, with standard loss correlations, and an efficient optimization algorithm, a new design system has been developed, which optimizes inlet and exit flow-field parameters for each blade row of a multistage axial flow turbine. The compressible steady state inviscid throughflow code based on streamline curvature method is suitable for fast and accurate flow calculation and performance prediction of a multistage axial flow turbine. A general purpose hybrid constrained optimization package, iSIGHT has been used, which includes the following modules: genetic algorithm, simulated annealing, modified method of feasible directions. The design system has been demonstrated using an example of a 5-stage low pressure steam turbine for 800MW thermal power plant previously designed by HANJUNG. The comparison of computed performance of initial and optimized design shows significant improvement in the turbine efficiency.

An Altitude Engine Test Facility(AETF) was built at the Korea Aerospace Research Institute in October 1999 and has been being operated for altitude testing of the gas turbine engines of 3,000 Ibf class or less. The AETF has been calibrated using several engines such as J69 engine of Teledyne Co. as a facility checkout engine. Based on the test results, uncertainty analyses on the air flow rate and thrust were performed according to ASME PTC 19.1-1998. As the analyses showed that the level of uncertainty was not satisfactory over the whole operating envelop, several modifications of the facility and testing method were made in order to improve the measurement uncertainty. As a result, the uncertainty of the air flow measurement was improved by 0.1 % over all the test conditions, and the net thrust measurement by upto 3%. The improved measurement uncertainties of air flow and thrust are 0.68-0.73% and 0.4-1.3%, respectively.

Local heat/mass transfer and friction loss in a square duct roughened with various types of continuous and discrete rib turbulators are investigated. The combined effects of the gap flows of the discrete ribs and the secondary flows are examined for the purpose of the reduction of thermally weak regions and the promotion of the uniformity of heat/mass transfer distributions as well as the augmentation of average heat/mass transfer. The rib-to-rib pitch to the rib height ratio (p/e) of 8 and the rib angles of 90 and 60 deg are selected with $e/D_{h}=0.08$. The vortical structure of the secondary flows induced by the parallel angled arrays are quite distinct from that induced by the cross angled arrays. This distinction influences on heat/mass transfer and friction loss in all the tested cases. The gap flows of the discrete ribs reduce the strength of the secondary flows but promote local turbulence and flow mixing. As a result, the fairly uniform heat/mass transfer distributions are obtained with two row gaps.

Since the characteristics of combustion and pollutant in Diesel engines were mainly affected by the characteristics of in-cylinder gas flow and fuel spray, an understanding of those was essential to the design of the D.I. Diesel engines. The improvement of volumetric efficiency of air charging into combustion chamber is a primary requirement to obtain better mean effective pressure of an engine. Since the air resistances in intake and exhaust flow passages, valve lift and valve shape influence greatly to the volumetric efficiency, it is very important to investigate the flow characteristics of intake and exhaust port which develops air motion in the combustion chamber. This paper presents the results of an experimental investigation of steady flow through the various kinds of commercial cylinder head ports, and the development procedures of HHI's H21/32 prototype cylinder head ports.

During cold operation period, fuel injection system directly contributes the unburned hydrocarbon formation in spark ignition engines. The relationship between injection parameters and HC emission behavior was investigated through a series of experiments. Spray behavior of port fuel injectors was characterized through a quantitative evaluation of mass concentration of liquid fuel by a patternator and PDA. 6-hole injector was found to produce finer spray than single hole one. Using a purpose-built test rig, the wall wetting fuel was measured, which was mostly affected by wall temperature. Varying coolant temperature($20{\sim}80^{\circ}C$), HC emissions were measured in a production engine. With respect to the different types of injectors, HC emission was also measured. In the 6-hole injector application, the engine produced less HC emission in low coolant temperature region. Though it produces much more amount of wetting fuel, it has the advantages of finer atomization quality. In high coolant temperature region, there was little effect between different types of injectors. The control schemes to reduce HC emissions during cold start could be suggested from the findings that the amount of fuel supply and HC emission could be reduced by utilizing fine spray and high intake wall temperature.

Diesel particulate trap is the most reliable system to reduce the particulate matters from diesel engine. Filter is the core component of DPF and ceramic monolith type is dominantly used, which is expensive and fragile relatively at thermal shock. Porous metal filter, which has superior thermal characteristics and low cost, was tested in order to analyze the regeneration performance by using with ferrocene additive. This filter showed the 72% filtration efficiency, additives itself diminished 48% of PM from engine out emission, and final PM reduction ratio of 89% was achieved by DPF system with D-13 test mode.

This paper represents the wall thinning analysis results for secondary side piping of two types of domestic nuclear power plants based on the DB establishment and F AC analysis study for NPP secondary system piping. CHECWORKS code utilized in this study has been applied world widely to wall thinning analyses for secondary side piping and its reliability has also been proved. The predicted wear rates for several piping systems of a pressurized water reactor NPP are compared with those of a pressurized heavy water reactor NPP and with the measured wear rates. On the basis of comparison results of the predicted and measured wear rates, the analysis results can be effectively applied to the development of a standard thinned pipe management program targeted all domestic nuclear power plants.

The Safety depressurization System(SDS) of KNGR prevents RCS from overpressurization by discharging high pressure and temperature coolant through the I-sparger into the IRWST during an accident. If IRWST water temperature rise locally, around the sparger, beyond $200_{\circ}$2000 F by the discharged coolant, unstable steam condensation can cause large pressure load on the IRWST wall. To investigate whether this condition can be avoided for the design basis event IOPOSRV(Inadvertent Opening of one Pilot Operated Safety Relief Valve), the flow and temperature distribution of water in the IRWST is calculated by using CFX 4.3 computational fluid dynamic code. According to the results, since pool water temperature does not exceeds temperature limit within 50 seconds after the opening of one POSRV, it can be assured that the integrity of IRWST wall is maintained.

As one of the advanced design features of the Korea next generation reactor, direct vessel injection (DVI) system is being considered instead of conventional cold leg injection (CLl) system. It is known that the DVI system greatly enhances the reliability of the emergency core cooling (ECC) system. However, there is still a dispute on its performance in terms of water delivery to the reactor core during the reflood period of a large-break loss-of-coolant accident (LOCA). Thus, experimental validation is under progress. In this paper, a new scaling method, using time and velocity reduced linear scaling law, is suggested for the design of a scaled-down experimental facility to investigate the direct ECC bypass phenomena in PWR downcomer.

In this study, an analysis tool to assess the susceptibility of steam generator tubes due to the flow-induced vibration was developed. The fluid-elastic instability analysis of the U-tube bundle for CE-type steam generator was accomplished. The effective mass distribution along the U-tube was obtained to calculate the natural frequency and dynamic mode shape. Finally, stability ratios for selected tubes are obtained.

Calandria tube wrapping each pressure tube is one of the key structural components of CANDU reactor(Calandria) which is consisted of many pressure tubes containing nuclear fuel assemblies. As the Calandria tube(made of zirconium alloy) is sagging due to its thermal and irradiation creep during the plant operation, it possibly contacts with liquid injection nozzle crossing beneath the Calandria tube, which subsequently results in difficulties on the safe operation. It is therefore necessary to check the gap for the confirmation of no contacts between the two tubes, Calandria tube and liquid injection tube, with a proper measure during the life of plant. In this study, an ultrasonic measurement method was selected among several methods investigated. The ultrasonic device being developed for the measurement of the gap was introduced and its preliminary performance test results were presented here. The gap between LIN and CT at site was measured using by this ultrasonic device at site.

Development of a small gas-turbine combustor for 100kW class APU(Auxiliary Power Unit) has been performed. This combustor is a reverse-annular type and has a tangential swirler in the liner head to improve the fuel/air mixing and flame stability. Three main and three pilot fuel injectors of the simplex pressure-swirl type are used. The performance target at the design condition includes a turbine inlet temperature of 1170K, a combustion efficiency of 99%, a pattern factor of 30%, and an engine durability of 3000 hours. Under developing the combustor, we conducted performance test of our first prototype(TS1) with some variants. As a result of the test, the performance targets of the combustor are satisfied except that the pattern factor is about 4% higher than target value. So, we redesigned the second prototype(TS2) and conduct performance test with the critical focus on pattern factor and exit mean temperature. We adopted TS2 four variant to check the improvement of pattern factor. As the result, the pattern factors of several variants were satisfied with the performance target. Finally, We chose the TS2A variant as a final combustor for our APU model.

For the successful development of the main engine of KSR(Korea Sounding Rocket)-III, Korea Aerospace Research Institute(KARI) carried out the experimental study on the subscale model engines. Several types of engines were tested on the Small Liquid Rocket Engine Test Facility. One of the typical test results of a Sub. engine(Sub. Mod.3) is presented here. It uses the Jet A-1 as fuel, liquid oxygen as oxidizer, and Tri-Ethyl Aluminium(TEA1) as ignition agent. The gas pressure feed system is adopted as a feeding mechanism and the design chamber pressure is 200psia. The physical phenomena are described in three regimes(ignition, transient, and steady state) with the pressure, thrust, flowrate and image data. And the pressure oscillation is analyzed in Fourier domain (<500Hz). Then we conclude that in this experiment, the engine shows the characteristic low frequency of 80Hz and it is stable for that frequency of pressure oscillation.

This paper discusses the combustion characteristics of water emulsified fuel produced by swirl impinging mixer. The experiments are given in a test boiler and a commercial boiler. Flame temperatures, flame shapes and exhaust emissions are measured in the test boiler, and thermal efficiency is tested in the commercial boiler. The test data show the water in oil makes the flame wider and shorter, the flame temperature lower and the NOx and smoke lower. Also, the commercial boiler test shows the thermal efficiency increases as well as the emissions decrease.

In this paper, the test and result of flow and combustion for 21AFR lean fuel models are described. The necessity to develop the low emission combustor has been issued from the concern on the increase of green house and the destruction of ozone layer. To evaluate the flow and combustion performance of new designed 21AFR lean modules, the hydraulic tests in stereo lithographic airflows models, the low pressure combustion tests in three injectors model for weak extinction and ignition and the high pressure combustion tests in single sector for NOx, SAE and efficiency are performed. The low pressure tests reveal that the governing parameters in weak extinction and ignition at atmospheric condition are prefilmer length, swirl flow rotation direction, secondary swirl angle and flow split. As a results of combustion test at high pressure, the efficiency and smoke level are satisfied with performance targets, but EINOx of 17.8 is higher than target value of 13.1. The high pressure tests show that the main parameters influenced on NOx are primary swirl angle, swirl flow rotation direction, heatshield exit angle and liner mixing hole location.

Flame propagation along vortex tube was experimentally investigated. The vortex tube was generated by the ejection of propane from a nozzle through a single stroke motion of a speaker and the ignition was induced from a single pulse laser. Non-reactive flow fields were visualized using shadow technique. From these images, vortex ring size and translational velocity were measured in order to determine the ignition time and position. Flame structure and flame speed were measured using high speed CCD camera. Flame speed was accelerated during the initial stage of flame kernel growth, and reached near constant value during steady propagation period. Near the completion of propagation, flame speed was decelerated and then extinguished. Flame speed along the non-premixed vortex tube was found to be linearly proportional to circulation, which was similar to that of the flame propagation along premixed vortex ring. Ignition position minimally affects the propagation characteristics. These imply that flame is propagating along the maximum speed locus expected to be along stoichiometric contour and also support the existence of tribrachial flames.

Experimental study was conducted to investigate the effects of axial forcing on the flame structures near the nozzle exit of non-premixed flame. The most notable observation is that the direction of vortical motions is changed at some ranges, according to the increase of excitation amplitude. Especially, the elongation flame and the phenomenon of In-burning are always occurred when the vortical motion turnabout. In an analysis of the flame/flow visualization by means of direct photography and RMS technique, a plausible explanation can be made that above phenomena are related only to the amplitude of phase average velocity between the instantaneous velocity elements of excited flow.

The effects of external excitation with various frequencies and amplitudes on the flame behavior and pollution emission characteristics from a laminar jet flame are experimentally investigated. Measurements of $NO_{x}$ emission indices($EINO_{x}$), performed in unconformed and vertical lifted flame at resonance frequency by strong excitation, have been conducted experimentally. It was also conducted to investigate the effects of excited frequency on $NO_{x}$ emissions with a various frequency ranged 0Hz to 2kHz. From the vertical lifted flame like turbulent of the excited jet with resonance frequency was shown that the dependence of $NO_{x}$ emission could be categorized into three groups: Group 1 of intermediate flame length and relative narrow flame volume yielding low $NO_{x}$ emission, Group 2 of short flame length but large flame volume yielding high $NO_{x}$ emission and Group 3 of long flame length with low temperature contours yielding high $NO_{x}$ emission.

Combustion chamber crevices in SI engines are identified as the largest contributor to the engine-out hydrocarbon emissions. The largest of crevice region is the piston ring pack crevice. To predict and understand the oxidation process of piston crevice hydrocarbons, a 3-dimensional numerical simulation method was developed. A engine shaped computational mesh with moving grid for piston and valve motions was constructed. And a 4-step oxidation model involving 7 species was used and the 16 coefficients in the rate expressions were optimized based on the results from a detailed chemical kinetic mechanism for the oxidation condition of engine combustion chamber. Propane was used as a fuel in order to eliminate oil layer absorption and liquid fuel effect.

Pressure resonance frequency that is caused in the combustion chamber can be interpreted to acoustic analysis. Until now the pressure resonance has been assumed and calculated to a disc type combustion chamber that neglected the combustion chamber height because the knock occurs near the TDC(top dead center). In this research FEM(fine element method) has been used to calculate the pressure resonance frequency inside the experimental engine combustion. The reduce error of the resonance frequency obtained by FEM has decreased about 50% compared to the calculation of Draper's equation. Due to the asymmetry in the shape of the combustion chamber that was neglected in Draper's equation we could find out that a new resonance frequency could be generated. To make the experimental results equal we could know that the speed of sound that satisfies Draper's equation was selected 13% higher than all the pent-roof type combustion considered.

To meet strict emission regulation while improving engine performances, common rail injection system which is suitable for electronic control, and capable of controlling injection quantity, timing, rate and pressure individually as well as realizing high pressure has been developed. At present study, a 8L DI diesel engine was converted to a single-cylinder experimental engine allowing optical access through an extended piston and a prototype of common rail injector in progress was applied to the engine. The combustion characteristics of the engine were analysed by using direct images and characteristics of the injector were analysed. We can not say that the results are always the same to general common rail injection system but that they are just characteristics of specific prototype injector.

A web-based client/server program, MOVIDIK(Motor-Operated-Valve Integrated Database Information of KEPCO) has been developed to perform a design basis safety evaluation for a motor-operated-valve(MOV) in the nuclear power plant. The MOVIDIK consists of seven analysis modules and one administrative module. The analysis module calculates a differential pressure on the valve disk, thrust/torque acting at a valve stem, maximum allowable stress, thermal-overload-relay selection, voltage degradation, actuator output and margin. In addition, the administrative module manages user information, approval system and code information. MOVIDIK controls a huge amount of evaluation data and piles up the safety information of safety-related MOV. The MOVIDIK will improve the efficiency of safety evaluation work and standardize the analysis process for the MOV.

Fluid properties which are most commonly used to evaluate spray atomization characteristics, are important because they affect velocity and size distribution of droplets. The purpose of this study was to incorporate the significant characteristics in atomization process of industrial etching spray and how each of them affects the design of precise pressure-swirl nozzles. The experiment was carried out with different viscosity and density of fluid. The macro characteristics of liquid spray, such as the spray angle and shape were captured by PMAS and the micro characteristics of liquid spray, such as droplet size and velocity were obtained by PDA. The mean velocity and SMD of droplets were measured along axial and radial direction. It was found that the higher viscosity and density resulted in the larger SMD and the lower mean velocity of droplets.

A spent fuel problem has prevented the nuclear power from claiming to be a completely clean energy source. The nuclear transmutation technology to incinerate the long lived radioactive nuclides and produce energy during the incineration process is believed to be one or the best solutions. HYPER(Hybrid Power Extraction Reactor) is the accelerator driven transmutation system which is being developed by KAERI(Korea Atomic Energy Research Institute). Some major feature of HYPER have been developed and employed. On-power fueling concepts are employed to keep system power constant with minimum variation of accelerator power. A hollow cylinder-type metal fuel is designed for the on-line refueling concept. Lead-bismuth(Pb-Bi) is adopted as a coolant and Spallation target material. HYPER is a subcritical reactor which needs an external neutron source. 1GeV proton beam is irradiated to Lead-bismuth(Pb-Bi) target inside HYPER, and spallation neutrons are produced. When proton beams are irradiated, much heat is also deposited in the Pb-Bi target and beam window which separates Pb-Bi and accelerator vacuum. Therfore, an effective cooling is needed for HYPER target. In this paper, we performed the thermal-hydraulic analysis of HYPER target using FLUENT code, and also calculated thermal and mechanical stress of the beam window using ANSYS code.

The reliability and accuracy of the information on control rod position are very important to the reactor safety and the design of the core protection system. In this study, a thorough investigation on the RSPT(Reed Switch Position Transmitter) type control rod position indication system and its actual implementation in the exiting nuclear power plants in Korea was performed first. A design of the control rod position indication system using reed switch for the CEDM on the system-integrated reactor SMART was developed based on the position indicator technology identified through the investigation. The feasibility of the design was evaluated by test of manufactured control rod position indicator using reed switch for SMART.

It is intended to develope an algorithm for dynamic simulation of fast-acting solenoid valves. The coupled equations of the electric, magnetic, and mechanical systems should be solved simultaneously in a transient nonlinear manner. The transient nonlinear electromagnetic field is analyzed by the Finite Element Method (FEM), which is coupled with nonlinear electronic circuitry. The dynamic movement of the solenoid valve is analyzed at every time step from the force balances acting on the plunger, which include the electromagnetic force calculated from the Finite Element analysis as well as the elastic force by a spring and the hydrodynamic pressure force along the flow passage. Dynamic responses of the solenoid valves predicted by this algorithm agree well with the experimental results including bouncing effects.

A 5-valve(intake 3-valve) engine has been developed to increase engine performance. These vehicles have a high power caused by the decrease of inertia mass of an intake valve and the increase of intake effective area. In this study, velocity profiles at near intake valves were inspected by using a two-color PIV and laser sheet method with tumble control valve(TCV). In addition, steady flow tests were performed to quantify tumble ratio on flow-fields generated with a TCV. These experimental results show that the tendency of the tunble ratio in intake 3-valve engine is different from the one in intake 2-valve engine. From this results, the intake flow characteristics around intake valves were made clear.

For the application of simultaneous $DeSO_{2}\;&\;DeNO_{x}$ equipment using non-thermal plasma process to the industrial and power plants, the many types of plasma device and process were studied. The e-beam and pulsed plasma corona discharge process are outstanding for the study to apply commercial large-scale plant from among these. In this paper, non-thermal plasma of technical trends and the characteristics of system developed by Doosan heavy industries & construction Co., Ltd. are explained. We have researched pulsed plasma corona discharge process since 1994. At the basis of reasonable results for the pilot plant, we constructed the demonstration plant at a domestic coal-fired power plant in 1999, as the previous step for commercial use. In near future, enough information about designs and costs of commercial-size system will be obtained.

This study is an experimental study on the characteristics of emission by changing catalytic converter position for cold-start. The measurements are done a changing of the distance between exhaust manifold and catalytic converter. It measured temperature of exhaust manifold, before and after catalytic converter at each position of experimental condition. and measured the characteristics of emission which is HC, CO, $CO_{2}$ and lambda at each position of experimental condition. The results show a few advantage about reduction of HC and CO as catalytic converter's temperature is raised quickly as closed exhaust manifold. but $CO_{2}$ has not the same trend of HC and CO. From measurement value of lambda, reduction effects of $NO_{x}$ are known a few advantage as increase of the distance between exhaust manifold and catalytic converter.

Gasoline engine manufacturers are currently considering designs that will result in low combustion air temperature for improvement of fuel consumption and emission levels. There are a variety of cooling systems that can be used to accomplish this goal. Coolong is therefore normally achieved through a balance of ram and fan action. This paper studies the various systems and compare the cooling performance for several conditions, based on a automotive engine. An experimental analysis was developed to predict the interaction of the fan system and the heat exchangers of the engine cooling system. The local temperature induced by the fan on the cooling system is measured. These experimental result were accomplished using air flow management techniques.

In order to develop an efficient turbine for wave energy conversion suitable for actual ocean conditions, a new type of the air turbine with staggered blades has been investigated experimentally. Experiments have been carried out under steady flow conditions. Both the running and starting characteristics under sinusoidally oscillating flow conditions are simulated by a CFD method using a quasi-steady analysis. It is known that the air turbine with staggered blades gives a better performance compared with conventional Wells turbine, and a proper design factor of the air turbine is discussed for the setting angle of the rotor.

The Wells turbine is one of the simplest and most promising self-rectifying air turbines which are useful for the systems of alternative energy development in near future, and it is economically desirable from the point of view of the practical use, as well. To investigate the effect of blade sweep on the performance of the Wells turbine, computations of a fully 3-D Navier-Stokes are carried out under steady flow conditions of NACA0020 blade. It is known that the performance of the Wells turbine is considerably influenced by the blade sweep. An optimum blade sweep ratio(f=0.35) for the NACA0020 is found to be the most promising for the practical use, and this value is in good agreement with the previous experiments. It is also found that the overall turbine performance for the NACA0020 is better than that for the CA9.

The emission of $NO_{x}$ during coal combustion is a major reason of environment impact. $NO_{x}$ is an acid rain precursor and participates in the generation of smog through ozone production. $NO_{x}$ can be divided into thermal $NO_{x}$, fuel $NO_{x}$ and prompt $NO_{x}$. Thermal $NO_{x}$ is formed in a highly temperature condition dependent. Fuel $NO_{x}$ is dependent on the local combustion characteristics and initial concentration of nitrogen bound compound, while prompt $NO_{x}$ is formed in a significant quantity in some combustion environments, such as low temperature and short residence times. This paper presents numerical simulation of the flow and combustion characteristics in the furnace of a tangential firing boiler of 500MW with burners installed at the every comer of the furnace. The purpose of this paper is to investigate the reduction of $NO_{x}$ emission in a 500MW pulverized coal tangential firing boiler with different OFA's and burner angles. Calculations with different air flow rates of over fired air(OFA) and burner angles are performed.

Numerical Analysis for liquid fuel combustion of horizontal firing boiler is performed. The mixture-fraction/PDF equilibrium chemistry model is used to predict the combustion of the vaporized fuel. P1 model for radiation effect is used. Superheater, reheater and economizer is modeled using porous with heat sink. Flow and temperature field is investigated, and distribution of thermal $NO_{x}$ and CO is investigated. Computation as the change of excess air and swirling is performed to investigate the change of thermal $NO_{x}$.

As a preceding process for developing design technology and establishing operation technology, the design procedure of the SCR(Selective Catalytic Reduction) pilot plant that can handle $1,000Nm^{3}/hr$ of flue gas was reported in this paper. And we also considered several factors that might cause abnormality of the plant in the designing process. The plant was designed and fabricated to test the $DeNO_{x}$ performances in variable operating conditions in the range of $3,000{\sim}36,000hr^{-1}/hr$ in space velocities, $1.67{\sim}6\;m/s$ in linear velocities, $200{\sim}500^{\circ}C$ temperatures, $300{\sim}1,000Nm^{3}/hr$ flow rates, and $0{\sim}1.4:1\;NH_{3}/NO$ ratios. In order to maintain the flow uniformity, the guide vanes and flow straightener were designed and constructed in the plant. The SCR pilot plant can be operated by the automatic control system, which enable to obtain performance data in real time and to set up the operating technology. The catalyst reactor consists of 4 catalyst layers and surface area of each layer can be adjusted to be of small size. Arrangement of catalysts per layer is $3{\times}6$ with the catalyst dimensions of $150{\times}150{\times}500mm(L{\times}W{\times}H)$.

Revision and establish of KS(Korean Standards) is currently actively discussed. It is just the time for a new world class standards under the new system with WTO(World Trade Organization). This paper is a part of "Researchs on the Standards in Railway rolling stock Field", as one of KS establish projects. The aim of this study is to define the requirements of railway rolling stock-test methods for running resistance. In the former KS, there is no items matched with this purpose. Therefore a new part of KS is proposed.

Numerical simulations of high temperature catalytic combustion have been performed for the application to a gas turbine combustor. Dependences of inlet temperature and pressure on the distributions of temperature and species concentrations were investigated using plug flow model with detailed homogeneous and heterogeneous chemistries of methane-air mixtures. Honeycomb typecombustor deposited with Pt catalyst of 100mm in length and 26mm in diameter is used. The results show that rapid increase of temperature profile occurs earlier with the increase of inlet temperature and the decrease of inlet pressure. The condition which catalytic combustion is stabilized exists at certain range of inlet temperature and pressure. The state of catalytic combustion is also confirmed by the distributions of species concentration.

Distributions of the parameters in proton exchange membrane fuel cell (PEMFC) has been analyzed numerically under steady-state and isothermal conditions. The distributions of the crucial parameters (e.g., temperature and pressure) in a PEMFC have a major impact on its safe and efficient operation. This paper predicts the performance of the model electrode plates by calculating the pressure and temperature distributions of working fluid. The calculated results of pressure and temperature at exit condition shows good agreement to experiments and gives details of flow pattern inside of electrode plates.

Hydrogen fueled engine with dual injection can achieve high power and high efficiency simultaneously. In this study, the suitable compression ratio of hydrogen fueled engine with dual injection were investigated including performance of this engine according to variation of compression ratio. As results, it was found that the suitable compression ratio of that was about CR=11, and torque and thermal efficiency increased by 6% and 7% respectively.

KSR-III propulsion system designed in KARI has a gas-pressurization system for propellant feeding system. This system uses a regulator for the control of the ullage pressure of propellant tank and a venturi for passive control of propellant flowrate. This system seems to be very reliable, but the flowrate of propellant varies according to the change of acceleration with the rocket flight. In this paper, dynamic characteristic of KSR-III propulsion feeding system was analyzed in flight condition. The purpose of this research is to find the variation of off ratio and propellant flowrate change for certification condition of engine reliability test.

A numerical investigation was performed to determine the effect of airfoil on the optimum flap height using NACA 00XX and 44XX airfoils. The six flaps which have 0.5% chord height difference were selected. A Navier-Stokes code, FLUENT, was used to calculate the flow field of the airfoil. The code was first tested as a benchmark by modelling flow around a NACA 4412 airfoil. Predictions of local pressure coefficients are found to be in good agreement with the result of the experimental result. For every NACA 00XX and 44XX airfoil, flap heights ranging from 0.0% to 2.5% chord were changed by 0.5% chord interval and their effects were also studied. Representative results from each case are presented graphically and discussed. It is concluded that this initial approach gives an idea for the future development of the wind turbine optimum design.