In this study, the intrinsic static/dynamic fracture toughness of Al 7175=T74 is evaluated from the apparent static/ dynamic toughness of notched specimen, The critical average stress fracture model is suggested to establish the relationship to predict the intrinsic fracture toughness from the apparent fracture toughness of a notched specimen. The critical average stress fracture model is established using the relation between the notch root radius and the effective distance calculated by finite element analysis. Also, effective distance is applied to estimate the failure criterion for the combustion pipe with notch. It is conclude that the true fracture toughness can be estimated from test results of apparent fracture toughness measured by using a notched specimen. Also, the effective In this study, the intrinsic static/dynamic fracture toughness of Al 7175=T74 is evaluated from the apparent static/ dynamic toughness of notched specimen, The critical average stress fracture model is suggested to establish the relationship to predict the intrinsic fracture toughness from the apparent fracture toughness of a notched specimen. The critical average stress fracture model is established using the relation between the notch root radius and the effective distance calculated by finite element analysis. Also, effective distance is applied to estimate the failure criterion for the combustion pipe with notch. It is conclude that the true fracture toughness can be estimated from test results of apparent fracture toughness measured by using a notched specimen. Also, the effective distance can be used to evaluate the failure criterion of structure with notch.

Large-sized pins are usually used to lift and handle large low speed diesel engine crankshaft. There has then been a need to reduce and optimize the weight of the traditionally used pins. Making an hole by cutting the inside of the pin out was investigated in view of static and fracture strength. To compensate the stress increase caused by the introduction of the inner hole, the groove in the circumferential direction pre-existing on the pin is to be removed. Finite element analysis was carried out for both the original model and weight reduced model. Stress intensity factors for semi-elliptical defects assumed on the pin for the original model and weight reduced model was calculated using the ASME method and compared with the fracture toughness test result of the pin material. The diameter of the cutting hole for the revised model was determined based on the analysis results.

The objective of this study is to investigate fracture toughness and fatigue crack propagation behavior in the Reduced Activation Ferritic Steel (RAFs) JLF-I. The fracture toughness tests were performed with various size(plane size and thickness) and various side groove of specimens. The fatigue crack propagation behavior of the JLF-I steel was investigated by the constant-amplitude loading test for the stress ratios R=O.I, 0.3 and 0.5 respectively. The effects of stress ratios and specimen size on the fatigue crack growth behaviors for JLF-I steel were discussed within the Paris law. The test results showed the standard CT 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 fatigue crack propagation rate of a half size specimen was similar to that of a full size specimen at the stress ratios of 0.1, 0.3 and 0.5 respectively. The fatigue crack propagation behavior of this material were evaluated by using a half size specimen.

A number of creep data were collected and filed for type 316LN stainless steels through literature survey and experimental data produced in KAERI. Using these data, polynomial equations for predicting creep life were obtained for Larson Miller (L-M), Qrr-Sherby-Dorn (O-S-D) and Manson-Haferd (M-H) parametric methods. In order to find out the suitability for them, the relative standard error (RSE) and standard error of estimate (SEE) values were obtained by statistical process of creep data. The O-S-D parameter showed better fitting to creep-rupture data than the L-M or the M-H parameters, and the three parametric methods did not generate the large difference in the SEE and the RSE values.

In this paper, the compensation of interfacial scatter due to adhesive layer and adherend thickness ratio variation was applied to improve measuring precision by calculating ultrasonic attenuation coefficient in the Al/Epoxy dissimilar bonded components. The optimum condition of theoretical value and experimental measuring accuracy by the ultrasonic method in the Al/Epoxy dissimilar bonded components have been investigated. From the experimental results, we proposed a measurement method of the interfacial crack lengths by the ultrasonic attenuation coefficient and discussed it.

The speed competition of optical disk drive has been accelerated with the fast advancement of its storage density and data transmission technology. The continuous increase of the spinning speed of CD meets the unexpected and catastrophic failure of disk during the operation. The effect of its thickness and outer radius of disk were investigated to reduce stresses and J-integral around the crack tip. The effect of its thickness was considered ahead of the crack tip. In the effect of outer radius of disk, linear elastic fracture mechanics was used to obtain the critical crack length, which indicates the onset length for unstable crack growth. This approach is so significant as to detect the growing crack by disk drive before the catastrophic failure, which will provide the standard size of its safety for high-speed disk drive.

A method to measure the interfacial toughness of film/substrate by nanoindenter is proposed. As the thickness of the film decreases, the measurement of the interfacial toughness requires the more sophisticated equipment such as nanoindenter. In this study, the nanoindenter is applied to the substrate near the interface of film/substrate in the direction perpendicular to the normal of the interface, causing the cohesive fracture of the substrate, followed by the interfacial cracking. The specimen of Cu($0.56 {\mu}m$)/Si(530 ${\mu}$) are made by sputtering the copper onto the silicon wafer. By scratching the copper surface, we can make the easy interfacial cracking during the nanoindentation. It is found that the averaged values of the interfacial toughness of the Cu/Si is $0.664{\pm}0.3\;J/m^2$ . The phase angle of the specimen in this study is ${\psi}{\simeq}-36.8^{\circ}$, computed by the method of Suo and Hutchinson.[1]

Metallic glasses are amorphous meta-stable solids and are now being processed in bulk form suitable for structural applications including impact. Bulk metallic glasses have many unique mechanical properties such as high yield strength and fracture toughness, good corrosion and wear resistance that distinguish them from crystalline metals and alloys. However, only a few studies could be found mentioning the dynamic response and damage of metallic glasses under impact or shock loading. In this study, we employed a small explosive detonator for the dynamic indentation on a Zr-based bulk amorphous metal in order to evaluate the damage behavior of bulk amorphous metal under impact loading. These results were compared with those of spherical indentation under quasi-static and impact loading. The interface bonded specimens were adopted to observe the appearances of subsurface damage induced during indentation under different loading conditions.

This paper suggested the method determing the stress intensity factor (SIF) for functionally gradient materials (FGMs) by photo elastic experimental method. The SIF for the center crack in a finite rectangulat plate with a linear variation of shear modulus with constant density and Poisson's ratio along the direction of the crack under mode I static loading is obtained. The exponential and linear variation of stress fields are used for obtaining the SIF. The greater crack length, the increaser the difference of the SIF between right and left side crack tip.

In this paper probabilistic fracture mechanics(PFM) approach is employed to evaluate the integrity of CANDU Zr-2.5Nb pressure tubes. Modified failure assessment diagram(Jr-FAD), plastic collapse, and critical crack length(CCL) approach are used for evaluating failure probability of the tubes. Jr-FAD was extended from the Kr-FAD because fracture of pressure tubes occurs in brittle manner due to hydrogen embrittlement of material by deuterium fluence. For developing the probabilistic integrity evaluation module, AECL procedures and fracture toughness parameters of EPRI were used.

Kim et al. described and compared other methods of measuring stress triaxiality using the displacements near the side necking, proved the validities of these models and explored the effect of location where the displacements are measured using three-dimensional finite element analysis for a standard CT specimen with 20% side-grooves. In this paper, the applicability of these models to various specimen and materials are examined in detail. To consider the effects of side groove, thickness of specimen, crack length, specimen geometry and strain hardening exponents, three-dimensional finite element simulation has been performed for various specimen geometries. For a case without a side groove, in the whole the difference between the stress triaxilaity analytically evaluated and directly determined is similar. For a case with a 20% side groove the stress triaxiality is measured at the area where ${\theta}$ is smaller than $60^{\circ}$, which excludes a side grooved area.

Estimation of the CANDU pressure tube deflection is important since the deflection may cause significant structural failure due to hydrogen diffusion and blister. However, there is no appropriate engineering model to estimate it exactly. The purpose of this paper is to propose a new analysis method and program to resolve this issue. For development of proper analysis method, a series of finite element analyses has been carried under elastic-creep condition. In addition, for effective estimation of the creep deflection, an analysis program named PC-DAS was developed based on the proposed method. Comparison of simple case study results with corresponding reference ones showed good agreement. Therefore, the proposed method and program can be utilized as one of valuable toolkit for integrity assessment of CANDU pressure tube.

In a primary reactor cooling system(RCS), a dissimilar weld zone exists between cast stainless steel(CF8M) in a pipe and low-alloy steel(SA508 cl.3) in a nozzle. Thermal aging is observed in CF8M as the RCS is exposed for a long period of time to a reactor operating temperature between 290 and $330^{\circ}C$, while no effect is observed in SA508 cl.3. The specimens are prepared by an artificially accelerated aging technique maintained for 300, 1800 and 3600 hrs at $430^{\circ}C$, respectively. The specimens for elastic-plastic fracture toughness tests are prepared one type, which notch is created in the heat affected zone(HAZ) of CF8M. And, the specimens for fatigue crack growth tests are prepared in three classes, which notches are created at the center of deposited zone, the HAZ of CF8M, and the HAZ of SA508 cl.3. From the experiments, the J-integral values with the increase of aging time decrease, and the differences of the fatigue crack growth behaviors are relatively small in the three classes specimens.

Environmental fatigue crack propagation of CF8M and CF8A steels used in the domestic PWR were investigated on the simulated PWR condition(Temperature: $316^{\circ}C$, Pressure: 15MPa). The test equipment for environmental fatigue(high temperature-high pressure loop, autoclave, load frame, measurement system) were designed. As-received and 60-year aged specimens were used in the test. To compare with environmental fatigue test, another test was performed in the air condition. The fracture surface of specimens were difficult to verify the fracture modes such as striation, intergranular crack and cleavage and so on. As the ferrite content of CF8M is increased, more particles covered fracture surface were peeled.

When the tube contacted to support, anti-vibration bar of the steam generator in nuclear power plant, the contact area is worn out by their relative displacement and contact force. Connors and Au-Yang found the relation between tube worn displacement and volume, or normal work rate at given gap size. The present analysis is obtained the relation between tube worn displacement and normal work rate at various gap size modifying Au-Yang's result. The results are compared with Connors and Yettisir and Pettigrew's results. The comparison shows that Yettisir and Pettigrew result is fairly good agreement with Connors and present results with gap clearance, 0.015in.

In this study, we proposed new method to mitigate tensile welding residual stress for preventing PWSCC in CRDM nozzle. Residual stress analysis using finite element method is performed to confirm benefit of the new method. In case of applying existing method, tensile axial residual stress decrease by about 28% and tensile hoop residual stress decrease by about 33%. In case of applying the new method, tensile axial residual stress decrease by about 32% and tensile hoop residual stress decrease by about 43%. Therefore, we conclude the new proposed method is more effective to prevent PWSCC than existing method.

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 embrittlement at the reactor operating temperature. The objective of this study is to summarize the method of estimating ferrite content, Charpy impact energy and J-R curve and to evaluate the thermal embrittlement of the cast austenitic stainless steel piping used in the domestic nuclear power plants. The result of evaluation, two domestic nuclear power plants used CF-8M and CF-8A material has adequate fracture toughness after saturation.

Suspension part should have enough endurance during its lifetime to protect passenger. Therefore, the coil spring is one of the major suspension part of an automobile. Corrosion fatigue strength of the coil spring depends on many factors including mechanical and environmental properties. In this paper, residual stresses by shot peening was analyzed using finite element analysis and evaluated its effect on corrosion fatigue strength.

According to our fatigue tests carried out at 20 Hz, R=0.1 on transversely butt~welded joints, fatigue strengths of as-welded specimens, that is, specimens having residual stress are higher than those of annealed specimens in short life range, but vice verse in long life range. This behavior seems to be concerned mainly with residual stress relaxation by applied loading. After analyzing the welding process, we conducted finite element analysis to quantify the degree of residual stress relaxation. By taking into account residual stress relaxation, modified Goodman diagram, and nominal stress, we evaluated the fatigue life of the welded joint from the S-N curve for the parent material. The estimated results are in a good agreement with the experimental results.

Endbeam is an important structural member of freight bogie for the support of service loading. In general, more than 25 years' durability is necessary. However, endbeam occur fatigue fracture in dynamic stress concentration location because comparatively strength and stiffness are low. Therefore, structure analysis is performed to evaluate structural problem of endbeam and local strain range as durability analysis. The number of cycles is extracted concerning the bogie in operation by measurement dynamic stress time history on critical part which is crack initiation in actual fact. At this time rainflow cycle counting is used to consider change of stress for operating condition. Based on the fatigue life curves and the stress analysis, the fatigue life of the endbeam is predicted and compared with the experimentally determined fatigue life, resulting in a fairly good correlation.

Fatigue life of welded joints are governed by the propagation of multiple collinear surface cracks distributed randomly along weld bead. These cracks propagate in mutual interaction and coalescence of them. To estimate the fatigue life, the influences of above two mechanisms on the fatigue life should be taken into account. These two mechanisms appear through the stress intensity factors disturbed mutually. However, it is difficult to calculate the stress intensity factors of multiple surface cracks located in vicinity of weld toe. The stress intensity factors are calculated normally by using the Mk-factors, but such Mk-factors are very rare in literature. In this study, the Mk-factors were obtained from a parametric study on crack length and depth, in which a finite element method is used. A fatigue test for a cruciform welded joint was conducted. The fatigue life of the tested specimen was estimated through present method with the informations obtained from the test, e.g. the number, size and locations of the cracks. The estimated and measured fatigue life showed a good agreement.

Chemically modified multiwalled carbon nanotubes with acids were incorporated into a epoxy matrix by in situ polymerization process, to improve the transfer of mechanical load through chemical bonds, which were demonstrated by infrared spectroscopy. And the mechanical properties of epoxy/carbon nanotube composites were measured to investigate the role of carbon nanotubes. The epoxy/carbon nanotube composites shows higher tensile strength and wear resistance than existing epoxy, with 1 or 2 wt. % addition of functionalized carbon nanotubes. The tensile strength with 7 wt. % carbon nanotibes is increased by a 28% and the wear resistance in exceptionally increased by an outstanding 100 times.

In this work, effects of hyper-elastic rubber material properties on the indentation load-deflection curve and subindenter deformation are first examined via [mite element (FE) analyses. An optimal data acquisition spot is selected, which features maximum strain energy density and negligible frictional effect. We then contrive two normalized functions. which map an indentation load vs. deflection curve into a strain energy density vs. first invariant curve. From the strain energy density vs. first invariant curve, we can extract the rubber material properties. This new spherical indentation approach produces the rubber material properties in a manner more effective than the common uniaxial tensile/compression tests. The indentation approach successfully measures the rubber material properties and the corresponding nominal stress.strain curve with an average error less than 3%.

Recently there has been a great world-wide interest in developing and characterizing new nano-structured materials. These newly developed materials are often prepared in limited quantities and shapes unsuitable for the extensive mechanical testing. The development of depth sensing indentation methods have introduced the advantage of load and depth measurement during the indentation cycle. In the present work, ZnO thin films are prepared on Si(111), Si(100) substrates at different temperatures by pulsed laser deposition(PLD) method. Because the potential energy in c-axis is low, the films always show c-axis orientation at the optimized conditions in spite of the different substrates. Thin films are investigated by X-ray diffractometer and Nano indentation equipment. From these measurements it is possible to get elastic modulus and hardness of ZnO thin films on Si substrates.

Porous anodic alumina has been used widely for corrosion protection of aluminum surfaces or as dielectric material in micro-electronics applications. It exhibits a homogeneous morphology of parallel pores which can easily be controlled between 10 and 400nm. It has been applied as a template for fabrication of the nanometerscale composite. In this study, mechanical properties of the AAO structures are measured by the nano indentation method. Nano indentation technique is one of the most effective method to measure the mechanical properties of nano-structures. Basically, hardness and elastic modulus can be obtained by the nano-indentation. Using the nano-indentation method, we investigated the mechanical properties of the AAO structure with different size of nano-holes. In results, we find the hole effect that changes the mechanical properties as size of nano hole.

The base material of optical disc is the amorphous polymer, polycarbonate. So crack growth on the rotating optical disc could happen by the static centrifugal force. This phenomenon is called static fatigue. Today's CD-R disc drive is being operated over 10,000 RPM. This is increasing the possibility of the disc fracture when operating. In this reason, new method to measure the static fatigue threshold quality and the way to calculate the threshold J-integral value and the safe crack length of the optical disc are studied. Finally the environmental effect to optical disc is also studied in this paper.

To predict the direction of the fatigue crack initiated from a hole under various types of biaxial fatigue loads with different phase difference and biaxiality, fatigue parameters were investigated. Axial and torsional biaxial fatigue loads were selected with the respective combination of five different phase differences of 0, 45, 90, 145 and 180 degrees and five biaxialities of 0, $1/{\sqrt{3}}$, 1, ${\sqrt{3}}$, ${\infty}$. Directions of the fatigue crack initiation around the hole were found to approach to the circumferential direction of the specimen with increment of the phase difference for fatigue tests with phase differences less than $90^{\circ}$. Whereas directions for tests with phase differences greater than $90^{\circ}$ went away from the circumferential direction and those were symmetric to the directions for tests with phase difference less than $90^{\circ}$. With increase of biaxilities, the fatigue crack initiated more apart from the circumferential direction of the specimen. These crack initiation direction were predicted using maximum tangential stress range and maximum shear stress range obtained at far-field and around the hole. Comparing these two stress parameters, The crack initiation direction can be successfully explained by using the direction of the maximum tangential stress range obtained around the hole and at far-field.

The low velocity impact characteristics of composite laminate curved beams are investigated to increase damage tolerance and reduce the deflection. Drop weight impact tests of the composite curved beam were performed with respect to pre-load, then the damage after impact was measured by macrography. Also, finite element analyses were performed using ABAQUS to investigate the stress state of composite curved beam with respect to pre-load and impact. From the investigation, it was found that pre-load of the composite curved beams had much influence on impact damage of the curved beam, which showed good agreement with the experiment results.

Fretting is a kind of surface degradation mechanism observed in mechanical components and structures. The fretting damage decrease in 50-70% of the plain fatigue strength. This may be observed in aircraft, automobile and nuclear power plant used in special environment and various loading conditions. In the present study, the characteristics of the fretting fatigue are investigated using the two aluminum alloy(Al2024-T3511 and Al7050-T7451). Through the experiment, it is found that the fretting fatigue strength of the Al7050-T7451 alloy decreased about 50% from the plain fatigue strength, while the fretting fatigue strength of the Al2024-T3511 alloy decreased about 45%. The tire track was widely observed in fracture surface area of oblique crack which was induced by contact pressure. These results can be the basic data to the structural integrity evaluation of aluminum alloy subjected to fretting damage.

The growth of the critical size bubble by diffusion process in viscoelastic medium was treated by an integral method for the concentration boundary layer adjacent to the bubble wall. In this study, we obtained a set of the first order time dependent equations to obtain bubble radius and gas pressure inside the bubble simultaneously. The calculated final cell sizes depending on the initial saturation pressure are in close agreement with the observed ones. The governing equations developed in this study may be used in polymer processing of microcellular foams.

This paper investigates asimply supported orthotropic rectangular laminate with viscous interfaces subjected to bending. Additional mathematical difficulty is involved due to the presence of viscous interfaces because the behavior of the laminate depends on time. A step-by-step state-space approach is suggested, which is directly based on the threedimensional theory of elasticity. In particular, Taylor's expansion theorem is employed to model the variations of field variables with time. The proposed method is suitable for analyzing laminated plate of arbitrary thickness. Numerical calculations are performed and it is shown that the viscous interfaces have a significant fluence on the response.

we can say that the increasing range of the value of GMT Sheet's tensile strength in the direction of fiber orientation is getting wider as the fiber content increases. It shows that the value of GMT Sheet's tensile strength in the direction of fiber orientation 90 is similar with the value of pp's intensity when fiber orientation function is J= 0.7, regardless of the fiber content. Tensile strength of GMT Sheet is affected by the fiber orientation distribution more than by the fiber content.

The anisotropy in coefficient of thermal expansion (CTE) between the in-plane and out-of-plane of 3-dimensional thick composite structures induces residual stresses and the large void content due to insufficient compaction of fabric composites, which results in low interlaminar strengths. In order to reduce the through thickness CTE and the void content, in this work, carbon fabric phenolic laminates were compacted by pressure generated by autoclave and a compressive jig, from which the through-thickness CTEs and the void contents were measured. From the measurement, it was found that the through-thickness CTE and the void content had different characteristics from ordinary composites due to gas produced during the cure reaction of phenolic resin.

In these days, the orthodontic surgery including lingual orthodontics has attracted a person' attention due to its functional and esthetic appreciation. The delivery of the optimal orthodontic treatment is greatly influenced by clinician' ability to predict and control the tooth movement by applying force system to dentition. The skeletal anchorage system with the miniscrew has been used recently in the lingual orthodontics to assist the anchorage control. Precise understanding of the force system produced from the various orthodontic appliances is necessary. However, the qualitative and quantitative effect of the miniscrew has not been identified well. In this paper, three dimensional finite element analysis is introduced on the lingual orthodontics to investigate the effect of anterior retraction force on the miniscrew and transpalatal arch wire. The purpose of this study is to determine the location of the miniscrew and the point of force application of the anchorage system in the lingual orthodontics. The analysis results indicate the efficient position of the miniscrew and the transpalatal arch wire in the lingual orthodontics.

SiC materials have been extensively studied for high temperature components in advanced energy system and advanced gas turbine. However, the brittle characteristics of SiC such as low fracture toughness and low strain-to fracture still impose a severe limitation on practical applications of SiC materials. For these reasons, $SiC_f/SiC$ composites can be considered as a promising for various structural materials, because of their good fracture toughness compared with monolithic SiC ceramics. But, high temperature and pressure lead to the degradation of the reinforcing fiber during the hot pressing. Therefore, reduction of sintering temperature and pressure is key requirements for the fabrication of $SiC_f/SiC$ composites by hot pressing method. In the present work, Monolithic LPS-SiC was fabricated by hot pressing method in Ar atmosphere at 1760 $^{\circ}C$, 1780 $^{\circ}C$, 1800 $^{\circ}C$ and 1820 $^{\circ}C$ under 20 MPa using $Al_2O_3-Y_2O_3$ system as sintering additives in order to low sintering temperature. The starting powder was high purity ${\beta}-SiC$ nano-powder with an average particle size of 30 nm. Monolithic LPS-SiC was evaluated in terms of sintering density, micro-structure, flexural strength, elastic modulus and so on. Sintered density, flexural strength and elastic modulus of fabricated LPS-SiC increased with increasing the sintering temperature. In the micro-structure of this specimen, it was found that grain of sintered body was grown from 30 nm to 200 nm.

: A DNA analysis system based on fluorescence analysis has to have a DNA amplifying thermal cycle system. DNA amplification is executed by the temperature control. Accuracy of fluorescence analysis is influenced by the temperature control technology. For that reason, the temperature control is core technology in developing the DNA analysis system. Therefore, the objective of this paper is to develop the hardware to apply thermoelectric module to the DNA amplifying thermal cycle system. In order to verify the developed hardware for controlling the temperature of thermoelectric module, a DNA amplifying thermal cycle test was performed. From the test, the developed hardware controlled the temperature of thermoelectric module successfully. Therefore, it is expected that the developed hardware can be applied to the DNA amplifying thermal cycle system.

Linepipe steels with a low carbon acicular ferrite microstructure have been recently developed to accommodate the current transportation condition of the gas and oil industry, and they are finally applied to West- East pipeline project in China. By adopting acicular microstructure, both better formability and better toughness could be obtained due to low yield ratio and fine grained microstructure. Mechanical properties of pipe are not greatly different from those of base plates or hot coils with a microstructure of acicular ferrite. Merits of introducing higher strength steels are well known, i.e., reducing the gauge of pipe and the material cost, increasing the welding speed and decreasing construction cost because of reducing the construction period. Threfore, gas and oil industry has required higher strength steel than APIX70 grade steel. Under this background, API-X80 steel has been developed and shall be applied to the several projects. In this paper, developing stage of API-X80 steel is also presented and discussed.

The failure assessment for corroded pipeline has been considered with the full scale burst test and the finite element analysis. The burst tests were conducted on 762 mm diameter, 17.5 mm wall thickness and API 5L X65 pipe that contained specially manufactured rectangular corrosion defect. The failure pressure for corroded pipeline was measured by burst testing and classified with respect to corrosion sizes and corroded regions - the body, the girth weld and the seam weld of pipe. Finite element analysis was carried out to derive failure criteria of corrosion defect on the pipe. A series of finite element analyses were performed to obtain a limit load solution for corrosion defects on the basis of burst test. As a result, the criteria for failure assessment of corrosion defect within the body, the girth weld and the seam weld of API 5L X65 gas pipeline were proposed.

The length of city gas pipeline is increasing with expansion of natural gas transmission rapidly. A lot of the expense was paid for repair and maintenance with increasing of pipeline length and the cost of repair and maintenance by the corrosion was the highest. It is necessary to evaluate integrity in case of thickness reduction by corrosion. There are a lot of assessment criteria for corrosion defect in foreign countries but they are not suitable for application in the country directly. In this work, we performed the burst test and the finite element analysis for city gas pipeline, KS D3507 and KS D3631 for city gas transmission, and developed the assessment method of corrosion defect, which is suitable for domestic condition.

A study on the probabilistic methodology for the estimation of the remaining life of Pressurized pipelines containing active corrosion defects is presented. This reliability assessment is earned out using extream value distribution of the corroded defects instead of already published failure perssure moded like NG18 or ASME B31G. The failure probability of pipelines depends on the number of corroded defects. and it could be calculated directly as the area exceeded a defined L V(Limited Value of corrosion depth). The remaining life of pressurized pipelines can also be estimated by the PDF of extream value distribution as calculating the exceeded area with a defined failure probability.

Tensile properties of hard coating material, TiN, were evaluated using micro-tensile testing system. TiN has been known as a hard coating material commonly used today. Micro-tensile testing system consisted of a micro tensile loading system and a micro-ESPI(Electronic Speckle Pattern Interferometry) system. Micro-tensile loading system had a maximum load capacity of 500mN and a resolution of 4.5 nm in stroke. TiN thin film $1{\mu}m$ thick was deposited on the Si wafer pre-deposited of $Si_3N_4$ film substrate by the closed field unbalanced magnetron sputtering (CFUBMS) process. Three kinds of micro-tensile specimen with the respective width of $50{\mu}m$, $100{\mu}m$ and $500{\mu}m$ were fabricated by MEMS process. The mechanical properties including tensile strength and elastic modulus were determined using the micro-tensile testing system and compared by those obtained by nano-indentation

Nano-indentation is used for measuring mechanical properties of thin films such as elastic modulus and hardness. For ductile materials, pile-up around the indenter causes the calculation of inaccurate projected contact area. This phenomenon was found by measurement of indentation shape using an atomic force microscope. In present study finite element analysis of nano-indentation was performed to compensate the effects of pile-up on the contact area. The result of finite element analysis was compared with that of nano-indentation for a ductile material. The analysis has demonstrated that the true contact area is greater than that calculated by nano-indentation. It is verified that the consideration of the effects of pile-up in nanoindentation for ductile materials using the finite element method is reasonable.

Nanometer-sized structures are being applied to many devices including micro/nano electronics, optoelectronics, quantum devices, MEMS/NEMS, biosensors, etc. Especially, the thin film with submicron thickness is a basic structure for fabricating these devices, but its mechanical behaviors are not well understood. The mechanical properties of the thin film are different from those of the bulk structure and are difficult to measure because of its handling inconvenience. Several techniques have been applied to mechanical characterization of the thin film, such as nanoindentation test, micro/nano tensile test, strip bending test, etc. In this study, we focus on the strip bending test because of its high accuracy and moderate specimen preparation efforts, and measure Au thin film, which is a very popular material in micro/nano electronic devices. Au film is deposited on Si substrate by evaporation process, of which thickness is 100nm. Using the strip bending test, we obtain elastic modulus, yield and ultimate tensile strength, and residual stress of Au thin film.

Mechanical property evaluation of micrometer-sized structures is necessary to help design reliable microelectromechanical systems(MEMS) devices. Most material properties are known to exhibit dependence on specimen size and such properties of microscale structures are not well characterized. This paper describes techniques developed for tensile testing of materials used in MEMS. Epi-polycrystalline silicon is currently the most widely used material, and its tensile strength has been measured as 1.52GPa. We have developed an uniaxial testing machine for testing microscale specimen using electro-magnetic actuator. The field magnet and the moving coil taken from an audio-speaker were utilized as the components of the actuator. Structure of specimen was designed and manufactured for easy handling and alignment. In addition to the static tensile tests, new techniques and procedures for measuring strength are described.

Wall thinning of carbon steel pipe components due to Flow-Accelerated Corrosion (FAC) is one of the most serious threats to the integrity of steam cycle piping systems in Nuclear Power Plants (NPP). Since the mid-1990s, secondary side piping systems in Korean NPPs have experienced wall thinning, leakages and ruptures caused by FAC. Korea Electric power Research Institute (KEPRI) and Korea Hydro & Nuclear Power Co., LTD. (KHNP) have conducted a study to develop the methodology for systematic pipe management and established the Korean Thinned Pipe Management Program (TPMP). To effectively maintain the integrity of piping system, FAC engineer should understand the criterions of the structural integrity evaluation and the safety margin assessment for the thinned pipe component. This paper describes the technical items of TPMP, and shows the example of the integrity evaluation and safety margin assessment for three thinned pipe component of a NPP.

This study is conducted to clarify the effect of internal pressure on the deformation and collapse behaviors of wall thinned elbow under in-plane bending moment. Thus the nonlinear three-dmensional finite element analyses were performed to obtain the moment-rotation curve of elbow contatining various wall thinning defects located at intrados and extrados under in-plane bending (closing and opening modes) with internal pressure of $0{\sim}15MPa.$ From the results of analysis, the effect of internal of collapse moment of elbow on the global deformation behavior of wall thinned elbow was discussed, and the dependence of collapse moment of elbow on the magnitude of internal pressure was investigated under different loading mode, defect location, and defect shape.

The most important mode of subcritical crack growth is primary water stress corrosion crack, which was the reported mechanism from the root cause analysis of the crack in the bimetallic welds. Stress corrosion crack growth evaluations was carried out for several flaw shapes of both axial and circumferential flaws, using the steady-state stresses including residual stresses. This evaluation considered the possibility of additional flaws in the primary loops of nuclear power plant, even though no such flaws have been identified by Ultrasonic Test. Consequently, Results show that the predicted flaw sizes will determine acceptability for continued service and maintenance.

Many analysis methods, including finite element method, have been suggested and used for assessing the integrity of cracked structures. In the paper, in order to analyze arbitrarily shaped three dimensional cracks in an infinite body, the finite element alternating method is extended. The cracks are modeled as a distribution of displacement discontinuities by the displacement discontinuity method and the symmetric Galerkin boundary element method. Applied the proposed method to several example problems for planner cracks in finite bodies, the accuracy and efficiency of the method were demonstrated.

The objective this study is to estimate the feasibility of acoustic emission method for the internal leak from the valves in nuclear power plants. From the experimental results, it was suggested that the acoustic emission method for monitoring of leak was feasible. When the background levels are higher than the acoustic signals from leak, we can detect the leak analyzing the spectrum of the remainders which take the background noise from the acoustic signals.

In this study, the pressure vessel piping with corrosion used during long term were investigated from the time-frequency analysis method. The damage of piping could be evaluated the attenuation factor by ultrasonic parameters such as center frequency and echo waveform. Based on NDE analysis by time-frequency analysis method, it should also be possible to evaluate from various damages and defects in piping members.

By the purpose to investigate the change of internal stress and mobile dislocation density in creep, the stress relaxation test is carried out in the condition of each strain. Mobile dislocation density increased until it reached minimum creep rate and after that, it decreased and internal stress didn't have the change approximately until it reached minimum creep rate and after that, it decreased. The stress relaxation rate is fast and approached zero after 1.5 seconds after the beginning of the stress relaxation. And the larger the applied stress is, the larger the internal stress is. By the evaluation of mobility of dislocation, the dislocations glide viscously in STS31OJlTB but it is the dislocations glide viscously which N passes by cutting Cr atom rather than typical viscosity movement after calculating mobility of dislocation.

This paper reviewed characteristics of fatigue crack behavior observed by changing various shapes of initial crack and magnitudes of loading in compact tension shear(CTS) specimen subjected to shear loading. In the high-loading condition, fatigue crack under shear loading propagated branching from the pre-crack tip. Meanwhile, the secondary fatigue crack in the low-loading condition which was created in the notch root due to friction on the pre-crack face grew to a main crack. Influenced by the mode II loading condition, fatigue crack propagation retardation appeared in the initial propagation region due to the reduction of crack driving force and friction on crack face. In both cases, however, fatigue cracks grew in tensile mode type. Propagation path of fatigue crack under the shear loading was 70 degree angle from the initial crack regardless of its shape and load magnitude.

Recently, the steel parts used for automobiles and trains are required to be used under higher stress than ever before in need of the weight down. In this study, high strength and superior toughness spring steels as the suspension material, used for automobile and railroad industries were utilized to carry out the following in vestigations; 1) To evaluate the characteristics of fatigue crack propagation, the experiments of fatigue crack growth were respectively carried out at the room temperature(RT), $100^{\circ}C$, $200^{\circ}C$ 2) Peening and unpeening materials at the each temperature were investigated for the effect on fracture toughness by compressive residual stress generated from the shot peeing.

