• Nuclear Engineering and Technology
• /
• v.56 no.1
• /
• pp.233-240
• /
• 2024

The wall friction correlations of oscillatory natural circulation loops are highly loop-specific, making it difficult to perform 1-D system simulations before obtaining specific experimental data. To better predict the friction characteristics, the nonlinear dynamics of a toroidal single-phase natural circulation loop were numerically investigated, and the transition effect was considered. The k-k_L-ω transitional turbulence and k-ω SST turbulence models were used to compute the flow characteristics of the loop under different heating powers varying from 0.48 to 1.0 W/cm², and the results of both models were compared with previous experiments. The mass flow rates and friction factors predicted by the k-k_L-ω model showed a better agreement with the experimental data than the results of the k-ω SST model. The oscillation frequencies calculated using both models agreed well with the experimental data. The k-k_L-ω transitional turbulence model provided better friction-factor predictions in oscillatory natural circulation loops because it can reproduce the temporal and spatial variation of the wall shear stress more accurately by capturing the movement of laminar, transition turbulent zones inside unstable natural circulation loops. This study shows that transition effects are a possible explanation for the highly loop-specific friction correlations observed in various oscillatory natural circulation loops.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.44 no.3
• /
• pp.329-336
• /
• 2024

In the case of multi-function weirs installed in Korea, the free hydraulic jump or the submerged hydraulic jump is occurred depending on the height of the gate opening and the tailwater level when the sluice gate of the movable weir is partially opened. In this study, the submerged hydraulic jump for the flows under the sluice gate were simulated and the mean flow, turbulence statistics, and relative water depth are investigated using numerical simulation. For numerical simulation, the unsteady Reynolds-averaged Navier-Stokes equation, volume of fluid method, and k-ωSST turbulence model were used. The numerical model was validated using the results of other researchers' previously performed experiments, and it was investigated that the numerical model appropriately simulates the two-phase flow in the hydraulic jump. In addition, the distribution of mean flow, turbulence statistics, and the length of recirculation region was investigated.

• 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 computational fluids engineering
• /
• v.9 no.2
• /
• pp.30-35
• /
• 2004

A numerical study of underexpanded jet and impingement on a wall mounted at various distances from the nozzle exit is presented. The 3-dimensional Wavier-Stokes equations and κ-ω turbulence equations are solved. The grids are constructed as overlapped grid systems to examine the distance effect. The DADI method is applied to obtain steady-state solutions. To avoid numerical instability such as the carbuncle phenomena that sometimes accompany approximate Riemann solver, the HLLE+ scheme is employed for the inviscid flux at the cell interfaces. A goal of this work is to apply a number of two-equation turbulence models based on the w equation to the impinging jet problem.

• Wind and Structures
• /
• v.37 no.4
• /
• pp.303-314
• /
• 2023

In this study, two steady RANS turbulence models (SST k-ω and Realizable k-ε) and four unsteady turbulence models (URANS SST k-ω and Realizable k-ε, SST-SAS, and SST-IDDES) are evaluated with respect to their capacity to predict crosswind characteristics on high-speed trains (HSTs). All of the numerical simulations are compared with the wind tunnel values and LES results to ensure the accuracy of each turbulence model. Specifically, the surface pressure distributions, time-averaged aerodynamic coefficients, flow fields, and computational cost are studied to determine the suitability of different models. Results suggest that the predictions of the pressure distributions and aerodynamic forces obtained from the steady and transient RANS models are almost the same. In particular, both SAS and IDDES exhibits similar predictions with wind tunnel test and LES, therefore, the SAS model is considered an attractive alternative for IDDES or LES in the crosswind study of trains. In addition, if the computational cost needs to be significantly reduced, the RANS SST k-ω model is shown to provide relatively reasonable results for the surface pressures and aerodynamic forces. As a result, the RANS SST k-ω model might be the most appropriate option for the expensive aerodynamic optimizations of trains using machine learning (ML) techniques because it balances solution accuracy and resource consumption.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.26 no.4
• /
• pp.419-428
• /
• 2020

