• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.63-66
• /
• 2005

A computational fluid dynamics(CFD) analysis has been made for fully developed turbulent flow in a triangular bare rod bundle with a pitch to diameter ratio (P/D) of 1.123. The nonlinear turbulence models predicted the turbulence-driven secondary flow in the triangular subchannel. The nonlinear quadratic κ-ε models by Speziale[1] and Myong-Kasagi[2] predicted turbulence structure in the rod bundle fairly well. The nonlinear quadratic and cubic k-ε models by Shih et al.[3] and Craft et al.[4] showed somewhat weaker anisotropic turbulence. The differential Reynolds stress model by Launder et al.[5} appeared to over predict the turbulence anisotropy in the rod bundle.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.7
• /
• pp.911-918
• /
• 1997

The turbulent viscous wake flows behind a single airfoil, two-dimensional stationary blade row and three-dimensional rotating blade row were calculated, and the numerical results were compared with experimental ones. The numerical technique was based on the SIMPLE algorithm using three turbulent closure models, standard k-.epsilon. model(WFM), low Reynolds number k-.epsilon. model(LRN) and Reynolds stress model (RSM). In the case of a single airfoil, WFM, LRN and RSM presented fairly good velocity distributions in the wake compared with experimental data. In the case of the stationary blade row, LRN and RSM presented better results than WFM for wake velocity distribution, and especially LRN showed best results among these three turbulent models. In the case of the rotating blade row, WFM and LRN showed fairly good agreement with experimental data of the three-dimensional velocity component distributions in the range from hub to mid span region. LRN was also superior to WFM in accuracy of prediction for the wake velocity distribution as same with the cases of a airfoil and the stationary blade row.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.12 no.6
• /
• pp.704-710
• /
• 2009

In this study, a comparison study of flutter analysis for the AGARD 445.6 wing with wind turnnel test data has been conducted in the subsonic, transonic and supersonic flow regions. Nonlinear aeroelastic using FSIPRO3D which is a generalized user-friendly fluid-structure analyses have been conducted for a 3D wing configuration considering shockwave and turbulent viscosity effects. The developed fluid-structure coupled analysis system is applied for aeroelastic computations combining computational structure dynamics(CSD), finite element method(FEM) and computations fluid dynamics(CFD) in the time domain. MSC/NASTRAN is used for the vibration analysis of a wing model, and then the result is applied to the FSIPRO3D module. the results for dynamic aeroelastic response using different turbulent models are presented for several Mach numbers. Calculated flutter boundary are compared with the wind-tunnel experimental and the results show very good agreements.

• Coupled systems mechanics
• /
• v.7 no.1
• /
• pp.95-109
• /
• 2018

The relevance of turbulent mixing in estuarine numerical models for stratified two-layer shallow water flows is analysed in this paper. A one-dimensional numerical model was developed for this purpose by extending an immiscible two-layer model with an additional source term, which accounts for turbulent mixing effects, namely the entrainment of fluid from the lower to the upper layer. The entrainment rate is quantified by an empirical equation as a function of the bulk Richardson number. A finite volume method based on an approximated Roe solver was used to solve the governing coupled system of partial differential equations. A comparison of numerical results with and without entrainment is presented to illustrate the influence of entrainment on both the salt-water intrusion length and lower layer dynamics. Furthermore, one example is given to demonstrate how entrainment terms may help to stabilize the numerical scheme and prevent a possible loss of hyperbolicity. Finally, the model with entrainment is validated by comparing the numerical results to field measurements.

• Wind and Structures
• /
• v.19 no.6
• /
• pp.687-708
• /
• 2014

The performance of the $k-{\varepsilon}$ and $k-{\omega}$ two-equation turbulence models was investigated in computational simulations of the neutrally stratified atmospheric boundary layer developing above various terrain types. This was achieved by using a proposed methodology that mimics the experimental setup in the boundary layer wind tunnel and accounts for a decrease in turbulence parameters with height, as observed in the atmosphere. An important feature of this approach is pressure regulation along the computational domain that is additionally supported by the nearly constant turbulent kinetic energy to Reynolds shear stress ratio at all heights. In addition to the mean velocity and turbulent kinetic energy commonly simulated in previous relevant studies, this approach focuses on the appropriate prediction of Reynolds shear stress as well. The computational results agree very well with experimental results. In particular, the difference between the calculated and measured mean velocity, turbulent kinetic energy and Reynolds shear stress profiles is less than ${\pm}10%$ in most parts of the computational domain.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.16 no.9 s.114
• /
• pp.937-948
• /
• 2006

