• Nuclear Engineering and Technology
• /
• v.56 no.7
• /
• pp.2524-2533
• /
• 2024

When analyzing containment thermal-hydraulics, computational fluid dynamics (CFD) is a powerful tool because multi-dimensional and local analysis is required for some accident scenarios. According to the previous study, neglecting steam bulk condensation in the CFD analysis leads to a significant error in boundary layer profiles. Validating the condensation model requires the experimental data near the condensing surface, however, available boundary layer data is quite limited. It is also important to confirm whether the heat and mass transfer analogy (HMTA) is still valid in the presence of bulk condensation. In this study, the boundary layer measurements on the vertical condensing surface in the presence of air were performed with the rectangular channel facility WINCS, which was designed to measure the velocity, temperature, and concentration boundary layers. We set the laminar flow condition and varied the Richardson number (1.0-23) and the steam volume fraction (0.35-0.57). The experimental results were used to validate CFD analysis and HMTA models. For the former, we implemented a bulk condensation model assuming local thermal equilibrium into the CFD code and confirmed its validity. For the latter, we validated the HMTA-based correlations, confirming that the mixed convection correlation reasonably predicted the sum of wall and bulk condensation rates.

• Journal of the Korean Society of Visualization
• /
• v.21 no.1
• /
• pp.47-56
• /
• 2023

Bubble condensation, which involves the interaction of bubbles within the subcooled liquid flow, plays an important role in the effective control of thermal devices. In this study, numerical simulations are performed using a VOF (Volume of Fluid) model to investigate the effect of tube diameter on bubble condensation. As the tube diameter decreases, condensation bubbles persist for a long time and disappear at a higher position. It is observed that for small tube diameters, the heat transfer coefficients of condensation bubbles, which is a quantitative parameter of condensation rate, are smaller than those for large tube diameters. When the tube diameter is small, the subcooled liquid around the condensing bubble is locally participated in the condensation of the bubble to fill the reduced volume of the bubble due to the generation of a backflow in the narrow space between the bubble and the wall, so that the heat transfer coefficient decreases.

• Proceedings of the KSME Conference
• /
• 2001.06d
• /
• pp.813-820
• /
• 2001

The Safety depressurization System(SDS) of KNGR prevents RCS from overpressurization by discharging high pressure and temperature coolant through the I-sparger into the IRWST during an accident. If IRWST water temperature rise locally, around the sparger, beyond $200_{\circ}$2000 F by the discharged coolant, unstable steam condensation can cause large pressure load on the IRWST wall. To investigate whether this condition can be avoided for the design basis event IOPOSRV(Inadvertent Opening of one Pilot Operated Safety Relief Valve), the flow and temperature distribution of water in the IRWST is calculated by using CFX 4.3 computational fluid dynamic code. According to the results, since pool water temperature does not exceeds temperature limit within 50 seconds after the opening of one POSRV, it can be assured that the integrity of IRWST wall is maintained.

• Nuclear Engineering and Technology
• /
• v.55 no.4
• /
• pp.1400-1409
• /
• 2023

In this study, three condensation models of the CUPID code, i.e., the resolved boundary layer approach (RBLA), heat and mass transfer analogy (HMTA) model, and an empirical correlation, were tested and validated against the COPAIN and CAU tests. An improvement on HMTA model was also made to use well-known heat transfer correlations and to take geometrical effect into consideration. The RBLA was a best option for simulating the COPAIN test, having mean relative error (MRE) about 0.072, followed by the modified HMTA model (MRE about 0.18). On the other hand, benchmark against CAU test (under natural convection and occurred on a slender tube) indicated that the modified HMTA model had better accuracy (MRE about 0.149) than the RBLA (MRE about 0.314). The HMTA model with wall function and the empirical correlation underestimated significantly, having MRE about 0.787 and 0.55 respectively. When using the HMTA model, consideration of geometrical effect such as tube curvature was essential; ignoring such effect leads to significant underestimation. The HMTA and the empirical correlation required significantly less computational resources than the RBLA model. Considering that the HMTA model was reasonable accurate, it may be preferable for large-scale simulations of containment.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1998.05a
• /
• pp.547-552
• /
• 1998

The standard RELAP5/MOD3.2 code were assessed with the condensation experiment in the presence of noncondensable gas in a vortical tube of PCCS of CP-1300. There are two wall film condensation models, the default model and the alternative model, in RELAP5/MOD3.2. The experimental apparatus was modeled with the two models, md simulations were performed for several sub-tests to be compared with the experimental results. In overall sense the simulation results showed that the default model of RELAP5/MOD3.2 under-predicts the heat transfer coefficients, while the alternative model over-predicts them throughout the condensing tube.

• Nuclear Engineering and Technology
• /
• v.33 no.2
• /
• pp.241-253
• /
• 2001

