In this study, we compare thermal performance between four different types of cold plates for humanoid robot cooling. Two commercially available cold plates made of copper have different dimensions and internal flow paths: One has $20{\times}20$ $mm^2$ base area with micro-channels and the other has $62.5{\times}62.5$ $mm^2$ base area with 85 round pin-fins. And two different types of cold plates of $20{\times}20$ $mm^2$ base area with 7 mm high are made of PC (polycarbonate), which aims to reduce the weight of cooling system. All cold plates are mounted on a $20{\times}20$ $mm^2$ copper block with two cartridge heaters of 30 $W/cm^2$. The overall heat transfer coefficient and thermal resistances for the liquid-cooled cold plates are obtained. The copper cold plate with micro-channels showed the best performance. Polycarbonate cold plates display fairly good thermal performance with more reduced system weight.

The microchannel heat sink is promising heat dissipation method for high heat flux source. Contrary to conventional circular channel, MEMS based microchannel had rectangular or trapezoidal cross-sectional shape. In our study, we conducted three dimensional conjugate heat transfer calculation for rectangular shape microchannel. First, we simulated that channel was completely drained with known heating power. As a result we obtained calibration line, which indicates heat loss was function of temperature. Second, we simulated single phase heat transfer with various mass flux, 100-400 $kg/m^2s$. In conclusion, the single phase test verified that the present heat loss evaluation method is applicable to micro scale heat transfer devices. Heat fluxes from each side wall shows difference due to non-uniform heating. However those ratios were correlated with supplied total heat. Finally, we proposed effective area correction factor to evaluate appropriate heat flux.

The micro plated heat exchangers were designed to transfer more heat/volume or mass than previous heat exchangers within the context of the design constraints specified. The increase of the surface-to-volume ratio results in an increase of the interfacial area. This enhances considerably the performance of a heat exchanger. This can be an important component in a wide range of applications fuel cell, aerospace, automotive, electronic system and home heating, etc). In this study, the performance evaluation of micro plated heat exchangers under the counter flows with straight and S-shaped channel are carried out. The pressure drop as well as inlet and outlet fluid temperature were measured at steady state under various operating conditions and the total heat transfer rate were also calculated.

In the present experimental study, the effect of interfacial tensions on pressure drop of air-water two-phase flow in round mini-channels was investigated. A glass (highly wettable) tube and a Teflon (poorly wettable) tube, both in 350 mm length but 1.8 mm and 1.59 mm in inner diameters each, were used for the tests. All the experiments were performed only in the plug flow regime, confirmed by visualization. In the glass tube, the gas plugs were surrounded by the liquid film along the inner periphery. On the other hand, the inner wall remained dry at the gas portion in the Teflon tube. The pressure drop of the plug flow in the Teflon tube without the liquid film) appeared much larger than in the glass tube (with the liquid film) due to dissipation of energy by movement of the wetting lines. In this paper, various correlations on the two-phase pressure drop of plug flows were compared and a modified correlation was proposed, taking account of the surface wettability.

In this study, the microchannel plated heat exchanger were numerically studied for the enhancement of heat transfer in the channel configuration. Unit cold and hot fluid region with the microchannel were modeled and periodic boundary condition at the side wall was applied to continuously repeating geometry. The material of micro-structured plate is STS304 and working fluid is water. Triangular obstacles were placed in micro channel to enhance heat transfer. The performance of microchannel plated heat exchangers were numerically investigated with various obstacle configuration and Reynolds number under the parallel and counter flows. Heat transfer rate has increased about 18% compared with straight channel, but pressure drop also increased about 3.5 times. The main factor of increasing of pressure drop and heat transfer rate is considered that the momentum was lost to collide against obstacles, generation of secondary flow and boundary layer separation, wake and vortex forming phenomena.

This paper presents the thermal characteristics for micro heat exchanger with different micro-channel shapes. The shapes of micro-channel has been manufactured sheet metal by chemical etching for the I shape of straight channel and V and W shapes of chevron feature and fabricated micro plated heat exchangers using the vacuum brazing of bonding technology. The experimental study has been performed on heat transfer and pressure drop characteristics with various Reynolds number for water to water at the counter flows. The average heat transfer rate of V and W shapes has been showed about 1.5${\sim}$1.6 times large than those of I shape. For the comparison of Nusselt number, it is known that the convective heat transfer of V and W shapes represent more effect than I shape. The pressure drops of V and W shapes are about 1.2${\sim}$1.7 times lager than those of I shape.

A micro cyclone combustor was developed to be used as a heat source of thermoelectric power generator (TPG). The cyclone combustor was designed so that fuel and air were supplied to the combustion chamber separately. The mixing and flow characteristics in the combustor were investigated numerically. The global equivalence ratio (${\Phi}$), defined using the fuel and air flow rates, was introduced to examine the flow features of the combustor. The mixing of fuel and air inside the combustor could be well understood using the fuel concentration distribution. It was found that the weak recirculating zone was formed upper the fuel-supplying tube in case of ${\Phi}$ < 1.0. In addition, it was found that small regions that have a negative axial velocity exist near the fuel injection ports. It is assumed that these negative axial velocity regions can stabilize a flame inside the micro cyclone combustor.

It has been well known that pool boiling CHF in nanofluids compared to pure water significantly increase due to the deposition of nanoparticles on heater surface. This study concerns the characteristics of the nanoparticle deposition layer and its influence on CHF. Pool boiling experiments were carried out with 0.01vol.% water-$TiO_2$ nanofluids to obtain various nanoparticle-deposited heaters. CHF on the prepared heaters was measured during pool boiling in pure water. The heater surfaces were visualized using scanning electron microscope (SEM) and also characterized using contact angle and capillarity. The results showed that the CHF enhancement in nanofluids was completely dependent upon the structural and physicochemical characteristics of the nanoparticle deposition layer.

The disjoining pressure is critical in modeling the transport phenomena in small scales. They are very useful in characterizing the non-continuum effects that are not negligible in heat and mass transports in the film of less than submicro-scales. We present he disjoining pressure of thin film absorbed on solid substrate using Molecular Dynamics Simulation (MD). The disjoining pressure with respect to the film thickness is accurately calculated in the resolution of a molecular scale. The characteristics of the pressure are discussed regarding the molecular nature of the fluid system like molecular diameter and intermolecular interaction. Also, the MD results are compared with those based on the macroscopic approximation of the slab-like density profile. Significant discrepancy is observed when the effective film thickness is less than several molecular diameter

We have been setting up experiments on propagation of shock waves generated by the pulsed laser ablation. One side of a thin metal foil is subjected to laser ablation as a shock wave propagates through the foil. The shock wave, which penetrates through the foil is reflected by an acoustic impedance which causes the metal foil to high-strain rate deform. This short time physics is captured on an ICCD camera. The focus of our research is applying shock wave and deformation of the thin foil from the ablation to accelerating micro-particles to a very high speed.

A transformer is a device that changes the current and voltage by electricity induced between coil and core steel, and it is composed of metals and insulating materials. In the core of the transformer, the thermal load is generated by electric loss and the high temperature can make the break of insulating. So we must cool down the temperature of transformer by external radiators. According to cooling fan's usage, there are two cooling types, OA(Oil Natural Air Natural) and FA(Oil Natural Air Forced). For this study, we used Fluent 6.2 and analyzed the cooling characteristic of radiator. we calculated 1-fin of detail modeling that is similar to honeycomb structure and multi-fin(18-fin) calculation for OA and FA types. For the sensitivity study, we have different positions(side, under) of cooling fans for forced convection of FA type. The calculation results were compared with the measurement data which obtained from 135.45/69kV ultra transformer flowrate and temperature test. The aim of the study is to assess the Fluent code prediction on the radiator calculation and to use the data for optimizing transformer radiator design.

In this study, we performed the CFD analysis to predict the cooling performance of the motor for golf-car application. The motor for the golf-car is different from the general motors in the cooling method and the operating condition. So, it is important to predict the motor's temperature generated by the electric losses. We investigated the inner temperature rise using the unsteady state analysis and compared with the experiment for the reliability evaluation of the CFD results through the present study.

The application of high velocity air jets to heat transfer surface has tremendous engineering potential in various industries. The impinging jets are therefore widely used for their enhanced transport characteristics, especially for drying of continuous materials such as sheet, film, carpets, forming materials and pallets. This paper presents the drying characteristics of plate type material and performs the test with the change of operating conditions in conveyer band dryer using the impinging jet nozzles. The factors influencing drying rate were confirmed, also the design data of band dryer using the impinging jet was obtained.

During a process of a nanoimprint for manufacturing LCD, a small temperature variation on the LCD glass can cause thermal stress and generate unexpected displacement. To avoid this trouble, a precision temperature control unit using thermoelectric modules is appropriate for nanoimprint processes. The unit consists of an air control system, a cooling water control system, and a power control system. The air control system includes a thermoelectric module, thermocouples measuring temperatures of air and a duct-stale fin, and two air fans. The heat generated by the thermoelectric module is absorbed by the cooling water control system. The power control system catches the temperature of the thermoelectric module, and a PID controller with SCR controls the input power of the thermoelectric module. Temperature control performance is evaluated by experiment and simulation. The temperature control unit is able to control the exit temperature about ${\pm}2^{\circ}C$ from the incoming fluid temperature, and the error range is ${\pm}0.1^{\circ}C$. However, the control time is approximately 30minute, which needs further study of active control

A thermoelectric module can be used for cooling or power generation. The basic requirements to achieve a significant thermoelectric performance are the same for both generators and coolers. Thermoelectric modules with $Bi_2Te_3$ materials are usually employed in the cooling applications below room temperature but it can also be used for the power generation in the similar temperature range. In the present study, the power generation with a $Bi_2Te_3$ thermoelectric module has been investigated. The temperature difference between the hot and cold sides of the module is maintained with electric heater and cold water from the circulating water bath. The result shows that the electric current generated increases with temperature difference and decreases with the load resistance. However, the voltage increases with both the temperature difference and load resistance. The electric power increases with temperature difference and it has the maximum value when the load resistance is about 4 ${\Omega}$ for a given device.

Heat transfer augmentation based on the process intensification concept in heat exchangers and thermal reactors has received much attention in recent years, mainly due to energy efficiency and environmental considerations. The concept consists of the development of novel apparatuses and techniques that, compared to those commonly used today, are expected to bring dramatic improvements in manufacturing and processing, substantially decreasing equipment size, energy consumption, and ultimately resulting in cheaper, sustainable technologies. The objective of this paper was to investigate the heat transfer characteristics of tubular thermal reactor using static mixing technology. Glycerin and water were used as the test fluids and water was used as the heating source. The results for heat transfer rate were strongly influenced by tube geometry and flow conditions.

This paper presents the results of a series of performance tests for the integral Stirling cryocooler. Infrared sensor systems incorporating cryocoolers are required to be qualified to the appropriate environmental specification. Integral Stirling cryocooler for thermal imaging system have matured to the stage of undergoing formal qualification test program. The thermal environmental test of the Stirling cryocooler is presented in this paper. We performed that low and high temperature keeping test from $-40^{\circ}C$ to $+71^{\circ}C$ and operating test at high and low temperature cyclic range with acceptance tests performed at scheduled intervals. Cooldown time to 80K and steady state input power at 80K were determined as a function of cooler components temperatures at the compressor, hot end and cold tip. Tests performed on this cooler have been successful with a measured cooldown time to 80K of less than 5 minutes 24seconds for $71^{\circ}C$ ambient temperature with input power of 11W

SMES systems need cryogenic cooling systems. Conduction cooling system has more effective, compact structure than cryogen. In general, 2 stage GM cryocoolers are used for conduction cooling of HTS SMES system. 1st stages of cryocoolers are used for the cooling of current leads and radiation shields, and 2nd stages of cryocoolers for HTS coil. For the effective conduction cooling of the HTS SMES system, the temperature difference between the cryocooler and HTS coil should be minimized. In this paper, a cryogenic conduction cooling system for HTS SMES is analyzed to evaluate the performance of the cooling system. The analysis is carried out for the steady state with the heat generation of the HTS coil and effects of the thermal contact resistance. The results show the effects of the heat generation and thermal contact resistance on the temperature distribution.

The river water heat source heat pump has the advantage in the performance compared to air source heat pump. Although its better performance, the large temperature difference between load and source makes system performance worse by nature. In this study, 2-stage compression is considered as the solution of this problem. Generally, heat pump is designed for maximum capacity rate, but it actually operates at part load condition in many cases. Therefore, an information on the part-load character is very important in view of the system overall performance. In this study, part-load performance tests of a R134a 2-stage compression heat pump were carried out over the river water and supply heating water temperature changes. The experimental results show that the system performance is influenced by the part load rates, river water temperature, load temperature, etc.

The purpose of this study is to investigate the performance of a sea water heat source cascade heat pump system. R717(Ammonia) is used for a low-stage working fluid while R134a is for a high-stage. In order to gain a high temperature supply water in winter season, the system is designed to perform a cascade cycle. In this study, two experiments were carried out. One is a system starting test from the low load temperature of $10^{\circ}C$. The other is a system performance investigation over the R717 compressor capacity changes. Experimental results show that when it starts from the low load temperature, the suction temperature of the low-stage compressor is higher than that of a high-stage. The system performance increases when a water source temperature or a low-stage compressor rotational frequency goes higher.

In this study, Ice slurry generator heat transfer characteristics are experimentally investigated for the ice slurry generating system with scraper which is pneumatically operated. The ice slurry generator has the same shape as the vertical double tube type heat exchanger. Refrigerant is flowing in the outside tube and ethylene glycol solution in the inside tube. Refrigerant and solution water are parallel flow type which entering bottom of generator and leaving top of generator. The experimentations are conducted under the various test conditions such as compressor speed and cooling water temperature. For the above experimental conditions, heat transfer characteristics of the ice slurry generating system are evaluated in terms of the overall heat transfer coefficient and the heat transfer rate. And the experimental results show that the heat transfer rate of the system is increased as the compressor speed increases and the cooling water temperature decreases.

The purpose of this study is to investigate the performance of an autocascade refrigeration system using the refrigerant mixtures of R-134a(64${\sim}$72wt%), R-23(20${\sim}$25wt%) and R-14(8${\sim}$12wt%) as working fluids by experiment. The experimental apparatus consisted of a autocascade system with a compressor, condenser, evaporator, gas-liquid separator and cascade condenser, etc. Two main causes to this failure were considered as follows. First, the pressure ratio in the compressor was extremely increased at the beginning of the experiment. Second, the outlet temperature in compressor was higher than its limited temperature. For a stable operation of the compressor, a portion of the mixed refrigerant was by-passed at the inlet of the heat exchanger and transferred directly to the suction of the compressor in the modified system. The experimental parameters were changed weight of R-134a(1800g, 2800g), R-23(700, 800, 900g) and R-14(300, 400, 500g). The results showed that the best performance in case of R-134a of 2800g, R-23 of 900g and R-14 of 500g.

KEPRI has studied anode-supported planar SOFCs and kW class stacks operated at intermediate temperature for development of a combined heat and power unit. A single cell composed of Ni-YSZ/FL/ScSZ/LSCF showed the maximum power density of 0.55 W/$cm^2$ at $650^{\circ}C$ and 1.8 W/$cm^2$ at $750^{\circ}C$. With 37 cells of 10${\times}10cm^2$ and stainless steel interconnects, a 1kW class SOFC stack was manufactured. When a 1kW class SOFC system was operated at $750^{\circ}C$ with city gas, it showed the power output of 1.3 kWe at 50 A. It also recuperated heat of 0.57-1.2 kWth according to the loaded current through combustion of unreacted anode off-gas. Recently, KEPRI is developing a new kW class SOFC stack and system to increase efficiency and durability at intermediate temperature.

This study focuses on the hydrogen embrittlement of iron tube for fuel line of PEMFC (Proton Exchange Membrane Fuel Cell). PEMFC is operated by feed of hydrogen as a reactant and steam for proton conductivity of membrane. However, the environment with hydrogen and steam occur the hydrogen-induced degradation in BOP system. When iron tube was exposed to hydrogen and steam condition for 24 hours, the oxide layer on the surface was decreased by reduction. When the ambient temperature was 90$^{\circ}C$ micro cracks were found on the surface than any other temperature. The mechanical strength of iron tube was 3% lower than that of non-experiment tube. Maximum tensile stress was decreased 8%.

The reformer is one of the most important chemical processes for the production of high purity hydrogen from fossil fuel. This study compares zero-dimensional model with CFD models for reaction analysis of methane-steam reformer. The zero-dimensional model is an empirical equation, however CFD model uses reactions of Arrhenius type. Because the reaction coefficients of the steam-methane catalytic reforming have not been reported before in the form of Arrhenius type, the present study aims to find the appropriate reaction coefficients. The used CFD code is Fluent 6.2 version. Several models are compared for the case of various operating temperature, mass of catalyst and steam to methane ratio.

The transient power characteristics of a PEM fuel cell stack was experimentally studied using a commercial 1.2kW PEM fuel cell ($Nexa^{TM}$ Power Module, Ballard Power System Inc.). The conditions in PEM fuel cell stack such as temperature and water content change rather slowly because of their large heat capacity and long channel length, which results in long transient time to converge to a steady state. The steady characteristics of the PEM fuel cell module was determined first, followed by the measurement of its transient characteristics upon stepwise and continuous load current changes. During the stepwise current change from 5A to 25A, the output voltage initially decreased below the steady voltage and then increased gradually. Similar behavior was also observed for the stepwise current change from 25A to 5A. This transient behavior is explained with reference to the evolution of the temperature and water content of the PEM fuel cell stack.

Hydrogen is a fuel of fuel cell system, which has powerful explosion possibility. Hence, the fuel cell system needs safety evaluation to prevent risk of hydrogen leakage. We use a actual size chamber of a common fuel cell module to analyze hydrogen. Hydrogen injection holes are located in lower part of the chamber in order to simulated hydrogen leakage. The hydrogen sensor can detect range of 0${\sim}$4%. Since the hydrogen gas, of which leaked amount is controled by MFC, are injected at the bottom holes, the transient sensor signals are measured. At a condition of 10cc/s of hydrogen leakage, the sensor detects hydrogen leakage after 22sec and there is also several seconds of time delay depending on the position of the sensor. This experimental data can be applied for the design of the hydrogen detection system and ventilation system of a residential fuel cell system.

Diesel is an excellent candidate fuel for fuel cell applications due to its high hydrogen density and well-established infrastructure. But, it is hard to guarantee desirable performance of diesel reformer because diesel reforming has several problems such as sulfur poisoning of catalyst and carbon deposition. We have been focusing on diesel autothermal reforming(ATR) for substantial period. It is reported that ATR of diesel has several technical advantages such as relatively high efficiency and fuel conversion compared to steam reforming(SR) and partial oxidation(POX). In this paper, we investigate characteristics of diesel reforming under various ratios of reactants(oxygen to carbon ratio, steam to carbon ratio) for improvement of reforming performances(high reforming efficiency, high fuel conversion, low carbon deposition). We also exhibit calculated heat balance of autothermal reformer at each condition to help thermal management of SOFC system.

Performance of single cell at solid oxide fuel cell (SOFC) system is largely affected by electrocatalytic and thermal properties of cathode. Samarium-based perovskite oxide material is recently recognized as promising cathode material for intermediate temperature-operating SOFC due to its high electrocatalytic property. Perovskite structured $Sm_{0.5}Sr_{0.5}CoO_{3-\delta}$ and its composite material, $Sm_{0.5}Sr_{0.5}CoO_{3-\delta}/Sm_{0.2}Ce_{0.8}O_{1.9}$ were investigated in terms of area specific resistance (ASR), thermal expansion coefficient (TEC), thermal cycling and long term performance. $Sm_{0.2}Ce_{0.8}O_{1.9}$ was used as electrolyte material. Electrochemical ac impedance spectroscopy (EIS) and dilatometer were used to measure the cathodic properties. Composite cathode ($Sm_{0.5}Sr_{0.5}CoO_{3-\delta}$: $Sm_{0.2}Ce_{0.8}O_{1.9}$ = 6:4) showed a good ASR of 0.13${\Omega}$ $cm^2$ at 650$^{\circ}C$ and its TEC value was 12.3${\times}$10-6/K at 600$^{\circ}C$ which is similar to the value of ceria-based electrolyte of 11.9${\times}$10-6/K. Performance of composite cathode was maintained with no degradation even after 13 times thermal cycle test.

Hydrogen is receiving much research attention as an alternative substitute for hydrocarbon fuels these days due to its cleanliness and renewability. However, hydrogen should be used with caution because of its high propensity for leak and wide flammable range. This study deals with a situation that hydrogen leaks and then forms a flammable mixture inside 1kW class residential fuel cell. The residential fuel cell was modeled as a box-shaped chamber with vent openings at the top and bottom, filled with various components such as reformer, desulfurizer, fuel cell stack and humidifier. Computational fluid dynamics (CFD) was used to simulate the diffusion, buoyant flow and accumulation of leaked hydrogen in the modeled chamber. From the simulation, the risk region vulnerable to flame was identified and the methods to minimize such hazardous region was discussed. When the vent openings are 1% of the total surface, as the quantity of hydrogen leakage increases the risk regions increases accordingly. As the vent openings of the total surface increased from 1% to 2.3%, averaged hydrogen mole fraction is under 1% in the system.

Numerical analysis was conducted to investigate the atomic layer deposition(ALD) of silicon nitride using silane and ammonia as precursors. The present study simulated the surface reactions for as-deposited $Si_3N_4$ as well as the kinetics for the reactions of $SiH_4$ and $NH_3$on the semiconductor wafer. The present numerical results showed that the ALD process is dependent on the activation constant. It was also shown that the low activation constant leads to the self-limiting reaction required for the ALD process. The inlet and wafer temperatures were 473 K and 823 K, respectively. The system pressure is 2 Torr.

Analysis of the internal state of the blast furnace is needed to predict and control the operating condition. Especially, it is important to develop modeling of blast furnace for predicting cohesive zone because shape of cohesive zone influences on overall operating condition of blast furnace such as gas flow, temperature distribution and chemical reactions. Because many previous blast furnace models assumed cohesive zone to be fixed, they can't evaluate change of cohesive zone shape by operation condition such as PCR, blast condition and production rate. In this study, an axi-symmetric 2-dimensional steady state model is proposed to simulate blast furnace process using the general purpose-simulation code. And Porous media is assumed for the gas flow and the potential flow for the solid flow. Velocity, pressure and temperature distribution for gas and solid are displayed as the simulation results. The cohesive zones are figured in 3 different operating conditions.

In order to predict the thermal radiation induced from alumina particle cloud in the plume of solid propellant motor, view factor method is applied to space shuttle SRB and the result is compared with that of monte carlo method. For this purpose, radiative characteristics, such as particle cloud temperature distribution, effective emissivity or emissive power of particle cloud are studied. In the case of effective emissivity, inverse wavelength method is applied and plume reduction characteristic length is used for emissive power distribution. As a result, thermal radiation using view factor method gives more conservative results than that using monte carlo method. So it can be used for preliminary design of thermal protection system.

This study numerically and experimentally investigated on thermal strain analysis of aluminum alloy casting mold using metal foundry. To predict the numerical result of thermal strain in Al alloy casting mold during the cooling process, it is performed the investigation of temperature distribution, stress and displacement based on the physical properties of Al alloy. In results of this study, Al alloy casting mold represented rapidly cooling graph during initial 20minutes after beginning cooling process, therefore value of stress and displacement is rapidly changed during initial 20minutes after beginning cooling process. In addition to, temperature distribution obtained by experiment confirmed corresponding pattern then compared numerical analysis with experiment. These results are distribute to make the effective and the high precision casting mold.

Microfin arrays with fin heights of 100 ${\mu}$m and 200 ${\mu}$m and six different spacings from 30 ${\mu}m$to 360 ${\mu}m$ are fabricated using the DRIE process. Natural convective heat transfer around the microfin arrays on both vertical and horizontal surfaces is experimentally examined. It turns out that the orientation effect of microfin arrays is negligible compared with macrofin arrays. The obtained heat transfer coefficients are compared with the existing heat transfer correlation for the macrofin arrays. It is concluded that the existing macrocorrelation is no longer valid for the microfin arrays. Relevant empirical correlations for microfin arrays on the vertical and horizontal surfaces are presented based on the present experimental data.

Recently, ultra-lightweight materials with open, periodic cell structures take much attention owing to its potential for multi-functionality such as load bearing, thermal dissipation, and actuation. This paper presents experimental results on the hydraulic and heat transfer characteristics for the Wire-woven Bulk Kagome(WBK) composed of aluminum 1100 wires. The overall pressure drop and heat transfer of the WBK specimen have been experimentally investigated under forced air convection condition. The pressure loss and heat transfer performance of the aluminum WBK are compared with other heat dissipation media. It was shown that heat transfer depended on relative density and surface area density. Comparison with metal foams and other heat dissipation media such as packed beds, lattice frame materials, louvered fins, and other materials suggests that the aluminum WBK competes favorably with the best available heat dissipation media in heat transfer performance.

