• 에너지공학
• /
• 제23권4호
• /
• pp.95-105
• /
• 2014

CUPID 코드는 기기 스케일(Component scale)의 2상 유동(Two-phase flow) 해석 코드로서 다양한 2상 유동 조건의 실험 자료를 이용하여 검증되어 왔다. 특히, CUPID 코드의 CANDU형 원자로 감속재 탱크 내부 유동 해석능력을 평가하기 위해 1/4 규모 축소 실험장치의 실험결과를 이용하여 검증한 바가 있다. 본 연구에서는 이전 연구를 바탕으로 CUPID 코드를 사용하여 실제 원자로 감속재탱크 내부의 열수력 거동을 해석하였다. 감속재 탱크의 내부 구조는 아주 복잡하기 때문에 다공질 매질 방법을 적용하였으며 탱크 입구노즐 또한 기기 스케일 코드의 취지에 부합하게 아주 단순화하여 모델하였다. 해석결과의 정확성을 결정하는 가장 중요한 요소는 입구노즐의 모델 방법에 있는 것으로 나타났다. 입구노즐을 단순하게 모델하여 입구유량을 경계조건으로 부여하고 발전소 정상운전조건으로 계산한 결과, 부력에 의한 열성층화 현상이 발생하였다. 이는 전혀 타당하지 않은 것으로 입구 유동의 모멘텀을 정확하게 모의하지 않아 발생한 것이 나타났다. 이를 개선하고자 입구 유량과 운동량을 동시에 보존시킬 수 있도록 입구 노즐 면적을 축소하고 속도는 증가시켜서 계산한 결과, 사실적인 내부 유동장을 얻을 수 있었다. 결론적으로 계산 비용효과가 뛰어난 다공질 매질 방법에 입각하여 CUPID 코드를 실규모 감속재 탱크 열유동 해석에 적용할 수 있음을 보였고, 입구노즐의 적절한 모델이 가장 중요한 요소임을 확인하였다.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1996년도 춘계학술발표회논문집(2)
• /
• pp.233-240
• /
• 1996

Analyses of the viscous and end effects on electromagnetic (EM) pumps of annular linear induction type for the sodium coolant circulation in Liquid Metal Fast Breeder Reactors have been carried out based on the MHD laminar flow analysis and the electromagnetic field theory. A one-dimensional MHD analysis for the liquid metal flowing through an annular channel has been performed on the basis of a simplified model of equivalent current sheets instead of three-phase currents in the discrete primary windings. The calculations show that the developed pressure difference resulted from electromagnetic and viscous forces in the liquid metal is expressed in terms of the slip, and that the viscous loss effects are negligible compared with electromagnetic driving forces except in the low-slip region where the pumps operate with very high flow velocities comparable with the synchronous velocity of the electromagnetic fields, which is not applicable to the practical EM pumps. A two-dimensional electromagnetic field analysis based on an equivalent current sheet model has found the vector potentials in closed form by means of the Fourier transform method. The resultant magnetic fields and driving forces exerted on the liquid metal reveal that the end effects due to finiteness of the pump length are formidable. In addition, a two-dimensional numerical analysis for vector potentials has been performed by the SOR iterative method on a realistic EM pump model with discretely-distributed currents in the primary windings. The numerical computations for the distributions of magnetic fields and developed pressure differences along the pump axial length also show considerable end effects at both inlet and outlet ends, especially at high flow velocities. Calculations of each magnetic force contribution indicate that the end effects are originated from the magnetic force caused by the induced current ( u x B ) generated by the liquid metal movement across the magnetic field rather than the one (E) produced by externally applied magnetic fields by three-phase winding currents. It is concluded that since the influences of the end effects in addition to viscous losses are extensive particularly in high-velocity operations of the EM pumps, it is necessary to find ways to suppress them, such as proper selection of the pump parameters and compensation of the end effects.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.320-325
• /
• 2010

