• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.8
• /
• pp.677-684
• /
• 2007

In the Present Study, a 3-D viscous flow solver, based on unstructured hybrid meshses containing tetrahedra, prisms and pyramids, has been developed. A finite-volume discretization scheme is used for solving the compressible Navier-Stokes equations. A cell-vertex median dual volume is used for spatial discretization. The one-equation Spalart-Allmaras turbulence model has been adopted to evaluate the eddy viscosity. Validation were made by computing laminar and turbulent flows around a 3-D wing for steady flows and turbulent flows around an oscillating 3-D wing in harmonic motion for unsteady flows.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.51-57
• /
• 2001

The eddy viscosity turbulence models were applied to predict the flows through a cascade, and the prediction performances of turbulence models were assessed by comparing with the experimental results for a controlled diffusion(CD) compressor blade. The original $\kappa-\omega$ turbulence model and $\kappa-\omega$ shear stress transport(SST) turbulence model were used as two-equation turbulence model which were enhanced for a low Reynolds number flow and the Baldwin-Lomax turbulence model was used as algebraic turbulence model. Farve averaged Wavier-Stokes equations in a two-dimensional, curvilinear coordinate system were solved by an implicit, cell-centered finite-volume computer code. The turbulence quantities are obtained by lagging when the men flow equations have been updated. The numerical analysis was made to the flows of CD compressor blade in a cascade at three different incidence angles (40. 43.4. 46 degrees). We found the reversion in the prediction performance of original $\kappa-\omega$ turbulence model and $\kappa-\omega$ SST turbulence model when the incidence angie increased. And the algebraic Baldwin-Lomax turbulence model showed inferiority to two-equation turbulence models.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.6
• /
• pp.714-724
• /
• 1985

Two-phase(air-solid, air-liquid droplet) turbulent round jet has been analyzed numerically using two equation turbulence model. The mean motion of suspending particles in air has been treated as the secondary fluid with virtual density and eddy viscosity. In this paper, the local mean velocity of secondary fluid is not assumed to be the same as that of the primary one. Dissipation rate of turbulent kinetic energy which arises because the particles can not catch up with the turbulent fluctuations of the primary fluid has been modelled by using the concept of Kolmogorov's spectral energy transfer. Numerical computations were performed for flows with different volume fraction of the dispersed phase and the diameter of particle. Results show that the total rate of turbulent energy dissipation, turbulent intensities and spreading rate of jets are reduced by the increase of volume fraction of dispersed phase. However it does not show consistent tendency with increasing the particle diameter. This investigation also shows that presence of particles in the fluid modifies the structure of the primary fluid flow significantly. Predicted velocity profiles and turbulence properties qualitatively agree with available data.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.9
• /
• pp.1073-1084
• /
• 1999

The effects of secondary flow on the structure of a turbulent wake generated by a flat plate was investigated experimentally. The secondary flow was induced In a $90^{\circ}$ curved duct in which the flat plate wake generator was installed. The wake generator was installed in such a way that the wake velocity gradient exists in the span wise direction of the curved duct. Measurements were made in the plane containing the mean radius of curvature where pressure gradient and curvature effects were small compared with the secondary flow effect. All six components of the Reynolds stresses were measured in the curved duct. Turbulence intensities in the curved wake are higher than those in the straight wake due to an increase of the turbulent kinetic energy production by the secondary flow. In the inner wake region, shear stress and strain in the plane containing the velocity gradient of the wake show opposite signs with respect to each other, so that eddy viscosity Is negative in this region. This indicates that gradient-diffusion type turbulence models are not appropriate to simulate this type of flow.

• 한국신재생에너지학회:학술대회논문집
• /
• 2011.05a
• /
• pp.59.2-59.2
• /
• 2011

Although there are many activities on the construction of wind farm to produce amount of power from the wind, in practice power productions are not as much as its expected capabilities. This is because a lack of both the prediction of wind resources and the aerodynamic analysis on turbines with far wake effects. In far wake region, there are velocity deficits and increases of the turbulence intensity which lead to the power losses of the next turbine and the increases of dynamic loadings which could reduce system's life. The analysis on power losses and the increases of fatigue loadings in the wind farm is needed to prevent these unwanted consequences. Therefore, in this study velocity deficits have been predicted and aerodynamic analysis on turbines in the far wake is carried out from these velocity profiles. Ainslie's eddy viscosity wake model is adopted to determine a wake velocity and aerodynamic analysis on wind turbines is predicted by the numerical methods such as blade element momentum theory(BEMT) and vortex lattice method(VLM). The results show that velocity recovery is more rapid in the wake region with higher turbulence intensity. Since the velocity deficit is larger when the turbine has higher thrust coefficient, there is a huge aerodynamic power loss at the downstream turbine.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.33 no.4
• /
• pp.73-83
• /
• 1991