High temperature fracture toughness characteristics of shot peened spring steel(SUP-9), which is used for automobile suspension system and railroad, was investigated in this paper. Fracture tougness test for room temperature, $100^{\circ}C$ , and $200^{\circ}C$ were evaluated by material test system(MTS). The experimental results show that the fracture toughness was improved by peened and unpeened. The fracture toughness for high temperature were also improved by peened and unpeened.

This work investigated fatigue characteristics of aluminum repaired by CFRP composites. Three specimens, cracked aluminum, cracked aluminum patched by CFRP, and plasma-treated aluminum patched by CFRP were used for the fatigue tests. The results showed that the fatigue crack growth behavior of cracked aluminum was improved by repairing the cracked area with composite patch. Specifically, the specimen repaired by composite patch showed about three times more fatigue life than the cracked aluminum. The plasma-treated aluminum repaired by composite patch showed about five times more fatigue life than the cracked aluminum.

In this study, an application of crack morphology variables in the Leak-Before-Break(LBB) evaluation for nuclear piping systems is investigated, including influences on the leakage crack size and crack instability loads. The crack surface roughness and the number of flow turns as a function of the crack opening displacement are applied to LBB evaluations for KSNP pressurizer surge line, for which fatigue and stress corrosion cracking are considered as failure mechanisms. As a result, there would be a significant impact on safety margins to acceptance criteria for the surge line if crack morphology variables are applied additionally to the current regulatory guide without re-analyses for justification of safety factors being applied on the leakage crack size and piping loads for evaluations.

The needs of walking assistive device such as the Ankle Foot Orthosis (A.F.O.) are getting greater than before. However, most of the A.F.O. are generally imported rather than domestic manufacturing. The major reason of high import reliability is the rack of impact properties of domestic commercial products. Therefore, this research is going to focus on the evaluation of impact properties of the A.F.O. which has the high import reliability. Unfortunately, these kinds of researches are not performed sufficiently. This research is going to evaluate impact energy behavior in composite materials such as the glass/epoxy (S-glass, $[0/90]_{2S}$) and the aramid/epoxy (Kevlar-29, woven type, 8 ply) of ankle foot orthosis. The approach methods were as follows. 1) The history of impact load and impact energy due to the various velocities. 2) Relationship between the deflection and damage shape according to the impact velocities. 3) The behavior of absorbed energy and residual strength rate due to the various impact velocities.

This research investigates the feasibility of ultrasonic microscope for nondestructive assessment of thermal degradation in artificially aged commercial Co-base superalloy, FSX414. This alloy has been used for high temperature structure applications such as stationary gas turbine blade and nozzle chamber in fossil plant. Microstructural change was found that the fine carbides became coarser and spheroidized in matrix as aging time increased. The leaky surface acoustic wave velocity gradually decreases by a maximum of 4.7% with increasing aging time up to 4,000hours. However, the longitudinal wave velocity has a little change. Also, it has a good correlation between leaky surface acoustic wave velocity and Vickers hardness. Consequently, LSAW can be used to examine the degree of degradation in thermally aged Co-base superalloy.

The deformation behavior of a $Zr_{55}Al_{10}Ni_5Cu_{30}$ bulk metallic glass under tensile loading at different range of strain rates and temperatures between 680 K and 740 K were investigated. The variation in the deformation behavior of $Zr_{55}Al_{10}Ni_5Cu_{30}$ bulk metallic glass which resulted from the crystallization induced during testing was reported. The$Zr_{55}Al_{10}Ni_5Cu_{30}$ bulk metallic glass has showed either homogeneous or inhomogeneous deformation depending on test condition. It exhibited a maximum elongation of about 560 % at the condition of $407^{\circ}C{\times}\;10^{-4}/s$. The flow behavior exhibited three different types and the flow stress became lower at higher temperatures and lower strain rates. The increase of the time elapsed during heating resulted in the partial crystallization of bulk metallic glass phase and eventually strain hardening and brittle fracture.

Recently the steel parts used for automiles and trains are required to be used under higher stress than ever before in need of the weight down. However, threr are a lot of problems with developing such of fatigue strength and fatigue life are mainly focused on by adopting residual stress. And got the following characteristics from crack growth test carried out stress ratio. Fatigue life shows more improvement in the Un-peening material. And Compressive residual stress of surface on the Shot-peening processed operate resistance force of fatigue. So we cam obtain fallowings. (1) The fatigue crack growth rate on stage II is conspicuous with the size of compressive residual stress and is dependent of Paris equation. (2) Although the maximum compressive residual stress is deeply and widely formed from surface, fatigue life does not improve than when maximum compressive residual stress is formed in surface. (3) The threshold stress intensity factor range is increased with increasing compressive residual stress.

We got the following characteristics from fatigue crack growth test carried out in the environment of room temperature and low temperature at $25^{circ}C$, $-60^{circ}C$, $-80^{circ}C$, and $-100^{circ}C$ in the range of stress ratio of 0.3 by means of opening mode displacement. And there is a difference between shot peened specimen and unpeened specimen. The purpose of this study is to predict the behavior of fatigue crack propagation as one of fracture mechanics on the compressive residual stress. Fatigue crack growth rate of shot peened metal was lower than that of unpeened metal. The compressive residual stress made an impact on tension and compression of the plasticity deformation in fatigue crack plasticity zone. That is. the constrained force about plasticity deformation was strengthened by resultant stress, which resulted from plasticity deformation and compressive residual stress in the process of fatigue crack propagation.

We investigated the effect of weld details on fatigue behavior of a material, JIS SM 490 A, with yielding strength of about 350 MPa and tensile strength of about 520 MPa. Tensile tests, instrumented indentation tests and fatigue tests were carried out on double V-grooved butt weld plates such as reinforcement removed, as-welded and weld toe ground. In addition plates with transverse fillet welded web, load carrying cruciform fillet welded plates, non-load-carrying cruciform fillet welded plates and longitudinal butt welded plates were tested. S-N curves for the above specimens were obtained and analyzed

Generally, component and FR-4 board are connected by solder joint. Because material properties of components and FR-4 board are different, component and FR-4 board show different coefficients of thermal expansion (CTE) and thus strains in component and board are different when they are heated. That is, the differences in CTE of component and FR-4 board cause the dissimilarity in shear strain and solder joint' failure. The first order Taylor series expansion of the limit state function incorporating with Tresca failure criterion is used in order to estimate the failure probability of solder joints under heated condition. Using shear stresses and shear strains appeared on the solder joint, we estimate the failure probability of solder joints with the Tresca failure criterion. The effects of random variables such as CTE, distance of the solder joint from the neutral point(DNP), temperature variation and height of solder on the failure probability of the solder joint are systematically studied by using the failure probability model with first order reliability method(FORM).

It is well known that a specific experimental method such as the Split Hopkinson Pressure Bar (SHPB) technique is the simplest experimental technique to determine the dynamic material properties under the impact compressive loading conditions with strain-rate of the order of $10^3/s{\sim}10^4/s$. This type of experimental procedure has been widely used with proper modification on the test setups to determine the varying dynamic response of materials for the dynamic boundary conditions such as tensile and fracture as well. In this paper, dynamic compressive deformation behaviors of an Ethylene Copolymer materials widely used for the isolation of vibration from varying structures under dynamic loading are estimated using the SHPB technique.

Electrical conductivity of Al-Si-Cu-Mg alloy foams of various density produced in powder metallurgical method has been measured using two probe electrical conductivity measurement method. Compressive mechanical properties such as elastic modulus and plastic plateau stress of foams were evaluated from electrical conductivity using power law relation and scaling laws of foam properties. Uni-axial compression test was also performed. Experimentally measured elastic modulus and plastic plateau stress were compared with the values evaluated from electrical conductivity. The computed values were in good agreement with the experimental result.

Nanoindentation technique has been used to measure the mechanical properties of aluminium alloy foam cell walls. Al-Si-Cu-Mg alloy foams of different compositions and different cell morphologies were produced using powder metallurgical method. Cell morphology of the foam was controlled during production by varying foaming time and temperature. Mechanical properties such as hardness and Young's modulus were calculated using two different methods: a continuous stiffness measurement (CSM) and an unloading stiffness measurement (USM) method. Experimental results showed that hardness and Young's modulus of Al-5%(wt.)Si-4%Cu-4%Mg (544 alloy) precursor and foam walls are higher than those of Al-3%Si-2%Cu-2%Mg (322 alloy) precursor and foam walls. It was noticed that mechanical properties of cell wall are different from those of precursor materials.

A comprehensive numerical study on the three point bending behavior of Aluminum foam-filled stainless steel tube has been performed. Aluminium alloy foams with various densities were produced and their mechanical properites were evaluated. Finite element(FE) analysis of three point bending test was performed to evaluate bending behavior of foam filled cylindrical structures. Results showed that foam filling offered remarkable increase of bending resistance and enhanced the crashworthiness of the structure. It turned out to prevent the inward fold formation at the compression flange, resulted into the multiple propagating folds and increased the load carrying capacity.

In order to secure the safety of power steering system, it is necessary to perform the strength and the fatigue analysis of pulley in this paper. The applied stress distribution of the pulley subject to combined loads condition was obtained using finite element analysis. Based on these results, the fatigue life of the pulley with the variation of the fatigue strength was evaluated using durability analysis simulator. The optimal hole size to improved the safety of the pulley was investigated using parametric study. Moreover, the predicted fatigue life cycle with the simulator was verified by experimental tests.

In this study, the optimal welding condition of the input power was selected experimentally through the ERW simulator, which is equal to welding status of ERW part in pipe. This condition is the input power 250kW in the heat treatment of the $900^{\circ}C$ normalizing derived from the nondestructive technique and impact energy. In order to evaluate the variation of the fatigue life in the pipe, fatigue surface crack growth test of base and optimal welded metal were performed statistically. As stress intensity factor range (${\Delta}K_s$) increases, the fatigue crack propagation rate (da/aN) of the base metal is faster than that of the welded joint. The variation of the fatigue life in the ERW pipe was estimated statistically using Monte-Carlo simulation with the standard deviation of material constants (C and m) of the paris law in the specimen.

In this study, we investigated compressive characteristics of seawater-absorbed carbon-epoxy composite under hydrostatic pressure environment. The hydrostatic pressures applied were 0.1 MPa, 100 MPa, 200 MPa, and 270 MPa. The results showed that the compressive elastic modulus increased about 10 % as the hydrostatic pressure increased from 0.1 MPa to 200 MPa. The modulus increased 2.3 % more as the pressure increased to 270 MPa. Fracture strength and fracture strain increased with pressure in a linear fashion. Fracture strength increased 28 % and fracture strain increased 8.5 % as the hydrostatic pressure increased from 0.1 MPa to 270 MPa.

Under longitudinal loading continuous fiber reinforced metal matrix composite(MMC) have interpreted an outstanding performance. However, the applicability of continuous fiber reinforced MMCs is somewhat limited due to their relatively poor transverse properties. Therefore, the transverse properties of MMCs are significantly influenced by the properties of the fiber/matrix interface. In this study, elastic-plastic behavior of transversely loaded unidirectional fiber reinforced metal matrix composites investigated by using elastic-plastic finite element analysis. Different fiber placement(square and hexagon) and fiber volume fractions were studied numerically. The interface was treated as three thin layer (with different properties) with a finite thickness between the fiber and the matrix. The analyses were based on a two-dimensional generalized plane strain model of a cross-section of an unidirectional composite by the ANSYS finite element analysis code.

Thermal barrier systems are susceptible to instability of the thermally grown oxide(TGO) at the interface between the bond coat(BC) and the thermal barrier coating(TBC). The instabilities have been linked to thermal cycling and initial geometrical imperfections, as well as to misfit strains due to oxide growth and expansion misfit. In this work, deformation of TGO near a surface groove due to thermal cycling has been observed at high temperatures, $1100^{circ}C$, $1150^{circ}C$, $1200^{circ}C$. The effect of peak temperature and the thickness of substrate are presented.

Lee and Kang measured side-necking deformation near a crack-tip for CT specimen using Stereoscopic Digital Speckle Photography and Digital Image Correlation. In this work the same technique was applied to SENB specimen. We happened to find that the deformation shape of the side-necking is similar to the one of plastic region estimated by McClictock using slip line theory. Based on volume constancy of plastic deformation as well as this finding, it is expected that a linear relationship holds between the volume of plastic deformation region and the one of side-necking upon the lateral surface of a specimen. To prove the idea, a preliminary study has been performed using 3-D finite element method on a model with modified boundary layer formulation. As the result, it is shown that the idea works well with acceptable error.

Leaking water coolant into stator electrical insulation is a growing concern for the aging water-cooled generator since leaks in the generator water-cooled stator winding can affect machine availability and insulation life. But a domestic techniques of such field are insufficient and depend wholly on GE or TOSHIBA technique. Therefore this paper introduces measuring principle and developed measuring system, which has been used to detecting wet absorption. We accomplished the experiment with a stator promotion of virtue which is used in actual power plant. Also, Experimental method of generator stator winding, which is investigated into wet absorption test.

KAERI has contrived 16 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 11 kinds of spacer grid shapes among them and the others are under review in USA, EC, China, and ROK. In this study, detailed 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.

Court sport shoes is consisted of several functional parts such as soles, upper and midfoot reinforcements. Currently, intensive research for court sport shoes considering functional parts is in progress world widely, but the shoes design relies only on the view point of kinesilogy and biomechanics until now. Thus, more scientific and reliable evaluation of shoes characteristics is definitely required. In this paper, we evaluate the landing impact of court sport shoes by using finite element method. We construct a shoes-leg coupled FEM model which can simulate effectively impact in running mode. From the numerical analysis results, the designer can establish the advanced design concepts and build up the detailed design standard for the specific court sport shoes under consideration.

The design modification of bottom guide structures of the instrumented capsule, which is used for the irradiation test in the HANARO reactor, was required because of the trouble of the bottom guide arm's pin during irradiation. The previous structure with 3-pin arms was changed into the cone shape of one body. The specimens of the bottom end cap ring with three different sizes (${\Phi}68/70/72mm$) were designed and manufactured. The out-pile tests for the capsule with previous and new three bottom guide structures were performed in the one-channel flow test facilities. In order to evaluate the compatibility with HANARO and the structural integrity of the capsule, a loading/unloading, a pressure drop, a thermal performance, a vibration, and an endurance test were conducted. From out-pile test results, the capsule with the cone shape bottom guide structures was found to be more stable than the previous structure and the optimized size of the bottom guide structure selected was 70mm in diameter. It is expected that the new bottom guide structures will be applicable to all material and special capsules which will be designed and manufactured for the irradiation tests in the future.

There is few research about contact problem for a rigid surface with an arbitrary shape in SPH. The variational equation based on the virtual work principle is derived and its solution is obtained by the penalty method. It is proposed a new method that can determine the parameters for a penetration and a penetration rate used in the penalty method. The reproducing condition is adopted to correct the deficiency of kernel on the boundary. In order to calculate a penetration of particles, after checking boundary particles for deformable body boundary normal vectors were determined on the rigid surface. Numerical simulations for models which have rigid surface with an arbitrary shape were conducted to validate the proposed method in 2D. The results of those analysis represent that the contact algorithm proposed in this study works properly.

The dynamic behavior of sheet metals must be examined to ensure the impact characteristics of auto-body by a finite element method. An appropriate experimental method has to be developed to acquire the material properties at the intermediate strain rate which is under 500/s in the crash analysis of auto-body. In this paper, tensile tests of various different steel sheets for an auto-body were performed to obtain the dynamic material properties with respect to the strain rate which is ranged from 0.003/sec to 200/sec. A high speed material testing machine was made for tension tests at the intermediate strain rate and the dimensions of specimens that can provide the reasonable results were determined by the finite element analysis. Stress-strain curves were obtained for each steel sheet from the dynamic tensile test and used to deduce the relationship of the yield stress and the elongation to the strain rate. These results are significant not only in the crashworthiness evaluation under car crash but also in the high speed metal forming.

In order to determine the most suitable material system which can achieve the lightweight design and fulfill the design requirements of carbody structures of Korean Tilting Train eXpress (TTX), aluminum carbody, composite carbody, and hybrid carbody combined with aluminum and composite structures were considered in present study. The finite-element analysis was used to verity the design requirements of the TTX carbody structures with the material system being considered in the design stages. The stresses in the carbody structures and deflections of underframe against static load cases were checked as design criteria. The results show that the hybrid carbody structures are beneficial with regard to weight savings and structural integrity when compared to aluminum and composite carbody structures.

New reliability-based topology optimization method is proposed by utilizing single-loop single vector approach, which approximate searching the most probable point in the probabilistic design domain analytically, to reduce the time cost and dealing with several constraints to handle practical design requirements. To examine uncertainties in the topology design of a structure, the modulus of elasticity of the material and applied loadings are considered as probabilistic design variables. The results of design examples show that the proposed method provides efficiency curtailing the time for the optimization process and accuracy satisfying the specified reliability.

A bushing is a device used in automotive suspension systems to reduce the load transmitted from the wheel to the frame of the vehicle. A bushing is a hollow cylinder, which is bonded to a solid steel shaft at its inner surface and a steel sleeve at its outer surface. The relation between the force applied to the shaft and the relative deformation of a bushing is nonlinear and exhibits features of viscoelasticity. A force-displacement relation for bushings is important for multibody dynamics numerical simulations. For the nonlinear viscoelastic axial response, Pipkin-Rogers model, the direct relation of force and displacement, has been derived from Lianis model and the sinusoidal input was used for Pipkin-Rogers model, and the affection of displacement with frequency change was studied with Pipkin-Rogers model.

For the rational design of welded joint, it is needed to assess repeatedly the fatigue life of the joint with various dimensions and welding conditions. In this paper, an automation of repeated process of fatigue life assessment for welded cruciform joint was studied. The process is related to stress analysis in vicinity of weld-toe and fatigue life assessment based on analyzed stress distribution. With the change of design condition including dimensions and/or welding heat input, the above two works have to be performed. Using the commercial tool for system integration, ModelCenter, an automation of the repeated process for welded cruciform joint based on 2D modeling was achieved. In this automation system, data exchanges between programs regardless of commercial or in-house one work well, and parametric studies for optimal design can be performed.

Flywheel Energy Storage System(FESS) consists of a high speed flywheel with an integral motor/generator suspended on non contact bearings and in an evacuated housing. Permanent magnet machines as the FESS motor/generator are a popular choice, since there are no excitation losses which means substantial increase in the efficiency. In this paper, the structural design method of rotor retainer for a high speed motor/generator are presented.

The oscillation of the fluid caused by external forces is called sloshing, which occurs in moving vehicles with contained liquid masses, such as aircraft. cars and liquid rocket and so on. This sloshing effect could be a severe problem in vehicle stability and control. So, various baffles are used in order to reduce the sloshing. The Lagrangian, Eulerian and ALE numerical method is widely used on the analysis of sloshing presently. But, these numerical methods are needed so many CPU time. In this study, for the reduction of the sloshing analysis time, me multi.layer perceptron artificial neural network is introduced and analysis results are presented.

This paper presents applications of the objective stress rates to stress update algorithms for transient shell dynamic analysis within the context of explicit time integration. The hypo elasto-plastic materials are assumed in establishing constitutive equations. The derivation of the objective stress rates are investigated by use of the Lie derivative. Comparison results are given between the Kirchhoff and Cauchy stress formulation. The Jacobian determination algorithm proposed in this paper is presented in association with the Belytschko-Lin-Tsay shell theory. Several numerical examples are demonstrated including contact and non-contact examples, by which proposed algorithms are compared with respect to the accuracy and effectiveness.

The radiated sound pressure induced by low-velocity impact is obtained by solving the Rayleigh integral equation. This paper established the sound analysis procedure using impact analysis model. For structurally radiated noise, the sound field is directly coupled to the structural motion. Therefore the impact response should be analyzed. The impact response is computed using the spring-mass model. And the influence of damage on the sound pressure and impacted force history of laminated were investigated. The results show that both radiated sound pressure and impact force history are strongly influenced by damage on laminated.

An aerospace vehicle in flight can be exposed to random gust which may cause critical structural failure. In this paper, the failure analysis is conducted for composite wing subjected to random gust. For this, the profile of random gust is idealized as a stationary Gaussian random process and the power spectral density (PSD) of wing bending moment induced by gust is obtained. The PSD function is converted to probabilistic distributions and the failure probability during total flight time is calculated by Monte Carlo simulation.