Understanding of the impinging behavior of an electrically charged spray is essential in determining appropriate operating conditions for electro-spraying of paints, surface coating materials and insecticides. In the present work, as an initial step, the wall impact of an electrically charged droplet has been experimentally investigated. The charged drops were directed on the surface of a paraffin wax, and the impinging behavior was visualized and recorded using a CCD camera to identify the impingement regime. The spread-rebound boundary for the charged drop turned out to be smaller compared to that for an electrically neutral droplet under the same surface condition. The shift of the transition criterion is considered to be due to the discrepancy between the maximum spread ratio of the electrically charged droplet and that of the neutral droplet.

This study was performed to develop a dryer for home clothes using thermoelectric module. The thermoelectric module was used as a heat source and a dehumidification device because it has heating part and cooling part at once. To design for maximizing the energy efficiency and the rate of dehumidification, the parameters of the dryer using thermoelectric module are heat capacity and air flow rate. This study showed that the thermoelectric module can be used in the clothes dryer and energy efficiency of clothes dryer be better than that of electric heating dryer.

A small-scale thermoelectric cooling system was built in an effort to enhance the performance of the refrigeration system by utilizing the water-cooled jacket which was attached on the hot side of the thermoelectric module. Considered design parameters for the water-cooled jacket included the geometry of the flow passage inside the jacket and the flow rate of cooling water. The higher flow rate of cooling water in the jacket resulted in a better performance of the refrigeration system. The introduction of geometrical complexity of the cooling water flow passage to the cooling jacket also showed significant improvement on the performance of the thermoelectric refrigeration system such as the cooling capacity and the COP of the refrigeration system.

This paper describes the possible way to improve the capacity, the efficiency and the pressure drop of $CO_2$ system. It is considered the use of an internal heat exchanger (IHX) to improve the performance of the $CO_2$ system. Experiment was performed by changing a tube shape, a tube number and a tube length of IHX to investigate the performance of IHX for $CO_2$ system. The focus of the present study is to obtain the concept on IHX optimal design. Experimental results show that design parameters are closely related with the capacity and the pressure drop of $CO_2$ system. In the transcritical $CO_2$ cycle, the system performance is very sensitive to the IHX design. System performance on operation condition and shape of IHX is also introduced.

This paper presents a mathematical model for predicting the frost behavior formed on heat exchanger fins, considering fin heat conduction under frosting condition. The model is composed of air-side, the frost layer, and fin region, and they are coupled to the frost layer. The frost behavior is more accurately predicted with fin heat conduction considered (Case A) than with a constant fin surface temperature assumed (Case B). The results indicate that the frost thickness and heat transfer rate for Case B are over-predicted in most regions of the fin, as compared to those for Case A. Also, for Case A, the maximum frost thickness varies little with the fin length variations, and the extension of the fin length over 30 mm contributes insignificantly to heat transfer.

This paper presents an optimal design of a micro evaporator which maximizes the heat transfer coefficient. Number of gaps, spanwise distance and streamwise distance are selected as the geometric design parameters. Mass flow rate of the refrigerant is selected as the non-geometric design parameter. Temperature at the surface of the heater is measured to valuate the heat transfer coefficient. Nine experiments are conducted using $L_9(3^4)$ orthogonal array. Maximum heat transfer coefficient is 640 W/$m^2K$ at the parameters of 2 gaps, 0.2 mm spanwise distance, 1.0 mm streamwise distance and 0.72 g/s mass flow rate. Among the 3 geometric parameters, the spanwise distance is the most sensitive parameter influencing the heat transfer coefficient. We conduct a second stage of experiment to increase the heat transfer coefficient by reselecting the mass flow rate. We concluded that 0.87 g/s is the optimized flow rate for an active micro cooler resulting in a heat transfer coefficient of 651 W/$m^2K$.

GSHP systems are used for air-conditioning systems in commercial buildings, schools, and factories because of low operating and maintenance costs. These systems use the earth as a heat source in heating and a heat sink in cooling mode. Ground heat exchangers are classified by a horizontal and vertical type according to the installation method. Vertical type is usually constructed by placing small diameter high density polyethylene tube in a vertical borehole. Vertical tube sizes range from 20 to 40 mm nominal diameter. Borehole depth range between 100 and 200 m depending on local drilling conditions and available equipment. In this study, to evaluate the performance of single u-tube with bentonite grouting, single u-tube with broken stone grouting and double u-tube bentonite grouting of vertical ground heat exchangers, test sections are buried on the earth and experimental apparatus is installed. Therefore the heat transfer performance and pressure loss of these are estimated.

The air and water flow distribution are experimentally studied for a round header-ten microchannel tube configuration. Three different inlet orientations (parallel, side, normal) were investigated. Tests were conducted with downward flow configuration for the mass flux from 70 to 130 kg/$m^2s$, quality from 0.2 to 0.6, non-dimensional protrusion depth (h/D) from 0.0 to 0.5. It is shown that, for almost all the test conditions, normal inlet yielded the best flow distribution, followed by side and parallel inlet. Possible reasoning is provided using flow visualization results.

In this study, a mathematical model for a thermal analysis of a flat heat pipe with a grooved wick structure is presented. The effects of the liquid-vapor interfacial shear stress, the contact angle, and the amount of liquid charge have been included in the proposed model. In particular, the axial variations of the wall temperature and the evaporation/condensation rates are considered by solving the one-dimensional conduction and the augmented Young-Laplace equations, respectively. In order to verify the model, the results obtained from the model are compared to existing experimental data.

This study deals with performance characteristics of heat pipe combined with a solid-liquid phase change material(PCM). The outer diameter of the heat pipe was 9.5 mm and the total length was 600 mm, where the evaporator, the adiabatic section and the condenser lengths were equally 200 mm. A paraffin wax having a melting point of 58.5$^{\circ}C$ was used as PCM. The paraffin container was attached to the adiabatic section of the heat pipe. The paraffin container had outer diameter of 18 mm, wall thickness of 1.2 mm and the total length of 100 mm. The heat pipe was tested with tilt angle of horizontal degree and favorite angle 10 degree, with evaporator lower position to provide stable operation of the heat pipe. Input thermal load was varied from 40 W, with increment of 40 W, to above 100 W until the maximum temperature of the heat pipe wall reached 200$^{\circ}C$. Test results of the PCM heat pipe were presented in comparison with conventional heat pipe of the same basic dimensions. The performance was analyzed in terms of temperature distribution, thermal resistance and heat transport capability.

A VCHP was fabricated and tested for its thermal performance. The container was made of copper, and the working fluid was water. STS-316 screen of mesh number 100 was inserted as a capillary structure. As a baseline performance, a normal heat pipe of the same dimensions was tested in advance to compare with VCHP, where an inert gas container was attached. The outer diameter of the heat pipe was 12.8 mm and the total length was 600 mm. The evaporator and the condenser lengths were both 200 mm. The thermal load ranged from 20 to 300W. Typical result revealed that the operating temperature of the VCHP stayed almost constant, while that of the normal heat pipe varied as much as 40$^{\circ}C$. Therefore, it was demonstrated that the VCHP is very effective for temperature control of heat-dissipating devices.

Operating characteristics of a loop heat pipe (LHP) having a bypass line was investigated experimentally. The LHP had a sintered metal wick as a capillary structure and methanol as a working fluid. The sintered metal wick was made of stainless steel of which the average pore size was 5 ${\mu}m$and porosity of 47%. A bypass line of a small diameter was attached between the vapor escape passage and the liquid reservoir. The dimension of the flat evaporator was $50(L){\times}40(W){\times}30(H)$ mm and that of the condenser was $50(L){\times}40(W){\times}11(H)$ mm. Wall and pipe material of the LHP was stainless steel and heating area was 35(W) mm${\times}$35(L) mm. The inner diameters of vapor and liquid transport lines were 4.0 mm and 2.0 mm, and the lengths of the two lines were both 0.5 m. The LHP was tested for three different tilt angles of horizontal, favorite tilt, and adverse tilt. The thermal load range was up to 290 W at the condenser above evaporation position. Furthermore, the effect of a bypass line on the start-up transient as well as steady-state operation was presented and discussed.

The Loop Heat Pipe(LHP) system uses capillary forces so as to pump the working fluid from heat acquisition to heat rejecting systems. The performance of the LHP systems depends mainly upon the operating performance of the wick structure. The capillary pressure increases with decreasing the pore size of the wick structure. By the way, the wick structure's permeability decreases with decreasing the pore size and the porosity. To obtain an ideal wick, the wick structure should possess several characteristics such as the small pore size, high porosity and chemical compatibility with working fluid. Sintered metal wicks have been mainly used as the capillary wick structure mounted in LHP because of the fact that the sintered metal wick has some advantages like convenient selection of wick material, smaller pore size and so on as well as high reliability. In this study, sintered metal wicks were developed to meet required several parameters to design the high performance LHP systems for obtaining even more effective cooling technologies.

Experimental study was conducted to evaluate the performance of a miniature loop heat pipe (MLHP) with non-inverted meniscus type capillary structure. All parts of MLHP in this study were made of copper including the capillary structure and the distilled water was used as a working fluid of MLHP. The outer diameter of evaporator was 9 mm and its length was 119 mm. The effective pore size of the capillary structure was 30 micron and its porosity was 60%. The vapor transport line, the liquid transport line and the condenser were consisted of single 4.0 mm copper tube. The distance between the evaporator and the condenser region was 200 mm and the length of the loop was 969 mm. This MLHP was operated successfully at any orientation but the gravity highly influenced the thermal performance of the MLHP. The maximum thermal load was 130 watts at the bottom heat mode and the 20 watts at the top heat mode.

The electricity conversion-efficiency of solar cell for commercial application is about 6-15%. More than 85% of the incoming solar energy is either reflected or absorbed as heat energy. Consequently, the working temperature of the photovoltaic cells increases considerably after prolonged operations and the cell's efficiency drops significantly. PV/T refers to the integration of a PV module and a solar thermal collector in a single piece of equipment. By cooling the PV module with a fluid steam like air or water, the electricity yield can be improved. At the same time, the heat pick-up by the fluid can be to support space heating or service hot-water systems. In this study, a pulsating heat pipe solar heat collector was combined with single-crystal silicon photovoltaic cell in hybrid energy-generating unit that simultaneously produced low temperature heat and heat and electricity. This experiment was investigating thermal and electrical efficiency for evaluation of a PV/T system.

The potential of using a double-faced wall stagnation flow burner in mass production of carbon nanotubes was evaluated experimentally and computationally. With nitrogen-diluted premixed ethylene-air flames established on the Nickel-coated stainless steel double-faced wall, the propensities of carbon nanotube formation were experimentally determined using SEM and FE-TEM images and Raman spectroscopy, while the flame structure was computationally predicted using a 3-dimensional CFD code with a reduced reaction mechanism. The uniformity and yields of synthesized carbon nanotubes were evaluated in terms of the flame stretch rates. Results show substantial increase of area on the wall surface where uniform carbon nanotubes are synthesized with using the double-faced wall stagnation flow burner due to enhanced uniformity of temperature distribution along the wall surface and support the potential of using a double-faced wall stagnation flow burner in mass production of carbon nanotubes.

In this study, the mechanism of enhanced thermal conductivity is elucidated on the bases of both electric double layer (EDL) and kinetic theory. A novel expression for the thermal conductivity of nanofluids is proposed and verified by applying to $Al_2O_3$ nanofluids with regard to various temperatures, volume fractions and particle sizes. In dilute nanofluids, the effects of Brownian motion and particle interaction on enhancing the thermal conductivity of nanofluids are quite comparable while the effect of particle interaction due to EDL is more prominent in dense nanofluids. The model presented in this paper shows that particle interaction due to the electrical double layer is the most responsible for the enhancement of thermal conductivity of nanofluids.

At nanoscales, the Boltzmann transport equation (BTE) can best describe the behavior of phonons which are energy carriers in crystalline materials. Through this study, the phonon transport in some micro/nanoscale problems was simulated with the Monte Carlo method which is a kind of the stochastic approach to the BTE. In the Monte Carlo method, the superparticles of which the number is the weighted value to the actual number of phonons are allowed to drift and be scattered by other ones based on the scattering probability. Accounting for the phonon dispersion relation and polarizations, we have confirmed the one-dimensional transient phonon transport in ballistic and diffusion limits, respectively. The thermal conductivity for GaAs was also calculated from the kinetic theory by using the proposed model. Besides, we simulated the electrostatic discharge event in the NMOS transistor as a two-dimensional problem by applying the Monte Carlo method.

The 3-omega (3-${\omega}$) method is utilized to measure the thermal conductivity of nanofluids. A metal line heater on a silicon nitride membrane bridge structure is microfabricated by a bulk silicon etching method. Localized measurement of the thermal conductivity within the nanofluids droplet is possible by the fabricated 3-${\omega}$ sensor. Time varying AC temperature amplitudes and thermal conductivities are measured to check the stability of the nanofluids containing multi-wall carbon nanotubes (MWCNTs). Stabilities of MWCNT nanofluids prepared with different chemical treatments are compared. Acid treated MWCNT showed best dispersion stability in water while MWCNTs dispersed in water with surfactants such as Gum Arabic and Sodium dodecyl benzene sulfate showed clear sign of gravity dependence.

In this paper, effective thermal conductivities of water-based $Al_2O_3$-nanofluids with low concentration from 0.01 vol. % to 0.3 vol. % are experimentally obtained by transient hot wire method (THWM). The water-based $Al_2O_3$-nanofluids are manufactured by two-step method which is widely used. To examine suspension and dispersion characteristics of the water-based $Al_2O_3$-nanofluids, Zeta potential as well as transmission electron micrograph (TEM) is observed. We confirm the manufactured $Al_2O_3$-nanofluids have good suspension and dispersion. The effective thermal conductivities of the water-based $Al_2O_3$-nanofluids with low concentration are enhanced up to 1.64% compared with that of DI water at $21^{\circ}C$. In addition, experimental results are compared with theoretical results from Jang and Choi model.

The relative contributions of phonon and electron to the thermal conductivity of silicon film with varying doping density are evaluated from the modified electron-phonon interaction model, which is applicable to the micro/nanoscale simulation of energy transport between energy carriers. The thermal conductivities of intrinsic silicon layer thicknesses from 20 nm to 500 nm are calculated and extended to the variation in n-type doping densities from 1.0 ${\times}$ $10^{18}$ to 5.0 ${\times}$ $10^{20}$ $cm^{-3}$, which agree well with the experimental data and theoretical model. From simulation results, the phonon and electron contributions to thermal conductivity are extracted. The electron contribution in the silicon is found to be not negligible above $10^{19}$ $cm^{-3}$, which can be classified as semimetal or metal by the value of its electrical resistivity at room temperature. The thermal conductivity due to electron is about 57.2% of the total thermal conductivity at doping concentration 5.0 ${\times}$ $10^{20}$ $cm^{-3}$ and silicon film thickness 100 nm.

Since the use of CFCs and HCFCs refrigerants are to be restricted due to the depletion of ozone layer, this experiment applies the R134a gas to study the performance characteristic from the superheat control for improving the energy efficiency. The experiments are carried out for the condensing pressure of refrigeration system from 1.5 MPa to 1.6 MPa by 0.05 MPa and for superheat temperature from 0$^{\circ}C$ to 5$^{\circ}C$ by 1$^{\circ}C$ at each condensing pressure. As a result of experiment, when the superheat temperature is 1$^{\circ}C$ at each condensing pressure, the refrigeration system has the highest performance.

An experimental study is performed to investigate the thermal and flow characteristics of subcooled liquid nitrogen in a natural circulation loop. Experimental apparatus is designed and constructed such that a closed loop is cooled at the top by a cryocooler and heated nearly at the bottom by cartridge heaters. Steady state is obtained by controlling the heating power to the cartridge heaters and a thin-film heater to reduce the cooling power of the cryocooler. Temperature is measured at several locations of the loop and the mass flow rate through the loop is estimated from the energy balance in terms of the measured temperatures. Experiment is repeated for various values of the vertical height between the cooling and heating parts. The results show that the heat transfer capability of the loop has a maximum at a certain value of height. The optimal height to maximize the heat transfer is in a good agreement with analytical prediction to take into account the buoyancy and viscous forces in the loop.

In order to investigate the variation on a heat transfer coefficient during evaporation of $CO_2$, basic experiment on the evaporation heat transfer characteristics in a horizontal micro-channel tube was performed. Hydraulic diameters of micro-channels were 0.68 and 1.46 mm. The experiment apparatus consisted of a test section, a DC power supply, a heater, a chiller, a mass flow meter, a pump and a measurement system. Experiments were conducted for various mass fluxes of 300 to 800 kg/$m^2s$, heat fluxes of 10 to 40 kW/$m^2$ and saturation temperatures of -5 to 5$^{\circ}C$. With the increase heat flux, the evaporation heat transfer coefficient increased. And the significantly change of the heat transfer coefficient was observed at any heat flux and mass flux. As the saturation temperature increased and the hydraulic diameter decreased, the heat transfer coefficient increased.

As a power transmission line supplying power to a densely populated city, the high temperature superconducting (HTS) cable is expected to one of the most effective cables with a compact size because of its high current density. The verification of HTS power cable system have been progressed by KEPRI. A cooling system for a 3-phase 100m HTS power cable with 22.9kV/1.25kA was installed and tested at KEPCO's Gochang power testing center in Korea. The system consists of a liquid nitrogen decompression cooling system with a cooling capacity of 3kW and a closed circulation system of subcooled liquid nitrogen. Several performance tests of the cable system with respect to the cooling such as cooling capacity, heat load and temperature stability, were performed at several temperatures. Thermal cycle test, cool-down to liquid nitrogen temperature and warm-up to room temperature, was also performed to investigate thermal cycle influences. The outline of the installed cooling system and performance test results are presented in this paper.

Liquefied gas vaporizer means machine to vaporize the liquefied gas as liquid nitrogen($LN_2$), liquefied natural gas(LNG), liquid oxygen($LO_2$) etc. In the air type vaporizer, the frozen dew is also created by temperature drop (below 273 K) on vaporizer surface. This problem increases as the time progresses and humidity increases. In addition, the frozen dew gradually becomes frost deposit consequently, heat transfer through vaporizer decreases because frost deposit form adiabatic sheet. Because of this reason, recent vaporizer system is installed as parallel type, this vaporizer system needs more expensive installation costs and more space. This paper was investigated on the heat transfer characteristics of liquefied gas vaporizer with super low temperature and this paper was carried out the numerical about air heating vaporizer with super low temperature. The numerical analysis on the heat transfer was studied on the effect of geometric parameters of the vaporizer, which are length 1000 mm of 4fin75le type vaporizer. 4fin75le means number of fin is 4 and height of fin is 75 mm.

Because of the ozone layer depletion and global warming, new alternative refrigerants are being developed. In this study, evaporation heat transfer characteristics of carbon dioxide flowing upward in a vertical tube have been investigated by experiment. Before the test section, a pre-heater is installed to adjust the inlet quality of the refrigerant to a desired value. A smooth tube with outer diameter of 5 mm and length of 1.44 m was selected as a test tube. The test was conducted at mass fluxes of 212 to 530 kg/$m^2s$, saturation temperature of -5 to 20$^{\circ}C$, and heat fluxes of 20 to 45 kW/$m^2$. As the vapor quality and mass fluxes increase, the heat transfer coefficients of carbon dioxide are decreased, and the heat transfer coefficients increase when the heat fluxes and saturation temperatures increase.

Newly structured metal-supported solid oxide fuel cell was fabricated and characterized by impedance analysis and galvanodynamic experiment. Using a cermet adhesive, thin ceramic layer composed of anode(Ni/YSZ) and electrolyte(YSZ) was joined with STS430 metal support of which flow channel was fabricated. $La_{0.8}Sr_{0.2}Co_{0.4}Mn_{0.6}O_3$ perovskite oxide was used as cathode material. Single cell performance was increased and saturated at operating time to 300hours at 800$^{\circ}C$ because of cathode sintering effect. The sintering effect was reinvestigated by half cell test and exchange current density was measured as 0.005A/$cm^2$. Maximum power density of the cell was 0.09W/$cm^2$ at 800$^{\circ}C$. Numerical analysis was carried out to classify main factors influencing the single cell performances. Compared to experimental IV curve, simulated curve based on experimental parameters such as exchange current density was in good agreement.

Water gas shift reactor in fuel processing is an important part that converts carbon monoxide into hydrogen. Fuel processing system for PEMFC usually has two stages of WGS reactors, which are high temperature and low temperature shifter. In this study we prepared noble metal catalysts and compared their performances with that of a commercial iron chromium oxide catalyst. Noble metal catalysts and the commercial catalyst showed quite different temperature dependence of carbon monoxide conversion. The conversion of carbon monoxide at the commercial catalyst was very low at medium temperature(${\sim}300^{\circ}C$) and increased rapidly as temperature increased while the conversion at noble metal catalysts was high in the medium temperature range and decreased as temperature increased, which is thermodynamically expected. Their characteristics agreed well with the literature published, and we are accomplishing further study for improvement of the noble metal catalysts.

Steam reforming and catalytic reforming of $CH_4$ conversion to produce synthesis gas require both high temperatures and high pressure. Non-thermal plasma is considered to be a promising technology for the hydrogen rich gas production from methane. In this study, three phase AC GlidArc plasma system was employed to investigate the effects of gas composition, gas flow rate, catalyst reactor temperature and applied electric power on the $CH_4$ and $H_2$ yield and the product distribution. The studied system consisted of three electrode and it connected AC generate power system different voltages. In this study, air was used for the partial oxidation of methane. The results showed that increasing gas flow rate, catalyst reactor temperature, or electric power enhanced $CH_4$ conversion and $H_2$ concentration. The reference conditions were found at a $O_2$/C molar ratio of 0.45, a feed flow rate of 4.9 ${\ell}$/min, and input power of 1kW for the maximum conversions of $CH_4$ with a high selectivity of $H_2$ and a low reactor energy density.

A new microcombustor configuration for a micro fuel-cell reformer integrated with a micro evaporator was studied experimentally and computationally. The present microcombustor is simply cylindrical to be easily fabricated but two-staged, expending downstream, to feasibly control ignition and stable burning. Results show that the aspect ratio of the first stage and the wall thickness of the microcombustors substantially affect ignition and thermal characteristics. For the optimized design conditions, a premixed microflame was easily ignited in the expanded second stage combustor, moved into the smaller first stage combustor, and finally stabilized therein. The measured and predicted temperature distributions across the microcombustor walls indicated that heat generated in the microcombustor is well transferred. Thus, the present microcombustor configuration could be applied to the practical micro reformers integrated with a micro evaporator for use of fuel cells.

A directly heated $SO_3$ decomposer for the sulfur-iodine and hybrid-sulfur processes has been introduced and analyzed by using a computational fluid dynamics code(CFD) with the CFX 5.7.1. The use of a directly heated decomposition reactor in conjunction with a VHTR allows higher decomposition reactor operating temperature. However, the thermochemical and hybrid hydrogen production processes accompanied with the high temperature and strongly corrosive operating conditions basically have material problems. In order to resolve these problems, we carried out the development of a structural material and equipment design technologies. The results show that the maximum temperature of the structural material (RA330) could be maintained at 800$^{\circ}C$ or less. Also, it can be seen that the mean temperature of the reaction region packed with catalysts in the $SO_3$ decomposition reactor could satisfy the temperature condition of around 850 $^{\circ}C$ which is the target temperature in this study.

Comprehensive analytical models focusing on the anode water loss, the cathode flooding, water equilibrium, and water management strategy are developed for polymer electrolyte fuel cells. Analytical solutions presented in this study are compared with two-dimensional computational results and shows a good agreement in predicting those critical characteristics of water. General features of water concentration profile as a function of membrane thickness and current density are presented to illustrate the net effect of the back-diffusion of water from the cathode to anode and the water production by the cathode catalytic reaction on water transport over a fuel cell domain. As one of practical applications, the required humidity level of feed streams for full saturation at the channel outlets are investigated as a function of the physical operating condition. These analytical models can provide good understanding on the characteristic water

This Paper focuses on modeling and simulation to analyze the characteristic of hybrid vehicle. The system includes proton exchange membrane fuel cell(PEMFC), photovoltaic generator(PV), lead-acid battery, motor, vehicle and controller. Main electricity is produced by the PEMFC and battery to meet the requirements of a user load. When vehicle is parked in a sunny place, extra power is generated by the photovotaics and is charged in a battery for next drive. Further we evaluate usefulness of this hybrid vehicle by using ADVISOR - the advanced vehicle simulator written in the Matlab/Simulink environment. According to simulation results, the extra power obtained by photovoltaics which have been explored in nature conditions can help to reduce the electrical load of PEMFC and increase the efficiency (over 30%).

In this study, the numerical analysis to estimate condensation heat and mass transfer of the condenser was carried out using the PMA (porous medium approach). In the PMA, the details of tube bundle in the condenser are replaced by the porous medium, and the flow resistance term is added in the momentum equation. In this regard, the PMA is quite helpful for the study of tube bundle in the large condenser. The pressure loss through tube bundle can be compensated by viscous and inertial momentum sink terms, which was validated numerically. Value of the pressure drop was compared to that of Butterworth correlation. Three dimensional analysis of condensation for McAllister condenser with the PMA was conducted using Fluent 6.2 and UDFs (use-defined functions). The result of condensation rate was analogous to previous results (experimental and numerical data).

