• 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.

• 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.

• Journal of the Korean Society of Visualization
• /
• v.19 no.1
• /
• pp.88-93
• /
• 2021

For a flow normal to a circular disk, the flow separation occurs from the edge of the disk and the flow recirculation zone exists behind the disk. Many existing studies conducted simulations of flow normal to a circular disk under low Reynolds numbers. Some studies performed LES or DES simulations under high Reynolds numbers. However, comparative study for different RANS models for high Reynolds numbers is very limited. This study presents numerical simulations of a flow normal to a circular disk using Realizable k-ε model and SST k-ω model. The recirculation bubble length and drag coefficient were compared with the experimental data. The SST k-ω model showed the excellent predictions for the recirculation bubble length and drag coefficient.

• Wind and Structures
• /
• v.17 no.1
• /
• pp.87-105
• /
• 2013

Modeling an equilibrium atmospheric boundary layer (ABL) in an empty computational domain has routinely been performed with the k-${\varepsilon}$ turbulence model. However, the research objects of structural wind engineering are bluff bodies, and the SST k-${\omega}$ turbulence model is more widely used in the numerical simulation of flow around bluff bodies than the k-${\varepsilon}$ turbulence model. Therefore, to simulate an equilibrium ABL based on the SST k-${\omega}$ turbulence model, the inlet profiles of the mean wind speed U, turbulence kinetic energy k, and specific dissipation rate ${\omega}$ are proposed, and the source terms for the U, k and ${\omega}$ are derived by satisfying their corresponding transport equations. Based on the proposed inlet profiles, numerical comparative studies with and without considering the source terms are carried out in an empty computational domain, and an actual numerical simulation with a trapezoidal hill is further conducted. It shows that when the source terms are considered, the profiles of U, k and ${\omega}$ are all maintained well along the empty computational domain and the accuracy of the actual numerical simulation is greatly improved. The present study could provide a new methodology for modeling the equilibrium ABL problem and for further CFD simulations with practical value.