KARI(Korea Aerospace Research Institute) has developed smart unmaned aerial vehicle(UAV) since 2002. Smart UAV has tilt rotor configuration which can take off and land vertically. For designing and developing smart UAV, it is necessary to obtain design loads. ARGON which use the panel method is multidisciplinary aircraft design program developed and modified by KARI and TsAGI. Panel method is very useful to obtain aerodynamic loads, so it have been used widely for aircraft loads analysis. For flight loads analysis, we have to prepare regulations and load conditions, and then design aerodynamic panel model, mass model and structure model. In this paper, we introduce the flight loads analysis procedure briefly, and show the smart UAV loads analysis procedure and result using ARGON.

This study is concerned with the braking distance estimation using frictional energy rate. First, steady state rolling analysis is performed, and using this result, the braking distance is estimated. Dynamic rolling analysis during entire braking time period is impratical, so that this study divides the vehicle velocity by 10km/h to reduce the analysis time. The multiplication of the slip rate and the shear stress provides the frictional energy rate. Using frictional energy rate, total braking distance is estimated, In addition, ABS(Anti-lock Brake System) is considered, and two type of slip ratios are compared, One is 15% slip ratio for the ABS condition, and the other is 100% slip ratio which leads lo the almost same braking distance as the elementary kinematic theory. A slip ratio is controlled by angular velocity in ABAQUS/Explicit, A 15% slip ratio gives the real vehicle's braking distance when the frictional energy occurred al disk pad is included. Disk pad's frictional energy rate is calculated by the theoretical approach.

In this paper, the optimum design of tire sidewall contour consisted of double plies for improving automobile maneuverability and tire durability is considered and a GUI program is developed for the purpose of the practical design. Each improvement of maneuverability and durability depends on the cord tension and strain energy distribution of tire sidewall. Satisfing trade-off method, which requires the judgment of aspiration levels, is used for the multi-objective optimization problem. Also, this paper presents the application to the practical sidewall contour design with the GUI program developed using visual Fortran.

In-arm type hydropeumatic suspension unit(ISU) is an equipment of armed tracked vehicle to absorb impact load and vibration from the irregular ground. During the operation of ISU, main piston moves forward and backward and oil flowing through damper transmits the external impact load to floating piston. Heat is generated in ISU by the oil pressure drop through the damper orifice and the friction between cylinder wall and two pistons. On the other hand, internal heat dissipatis outside via heat convection. Occurrence of high temperature can deteriorate durability of major components and basic function of ISU. And, it can cause fatal problem in the ISU life time and the sealing performance of piston rings. As well, the spring constant change of nitrogen gas that is caused by the temperature rise exerts the negative effect to the vehicle stability. Therefore, in this paper, we analyze the heat transfer analysis of the entire ISU unit, by finite element method, with the outside flow velocities 8m/s and 10m/s.

Molecular dynamics simulation of nanolithography by AFM is conducted to study nucleation of various defects, and their subsequent development and interactions as well. During nanolithography via AFM, dislocation loops are emitted along the top surface, and resourceful defect interactions such as, formation of voids chain via the motion of a jog, and creations of extended nodes and Lomer-Cottrell Lock are observed.

Fully non-local Quasicontinuum method using sub-divided region with Hermite interpolation function is proposed for simulation of carbon nanotube. Tersoff-Brenner potential is adopted for interaction of bonded atoms and also van der Waals force for non-bonded interaction. Bending of single wall carbon nanotube with chirality (20,0) and 15nm length is simulated up to 90 degree. Bending of double wall carbon nanotube with chirality (20,0) and (12,0) is simulated up to 65 degree. Bending of four wall carbon nanotube is simulated up to 45 degree.

We present a multiscale scheme which describes the dynamic pictures of atoms in the multiple length-scale systems. Large-scale atomic systems are reduced to coarse grained system by the quasicontinuum, of which the dynamic pathways are rendered by the action-derived molecular dynamics proved effective for multiple time-scale problems such as rare events. Adatom diffusions on the metal (001) surface are selected for our numerical examples. The energy barriers of the diffusions and the real dynamic trajectories of the adatoms are calculated.

We propose a 2D 'crack' element for the simulation of propagating crack with minimal remeshing. A regular finite element containing the crack tip is replaced with this novel crack element, while the elements which the crack has passed are split into two transition elements. Singular elements can easily be implemented into this crack element to represent the crack-tip singularity without enrichment. Both crack element and transition element proposed in our formulation are mapped from corresponding master elements which are commonly built using the moving least-square (MLS) approximation only in the natural coordinate. In numerical examples, the accuracy of stress intensity factor $K_I$ is demonstrated and the crack propagation in a plate is simulated.

For constructing virtual environment it is more natural to model object as deformable body than as rigid body. High accuracy of simulation of model and low-latency computational cost for real-time simulation should be guaranteed. We pre-compute Green function through finite element analysis of deformable body and it is possible to simulate deformation of body in real-time environment using Capacitance Matrix Algorithm. Also, the capacitance matrix algorithm enables to construct the haptic rendering which serves the reaction force through a haptic device. The Newmark scheme is used for the more realistic haptic rendering and dynamic simulation in real-time.

The purpose of this study is to estimate the availability of acoustic emission method to the internal leak of the valves at nuclear power plants. The acoustic emission method was applied to the valves at the site, and the background noise was measured for the abnormal plant condition. From the comparison of the background noise data with the experimental results as to relation between leak flow and acoustic signal, the minimum leak flow rates that can be detected by acoustic signal was suggested. When the background levels are higher than the acoustic signal, the method described below was considered that the analysis the remainder among the background noise frequency spectrum and the acoustic signal spectrum become a very useful leak detection method. A few experimental examples of the spectrum analysis that varied the background noise characteristic were given.

Shell finite elements are widely used for the analysis of thin section objects such as sheet metal parts, automobile bodies and et al. due to their computational efficiency. Since many of input data for finite element analysis are given as solid models or triangulated surface models, one should extract midsurface information from these input data initially and then construct shell meshes on the extracted midsurfaces. In this paper, a method of generating shell elements on midsurfaces directly from input models have been proposed. In order to construct shell meshes, the input models should be triangulated on surfaces first, and then tetrahedral elements are generated by using an advancing front method, and finally mid shell surfaces are obtained from tetrahedral meshes. Some examples are given to demonstrate the efficiency of the proposed method.

A numerical method was developed that imposes constraint condition on the order parameters in martensitic phase transformation. In the method, an amplitude function having values of 1 or 0 was multiplied to transformation rates. The merit of the method is that the imposition of the constraint condition is more straightforward than a method with Lagrangian multiplier and easy to implement in the tangent modulus method. The developed method is applied to three-dimensional finite element analyses of single and poly crystalline shape memory alloys.

There are a lot of manholes such as for water supply, sewage, telecommunication cable, traffic sign, electricity supply, and rainwater, etc. Conventional manholes installed on a road are impossible to adjust height, so that they should be entirely excavated to reinstall or repair. This entire excavation of a manhole causes too much time-consuming work, waste of resources, and obstruction of traffic. In this study, in order to solve the above mentioned problems, a cover, outer and inner parts of a manhole are integrated by gear-shaped parts located between outer and inner parts of a manhole. Mechanical design is performed to determine dimension of gear-shaped parts by Taguchi orthogonal array table. Cast molds for a gear-shaped manhole are also manufactured.

In this study, dome port bringing up for discussion where the ramjet occurs in flying it presents the tendency of distribution of thermal contour due to temperature and pressure. It is assumed that the material of ramjet is steel for the ease of result analysis. It applied matrial property which it follows by temperature and input boundary condition that changing temperature and pressure on each region by time difference for transient analysis. Thermal analysis region is decided until dome port part is separated and operate analysis in 0.5 second. Finally we draw tendency of thermal contour in ramjet dome port part by temperature and pressure.

This paper includes the comparative study between Cr plating and nitriding process with an aim at improving corrosion, wear and maintainability for KNR(Korean Next Generation Rifle) 5.56mm barrel. The endurance test was conducted to compare the performance of standard barrel, Cr plating barrel and nitriding barrel. Main activities are described as follows; optimal Cr plating and nitriding process set-up for KNR 5.56mm barrel; durability test of each barrel(20,000 rounds); salt water immersion test; dispersion, initial velocity, inner diameter data acquisition. According to the results of this firing test, Cr plating barrel is superior to standard barrel and nitriding barrel in view of corrosion, wear and maintainability

In the FE analysis of sheet metal forming, efficient results can be obtained by using shell elements rather than using solid elements. However, shell elements have some limitations to describe three-dimensional material laws. In the recent years, solid-shell element, which has only translational degree of freedom like solid element, has been presented. The assumed nature strain (ANS) and enhanced assumed strain (EAS) methods can be used to remove several solid-shell locking problems. In this paper, ANS method was used for diminish transverse shear locking and EAS method for thickness locking. Using the element, the steel pipe making process from flat plate analyzed effectively, which is including bending and welding.

This paper deals with feedback control of a hydraulic unit for direct yaw moment control, a method used to actively maintain the dynamic stability of an automobile. The uncertain parameters and complex structure naturally call for empirical modeling of the hydraulic unit, which readily results in a control-oriented model with high fidelity. The identified model is cross-validated against experimental data under various conditions, which helps to establish model uncertainty. Then, the $H_{\infty}$ optimization technique is employed to synthesize a controller with guaranteed robust stability and performance against the model uncertainty. The performance of the synthesized controller is verified using experimental results, which shows the viability of the proposed approach in a real-world application.

This paper presents an adaptive approach to control the amount of slip of the torque converter bypass clutch using its estimated friction coefficient. The proposed approach can be readily implemented using the inexpensive speed sensors currently installed in an automobile. A measurement feedback control law to drive the slip error to zero together with an adaptation law to identify the unknown friction coefficient is developed using the Lyapunov control design method. The robustness of the control and adaptation laws to parametric and/or torque uncertainties as well as the convergence of the friction coefficient are investigated. Simulation results verify the viability of the proposed control algorithm in real-world vehicle control applications.

For the safety of railway, it should be evaluated for the running safety by measuring the derailment coefficient. Although railway has run the fixed and maintained rail, some of railway is derailed. This report shows the results that performed the static load test, wheelset manufacturing for test, main line running test on the basis of the derailment theory and experience. It is executed main line test into more than 80km/h for estimating the curving performance and running safety of Gwangju EMU. As the test results, could confirm the curving performance and running safety of Gwangju EMU from the results of the wheel unloading, lateral force, derailment coefficient etc. Derailment coefficient was less than 0.8, and lateral force allowance limit and wheel load reduction ratio were enough safe.

This paper presents a simplified dynamics of railway vehicle for a tilting train simulator. The tilting train simulator has 6 electric-driven actuators and a visualization system with 1600mm-diameter dome screen. The each system shares the data by ethernet. In order to analyze the dynamics of railway vehicle and transfer the results of the analysis to the other system of the tilting train simulator in realtime base, We assumed the tilting train as a simplified rigid body model with primary and secondary suspensions. The simplified vehicle model has a 17-DOF. Through the running analysis on the tight curve with various radius, we verified the simplified vehicle model.

The effects of two important numerical procedures on the high precision structural analysis are investigated in this study. The two numerical procedures include continuous variable approximation and time integration. For the continuous variable approximation, polynomial mode functions generated by the Gram-Schmidt process are introduced and the numerical results obtained by employing the polynomial mode functions are compared to those obtained by classical beam mode functions. The effect of the time integration procedure on the analysis precision is also investigated. It is found that the two procedures affect the precision of structural analysis significantly.

In this study, a position synchronous control algorithm being applied to two-axes pneumatic cylinder driving apparatus is proposed. The position synchronous control algorithm is composed of position controller and synchronous controller. The position controller is designed to minimize the effect of several nonlinear characteristics of the driving apparatus. The synchronous controller is designed to reduce the synchronous error. The effectiveness of the proposed controller is proved by simulation results.

This work deals with the synthesis of discrete-time nonlinear controller for input time-delay existing nonlinear system and proposes a new effective method to compensate the influence of input time-delay. The controller is synthesised by using input/output linearization. Under the circumstance that input time-delay exist, controller have to produce future value that will be needed for system. On account of this reason described, a weighted average predictor of combined states is adopted. Using the discretization via Euler method, numerical simulations about Van der Pol system are performed to evaluate performance of the proposed method.

This paper describes the fault detection and failsafe logic to be used in the Electronic Stability Program (ESP). The Aim of this paper is prevention of erroneous control in the ESP. This paper introduces the fault detection logic and evaluation of residual signals. Failsafe logic consist of four redundant sub-models and they can be used for the detection of faults in each sensor (yaw rate, lateral acceleration, steering wheel angle). We presents two mathematical residual generation method ; one is the method by the average value, and the other is the method by the minimum value of the each residual. We verify a failsafe logic using vehicle test results, also we compare vehicle model based simulation results with test vehicle results.

The work presented here deals with controller design by graphical method based on proportional control ramp response. The design aims at the improvement of transient response, disturbance rejection capability, steady-state error reduction with stability preservation. The first step is to generate tracking-error curve with proportional control only and decide the added error signal shape on the error curve. The effectiveness of the proposed controller is confirmed through the simulation and experiment.

This study proposed the analysis of mass position detection and modified stiffness due to the change of the mass and stiffness of structure by using the original and modified dynamic characteristics. The method is applied to examples of a cantilever and 3 degree of freedom by modifying the mass. The predicted detection of mass positions and magnitudes are in good agreement with these from the structural reanalysis using the modified mass.

Dynamic modeling of an eccentric MEMS gyroscope is presented and the dynamic characteristics of the gyroscope are investigated with the modeling method. It is found that the eccentricity of the MEMS gyroscope affects the dynamic characteristics significantly. Different from conventional MEMS gyroscopes, the zero-rate output is significantly reduced in this gyroscope. To obtain general guidelines of the gyroscope design, dimensionless parameters are first identified and the effects of the parameters on the gyroscope performance measures are investigated.

Leak noise is a good source to identify the exact location of a leak point of underground water pipelines. Water leak generates broadband sound from a leak location and this sound propagation due to leak in water pipelines is not a non-dispersive wave any more because of the surrounding pipes and soil. However, the necessity of long-range detection of this leak location makes to identify low-frequency acoustic waves rather than high frequency ones. Acoustic wave propagation coupled with surrounding boundaries including cast iron pipes is theoretically analyzed and the wave velocity was confirmed with experiment. The leak locations were identified both by the acoustic emission (AE) method and the cross-correlation method. In a short-range distance, both the AE method and cross-correlation method are effective to detect leak position. However, the detection for a long-range distance required a lower frequency range accelerometers only because higher frequency waves were attenuated very quickly with the increase of propagation paths. Two algorithms for the cross-correlation function were suggested, and a long-range detection has been achieved at real underground water pipelines longer than 300m.

This paper introduces a link-type tracked vehicle which is developed for demining operations. The vehicle consists of three parts - front frame, rear frame and body. The front frame is connected to the rear frame by a rotational passive adaptation mechanism which is a driving mechanism of the vehicle. Additionally, the demining system which is adaptable to mobile robot is developed to clear small Anti-Personnel(AP) mines with inexplosive method. In other words, assembled rakes unearth mines by their opposite rotation to the direction of the robot. Finally, the motions of demining rakes and design parameters of the demining system are analyzed.

In this paper, we propose new parallel mechanism of a 3 dimensional biped robot whose each leg is composed of two 3-dof parallel platforms linked serially. This proposed parallel mechanism is able to move freely in the man-made environment and is applied to various fields, such as medical, welfare, and so on. And a total weight of each leg is expected to be lighter than serial linked leg. One side leg consists of a 3-dof orientation platform and 3-dof asymmetric parallel platform. The former consists of three active linear actuators and seven passive joints, and the latter of two active linear actuators, one active rotational actuator and eight passive joints. Thus, there are two kinds of parallel platforms each chain's elements and active joint's positions are different for the biped robot to move freely like a serial link without the kinematics constraints. The effectiveness and the performance of the proposed parallel mechanism and locomotion trajectory are shown in computer simulations with a 12-DOF parallel biped robot.

The outlet feeder pipe thinning in a PHWR (Pressured Heavy Water Reactor) is caused by high pressure steam flow inside the pipe, which is a well known degradation mechanism called FAC (Flow Assisted Corrosion). In order to monitor the degradation, the thickness of the outlet bends closed to the exit of the pressure tube should be measured and analyzed at every official overhaul. This paper develops an automatic feeder pipe inspection system that can minimize the irradiation dose by automating the measurement process. The robot can move by itself on the feeder pipe by using an inch worm mechanism, which is constructed by two gripper bodies that can fix their body on the pipe and one extendable and retractable body connected the two gripper bodies to move forward and backward.

In this paper, an intelligent fault diagnosis system is proposed for induction motors through the combination of feature extraction, genetic algorithm (GA) and neural network (ANN) techniques. Features are extracted from motor vibration signals, while reducing data transfers and making on-line application available. GA is used to select most significant features from whole feature database and optimize the ANN structure parameter. Optimized ANN diagnoses the condition of induction motors online after trained by the selected features. The combination of advanced techniques reduces the learning time and increases the diagnosis accuracy. The efficiency of the proposed system is demonstrated through motor faults of electrical and mechanical origin on the induction motors. The results of the test indicate that the proposed system is promising for real time application.

The performance of a mixed $H_{\infty}/H_2$ design with pole placement constraints based on robust vibration control for a piezo/beam system is investigated. The governing equation of motion for the piezo/beam system is derived by Hamilton's principle. The assumed mode method is used to discretize the governing equation into a set of ordinary differential equation. A robust controller is designed by $H_{\infty}/H_2$ feedback control law that satisfies additional constraints on the closed-loop pole location in the face of model uncertainties, which are derived for a general class of convex regions of the complex plane. These constraints are expressed in terms of linear matrix inequalities (LMIs) approach for the multiobjective synthesis. The validity and applicability of this approach for vibration suppressions of SMART structural systems are discussed by damping out the multiple vibrational modes of the piezo/beam system.

In this paper, modeling methods for the structural dynamic analysis employing single reference frame are presented and their modal and transient analysis results are compared. The geometric stiffening effects often occur when structures undergo large overall motion. These effects were considered in several structural previous modeling methods but the role of reference frame has never been scrutinized. In this study, modeling methods employing single reference frame are presented, and their numerical results are compared. The results show that discrepancy between the two modeling methods increases as the eccentricity of the structural system and the magnitude of the large overall motion increase.

In many applications, piezoelectric motor has been used in the areas where excellent controllability, fine position resolution are required or magnetic field noise should be eliminated. However the piezoelectric motor has two major demerits. One is difficulty in maintaining vibration amplitude constantly with temperature rise and wear, and the other is heat generation induced by dielectric and mechanical loss. In this thesis, piezoelectric motor to overcome these problems is proposed. And proposed piezoelectric motor is operated using momentum exchange between bimorph and rotor. To maximize steady state velocity and static torque of proposed motor, the guideline is established using two bimorph models. Then the guideline is partially proved by comparison between simulation and experiment. There was no heat generation in the few hours of experiments.

In this paper, a sliding mode control based on disturbance observer is proposed to attenuate disturbance responses in an active magnetic bearing system, which is subject to base motion. An algorithm for exactly decoupling the disturbance estimation dynamics from the sliding mode dynamics is developed. It is also shown that the proposed method preserves the robustness of the sliding mode and asymtotically achieves zero regulation error, in the presence of external disturbances and parametric uncertainties. The proposed control is applied to a 2-DOF active magnetic bearing system subject to base motion. The feasibility of the proposed technique is illustrated, and the results of an experimental demonstration are shown.

In this paper, an angular self-sensing algorithm is proposed and implemented to a Lorentz force type integrated motor-bearing system. It is based on the principle that the flux linkages of stator windings, calculated from the voltage and torque control current, are the functions of the rotor angle. The tracking angular position error is proven to vanish using the Lyapunov stability method, and the experimental results show that the initial error decays within about 5 seconds. It is found that the resolution of the algorithm remains about 1º over the speed range of 100 to 1000 rpm.

In this study, the new prosthetic hand which can perform sensory feedback is developed. Although the continuous development of current prosthetic hands, they have some drawbacks. The drawbacks are the low grasping capabilities, the lack of sensory information given to the user, and so on. Despite the improvement of the function of fingers, the sensory information problem still remains undeveloped. By using differential gear and underactuated mechanical design, it can be solved with the minimum additional weight. And it will be applied to all of the prosthetic hands.

In this paper, a blackbox approach is carried out to model the nonlinear dynamic bushing model. One-axis durability test is performed to describe the mechanical behavior of typical vehicle elastomeric components. The results of the tests are used to develop an empirical bushing model with an artificial neural network. The back propagation algorithm is used to obtain the weighting factor of the neural network. Since the output for a dynamic system depends on the histories of inputs and outputs, Narendra's algorithm of 'NARMAX' form is employed in the neural network bushing module. A numerical example is carried out to verify the developed bushing model.

This paper presents a three dimensional modeling and simulations of operation and running of a wheel loader using the ADAMS program. A wheel loader consists of a bucket, a boom, a crank, a front frame, a rear frame, a bucket cylinder, two boom cylinders, two steering cylinders, nine spherical joints, six universal joints, five translation joints, three inline joints, a revolute and a fixed joint. Judging from the actual degrees of freedom of the wheel loader, proper kinematic joints are selected to exclude redundant constraints in the modeling. Through the running simulation over a bump with the three dimensional modeling, the joint reaction forces are calculated.

The DC Motor of a vehicle may cause noise and vibration due to high-speed revolution, which can make a driver feel uncomfortable. There have been various studies that attempted to solve these problems, mostly focusing on the causes of noise and vibration and the means of preventing them. The CAE methodology is more efficient than a real test for the purpose of looking for various design parameters to reduce the noise and vibration of the DC motor. In this study, a design process for reducing brush noise is presented with the use of a computer model, which is made by using a multi-body dynamics program (DADS). The design parameters to reduce the brush noise and vibration were proposed using a computer model. They were used to reduce the noise and vibration of a DC motor and verified by the test results of a fan DC motor in a vehicle. This method may be applicable to various DC motors.

This paper presents the inverse dynamic analysis of the wheel loader manipulator based on the experimental data. A three dimensional rigid multi-body model of the wheel loader manipulator was built up. The inverse dynamic analysis for the typical operation mode was carried out by the ADAMS program. In order to verify the analysis result with the measured one, the hydraulic pressure and displacements of the cylinders were measured and the inverse dynamic analysis was carried out using experimental data. From the results of the analysis and measurement, it was concluded that the computational driving force showed good agreement with the measured one.

In this paper, the activity of Japanese researchers relating multibody dynamics is introduced. The author stayed at the CCR(Center for Collaborative Research) in the university of Tokyo as a visiting scholar during December 10, 2003 to March 9, 2004. The information obtained from several universities (University of Tokyo, Sophia university, Nihon university, Tsukuba university), several industries and research institutes (JR Central in Nagoya, Toyota Research Center in Nagoya., Subway department in Tokyo, and JARI at Tsukuba), several research groups (JSCM, JSME, JSAE) was summarized with several photos. Professor Yoshihiro Suda's laboratory at the university of Tokyo, professor Yoshiaki Terumichi's laboratory at the Sophia university, and professor Ichiro Kageyama's laboratory at the Nihon university are introduced.

This research proposes an implementation method of linearized equations of motion for multibody systems with closed loops. The null space of the constraint Jacobian is first pre multiplied to the equations of motion to eliminate the Lagrange multiplier and the equations of motion are reduced down to a minimum set of ordinary differential equations. The resulting differential equations are functions of all relative coordinates, velocities, and accelerations. Since the coordinates, velocities, and accelerations are tightly coupled by the position, velocity, and acceleration level constraints, direct substitution of the relationships among these variables yields very complicated equations to be implemented. As a consequence, the reduced equations of motion are perturbed with respect to the variations of all coordinates, velocities, and accelerations, which are coupled by the constraints. The position, velocity and acceleration level constraints are also perturbed to obtain the relationships between the variations of all relative coordinates, velocities, and accelerations and variations of the independent ones. The perturbed constraint equations are then simultaneously solved for variations of all coordinates, velocities, and accelerations only in terms of the variations of the independent coordinates, velocities, and accelerations. Finally, the relationships between the variations of all coordinates, velocities, accelerations and these of the independent ones are substituted into the variational equations of motion to obtain the linearized equations of motion only in terms of the independent coordinate, velocity, and acceleration variations.

Multibody dynamics formulation has been developed based on relative cartesian coordinates for subsystem analysis. Relative cartesian coordinates are defined with respect to a reference body of a subsystem. Relative cartesian formulation inherits the same merits of absolute cartesian formulation, such as generality and easy implementation. Two methods have been applied. One is Largrange Multiplier Elimination method and the other is independent coordinate method. A 1/4 car simulation has been carried out to verify the formulations. Since both methods provide identical results, it proves the validity of the formulation.

Vibration characteristics of the FIV test loop for the Flow-Induced Vibration(FIV) study of a PWR partial(5x5) fuel assembly are investigated by the Finite Element(FE) analysis and the modal test. For the FE analysis, 3-D beam element is used for the pipes and the test section and mass element used for the valves and flanges. The 'U' restrainer stiffness determined by numerical simulation is used for the FE model. The result of the FE analysis is compared with that of the modal test. The higher mode similarity between the test and analysis is observed in a few low modes. After that, the mode similarity reduce as the mode goes high. It is concluded that the first to the third vibration modes are observed at the lower parts of the 6 inches restoring line, followed by a local mode at the test section, and the natural frequencies of the modes are 22.4 Hz, 26.0 Hz, 27.5 Hz and 31.4 Hz.

The paper deals with the dynamic stability of a cantilevered Timoshenko beam on partial elastic foundations subjected to a follower force. The beam is assumed to be a Timoshenko beam with a concentrated mass taking into account its rotary inertia and shear deformation. Governing equations are derived by extended Hamilton's principle, and FEM is applied to solve the discretized equation. Critical follower force depending on the attachment ratios of partial elastic foundations, concentrated mass and rotary inertia of the beam is fully investigated.

As the tasks of robots become more diverse, some complicated tasks have come to require force and position hybrid control. A compliant device can be used to control force and position simultaneously by separating the twist of the robot's end effector from the twist of compliance and freedom by using stiffness mapping of the compliant device. The development of a fuzzy gain scheduling scheme of control for a robot with a compliant device is described in this paper. Fuzzy rules and reasoning are performed on-line to determine the gain of twists based on wrench error and twist error and twist of compliance and twist of freedom ratio. Simulation results demonstrate that better control performance can be achieved in comparison with constant gain control.

AFM tip has been used for surface profiling with a fine resolution, but there is a barrier to improve its performance because of the low aspect ratio. Many researchers have solved this problem with attaching carbon nanotube (CNT) to Si-tip. In this paper, we proposed the aligner system that composed of dual type stage system, and these stages could attach a carbon nanotube to tungsten-tip in vacuum condition. We used tungsten tip instead of Si-tip because of its conductivity. The aligner system proposed in this paper has 10 degree-of-freedom that 3 in the first stage and 7 in the second stage. With picomotors and piezotube, the first stage has the resolution about several tens of nm and the second stage has a resolution about a nm. We experimented on characterization of Nano Aligner System and operated picomotors in SEM environment.