Under severe accidents, the pressure and temperature response has an important role for the integrity of a nuclear power plant containment. The history of the pressure and temperature is characterized by the amount and state of steam/air mixture in a containment. Recently, the heat transfer rate to the structure surface is supposed to be increased by the wavy interface formed on condensate film. However, in the calculation by using CONTAIN code, the condensation heat transfer on a containment wall is calculated by assuming the smooth interface and has a tendency to be underestimated for safety. In order to obtain the best- estimate heat transfer calculation, we investigated the condensation heat transfer model in CONTAIN 1.2 code and adopted the new forced convection correlation which is considering wavy interface. By using the film tracking model in CONTAIN 1.2 code, the condensate film is treated to consider the effect of wavy interface. And also, it was carried out to investigate the effect of the different cell modelings - 5-cell and 10-cell modeling - for KNGR(Korean Next Generation Reactor) containment phenomena during a severe accident. The effect of wavy interface on condensate film appears to cause the decrease of peak temperature and pressure response . In order to obtain more adequate results, the proper cell modeling was required to consider the proper flow of steam/air mixture.

• Particle and aerosol research
• /
• v.9 no.3
• /
• pp.163-171
• /
• 2013

Soot particles emitted from combustion processes are often coated by non-absorbing organic materials, which enhance the global warming effect of soot particles. It is of importance to study the condensation characteristics of soot particles experimentally and theoretically to reduce the uncertainty of the climate impact of soot particles. In this study, the condensational growth of soot particles in a tubular coater was modeled by a one-dimensional (1D) plug flow model and a two-dimensional (2D) laminar flow model. The effects of 2D heat and mass transports on the predicted particle growth were investigated. The temperature and coating material vapor concentration distributions in radial direction, which the 1D model could not accounted for, affected substantially the particle growth in the coater. Under the simulated conditions, the differences between the temperatures and vapor concentrations near the wall and at the tube center were large. The neglect of these variations by the 1D model resulted in a large error in modeling the mass transfer and aerosol dynamics occurring in the coater. The 1D model predicted the average temperature and vapor concentration quite accurately but overestimated the average diameter of the growing particles considerably. At the outermost grid, at which condensation begins earliest due to the lowest temperature and saturation vapor concentration, condensing vapor was exhausted rapidly because of the competition between condensations on the wall and on the particle surface, decreasing the growth rate. At the center of the tube, on the other hand, the growth rate was low due to high temperature and saturation vapor concentration. The effects of Brownian diffusion and thermophoresis were not high enough to transport the coating material vapor quickly from the tube center to the wall. The 1D model based on perfect radial mixing could not take into account this phenomenon, resulting in a much higher growth rate than what the 2D model predicted. The result of this study indicates that contrary to a previous report for a thermodenuder, 2D heat and mass transports must be taken into account to model accurately the condensational particle growth in a coater.

• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.80-86
• /
• 2005

In order to control the transonic flow field with a shock wave, a condensing flow was produced by an expansion of moist air on a circular bump model and shock waves were occurred in the supersonic parts of the fields. Furthermore, the additional passive technique of shock-boundary layer interaction using the porous wall with a cavity underneath was adopted in this flow field. The effects of these methods on the shock wave characteristics were investigated numerically. The result showed that the flow fields might be effectively controlled by the suitable combination between non-equilibrium condensation and the position of porous wall.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.187-188
• /
• 2003

In order to control the transonic flow field with shock wave, a condensing flow was produced by an expansion of moist air on a circular bump model and shock waves were occurred in the supersonic parts of the fields. Furthermore, the additional passive technique of shock - boundary layer interaction using the porous wall with a cavity underneath was adopted in this flow field. The effects of these methods on the shock wave characteristics were investigated numerically. The result showed that the flow fields might be effectively controlled by the suitable combination between non-equilibrium condensation and the position of porous wall.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.8 no.6 s.40
• /
• pp.55-66
• /
• 2004

Recently, many high-rise apartment buildings using the box system, composed of only reinforced concrete walls and slabs, have been constructed. In residential buildings such as apartments, vibrations occur from various sources and these vibrations transfer to neighboring residential units through walls and slabs. It is necessary to use a refined finite element model for an accurate vibration analysis of shear wall building structures. But it would take significant amount of computational time and memory if the entire building structure were subdivided into a finer mesh. Therefore, an efficient analytical method, which has only translational DOFs perpendicular to walls or slabs by the matrix condensation technique, is proposed in this study to obtain accurate results in significantly reduced computational time. If all of the DOFs except those perpendicular to walls or slabs in the shear wall structure eliminated using the matrix condensation technique at a time, the computational time for the matrix condensation would be significant. Thus, the modeling technique using super elements and substructuring technique is proposed to reduce the computational time for the matrix condensation. Dynamic analysis of 3-story and 5-story shear wall example structures were performed to verify the efficiency and accuracy of the proposed method. It was confirmed that the proposed method can provide the results with outstanding accuracy requiring significantly reduced computational time and memory.