• Nuclear Engineering and Technology
• /
• v.55 no.8
• /
• pp.3017-3029
• /
• 2023

In the late in-vessel phase of a nuclear reactor severe accident, the internal heat transfer and crust evolution during the debris bed melting process have important effects on the thermal load distribution along the vessel wall, and further affect the reactor pressure vessel (RPV) failure mode and the state of melt during leakage. This study coupled the phase change model and large eddy simulation to investigate the variations of the temperature, melt liquid fraction, crust and heat flux distributions during the debris bed melting process in the hypothetical severe accident of HPR1000. The results indicated that the heat flow towards the vessel wall and upper surface were similar at the beginning stage of debris melting, but the upward heat flow increased significantly as the development of the molten pool. The maximum heat flux towards the vessel wall reached 0.4 MW/m². The thickness of lower crust decreased as the debris melting. It was much thicker at the bottom region with the azimuthal angle below 20° and decreased rapidly at the azimuthal angle around 20-50°. The maximum and minimum thicknesses were 2 and 90 mm, respectively. By contrast, the distribution of upper crust was uniform and reached stable state much earlier than the lower crust, with the thickness of about 10 mm. Moreover, the sensitivity analysis of initial condition indicated that as the decrease of time interval from reactor scram to debris bed dried-out, the maximum debris temperature and melt fraction became larger, the lower crust thickness became thinner, but the upper crust had no significant change. The sensitivity analysis of in-vessel retention (IVR) strategies indicated that the passive and active external reactor vessel cooling (ERVC) had little effect on the internal heat transfer and crust evolution. In the case not considering the internal reactor vessel cooling (IRVC), the upper crust was not obvious.

• Journal of Advanced Marine Engineering and Technology
• /
• v.27 no.5
• /
• pp.617-623
• /
• 2003

In the present study, the characteristics of upward bubble flow were experimentally investigated in a liquid bath. The velocity of upward bubble flow was calculated for two different experimental conditions：1) bubble flow without kinetic energy 2) bubble flow with kinetic energy. Bubble flow without kinetic energy starts to undergo the effect of buoyancy l0cm away from the nozzle. Whereas. kinetic energy is dominant before 30 cm away from the nozzle in bubble flow but after this point kinetic energy and inertial force are applied on bubble flow at the same time In addition, as the flow rate increases the maximum velocity point moves to the nozzle. The velocity Profiles near free surface is extremely irregular due to surface flow. Gas volume fraction is high near the nozzle due to gas concentration. but decreases with the increasement of axial position. Gas volume fraction does not vary after the axial position, z=60 in spite of the increasement of flow.

• 유체기계공업학회:학술대회논문집
• /
• 1999.12a
• /
• pp.118-122
• /
• 1999

Thermal mass flow meter was developed using principle of convective heat transfer. The advantage of thermal mass flow meter is measuring mass flow directly, therefore, it is not required to use densitometer or temperature/pressure and DP gages. The final accuracy of this thermal mass flow meter is $\pm1.0{\%}$ or better, reproducibility is $\pm0.2{\%}$, and the response time is 600 ms. The thermal mass flow meter was developed from a single point to multi-points (maximum is 9 points), and the number of points is determined according to desired accuracy and size of piping/duct. Since this thermal mass flow meter adopted microprocessor based design, it is intrinsically accurate, self-error detectable, and has self-diagnosis function. The applications of this thermal mass flow meter are for measurement and control of HVAC air flow, other gas flow, and liquid flow.

• Journal of Korean Society of Industrial and Systems Engineering
• /
• v.20 no.44
• /
• pp.81-92
• /
• 1997

The purpose of this paper is to develop a method of paramatric analysis that can be applied to an optimal solution of a maximum flow problem. We first define the transformed network corresponding to a given network. In such a network, we conduct paramatric analysis by determining changes in the optimal solution precipitated by changes in the capacity as the arc capacity varies from 0 to infinite. By this method we can easily calculate not only the characteristic region where the given optimal solution remains unchanged, but also the characteristic region where the value of the maximal flow gradually increases or decreases. The proposed method is demonstrated by numerical example.

• Journal of Korean Institute of Industrial Engineers
• /
• v.23 no.4
• /
• pp.825-833
• /
• 1997