The present study investigates the two dimensional flow and heat/mass transfer characteristics of wavy duct with various corrugation angles. For the heat/mass transfer coefficients, a naphthalene sublimation technique is used. Numerical analysis and wall pressure measurement show detailed two dimensional flow features. The corrugation angles change from 145$^{\circ}$ to 100$^{\circ}$. The operating Reynolds numbers based on the duct hydraulic diameter vary from 700 to 3,000. The duct aspect ratio maintains 7.3. On the pressure wall, strong flow mixing enhances heat/mass transfer coefficients at the front position. In addition, the rear side of pressure wall, the near of peak, is affected by the acceleration and the shedding of main flow. On the suction wall, however, flow separation and reattachment lead to the valley and the peak of heat/mass transfer coefficient. Also, highly increasing boundary layer at the suction wall affects the decrease of heat/masst transfer. As decreasing corrugation angles, the spanwise average Sherwood number increases and the peak or the valley positions of the local Sherwood number are varied.

Brazed Plate Heat Exchanger (BPHE) is a type of compact plate heat exchanger with parallel corrugated plates which are brazed together in series. Each plate hascorrugation called herringbone pattern. Inside a BPHE, hot fluid and cold fluid alternate its flow direction to establish counter current flow configuration. Two-phase flow heat transfer and pressure drop of R-22 in BPHE were experimentally measured in this study. In the present experiments, single-phase region and two-phase region coexist in a BPHE. Therefore, the inside of a BPHE have to be divided into single phase region and two phase region and analyzed accordingly. The results from the single phase flow analysis are then extended to the two phase flow analysis to correlate the condensation and evaporation heat transfer and pressure drop for the refrigerant R-22 in the BPHEs. Previous models for two- phase friction factor have been compared with the present experimental results.

This study is intended to analyze the effect of thin ink-film thickness around rotating printing roll on the printing quality in the gravure printing process which is used for making electronics circuit like a RFID tag with a conductive ink. The present work numerically estimates the film thickness around rotating roller partially immersed in ink, for which the volume of fluid (VOF) method was adopted to figure out the film formation process around rotating roller. Parameter studies were performed to compare the effect of ink viscosity, surface tension, roller rotating speed, immersed angle on the film thickness. The result indicates that the film thickness has a strong dependency on the fluid viscosity, while the surface tension has negligible effect.

The main objective of this study is to analyze the heat transfer characteristics in the electric melting furnace. Local temperatures are measured at various location in the furnace using the B-type thermocouples. In this paper, the numerical simulation was performed using the ANSYS software, and compared with experimental data. Mathematical heat transfer model for the prediction of temperature distribution has been developed by considering the thermal radiation among heating element, crucible and insulating materials. The results show that the temperature distributions predicted by the numerical simulation agree with experimental results comparatively.

Thermal transport from vertical heated surface to falling liquid film in a channel has been investigated experimentally. Air-flow is introduced into channel to make a counter flow against falling liquid film. This problem is of particular interest in the design of direct contact heat exchange system, such as cooling tower, evaporative cooling system, absorption cooling system, and distillation system. The effects of channel width and air flow rate on the heat transfer to falling liquid film are studied in detail. The results obtained indicate that heat transfer rate is gradually decreased with an increase in the channel width without air flow as well as with air flow in a channel. It is also found that heat transfer rate of air-flow is increased while heat transfer rate of falling liquid film is decreased with an increase in the air flow rate at a given channel width. However, total heat transfer rate form the heated surface is increased as the air flow rate is increased.

Numerical simulations are conducted to analyze conjugate heat transfer characteristics in a ribbed channel. In this simulation, the effects of Reynolds number and heat capacity of the solid channel wall on convective heat transfer are observed in the turbulent flow regime. In the case of the conducting wall against isothermal wall, the relative ratio of the thermal resistance between the solid wall and the flow field varies with Reynolds number. Thus the characteristics of the conjugate heat transfer are changed with the Reynolds number. Heat capacity ratio affects the temperature fluctuation inside solid wall. The temperature fluctuation inside the solid wall decreases with increasing the heat capacity of the solid wall so that the convective heat transfer increases. When the thermal conductivity ratio is smaller than 10, the effects of flow characteristics on heat transfer are changed.

Experiments were performed to determine the thermal (or turbulent) diffusion coefficient (TDC) and to investigate the critical heat flux (CHF) performance in the 5${\times}$5 rod bundle with 5 unheated rods which are supported by Hybrid Mixing Vane. In this study, HFC-134a fluid was used as working fluid and the fluid temperature were measured in the important subchannels. To determine the TDC value, the measured fluid temperatures were compared with the predicted values obtained from the MATRA code. The best optimized value of ${\beta}$ was found to be 0.02 by considering prediction statistics, i.e., average and standard deviations of the differences between the experimental results and code calculations. Using the best optimized value of ${\beta}$ as 0.02, the MATRA code predicts the test results of the fluid temperature within ${\pm}$1.0 % of error. According to the experimental results on CHF of 5 non-heating guide tubes, the case with non-heating guide tube showed a little good performance in terms of CHF.

The application analysis of bubble pump on the domestic solar water heater system is presented. The system investigated in this study is a passive device, self pumping and self regulating. It was test to use the bubble pump on solar water heater system. The test experiment has been taken on the existed vacuum tube about the efficiency, working fluid temperature and pressure and circulated power. In order to check the working temperature and working pressure effectively, the bubble pump was test separated from the solar water heater. The equipment consists of the bubble pump, heater and heat exchanger. The main structure of bubble pump was design depend on the character of two phase flow. The complete system was instrumented to measure pressures, temperatures and flow-rates at various locations. The theory analysis of design bubble pump has been given and the experiment design has been included in the paper.

This paper presents an investigative study on the efficacy of a new physical water treatment (PWT) technology using an oscillating electric field to mitigate mineral fouling in heat exchangers. Parallel graphite electrode plates immersed in water were used to generate the electric field directly in water. Artificial hard water at 500 ppm hardness was used in all fouling tests. The inlet temperatures were maintained at 23.5${\pm}$0.5$^{\circ}C$ and 85${\pm}$0.5$^{\circ}C$ for cold and hot water sides, respectively. The results at a cold water-side velocity of 0.3 m/s showed a 16-60% drop in fouling resistances from the baseline test depending on the frequency of the electric field for the PWT-treated cases.

To return LNG to a gaseous state, it is fed into a regasification plant. For inland area, where pipelines do not exist or difficult to construct, the LNG is delivered to the inland-receiving terminal available at that area, regasified and delivered to consumers. At inland-receiving terminal, air-heating vaporizer type is usually used. To contribute in developing an efficient air-heating vaporizer, experiment on finned type air-heating vaporizer using 8 fins, 50 mm fin length (8fin50le) and 4 fins, 75 mm fin length (4fin75le) were conducted. The experiments were conducted by varying the ambient condition and the length of the vaporizer. The ambient air was controlled so that it has the same condition with air condition in every season available. $LN_2$ is used to substitute LNG because of safety reason. The results show that characteristics of the finned type 4fin75le vaporizer are comparable to finned type 8fin50le vaporizer at all season.

This study investigates on the inhalation force improvements of hood consisted one of the local ventilation systems attached the new device named as gas-guidance-device for removed a pollution source. The numerical method applying finite element method is calculated the velocity and pressure distributions of a moving fluid at the beginning and the inside of a hood with and without the gas-guidance-device in hood. And, the experimental study is measured the wind velocity using the anemometer at the same condition of numerical study. Also, the optimum shape of gas-guidance-device which is suitable for hood shape derived from the numerical and experimental results. The results of this study is supplied the important data to an industrial field for control of a pollution source in the engineering aspect. Moreover, the introduced technique of hood attached the gas-guidance-device is very useful to remove the harmful materials such as dust and waste happened in the manufacturing factory.

The practical studies on the method of particle contamination control for yield enhancement in the cleanroom were carried out. The method of the contamination control was considered, which is composed of data collection, data analysis, improvement action, verification, and implement control. The composition analysis for data collection and data analysis was used in the cellular phone module packaging lines. And this method was evaluated by the variation of yield loss between before and after improvement action. In case that the composition analysis was applied, the critical sources were selected and yield loss reduction through improvement actions was also investigated. From these results, it is concluded that the composition analysis is effective solutions for particle contamination control in the cleanroom.

The practical studies on the method of particle contamination control for yield enhancement in the cleanroom were carried out. The method of the contamination control was proposed, which are composed of data collection, data analysis, improvement action, verification, and implement control. The partition check method for data collection and data analysis was used in the cellular phone module production lines. And this method was evaluated by the variation of yield loss between before and after improvement action. In case that the partition check method was applied, the critical process step was selected and yield loss reduction through improvement actions was observed. From these results, it is concluded that the partition check method is effective solution for particle contamination control in the cleanroom production lines.

To generate negative ion, a small dielectric barrier discharge (DBD) plasma reactor was used in this study and operated by high AC voltage. With increasing of voltage, we can get more negative ions. However unfortunately, if the input voltage is too high, it will also cause formation of ozone which is very harmful to human being health. So the work of finding out the best condition of Voltage and frequency was carried out firstly. After several times of measurement, operating at 20 kHz frequency is the best condition generating high ion concentration without ozone. For the purpose of finding out the best reactor structure, two types of surface dielectric barrier discharge (DBD) reactors were examined to produce negative oxygen ions at the conditions of 20 kHz frequency. The results indicated that the surface DBD reactor with several small tips showed better characteristics for generation of negative oxygen ions at the same condition.

In this work, combustion-synthesized $TiO_2$ nanoparticles were used for the photocatalytic degradation of methylene blue with UV light irradiation. Also the results were compared with those of commercial $TiO_2$ nanoparticles (Degussa, P-25). Particle characteristics of the two were analyzed thru the SEM, TEM, and XRD. In spite of the lower specific surface area than that of P-25, the $TiO_2$ nanoparticles formed in this study showed the relatively good ability to degrade the concentration of the organics.

Influence of changing combustor pressure on flame stabilization and emission index in the swirl-stabilized flame was investigated. The combustor pressure was controlled by suction fan at combustor exit. Pressure index ($P^{\ast}$=Pabs/Patm), where Pabs and Patm indicated the absolute pressure and atmosphere pressure, respectively, was controlled in the range of 0.7${\sim}$1.3 for each equivalence ratio conditions. The flammable limits of swirl flames were largely influenced by changing combustor pressure and they showed similar tendency with laminar flames. NOx emission index decreased with decreasing pressure index for overall equivalence ratio conditions. R.m.s. of pressure fluctuations is increased with decreasing combustor pressure. This flame fluctuation caused incomplete combustion, hence CO emission index increased. These oscillating flames were measured by simultaneous $CH^{\ast}$ chemiluminescence time-series visualization and pressure fluctuation measurement.

A new microemitter (microcombustor) configuration for a micro thermophotovoltaic system in which thermal energy is directly converted into electrical energy through thermal radiation was investigated experimentally and computationally. The microemitter as a thermal heat source was designed for a few watt power-generating micro thermophotovoltaic system. In order to satisfy the primary requirements for designing the microemitter, i.e., stable burning in the small confinement and maximum heat transfer through the emitting walls but uniform distribution of temperature along the walls, the present microemitter is cylindrical with an annular-type shield for heat recirculation to apply for the excessive enthalpy concept. Results show that the heat recirculation substantially improves the performance of the microemitter: the observed and predicted thermal radiation from the microemitter walls indicated that heat generated in the microemitter is uniformly emitted.

Spray combustion characteristics of a conducting fuel electrospray have been studied for clean combustion technology. The multiplexing system which can retain the characteristics of the cone-jet mode is inevitable for the electrospray application. Charged micro droplets can be obtained in almost uniform size during operating the electrospray in the cone-jet mode. This experiment device set up the multiplexed grooved nozzle system with the extractor. Using the grooved nozzle, the stable cone-jet mode can be achieved at the each groove in the grooved mode. This electrospray system was applied to the diffusion combustion. It is the first step to discover the diffusion combustion characteristics of the electrospray. In case of the single grooved nozzle electrospray, the diffusion flames are occurred at each jet of grooved mode and they are quite stable. The exhaust gas analysis was indicated that there is the critical point which can make very stable diffusion combustion.

Characteristics of turbulent lifted flames in coflow jets with the varying initial temperature have recently been investigated about only propane case diluted by nitrogen. The investigation has firstly improved a premixed flame model and a large scale mixing model among competing theories on the stabilization mechanism of turbulent flame to be suitable for a high temperature condition. In this research, about methane with good availability to apply for a practical combustor as clean fuel, its characteristics of turbulent nonpremixed flame have been studied experimentally. The results have shown an effectiveness of the premixed flame model and the large scale mixing model considered initial temperature variation. Additionally, considering the axial distance where the mean fuel concentration falls below the stoichiometric level along the center line of the jet according to diluting nitrogen, the premixed flame model have more accurately been improved.

Direct Numerical Simulation was carried out to predict mass transfer in turbulent flow around a rotating stepped cylinder. This investigation is a follow-up study of Nesic et al. [Corrosion, Vol. 56, No. 10, pp. 1005 - 1014] The original motivation of this work stemmed from the efforts to design a simple device which can generate flows of high turbulence intensity at low cost for corrosion researchers. Two cases were considered; Sc=1 and 10 both at Re=335. Here, Sc and Re stand for Schmidt number and Reynolds number, respectively, based on the step height and the surface speed of the cylinder upstream the step. Main focus was placed on the correlation between turbulent fluctuation and concentration field. The spatio-temporal evolution of concentration field is discussed. The numerical results are qualitatively compared with those of the experiment conducted with the same flow configuration.

The 2nd order response surface method (RSM) has been carried out to get optimum thermal design for enhanced heat transfer on square channel with bleed holes. The RSM was used as an optimization technique with Reynolds-averaged navier-stokes equation. Turbulence model for heat transfer analysis used RNG k-epsilon model. The wall function used enhanced wall function. Numerical local heat transfer coefficients were similar to the experimental tendency. Two design variables such as attack angle of rib (${\alpha}$), rib pitch-to-rib height ratio (p/e) were chosen. Operation condition considered bleeding ratio per bleed hole ($BR_{hole}$). A response surface were constructed by the design variables and operation condition. As a result, adjusted $R^2$ was more than 0.9. Optimization results of various objective function were similar to heat transfer in channel with and without bleed flow. But friction factor was lower than channel without bleed flow.

The present research was experimentally performed to analyze the effect of delta-wing vortex generators(DWVG) on the heat transfer of fin surface of the plate fin-oval tube. The local heat transfer coefficient of the fin surface for four kinds of DWVG's arrangement was measured by the naphthalene sublimation technique for Reynolds numbers ranging from 2000 to 3200. The results showed that the heat transfer of the plate fin-oval tube can be significantly enhanced by DWVG for relatively low Reynolds numbers.

Freon CHF experiments are carried out to investigate the CHF enhancements by mixing vane shapes of spacer grids in nuclear fuel assembly. The experiments were performed for a wide range mass flux, 50$\sim}$3000 kg/$m^2s$. Three kinds of spacer grids in 5${\times}$5 rod bundles are tested: no mixing vane grids, hybrid mixing vane grids, and split mixing vane grids. The CHF performances are compared along with the data belong to the PWR operating conditions based on a water equivalence through a fluid-to-fluid modeling method. The average of the data in this range is 16.4% for 37 data of hybrid vane grid and 12.5% for 24 data of split vane. In the lower mass flux, however, the split vane grid shows slightly higher performance than the hybrid vane grid.

Heat transfer distributions and friction factors in square channels (3.5 ${\times}$ 3.5 cm) with twisted tape inserts and with twisted tape inserts plus interrupted ribs are respectively investigated. Tests are performed for Reynolds numbers ranging from 8,900 to 29,000. The rib height-to-channel hydraulic diameter, e/Dh, is kept at 0.057 and test section length-to-hydraulic diameter, L/Dh is 30. The twisted tape is 0.1 mm thick carbon steel sheet with diameter of 3.3cm, length of 90cm, and 2.5 turns. The square ribs are arranged to follow the trace of the twisted tape and along the flow direction defined as axial interrupted ribs. Each wall of the square channel is composed of isolated aluminum sections. The following conclusions from the experimental study were drawn as: 1) In the 4 heating wall channel with twisted tape inserts, Nusselt number based on bottom wall temperature is enhanced by 1.2 - 1.6 times if adding the axial interrupted ribs on the bottom wall only. 2) The twisted tape with interrupted ribs under the two-sided heating condition produces the highest heat transfer performance. 3) Friction factor data obtained for the square channel with twisted tape inserts plus axial interrupted ribs are less than those in the past publications for circular tubes with axial interrupted ribs and twisted tape inserts.

Since the heat generation in a chip increases as the components are miniaturized and the computing speed becomes faster, suitable heat dissipation has become one of the primary limiting factors to ensure the reliable operation of the electronic devices. A pin-fin array could be used as an alterative cooling system of the electronic equipment. In this study, convective heat transfer through the pin-fin array is analyzed experimentally based on porous medium approach. The influence of the structure of the pin-fin array including the pin-fin spacing, the pin diameter and plate length on heat transfer characteristics is investigated and compared with the previous analytical results and existing correlation equations. Nowadays, electronic and mechanical devices become smaller and smaller. In this sense, the main purpose of this study is to decide the optimum pin-fin arrangement to get similar heat transfer performance when the length of the existing cooling system is reduced as a half.

A rectangular fin with variable fin height, fin length and surrounding temperature is analyzed using a one-dimensional analytical method. Both the heat loss from a rectangular fin with non-insulated fin tip and that with insulated fin tip are presented as a function of the fin height, fin tip length and the convection characteristic number. The relative error in the heat loss of these two cases is also given as a function of the same variables. One of the results shows that the trend of heat loss for both cases with the variation of given variables is similar even though the relative error increases as the shape of the fin becomes shorter and fatter.

Heat transfer from three-dimensional heat-generating modules was investigated. A simulated electronic module in an array configured with dummy module elements was used to measure the average heat transfer coefficients. Various module arrangements were tested using module spacings of 0.85 and 1.15 cm for six Reynolds numbers ranging from 500 to 975. The results show that a module placed in-line with and upstream of a heated module results in the heat transfer enhancement due to a high level in turbulence prompted by upstream modules. The highest enhancement occurs when the separation distance between modules is close to the module length in the flow direction. Flow visualization reveals laminar flow on the front of the first module, slow recirculation regions on the sides parallel to the air stream, and turbulence on the back side. It appears that the first module serves to trip the air stream and produce a high level of turbulence, which enhances the heat transfer rate downstream.

The microchannel waterblock has a good capability in the cooling of electronic devices. The object of this paper is to estiblish the scheme of design for the microchannel waterblock. The effects of flow rate and channel size on the cooling performances are investigated. It was found that the optimum flow rates were ragned from 0.7 lpm to 1.4 lpm. The thermal resistance at 2.0 lpm and 100 W was 0.13 $^{\circ}C$/W. Decrease in the width of channels is more effective for the improvement in the cooling performances of microchannel waterblock than increase in the height of channels. The increase of pressure drop resulted from decrease in the width of channels can be decreased by increasing the hight of channels.

The demand of high speed and miniaturization of electronic devices results in increased power dissipation requirement for thermal management. In this work, the effects of microchannel width, height and liquid flowrate on the cooling performances of microchannel waterblock are investigated experimentally. The microchannel waterblock considered ranged in width from 0.5 to 0.9 mm, with the channel height being nominally 1.7 to 9 times the width in each case. The experiments were conducted using water, over a liquid flow rate ranging from 0.2 to 2.0 lpm. The base temperature, thermal resistance and pressure drop increase with increasing of liquid flow rate. The measured thermal resistances ranged from 0.10 to 0.23 $^{\circ}C$/W for the channel 5.

In order to guarantee the performance of electronic products. It is needed to improve the cooling performance of heat sink. So this paper has been made to investigate the cooling performance for the aluminum heat sink using pulsating heat pipe(PHP). The pulsating heat pipe was used as a heat spreader. Working fluid was R-22. Heater (50 mm ${\times}$ 50 mm ${\times}$ 3mm) was attached to heat sink and it generated 30W, 60W, 80W, 100W. Heat sink was tested for forced convection with 1${\sim}$4m/s of inlet air velocity. And both type heat sinks were carried out by using CFD simulation. This study showed that pulsating heat pipe can be a good tool to improve cooling performance of heat sink.

The multi channel of DDI which is the core part of the LCD-TV has been propelled. When multi channel in DDI is introduced, it brings a thermal problem because of the increased power. To solve the thermal problem of the DDI it needs to be investigated each at the package level and module level. It is important to extract the junction temperature(Tj) of DDI clearly from the system level. The objective of this research is to construct a compact model. The compact model is to reduce LCD module including DDI. When the compact model is used, it will be able to easily handle the boundary condition and accurately predict the temperature. Consequently, the temperature of DDI can be calculated easily at the system level. Through this research,we also proposed the cooling plan of DDI for a protection of overheating. The cooling plan was utilized in DDI design.

Hot spots on thin wafers of IC packages are becoming important issues in thermal and electrical engineering fields. To investigate these hot spots, we developed a Diode Temperature Sensor Array (DTSA) that consists of an array of 32 ${\times}$32 diodes (1,024 diodes) in a 8 mm ${\times}$ 8 mm surface area. To know specifically the hot spot temperature which is affected by the chip thickness and a generated power, we made the DTSA chips, which have 21.5 ${\mu}m$, 31 ${\mu}m$, 42 ${\mu}m$, 100 ${\mu}m$, 200 ${\mu}m$, and 400 ${\mu}m$ thickness using the CMP process. And we conducted the experiment using various heater power conditions (0.2 W, 0.3 W, 0.4 W, 0.5 W). In order to validate experimental results, we performed a numerical simulation. Errors between experimental results and numerical data are less than 4%. Finally, we proposed a correlation for the hot spot temperature as a function of the generated power and the wafer thickness based on the results of the experiment. This correlation can give an easy estimate of the hot spot temperature for flip chip packaging when the wafer thickness and the generated power are given.

The performance of polysilicon thin film transistor (p-Si TFT) has an important role in the operation of active matrix liquid crystal displays. To fabricate the p-Si TFTs that have uniform characteristics, understanding of the recrystallization mechanism of silicon is crucial. Especially, the analysis of the transient temperature variation and the liquid-solid interface motion is required to find the mechanism. The thermal conductivity is one of the most important parameters to understand the mechanism. In this work, a KrF eximer laser beam was irradiated to amorphous silicon thin films. We measured the transient reflectivity at the wavelength of 633 nm. We carried out the numerical simulation of one dimension conduction equation so that we determined the most well-fitted thermal conductivity by comparing the numerically obtained transient reflectivity with the experimentally measured one. The experimentally determined thermal conductivity of amorphous silicon thin films is 1.5 W/mK.

In present study, a temperature field of specimens which was coated with fluorescence dye such as Rhodamine-B(Rh-B) has been measured, based on the fluorescence intensity. Silica(SiO2) nano porous structure with 1um thickness was constructed on a cover glass, and fluorescence dye was digested into these porous thin films. To optimize manufacturing coating process, various solvents, Rh-B concentration, and other chemical materials were applied to fabricate the specimen and all specimens were measured on the various temperature conditions. For the measurement, a 14 bit cooled CCD camera with 1600 by 1200 spatial resolution is equipped with epifluorescence microscope to obtain only fluorescence intensity from 1.2 mm by 0.9 mm field of view of the illuminated coated specimen.

We have analyzed the three-dimensional thermal conduction in anisotropic materials using nonsymmetric-Fourier transforms. And a complete theoretical treatment of the photothermal deflection spectroscopy has been performed for thermal conductivity measurement in anisotropic medium. Thermal conductivity tensor was determined by the deflection angle and phase angle with the relative position between the heating and probe beams. The influence of the parameters, such as modulation frequency of the heating beam, the thermal conductivity tensor, was investigated.

Clathrate compound is the material that host in hydrogen bond forms cage and guest is included into it and combined. Crystallization of hydrate is generated at higher temperature than that of ice from pure water. And physical properties according to temperature are stable and congruent melting phenomenon is occurred without phase separation. But clathrate compound still had supercooling problem occurred in the course of phase change and supercooling should be minimized because it affects efficiency of equipment very much. Therefore, various studies on additives to restrain this or heat storage methods are needed. In this study was investigated the cooling characteristics of the TMA-water clathrate compound including TMA (Tri-methyl-amine, $(CH_3)_3N)$ of 20${\sim}$25 wt% as a low temperature latent heat storage material. And ethanol$(CH_3CH_2OH)$ was added and its cooling characteristics were studied experimentally to restrain supercooling of TMA-water clathrate compound.