A robot manipulator is usually operated in two modes: free motion and constraint motion according to the fact whether the robot comes into contact with the environment or not. At the moment of contact, impact occurs, and sometimes, it can possibly degrade the robot's performance such as vibration and at worst, shortens its lifetime. In this article, a new proposed algorithm is described by introducing a command signal modification method on the basis of impedance control and a validity of the proposed algorithm is demonstrated by showing the simulation study.

Railway wheel/rail contact conditions have an influence on dynamic behavior of rolling stock. If there are problems of incompatibility between wheel and rail, damages like wheel wear, wheel spalling, rail wear, etc are occurred. Especially wheel and rail profiles are important factor of vehicle curving performance, so compatibility study between wheel and rail has to be carried out preferentially, In this study, we have analyzed the compatibility between wheel and rail of KNR conventional line to improve the maintenance efficiency of wheel and rail. Thus we showed the results relating to wheel/rail geometric contact, vehicle running performances as the change of wheel/rail combination.

Recently, the trend in architectural forms has been toward taller and larger building. The building materials, therefore, are getting larger and heavier as wall. Most of the construction projects are, however, dependent on outdated equipment and human resources. Construction processes up to now face a number of problems, including dangerous work, high construction cost and heterogeneous construction quality. In various construction sites, automation in construction has been introduced to address these problems. This paper proposes a human-machine cooperative system in the construction site; the system utilizes construction of a curtain wall in tall buildings. The use of automation system at construction sites can reduce the need of human involvement. Construction time and cost can be reduced as well. The most important aspect of the use of automation system at construction sites is prevention of accidents.

It is necessary to detect the tilt angle for control or monitoring of various systems such as ship, aircraft recreation facilities etc as well as bridge or building systems. However the electronic liquid charge type tilt sensor which is one of typical tilt sensor has many problems. Those are remarkably slow response time and limited mounting condition because or liquid viscosity coefficient and inertia etc. In this study we propose a tilt angle instrumentation method using piezoelectricity acceleration sensor. The method can he applied on moving mount We verified the validity of the method through experiment.

Slip displacement is brought into focus to study the tube fretting wear that occurs on the contact by the springs. An oscillating tube was in contact with plate support springs. The contact condition was varied as normal force 5 N, and gaps of 0.1 and 0.2 mm in the experiment. The oscillation range of the tube was also varied as 0.2, 0.3, 0.4 and 0.7 mm. Formulas for predicting the slip displacement range were derived in terms of the vibration amplitudes measured during the tube oscillation. It was found that the slip displacement in transverse direction was much higher ($720{\sim}33000$ times) than that in axial one. This resulted in the severer wear on the contact suffered from transverse slip.

Effects of ion implantation and deposition on the tribological properties of DLC film as a function of implanted energies and process times were investigated. TiC ions were implanted and deposited on the Si-wafer substrates followed by DLC coating using ion beam deposition method. In order to study tribological properties such as friction coefficient and behavior of DLC film on the modified surface as a function of implanted energies and process times, we used a ball-on-disc type apparatus in the atmospheric environment. From results of wear test, as the implanted energy was increased, the friction coefficient was more stable below 0.1.

The cylinder without seal has a piston with air bearing which is partly cylindrical and conical shape. The description of system geometry is deviation by the flow rate equation. 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. As a result of leakage flow test, it is evaluation to air bearing in sealless cylinder.

The role of transfer films formed during sliding of polymer composites against steel counterfaces was studied in terms of the tribological behaviors of composites. Four kinds of composites were included in this study: (1) unfilled PPS, (2) PPS+2%CuO, (3) PPS+2%CuO+5% carbon fiber (CF), and (4) PPS+2%CuO+15%Kevlar. The filler material CuO was in nanoscale particulate form and the reinforcing material was in the form of short fibers. The composites were prepared by compression molding at $310^{\circ}C$ and sliding tests were run in the pin-on-disk sliding configuration. The counterface was made of tool steel hardened to 55-60 HRC and finished to a surface roughness of 0.09-0.10 ${\mu}m$ Ra. Wear tests were run for 6 hrs at the sliding speed of 1 m/s and contact pressure of 0.65 MPa. Transfer films formed on the counterfaces during sliding were investigated using AFM and SEM. The results showed that as the transfer film became smooth and uniform, wear rate decreased. PPS+2%CuO+15%Kevlar composite showed the lowest steady state wear rate in this study and its transfer film showed the smoothest and the most uniform characteristics. The examination of worn surfaces of PPS+2%CuO composite using X-ray area scanning (dot mapping) showed back-transfer of steel counterface material to the polymer pin surface. This behavior is believed to strengthen the polymer pin surface during sliding thereby contributing to the decrease in wear rate.

A pneumatic cylinder used to sliding seal which seal the element one to another in relative motion. It is difficult to accuracy control because of adhesion phenomenon. To confirm this phenomenon, it is carried to friction force test and analysis for bellows type rodless cylinder. In the rodless cylinder of this type, friction force test is very important. Through the theoretical analysis of shock absorber obtained the optimal design of bellows type rodless cylinder. As the result of the friction force test, LM Guide type is suitable for work under low friction force.

We present high aspect ratio 100nm-scale nickel stamper using e-beam lithography process and Cr/Qz mask for the injection molding process of nano grating patterns. Conventional photolithography blank mask (CrON/Cr/Qz) consists of quartz substrate, Cr layer of UV protection and CrON of anti-reflection layer. We have used Cr/Qz blank mask without anti-reflection layer of CrON which is non-conductive material and ebeam lithography process in order to simplify the nickel electroplating process. In nickel electroplating process, we have used Cr layer of UV protection as seed layer of nickel electroplating. Fabrication conditions of photolithography mask using e-beam lithography are optimized with respect to CrON/Cr/Qz blank mask. In this paper, we have optimized e-beam lithography process using Cr/Qz blank mask and fabricated nickel stamper using Cr seed layer. CrON/Cr/Qz blank mask and Cr/Qz blank mask require optimal e-beam dosage of $10.0{\mu}C/cm^2$ and $8.5{\mu}C/cm^2$, respectively. Finally, we have fabricated $116nm{\pm}6nm-width$ and $240nm{\pm}20nm-height$ nickel grating stamper for the injection molding pattern.

In the present study, the effects of orthogonal ribs on structural warpage during injection molding process were investigated. Basic ribbed models for the evaluation of degree of warpage were introduced and designed. Injection molds for these models are manufactured based on the full 3D CAD/CAM technology and specimens are prepared for experiment. Numerical analysis using commercial plastic injection molding analysis software was also performed to compare the results with experimental ones. The variations of materials and parameters such as injection time, mold temperature, melt temperature, holding time were considered in the present work. It was shown that orthogonal ribs have significant effect on the reduction of warpage during the injection molding process.

A new method is proposed to fabricate a reusable qualtz master with order of 100 nm dot pattern on its surface. Some fabrication conditions such as dose are investigated to find optimal condition. This reusable qualtz master is used directly as a stamper to injection mold the dot patterns. Polycarbonate and Polyoxymethylene are used as molding materials and the effect of the mold temperature is also investigated to see the moldabilty of the injection molding for very fine dot features.

In this paper, the effects of grinding conditions on the surface roughness of $ZrO_2$ ferrule applying to the optical fiber connector were investigated. The mesh number of grindstone, the grinding speed and the time schedule during grinding were regarded as main cutting parameters that have an effect on the surface roughness. And then, optimal combinations for chamfering grinding were obtained by using the Taguchi method. In addition, analytic values for maximum surface roughness ($R_{max}$) estimated by the theoretical equations which were derived from the formative model of surface roughness on the chamfering grinding were compared with those of the experiments.

A numerical method and algorithms is proposed to perform optimization of non-linear response structures. An analytical and numerical method based finite element method is also proposed for the transformation of non-linear response into linear response. Loads transformed from this method are defined as the equivalent linear loads. With the loads and the transformed response, linear static optimization is performed for nonlinear response structure with geometric and/or material non-linearity. The results of the optimization are compared with them of typical non-linear response optimization using finite difference method. The proposed method is very efficient and derives good solution.

The structural optimization has been carried out in the continuous design space or in the discrete design space. Generally, available designs are discrete in design practice. But methods for discrete variables are extremely expensive in computational cost. In order to overcome this weakness, an iterative optimization algorithm was proposed for design in the discrete space, which is called as a sequential algorithm using orthogonal arrays (SOA). We focus to verify the fact that the local solution can be obtained throughout the optimization with this algorithm. The local solution is defined in discrete design space. Then the search space, which is the set of candidate values of each design variables formed by the neighborhood of current design point, is defined. It is verified that a local solution can be founded by moving sequentially the search space. The SOA algorithm has been applied to problems such as truss type structures. Then it is confirmed that a local solution can be obtained using the SOA algorithm

All the loads in the real world are dynamic loads and it is well known that structural optimization under dynamic loads is very difficult. Thus the dynamic loads are often transformed to the static loads using dynamic factors. However, due to the difference of load characters, there can be considerable differences between the results from static and dynamic analyses. When the natural frequency of a structure is high, the dynamic analysis result is similar to that of static analysis due to the small inertia effect on the behavior of the structure. However, if the natural frequency is low, the inertia effect should not be ignored. Then, the behavior of the dynamic system is different from that of the static system. The difference of the two cases can be explained from the relationship between the homogeneous and the particular solutions of the differential equation that governs the behavior of the structure. Through various examples, the difference between the dynamic analysis and the static analysis are shown. Also the optimization results considering dynamic loads are compared with static loads.

Though the standard element density-based topology optimization method has been applied for the optimal design of multiphysics systems, some theoretical problems, such as material interpolation, undershoot temperature prediction, and unstable elements, still remain to be overcome. The objective of this investigation is to present a new topology optimization formulation based on the element connectivity parameterization (ECP) in order to avoid the numerical problems in multiphysics system design and improve optimization results. To show the validity of the proposed approach, the designs of an optimal thermal dissipation and an electro-thermal-compliant actuator were considered.

A micro-pressure wave is generated by the high-speed train which enters a tunnel, and it causes explosive noise and vibration at the exit. It is known that train speed, train-tunnel area ratio, nose slenderness and nose shape mainly influence on generating micro-pressure wave. So it is required to minimize it by searching optimal values of such train shape factors and tunnel condition. In this study, response surface model, one of approximation models, is used to perform optimization effectively and analyze sensitivity of design variables. Owen's randomized orthogonal array and D-optimal Design are used to construct response surface model. In order to increase accuracy of model, stepwise regression is selected. Finally SQP(Sequential Quadratic Programming) optimization algorithm is used to minimize the maximum micro-pressure wave by using built approximation model.

Recent engineering process requires fast development and manufacturing of the products. This paper mainly discusses the process of rapid product development (RPD) from the reverse engineering to the optimal design. A laser scanning system scans a product and the efficient data processing method reduces the scanned point data. The reduced (scanned) points model is transformed to a finite element model without the construction of a CAD model. Since CAD modeling is a time-consuming work, skipping this step can save much time. This FE model is updated from the result based on the structural characteristics from modal test of the real model. For FE model updating, Response Surface Method is adopted. Finally, the updated FE model is optimized using the reliability-based topology optimization, which is developed recently. All these processes are applied to the design of an upper part model of a cellular phone.

In the proto design stage of a new car, the performance of an occupant protection system is often evaluated by CAE instead of the real test. CAE predicts and recommends the appropriate design values hence reducing the number of the real tests. However, the existing researches using CAE in predicting the performances do not consider the uncertainties of parameters, in which inconsistency between the actual test results and CAE exists. In this research, the optimization procedure of a protection system such as airbag and load limiter is suggested for the frontal collision. The DACE modeling known as Kriging interpolation is introduced to obtain the meta model of the system followed by the tabu search method to determine a global optimum. Finally, the distribution of a suggested design is determined through the Monte-Carlo Simulation.

A proportional pressure control valve has a nonmagnetic ring which is inserted in between a coil and plunger and it can get attraction force in proportion to input current by an influence of control cone. Therefore, a proportional pressure control valve is applied to a servo system as substitution of servo valve and an on-off solenoid valve widely because control of a high level is possible and pollution level is low. The purpose of this study is to develop domestic model of a proportional pressure control valve, and a test model was designed and manufactured through valve system analysis and finite element analysis. And comparison between results of theoretical analysis and static / dynamic characteristics test were carried out on a manufactured test model, and it was confirmed that it has excellent performance in comparison with other foreign products.

An efficient boundary-based technique is developed for addressing shape design sensitivity analysis in supercavitating flow problem. An analytical sensitivity formula in the form of a boundary integral is derived based on the continuum formulation for a general functional defined in potential flow problems. The formula, which is expressed in terms of the boundary solutions and shape variation vectors, can be conveniently used for gradient computation in a variety of shape design in potential flow problems. While the sensitivity can be calculated independent of the analysis means, such as the finite element method (FEM) or the boundary element method (BEM), the FEM is used for the analysis in this study because of its popularity and easy-touse features. The advantage of using a boundary-based method is that the shape variation vectors are needed only on the boundary, not over the whole domain. The boundary shape variation vectors are conveniently computed by using finite perturbations of the shape geometry instead of complex analytical differentiation of the geometry functions. The supercavitating flow problem is chosen to illustrate the efficiency of the proposed methodology. Implementation issues for and optimization procedure are addressed in this flow problem.

This paper presents the adjoint variable design sensitivity analysis for thermal systems considering both conduction and convection heat transfer. Both nodal temperature and total heat flow are considered to be objective functions and design sensitivity formulas are derived for each case. For the case of convection heat transfer, the adjoint analysis is carefully proceeded to obtain a precise result. A topology optimization example is examined for a simple planar square plate in order to design a heat exchanger as verification.

This paper presents a new approach regarding thermal characteristics associated with a design of the high efficiency motor. Electrical conduction materials, such as coil and aluminum embedded in the core generate high heat exerting negative influence on both lifetime and performance of machine. Thus, it is necessary to design high efficiency motor considering heat transfer in order to improve motor performance and to be protected against overheating. In this paper, firstly, numerical analysis of electromagnetic field is carried out by the nonlinear transient finite element method (FEM). Secondly, the linear static FEA of magneto-thermal field is implemented by applying source current computed by the nonlinear transient analysis. FE results are validated in terms of electromagnetics and heat transfer by experiments. And then, the pseudo-transient topology optimization using a multi-objective function is performed. The proposed method is applied to a squirrel cage single-phase induction motor of the scroll compressor.

Ceramic ferrules which are major parts of the fiber optic connectors are requiring high precisions in grinding. After grinding, it is the problem that subsurface damages cannot be removed. The objective of this study was to analyze the grinding force and the associated stress generated in a ceramic ferrule during cylindrically external griding using finite element analysis(FEA). A two-dimensional finite element model was constructed with the grinding parameters and the mechanical properties of the ferrule as input variables. The size of the geometric model was the same with ceramic ferrule. The experimental results are investigated by SEM photograph and compared with the results from FEM. The result of FEA showed a good agreement with that of experiment.

Many mid-sized companies in Korean automotive industry have attempted to solve the lack of human power, to control the quality of components, to improve the production rate, and to innovate the manufacturing line. The goals of this study are to analyze the production rate of an automotive component manufacturing line using simulation software, to construct a Flexible Automation Subassembly (FAS) system and to suggest an optimized layout design using FAS line. In this research, the simulation model for manufacturing line was developed and used the realistic data (production planning, component type, working order, process time, queue time, line rules, etc) of a medium sized company in Korean automotive industry. To complete this research, a simulation software 'ARENA' was used.

In this study, we have developed the simulation tool in order to investigate driving trajectory of AGV for container transport. AGV for container transport is different from the indoor AGV in that it is a large size structure at being loaded the weight of 40 ton. And AGV for container transport is applied to front wheel steering, rear wheel steering, all wheel steering, and crap steering. Therefore, we have developed the simulation tool considering dynamic problems and center of turning in accordance with four way of steering modes. Throughout some computer simulations, we have confirmed that this tool is useful to analysis dynamic problems and to calculate minimum radius of turning for large size vehicles.

In silicon wafer manufacturing process, the grinding process has been adopted to improve the flatness of wafer. The grinding of wafer is usually used by the infeed grinding machine. The infeed grinding machine has been depended on imports. Therefore, it is necessary to develop the infeed grinding machine because the demand of the infeed grinding machine is increasing more and more. This paper describes the technologies of infeed grinding machine and intend to introduce the studies in the development of the intelligent grinding system for grinding of wafer. The air bearing spindle for the infeed grinding machine was developed by domestic technologies and the grinding part design of the intelligent grinding system for wafer grinding was completed.

In this study, the ultra precision centerless grinder for ferrule grinding was designed. The centerless grinder is composed of the high damping concrete bed, grinding wheel spindle unit, regulating wheel spindle unit, feeding table and dressing unit. For a newly developed centerless grinder, hydrostatic system with high precision feeding and high stiffness was proposed. The grinding and regulating wheel spindle units were composed of hydrostatic spindle and feeding table was hydrostatic table. The prototype of hydrostatic table was manufactured and tested.

Today the demand of the optical communication components has been increased. Zirconia Ferrule has become the one of the most important elements because it determines transmission efficiency and quality of information in the optical communication system. Grinding is the major process in the ferrule manufacturing process which require high processing precision. In this reseach, specially designed spindle, chucking system, loading & unloading system and cooling system, as a supporting experimental equipment for development of an Intelligent Coaxial Grinding System (ICGS) for Zirconia Ferrule processing, is developed. We are also analized the adaptability of ICGS in practical use, through the way of evaluation for the performance of the each systems above.

As consumer in optics, electronics, aerospace and electronics industry grow, the demand for ultra precision aspherical surface lens increases higher. To enhance the precision and productivity of ultra precision aspherical surface micro lens, the following specification of ultra precision grinding system is required: the highest rotational speed of the grinder is 100,000rpm and its turning accuracy is $0.1{\mu}m$, positioning accuracy is $0.1{\mu}m$. The development process of the grinding system for the ultra precision aspherical surface micro lens for optoelectronics industry is introduced. In the work reported in this paper, an intelligent grinding system for ultra precision aspherical surface machining was designed by considering the factors affecting the surface roughness and profiles accuracy. An aerostatic form was adopted to build the spindle of the workpiece and the spindle of grinder and ultra precision LM guide way was adopted in this system.

As recently optical communication industry is developed, request of optical communication part is increased. Ferrule is very important part which determines transmission efficiency and quality of information in the optical communication part. Most of ferrule processes are grinding which request high processing precision. The ultra precision centerless grinding machine for ferrule grinding was designed. The centerless grinding machine is composed of the high damping bed, grinding wheel spindle unit, regulating wheel spindle unit, feeding table and dressing unit. Reliability prediction was very important for the high quality design. In this study, centerless grinding machine was predicted reliability.

RMQC, rail mounted quay crane, lifts and moves large container by means of hoist and trolley motion on the dock. Our company is trying to develop advanced RMQC applying the concepts of elevator hoist and container conveyer for the automation and high efficiency in handling the boxes. Prior to the development of new products, it is necessary to check the capacity of the new systems using simulation program. The program simulates container-handling rates and gives some design factors for the new-type cranes.

In a foaming process of microcellular plastics (MCPs) with a injection molding, research on the viscosity change that occurs when the gas is injected to the polymer has received little attention despite its importance. The purpose of this paper is to provide the basic data required to determine the processing condition by measuring viscosity changes against the gas injection rates of the blowing agent, and to verify the influence of the viscosity change on the flow condition of polymer inside the mold at the injection process.

Welding characteristics of STS304 stainless steel and Inconel 600 using a continuous wave Nd:YAG laser are experimentally investigated. Alloy 600 being used in steam generator tubing of pressurized water reactor(PWR) exposed to some corrosion environment, stress corrosion cracking can occur on this material. Presented here are the results from a series of experiments in which dissimilar metal welds were made using the gas tungsten arc welding process with pure argon shielding gas. But It is well known that solidification cracking susceptibility of austenitic alloys depends on the solidification temperature range and amount/distribution of solute rich liquid that exists at the terminal stages of solidification. An experimental study was conducted to determine effects of welding parameters and to optimize those parameters that have the most influence on eliminating or reducing the extent welding zone formation at dissimilar metal welds.

Because of the various length condition of carbon nanotube, it is very confined the application area of the single CNT tip. In this paper, it was proposed the cutting technique of single CNT tip and the relationship between the etched volume and the amount of the applied charge. It is possible to control the length of single CNT tip arbitrary using this technique. The etched volume and length in the single CNT tip can be predicted with the amount of the applied charge. It is very easy to make a single CNT tip with proper length using this technique and to make nanotweezer that was composed two single CNT tips.

The vertical roller mill is the important machine grinding and mixing various crude materials in the manufacturing process of portland cement. Vertical roller mill is subjected to the cyclic bending stress by rollers load. It demands $4{\times}10^7$ cycle but has $4{\times}10^6{\sim}8{\times}10^6$ cycle. It fractures at the edge of grinding path of outside roller. The repair expense for it amounts to 30% of total maintenance. Therefore, this study shows optimal design for vertical roller mill using DOE and neural network.

In recent years, demands on microelectrode have been greatly enhanced because of its potential applications and mass production of microelectrodes is needed. An electrochemical fabrication is used as an method for the simple and cheap fabrication of multi microelectrodes. In this paper, one dimensional microelectrode array is used for fabricating of multi electrodes. A diffusion layer which is formed near the electrode surface has an effect on the shape error of multi microelectrodes. The optimal distance between electrodes to minimize shape errors of multi electrodes is investigated. Multi microelectrodes which has several tens of and hundreds of micrometer in diameter are fabricated at a time.

The tooth profile of helical gear that is manufactured by the rack cutter consists of the involute curve and the trochoid curve. However, gear designers are very hard to calculate the exact profile because it needs very complex information about the gear manufacturing. Therefore, the purpose of this study is to develop the automatic generation program for the helical gear using the Solid Works API. First, involute and trochoide coordinates by the rack cutter are calculated. Using the data, Visual Basic program for the helical gear model is coded. This work gives us the quick helical gear modeling and can be used as the modeling for the finite element analysis.

Recently the high speed end-milling processing is demanded the high-precise technique with good surface rougj1ness and rapid time in aircraft, automobile part and molding industry. The working factors of high speed end-milling has an effect on surface roughness of cutting surface. Therefore this study was carried out to analyze the working factors to get the optimum surface roughness by design of experiment. From this study, surface roughness have an much effect according to priority on Spindle speed, feed rale, hardness and axial depth of cut By design of experiment, it is effectively represented shape characteristics of surface roughness in high speed end-milling And determination($R^2$) coefficient of regression equation had a satisfactory reliability of 89.7% and regression equation of surface roughness is made by regression analysis.

Laser beam welding is increasingly being used in welding of structural steels. The laser welding process is one of the most advanced manufacturing technologies owing to its high speed and deep penetration. The thermal cycles associated with laser welding are generally much faster than those involved in conventional arc welding processes, leading to a rather small weld zone. Experiments are performed for Inconel 600 plates changing several process parameter such as laser power, welding speed, shielding gas flow rate, presence of surface pollution, with fixed or variable gap and misalignment between plate and plate, etc. The follow conclusions can be drawn that laser power and welding speed have a pronounced effect on size and shape of the fusion zone. Increase in welding speed resulted in an increase in weld depth/ aspect ratio and hence a decrease in the fusion zone size. The penetration depth increased with the increase in laser power . Welding characteristics of austienite Inconel 600 using a continuous wave Nd:YAG laser are experimentally investigated. This paper describes the weld ability of inconel 600 for machine structural use by Nd:YAG laser.

Generally in the case of parts used for precision products, tolerance of parts is very small. So inaccuate size of molding parts generates serious problems. Therefore, it's necessary to secure data about shrinkage on each condition or study about manufacturing process which reduces shrinkage. To apply MCPs to manufature of plastic product, this paper verifies how the amount of gas and Talc can affect to cell-morphology, and examines the relation between shrinkage and cell-morphology by using ASTM specimen formed by MCPs process.

The purpose of this study is to investigate the supercooling improvement of TMA30wt% clathrate when the chloroform is added to it. For this purpose, phase change temperature and supercooling are measured and evaluated experimentally in heat source of $-7^{/circ}C$. The results show that phase change temperature and supercooling are improved. From the results, this research can provide and important data for the low temperature thermal storage.

To develop and improve a switchgear, the prediction of the pressure rising within the switchgear is very important. This study investigates the pressure rising characteristics of switchgear in order to evaluate the result of arc fault test. The pressure rising time at the four points of measurement calculated by CFD is well accord with the experimental results. The maximum pressure within the switchgear estimated by CFD is about 1.0bar, the pressure from experiment is 0.7 bar. The results of this study are able to be used to improve the performance of existing switchgear and to develop a new type switchgear.

The physical model considered here is a horizontal layer of fluid heated below and cooldabove with a heat-generating conducting body placed at the center of the layer. The body genrates a constant amount of heat as initial condition. Two-dimensional solution for unsteady natural convection is obtained using an accurate and efficient Chebyshev spectral methodology for various of Rayleigh number from $10^3$ to $10^6$. Multi-domain Technique is used to handle heat-generating conducting body. The results for the case of heat-generating body are also compared to those of adaibatic body.

The inverse heat conduction problem (IHCP) is a problem of estimating boundary condition from temperature measurement at one or more interior points. Neural networks are general information processing systems inspired by the connectionist theory of human brain. By properly training the network by the learning rule, the neural network method can handle many non-linear or other complex problems. In this work, neural network is applied to complicated inverse heat conduction problems. Efficiency of the procedure is enhanced by incorporating the radial basis functions (RBF). The RBF is trained faster than other neural network and can find smooth solution. In order to demonstrate the effectiveness of the current scheme, a typical one-dimensional IHCP is considered. At one surface, the temperature as well as the heat flux is known. The unknown temperature of interest is estimated on the other side of the slab. The results from the proposed method based on RBF neural network are compared with the conventional method.

In this paper, the effect of tip clearance on the cooling performance of the microchannel heat sink is presented under the fixed pumping power condition. For the various types of microchannel heat sink having different size of fin width and channel width, experimental study is conducted. Through the experiment, the tip clearance effect is investigated by increasing tip clearance from zero. As a result, it is shown that cooling performance of heat sink with tip clearance is better than that of heat sink without tip clearance. For the microchannel heat sink with tip clearance, the optimum conditions for cooling performance is also studied.

A fully coupled fluid flow, heat, and solute transport model was developed to investigate turbulent flow, solidification, and macrosegregation in a continuous casting process of steel slab with EMBR. Transport equations of mass, momentum, energy, and species for a binary iron-carbon alloy system were solved using a continuum model. The electromagnetic field was described by the Maxwell equations. A finite-volume method was employed to solve the conservation equations associated with appropriate boundary conditions. The effects of intensity of magnetic field and carbon segregation were investigated. The electromagnetic field reduces the velocity of molten flow in the mold and an increase in the percentage of C in steel results in a decrease of carbon segregation ratio.