The purpose of this paper is to develop a method of postoptimality analysis that can be applied to an optimal solution of a maximum flow problem. We first use the transformed network corresponding to a given network. In such a network we conduct postoptimality analysis by determining changes in the optimal solution precipitated by changes in the capacity as the arc capacity varies from 0 to infinite. By this method we can easily calculate not only the characteristic region where the given optimal solution remains unchanged, but also the characteristic region where the value of the maximal flow gradually increases or decreases. The proposed method is demonstrated by numerical example.

• Journal of the Korean Society for Precision Engineering
• /
• v.22 no.3 s.168
• /
• pp.154-161
• /
• 2005

In this study, the leaflet motion of a mechanical heart valve and the characteristics of two-dimensional transient blood flow in an elastic blood vessel have been numerically investigated by using fluid-structure interaction method. Here, blood has been assumed as a Newtonian, incompressible fluid. Pressure profiles have been used as boundary conditions at the ventricle and the aorta. As a result, closing motion of the leaflet is faster than opening one. While opening angles of leaflet grow up, vortex is detected at the sinus and backward of the leaflets. When the leaflet is fully closed, vortex is detected at the ventricle and at that moment maximum displacement of the elastic blood vessel is observed in the vicinity of the sinus region. Maximum displacement is caused in association with the blood flow that is oriented toward the elastic blood vessel.

• Nuclear Engineering and Technology
• /
• v.28 no.4
• /
• pp.405-414
• /
• 1996

In this paper, the unsteady 2-dimensional turbulent flow model for thermal stratification in a pressurizer surge line of PWR plant is proposed to numerically investigate the heat transfer and flow characteristics. The turbulence model is adapted to the low Reynolds number K-$\varepsilon$ model (Davidson model). The dimensionless governing equations are solved by using the SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm. The results are compared with simulated experimental results of TEMR Test. The time-dependent temperature profiles in the fluid and pipe nil are shown with the thermal stratification occurring in the horizontal section of the pipe. The corresponding thermal stresses are also presented. The numerical result for thermal stratification by the outsurge during heatup operation of PWR shows that the maximum dimensionless temperature difference is about 0.83 between hot and cold sections of pipe well and the maximum thermal stress is calculated about 322MPa at the dimensionless time 28.5 under given conditions.

• Journal of the Korean Society of Combustion
• /
• v.15 no.3
• /
• pp.67-73
• /
• 2010

The temperature distribution of diesel particulate filter with five partitioned electric heaters is numerically analyzed to investigate the condition of regenerating ceramic filter. The commercial code STAR-$CCM+^{(R)}$ is utilized to simulate multi-dimensional steady hot air flow in DPF. In order to verify the computational results, thermocouples are used to measure the temperature distribution in DPF. Computational results agree well with experimental ones. The results show that the maximum temperature in DPF is lowered as the mass flow rate of exhaust gas increases, which means that the more power in heater will be necessary as the engine speed increases. Compared with heater placed at center, heater at circumference has the higher maximum temperature in DPF. The maldistribution of flow field in front of heater has the main influence on the temperature distribution in DPF.

• Journal of Advanced Marine Engineering and Technology
• /
• v.18 no.3
• /
• pp.17-25
• /
• 1994

Fully developed turbulent flow through three concentric annuli with both the rough inner and outer walls was experimentally investigated for a Reynolds number range Re = 15, 000 - 85, 000. Measurements were carried out for the pressure drop, the positions of zero shear stress and maximum velocity, and the velocity distributions in annuli of radius ratios, ${\alpha}$=0.26, 0.4 and 0.56, respectively. The experimental results showed that the positions of zero shear stress and maximum velocity were only weakly dependent on the Reynolds number.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.10 no.1
• /
• pp.65-70
• /
• 2001

The hear transfer mechanism initiating the friction welding is examined and a transient three dimensional heat conduc-tion model for the welding of two dissimilar cylindrical metal bars is investigated. The cylindrical metal bars are made of materials made of A2024 and SM 45C. Numerical simulations of heat flow are performed using the finite volume method. Respectively. Commercial FLUENT code is used in the heat flow simulation and maximum temperature and distribution of temperature are calculated. Temperature of friction welded joining face is compared with the temperature distribution measured by experiment and numerical simulation. The maximum temperature of friction welded joining face is lower than melting point of A2024-T6 aluminum alloy using insert metal. The temperature distribution of friction welded join- ing face with insert metal is more uniform than that of without inset metal.