An inverse radiation analysis is presented for the estimation of the wall emissivities for an absorbing, emitting, and scattering media with diffusely emitting and reflecting opaque boundaries. In this study, a repulsive particle swarm optimization(RPSO) algorithm which is a relatively recent heuristic search method is proposed as an effective method for improving the search efficiency for unknown parameters. To verify the performance of the proposed RPSO algorithm, it is compared with a basic particle swarm optimization(PSO) algorithm and a hybrid genetic algorithm(HGA) for the inverse radiation problem with estimating the wall emissivities in a two-dimensional irregular medium when the measured temperatures are given at only four data positions. A finite-volume method is applied to solve the radiative transfer equation of a direct problem to obtain measured temperatures.

Numerical simulations are conducted for the analysis of a thermal mass air flow sensor with periodic heating pulses on silicon-nitride ($Si_3N_4$) thin membrane structure. This study aims to find the locations of temperature sensors on the thin membrane and the heating pulse conditions, that the higher sensitivity can be achieved, for the development of a MEMS fabricated mass air flow sensor which is driven in periodic heating pulse. The simulations, thus, focus on the membrane temperature profile according to variation of the flow velocity, heating duration time and imposed power. The flow velocity of the simulations is ranging from 3 m/s to 35 m/s, heating duration time from 1 ms to 3 ms and imposed power from 50 mW to 90 mW. The corresponding Reynolds numbers vary from 1000 to 10000.

A gap conductance is very important factor which can affect nuclear fuel temperature. Especially, in case of an annular fuel, a gap conductance effect can lead an unexpected heat split phenomena which is caused by a large difference of an inner and outer gap conductance. The gap conductance mechanism is very complicated behavior due to the its strong dependency on microscopic factors such as a contact surface roughness, local contact pressure and local temperature. In this paper, for the decision of test temperature and pressure range, a procedure and calculation results of in-reactor fuel temperature and pressure analysis are summarized which can be applied to test equipment design and determination of test matrix. Based upon analysis results, it is concluded that the minimum and maximum test temperature are $300^{\circ}C$ and $530^{\circ}C$ respectively, and the maximum pellet/cladding interfacial contact pressure should be observed up to 45MPa.

Self-blast circuit breakers utilize the energy dissipated by the arc itself to create the required conditions for arc quenching during the current zero. The high-current simulation provides information about the mixing process of the hot PTFE cloud with $SF_6$ gas which is difficult to access for measurement. But it is also hard to simulate flow phenomenon because the flow in interrupter with high current, $SF_6$-PTFE mixture vapor and complex physical behavior including radiation, calculation of electric field. Using a commercial computational fluid dynamics(CFD) package, the conservation equation for the gas and temperature, velocity and electric fields within breaker can be solved. Results show good agreement between the predicted and measured pressure rise in the thermal chamber.

The vacuum interrupter (VI) is used for medium-voltage switching circuits due to its abilities and advantages as a compact and environmental friendly circuit breaker. In general, the application of a sufficiently strong axial magnetic field (AMF) permits the arc to be maintained in a diffused mode to a high-current vacuum arc. A full understanding of the vacuum arc physics is very important since it can aid to improve the performance of vacuum interrupter. In order to closely examine the vacuum arc phenomena, it is necessary to predict the magnetohydrodynamic (MHD) characteristics by the multidisciplinary numerical modeling, which is coupled with the electromagnetic and hydrodynamic fields, simultaneously. In this study, we have investigated the electromagnetic behaviors of high-current vacuum arcs for two different types of AMF contacts, which are coil-type and cup-type, using a commercial finite element analysis (FEA) package, ANSYS. The present results are compared with those of MAXWELL 3D, a reliable electromagnetic analysis software, for verification.

The coupled simulation is performed to find out the interaction of arc plasmas with surrounding materials in a thermal puffer plasma chamber. In order to be more realistic, PTFE nozzle ablation and Cu electrode evaporation, which are caused by high temperature of arc plasmas, are considered together. Pressure rise and temperature inside the chamber generated during the whole arcing history are investigated with the applied currents. It is very important to define how thermal flow and mass transfer are processing between the arc plasma and surrounding materials for further understanding complex physics inside the chamber. It is concluded that the result might be very useful to understand the mechanism happened inside and to design thermal puffer plasma chambers, but further experimental studies are required to verify the results for the more practical applications.

Effect of natural convention for hot plate surface temperature uniformity was studied by experiments that were adjusted height of chamber and temperature difference. The hot plate chamber is composed of the hot plate and the upper heater and adiabatic vertical wall. The hot plate diameter is 220mm and maintains temperature at $150^{\circ}C$. Flow pattern compares with surface temperature and confirms that natural convection affects on temperature uniformity of hot plate surface. In case, temperature non-uniformity of hot plate surface is due to heater pattern, lots of weak and small flow cells more improve temperature uniformity than stronger flow cells or non-developing flow cell. Improve temperature uniformity $1.2^{\circ}C$ when developing weak and small flow cells.

In this paper, thermo-flow characteristics of an outer-rotor type of a BLDC motor are numerically analyzed using three-dimensional turbulence modeling. In an advance design of BLDC motor, cooling blades and holes are preferred for the enhanced cooling performances. Rotating the blades and holes generates axial air flow passing through stator slots, which cools down stator by forced convection. For the present study, a new design of the BLDC motor has been developed and major design parameters such as the arrangement of cooling holes, the area of cooling holes, and cooling blades and the cooling blade angle, are analyzed for the enhanced convective heat transfer rate. It is found that the convective heat transfer rate of the new BLDC motor model is increased by about 8.1%, compared to that of the reference model.

Solid propellant gas generators (SPGG) play a role as a turbopump starter in liquid propellant propulsion systems by supplying pressurized gas to power turbines for engine start. For such a purpose, the propellants should burn with a relative low flame temperature and the combustion gas should not contain corrosive constituents such as chlorine compounds. In accordance with these requirements, stabilized AN-based propellants have been usually used as the most appropriate oxidizer for propellant compositions. However, the burning area of the propellant intends to increase to satisfy the required mass flux because of its low burning rate. Consequently the burning area incensement brings on the SPGG size augmentation. A flow restriction such as filters is applied to decrease the SPGG size by rising up the combustion pressure resulting in increasing the burning rate. The feasibility of the size reduction of SPGG by the employment of filters have been studied. The preliminary results of this study show that the considerable reduction of SPGG size would be achievable just by installing a filter with relatively high pressure loss coefficient.

Thermoforming is one of the most versatile and economical processes available for the manufacturing polymer products. The drawback of thermoforming is difficult to get uniform thickness of final products. For the distribution of thickness strongly depends on the temperature distribution of sheet, the adjustment of heater power is very important In this paper, an optimization study for getting uniform temperature distribution was carried out using dual optimization steps. At first, the steady state optimal distribution of heater power is searched by numerical optimization to get uniform temperature of sheet surface. In the second step, time-dependent optimal heater inputs have been found out to decrease the temperature difference through the direction of thickness using Rseponse Surface Method and D-optimal method. The optimization results show that the time-dependent optimal heater power distribution gives acceptable uniform sheet temperature in the field of forming temperature..

Launch vehicles are exposed to aerodynamic heating conditions while flying at high Mach numbers in the atmosphere. In this study aerodynamic heating test for fairing nose-cone was done using ATSF(Aerodynamic Thermal Simulation Facility) and Engineering Model for fairing. ATSF is a facility that can simulate given temperature profile using about 4,000 halogen heaters on fairing model. Aerodynamic heating profile is got from result of thermal analysis using MINIVER, Thermal Desktop and SINDA/FLUINT. After aerodynamic heat test, it is found that initial temperature of fairing inner surface and thickness of BMS has important effects on temperature of fairing inner surface. Also it is confirmed that maximum temperature of fairing nose-cone inner surface during flight is lower than allowable temperature limit. Later, thermal correlation between thermal analysis and experimental results will be done using aerodynamic heating test result

A computational study of BLDC motor is presented to elucidate thermo-flow characteristics in winding and bearing with heat generation. Rotation of rotor and blades drives influx of ambient air into the rotor inlet and the inflow rates are predicted more at the front-side inlet than at the rear-side, which can be ascribed to the different pressure distribution. Recirculation zone appears in the tiny interfaces between windings, however, showing the enhanced cooling performance due to the higher velocity distribution near the rotor wall. In contrast, flow separation and incline angle of bearing groove, and relatively slower velocity distribution cause poor cooling performance and therefore the redesign of the bearing groove is significantly required.

Weather Strip(W/S) is a rubber part to proof water, sound and dust for opening and shutting devices including vehicle doors. And it requires high dimension precision and durability to proof water, noise, vibration and etc. But ironically it itself makes some wind noise because of some protuberance with glasses. The air flow analysis of door part of vehicle makes it possible to calculate and find out the cause of wind noise. In previous analysis, we focus on the numerical air flow analysis of the automobile side part. We do 2D-C.F.D first and 3D second. Through simulations, we can calculate the amount of sound pressure level at the glass run and find out the effects of glass run to make wind noise. Finally we can improve shape of glass run to reduce wind noise although it is small amounts of sound pressure reduction compared with total vehicle noise level.

Shape optimization of an internal cooling passage with staggered dimples on single surface is performed and performances of surrogates are evaluated in this paper. Optimizations are performed so that turbulent heat transfer can be enhanced compromising with pressure loss due to friction. The three-dimensional governing differential equations have been solved to find the overall Nusselt number and friction factor which are related to the objective functions of this problem. Three design variables were selected among the dimensionless geometric variables. Basic surrogate models such as second order polynomial response surface approximation (RSA), Kriging meta-modeling technique, radial basis neural network (RBNN), and derived press based averaged (PBA) surrogate model are constructed. The optimal points are searched from the above constructed surrogates by sequential quadratic programming (SQP). It is shown that use of multiple surrogates can increase the robustness in prediction of better design with minimum computational cost.

In the present work, flow characteristics analysis has been performed for steam turbine bypass control valve (single-path type). The numerical analysis is performed by solving three-dimensional Reynolds-averaged Navier-Stokes (RANS) equations. Shear stress transport (SST) model is used as turbulence closure. Symmetry condition is applied at the mid plane of the valve while adiabatic condition is used at the outer wall of the cage. Grid independency test is performed to find the optimal number of grid points. The pressure and temperature distributions on the outer wall of the cage are analyzed. Mass flow rate at maximum plug opening condition is compared with the designed mass flow rate.

TFT LCD industries have endeavored to adopt the piezo Drop-On-Demand (DOD) inkjet printing technology to production lines of TFT LCD color filters and these efforts have significantly been based on experimental works. Because the degree of complexity in the piezo DOD inkjet printing technology results in too many combination of parameters, the matrix experimental method to investigate all possible sets of parameters becomes ineffective and hence the basic understanding of the piezo Drop-On-Demand inkjet print technology becomes important. In this study, one possible cause of the droplet volume variation across nozzles, which might cause visible swathe marks on an inkjet patterned TFT LCD color filter, is theoretically investiaged and new R&D directions are suggested.

The suppression of fluid force acting on a square prism near plane wall was studied by attaching fences on the corners of the prism. The height of the fence was 10% of the square width and the range of Reynolds number considered was Re=$2.0{\times}10^4$. The experimental parameters were the attaching position and numbers of fences, the space ratios G/B(G/B=0.1${\sim}$1.2) between prism and plane wall. The average drag coefficients were increased and the average lift coefficients were decreased and increased with the space ratios toward plane wall. The drag of the prism was reduced average 7.6% with the space ratios by attaching the normal fence at the rear and upper corner and the horizontal normal fence at the rear and lower corner on the prism. In this case, the separated flow at the front and upper corner was reattached on the upper side of the prism and the vortex streets between the prism and plane wall were appeared more slowly than that of prism without fences.

In order to investigate the vortex body frame interaction around the side mirror of a passenger car, velocity vector fields in the wake, pressure distributions and boundary layer flows over both the mirror surface and the mirror housing, have been measured by several experimental tools. It was resulted that only within an half downstream distance of the mirror span there appears the recirculation zone, and also found that vortex trail towards to the driver side window between A and B pillars, making the acoustic noise and vibration. Wake vortex rolls up after this recirculating zone and makes the trail of the vortex center towards the driver side window, which was also confirmed by measurements of wake velocity vectors in the vertical sections of the trail and visualization over the side mirror surfaces as well. It was also observed that total pressure distribution over the mirror surface has the minimum peak near the lower tip region which can be considered as the origin of the vortex center. It can be concluded that the geometrical modification of the lower tip and the upper root area of the mirror housing is the key to control the wake vortex.

A digital micro holographic particle tracking velocimetry (HPTV) system consisting of a high-speed camera and a single laser with acoustic optical modulator (AOM) chopper was established. The digital micro HPTV system was applied to water flow in a micro curved-tube for measuring instantaneous 3-D velocity field data consecutively. The micro curved-tube is using to reproduce the dorsal aorta or utilize in various lap-on-a-chip. The temporal evolution of a three-dimensional water flow in the micro curved-tube (the curvature, ${\kappa}$=1/${\phi}$, 2/${\phi}$, 4/${\phi}$, 8/${\phi}$) of 100 ${\mu}m$ and 300 ${\mu}m$ inner diameters was obtained and mean velocity field distribution was obtained by statistical-averaging the instantaneous velocity fields.

In this paper, a technique of simultaneously measuring the velocity and the temperature in micro-scale flow is proposed. This method uses particle tracking velocimetry (PTV) for measuring the velocity and laser induced fluorescence (LIF) for measuring the temperature. To measure the accurate velocity and temperature, images for PTV and for LIF are separated by using two light sources and a shutter which is synchronized with a camera. By using only one camera, measurement system can be simplified and the error from complicate optical system can be minimized. Error analyses regarding the concentrations of fluorescent dye and particle and the light source fluctuation are also conducted. It is found that the error of the temperature and the velocity highly depends on the concentration of fluorescent particles which are used for PTV. This technique is applied to the simultaneous measurement of the velocity and the temperature in the electrokinetic flow. It is found that the velocity and temperature vary with the electric field strength and the concentration of electrolyte.

The effect of incidence angle on the three-dimensional flow and aerodynamic loss in the downstream region of a high-turning turbine rotor blade has been investigated with a straight miniature five-hole probe. The incidence angle is changed to be +10, +5, 0, -10, -20, -30 and -40 degrees. The results show that the positive incidence reinforces the three-dimensional vortical flows within the turbine passage including the passage vortex, but the negative incidence weaken them significantly. A small increment in the positive incidence angle results in a remarkable aerodynamic loss increase, while increasing the incidence angle in the negative range leads to a very small change in the aerodynamic loss.

In this paper, the performance analysis is experimentally carried out in order to select the best cooling fan and shroud considering both cooling performance and noise reduction. 4 cooling fans have been tested in the fan tester and the real excavator. In order to obtain the performance of the cooling fans, flow capacity has been estimated by measuring flow velocity using a hot wire anemometer, and noise radiation has been also measured to estimate the fan noise. Characteristics of a box-type and a streamlined shroud have been examined by changing the immersion depth of cooling fans. Based upon the results, the best cooling fan is selected. Finally, the criterion to select the best cooling fan has been set up.

Flows in the centrifugal compressor volute with circular cross section are numerically investigated. The computational grid for the calculation utilized a multi-block arrangement to form a butterfly grid and flow calculations are performed using commercial CFD software, CFX-TASCflow. The centrifugal compressor of this study has axial diffuser after radial diffuser because of the shape of inlet duct and installation constraints. Due to this feature the swirling flow pattern is different from the other investigations. The flow inside volute is very complex and three dimensional with strong vortex and recirculation through volute tongue. The calculation results show circumferential variations of the swirl and through flow velocity and pressure distribution. The mechanism deciding flow structure is explained by considering the force balance in volute cross section. And static pressure recovery and total pressure loss are estimated from the calculated results and compared with Japikse model.

Shell & Tube type Oil Cooler is widely used for hydraulic presses, die casting machines, generation equipments, machine tools and construction heavy machinery. Temperature of oil in the hydraulic system changes viscosity and thickness of oil film. They have a bad effect to performance and lubrication of hydraulic machinery, so it is important to know exactly the heat exchanging efficiency of oil cooler for controlling oil temperature. But most Korean manufacturers do not have test equipment for oil cooler, so they cannot carry out the efficiency test of oil cooler and it is impossible to verify its performance. This paper includes information of construction of necessary utilities for oil cooler test and design and manufacture of test equipment. One can select the optimum product by obtaining performance data through tests of various kinds of oil coolers. And also the paper developed a program which can be easily used for design of 2D and 3D drawings of oil cooler.

This paper describes on pressure drop in a circular pipe of refuse collecting system. The flow characteristics inside the circular pipe are analyzed by three-dimensional Navier-Stokes analysis. In numerical analysis, an organic waste is modeled using the data obtained by site survey. Pressure drop obtained by numerical simulation is compared to the value obtained by experimental measurements for the two kinds of pipe; straight and bended type. The pressure drop obtained by numerical simulation has a good agreement with that of experiments. It is noted that the accurate prediction of pressure drop in the waste pipe is very important to determine the performance of turbo blower used in making a suction pressure in the waste pipe. Especially, the pressure drop for an organic waste is analyzed according to the mass flow rate of waste.

This study aims at analyzing the flow characteristics of the electromagnetic pump using a linear induction motor (LIM) for transferring molten metals. The flow characteristics of the pump are simulated by magnetohydrodynamic(MHD) program. In this system, the LIM is used for transferring molten metal by electromagnetic force. The molten metal is treated as the secondary part of the LIM. Since the LIM produces an electromagnetic force in the duct, the molten metal can flow from the furnace to the reservoir. The flow characteristics of the pump are analyzed using MHD program for magnetic field of 0.1[T] in duct. In order to prove the analysis, we made a prototype electromagnetic pump using LIM.

Air-water two phase natural circulation flow in the T-HERMES (Thermo-Hydraulic Evaluation of Reactor cooling Mechanism by External Self-induced flow)-1D experiment has been evaluated to verify and evaluate the experimental results by using the RELAP5/MOD3 computer code. The RELAP5 results have shown that an increase in the coolant inlet area leads to an increase in the water circulation mass flow rate. However, the water outlet area does not effective on the water circulation mass flow rate. As the coolant outlet moves to a lower position, the water circulation mass flow rate decreases. The water level is not effective on the water circulation mass flow rate. As the height increases in the air injection part, the void fraction increases. However, the void fraction in the upper part of the air injector maintains a constant value. An increase in the air injection mass flow rate leads to an increase in the local void fraction, but it is not effective on the local pressure.

Rotary performance and flow resistance induced by each element in suction nozzle of a vacuum cleaner with a double-blade rotary fan are investigated numerically and its relation with flow-induced noise and rotary performance is examined. Flow resistance and vorticity in suction nozzle are calculated and it is found that they are closely related with flow-induced noise. Gap between double blades, adoption of cross-flow fan, enlargement of flow inlet area, and optimization of blade number are tested for noise reduction. Finally, the effects of each method are verified experimentally. It is found that several combinations of the proposed methods can be adopted for noise reduction although the degree of reduction is not much.

The magnetohydrodynamic(MHD) pressure drop along a liquid sodium flow was measured in a rectangular duct under a transverse magnetic field. The test section was made of a 3 mm thick stainless steel SUS304 with a $74{\times}5mm^2$ rectangular flow channel. The range of experimental parameters was roughly B=0${\sim}$0.18T and U=0${\sim}$0.9m/s at around $200^{\circ}C$. The differential pressure was measured by a diaphragm seal-type pressure transmitter filled with a high temperature silicon oil within 0.1MPa. The experimental results show a similar pressure drop with the theoretical estimation according to a change of the flow velocity and the magnetic field.

This study represents numerical study on the thermal and fluid flow characteristics of exhaust gas in a motorcycle muffler. The reference engine was used 124.cc small displacement. Numerical analysis with computational fluid dynamics(CFD) was carried out to investigate the exhaust gas that flow into a motorcycle muffler. The STAR-CD S/W used to three dimensional steady state CFD analysis in a muffler. And than We got the information of static pressure it is used to structural analysis ant the first baffle plate using the commercial CAE code ANSYS workbench. Exhaust gas flow third chamber from frist chamber and running second chamber. A simulation result shows that each chamber of muffler temperature is about 460 K, 445 and 463K and pressure is about 22,000 Pa, 16,000 Pa and 10,000 Pa.

Recently, increasing demand on not only lighter but also extremely mobile battery make micro fuel cell device very attractive alternative. By reducing the size of fuel cell, surface tension becomes dominant factor with minor gravitational effect. Therefore, it is very difficult to detach the $CO_2$ bubble generating on a cathode side in ${\mu}DMFC$ (micro direct methanol fuel cell). The degassing of a $CO_2$ bubble has drawn quite attention especially for ${\mu}DMFC$ due to its considerable effect on overall machine performance. Our attention has been paid to the dynamic behavior of immiscible bubble attached to the one side of the wall on 2D rectangular channel subject to external shear flow. We use Level Contour Reconstruction Method (LCRM) which is simplified version of front tracking method to track the bubble interface motion. Effects of Reynolds number, Weber number, advancing/receding contact angle and property ratio on bubble detachment characteristic has been numerically identified.

Fluid flow and temperature distribution inside a molten zinc bath were investigated by computational fluid dynamics method. Modeling the channel inductor where alternating current of 60Hz was applied, Lorentz force and generated heat were obtained and later supplemented as source terms to momentum and heat equations. The present work validates CFD technique is effectively adopted when the inductor hardware modification or its configuration is considered for the optimum flows.

