The purpose of the present study is to investigate the capillary heat transportation limitation in heat pipe according to the change of screen mesh wick porosity. Diameter of pipe was 6 mm, and mesh numbers are 100, 150, 200 and 250 and water was selected as a working fluid. According to the change of wick porosity and mesh number, the capillary pressure, pumping pressure, liquid friction coefficient in wick, vapor friction coefficient, and capillary heat transportation limitation are analyzed by theoretical design method of a heat pipe. As some results, the capillary heat transportation limitation in screen mesh wick heat pipe is largely affected by wick porosity and mesh number.

heat pipe with screen meshed wick. The heat pipe was designed in 200 screen meshs, 500 mm length and 12mm O.D tube of copper, water as working fluid and nitrogen as non-condensible gas. Heat pipe used in this study has evaporator, condenser and adiabatic section, respectively. Experimental data of wall temperature distribution along axial length is presented for heat transport capacity, condensor cooling water temperature change, degrees of an inclination angle, and operating temperature.

Heat pipes with binary mixture fabricated and tested for applications where condenser temperature is in a range of $10^{\circ}C$ to $130^{\circ}C$. The pipe materials 8.0 mm O.D. cupper tube and the working fluids are ethanol-water mixtures. The total length of test of the heat pipe was 1710mm in which evaporator section was 1570mm, adiabatic section was 50mm and condenser section was 90mm. Mixing ratios of ethanol and water could be variable in mole fraction. Temperature of condenser section was $10^{\circ}C$, $80^{\circ}C$ and $130^{\circ}C$. Heat pipe performance experimental study was accomplished with change of mixing ratio in these temperatures. The fill charge ratio was 20% of the heat pipe volume. Wick structure was woven-wire and method of experimental work was that thermal load was increased 20W step until the heat pipe wall temperature reached at $150^{\circ}C$. Results were following: At coolant $10^{\circ}C$ and $130^{\circ}C$, mixing ratio that have beat thermal performance was 0.8M+ and at coolant $80^{\circ}C$, was 0.3 ${\sim}$ 0.5 M+.

This paper is to research the heat transfer characteristic in copper-water variable conductance heat pipes(VCHP) with a non-condensable gas and gas reservoir. The heat transfer characteristics in the VCHP have not yet been studied much researches. VCHP are used in many applications. These applications range from thermal control of components and systems on satellites, to precise temperature calibration duties, conventional electronics temperature control and thermal diodes. The practical use of VCHP is a simple way to control the temperature of satellites. As the quantity of NCG was increased, there was an increase in the saturation vapor temperatures. As the input heat has loaded from 90 W to 110 W, the difference of the evaporator surface is lower than $10^{\circ}C$.

Several desiccant cooling systems have been developed in terms of cost and performance. In this study a fin-tube exchanger has been used for liquid desiccant dehumidification system. This dehumidifier has been designed to study the absorption characteristic of the aqueous triethylene glycol(TEG) solution which has the flow range from 20 to 50 LPM. The dehumidifier performance characteristic of working factor variables such as inlet solution flow rate, air flow rate, solution concentration, solution temperature, brine temperature, air temperature and inlet air relative humidity has been analyzed. The result of this experiment can provide useful data for hybrid air conditioning system.

Recently, the power consumption and heat generation in an electronic equipment increase as the components are miniaturized and the computing speed becomes faster. Effective cooling method is required to ensure the guaranteed performance and reliable operation of the electronic devices. The aim of the present study is to investigate the cooling characteristics of a pin-fin heat exchanger as a candidate for cooling system of the electronic devices. Various configuration of the pin-fin array is selected in order to find out the effect of spacing and diameter of the pin-fin on the heat transfer characteristics. The results are compared with the experimental data or correlations of several researchers for the channel flow with pin-fin arrays.

An experiment was performed to study heat transfer characteristics between corrugated heat exchanger and flat plate type one. While heat capacity(13.86kW) was provided constantly and the flow speed was varied from 2.8 to 17.9m/s, The temperature and the pressure drop were measured. Furthermore, Heat transfer coefficient, Colburn factor and Nusselt number were calculated using them. With increase of the flow speed for both exchangers, the coefficient and the pressure drop increased, but Colburn factor decreased. The coefficient, pressure drop and Colburn factor of the corrugated type were all higher than those of the flat one, which is due to the flow interruption with recirculation and reattachment of the corrugated type. The empirical correlations of Nusselt number were suggested for the tested two heat exchangers.

In this study the optimization of plate-fin type heat sink for the thermal stability is performed numerically. The optimum design variables are obtained when the temperature rise and the pressure drop are minimized simultaneously. The flow and thermal fields are predicted using the finite volume method and the optimization is carried out by using the sequential quadratic programming (SQP) method which is widely used in the constrained nonlinear optimization problem. The results show that when the temperature rise is less than 34.6 K, the optimal design variables are as follows; $B_{1}$ = 2.468 mm, $B_{2}$ = 1.365 mm, and t = 10.962 mm. The Pareto optimal solutions are also presented for the pressure drop and the temperature rise.

Scale is formed when hard water is heated or cooled in heat transfer equipments such as heat exchangers, condensers, evaporators, cooling towers, boilers, and pipe walls. When scale deposits in a heat exchanging surface, it is traditionally called fouling. The objective of the present study was to compare the fouling characteristics of river and tap water in a heat exchanging model. FromtheSEM analyses for tap water the $calciteformofCaCO3_{3}$ was formed. For river water, however, the $aragoniteCaCO_{3}$ wasformed.In order to investigate velocity effects on the fouling characteristics in the heat exchanging model, the inlet velocity was varied with 0.5, 1.0 and 1.5 m/s, respectively. The fouling characteristics of river water were quite different from those of tap water. For the case of the 'velocity of 1.5m/s', the overall heat transfer coefficient was reduced up to 26% than that of the 'velocity of 0.5m/s'

As operating time of heat exchangers progresses, fouling generated by water-borne deposits increases and thermal performance decreases. The fouling is known to interfere with normal flow characteristics and reduce thermal efficiencies of heat exchangers. The heat exchangers of nuclear power plants have been analyzed in terms of the heat flux and heat transfer coefficient at test conditions based on the ASME OM-S/G-Part 2 as a means of heat exchanger management. It is hard to estimate the heat performance trend and to establish the future management plan. This paper describes the fouling evaluation method which can evaluate the thermal performance for heat exchangers and estimate the future fouling variations and the plugging margin evaluation method which can reflect the current fouling level developed in this study. To develop the fouling and plugging margin evaluation methods for heat exchangers, fouling factor was introduced based on the ASME O&M codes and TEMA standards. For the purpose of verifying the two evaluation methods, the fouling and plugging margin evaluations were performed for a component cooling heat exchanger in a nuclear power plant.

Fouling is very serious problem in heat exchanger because it rapidly deteriorates the performance of heat exchanger. Cross-flow heat exchanger with vortex generators is developed, which enhance heat transfer and reduce fouling. In the present heat exchanger, shell and baffle are removed from the conventional shell-and-tube heat exchanger. The naphthalene sublimation technique is employed to measure the local heat transfer coefficients. The experiments are performed for single circular tube, staggered array tube bank and in-line array tube bank with and without vortex generators. Local and average Nusselt numbers of single tube and tube bank with vortex generator are investigated and compared to those of without vortex generator.

Numerical calculation has been performed to investigate the effect of inlet configuration on the growth rate of GaN layer on the heated susceptor. The conventional single inlet, where the gas is mixed by force in the inlet, is compared with separated flow inlet. Two-parallel gas flow $H_{2}$ and $NH_{3}$ are separated by a plate with finite length which are also parallel to the susceptor. The effect of separated plate length, carrier gas and flow rate of each precursor on the mixing of reactant gases and growth rate were investigated. Furthermore the three dimensional model is employed to predict the transverse variation of growth rate.

The InP thin films grown by metalorganic chemical vapor deposition (MOCVD) are widely used to optoelectronic devices such as laser diodes, wave-guides and optical modulators. Effects of various parameters controlling film growth rate such as gas-phase reaction rate constant, surface reaction rate constant and mass diffusivity are numerically investigated. Results show that at the upstream region where film growth rate increases with the flow direction, diffusion including thermal diffusion plays an important role. At the downstream region where the growth rate decreases with flow direction, film deposition mechanism is revealed as a mass-transport limited. Mass transport characteristics are also studied using systematic analyses.

In this study we propose an unsteady I-dimensional model of an iron ore sintering bed with multiple solid phases, which confers a phase on each solid material. This model contains coke combustion, limestone decomposition, gaseous reaction, heat transfers between each phase, and geometric changes of the solid particles. Simulation results are compared with the limited experimental data set of various coke contents and air supply rates. Effect of the coke diameter is also evaluated. They predict the experimental results well and show applicabilities to the various system of the fuel bed with various solid materials.

In a closed square cavity filled with a liquid, a cooled the upper horizontal wall and a heated the lower horizontal wall, the flow isn't generated under the ground-based condition when Rayleigh number is lower than 1700. In such case the flow phenomena near an air bubble under a cooled horizontal wall were investigated. The temperature and the flow fields were studied by using the Thermo-sensitive Liquid-Crystal and the image processing. The qualitative analysis for the temperature and the flow fields were carried out by applying the image processing technique to the original data. Injecting bubble at the center point of upper cooled wall, the symmetry shape of two vortexes near an air bubble was observed. The bubble size increased, the size of velocity and the magnitude of velocity increased. In spite of elapsed time, a pair of two vortexes was the unique and steady-state flow in a square cavity and wasn't induce to the other flow in the surround region.

Numerical analyses have been performed to obtain the absorption heat and mass transfer coefficients and the absorption mass flux from a falling film of LiBr solution. In the present study, the behavior of laminar-wavy falling film in the vertical absorber was studied analytically and experimentally. The change of absorption performance on mean film thickness, wave amplitude, wave celerity was analysed. The heat and mass transfer equations are solved simultaneously to give the temperature and concentration variations at the LiBr solution/refrigeration vapor interface and at the wall. Effects of uniform film, wavy film and film Reynolds number on the heat and mass transfer coefficients have been estimated. The analytical results of the uniform and wavy falling film in the bare tube was higher than experimental result for $Rd_{t}<100$. The absorption performance showed the maximum at the wavy film by the insert device(spring).

An experimental study is conducted to investigate the effects of duct corrugation angle on heat/mass transfer characteristics in wavy ducts by using a naphthalene sublimation technique. The corrugation angles of the wavy ducts are $145^{\circ}$ , $130^{\circ}$ and $115^{\circ}$ . and the Reynolds numbers based on the duct hydraulic diameter vary from 300 to 3,000. At the low $Re(Re{\leq}1000)$, high heat/mass transfer regions are formed by the secondary vortex flows called Taylor-Gortler vortices on both pressure-side and suction-side walls. At the high $Re(Re{\geq}1000)$, the effects of these secondary flows are vanished. As corrugation angle decreases, the local peak Sh induced by Taylor-Gertler vortices are increased and average Sh also enhanced. More pumping power (pressure loss) is required with the smaller corrugation angle due to the stronger secondary vortex flows.

In this paper, the current lead for superconducting device is studied by numerical method. The current lead is cooled by surrounded $N_{2}$ gas by natural convection. The heat conduction equation for current lead and boundary layer equation for $N_{2}$ gas must be solved simultaneously. The boundary layer equation for $N_{2}$ gas is highly nonlinear for varied temperature of current lead. So the linearization method is adopted for simplicity. Numerical results using natural convection cooling are compared with the conventional cooling methods such as conduction cooling and vapor cooling methods. The main difference of natural convection cooing is the non-zero temperature gradient at the top of current lead for the minimum heat dissipation into superconducting devices. For the optimized conduction-cooling and vapor-cooling current leads, the temperature gradient at the top of current lead is zero. Also, the heat flow at the cold end is much smaller than conduction cooling case.

The physical model considered here is a horizontal layer of fluid heated below and cooled above with a periodic array of evenly spaced square cylinders placed at the center of the layer, whose aspect ratio here varies from unity to twelve. Periodic boundary condition is employed along the horizontal direction to allow for lateral freedom for the convection cells. Two-dimensional solution for unsteady natural convection is obtained using an accurate and efficient Chebyshev spectral multi-domain methodology for a given Rayleigh numbers of $10^{6}$.

In this study, smoke movement in tunnel fires was investigated with various aspect ratio(0.5, 0.667, 1.0, 1.5, 2.0) of tunnel cross section. Reduced-scale experiments were carried out under the Froude scaling using 8.27 kW ethanol pool fire. Temperatures were measured under the ceiling and vertical direction along the center of the tunnel. Smoke front velocity and temperature decrease rate were reduced as higher aspect ratio of the tunnel cross-section. Smoke movement was evaluated by analysis of vertical temperature distribution 3 m downstream from the fire source. Elevation of smoke interface according to N percent rule was under about 60% of tunnel height.

A study of argon droplet vaporization is conducted using molecular dynamics, instead of using traditional methods such as the Navier-Stokes equation. Molecular dynamics uses Lagrangian frame to describe molecular behavior in a system and uses only momentum and position data of all molecules in the system. So every property is not a hypothetical input but a statistical result calculated from the momentum and position data. This work performed a simulation of the complete vaporization of a three dimensional submicron argon droplet within quiescent environment. Lennard-Jones 12-6 potential function is used as a intermolecular potential function. The molecular configuration is examined while an initially non-spherical droplet is changed into the spherical shape and droplet evaporates. And the droplet radius versus time is calculated with temperature and pressure profile.

A study of high-pressure n-heptane droplet vaporization is conducted with emphasis placed on equilibrium at vapor-liquid interface. General frame of previous rigorous model[1] is retained but tailored for flash equilibrium calculation of vapor-liquid interfacial thermodynamics. The model is based on complete time-dependent conservation equations with a full account of variable properties and vapor-liquid interfacial thermodynamics. The influences of high-pressure phenomena, including ambient gas solubility, thermodynamic non-ideality, and property variation on the droplet evaporation are investigated. The governing equations and associated moving interfacial boundary conditions are solved numerically using a implicit scheme with the preconditioning method and the dual time integration technique. And a parametric study of entire droplet vaporization history as a function of ambient pressure, temperature has been conducted. Some computational results are compared with Sato's experimental data for the validation of calculations. For low ambient temperatures, the droplet lifetime first increases with pressures, then decreases for high pressures. For higher ambient temperatures, the droplet lifetime increase with less amplitude than that of low ambient temperatures, which then decreases with more amplitude than that of low temperatures. The solubility of nitrogen can not be neglected in the high pressure and it becomes higher as the pressure goes up.

In this study, numerical investigation has been performed on the evolution of weld pool geometry with moving free surface during low-energy density laser welding process. The free surface elevates near the weld pool edge if ${\partial}{\sigma}/dT$ is dominantly negative. On the contrary, the free surface rise at the center of weld pool in case of mainly positive ${\partial}{\sigma}/dT$. The predicted weld pool width and depth with moving free surface are 5∼15%$5{\sim}15%$ greater than those with flat weld pool surface. It is considered that weld pool surface oscillation during melting process augments convective heat transfer rate in the weld pool.

Laser melting problems with deformed substrates are investigated by axisymmetric numerical simulations. Source-based method is used to solve the energy equation, and the momentum equations are solved in the liquid domain with SIMPLER algorithm. Using a laser beam with a top-hat heat flux distribution, this study is performed to examine the effect of surface deformation, beam power density and surface tension force on the melt pool during laser melting. Surface temperature decreases with increasing surface deformation, while surface velocity increases. It is found that surface deformation, beam power density and surface tension force have a very significant effect on heat transfer and fluid flow during laser melting.

In this paper we have studied the characteristics of temperature fields in a Molded Case Circuit Breaker (MCCB). A switch and a trip device are arranged in narrow space of the MCCB. Thus, thermal factors have been risen to important problem. The temperature rise means an energy loss and it becomes the reason of fatal fault in devices. Also, the temperature rise on the connection and the contact parts is regulated up to $115^{\circ}C$ and $105^{\circ}C$, respectively. Therefore, a study for preventing the temperature rise should be investigated. A numerical analysis method that has confidence might be preceded for this purpose. In order to verify the confidence of temperature field analysis, the results of the numerical analysis are compared with those of experimental one for the same model. Comparison results show a qualitatively good agreement within ${\pm}5%$ errors.

We performed the numerical study on the characteristics of the airflow and cross contamination in the photolithography process cleanroom. The nonuniformity, the deflection angle and the global cross contamination were used for analyzing the characteristics and performances of cleanroom. We knew that the airflow characteristics of the cleanrooms are largely affected by the porosity of panel and the adjustment of dampers. And the numerical result showed that the global cross contamination varies with the location of source and the passage of time.

This study is analysis on the heat transfer of exhaust gas into the muffler at passenger Car. Numerical analysis with Computational fluid Dynamics(CFD) was carried out to investigate exhaust gas flow. The STAR-CD S/W used for the three dimensional steady state CFD analysis in a muffler. The Navier-Stokes Equation is solved with the SIMPLE method in a general cartesian coordinates system. Result of numerical simulation; Inlet and outlet temperature shown about ${\Delta}T=239K$, 216K, 202K at in the muffler. Heat transfer was progressed quickly by atmospheric temperature of muffler external at in the near wall.

Analytical approach is applied to predict temperature of satellite box under worst hot condition from fairing jettison to separation. The box is tried to solve analytically which is exposed to known environmental heating condition and has internal heating and irradiation to space. For a single thermal mass, governing equation for temperature is simplified to 1st order ordinary differential equation(ODE) by several assumptions. Two cases are considered. The one is for constant mass box and the other is for mass-varying box. Each case has three different analytical solution by sign of constant term in ODE. One analytical solution for constant mass is applied to physical launch stage condition. It is concluded that the present analytical method can be used to quick predict the temperature of a satellite box and check whether a satellite is safe against space environment during launch stage.

The thermal design of the bus bar of a Gas-Insulated Switchgear(GIS) becomes important since the current-carrying capacity of the GIS is limited by maximum operating temperature. In order to predict temperature rise of the bus bar, a program has been developed. Various heat sources possibly generated in the bus bar are calculated in the program. To estimate temperature rises at the bus bar caused by the heat balance between the heat generation and heat transfer, the finite volume method as well as the $4^{th}$ order Runge-Kutta method has been employed. In the experiments, temperature rises at conductor, contact part and external tank are measured for full-scale gas-insulated bus bars. The comparisons of the predicted values of the heat balance calculation to those of the experiments are made. From the comparisons, it is concluded that the developed program can predict the temperature rise of the bus bar quite well.

An inverse problem to determine two-dimensional total heat exchange factor is studied for the prediction of the billet temperature in the reheating furnace. Temperature measurements by the experiment are used in the inverse analysis. This inverse analysis employs the conjugate gradient method. The total heat exchange factors for 12-zones of the cross-section of the billet are estimated. The estimated temperatures at measurement locations are in good agreements with the measured temperatures.

An experimental study has been conducted to measure the temperature fields for two and three staggered rows of the rectangular shaped-holes with high turbulence intensity. 10 % turbulence intensity is obtained by installation of two kinds of grids which have different shapes. One grid which is installed at 30d upstream from center of 1st row of holes is composed of vertical cylinders of which diameter is 10 mm and center to center distance is 18 mm. The other installed 15d apart to upstream from center of 1st row of holes which has square pattern is constructed of 3 rum square bars and bar spacing is 25 mm. Temperature fields are measured by using a thermocouple rake which is attached on three-axis traversing system. The results show that the overall values are decreased and the thicker film of coolant is fanned downstream of rows of holes for high mainstream turbulence intensity.

The heat (mass) transfer characteristics in the tip-leakage flow region of a high-turning first-stage turbine rotor blade has been investigated by employing the naphthalene sublimation technique. The heat transfer data in the tip-leakage flow area for the tip clearance-to-span ratio, h/s, of 2.0% are compared with those in endwall three-dimensional flow region without tip clearance (h/s = 0.0 %). The result shows that the thermal load in the tip-leakage flow region for h/s = 2.0% is more severe than that in the endwall flow region for h/s = 0.0%. The thermal loads even at the leading and trailing edges for h/s = 2.0% are found larger than those for h/s = 0.0%. The tip-leakage flow results in heat transfer augmentations near the tip on both pressure and suction sides in comparison with the mid-span results.

A gas turbine blade has an internal cooling passage equipped with ribs, which can be modeled as a ribbed channel. We have studied a flow inside a ribbed channel using large eddy simulaton (LES) with a dynamic subgrid-scale model. The simulation results are compared with the experimental ones. The turbulence intensity and local heat transfer near the rib have not been well captured by the conventional Reynolds averaged Navier-Stokes simulation (RANS). However, these variables obtained by the present LES agree well with those from experiments. From the instantaneous velocity and temperature fields, we explain the mechanisms responsible for the local peaks in the heat transfer distribution along the channel wall. We have also investigated the effect of rotation on the flow and heat transfer in the ribbed channel.

In the present study, the compression process in scroll compressor was simulated in consideration of flow leakage and heat transfer. Tangential and radial leakages of the refrigerant between the scrolls were considered as nozzle flow. The experiment was first conducted with a scroll compressor for automobile air conditioning system and R134a as a refrigerant. Temperature and pressure were measured at the suction and discharge ports of the compressor to determine the thermodynamic states of the refrigerant flow. Temperature distribution of the scroll with the involute angle was also measured by thermocouples that were installed inside the scroll. Measured temperature distribution was compared with the numerical results. From this result, the thermal effect of mechanical contact was found to be important in heat transfer of the compression process.

The condensation water temperature $(T_{CW})$ distributions on the combi-drum washer of condensation type which the cloth washes, rinses, spins and dry cloth over a heater were measured using the thermocouples. During the drying process, in order to find out the relationship between the $T_{CW}$ distribution and the cloth dryness level $(C_{D1})$, the experiment was made for the dehydration rate of the cloth $(C_{D0})$ is 63%, the weight of cloth $(W_{C0})$ from 1 to 6kg. In result, according to $W_{C0}$ and $C_{D0}$, the distributions of $T_{CW}$ are very different. Using this tendency, it was appropriate and it developed the automatic drying algorithm and it was gotten dryness level regardless of $W_{C0}$, $C_{D0}$ and the other factor.

A fine hot-wire is used both as a heating element and a temperature sensor in transient hot-wire method. The traditional sensor system is unnecessarily big so that it takes large fluid volume to measure the thermal conductivity. To dramatically reduce this fluid volume, a new sensor fabrication and a data processing method are proposed in this article. Contrast to the conventional and most popular two wire sensor, the new sensor system is made up of divided multiple long and short wires. Through validation experiments, it is found that the measured thermal conductivities of the glycerin are exactly same each other between the conventional and proposed new method. Also some technical considerations in arranging the multiple wires are briefly discussed.

A thermal vacuum chamber is mainly used to simulate thermal environments of a test satellite in satellite orbits in which daily temperature variations range from 80K to above 400K depending on solar radiation under the vacuum below $10^{-4}$ torr. The test facility is quite complex and consists of expensive parts. So any modification of control software is discouraged in fear of unexpected system failure. The purpose of this study is to develop a realtime dynamics model of the thermal vacuum chamber in view of controller design and simulate its electrical inputs and outputs for interface with a PLC (programmable logic controller). A PLC program that was used in the thermal vacuum chamber is applied to the realtime simulator. The realized simulator dynamics is found to be quite similar to that of the thermal vacuum chamber and serve to an appropriate plant to verify the control performance of a programmed PLC.

Recent advances in room-temperature, visible and near-IR diode laser sources for telecommunication, optical data storage applications are enabling combustion diagnostics system based on diode laser absorption spectroscopy. In contrast to some traditional sampling-based gas-sensing instruments, tunable diode laser absorption spectroscopy system is advantageous because of their non-invasive nature, high sensitivity, fast response time and real-time measurement capability. So, combined with fiber-optics and high sensitive detection strategies, compact and portable sensor system arc now appearing for a variety of applications. The objective of this research is to take advantage of distributed feed-back diode laser and measure the $CO_{2}$ concentration (by using direct absorption and wavelength modulation spectroscopy methods). In addition to survey spectra of $CO_{2}$ bands and spectroscopic parameters between 1565 and 1579 run were computed at temperatures between 296 and 1200 K (by using HITRAN 2000 database). It experimentally found out that the features of direct absorption and wavelength modulation spectroscopy methods.

It has recently been shown that the instability of thin film of a nanoscale can be used in the processes of building nano-size structures, which have potential practical importance in nanotechnology. Molecular dynamics simulation is conducted to probe the thin fluid film of a nano-size and its dynamic behavior during destabilization and structure formation. Non-continuum characteristics are shown in the properties like pressure tensor, viscosity, and thermal conductivity. The thermocapillary force induces a slow growth of long waves in the scale considered. A long-range interaction with the solid wall induces vertical structures, whose formation time and space between neighbors are proportional to the strength of the interaction.

A brand-new micro actuator is introduced in this paper. This device is one of thermo-pneumatic actuators, and based on two distinct principles of snap-through buckling and phase change. These coupled phenomena affect each other positively and will outrun the performance of an ordinary thermo-pneumatic actuator. Our efforts are focused on comprehensive analysis on the driving force of the actuator. For the analysis, we explain each principle and offer approximated models for the buckling and phase change. The calculation results from each model are compared to experimental data. The comparison between prediction from models and data from experiments is within the satisfaction in spite of a lot of approximations.