This paper deals with the numerical analysis of natural convection flow induced by the density inversion effect of water inside horizontal pipe. The numerical method is based on SIMPLE/PWIM in general coordinate for its wide applicabilities. The numerical tool was validated through the comparison with the previous results concerning the density inversion effect of water It is shown that the developed numerical tool could predict the flow pattern and the heat transfer phenomena qualitatively And it is also found that the density inversion effect of water has significant effects on the flow pattern.

For the design of Aluminium pool furnace, position of burner and pool depth effects on flow and temperature field in Aluminium pool furnace are examined by the commercial computational code, CFD-ACE+. From the results, position of burner which is on the same face in side wall is better to distribute the flow field in Al furnace. That yields temperature to distribute more uniformly. And the burner position is on upper wall, fire frame reach pool surface. Customer must consider that, because it make Aluminium to oxidize.

Film growth rate of InP and GaAs using TMI, TMG, TBA and TBP is numerically predicted and compared to the experimental results. Obtained results show that the film growth rate is very sensitive to the thermal condition in the reactor. To obtain exact thermal boundary conditions at the reactor walls, we analyzed the gas flow and heat transfer in the reactor including outer tube as well as the inner reactor parts using a full three-dimensional model. The results indicate that the exact thermal boundary conditions are important to get precise film growth rate prediction.

In this study the optimization of plate-fin type heat sink with vortex generator for thermal stability is conducted numerically. To acquire the optimal design variables, the CFD and mathematical optimization are integrated. The flow and thermal fields are predicted using the finite volume method. The optimization is carried out by means of the sequential quadratic programming (SQP) method. The results show that when the temperature rise is less than 40 K, the optimal design variables are as follows; $B_1=2.584mm$, $B_2=1.741mm$, and t = 7.914 mm. Comparing with the initial design, the temperature rise is reduced by 4.2 K, while the pressure drop is increased by 9.43 Pa. The Pareto optimal solutions are also presented between the pressure drop and the temperature rise.

In this paper, the pressure drops were investigated according to the sintered porous wick structure in loop heat pipe(LHP) by theoretical analysis. LHP has the wick only in evaporator for the circulation of working fluid, so utilizes porous wick structure which pore diameter is very small for large capillary force. This paper investigates the effects of different parameters on the pressure drops of the LHP such as particle diameter of sintered porous wick, wick porosity, vapor line diameter, thickness of wick and heating capacity. Working fluid is water and the material of sintered porous wick is copper. According to the these different parameters, capillary pressure, pressure drop in wick were analized by theoretical design method of LHP.

In the present work, the practical thermal design process of 1 DIN car DVD receiver described. In the course of its efficient design, CAE technique was essentially used. CAE technique has reduced research period, man power and material cost but has increased research convenience, organized results and persuasive power. CAE technique helped to study parameters such as vent, fan and heat sink. Using these elements, it tried to meet optimal thermal solution. But safety standard, printed circuit board and framework mechanism should be considered as the constraint. To overcome these constraints, we tried to communicate and compromise with projectors in charge. After all, the price of those efforts has made the most competitive heat sink for heat dissipation in the 1 DIN car DVD receiver market. Moreover, we are trying to save $3 per product by removing fan. This paper is supposed to show an example of the CAE technique and help thermal designers to make electronic packaging goods.

In this study, single-phase heat transfer experiments were conducted with oil cooler with offset strip fin using water. An experimental water loop has been developed to measure the single-phase heat transfer coefficient in a vertical oil cooler. Downflow of hot water in one channel receives heal from the cold water upflow of water in the other channel. Similar to the case of a plate heat exchanger, even at a very low Reynolds number, the flow in the on cooler with offset strip fin remains turbulent. The present data show that the heat transfer coefficient increases with the Reynolds number. Based. On the present data, empirical correlation of the heat transfer coefficient was proposed. Also, performance prediction analysis for oil cooler were executed and compared with experiments. ${\varepsilon}-NTU$ method was used in this prediction program. Independent variables are flow rates and inlet temperature. Compared with experimental data, the accuracy of the program is within the error bounds of ${\pm}5$% in the heat transfer rate.

This study is focused on the comparison between the internal saturation temperature of the working fluid and the surface temperature of adiabatic zone of two-phase closed thermosyphons with various helical grooves. Distilled water, methanol and 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 results show that the numbers of grooves and the type of working fluids are very important factors for the operation of thermosyphons. A good agreement between the internal saturation temperature of working fluid and the surface temperature of adiabatic zone of two-phase closed thermosyphons with various helical grooves is obtained.

To develop more compact and light generators which have high capacity, the most important thing that should be considered is the inner cooling system. Under all circumstances, the temperature of rotor and stator windings must be kept below the maximum temperature of insulation to maintain reliability and prolong durability of the machine. Therefore, the development of more effective cooling system and the exact prediction of windings are essential to produce our unique generator model which is reliable and competitive in international market. In this study, the flow of cooling air and the temperature distribution of winding is analyzed by using computational fluid dynamics. This analysis can lead to optimize the structure of cooling system and predict a local temperature rise.

Decomposition of ozone at room temperature was investigated comparatively with commercial monolithic ozone decomposition catalyst (ODC, $MnO_2$) and monolithic photo catalyst ($TiO_2$). The effects of residence time, UV (ultraviolet) light dependence and ozone concentration on the conversion was presented. UV ray was irradiated using BLB (black light blue) lamp ($315{\sim}400$ nm), supplied with a constant intensity in the reactor. The concentration of ozone in the square-shape reactor can be controlled by combining the DBD (dielectric barrier discharge) reactor with an AC high voltage supply system. The catalytic performance, in presence of UV irradiation did not show significant changes for $MnO_2$ catalyst. $TiO_2$ catalyst was the different case, which showed higher decomposition activity in presence of UV irradiation. Deactivation of catalyst detected by real-time ozone monitor for 120 hours with a constant inlet ozone concentration.

Dielectric Barrier Discharge (DBD) in air, which has been established for the production of large quantities of ozone, is more recently being applied to a wider range of aftertreatment processes for HAPs (Hazardous Air Pollutants). Although DBD has high electron density and energy, its potential use as precharging nano and submicron particles are not well known. In this work, we measured I-V characteristics of DBD and estimated collection efficiency of the particles by DBD type 2-stage ESP. To examine the particle collection with various applied voltage waveforms of DBD for nano and submicron sized, bimodal particles were generated by a electrical tube furnace and an atomizer.

The temperature and soot particle measurement technique in a laminar diffusion flame have been studied to investigate the characteristics of soot particle with temperature using a co-flow burner. The temperature distributions in the flame were measured by rapid insertion of a R-type thermocouple and the soot particles were detected were detected by LEM/LIS techniques. In these measurement, soot volume fraction, number density and soot diameters were analyzed experimentally. As a results, the spacial distributions of particle volume fraction, soot diameter, and number density are mapped throughout the flame using the Rayleigh theory for the scattering of light by absorbing particles. A laser extinction method was used to measure the soot volume fraction and Laser induced scattering method was used to measure the soot particle diameter and number density. Also, we measured temperature without the effect of soot particles attached to the thermocouple junction, which is close to the nozzle. In this result, we found that upstream zone has a unstable flowing in co-flow diffusion flame and the y-axis temperature of flame has a uniform temperature distribution in the most soot volume fraction zone.

We performed the numerical study on the comparison of thermal comport performance indices for cooling loads in the classroom when the 4-way cassette air-conditioner is mounted on the ceiling. We investigated the velocity and the temperature distribution of the classroom as with respect to the variation of the air diffusion angle of the air-conditioner. Air diffusion performance index and Predicted mean vote were used for analyzing the characteristics of the thermal comport in the classroom and comparing their values each other. From the numerical results, we knew that the thermal comport is largely affected by the air diffusion angle and velocity of the air-conditioner. And we also found out that the qualitative tendency of the distribution between Air diffusion performance index and Predicted mean vote is very similar in all occupied zone.

Microsystem technology has been applied to space technology and became one of the enabling technology by which low cost and high efficiency are achievable. Micro propulsion system is a key technology in the miniature satellite because micro satellite requires very small and precise thrust force for maneuvering and attitude control. In this paper research on micro solid propellant thruster is reported. Micro solid propellant thruster has four basic components; micro combustion chamber, micro nozzle, solid propellant and micro igniter. In this research igniter, solid propellant and combustion chamber are focused. Micro igniter was fabricated through typical micromachining and evaluated. The characteristic of solid propellant was investigated to observe burning characteristic and to obtain burning velocity. Change of thrust force and the amount of energy loss following scale down at micro combustion chamber were estimated by numerical simulation based on empirical data and through the calculation normalized specific impulses were compared to figure out the efficiency of combustion chamber.

In order for studying pressure-coupled dynamic responses of droplet vaporization, open-loop experiment of an isolated droplet vaporization exposed to pressure perturbations in stagnant gaseous environment is numerically conducted. Governing equations are solved for flow parameters at gas and liquid phases separately and thermodynamic parameters at the interfacial boundary are matched for problem closure. For high-pressure effects, vapor-liquid interfacial thermodynamics is rigorously treated. A series of parametric calculations in terms of mean pressure level and wave frequencies are carried out employing a n-pentane droplet in stagnant gaseous nitrogen. Results show that wave instability in view of pressure-coupled vaporization response seems more susceptible at higher pressures and higher wave frequencies. Mass evaporation rate responding to pressure waves is amplified with increase in pressure due to substantial reduction in latent heat of vaporization. Augmentation of perturbation frequency also enhances amplification due to the reduction of phase differences between pressure perturbation and surface temperature fluctuation.

Acoustic characteristics in an industrial gas-turbine combustor are numerically investigated by adopting linear acoustic analysis. Spatially non-homogeneous temperature field in the combustor is considered in the numerical calculation and the characteristics are analyzed in view of acoustic instability. Acoustic analysis are conducted in the combustors without and with acoustic resonator, which is one of combustion stabilization devices. It has been reported that severe pressure fluctuation frequently occurs in the adopted combustor, and the measured signal of pressure oscillation is compared with the acoustic-pressure response from the numerical calculation. The numerical results are in a good agreement with the measurement data. In this regard, the phenomenon of pressure fluctuation in the combustor could be caused by acoustic instability. The acoustic effects of the resonators are analyzed in the viewpoints of both the frequency tuning and the damping capacity.

The characteristics of soot near the soot inception point for an ethene-air flame was carried out in a WSR (well-stirred reactor). The new sampling tool like the temperature controlled filter system was introduced to minimize the condensation during sampling. The new analysis tools applied include the real time size distribution analysis with the Nano-DMA, particle size by transmission electron microscopy, C/H analysis, g filter analysis, and thermogravimetric analysis using both non-oxidative and oxidative pyrolysis. The WSR can generate young soot particles that can be collected and examined to gain insight into inception. For the current conditions, soot does not form for ${\Phi}=1.9$, inception occurs at or before ${\Phi}=2.0$, and inception combined with soot surface growth and/or coagulation occurs for ${\Phi=2.1}$. The filter samples for ${\Phi}$=1.9 are composed of volatile compounds that evolve at relatively low temperatures when heated in the presence or absence of $O_2$. The samples collected from the WSR at ${\Phi}=2.0$ and ${\Phi}=2.1$ are precursor-like in morphology and size. They have higher C/H ratios and lower organic percentages than precursor particles, but they are clearly not fully carbonized soot. The WSR PAH distribution is similar to that in young soot from inverse flames.

Synthesis of carbon nanotubes and nanofibers on a catalytic metal substrate, using an ethylene fueled inverse diffusion flame, was investigated. Multi-walled carbon nanotubes, with diameters of 20 - 60nm, were formed on the substrate coated with nickel-nitrate in the region of 5 - 6mm from the flame center along the radial direction. The gas temperature for this region was ranging from about 1400 to 900K. Nickel particles originated from the coated nickel-nitrate on the substrate were the major catalyst for the formation of the nanomaterials. HR-TEM and Raman spectrum revealed that synthesized carbon nanotubes had multi-walled structures with some defective graphite layers at walls.

The present study deals with the unique characteristics of hydrogen jet diffusion flames, such as split flames and reignition phenomenon. The split flames are composed of a small flamelet on the nozzle rim and a lifted main flame at downstream. When mass flow rates of fuel reach a critical point, a small-sized flamelet is found to remain in the vicinity of the nozzle exit and the flame reignition subsequent to blowout of main flame occurs repeatedly. In this study, the non-luminous hydrogen jet diffusion flames are visualized by using schlieren technique in order to analyze the combustion characteristics of hydrogen jet diffusion flames with focus on the flame reignition phenomenon.

An experimental study is performed on atomization characteristics and stable operating condition for injection of high viscous waste vegetable oil using effervescent atomizer with two aerator tubes. Consideration is given to the effects of ALR and liquid viscosity on the velocity and mean diameter of the injected droplet. It is found that (i) as ALR increases, the axial velocity of the droplet is increased, while half-velocity width and SMD are decreased regardless of the change in liquid viscosities, (ii) the rate of fine drop distribution occupied in the total spray field is increased with an increasing in ALR, and (iii) the effect of viscosity on atomization characteristics is minor.

In the present work, a series of experiments have been performed on electro-hydrodynamic atomization of non-conducting liquid using a charge injection type nozzle. Effects of liquid flow rate, input voltage, and distance between the needle and the ground electrode (nozzle-embedded metal plate) have been examined. For fixed electrode distances, total and spray currents increase with increase of liquid flow rate and input voltage. When the distance between the needle and the ground electrode becomes closer, total, leakage and spray current increase, but the onset voltage for dielectric breakdown decreases. When the electric field strength of the liquid jet exceeds that for the air breakdown, a portion of the electric charges in the liquid jet is dissipated into the ambient air, and the charge density shows a limiting value. Atomization quality can be improved by increasing the flow rate because the higher charge density is achieved with the larger liquid velocity in addition to the enhanced aerodynamic effect.

The Flame length of $CH_4$ with the Oxidizer of air and $O_2$ has been measured respectively for the nozzle diameter of 1.6mm, 2.7mm, 4.4mm and 7.7mm. In all $CH_4$ flame on oxidizer of air and $O_2$. the flame length was independent of the initial jet diameter, dependent only on the flowrate in laminar flame regime, and in turbulent flame dependent on the initial jet diameter. Using correlation equation of Delichatsios, the flame length has been expected exactly for $CH_4/air$ flame, but has been underestimated for $CH_4/O_2$ flame. This paper has proposed correlation equation of $CH_4/O_2$ flame.

Nowadays CFCs and HCFCs refrigerants are restricted because it cause depletion of ozone layer. Accordingly, this experiment apply the $NH_3$ gas and not CFCs and HCFCs for refrigerant to study the performance characteristic from the superheat control and improve the energy efficiency from the high performance. The condensing pressure of refrigeration system is increased from 14.5bar to 16bar by 0.5bar and superheat temperature is increased from $0^{\circ}C$ to $10^{\circ}C$ by $1^{\circ}C$ at each condensing pressure. As the result of experiment, when the superheat temperature is $1^{\circ}C$ at each condensing pressure, the refrigeration system has the high performance.

Convective boiling heat transfer coefficients of R-22 were obtained in a flat extruded aluminum tube with $D_h=1.41mm$ . The test range covered mass flux from 200 to 600 $kg/m^2s$, heat flux from 5 to 15 $kW/m^2$ and saturation temperature from $5^{\circ}C$ to $15^{\circ}C$ . The heat transfer coefficient curve shows a decreasing trend after a certain quality(critical quality). The critical quality decreases as the heat flux increases, and as the mass flux decreases. The early dryout at a high heat flux results in a unique 'cross-over' of the heat transfer coefficient curves. The heat transfer coefficient increases as the mass flux increases. At a low quality region, however, the effect of mass flux is not prominent. The heat transfer coefficient increases as the saturation temperature increases. The effect of saturation temperature, however, diminishes as the heat flux decreases. Both the Shah and the Kandlikar correlations underpredict the low mass flux and overpredict the high mass flux data.

This study is concerned with the thermal analysis during the cool-down period of 135,000 $m^3$ class GT-96 membrane type LNG carrier under IMO design condition. The cool-down is performed to cool the insulation wall and the natural gas in cargo tank for six hours to avoid the thermal shock at the start of loading of $-163^{\circ}C$ LNG. During the cool-down period, the spraying rate for the NG cooling decreases as the temperature of NG falls clown from $-40^{\circ}C$ to $-130^{\circ}C$ and the spraying rate for the insulation wall cooling increases as the temperature gradient of the insulation wall is large. It was confirmed that there existed the largest temperature decrease at the 1 st barrier and 1st insulation, which are among the insulation wall, especially in the top side of the insulation wall. By the 3-D numerical calculation about the cargo tank and the cofferdam during the cool-down period, the temperature variation in hulls and insulations is precisely predicted.

Thermal error compensation has been developed for CNC (Computer Numerical Control) machining center with moving heat sources. The thermal error in CNC machining center has an effect on machining accuracy more than the geometric error does. Thus, temperature distributions of a spindle unit have been analyzed numerically by a Finite Differential Method and experimentally by an infrared (IR) camera in this study. A multiple variable method has been derived to estimate the thermal deformation of the machine origin stably and effectively after measuring deformation and temperature data. The experimental results for a vertical machining center have shown that the thermal errors of the machine origins were reduced more than 30% by the developed method.

The objective of this work is to propose calibration facility in which a thin film type heat flux sensor can be calibrated under convective flow condition by using a small wind tunnel with the constant temperature plate condition. A small wind tunnel has been built to produce a boundary layer shear flow above a constant temperature copper plate. 12-independent copper blocks, thin film heaters, insulators and temperature controllers were used to keep the temperature of flat plate constant at a specified temperature. Three commercial thin film-type heat flux sensors were tested. Convective calibrations of these gages were performed over the available heat flux range of $1.4{\sim}2.5kW/m^2$. The uncertainty in the heat flux measurements in the convective-type heat flux calibration facility was ${\pm}2.07%$. Non-dimensional sensitivity is proposed to compare the sensitivity calibrated by manufacturer and that of experiment conducted in this study.

New method for the design, fabrication, and calibration of micro-machined heat flux sensor has been developed. Two types of micro-machined heat flux sensor having different thicknesses of the thermal-resistance layer are fabricated using the MEMS technique. Photo-resist patterning using a chrome mask, bulk-etching and copper-nickel sputtering using a shadow mask are applied to make heat flux sensors, which are calibrated in the convection-type heat flux calibration facility. The sensitivity of the device varies with thermal-resistance layer, and hence can be used to measure the heat flux in heat-transfer phenomena.

This article investigates the heat transfer characteristics in a RTA system for LCD manufacturing and suggests a way to evaluate the quality of a poly-Si film from the thin film optics analysis. The transient and one-dimensional conductive/radiative heat transfer equation considering wave interference effect is solved to predict surface temperatures of thin films. In dealing with radiative heat transfer, a one-dimensional two-flux method is used and the ray tracing method is also utilized to account for the wave interference effects. It is assumed that each interface is assumed diffusive but the spectral radiative properties are included. It is found that the selective heating region exists for various wavelengths and consequently may contribute to heat the poly-Si film. Using the formalism of the characteristic transmission matrix, the lumped structure reflectance, transmittance, and absorptance are calculated and they are compared with experimental data of the poly-Si film during the SPC process via the FE-RTA (Field-Enhanced RTA) technology.

Understanding the thermal conductivity and heat transfer processes in superlattice structures is critical for the development of thermoelectric materials and optoelectronic devices based on quantum structures. $Chen^{(1)}$ developed ballistic diffusive equation(BDE) for alternatives of the Boltzmann equation that can be applied to the complex geometrical situation. In this study, a simulation code based on BDE is developed and applied to the 1-dimensional transient heat conduction across a thin film and transient 2-dimensional heat conduction across the film with heater. The obtained results are compared to the results of the $Chen^{(1)}$ and Yang and $Chen^{(1)}$. Finally, steady 2-dimensional heat conduction in the quantum dot superlattice are solved to obtain the equivalent thermal conductivity of the lattice and also compared with the experimental data from $Borca-Tasciuc^{(2)}$.

From the viewpoint of a macro state, there is no thermal boundary resistance (TBR) at an interface if both surfaces at an interface are perfectly contacted. However, recent molecular dynamics (MD) studies reported that there still exists the TDR at the interface in an ideal epitaxial superlttice. Our previous studies suggested the model to predict the TBR not only quantitatively also qualitatively in superlattices. The suggested model was based on the classical theory of a wave reflection, and provided highly satisfactory results for an engineering purpose. However, it was not the complete model because our previous model was derived by considering only the effects from a mass ratio and a potential ratio of two species. The interaction of two species presented by the Lennard-Jones (L-J) potential is governed by the mutual ratio of the masses, the potential well depths, and the diameters. In this study, we performed the preliminary simulations to investigate the effect resulting from the diameter ratio of two species for the completion of our model and confirmed that it was also a ruling factor to the TBR at an interface in superlattices.

Investigators have been perplexed with the thermal phenomena behind the recently discovered nanofluids, fluids with unprecedented stability of suspended nanoparticles although huge difference in the density of nanoparticles and fluid. For example, nanofluids have anomalously high thermal conductivities at very low fraction, strongly temperature-dependent and size-dependent conductivities, and three-fold higher critical heat flux than that of base fluids. Traditional conductivity theories such as the Maxwell or other macroscale approaches cannot explain why nanofluids have these intriguing features. So in this paper, we devise a theoretical model that accounts for the fundamental role of dynamic nanoparticles in nanofluids. The proposed model not only captures the concentration and temperature-dependent conductivity, but also predicts strongly size-dependent conductivity. Furthermore, we physically explain the new phenomena for nanofluids. In addition, based on a proposed model, the effects of various parameters such as the ratio of thermal conductivity of nanofluids to that of a base fluid, volume fraction, nanoparticle size, and temperature on the thermal conductivities of nanofluids are investigated.

Recently solar energy is representative in the technology development and spread of alternative energy. Specially in condition of solar collectors, they have had very various shape. This paper reports experimental study about the application of Oscillating Capillary Tube Heat Pipe to flat-plate solar collector. In conclusion, overall temperature distribution of OCHP was investigated by charging ratio and inclination angle. Respective charging ratio is 15%, 20%, 40% and respective inclination angle is horizontal, $30^{\circ}$, $45^{\circ}$, $60^{\circ}$, perpendicular. As a result of experiment, charging ratio 20% heat pipe has shown the most uniform temperature distribution and also performance of heat transfer has been the best.

Experimental study is performed to investigate the effect of heat load and operating temperature on the thermal performance of a heat pipe with screen mesh wick. The heat pipe was designed in 200 screen meshes, 500mm length and 12.7mm O.D tube of copper, water as working fluid(4.8g) and nitrogen as non-condensible gas(NCG). The heat pipe used in this study has evaporator, condenser and adiabatic section, respectively. Experimental data of axial wall temperature distribution is presented for heat transport capacity, the temperature of cooling water of condenser, inclination angle, and operating temperature. For the results from this study, it is found that, for the same charging mass of working fluid, the initial operating temperature and the overall wall temperatures of heat pipe are higher for NCG charging mass of $5.0{\times}10^{-6}kg$ and $3.4{\times}10^{-6}kg$, than that of $1.0{\times}10^{-6}kg$.

A series of experiments was conducted to examine characteristics of a grooved flat-strip heat pipe having multiple heat sources. The inner grooves of the heat pipe have the aspect ratio of 1 to $2.5(0.42{\times}1.05$ mm) whose pitch was 0.6 mm. Four block heaters ($10{\times}20$ mm) were placed in the evaporator section at intervals of 20 mm and six different heating modes were tested. The maximum surface heat flux of 80 $W/cm^2$ was achieved while the operating temperature was kept below $100^{\circ}C$, In the nearest heating mode (from the condenser location), the heat pipe exhibited more stable temperature distribution than the far heating mode where the heaters is located furthest from the condenser.

Concentric annular heat pipes having the length of 200 mm and the outside diameter of 25.4 mm were manufactured and tested. The inside diameters of the heat pipes were 11.3 mm and 8.1 mm and the material of the container was copper. To compare with thermal performances between the concentric annular heat pipe and a copper block with same shape, start-up and isothermal performance were tested. In the result, the start-up of the concentric annular heat pipe was influenced by the start-up of their heat source. The concentric annular heat pipe reached at steady state faster then the same shape of copper block. The maximum wall temperature difference of the concentric annular heat pipe on whole surface was $4.6^{\circ}C$, and the case of the copper block was $16.5^{\circ}C$.

Optimization of a heat exchanger with internally finned circular tubes has been performed for three-dimensional periodically fully developed turbulent flow and heat transfer. The design variables of fin number N, fin width ($d_1,d_2$) and fin height(H), are numerically optimized for the limiting conditions of $N=22{\sim}37$, $d_1=0.5{\sim}1.5$ mm, $d_2=0.5{\sim}1.5$ mm, $H=0.1{\sim}1.5$. Due to the periodic boundary conditions along main flow direction, the three layers of meshes are considered. The CFD and the mathematical optimization are coupled to optimize the heat exchanger. 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.

Distillation is one of many processes that can be used for the purification of water. It requires an energy input such as heat and solar radiation is a possible source of energy. In this process, water is evaporated, thus separating water vapor from the solute. The vapor is then condensed to pure water. The temperature of the water, the cavity of the enclosure and the surface of the glass was measured everyday. Fifteen points were chosen for the temperature measure using a thermocouple. The inner wall and the bottom of each still was painted black for good absorption of heat. The enforced glass was used for the cover for the entering of solar energy. The size of all of the water baths was the same, but the glass of the rectangular form had a tilted angle. In the case of fine or general weather, the volume of condensed water produced by $45^{\circ}$ and $60^{\circ}$ is very irregular compared to that of $15^{\circ}$ and $30^{\circ}$, In case of a rainy day, the occurrence of the volume of condensed water was similar to that of $15^{\circ}$ and $30^{\circ}$, but the volume of condensed water produced by $30^{\circ}$ was the highest.

The present paper is devoted to the modeling based on an averaging approach for microchannel heat sinks. Firstly, analytic solutions for velocity and temperature distributions for low-aspect-ratio microchannel heat sinks are presented by using the averaging approach. When the aspect ratio of the microchannel is smaller than 1, analytic solutions accurately evaluate thermal resistances of heat sinks while the previous model cannot predict thermal resistances. Secondly, asymptotic solutions for velocity and temperature distributions at low-aspect-ratio limit and at high-aspect-ratio limit are presented by using the scale analysis. Asymptotic solutions are very simple, but shown to predict thermal resistances accurately.

The internal cooling passage of a gas turbine blade can be modeled as a ribbed channel. Most studies have considered square ribs. However, the ribs can be rounded due to improper manufacturing or wear during the operation. Hence, we have studied two different rib geometries in this study, i.e. square and semicircle ribs. We have performed large eddy simulations (LES) and experiments to validate the results from the simulations. LES predicts the detailed flow and thermal features, which have not been captured by simulations using turbulence models. By investigating the instantaneous flow and thermal fields, we propose the mechanisms for the local heat transfer distribution between ribs. For both the geometries, heat transfer is enhanced by the entrainment of the cold fluid by the vortical motions and impingement of the entrained cold fluid on the ribs.

