• Journal of computational fluids engineering
• /
• v.21 no.2
• /
• pp.54-61
• /
• 2016

When spent fuel assemblies from the reactor of nuclear power plants(NPPs) are transported, the assemblies are exposed to short-term operations that can affect the peak cladding temperature of spent fuel assemblies. Therefore, it needs to perform the analysis of heat transfer on spent fuel dry cask during the operation. For 3 dimensional computational fluid dynamnics(CFD) simulation, it is proposed that the short-term operation is divided into three processes: Wet, dry, and vacuum drying condition. The three processes have different heat transfer mode and medium. Metal transportation cask, which is Korea Radioactive Waste Agency(KORAD)'s developing cask, is evaluated by the methods proposed in this work. During working hours, the boiling at wet process does not occur in the cask and the peak cladding temperatures of all processes remain below $400^{\circ}C$. The maximum peak cladding temperature is $173.8^{\circ}C$ at vacuum drying process and the temperature rise of dry, and vacuum drying process occurs steeply.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.2
• /
• pp.193-203
• /
• 2014

Temperature-dependent stress analysis and heat transfer analysis of a rotating gas turbine blade made of functionally graded materials (FGMs) are presented considering turbine operating temperature and ceramic particle size. The material properties of functionally graded materials are assumed to vary continuously and smoothly across the thickness of the thin-walled blade. For obtaining system stiffness reflecting these characteristics, the one-dimensional heat transfer equation is applied along the thickness of the thin-walled blade for determining the temperature distribution. Using the results of the temperature analysis, the equations of motion of a rotating blade are derived with hybrid deformation variable modeling method along with the Rayleigh-Ritz assumed mode methods. The validity of the derived rotating blade model is evaluated by comparing its transient responses and temperature distribution with the results obtained using a commercial finite element code. The maximum tensile stress with operating speed and gradient index are obtained. Furthermore, the gradient index that minimizes blade temperature was investigated.

• 한국연소학회:학술대회논문집
• /
• 2004.11a
• /
• pp.197-205
• /
• 2004

Most of gas turbines is operated by the type of dry premixed combustion to reduce NOx emission and economize fuel consumption. However this type operation, combustion induced instability brought failure problems cause by high pressure and heat release fluctuations. Though there has been lots of studies since Lord Rayleigh to understand this instability mechanism and control the instabilities, none of them made matters clear. In order to understand the instability phenomena, a simple experimental study with dump combustor was conducted at the moderate pressure and ambient temperature conditions. From this model gas turbine combustor self-excited instabilities at the resonance mode(200Hz) and bulk mode(10Hz) were occurred and observed at the three points of view; pressure, heat release and equivalence ratio which are acquired by peizo-electric transducer, HICCD camera and acetone LIF respectively. From this results we could see the instability mechanism clear with the account of time scale analysis which explained by the propagation of pressure wave to the upward of mixture stream and convectional transfer of the equivalence ratio fluctuation by this pressure fluctuation.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1997.04a
• /
• pp.552-556
• /
• 1997

The Direct Metal Prototyping(DMP), one of the rapid prototyping technologies, allows the manufacturing of three-dimensional metallic parts using metal powders directly from the CAD data. Laser power and scanning speed are the most important variables of the process. The objective of this study is to obtain the design data for laser power and scanning speed to bond metal powders effectively using the finite element method. To obtain the design values, a numerical analysis considering two-dimensional heat transfer during the sintering of metal powder layers of the process was performed. The laser beam has been modeled to have directionality in its heat flux distribution, i. e., in the scanning direction a Gaussian beam mode distribution has been assumed and in the thickness direction a square beam mode distribution. The three-dimensional irregular distribution of metal powders of the powder layer is idealized as two-dimensional distribution in which metal powders are located regularly and periodically on the plate. In this study the design values of laser power vs scanning speed have been obtained. Temperature distribution and temperature variation of the powder layers with respect to time have been predicted. The commputed dsign data will be useful in determining the initial conditions of the process.

• Nuclear Engineering and Technology
• /
• v.44 no.3
• /
• pp.283-296
• /
• 2012

The natural convection in a horizontal enclosure heated from the bottom wall, cooled at the top wall, and having a square adiabatic body in the center is studied. Three different Prandtl numbers (0.01, 0.7 and 7) are considered for the investigation of the effect of the Prandtl number on natural convection. Adiabatic boundary conditions are employed for the side walls. A two-dimensional solution for unsteady natural convection is obtained, using an accurate and efficient Chebyshev spectral methodology for different Rayleigh numbers varying over the range of $10_3$ to $10_6$. It had been experimentally reported that the heat transfer mode becomes oscillatory when Pr is out of a specific Pr band beyond the critical Ra. In this study, we reproduced this phenomenon numerically. It was found that when Ra=$10_6$, only the case for intermediate Pr (=0.7) reached a non-changing steady state and the low and high Pr number cases (Pr=0.01 and 7) showed a periodically oscillatory fashion hydrodynamically and thermally. The variation of time- and surface-averaged Nusselt numbers on the hot and cold walls for different Rayleigh numbers and Prandtl numbers are presented to show the overall heat transfer characteristics in the system. Further, the isotherms and streamline distributions are presented in detail to compare the physics related to their thermal behavior.