This study concerns the performance of condensing heat transfer in two-phase closed thermosyphons with various helical grooves. Distilled water, methanol, ethanol have been used as the working fluid. In the present work, a copper tube of the length of 1200mm and 14.28mm of inside diameter is used as the container of the thermosyphon. Each of the evaporator and the condenser section has a length of 550mm, while the remaining part of the thermosyphon tube is adiabatic section. A experimental study was carried out for analyzing the performances of having 50, 60, 70, 80, 90 helical grooves. A plain thermosyphon having the same inner and outer diameter as the grooved thermosyphons is also tested for the comparison. The type of working fluid and the numbers of grooves of the thermosyphons with various helical grooves have been used as the experimental parameters. The experimental results have been assessed and compared with existing theories. The results show that the type of working fluids are very important factors for the operation of thermosyphons. And the maximum enhancement (i.e. the ratio of the heat transfer coefficients the helical thermosyphons to plain thermosyphons) is $1.5{\sim}2$ for condensation.

Vapor-liquid equilibrium data were obtained for system of propane + R227ea (Heptafluoropropane) over the temperature range from 253.15 K to 323.15 K at 10 K intervals. Experiments were performed in a circulation type apparatus by injecting vapor through liquid pool using a magnetic pump. This system forms azeotrope in the temperature range of this study. The experimental results were correlated with the Peng-Robinson (PR) equation of state and Redlich-Kwong-Soave (RKS) equation of state using the van der Walls one-fluid mixing rule and were compared with each other.

TMA clathrate that is used as PCM of low temperature thermal storage system in this research creates hydrate crystallization at higher temperature then pure water, and its application is expected as PCM because of comparatively big latent heat without phase separation phenomenon. Acetone, Ethylen Glycol, and Ethanol is used as additive and evaluated experimentally for the purpose of the improvement in subcooling of TMA clathrate. In view of the results so far achieved subcooling is improved, the running time of the refrigerator is reduced. Thus the results are expected to use for the increase of coefficient of performance of low temperature thermal storage system in the building.

A computer program to calculate thermodynamic properties of oxygen is developed. Procedures for the calculation is briefly discussed. The program calculates unknown thermodynamic properties fixing the state with two independent input properties. If input value by user is inappropriate, it displays an error message. In addition user can change units with easy. The program developed in this work can be utilized to calculate parameters required for the simulation and design of an equipment using oxygen.

As motor performance enhancement by improving electric design has reached its limit and downsizing issue has risen, the importance of thermal design is increasing. In this study, the flow and temperature distribution were reviewed with the help of CFD analysis and this result was compared with the experimental results. Furthermore, parametric analysis with thermal design structure showed that axial duct width but fan capacity is a critical factor to lower the hot spot temperature in electric motor.

The numerical study was carried out to effectively cool Induction motor applied to a washing machine. The outer rotor made of steel periodically spins up and down. The stator consists of the thin layered iron plates and copper coil. The effective cooling system is necessary to solve the reliability problem caused by the electric losses at the coil and the iron plate. Because the heat transfer rate of the natural convection in partially open space is generally low, thus it is necessary to enhance the heat transfer using rotating perforated plate. The flow phenomena around the motor are very complex due to the motor geometry and the outer rotor motion. The mixed convection takes place due to the slow rotation speed. The three dimensional flow simulation was performed using rotating reference frame technique and Boussinesq approximation but the radiation effect was neglected. It was found that the angle and direction of the cooling blades play an important role in the stator temperature.

This paper reports the fluid flow and heat transfer around a module cooled by forced air flow generated by a piezoelectric(PZT) cooling fan. A flexible PZT fan with distortion in a fluid transport system of comparatively simple structure which was mounted on a PCB in a parallel-plate channel($450{\times}80{\times}700mm^3$) accelerates surrounding fluid locally. Input voltages of 20-100V and a resonance frequency of 23Hz were used to vibrate the cooling fan. Input power to the module was 4W. The cooling effect using a PZT fan was larger than that of free convection. Fluid flow around the module were visualized by using PIV system. The temperature distribution around heated module were visualized by using liquid crystal film(LCF). We found that the flow type was y-shaped and the cooling effect was increased by the wake generated by a piezoelectric cooling fan.

Analysis has been carried out to investigate the temperature variation and the uniformity of the temperature distribution of the glass panel by infrared radiant heating. Halogen lamps are used to heat the panel and located near the top and bottom of the rectangular chamber. The thermal energy is transfered only by radiation and the radiation exchange occurs only on the solid surfaces and is considered by using the view factor. The results show that the uniformity of the temperature distribution of the panel is improved but the time for heating increases as the wall reflectivity is large. The temperature difference reaches a maximum in the early stage of the heating process and then decreases until it reaches the uniform steady-state value.

The narrow band-averaged transmissivity of $CO_2-H_2O$ mixtures is expressed by multiplying the transmissivities of $CO_2$ and $H_2O$. Applying the multiplication property of narrow band transmissivities for gas mixtures of $CO_2-H_2O$, the number of gray gases, required for accurate representation of the absorption characteristics by using the narrow band based WSGGM, is significantly increased. To reduce the computational loads by reducing the number of gray gases, we propose a gray gas regrouping process where the gray gases used for the WSGGM are regrouped into a specified number of groups according to the magnitudes of absorption coefficients. To evaluate the proposed WSGGM for gas mixtures, the radiative transfer problems through three-dimensional gas media are considered. The radiative source terms and the radiative heat fluxes obtained by using the proposed method are fairly well compared to those obtained by using the SNB model. The regrouping technique results in an excellent computational efficiency with minor loss of accuracy.

Among the thermal analysis methods for Multi-Layer Insulation(MLI), effective emittance, diffusion MLI node and arithmetic MLI node methods are compared. The methods have been applied to the aluminum panel under the low earth orbit environment. TRASYS program is used for geometrical math modeling and SINDA program for thermal math modeling and temperature calculation. Test cases are selected according to MLI area on the panel. Temperature results are calculated and compared under the ratio of absorptivity and emissivity.

A study on the engineering design and numerical thermo-hydraulic analysis for KSTAR TF coil structure cooling system has been conducted. The numerical analyses have been done to verify the engineering design of cooling using the commercial code, FLUENT and in-house code for calculating helium properties which varies with cooling tube's heat transfer. Through the engineering design process based on the steady heat balance concepts, the circular stainless steel tube with inner diameter of 4 mm for TF coil has been selected as cooling tube. From normal operation mode analysis results, total 28 cooling tubes were finally chosen. Also, three dimensional cool down analysis for TF coil with designed cooling tube was satisfied with next three design criteria. First is cooling work termination within a month, second is maximum temperature difference within 50 K in TF coil structure and third is exit helium pressure above 2 bar. Consequently, these cool down scenario results can afford to adopt as operating scenario data when KSTAR facilities operate.

Demand of the LCD-TV is remarkably increasing with development of the LCD technology in these days. Thus, this research has analyzed thermal problems such as heat transfer characteristics inside and outside the LCD-TV using numerical simulation and experiments. The simulated results have been compared with the experimental results using an infrared (IR) camera and T-type thermocouples. The optimal design of structure has been proposed to improve the thermal efficiency of radiation from the comparison.

Thermal mass flow meters (TMFMs) are most widely used for measuring mass flow rates in the semiconductor industry. A TMFM should have a short response time in order to measure the time-varying flow rate rapidly and accurately. Therefore it is important to study transient heat transfer phenomena in the sensor tube of a TMFM that is the most critical part in the TMFM. In the present work, a simple numerical model for transient heat transfer phenomena of the sensor tube of a TMFM is presented. Numerical solutions for the tube and fluid temperatures in a transient state are obtained using the proposed model and compared with experimental results to validate the proposed model. Based on numerical solutions, heat transfer mechanism in a transient state in the sensor tube is explained. Finally, a correlation for predicting the response time of a sensor tube is presented. The correlation is verified by experimental results.

The concept of lean-premixed combustion in gas turbine combustor operation has become a standard in recent years as an effective means to meet stringent enviromental standards on NOx emissions. The combustion characteristics of 75 kW class lean premixed combustor were investigated at the conditions of high temperature and ambient pressure. The exit temperature and emissions of CO and NOx were measured at the center of exit plane. The high temperature air of $550K{\sim}650K$ was supplied through air preheater. As expected, experimental results indicate that NOx emission was increased and CO emission was decreased by increasing inlet air temperature. But CO emission measured at the center of exit plane was increased because of the non-uniform radial direction profiles. The Semi-Empirical Correlation method was applied to obtain the design point emissions of NOx and CO. Also the flame temperature, CO and NOx emissions were measured along the centerline of liner at 650K inlet air temperature to determine the position of dilution holes.

Acoustic-Pressure Response of diluted hydrogen-air diffusion flames is investigated numerically by adopting a fully unsteady analysis of flame structures. In the low-pressure regime, the amplification index remains low and constant at low frequencies. As acoustic frequency increases, finite-rate chemistry is enhanced through a nonlinear accumulation of heat release rate, leading to a high amplification index. Finally, the flame responses decrease at high frequency due to the response lag of the transport zone. For a medium-pressure operation and low-frequency excitation, the amplification index is low and constant. It then decreases at moderate frequencies. As frequency increases further, the amplification index increases appreciably due to an intense accumulation effect.

Combustion using oxygen enriched air is known as a technology which can increase flame stability as well as thermal efficiency due to improvement of the burning rate. Lift-off, blowout limit and flame length were examined as a function of jet velocity, coflow velocity and OEC(Oxygen Enriched Concentration). Blowout limit of the flame below OEC 25% decreased with coflow velocity, but the limit above OEC 25% increased inversely. Lift-off height decreased with increase of OEC. Especially lift-off hardly occurred in the condition above OEC 40%. Flame length of the flames above OEC 40% was increased until the blowout occurred. Flame stability became improved since lift-off and blowout limit increased much with increase of OEC.

Oxygen-enriched non-premixed flame characteristics was investigated numerically with variation of dilution methods, which are Flue gas recirculation(FGR) and fuel injection recirculation(FIR). In order to compare flamelets in various oxygen-enrichment conditions reasonably, the adiabatic flame temperature and Damkoller number was held fixed by modulating amount of diluents to fuel and oxidizer stream and by varying global strain rate of flame respectively. Also modified GRI 3.0 reaction mechanism was utilized, which was able to predict oxygen-enriched methane flame correctly. Fundamental flame characteristics such as structure, heat release rate and extinction with FGR and FIR were compared in various oxygen enrichment conditions.

Computations were performed to investigate the spray characteristics of the twin fluid nozzle in three stage heavy-oil combustion burner. The burner geometry and flow conditions were provided by a burner company. The goal of the study is to estimate mean droplet size, initial velocity and spread factor of the nozzle through comparison between experiments and numerical analyses. Air stage ratio is 2:4:4 by mass, and O2 in exhaust gas is about 4 % by volume. Here, the agreement between the experiment and numerical analyses is evaluated by NOx generation. Spray characteristics will be linearly interpolated between fuel consumption rate l20L/h and 240 L/h.

The purpose of this study is to ascertain the ability of New type cutter using Hybrid Rocket Propulsion System to cut normal carbon steel and also compound metal like stainless steel which cannot be cut by regular oxygen-acetylene cutter. To compare cutting performance, Two different types of experiment with oxygen-acetylene and Hybrid Combustion cutters were performed. As a result, Hybrid Combustion cutter is used to cut both carbon steel and stainless steel with cutting speed of 400mm/min(carbon steel) and 250mm/min(stainless steel). Otherwise, oxygen-acetylene cutter can be used to cut only carbon steel with cutting speed of 500 $^{\sim}$ 700mm/min. The possibility of Hybrid Combustion cutter as a cutting machine was confirmed.

Flame onset and propagation within hydrogen premixed gas mixture are numerically investigated in an rectangular enclosure. A detailed chemistry for hydrogen reaction is applied to anticipate the thermochemical behavior of intermediate species appropriately. To facilitate computation, 10 species and 16 elementary reaction steps for hydrogen combustion are taken into account. On the basis of 30 % of hydrogen concentration in hydrogen-air mixture, the effects of position and quantity of ignition sources on the flame evolution are analyzed. From the simulation results, the means that can lessen the possible hazard caused by flame propagation are suggested.

The objective of this study is a qualitative comparison between line-integrated OH chemiluminescence ($OH{\ast}$) image and its Abel inversion image at different phase of the oscillating pressure field. PIV(Particle Image Velocimetry) measurements were conducted under non-reacting conditions to see the global flow structure. Also NOx emission was measured to investigate the effect of fuel-air premixing on combustion instability and emission characteristics. Experiments were carried out in an atmospheric pressure, laboratory-scale dump combustor operating on natural gas. Combustion instabilities in present study exhibited a longitudinal mode with a dominant frequency of ${\sim}341.8$ Hz, which corresponded to a quarter wave mode of combustor. Heat release and pressure waves were in-phase when instability occurred. Results gave an insight about the location where the strong coherence of pressure and heat release existed. Also an additional information on active control to suppress the combustion instabilities was obtained. For lean premixed combustion, strong correlation between $OH{\ast}$ and NOx emissions was expected largely due to the exponential dependence of thermal NOx mechanism on flame temperature.