The heat (mass) transfer characteristics on the blade surface of a first-stage turbine rotor cascade has been investigated by employing the naphthalene sublimation technique. A four-axis profile measurement system is employed for the measurements of the local heat (mass) transfer coefficient on the curved blade surface. The experiments are carried out for two free-stream turbulence intensities of 1.2% and 14.7%. The high free-stream turbulence results in more uniform distributions of heat load on the both pressure and suction surfaces and in an early boundary-layer separation on the suction surface. The heat (mass) transfer enhancement on the suction surface due to the endwall vortices is found to be relatively small under the high free-stream turbulence.

The laminar impinging jet thermal fields were investigated with or without magnetic fields. The transient phenomenon from steady to unsteady flow was founded at specific Reynolds number ranges. In unsteady flow region, the magnetic fields make flow stable. So the characteristics of heat transfer at impingement wall are changed

The water jet impingement cooling is one of the techniques to remove heat from high heat flux equipments. We investigate the local heat transfer of the confined water impinging jet and the effect of nozzle collar to enhance the heat transfer in the free surface jet and submerged jet. Boiling is initiated from the furthest downstream and the wall temperature increase is reduced with developing boiling, forming the flat temperature distributions. The reduction in the nozzle-to-surface distance for $H/W{\leq}1$ causes the significant increases and distribution changes in heat transfer. Developed boiling reduces the differences in heat transfer for various conditions. The nozzle collar is employed at the nozzle exit. The distances from heated surface to guide plate, $H_c$ are 0.25W, 0.5W and 1.0W. The liquid film thickness is reduced and the velocity of wall jet increase as decreased spacing of collar to heated surface. Heat transfer is enhanced for region from the stagnation to $x/W{\sim}8$ in the free surface jet and to $x/W{\sim}5$ in the submerged jet. For nucleate boiling region of further downstream, the heat transfer by the nozzle collar is decreased in submerged jet compare with higher velocity condition. It is because the increased velocity by collar is de-accelerated at downstream.

By using unique experimental techniques and careful construction of the experimental apparatus, the characteristics of the local heat transfer were investigated using the condensing R134a two-phase flow, in horizontal single mini-channels. The circular channels ($D_h=0.493$, 0.691, and 1.067 mm) and rectangular channels ($D_h=0.494$, 0.658, and 0.972 mm) were tested and compared. Tests were performed for a mass flux of 100, 200, 400, and 600 $kg/m^2s$, a heat flux of 5 to 20 $kW/m^2$, and a saturation temperature of $40^{\circ}C$. In this study, effect of heat flux, mass flux, vapor qualities, hydraulic diameter, and channel geometry on flow condensation were investigated and the experimental local condensation heat transfer coefficients are shown. The experimental data of condensation Nusselt number are compared with existing correlations.

In this study, micro square heaters having dimensions of $65{\times}65{\mu}m^2$and $100{\times}100{\mu}m^2$ were fabricated and bubble nucleation experiments on the heaters were performed. Bubble nucleation temperature was also measured using a bridge circuit and the photographs of bubble nucleation and subsequent growth were taken by a camera with a flash unit. Measured bubble nucleation temperatures were found to be closer to the superheat limit of working fluid (FC-72). Also quasi-1D analyses for the square heaters were performed. The quasi-1D analysis yielded proper temperature distribution of the square heater at steady state, however failed to predict the temperature rise up to the steady state. Similar time dependent temperature can be obtained with proper value of thermal diffusivity. For the $100{\times}100{\mu}m^2$ square heater, nucleation of several bubbles was observed while only one bubble was observed to be nucleated on $65{\times}65{\mu}m^2$ heater.

'Nanofluids' means suspension of common fluids with particles of the order of nanometers in size. The present research is an experimental study of critical heat flux (CHF) behavior in pool boiling of $water-TiO_2$ nanofluids under atmospheric pressure. CHF for pure water and $water-TiO_2$ nanofluids were respectively measured using disk-type copper block heater with 10mm diameter, and CHF of water with surfactant was also measured to consider the effect of surfactant used to disperse nanoparticle. The results show a large increase in CHF for $water-TiO_2$ nanofluids compared to pure water. After CHF occurred, heat flux in pool boiling for $water-TiO_2$ nanofluids was maintained in considerable value, but not for pure water.

Plasma Enhanced Chemical Vapor Deposition (PECVD) process uses unique property of plasma to modify surfaces and to achieve the high deposition rates. In this study, a vertical thermal RF-PECVD (Radio Frequency-PECVD) reactor is modeled to investigate thermal flow and the deposition rates with various shapes of the showerhead. The showerhead in the CVD reactor has the shape of a ring and gases are injected in parallel with the susceptor, which is a rotating disk. In order to achieve the high deposition rates, we have simulated the thermal flow fields in the reactor with several showerhead models. Especially the effects of the number of injection holes and the rotating speed of the susceptor are studied. Using a commercial code, CFDACE, which uses FVM (Finite Volume Method) and SIMPLE algorithm, governing equations have been solved for the pressure, mass-flow rates and temperature distributions in the CVD reactor. With the help of the Nusselt number and Sherwood number, the heat and mass transfers on the susceptor are investigated. In order to characteristics of measure the flatness of the layer, furthermore, the relative growth rate (RGR) is considered.

This report presents experimental results on the heat transfer coefficients in the boiling region of spray cooling for actual metallurgical process. In this study, the heat flux distribution of a two dimensional dilute spray impinging on hot plate was experimentally investigated. Based on the experimental results, they classified the heat transfer area into the stagnation region and wall-flow region. In the stagnation region, the local heat transfer coefficient relates mainly to the droplet-flow-rate supplied from spray nozzle directly, so the local heat transfer coefficients is good agreement with the predicted values from correlation for spray cooling proposed by former report However, the local heat transfer coefficient in wall-flow region is larger than predicted values, and it is found that the rebounding droplets-flow-rate must be accurately evaluated to predict the local heat transfer coefficient in this region.

For the present study, three heat pipes with different thickness of sintered metal wick were manufactured, and their operational performances, such as capillary limit and thermal resistance were tested and compared with theoretical predictions. Experimental results showed good agreement with those by the theoretical model, and that seemed to present that the sintering process we had developed in the present study was valid.

WSGGM with gray gas regrouping is successfully applied to study the flame structure of counter flow flames including effect of radiative transfer. The statistical narrow band model is used to obtain the benchmark solutions. Results obtained by using the optically thin model are shown to overestimate the emission and to predict the flame structures inadequately especially for optically thick and low stretch rate flames. Computed results by using the WSGGM with 10 gray gases and SNB model show reasonable agreements with each other, and the required calculation time for the WSGGM is acceptable for engineering applications.

Flue gas recirculation (FGR) is a method used to control oxides of nitrogen ($NO_x$) in combustion system. The recirculated flue gases resulted in slow reaction and low flame temperatures, which in turn resulted in decreased thermal NO production. Recently, it has been demonstrated that introducing the recirculated flue gas in the fuel stream, that is, the fuel induced recirculation (FIR), resulted in a much greater reduction in $NO_x$ per unit mass of recirculated gas, as compared to introducing the flue gases in air. In the present study, the effect on $NO_x$ reduction in turbulent swirl flame in laboratory scale using FGR/FIR methods through the dilution using $N_2$ and $CO_2$. Results. show the $CO_2$ dilution is more effective $NO_x$ reduction methods because of large temperature drop due to the larger specific heat $CO_2$ compared to $N_2$. FIR is more effective to reduce $NO_x$ emission than FGR when the same recirculation ratio of dilution gas.

A novel technique for fine particle beam focusing under the atmospheric pressure is introduced using a radiation pressure assisted aerodynamic lens. To introduce the radiation pressure in the aerodynamic focusing system, a 25 mm plano-convex lens having 2.5 mm hole at its center is used as an orifice. The particle beam width is measured for various laser power, particle size, and flow velocity. In addition, the effect of the laser characteristics on the beam focusing is evaluated comparing an Ar-Ion continuous wave laser and a pulsed Nd-YAG laser. For the pure aerodynamic focusing system, the particle beam width was decreased as increasing particle size and Reynolds number. For the particle diameter of 0.5 ${\mu}m$, the particle beam was broken due to the secondary flow at Reynolds number of 694. Using the Ar-Ion CW laser, the particle beam width becomes smaller than that of the pure aerodynamic focusing system about 16 %, 11.4 % and 9.6 % for PSL particle size of 2.5 ${\mu}m$, 1.0 ${\mu}m$, and 0.5 ${\mu}m$ respectively at the Reynolds number of 320. Particle beam width was minimized around the laser power of 0.2 W. However, as increasing the laser power higher than 0.4 W, the particle beam width was increased a little and it approached almost a constant value which is still smaller than that of the pure aerodynamic focusing system. The radiation pressure effect on the particle beam width is intensified as Reynolds number decreases or particle size increases relatively. On the other hand, using 30 Hz pulsed Nd-YAG laser, the effect of the radiation pressure on the particle beam width was not distinct unlike Ar-Ion CW laser.

The effects of additive ions on the generation of metallic nanoparticles were evaluated using a corona induced supersonic nozzle. Applying the corona discharge to the nanoparticle generator, a tungsten needle and the supersonic nozzle are used as an anode electrode and a cathode electrode respectively. The corona ions act as nuclei for the silver vapor condensation. The ion density was controlled precisely as varying the applied voltage between electrode and nozzle. The mean diameter of the silver particle decreases as the ion density increases. However, the number concentration of the silver particle tended to increase with the ion density. The size distribution is more uniform as the ion density increases.

An inverse boundary analysis of surface radiation in an axisymmetric cylindrical enclosure has been conducted in this study. Net energy exchange method was used to calculate the radiative heat flux on each surface, and a hybrid genetic algorithm was adopted to minimize an objective function, which is expressed by sum of square errors between estimated and measured heat fluxes on the design surface. We have examined the effects of the measurement error as well as the number of measurement points on the estimation accuracy.

The purpose of this study is to investigate the effect of swirler angle and swirl chamber aspect ration of nozzle on the characteristics of single and twin spray. The performances of nozzle has been investigated by measurements of spray angle, droplet size, velocity and Weber number at a water pressure 0.4MHz. Visualization of spray was conducted to obtain the spray angle and breakup process. The spray characteristics such as droplet size and velocity were measured by Phase Doppler Anemometry(PDA). It was found that the smaller swirler angle, the larger axial velocity became. It was also shown that the larger aspect ratio, the smaller droplet diameter became.

In this study, we have manufactured the LHP(Loop Heat Pipe) with sintered metal wick and investigated the working characteristics of LHP experimentally.Water was used as a working fluid and fill charge rate was changed. LHP basically consist of the separated vapor/liquid channels, evaporator having sintered metal wick(effective pore diameter :$16{\sim}19{\mu}m$), and condenser cooled by water. The diameter of vapor/liquid line tube are 3.2mm/6.35mm, respectively. Heat transfer rate and thermal resistance was represented to study the basic characteristics of LHP at each conditions

The objective of this study is to simulate the etching characteristics with oscillation angle for the optimization of etching system. The etching characteristics were analyzed under different etching conditions. The spray characteristics were measured by Phase Doppler Anemometer (PDA). The correlation between the spray characteristics and the etching characteristics was investigated and used for fundamental data to simulate the etching characteristics with oscillation angle. The smaller coefficient of variation, the more uniform etching characteristic distribution became. It was found that numerical predictions of etching factor generally agreed well with the measured results with distance from nozzle tip. Oscillation leads to decrease of etching factor and increase of uniformity.

This paper presents a semi-empirical model to predict the frost growth formed on the cold cylinder surface. The model is composed of the correlations for frost properties including the various frosting parameters and local heat transfer coefficient. The effects of varying the correlations for local heat transfer coefficient on the frost growth are examined to establish the model. The numerical results are compared with experimental data obtained by the previous researchers. The results agree well with the experimental data within a maximum error of 13%. As the results, the frost thickness decreases with changing angular position from front stagnation to separation point. Also the effects of air velocity on the frost growth are negligible, as compared to the other frosting parameters.

This paper presents a mathematical model to predict the frost properties and heal and mass transfer within the frost layer formed on a cold plate. The laminar flow equations for the air-side are analyzed. and the empirical correlations of local frost properties are employed in order to predict the frost layer growth. The correlations of local frost density and effective thermal conductivity of frost layer, obtained from various experimental conditions, are derived as functions of various frosting parameters (Reynolds number, frost surface temperature, absolute humidity and temperature of moist air, cooling plate temperature, and frost density). The numerical results are compared with experimental data and the results of various models to validate the present model, and agree well with experimental data within a maximum error of 10%. The heat and mass transfer coefficients obtained from the numerical analyses are presented, as the results, it is found that the model for frost growth using the correlation of heat transfer coefficient without solving air flow have a limitation in its application.

This paper introduced about characteristics on the non-steady heat transfer of STS 304 hollow cylinder, In the non-steady state, the specific heat and conductivity are depended on the temperature variations, and these properties affect to the governing equation on heat conduction. But the most of numerical analysis on heat conduction is assumed to constant properties which is conductivity and specific heat. Assuming that conduction is assumed to constant properties which is conductivity and specific heat. Assuming that the properties are reacted sensitively, the numerical results can have the difference of between constant properties with non-constant properties. The main parameters are specific heat and conductivity. The temperature distributions of the STS 304 hollow cylinder became in steady state after 4 minutes in case of the constant properties. As the conductivity in varied with temperature, the temperature distributions became in steady state after 15 minutes. Therefore, a numerical analysis of the non steady state heat transfer will has to apply that conductivity varied with temperature.

This work applies the nonequilibrium molecular dynamics simulation method to study a Lennard-Jones liquid thin film suspended in the vapor and calculates the thermal conductivity by linear response function. As a preliminary test, the thermal conductivity of pure argon fluid are calculated by nonequilibrium molecular dynamics simulation. It is found that the thermal conductivity decrease with decreasing the density. When both argon liquid and vapor phase are present, the effects of the system temperature on the thermal conductivity are investigated. It can be seen that the thermal conductivity of liquid-vapor interface is constant with increasing the temperature

Insulating a heat from the environment is the most important in cryogenic applications like HTS cable system. Vacuum and MLI(Multi-Layer Insulation) have been widely used to get highly efficient cryogenic insulation. In this study, experiment on annular cylinder regarded as basic model of HTS cable cryostat has been performed to measure the heat loads. To investigate the effectiveness of radiation shields in cryogenic insulation system, the experiments are carried out to various number of the shields. The measured values are compared with the results estimated theoretically. The heat loads invaded from the environment was determined by liquid nitrogen boil-off calorimetry.

In order to evaluate the performance of fire-fighting agents used to protect structures from heat and fire damages, the thermal characteristics of fire-protection foams are experimentally investigated. The current research focuses on the destruction of a fire-fighting foam subjected to heat radiation. A simple repeatable test for fire-protection foams subjected to fire radiation is developed. This test involves foam generation equipment, a fire source for heat generation, repeatable test procedures, and data acquisition techniques. Results of the experimental procedure indicated that each thermocouple within the foam responded in a similar manner and gradually to a temperature of $15^{\circ}C{\sim}20^{\circ}C$. At this point, each trace generally rises to a temperature of approximately $90^{\circ}C$. The temperature gradient in the foam as time passes increases with increased foam expansion ratio. In addition, it is determined that the temperature gradient along the foam for depth decreases with increased foam expansion ratio.

An experimental study has been carried out to investigate the frosting behaviors of thermally conductive plastic(PBT based resin) resin by comparing with those of aluminum and some plastic(PTFE based resin) test specimens. It is found that the frosting behavior of plastic specimens with 1 mm thickness show similar trend with aluminum except PTFE. The properties of frost formed on the specimens are affected by both thermal conductivityand surface characteristics (hydrophilic/hydrophobic) of the materials. It can be said that the heat and mass transfer rate of plastic materials are almost equivalent with those of aluminum.

The present experimental study aimed to investigate the feasibility of physical water treatment (PWT) devices using catalytic materials to mitigation of crystalization fouling. Two PWT devices having different shape and material were used. The results from microscopic observation and SEM photographs demonstrate that the crystal structure of $CaCO_3$ is Argonite type without water treatment while Calcite type with water treatment, which clearly shows the feasibility of PWT techinque using catalytic materials to mitigate crystalization fouling. For 500 ppm solution, the fouling resistance decreased up to about 23% due to physical water treatment using catalytic materials. The results also reveal that if two more techniques are used simultaneously the additional effects could be expected.

A numerical analysis has been performed to give an understanding of the physics of a compressible base flow with mass bleed in a Mach 2.47 freestream. Axisymmetric, compressible mass-averaged Navier-Stokes equations are computed using a two-equation turbulence model, standard ${\kappa}-{\omega}$, and a fully implicit finite volume scheme. The mass bleed is characterized by the change in the mass flow rate of the bleed jet non-dimensionalized by the product of the base area and freestream mass flux. The result showing that there is an optimum bleed condition with maximum base pressure, leading to a minimum base drag, is clearly predicted and the validation with experimental data shows reasonable agreement.

When a shock wave arrives at a duct, an impulsive wave is discharged from the duct exit and causes serious noise and vibration problems. In the current study, the characteristics of the impulsive wave discharged from a partial closed duct exit is numerically investigated using a CFD method. The Yee-Roe- Davis's total variation diminishing(TVD) scheme is used to solve the axisymmetric, unsteady, compressible Euler equations. With several partial closed duct exits, the Mach number of the incident shock wave $M_s$ and the distance L/D between the duct exit and a flat plate are varied in the range of $M_s$ = 1.01 ${\sim}$ 1.50 and L/D = 1.0 ${\sim}$ 4.0, respectively. The results obtained show that the magnitude of the impulsive wave impinging upon the flat plate strongly depends upon $M_s$, L/D and the partial closure of duct exit. The impulsive wave on the flat plate can be considerably alleviated by the partial closure of duct exit and, thus, the present method can be a passive control for the impulsive wave.

In this paper, the unsteady behavior of the viscous flow field past an impulsively started elliptic cylinder is studied numerically. In order to analyze flow field, we introduce vortex particle method. The vorticity transport equation is solved by fractional step algorithm which splits into convection term and diffusion term. The convection term is calculated with Biot-Savart law, the no-through boundary condition is employed on solid boundaries. The diffusion term is modified based on the scheme of particle strength exchange. The particle redistributed scheme for general geometry is adapted. The flows around an elliptic cylinder are investigated for various attack angles at Reynolds number 200. The comparison between numerical results of present study and experimental data shows good agreements.

Numerical simulations are performed to investigate the characteristics of flow past a sphere in uniform shear. The Reynolds numbers considered are Re=300, 425 and 480 based on the inlet center velocity and sphere diameter. The non-dimensional shear rate K of the inlet uniform shear is varied from 0 to 0.15. At Re=300, the head of the hairpin vortex loop always locates on the high-velocity side in uniform shear, and the flow maintains the planar symmetry. At Re=425 and 480, the irregularity in the location and strength of the hairpin vortex appearing in uniform inlet flow is much reduced in uniform shear, but the flows still keep the asymmetry for most inlet shear rates. However, in the cases of K=0.075 and 0.1 at Re=425, the flows become planar symmetric and their characteristics of the evolution of the hairpin vortex loops are different from those of asymmetric flows. A hysteresis phenomenon switching from the planar symmetry to the asymmetry (or vice versa) depending on the initial condition is also observed at Re=425.

Diffuser is an important fluid-mechanical equipment to convert kinetic energy into pressure energy. Many of the experimental and theoretical researches have been done in a diffuser but the understanding of energy transfer and detailed mechanism of energy dissipation is unclear. In this study, computations were performed using a numerical method with SIMPLE algorithm for conical diffuser with various diffuser angles and diffuser lengths. Also, we investigated the pressure recovery coefficient in conical diffuser by inserting strut. In this paper, we showed that the strut can cause a rising pressure recovery.

Experimental study of separated flow over a NACA0012 airfoil is conducted at $Re=2{\times}10^5$ when periodic wakes pass over the airfoil. The wakes are periodically generated by circular cylinders upstream of the airfoil. The measurement of surface pressure and surface visualization at various angles of attack are carried out without and with passing wakes. Without passing wakes, a separation bubble at the leading edge of the suction surface is formed at an angle of attack, found from a local plateau in the streamwise pressure distribution and two distinct lines in the surface flow visualization. With passing wakes, however, the bubble disappears. Owing to passing wakes, the lift increases at high angle of attack and the angle of stall also increases.

An experimental study was carried out to investigate influence of flow conditions on a boundary layer to the near-wake of a flat plate. The flow condition in the vicinity of trailing edge that is influenced by upstream condition history is an essential factor that determines the physical characteristics of a near-wake. Various tripping wires were used to change boundary layer flow condition of upstream at the freestream velocity of 6.0 m/sec. Measurements of the boundary layer and near-wake according to the change of upstream conditions were conducted by using both I-probe(55P14 for boundary layer) and X-probe(55P61 for wake). Normalized velocity profiles of the boundary layer were shown the flow types such as laminar boundary layer, transition, and turbulent boundary layer at 0.95C from the leading edge. The velocity and turbulence intensity profiles of the near-wake for the case of laminar boundary layer at the flat plate surface exhibited a defect and a double peak showing perfect symmetry, respectively.

The supersonic impinging jet has been extensively applied to rocket launching system, gas jet cutting control, gas turbine blade cooling, etc. In such applications, wall temperature of an object on which supersonic jet impinges is a very important factor to determine the performance and life of the device. However, wall temperature data of supersonic impinging jets are not enough to data. The present study describes an experimental work to measure the wall temperatures of a vertical flat plate on which supersonic, dual, coaxial jet impinges. An Infrared camera is employed to measure the wall temperature distribution on the impinging plate. The pressure ratio of the jet is varied to obtain the supersonic jets in the range of over-expanded to moderately under-expanded conditions at the exit of coaxial nozzle. The distance between the coaxial nozzle and the flat plate was also varied. The coaxial jet flows are visualized using a Shadow optical method. The results show that the wall temperature distribution of the impinging plate is strongly dependent on the jet pressure ratio and the distance between the nozzle and plate.

Experimental and computational studies were performed to investigate the effectiveness of a thrust vectoring method using a counterflow concept. A shadowgraph method was used to visualize the supersonic jet expanded from a two-dimensional convergent-divergent nozzle and deflected by a now suction. The primary nozzle pressure and suction nozzle pressure ratios are varied between 3.0 and 5.0, and between 0.2 and 1.0 respectively. The present experimental and computational results showed that, for a given primary nozzle pressure ratio, a decrease in the suction nozzle pressure ratio produced an increased thrust vector angle, and during the change processes of the suction pressure, a hysteresis effect of the thrust vectoring was found through the wall pressure distributions.

When a shock wave propagates into a Helmholtz resonator, very complicated wave phenomena are formed both inside and outside the resonator tube. Shock wave reflection, shock focusing phenomena and shock-vortex interactions cause strong pressure fluctuations inside the resonator, consequently leading to powerful sound emission. In the present study, the wave phenomena inside and outside the Helmholtz resonator are, in detail, investigated with a help of CFD. The Mach number of the incident shock wave is varied below 2.0 and several types of resonators are tested to investigate the influence of resonator geometry on the wave phenomena. A TVD scheme is employed to solve the axisymmetric, compressible, Euler equations. The results obtained show that the configuration of the Helmholtz resonator significantly affects the peak pressure of shock wave focusing, its location, the amplitude of the discharged wave and resonance frequency.

The near field structures of an under-expanded, dual, coaxial, jets issuing from the coaxial nozzles with four different geometries are visualized by using a shadowgraph optical method. Experiments are conducted to investigate the effects of the nozzle-lip thickness, secondary stream thickness, the nozzle pressure ratio and the secondary swirl stream on the characteristics of under-expanded jets. The results show that the presence of secondary annular swirling stream causes the Mach disk to move further downstream and increases its diameter, which decreases with a decrease in the nozzle-lip thickness. The secondary stream thickness has an influence on the location of an annular shock wave.

Theoretical study has been conducted to clarify pressure characteristics of KTX(Korea Train eXpress) in tunnel. The external and internal pressure of rolling stock have been measured by using the atmospheric pressure sensors and portable data acquisition system on Seoul-Busan high speed railroad line. These pressure change may give rise to the ear-discomfort for passenger and fatigue for car body. In this study, the tunnels from 200m to 4000m in length have been chosen for the investigation of tunnel length effects. From the results of experiment, the pattern of pressure change generally agrees to RTRI's experimental result for Shinkansen. We found that there are similar patterns of external pressure variation for each critical tunnel length. The critical tunnel lengths are governed by train speed, train length and sonic velocity.

Hard disk drives (HDD) in computer are used extensively as data storage capacity. The trend in the computer industry to produce smaller disk drives rotating at higher speeds requires an improved understanding of fluid motion in the space between disks. Laser sheet and digital camera was used for 2-dimensional visualization of the unsteady flow between co-rotating disks in air with a cylindrical enclosure (or shroud). Geometric parameters are gap height (H) between disks, and gap distance (G) between disk tip and shroud. The lobe-structured boundary between inner region and outer region was detected by inserted particles, and the number of dominant vortices was determined clearly It is found from flow visualization that the number of vortex cells can be correlated with Reynolds number based on H which is defined as $Re_H={\Omega}RH/v$ ranging from $7.96{\times}10^2$ to $1.43{\times}10^4$, and decreases as the disk speed increases. The lobe pattern by vortex cells is changed to a circular pattern for the wide gap than narrow one.

In this study, a experimental method has been introduced for the various exhaust pipe geometry of 4-stroke single engine. The main experimental parameters are the variation of exhaust pipe diameters and lengths, to measuring the pulsating flow when the intake and exhaust valves are working, As the results of experimental test, the various exhaust geometry were influenced strongly on the exhaust pressure. As the exhaust pipe diameter was decreased, the amplitude and the number of compression wave in exhaust pressure was increased. According to decreasing pipe diameter, the number of compression wave in exhaust pressure was decreased. When the pipe diameter was increase, the second amplitude was increased.

Leading edge extension(LEX) in a highly swept shape applied to a delta wing features the modem air-fighters. The LEX vortices generated upon the upper surface of the wing at high angle of attack enhance the lift force of the delta wing by way of increased negative suction pressure over the surfaces. The present 3-D stereo PIV includes the Identification of 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 criterion and so on. A delta wing model with or without LEX was immersed in a circulating water channel. Two high-resolution, high-speed digital cameras($1280pixel{\times}1024pixel$) were used to allow the time-resolved animation work. The present dynamic stereo PIV represents the complicated vortex behavior, especially, in terms of time-dependent characteristics of the vortices at given measuring sections. Quantities such as three velocity vector components, vorticity and other flow information can be easily visualized via the 3D time-resolved post-processing to make the easy understanding of the LEX effect or vortex emerging and collapse which are important phenomena occurring in the field of delta wing aerodynamics.

