• 설비공학논문집
• /
• 제1권2호
• /
• pp.128-137
• /
• 1989

Numerical solutions for two-dimensional, steady, free convection are presented for a cylinder filled with saturated porous media. An annulus is bounded by inner wall with constant heat flux and two adiabatic horizontal walls with outer wall isothermally cooled. Governing equations are numerically solved for the range of Aspect Ratio 1 to 20, Radius Ratio, 1 to 20, and Rayleigh number, 50 to $10^4$ by Finite Difference method utilizing upwind scheme. Results are presented in terms of stream lines and isotherms, temperature distributions and local Nusselt numbers at the heated wall. Average Nusselt numbers are also presented for the comparisons.

• 대한기계학회논문집B
• /
• 제20권1호
• /
• pp.266-274
• /
• 1996

An experimental work was conducted to investigate the characteristics of the onset of flooding in vertical annuli with a direct injection mode using air and water. The onset of flooding was determined by means of pressure drop measurement while the air velocity was increased gradually under fixed liquid flow rates. Data of the onset of flooding were collected for various combinations of the tube size and the nozzle number. A theoretical analysis of the onset of flooding was also performed based on an envelope theory. The result shows that the onset of flooding in small-scale annuli can be predicted relatively well by the theory. A modified Wallis parameter was used to investigate the scaling effect of flooding phenomena in the annuli, indicating a relatively reasonable result. The number of nozzle has no effect on the flooding velocity when liquid was injected through 2, 3, 4 and 6 nozzles but the initiation of flooding was significantly expedited when 12 nozzles were employed for liquid injection.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2007년도 춘계학술대회B
• /
• pp.2696-2701
• /
• 2007

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

• Journal of Advanced Marine Engineering and Technology
• /
• 제10권2호
• /
• pp.156-164
• /
• 1986

In the case of boiling on high temperature wall, vapor film covers fully or parcially the surface. This phenomenon, film boiling or transition boiling, is very important in the surface heat treatment of metal, design of cryogenic heat exchanger and emergency cooling of nuclear reactor. Mainly supposed hydraulic-thermal accidents in nuclear reactor are LCCA (Loss of Coolant Accident) and PCM (Power-Cooling Mismatch). Recently, world-wide studies on reflooding of high temperature rod bundles after the occurrence of the above accidents focus attention on wall temperature history and required time in transient cooling process, wall superheat at rewet point, heat flux-wall superheat relationship beyond the transition boiling region, and two-phase flow state near the surface. It is considered that the further systematical study in this field will be in need in spite of the previous results in ref. (2), (3), (4). The paper is the study about the fast transient cooling process following the wall temperature excursion under the CHF (Critical Heat Flux) condition in a forced convective subcooled boiling system. The test section is a vertically arranged concentric annulus of 800 mm long and 10 mm hydraulic diameter. The inner tube, SUS 304 of 400 mm long, 8 mm I.D, and 7 mm O.D., is heated uniformly by the low voltage AC power. The wall temperature measurements were performed at the axial distance from the inlet of the heating tube, z=390 mm. 6 chromel- alumel thermocouples of 76 .mu.m were press fitted to the inner surface of the heating tube periphery. To investigate the heat transfer characteristics during the fast transient cooling process, the outer surface (fluid side) temperature and the surface heat flux are computed from the measured inner surface temperature history by means of a numerical method for inverse problems of transient heat conduction. Present cooling (boiling) curve is sufficiently compared with the previous results.

• 한국결정성장학회:학술대회논문집
• /
• 한국결정성장학회 1996년도 The 9th KACG Technical Annual Meeting and the 3rd Korea-Japan EMGS (Electronic Materials Growth Symposium)
• /
• pp.179-200
• /
• 1996

The intrinsic instabilities of fluid flow occurred in the melt of the Czochralski crystal growth system Czochralski method, asymmetric flow patterns and temperature profiles in the melt have been studied by many researchers. The idea that the non-symmetric structure of the growing equipment is responsible for the asymmetric profiles is usually accepted at the first time. However further researches revealed that some intrinsic instabilities not related to the non-symmetric equipment structure in the melt could also appear. Ristorcelli had pointed out that there are many possible causes of instabilities in the melt. The instabilities appears because of the coupling effects of fluid flow and temperature profiles in the melt. Among the instabilities, the B nard type instabilities with no or low crucible rotation rates are analyzed by the visualizing experiments using X-ray radiography and the 3-D numerical simulation in this study. The velocity profiles in the Silicon melt at different crucible rotation rates were measured using X-ray radiography method using tungsten tracers in the melt. The results showed that there exits two types of fluid flow mode. One is axisymmetric flow, the other is asymmetric flow. In the axisymmetric flow, the trajectory of the tracers show torus pattern. However, more exact measurement of the axisymmetrc case shows that this flow field has small non-axisymmetric components of the velocity. When fluid flow is asymmetric, the tracers show random motion from the fixed view point. On the other hand, when the observer rotates to the same velocity of the crucible, the trajectory of the tracer show a rotating motion, the center of the motion is not same the center of the melt. The temperature of a point in the melt were measured using thermocouples with different rotating rates. Measured temperatures oscillated. Such kind of oscillations are also measured by the other researchers. The behavior of temperature oscillations were quite different between at low rotations and at high rotations. Above experimental results means that the fluid flow and temperature profiles in the melt is not symmetric, and then the mode of the asymmetric is changed when rotation rates are changed. To compare with these experimental results, the fluid flow and temperature profiles at no rotation and 8 rpm of crucible rotation rates on the same size of crucible is calculated using a 3-dimensional numerical simulation. A finite different method is adopted for this simulation. 50×30×30 grids are used. The numerical simulation also showed that the velocity and flow profiles are changed when rotation rates change. Futhermore, the flow patterns and temperature profiles of both cases are not axisymmetric even though axisymmetric boundary conditions are used. Several cells appear at no rotation. The cells are formed by the unstable vertical temperature profiles (upper region is colder than lower part) beneath the free surface of the melt. When the temperature profile is combined with density difference (Rayleigh-B nard instability) or surface tension difference (Marangoni-B nard instability) on temperature, cell structures are naturally formed. Both sources of instabilities are coupled to the cell structures in the melt of the Czochralski process. With high rotation rates, the shape of the fluid field is changed to another type of asymmetric profile. Because of the velocity profile, isothermal lines on the plane vertical to the centerline change to elliptic. When the velocity profiles are plotted at the rotating view point, two vortices appear at the both sides of centerline. These vortices seem to be the main reason of the tracer behavior shown in the asymmetric velocity experiment. This profile is quite similar to the profiles created by the baroclinic instability on the rotating annulus. The temperature profiles obtained from the numerical calculations and Fourier transforms of it are quite similar to the results of the experiment. bove esults intend that at least two types of intrinsic instabilities can occur in the melt of Czochralski growing systems. Because the instabilities cause temperature fluctuations in the melt and near the crystal-melt interface, some defects may be generated by them. When the crucible size becomes large, the intensity of the instabilities should increase. Therefore, to produce large single crystals with good quality, the behavior of the intrinsic instabilities in the melt as well as the effects of the instabilities on the defects in the ingot should be studied. As one of the cause of the defects in the large diameter Silicon single crystal grown by the