It is shown that the effect of variable parameters on flame structures and NOx emissions in the laminar partially premixed methane-air flames with a co-axial Bunsen burner. Objectives of this paper is to understand the effects of flow variables on NOx emissions and the flame structure with OH chemiluminescence, including reconstructed image by abel inversion processing at each conditions. A fuel flowrate of 200 [cc/min] was fixed and the amount of air was varied from 400 to 1200 [cc/min]. The experimental variables were equivalence ratio(${\Phi}$ fuel split percentage(${\sigma}$ and inner tube recess(x/D). Flow conditions were ranged from $1.36{\sim}4.76$(equivalence ratio), $50{\sim}100$(fuel split percentage) and $0{\sim}20$(inner tube recess). NOx analyzer and ICCD camera with a OH filter were used as a main experimental apparatus. In addition, Abel inversion, which is a kind of tomography and valuable to estimate a two-dimensional structure of co-axial flames from cubical information, was employed for combustion diagnostics. Results from this study indicate that the main effects depend on equivalence ratio and next sigma, x/D for NOx production and OH formation. Throughout Abel inversion, we could affirm the maximum position and the tendency of OH radical intensity by variants at five axial heights above the burner exit.

There arc two methods to evaluate the combustion product, that is analysis method and measurement method. Properties of the combustion products, constituents, and quanties, can be used to calculate boiler efficiency and beat losses. Using combustion evaluation methods combustion calculation were performed at certain condition and analyze the results. In this paper by comparing the results between two methods, deviation, correction and possibility of site application were confirmed

It is important that overheat protection of super heated tube in boiler operation and maintenance. The overheat of super heat tube can make damage and rupture of tube material, which causes accidental shutdown of boiler. The super heated tube overheat is almost due to the lack of uniformity of gas temperature distribution. There are two ways to protect overheat of super heated tube. The one is to control hot gas operation pattern which is temperature or flow distribution. the other is to control super heated steam flow distribution. The former is difficult than the later, because of control device design. In this paper steam flow control method which uses orifices is proposed to protect overheat of super heat tube.

Coal-burning utilities are facing a major NOx control compliance challenge due to the heavy emission regulation. In response to this challenge, some applicative technologies to effectively reduce NOx are developed and applied in the pulverized coal power plants. One of these is low NOx burner(LNB) equipped with multi-staged air register. In this study, NOx emission rate and flame shapes are investigated with secondary and tertiary air flow conditions in air staged coal burner, and the optimal windows of flow conditions to minimize NOx emission rate are found out. The test conditions treated in this study are the flow rate, swirl direction and intensity and throat injection velocity of secondary and tertiary air.

Three-dimensional experimental analysis was conducted in the pulverizer simplified isothermal model. The experiment model was constructed on a 1/3.5 scale of 500MW pulverizer. The purpose of this study is to investigate the effect of design parameters on the pulverized coal separator efficiency. Where used pulverized coal separator design parameters are guide vane angle, static classifier angle, dynamic classifier rpm. Taguchi method was used to find the effective design parameters related to pulverized coal separator efficiency. The results of the experiment showed that guide vane angle and dynamic classifier rpm were the design key parameters. In addition to the total number of experiment cases were reduced by Taguchi method.

The effect of recirculated exhaust gas on exhaust emissions under four kinds of nozzle tip with the different fuel consumption rate are experimentally investigated by using an once-through boiler with FGR system. The purpose of this study is to develop the FGR control system for reducing NOx in a boiler. Intake and exhaust oxygen concentrations, and equivalence ratio are applied to discuss the effect of FGR rate on exhaust emissions at various fuel consumption rates. It is found that NOx emissions are decreased, while soot emissions are increased owing to the drop of intake and exhaust oxygen concentrations, and the rise of equivalence ratio as FGR rates are elevated.

In the present study, preliminary investigation were carried out by analysis of energy system(heat and electricity) based on phosphoric acid fuel-cell of 50 kW for eco-apartment houses. Analysis model were consisted of fuel cell energy system, secondary energy unit and residential building of 5 stories with 20 and 40 households. And the investigation results reviewed under load pattern of heat and electric power of the apartment houses. The results showed mismatch between the needed heat load pattern and output of fuel cell energy system. The mismatch rate were assessed about 10-180% of heat load for apartment houses with season. We found that secondary energy unit are needed in order to supply insufficient heat.

The main reason to select the maxium plant efficiency through the performance test in fossil power plant is to increase the efficiency of power plant as well as saving energy collated with the policy of government. This study is aimed at unerstanding the variantion trend of efficiency and analyzing the efficiency of boiler and turbine through each of the performance test. Ultimately, the maxium efficiency of power plant will be presented in super-critical pressure type power plant.

Naphthalene sublimation technique is used to investigate the average and local heat transfer from the circular rod, and to determine the average and local heat transfer from the circular rod with and without square wing type mixing vane in axial flow. The experiments are performed for a circular rod and flat plate with and without mixing vane in wind tunnel. In comparison with flat plate and circular rod in axial flow, averaged Nusselts number is increased 2 times as the increase of Reynolds number with mixing vane. Longitudinal vortex induced by square wing type has the stronger vortex strength, so square wing type vortex generator shows an effect further in downstream.

In the recent TROI experiments, melts of zirconia and two different compositions of corium were used to observe the occurrence of a steam explosion when it came into contact with water at two different depths. The compositions of the corium were 70 : 30 and 80 : 20 in weight percent of $UO_{2}$ and $ZrO_{2}$, and the mass of the corium was about 10kg. The depth of water in the interaction vessel was 67cm and 130cm. A steam explosion did not occur in the interaction between 80 : 20 corium melt and water at 130cm depth, while steam spikes were observed in the interactions between corium melts of two different compositions and water at 67cm depth. A strong steam explosion occurred in the interaction between 5.43kg of zirconia melt and water at 67cm depth. This fact shows that the explosivity of zirconia is much greater than that of corium.

A low temperature difference model Stirling engine is a small Stirling engine running with several degree of temperature difference without power output. In this study, the design parameters to give maximum power are discussed. As results, the phase angle is about 100 degree, and compression ratio is 1.5% of the ratio of heat source temperatures at maximum power condition.

Hot gas path components, which are made of nickel based superalloys, are subject to periodic replacement due to degradation of thermomechanical properties that might bring catastrophic failure during normal operation of gas turbine units. In order to rejuvenate the metallurgical condition of the serviced components, heat treating techniques such as solution annealing and aging heat treatments have widely been employed. However, the effectiveness of those typical heat treatments is not apparent enough in terms of quantitative grounds. On the other hand the demand of the rejuvenation heat treatment and hot isostatic pressing (HIP) have constantly been raised by the end users. Therefore it is necessary to verify how the typical heat treating techniques affect to the aged and degraded material. As the result of experimental work in this study, GTD-111 and GTD-222 Ni-based superalloys were collected and analyzed quantitatively through microscopic observation, microhardness evaluation and creep test.

Detailed heat transfer coefficient distributions on a gas turbine blade tip were measured using a hue-detection base transient liquid crystals technique. The heat transfer coefficients on the shroud and near tip regions of the pressure and suction sides of a blade were also measured. Both plane tip and squealer tip blade were considered. The heat transfer measurements were taken at the three different tip gap clearance of 1.0%, 1.5%, and 2.5% of blade span. Results show the overall heat transfer coefficients on the tip and shroud with squealer tip blade were lower than those with plane tip blade. However, the reductions of heat transfer coefficients near the tip regions of the pressure and suction sides were not remarkable.

The present study investigates heat/mass transfer characteristics in rotating two-pass ducts of three different aspect ratios with 90-ribbed surfaces. The results show that the flows generated by the 180-turn, rib turbulators, and duct rotation. The curvature of the 180-turn produces Dean vortices causing high heat/mass transfer in the turning region and in the upstream region of the second-pass. The rib turbulators disturb the main flow by producing reattachment and separation near the ribbed surfaces, and increase heat/mass transfer in the region between the ribs. As the rotation number increases, the heat/mass transfer discrepancy between the leading and the trailing surfaces become larger.

The characteristics of instantaneous wall-surface temperature of impinging plate in case of ultra high pressure injection have been measured and analyzed by using thin film instantaneous temperature probe and ultra high pressure injection equipment. The decreasing rate of temperature was greater in case of higher temperature of impinging plate. Temperature drop was largest at center of piston and it was slight for others. Instantaneous temperature decreases rapidly with increasing injection pressure. But above 2,500bar of injection pressure, the decreasing rates are slightly affected by increasing injection pressure.

In this study, the effects of oxygen component in fuel and exhaust gas recirculation(EGR) method on the exhaust emissions has been investigated for a D.I. diesel engine. It was tested to estimate change of exhaust emission characteristics for the commercial diesel fuel and oxygenate blended fuel which has five kinds of blending ratio. Dimethoxy methane(DMM) contains oxygen component 42.5% in itself. and it is a kind of effective oxygenated fuel for reduction of smoke emission. It was affirmed that smoke emission was decreased with increasing of DMM blending ratio. But, NOx emission was increased compared with commercial diesel fuel. It was needed a NOx reduction countermeasure that EGR method was used as a countermeasure for NOx reduction. It was found that simultaneous reduction of smoke and NOx emission was achieved with DMM blended fuel and cooled EGR method($10{\sim}15%$).

The one of the most important subject to develop a LPDi engine is to suppress the generation of bubble inside LPG direct injector. For the purpose of this, in this study, the analogy visualization injector to visualize the generation and behavior of bubble, is manufactured and the bubbling phenomenon and behaviors are visualized and studied. The bubble inside the injector is generated at injection hole and after rising by buoyancy, it disappear around the top of a nozzle. The number of bubble generated is little changed regardless of the lapse of time but it is increased remarkably as the temperature around the injector is increased. With injection, the temperature around the injector at which the bubble is generated in_cylinder is much lower than that without injection because the transient pressure drop of fuel by injection.

The modeling of PEM (Proton Exchange Membrane) fuel cell system consisting of fuel cell stack and BOP (Balance of Plant) is presented in this paper. The effects of temperature, pressure (air, hydrogen), and humidity on the fuel cell system performance were mainly investigated using thermo-dynamical and electro-chemical equations. To understand the power distribution characteristics of fuel cell system, the effects of operating temperature and air pressure on maximum power and system power were also demonstrated. Through this study, we can get the basic insight into the fuel cell stack and BOP component sizing and it can be used effectively for the optimization of the practical fuel cell systems in purpose.

The purpose of this study is to obtain the accurate prediction for the atomization and vaporization processes of GDI spray. Atomization process is modeled using hybrid model that is composed of Linearized Instability Sheet Atomization (LISA) model and Aerodynamically Progressed TAB (APTAB) model. Vaporization process is modeled using Spalding model and Abramzon & Sirignano model. To obtain the experimental results for comparing with calculated results, the cross-sectional images of liquid and vapor phases and SMD distribution were acquired by exciplex fluorescence method and Phase Doppler Analyzer respectively. The experiment and computation was performed at the ambient pressure of 0.1 MPa, 0.5 MPa, 1.0 MPa and the ambient temperature of 293K, 473K. The calculated results by modified KIVA-II code show good agreement with experimental results.

This study was performed to investigate the behavior of vapor phase of fuel mixtures with different piston cavity diameters in a optically accessible engine. The images of vapor phases were measured in the motoring engine using exciplex fluorescence method. The conventional engine was modified as GDI engine with swirl flow. Fuel was injected into atmospheric nitrogen to prevent quenching phenomenon by oxygen. Injection pressure is 5.1MPa. Two dimensional spray fluorescence image of vapor phases was acquired to analyze spray behavior and fuel distribution inside of cylinder. Three injection timings were set at BTDC $180^{\circ}$, $60^{\circ}$and $60^{\circ}$. With a fuel injection timing of BTDC $60^{\circ}$, fuel-rich mixture was concentrated in near the cavity center. With a fuel injection timing of BTDC $60^{\circ}$, fuel-rich mixture level in the center region was highest in the S-type during the late compression stroke. With a fuel injection timing of BTDC $180^{\circ}$, fuel was not affected in a piston cavity and generally distributed as homogeneous mixture.

This paper describes spray characteristics of a swirl injector which is intended for use in a HCCI engine. Many optical diagnostics such as laser diffraction methods, and high speed camera photography are applied to measure the spray drop diameter and to investigate the spray development process. The effect of fuel properties on the spray characteristics was investigated using three different fuels because HCCI combustion is tolerant of the chemical composition of various fuels. From these results, the HCCI injector formed a hollow cone sheet spray rather than a liquid jet and the atomization efficiency is high for the low-pressure injector. The SMD of test injector was ranged from $15{\mu}m$ ${\mu}m$ We also found that the spray breakup characteristics were dependent on the fuel properties such as density, viscosity, and surface tension.

The mass flow rate of gas flow through critical nozzle depends on the nozzle supply conditions and the cross-sectional area at the nozzle throat. In order that the critical nozzle can be operated at a wide range of supply conditions, the nozzle throat diameter should be controlled to change the flow passage area. This can be achieved by means of a variable critical nozzle. In the present study, both experimental and computational works are performed to develop variable critical nozzle. A cone-cylinder with a diameter of d is inserted into conventional critical nozzle. It can move both upstream and downstream, thereby changing the cross-sectional area of the nozzle throat. Computational work using the axisymmetric, compressible Navier-Stokes equations is carried out to simulate the variable critical nozzle flow. An experiment is performed to measure the mass flow rate through variable critical nozzle. The present computational results are in close agreement with measured ones. The boundary layer displacement and momentum thickness are given as a function of Reynolds number. An empirical equation is obtained to predict the discharge coefficient of variable critical nozzle.

In steel-making process of iron and steel industry, the purity and quality of steel can be dependent on the amount of CO contained in the molten metal. Recently, the supersonic oxygen jet is being applied to the molten metal in the electric furnace and thus reduces the CO amount through the chemical reactions between the oxygen jet and molten metal, leading to a better quality of steel. In this application, the supersonic oxygen jet is limited in the distance over which the supersonic velocity is maintained. In order to get longer supersonic jet propagation into the molten metal, a supersonic coherent jet is suggested as one of the alternatives which are applicable to the electric furnace system. It has a flame around the conventional supersonic jet and thus the entrainment effect of the surrounding gas into the supersonic jet is reduced, leading to a longer propagation of the supersonic jet. In this regard, gasdynamics mechanism about why the combustion phenomenon surrounding the supersonic jet causes the jet core length to be longer is not yet clarified. The present study investigates the major characteristics of the supersonic coherent jet, compared with the conventional supersonic jet. A computational study is carried out to solve the compressible, axisymmetric Navier-Stokes equations. The computational results of the supersonic coherent jet are compared with the conventional supersonic jets.

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 of 2.0, and the operation pressure ratio is varied to obtain the different flow features in the nozzle flow. The injection flow rate is varied by means of the injection port pressure. Test conditions are in the range of the nozzle pressure ratio from 3.0 to 8.26 and the injection pressure ratio from 0 to 1.0. The present computational results show that, for a given nozzle pressure ratio, an increase of the injection pressure ratio produces increased thrust vector angle, but decreases the thrust efficiency.

Spiral jet is characterized by a wide region of the free vortex flow with a steep axial velocity gradient, while swirl jet is largely governed by the forced vortex flow and has a very low axial velocity at the jet axis. However, detailed generation mechanism of spiral flow components is not well understood, although the spiral jet is extensively applied in a variety of industrial field. In general, it is known that spiral jet is generated by the radial flow injection through an annular slit which is installed at the inlet of convergent nozzle. The objective of the present study is to understand the flow characteristics of the spiral jet, using a computational method. A finite volume scheme is used to solve 3-dimensional Navier-Stokes equations with RNG ${\kappa}-{\varepsilon}$ turbulent model. The computational results are validated by the previous experimental data. It is found that the spiral jet is generated by coanda effect at the inlet of the convergent nozzle and its fundamental features are dependent the pressure ratio of the radial flow through the annular slit and the coanda wall curvature.

Experiments are performed to investigate the detailed structure of underexpanded twin jet impinging on a perpendicular flat plate. The major parameters, such as nozzle operating pressure and nozzle spacing, are varied to create different jet flow fields resulted from the complicated interactions of the twin jets. From the surface pressure measurements and shadowgraphs taken by schlieren optical system, the jet structure is strongly dependent on the nozzle operation pressure and the spacing. The results obtained show that the closer nozzle spacing may induce to decrease the diameter of the Mach disk within the first shock cell in the underexpanded twin jet. With the increasing nozzle operating pressure and decreasing the nozzle spacing, a new shock wave appears at the entrainment region between the two jets, due to the enhancement of mixing effect of the both jets. The closer nozzle spacing makes the overall impinging pressure level higher, while severe pressure oscillation along the axis of symmetry. Furthermore it is recommended the wider spacing to obtain higher thrust under the present experimental conditions.

A computational study is performed to clarify the characteristics of supersonic moist air jet issuing from a simple sonic nozzle. The effects of the initial supersaturation on the Mach disk diameter and location, the barrel shock wave and jet boundary structures are investigated in details. The axisymmetric, compressible, Navier-Stokes equations, coupled with droplet growth equation, are solved using a third-order MUSCL type TVD finite-difference scheme. It is found that the Mach disk diameter increases with an increase in relative humidity of moist air. while its location is not significantly dependent on the relative humidity. As the relative humidity increases, the barrel shock wave and jet boundary are more expanded due to the local static pressure rise of nonequilibrium condensation.

Supersonic jet issuing from a nozzle invariably cause high-frequency noises. These consist of three principal components ; the turbulent mixing noise, the broadband shock-associated noise, and the screech tone. In present study, it was experimentally investigated to the effect of nozzle lip thickness on the characteristics of supersonic jet noise. The convergent-divergent nozzle of a design Mach number 2.0 was used in experiment. With three different nozzle-lip thicknesses, the jet pressure ratio was varied in the range between 2.0 and 12.0. Acoustic measurements were conducted by microphones in an anechoic room, and the major structures of the supersonic jets were visualized by a Schlieren optical system to investigate the effect of nozzle lip thickness. The measured results show that the characteristics of supersonic jet noise, such as overall sound pressure level (OASPL) and screech frequency, strongly depend upon the thickness of nozzle-lip.

A cone cylinder is used to obtain variable operation conditions of a sonic ejector-diffuser system. The cone cylinder is movable to change the ejector area ratio, thus obtaining variable mass flow rates. The present study investigates the effects of ejector throat area ratio and operating pressure ratio on the entrainment of secondary stream. The numerical simulations are based on a fully implicit finite volume scheme of the compressible, Reynolds-Averaged, Navier-Stokes equations. The ejector throat area is varied between 3.94 and 8.05, and the operating pressure ratio is changed from 3.0 to 9.0. The results show that the entrainment ratio and mass flux ratio become more dependent on the ejector throat area ratio, when the pressure operating ratio is low. The total pressure losses produced in the present ejector system increase with the operating pressure ratio and the ejector area ratio, but for a given operating pressure ratio, the losses are not significantly dependent on the ejector area ratio when it is larger than about 5.0.

Gas flow through orifice is encountered in many diverse fields of engineering applications. In order to investigate the critical gas flow through an orifice system, a computational analysis is performed using axisymmetric, compressible, Navier-Stokes equations which are numerically solved by a fully implicit finite volume method. In the present study, the discharge coefficients of two different types of orifices which are a straight-bore orifice and a sharp-edged orifice, are predicted to obtain the critical flow conditions. The present CFD data are compared with the previous experimental results. The present computational results show that the critical mass flow rate through orifice is well predicted and it is a strong function of Reynolds number. The discharge coefficient increases with the orifice diameter.

Perforated wall has long been employed to control a variety of flow phenomena. It has been, in general, characterized by a porosity of the perforated wall. However, this porosity value does not take account of the number and detailed shape of porous holes, but is defined by only the ratio of the perforated area to total wall surface area. In order to quantify the porous wall effects on the flow control performance, an effective porosity should be known with the detailed flow properties inside the porous holes. In the present study, a theoretical analysis using a small disturbance method is performed to investigate detailed flow information through porous hole and a computational work is also carried out using the two-dimensional, compressible Navier-Stokes equations. Both the results are compared with existing experimental data. The gasdynamical porosity is defined to elucidate the effect of perforated wall.

When a shock wave arrives at an open end of tube, an impulse wave is discharged from the tube exit and causes serious noise and vibration problems. In the current study, the effect of the cross-sectional area of tube on the impulse wave is numerically investigated using a CFD method. The Harten-Yee's total variation diminishing(TVD) scheme is used to solve the axisymmetric, two-dimensional, unsteady, compressible Euler equations. With three different cross-sectional areas of tube, the Mach number of the incident shock wave $M_{s}$ is varied between 1.01 and 1.5. The results obtained show that the directivity and magnitude of impulse wave strongly depend upon the Mach number of incident shock wave and are influenced by the tube area. It is also known that the tube cross-sectional area significantly affects the magnitude of impulse wave at or near the tube axis.

When a shock wave arrives at an open end of tube, an impulse wave is discharged from the tube exit and complicated flow is formed near tube exit. The flow field is influenced by the cross-sectional geometry of tube exit, such as circular, square, rectangular, trapezoid and etc. In the current study, three-dimensional propagation characteristics of impulse wave discharged from the tube exit with non-circular cross section are numerically investigated using a CFD method. Total variation diminishing (TVD) scheme is used to solve the three-dimensional, unsteady, compressible Euler equations. Computations are performed for the Mach numbers of the incident shock wave $M_{s}$ below 1.5. The results obtained show that the peak pressure of the impulse wave and propagation directivity depends on the cross-sectional geometry of tube exit and the Mach number of incident shock wave.

Flight bodies are subject to highly unstable and severe flow conditions during taking-off and landing periods. In this situation, the flight bodies essentially experience accelerating or decelerating flows, and the aerodynamic characteristics can be completely different from those of steady flows. In the present study, unsteady aerodynamic characteristics of an aerofoil accelerating at subsonic speeds are investigated using a computational method. Two-dimensional, unsteady, compressible Navier-Stokes simulations are conducted with a one-equation turbulence model, Spalart-Allmaras, and a fully implicit finite volume scheme. An acceleration factor is defined to specify the unsteady aerodynamics of the aerofoil. The results show that the acceleration of the subsonic aerofoil generally leads to a variation in aerodynamic characteristics and it is more significant at angles of attack.

A computational study is performed to clarify the characteristics of supersonic moist air jet issuing from a simple sonic nozzle. The effects of the initial supersaturation on the Mach disk diameter and location, the barrel shock wave and jet boundary structures are investigated in details. The axisymmetric, compressible, Navier-Stokes equations, coupled with droplet growth equation, are solved using a third-order MUSCL type TVD finite-difference scheme. It is found that the Mach disk diameter increases with an increase in relative humidity of moist air. while its location is not significantly dependent on the relative humidity. As the relative humidity increases, the barrel shock wave and jet boundary are more expanded due to the local static pressure rise of nonequilibrium condensation.

When a safety valve equipped in a nuclear power plant opens in an instant by an accident, a moving shock wave propagates downstream the valve, inducing a complicated unsteady flow field. The moving shock wave may exert severe load to the structure. So, to reduce the load acting on the wall of POSRV, a gradual opening of POSRV is adopted in general. In theses connections, a numerical work is performed to investigate the effect of valve opening time on the unsteady flow fields downstream of the valve. Compressible, two-dimensional Navier-Stokes equations are used with the finite volume method. The obtained results show that sharp pressure rise through moving shock tor the case of instant opening is attenuated by employing the gradual opening of valve. It is turned that the flows for the two cases of gradual valve opening time show the similar to that of highly under-expanded one in jet structure having expansion and compression waves and Mach stem. Also, comparing with the results for the two cases of opening time, the shorter the valve opening is, the pressure gradient at the downstream of the valve becomes softly.

The present study addresses an analytical investigation to understand the characteristics of gas flow in the High-Velocity Oxy-Fuel(HVOF) thermal spray gun. One-dimensional analysis is extended to involve the effects of the wall friction and powder particle diameter. From the present analysis it is well known that the flow characteristics inside and outside the thermal spray gun is varied depending on the combustion chamber pressure. The thermal spray gun flow is characterized by six different patterns. The powder particle size and wall friction significantly influence the powder particle velocity. The particle velocity decreases with an increase in the powder particle size. This implies that the combustion chamber pressure should be increased to achieve a higher velocity of the powder particle.

The present study addresses experimental results to investigate the details of the near field flow structures produced in the under-expanded, dual, coaxial, swirling, jet. The sonic/supersonic swirling jets are emitted from the sonic inner nozzle and the outer annular nozzle produce the co-swirling and counter swirling against the primary swirling jet, respectively. The interactions between both the secondary annular swirling and primary inner supersonic swirling jets are quantified by the pitot impact and static pressure measurements and visualized by using the Schliern optical method. The experiment is performed for different swirl intensity and pressure ratio. The results obtained show that the secondary co-swirling jet significantly changes the inner under-expanded swirling jet, such as the recirculation zone, pressure distribution, through strong interactions between both the swirling jets and the effects of the secondary counter-swirling jet is similar to the secondary co-swirl jet case.

The present study addresses experimental results to investigate the effect of the jet supply chamber configuration on the sonic/supersonic swirling jets, as the case study. The experiment is carried out using the convergent nozzle with a various different chamber configurations upstream the nozzle throat, which is composed of four tangential inlet holes for the swirling flows. The jet pressure ratio is varied between 3.0 and 7.0. The sonic/supersonic swirling jet flows are specified by the pitot impact and static pressure measurements and visualized using the Shadowgraph method. The results show that the major structures of the sonic/supersonic swirling jet are strongly influenced by the jet supply chamber.

Thrust vector control using a coflow-counterflow concept is achieved by suction and blowing through a slot adjacent to a primary jet which is shrouded by a suction collar. In the present study, the flow characteristics of thrust vectoring is investigated using a numerical method. The nozzle has a design Mach number of 2.0, and the operation pressure ratio is varied to obtain various flow features of the nozzle flow. Test conditions are in the range of the nozzle pressure ratio from 6.0 to 10.0, and a suction pressure from 90kPa to 35kPa. Two-dimensional, compressible Navier-Stokes computations are conducted with RNG ${\kappa}-{\varepsilon}$ turbulence model. The computational results provide an understanding of the detailed physics of the thrust vectoring process. It is found that an increase in the nozzle pressure ratio leads to increased thrust efficiency but reduces the thrust vector angle.

A CFD analysis has been performed to investigate turbulent heat transfer in a triangular rod bundle with a pitch-to-diameter ratio(P/D) of 1.06. Anisotropic turbulence models predicted the turbulence-driven secondary flow in the triangular subchannel and the distributions of time mean velocity and temperature showing significantly improved agreement with the measurements over the linear standard ${\kappa}-{\varepsilon}$. The anisotropic turbulence models predicted turbulence structure in large flow region fairly well but could not predict the very high turbulent intensity of azimuthal velocity observed in narrow flow region(gap).

Turbulent temperature field in a channel with wall injection has been investigated using dynamic mixed model(DMM). This flow is pertinent to internal flows inside the hybrid rocket motors. In general, the results obtained with DMM are in better agreement with DNS results compared to those of dynamic Smagorinsky model(DSM). Such favorable features of DMM are attributed to the fact that it explicitly calculates the modified Leonard stress term which takes care of the local interaction between resolved and SGS stresses and only models the remaining cross and SGS Reynolds stress terms.

This experiment has been carried out to find the structure of turbulent boundary layer with instantaneous velocity fields obtained in stream-wall-normal planes using a stereo-PIV (Particle Image Velocimetry) method. And it has been measured perpendicular plane and horizontal plane with hairpin vortex structure by Reynolds number change and made third dimension shape for section of horizontal plane through stereo-PIV. In the outer layer hairpin vortices occur in streamwise-aligned packets that propagate with small velocity dispersion. A streaky structure is composed of counter-rotating vortex. According as y+ increases, streaky structure's interval space decrease, and it shows that hairpin shape of prior research is vertified. The objective of the present research is to gain a better understanding of coherent structures in the outer of wall turbulence by experimentally examining coherent structures.

Fluid-elastic instability is believed to be a cause of the large-amplitude vibration and resulting rapid wear of heat exchanger tubes when the flow velocity exceeds a critical value. For sub-critical flow velocities, the random turbulence excitation is the main mechanism to be considered in predicting the long-term wear of steam generator tubes. Since flow-induced interactions of the tubes with tube supports in the sub-critical flow velocity can cause a localized tube wear, tube movement in the clearance between the tube and tube support as well as the normal contact force on the tubes by fluid should be maintained as low as possible. A simplified method is used for predicting fretting-wear damage of the double $90^{\circ}$U-bend tubes. The approach employed is based on the straight single-span tube analytical model proposed by Connors, the linear structural dynamic theory of Appendix N-1300 to ASME Section III and the Archard's equation for adhesive wear. Results from the presented method show a similar trend compared with the field data. This method can be utilized to predict the fretting-wear of the double $90^{\circ}$U-bend tubes in steam generators.

A numerical study is made of a flow in a cylinder with a rotating grooved endwall disk. The aim is to describe differences in the flow fields when there is concentrically-grooved obstacle characterized by amplitude(a) and wave number(N). The Reynolds number(Re) is varied from $10^{3}$ to $10^{4}$ and the aspect ratio(Ar) fixed to 1.0 for the most part of the simulation. For the various cases of amplitude(a) and wave number(N), numerical results are acquired. As the endwall groove roughness increases until certain limit, the interior azimuthal velocity component(v) increases drastically. But over the limit, the swirl motion chararcterized by velocity v decreases and finally it approaches much alike Ar=1.0-a case. The reason of activating swirl motion is based on increasing of torque transported by endwall disk. Torque coefficients($C_{T}$) are aquired for the various (a,N,Re) combinations and the limiting phenomena of swirl motion activation is explained.

An analysis on the steady-state has been made of flow of a compressible fluid rapidly-rotating in a pipe. The flow is induced by an small arbitrary azimuthally-varying thermal forcing added on the basic state of rigid body isothermal rotation. The system Ekman number is assumed to be very small value. Analytic solutions have been obtained for axisymmetric and non-axisymmetric types, in which the axisymmetric solution comes from the azimuthally-averaged wall boundary condition and the non-axisymmetric solution from fluctuating wall boundary condition.

It has been studied the flow near a rotating disk with surface topography. The system Ekman number is assumed very small, i.e., $E[{\equiv}\frac{\nu}{{\Omega}^{\ast}L^{\ast2}}]<<1$ in which $L^{\ast}$ denotes a disk radius, ${\nu}$ kinematic viscosity of the fluid and ${\Omega}^{\ast}$ angular velocity of the basic state. Disk surface has a sinusoidal topographic variation along radial coordinate, i.e., $z={\delta}cos(2{\pi}{\omega}r)$, where ${\delta}$ and ${\omega}$ are, respectively, nondimensional amplitude and wave number of the disk surface. Analytic solutions, being useful over the parametric ranges of ${\delta}{\sim}O$( $E^{1/2}$ ) and ${\omega}{\leq}O$ ( $E^{1/2}$ ), are secured in a series-function form of Fourier-Bessel type. An asymptotic behavior, when $E{\rightarrow}0$, is clarified as : for a disk with surface roughness, in contrast to the case of a flat disk, the azimuthal velocity increases in magnitude, together with the thickening boundary layer. The radial velocity, however, decreases in magnitude as the amplitude of surface waviness increases. Consequently, the overall Ekman pumping at the edge of the boundary layer remains unchanged, maintaining the constant value equal to that of the flat disk.

In this paper, we report the numerical and experimental solutions of the vortical flows driven by an impingement of fluid from the bottom wall of a circular cylinder. We managed to visualize successfully the flow pattern shown on the vertical plane through the container axis. The numerical results are shown to compare well with the experimental results for the case of infinity Rossby number. The satisfactory agreement between the two results was possible when in the numerics the free surface was treated as a solid wall so that a no-slip condition was applied on the surface. The numerical solutions reveal that inertial oscillation plays an important role at small Rossby numbers, or at a large background rotation.

Quantitative performance test on the conventional 2D-PIV and the hybrid angular 3D-PIV (Stereoscopic PIV) was carried out. LES Data sets on an impinging jet which are provided on the webpage(http://www.vsj.or.jp/piv) for the PIV Standard Project were used for the generation of virtual images. The generated virtual images were used for the 2D-PIV and 3D-PIV measurements. The measurement results showed that the results obtained by 2D-PIV on average values are closer to the LES data than those obtained by 3D-PIV, but the turbulent properties obtained by 2D-PIV are largely underestimated than those obtained by 3D-PIV.

A wind-tunnel experiment is carried out to examine the applicability of a new passive device, wake disrupter, to flow over a model vehicle for drag reduction. The wake disrupter is a small-size rectangular body attached to a part of the trailing edge of the model vehicle, designed to perturb an essentially two-dimensional nature of wake. A pair of wake disrupter is mounted on the mid-span at the upper and lower trailing edges. From a parametric study about the size of wake disrupter, it is found that the optimum disrupter increases the base pressure by about 20%. Large eddy simulation is also conducted to confirm the experimental result, and shows that the wake is indeed disrupted by the present device.

A process of 3-D stereo particle image velocimetry(PIV)was developed for the measurement of an illuminated sliced section field of 3-D complex flows. The present method includes modeling of camera by a calibrator based on the homogeneous coordinate system, transformation of the oblique-angled image to the right-angled image, identification of 2-D velocity vectors by 2-D cross-correlation equation, stereo matching of 2-D velocity vectors of two cameras, accurate calculation of 3-D velocity vectors by homogeneous coordinate system, removal of error vectors by a statistical method followed by a continuity equation criteria, and finally 3-D display as the post processing. An experimental system was also used for the application of the proposed method. Two high speed digital CCD cameras and an Argon-Ion Laser for the illumination were adopted to clarify the time-dependent characteristics of the leading edge extension(LEX) in a highly swept shape applied to a delta wing found in modern air-fighters.

Recently, the recording density of hard disk drives has improved at an annual percentage rate of 100%. Therefore for faster access, higher disk rotational speeds will be required. The influence of the airflow produced by the rotation of a disk on the positioning accuracy has become a serious topic of research and the aerodynamic aspect of hard disk drives is now quite considerable with the increases in recording density and higher rotational speeds. Unsteady airflow in an actual hard disk drive is numerically simulated by using LES(Large Eddy Simulation) technique, we could predicted and aerodynamic mechanism that was related actuators' surroundings in HDD. At a result, with modifying the various shapes of the E-block and Damper, we estimated the characteristic of the influence of airflow in HDDs.

With the increase in recording density and higher rotation speeds, the aerodynamic aspect of hard disk drives (HDDs) is now quite significant. To achieve high TPI (tracks per inch) for hard disk drives, the actuator must induce less air turbulence. Also, alternative cross sections for the actuator arms are proposed for the reduction of air turbulence. In this study, off-track vibration of actuator arms used in of 3.5" hard disk drives with different cross sections have been measured by using Laser Doppler Vibrometer (LDV). It has been found that the vibration of the actuator arms with modified cross sections is lower than the conventional one.

In this paper, we investigated the disc warping in high-speed slim-type optical disc drive. Recently, the information storage devices are increasing track density and higher rotation speed to enhance their recording capacity and their data transfer rate. Generally, ODD used in the Lap-top Computer has small inner space. So, the flow instability of inner space is rapidly increased as its higher rotation speed. An extreme case, the flow instability causes the malfunction of readout at pick-up of drives. The experiments and numerical analysis were carried out for several cases, the result shows the influence of airflow to the disc warping.

The variations of temperature and velocity fields in a Hele-Shaw convection cell (HSC) were investigated using a holographic interferometry and 2-D PIV system with varying Rayleigh number. To measure quasisteady changes of temperature field, two different measurement methods of holographic interferometry; double-exposure method and real-time method, were employed. In the double-exposure method, unwanted waves can be eliminated effectively using digital image processing technique and the reconstruction images are clear, but transient flow structure cannot be reconstructed clearly. On the other hand, transient convective flow can be reconstructed well using the real-time method. However, the fringe patterns reconstructed by the real-time method contain more noises, compared with the double-exposure method. Experimental results show a steady flow pattern at low Rayleigh numbers and a time-dependent periodic flow structure at high Rayleigh numbers. The periodic flow pattern at high Rayleigh numbers obtained by the real-time holographic interferometer method is in a good agreement with the PIV results.

The PSP technique has been used to measure pressure distribution on a model surface quantitatively. The objective of this study is to develop PSP technique which can be applied to low-speed flows. Four different PSP formulations including two porphyrins (PtOEP and PtTFPP) and two polymers (Poly(TMSP) and RTV-118) were tested and the performance of each combination was checked. In the static calibration, the luminescent intensity of the PSP coatings was measured from 0kPa to 11kPa with 0.5, 1, 2kPa increments. Among 4 PSP formulations tested, the combination of PtOEP and RTV-118 shows the best performance. The PSP technique was applied to an oblique impinging jet to measure the pressure field distribution on the impinging plate.

Simultaneous measurement with PLIF(Planar Laser-Induced Fluorescence) and Stereo-PIV(Stereoscopic Particle Image Velocimetry) was performed to investigate the structural characteristics of flow field in Rushton Turbine Mixer. Instantaneous 3D velocity fields are measured by two 2K ${\times}$ 2K CCD cameras focused on an object plane with the angular displacement methods while the concentration fields are obtained through the measurement of the fluorescence intensity of Rhodamine B tracer excited by the second pulse of Nd:Yag laser light. Image distortion due to the camera view-angle is compensated by a mapping function. Finally, the spatial structures of turbulent mixing around Rushton turbine were identified by the calculation of cross-correlation fields between the velocity and concentration field.

The general cooling tower is a device for making a cooling water in refrigerant condensers or industrial process heat exchangers. The present cooling tower have defects with noises, complicated structure and environmental problems. In this paper, we constituted a new water cooling system by using a evaporating latent heat in an enclosed tank, and this system is consisted of an enclosed vacuum tank and water driving ejector system. Several experimental cases were carried out for improvement methods of high vacuum pressure and water cooling characteristics. The ejector performance was tested in case of water temperature variations that flows in the ejector. Based on the vacuum pressure by water driving ejector, the water cooling characteristics were investigated for the vaporized air condensing effects.

A numerical method is presented for computing unsteady incompressible two-phase flows with immersed solids. The method is based on a level set technique for capturing the phase interface, which is modified to satisfy a contact angle condition at the solid-fluid interface as well as to achieve mass conservation during the whole calculation procedure. The modified level set method is applied for numerical simulation of bubble deformation in a micro channel with a cylindrical solid block and liquid jet from a micro nozzle.

In the present study, the gas injection system based on air-water model was designed to investigate the behavior characteristics of bubbles injected into a ladle. The parameters such as gas volume fraction and bubble rise velocity were exprementally measured in a gas-liquid flow region. To measure gas volume fraction, an electo-conductivity probe was used and bubble rise velocity was obtained by a high speed CCD camera. Gas volume fraction was symmetric to the axis of nozzle secured on the bottom of a ladle. The bubble rise velocity was calculated for two different experimental conditions. That is, gas flow conditions were following two case: 1) Q = $0.63{\times}10^{-4}$ $m^{3}/s$, 2) $1.26{\times}10^{-4}$ $m^{3}/s$. As a gas injected into the liquid ladle, the liquid-phase region is circulated by bubbles' behavior. The bubble rise velocity was influenced of the circulation flow of liquid phase. As a result, the bubble rise velocity was appeared higher middle region of ladle than near the nozzle.

The aspect ratio is the main parameter which governs the outer flow pattern and nozzle performance. And in this study, some flow characteristics with the variation of nozzle aspect ratios such as mean pressure distributions along the center line of the outer flow, flow coefficients and the diffusion angles have been experimentally investigated. Through the experimental analysis, the higher aspect ratio was known to decrease the jet kinetic energy because of the friction losses at the outer of nozzle. As the result, it is found that the nozzle performance depends mainly on the aspect ratio of nozzle.

The Ejector is used to get low pressure, and it has been applied to a lot of industry field like the heat engine, the fluid instrument power plant, the food industry, environment industry etc... because there are not any problem even it is mixed with a any kind of liquid, gas, and solid. The flow characteristics in ejector are investigated by PIV and CFD. The experiment using PIV measurement for mixing pipe's flow characteristics acquired velocity distribution, kinetic energy distribution, and whirlpool . (Condition : when mixing pipe's diameter ratio is 1:1.9, and the flux is $Q_{1}=1.136{\imath}/s$, $Q_{2}=1.706{\imath}/s$, $Q_{3}=2.276{\imath}/s$. Based on the PIV and the CFD results, the flow characteristics in ejector are discussed, and it shows the validity of this study.

This paper reports the characteristics of the three dimensional turbulent flow in the rectangular-sectioned 180 degree bends by Hot-wire anemometer. Grande and Kool proposed a cooling law for the measurements of the flow through the narrow passage. The authors noticed that the calibration coefficients of original method are not constant and fairly sensitive to the flow approaching angle. Measured voltages are converted to three velocity and six Reynolds stress components using the modified method in which the coefficients are treated as a function of approaching angle.

In this study, a numerical optimization to find the optimal shape of streamwise periodic ribs mounted on both of the principal walls is performed to enhance turbulent heat transfer in a rectangular channel. The optimization is based on Navier-Stokes analysis of flow and heat transfer with $k-{\varepsilon}$ turbulence model and is implemented using response surface method. The width-to-height ratio of a rib, rib height-to-channel height ratio, rib pitch to rib height ratio and distance between opposite ribs to rib height ratio are chosen as design variables. The object function is defined as a function of heat transfer coefficient and friction drag coefficient with weighting factor. Optimum shapes of the rib have been investigated for the range of 0.0 to 0.1 of weighting factor.

The main objective of present study is to obtain information for flow characteristics, such as velocity and wall shear stress, of bifurcation in blood vessel. Branch flows for Newtonian fluids are simulated by using Fluent V.6.0. The numerical simulations are carried out for five cases divided by different values of bifurcation angle and area ratio. As a result of simulation, high wall shear stress is appeared at the bifurcated region. As increasing bifurcation angle, pressure drop is increasing. In addition, as the area is decreasing, pressure drop and wall shear stress is increasing.

Many clinical studies have suggested that the blood flow in ideal geometry is involved in the development of atherosclerosis. This study simulated blood flow in the abdominal artery with real geometry to investigate MWSS(mean wall shear stress), AWSS(amplitude of wall shear stress) and OSI(oscillator shear index). The calculation grid for the real geometry was constructed by extracting the surface of arterial wall from CT(Computed Tomography) or MRI(Magnetic Resonance Imaging) sheets called as DICOM (Digital Imaging and Communications in Medicines). The calculated MWSS, AWSS and OSI are much different from those of ideal geometry calculation. The MWSS increased while the AWSS decreased. Many shear forces are related to shapes of gradient. This paper will give clinical datum where the MWSS, AWSS and OSI are strong or weak. The hemodynamic analysis based on real geometry can provide surgeons with more reliable information about the effect of blood flow.

The purpose of the present study was to establish the mechanism of the generation of atherosclerosis by analyzing the hemodynamic variables in the coronary artery where atherosclerosis occurs frequently. From the previous results, the stenosis phenomena due to atherosclerosis were related to not only biochemical reaction between blood and blood vessel but also the hemodynamic factors like flow separation and oscillatory wall shear stress. The present study aimed to investigate the causes of the generation and progression of atherosclerosis in the coronary artery. This study also aimed to develop the softwares which generate automatically three dimensional vascular models obtained by the angiogram images and the computer vision techniques. In the present study, the flow patterns for full three-dimensional hemodynamic characteristics were analyzed. To understand the three-dimensional hemodynamic characteristics, the wall shear stress distributions and secondary flows were investigated quantitatively.

Numerical simulations are conducted to investigate the mechanism of hovering flight by single flapping wing, and to examine the effect of the phase difference between the fore- and hindwings in hovering flight by two flapping wings. The numerical method used is based on an immersed boundary method in Cartesian coordinates. The Reynolds number considered is Re=150 based on the maximum translational velocity and chord length of the wing. For single flapping wing, the stroke plane angles are $0^{\circ}$, $30^{\circ}$, $60^{\circ}$, $75^{\circ}$ and $90^{\circ}$ and the downstroke angles of attack are varied for each stroke angle. Results show that for each stroke plane angle, there is an optimal angle of attack to maximize the vertical force. Below the stroke angle of $60^{\circ}$, wake capturing reduces the negative vertical force during the upstroke. For two flapping wings, The phase lags of the hindwing are $0^{\circ}$, $90^{\circ}$, $180^{\circ}$ and $270^{\circ}$. The amplitudes of the stroke are 2.5 and 4.0 times the chord length at each phase lag. The results show that maximum vertical force is generated when the phase lag is zero, and the amplitude of the vertical force is minimum at the phase lag of $180^{\circ}$.

The objective of the present study was to visualize the pulsatile flow field in a branching model by using the high-resolution PIV system. A bifurcated flow system was built for the experiments in the pulsatile flow. Harvard pulsatile pump was used to generate the pulsatile velocity waveforms. Conifer powder as the tracing particles was added to water to visualize the flow fields. Two consecutive particle images at several cross sections of the flow filed were captured by the CCD cameras ($1K{\ast}1K$ and $640{\ast}480$). The results after the image processing clearly showed the recirculation zones and the formation of the paired secondary flows from the distal to the apex in the bifurcated model. The results also indicated that the flow velocities in the inner wall moved faster than those in the outer wall due to the inertial force effects and the helical motions generated in the branch flows as the flow proceeded toward the outer wall. While the PIV images from the $1K{\ast}1K$ camera were closer to the simulation results thantheimagesfromthe640${\ast}$480camera,bothresultsofthePIVexperimentsusingthetwocamerasgenerallyagreed quitewellwiththeresultsfromthenumericalsimulation.

In this paper, the performance of a domestic Kim-Chi refrigerator is predicted by using a calculation model & experiment. The objectives of this study are to find out the best design points of a refrigeration system and calculate an adiabatic characteristic to change outdoor temperature. The best design points such as refrigerant charge and capillary length were experimentally investigated. And the calculation model is conducted as a function of calculation parameters and outdoor temperature. According to this study results, the best design points of a refrigeration system are each 95g of a refrigerant charge and 3500/3500mm of capillary lengths. And the power consumption is 13.578 Kwh/month. And a part of the worst heat loss is a front side of a domestic Kim-Chi refrigerator body.

In this study, a computer analysis has been developed for predicting the pipe pressure of the intake and exhaust manifold in a single cylinder engine. To get the boundary conditions for a numerical analysis, one dimensional and unsteady gas dynamic calculation is performed by using the MOC(Method Of Characteristic). The main numerical parameters are the variation of the exhaust pipe diameters to calculate the pulsating flow when the intake and exhaust valves are working. As the results of numerical analysis, the shapes and distributions of the exhaust pipe pressures were influenced strongly on the cylinder pressure. As the exhaust pipe diameter is decreased, the amplitude of exhaust pressure is large and the cylinder pressure was showed low in the region of intake valve opening time.

This work presents a numerical analysis of two-phase natural circulation flow in reactor cavity under external vessel cooling. Steady, incompressible, three-dimensional Reynolds-averaged Navier-Stokes equations for multiphase flows with zero equation turbulence model are solved to predict the shear key effect on the circulation rate of cooling water and the distribution of void fraction according to the different mass flow of inlet air. Results show that shear key has a positive effect on the circulation rate of cooling water and induce a local increase of void fraction below the shear key, but not remarkably.

We have studied on the main causes of pressure drop in the paint hose under the spraying with airless spray pump and suggest calculation method for selecting the optimum pressure drop during the spraying. Among many factors contributing the pressure drop during the airless spraying procedure of paint, some of the key factors and the correct method for checking pressure loss between airless pump and spray gun are addressed. We have developed pressure loss calculation method which depending on hose length and diameter, viscosity and flow rate in the painting hose during the spraying. also we have developed calculation equation for the expected spray tip flow rate which depending on pressure and specific gravity and tip size.

A helical vane is applied to reduce the magnitude of the impulse wave discharged from the exit of a duct. A shock tube with an open end is used to investigate the effect of the helical vanes on the impulse wave magnitude. Four different types of helical vanes are installed into the low-pressure tube of shock tube. The magnitude of the incident shock wave is varied below 1.25, and the magnitude of impulse wave is measured using a pressure transducer mounted on a wedge probe. Instant images of the impulse wave are obtained by means of the Schlieren optical method. The present experimental results show that the helical vane considerably reduces the magnitude of the impulse wave and the vane effects are more remarkable for stronger incident shock wave.

The objective of the present study is to investigate the changes in the acoustic source characteristics and far-field noise propagation in an incompressible round jet at Re=10000 for single-frequency excitations using large eddy simulation and Lighthill acoustic analogy. We apply excitations at a frequency corresponding to the jet-column mode ($St_{D}=0.85$) or maximum growth rate in the shear layer ( $St_{\theta}=0.017$ ). The acoustic source derived from the Lighthill acoustic analogy is the second spatial derivative of the Reynolds stresses. In the case of $St_{D}=0.85$, vortex ring and large scale structures are dominant sources, whereas in the case of $St_{\theta}=0.017$, the main sources are located at an upstream position along the shear layer than in the uncontrolled case. Also, the far-field noise propagates along the axial direction due to excitation.

In this study, 2D computations of the Aeolian tones for some obstacles (circular cylinder, square cylinder and NACA0012 airfoil) are simulated. First of all, we calculate the flow noise generated by a uniform flow around a two-dimensional circular cylinder at Re=150 are simulated by applying the finite difference lattice Boltzmann method (FDLBM). The third-order-accurate up-wind scheme (UTOPIA) is used for the spatial derivatives, and the second-order-accurate Runge-Kutta scheme is applied for the time marching. The results show that we successively capture very small acoustic pressure fluctuation with the same frequency of the Karman vortex street compared with the pressure fluctuation around a circular cylinder. The propagation velocity of the acoustic waves shows that the points of peak pressure are biased upstream due to the Doppler effect in the uniform flow. For the downstream, on the other hand, it is faster. To investigate the effect of the lattice dependence, furthermore, simulations of the Aeolian tones at the low Reynolds number radiated by a square cylinder and a NACA0012 airfoil with a blunt trailing edge at high incidence are also investigated.

Flow-induced noise propagated from flow over a sphere is numerically investigated for laminar flow at Re = 300 and 425, and for turbulent flow at Re = 3700 and $10^4$, where the Reynolds number is based on the freestream velocity and the sphere diameter. The numerical method used for obtaining the flow over a sphere is based on an immersed boundary method in a cylindrical coordinate system. The Curle’s solutions of the Lighthill’s acoustic analogy with and without the far-field and compact-source approximation are used in order to investigate the noise field from flow over a sphere. Since the drag and lift forces change irregularly in time at Re = 425, 3700 and $10^{4}$, the noise propagates in a complicated manner. At Re = 300, 425 and $10^{4}$, the noise from dipole sources is much larger than that from quadrupole sources. On the other hand, at Re = 3700, the quadrupole source becomes dominant. The temporal variation of the flow-induced noise around a sphere is obtained at some observation points, which shows that the peak frequency corresponds to the Strouhal number associated with the wake instability.

The pumping characteristics of a disk-type drag pump (DTDP) from free molecular flow region to the slip flow region are calculated by the direct simulation Monte Carlo (DSMC) method. In this study, the pumping performance is studied numerically for several channel depths. The interaction between molecules is modeled by variable hard-sphere (VHS). The no time counter method is used as a collision sampling technique. The clearance between rotor and stator is considered an effect on performance. Spiral channels are cut on both upper and lower sides of rotating disks, and stationary disks are planar. A three-dimensional DSMC method for the analysis of steady rarefied flows in a single-stage DTDP has been developed. Velocity and density fields were obtained by the DSMC simulation in the rotor. The present experimental data in the outlet pressure range of $7.5{\times}10^{-3}{\sim}4$ Torr were compared with the DSMC results in the single-stage DTDP. Comparison between the experimental data and DSMC results showed good agreement.

Experiments were done for the three dimensional unsteady flow in a counter rotating axial flow fan under stable operating condition. Flow fields in a counter rotating axial flow fan were measured at cross-sectional planes of the upstream and downstream of each rotor. Cross sectional flow patterns were investigated through the acquired data by the $45^{\circ}$ inclined hot-wire. Flow characteristics such as tip vortex, secondary flow and tip leakage flow were confirmed through axial, radial and tangential velocity vector plot. Swirl velocity, which was generated by the front rotor, was recovered in the form of static pressure rise by the rear rotor except for hub and tip regions.

This work analyzes the effects of the independent variation of different geometric dimensions of compressor valves on the effective flow and force areas using a circular valve plate, such as different geometry of the valve seat, and the valve reed is opened and closed by pressure pulsation, the flow characteristic of the refrigerant passing the valve is very important. In the present study, a circular disk with inclination is assumed to be the valve reed of a reciprocating compressor and numerical analysis of three dimensional velocity fields are performed for theradial flow through the valve model. The effective flow and force area which are required to predict the efficiency of the valve are required to predict the efficiency of the valve are measured and compared with the numerical analysis in this research.

In the present study, the performance characteristics and the number of stall cells during rotating stall of a centrifugal air compressor were experimentally investigated. If the flow coefficient is lower than 0.150, the static pressure at impeller inlet is higher than that at inlet duct And reverse flow is observed under these flow coefficient region. Maximum adiabatic efficiency is obtained for the tested compressor around flow coefficient 0.128, and it is independent of compressor rotating speed. The number of stall cells and their rotational speeds are distinctive features of the rotating stall phenomenon. The present study is mainly concerned with the number of stall cells and their rotational speeds. The interpretation method of visualization is based on the pressure distribution in the circumference pressure fields while plotting the pressure and its harmonics variations in time in polar coordinates.

A study for increase of fan flow rate by geometric modification has been conducted to decrease temperature rise of marine generator inner part. Through experiment of a real product, a performance curve for various flow resistances was obtained. Flow analyses for each cases were done by using commercial code-FLUENT and the results were very similar to experimental data (0.7% deviation at normal operating condition). Through flow analysis results for various design geometric modification, a scroll type fan was adopted as a best design geometry with 100Pa more pressure and 22% more flow rate than original fan.

Unsteady nature of a tip leakage vortex in an axial flow fan operating at a design and off-design operating conditions has been investigated by measuring the velocity fluctuation in a blade passage with a rotating hotwire probe sensor. Two hot-wire probe sensors rotating with the fan rotor were also introduced to obtain the cross-correlation coefficient between the two sensors located in the vortical flow as well as the fluctuating velocity. The results show that the vortical flow structure near the rotor tip can be clearly observed at the quasi-orthogonal planes to a tip leakage vortex. The leakage vortex is enlarged as the flow rate is decreased, thus resulting in the high blockage to main flow. The spectral peaks due to the fluctuating velocity near the rotor tip are mainly observed in the reverse flow region at higher flow rates than the peak pressure operating condition. However, no peak frequency presents near the rotor tip for near stall condition.

A cross-flow fan consists of an impeller, a stabilizer and a rearguider. When it applied for an air conditioner, an evaporator should be added. It relatively makes high dynamic pressure at low speed because a working fluid passes through an impeller blade twice and blades have a forward curved shape. Therefore, the performance of a cross-flow fan is influenced 25% by the impeller, 60% by the rearguider and the stabilizer, 15% by the heat exchanger. At the low flow rate, there are a rapid pressure head reduction, a noise increase and an unsteady flow against a stabilizer and a rearguider. Moreover, the reciprocal relation between the impeller and the flow passage is the important factor for performance improvement of the cross-flow tan because each parameter is independent. The performance characteristics in the cross-flow fan are graphically depicted with various impeller outlet angles and rearguiders.

Numerical investigations have been performed to examine the effects of the computational grids on the prediction of the flow characteristics inside the turbine cascades. Three kinds of grid system based on H-type grid are applied to the high-turning transonic turbine rotor blades and comparisons with the experimental data and the numerical results of each grid structure have been done. In addition, the grid sensitivity on the estimation of the blade performances has been investigated.

Mixers are used in various industrial fields where it is necessary to intimately mix two reactants in a short period of time. However, despite their widespread use, complex unsteady flow characteristics of industrial mixers are not systematic investigated. The present study aimed to clarify unsteady flow characteristics induced by various impellers in a tank. Impellers arc hydrofoil turbine and neo-hydrofoil turbine types. A high speed CCD camera and an Ar-Ion laser for illumination were adopted to clarify the time-dependent flow characteristics of the mixers. The rotating speed of impellers increased from 6Hz to 60Hz by 6Hz, The maximum velocity around neo-hydrofoil impeller is higher than the hydrofoil type impeller. These two types of turbine shows that typical flow characteristics of axial turbine and suitable for mixing high-viscosity materials.

The atomization characteristics were investigated through the influence of the change of GLR and the change of working fluid on droplet size distribution and mean diameter of drop produced by effervescent atomizer. For simultaneous injection of water and high viscous waste vegetable oil, effervescent atomizer with two aerator tubes was specially designed. From the experimental results, regardless of mass fraction of vegetable oil in working fluids, it is expected that effervescent atomizer will exhibit excellent atomization performance at the high GLR conditions.

Stall inception means the phenomena of rotating stall initiation. The initiation mechanism of rotating stall, the existence of stall precursor, the behavior of stall precursor, stall warning scheme and control scheme are the main interests in stall inception research. Compared to the studies on rotating stall which has long history, the stall inception has been studied for about recent 20 years. After the first discovery of stall precursor in about 20 years ago, many studies were reported on stall inception phenomena. The inception pattern of "mode" and "spike" were found, and some of its characteristics are known. And now the stall inception has become one of the fascinating fields in turbomachinery. This paper reviews the results and analysis methods on stall inception studies.

In spite of its simple design, structure and operating mechanism, swing check valves are one of the critical components which adversely affect the safety of the nuclear power plants if they fail to function properly. Therefore, it is important to evaluate the performance condition of the swing check valves in safety-related systems. The performance characteristics of swing check valves include opening characteristics, the minimum required flow velocity, the pressure drop at design flow, the disc stability, and the effect of the upstream disturbances. Among factors to identify the performance of a swing check valve, a method to evaluate the opening characteristics and the minimum required flow velocity, which guarantees to fully open the disc and hold the disc without motion, are presented to determine the operating region of the swing check valve, such as stable, tapping, or oscillation. Based on the determined operating region and opening characteristics, the simple methods of wear and fatigue analyses of the specific parts of the valve are also described.

A front-end cooling fan is designed and tested in the present study. The design technique is developed using the one-dimensional inviscid flow through the fan blade, the empirical equations, and the experimental correlations. Design data for the blade can be obtained for a given flow rate and a pressure rise. A parabolic function is used to generate a sweep of the fan. Characteristics of the blade geometry are discussed between the huh and the tip. The fan is tested in the fan test unit. The measured volume flow rate at the operating point is in good agreement with that of the design specifications. Sound pressure levels of the noise are predicted with the Ffowcs Williams-Hawkings equations. Calculation results of the sound pressure level(SPL) 1m away from the fan are obtained and cpmpared with the measured data.

An experiment was conducted to characterize water droplet charging performance of an electrostatic spraying nozzle for an electrostatic wet scrubber. Charge-to-mass ratios, the nozzle currents divided by the mass flow rate of water were obtained with respect to the applied voltage to the ring-electrode for 2 different flow conditions. It was shown that the charge-to-mass ratio increased in proportion to the applied voltage and tended to saturate at a certain higher voltage.

A new model for dense gas dispersion is formulated within the Lagrangian framework. In several accidental released situations, denser-than-air vapour clouds are formed which exhibit dispersion behavior markedly different from that observed for passive atmospheric pollutants. For relevant prediction of dense gas dispersion, the gravity and entrainment effects need to implemented. The model deals with negative buoyancy which is affected by gravity. Also, the model is subjected to entrainment. The mean downward motion of each particle was accounted for by considering the Langevin equation with buoyancy correction term.

This paper describes a shower type ultrasonic cleaning system using particle acceleration, developed for better effectiveness in washing out agricultural chemicals from fruits or vegetables, compared to the existing bath type ultrasonic cleaning system. The shower type system consists of a signal generator, a transducer, and a spouting nozzle, mainly. The system has been tested based on the method recommended by Korea Food and Drug Administration. It is found that shower type system shows better performance and cost-effectiveness than the bath type system that is commercially available .

A Navier-Stokes equation solver for incompressible viscous flows with free surface is developed and tested. This is based upon a fractional time step method and a non-staggered finite volume formulation for unstructured meshes. For time advancement scheme, Adams -Bashforth method for convective term and Crank-Nicolson method for diffusive term are applied. The interface between two fluids with different fluid properties is tracked with Piecewise Linear Interface Calculation(PLIC) Volume-of-Fluid(VOF) methods. Computational results are presented for some test problems: the broken dam, the sloshing in a rectangular tank, the filling of a cylindrical tank.

In most industrial applications, the geometrical complexity is combined with the moving boundaries. These problems considerably increase the computational difficulties since they require, respectively, regeneration and deformation of the grid. As a result, engineering flow simulation is restricted. In order to solve this kind of problems the immersed boundary method was developed. In this study, the immersed boundary method is applied to the numerical simulation of stationary, rotating and oscillating cylinders in the 2-dimensional square cavity. No-slip velocity boundary conditions are given by imposing feedback forcing term to the momentum equation. Besides, this technique is used with a second-order accurate interpolation scheme in order to improve the accuracy of flow near the immersed boundaries. The governing equations for the mass and momentum using the immersed boundary method are discretized on the non-staggered grid by using the finite volume method(FVM). This study presents the possibility of the immersed boundary method to apply to the complex flow experienced in the industrial applications.

In this paper, slurry fluid motion, abrasive particle motion, and roles of groove patterns on the pads are numerically investigated in the 2D and 3D geometries. The simulation results are analyzed in terms of experimental removal rate and WIWNU (within wafer non-uniformity) for ILD (inter level dielectric) CMP process. Numerical investigations reveal that the grooves in the pad behave as uniform distributor of abrasive particles and enhance the removal rate by increasing shear stress. Higher removal rate and desirable uniformity are numerically and experimentally observed at the pad with grooves. Numerical analysis is very well matched with experimental results and helpful for understanding polishing mechanism and local physics.

For simulation of three-dimensional turbulent flow with LES and DNS takes much time and expense with current available computing resources. It is nearly impossible to simulate turbulent flow with high Reynolds number. So, the emerging alternative is the Grid computing for needed computation power and working environment. In this study, the CFD code was parallelized to adapt it for the parallel computing under the Grid environment. In the first place, the Grid environment was built to connect the PC-Cluster facilities belong to the different institutions using communication network system. And CFD applications were calculated to check the performance of the parallel code developed for the Grid environment. Although it is a fundamental study, it brings about a important meaning as first step in research of the Grid.

Flow analyses using FLUENT 5.4 code were performed to validate the application of the local bidirectional flow tube in case of water and air flow. In this study, sensitivity studies have been performed to optimize the design features of the bi-directional flow tube which can be applied for the various experimental conditions. 2-dimensional axisymmetric steady state flow analyses have been performed. By calculating the differential pressure at both the front and the rear hole of the flow tube, K values were evaluated. The K values show good linearity regardless of the tube sizes and the Re numbers in both water and air flow. And system pressure and water subcooling didn’'t affect the K values. Under the elevated pressure of 80bar with 80K water subcooling, the K value indicates a similar trend with the case of 2bar with 80K water subcooling.

The analysis involves an adaptive grid that is created under a criterion of element categorization of filling states and locations in the total region at each time step. By using an adaptive grid wherein the elements, finer than those in internal and external regions, are distributed at the surface region through refinement and coarsening procedures, a more efficient analysis of transient fluid flow with free surface is achieved. Using the proposed numerical technique, the collapse of a dam is analyzed. The numerical results agree well with the theoretical solutions as well as with the experimental results. Through comparisons with the numerical results of several cases using different types of grids, the efficiency of the proposed technique is verified.

In this paper, we present that now and heat transfer characteristics of oscillatory impinging jet which have been numerically investigated using parallel computations. Numerical value were obtained for dimensionless distance H=4, dimensionless outlet length L=45 and Reynolds number Re= 1500. It was found that the oscillatory impinging jet generated the regular heat transfer area even though the maximum heat transfer is lower than nonoscillatory impinging jet. We also found that heat transfer depends on the period of nozzle for the oscillatory impinging jet.

The fluid propulsion mechanism of two oscillating flat plates is studied numerically using a discrete vortex method. Presently, the flat plates are assumed to be rigid. To analyze the closely coupled aerodynamic interference between the flat plates, a core addition scheme and a vortex core model are combined together. A calculated wake pattern for a flat plate in heaving oscillation motion is compared with the flow visualization. The effect of wake shapes on the aerodynamic characteristics of the flat plate in pitching oscillation is investigated. The velocity profiles behind the flat plates in pitching oscillations are plotted to investigate the possible thrust generation mechanism.

Direct numerical simulations of flow over a circular cylinder are performed at two different Reynolds numbers (Re=220 and 300) that correspond to three-dimensional instabilities of mode A and mode B, respectively, to investigate the characteristics of drag and lift at these Reynolds numbers. The drag and lift coefficients are measured locally along the spanwise direction and their characteristics are studied in detail. The variation of total drag in time is large at Re=220, and the total drag becomes minimum when vortex dislocation occurs in the wake. The drag and lift variations in space are also closely associated with the evolution of vortex dislocation at this Reynolds number. At Re=300, vortex dislocation is not found in the wake and temporal variations of drag and lift are much smaller than those at Re=220, but their spatial variations are quite large due to the near-wake secondary vortices existing in the mode B instability.

CFD data compression method based on supercompact multiwavelets is presented. High data compression can be achieved when taking advantage of the compact nature of multiwavelets. Thresholding technique is also a matter of primary concern in determining pressure ratio. In this paper, we apply thresholding for multiwavelets that considers the coefficient vector as a whole rather than thresholding individual elements. Various thresholding methods are described briefly. CFD data compression suggests that the multivariate thresholding method is suitable for supercompact multiwavelets.

Dynamic L/UL system has many merits, but it can develop an undesirable collision during dynamic loading process. In this paper, the dynamics of negative pressure pico slider during the loading process was investigated by numerical simulation. A simplified L/UL model for the suspension system was presented, and a simulation code was built to analyze the motion of the slider. A slider deigns have been simulated at various disk rotating speeds, skew angles of slider. By selection an optimal RPM and pre-skew angle, we can decrease the amount of collision and smoothen the loading process for a given slider-suspension design.

This study aimed to analyze spray characteristics and the ambient flow field in the mixture preparation state of the premixed combustion stage. It is very important to understand the spray characteristics and the fuel injection conditions in direct injection diesel engine because the emission gas compositions from diesel engines are related to spray formation processes of the premixed combustion stage. The numerical simulation was performed using the STAR-CD which is a commercial CFD code. Computed results of the transient high pressure diesel spray were compared with experimental results of the same spray injection condition in the constant volume chamber. The results show that spray patterns of numerical simulation agree with this experimental results comparatively.

The present study addresses a computational work to investigate the influence of a turbulent wake flow on the pressure recovery of a subsonic diffuser. The turbulent wake is generated by a cylinder with a small diameter, which is installed at the inlet of a 2-dimensional diffuser. Computation are applied to three-dimensional steady Navier-Stokes equations. The fully implicit finite volume scheme is used to discretize the governing equations. The computational results are qualitatively well compared to the experimental results. The results show that the pressure recovery of the subsonic diffuser is dependent on the diameter and location of cylinder. It is found that a certain diameter and location of the cylinder to generate the turbulent wake give a better pressure recovery, compared with no cylinder flow.

A turbulent wake generated by a cylinder in nozzle contraction affects to the jet flow characteristics. In this study, a computational work to investigate the effect of the turbulent wake on two-dimensional subsonic jet was carried out with three different kinds of nozzle. Computations are applied to the two-dimensional unsteady, Navier-Stokes equations. Several kinds of turbulent models and wall functions are employed to validate the computational predictions. It was known that the wake flow enhanced the spread of the jet flow, compared with no wake flow condition. It was also found that the jet core is shortened by the wake flow developed from a control cylinder.

The flow characteristics along a intake/compression process are very important for the combustion process. The intake/compression flow fields are related to the piston shape of engine. The flow fields are analysed by using the ICEM-CFD IC3M code for the rapid mesh-generation and by using the STAR-CD code for the calculations. The influences of the piston bowl shapes were investigated. The results showed that piston shapes had influences on a intake/compression flow and offered the definite basic data in a design side.

Ultrasonic flow metering(UFM) technology is being received much attention from a variety of industrial fields to exactly measure the flow rate. The UFM has much advantage over other conventional flow meter systems, since it has no moving parts, and offers good accuracy and reliability without giving any disturbances to measure the flow rate, thereby not causing pressure losses in the flow fields. In the present study, 3-dimensional, unsteady, compressible Navier-Stokes equations are solved by a finite volume scheme, based upon the second order upwind scheme for spatial derivatives and the multi-stage Runge-Kutta integral method for time derivatives. In order to simulate multi-path ultrasonic flow meter, an excited pressure signal is applied to three different locations upstream, and the pressure signals are received at three different locations downstream. The mean flow velocities are calculated by the time difference between upstream and downstream propagating pressure signals. The obtained results show that the present CFD method simulates successfully ultrasonic meter gas flow and the mean velocity measured along the chord near the wall is considerably influenced by the boundary layers.

Electroosmotic flow induced by an applied electrostatic potential field in microchannel is analyzed in this study. The electroosmotic flow is an alternative to pressure driven flow in microchannels, but the usage has been limited to the simple cases. In this study, We analyze electroosmotic flow driven by inhomogeneous surface charge on the channel wall. The surface charge varies along a direction perpendicular to the electric field in order to generate the electroosmotic flow. A numerical results substantiate the highly efficient mixing performance. It is highly the beneficial to fabrication process since only straight microchannel rather than complex geometry is enough to yield efficient mixing.

Bubble growth on microheater has been experimentally investigated in this study. The experiment was performed using platinum microheaters having dimensions of 300 ${\mu}m$ or 50 ${\mu}m$ in length, 20 ${\mu}m$ or 5 ${\mu}m$ in width, and $0.2{\pm}0.01$ ${\mu}m$ in thickness. A high speed video camera was used to observe bubble growth at 2,000 frames per second. Microheater temperature was measured at the rate of 300 Hz. with a data acquisition system. Bubble nucleation frequency increased with working fluid temperature. Although the slope of temperature drop was similar in all cases, the magnitude of temperature drop was different. The temperature profiles and the high speed camera images were combined to explain temperature drop.

In this paper, we present a technology of producing anew chaotic micromixer, named Micromixer with Arranged Blocks(MAB), and the experimental result of the mixing performance. Chaotic mixing was successfully achieved by introducing periodic perturbation in the field of the channel flow by means of slanted blocks. The MAB was made by an RP(Rapid Prototyping) technology. We performed flow visualization experiments for the quantification of the mixing performance with the MAB. Lyapunov exponent was measured to be 0.3557 and 0.1305 for the block height 0.8 and 0.2 times the channel width.

The problem of predicting the fracture strength behavior in orthotropic plate with a crack inclined with respect to the principal material axes is analyzed. Both the load to cause fracture and the crack direction of crack growth arc of interest. The theoretical results based on the normal stress ration theory show significant effects of biaxial loading and the fiber orientation on the crack growth angle and the critical stress. The additional term in the asymptotic expansion of the crack tip stress field appears to provide more accurate critical stress prediction.

Thermo-mechanical behavior of flip-chip plastic ball grid array (FC-PBGA) packages are characterized by high sensitive $Moir{\acute{e}}$ interferometry. $Moir{\acute{e}}$ fringe patterns are recorded and analyzed for several temperatures. Deformation analysis of bending displacements of the packages and average strains in the solder balls for a single-sided package assembly and a double-sided package assembly are presented. The bending displacement of the double-sided package assembly is smaller than that of the single-sided one. The largest of effective strain occurred in the solder ball located at the edge of the chip and its magnitude of the double-sided package assembly is greater than that of single-sided one.

A wind turbine obtains its power input by converting the force of the wind into a torque (turning force) acting on the rotor blades. The amount of energy which the wind transfers to the rotor depends on the density of the air, the rotor area, and the wind speed. Because it has long term operating life and very complex load condition, the fatigue strength of each component must be considered. In this paper, we calculated the load condition by wind using a combined blade elemental theory and a FEM based analytical approach was use to evaluate the fatigue strength of a Hub of wind turbine. The effect of tensile mean stress was taken into account by the modified Goodman diagram. Using this approaches, we evaluated the fatigue strength of hub and main shaft and improved the design.

Adhesive-bonded joints are widely used in the industry. Recently aircraft applications of adhesive bonding joints have been increased extensively in automobile and air industry. Because adhesives which are available for structural applications have been developed a lot and understanding of adhesive bonding has been improved so much. In this study, as the fundamental research of design of adhesive bonding joints, this study considers specimen shape are affect strength and durability of Al/Polymer lap joints. In this research, cross head speed difference were concerned to evaluate their effects on the adhesive strength. Cross head speed makes a change 0.05mm/min, 0.5mm/min, 5mm/min. The result is load-displacement diagram showed brittleness fracture tendency. Fracture tendency that is shown enough on stress distribution of trigonal single lap joint and trigonal edged single lap joint occur the inside of adhesive.

In this study, to investigate the effects of omitting low-amplitude cycles from a flight-simulation loading, crack growth tests are conducted on 2124-T851 aluminum alloy specimens. Three test spectra are generated by omitting small load ranges as counted by the rain-flow count method. The crack growth test results are compared with the data obtained from the flight-simulation loading. The experimental results show that omission of the load ranges below 5% of the maximum load does not significantly affect crack growth behavior, because these are below the initial stress intensity factor range. However, in the case of omitting the load ranges below 15% of the maximum load, crack growth rates decrease, and therefore crack growth curve deviates from the crack growth data under the flight-simulation loading. To optimize the load range that can be omitted, crack growth curves are simulated by the stochastic crack growth model. The prediction shows that the omission level can be extended to 8% of the maximum load and test time can be reduced by 59%.

Aluminum carbody for rolling stocks is light and perfectly recycled, but includes severe defects which are very dangerous to fatigue strength. Structural integrity assessment for the carbody by static load test has been performed up to date. In this study, to evaluate fatigue strength of the aluminum carbody of urban transit unit. a testing method to simulate dynamic loading condition was proposed and the fatigue strength of the carbody was evaluated. The dynamic load test results showed that the alternating stress ranges were different from the estimated ranges based on the static test results. Excessive stress ranges at the center are thought to come from the flexible motion of the carbody. published fatigue test data for aluminum components, but variation of alternating acceleration along the length due to flexibility of carbody yielded unexpected results. Because fatigue strength based on the static test results may be overestimated at the center, modification of testing method is necessary.

CD-R disc drive is one of the basic options of personal computer today. In this trend, various CD-R discs can be purchased in every market. Even very low quality discs are in market too. Today's CD-R disc drive is being operated over 10,000 RPM. So the possibility of the disc fracture is growing faster. Sometimes, during the test of a new disc drive by various methods to confirm the quality of it, fracture of the disc happens. And it happens in end user's personal computer scarcely. In this reason, new methods to confirm the quality and the failure mechanism of the optical disc are studied in this paper.

In this study, corrosion fatigue characteristics of CF8M and CF8A steel were investigated on the simulated PWR condition(Temp.:$316^{\circ}C$, Pres.: 15:MPa). To make the simulated PWR condition. the special test machine consisted of INSTRON, Autoclave, LOOP and Measurement system was developed. As ${\Delta}K$ is ranged from 11 to $20MPa{\sqrt{m}}$, Crack growth rate of PWR condition is faster than air condition. Above $20MPa{\sqrt{m}}$, the crack growth rate of PWR and air condition is similar. Corrosion fatigue characteristics regardless of the ferrite contents($10{\sim}25wt.%$) is not different. After the test, the fracture surface of specimens was examined. It was difficult to verify the fracture modes such as striation, intergranular crack and cleavage and so on. As the ferrite content of CF8M is increased, the more particles covered fracture surface were peeled.

Method for contact failure mitigation is studied in this paper. The focus is laid on the contact shape that eventually influences the internal stresses. Contact mechanics is consulted within the frame of plane problem. Hertzian contact, rounded punch and uniform traction profiles are considered. Frictional as well as frictionless contact is also considered. As results, the higher traction profile induced by the rounded punch reveals the greatest among the considered shapes. Therefore, it is suggested to increase the edge radius as large as possible if a contact body of punch shape needs to be designed. It is also found that uniform traction cannot always provide the solution of contact failure mitigation.

In this paper, we conducted finite element analysis to investigate the residual stress redistributions of weldment due to cutting. To evaluate the effect of the residual stress on the fatigue behavior of weldment, test specimens are commonly cut from the weldment, but the distributions of the residual stress in the cut specimen should be different from those in the original weldment. Our work is to evaluate the difference between the residual stresses before and after weldment-cutting to understand the effect of cutting on the residual stress. Transient heat analysis, elastic-plastic mechanical analysis and element removal technique are used to simulate the welding and cutting procedures on the commercial finite element code ABAQUS.

The flow of electrical current through a microscopic actual contact spot between two conductors is influenced by the flow through adjacent contact spots. A smoothed version of this interaction effect is developed and used to predict the contact resistance when the statistical size and spatial distribution of contact spots is known. To illustrate the use of the method, an idealized fractal rough surface is defined using the random midpoint displacement algorithm and the size distribution of contact spots is assumed to be given by the intersection of this surface with a constant height plane. With these assumptions, it is shown that including finer scale detail in the fractal surface, equivalent to reducing the sampling length in the measurement of the surface, causes the predicted resistance to approach the perfect contact limit.

The problem of residual stresses and fatigue behavior in welded structures is the main concern of welding research fields. The residual stresses and distortions of structures by welding exert negative effect on the safety of mechanical structures. Postweld heat treatment is usually carried out to relieve this residual stresses of welded joints. In this paper the influence of postweld heat treatment on fatigue life of butt-welded joint was investigated. To predict the effect of PWHT, an analytical model is developed by finite element and local strain approach and the result of fatigue life analysis is compared to experimental results. It is demonstrated that fatigue life estimates closely approximate the experimental results and PWHT provides some increase of fatigue lives in long-life fatigue region and no increase in short-life fatigue region because of the residual stress relaxation under tensile loads.

S-N fatigue tests were conducted to investigate the fatigue strength of small diameter socket and butt welded joints made of carbon steels. Experimental parameters were pipe diameter, throat depth, shape of socket welds and welding procedure. Filler metals used in SMAW and GTAW procedure were E9016-G with diameter of 4.0 mm and ER70S-G with diameter of 2.4 mm. API 5L Gr.B pipes were adopted as a small diameter branch pipes. All socket fittings were machined from ASTM A105 carbon steel. Tensile strength was not affected by the welding procedure. Fatigue strength in socket weld joints increased with increasing pipe diameter, area of weld metal and weld leg length of pipe side.

Vacuum vessel of the KSTAR(Korea Superconducting Tokamak Advanced Research) tokamak is a fully welded structure with D-shaped cross-section. According to the requirements of the physics design, sixteen horizontal ports, sixteen slanted ports, sixteen baking and cooling ports, and twenty-four top and bottom vertical ports are designed for the diagnostics, plasma heating, vacuum pumping, and baking and cooling. Bellows on these ports are used for flexible components to absorb the relative displacement due to the vacuum vessel thermal expansion and the electromagnetic force between the vacuum vessel and the cryostat ports. Fatigue strength evaluation was performed to decide the dimension of the bellows. In order to assure the quality of the bellows, a prototype bellows for the neutral beam injection port has been fabricated and tested prior to main fabrication. It was conformed that the prototype bellows has sufficient fatigue strength and vacuum reliability in the expected load conditions.

The fatigue strength of welded joint is influenced by the welding residual stress which is relaxed depending on local stress distributed in vicinity of stress raisers, eg. under cut, overlap and blow hole. To evaluate its fatigue life the geometry of the stress raisers and the welding residual stress should be taken into account. The several methods based on notch strain approach have been proposed in order to consider the two factors above mentioned. These methods, however, have shown considerable differences between analytical and experimental results. It is due to the fact that the amount of the relaxed welding residual stress evaluated by the cyclic stress-strain relationship do not correspond with that occurred in reality. In this paper the residual stress relaxation model based on experimental results was used in order to reduce the discrepancy of the estimated amount of the relaxed welding residual stress. Under an assumption of the superimposition of the relaxed welding residual stress and the local stress, a modified notch strain approach was proposed and verified to the cruciform welded joint.

The initial crack often occurs on the bonded interface and it is the general cause of the interface fracture. It is very significant to establish the measurement method of interfacial crack by applying the ultrasonic technology into the interface of bonded dissimilar materials. In this paper, the interfacial crack length was measured by ultrasonic attenuation coefficient in the Al/Epoxy bonded dissimilar materials of double-cantilever beam(DCB). The energy release rate, G, was obtained by the experimental and Ripling's equation measurement of compliance. The experimental results represent that the relation between interfacial crack length for the ultrasonic attenuation coefficient and energy release rate is increased proportionally. From the experimental results, a measurement method of the interfacial crack length by the ultrasonic attenuation coefficient was proposed and discussed.

The recent tendency in the automobile industries is toward light weighting vehicle body to improve the problems by environmental pollution as well as improving fuel cost. The effective way to reduce the weight of vehicle body seems to be application of new materials for body structure and such trend is remarkable. Among the various materials for vehicle body, stainless steel sheet (for example, 301L and 304L), TRIP steel and cold rolled steel sheets are under the interests. However, in order to guarantee reliability of new material and to establish the long life design criteria of body structure, it is important and require condition to assess spot weldability of them and fatigue strength of spot welded lap joints which were fabricated under optimized spot welding condition. And, recently, a new issue in the design of the spot welded structure is to predict economically fatigue design criterion without additional fatigue tests. In general, for fatigue design of the spot-welded thin sheet structure, additional fatigue tests according to the welding condition, material, joint type, and fatigue loading condition are generally required. This indicates that much cost and time for it should be consumed. Therefore, in this paper, the maximum stresses at nugget edge of spot weld were calculated through nonlinear finite element analysis first. And next, obtained the ${\Delta}P-N_{f}$ relation through the actual fatigue tests on spot welded lap joints of similar and dissimilar high strength steel sheets. And then, the ${\Delta}P-N_{f}$ relation was rearranged in the ${\Delta}{\sigma}-N_{f}$ relation. From this ${\Delta}{\sigma}-N_{f}$ relation, developed the fatigue design technology for spot welded lap joints of them welded using the optimized welding conditions.

Strength evaluation was carried out for the fillet and oil hole of crankshaft of medium speed diesel HiMSEN engine to verify initial concept design. Alternating torque obtained from torsional vibration analysis and radial force due to firing pressure were applied. It was assumed that the maximum alternating torque and radial force occur simultaneously. Weak points in view of fatigue are proceeding fillet and crank pin fillet area and the minimum normalized fatigue safety factor is 1.17 at crank pin fillet. The fatigue strength of the oil hole was evaluated to verify the effect of the surface roughness of the oil hole. As results, the specific level of the inner surface roughness and the polishing depth of the oil hole for sufficient fatigue strength was suggested. The maximum stress value and stress distribution at the inner surface of the oil hole can be easily estimated at initial design stage by the newly developed method.

Suspension system of vehicle have enough endurance during its life time to protect passenger. Spring is one of major part of vehicle. Thus, a fatigue strength evaluation for leaf spring based on road load response was carried out. At first, strain of leaf spring is measured on the city condition and proving ground condition. And next, the damage analysis of road load response data was carried out. And fatigue test of leaf spring were also carried out. Based on -N life relation, fatigue life of leaf spring was evaluated at belgian mode, city mode and drawing test specification. After that, it is compared the design life of leaf spring and evaluated fatigue life by belgian mode, city mode and drawing test specification. From the above, the maximum load-fatigue life relation of leaf spring was defined by test. and new test target of belgian mode and city mode was proposed to accept design specification of leaf spring. It is expect that proposed test target can verify leaf spring fatigue endurance at specific road condition.

The pressure tubes, which contain high temperature heavy water and fuel, are within the core of a CANDU nuclear reactor, and are thus subjected to high stresses, temperature gradient, and neutron flux. Further, it is well known that pressure tubes of cold-worked Zr-2.5Nb materials result in hydrogen diffusion, which create fully-hydrided regions (frequently called Blister). Thus a proper investigation of hydrogen diffusion within zirconium-alloy nuclear components, such as CANDU pressure tube and fuel channels is essential to predict the structural integrity of these components. In this respect, this paper presents numerical investigation of hydrogen diffusion to quantify the hydrogen concentration for blister growth of CANDU pressure tube. For this purpose, coupled temperature-hydrogen diffusion analyses are performed by means of two-dimensional finite element analysis. Comparison of predicted temperature field and blister with published test data shows good agreement.

This paper addresses a theoretical study to calculate the amount of the stored energy due to vacancies during high-strain-rate deformation. The study concerns the role of excess vacancies, which can play an important role to increase the amount of stored energy. Molecular dynamics simulation using a 3D model is carried out and the result clearly shows that the excess vacancies are credited to generation of the stored energy.

Zr-based bulk metallic glasses have a significant mechanical properties such as high strength and elastic strain limit, and a good processing ability due to the deformation behavior such as superplasticity under supercooled liquid region. Recently, many researches on the determination of optimum working condition in various bulk metallic glasses have been carried out. In this study, the deformation behavior and forming conditions of $Zr_{55}Cu_{30}Al_{10}Ni_{5}$ bulk metallic glass were investigated under three different strain rates and at various temperatures between 627K and 727K. The glass transition temperature, crystallization temperature and supercooled liquid region of $Zr_{55}Cu_{30}Al_{10}Ni_{5}$ bulk metallic glass are 680K, 762K and 82K, respectively.

The dynamic properties of short-fiber reinforced Chloroprene rubber for vibration isolators have been studied as functions of interphase conditions and fiber content. The loss factor showed the maximum at strain amplitude 2%, and increased 0.09 for matrix, 0.05 for reinforced rubber with increasing frequency respectively. The dynamic ratio rapidly decreased with increasing strain amplitude, and some increased with increasing frequency. The better interphase condition showed the lower dynamic ratio. Therefore, the short-fiber reinforced rubber could have the better isolation in frequency ratio(${\sqrt{2}}min$.) compared to frequency ratio(${\sqrt{2}}max$.). And we have investigate the possibility of applying short-fiber reinforced rubber to automotive engine mount.

Behavior of aluminum alloy embedding a particle was investigated at room temperature under ECAP. Finite element analysis by using ABAQUS shows that ECAP is a useful tool for eliminating residual porosity in the specimen, and much more effective under friction condition. The simulation, however, shows considerably low density distributions for matrix near a particle at which rich defects may occur during severe deformation. Finite element results of effective strains and deformed shapes for matrix with a particle were compared with theoretical calculations under simple shear stress. Also, based on the distribution of the maximum principal stress in the specimen, Weibull fracture probability was obtained for particle sizes and particle-coating layer materials. The probability was useful to predict the trend of more susceptible failure of a brittle coating layer than a particle without an interphase in metal matrix composites.

Effective properties of ceramic matrix plain woven textile composites were calculated using finite element analysis. The considered geometry is a unit cell of plain weave and the analysis was performed by commercial finite element program, ANSYS. The materials for analysis are 3 types for matrix, 1 type for fiber with various volume fraction. The result indicates that the effective properties of ceramic matrix composites can be controlled by the volume fraction. The result can be used for numerical analysis using ceramic matrix composites.

When the tube contacted to support and antivibration bar of the steam generator in unclear power plant, the contact area is worn out by their relative displacement. In the study, wear depths of the tube inclined to tube support and antivibration bar are approximately predicted by a method using the contact load and relative displacement. In the case of the inclined contact, the results show wear depths of the steam generator tube predicted by the impact model are larger than those by the sliding model.

Cast austenitic stainless steel is used for several components, such as primary coolant piping, elbow, pump casing and valve bodies in light water reactors. These components are subject to thermal aging at the reactor operating temperature. Thermal aging results in spinodal decomposition of the delta-ferrite leading to increased strength and decreased toughness. This study shows that ferrite content can be predicted by use of the artificial neural network. The neural network has trained learning data of chemical components and ferrite contents using backpropagation learning process. The predicted results of the ferrite content using trained neural network are in good agreement with experimental ones.

One of the main degradation of steam generator tubes is stress corrosion cracking induced by residual stress. The resulting damages can cause tube bursting or leakage of the primary water which contained radioactivity. Primary water stress corrosion crack occurs at the location of tube/tubesheet hard rolled transition zone. In order to investigate the effect of shot peening on stress corrosion cracking, stress intensity factors are calculated for the crack which is located in the induced residual stress field.

The 40% of wall criterion, which is generally used for the plugging of steam generator tubes, may be applied only to a single crack. In the previous study, a total of 9 failure models were introduced to estimate the local failure of the ligament between cracks and the optimum coalescence model of multiple collinear cracks was determined among these models. It is, however, known that parallel axial cracks are more frequently detected during an in-service inspection than collinear axial cracks. The objective of this study is to determine the plastic collapse model which can be applied to the steam generator tube containing two parallel axial through-wall cracks. Nine previously proposed local failure models were selected as the candidates. Subsequently interaction effects between two adjacent cracks were evaluated to screen them. Plastic collapse tests for the plate with two parallel through-wall cracks and finite element analyses were performed for the determination of the optimum plastic collapse model. By comparing the test results with the prediction results obtained from the candidate models, a plastic zone contact model was selected as an optimum model.

The aim of this study is to investigate the hydride embrittlement when the LBB evaluation is carried out for the integrity of PHWR Pressure Tubes. The transverse tensile and CCT tests were performed at three hydrogen concentrations while the test temperatures were changed (RT to $300^{\circ}C$). The specimens were directly machined from the pressure tube retaining original curvature. Both the transverse tensile and the fracture toughness tests showed the hydrogen embrittlement clearly at RT but this phenomenon was disappeared while the test temperature arrived over $250^{\circ}C$. Using the DHC test results, the CCL and LBB time were calculated and compared. The hydride embrittlement behavior at the LBB evaluation was definitely showed.

Most domestic fossil power plant have exceeded 100,000 hours of operation with the severe operating condition. Among the critical components of fossil power plant, high temperature steam pipe system have had a many problems and damage from unstable displacement behavior because of frequent start up and shut down. In order to prevent the serious damage and failure of the critical pipe system in fossil power plant, 3-dimensional displacement measurement system were developed for the on-line monitoring system. 3-D Measurement system was developed with using the LVDT type sensor and rotary encoder type sensor, this system was installed and operated on the real power plant successfully. In the future time, network system of on-line diagnosis for critical pipe will be designed.

Effect of transverse electric field on crack kinking in ferroelectric ceramics subjected to purely electric loading is investigated. It is shown that the shape and size of the domain switching zone depends strongly on the direction of the applied electric field as well as the ratio of the transverse electric field to the coercive electric field. Under small-scale conditions, mode I and II stress intensity factors induced by ferroelectric domain switching are numerically obtained. The crack kinking in ferroelectrics is also discussed.

Pressure tubes are major component of nuclear reactor, but only selected samples are periodically examined due to numerous numbers of tubes. Current in-service inspection result show there is high probability of flaw existence at un-inspected pressure tube. Probabilistic analysis is applied in this study for the integrity assessment of un-inspected pressure tube. But all the current integrity evaluations procedures are based on conventional deterministic approaches. So many integrity evaluation parameters are not directly apply to probabilistic analysis. As a result of this study failure assessment diagram are proposed based on test data.

The specimen materials used in this research is bimaterial. The static photoelastic experiment was applied to them. And then the specimens used in photoelastic experiment were fractured under static load. The static photoelastic hybrid method was introduced and it's validity had been assured. The static photoelastic hybrid method was applied to the Minimum Strain Energy Density Criterion, the Maximum Tangential Stress Criterion and Mode Mixity. Crack propagation criterion by the static photoelastic hybrid method was introduced and it was applied to the above various failure theories. Comparing the experimental initial angle of crack propagation with the theoretical initial angle of crack propagation from the various failure criterions. And then the optimal crack propagation criterion was suggested and it's validity was assured.

Pitting wear is a dominant form of polyethylene surface damage in total knee replacements, and may originate from surface cracks that propagate under repeated tribological contact. In this study, stress intensity factors, $K_{I}$ and $K_{II}$, were calculated for a surface crack in a polyethylene - CoCr - bone system under the rolling and/or sliding contact pressures. Crack length and load location were considered in determination of probable crack propagation mechanisms and fracture modes. Positive $K_{I}$ values were obtained for shorter cracks in rolling contact and for all crack lengths when the sliding load was apart from the crack. $K_{II}$, was the greatest when the load was directly adjacent to the crack $(g/a={\pm}1)$. Sliding friction caused a substantial increase of both $K_{I}^{max}$ and $K_{II}^{max}$. The effective Mode I stress intensity factors, $K_{eff}$, were the greatest at $g/a={\pm}1$, showing the significance of high shear stresses generated by loads adjacent to surface cracks. Such behavior of $K_{eff}$ suggests mechanisms for surface pitting by which surface cracks may propagate along their original plane under repeated rolling or sliding contact.

The molecular cluster model for the homogeneous bubble nucleation was extended to predict the bubble nucleation events in elastomers, polymers and polymer solutions. For the bubble nucleation in elastomers and polymers, the strain energy overcome by a critical bubble was also considered. The calculation results for the number of bubbles nucleated are in good agreement with observed ones.

A failure analysis of holddown spring screw was performed using fracture mechanics approach. The spring screw was designed such that it was capable of sustaining the loads imposed by the initial tensile preload and operational loads. In order to investigate the cause of failure, a stress analysis of the top nozzle spring assembly was done using finite element analysis and a life prediction of the screw was made using a fracture mechanics approach. The elastic-plastic finite element analysis showed that the local stresses at the critical regions of head-shank fillet and thread root significantly exceeded than the yield strength of the screw material, resulting in local plastic deformation. Primary water stress corrosion cracking life of the Inconel 600 screw was predicted by using integration of the Scott model and resulted in 1.42 years, which was fairly close to the actual service life of the holddown spring screw.

Reduced activation ferritic steel (JLF-1) is considered as a promising candidate material for blanket or first-wall structure of D-T fusion reactors. The fracture tests of fracture resistance curve (J-R curve) and $J_{IC}$ are desirable to investigate the exact fracture toughness of JLF-1 steel, since it has a high ductility. The fracture toughness of JLF-1 steel is affected by the configuration of test specimen such side groove, specimen thickness or specimen size. In this study, the fracture toughness tests were performed with various size(plane size and thickness) and various side groove of specimens. The test results showed the standard specimen with the side groove of 40 % represented a valid fracture toughness. The fracture resistance curve increased with increasing plane size and decreased with increasing thickness. However, the fracture resistance curve of half size specimen was similar to that of the standard specimen.

The fracture behaviors of Zr-based bulk amorphous metals(BAMs) having compositions of $Zr_{55}Al_{10}Ni_{5}Cu_{30}$, were investigated under impact loading and quasi-static conditions. For experiments, a newly devised instrumented impact testing apparatus and the subsize Charpy specimens were used. The influences of loading rate and the notch shape on the fracture behavior of the Zr-based BAM were examined. The Zr-based BAMs showed an elastic deformation behavior without any plastic deformation on it before fracture. Most fracture energies were absorbed in the process of the crack initiation. The maximum load and fracture absorbed energy under quasi-static condition were larger than those under impact condition. However, there existed relatively insignificant notch shape effect. Fracture surfaces under impact loading were smoother than those under quasi-static loading. The absorbed fracture energy appeared differently depending on the extent of the vein-like pattern region due to the shear bands developed at the notch tip. It can be found that the fracture energy of the Zr-Al-Ni-Cu alloy is closely related with the development of shear bands during fracture.

This paper present a procedure of meter-out flow control method for dump valve in full-scale airframe test. Emergency stop, which results in dump state, can be happened during full-scale airframe test by several causes. Because servo valve can't control hydraulics actuator in the dump state, pressure in cylinder chamber may rise abruptly and overload can be acted to the test article. In this paper, the procedure and technology of orifice setting are investigated to protect the test article from unexpected loads by dump. The test results show that the presented methods decrease peak loads and improve unloading characteristics of hydraulic actuators in the dump state.

Oil supply system is one of the most important part of Turbine. Lubricating oil of bearings supplied by oil pump. Failure of Oil supply pump critical damaged parts of Turbine, especially bearings. In this paper we have discussed the serious damage of turbine bearings due to failure of Oil supply pump.

TF magnet structures are the main structural components in the KSTAR magnet systems to protect the superconducting coils from mechanical, electrical, and thermal loads. TF coil structure supports CS and PF coil system. The inter-coil structure contains adjustable shear keys and conical bolts to provide pre-loading in toroidal direction and to resist against in-plane and out-of-plane forces that are the most critical loads on the TF magnet system. The conical bolts and shear keys are specially designed to assemble easily and to provide a convenient accommodation for a good alignment. The connection plate that is one of the prototype fabrications had been manufactured to study adjustability of conical bolts and shear keys for assembly of TF coil structure. We could measure the misalignments at the keyways and conical holes with the misalignment measuring instrument.

This paper presents a case history of piping failures on power plant. The root cause of the failure was defined to set the optimal countermeasures. The failure comes from transient vibration and the 1st stress increased at the hydraulic oil supply system of control valves for high pressure steam turbine.

In this paper, a new meshfree technique which improves the numerical integration accuracy is introduced. This new method called the Petrov-Galerkin natural element(PG-NE) is based on the Voronoi diagram and the Delaunay triangulation which is based on the same concept used for conventional natural element method called the Bubnov-Galerkin natural element(BG-NE). But, unlike BG-NE method, the test shape function is differently chosen from the trial shape function. The proposed technique ensures that the numerical integration error is remarkably reduced.

In the developement of LNG cargo, the current concern focuses on the slim design of insulation layer to increase the LNG carrying capacity. Not only thermal stability with BOR(Boil-Off Rate) but structual stability against the LNG weight and the sloshing phenomenon must be also considered. In this paper, we applied the stitched sandwitch composite called the 3X-Board which is stitched through the core thickness direction using glass fiber to the LNG cargo. We evaluated the thermal and structural characteristics of 3X-Board by changing the core thickness and the material, in order to explore a validity for the slim design through the finite element analysis.

This paper presents a newly implemented parallel finite element procedure for contact-impact problems. Three sub-algorithms are includes in the proposed parallel contact-impact procedure, such as a parallel Belytschko-Lin-Tsay (BLT) shell element generation, a parallel explicit time integration scheme, and a parallel contact search algorithm based on the master slave slide-line algorithm. The underlying focus of the algorithms is on its effectiveness and efficiency for inclusion in future finite element systems on parallel computers. Throughout this research, a prototype code, named GT-PARADYN, is developed on the IBM SP2, a distributed-memory computer. Some numerical examples are provided to demonstrate the timing results of the procedure, discussing the accuracy and efficiency of the code.

The objective of the present research is to develop a wavelet-based multiscale adaptive Galerkin method for membrane eigenvalue analysis. Since approximate eigensolutions at a certain resolution level can be good guesses, which play an important role in typical iterative solvers, at the next resolution level, the multiresolution iterative solution approach by wavelets can improve the solutionconvergence rate substantially. The intrinsic difference checking nature of wavelets can be also utilized effectively to develop an adaptive strategy. The present wavelet-based approach will be implemented for the simplest vector iteration method, but some important aspects, such as convergence speedup, and the reduction in the number of nodes can be clearly demonstrated.

KAERI has contrived 15 kinds of spacer grid shapes of its own since 1997 and applied for domestic and foreign patents. To date, KAERI has obtained US and ROK patents for 6 kinds of spacer grid shapes among them and the others are under review in USA, EC, China, and ROK. In this study, mechanical/structural performance analysis and test on two spacer grid shapes that are assumed to be the most effective candidates for the spacer grid of the next generation nuclear fuel in Korea was carried out. The result has shown that the performances of the candidates are better or not worse than those of the current spacer grid.

This paper presents the static and dynamic measurements for the strength and motion characteristics as well as the improved procedures to assess strength of wheel loaders. Two scenarios for static measurement were decided by which cylinder was actuating. The dynamic measurement was performed for two types of motion, that is, simple reciprocation of the working devices and actual working motion including traveling, digging and dumping. The measured items were stresses, cylinder pressures and strokes. Stress induced by bucket working showed higher level than that by boom working. The measured cylinder speeds were relatively superior to the design speeds. Working stress histories were thought to be closer to static rather than dynamic. A fully assembled FE model was prepared for structural analysis. In this paper, a more simple method was suggested to avoid nonlinearity caused by heave of rear frame under digging forces. Also how brake affected on structural behavior and digging force was examined closely in relation with tire pressure. It was confirmed that the overall stress level of wheel loader during turning traveling with loaded bucket was far lower than the yield stress of material.

In many practical engineering design problems, there are some design and manufacturing considerations that are difficult or infeasible to express in terms of an objective function or a constraint. In this situation, a set of optimal candidate designs having different topological complexities, not just a single optimal design, is preferred. To generate systematically such design candidates, we propose a hierarchical multiscale design resolution control scheme. In order to adjust its topological complexity by choosing a different starting resolution level in the hierarchical design space, we propose to employ a general M-band wavelet transform in transforming the original design space into the multiscale design space.

The systematic method to construct equivalent static load from the given dynamic load is proposed in the present study. Previously reported works to construct equivalent static load were based on ad hoc methods. They may results in unreliable structural design. The present study proposes a selection scheme of degrees of freedom(d.o.f) for imposing the equivalent static loads. The d.o.fs are selected by Two-level condensation scheme(TLCS). TLCS consists of two two-steps. The first step is the energy estimation in element-level and the second step consists of the traditional sequential elimination precudure. Through several numerical examples, the efficiency and reliability of proposed scheme is verified.

It is important to determine supporting locations for structural stability of a structure under non-uniform load in space interfered by other parts. In this case, There are many local optima with discontinuous design space. Therefore, The traditional optimization methods based on derivative are not suitable. Whereas, Genetic algorithm(GA) based on stochastic search technique is a very robust and general method. This paper has been presented to determine supporting locations of the vertical supports for reducing stress of the KSTAR(Korea super Superconducting Tokamak Advanced Research) IVCC(In-vessel control coil) under non-uniform electromagnetic load and space interfered by other parts using genetic algorithm. For this study, we develop a program combining finite element analysis with a genetic algorithm to perform structural analysis of IVCC. In addition, this paper presents a technique to perform optimization with FEM when design variables are trapped in an incongruent design space.

In the present study, an integrated framework of geometric modeling, analysis, and design optimization is proposed. Geometric modeling is based on B-spline surface representation. Geometrically-exact shell finite element is implemented in analysis module. Control points of the surface are selected as design variables for optimization, which can make the interaction easier between analysis and surface representation. Sequential linear programming(SLP) is adopted for the shape optimization of surfaces. For the computation of shape sensitivities, semi-analytical method is used. The developed integrated framework should serve as a powerful tool for the geometric modeling, analysis, and shape design of surfaces.

An efficient method is developed for the shape optimization of 2-D structures. The sequential linear programming is used for minimization problems. Selected set of master nodes are employed as design variables and assigned to move towards the normal direction. After adapting the nodes on the design boundary, the B-spline curves and mesh smoothing schemes are used to maintain the finite element in good quality. Finally, a numerical implementation of optimum design of an automobile torque converter piston subjected to pressure and centrifugal loads is presented. The results shows additional weight up to 13% may be saved after the shape optimization.

In this paper, an integrated optimal design software system for structural components has been developed which interfaces existing commercial codes for CAD, CAE and Optimization. They include specialized optimal design software codes such as iSIGHT and VisualDOC, optimization module imbedded in CAD software developed by CAD developers, and optimal design software systems based on API of commercial CAD software. The advantages of the CAD imbedded optimal design approach and those of specialized optimal design software are taken to develop the system. The user defines optimal design formulation in the user interface for problem definition in the CAD control stage, where design variables are directly selectable from the CAD model and various properties and performance functions defined. The commercial CAD codes, Open I-DEAS are used for the development. The resulting software is minimally connected to CAD and CAE systems while keeping maximum independence from each other. This assures flexibility and freedom for problem definition. Fatigue life optimization is taken as a nontrivial application area. As a specific example, the shape design of a knuckle part of an automobile is performed, where the minimum fatigue life over the material domain in terms of the number of cycles of a curb strike are maximized under the constraint of not exceeding the current mass. The fatigue life has been improved by four times of the initial life. The developed software is illustrated to maintain the advantages of existing optimal design software systems while improving independency and flexibility.

There is few research about contact problem in SPH because it is primarily suitable to analyze the large deformation problem. However, an elasto-plastic problem with small deformation need to be considered about contact characteristics. The numerical formulating methods for SPH is induced to be able to obtain solutions based on a variational method in contact problem. The contact algorithm presented is applied to the elastic impact problem in 1D and 2D. The results show thai an imaginary tension and a numerical instability which happen in impacting between different materials can be removed and contact forces which could not have been calculated are able to obtain.

The effect of the metal mold on densification behavior of stainless steel 316L powder was investigated under warm isostatic pressing with metal mold. We use lead as metal mold and obtain experimental data of metal mold property. To simulate densification of metal powder, the elastoplastic constitutive equation proposed by Shima and Oyane was implemented into a finite element program (ABAQUS) under warm die pressing and warm isostatic pressing with metal mold. Finite element results were compared with experimental data for densification and deformation of metal powder under warm isostatic pressing and warm die pressing.

Due to the application of Integrated Head Assembly (IHA) in KSNP+ reactor design, an investigation of reactor internals integrity is carried out to assure that the adoption of IHA does not affect the safety of reactor operation. One of the postulated accident events is the R.V. closure head fall from 5.5m high directly above the reactor vessel that may occur during the refueling operation. The analysis model consists of lumped mass elements of the entire reactor vessel and internals. Because of extreme load, separate elastic-plastic analyses are done for the members that undergo plastic deformation. The analysis verified that the stresses of the reactor internals and the fuel assemblies are within the bound of allowable stress limits and the integrity of the fuel assemblies is maintained.

These days, spinal interbody arthrodesis using fusion cage is very popular. The cage used for the spinal interbody arthrodesis is mainly inserted from the posterior of the spine. Accordingly, there could possibly occur damages at posterior and results in instability of structure. Moreover, one or two cages are inserted depending on the patients. In this study, it is attempted to evaluate the stability quantitatively by comparing two cases where one and two cages are inserted. For this purpose, a very fine 3-dimensional finite element model of vertebra is generated from the MRI data. From this vertebra model, two models are made: one with one cage and the other with two cages. Finally, finite element analys is performed for these two models and both of the mechanical behaviors are examined In addition, the effect on the stability is evaluated and compared quantitatively.

PMMA which is used as the bone cement for vertebroplasty is able to be a supporter, as a fixing supporter role, for broken trabecular structure, caused by the compressed fracture of spine on aged osteoporosis. In this thesis, as experimenting apparent density of bone pieces, we have figured out support extent of Young's modulus as classifying the bone pieces injected PMMA and the others which are not. In case of low apparent density of PMMA in some bone, Young's modulus seems to be more supportable to bone. On the other hand, if apparent density of bones is normal, injection of PMMA is not very effective on improvement in Young's modulus of bone cement injection.

Current prostheses for amputees are generally extrinsic wearing socket type that the coupling between body stump and appliance wraps the soft tissue and this structure causes several problems :applying direct weight to soft tissue such as skin and muscle, skin trouble of contacting area and pain. In this study, osseointegration implant is a method to directly connect prosthesis to the residual stump skeletal tissue of arm, finger and leg through surgical operation. Technology presented in this paper essentially solves the problems of pain and abnormal weight transfer system indicated above and recovers the functions of the amputated arm and leg. In this paper, implant shape was designed for the first step for the development of osseointegration implant and then we studied the possibility to apply this osseointegration implant to human body by performing implant insertion operation to beagle tibia for the clinical animal test and normal beagle's gait analysis was executed in order to quantitatively verify the beagle's skeletal functions after the implant insertion.

A ferroelectric (called piezoelectric afterwards) wafer has been widely used as a key component of actuators or sensors of a layer type. According to recent researches, the piezoelectric wafer behaves in a nonlinear way under excessive electro-mechanical loadings. In the present paper, one-dimensional constitutive equations for the nonlinear behavior of a piezoelectric wafer are proposed based on the principles of thermodynamics and a simple viscoplasticity theory. The predictions of the developed model are compared with experimental observations.

Pressurization test is usually required in aircraft full-scale static test. There are several test conditions including pressurization of cockpit, fuselage fuel tank, air inlet duct for T-50 full-scale static test. In this paper, the test conditions, equipment, piping analysis for the pressurization test are introduced. Tank simulation test is performed to verify the validity of piping analysis and to find good tuning parameters for the pressurization channel in the servo controller. Several test setup for pressurization of T-50 test is introduced. Test article volume is filled by form, $60%{\sim}80%$ volume is reduced for the T-50 full-scale static test. Pressurization system is connected to servo controller which also controls hydraulic actuator. Load and pressure control is synchronized by using the same servo controller during T-50 test. Typical control result for pressurization test condition is shown. Pressurization tests of T-50 full-scale static test was completed successfully.

In this study, analytic stress-displacement solutions are obtained by using a shear lag modeling constructed for the spliced joint area with a splicing gap filled with adhesive material of elastic modulus $E_{a}$ in the fiber metal laminate (FML) which is known to have excellent fatigue, corrosion and fire-flame resistant characteristics while with relatively low densities compared to the conventional aluminum alloys for lightweight structures.

Tire subjected to inflation pressure and the load of passenger car has a important roll to dynamic performance characteristics of radial tires such as cornering and braking. These performance characteristics are directly related to the belt cord durability of tires. In this paper, we analyze the mechanism of belt cord broken due to cord tension at braking theoretically.

Mechanical behavior of a piezoelectric actuator is studied as a preceding research for the manufacturing of three-dimensional micro-structures. It is needed to examine the simulation of a piezoelectric actuator according to applied direction of voltage, by researching displacement characterization of piezoelectric material through piezoelectric theory. To this end, finite element modeling is employed to study the response of a piezoelectric material under the various input voltages. Where the actuator is simulated by use of ANSYS. To avoid direct contact piezoelectric material with working fluid, silicon, polymer, etc., the actuator is modeled with nickel fixed diaphragm.

In MEMS devices, packaging induced stress or stress induced structure deformation become increasing concerns since it directly affects the performance of the device. In this paper, deformation behavior of MEMS gyroscope package subjected to temparature change is investigated using high-sensitivity $Moir{\acute{e}}$ interferometry. Using the real-time $Moir{\acute{e}}$ setup, fringe patterns are recorded and analyzed at several temperatures. Temperature dependent analyses of warpages and extensions/contractions of the package are presented. Linear elastic behavior is documented in the temperature region of room temperature to $125^{\circ}C$. Analysis of the package reveals that global bending occurs due to the mismatch of thermal expansion coefficient between the chip, the molding compond and the PCB.

In this study, the topology optimization is applied to the design of a piezoelectric microactuator satisfying the specific mean transduction ratio(MTR). The optimization problem is formulated to minimize the difference between the specified and the current mean transduction ratio. In order to analyze the response of the piezoelectric-structure coupled system, both the structural and the electric potential are considered in the finite element method. The optimization problem is resolved by using Sequential Linear Programming(SLP) and the results of test problems show that the design of a piezoelectric microactuator with specified mean transduction ratio can be obtained.

Impact analysis of TFT-LCD module is very complicated because the structure is assisted with thin, small and non-uniform geometry. Especially, finite element modeling is more difficult and need time-consuming efforts. In this study, we developed LCD Impact Analysis System (LIAS) for the purpose of reducing the analysis time without accuracy reduction. This system contains pre-meshing data, material database, shock condition, auto-reporting etc. PATRAN and DYNA3D is used for meshing and solving. Previously, we performed impact test and reviewed the accuracy of analysis results. Simply we can control design parameters, the procedure such as meshing, running and reporting which are partially auto-prepared. By adopting proposed system, it is expected to achieve efficient impact analysis of LCD module.

In this study, when structure which is combined by welding is receiving internal pressure, finite element analysis to confirm stability of structure and reliability of welding part is achieved. And we analyze the results. Also, if stability of the structure and reliability of the welding part are not defined, research that look for method to change design to receive stability and reliability is achieved.

Satellite upper platform is optimized by response surface method which has non-gradient, semi-glogal, discrete and fast convergency characteristics. Sampling points are extracted by design of experiments using Central Composite Method and Factorial Design. Also response surface is generated by the various regression functions. Structure analysis is execuated with regard for static and dynamic environment in launching stage. As a result response surface method is superior to other optimization method with respect to optimum value and cost of computation time. Also a confidence is varified in the various regression models.

Recently, simulated annealing algorithms have widely been applied to many structural optimization problems. In this paper, simulated annealing, boltzmann annealing, fast annealing and adaptive simulated annealing are applied to optimization of truss structures for improvement quality of objective function and number of function evaluation. These algorithms are classified by cooling schedule. The authors have changed parameters of ASA's cooling schedule and the influence of cooling schedule parameters on structural optimization obtained is discussed. In addition, cooling schedule of BA and ASA mixed is applied to 10 bar-truss structure.

Computer is a very powerful machine which is widely using for data processing, DB construction, peripheral device control, image processing etc. Consequently, many researches and developments have progressed for high performance processing unit, and other devices. Especially, the core units such as semiconductor parts are rapidly growing so that high-integration, high-performance, microminiat turization is possible. The packaging in the semiconductor industry is very important technique to de determine the performance of the system that the semiconductor is used. In this paper, the inspection of the inner defects such as delamination, void, crack, etc. in the semiconductor packages is studied. ESPI which is a non-contact, non-destructive, and full-field inspection method is used for the inner defect inspection and its results are compared with that of C-Scan method.

A scratch tester was developed to evaluate the adhesive strength at interface between thin-film and substrate(silicon wafer). Under force control, the scratch tester can measure the normal and the tangential forces simultaneously as the probe tip of the equipment approaches to the interface between thin-film and substrate of wafer. The capacity of each component of force sensor is 0.1 N ${\sim}$ 100 N. In addition, the tester can detect the signal of elastic wave from AE sensor(frequency range of 900 kHz) attached to the probe tip and evaluate the bonding strength of interface. Using the developed scratch tester, the feasibility test was performed to evaluate the adhesive strength of thin-film.

The pantograph for Korean High Speed Train was developed by home-grown technology. In this study, a system to measure the contact force of pantograph is developed and installed on the prototype high speed train, Contact force prevents the pantograph from separating from the catenary. However, excessive contact force causes rapid erosion of catenary. The contact force can be divided into lift force and spring force. Contact force measurement is conducted while the train runs on the test track. The lift force is measured by the load cell on the roof separately and combined with the spring force of pan-head to form the contact force. Measured results show that the contact force of the pantograph of Korean High Speed Train is below the upper limit regulated by the high speed train standards. The contact force measuring system provides data to evaluate safety of the catenary system.

This paper presents an improved identification algorithm of active magnetic bearing rotor systems considering sensor and actuator dynamics. An AMB rotor system has both real and complex poles so that it is very hard to identify them together. In previous research, a linear transformation through a fictitious proportional feedback was used in order to shift the real poles close to the imaginary axis. However, the identification result highly depends on the fictitious feedback gain, and it is not easy to identify the additional dynamics including sensor and actuators at the same time. First, this paper discusses the necessity and a selection criterion of the fictitious feedback gain. An appropriate feedback gain minimizes dominant SVD(Singular Value Decomposition) error through maximizing rank deficiency. Second, more improvement in the identification is achieved through separating the common additional dynamics in all elements of frequency response matrix. The feasibility of the proposed identification algorithm is proved with two theoretical AMB rotor models. Finally, the proposed scheme is compared with previous identification methods using experimental data, and a great improvement in model quality and large amount of time saving can be achieved with the proposed method.

In this paper, fault tolerant mechanisms are presented for a servo manipulator system designed to operate in a hot cell. A hot cell is a sealed and shielded room to handle radioactive materials, and it is dangerous for people to work in the hot cell. So, remote operations are necessary to handle the radioactive materials in the hot cell. KAERI has developed a servo manipulator system to perform such remote operations. However, since electric components such as servo motors are weakened with radiation, fault tolerant mechanisms have to be considered. For fault tolerance of the servo manipulator system, hardware and software redundancy has been considered. In the case of hardware, radioactive resistant electric components such as cables and connectors have been adopted and motors driving a transport have been duplicated. In case of software, a reconfiguration algorithm accommodating one motor's failure has been developed. The algorithm uses redundant axes to recover the end effector's motion in spite of one motor's failure.

The senses that a doctor can feel is limited in MIS(Minimal Invasive Surgery) which guarantees the fast recovery of the patient and minimal incision for going in and out of instruments through the tissue of the patient. In particular, the surgical robotic teleoperation system developed recently serves with only the information of eyesight and auditory sense. Therefore force-reflection is the most demanded element of the senses in manipulating surgical instruments. In this paper, we designed the Master system and the 2 D.O.F grasper for the robotic teleoperation system(Slave) that has two force sensors on the grasper. Particularly, we focused on serve to master's handle with the contact force between tissue and the grasper of Slave.

Kinematic calibration enhances absolute accuracy by compensating for the fabrication tolerances and installation errors. Effectiveness of calibration procedures depends greatly on the measurements performed. This paper investigates identifiable parameters and optimum postures for four different calibration procedures - measuring postures completely with inverse kinematic residuals, measuring postures completely with forward kinematics residuals, measuring only the three position components, and restraining the mobility of the end-effector using a constraint link. The study is performed for a six degree-of-freedom fully parallel HexaSlide type parallel manipulator, HSM. Results verify that all parameters are identifiable with complete posture measurements. For the case of position measurements, one and for the case of constraint link, three parameters were found non-identifiable. Selecting postures for measurement is also an important issue for efficient calibration procedure. Typically, the condition number of the identification Jacobian is minimized to find optimum postures. Optimal postures showed the same trend of orienting themselves on the boundaries of the search space.

This paper presents a new family of 4-DoF parallel mechanism with two platforms. The new mechanism is composed of front and rear platforms, and three limbs. Two limbs with 6dof joint (P-P-S-P) are attached to the each platform and are perpendicular to baseplate, while the middle limb with 4-Dof joints (R-R-R-P or R-R-P-P) is attached to the revolute joint that connect front and rear platform. The two-DoF-driving mechanism at the middle limb with two base-fixed prismatic actuators can generate the heaving and roll motions or two translational motions. Therefore, Therefore, the new 4-Dof parallel mechanism (1T-3R) can generate pitch motions at each platforms, roll, and heaving motions, while another type of new 4-Dof parallel mechanism (2T-2R) can generate pitch motions at each platforms, x and z translational motions. For 1T-3R mechanism, kinematic analyses including inverse, forward kinematics, and Jacobian are performed.

Recently, higher and bigger building is trend of the construction. Accordingly the building material is getting bigger and construction equipments are developed. But operation is still depends on human resource. Therefore there have several problems that are safety, laborious operation, and shortage of worker. In the various construction sites, the automation in construction is introduced to solve these problems. This paper proposed the automation system in construction that installs curtain wall in a high building. The system is expected effects that are reduction of a construction period, retrenchment of the cost and assurance of safety.

Kinematic calibration is a process whereby the actual values of geometric parameters are estimated so as to minimize the error in absolute positioning. Measuring all components of Cartesian posture, particularly the orientation, can be difficult. With partial pose measurements, all parameters may not be identifiable. This paper proposes a new device that can be used to identify all kinematic parameters with partial pose measurements. Study is performed for a six degree-of-freedom fully parallel Hexa Slide manipulator. The device, however, is general and can be used for other parallel manipulators. The proposed device consists of a link with U joints on both sides and is equipped with a rotary sensor and a biaxial inclinometer. When attached between the base and the mobile platform, the device restricts the end-effector's motion to five degree-of-freedom and can measure position of the end-effector and one of its rotations. Numerical analyses of the identification Jacobian reveal that all parameters are identifiable. Computer simulations show that the identification is robust for the errors in the initial guess and the measurement noise.

In order to avoid the unit inconsistency problem in the conventional Jacobian matrix, previously we presented new formulation of a dimensionally homogeneous inverse Jacobian matrix for parallel manipulators with a planar mobile platform by using three end-effector points based on the velocity relationship [1]. This paper presents force relationships between joint forces and Cartesian forces at the three End-Effector points. The derived force relationships can then be used for analyses of the input/output force transmission. These analyses, forward and inverse force transmission analyses, depend on the singular values of the derived dimensionally homogeneous Jacobian matrix. Using the proposed force relationship, a numerical example is presented for actuator size design of a 3-RRR planar parallel manipulator.

This paper describes an automatic weld seam tracking method of plasma arc welding system designed for the corrugation panel that consists of a linear section and a curved section with various curvatures. Due to the complexity of the panel shape, it is difficult to find a seam and operate a torch manually in the welding process. So, the laser vision sensor for seam tracking is designed for sensing the seam position and controlling a torch automatically. To achieve precise seam tracking, the design of sensor head, image simulation, and calibration are carried out. Through a series of experiment result, compensation algorithm is added and real time error compensation is achieved. The experiment result shows that this vision sensor works effectively. It will provide more precise welding performance and convenience to the operator.

This paper presents the calibration for the parallel typed tilting table. The calibration system needs only simple sensing device which is a digital indicator to measure the orientation of a table. The calibration algorithm is developed by a measurement operator. It eliminates the concern about the poor parameter observability due to a large number of parameters of parallel-mechanism. This paper uses the QR-decomposition to find the optimal calibration configurations maximizing the linear independence of rows of a observation matrix. The number of identifiable parameters is examined by the rank of the observation matrix, which represents the parameter observability. The method is applied to a Parallel-typed Tilting Table and all the necessary kinematic parameters are identifiable.

Automation of roadway sign painting offers more safety for workers, shortening of work period, etc., compared with manual painting. In this study, an automated system using a gantry-type robot was developed for roadway sign painting which has been done manually. Any characters (Korean and English) as well as symbols can be painted by the system. A simulator was also developed, which can show the painted results in advance. The developed system performed well, and the signs painted by the system were found to be as accurate as those made by the simulator.

In this paper deals with a unique method for measuring vehicle states such as body sideslip angle and tire sideslip angle using GPS velocity information in conjunction with other sensors. A method for integrating Inertial Navigation System (INS) sensors with GPS measurements to provide higher update rate estimates of the vehicle states is presented, and the method can be used to estimate the tire cornering stiffness. The experimental results for the GPS velocity-based sideslip angle measurement and cornering stiffness estimates are compared with the theoretical predictions. From the experimental results, it can be concluded that the proposed method has an advantage for future implementation in a vehicle safety system.

The automatic load transfer switch (ALTS) typically automatically transfers electrical loads from a normal electrical power source to an emergency electrical power source upon reduction or loss of normal power source voltage. It can also automatically re-transfer the load to the normal power source when the normal voltage has been restored within acceptable limits. The transfer operation of ALTS is accomplished by a spring-driven linkage mechanism. In this paper we build a dynamic model of driving mechanism for ALTS using ADAMS and checked the characteristics of the transfer operation. Finally we performed a detailed design of the driving mechanism through results of analysis and confirmed it to satisfy design requirements.

The model of azimuth driving servo system with a flexible antenna in a microwave seeker has been derived in this paper. The validity of the model is verified by comparing the result of the model with that of experiment. It is found that one should consider an antenna as a flexible body in case of modeling the dynamics of the microwave seeker. It is also known that the effect of reducing backlash magnitude for extending the bandwidth in the system with a flexible antenna is smaller than that in the system with a stiff antenna. It is thought that the model-based design optimization of the microwave seeker will become possible by virtue of the derived model, when a weight reduction and a bandwidth extension are required.

Air entrainment can become a significant problem in a web handling process. The development of air film between a web and an idle roller can cause a reduction of traction and traction coefficient, by which a slip is occurrred. Computational and experimental study was carried out to describe the slippage of an idle roller for given operating conditions, tension and web velocity. An extended mathematical model to find out a slip condition was developed by using the models of air film height, dynamic traction coefficient, and torque balance of a rotational roller. And by using the extended model, a mechanism to define the slippage between the roller and the moving web was suggested. The results of simulation and experiment showed that the extended dynamic model could properly characterize the rotational motion of the idle roller by considering dynamic traction coefficient. By examining the rotational motion of the idle roller with web dynamics(speed), the mechanism to define al slip condition between the roller and the web was found to be effective.

The design of a high speed axial piston pump for hydrostatic transmission systems requires specific understanding on where and how much its internal frictional and flow losses are generated. In this study, the frictional loss of a bentaxis type hydraulic piston pump was analyzed in order to find out which design factors influence the mechanical efficiency most significantly. To this end, the friction coefficients of the sliding components were experimentally identified by a specially constructed tribometer. Applying them to the three-dimensional dynamic model of the pump presented by Doh and Hong [1], the friction torques generated by the sliding components such as piston head , bearing and valve plate were theoretically computed. The accuracy of the computed results was confirmed by the comparison with the experimentally measured mechanical efficiency. In this paper, it is shown that the viscous friction on the valve plate and the drive shaft bearing is the primary sources of the frictional losses of the bent-axis type pump, while the friction forces on the piston contribute to them only slightly.

This study validated the musculoskeletal model of the human lower extremity by comparative study between calculated muscle parameters through simulation using modified hill-type model and measured them through isokinetic exercise. And the relationship between muscle forces and moments participated in motion was quantified from the results of simulation. For simulation of isokinetic motion, a three-dimensional anatomical knee model was constructed using trials of gait analysis and the EMG-force model was used to determine muscle activation level exciting muscles. The modified Hill-type model was used to calculate individual muscle forces and moments in dynmaic analysis and the results were validated by comparing them of experiments on BIODEX. The results showed that there was a high correlation between calculated torques from simulation and measured them from experiments for isokinetic motion(R=0.97). Therefore we concluded that the simulation by using musculoskeletal model was so useful means to predict and convalesce musculoskeletal-related diseases, and analyze unrealizable experiment such as clash condition.

Tapping mode atomic force microscopy (TM-AFM) utilizes the dynamic response of a resonating probe tip as it approaches and retracts from a sample to measure the topography and material properties of a nanostructure. We present recent results based on nonlinear dynamical systems theory, computational continuation techniques and detailed experiments that yield new perspectives and insight into AFM. A dynamic model including van der Waals and Derjaguin-Muller-Toporov (DMT) contact forces demonstrates that periodic solutions can be represented with respect to the approach distance and excitation frequency. Turning points on the surface lead to hysteretic amplitude jumps as the tip nears/retracts from the sample. Experiments are performed using a tapping mode tip on a graphite sample to verify the predictions.

PGO(Powered Gait Orthosis) mounted with pneumatic muscle as an actuator is upgraded model from RGO(Replicate Gait Orthosis) for paraplegia patients to walk easy and safe. Pneumatic muscles supply powers to both hip joint during PGO gait. The objective of this research is to develop the PGO gait simulation model. Dynamic model of PGO linkage system is processed. Mathematical model of pneumatic muscle was developed and combined it with PGO linkage system. Developed simulation model will be used as a tool for evaluation of the efficiency of pneumatic muscle and for analysis the PGO system.

The powershift transmission have the advantage of easier operation and higher efficiency by using the hydraulic clutch and mechanical power transfer system. It is important to control the engaging pressure and time. The hydraulic control system is used for these controls of the modulator valve, the accumulator, the sump valve and etc. This study have made a simulator for verifying the pressure characteristic of the shuttle powershift transmission and developed the computer simulation model of the hydraulic components and system by using 'AMESim'. As a result, the design parameters which have an effect on the pressure modulation are verified to the spring stiffness of the modulator valve and the volume of the accumulator.

Time-Optimal solutions for attitude control with reaction wheels as well as with thrusters are studied. The suggested varying-time-sharing ratio thrusting is found to reduce the maneuvering time enormously. The hybrid control such as sequential hybrid and simultaneous hybrid with reaction wheels and thrusters are considered. The results show that simultaneous hybrid method reduces the maneuver time very much. Spacecraft model is KOrea Multi-Purpose SATellite(KOMPSAT)-II, which is being developed by KARI in KOREA as an agile maneuvering satellite.

Simple and effectively developed learning control logic is used to control vibration of U type Tuned Liquid Damper system. The purpose of this paper is design optimal control system to deal with unknown errors from nonlinearity and variation that cost modeling difficulty in complex structure and is followed with the desired behavior. Finally this hybrid control method applied to U type Tuned Liquid Damper structure gives the benefit from better performance of precision and stability of the structure by reducing vibration effect. This research leads to safety design in various structure to robust unspecified foreign disturbances such as earthquake.

Direct drive servovalve(DDV) is a kind of one-stage valve since the rotary motion of DC motor is directly transferred to the linear motion of valve spool through the link. Since the structure of DDV is simple, it is less expensive, more reliable and offers reduced internal leakage and reduced sensitivity to fluid contamination. However, the flow force effect on the spool motion is significant such that it induces large steady-state error in a step response. If the proportional control gain is increased to reduce the steady-state error, the system becomes unstable. In order to satisfy the system design requirements, the lead-lag controller is designed using the complex method that is one kind of constrained direct search method.

As the internet communication is prevalent in recent years, it becomes quite possible to monitor and control some mechanical plants from the remote place through the TCP/IP communication. Such a concept is expected to be applied to many industrial systems for easy maintenance and trouble shooting as well as various kinds of expensive test equipments for sharing. In this research, remote data monitoring and speed control for a DC-motor is implemented and tested through TCP/IP communication with embedded micro-controllers. It showed the possibility of reliable remote control system design utilizing the internet communication.

Vibration causes noise and sometimes makes structure unstable. Especially, due to the efforts of lightening, deformation of flexible structure is increased in its shape. Just a little disturbance can cause vibration and low damping ratio makes residual vibration last long time. This research is concerned with the model reference neuro-controller design for the vibration suppression of smart structures. By using a model reference neurocontroller, which is one of the algorithms of adaptive control, we performed an adaptive control of flexible cantilever plate and opened box structure with piezoelectric materials. The proposed adaptive vibration control algorithm, a model reference neuro-controller, was proved in its effectiveness by applying to an opened box structure. The model reference neuro-controller is implemented with DSP, and the real-time adaptive vibration control experiment results confirm that the model reference neuro-controller is reliable.

In the semiconductor and the optical industry, a new transport system which can replace the conventional sliding system is required. These systems are driven by the magnetic field and the conveyer belt. The magnetic field damages semiconductor and contact force scratches the optical lens. The ultrasonic wave driven system can solve these problems. In this paper, the object transport system using the excitation of ultrasonic wave is proposed. The experiments for finding the optimal excitation frequency, finding phasedifference between two ultrasonic wave generators are performed. The relationship of transporting speed according to the change of flexural beam shape is verified. In addition, the system performance for practical use is evaluated.

In this study, position and force simultaneous trajectory tracking control system with pneumatic cylinder driving apparatus is proposed. The pneumatic cylinder driving apparatus that consists of two pneumatic cylinders constrained in series and two proportional flow control valves offers a considerable advantage as to non-interaction of the actuators because of the low stiffness of the pneumatic cylinders. The controller applied to the driving system is composed of a non-interaction controller to compensate for interaction of two cylinders and a disturbance observer to reduce the effect of model discrepancy of the driving system in the low frequency range that cannot be suppressed by the non-interaction controller. The experimental results with the proposed control system show that the interacting effects of two cylinders are eliminated remarkably and the proposed control system tracks the given position and force trajectories accurately.

Recently, the height of building is more high and the number of moving people is increasing in the building. So we want to be necessary more effective vertical moving tool. Most of high intelligent building can satisfy tills need using many elevators. Many elevators system should need to distribute and distribute many elevator effectively. This paper effectively use many elevators to reduce customer' waiting time and propose the model of mathematical group control system.

An iterative modal analysis approach is developed to determine the effect of the transverse open cracks and the moving mass on the dynamic behavior of simply supported pipe conveying fluid. The equation of motion is derived by using Lagrange's equation. The influences of the velocity of moving mass, the velocity of fluid flow and a crack have been studied on the dynamic behavior of a simply supported pipe system by numerical method. The crack section is represented by a local flexibility matrix connecting two undamaged beam segments. that is, the crack is modelled as a rotational spring. Totally, as the velocity of fluid flow is increased, the mid-span deflection of simply supported pipe conveying fluid is increased. The position of the crack is middle point of the pipe, the mid-span deflection of simply supported pipe presents maximum deflection.

The purpose of this study is to analyze the phenomena of the thermally-induced vibration for the flexible space structure due to abrupt change of radiation heating circumstance using the numerical analyze and experiment test. In order to verify this structure, numerical approaches on the simplified flexible tube were compared with experimental test results at the ground experimental facility In this analyze, it was found that the thermal deformation occurs firstly due to fast radiation heating of flexible structure and then the thermally-induced vibration would be induced due to small periodic change of temperature. According to comparison of numerical and experimental result, in case of no tip mass, the first mode vibration by the numerical analyze was O.78Hz same as that of the experimental result However in case of increase tip-masses of 8g l6g, 50g and 100g, the first modes vibration theoretical analyze were 1.75Hz, 1.3Hz, 0.87Hz and O.73Hz, in decrease trend respectively and those by experimental test were 234Hz, 1.5Hz, O.78Hz and O.78Hz in decrease trend respectively Although using the simpled equation for the estimation, the estimation results were similar to experimental results.

Korea Gas Corporation(KOGAS) is a Liquified Natural Gas(LNG) supplier through out the Korea. LNG, which is imported wholly from foreign countries, is compressed 1/600 for easy transportation and is stored in a liquid state in the storage tanks at Incheon, Pyeongtaek and Tongyeong. At LNG receiving terminals, LNG is vaporized to natural gas before supplying to City Gas Consumer of Power Plant. The secondary pump is a equipment which compress LNG from 1- kgf/cm2 to 70 kgf/cm2. The secondary pump at Tongyeong LNG receiving terminal is consisted of two pumps in one underground PIT, and is connected to supporting structures. It is therefore expected that there is a vibration problem whit the pump and was found that high level vibration was occurred in a low frequency band($5^{\sim}10Hz$). In this paper, the vibration of secondary pump was analyzed, and the main cause of vibration was found out.

The market of the horizontal axis washing machine (drum washing machine) has been growing drastically in Korea by about 80% annually since 2000. As market grows fast, the customerTs demands concerning quality becomes more strict and various. Imbalance sensing is a key technology to reduce the NVH problem in a washing machine, because the laundry is time-variant and uncontrollable source of imbalance, which can cause more than 200kgf exciting force. In this paper, imbalance-sensing methods are briefly reviewed, new acceleration sensing circuits are examined, and finally the control algorithm of spinning process is proposed and validated.

This paper predicts the modified proportional damping structural eigenvectors and eigenvalues due to the change in the mass and stiffness of a proportional damping structure by iterative calculation of the sensitivity coefficient using the original dynamic characteristics. The method is applied to examples of a cantilever and 3 degree of freedom lumped mass model by modifing the mass and stiffness. The predicted dynamic characteristics are in good agreement with these from the structural reanalysis using the modified mass and stiffness.

This paper introduces a fast Fourier transform (FFT)-based spectral analysis method for the transient responses as well as the steady-state responses of linear discrete systems. The force vibration of a viscously damped three-DOF system is considered as the illustrative numerical example. The proposed spectral analysis method is evaluated by comparing with the exact analytical solutions as well as with the numerical solutions obtained by the Runge-Kutta method.

On this study, we improved the efficiency applying algorithm that is repeatedly using orthogonal array in discrete design space and filling a defect of gradient method in continuous design space. we showed optimal ply angle that maximized 1st natural frequency of CFRP laminated composite cantilever plate by each aspect ratio. A finite element analysis on the CFRP laminated composite cantilever plate using orthogonal array is carried out, and the results are compared with those obtained by modal testing.

Exact dynamic stiffness model for a uniform straight pipeline conveying unsteady fluid is formulated from a set of fully coupled pipe-dynamic equations of motion, in which the fluid pressure and velocity of internal flow as well as the transverse and axial displacements of the pipeline are all treated as dependent variables. The accuracy of the dynamic stiffness model formulated herein is first verified by comparing its solutions with those obtained by the conventional finite element model. The spectral element analysis based on the present dynamic stiffness model is then conducted to investigate the effects of fluid parameters on the dynamics and stability of an example pipeline problem.

In this paper, the spectral element model is derived for the vibration and stability analyses of an axially moving viscoelastic beam subjected to axial tension. The viscoelastic material is represented by using a one-dimensional constitutive equation of hereditary integral type. The accuracy of the present spectral element model is first verified by comparing the eigenvalues obtained by the present spectral element model-based SEM with those obtained by the exact theory and the conventional FEM. The effects of viscoelasticity on the vibration and stability of an example moving viscoelastic beam are numerically investigated.

Some characteristics of the sensitivities of the eigenvalues for beams have been found in the paper. For cantilever beams and simply supported beams, the sensitivities of the eigenvalues to the stiffness correction factor of one element are equal and opposite to the sensitivities to the mass correction factor of the symmetrically positioned element. The relationship means that to increase stiffness in one element has the same effects on the eigenvalues as to decrease mass by the same proportion in the symmetrically positioned element. For beams with other boundary conditions, however, the relationship does not hold.

This paper concerns the analytical modeling and dynamic analysis of advanced cantilevered blade structure implemented by a dual approach based on structural tailoring and viscoelastic materials technology. Whereas structural tailoring uses the directionality properties of advanced composite materials, the passive materials technology exploits the damping capabilities of viscoelastic material(VEM) embedded into the host structure. The structure is modeled as a composite thin-walled beam incorporating a number of nonclassical features such as transverse shear, secondary warping, anisotropy of constituent materials, and rotary inertias. The case of VEM spreaded over the entire span of the structure is considered. The displayed numerical results provide a comprehensive picture of the synergisitic implications of the application of both techniques, namely, the tailoring and damping technology on vibration response of thin-walled beam structure exposed to external time-dependent excitations.

In order to investigate the mechanism of brake squeal noise, the sound and vibration of an actua1 brake system were measured using a brake dynamometer. The experimental results show that disc run-out varies with brake line pressure and the factor of squeal generation is the run-out due to the misalignment of brake disk. A three degrees of freedom friction model is developed for the disk brake system where the run-out effect and nonlinear friction characteristic are considered. The results of numerical analysis of the model agree well with the experimental results. Also, the stability analysis of the model was performed to predict the generation of brake squeal due to the design parameter modification of brake systems. The results show that the squeal generation depends on the nm-out rather than the friction characteristic between the pad and the disk of brake.

This paper deals with a vibration control analysis of a rotating composite blade, modeled as a tapered thinwalled beam induced by heat flux. The displayed results reveal that the thermally induced vibration yields a detrimental repercussions upon their dynamic responses. The blade consists of host graphite epoxy laminate with surface and spanwise distributed transversely isotropic (PZT-4) sensors and actuators. The controller is implemented via the negative velocity and displacement feedback control methodology, which prove to overcome the deleterious effect associated with the thermally induced vibration. The structure is modeled as a composite thin-walled beam incorporating a number of nonclassical features such as transverse shear, secondary warping, anisotropy of constituent materials, and rotary inertias.

Absorptive material arrangement method for effective interior noise control is proposed. Sound field with arbitrary boundary condition is formulated by Kirchhoff-Helmholtz integral equation. A simple example such as a rectangular cavity will present physical meaning between changing boundary condition and control of sound field. The effect of changing boundary condition is expressed in modal admittance. From this formulation, an admittance map is presented. The admittance map is the figure to represent position where absorptive material is attached. The admittance map can be assigned to each resonant frequency. There, however, may be common area of those maps. Then, frequency robust arrangement of absorptive material in noise control will be presented.

The acceptance test of KTX has been performed in Korea. During the test, lateral vibration of carbody over the accepted value called swat was found. KTX has 20 car trainsed formation whose trailer cars are linked by articulate bogies. So this study is performed to see the effects of long trainsed formation on vehicle dynamics and the train stability by 20 car vehicle model. Firstly the reliable vehicle model which shows well the tendencies appeared in the tests on the high speed test line is required to find the cause of lateral vibration and the countermeasure. Vehicle model was made for the analysis with VAMPIRE. The analysis results show that secondary air spring lateral stiffness is the most significant parameter to cause carbody lateral vibration. Mode analysis results show that he least damped mode shape is similar to the vibration pattern shown in the tests that the amplitude of the motion increases along the train set and decreases in the tail part. For the case of short train formation with 7 or 10cars, sway does not happen. But in the case of longer train formation with 16 or 20 cars, sway was found.

ECCOMAS Thematic Conference Multibody 2003 was held at IST (Instituto Superior Technico), Lisbon, Portugal from July 1 to July 4. 2003. And MBDV(Multibody Dynamics and Vibration) in the 2003 ASME DETC was held at Chicago, U.S.A. from September 2 to September 6. In this paper, the presented papers in these conferences were reviewed and the trends in the multibody dynamics are summarized. The session titles in these conferences include Flexible Multibody Dynamics, Vehicle Dynamics, Contact, Biomechanics, Real-time Challenges, Spatial manipulator and Control, Multidisciplinary Applications, and Advanced Education. The poster session was also organized for more discussions in the Multibody2003 conference.

This paper presents a computational method for deformation modes of a flexible body in multibody system from the experimental modal analysis and an efficient method for flexible multibody dynamic analysis by use of the modes. It is difficult to directly use experimental modal parameters in flexible multibody dynamic analysis. The major reasons are that there are many inconsistencies between experimental and analytical modal data and experimental noises are inherent in the experimental data. So two methods, such as, a method for ortho-normalization of experimental modes and the other one for mode expansion, are suggested to gain deformation modes of a flexible body from the experimental modal parameters. Using the virtual work principle, the equation of motion of a flexible body is derived. The effectiveness of the proposed method will be verified in the numerical example of cab vibration of a truck by comparing analysis and test results.

Two different subsystem synthesis methods with independent generalized coordinates have been developed and compared. In each formulation, the subsystem equations of motion are generated in terms of independent generalized coordinates. The first formulation is based on the relative Cartesian coordinates with respect to moving subsystem base (virtual) body. The second formulation is based on the relative joint coordinates using recursive formulation. Computational efficiency of the formulations has been compared theoretically by the operational counting method.

In this paper, the absolute nodal coordinate formulation was introduced to describe the large deformation problems. And also, the modal coordinates were employed to represent the small elastic deformation. A new hybrid formulation was developed to combine the modal coordinates and the absolute nodal coordinates. A spherical joint and the DOT1 constraint were developed to carry out the numerical simulation of mechanical systems with kinematic joints. A beam example was suggested to show the new formulation. The simulation results using the modal coordinates and the absolute nodal coordinates show a good agreement to the experiments.

There have been various studies to find out the correlation between subjective and objective evaluation ride comfort in vehicle. To be proper, subjective evaluation considers the signal pattern of vibration generated by vehicle. In this study, driving test is carried out to extract the correlation between subjective evaluation by subjects and vibration signals measured on the seat when a test vehicle is passing over test bumpers, which have two sorts of width change and height change respectively. A unique reference bumper is used to be compared with the relative comfort of test bumpers.

Effects of electromagnetic force which is one of the most important factor of metal transfer that affects bead geometry and microstructure of weld metal in GMAW(gas metal arc welding). In this paper, different ways of external electromagnetic forces were applied on GMAW process and their effects on the welding were studied. On certain conditions, better bead geometry, better influence on the arc and metal transfer mode and higher welding efficiency could be obtained. Experimental methods and their results will be presented.

The paper shows properties such as vickers hardness, yielding and ultimate stresses for the weld zone of the butt and the lap jointed specimens, and compare maximum loads, stress-strain curves, deformation appearance after guided bending test and fracture appearance for butt and lap jointed specimens. The research in this experiment also shows the weldability of the butt joint specimen is better than that of the lap joint specimen using FSW with $2mm^{t}$ aluminum alloy sheet in milling machine.

Welding automation in shipbuilding process, especially in the assembly line is considered to be a difficult job because welding part is too huge , various and unstructured for a welding robot to weld the whole part automatically. We developed an automatic welding robot to improve those difficult process. This paper show how to systematically operate the integrated automation system which consists of several robots. We introduce our software and system integration method. Specially we focus that network communication and operating process. The developed system visualizes the operation environment using Open Inventor and communicates with the entire system via TCP/IP and FTP.

For coping with various ship types easily, welding automation using CAD/CAM is demanded and and developed. In this paper we propose a dedicated CAM system which generates MACRO files to control robot manipulators. The paper contains CAD interface, virtual simulation, macro generation and job schedulling. At first, it defines and extracts weldline from CAD data, and generates proper MACRO programs. And it obtains optimum job schedules. This system removes the manual work, and consequently reduces the overall lead time. And it reduces costs and time for developing robot welding programs. This can be expanded to virtual factory simulation technology. Moreover, it is possible to apply this system for automation of Cutting and Painting

Current trend of design technologies shows engineers to objectify or automate the given decision-making process. The numerical optimization is an example of such technologies. However, in numerical optimization, the uncertainties are uncontrollable to efficiently objectify or automate the process. To better manage these uncertainties, Taguchi method, reliability-based optimization and robust optimization are being used. Based on the independence axiom of axiomatic design theory that illustrates the relationship between desired specifications and design parameters, the designs can be classified into three types: uncoupled, decoupled and coupled. To best approach the target performance with the maximum robustness is one of the main functional requirements of a mechanical system. Most engineering designs are pertaining to either coupled or decoupled ones, but these designs cannot currently accomplish a real robustness thus a trade-off between performance and robustness has to be made. In this research, the game theory will be applied to optimize the trade-off.

Reverse Engineering has wasted high cost and much of time because there were few softwares for supporting post-process of 3D scanning efficiently. However today Reverse Engineering is attracted by various researchers on a count of its practical use and gradually improved quality. Accordingly, many industries induce the technology of Reverse Engineering. For keeping with this general trend, the algorithms of generating meshes from 3D scanning data and visualizing their cross section are suggested in this research. It is expected that output data from the algorithms can be used in inspecting errors between product designs and their final products, and in finding a way to improve quality by considering ideal model data.

A hybrid finite element analysis was used to analyze the influence of ring gear rim thickness and spline number on the static properties of an epicyclic gear system with manufacturing errors. Both of these parameters affected the bearing force and critical stress. The effect of changes in the rim thickness on the load sharing between the gears depended on the type of manufacturing error. Ring flexibility improved the load sharing between planetary gears only in systems with planet tooth thickness or planet tangential errors; for other types of error, ring flexibility worsened the load sharing. To improve load sharing, rim thickness and spline number should be controlled within a specific range. The effect of the ring gear boundary condition was more apparent in a system with errors than in a normal system.

This research develops a robot as the device which constructs underwater harbour. This construction is to build a breakwater, which is dangerous and difficult. The hydraulic parallel mechanism-typed robot is developed to mechanize the construction by operating of a stoneworker (or diver) through a joystick. The six-dof robot is able to carry 2-3 tons' heavy stone and put it on the surface of breakwater. This paper presents the mechanical design of a miniature robot, its control and application for the breakwater construction.

Recently with the development of semiconductor technology, the miniaturization of parts and products as well as their high precision is required. In addition, as the national competitiveness is increasingly affected by the development of the micro parts through micro machining technology, the study of the micro machining technology is being conducted in many countries. The goal of this study is to fabricate micro tools under the size of $20{\mu}m$ and to machine micro holes using them. The fabrication is done by grinding and the application of ELID to the grinding wheel. The surface roughness of the micro tools is measured to evaluate the study.

After the micro turning lathe was developed in the last year by AMR Laboratory in Yeungnam university, a micromilling machine is developed for micro/meso machining. This machine is integrated with PZT-driven micro-sliders, micro-linear encoders, aerostatic spindle which has maximum 150,000 rpm. It is applicable to milling and drilling of micro scale. This paper presents the possibility of micro/meso machining and characteristics of micro end milling process by using micro machine. A machining of micro parts using 0.2 mm flat end mill was achieved by micro-milling machine. Experimental results show the machining capability and positional accuracy of this machine is good enough for machining micro parts.

Surface and edge finishing processes are important technological operations of in parts machining. Quality of the finished parts directly affect the performance of the whole product. Especially, edge quality, which depends on burr control, is extremely important. Burrs are undesirable projections of the material beyond the edge of the workpiece. A number of deburring processes have been developed such as barreling, brushing, chemical methods etc. But, there are only few publications in the area of applying ultrasonics to deburring. When ultrasonic vibration propagates in the liquid medium, a large number of bubbles are formed. These bubbles generate an extremely strong force, which can be used to remove burrs. Cavitation is used as a term to describe the erosion of parts caused by the action of cavities in liquid. The object of this study is to analyze the effects of ultrasonic cavitation in the deburring process. For this purpose, we introduce a new ultrasonic cavitation method, which efficiently removes the burrs. Experimental parameters to verify the deburring effects of ultrasonic cavitations are ultrasonic power, amplitude, distant of the transducer from the workpiece, deburring time and abrasive. It has been shown that deburring with ultrasonic cavitation in water is effective to burrs.

The purpose of this study is to clarify the bulk/sheet forming characteristics of bulk amorphous alloys in the supercooled liquid state. The temperature dependences of Newtonian viscosities of amorphous materials are obtained based on the previous experimental works. Finite element analyses for compression forming and sheet deep drawing of amorphous materials are performed. Effects of friction coefficients and temperature are examined and formability of amorphous material is explained in detail.

Tube hydroforming provides a number of advantages over conventional stamping process, including fewer secondary operation, weight reduction, assembly simplification, adaptability to forming of complex structural components and improved structural strength and stiffness. In this study, the effect of the heat treatment on the hydro-formability has been investigated. By using the mild steel tube bulging test is performed at various heat treatment conditions to evaluate the hydro-formability.

The sapphire wafer for blue light emitting device was manufactured by the implementation of the chemical and mechanical polishing process. The surface polishing of crystalline sapphire wafer was characterized by zeta potential measurement. The reduction process with the alkali slurry provides the surface chemical reaction with sapphire atoms. The poly-urethane pad also provides the frictional force to take out the chemically-reacted surface layers. The surface roughness was measured by the atomic force microscope and the crystalline quality was characterized by the double crystal X -ray diffraction analysis.

End-milling have been used widely in industrial system because it is effective to a material manufacturing with various shape. Recently the end-milling processing is needed the high-precise technique with good surface roughness and rapid time in precision machine part and electronic part. The optimum surface roughness has an effect on end-milling condition such as, cutting direction, spindle speed, feed rate and depth of cut, etc. Therefore this study was carried out to presume for mutual relation of end-milling condition to get the optimum surface roughness by regression analysis. The results shown that coefficient of determination($R^{2}m$) of regression equation has a fine reliability over 80% and regression equation of surface rough is made by regression analysis.

In a foaming process of microcellular plastics (MCPs) with a batch process, amorphous plastics and crystalline plastics have different characteristics for a foaming temperature. It is known that a foaming of amorphous plastics occurs at the temperature above a glass transition temperature, however, it is discovered that crystalline plastics do not take place above a glass transition temperature without exception, and even though the foaming occurs, it does not in all the range. In this research, to measure foaming temperature of crystalline polymer, a foaming experiment was performed using one of the typical crystalline polymer, polypropylene. To analyze whether the foaming occurs both at amorphous and crystalline fields, SEM was applied

Temperature and velocity distributions of hot air flows in rotational molding machines with two different shapes and structures of oven and inlet were investigated by using FLUENT, a commercial computational fluid dynamics code. The shape and structure of oven and inlet in current rotational molding machine were improved. Two different sizes of mold inside each oven were considered in the analysis. Temperature and velocity distributions of hot air flows in two different rotational molding machines were compared to each other. In order to reduce cycle time and improve product quality in current rotational molding machine, the improved shape and structure of oven and inlet were proposed.

Collaborative product commerce (CPC) supports a collaboration that a global enterprise and customer related to life cycle of product share product information and a collaboration process for the collaboration, and integrating applications. In this paper, we use common data schema in order to solve a interoperability problem about shared product information between enterprises. And we map to common data schema from each other different data format. Therefore we implement CPC Adaptor in order to integrate distributed product information.

Approximate optimization has become popular in engineering field such as MDO and Crash analysis which is time consuming. To accomplish efficient approximate optimization, accuracy of approximate model is very important. As surrogate model, Kriging have been widely used approximating highly nonlinear system . Because Kriging employs interpolation method, it is adequate for deterministic computer simulation. Because there are no random errors and measurement errors in deterministic computer simulation, instead of classical DOE ,space filling experiment design which fills uniformly design space should be applied. In this work, various space filling designs such as maximin distance design, maximum entropy design are reviewed. And new design improving maximum entropy design is suggested and compared.

This study is concerned about the development of three-dimensional bending machine for heat exchanger. Recently, three-dimensional bending is required for various heat exchanger. The purpose of this study is design of three-dimensional bending machine by analysis of bending process and structural analysis simulation. The analysis is carried out by FEM simulation using DEFORM and CATIA V5 software. The copper-tube is modeled by shell elements and the machine is modeled by placing proper shell and solid finite elements and fictitious mass properties to represent the real one. The final results of analysis are applied to the design of three-dimensional bending machine and the machine is successfully developed.

The properties of a machine tool greatly affect machining quality since a machine tool has large variance in its features. Machine tool makers want to find best machining condition with the one that they have built. Machine builders need to develop test specimen since it helps finding characteristics of machine tools when the machining properties of the specimen are analyzed. This paper develops test specimen to identify features of the main spindle, the feeding device, and the frame of a machine tool. The specimen is machined with a high speed machine and the features of the machine are analyzed with test items. They are surface roughness, overshoot in axial movement, errors in circular movement, feeding with small movement, and compensational error. This work can improve usability for a machine tool in machining practice.

A novel piezoelectric micropump using active check valves in place of conventional passive check valves in inlet and outlet has been proposed and investigated. It actively controls open/close motion of check valves using piezoelectric actuator for expansion/contraction of pump chamber. In this paper, bi-directional flow characteristics and load characteristics are experimentally investigated using an adequate timing control for valve closing motion with a prototype micropump fabricated with the effective size of $17{\times}8{\times}11mm^{3}$. From the experimental results, it is ascertained that optimal values of phase shift against voltage to drive pump chamber for realization of a miniaturized but powerful micropump, are $15^{\circ}$ in inlet check valve and $195^{\circ}$ in outlet. Based on the obtained results, a sheet-type active shuttle valve that has a unified valve-body for inlet and outlet check valves is proposed. A micropump with the effective size of $10{\times}10{\times}10mm^{3}$ is fabricated and basic characteristics are experimentally investigated.

Optimization has been used in many engineering problems and must be one of the essential components during design process now. It is the process of maximizing the performance called an objective function of a system while satisfying some constraints, so finite element method is generally required in order to obtain these values during optimization. However, there are some difficulties to obtain them by means of FEM, where the changes of design variables cause the distortion and the regeneration of mesh that may result in inaccuracy and inefficiency. In order to overcome this problem, this paper proposed an alternative that the eXtended FEM introduced and developed by Ted Belytschko was applied to the optimization process because the key points of the X-FEM lie in that the discontinuity can be represented independently on the mesh by a function called in an enrichment function.

When a projectile travels at high speed underwater, supercavitating flow arises, in which a huge cavity is generated behind the projectile so that only the nose, i.e., the cavitator, of the projectile is wetted, while the rest of it should be surrounded by the cavity. In that case, the projectile can achieve very high speed due to the reduced drag. Furthermore if the nose of the body is shaped properly, the attendant pressure drag can be maintained at a very low value, so that the overall drag is also reduced dramatically. In this study, shape optimization technique is employed to determine the optimum cavitator shape for minimum drag, given certain operating conditions. Shape optimization technique is also used to solve the potential flow problem for any given cavitator, which is a free boundary value problem having the cavity shape as unknown a priori. Analytical sensitivities are derived for various shape parameters in order to implement a gradient-based optimization algorithm. Simultaneous optimization technique is proposed for efficient cavitator shape optimization, in which the cavity and cavitator shape are determined in a single optimization routine.

The response surface method (RSM) became one of famous meta modeling techniques, however its approximation errors give designers several restrictions. Classical RSM uses the least squares method (LSM) to find the best fitting approximation models from the all given data. This paper discusses how to construct RSM efficiently and accurately using moving least squares method (MLSM) with sensitivity information. In this method, several parameters should be determined during the construction of RSM. Parametric study and optimization for these parameters are performed. Several difficulties during approximation processes are described and numerical examples are demonstrated to verify the efficiency of this method.

In this paper, we proposed a new method to control the length of carbon nano tube in the single CNT probe. A single CNT probe was composed of a tungsten tip made by the electrochemical etching and carbon nano tube which was grown by CVD and prepared through the sonication. The two components were attached with the carbon tape. Since the length of CNT can not be controlled during the manufacturing, the post process is needed to shorten the CNT. In this paper, we proposed the method of electrochemical process. The process was done under the optical microscope and the results were checked by SEM. The diameter of the carbon nano tube used in this paper was about 130nm because the above process had to be done with the optical microscope. Using the method proposed in this paper, we can control the length of the nano tube tip.

A new thermal inkjet printer head on SOI wafer with virtual valve was proposed. It was composed of two rectangular heaters with same size. So we could call it T-jet(Twin jet). T-jet has a lot of merits. It has the advantage of being fabricated with one wafer and is easy to change the size of chamber, nozzle, restrictor and so on. However, above all, It is the best point that T-jet has a virtual valve. And it was manufactured on SOI wafer. The chamber was formed in its upper silicon whose thickness was 40um. The chamber's bottom layer was silicon dioxide of SOI wafer and two heaters were located underneath the chamber's ceiling. And the restirctor was made beside the chamber. Nozzle was molded by process of Ni plating. Ni was 30um thick. Nozzle ejection test was performed by printer head having 56 nozzles in 2 columns with 600NPI(nozzle per inch) and black ink. It measured a drop velocity of 12m/s, a drop volume of 30pl, and a maximum firing frequency of 12KHz for single nozzle ejection. Throwing out the ink drop in whole nozzles at the same time, it was observed that the uniformity of the drop velocity and volume was less than 4%.

Laser heat treatment is an effective technique used to improve the tribological properties and also to increase the service life of automobile components such as camshafts, crankshatfs, lorry brake drums and gears. High power $CO_{2}$ lasers and Nd:YAG lasers are employed for localized hardening of materials and hence are of potential application in the automobile industries. The heat is conducted rapidly into the bulk of the specimen causing self-quenching to occur and the formation of matrensitic structure. In this investigation, the microstructrual features occurring in Nd:YAG laser hardening SM45C and $STD_{11}$ steel are discussed with the use of optical microscopic and scanning electron microscopic analysis. Moreover, This paper describes the optimisation of the processing parameters for maximum hardened depth of SM45C and $STD_{11}$ steel specimens of 10mm thickness by using CW Nd:YAG laser.

The study of Sealless Cylinder is presented. The cylinder has a piston with air bearing. The piston has a partly cylindrical and partly conical shape. The description of system geometry is follows by the flow rate equations. Then pressure distribution and Bearing force equations are derived. Several non dimensional parameters are suggested. The relationship among bearing force, leakage flow and geometry of the bearing is investigated by simulation. And determination method for optimal design of sealless cylinder is given. A prototype of sealless cylinder which had rod bearing with four pockets, five pockets, and six pockets was built respectively.

By development of dissolvable technique and equipment for warhead, empty cartridge and ammunition in small caliver, pollution of environment and waste of resources problems brought by existing incinerative abrogation can be fundamentally prevented. In addition, Automatic high-speed mechanically dissolving technique using indexing equipment developed in this study makes possible curtailment of manpower and recycling treatment of recources.

Nano-imprint technology has been vigorously studied by many researchers for it is one of the most promising technologies for manufacturing the pattern with its critical dimension below 100nm. In the nano-imprint technology, nano patterns are transferred on a polymer film and the transferred patterns are used as an etch mask to define the designed patterns on a substrate or a metal layer. To this end, it is important to keep the residual thickness of the imprinted polymer film uniform. In this study, a novel measurement technique to measure the residual thickness of films is proposed based on nanoindentation theory. This technique has advantages of saving time and measuring the residual thickness of highly-localized portions in comparison with other techniques, but has limitation of requiring calibration process.

Nanometer-sized structures are being applied to many fields including micro/nano electronics, optoelectronics, quantum computing, biosensors, etc. Micro contact printing is one of the most promising methods for manufacturing the nanometer-sized structures. The crucial element for the micro contact printing is the nano-resolution printing technique using polymeric stamps. In this study, a multi-scale analysis scheme for simulating the micro contact printing process is proposed and some useful analysis results are presented. Using the slip-link model [1], the dependency of viscoelasticity on molecular weight of polymer stamp is predicted. Deformation behaviors of polymeric stamps are analyzed using finite element method based upon the predicted viscoelastic properties.

We have applied spectroscopic ellipsometry to investigate $high-{\kappa}$ dielectric thin films and correlate their optical properties with fabrication processes, in particular, with high temperature annealing. The use of high-k dielectrics such as $HfO_{2}$, $Ta_{2}O_{5}$, $TiO_{2}$, and $ZrO_{2}$ as the replacement for $SiO_{2}$ as the gate dielectric in CMOS devices has received much attention recently due to its high dielectric constant. From the characteristics found in the pseudo-dielectric functions or the Tauc-Lorentz dispersions, the optical properties such as optical band gap, polycrystallization, and optical density will be discussed.

The precise structural analysis of $SnO_{2}$ thin film, which was prepared by PECVD and thickness 2400 ${\AA}$, was tried to do the structural refinement using X -ray diffraction data. The observed diffraction patterns of $SnO_{2}$ thin film had the strongly preferred orientation effect. WIMV method was used to correct the preferred orientation effect. The final weighted R-factor, $R_{WD}$ was 7.92 %. The lattice parameters, a = b == 4.7366(1) ${\AA}$ and c = 3.1937(1) ${\AA}$, were almost in accordance with ones of $SnO_{2}$ powder.