Turbulent boundary-layer over a micro-riblet film(MRF) was investigated experimentally. The MRF has sharp V-shaped micro scale grooves of $300{\mu}m$ in width and $176.8{\mu}m$ in height. Particle image velocimetry(PIV) system was employed to measure velocity fields of flow over the MRF coated plate. Flow over a smooth plate was also measured for comparison. The PIV measurements were taken in the streamwise wall-normal planes at Re$\theta$= 985 and 2342. Vortex structures of the flow were analyzed by extracting the swirling strength as an unambiguous vortex-identification criterion. As a result the number of spanwise vortices with clockwise(negative) rotation decreases rapidly in the near-wall region(y<0.2h), but decreases slowly in the outer region(0.2h

This work presents microchannel heat sink shape optimization procedure using Kriging method. Design variables relating to microchannel width, depth and fin width are selected, and thermal resistance has been taken as objective function. Design points are selected through a three-level fractional factorial design of sampling method. Navier-Stokes and energy equations for laminar flow and conjugate heat transfer are solved at these design points using a finite volume solver. Solutions are carefully validated with experimental results. Using the numerically evaluated objective function, a surrogate model (Kriging) is constructed and optimum point is searched by sequential quadratic programming. The process of shape optimization greatly improves the thermal performance of microchannel heat sink under constant pumping power.

An investigation of low Reynolds number flow in nozzles and diffusers which are widely used in the valveless micropump is presented. Flow characteristics in the nozzle and diffuser are explained in view of viscous effect and flow oscillation induced by pumping membrane. These calculation results show that the rectification property of valveless micropump is due to a flow separation in the diffuser and the separation is largely originated from the flow oscillation. Under the assumptions of steady flow velocity profile and flow separation in the diffuser, simplified analytical models are provided to see the dependency of rectification on the micropump geometry. Geometric parameters of channel length, nozzle throat, chamber size, and converging/diverging angle are depicted through the analytical models in low Reynolds number flow, and the prediction and experimental results are compared. This theoretical study can be used to determine the optimum geometry of valveless micropump.

The interaction between metal molecules and liquid metal molecules was modeled in the molecular scale and simulated by the molecular dynamics method in order to understand behaviors of the cluster on metallic surface in collision process. Lennard-Jones potential had been used as intermolecular potential, and only attraction 때 d repulsion had been used for the behavior of the cluster on the metal surface. As results, the behavior of the cluster was so much influenced by the cluster of liquid metal temperature and function of molecules forces, such as attraction and repulsion, in the collision progress. These results of simulation will be the foundation for the micro fabrication manufacturing by using spray application technology.

Non-equilibrium first order extrapolation boundary condition proposed by Guo et $al.^{(9)}$ proposed has a good application for complex geometries, a second order accuracy and a treatment on non-slip wall boundary condition easily. However it has a lack of the numerical stability from high Reynolds number. Guo et $al.^{(9)}$ substituted the density value of adjacent nodes for the density of boundary nodes. This procedure causes the numerical instability on the boundary. In this paper, we derived a procedure of density extrapolation and compared to previous results.

This work presents numerical optimization for design of a blade stacking line of a low speed axial flow fan with a fast and elitist Non-Dominated Sorting of Genetic Algorithm (NSGA-II) of multi-objective optimization using three-dimensional Navier-Stokes analysis. Reynolds-averaged Navier-Stokes (RANS) equations with ${\kappa}-{\varepsilon}$ turbulence model are discretized with finite volume approximations and solved on unstructured grids. Regression analysis is performed to get second order polynomial response which is used to generate Pareto optimal front with help of NSGA-II and local search strategy with weighted sum approach to refine the result obtained by NSGA-II to get better Pareto optimal front. Four geometric variables related to spanwise distributions of sweep and lean of blade stacking line are chosen as design variables to find higher performed fan blade. The performance is measured in terms of the objectives; total efficiency, total pressure and torque. Hence the motive of the optimization is to enhance total efficiency and total pressure and to reduce torque.

For the analysis of a two-phase flow, the interaction between two phases such as the interfacial momentum or heat transfer is proportional to the interfacial area. So the interfacial area concentration (IAC) is one of the most important parameters governing the behavior of each phase. This study focuses on the development of a computational fluid dynamics (CFD) code for investigating a boiling flow with a one-group IAC transport equation. It was based on the two-fluid model and governing equations were calculated by SMAC algorithm. For checking the robustness of the developed code, the experiment of a subcooled boiling in a vertical annulus channel was analyzed to validate the capability of the IAC transport equation. As the results, the developed code was confirmed to have the capability in predicting multi-dimensional phenomena of vapor generation and propagation in a subcooled boiling.

The purpose of this study is to calculate the behavior of molecules for the generation of homogeneous thin-films in the process of spray deposition. The calculation system was composed of a suface molecular region and droplet molecular region. The thin-film was generated when droplet molecules fell to surface molecules. Lennard-Jones potential had been used as intermolecular potential, and only attraction 때 d repulsion had been used for the behavior of the droplet on the solid surface. As results, the behavior of the droplet was so much influenced by the surface temperature in the spray deposition process. High temperature of surface has higher porosity and larger spread area. It was found that simulation results generally agreed well with previous the experimental results. This simulation result will be the foundation for the deposition processes of industry.

The aim of this paper is further studies to achieve deeper understanding in this field. First investigate the influence of operating conditions and design parameters on the hydrodynamics and the mass transfer properties of a loop reactor. This paper provides a literature review on the ejectors applications in the mixing system. A number of studies are grouped and discussed in several topics such as the background, theory of ejector, mixing characteristics, optimization of the system. Since the high efficiencies reactor using ejector widely used in gas-liquid system, especially in a number of chemical and biochemical processes. This is due to their high efficiency in gas dispersion resulting in high mass transfer rate and low power requirements. Thus ejector has been applied to the mixing system. An investigation on hydrodynamics and mass transfer characteristics of gas-liquid ejector has been carried out using three-dimensional CFD modeling.

HANARO, an open-tank-in-pool type research reactor of 30 MWth power in Korea, has been operating normally since its initial criticality in February, 1995. For the last about ten years, A cooling tower of a secondary cooling system has been operated normally in HANARO. Last year, the cooling tower has been overhauled for preservative maintenance including fills, eliminators, wood support, water distribution system, motors, driving shafts, gear reducers, basements, blades and etc. This paper describes the results of the overhaul. As results, it is confirmed that the cooling tower maintains a good operability through a filed test. And a cooling capability will be tested when a wet bulb temperature is maintained about 28 $^{\circ}C$ in summer and the reactor is operated with the full power.

Field measurements have revealed that the pressure drop over a borehole during drilling of a slim oil well or a well with a long reach can depend significantly on the rotation speed of the drill pipe. An accurate prediction of the annular frictional pressure drop is therefore important for conditions where the annular clearance is small. An experimental study was carried out to study solid-liquid two phase flow in a slim hole annulus. Annular velocities of carrier fluids varied from 0.2 m/s to 1.5 m/s. The carrier fluids which were utilized included tap water and CMC water solutions. Pressure drops and average flow rates were measured for the parameters such as inner-pipe rotary speed, carrier fluid velocity, hole inclination and particle injection rate. For both water and CMC solutions, the higher the concentration of the solid particles is, the larger the pressure gradients become.

Turbo-type molecular drag pumps ( MDPs ) are used in the liquid crystal display ( LCD ), semiconductor and other thin film industries. Siegbahn ( disk-type ) molecular drag pumps are used as high-pressure stages in the hybrid-type turbomolecular pumps, where they can operate in the viscous, the transition and the free molecular flow regime. In this study is performed to investigate the pumping characteristics of three-stage disk-type molecular drag pump ( DTDP ) in the molecular transition flow region. The experiments are measured using five vacuum pressure gauges in the positions for rotors of DTDP. The test is performed with nitrogen gas ( $N_2$ ).

In this paper, incompressible flow over a cylinder near a plane wall using the Immersed Boundary. Finite Difference Lattice Boltzmann Method (IB-FDLBM) is implemented. In this present method, FDLBM is mixed with IBM by using the equilibrium velocity. We introduce IBM so that we can easy to simulate bluff-bodies. With this numerical procedure, the flow past a circular cylinder near a wall is simulated. We calculated the flow patterns about various Reynolds numbers and gap ratios between a circular cylinder and plane wall. So these are enabled to observe for vortex shedding. The numerical results are found to be in good agreement with those of previous studies.

Strong wind flow around a building complex was numerically studied by LES. The original motivation of this work stemmed from the efforts to develop a risk assessment technique for windstorm hazards. Lagrangian-averaged scale-invariant dynamic subgrid-scale model was used for turbulence modeling, and a log-law-based wall model was employed on all the solid surfaces including the ground and the surface of buildings to replace the no-slip condition. The shape of buildings was implemented on the Cartesian grid system by an immersed boundary method. Key flow quantities for the risk assessment such as mean and RMS values of pressure on the surface of the selected buildings are presented. In addition, characteristics of the velocity field at some selected locations vital to safety of human beings is also reported.

When a small tube is dipped into a liquid surface, surface tension forces cause the liquid in the tube to rise vertically against the gravity. When the tube is flexible, hydrostatic pressure difference caused by the capillary flow deforms the tube and the deformation which narrows the flow route changes the rising velocity. We study a simple model of this elastocapillary interaction in the context of the surface-tension-driven vertical rise of a liquid between two long flexible hydrophilic sheets that are held a small distance apart at one end. We provide an analytical theory for the rise rate of the liquid and show that our experiments are consistent with the theory.

A numerical study has been performed on the flow past a two-staged orifice in a rectangular duct. The flow field including the recirculation region behind the orifice was investigated and the pressure drop was calculated. Water was used as a working fluid and the flow was treated as the turbulent flow, of which the Raynolds number was 6000. The main parameters for the pressure drop and the recirculation region were the orifice's inclined angle against the duct, the interval between two orifices, the shape of the orifice's hole having the same area, and the change of the hole position at the same interval. The variation of the flow field was investigated with each parameter. Consequently, it was found that the most dominant parameter influencing the drop of the pressure was the change of the hole position at the same interval between orifices. Especially when the interval between orifices was narrow and the relative position the holes was changed, its effect to the flow field was shown most drastically as a result of this study. The SIMPLER algorithm with FLUENT code was employed to analyze the flow field.

In the present study, we reduce the drag and lift fluctuations of the sphere by providing a linear proportional control. For this purpose, we measure the radial velocity along the centerline in the wake and provide blowing and suction at a part of sphere surface based on the measured velocity. Zero-net mass flow rate is satisfied during the control. This control is applied to the flow over a sphere at Re=300 and 425. We vary the sensing location at $0.8d{\leq}X_s{\leq}1.3d$ and find that the most effective sensing region coincides with the location at which minimum correlation between the lift and sensing-velocity directions occurs. As a result, the lift and drag fluctuations are significantly reduced.

A cyclone design concept named Convex cyclone was developed to reduce pressure losses. Contrary to conventional cylinder-on-con type cyclone, inner wall of Convex cyclone are defined with a continuous curve and it has convex shape body. The discontinuity of inner diameter variation rate of cylinder-on-con type cyclone cause additional pressure loss. Continuous wall of Convex cyclone prevent additional pressure loss. In order to verify Convex cyclone design concept, we make a comparative experiments between Stairmand HE and Convex cyclone. Experimental Convex cyclone designed based on Stairmand HE model, and inner wall are defined with circular arch. The experimental result clearly shows that Convex cyclone can achieve maximum 50% pressure loss reduction with a few percent of collection efficiency drop. In addition, the experimental results indicated the existence of optimum convexity, minimum pressure loss, of cyclone wall.

Recently, as the semiconductor production technology develops, there has been growing interest in the cooling system using micro fluid pump. Among the various types of micro fluid pump, the valve-less diffuser/nozzle has been extensively studied in recent years. However, the flat-walled diffuser/nozzle flow has not been clearly looked into due to its non-linear characteristics. In this paper, the flow characteristics of the flat-walled diffuser/nozzle have been analyzed using similitude model and simulations. Similitude models are designed so that the flow pattern is same as that of 1/10 scale flow by using high viscous fluid as working fluid. The results are compared to the simulations. It is shown that the flow characteristics of 2D simulation are different from 3D simulations at high Re region, and the measured pump efficiency is highly dependent on the pressure difference as well as the channel geometry. From these results, the desirable conditions for the efficient pump is discussed.

In vacuum ejector-diffuser systems where a finite volume secondary chamber is used, the secondary jet exhibits transient characteristics during start-up. A steady state is achieved after some time in which mass entrainment prevails indefinitely inside the ejector, though there is no flow from the secondary chamber. An attempt is made in this work to study the infinite entrainment of secondary jet into the primary jet from a finite secondary chamber, with the help of a computational fluid dynamics method. The present study is also intended to identify the operating range of vacuum ejector-diffuser systems where the steady flow assumption can be applied without uncertainty. The results obtained show that the only condition in which an infinite mass entrainment is possible is the generation of a re-circulation zone near the primary nozzle exit. The flow in the secondary chamber attains a state of dynamic equilibrium at this point. Steady flow assumption is valid only after this point.

The ballistic range has long been employed in a variety of engineering fields such as high-speed impact engineering, projectile aerodynamics and aeroballistics, since it can create very high-pressure states in a short time. Since the operation of the ballistic range includes many complicated processes, each should be studied in detail for the best operation of the device. One of the main processes which have a major influence in its operation is the compression of the driver gas. Most of the studies available in this field hardly discuss this process in detail and thus lack a proper understanding of its effect. In the present study, a computational analysis has been made to investigate the compression process in the pump tube of a ballistic range. The results obtained are validated with some experimental data. It is seen that the pump tube parameters and the piston mass significantly affect the compression process and the time to build up the required diaphragm rupture pressure.

Early detection of a disease is important to tackle treatment issues in a better manner. Several diagnostic techniques are in use, these days; for such purpose and tremendous research is going on to develop newer and newer methods. However, more work is required to be done to develop cheap and reliable early detection techniques. Micro-fluidic chips are also playing key role to deliver new devices for better health care. The present study focuses on a review of recent developments in the interrogation of different techniques and present state-of-the-art of microfluidic sensor for better, quick, easy, rapid, early, inexpensive and portable POCT (Point of Care testing device) device for a particular study, in this case, bone disease called osteoporosis. Some simulations of the microchip are also made to enable feasibility of the development of a blood-chip-based system. The proposed device will assist in early detection of diseases in an effective and successful manner.

This study presents an impedance method of in-vitro characterization of the thrombotic potential of whole blood. Whole blood samples of 0.2 cc were put into a micro-cell with embedded three electrodes immediately after venepuncture at $37^{\circ}C$. Anti-coagulated blood samples were also collected for hematocrit and blood viscosity analyses. The rate of change of electron flow was measured, which indicates the inverse of the thrombotic potential. A sudden decrease in the rate of change of electron flow was observed at a time equal to approximately 110 seconds. This sudden decrease was significantly delayed in anti-coagulated samples. After the sudden decrease, the rate continued to decrease, reaching a minimum value in unadulterated samples while the change in the rate in the anti-coagulated ones was found rather moderate. Based on these preliminary findings, the present method may be of used as a new tool for the diagnosis of the thrombotic potential of whole blood.

In their long journey through the cardiovascular circuit, erythrocytes are coerced to shape transform and assume different shapes on account of varying flow conditions in different blood vessels. The present work was aimed to visualize these erythrocyte shape transformations by an invitro microcirculatory model, and assess multi-shape erythrocyte deformability. The model uses an in-house fabricated, inexpensive disposable micro flow channel to mimic certain invivo conditions and a fast frame video microscopic system for imaging the shape changes in erythrocytes. Results show the multi-shape transformation of erythrocyte christened as discoidal shape, the asymmetrically deformed 'hat' and 'bullet-like' shapes, and the axially deformed 'slipper' and 'spindle-like' shapes. Specific erythrocyte showed the shape transition and transformation while passing through the observed window. The obtained erythrocyte shapes very analyzed for deformability index using image processing techniques that varied significantly (p <0.001) for different shapes as compared with the resting shape.

It has been confirmed that implementation of the no-slip boundary conditions for the lattice-Boltzmann method play an important role in the overall accuracy of the numerical solutions as well as the stability of the solution procedure. We in this paper propose a new algorithm, i.e. the method of the dynamic boundary condition for no-slip boundary condition. The distribution functions on the wall along each of the links across the physical boundary are assumed to be composed of equilibrium and nonequilibrium parts which inherit the idea of Guo's extrapolation method. In the proposed algorithm, we apply a dynamic equation to reflect the computational slip velocity error occurred on the actual wall boundary to the correction; the calculated slip velocity error dynamically corrects the fictitious velocity on the wall nodes which are subsequently employed to the computation of equilibrium distribution functions on the wall nodes. Along with the dynamic selfcorrecting process, the calculation efficiently approaches the steady state. Numerical results show that the dynamic boundary method is featured with high accuracy and simplicity.

Three dimensional numerical studies were carried out to investigate the effect of aspect ratio on gas slip flow in rectangular microchannels. We focused on aspect ratio effect on slip velocity, pressure distribution and mass flow rate. As aspect ratio decreases the wall slip velocity also decreases. As a result nonlinearity of pressure distribution increases. The slip velocities on sides and top/bottom walls are different and this difference decreases with increasing aspect ratio. These two velocities are equal when aspect ratio is 1. The ratios of slip mass flow rate over noslip mass flow rate increases with increasing aspect ratios.

We conducted a numerical study of AC-electroosmotic (alternating current) effect on the fluid flow and mixing in a 3-D microchannel. The microchannel used as an efficient micro-mixer is composed of a channel and a series of pairs of electrodes attached in zigzag pattern on the bottom wall. The AC electric field is applied to the electrodes so that a steady flow current takes place around the electrodes. This current is flowing across the channel and thus contributing to the mixing of the fluid within the channel. We performed numerical simulations by using a commercial code to obtain a steady flow field. This steady flow is then used in evaluation of the mixing performance via the concept of mixing index. It was found that good combination of two kinds of electrode, which gave us a good mixing, is not simple harmonic. And when the length ratio of these two kinds of electrode is 2:1, we can get the best mixing effect.

The purpose of this paper is to study the influence of operating condition and ejector geometries on the hydrodynamics and on the mass transfer characteristic of ejector. The CFD results were validated with available experimental data. Flow field analyses and predictions of ejector performance were also carried out. Variation on the operating conditions was made by varying the gas-liquid flow rate ratio in the range of 0.2 to 1.2. The ejector configuration was also varied on the length to diameter ratio of mixing tube ($L_M/D_M$) in the range of 4 to 10. CFD studies show that at $L_M/D_M$ 5.5, the volumetric mass transfer coefficient increases with respect to gas flow rates. Meanwhile, at $L_M/D_M$ 4, the plot of volumetric mass transfer coefficient to gas-liquid flow rates ratio reach maximum at gas-liquid flow rates ratio of 0.6. This study also shows that volumetric mass transfer coefficient decrease with respect to the increase of mixing tube length.

Numerical simulation on the two-stage centrifugal compressor with Low Solidity Vaned Diffuser LSVD) for HFC-134a Turbo-chiller was performed using a commercial code. The comparative study with experimental results from other compressor was also investigated to testify the simulation schemes. The numerical analysis was separately simulated for each stage of the compressor and the effect of impeller-diffuser flow interaction was considered. Setting angle of the diffuser vane changed in the range of 15 deg. and the effects on its variation were discussed in detail including the flow analysis in the passage of the compressor. The vane setting angle obtained from the preliminary design was slightly adjusted to the optimal value by the performance enhancement in terms of pressure recovery and flow characteristics.

In this study, one-way fluid structure interaction analysis(FSI) on wind turbine blade was performed. Both a quantitative fluid analysis on 3-bladed wind turbine and a structural analysis using the surface pressure data resulting from fluid analysis were carried out. Streamlines and angle of attack was easily acquired from analysis results, we showed the inlet velocity that the stall begins to occur. In the structural analysis, structural displacement and maximum stress of the two comparative models was calculated. The location that has maximum stress was found. The pressure difference between back and front part of the blade increases as the inlet velocity increase. The torque and maximum with regard to inlet velocity was also presented.

Recently, small hydropower attracts attention because of its clean, renewable and abundant energy resources to develop. However, suitable turbine type is not determined yet in the range of small hydropower and it is necessary to study for the effective turbine type. Moreover, relatively high manufacturing cost by the complex structure of the turbine is the highest barrier for developing the small hydropower turbine. Therefore, a cross-flow turbine is adopted because of its simple structure and high possibility of applying to small hydropower. The purpose of this study is to examine the optimum configuration of nozzle shape to further optimize the cross-flow hydraulic turbine structure and to improve the performance. The results show that pressure on the runner blade in Stage 1 and velocity at nozzle outlet have close relation to the turbine performance.

Micro hydraulic power generation of which the output is less or equal to a 100kW is attracting considerable attention. This is because of its small, simple, renewable, and abundant energy resources. By using a small hydropower generator of which main concept is based on using the different water pressure levels in pipe lines, energy which was initially wasted by use of a reducing valve at the end of the pipeline, is collected by turbine in the hydropower generator. A propeller shaped hydroturbine has been used in order to use this renewable pressure energy. In this study, in order to acquire basic design data of tubular type hydraulic turbine, output power, head, efficiency characteristics due to the flow coefficient are examined in detail. Moreover influences of pressure and velocity distributions with the variations of runner vane angle on turbine performance are investigated by using a commercial CFD code.

Flow-induced forces on two identical nearby circular cylinders immersed in the cross flow at Re =100 were numerically studied. We consider all possible arrangements of the two circular cylinders in terms of the distance between the two cylinders and the inclination angle with respect to the direction of the main flow. It turns out that significant changes in the characteristics of flow-induced forces are noticed depending on how the two circular cylinders are positioned, resulting in quantitative changes of force coefficients on both cylinders. Collecting all the numerical results obtained, we propose a contour diagram for drag coefficient and lift coefficient for each of the two cylinders. The perfect geometrical symmetry implied in the flow configuration allows one to use those diagrams to estimate flow-induced forces on two identical circular cylinders arbitrarily positioned in physical space with respect to the main flow direction.

Flow patterns in the presence of two identical nearby circular cylinders at =100 were numerically studied. We considered all possible arrangements of the two circular cylinders in terms of the distance between the two cylinders and the inclination angle with respect to the direction of the main flow. Eight distinct flow patterns were identified based on vorticity contours and streamlines, which are Base-Bleed, Biased-Base-Bleed, Shear- Layer-Reattachment, Induced-Separation, Vortex-Impingement, Flip-Flopping, Modulated Periodic, and Synchronized-Vortex-Shedding. Collecting all the numerical results, we propose a general flow pattern diagram for flows past the two cylinders. The perfect geometrical symmetry implied in the flow configuration allows one to use this diagram to distinguish flow patterns in the presence of two identical circular cylinders arbitrarily positioned in physical space with respect to the main flow direction.

Flow and aerodynamic characteristics were analyzed numerically for a commercial passenger airplane, Boeing 747-400, flying in the cruising condition. The model geometry with 100:1 in scale was obtained by the photo scanning measurement with the maximum error of 1.4% comparing with the real airplane dimension. The three-dimensional inviscid steady compressible governing equations were solved by the finite volume method in the unstructured grid system. The convective terms were treated by the Crank-Nicholson and first-order upwind schemes. In the computational results, the strong wing-tip vortices were clearly observed and the pressure contours on the airplane surface were suggested. The lift and drag forces in the wing with engines increase by 1.49% and 3.9%, respectively compared with the case without engines. The aerodynamic forces were estimated quantitatively for each element which consists of the airplane.

The steady non-reacted compressible flow field in a symmetric micro-thruster, which is used for the accurate attitude control of a satellite, is analyzed varying the nozzle pressure ratio (NPR) to investigate the plume characteristics. The nozzle throat diameter is 0.06 inch and the area ratio is 56. The recirculation region is found just behind the normal shock at the several NPRs due to the locally adverse pressure gradient along the nozzle centerline when the environmental pressure is atmospheric. This phenomenon, the cause of flow loss, is similar to the flow behind a blunt body. As NPR increases the location of Mach disk, characteristics of the normal shock, moves downstream and its strength increases. The Mach number distribution appears in a wave-type patter after the normal shock because oblique shocks are reflected on the shock boundaries especially when NPRs are very high.

In this study, a passive control using a boat-tail device is conducted for a three-dimensional car model in ground proximity. We consider various boat-tails and investigate the mechanism of drag reduction by them. By varying the length and slant angle of boat-tail, we obtain drag reductions up to 40%. From the oil-surface flow visualization and hot-wire measurement, the drag reduction by the boat-tail is characterized by the shear-layer instability and reattachment on the boat-tail, forming a small separation bubble at the upstream part of boat-tail surface, resulting in the delay of main separation and drag reduction. At high slant angles, the flow fully separates and drag is nearly same as that of no control.

The mixing vanes attached to the spacer grid of rod bundles are used to improve the heat transfer in heat exchanger devices by controlling the characteristics of the flow structures and turbulence. In this study, velocity patterns induced by two types of mixing vane(split and swirl vane) are measured by the PIV technique to better understand how to effect on the cross and secondary vortex flow patterns in $5{\times}$ rod bundle simulating the fuel assembly of the nuclear reactor. A successful measurement of the lateral velocity patterns was conducted using a specially designed beam sheet generator and experimental loop at KAERI. As the result, we found that for the cross flow between subchannels, the split vane is more effective than the swirl vane, while for the secondary vortex flow in each subchannel, the swirl vane's one is larger and longer than split vane's one.

Dynamic responses of electrorheological (ER) fluids in steady pressure flow to stepwise electric field excitations are investigated experimentally. The transient periods under various applied electric fields and flow velocities were determined from the pressure behavior of the ER fluid in the flow channel with two parallel-plate electrodes. The pressure response times were exponentially decreased with the increase of the flow velocity, but increased with the increase of the applied electric field strength. In order to investigate the cluster structure formation of the ER particles, it was verified using the flow visualization technique that the transient response of ER fluids in the flow mode is assigned to the densification process in the competition of the electric field-induced particle attractive interaction forces and the hydrodynamic forces, unlike that in the shear mode determined by the aggregation process.

In this study, the effect of particle aggregation on dynamic response time of Electrorheological (ER) fluid is investigated. The particle aggregation time is defined as the time interval between the application of the field and the formation of the first chain bridging the two electrodes. The dynamic response times of an ER fluid sheared between two concentric cylinders have been obtained under two different experimental conditions: the one is that the electric field is induced before shearing, and the other is that the electric field is induced after shearing. From the difference between two response times, the particle aggregation times are determined under various electric fields and shear rates. The experimental results show that the aggregation rate is decreased with an increase of shear rate, while electric field has little effect on it. Therefore, it is verified that the hydrodynamic force hinders the formation of chain-like structures.

Heat transfer on two identical nearby circular cylinders immersed in the uniform cross flow at Re = 100 and Pr = 7.0 was numerically studied. We consider all possible arrangements of the two circular cylinders in terms of the distance between the two cylinders and the inclination angle with respect to the direction of the main flow. It turns out that significant changes in the characteristics of heat transfer are noticed depending on how the two circular cylinders are positioned, resulting in quantitative changes of heat transfer coefficients on both cylinders. Collecting all the numerical results obtained, we propose a contour diagram for averaged Nusselt number for each of the two cylinders. The perfect geometrical symmetry implied in the flow configuration allows one to use those diagrams to estimate heat transfer rates on two identical circular cylinders arbitrarily positioned in physical space with respect to the main flow direction.

A numerical analysis was carried out to investigate the degradation of Volatile Organic Compounds (VOC) in photocatalytic microreactors with different inner geometries. Two different cases of microreactor were considered, namely, one microreactor has bump on the channel and the other has no bump on the channel. The removal efficiency of VOC has been calculated by the Langmuir-Hinshelwood reaction rate equation that was obtained from the experimental results. From the numerical calculations, it was observed that the conversion ratio of VOC for the microchannel with bump is about 4.5% greater than the microchannel without bump. And the mass transfer characteristics in the microreactor are also shown in numerical results. These results can be used effectively for the photocatalytic numerical analysis.

A solver for icing and condensation of water has been developed. The phase change process was solved by the enthalpy method. For the code validation, the temperature and the phase change from water to ice of the driven cavity were calculated. Also, the melting process of the frost on the windshield glass of an automobile has been simulated. The calculation showed a good agreement with analytical solution and other numerical results. Using the present validated code, the condensation of water vapor has been first tried. The computed results provided some physical features of condensation phenomena even though experimental data and other numerical data were not available. For future work, it is recommended to throughly investigate the effects of boundary conditions on the solution.