A multiphase CFD analysis is performed to investigate the effect of near-wall grid for simulating a subcooled boiling flow in vertical tube. The multiphase flow model used in this CFD analysis is the two-fluid model in which liquid(water) and vapor(steam) are considered as continuous and dispersed fluids, respectively. A wall boiling model is also used to simulate the subcooled boiling heat transfer at the heated wall boundary. The diameter and heated length of tube are 0.0154 m and 2 m, respectively. The system pressure in tube is 4.5 MPa and the inlet subcooling is 60 K. The near-wall grid size in the non-dimensional wall unit ($y_{w}^{+}$) was examined from 64 to 172 at the outlet boundary. The CFD calculations predicted the void distributions as well as the liquid and wall temperatures in tube. The predicted axial variations of the void fraction and the wall temperature are compared with the measured ones. The CFD prediction of the wall temperature is shown to slightly depend on the near-wall grid size but the axial void prediction has somewhat large dependency. The CFD prediction was found to show a better agreement with the measured one for the large near-wall grid, e.g., $y_{w}^{+}$ > 100.

• 대한기계학회논문집B
• /
• 제29권11호
• /
• pp.1173-1181
• /
• 2005

The cooling effect of a mobile phone using PCM(Phase Change Material) module has been numerically investigated. A transient three-dimensional numerical analysis of heat and fluid flow with natural convection is performed in this study. Governing conservation equations for mass, momentum and energy are solved by an implicit finite volume method. An enthalpy-porosity technique has been used for modeling of the melting process. Two different ways of placing the PCM module are considered. One is to place a PCM module between the substrate and battery pack, and the other is to place a PCM module between MCM(multichip module) and battery pack. Three different types of PCMs are used to predict the performance of PCM. The results show that passive cooling with PCM can reduce the temperature rise and the effect of natural convection in PCM module considered in this study is negligible.

• 한국초전도ㆍ저온공학회논문지
• /
• 제15권3호
• /
• pp.49-52
• /
• 2013

An experiment to investigate the natural circulation of a cryogen has been performed. The study is motivated mainly by our recent development of cryogenic cooling system for prototype superconducting cyclotron without any circulating pump. In the natural circulation loop system, a cooling channel is attached on the outer surface of the aluminium block and the liquid nitrogen passes through inside of the channel to cool the block indirectly. A cryocooler as a heat sink is located at the top to re-condense cryogenic vapor coming from the aluminium block in which electrical heater is installed as a heat source. The main dimensions are determined using the relevant analysis and the natural circulation loop is successfully fabricated. The temperature distributions in the loop are measured during initial cool-down process and in steady state, from which the modified Grashof numbers are calculated and compared with the existing correlation estimated with one-dimensional analysis for steady state flow.

• 한국자동차공학회논문집
• /
• 제18권5호
• /
• pp.108-114
• /
• 2010

The injection nozzle of an electro-hydraulic injector is being opened and closed by movement of a injector's needle which is balanced by pressure at the nozzle seat and at the needle control chamber, at the opposite end of the needle. In this study, the effects of needle movement in a piezo-driven injector on unsteady cavitating flows behavior inside nozzle were investigated by cavitation numerical model based on the Eulerian-Lagrangian approach. Aimed at simulating the 3-D two-phase flow behavior, the three dimensional geometry model along the central cross-section regarding of one injection hole with real design data of a piezo-driven diesel injector has been used to simulate the cavitating flows for injection time by at fully transient simulation with cavitation model. The cavitation model incorporates many of the fundamental physical processes assumed to take place in cavitating flows. The simulations performed were both fully transient and 'pseudo' steady state, even if under steady state boundary conditions. As this research results, we found that it could analyze the effect the pressure drop to the sudden acceleration of fuel, which is due to the fastest response of needle, on the degree of cavitation existed in piezo-driven injector nozzle.

• 한국연소학회:학술대회논문집
• /
• 한국연소학회 2004년도 제29회 KOSCI SYMPOSIUM 논문집
• /
• pp.17-23
• /
• 2004