The light-transmission technique has been applied to a slit rheometer for measuring red blood cell aggregation as well as blood viscosity over a range of shear rates. For measurement of blood viscosity and aggregation, instantaneous pressure and transmit-light intensity are measured with time. Using a precision pressure measurement, one can determine the shear stress and shear rate. In addition, a transmitted light through a blood sample indicates degree of RBC aggregation. With abruptly flowing with high shear rate, RBCs rapidly disaggregate and the intensity of the transmitted light becomes low. When continuously flowing with decreasing shear rate, RBCs tend to re-aggregate and the corresponding transmit-intensity gradually increases with time. The light intensity as a degree of RBC aggregation is plotted against shear rate and compared with blood viscosity. The advantages of this design are dual measurement at a time, simplicity, i.e., ease of operation and no moving parts, low cost, short operating time, and the disposable kit which is contacted with blood sample.

A slit-flow apparatus with laser diffraction method has been developed with significant advances in ektacytometry design, operation and data analysis. In the slit-flow ektacytometry (or laser-diffractometry), the deformation of red blood cells subjected to continuously decreasing shear stress in slit flow is measured. A laser beam traverses a diluted blood suspension flowing through a slit and is diffracted by RBCs in the volume. The diffraction patterns are captured by a CCD-video camera, linked to a frame grabber integrated with a computer, while the differential pressure variation is measured by a pressure transducer. Both measurements of laser-diffraction image and pressure with respect to time enable to determine deformation index and the shear stress. The range of shear stress of 0 ${\sim}$ 35 Pa and measuring time is less than 2 min. When deforming under decreasing shear stress, RBCs change gradually from the prolate ellipsoid towards a circular biconcave morphology. The Deformation Index (DI) as a measure of RBC deformability is determined from an isointensity curve in the diffraction pattern using an ellipse-fitting program. The advantages of this design are simplicity, i.e., ease of operation and no moving parts, low cost, short operating time, and the disposable kit which is contacted with blood sample.

CT 사진을 이용하여 살아있는 사람의 비강 형상을 얻고 RP 를 이용하여 모형을 만들고 이 모형 으로 정교한 비강 유로 모형을 제작하였고, 호흡을 정확하게 모사하는 펌프를 만들었다. 사람의 호흡 데이터를 이용하여 캠을 제작하고 대형 피스 톤 펌프를 만들어 사람의 호흡을 정확 하게 모사 하였다. 이를 이용하여 생리적 주기를 갖는 비강 내 유동에 대한 결과를 획득하고 기존의 일정유량 하의 실험결과와 비교 분석하였다.

Synchrotron X-ray micro-imaging technique was employed to monitor non-invasively the refilling process of water inside the xylem vessels in bamboo leaves. The consecutive phase-contrast X-ray images clearly show both plant anatomy and the transport of water inside the xylem vessels. Traces of water-rise, vapor bubbles and variations of contact angle between the water front and the xylem wall were measured in real time. During the refilling process, air bubbles are removed when the rising water front halts at a vessel end for a while. Subsequently, it starts rising again at a higher velocity than the normal refilling speed. Repeated cavitation seems to deteriorate the refilling ability in xylem vessels. In dark environment, the water refilling process in xylem vessels is facilitated more effectively than in bright illuminated conditions. Finally, X-ray micro-imaging was famed to be a powerful, high resolution, real time imaging tool to investigate the water refilling process in xylem vessels.

The waterhammer has recently become more important because the pumping stations were big and the systems conveying the fluid through the large and long transmission pipelines were complex. When the pumps are started or stopped for the operation or tripped due to the power failure. the hydraulic transients occur as a result of the sudden change in velocity. In this paper, the field tests on the waterhammer by the startup, stoppage, and power failure of a centrifugal pump were carried out for Yongma transmission pumping station in Seoul. The experimental results were compared with that of the numerical calculations. in which results the procedure of controlled pump normal shut-down and the two-step closing mode of controlling the ball valve for pump emergency stop are proposed to reduce the pressure surge.

In this study, we are investigated experimentally the pumping characteristics about the pumping channel shapes of disk-type drag pump (DTDP). We are experimented the pumping performance about the rotors which have channel or do not exist. The channel disk-type rotor has spiral channels both upper and lower part, and stator is planar. The planar disk-type rotor hasn't channel and stator has spiral channels both upper and lower part. The flow-meter method is adopted to calculate the pumping speed. Compression ratio and pumping speeds for the nitrogen gas are measured under the inlet pressure range of 0.001 ${\sim}$ 4 Torr. The maximum of compression ratio was about 3300 for three-stage DTDP (channel disk-type rotor), 1000 for four-stage (planar disk-type rotor) and two-stage DTDP (channel disk-type rotor) at zero throughput. The ultimate pressure was $1.6{\times}10^{-6}$ Torr for three-stage DTDP (channel disk-type rotor), $2.5{\times}10^{-6}$ Torr for four-stage DTDP (planar disk-type rotor).

Pressurization system of a 50m unmanned airship was manufactured according to the detailed design. In this paper. the whole procedures including system design. simulation and fabrication were described. The fabricated part were ground tested to check compliances with design requirements. Ground tests include operational tests, leakage tests, endurace tests and low temperature environmental tests. Results shows pressurization system of a 50m unmanned airship meets design requirements. Currently. pressurization system is installed to the KARI airship, Via 50m and performance verification through the flight tests are being conducted.

An in-house supersonic ejector was designed to ensure low pressure and high speed scavenging of resonating cavity of chemical lasers. For given primary flow condition, 100g/s secondary mass flow rate was observed at the design pressure. Performance validation of a supersonic ejector system along with an investigation of effects of supersonic diffuser was conducted. Placement of diffuser at the secondary inlet further reduced diffuser upstream pressure to 1/4-1/5 relieving the local to the primary supply unit. In order to increase the secondary flow, we put two ejectors capable of removing 50g/s each of secondary flows together to deal with higher mass flow. Test of the parallel unit demonstrated the secondary flow rate was proportional to the numbers of individual units that were brought together. Additionally, flow calculations with a commercial code were carried out in every case of experiment and compared with results.

The present study addresses a computational work to investigate the influence of a turbulent wake flow on the pressure recovery of a two-dimensional subsonic diffuser. The turbulent wake is generated by a cylinder with a small diameter, which is installed at the diffuser inlet. Computation is applied to two-dimensional steady Navier-Stokes equations. The computational results are qualitatively well compared to existing experimental data. The results show that the diffuser pressure recovery is strongly dependent on the diameter and location of the cylinder. It is found that there is a certain diameter and location of cylinder for the diffuser pressure recovery to be most enhanced. Compared with no cylinder case, the diffuser performance increases up 24%.

This paper presents the response surface optimization method using three-dimensional Navier-Stokes analysis to optimize the shape of a forward-curved blades centrifugal fan. For numerical analysis, Reynolds-averaged Navier-Stokes equations with $k-{\varepsilon}$ turbulence model are discretized with finite volume approximations. In order to reduce huge computing time due to a large number of blades in forward-curved blades centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models. Three geometric variables, i.e., location of cut off, radius of cut off, and width of impeller, and one operating variable, i.e., flow rate, were selected as design variables. As a main result of the optimization, the efficiency was successfully improved. And, optimum design flow rate was found by using flow rate as one of design variables. It was found that the optimization process provides reliable design of this kind of fans with reasonable computing time.

In this paper, a small supersonic wind tunnel is designed and built to study the flow characteristics of a supersonic impulse turbine cascade. The flow is visualized by means of a single pass Schlieren system. The supersonic cascade with 3-dimensional supersonic nozzle was tested over a wide range of nozzle installation angle. Highly complicated flow patterns including shocks, nozzle-cascade interaction and shock boundary layer interactions are observed.

Turbo-pump system, an essential component of liquid rockets and induced weapons, adopts a partial admission axial turbine which drives pump. And the turbine of a turbo-pump system is usually operated at supersonic condition due to its high loading chracteristics. Therefore, reseaches about flow and performance characteristics of a partial admission supersonic turbine must be preceeded to progress the aerospace and defense industries as well as the development of turbo-pump systems. In this study, flow characterisitics within blades of the partial admission supersonic turbine are numerically investigated by using Fine Turbo, a commercial CFD Code. Before performing the numercial analyses, to verify accuracy of the numerical result computed by Fine Turbo, I performed the comparison between the numerical results with J.J.Cho' experimental results. It is found that the numerical results show good agreement with the experimental results. Computations about the partial admission supersonic turbine have been performed to investigate flow characteristics including shock patterns. It is also found that the flow and performance of partial admission supersonic turbine are largely depend on shocks ocurred in the nozzle and at the leading edge of blades, expansion or compression at exit of nozzle and separations occurred in passage.

It is important to measure precisely the size and velocity of micro-bubbles used in various field. The synchrotron X-ray micro-imaging technique was employed to measure the size and velocity of micro-bubbles moving in an opaque tube simultaneously. Phase contrast images were obtained at interfaces of micro-bubbles between water and air due to their different refractive indices. The X-ray micro-imaging technique was found to measure an optical fiber with an accuracy of 0.2%. Micro-bubbles of $10{\sim}60{\mu}m$ diameter moving upward in an opaque tube (${\phi}=2.7mm$) were tested to measure bubble size and up-rising velocity. For DI water, the measured velocity of micro-bubbles is nearly proportional to the square of bubble size, agreed well with the theoretical result. In addition, the synchrotron X-ray micro-imaging technique can measure accurately the size and velocity of several overlapped micro-bubbles.

In order to measure non-intrusively velocity profile in liquid metal flow, a modified electromagnetic flowmeter was designed, which was based on electromagnetic tomography technique. Under the assumption that flow is fully-developed, axisymmetric and rectilinear, the velocity profile was reconstructed after the flowmeter equation, the first kind of Fredholm integration equation, was linearized. In reconstruction process Tikhonov regularization method with regularization parameter was used. The reconstructed velocity profile had the nearly same as turbulent flow profile which was approximately represented as log law. In addition, flowmeter output for a fixed magnet rotation angle was linearly proportional to flow rate. When magnet rotation angle was $54^{\circ}$, axisymmetric weight function was nearly uniform so that the flowmeter gives a constant signal for any fully-developed, axisymmetric and rectilinear profile with a constant flow rate.

The holographic velocimetry system has a significant potential for the measurements of three dimensional velocities of particles. In this study, orthogonal two-side holography system was developed to obtain three dimensional velocities and sizes of spray droplets. To get high quality of reconstructed images, singe-exposure holography at two time moments and two orthogonal sides was adopted instead of multi-exposure, single-side holography. From three dimensional positions of droplets determined by reconstruction and image processing system, the three dimensional velocities and sizes of each droplet was extracted using the PTV algorithm. To determine the position of particles in the optical axis, a new focusing parameter was introduced based on the correlation between two droplet images at the same distance. The measured results by holography system were compared with those by the PDPA.

Recently, Plastic optical products are widely used. Injection molding process has advantages of low cost and high productivity. However, it remains a residual birefringence and residual stress by difference cooling. The present study focused on measuring birefringence in optical plastic parts using interference color pattern. The main idea of an analysis comes from chromatic aberration which is caused by difference light wavelengths. As a result, a complete system measuring the high order 2-D birefringence pattern was built. Further investigation is under way to improve the accuracy of birefringence measurement system by diode laser.

A 4D-PTV system was constructed. The measurement system consists of three high-speed high-definition cameras, Nd-Yag laser(10mJ, 2000fps) and a host computer. The GA-3D-PTV algorithm was used to extract three-dimensional velocity vectors in the measurement volume. A horizontal impinged jet flow was measured with the constructed system. The Reynolds number is about 40,000. Spatial temporal evolution of the jet flow was examined in detail and physical properties such as spatial distributions of vorticity and turbulent kinetic energy were obtained with the constructed system.

The performances of the newly constructed Stereoscopic PTV and Stereoscopic PIV which had been completed based on a 3D-PTV principle are tested using the Standard Images. Virtual images were produced for the benchmark tests of the constructed two Stereoscopic techniques. The arrangement of the two cameras was based on angular position. The calibration of cameras and the pair-matching of the three-dimensional velocity vectors were based on Genetic Algorithm based 3D-PTV technique. The Standard data of LES are on the impinged jet proposed by VSJ. It is shown that the results obtained by Stereoscopic PTV have better reliability than those by Stereoscopic PTV.

A new algorithm, CBIR (Correlation-Based Image Registration) was proposed to improve the resolution of image registration for PSP (Pressure-Sensitive Paint). The local displacement vectors were obtained by finding the displacement which maximizes the cross-correlation between two interrogation windows of 'wind-off' and 'wind-on' images. A recursive multigrid processing was employed to increase the non-linear spatial resolutions. The variations of image were precisely measured without identifying the control points.

In the present paper, a micro flow sensor, which can be used at bio-delivery systems and micro heat pumps, is developed. For this, the micro flow sensor is integrated on a quartz wafer ($SiO_2$) and is manufactured by simple and convenient microfabrication processes. The micro flow sensor aims for measuring mass flow rates in the low range of about $0{\sim}20$ SCCM. The micro flow sensor is composed of temperature sensors, a heater, and a flow microchannel. The temperature sensors and the heater are manufactured by the sputtering processes in this study. In the microfabrication processes, stainless steel masks with different patterns are used to deposit alumel and chromel for temperature sensors and nichrome for the heater on the quartz wafer. The microchannel is made of Polydimethylsiloxane(PDMS) easily. A deposited quartz wafer is bonded to the PDMS microchannel by using the air plasma. Finally, we confirmed the good operation of the present micro flow sensor by measuring flow rate.

The present study investigates flow characteristics in an optical disc drive. Detailed knowledge of the flow characteristics is essential to analyze flow-induced noise and vibration, forced convection and flow friction loss. The ODD used in the personal computer is used for the experiment and rotating velocity of disc is under the 4500 rpm. Time-resolved velocity components and velocity spectrum are obtained using the laser Doppler anemometry (LDA). The results show that the front holes reduce now-induced noise and the position of pickup body affects flow near the window. In addition, il is possible for cooling of heat sources in an optical disc drive through measuring the flow fields under the tray.

DOS and NaCl aerosol particles with geometric mean diameter of $0.1{\sim}3.0{\mu}m$ geometric standard deviation of $1.1{\sim}1.8$ and total number concentration of $450{\sim}400$ $particles/cm^3$ were used to determine collection efficiencies of a duct-type wet scrubber with respect to particle size. The tested operating variables included air velocity and water injection rate. It was shown from the experimental results that the collection efficiencies increased with increasing water injection rate and decreasing air velocity. It was also seen that the collection efficiency of the Duct-type wet scrubber is mainly governed by the mechanism of inertial impaction.

Stochastic nature of subgrid-scale stress causes the predictability problem in large eddy simulation (LES) by which the LES solution field decorrelates with field from filtered directnumerical simulation (DNS). In order to evaluate the predictability limit in a priori sense, the information on the interplay between resolved scale and subgrid-scale (SGS) is required. In this study, the analysis on the inter-scale interaction is performed by applying tophat and cutoff filters to DNS database of flow over a circular cylinder at Reynolds number of 3900. The effect of filter shape is investigated on the interpretation of correlation between scales. A critique is given on the use of tophat filter for SGS analysis using DNS database. It is shown that correlations between Karman vortex and SGS kinetic energy drastically decrease when the cutoff filter is used, which implies that the small scale universality holds even in the presence of the large scale coherent structure.

Sinusoidal wavy channel is one of the most commonly used devices in the industry for achieving mixing and heat transfer. Here we report on results obtained from the DNS of flow inside the wavy channel performed using the finite volume technique. As a primary stage to obtain the optimal design for heat transfer and mixing, this study observed the basic flow structures in a wavy channel. The mass flow rate is kept constant with friction Reynolds number of $Re_{\tau}$ = 140 . Time- and space-averaged and instantaneous flow fields are illustrated to observe the flow structures. Although the direct comparison of results between turbulent wavy and flat channel is somehow difficult due to the different flow phenomena derived from different configuration, here the mean streamwise velocity and RMS of velocities at same $Re_{\tau}$ of two different channels are compared. The basic difference between wavy and flat channel flow is the existence of small scale wall vortices along the walls in a wavy channel. These vortices make flow more complex, which will accompany the increase of heat transfer, pressure drop and drag.

Discrete Wavelet Transform (DWT) has been applied to the Direct Numerical Simulation (DNS) data of turbulent channel flow. DWT splits the turbulent flow into two orthogonal parts, one corresponding to coherent structures and the other to incoherent background flow. The coherent structure is extracted from not vorticity field but velocity's since the channel flow is not isotropic. By comparing DWT's result of channel flow with that of isotropic flow, it is shown that coherent structure maintains the properties of original channel flow. The velocity field of coherent structures can be represented by few wavelet modes and that these modes are sufficient to reproduce the velocity probability distribution function (PDF) and the energy spectrum over the entire inertial range. The remaining incoherent background flow is homogeneous, has small amplitude, and is uncorrelated. These results are compared with those obtained for the same compression rate using large eddy simulation (LES) filtering. In contrast to the incoherent background flow of DWT, the LES subgrid scales have a much larger amplitude and are correlated, which makes their statistical modeling more difficult.

The common method to improve heat transfer in Nuclear fuel rod bundle is install a mixing vane in space grid. The previous split mixing vane is guides cooling water to swirl flow in sub-channel of fuel assembly. But, this swirl flow decade rapidly after mixing vane and the effect of enhancing the heat transfer vanish behind this short region. The large scale secondary vortex flow was generated by rearranging the inclined angle direction of mixing vanes to the coordinated directions. This LSVF mixing vanes generate the most strong secondary flow vortices which maintain about 35 $D_H$ after the spacer grid and the streamwise vorticity in subchannel with LSVF mixing vane sustain two times more than that in subchannel with split mixing vane. The turbulent kinetic energy and the Reynolds stresses generated by the mixing vanes have nearly same scales but maintain twice more than previous type.

A numerical simulation was performed of flow behind a squareback car with a rear spoiler. Influence of the rear spoiler on drag force has been studied. A lattice Boltzmann method was utilized to portray the unsteady aerodynamics of wake flows. The pressure distributions were employed to examine the vortex formation mode against the rear spoiler. It was found that the separation flow at roof end and c-pillar makes three dimensional vortex structures and the rear spoiler increases pressure on the rear glass surface.

In this paper, we present scavenge characteristics in a small subchamber of HCCI. It is very important to enhance scavenge rate because ignition in a chamber sometimes does not happen. To understand this phenomenon numerical tool was performed using the FLUENT which is a commercial code. Focus is given to the effect on the scavenge rate of the geometric factor that is the angle of nozzle injection. The numerical results show that the scavenge ones in the subchamber heavily depend on the nozzle angle. It was found that the scavenge rate is more effective at angled nozzle.

The HANARO, a multi-purpose research reactor of 30 MWth, open-tank-in-pool type, has been under normal operation since its initial criticality in February, 1995. The HANARO is composed of inlet plenum, grid plate, core channel with flow tubes and chimney. The reactor core channel is located at about twelve meters (12 m) depth of the reactor pool and cooled by the upward flow that the coolant enters the lower inlet of the plenum,. rises up through the grid plate and the core channel and comes out from the outlet of chimney. A guide tube is extended from the reactor core to the top of the reactor chimney for easily un/loading a target under the reactor normal operation. But active coolant through the core can be quickly raised up to the top of the chimney through the guide tube by a jet flow. This paper describes an analytical analysis that is the study of the flow behavior through the guide tube under reactor normal operation and unloading the target. As results, it was conformed through the analysis results that the guide jet is suppressed under the top of the chimney after modifying the orifice diameter of 37.5 mm to 31 mm.

'Mixing Index($D_I$)'s generally used to measure the degree of mixing. A new method to calculate $D_I$ was proposed, when insoluble solution flows in micromixer. 'Vortex Index (${\Omega}_I$)'which indicate the degree of chaotic advection, is defined and formulated. A lots of arbitrary shaped microchannels were tested to calculate the $D_I$ and ${\Omega}_I$. And then a simple algebraic equation, $D_I=A{\Omega}_I+B$, was obtained. This equation may be used instead of partial differential equation, concentration equation.

The laminar impinging jet flow fields were investigated with or without magnetic fields. The transient phenomenon from steady to unsteady flow was founded at specific Reynolds number ranges. In unsteady flow region, the magnetic fields make flow stable. So the characteristics of fluid flow at impingement wall are changed

Sedimentation phenomenon of suspended solids occurs by the gravitational force. Pollution particles are separated from slowly flowing wastewater in clarifier. Recently, the sludge suction collector is preferred rather than the scraper type sludge collector due to enhancement of the clarifier efficiency. The sludge suction collector is usually operated by the user's experiences without any scientific and/or technical consideration. To evaluate the performance of sludge suction collector, the three dimensional computer simulation was conducted by the finite volume method. To analyze the performance, the velocity vectors and the suction flow rates of the orifices were investigated. The result showed that each suction flow rate through out the collector was equivalent in the sludge suction collector and the efficiency of suction collector was good to remove high concentrated sludge in clarifier.

EK pumps packed with particles inside capillaries are involved in the mixed electroosmotic flow and pressure driven flow. For analysis in the porous EK pumps, the volume-averaging technique is applied to derive the volume-averaged equations for momentum and electrical potential. By using the volume-averaged equations, analytical solutions for electric potential and velocity distribution due to the mixed electroosmotic and pressure driven flows are obtained. The present analysis is validated by comparison with numerical and experimental results for the case of microchannel EK pumps.

The present study concerns a computational study of fully developed laminar flow of a Newtonian fluid through an eccentric annulus with a combined bulk axial flow and inner cylinder rotation. This study considers the identical flow geometry as in the calculation of Escudier et $al.^{(3)}$ An unexpected feature of the calculations for eccentricity ${\varepsilon}$)0.7 is the appearance of a second peak in the axial velocity, located in the narrowing gap. The distribution of the axial component of the surface shear stress has a maximum in the narrowing gap and a minimum in the widening gap.

Acoustic pressure field around the cross-flow fan is predicted by a hydrodynamic-acoustic splitting method. Unsteady flow field is obtained by solving the incompressible Navier-Stokes equations using an unstructured finite-volume method on the triangular meshes, while the acoustic waves generated inside the cross-flow fan are predicted by solving the perturbed compressible equations(PCE) with a 6th-order compact finite difference method. Computational results show that the acoustic waves of BPF tone are generated by interactions of the blades wakes with the stabilizer, which then are reflected from the rear-guider and mainly propagate towards the fan inlet. Also, a directivity of BPF noise predicted by the PCE is noticeably different from that of the FW-H equations, in which a fan casing effect cannot be included.

This paper deals with the evaluation of several boundary conditions which are commonly used in the lattice Boltzmann equation method. 2-D channel flow(poiseui1le flow) and lid-driven cavity flow was selected as a test problem of this study, because there exist an analytic solution and previous study which could be used for a benchmarking test. It was found that lattice Boltzmann method still needs more considerations of stability and physical consistency, though it could predict the flow patterns both qualitatively and quantitatively.

To analyze the fluid flow and thermal characteristics in a nanoscale system, the planar Poiseuille flow of a Lennar-Jones liquid through parallel plates formed by fixed atoms is studied using nonequilibrium molecular dynamics simulations. The role of important simulation parameters such as the channel width, the magnitude of external field, the temperatures of the top and bottom plates, and the interaction potential parameter between fluid and wall atoms, which affect flow patterns and heat transfer rate inside the channel, are investigated. Under the various simulation conditions, interesting phenomena deviated from the continuum predictions have found.

In this paper the numerical method with a zonal embedded grid system was applied to the spin-up flow within a semi-circular container. Flow visualization was also performed on a rotating table. The results show that at a low level of damping (i.e. low viscosity and high liquid depth) a cyclonic cell produced initially near the left-hand-side corner of the container moves along a wall and merged with the cell on the right-hand-side.

The present article proposes a new droplet collision model including the stretching separation regime and the formation of satellite droplets. The new model consists of a several equations to calculate the post-collision characteristics of colliding droplets and satellite droplets. These equations are derived from the energy balance of droplets between before and after collision. For binary collision of water droplets, the new model shows good agreement with experimental data for the number of satellite droplets. Nevertheless, it is thought that, in order to guarantee the generality of the new model, the improvements should be performed to consider the effects of the bouncing and the reflexive separation, which is essential process in the collision of hydrocarbon droplets.

Numerical analysis of the heat transfer associated with droplet impact on a hot solid surface is performed by solving the mass, momentum and energy equations for the liquid-gas region. The deformed droplet shape is tracked by a level set method which is modified to achieve volume conservation during the whole calculation procedure and to include the effect of contact angle at the wall. The numerical method is validated through test calculations for the cases reported in the literature. Based on the numerical results, the effects of advancing/receding contact angle, impact velocity and droplet size on the heat transfer during droplet impact are quantified.

A preconditioned numerical method for gas-liquid two-phase flows is applied to solve cavitating flow. The present method employs a finite-difference dual time-stepping integration procedure and the MUSCLTVD scheme. A homogeneous equilibrium cavitation model is used. The present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation, large density changes and incompressible flow characteristics at low Mach number. Some internal flows such as convergent-divergent nozzles are computed using this method. Comparisons of predicted and experimental results are provided and discussed.

An immersed boundary method for simulation of density-stratified flows is developed and applied to computation of viscous flows over two-dimensional obstacles in a bounded domain under stable density stratification. Density sources/sinks are introduced on the body surface. Two obstacle shapes are used, a vertical barrier and a smooth cosine-shaped hill; weak stratification, defined by $K=ND/{\pi}U{\leq}1$, where U, N, and D are the upstream velocity, buoyancy frequency, and domain height, respectively, is considered. The results are consistent with other authors' calculations, and shed light on computation of density-stratified flows in complex geometries.

An edge tone is the discrete tone or narrow-band sound produced by an oscillating free shear layer impinging on a rigid surface. In this paper we present a two-dimensional edge tone to predict the frequency characteristics of the discrete oscillations of a jet-edge feedback cycle by the finite difference lattice Boltzmann method. We use a new lattice BGK compressible fluid model that has an additional term and allow larger time increment comparing a conventional FDLB model, and also use a boundary fitted coordinates. The jet is chosen long enough in order to guarantee the parabolic velocity profile of the jet at the outlet, and the edge consists of a wedge with an angle of ${\alpha}=23^{\circ}$ . At a stand-off distance ${\omega}$ , the edge is inserted along the centreline of the jet, and a sinuous instability wave with real frequency f is assumed to be created in the vicinity of the nozzle and to propagate towards the downstream. We have succeeded in capturing very small pressure fluctuations result from periodically oscillation of jet around the edge. That pressure fluctuations propagate with the sound speed. Its interaction with the wedge produces an irrotational feedback field which, near the nozzle exit, is a periodic transverse flow producing the singularities at the nozzle lips.