The cooling water system for the PEFP 20 MeV proton accelerator was established and tested to obtain the precise resonance frequency of DTL through the temperature control of cooling water. The water temperature in the main flow loop was manipulated by adjusting the proportion of hot water returning from the DTL structures through the heat exchanger loop. Due to low duty factor operation and insufficient cooling loop installation of the DTL tanks, the manual mode operation was applied to maintain the DTL temperatures close to their resonance temperatures. The optimized process conditions with flow balancing and pressure drop in the DTL cooling systems are reported.

Turbocharger has been widely used in many passenger cars in application with diesel engines because of high power and fuel efficiency. However, flow-induced noise (whoosh or hissing noise) which is generated within the compressor during its operation at marginal surge line can deteriorate noise characteristics. Hissing noise excitation was associated with the generation of turbulence within the turbocharger compressor and radiated through the transmission path in turbocharger system. In this study, a sharp-edged reactive-type muffler was devised and installed in the transmission path to reduce the hissing noise. Acoustic and fluid dynamic characteristics for the muffler were investigated which is related to the unsteadiness of turbulence and pressure in turbocharger system. A transfer matrix method was used to analyze the transmission loss of the muffler. Simple expansion muffler with extended tube of the reactive type is proposed for the reduction of high frequency component noise. Turbulence computation was carried out by a standard ${\kappa}-{\varepsilon}$ model. An optimal design condition of the muffler was obtained by extensive acoustic and fluid dynamic analysis on the engine dynamometer with anechoic chamber. A significant reduction of the hissing noise was achieved at the optimal design of the muffler as compared with the conventional turbocharger system.

This is a report of a feasibility study on the reduction of harmful substances such as particulate matters and nitric oxides emitted from diesel engines by using a plasma reforming system that can generate hydrogen-rich gas. In this paper, an exhaust reduction mechanism of the non-thermal plasma reaction was investigated to perform its efficiency and characteristics on producing hydrogen-rich gas. Firstly, we explain briefly the chemistry of hydrocarbon reforming. The experimental system is showed in the second part. Finally, we demonstrate the feasibility of producing hydrogen using non-thermal plasma. The experimental results are focused on the influence of the different operating parameters (air ratio, inlet flow rates, voltage) on the reformer efficiency and the composition of the produced gas.

The study of model for velocity and turbulence within the urban canopy was carried out. To evaluate existing urban model we conducted wind tunnel experiment and large-eddy simulation (LES). Mean velocity profile and turbulence are measured within simple three different obstacle arrays. To obtain supplemental data and to verify morphological model large-eddy simulation was performed. Several methods have been used to achieve embodying the flow field in urban area. Recently, morphological method obtaining flow parameters from the statistical or physical representation of obstacle elements is a arising method. It was found that all morphological model, evaluated in this study, over predict the friction velocity, most sensitive one among the flow parameters. Velocity and turbulence in the urban canopy layer were improved by the correction using 'true' friction velocity.

Prediction methods for cavitation noise are presented. At first, direct numerical simulation of cavitating flow noise has been performed, and acoustic analogy equation based on the cavitation noise modeling is derived. For the direct numerical simulation, a density based homogenous equilibrium model is employed to simulate cavitating two-phase flow and the governing equations are solved with high-order numerical schemes to resolve cavitation noise. The compressible Navier-Stokes equations for mixture fluids are discretized with a sixth-order central compact scheme, and the steep gradient of flow variables and supersonic regions are treated with the selective spatial filtering technique. The direct simulation of cavitating flow noise is performed for a 2D circular cylinder at cavitation number 0.7 and 1. The far-field noise is also predicted with the derived analogy equation. Noise spectrum predicted with the equation is well compared with the result of direct numerical simulation and also agree well with the theory.

To diagnose circulatory diseases in the viewpoint of hemodynamics, we need to get quantitative hemodynamic information of blood flows related with the vascular diseases with high spatial resolution of tens micrometer and high temporal resolution in the order of millisecond. For investigating in-vivo hemodynamic phenomena, a new diagnosing technique combining medical radiography and PIV method was newly proposed and developed. This angiographic PIV technique consists of a medical X-ray tube, an X-ray CCD camera, a shutter module for double pulses of X-ray, and a synchronizer. The feasibility of the angiographic PIV technique was tested and quantitative flow velocity field distribution of a flow inside an opaque conduit was acquired by the developed system. It can be used for measuring flow phenomena of nontransparent fluids inside opaque conduits.

A model of the cardiovascular system coupling cell, hemodynamics and autonomic nervecontrol function is proposed for analyzing heart mechanics. We developed a comprehensive cardiovascular model with multi-physics and multi-scale characteristics that simulates the physiological events from membrane excitation of a cardiac cell to contraction of the human heart and systemic blood circulation and ultimately to autonomic nerve control. Using this model, we delineatedthe cellular mechanism of heart contractility mediated by nerve control function. To verify the integrated method, we simulated a 10% hemorrhage, which involves cardiac cell mechanics, circulatory hemodynamics, and nerve control function. The computed and experimental results were compared. Using this methodology, the state of cardiac contractility, influenced by diverse properties such as the afterload and nerve control systems, is easily assessed in an integrated manner.

A 3D human ventricular model is proposed to simulate an integrative analysis of heart physiology and blood hemodynamics. This consists of the models of electrophysiology of human cells, electric wave propagation of tissue, heart solid mechanics, and 3D blood hemodynamics. The 3D geometry of human heart is discretized to a finite element mesh for the simulation of electric wave propagation and mechanics of heart. In cellular level, excitations by action potential are simulated using the existing human model. Then the contraction mechanics of a whole cell is incorporated to the excitation model. The excitation propagation to ventricular cells are transiently computed in the 3D cardiac tissue using a mono-domain method of electric wave propagation in cardiac tissue. Blood hemodynamics in heart is also considered and incorporated with muscle contraction. We use a PISO type finite element method to simulate the blood hemodynmaics in the human ventricular model.

Tumor hemodynamics in vascular state is numerically simulated using pressure node solution. The tumor angiogenesis pattern in our previous study is used for the geometry of vessel networks. For tumor angiogenesis, the equation that governed angiogenesis comprises a tumor angiogenesis factor (TAF) conservation equation in time and space, which is solved numerically using the Galerkin finite element method. A stochastic process model is used to simulate vessel formation and vessel. In this study, we use a two-dimensional model with planar vessel structure. Hemodynamics in vessel is assumed as incompressible steady flow with Newtonian fluid properties. In parent vessel, arterial pressure is assigned as a boundary condition whereas a constant terminal pressure is specified in tumor inside. Kirchhoff's law is applied to each pressure node to simulate the pressure distribution in vessel networks. Transient pressure distribution along with angiogenesis pattern is presented to investigate the effect of tumor growth in tumor hemodynamics.

To control the coating thickness of zinc in the process of continuous hot-dip galvanizing, it is known from early day that the gas wiping through an air knife system is the most effective because of the obtainable of uniformity of coating thickness, possibility of thin coating, working ability in high speed and simplicity of control. But, the gas wiping using in the galvanizing process brings about a problem of splashing from the strip edge for a certain high speed of coating. And, it is known that the problem of splashing is caused mainly by the existence of separation bubble at the neighbor of the strip surface. In theses connections, in the present study, we proposed two kinds of air knife systems having the same expansion rate of nozzle, and the jet structures and coating thicknesses from a conventional and new proposed nozzles are compared. In numerical analysis, the governing equations consisted of two-dimensional time dependent Navier-Stokes equations, standard ${\kappa}-{\varepsilon}$ turbulence model to solve turbulence stress and so on are employed. As a result, it is found that it had better to use the constant rate nozzle from the point view of the energy saving to obtain the same coating thickness. Also, to reduce the size of separation bubble and to enhance the cutting ability at the strip, it is recommendable to use an air knife having the constant expansion rate nozzle.

When the pressure at the weak spot established at a certain part of a high pressure vessel or piping system exceeds a design pressure, this weak spot is burst, and the pressurized gas emitted through the weak spot will cause a compression wave system. In this connection, in the present study, an experimental study by using a conventional shock tube facility is performed to estimate the effects of the material of diaphragm, curvature radius and thickness of materials on the valve opening time in diaphragm. Pressure sensor having 500kHz in natural frequency is installed at 35mm downstream of the rupture diaphragm to measure the static pressure history of propagating and being accumulated compression wave. 4 kinds of materials are used as diaphragm that is aluminium, copper, stainless steel and zinc. The diaphragm radii of curvature R are ${\infty}$, 120mm and 60, respectively. And the depth for $90^{\circ}$ groove is 0.04mm. It is found that the smaller the tensile strength and elongation of the rupture diaphragm is, the smaller the radius of curvature of the rupture diaphragm is, and for the same conditions the thinner the thickness of the rupture diaphragm is, the shorter the valve opening time becomes. Also, the tensile strength, elongation and the radius of curvature of the rupture diaphragm for the same conditions are smaller, the maximum pressure rise caused by the coalescences of the compression wave is smaller. Finally the pressure ratio is higher, the valve opening time is shortened and gradient of pressure increment is more steepen.

In a usual painting process, a liquid drop spreads on canvas by being dragged along a paintbrush. To obtain the fundamental understanding of the painting process from the mechanical point of view, we experimentally investigate various dynamic behavior of a liquid drop that spreads between moving solid plates. It is shown that three distinct types of drop spreading take place, i.e. shearing, spreading, and intact dragging, depending on the liquid viscosity and surface tension, the plate speed, and the wettability. We suggest a regime map based on the capillary number and the receding contact angle, which indicates the boundaries between different types of spreading behavior in a dimensionless space.

The mode change from Taylor cone-jet to dripping in electrospraying has been analytically investigated. The change has been predicted by the dynamic behavior of a liquid drop at the tip of the cone-jet. Conservation laws are applied to determine the upward motion of the drop, and an instability model of electrified jets is used to determine the jet breakup. Finally, for the first time, the analysis enables prediction of the transition in terms of the Weber number and electric Bond number. The predictions are in good agreement with experimental data.

Feedwater heaters of many nuclear power plants have recently experienced severe wall thinning damage, which will increase as operating time progresses. Several nuclear power plants in Korea have experienced wall thinning damage in the area around the impingement baffle - installed downstream of the high pressure turbine extraction steam line - inside number 5A and 5B feedwater heaters. At that point, the extracted steam from the high pressure turbine is two phase fluid at high temperature, high pressure, and high speed. Since it flows in reverse direction after impinging the impingement baffle, the shell wall of the number 5 high pressure feedwater heater may be affected by flow-accelerated corrosion. This paper describes the comparisons between the numerical analysis results using the FLUENT code and the down scale experimental data which effect on disclosing of the shell wall thinning of the high pressure feedwater heaters by porous plate.

In the present study, the numerical investigation on the effects of water-mist characteristics has been carried out for the fire suppression mechanism. The FDS are used to simulate the interaction of fire plume and water mists, and program describes the fire-driven flows using LES turbulence model, the mixture fraction combustion model, the finite volume method of radiation transport for a non-scattering gray gas, and conjugate heat transfer between wall and gas flow. The numerical model is consisted of a rectangular enclosure of $L{\times}W{\times}H=1.5{\times}1.5{\times}2.0m$ and a water mist nozzle that be installed 1.8m from fire pool. In the study, the parameters of nozzle for simulation are the droplet size and the spray velocity. Finally, the droplet size influences to fire flume on fire suppression than spray velocity because of the effect of terminal velocity, and the optimal condition for fire suppression is that the droplet size and the spray velocity are $100{\mu}m$ and 20m/s, respectively.

Gas hydrates are ice-like crystalline compounds that form under low temperature and elevated pressure conditions. Although hydrate formation can pose serious flow-assurance problems in the gas pipelines or facilities, gas hydrates present a novel means for natural gas storage and transportation with potential applications in a wide variety of areas. An important property of hydrates that makes them attractive for use in gas storage and transportation is their very high gas-to-solid ratio. In addition to the high gas content, gas hydrates are remarkably stable. The main barrier to development of gas hydrate technology is the lack of an effective method to mass produce gas hydrate in solid form. The first objective of this study is investigating the characteristics of gas hydrate formation related to several factors such as pressure, temperature, water-to-storage volume ratio, concentration of SDS, heat transfer and whether stirred or not respectively. And the second objective is clarifying the relation between the formation efficiency and each factor in order to find the proper way or direction to improve the formation performance.

We found that there exists a flow inside a droplet in AC electrowetting, which is distinct from DC electrowetting. In order to investigate the origin of the flow inside the droplet, we performed an experiment and numerical simulation. It is conjecture, based on the results of the experiment and numerical simulation, the flow is caused by the so called induced-charge electroosmosis at high frequencies, and by droplet oscillation at low frequencies.

A manipulator is operated for the motion of mechanical hands or arms by mechanical mechanism. When a cup including liquid inside is shifted by a manipulator, it is important to know how a free surface of the liquid moves. In this study, non-dimensional parameters have been found that affect the rise of the free surface in a cup moving with constant acceleration. The non-dimensional parameters are the dimensionless time, the ratio of inertia effect to vicous effect (Reynolds number), aspect ratio of the liquid inside the cup and acceleration ratio (Froude number). Through this study, the height of the free surface rise in a cup has been predicted. Generally the maximum rise of the free surface is dependent on the Reynolds number and Froude number strongly, but on the aspect ratio weakly. But the influence of the aspect ratio on the maximum rise of the free surface in not negligible in the range 10 < Re < 100.

Critical nozzle has been frequently employed to measure the flow rate of various gases, but hydrogen gas, especially being at high-pressure condition, was not nearly dealt with the critical nozzle due to treatment danger. According to a few experimental data obtained recently, it was reported that the discharge coefficient of hydrogen gas through the critical nozzle exceeds unity in a specific range of Reynolds number. No detailed explanation on such an unreasonable value was made, but it was vaguely inferred as real gas effects. For the purpose of practical use of high-pressure hydrogen gas, systematic research is required to clarify the critical nozzle flow of high-pressure hydrogen gas. In the present study, a computational fluid dynamics(CFD) method has been applied to predict the critical nozzle flow of high-pressure hydrogen gas. Redlich-Kwong equation of state that take account for the forces and volume of molecules of hydrogen gas were incorporated into the axisymmetric, compressible Navier-Stokes equations. A fully implicit finite volume scheme was used to numerically solve the governing equations. The computational results were validated with some experimental data available. The results show that the coefficient of discharge coefficient is mainly influenced by the compressibility factor and the specific heat ratio, which appear more remarkable as the inlet total pressure of hydrogen gas increases.

The effectiveness of passive control techniques for alleviating the pressure oscillation generated in a supersonic cavity flow was investigated numerically and experimentally, respectively. The control device includes a sub-cavity installed near the leading edge of a rectangular cavity. Time-dependent supersonic cavity flow characteristics with turbulent features were examined by using the three-dimensional, mass-averaged Navier-Stokes computation based on a finite volume scheme and large eddy simulation. The results show that the pressure oscillation near the trailing edge dominates overall time-dependent cavity pressure variations. Such an oscillation can be attenuated more significantly in the presence of the sub-cavity compared with the cavity without sub-cavity, and a larger sub-cavity leads to better control performance.

An experimental study of the dynamic characteristics of the free rising oblate spherical bubble is investigated. As noted by Saffman, the characteristics of the spiral motion are defined with parameters of the spiral frequency, spiral radius, and rising velocity. High speed camera recorded every detail information of free rising bubble. The spiral number, the bubble rise velocity, and the angular velocities were measured for the bubble size between 1.0mm to 20.0mm in diameter. In order to make clear trajectory, we employed the Fast Fourier Transformation for the normal digital camera data to synchronize with the high speed camera data. It was found that the spiral number suggested here was monotonically decreased as the bubble size increases. The present observation, however tells us that previous Saffman's formulation shows a good agreement with the trend, but over estimated spiral number. Therefore, it is recommended that Saffman's theoretical study is needed to be improved.

In this paper, experimental data on flow pattern transition of inclination angles from 0-90 are presented. A test section is constructed 2 mm long and I.D 1inch using transparent material. The test section is supported by aluminum frame that can be placed with any arbitrary inclined angles. The air-water two-phase flow is observed at room temperature and atmospheric condition using both high speed camera and void impedance meter. The signal is sampled with sampling rate 1kHz and is analyzed under fully-developed condition. Based on experimental data, flow pattern maps are made for various inclination angles. As increasing the inclination angels from 0 to 90, the flow pattern transitions on the plane jg-jf are changed, such as stratified flow to plug flow or slug flow or plug flow to bubbly flow. The transition lines between pattern regimes are moved or sometimes disappeared due to its inclined angle.

To understand the complex phenomena performed in steam explosion, the fast and global measurement of the steam distribution is imperative for this extremely rapid transient stimulation of the bubble breakup and coalescence due to turbulent eddies and shock waves. TROI, the experimental facility requests more robust sensor system to meet this requirement. In Europe, researchers are prefer a X-ray method but this method is very expensive and has limited measurement range. There is an alternative technology such as ECT. Because of TROI's geometry, however, we need axial tomography method. This paper reviews image reconstruction algorethms for axial tomography, including Tikhonov regularization and iterative Tikhonov regularization. Axial tomography method is examined by simulation and experiment for typical permittivity distributions. Future works in axial tomography technology is discussed.

In this research, Rayleigh instability of gas-liquid flow in annular pipe is studied in film boiling using viscous potential flow. Viscous potential flow is a kind of approximation of gas-liquid interface considering velocity field as potential including viscosity. A dispersion relation is obtained including the effect of heat and mass transfer and viscosity. New expression for dispersion relation in film boiling and critical wave number is obtained. Viscosity and heat and mass transfer have a stabilizing effect on instability and its effect appears in maximum growth rate and critical wave number. And the existence of marginal stability region is shown.

A few fluid structure interaction analyses have been developed for turbomachinery blades in comparison with aircraft wings. Also, the existing aeroelastic analyses for turbomachinery blades have been mostly limited to cases with a steady freestream. In reality, however, the inflowing freestream is often pulsating. Therefore, this paper presents stability and forced response analyses of an isolated three-dimensional blade under pulsating freestream conditions. A new three-dimensional unsteady vortex lattice model under a pulsating freestream has been developed in discrete time domain to examine unsteady aerodynamic forces acting on a vibrating blade. The blade's structural behaviors have been analyzed by using a three-dimensional plate model. In the aeroelastic analysis, the flutter onset of a blade under pulsating freestream is predicted by the Floquet analysis. The new time domain method can predict aeroelastic stability as well as time history.

This paper presents a numerical study to enhance the mixing of polyurethane resin in a container. In general, the properties of polyurethane resin vary with the production environment, such as temperature. However in this study we assumed that the dynamic viscosity of the polyurethane is kept constant at 15 [Pa s]. We computed the flow solution and visualized the mixing pattern for different shapes of mixers by using commercial code, ANSYS CFX. In order to quantify mixing, we employed the concept of mixing index by volume integration.

A wind turbine is one of the most popular energy conversion systems to generate electricity from the natural renewable energy source and an axial-flow type wind turbine is the most popular system for the electricity generation in the wind farm nowadays. In this study, a cross-flow type turbine has been studied for the application of wind turbine for electricity generation. The target capacity of electric power generation of the model wind turbine developing on the project is 12 volts, 130A/H (about 1.56kW). The important design parameters of the model turbine impeller are the inlet and exit angle of the turbine blade, number of blade, hub/tip ratio and the exit flow angle of the casing. In this study, the radial equilibrium theorem was used to decide the inlet and exit angle of the impller blade and CFD technique was used to have the performance analysis of the designed model power turbine to find out the optimum geometry of the CPT impeller and casing. The designed CPT with 24 impeller blades at ${\alpha}=82^{\circ}$, ${\beta}=40^{\circ}$ of turbine blade angle was estimated to generate 284.6 N.m of indicated torque and 2.14kW of indicated power.

The technology of AC solenoid valves is now considered as a core technology in the fields of the production line of semi-conductor chips and the micro fluid chips for medical applications. And AC solenoid valves, which operate by compressed air, are characterized by high speed response, great repeatability and that the pressure on the cross sectional area of poppet is kept constant regardless of the fluctuation of the pressure exerted on the ports. In this study, AC solenoid valves that posses the high-speed responsibility and the high rate of flow have designed and analyzed through the law of equivalent magnetic circuit and Finite Element Method (FEM) respectively. In case of poppet, Flow field characteristic was analyzed by the variation of poppet and it was able to display flow field by changing the location of the poppet. Also, we verified possibility of the design through the static and dynamic pressure and the 3D distribution curve of the force by working the front poppet.

HVOF thermal spray guns are now being widely used to produce protective coatings, on the surfaces of engineering components. HVOF technology employs a combustion process to heat the gas flow and melt the coating materials which are particles of metals, alloys or cermets. Particle flow which is accelerated to high velocities and combustion gas stream are deposited on a substrate. In order to obtain good quality coatings, the analysis of torch design must be performed. The reason is that the design parameters of torch influence gas dynamic behaviors. In this study, numerical analysis is performed to predict the gas dynamic behaviors in a HVOF thermal spray gun with various torch shapes. The CFD model is used to deduce the effect of changes in nozzle geometry on gas dynamics. Using a commercial code, FLUENT which uses Finite Volume Method and SIMPLE algorithm, governing equations have been solved for the pressure, velocity and temperature distributions in the HVOF thermal spray torch.

Aerodynamic forces and moments have been used to control rocket propelled vehicles. If control is required at very low speed, Those systems only provide a limited capability because aerodynamic control force is proportional to the air density and low dynamic pressure. But thrust vector control(TVC) can overcome the disadvantages. TVC is the method which generates the side force and roll moment by controlling exhausted gas directly in a rocket nozzle. TVC is classified by mechanical and fluid dynamic methods. Mechanical methods can change the flow direction by several objects installed in a rocket nozzle exhaust such as tapered ramp tabs and jet vane. Fluid dynamic methods control the flight direction with the injection of secondary gaseous flows into the rocket nozzle. The tapered ramp tabs of mechanical methods are used in this paper. They installed at the rear in the rocket nozzle could be freely moved along axial and radial direction on the mounting ring to provide the mass flow rate which is injected from the rocket nozzle. In this paper, the conceptual design and the performance study on the tapered ramp tabs of the thurst vector control has been carried out using the supersonic cold flow system and shadow graph. Numerical simulation was also performed to study flow characteristics and interactions between ramp tabs. This paper provides to analyze the location of normal shock wave and distribution of surface pressure on the region enclosed by the tapered ramp tabs.

The conventional control valves have been used at the locations occurring high differential pressure and high temperature which causes cavitation, flashing, severe vibration due to abrupt flow change, and sudden pressure drop. Previous studies concerning control valves focused to prevent damage of valve trim due to the internal leak and low flow rate. The newly designed helical trim of control valve has been installed at the location of high pressure change and high temperature in a power plant, and operated for evaluation. It is confirmed that the new control valve developed in this study generates flow characteristics in comparison with previous helical trim of control valves.

Shape optimization of the 3-dimensional WIG airfoil with 3.0-aspect ratio has been performed by using the multi-objective genetic algorithm. The WIG ship effectively floating above the surface by the ram effect and the virtual additional aspect ratio by a ground is one of next-generation and cost-effective transportations. Unlike the airplane flying out of the ground effect, a WIG ship has possibility to capsize because of unsatisfying the static stability. The WIG ship should satisfy aerodynamic properties as well as a static stability. They tend to strong contradict and it is difficult to satisfy aerodynamic properties and static stability simultaneously. It is inevitable that lift force has to scarify to obtain a static stability. Multi-objective optimization technique that the individual objectives are considered separately instead of weighting can overcome the conflict. Due to handling individual objectives, the optimum cannot be unique but a set of nondominated potential solutions: pareto optimum. There are three objectives; lift coefficient, lift-to-drag ratio and static stability. After a few evolutions, the non-dominated pareto individuals can be obtained. Pareto sets are all the set of possible and excellent solution across the design space. At any selections of the pareto set, these are no better solutions in all design space

When a molten corium is relocated in a lower head of a reactor vessel, the ERVC (External Reactor Vessel Cooling) system is actuated as coolant is supplied into a reactor cavity to remove a decay heat from the molten corium during a severe accident. To achieve this severe accident mitigation strategy, the two-phase natural circulation flow in the annular gap between the external reactor vessel and the insulation should be formed sufficiently by designing the coolant inlet/outlet area and gap size adequately on the insulation device. For this reason, one-dimensional natural circulation flow tests were conducted to estimate the natural circulation flow under the ERVC condition of APR1400. The experimental facility is one-dimensional and scaled-down as the half height and 1/238 rectangular channel area of the APR1400 reactor vessel. As the water inlet area increased, the natural circulation mass flow rate asymptotically increased, that is, it converged at a specific value. And the circulation mass flow rate also increased as the outlet area, injected air flow rate, and outlet height increased. But the circulation mass flow rate was not changed along with the external water level variation if the water level was higher than the outlet height.

