• Wind and Structures
• /
• v.17 no.4
• /
• pp.361-377
• /
• 2013

Wall functions have been widely used in computational fluid dynamics (CFD) simulations and can save significant computational costs compared to other near-wall flow treatment strategies. However, most of the existing wall functions were based on the asymptotic characteristics of near-wall flow quantities, which are inapplicable in complex and non-equilibrium flows. A modified wall function is thus derived in this study based on flow over a plate at zero-pressure gradient, instead of on the basis of asymptotic formulations. Turbulent kinetic energy generation ($G_P$), dissipation rate (${\varepsilon}$) and shear stress (${\tau}_{\omega}$) are composed together as the near-wall expressions. Performances of the modified wall function combined with the nonlinear realizable k-${\varepsilon}$ turbulence model are investigated in homogeneous equilibrium atmosphere boundary layer (ABL) and flow around a 6 m cube. The computational results and associated comparisons to available full-scale measurements show a clear improvement over the standard wall function, especially in reproducing the boundary layer flow. It is demonstrated through the two case studies that the modified wall function is indeed adaptive and can yield accurate prediction results, in spite of its simplicity.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.12
• /
• pp.2056-2061
• /
• 2001

Non-proportional loading tests of Nylon 66 at room temperature exhibit path dependent behavior and plasticity-relaxation interactions. The uniaxial formulation of the viscoplasticity theory based on overstress (VBO), which has been used to reproduce the nonlinear strain rate sensitivity, relaxation, significant recovery and cyclic softening behaviors of Nylon 66, is extended to three-dimensions to predict the response in strain-controlled, comer-path tests. VBO consists of a flow law that is easily written for either the stress or the strain as the independent variable. The flow law depends on the overstress, the difference between the stress and the equilibrium stress that is a state variable in VBO. The evolution law of the equilibrium stress in turn contains two additional state variables, the kinematic stress and the isotropic stress. The simulations show that the constitutive model is competent at modeling the deformation behavior of Nylon 66 and other solid polymers.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.404-408
• /
• 2005

The Reynolds strss equation with turbulent energy-length scales was simplified in the nearly homogeneous turbulent equilibrium flow and a modified Reynolds stress model was proposed. Tn the model proposed in the present study, Reynolds stresses can be expressed in the form of algebraic equation, so that the turbulent stresses and related quantities are calculated through relatively simple procedures. The model predicted well the turbulent shear stresses of homogeneous flow in local equilibrium state obtained from experimental results published earlier Constants used In the model was determined universally and its validity was discussed briefly.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.19 no.4
• /
• pp.483-491
• /
• 2016

A simple aerodynamic calculation program for high Mach number flow is developed by combining the modified Newtonian method with Tannehill's curve fits for the thermodynamic properties of air in equilibrium state. Aerodynamic characteristics for a parabolic nose are predicted and compared with CFD(Computational Fluid Dynamics) analysis results. Comparison shows good agreements, and the developed program is expected to be a practical tool for slender body aerodynamic calculation for high Mach number flow.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.13-17
• /
• 2006

RoeM and AUSMPW+ schemes are two of the most accurate and efficient schemes which are recently developed for the analysis of single phase gas dynamics. In this paper, we developed two-phase versions of these schemes for the analysis of gas-liquid large density ratio two-phase flow. We adopt homogeneous equilibrium model (HEM) using mass fraction to describe different two phases. In the Eulerian-Eulerian framework, HEM assumes dynamic and thermal equilibrium of the two phases in the same computational mesh. From the mixture equation of state (EOS), we derived new shock-discontinuity sensing term (SDST), which is commonly used in RoeM and AUSMPW+ for the stable numerical flux calculation. The proposed two-phase versions of RoeM and AUSMPW+ schemes are applied on several air-water two-phase test problems. In spite of the large discrepancy of material properties such as density, enthalpy, and speed of sound, the numerical results show that both schemes provide very satisfactory solutions.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1903-1908
• /
• 2004