In the present work, the numerical model was refined to predict the thermal analysis of energy storage in a fixed beds during (charging ,storing, discharging) mode. The governing energy equations of both fluid and the solid particles along with their initial and boundary conditions are derived using a two-phase, one dimensional model. The refined model is carried out by taking into account change of (air density , air specific heat) with air temperature and also by taking into considerations heat losses from bed to surrounding. Finite difference method was used to obtain solution of two governing energy equations of both fluid and solid particles through a computer program especially constructed for this purpose. The temperature field for the air and the solid are obtained, also efficiency of energy stored inside the bed is computed. Finally using refined model the effect of air flow rate per unit area Ga (0.2, 0.3, and 0.4 kg/$m^2$-s), and inlet air temperature (200, 250, 300 $^{\circ}C$) on energy storage characteristics was studied in three mode ( charging ,storing, discharging). The rock particles of diameter 1 em is used as bed material in this research.

• 한국연소학회:학술대회논문집
• /
• 한국연소학회 2004년도 제28회 KOSCO SYMPOSIUM 논문집
• /
• pp.158-164
• /
• 2004

In the present work, the numerical model was refined to predict the thermal analysis of energy storage in a fixed beds during charging mode. The governing energy equations of both fluid and the solid particles along with their initial and boundary conditions are derived using a two-phase, one dimensional model. The refined model is carried out by taking into account change of (air density , air specific heat) with air temperature and also by taking into considerations heat losses from bed to surrounding. Finite difference method was used to obtain solution of two governing energy equations of both fluid and solid particles through a computer program especially constructed for this purpose. The temperature field for the air and the solid are obtained, also energy stored inside the bed is computed. A comparison between refined model and non refined model is done. Finally using refined model the effect of bed material (Glass, Fine clay ,and aluminum ), and air flow rate per unit area Ga (0.3, 0.4, and 0.5 kg/$m^2$-s) on energy storage characteristics was studied.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
• /
• pp.383-389
• /
• 2003

Clay minerals show a distinct induced-polarization phenomenon, which is one of the most important factors for predicting groundwater flow and contaminant transport. This paper presents a step-by-step process to estimate Cole-Cole parameters from spectral induced-polarization (IP) data measured on the surface of three-dimensional earth. First, the inversion of low-frequency resistivity survey data is made to identify the dc resistivity ${\rho}_dc$ of a volume having IP effects. The other parameters, chargeability m, time constant $\tau$, and frequency dependence c, are sought for the polarizable volume. Next, using multi-frequency data, c can be obtained as high or low asymptotes of the slope of log phase vs. log frequency. Further, for low m, intrinsic $\tau$ is approximated by apparent one, ${\tau}_a$, which is derived from the relation ${{\omega}{\tau}}_a$=1 at an angular frequency $\omega$, where the imaginary component of spectral IP data has an extreme value. Finally, to obtain intrinsic m a two-step linearized procedure has been derived. For a body of given $\tau$ and c, forward modeling with a progression of m values yields a plot of observed vs. intrinsic imaginary components for a frequency. Since this plot is essentially linear, to extract the intrinsic imaginary component is quite simple with an observed value. Using the plot of intrinsic imaginary component vs. m, intrinsic m is determined. We present a synthetic example to illustrate that the Cole-Cole parameters can be recovered from spectral IP data.

• 대한기계학회논문집
• /
• 제16권6호
• /
• pp.1205-1215
• /
• 1992

본 연구에서는 1차원 상변화문제를 확장하여 사각형 주형내에서의 접촉열저항 을 고려한 2차원 상변화문제를 해석하고자 한다.즉 고상/액상 경계면인 상변화면을 추적하기 위하여 온도수식법(temperature formulation method)을 적용시 나타나는 복 잡성을 극복하기 위하여 제시된 엔탈피법을 이용하여 순수물질의 경우 뿐만 아니라 2상영역이 존재하는 합금에서의 온도분포, 상변화면의 이동속도 및 총응고에 소요되는 시간등의 열적특성을 해석하며, 또한 이러한 열적특성에 영향을 미치는 매개변수를 도출하고자 주형의 형상, 주형 외벽에서의 열전달계수, 접촉 열전달계수, 주형의 열전 도계수 및 주물의 초기온도등의 변화에 따른 효과를 고찰한다.