• Transactions of the Korean hydrogen and new energy society
• /
• v.26 no.2
• /
• pp.184-191
• /
• 2015

As the power electronics system increases the frequency, the power loss and thermal management are paid more attention. This research presents a real time model of dissipation power with junction temperature response for 120kw IGBT inverter which is applied to the thermal management of high power IGBT inverter. Since the computational time is critical for real time simulation, look-up tables of IGBT module characteristic curve are implemented. The power loss from IGBT provides a clue to calculate the temperature of each module of IGBT. In this study, temperature of each layer in IGBT is predicted by lumped capacitance analysis of layers with convective heat transfer. The power loss and temperature of layers in IGBT is then communicated due to mutual dependence. In the dynamic model, PWM pulses are employed to calculation real time IGBT and diode power loss. Under Matlab/Simulink$^{(R)}$ environment, the dynamic model is validated with experiment. Results showed that the dynamic response of power loss is closely coupled with effective thermal management. The convective heat transfer is enough to achieve proper thermal management under guideline temperature.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.10
• /
• pp.2595-2606
• /
• 1995

Buoyancy-assisted melting of an unconstrained ice in an isothermally heated horizontal enclosure was numerically analyzed in a range of wall temperatures encompassing the density inversion point. The problem as posed here involves two physically distinct domains each of which has its own scales and respective heat transfer mode. These two domains join at the junction where the liquid squeezed out of the film region flushes into the lower melt pool. Both of these domains have been treated separately in the literature by a patching technique which invokes several, otherwise unnecessary, assumptions. The present study eliminates successfully such a superfluous procedure by treating the film and lower melt pool regions as a single domain. As a result of this efficient solution procedure, the interaction of the water stream ejected at the junction and the natural convection in the melt pool could be clarified for different wall temperatures. Though limited by two-dimensionality, the present results conformed indirectly the earlier reported transition of the flow pattern, as the wall temperature was increased over the density inversion point. The transient evolution of the melting surface, the time rate of change in melt volume fraction, the local and temporal variation of the heat transfer coefficients are analyzed and presented.

• Journal of computational fluids engineering
• /
• v.16 no.3
• /
• pp.29-36
• /
• 2011

The natural convection in a horizontal enclosure heated from the bottom wall, cooled at the top wall, and having a square adiabatic body at its centered area was studied. Three different Prandtl numbers (0.01, 0.7 and 7) were considered for an effect of the Prandtl number on natural convection. A two-dimensional solution for unsteady natural convection was obtained, using Chebyshev spectral methodology for different Rayleigh numbers varying over the range of $10^4$ to $10^6$. It had been experimentally and numerically reported [1,2] that the heat transfer mode becomes oscillatory when Pr is out of a specific Pr band beyond the critical Ra. In this study, we reproduced this phenomenon numerically. The variation of time- and surface-averaged Nusselt numbers on the hot and cold walls for different Rayleigh numbers and Prandtl numbers was presented to show the overall heat transfer characteristics in the system. And also, the isotherms and streamline distributions were presented in detail to compare the physics related to their thermal behavior.

• Progress in Superconductivity and Cryogenics
• /
• v.11 no.3
• /
• pp.40-45
• /
• 2009

Many superconductor applications such as MRI and SMES must be operated in persistent current mode to eliminate the electrical ohmic loss. This paper presents the characteristic analysis of the high temperature superconducting (HTS) power supply made of YBCO coated conductor (CC). In this research, we have manufactured the HTS power supply to charge the 0.73 mH HTS double-pancake magnet made of YBCO CC. Among the all design parameters, the heater triggerring time and magnet applying time were the most important factors for the best performance of the HTS power supply. In this paper, three-dimensional simulation through finite element method (FEM) was used to study the heat transfer in YBCO CC and the magnetic field of the magnetic circuit. Based upon these results, the final operational sequence could be determined to generate the pumping current. In the experiment, the maximum pumping current reached about 16 A.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.13 no.6
• /
• pp.142-153
• /
• 2005

The brake judder is a phenomenon that the steering wheel is abnormally vibrating when the car is braked at a high speed. It is classified by the cold and the hot judder. The former is generated due to the initial uneven disk surface and the latter is resulted from the uneven heat spots on disc surface by repeatedly braking. There are two ways to reduce the judder. One is to control vibration by modification of the disk shapes and pad ingredients. The other is to improve modal characteristics of the suspension system. The latter approach is used in this research. In this paper, the real vehicle test and computer simulation are considered to systematically understand the judder phenomenon of the vehicle. The Macpherson strut suspension is employed. Especially, the judder sensitivity is calculated based on design sensitivity analysis. A bush stiffness was reworked and braking test was done to verify the sensitivity result. The judder reduction by the mode control was verified.