In this study, a fluid/structure coupled analysis system for simulating complex flow-induced vibration (FIV) phenomenon of cascades has been developed. The flow is modeled using Euler and Wavier-Stokes equations with different turbulent models. The fluid domains are modeled using the unstructured grid system with dynamic deformations due to the motion of structural boundary. The Spalart-Allmaras (S-A) and the SST ${\kappa}-{\omega}$ turbulent models are used to predict the transonic turbulent flows. A fully implicit time marching scheme based on the Newmark direct integration method is used in order to solve the coupled governing equations for viscous flow-induced vibration phenomena. For the purpose of validation for the developed FIV analysis system, comparison results for computational analyses of steady and unsteady aerodynamics and flutter analyses are presented in the transonic flow region. In addition, flow-induced vibration analyses for the isolated cascade and multi-blades cascade models have been conducted to show the physical fluid-structure interaction effects in the time domain.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.05a
• /
• pp.793-802
• /
• 2006

In this study, a fluid/structure coupled analysis system for simulating complex flow-induced vibration (FIV) phenomenon of cascades has been developed. The flow is modeled using Euler and Wavier-Stokes equations with different turbulent models. The fluid domains are modeled using the unstructured grid system with dynamic deformations due to the motion of structural boundary. The Spalart-Allmaras (S-A) and the SST ${\kappa}-{\omega}$ turbulent models are used to predict the transonic turbulent flows. A fully implicit time marching scheme based on the Newmark direct integration method is used in order to solve the coupled governing equations for viscous flow-induced vibration phenomena. For the purpose of validation for the developed FIV analysis system, comparison results for computational analyses of steady and unsteady aerodynamics and flutter analyses are presented in the transonic flow region. In addition, flow-induced vibration analyses for the isolated cascade and multi-blades cascade models have been conducted to show the physical fluid-structure interaction effects in the time domain.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.17 no.1
• /
• pp.61-69
• /
• 2013

Numerical simulations have been performed for assessment of compressibility correction of two-equation turbulent models suitable for large scale separation flows perturbed by a pintle strokes. Two-equation turbulence models, the low Reynolds k-${\varepsilon}$ and the k-${\omega}$ SST models with or without compressibility correction proposed by Wilcox and Sarkar are evaluated. The detail flow structures are observed and static pressures along nozzle wall are compared with experimental results. Mach disk location and pressure recovery profiles in flow separation region are noticeably distinct between turbulent models of k-${\varepsilon}$ and k-${\omega}$ SST. The compressible effect corrections to those models improve resolving of separation flow behaviors. The compressibility corrections to k-${\varepsilon}$ model have provided very comparable results with test data.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.6 no.4
• /
• pp.331-340
• /
• 1982

Numerical analysis is made on the turbulent heat transfer with suspension of solid particles in circular tube with constant heat flux. The mean motion of suspending particles in mixture is treated as the secondary gas flow with virtual density and viscosity. Our modeling of turbulent transport phenomena of suspension flow is based on this assumption and conventional mixing length theory. This paper gives the evidence that the mixing length models can be extended to close the governing equations for two phase turbulent flow with solid boundary at a first order level. Results on Nusselt numbers obtained by analytical treatments are compared with available experimental data and discussed. They suggest that the most important parameters of two phase turbulent heat transfer phenomena are relative particle diameter to pipe diameter, gas-solid loading ratio, and specific heat of suspending material.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.8
• /
• pp.680-687
• /
• 2007

The paper presents a numerical simulation of flow of a turbulent boundary layer, representing a typical wind environment and matching a series of wind tunnel observations. The simulations are carried out at a Reynolds number of 20,000, based on the velocity U at a pseudo-height h, and large enough that the flow be effectively Reynolds number independent. Some wall models are proposed for the LES(Large Eddy Simulation) of the turbulent boundary layer over a rough surface. The Jenson number, $J=h/z_0$, based on the roughness length $z_0$, is 600 to match the wind tunnel data. The computational mesh is uniform with a spacing of h/32, as this aids rapid convergence of the multigrid solver, and the governing equations are discretised using second order finite differences within a parallel multiblock environment. The results presented include the comparison between wind tunnel measurements and LES computations of the turbulent boundary layer over rough surface.