• Nuclear Engineering and Technology
• /
• v.55 no.10
• /
• pp.3754-3767
• /
• 2023

Turbulent jet discharges in subcooled water pools are essential for safety systems in nuclear power plants, specifically in the pressure suppression pool of boiling water reactors and In-containment Refueling Water Storage Tank of advanced pressurized water reactors. The gas and liquid flow in these systems is investigated using multiphase flow analysis. This field has been extensively examined using a combination of experiments, theoretical models, and Computational Fluid Dynamics (CFD) simulations. ANSYS CFX offers two approaches to model multiphase flow behavior. The non-homogeneous Eulerian-Eulerian Model has been used in this work; it computes global information and is more convenient to study interpenetrated fluids. This study utilized the Large Eddy Simulation Model as the turbulence model, as it is better suited for non-stationary and buoyant flows. The CFD results of this study were validated with experimental data and theoretical results previously obtained. The figures of merit dimensionless penetration length and the dimensionless buoyancy length show good agreement with the experimental measurements. Correlations for these variables were obtained as a function of dimensionless numbers to give generality using only initial boundary conditions. CFD numerical model developed in this research has the capability to simulate the behavior of non-condensable gases discharged in water.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.1032-1043
• /
• 1989

The heat transfer characteristics of the drag reducing polymer solutions are investigated experimentally in the thermal entrance region of circular tube flows. Fluids used in experiments are the aqueous solutions of high molecular polymer, polyacrylamide Separan AP-273 and the range of polymer concentrations is from 20 to 1000 wppm. Two stainless steel tubes with inside diameter 8.5mm(L/D=712) and 10.3mm(L/D=1160) are used for the heat transfer flow loops. The flow loop is set up to measure friction factors and heat transfer coefficients of test sections in two different modes; the recirculating flow system and once-through flow system. The test tubes are heated directly by electricity to apply the constant heat flux boundary conditions to the wall. Three different types of adaptors are used to observe the effects of the upstream flow conditions of the heat transfer test sections. The viscosity and characteristic relaxation time of the test fluids circulating in the flow system are measured by the capillary tube viscometer and falling ball viscometer at regular time intervals. The installed adaptors exhibit slight effect on the entrance heat transfer of Newtonian fluid. However, no noticeable effects are observed for the entrance heat transfer of the drag reducing fluids. The order of magnitude of the thermal entrance lengths of the drag reducing fluids which follow the minimum friction asymptote is much longer than that of Newtonian fluids in turbulent flows. A new dimensionless parameter, the viscoelastic Graetz number, is defined and all the experimental data are recasted in terms of the viscoelastic Graetz number. The local Nusselt number of the viscoelastic fluids is represented as a function of flow behavior index n and the viscoelastic Graetz number. As degradation continues the viscosity and the characteristic relaxation time of the testing fluids decrease. Weissenberg number defined by the relaxation time and D/V appears to be a proper dimensionless parameter in describing degradation effects on heat transfer of the viscoelastic fluids.

• Journal of Korea Water Resources Association
• /
• v.57 no.8
• /
• pp.493-507
• /
• 2024

RANS-based CFD analysis is widely applied in various engineering fields, including practical hydraulic engineering, due to its high computational efficiency. However, problems of non-physical behavior in the analysis of two phase flow, such as free surfaces, have long been raised. The two-equation turbulence models used in general RANS-based analysis were developed for single phase flow and simulate unrealistically high turbulence energy at the interface where there are abrupt changes in fluid density. To solve this issue, one of the methods recently developed is the buoyancy-modified turbulence model, which has been partially validated in coastal engineering, but has not been applied to open channel flows. In this study, the applicability of the buoyancy-modified turbulence model is evaluated using the VOF method in the open-source program OpenFoam. The results of the uniform flow showed that both the buoyancy-modified k-𝜖 model and the buoyancy-modified k-ω SST model effectively simulated the reduction of turbulence energy near the free surface. Specifically, the buoyancy-modified k-ω SST model accurately simulated the vertical velocity distribution. Additionally, the model is applied to dam-break flows to examine cases with significant surface variation and cavity formation. The simulation results show that the buoyancy-modified turbulence models produce varying results depending on the VOF method and shows non-physical behavior different from experimental results. While the buoyancy-modified turbulence model is applicable in cases with stable surface shapes, it still has limitations in general application when there are rapid changes in the free surface. It is concluded that appropriate adjustments to the turbulence model are necessary for flows with rapid surface changes or cavity formation.

• Journal of Korea Water Resources Association
• /
• v.53 no.1
• /
• pp.1-13
• /
• 2020

The most common approach for computing engineering flow problems at high Reynolds number is still the Reynolds-averaged Navier-Stokes (RANS) computations based on turbulence models with wall functions. The recently developed generalized wall functions blending between the wall-limiting viscous and the outer logarithmic relations ensure a smooth transition of flow quantities across two regions. The performances and convergence properties of widely used turbulence models with wall functions that are applicable for turbulence kinetic energy (TKE), turbulent and specific dissipation rates, and eddy viscosity are presented through a series of near wall flow simulations. The present results show that RNG k-𝜖 model should be carefully applied with small tolerance to get the stable solution when the first grid lies in the buffer layer. The standard k-𝜖 and RNG k-𝜖 models are not sensitive to the selection of wall functions for both TKE and eddy viscosity, while the k-ω SST model should be applied together with kL-wall function for TKE and nutUB-wall functions for eddy viscosity to ensure accurate and stable boundary conditions. The applications to a backward-facing step flow at Re=155,000 reveal that the reattachment length is reasonably well predicted on appropriately refined mesh by all turbulence models, except the standard k-𝜖 model which about 13% underestimates the reattachment length regardless of the grid refinement.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.18 no.3
• /
• pp.673-680
• /
• 2017

