• Nuclear Engineering and Technology
• /
• v.56 no.4
• /
• pp.1310-1319
• /
• 2024

The reactor cavity cooling system (RCCS) is a passive reactor safety system commonly present in the designs of High-Temperature Gas-cooled Reactors (HTGR) that removes heat from the reactor pressure vessel by means of natural convection and radiation. It is one of the factors responsible for ensuring that the reactor does not melt down under any plausible accident scenario. For the simulation of accident scenarios, which are transient phenomena unfolding over a span of up to several days, intermediate fidelity methods and system codes must be employed to limit the models' execution time. These models can quantify radiation heat transfer well, but heat transfer caused by natural convection must be quantified with the use of correlations for the heat transfer coefficient. It is difficult to obtain reliable correlations for HTGR RCCS heat transfer coefficients experimentally due to such a system's size. They could, however, be obtained from high-fidelity steady-state simulations of RCCSs. The Rayleigh number in RCCSs is too high for using a Direct Numerical Simulation (DNS) technique; thus, a Reynolds-Averaged Navier-Stokes (RANS) approach must be employed. There are many RANS models, each performing best under different geometry and fluid flow conditions. To find the most suitable one for simulating an RCCS, the RANS models need to be validated. This work benchmarks various RANS models against three experiments performed on the HTTR RCCS Mockup by the Japanese Atomic Energy Agency (JAEA) in 1993. This facility is a 1/6 scale model of a vessel cooling system (VCS) for the High Temperature Engineering Test Reactor (HTTR), which is operated by JAEA. Multiple RANS models were evaluated on a simplified 2d-axisymmetric geometry. They were found to reproduce the experimental temperature profiles with errors of up to 22% for the lowest temperature benchmark and 15% for the higher temperature benchmarks. The results highlight that the pragmatic turbulence models need to be validated for high Rayleigh natural convection-driven flows and improved accordingly, more publicly available experimental data of RCCS resembling experiments is needed and indicate that a 2d-axisymmetric geometry approximation is likely insufficient to capture all the relevant phenomena in RCCS simulations.

• Wind and Structures
• /
• v.7 no.5
• /
• pp.317-332
• /
• 2004

Wind tunnel experiments were conducted under highly turbulent and disturbed flow conditions over a solid/perforated plate with a long splitter plate in its plane of symmetry. The effect of varied level of perforation of the normal plate on fluctuating velocities and fluctuating pressures measured across and along the separation bubble was studied. The different perforation levels of the normal plate; that is 0%, 10%, 20%, 30%, 40% and 50% are studied. The Reynolds number based on step height was varied from $4{\times}10^3$ to $1.2{\times}10^4$. The shape and size of the bubble vary with different perforation level of the normal plate that is to say the bubble is reduced both in height and length up to 30% perforation level. For higher perforation of the normal plate, bubble is completely swept out. The peak turbulence value occurs around 0.7 to 0.8 times the reattachment length. The turbulence intensity values are highest for the case of solid normal plate (bleed air is absent) and are lowest for the case of 50% perforation of the normal plate (bleed air is maximum in the present study). From the analysis of data it is observed that $\sqrt{\overline{u^{{\prime}2}}}/(\sqrt{\overline{u^{{\prime}2}}})_{max}$, (the ratio of RMS velocity fluctuation to maximum RMS velocity fluctuation), is uniquely related with dimensionless distance y/Y', (the ratio of distance normal to splitter plate to the distance where RMS velocity fluctuation is half its maximum value) for all the perforated normal plates. It is interesting to note that for 50% perforation of the normal plate, the RMS pressure fluctuation in the flow field gets reduced to around 60% as compared to that for solid normal plate. Analysis of the results show that the ratio [$C^{\prime}_p$ max/$-C_{pb}(1-{\eta})$], where $C^{\prime}_p$ max is the maximum coefficient of fluctuating pressure, $C_{pb}$ is the coefficient of base pressure and ${\eta}$ is the perforation level (ratio of open to total area), for surface RMS pressure fluctuation levels seems to be constant and has value of about 0.22. Similar analysis show that the ratio $[C^{\prime}_p$ max/$-C_{pb}(1-{\eta})]$ for flow field RMS pressure fluctuation levels seems to be constant and has a value of about 0.32.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.28 no.7
• /
• pp.1244-1251
• /
• 2022