The flow past a wall mounted square cylinder, a typical and basic shape of building, bridge or offshore structure, was simulated using URANS computation through adoption of three turbulence models, namely, the k-ε model, k-ω model, and the v²-f model. It is well known that this flow is naturally unstable due to the Karman vortex shedding and exhibits a complex flow structure in the wake region. The mean flow field including velocity profiles and the dominant frequency of flow oscillation that was from the simulations discussed earlier were compared with the experimental data observed by Wang et al. (2004; 2006). Based on these comparisons it was found that the v²-f model is most accurate for the URANS simulation; moreover, the k-ω model is also acceptable. However, the k-ε model was found to be unsuitable in this case. Therefore, v²-f model is proved to be an excellent choice for the analysis of flow with massive separation. Therefore, it is expected to be used in future by studies aiming to control the flow separation.

• Wind and Structures
• /
• v.31 no.4
• /
• pp.351-362
• /
• 2020

In this study, the impact of roof slope on the flow characteristics over low-sloped gable roofs was investigated using steady computational fluid dynamics (CFD) simulations based on a k-ω SST turbulence model. A measurement database of the flow field over a scaled model of 15° was created using particle image velocimetry (PIV). Sensitivity analyses for the grid resolutions and turbulence models were performed. Among the three common Reynolds-averaged Navier-Stokes equations (RANS) models, the k-ω SST model exhibited a better performance, followed by the RNG model and then the realizable k-ε model. Next, the flow properties over the differently sloped (0° to 25°) building models were determined. It was found that the effect of roof slope on the flow characteristics was identified by changing the position and size of the separation bubbles, 15° was found to be approximately the sensitive slope at which the distribution of the separation bubbles changed significantly. Additionally, it is suggested additional attention focused on the distributions of the negative pressure on the windward surfaces (especially 5° and 10° roofs) and the possible snow redistribution on the leeward surfaces.

• Nuclear Engineering and Technology
• /
• v.54 no.11
• /
• pp.4296-4309
• /
• 2022

Mesh and turbulence model sensitivity analyses have been performed on computational fluid dynamics simulations executed with Hydra and ANSYS Fluent for a single CANadian Deuterium Uranium (CANDU) 37M nuclear fuel bundle placed within a standard pressure tube. The goal of this work was to perform a methodical analysis to objectively determine an appropriate mesh and to gauge the sensitivity of different turbulence models for CANDU subchannel flow under isothermal conditions. The boundary conditions and material properties are representative of normal operating conditions in a high-powered channel of the Darlington Nuclear Generating Station. Four meshes were generated with ANSYS Workbench Meshing, ranging from 22 to 84 million cells, and analyzed here to determine an appropriate level of mesh resolution and quality. Five turbulence models were compared in the turbulence model sensitivity analysis: standard k - ε, RNG k - ε, realizable k - ε, SST k - ω, and the Reynolds Stress Model. The intent of this work was to gain confidence in mesh generation and turbulence model selection of a single bundle to inform the decision making of subsequent investigations of an entire fuel channel containing a string of twelve bundles.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.5
• /
• pp.375-384
• /
• 2012

This paper presents numerical results comparing the performance of the 2008 Wilcox $\mathcal{k}-{\omega}$ turbulence model to the one of the 1988 Wilcox $\mathcal{k}-{\omega}$ model for supersonic flows. A comparison with experimental data is offered for a shock wave/turbulent boundary layer interaction case and two ramp injector mixing cases. Furthermore, a comparison is performed with empirical correlations on the basis of skin friction for flow over a flat plate and shear layer growth for a free shear layer. It is found that the maximum injectant mass fraction of some ramp injector cases is better predicted using the 1988 Wilcox model. On the other hand, the 2008 model performs better in simulating shock-boundary layer cases.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.109-116
• /
• 1997

The objective of this paper is to develop a computational model for the prediction of the pollutant spread from a mass source over a complex terrain. The model comprises a two-dimensional, steady state flow model and a concentration model which employs the results of the computed flow field. The computational model is applied to predict the spread of pollutants for Sanbon city, and the two cases have been compard with the results of SF/sub 6/ trace experiments.