sump pump is a system that draws in water that is stored in a dam or reservoir. They are used to pump large amounts of water for cooling systems in large power plants, such as thermal and nuclear plants. However, if the flow and sump pump ratio are small, the flow rate increases around the inlet port. This causes a turbulent vortex or swirl flows. The turbulent flow reduces the performance and can cause failure. Various methods have been devised to solve the problem, but a correct solution has not been found for low water level. The most efficient solution is to install an anti-vortex device (AVD) or increase the length of the sump inlet, which makes the flow uniform. This paper presents a computational fluid dynamics (CFD) analysis of the flow characteristics in a sump pump for different sump inlet lengths and AVD types. Modeling was performed in three stages based on the pump intake, sump, and pump. For accurate analysis, the grid was made denser in the intake part, and the grid for the sump pump and AVD were also dense. 1.2-1.5 million grid elements were generated using ANSYS ICEM-CFD 14.5 with a mixture of tetra and prism elements. The analysis was done using the SST turbulence model of ANSYS CFX14.5, a commercial CFD program. The conditions were as follows: H.W.L 6.0 m, L.W.L 3.5, Qmax 4.000 kg/s, Qavg 3.500 kg/s Qmin 2.500 kg/s. The results of analysis by the vertex angle and velocity distribution are as follows. A sump pump with an Ext E-type AVD was accepted at a high water level. However, further studies are needed for a low water level using the Ext E-type AVD as a base.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.24 no.5
• /
• pp.619-627
• /
• 2018

As ships become faster, larger and are required to meet higher standards, the importance of flow noise is highlighted. However, unlike in the aeroacoustics field for airplanes and trains (where flow noise is considered in design), acoustics are not considered in the marine field. In this study, analysis procedures for hull-induced flow noise are established to investigate the flow noise characteristics of a wave-piercing hull form that can negate the effect of wave-breaking. The principal mechanisms behind hull-induced flow noise are fluid-structure interactions between complex flows underneath the turbulent boundary layer and the hull. Noise induced by the turbulent boundary layer was calculated using wall pressure fluctuation and energy flow analysis methods. The results obtained show that noise characteristics can be distinguished by frequency range and hull region. Also, the low-frequency range is affected by hull forms such that it is correlated with ship speed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.575-589
• /
• 2008

In the present work, we investigate the hydrodynamic behavior of a turbulent bore, such as tsunami bore and tidal bore, generated by the removal of a gate with water impounded on one side. The bore generation system is similar to that used in a general dam-break problem. In order to the numerical simulation of the formation and propagation of a bore, we consider the incompressible flows of two immiscible fluids, liquid and gas, governed by the Navier-Stokes equations. The interface tracking between two fluids is achieved by the volume-of-fluid (VOF) technique and the M-type cubic interpolated propagation (MCIP) scheme is used to solve the Navier-Stokes equations. The MCIP method is a low diffusive and stable scheme and is generally extended the original one-dimensional CIP to higher dimensions, using a fractional step technique. Further, large eddy simulation (LES) closure scheme, a cost-effective approach to turbulence simulation, is used to predict the evolution of quantities associated with turbulence. In order to verify the applicability of the developed numerical model to the bore simulation, laboratory experiments are performed in a wave tank. Comparisons are made between the numerical results by the present model and the experimental data and good agreement is achieved.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.4
• /
• pp.417-424
• /
• 2011

Natural convection heat transfer experiments from two parallel horizontal cylinders were performed varying the Pitch-to-Diameter ratio (P/D) of 1.02-9 at Sc of 2,014 to 8,334 and $Ra_D$ of $1.5{\times}10^8$ to $4.5{\times}10^{10}$. Mass transfer experiments that are analogous to the heat transfer experiments were performed using copper electroplating system. In all cases, the measured heat transfer rates for the lower cylinder agreed well with the existing heat transfer correlations developed from a single cylinder. For laminar flows, the measured heat transfer rates of the upper cylinder were less than those of the lower cylinder at P/D less than about 1.5. However, as the P/D increased, the heat transfer rates of the upper cylinder increased. For turbulent flows, the heat transfer rates of the upper cylinder were considerably similar to those of the lower cylinder when the P/D is approximately unity. In contrast, as the P/D increased, the heat transfer rates of the upper cylinder were always higher than those of the lower cylinder.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.9
• /
• pp.126-132
• /
• 2016

This study numerically investigates the characteristics of chemical reactions and thermal deformation in a steam reformer. These phenomena are significantly affected by the high-temperature burner gas and the process gas conditions. Because the high temperature of the burner gas ranges from 800 to 1000 K, the reformer tubes undergo substantial thermal deformation, eventually resulting in structural failure. Thus, it is necessary to understand the characteristics of the reaction and thermal deformation under the operating conditions to evaluate the reformer tubes for sustainable, stable operation. Extensive numerical simulations were carried out using commercial CFD code (ANSYS FLUENT/MECHANICA Ver. 13.0) while considering three-dimensional turbulent flows and combined heat transfer including conduction, convection, and radiation. Structural analysis considering conjugated heat transfer between solid tubes and fluid flows was conducted using the Fluid-Solid Interaction (FSI) method. The results show that when the injection temperature of the process gas and burner gas decreased, the hydrogen production rate decreased significantly, and thermal deformation decreased by at least 15 to 20%.

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.7-13
• /
• 2010

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 5.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}10^6$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental datum in the low aspect ratio cavity (L/D = ~4.5). In the high aspect ratio cavity, however, there are other low dominant frequencies of the leading edge shear layer with the dominant frequencies of the feedback mechanism.