Marine accidents due to a loss of stability, have been gradually increasing over the last decade. Measures must be taken on the roll reduction of a ship. Amongst the measures, building an anti-roll tank in a ship is recognized as the most simple and effective way to reduce the roll motion. Therefore, this study aims to develop a computational model for a U-tube type anti-roll tank and to validate it by experiment. In particular, to validate the developed computational model, the height of the free surface in the tank was measured in the experiment. To develop a computational model, the mesh dependency test was carried out. Further, the effects of a turbulence model, time step size, and the number of iterations on the numerical solution were analyzed. In summary, a U-tube type anti-roll tank simulation had to be performed accurately with conditions of a realizable k-𝜖 turbulence model, 10^-2s time step size, and 15 iterations. In validation, the two cases of measured data from the experiment were compared with the numerical results. In the present study, STAR-CCM+ (ver. 17.02), a RANS-based commercial solver was used.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.127-136
• /
• 2022

In this work a multi-fidelity non-intrusive polynomial chaos (MF-NIPC) has been applied to a structural wind engineering problem in architectural design for the first time. In architectural design it is important to design structures that are safe in a range of wind directions and speeds. For this reason, the computational models used to design buildings and bridges must account for the uncertainties associated with the interaction between the structure and wind. In order to use the numerical simulations for the design, the numerical models must be validated by experi-mental data, and uncertainties contained in the experiments should also be taken into account. Uncertainty Quantifi-cation has been increasingly used for CFD simulations to consider such uncertainties. Typically, CFD simulations are computationally expensive, motivating the increased interest in multi-fidelity methods due to their ability to lev-erage limited data sets of high-fidelity data with evaluations of more computationally inexpensive models. Previous-ly, the multi-fidelity framework has been applied to CFD simulations for the purposes of optimization, rather than for the statistical assessment of candidate design. In this paper MF-NIPC method is applied to flow around a rectan-gular 5:1 cylinder, which has been thoroughly investigated for architectural design. The purpose of UQ is validation of numerical simulation results with experimental data, therefore the radius of curvature of the rectangular cylinder corners and the angle of attack are considered to be random variables, which are known to contain uncertainties when wind tunnel tests are carried out. Computational Fluid Dynamics (CFD) simulations are solved by a solver that employs the Finite Element Method (FEM) for two turbulence modeling approaches of the incompressible Navier-Stokes equations: Unsteady Reynolds Averaged Navier Stokes (URANS) and the Large Eddy simulation (LES). The results of the uncertainty analysis with CFD are compared to experimental data in terms of time-averaged pressure coefficients and bulk parameters. In addition, the accuracy and efficiency of the multi-fidelity framework is demonstrated through a comparison with the results of the high-fidelity model.

• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1049-1058
• /
• 2023

The flow passing through river-crossing structures such as weirs and low-fall dams is dominated by rapidly varied flow including hydraulic jump. The intense unsteadiness of flow velocity and free surface profile affects the stability of such hydraulic structures. In particular, the steady hydraulic jump generated at high Froude number conditions includes remarkably air entrainment, making the flow characteristics more complicated. In this study, a large-eddy simulation was performed for turbulence effect and the hybrid VoF technique to simulate the steady hydraulic jump at the Froude number of 7.3 and the Reynolds number of 15,700. The results of the numerical simulation showed that the instantaneous maximum pressure and time-average pressure distribution calculated on the bottom surface downstream of the structure could be reasonably well reproduced being in good agreement with the experimental values. However, the instantaneous minimum pressure distribution in the direct downstream of the structure shows the opposite pattern to the target experimental measurement value. However, the numerical simulation performed in this study is considered to reasonably predict the minimum pressure distributions observed in various experiments conducted at similar conditions. The vertical distributions of flow velocity and air concentration computed in the center of the hydraulic jump were found to be in good agreement with the experimental results measured under similar conditions, showing self-similarity. These results show that the large eddy simulation and hybrid VoF techniques applied in this study can reproduce the hydraulic jump with strong air entrainment and the resulting intense free surface and pressure fluctuations at high Froude number conditions.

• Solar Energy
• /
• v.13 no.1
• /
• pp.11-21
• /
• 1993

