• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.45 no.3
• /
• pp.233-240
• /
• 2017

The flowfield around transonic wing-body configuration was simulated using in-house CFD code and compared with the experimental data to understand the influence of several features of CFD(Computational Fluid Dynamics) ; grid dependency, turbulence models, spatial discretization, and viscosity. The wing-body configuration consists of a simple planform RAE Wing 'A' with an RAE 101 airfoil section and an axisymmetric body. The in-house CFD code is a compressible Euler/Navier-Stokes solver based on unstructured grid. For the turbulence model, the $k-{\omega}$ model, the Spalart-Allmaras model, and the $k-{\omega}$ SST model were applied. For the spatial discretization method, the central differencing scheme with Jameson's artificial viscosity and Roe's upwind differencing scheme were applied. The results calculated were generally in good agreement with experimental data. However, it was shown that the pressure distribution and shock-wave position were slightly affected by the turbulence models and the spatial discretization methods. It was known that the turbulent viscous effect should be considered in order to predict the accurate shock wave position.

• 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 Ocean Engineering and Technology
• /
• v.11 no.2
• /
• pp.100-106
• /
• 1997

Viscous flow around an actual ship is calculated by an use of RANS(Reynolds-averaged Navier-Stokes) solver. Reynolds stress is modelled by using k-$\varepsilon$ turbulence model and the law of wall is applied near the body. Body fitted coordinates are introduced for the treatment of the complex boundary of the ship hull form. The transformed equations in the computational domain are numerically solved by an employment of FVM(Finite Volume Method). SIMPLE(Semi-Implcit Pressure Linked Equation) method is adopted in the calculation of pressure and the solution of the disssssssscretized equation is obtained by the line-by-line method with the use of TDMA(Tri-Diagonal Matrix Algorithme). The subject ship model of actual calculation is 4,410 TEU class container carrier. For 4 geosim models the calculated viscous resistancce values are compared with the model test results and analyzed on their componentss. The resistance performance of an actual ship is predicted very resonably, so this mothod may be utilized as a design tool of hull form.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.6
• /
• pp.772-778
• /
• 2002

A numerical analysis of turbulent gas-particle two-phase flow is performed in conjunction with the experiments of Fackrell ＆ Robins and Raupach & Legg that considered ground-level source and/or elevated source flat plate flow. K-$\omega$ turbulence model is used in order to analyze fully turbulent flow field and the concentration equation with settling velocity is adopted for the concentration field. The model of Einstein and Chien is applied that couples the velocity field and the concentration field. Turbulent eddy viscosity is re-evaluated in this model. The present numerical results have good agreement between the simulation and the experimental data for the mean flow velocities and particle concentrations. While the previous study shows about 27% error in the vicinity of the source of particle concentration, the .present study allows about 14% error. A new turbulent gas-particle flow model developed by this study is able to cut down error by 13% at a near source.

• Journal of computational fluids engineering
• /
• v.14 no.3
• /
• pp.115-122
• /
• 2009

Numerical methods are applied to simulate the smoke behavior in a ventilated tunnel using large eddy simulation (LES) which is incorporated in FDS (Fire Dynamics Simulator) with proper combustion and radiation model. In this study, present numerical results are compared with data obtained from experiments on pool fires in a ventilated tunnel. The model tunnel is $182m(L){\times}5.4m(W){\times}2.4m(H)$. Two fire scenarios with different ventilation rates are considered with two different fire strengths. The present results are analyzed with those from LES without combustion and radiation model and from RANS ($\kappa-\epsilon$) model as well. Temperature distributions caused by fire in tunnel are compared with each other. It is found that thermal stratification and smoke back-layer can be predicted by FDS and the temperature predictions by FDS show better results than LES without combustion and radiation model. The FDS solver, however, failed to predict correct flow pattern when the high ventilation rate is considered in tunnel because of the defects in the tunnel-inlet turbulence and the near-wall turbulence.

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

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.1
• /
• pp.67-77
• /
• 2014

In this study, Explicit Algebraic Reynolds Stress Model (EARSM) which is based on the existing ${\kappa}-{\omega}$ model has been applied to the flow field analysis around ship hulls. Existing transport equations for the turbulent kinetic energy and the dissipation rate are used in almost the same form and anisotropy terms of Reynolds stresses are newly considered. The well-known KVLCC2 and KCS hull forms are selected as validation cases, which were also used in 2010 Workshop on CFD in Ship Hydrodynamics. In case of KVLCC2 double model, comparison of mean velocity distribution, turbulent kinetic energy, and Reynolds stresses near the propeller plane has been carried out and wave elevation and wave profiles have been additionally studied for KCS and KVLCC2 with free surface models. Some improved results for mean velocity distribution at the propeller plane have been obtained while there is little change in free surface wave profiles.

• Journal of the Korean Society of Visualization
• /
• v.20 no.3
• /
• pp.102-108
• /
• 2022

In the fan module of the intercooling refrigerator, a drain hole structure was designed for stable drainage of defrost water. However, the airflow passing through the drain hole can disturb flow features around the evaporator. Since this backflow leads to an increase in flow loss, the accurate experimental and numerical analyses are important to understand the flow characteristics around the fan module. Considering the complex geometry around the fan module, three different turbulence models (Standard k-ε model, SST k-ω model, Reynolds stress model) were used in computational fluid dynamics (CFD) analysis. According to the quantitative and qualitative comparison results, the Standard k-ε model was most suitable for the research object. High-accuracy results well match with the experiment result and overcome the limitation of the experiment setup. The method used in this study can be applied to a similar research object with an orifice outflow driven by a rotating blade.

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