The tone generation mechanism and aeroacoustic characteristics have been investigated for flow over open cavities using direct acoustic numerical simulations. Physically the tone generation mechanism of open cavity is more complicated when flow instabilities are excited by the correlation effects of flow parameters. From non-dimensional parameter studies in very low Mach number range, it is shown that characteristics of cavity resonance inherently involve typical acoustic pattern at each discrete tone frequency, and especially in laminar flow the fundamental tone frequency is determined within flow instability criterion of laminar shear layer as well as cavity geometry, length to depth ratio.

A 1/21.6 scaled non-heating experimental facility was prepared utilizing the results of a scaling analysis to simulate the APR1400 reactor and insulation system. The behaviors of the air bubble-induced two-phase natural circulation flow in the insulation gap were observed, and the liquid mass flow rates driven by natural circulation loop were measured by varying the injected air flow rate and distribution. As the injected air flow rates increased, the natural circulation flow rates also increased. Both the longitudinal and the latitudinal distributions of the injected air affected the natural circulation flow rates, especially, the longitudinal effect is more larger.

Nucleate boiling experiments with constant wall temperature of heating surface were performed using R113 for almost saturated pool boiling conditions. A microscale heater array and Wheatstone bridge circuits were used to maintain a constant wall temperature condition and to measure the heat flow rate with high temporal and spatial resolutions. Bubble images during the bubble growth were taken as 5000 frames a sec using a high-speed CCD camera synchronized with the heat flow rate measurements. The geometry of the bubble during growth time could be obtained from the captured bubble images. The bubble growth behavior was analyzed using the new dimensionless parameters for each growth regions to permit comparisons with previous results at the same scale. We found that the new dimensionless parameters can describe the whole growth region as initial and later respectively. The comparisons showed good agreement in the initial and thermal growth regions. The required heat flow rate for the volume change of the observed bubble was estimated to be larger than the instantaneous heat flow rate measured at the wall. Heat, which is different from the instantaneous heat supplied through the heating wall, can be estimated as being transferred through the interface between bubble and liquid even with saturated pool conditions. This phenomenon under a saturated pool condition needs to be analyzed and the data from this study can supply the good experimental data with the precise boundary condition (constant wall temperature).

Recently, several kinds of experimental and computational studies are being carried out to investigate the WIG aerodynamic characteristics which are of practical importance to develop the new ground transportation vehicle system. These works are mainly based upon conventional wind tunnel tests, but many problems associated with the WIG aerodynamic characteristics cannot be satisfactorily resolved due to the wind tunnel blockage effects or string problems to support the test object. To do this, it is necessary to develop a novel simulator appropriate to the WIG aerodynamics. The objective of the present study is to clarify the aerodynamic characteristics of a new developed WIG simulator, which is able to imitate real WIG flow circumstances such as gradually decelerating and accelerating flows.

In this paper, The performance of Kim-Chi refrigerator with three evaporator and one compressor was investigated in employing 55% propane and 45% isobutane (R290/R600a) refrigerant mixture as an alternative refrigerant of R134a. The drop in test was performed by varying both refrigerant charge amount and capillary tube length in order to find both the performance and reliability of a small multi-refrigeration system. As a result, Both the power consumption and COP is increased by about 15% and 10%, respectively as compared to the baseline R134a system. In addition, the propane/isobutene refrigerant mixture system took advantage of the minimization of modification and redesigning of system components because of similar thermodynamic properties with R134a such as saturation pressure, temperature, normal boiling point(NBP) characteristics

The transient nature and complex flow geometries of two-phase gas-liquid flows cause fundamental difficulties when measuring flow velocity using an electromagnetic flowmeter. Recently, a current-sensing flowmeter was introduced to obtain measurements with high temporal resolution (Ahn et $al.^{(1)}$). In this study, current-sensing flowmeter theory was applied to measure the fast velocity transients in slug flows. To do this, the velocity fields of axisymmetric gas-liquid slug flow in a vertical pipe were obtained using Volume-of-Fluid (VOF) method and the virtual potential distributions for the electrodes of finite size were also computed using the finite volume method for the simulated slug flow. The output signal prediction for slug flow was carried out from the velocity and virtual potential (or weight function) fields. The flowmeter was numerically calibrated to obtain the cross-sectional liquid mean velocity at an electrode plane from the predicted output signal. Two calibration parameters are required for this procedure: a flow pattern coefficient and a localization parameter. The flow pattern coefficient was defined by the ratio of the liquid resistance between the electrodes for two-phase flow with respect to that for single-phase flow, and the localization parameter was introduced to avoid errors in the flowmeter readings caused by liquid acceleration or deceleration around the electrodes. These parameters were also calculated from the computed velocity and virtual potential fields. The results can be used to obtain the liquid mean velocity from the slug flow signal measured by a current-sensing flowmeter.

The objective of this study is to perform the numerical investigation of flow characteristics in static mixers. Simulations are carried out for mixers consisting of up to six Kenics and PPM elements placed end-to-end at an angle of $90^{\circ}$and for a range of Reynolds number($1{\leq}Re{\leq}100$). The pressure drop across a six-element Kenics mixer is computed and compared with the previous experimental correlations. The results are in good agreement with the previous correlations. The simulated flow field of Kenics mixer is extremely complex and contains regions of transverse flow that is dominated by the interaction of vortices produced by the mixer elements.

A collection efficiency of cyclone is influenced by cut-size and slope of grade efficiency curve. It has been recognized that the collection efficiency is improved when the cut-size is reduced. However, effects of the slope have never been studied so far. In this study, we analyze a relationship between slope of grade efficiency and collection efficiency in two-stage cyclones. In single stage cyclones, higher slope cyclones have high efficiency. On the contrary to single stage cyclone, collection efficiency of two-stage cyclone have the maximum value when the first cyclone has a lower slope and second cyclone is high.

The performance of small cyclone is analysed by an experiment for the purpose of developing a bag-less vacuum cleaner. For the high collection efficiency and low pressure loss cyclone, the effect of cyclone inlet feature must be well understood. Four types of the helical inlet are considered to compare with the normal tangential inlet, and also various inlet velocities are used to each inlet type. Based on the reference dimension, each type of inlet shows the changes of the grade efficiency and pressure loss which determine the cyclone quality. The results show that the helical inlet has the smaller cut-size but bigger pressure loss than the tangential inlet. And the degree of opening area influences factors of cyclone performance. As the inlet velocity is increased, the cut-size becomes smaller and the pressure loss becomes bigger of each cyclone. Further studies are required to understand the optimized helical inlet of cyclone.

The three-dimensional flow in a turbine nozzle guide vane passage causes large secondary loss through the passage and increased heat transfer on the blade surface. In order to reduce or control these secondary flows, a linear cascade with a contoured endwall configuration was used and changes in the three-dimensional flow field were analyzed and discussed. Measurements of secondary flow velocity and total pressure loss within the passage have been performed by means of five-hole probes. The investigation was carried out at fixed exit Reynolds number of $4.0{\times}10^5$. The objective of this study is to document the development of the three-dimensional flow in a turbine nozzle guide vane cascade with modified endwall. The results show that the development of passage vortex and cross flow in the cascade composed of one flat and one contoured endwalls are affected by the flow acceleration which occurs in contoured endwall side. The overall loss is reduced near the flat endwall rather than contoured endwall.

A cross-flow fan 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 exists a rapid pressure head reduction, a noise increase and an unsteady flow against a stabilizer and a rearguider. Moreover, it is difficult to analyze the reciprocal relations of the cross-flow fan because each parameter is independent. Numerical analyses are conducted with different starting angles of the rearguider. Two-dimensional, unsteady governing equations are solved, using FVM, PISO algorithm, sliding grid system and ${\kappa}-{\varepsilon}$ standard turbulence model.

Hemodynamics behavior of the blood flow is influenced by the presence of the arterial stenosis. If stenosis is present in an artery, normal blood flow is disturbed. In the present study, characteristics of steady and pulsatile flow of non-Newtonian fluid, the effects of stenosised geometry are analyzed by numerical simulation. One interesting point is that non-symmetric solutions were obtained at severity stenosis, although the stenosis and the boundary condition were all axisymmetric.

The present study was carried out to develop highly efficient RC ornithopter 'Songgolmae' powered by motor and battery. Designed electric ornithopter weighs 277 grams and has 3 channels radio control. 1t runs on an electric motor by a lithium polymer battery and has a gear ratio of about $75{\sim}95$ to 1 to flap its 88 cm wingspan. The aerodynamic performance of the ornithopter, applied to a flapping motion only, was validated by flight tests. Flight times have exceeded 23 minutes until the battery was used up. The flight test results indicate that the ornithopter developed here has sufficient thrust to propel itself in a forward flight. From the economical point of view and the handling of the RC ornithopter, it can be said that the developed robot ornithopter is an effective RC ornithopter. This radio controlled ornithopter flies its way high into the sky just like a real bird flies.

Experiments were conducted to investigate the flow characteristics of water through rectangular PDMS microchannels with a hydraulic diameter ranging from 66.67 to 200 ${\mu}m$. In the experiments, the flow rate and pressure drop across the microchannels were measured at steady states. The experimental results were compared with the predictions from the conventional laminar flow theory. A significant difference between the experimental data and the theoretical predictions was found. Experimental results indicate that the pressure gradient and flow friction in microchannels are higher than those from the conventional laminar flow theory. This may be attributed to the fact that there exists effect of surface roughness of the microchannels. In this study, a surface roughness model is implemented to interpret the experimental data. A good agreement between the experimental data and the numerical predictions with a surface roughness model were found.

Ultrasonic cleaning is performed by cavitation which is caused by the change of the sound pressure due to the vibration in a cleaning tank. In this study, experiments on electric power and sound pressure with various temperatures, dissolved oxygen and the level of the fluid was done in order to find out how the changes in a cleaning tank affect cavitation. As a result of a series of experiments, we found that transducer impedance changes periodically in response to the variances of fluid and have a direct influence on cleaning efficiency.

In order to develop a high efficiency fire pump, its performance characteristics with various operating conditions are investigated. The governing equations are derived from making using of three-dimensional Navier-Stokes equations with the standard ${\kappa}-{\varepsilon}$ turbulence model and SIMPLE algorithm. Using a commercial code, CFX, pressure distribution and flow fields in a fire pump are calculated with various ranges of rotating speed 800-2400 rpm. Particularly, calculations with multiple frames of reference method between the rotating and stationary parts of the domain are carried out. With the help of numerical results, correlation formula between the casing pressure and the efficiency is derived.

There are multistage preheaters in the power generation plan to improve the thermal efficiency of the plant and to prevent the components from the thermal shock. The energy source of these heaters comes from the extracted two phase fluid of working system. These two-phase fluid can cause the so-called Flow Accelerated Corrosion(FAC) in the extracting piping and the bubble plate of the heater for example, in case of point Beach Nuclear Power Plant and in the Wolsung Nuclear Power Plant. The FAC is due to the mass transport of the thin oxide layer by the convection. FAC is dependent on many parameters such as the operation temperature, void fraction, the fluid velocity and pH of fluid and so on. Therefore, in this paper velocity was calculated by FLUENT code in order to find out the root cause of the wall thinning of the feedwater heaters. It also includeed in the fluid mixing analysis model are around the number 5A feedwater heater shell including the extraction pipeline. To identify the relation between the local velocities and wall thinning, the local velocities according to the analysis results were compared with distribution of the shell wall thicknes by ultrasonic test.

The effects of nozzle-lip thickness on the relationship between screech tone and broadband shock-associated noise were experimentally investigated using a convergent-divergent nozzle with a design Mach number of 2.0. Overall sound pressure levels (OASPL) and noise spectra were obtained at far-field locations. Schlieren optical system was used to visualize the flow-fields of supersonic jets. A baffle plate was installed at the exit of the nozzle and its size was varied to obtain different nozzle-lip thicknesses. Experiment was carried out over a wide range of nozzle pressure ratios from 2.0 and 18.0, which corresponds to over- and under-expanded conditions. The results obtained clearly show that the screech tones are influenced by the nozzle-lip thickness. It is found that the screech tone and its peak amplitude are strongly dependent on whether the jet is over-expanded and under-expanded at the nozzle exit.

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 a conical convergent nozzle. The present study describes a computational work to investigate the effects of annular slit on the spiral jet. In the present computation, a finite volume scheme is used to solve three dimensional Naver-Stokes equations with RNG ${\kappa}-{\varepsilon}$ turbulent model. The annular slit width and the pressure ratio of the spiral jet are varied to obtain different spiral flows inside the conical convergent nozzle. The present computational results are compared with the previous experimental data. The results obtained obviously show that the annular slit width and the pressure ratio of the spiral jet strongly influence the characteristics of the spiral jets, such as tangential and axial velocities.

A cone cylinder is used to obtain variable operation conditions for the sonic ejector-diffuser system. The cone cylinder is designed to move upstream and downstream to change the ejector throat 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 for the variable sonic ejector system. In experiment, the ejector throat area is varied in the range from ${\psi}=11.88$ to 66.69, and the operating pressure ratio from $p_{0p}/p_a=1.25$ to 9.0. The results show that the variable sonic ejector system is suitable for a required entrainment ratio of secondary stream by altering the ejector throat area ratio and operating pressure ratio.

Recently series of steam explosion experiments have been performed in the TROI facility to identify the influence of corium compositions on the occurrence of a spontaneous steam explosion varying corium melt composition. The compositions of the corium were 0 : 100, 50 : 50, 70 : 30, 80 : 20 and 87 : 13 at weight percent of $UO_2$ to $ZrO_2$, and the mass of the corium was about 10kg. Corium melt at 0 : 100 weight percent (pure zirconia) caused a strong spontaneous steam explosion, and melt at 70 : 30 weight percent(eutectic corium) led to a weak steam spike, while melts at other compositions did not result in spontaneous steam explosions, when they came into contact with 67cm deep water pool at room temperature. It seems that the explosivity of pure zirconia is stronger than that of corium at other compositions and a steam explosion is not likely to occur with corium melts at non-eutectic compositions which are included in mushy zone region.

Inconel-600 alloy has been used as steam generator tube material for current pressurized water reactors (PWRs). The long-term operation of steam generators showed that the use of this material induced localized corrosion damages and increased tube wear of steam generator. To protect these problems, steam generator tube material is being changed to Inconel-690 alloy. Based on the current trend, we have chosen Inconel 690 as the Advanced Power Reactor 1400 (APR1400) steam generator(SG) tube material and performed the design optimization of preventive measure against tube fretting wear for the APR1400 steam generator. In this paper, we examined the technical consideration in this modification : the selection of material, wear characteristics, effect of the Egg-crate Flow Distribution Plate installation, and effect analysis of vertical strip installation.

The pressurized light water cooled, medium power (330 MWt) SMART (System-integrated Modular Advanced ReacTor) has been under development at KAERI for a dual purpose : seawater desalination and electricity generation. The SMART design verification phase was followed to conduct various separate effects tests and comprehensive integral effect tests. The high temperature / high pressure thermal-hydraulic test facility, VISTA(Experimental Verification by Integral Simulation of Transient and Accidents) has been constructed to simulate the SMART-P (the one fifth scaled pilot plant) by KAERI. Experimental tests have been performed to investigate the thermal-hydraulic dynamic characteristics of the primary and the secondary systems. Heat transfer characteristics and natural circulation performance of the PRHRS (Passive Residual Heat Removal System) of SMART-P were also investigated using the VISTA facility. The coolant flows steadily in the natural circulation loop which is composed of the steam generator (SG) primary side, the secondary system, and the PRHRS. The heat transfers through the PRHRS heat exchanger and ECT are sufficient enough to enable the natural circulation of the coolant.

Most power plants have operated many independent computerize systems for maintenance. Independence of systems have caused complexity of business process and inconvenience of computer system management. Because the equipment and material master data is not standardize and structurize, it is difficult to manage equipment maintenance history and material delivery. Especially equipment classification criterion is important for standardization of every maintenance information. It is necessary to integrate function of independent systems for business process simplification and rapid work flow. this paper provides equipment classification criterion design and system integration method with the case of live system development.

The main reason to analyze heat transfer in boiler inside through the performance test in fossil power plant is to increase plant high efficiency and energy saving movement in the government. Tins study intends to have trend and analyze the boiler heat transfer through the performance test, so it may give us the heat distribution in boiler inside in super-critical and sub-critica1 pressure type power plant

Performance of a solid oxide fuel cell (SOFC) can be enhanced by converting thermal energy of its high temperature exhaust gas to mechanical power using a micro gas turbine (MGT). A MGT plays also an important role to pressurize and warm up inlet gas streams of the SOFC. In this study, the influence of performance characteristics of the tubular SOFC on the hybrid power system is discussed. For this purpose, detailed heat and mass transfer with reforming and electrochemical reactions in the SOFC are mathematically modeled, and their results are reflected to the performance analysis. The analysis target is 220kWe SOFC/MGT hybrid system based on the tubular SOFC developed by Siemens-Westinghouse. Special attention is paid to the ohmic losses in the tubular SOFC counting not only current flow in radial direction, but also current flow in circumferential direction through the anode and cathode.

Future exhaust gas limits for diesel-driven passenger cars will force the automotive industry to significantly improve the performance of engine. Since modern common-rail injection systems deliver more degrees of freedom referring to the injection process, again the optimization of the injection process could offer a possibility to meet the exhaust gas limits. This study describes the characteristic the pilot spray structure of piezo-driven injector for a passenger car common-rail system to be applicable multiple injection caused by fast response rather than solenoid-driven injector. The piezo-driven injector is prototype injector with same needle chamber of solenoid injector and the solenoid-driven one is commercial injector. The pilot spray characteristic such as spray tip penetration, spray speed, spray angle were obtained by spray images, which is measured by the Mie scattering method with optical system for high-speed temporal photography. It was found that piezo-driven injector effected electric change as important factor and showed faster response than solenoid-driven injector.

In this study, the variation of spark plug location in the combustion chamber was investigated for the sake of emission characteristics from SI engine by using PDA valve. The swirl is ong of the important parameters that effects emission characteristics. PDA valve has been used to satisfy the requirements of sufficient swirl generation to improve combustion and emission reduction to effect on flow profile on a combustion chamber. Especially, the variation of spark plug location have an important effect to analyze exhaust gas and the early flame propagative process. Therefore, this test is forced that injection timing, spark timing and intake air motion govern the stable combustion. From the results, it showed that the variable spark plug location and PDA valve can be reduced exhaust gas.

Harmful elements from the exhaust gases are caused by incomplete combustion of mixture inside the engine cylinder and this abnormal combustion like misfire or partial burning is the direct cause of the air pollution and engine performance degradation. In this study, I obtain the shapes of spark, voltage and current generated when changing the experimental parameters such as grounded electrode shapes, electrode gap and the material of center electrodes. After that, I produce ignition energy by using the voltage and current and classify ignition energy into capacitive discharge energy and inductive discharge energy.

An experimental study on heat transfer characteristics near the critical pressure has been performed with an internally-heated vertical annular channel cooled by R-134a fluid. Two series of tests have been completed: (a) steady-state critical heat flux (CHF) and (b) heat transfer tests for pressure reduction transients through the critical pressure. In the present experimental range, the steady-state CHF decreases with the increase of the system pressure For a fixed inlet mass flux and subcooling, the CHF falls sharply at about 3.8 MPa and shows a trend toward converging to zero as the pressure approaches the critical point of 4.059 MPa. The CHF phenomenon near the critical pressure does not lead to an abrupt temperature rise of the heated wall because the CHF occurred at remarkably low power levels. In the pressure reduction transient experiments, as soon as the pressure passed through the critical pressure, the wall temperatures rise rapidly up to a very high value due to the occurrence of the departure from nucleate boiling. The wall temperature reaches a maximum at the saturation point of the outlet temperature, then tends to decrease gradually.

Amount of Electricity which product generator decide control valve at Turbine. Operating method of Control valve have two mode. First operating method is Partial Arc Admission, and second operating method is Full Arc Admission. Failure of Control Valve have on serious damage electricity lineage. In this Paper, We have investigated resonance that Control Valve spring casing.

Plate-fin type recuperators for the gas turbine/fuel cell hybrid power generating system were designed using commercial design software, MUSE. Heat transfer efficiency and total pressure drop in the recuperator were calculated to confirm required recuperator performance. Both counter flow and cross flow type plate-fin recuperators were designed. Results show that the counter flow type has higher efficiency and short core length, but the cross flow type is simpler to construct because the cross flow type does not need additional distributors. Two or three headers for the each recuperator core will be designed and tested to evaluate best header design. The designed recuperators and headers which will be designed later will be constructed, tested, and used in gas turbine/fuel cell hybrid power generating system.

This paper shows how to improve the efficiency and output and to reduce NOx emission of Seoinchon GE 7F gas turbine, Korea Western Power Co. by replacing the existing 7F gas turbine with new 7FA+e gas turbine because the performance of 7F gas turbine was degraded due to long term operation. In this paper, we will study gas turbine development trend and O&M technology. Finally, we will review for uprate of Seoinchon 7F gas turbine to help someone to improve their units in the future.

The main objective of the present study is to investigate the performance characteristics of a ground source heat pump (GSHP) system with a 130 m vertical 60.5 mm nominal diameter U-bend ground heat exchanger. In order to evaluate the performance analysis, the GSHP system connected to a test room with 90 $m^2$ floor area in the Korea Institute of Construction Technology ($37^{\circ}39'$ N, $126^{\circ}48'$ E) was designed and constructed. This GSHP system mainly consisted of ground heat exchanger, indoor heat pump and measuring devices. The cooling and heating loads of the test room were 5.5 and 7.2 kW at design conditions, respectively. The experimental results were obtained from July to January in cooling and heating season of $2003{\sim}2004$. The cooling and heating performance coefficients of the system were determined from the experimental results. The average cooling and heating COPs for the system were obtained to be 4.82 and 3.02, respectively. The temperature variations in ground and the ground heat exchanger surface at different depths were also measured.

In this study, analysis is made for the effects of groove shape on the thermal performance of a axial groove heat pipe. The mathematical models of two-phase flow in grooved heat pipe are presented for the capillary limitation in steady state. Generally, the heat pipe performance depends on the capillary pressure and liquid flow. The friction force of liquid flow through the groove increases with the groove width decreased, and then the capillary pressure is improved in the gas-liquid interface of groove. Therefore, the optimal groove width shaper exists for the maximum thermal performance of heat pipe. In this paper, the optimal groove shape and scale are presented by considering both capillary pressure and liquid flow.

The present study has been conducted to develop a heat pump system using river water of temperature energy which not only belongs to unutilized energy but is a kind of good heat source due to maintain its temperature in a certain degree regardless of seasonal variation. The system did not meet the proposed performance after setup. In this paper, the system performance affected by refrigerant Oil, by pressure drop, or by other factors has been discussed. The followings were obtained : (1) Refrigerant Oil mixture rate was 2.5 in weight percentage, (2) Pressure drop through evaporator was 29.1kPa($3.1^{\circ}C$ in saturated tempearture) (3) Pressure drop from the end of evaporator to compressor inlet was 39.8kPa($4.0^{\circ}C$ in saturated tempearture). (4) The system performance can to be improved by modifying a part of pipe line to compressor, and reducing pressure drop through heat exchangers.

This study represented experimental research on the flat plate solar collector. For the flat plate Solar system, it is sensitive of the Global radiation. In Actually, it suppose to be dependent on the direct radiation. Also, the existing method's factors are depend upon Global radiation in the flat plate collector system. therefore it needs which is depend upon direct radiation. In this experiment, the flat plate collector is used for obtaining the method's factors of the direct radiation. As a result, the correct $({\tau}{\alpha})_e$ is found out for practical value.

As part of study on thermal hydraulic behavior in the reactor cavity under external vessel cooling in the APR (Advanced Power Reactor) 1400, one dimensional two phase flow of steady state in the reactor cavity have been analyzed to investigate a coolant circulation mass flow rate in the annulus region between the reactor vessel and the insulation material using the RELAP5/MOD3 computer code. The RELAP5/MOD3 results have shown that a two phase natural circulation flow of 300 - 600 kg/s is generated in the annulus region between the reactor vessel and the insulation material when the external vessel cooling has been applied in the APR 1400. An increase in the heat flux of the inner vessel leads to an increase of the coolant mass flow rate. An increase in the coolant outlet area leads to an increase in the coolant circulation mass flow rate, but the coolant inlet area does not effective on the coolant circulation mass flow rate. The change of the lower coolant outlet to a lower position affects the coolant circulation mass flow rate, but the variation trend is not consistent.

To evaluate the variation of spring stiffness in nuclear plant operating condition, load-displacement tests ($P-{\delta}$ test) were performed using two kinds of space grid springs in high temperature and high pressure water. With increasing temperature, stiffness of each spring gradually decreased except $100{\sim}150^{\circ}C$. It is apparently showed that spring with convex shape had a relatively high stability of spring stiffness at high temperature compared with I-shaped spring. It is suggested that the variation of spring stiffness with temperature and spring shape should be considered as an important variable in the design and analysis of the fuel assembly.

In this study, five experiments were carried out, with an orifice located downstream of a fan in case I and upstream distance to the fan in the rest cases(case $II{\sim}$), so as to determine the optimal location of the orifice and reduce the size of airflow measurement device. The resulting flow rate-pressure drop correlations were found to satisfy the limitation of SMACNA standard, which specified an error of ${\pm}7.5%$ based on the real flow rate. The best outcome was achieved with the orifice located midway of the orifice duct 4 times its diameter long.

As the friction between two mated materials generates electrostatic voltage, vehicles especially in the cold climate has been accumulated an electrostatic voltage by the reason of fuel flow in the fuel line, air flows out side of vehicle, etc. In this study, investigations of electrostatic discharge characteristics has been carried out in the fuel systems of gasoline engines depending on the environment temperature. For the experiment, conventional fuel filters(paint coated steel case) and specially made testing equipment are prepared. The experimental results shows, an accumulation speed of electrostatic voltage depends on the environmental temperature and new methods of testing procedures are required compare to SAE standards.

For designing high Tc superconductor proceeding bearing(HTSJB) which is used on a flywheel energy storage system which requests the free of the bearing loss caused by the friction, it is necessary to understand the basic characteristics of the classical superconductor proceeding bearing because the mechanical characteristics of the HTSJB are identified by the magnetic relationships between the permanent magnet(PM) and the high Tc superconductor(HTS). In this paper, using the method, frozen image model, the force problems between the PM and the HTS were solved and then the dynamic characteristics of the rotor inside of the HTSJB can be expected in advance by using the basic characteristics between the PM and the HTS. The coefficient of friction of the HTSJB was measured in the vacuum environment. From the results, the mechanical characteristics of HTSJB can be identified using the numerical models.

Current common rail engines are equipped with cooled EGR systems by using an engine cooling water system. In this study, investigations of exhaust gas reduction characteristics have been carried out in the common rail engine system depending on the EGR rate variation. The experimental results shows that NOx reduces and smoke increases as the EGR rate increases.