The purpose of this study was to investigate the enhancement of heat transfer without additional external power in the case of rectangular impingement air jet vertically on the flat heating surface. The technique used in the present study was placement of square rod bundles as a turbulence promoter in front of the heat transfer surface. The results obtained through this study were summerized as follws. High heat transfer enhancement was achieved by inserting rods in front of the heating flat plate. According to visulaization, it was examined because of flow acceleration and separation and disturbance of boundary layer. The smaller clerance between rod and heating plate was, the larger heat transfer effect became at each H/B. Arverage Nusselt number reached maximum at H/B=10 and the local augmentation rate of heat transfer became maximum at H/B=2. The maximum average heat transfer enhancement rate increase about 43% for the case of X/B=2 and C=1mm, compared to a flat plate without rods. The correlating equation of average Nusselt number and Reynolds number was obatined. As follws : ${\overline{Nu}}_0=1.249Re^{0.465}(C/A)^{-0.033}(H/B)^{0.013}$.

• Korean Chemical Engineering Research
• /
• v.59 no.1
• /
• pp.94-99
• /
• 2021

It was found that mixing patterns suddenly changed at an impeller rotation speed in a dished bottom vessel with dual Rushton turbines. Two isolated mixing regions like doughnuts rings generated at a low rotational speed and three isolated mixing regions generated at a higher speed. This phenomenon was observed at the mixing condition in transition area, where the power number with baffle was the same as that without baffle. We found a phenomenon in which the flow state in a dish-bottom agitation tank equipped with a two-stage Rushton turbine blade changes at a certain rotational speed. In the laminar flow region, the isolated stable donut rings were formed even when the rotational speed was changed, and no specific variation in the mixing pattern was observed. In the transition region, the two isolated thick unmixed donut rings do not change even if the rotation speed is changed in the flat bottom vessel, whereas in the dished bottom vessel, when the rotation speed is 450 rpm, the two isolated thick unmixed donut rings are changed to three isolated thin donut rings and then improved mixing. In the dished bottom vessel, in the range of Re=138~178, the isolated ring-shaped unmixed region appeared in three places and the size was also large. But in the flat bottom vessel, the isolated thick ring-shaped unmixed region appeared in two places in Re=116~176 and the size was also small. It appeared in two places, and the size was also small. The condition in which this phenomenon is observed is a transition region, and it was found that when the baffle plate is attached, the power number starts to increase compared to when the baffle plate is not present. In addition, when the mixing Reynolds number exceeded 300 and a slight turbulence was mixed in the flow state, the disconnection of these flow pattern was resolved and the mixture was completely mixed.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.18 no.2
• /
• pp.235-242
• /
• 2017

A driving controller for active destruction protections can be applied to machinery, aerospace and military fields. In particular, this controller can be used to track and attack enemy flying objects through the active control. It is important to ensure reliability of the driving controller since its operation should be kept with precision to the target point. The temperature of the environment where the driving controller is used is about -32 C ~ 50 C (241~323 ). Heat generated in the driving controller should be maintained below a certain threshold (85 C (358 )) to ensure reliability; therefore, the study and analysis of the heat flow characteristics in the driving controller are required. In this research, commercial software Solid-Works Flow Simulation was used for the numerical simulation assuming a low Reynolds number turbulence model and an incompressible viscous flow. The goal of this paper is to design the driving controller safely by analyzing the characteristics of the heat flow inside of the controller composed of chips or boards. Our analysis shows temperature distributions for boards and chips below a certain threshold.

• Fire Science and Engineering
• /
• v.33 no.4
• /
• pp.35-40
• /
• 2019

In this study, the mass flow rate of the heat release rate equation, which is the major factor of the oxygen consumption method, was analyzed for the fundamental investigation of the cone-calorimeter (5 m length and 0.3 m diameter). The shapes of a completely empty inside, 3 mm pore diameter mesh and pore diameter 10 mm honeycomb with 0.76 porosity were constructed using the cone-calorimeter. To calculate the mass flow rate, four bi-directional probes and thermocouples were installed in a uniform position in the vertical direction of flow. The velocity gradient and flow perturbation were measured from the increase in Reynolds number. As the flow capacity increased, the speed gradient increased in all three shapes relative to the turbulence intensity. In addition, the deviation of extended uncertainty to the mass flow was completely low in the order of empty space, mesh (d_p = 3 mm) and honeycomb (d_p = 10 mm and 𝜖 = 0.76) at the 95% confidence level. The results can be used in designs to improve the flow stability of the cone calorimeter.