Although noise of a RAC can be reduced effectively by decreasing RPM, condensation problems can occur to reduce reliability of the RAC for low RPM. Thus, this research has been performed to propose a design guideline of the RAC for low-noise RPM with high reliability. The internal and external flows of the RAC have been visualized and analyzed by a PIV technique to solve the condensation problem at an outlet and impeller. Then, the design guideline has been proposed by the analyzed results and confirmed by wind-tunnel and noise tests to reduce the condensation problem. Finally the shapes of the outlet with reduced condensation problem and the impeller with low noise have been obtained in this study.

In this paper, natural convection of a magnetic fluids(W-40) in a half circular pipe enclosure are investigated by numerical and experimental method. One side wall is kept at a constant temperature(25$^{\circ}C$), and the opposite side wall is also kept at a constant temperature(20$^{\circ}C$). Under above conditions, various magnitudes of the magnetic fields were applied up. Theoretical study through the governing equation derived by Siliomis is carried out with numerical analysis by the GSMAC Method. And the thermo-sensitive liquid crystal film(R20C5A) is utilized in order to visualize wall-temperature distributions as an experimental method. This study has resulted in the following fact that the natural convection of a magnetic fluids are controlled by the direction and intensity of the magnetic fields.

The present study has been carried out to investigate the pumping characteristics for different figures of inlet and outlet in diffuser/nozzle based on piezoelectric micropump. Piezoelectric micropump system consists of several parts like a pumping chamber, diffuser/nozzle, piezoelectric element and tubes. Parts of the micropump connected with diffuser/nozzle and tubes have been analyzed.. The magnified parts have been classified into two different models based on their resistance. These models have been further classified into six models with each one having three different angles at the magnified parts. Each model has been compared and analyzed using the simulation tool, namely, CFD-ACE depending on their flow rates and characteristics.

The effect of Lorentz force(Electromagnetic force) on the liquid metal flow has been investigated. The flow velocity has been calculated by treating the Lorentz force as a source term in the Navier-Stokes equation. The liquid metal flow in the rectangular duct of an electromagnetic pump was analyzed with the Lorentz force varied.

In this paper, the inverse shape design is introduced using the permeable wall boundary condition. Inverse shape design defines the blade shape for the prescribed Mach numbers or pressure distribution on its surface. It calculates the normal mass flux from the difference between the calculated and prescribed pressure at the surface. A new geometry can be achieved after applying the quasi one-dimensional continuity equation from the leading edge to the trailing edge. For validation of this method, two test cases are studied. The first test case of inverse shape design illustrates the cosine bump with a strong shock. After seven geometry modifications, the shock-free bump geometry can be obtained. The second example concerns the redesign of a transonic turbine cascade. The initial isentropic Mach distribution has a peak on the upper surface. The target isentropic Mach number distribution was imposed smoothly. The peak of Mach distribution has disappeared at the final geometry. This proposed inverse design method has proven to be an efficient and robust tool in turbomachinery design fields.

Water resource can be examined using biological sensors. Algae has been one of the biological sensors used to evaluate and monitor the water pollution. The monitoring system, however, could determine whether the examined water was safe or not. It needs additional expensive chemical test to figure out the cause of the water pollution. In this study, an endeavor is given to identify the toxicant in the water using the shape of the chlorophyll fluorescence induction curve(FIC) from algae using monitoring system. Fundamental curves are obtained from the experiments with specified amount of toxicant. Baysian method is utilized to determine the unknown toxicant in the water by comparing it with the fundamental curves. The results shows that the proposed method works fairly well.

The mixing method of water and chemicals is significant in a small-scale water supply system because drinking water should be supplied with a certain quantity of remaining chemicals maintained. In the present study, the concentration distribution and the mixing index were obtained from four models, which were to find out the optimal mixing method of water and chemicals. The two models brought the good mixing effects out of the four, one for providing chemicals from the center of water supply pipe and the other for setting up the semicircle block at the downstream of the chemicals-providing pipe. As a result, the mixing effect was found out to be increased due to the diffusion and the disturbance of flows. In conclusion, these results are expected to contribute to designing the optimal mixing system.

The objective of this paper is to elucidate the characteristic performance of the mini-vacuum pump with various heights of eccentric shaft. The mini-diaphragm pump is composed of DC motor, eccentric shaft and diaphragm and it is operated by PCB panel. The height of eccentric shaft is changed that controls the quantities of air flow to improve the performance of vacuum pressure. This device is manufactured in order to embody the vacuum pressure with 200 mmHg. The heights of eccentric shaft which is used in present experiment are located in the range from 3 mm to 5.5 mm. The vacuum pressure distributions with each height of eccentric shaft was measured and the results were graphically depicted.

The objective of this paper was to investigate the mixing characteristics of an three-stage inline coagulant mixing system experimentally. Wastewater samples of pH 8.5 and initial turbidity 1,000NTU were taken from a site of tunneling work. At the constant dosage, 0.36mL/L, of polymer as coagulant aids, the coagulation efficiency with the dosage of PAC as coagulant was about 4${\sim}$6% at 10 minutes after sampling. In the case of 2 different velocity gradient conditions, the efficiency of turbidity removal was increased about 6.5${\sim}$8% with increasing the dosage of coagulant while, the efficiency was increased about 20${\sim}$21.5% with increasing the dosage of coagulant aids. The efficiency of turbidity removal with the settling time after sampling was about 90% after 1 minute, and the efficiency was about 95% after 5 minutes.

A numerical analysis for temperature distribution of four different materials such as iron, silicon, aluminum and PVC has been performed in this study to predict thermal behaviors of combined weapon systems in a large environmental tester. Thus, experimental conditions have been proposed using a calculating software (SolidWorks 2007 COSMOS FloWorks) to prepare for field tests and analyze heat flow inside the environmental tester and temperature distributions of materials. The boundary conditions of the analysis are composed of inlet and outlet conditions of the environmental tester with different pressures and the limit of low temperature of -30$^{\circ}C$. The soaking time of the system in the environmental tester has been calculated by this commercial program in this study to carry out the experiment.

This paper describes an experimental work to investigate the effect of mesh screen device on the jet structure and acoustic characteristics of over-expanded supersonic jet. The mesh screen device is placed into the supersonic jet stream. In order to perturb mainly the initial jet shear layer, the hole is perforated in the central part of the mesh screen. The diameter of the perforated hole and the location of mesh screen device are varied. A Schlieren optical system is used to visualize the flow fields of supersonic jet without and with the mesh screen device. Pitot pressure measurement is carried out to obtain the pressure distribution in the jet flow. Acoustic measurement also is performed to obtain the OASPL and noise spectra. The results obtained show that the jet structure and the jet noise control effectiveness is strongly dependent upon the diameter of the perforated hole and the location of the mesh screen device in the jet stream. Provided that the mesh screen device is placed at the location to perturb effectively the initial shear layer, the present control method is effective in suppressing the supersonic jet noise.

In the field of yarn dyeing, the most generally employed method is a type of package dyeing which uses a package of cheeses stacked on a spindle made of a perforated tube. In order to understand the process of level dyeing, it is essential to perform a study of the porous flow through the spindle for the cheese dyeing method. In this paper, the axisymmetric, incompressible, Navier-Stokes equations are solved for several spindle configurations using a fully implicit finite volume scheme. For investigating the flow patterns through the spindle, porous diameter and porosity is varied in the present study. The computational results show that the total pressure loss depends only on the velocity of inflow regardless of porous diameter and porosity and a large percentage of the mass flow rate through the spindle is discharged at the upside of the spindle. Therefore, it is required to design a new spindle to obtain the level dyeing.

The OPS (Optical Particle Sensor) using light scattering from the particles (real-time measurement without physical contact to the particles) can be used for cleanroom or atmospheric environment monitoring. For particles smaller than 300 nm, the detection efficiency becomes lower as scattered light decreases with particle size. To obtain higher detection efficiency with small particles, the flow field in particle chamber and the detection volume should be designed optimally to achieve maximum scattered light from the particles. In this study, a commercial computational fluid dynamics software FLUENT was used to simulate the gas flow field and particle trajectories with various optical chamber designs for 300 nm PSL particle. For estimation of laser viewing volume, we used a commercial computational optical design program ZEMAX. The results will be a great help in the development of OPS which can measure small particles with higher detection efficiency.

As the minimum feature size decreases, control of contamination by nanoparticles is getting more attention in semiconductor process. Cleaning technology which removes nanoparticles is essential to increase yield. A reference wafer on which particles with known size and number are deposited is needed to evaluate the cleaning process. We simulated particle trajectories in the chamber by using FLUENT and designed a particle deposition system which consists of scanning mobility particle sizer (SMPS) and deposition chamber. Charged monodisperse particles are generated using SMPS and deposited on the wafer by electrostatic force. The experimental results agreed with the simulation results well in terms of particle number and deposition area according to particle size, flow rate and deposition voltage.

Recently, we have many problems on the process of the sludge. In past, the sewage sludge was treated by reclaimed land or thrown away in the sea. But these methods caused environmental pollution. Today, many researchers are studying various methods for reducing its volume. One of these method, this study is to reduce the moisture of sewage sludge and to solidify it using a dryer and curing equipment. In this research, we investigated about design parameter and operation condition of the equipment. The curing equipment reduces the percentage of water content from 30% of dryer to 10%. So, we have to study the curing characteristics and performance of curing equipment. For example, there are internal flow characteristics and change of the percentage of water content. And we investigated the change of data at outlet along the initial condition, temperature, humidity and air flow. Using this data, we achieve the experimental results of curing efficiency by each geometry and operating condition. And we also investigated numerical analysis of internal flow using CFD code. This research is basic study for optimal design of the curing equipment.

The interferometer method with nano-scale spatial resolution has been developed in this study. To enhance the accuracy of the previous developed method, the 14 bit cooled CCD camera with 1280 by 980 spatial resolution was applied to the measurement. And optical alignment has been carried out on the highly accurate position sensors with 500nm resolution so as to be able to calibrate the detected interference image with the field of view. Also the measurements were applied to the ultra thin oil film between the Al coated cylinder mirror with 38.1mm radius and 0.5mm cover glass to verify the developed method. The measured result showed the good agreement with the used cylinder curvature with ${\pm}$5.18run uncertainty.

Taylor-Couette flow may appear when the angular velocity is different between two concentric rotating cylinders. This kind of Taylor-vortex flow can be easily seen in lots of engineering problems. In general the geometries of rotating cylinders are generally complex in these cases. In this study, we investigated Taylor-Couette flow when the outer cylinder has the slit along the annulus. The radius ratio and aspect ratio of the experimental model used was 0.825 and 48, respectively. The depth of slits is 5mm and total 18 slits are azimuthally located along the inner wall of outer cylinder. We used PIV method to measure the flow and applied index matching method to resolve the complex geometry effect. The results show the model with slit has no stable wavy vortex region above Re=143.

Elastin-Like polypeptide (ELP) composed of elastin-based repeating units is an artificial biomaterial which is biocompatible and non-immunogenic. ELP shows a characteristic inverse phase transition between hydrophobic and hydrophilic phase by external stimuli such as salt, pH and temperature. In this study, ELP coated PS (polystyrene) beads are packed in tubing and the thermo -responsive flow characteristics of the packed bed are investigated. Preliminary test results show that the control of the fluid flow can be achieved by using the temperature driven phase transition effect of the ELP coated beads in a microchannel.

The present study addresses a method to operate a fuel-cell system effectively using a recirculation ejector which recycles wasted hydrogen gas. Configuration of a recirculation ejector is changed to investigate the flow behavior through it under varying operating conditions, and how such conditions affect the fuel-cell hydrogen cycle. The numerical simulations are based on a fully implicit finite volume scheme of the axisymmetric, compressible, Reynolds-Averaged, Navier-Stokes equations for hydrogen gas, and are compared with available experimental data for validation. The results show that a hydrogen recirculation ratio is effectively controlled by a configurational alteration within the operational region in which the recirculation passage doesn't plugged by a sonic line.

An analysis program for pedestrian flow has been developed to investigate the flow patterns of passenger in railway stations. Analysis algorithms for pedestrian flow based on DEM(Discrete Element Method) are newly developed. There are lots of similarity between particle-laden two phase flow and passenger flow. The velocity component of 1st phase corresponds to the unit vector of calculation cell, each particle to passenger, volume fraction to population density and the particle velocity to the walking velocity, etc. And, the walking velocity of passenger is also represented by the function of population density. Key algorithms are developed to determine the position of passenger, population density and numbering to each passenger. By using the developed program, we compared the simulation results of the effects of the location and size of exit and elapsed time.

The performance study shows the result using two different methods which are used to control missile or aircraft. One is the Thrust Vector Control(TVC) method for the aviation of next generation and the other is the present effective Shroud Jet-vane System(SJVS) method for the satellite effector development. The research was done through the performance estimation using the numerical simulation analysis, the modelling, the performance measuring using the model, the investigation of the flow visualization and aerodynamic performance with the enforced power to the vane and the result comparison.

This study numerically investigates the unsteady flow and acoustic characteristics of a flapping wing using a hydrodynamic/acoustic splitting method. The Reynolds number based on the maximum translation velocity of the wing is Re=8800 and Mach number is M=0.0485. The flow around the flapping wing is predicted by solving the two-dimensional incompressible Navier-Stokes equations (INS) and the acoustic field is calculated by the linearized perturbed compressible equations (LPCE), both solved in moving coordinates. Numerical results show that the hovering sound is largely generated by wing translation (transverse and tangential), which have different dipole sources with different mechanisms. As a distinctive feature of the flapping sound, it is also shown that the dominant frequency varies around the wing.

Gas-liquid separator has been designed for the sake of reducing expenses associated with production operations. To date, a number of gas-liquid separators have been installed and put to use for various applications. Despite the advantages of simple and compact configuration of separator with no moving part, its efficient operation is limited in terms of total pressure losses, separation performance and flow-induced noise and vibration, which are closely associated with the very complicated flow phenomena involved. In the present study, a gas-liquid centrifugal separator with a swirl vane is investigated for the purpose of water separation from compressed moisture air. The 3D Navier-Stokes equations are numerically solved using a fully implicit finite volume scheme. Based upon the obtained solutions, tangential velocities, centrifugal forces, vortices and total pressure losses are analyzed to find out the best design parameters. From the present study, several attempts are made to improve the performance of conventional separators of centrifugal type.

Recently, nuclear power plant companies have been extending the turbine valve test interval to reduce the potential of the reactor trip accompanied with a turbine valve test and to improve the NPP's economy through the reduction of unexpected plant trip or decreased operation. In these regards, the extension of the test interval for turbine valves was reviewed in detail. The effect on the destructive overspeed probability due to the test interval change of turbine valves is evaluated by Fault Tree Analysis(FTA) method. Even though the test interval of turbine valves is changed from 1 month to 3 months, the analysis result shows that the reliability of turbine over speed protection system meets acceptance criteria of 1.0E-4/yr. This result will be used as the technical basis on the extension of the test interval for turbine valves. In this paper, the propriety of the turbine valve test interval extension is explained through the review on the turbine valve test interval status of turbine overspeed protection system, the analysis on the annual turbine missile frequency and the probability evaluation of the destructive overspeed due to the test interval extension.

This study aims to analyse the influence of steam injection on the performance of hybrid systems combining a solid oxide fuel cell and a gas turbine. The steam is generated by recovering heat from the exhaust gas. Two system configurations, with difference being the operating pressure of the SOFC, are examined and effects of steam injection on performances of the two systems are compared. Two representative gas turbine pressure ratios are simulated and a wide range of both the fuel cell temperature and the turbine inlet temperature is examined. Without steam injection, the pressurized system generally exhibits better system efficiency than the ambient pressure system. Steam injection increases system power capacity for all design cases. However, its effect on system efficiency varies much depending on design conditions. The pressurized system hardly takes advantage of the steam injection in terms of the system efficiency. On the other hand, steam injection contributes to the efficiency improvement of the ambient pressure system in some design conditions. A higher pressure ratio provides a better chance of efficiency increase due to steam injection.

Design performances of various configurations of the PEMFC/GT hybrid systems have been evaluated. Based on PEMFC adopting steam reforming, various system configurations (one ambient pressure configuration and three different pressurized configurations) were designed and their performances were compared. Their Performances are also compared with the reference PEMFC system. Influences of turbine inlet temperature, pressure ratio on the hybrid systems performance were investigated and design ranges exhibits better efficiency than the PEMFC system were presented. One of the pressurized system may have much higher efficiency than the PEMFC system, while other systems hardly provide efficiency upgrade.

This paper describes the results of steam turbine performance tests. The objectives of performance test is to exactly evaluate the degradation(decrease in performance) of the coal-fired steam turbine generator in order to provide plant information to help performance engineers identify problems, improve performance, and make economic decisions about scheduling maintenance and optimizing operation. To achieve these goals, the periodic thermal performance tests have been carried out since the initial operation period, 1997. We made the calculation program and guidelines for the tests and developed the performance index of the turbine cycle on the basis of the ASME PTC. By comparing the performance changes throughout the whole operation period, we confirmed the performance reliabilities of the turbine and its conditions.

Integrated Gasification Combined Cycle (IGCC) power plant converts coal to syngas, which is mainly composed with hydrogen and carbon monoxide, by the gasification process and produces electric power by the gas and steam turbine combined cycle power plant. The purpose of this study is to investigate the influence of the syngas to the performance of a gas turbine in a combined cycle power plant. For this purpose, a commercial gas turbine is selected and its performance characteristics are analyzed with three different fuels, i.e., natural gas ($CH_4$), syngas and hydrogen. It is found that different heating values of those fuels and chemical compositions in their combustion gases are the causes in the different performance characteristics.

Endwall losses contribute significantly to the overall losses in modern turbomachinery, especially when aerodynamic airfoil load and pressure ratio are increased. Hence, reducing the extend and intensity of the secondary flow structures helps to enhance overall efficiency. From the large range of viable approaches, a promising combination positioning and height of endwall contouring was chosen. The objective of this study is to document the three-dimensional flow in a turbine cascade in terms of streamwise vorticity, total pressure loss distribution and static pressure distribution on the endwall and blade surface and to propose an appropriate positioning and height of the endwall contouring which show best secondary, overall loss reduction among the simulated endwall. The flow through the gas turbine were numerically analyzed using three dimensional Navier-Stroke equations with a commercial CFD code ANSYS CFX-10. The result shows that the overall loss is reduced near the flat endwall rather than contoured endwall, and the case of contoured endwall installed at 30% from leading edge with height of 25% for span showed best performance.

A large eddy simulation (LES) has been conducted for the flow and heat transfer in a dimpled channel. Two dimple depths of 0.2 and 0.3 times of the dimple print diameter (= D) have been compared at the bulk Reynolds number of 20,000. Three Reynolds numbers of 5,000, 10,000 and 20,000 have been studied, while the dimple depth is kept as 0.2 D. With the deeper dimple, the flow reattachment occurs father downstream inside the dimple, so that the heat transfer is not as effectively enhanced as the case with shallow ones. At the low Reynolds number of 5,000, the Nusselt number ratio is as high as those for the higher Reynolds number, although the value of heat transfer coefficient decreases because of the weak shear layer vortices.

This paper describes the proceedings of creating countermeasures after analysis and maintenance be able to conduct operation safely in a power plant. in order to operate the power plant in a stable and reliable way, the best condition of the govemor system can be maintained through the response characteristic analysis of the control device for the pitch blade control hydraulic actuator. The fan pitch blade control hydraulic actuator of a 500MW large-scale boiler is frequently operated under normal operation conditions. Common problems or malfunctions of the pitch blade control hydraulic actuators leads to the decline of boiler thermal efficiency and unexpected power plant trip. The inlet and outlet gas can be controlled by using the fan pitch blade control hydraulic actuator in order to regulate the internal pressure of the furnace and control the frequency in the power plant facility which utilizes soft coals as a power source.

Combustion instability is a serious obstacle for the lean premixed combustion of gas turbines, and can even cause fatal damage to the combustor and the entire system. Thus, improved understanding of the mechanisms of combustion instability is necessary for designing and operating gas turbine combustors. In this study, in order to understand the instability phenomena, an experimental study was conducted in a rearwardstep dump combustor with LPG and air. The fluctuations of pressure and heat release were measured by piezoelectric pressure sensor and High speed Intensified Charge Coupled Device (ICCD) camera respectively. Various types of combustion modes occurred in accordance with the equivalence ratio and the fuel supplying conditions. The unmixedness of the fuel and air can be controlled by changing the mixing distance ($L_{fuel}$). It is found that the unmixedness of the fuel and air affects the characteristics of flame behavior and pressure fluctuations in a lean premixed flame.

DME (Dimethyl Ether, $CH_3OCH_3$) has highly attracted attention as an alternative fuel for transportation, power generation and LPG substitute owing to its easy transportation and cleanliness. This study was conducted to verify the combustion performance and to identify potential problems when DME is fuelled to a gas turbine. GE7EA gas turbine of Pyong-Tak power plant was selected as a target to apply the DME. Combustion tests were conducted by comparing DME with methane, which is a major component of natural gas, in terms of combustion instability, $NO_X$ and CO emissions, and the outlet temperature of the combustion chamber. The results of the performance tests show that DME is very clean but has a low combustion efficiency in low load condition. From the results of the fuel nozzle temperature we have ascertained that DME is easy to flash back, and this property should be considered when operating a gas turbine and retrofitting a burner.

A novel oxy-fuel burner for a boiler has been devised and composed into a 50 kW class boiler system. A series of test has been conducted to show the characteristics of combustion, exhaust gas and the boiler. Numerical simulations have been also performed and validated against the experimental data to discuss detailed physics. The oxy-fuel burner can effectively heat the combustion chamber with the significantly reduced combustion gas, which enables to realize the compactness of the system. The composition of exhaust gas reveals that the sealing of the system is crucial to achieve high $CO_2$ concentration and low $NO_X$ emission.

Lean combustion is one of the most promising method for increasing engine efficiency and reducing the exhaust emission from SI gas engines. Due to the possibility of partial burn and misfire, however, under lean burn operation, stable flame kernel formation and fast burn rate are needed to guarantee a successful subsequent combustion. Experiment data were obtained on a single-cylinder CNG fueled SI engine to investigate the effect of direct injection, spark timing and variation of injection timing. Experimental results show that lean burn limit is ${\lambda}$=1.3 with port injection, and expansion of lean burn limit ${\lambda}$=1.4 with direct injection method, due to increase of turbulence intensity in cylinder and stratified charge. Combustion duration in lean region is improved by using the variation of injection timing.

Compressed natural gas has good potential for alternative vehicle fuel due to its economical and clean characteristics. However, the composition of natural gas based on production location is known to affect performance and emissions of CNG engine. Thus, the objective of this paper is to clarify the effect of fuel composition on combustion and emissions of CNG engine. This paper presents combustion characteristics obtained from running a 2.5L, 4-cylinder CNG engine retrofitted IDI diesel engine with engine dynamometer. BSFC, emissions, fuel consumption and combustion pressure were measured under steady state operating conditions especially at partial load for CNG engine. Based on the experimental results, we found that CNG composition affects engine performance, fuel conversion efficiency and burning rate.

These days the exhaustion of petroleum resources and environmental problems are getting serious. Many researchers are focused on low emission and high performance vehicles. Therefore, we should concern about emission regulation when we design a new car. In this study, we investigated the characteristics of the traditional mechanical engine cooling systems which control the engine temperature using engine speed and wax type thermostat. This experiment used three components which are Radiator fan, water pump and water valve controlled by an electronic system based on the engine status (load, speed). We elucidated how different between traditional mechanical cooling system and electronic cooling system which control coolant temperature and coolant flow rate in a DI diesel engine in this paper. The results revealed a fuel saving and an emission (CO, HC) reduction on NEDC cycle.

Ignition delay of second injection of HSDI diesel engine was usually much shorter than that of first injection. It is due to the interaction between radicals generated during the combustion process, and mixed gas of second injection. In this paper, To analyze combustion phenomena of multiple injection mode in HSDI diesel engine effectively, two-dimensional flamelet combustion model was modified. To reduce calculation time, two-dimensional flamelet equations were only applied near stoichiometric region. If this region was ignited, species and temperature of other region were changed to the steady-state solutions of one dimensional flamelet equations. By this method calculation time for solving flamelet equations was reduced to 20 percents, thought the results were almost same. Modified flamelet combustion model was coupled to commercial CFD code interactively using user subroutine.

Many researches have been studied on in-cylinder flow as one of dominant effects for an engine combustion. Specially because the combustion flame speed is mainly determined by the turbulence at the end of compression process. Tumble and Turbulence ahead of combustion is very important phenomenon. As this phenomenon make research certainly, combustion condition will effectively be improve. This paper describes analytical results of the tumble flow, intensity, turbulence inside the cylinder of maritime engine. 3-D computation has been performed by using STAR-CD v3.26 solver and es-ice

To analyze the influence of valve overlap period on a backfire occurrence, the single cylinder research engine with MCVVT(Mechanical Continuous Variable Valve Timing) system is developed and backfire limit equivalence ratio defined as fuel-air ratio equivalence ratio at which backfire occurs is examined according to various valve overlap period. The MCVVT is the system to control valve overlap period by mechanical device. It is estimated that the lower valve overlap period has the higher backfire limit equivalence ratio though the same energy is supplied. When the valve overlap period is changed from 30$^{circ}$ CA to 0$^{circ}$ CA, backfire limit equivalence ratio is increased 74%, approximately. It means that valve overlap period is concern in backfire occurrence, and may be one of the methods for controlling back fire occurred in a $H_2$ engine.