A preconditioned numerical method for gas-liquid two-phase flows is applied to solve cavitating flow. The present method employs a finite-difference dual time-stepping integration procedure and the MUSCLTVD scheme. A homogeneous equilibrium cavitation model is used. The present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation, large density changes and incompressible flow characteristics at low Mach number. Some internal flows such as convergent-divergent nozzles are computed using this method. Comparisons of predicted and experimental results are provided and discussed.

• Journal of Korea Water Resources Association
• /
• v.45 no.2
• /
• pp.217-227
• /
• 2012

This paper presents a quasi-steady model for numerical simulations of reservoir sedimentation and reservoir flushing. The quasi-steady model is based on the assumption that the flow is steady with time-dependent stream morphology change. This is reasonable because stream morphology changes over a long period, while the flow changes rapidly. The proposed model is first applied to two laboratory experiments for reservoir sedimentation. The channel is shown to be adjusted to new sediment supply at the upstream by changing both the flow depth and slope. Simulated water surface and bed profiles compare favorably to measured data. The model is also applied to reservoir flushing. Good agreement between simulated and measured data is not obtained due to time variation of outflow generated to facilitate the flushing in the experiment. Finally, relationships for equilibrium flow depth and bed slope are proposed and tested through numerical experiments.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.6
• /
• pp.23-29
• /
• 2003

Transonic two-phase flow of Saturated humid air, in which relative humidity is 100%, with various condensation processes around thin airfoils is investigated. The study uses an extended transonic small-disturbance(TSD) model of Rusak and Lee [11, 12] which includes effects of heat addition to the flow due to condensation. Two possible limit types of condensation processes are considered. In the nonequilibrium and homogeneous process, the condensate mass fraction is calculated according to classical nucleation and droplet growth rate models. In the equilibrium process, the condensate mass fraction is calculated by assuming an isentropic process. The flow and condensation equations are solved numerical1y by iterative computations. Results under same upstream conditions describe the flow structure, field of condensate, and pressure distribution on airfoil's surfaces. It is found that flow characteristics, such as position and strength of shock waves and airfoil’s pressure distribution, are different for the two condensation processes. Yet, in each case, heat addition as a result of condensation causes significant changes in flow behavior and affects the aerodynamic performance of airfoils.

• Nuclear Engineering and Technology
• /
• v.31 no.6
• /
• pp.572-585
• /
• 1999

The characteristics of two-phase flow instability in a vertical boiling channel connecting with an unheated riser are investigated through the linear stability analysis model. Various two-phase flow models, including thermal non-equilibrium effects, are taken into account for establishing a physical model in the time domain. A classical approach to the frequency response method is adopted for the stability analysis by employing the D-partition method. The adequacy of the linear model is verified by evaluating experimental data at high quality conditions. It reveals that the flow-pattern-dependent drift velocity model enhances the prediction accuracy while the homogeneous equilibrium model shows the most conservative predictions. The characteristics of density wave oscillations under low-power and low-quality conditions are investigated by devising a simple model which accounts for the gravitational and frictional pressure losses along the channel. The necessary conditions for the occurrences of type-I instability and flow excursion are deduced from the one-dimensional D-partition analysis. The parametric effects of some design variables on low quality oscillations are also investigated.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.7 no.1
• /
• pp.57-65
• /
• 2003

As a compressed water is rapidly expanded through a nozzle, two-phase flow of vapor and liquid is formed in the nozzle due to the flash evaporation. In the present study, critical flow of two-phase fluids is analysized using an Isentropic-Homogeneous-Equilibrium model and a Leung model. Calculation results show that the choke of the two-phase flow can be two different types of continuous and discontinuous chokings. For the stagnation pressure below 10 Mpa it is found that the continuous choking, which is similar to the choking phenomenon of single-phase gas flow, is possible only when the degree of subcooling is less than 10K.