A numerical simulation of a 2-dimensional tidal flow in a shallow sea was performed using the frequency domain finite element method. In this study, to overcome the inherent problems of a time domain model which requires high eddy viscosity and small time steps to insure numerical stability, the harmonic function incorporated with the linearized function of governing equations was applied. Calculations were carried out using the developed tidal model(TIDE) in a rectangular channel of lOm(depth) X 4km (width) X 25km(length) under the condition of tidal waves entering the channel closed at one end for both with and without bottom friction damping. The predicted velocities and water levels at different points of the channel were in close agreement with less than 1 % error between the numerical and analytical solutions. The results showed that the characteristics of the tidal flow were greatly affected by the magnitude of tidal elevation forcing, and not by on surface friction, wind, or the linear bottom friction when the value was less than 0.01. For the optimum size of grid to obtain a consistent solution, the ratio between the length of the maximum grid and the tidal wave length should be less than 0.0018. It was concluded that the finite element tidal model(TIDE) developed in this study could handle the numerical simulation of tidal flows for more complex geometrical conditions.

• Nuclear Engineering and Technology
• /
• v.26 no.3
• /
• pp.345-354
• /
• 1994

Numerical Predictions including secondary flows have been Performed for fully developed turbulent single-phase rod bundle flows. The k-$\varepsilon$ turbulence model(two equation model) for the isotropic eddy viscosity, together with an algebraic stress model for generating secondary velocities, enabled the prediction of mean axial velocities, secondary velocities, and turbulent kinetic energy and turbulent stresses. Comparisons with experiment hate shown that the influence of secondary motion on mean flow and turbulence is dearly evident. The convective transport effects of secondary flow on the velocity field have been identified.

• Journal of Ocean Engineering and Technology
• /
• v.24 no.1
• /
• pp.116-122
• /
• 2010

The aim of this study is to examine the effects of shape and plane arrangement of submerged breakwaters on 3-D wave breaking characteristics over them. First, the numerical model, which is able to consider the flow through a porous medium with inertial, laminar, and turbulent resistance terms, i.e. simulate directly WAve Structure Seabed/Sandy beach interaction, and can determine the eddy viscosity with a LES turbulent model in a 3-Dimensional wave field (LES-WASS-3D), has been validated by a comparison with Goda's equation for breaking wave heights. And then, using the numerical results, the wave breaking points over the crest of submerged breakwaters have been examined in relation to the shape and plane arrangement of submerged breakwaters. Moreover, the wave height distribution and upper flow around submerged breakwaters have been also discussed, as well as the distribution of the wave breaking points over the beach.

• Journal of Ocean Engineering and Technology
• /
• v.21 no.6
• /
• pp.34-41
• /
• 2007

The present study investigates numerically damping characteristics of irregular waves passing over a permeable seabeds. At first, the numerical model, which is able to consider the flow through a porous medium with inertial, laminar and turbulent resistance terms and determine the eddy viscosity with LES turbulent model, is validated by comparing with existing experimental data. And then, the numerical test on irregular wave damping over a permeable seabeds is performed in case that wave and flume conditions are changed. It is revealed from the numerical results that the more porosity and mean grain are increased, the more wave damping is increased. Also, the effect of wave period on damping of irregular waves over a permeable seabed is discussed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.5
• /
• pp.1661-1670
• /
• 1996

Incompressible flow over a backward-facing step is computed by low Reynolds number turbulence models in order to compare with direct simulation results. In this study, selected low Reynolds number 1st and 2nd (Algebraic Stress Model : ASM) moment closure turbulence models are adopted and compared with each other. Each turbulence model predicts different flow characteristics, different re-attachment point, velocity profiles and Reynolds stress distribution etc. Results by .kappa.-.epsilon. turbulence models indicate that predicted re-attachment lengths are shorter than those by standard model. Turbulent intensity and eddy viscosity by low Reynolds number .kappa.-.epsilon. models are still greater than DNS results. The results by algebraic stress model (ASM) are more reasonable than those by .kappa.-.epsilon. models. The convective scheme is QUICK (Quadratic Upstream Interpolation for Convective Kinematics) and SIMPLE algorithm is adopted. Reynolds number based on step height and inlet free stream velocity is 5100.