THC(Total Hydrocarbon) emissions during cold start and warm-up phase constitute the majority of THC emissions during the FTP-75 mode. As the basic approach to improve the emission performance of Gasoline engine during transient phase, the effect of spark timing retard from MBT on THC emission characteristics is studied by engine test using a Fast response Flame Ionization Detector(FFID). A cyclic analysis of the combustion process shows that high THC emissions are produced first few cycles during the transient phase. This paper presents the results of engine performance and emission of Gasoline engine with various spark timing. consequently, This paper was focused on the combustion phenomena with various spark timing during transient phase which was analyzed by Fast response Flame Ionization Detector (FFID) equipment to measure the cyclic THC emission characteristics.

Conventional LPG pump for Liquified Petroleum injection(LPi) engine has been adopted vane type. But the BLDC type fuel pump for LPi system has complicated structure and its price is high. Therefore, as a alternative, this study has mainly focused on the development of turbine type LPG pump which has lower cost and simple structure than conventional BLDC type. To verify the possibility of substitute the performance tests were performed for each fuel pump. The comparative items were pressure settling time, variation of fuel outlet temperature and engine performance of hot restart ability. As a result, performances of turbine type LPG pump were equivalent or high comparing to the BLDC type all over the tests for different fuel composition.

The purpose of the research is to investigate of diesel spray combustion for simultaneously reduce way NOx and PM. The pressure diesel injection were done into intermediates that are generated by very lean DME HCCI combustion using a RCM. The concentration of intermediate could not be directly measured; we estimated it by CHEMKIN calculation. DME HCCI characteristic is surveyed. Validations of the CHEMKIN calculation were confirmed pressure rise of an experiment and pressure rise of a calculation. Using a framing streak camera captured two dimensional spontaneous luminescence images from chemical species at low temperature reaction(LTR) and high temperature reaction (HTR). Also, the combustion events were observed by high-speed direct photography, the ignition and combustion were analyzed by the combustion chamber pressure profiles.

The government have been tightening EM regulation gradually but the effect is not good because of rapid increase of vehicles. And medium & heavy duty diesel vehicles, even though the number is small, exhaust very large pollutants(about over 50%). Especially it is more severe about old trucks and buses. Accordingly, CNG vehicle and the retrofit of diesel to CNG must be an alternative in order to protect the atmospheric environment and improve the air quality in the metropolitan area. The main object of this study is to secure the retrofit technology of diesel to CNG vehicle and the management system of CNG engine. we passed the government emission certification test. In addition to this, the mass production for retrofit is also studied. Results of emission and durability test for certification are as follows; there was no problem during 30,000km vehicle durability test and good emission levels satisfying the regulation.

This paper presents effects of carburetor throttle angle on the performance characteristics of a small sparkignited, gasoline engine. The engine used in the test is a single cylinder, two-stroke, air-cooled 26cc SI engine for brush-cutter. We measured the rpm, torque, power, and fuel consumption according to the six different throttle angle conditions of the rotary-type carburetor and to the engine dynamometer loads. We had concluded that maximum power happened at 5000${\sim}$6000rpm and at the same condition was the minimum specific fuel consumption.

This paper describes a research for the performance improvement of the straight-bladed vertical axis wind turbine. To improve the performance of VAWT, the individual blade pitch control method is adopted. For the wind turbine, CFD analysis is carried out by changing blade pitch angle according to the change of wind speed and wind direction. By this method, capacity and power efficiency of VAWT are obtained according to the wind speed and rotating of rotor, and could predict the overall performance of VAWT. It was manufactured to verify performance of the experimental system that consists of rotor including four blades and base. Furthermore, torque sensor and power generator were installed. Also, active controller which can change the pitch angle of the individual blade according to the wind speed and direction was used.

The present paper is devoted to investigate dynamic effect and steady-state performance of methane autothermal reformer theoretically and numerically. In order to simplify the complicated phenomena in the system, axisymmetric heterogeneous reactor model is developed. As autothermal reaction takes places on catalyst surface between bulk gas and catalyst, volume averaging method is incorporated using porous medium approach. To understand the start-up process which occurs in the reactor is highly important. Therefore, in this paper we get various goverining equations to find out transient and steady solutions and time scale for start-up introducing dimensionless variables. Start-up is a significant issue in reforming reaction for automobile system and fueling of SOFC-based auxiliary power units. This paper deals with characteristics of heat and mass transfer and predicted light-off time in the reformer as oxygen to carbon ratio ($O_2$/C) and amount of feeding gas.

The smart energy system is the integrated power system in which the power components including central station generation, distributed generation, renewable power generation, energy storage, and communications and controls are complexly connected with each other. In smart energy system, it is very important how to configure the diverse power generations and how to determine the operation mode of the chosen components with economic feasibility. In this study, we introduce the optimal planning method based on both economic feasibility and load profiles and its applications for the smart energy system in apartment. This method was considered very useful to determine the configuration and to decide the optimal operation mode of the smart energy system.

We developed micro power generation system using piezoelectric materials. In our system, the ambient vibrating energy is converting to electric energy by deflection of piezoelectric beams. The system consists of energy generating parts, converting enhancement parts, electric regulation and charging parts, and interface with small-energy-consuming mobile devices. The geometry of piezoelectric beams, the source of vibrating energy, and the electric load of target application determine the characteristics of generating electric power, such as impedance, voltage, current and power density. Therefore, we made a model for analysis of generating power with given information such as piezoelectric materials, geometry, vibration type, and mass. With this model, we can calculate capacitance of piezoelectric beams, generating voltage, current, and power. To obtain maximum energy transfer efficiency, we approached this study in the view of material, electrical, and mechanical engineering

The design of a ground-source heat pump system includes specifications for a ground loop heat exchanger where the heat transfer rate depends on the thermal conductivity of the ground. To evaluate this heat transfer property, in-situ thermal response tests on four vertical test boreholes with different grouting materials were conducted by adding a monitored amount of heat to water over various test lengths. By measuring the water temperatures entering and exiting the loop, water flow rate, and heat load, effective thermal conductivity values of the ground were determined. The effect of increasing thermal conductivity of grouting materials from 0.82 to 1.05 W/m$^{\circ}C$ resulted in overall increases in effective ground thermal conductivity by 25.8% to 69.5%.

In proton exchange membrane fuel cell, the heat is generated at the catalyst layer as result of exothermic electrochemical reaction. This heat increases temperature of gas diffusion layer and membrane whose conductivity is very sensitive to humidity, function of temperature. So it is very important to analysis heat transfer through fuel cell to maintain temperature at specified range. In this paper numerical simulation was done including reversible, irreversible, ionic resistance, water formation loss to source term of energy equation. Results show that irreversible and water formation loss contributes mainly to energy source term and as current density increases, all of energy source terms become increased and Nusselt number is increased as results of more heat generation. Particularly irreversible loss is found to be predominant among the all energy source and water formation at cathode channel influences the temperature distribution of fuel cell greatly.

To compare the sensor performance of AE leak diagnosis system which can measure valve leak conditions, AE activities such as RMS voltage level, AE signal trend, leak rate degree according to AE database, FFT spectrum were measured on valve of the simulated test system for power plant. AE activities were recorded and analyzed from various operating conditions including different temperature, pressure difference, valve size and fluid using both piezoelectric acoustic emission sensor and Pb-Free acoustic emission sensor. The results of this study are utilized to select the type of sensors, the frequency band for filtering and thereby to improve the signal-to-noise ratio for diagnosis or monitoring of valves in operation. As the final result of application study above, portable type leak diagnosis system by AE was developed. The outcome of the study can be definitely applied as a means of the diagnosis or monitoring system for energy saving and prevention of accident for power plant valve.

Characteristics of heat transfer and pressure drops of fills for solar tower volumetric air receivers are experimentally investigated with the material and the thickness. The volumetric air receiver considered in this paper consists of a ceramic tube and fills are inserted in the ceramic tube. Air is used as the working fluid. Two materials, which are a honeycomb(diameter: 100mm, thickness: 30mm) and laminated mesh(diameter: 100mm, thickness: 1mm), are considered as the fills. In order to investigate the characteristics of heat transfer of fills, this volumetric air receiver is heated by an electric heater and air temperatures in ceramic tube are measured. Also, the radiative shields are installed to measure the only air temperature. In addition, the pressure losses are measured with the thickness of fills while the air goes through the fills inserted in an acrylic tube. The flow becomes turbulent and fully developed in front of the fills. The results show that the heat transfer and pressure drop characteristics of the laminated mesh are superior to those of the honeycomb.

Heat transfer experiments at a vertical annulus passage were carried out in the SPHINX(Supercritical Pressure Heat Transfer Investigation for NeXt Generation) to investigate the heat transfer behaviors of supercritical $CO_2$. The collected test data are to be used for the reactor core design of the SCWR (SuperCritical Water-cooled Reactor). The mass flux was in the range of 400${\sim}$1200 kg/$m^2$s and the heat flux was chosen up to 150 kW/$m^2$. The selected pressures were 7.75 and 8.12 MPa. The heat transfer data were analyzed and compared with the previous tube test data. The test results showed that the heat transfer characteristics were similar to those of the tube in case of a normal heat transfer mode and degree of heat transfer deterioration became smaller than that in the tube. Comparison of the experimental heat transfer coefficients with the predicted ones by the existing correlations showed that there was not a distinct difference between the correlations.

This study was performed to acquire the reliable in-situ thermal conductivity of closed type ground heat exchanger used in ground source heat pump. We selected four sites(Cheonan, Daejeon, Daegu, Gwangju) which are central area of South Korea. Test results show that the effective thermal conductivities are 2.33 W/m$^{\circ}C$, 2.50 W/m$^{\circ}C$, 2.75 W/m$^{\circ}C$ and 2.86 W/m$^{\circ}C$. From this data, we can see that thermal conductivity varies about the range of 23% with the sites. Also, thermal conductivity increases up to 20% by changing grouting material from low salica sand to high one.

The performance characteristics of the polymer electrolyte fuel cells (PEFCs) were investigated under various humidification conditions at steady-state and transient conditions. The PEFC studied in this study was characterized by I-V curves in potentiostatic mode. The I-V curves representing steady-state performance were obtained from OCV to 0.25V, and the dynamic performance responses were obtained at some points of voltages. The anodic external humidification was applied and the humidity was controlled from 20% to 100%. The effects of relative humidity of hydrogen were measured with the dry air at the cathode. At high voltage region, the performance at high temperature was higher, but at low voltage region, low temperature condition showed the higher performance. The dynamic responses were observed at the instant when the voltage of the PEFC was changed. It was observed that the performance reached steady-state earlier with the increase of temperature.

In recent semiconductor manufacturing clean rooms, air washers are used to remove airborne gaseous contaminants such as $NH_3$, SOx and organic gases from the outdoor air introduced into clean room. In the present study, an experiment was carried out to examine the improvement of removal efficiency for the gaseous contaminants. In order to improve the gas removal efficiency, a hot water contact heat exchanger was installed upstream of the air washer to heat and humidify the incoming outdoor air before entering the air washer.

The neutron imaging technique was used to investigate the water discharge characteristics at PEMFC. Prior to investigation of water discharge characteristics, the linear attenuation coefficient for water at Neutron Radiography Facility (NRF) was calibrated. The feasibility test apparatus was consisted of pressurized air and water in order to simulate the actual operating PEMFC. The feasibility tests have been performed at 1-parallel serpentine type with 100 $cm^2$ active area and different air flow rate (1, 2, and 4 lpm). The total water volume variations at each condition were calculated from the neutron images. The water at channel is well discharged as soon as supplying the pressurized air into the PEMFC. However, because the water at MEA isn't removed the total water volume is constant after 150. Therefore more effective method is needed in order to discharge water at MEA, and the neutron imaging technique is helpful for it.

This paper shows that inlet humidity condition at cathode side is one of dominant parameters affecting the performance of PEMFC. To investigate effects of inlet humidity condition, the performance measurements were conducted for a single PEMFC with two operating variables : cathode relative humidity and dry condition in anode dry. The fuel cell employed for the experiments is a unit PEMFC with a 25$Cm^2$, Nafion$^(R)$112 membrane. As a result of this study, the cell performance is getting higher by increasing inlet humidity condition at cathode side. The cell performance is different from each operating temperature an it has maximum30% higher than dry condition at 60$^{\circ}C$ operating temperature with 80% relative humidity.

An experimental study was carried out to investigate gasification process of wood sawdust in the I-dimensional downdraft fixed bed gasifier. The preheated air was used oxidizer and steam were used as a gasifying agent. The operating parameters, the supplied air temperature and steam were used. The oxidizer temperature was varied from 500K to 620K and vapor was added. The gasification process was monitored by measuring temperature at three position near the biomass using R-type thermocouples and the syngas composition was analyzed by gas chromatograph. The change of hydrogen and carbon monoxide, carbon dioxide, methane was observed. Overall, the volume fraction of hydrogen and methane were increased widely as increasing the oxidizer temperature and adding steam.

The aim of this study is to find out the condition that generates maximum $H_2$ through the calculation of equilibrium model with conditions of pyrolysis gases of Refuse Plastic Fuel(RPF). This study deals with the computational simulation of a RPF gasification using an equilibrium model based on minimization of the Gibbs free energy. An equilibrium analysis was carried out to determine species composition of Syngas in RPF gasification and reactions to variation of temperature, $O_2$/Fuel ratio and Steam/Fuel ratio. Calculated results showed that $O_2$/Fuel ratio, Steam/Fuel ratio and temperature affected on mole fraction of $H_2$, CO.

For proton exchange membrane fuel cell, it is very important to design the flow channel on separation plate optimally to maximize the current density at same electrochemical reaction surface and reduce the concentration polarization occurred at high current density. In this paper, three dimensional computation model including anode and cathode domain together was developed to examine effects of flow patterns and operation conditions such as humidity and operating temperature on performance of fuel cell. Results show that voltage at counter flow condition is higher than that at coflow condition in parallel and interdigitated flow pattern. And fuel cell with interdigitated flow pattern which has better mass transport by convection flow through gas diffusion layer has higher performance than with parallel flow pattern but its pressure drop is increased such that the trade off between performance and pressure drop should be considered for selection of flow pattern of fuel cell.

The thermal behavior of a building in response to heat input from an active solar space heating system is analysed to determine the effect of the variable storage tank temperature on the cycling rate, on and off temperature of a heating cycle and on the comfort characteristics of room air temperature. A computer simulation of the system behavior has been performed and verified by comparisons with various parameters. Especially, this study is focused on the effect of the system's performance when subjected to dynamic cooling loads. The heat input to the absorption system is provided by an array of solar collectors that coupled to a thermal storage tank.

Experiments were carried out in an atmopheric pressure, lab-scale gas turbine combustor to see the effect of partial premixing on unstable flame structure and $NO_X$ emission characteristics. The swirl angle is 45 deg., fuel-air mixing degrees were varied 0, 50, and 100% respectively at equivalence ration ranging from 0.53 to 0.79. The evaluation of phased-locked OH chemiluminescence images were acquired with an ICCD. $NO_X$ emission characteristics were also investigated at each experimental condition. The effect of the fuel-air mixing degree on the flame structure was obtained from phase-locked $OH^*$ images. And it was obtained from local heat release characteristics that the information about the region which the combustion instability was amplified or damped. It also could be confirmed that ${\sigma}$ has greatly influence on $NO_X $emission characteristics at lean regimes. It would be expected that it could provide invaluable data for understanding the mechanism of combustion instability

Leakage would happen because of the damage of high temperature and high-pressure valve in nuclear power plant. condition based prevention maintenance is essential by using the suitable method based on local condition. Energy loss prevention can prevent from an accurate test, Local actually and ability. The methods of test for high energy fluid leakage at present are analysis of ${\Delta}$T, AE(Acoustic Emission) analysis, and thermal image. The result for test of AC (Main steam) system in YNG 2 Unit reveals that the AE occurred clearly in leakage situation, but thermal image didn't occur. It is identified that leakage is occurred when the orifice located front and back of valve operates. It shows that making a impatient judgment by using the single method if it is leakage is containing uncertainty. So I think that using the Multi-Measuring method is more sound judgment than Single-Measuring method.

Power uprate is the process of increasing the maximum power level at which a commercial nuclear power plant may operate. Power uprate applications(113 units) for NPPs(Nuclear Power Plants) were recently approved in the United States. Utilities have been using power uprates since the 1970s as a way of increasing the power output of their nuclear plants. To increase the power output of a reactor, typically more highly enriched uranium fuel and/or more fresh fuel is used. This enables the reactor to produce more thermal energy and therefore more steam, driving a turbine generator to produce electricity. In this paper, the propriety of power uprate is explained through the review on the power uprate method and the changes of the physical parameters due to power uprate. The analysis results showed that the CDF(Core Damage Frequency) and LERF(Large Early Release Frequency) are affected in the current probabilistic safety assessment (PSA) model.

The effects of the wall geometry on the spray-wall impingement process of a hollow-cone fuel spray emerging from a high-pressure swirl injector of the Gasoline Direct Injection (GDI) engine were investigated by means of a numerical method. The ized Instability Sheet Atomization (LISA) & Aerodynamically Progressed Taylor Analogy Breakup (APTAB) model for spray atomization process and the Gosman model were applied to model the atomization and wall impingement process of the spray. The calculation results of spray characteristics, such as a spray development process and a radial distance after wall impingement, compared with the experimental ones by the Laser Induced Exciplex Fluorescence (LIEF) technique. It was found that the radial distance of the cavity angle of 90$^{circ]$ after wall impingement was the shortest and the ring shaped vortex was generated near the wall after spray-wall impingement process.

A dual side cooling annular fuel having internal and external coolant channels has many advantages basically due to low fuel temperature and high DNBR margin, which can make a significant increase of core power density possible. So recently a 12x12 square annular fuel array was proposed for the fuel assembly to be reloaded without structural interference with operating reactors of OPR-1000s. Even through the inherent potential of the annular fuel on the high power density, it may be seriously eroded in the case of a severe unbalanced mass flux split to the internal and external channels in standpoint of DNB. Mass flux split is determined pressure drop characteristics between inner and outer channels. The spacer grids binding fuel array influence greatly the pressure drop in outer channels and the mass flux split. As an important factor of DNB behavior, the enthalpy differences at both channel exits were evaluated using the mass flux splits.

Two steam explosion experiments were performed in the TROI facility by using metal-added molten corium (core material) which is produced during a postulated severe accident in the nuclear reactor. A triggered steam explosion occurred in a case, but no triggered steam explosion did in the other case. The dynamic pressure and the dynamic load measured in the former experiment show a stronger explosion that those performed previously with oxidic corium. A steam explosion is prohibited when the melt temperature is low, because the melt is easily solidified to prevent a liquid-liquid interaction.

The TEXAS-V code tuned for TROI-13 was used for analyzing the parametric findings in TROI experiments. The calculations on the melt composition are relatively similar to the TROI experimental results. The water depth effect in TEXAS-V code seems to be consistent with TROI experiments in some degree. The water area effect of TEXAS-V calculations seems not to be harmonious to that in TROI experiments. This seems to indicate that TEXAS-V as 1-dimensional code or as the numerical steam explosion has a limitation on estimating area effect. Thus, TEXAS-V tuned for TROI-13 seems to have an ability to estimate the parametric effect of TROI experiments. The evaluated TEXAS-V was used for estimating the ex-vessel steam explosion load. The calculated explosion pressure and load were about 40 MPa and 75 kPa.sec, which are not much threatening level for containment integrity, but are arguable value for the integrity.

If a molten core is released from a reactor vessel into a reactor cavity during a severe accident, an important safety issue of coolability of the molten core from top-flooding and concrete ablation due to a molten core concrete interaction (MCCI) is still unresolved. The released molten core debris would attack the concrete wall and basemat of the reactor cavity, which will lead to inevitable concrete decompositions and possible radiological releases. In a OECD/MCCI project scheduled for 4 years from 2002. 1 to 2005. 12, a series of tests were performed to secure the data for cooling the molten core spread out at the reactor cavity and for the 2-D long-term core concrete interaction (CCI). The tests included not only separate effect tests such as a melt eruption, water ingression, and crust failure tests with a prototypic material but also 2-D CCI tests with a prototypic material under dry and flooded cavity conditions. The paper deals with the transient simulations on the CCI-2 test by using a severe accident analysis code, CORQUENCH, which was developed at Argonne National Laboratory (ANL). Similar simulations had been already per for me d by using MELCOR 1.8.5 code. Unlike the MELCOR 1.8.5, the CORQUENCH includes a melt eruption mode I and a newly developed water ingression model based on the water ingression tests under the OECD/MCCI project. In order to adjust the geometrical differences between the CCI-2 test (rectangular geometry) and the simulations (cylindrical geometry), the same scaling methodology as used in the MELCOR simulation was applied. For the direct comparison of the simulation results, the same inputs for the MELCOR simulation were used. The simulation results were compared with the previous results by using MELCOR 1.8.5.

The internal structure is subjected to dynamic analysis due to the structural integrity. The internal structure shall be installed in the vertical hole call IR1 of reactor core. In order to verify the deflection of the internal structure, the mode and response spectrum analysis of the internal structure was performed. The natural frequency of the internal structure is 11.6 Hz(mode 1 and 2) and deflections of the internal structure are less than values of allowable design (3.2 mm).

As the power plant industry has been changed into competition structure, power generation companies do more with less by increasing capacity and lowering operation costs. In order to achieve this goals, an on-line real-time performance monitoring system has been needed to introduced to fossil power plant. The system represents a suite of related software modules which consist of on-line data, and on-line performance modules. This system can help the plant staff get the most out of their facilities by continuously monitoring deviations in equipment performance and the impact on those deviations on plant power, heat rate and operating cost. This paper shows the comparison of design value with acceptance test and current(measured) value.

The energy harvesting using smart materials has been extensively investigated to supply electric power to wireless sensor systems. In this paper, the energy harvesting using eddy current was studied with the integrated magnetic cantilever beam system. If a large conductive metal plate moves through a magnetic field which intersects perpendicularly to the sheet, the magnetic field will induce small rings of current which will actually create internal magnetic fields opposing the change. This eddy current that was induced in the coiled conductive sheet from the mechanical vibration was converted to chemical energy by charging batteries. The experimental results show that the eddy current generated the electric power up to max 31.2mW. Additionally the vibration reduction of the mechanical cantilever beam was observed by the energy dissipation in the electro-magnetic coupled system. The present result shows that the vibration level of the first natural frequency was reduced up to 7.7dB

Nuclear power plant has many external fire hydrants that have to operate in the state of emergency such as facility fire, forest fire. The valve stem of one among them was broken 3 times for 4 years. It had long valve stem and operated under high water pressure. The elongation and the tensile strength for the broken valve stem was measured to examine the defect of material property. And the vibration level and the natural frequencies was detected to check the resonance. As the result of a diagnosis, the cause of this fault is proven buckling of long valve stem.

An analysis program of specific impulse has been developed for a gas generator cycle rocket engine. The program has been verified by comparing the published performance data of the same cycle engine with RP-1 as fuel. A model for pressure drop of regenerative cooling and film cooling mass flow rate has been suggested to satisfy the necessary cooling condition with Jet-A1 as fuel. The engine mixture ratio is defined by the film cooling mass flow rate and the core mixture ratio. The optimal condition of the combustor pressure and engine mixture ratio has been found for maximum specific impulse.

In this study, the Peltier effect was applied to eliminate moistures in the air enclosed by a cabinet. We have developed the new electronic dehumidifier which has a new function of automatically evaporating the condensed water inner cabinet into the outside air. To obtain this function, the processes of dehumidification is that it condensed the moistures on the cold side heat sink and drained it into the hot side heat sink by the both gravitational and capillary forces and the droplets on the hot side heat sink surface was evaporated into the air outside the cabinet by the heat conducted through the hot side heat sink surface and the forced heat convection through the fan for cooling hot side heat sink. Compared to existing electronic dehumidifiers, this manufactured one showed a good performance that the electric power consumption for the same dehumidifying quantity was reduced by 50% compared with that of existing ones.

The advantages of the SITVC(secondary injection thrust vector control) technique over mechanical thrust vector systems include a reduction in both the nozzle weight and complexity due to the elimination of the mechanical actuators that are used in conventional vectoring. Computational study is performed to understand the fluidic thrust vectoring control of an axisymmetric nozzle, in which secondary gas injection is made in the divergent section of the nozzle. The nozzle has a design mach number 3. The effect of injection hole number and shape of secondary jet on the mach number distribution of SITVC were investigated. The standard ${\kappa}$ - ${\epsilon}$ turbulence model solved the complex three-dimensional nozzle flows perturbed by the secondary gas jet. The numerical code was validated by experiment. The results showed that the mach number distribution of circular and square nozzle are similar each other. As number of second injection hole increasing, a effect of deflection was decreased.