• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.1
• /
• pp.8-18
• /
• 2015

A numerical simulation for a nonslender BWB UCAV configuration with a rounded leading edge and span of 1.0 m was performed to analyze its aerodynamic characteristics. Numerical results were compared with experimental data obtained at a free stream velocity of 50 m/s and at angles of attack from -4 to $26^{\circ}$. The Reynolds number, based on the mean chord length, is $1.25{\times}106$. 3D multi-block hexahedral grids are used to guarantee good grid quality and to efficiently resolve the boundary layer. Menter's shear stress transport model and two transition models (${\gamma}-Re_{\theta}$ model and ${\gamma}$ model) were used to assess the effect of the laminar/turbulent transition on the flow characteristics. Aerodynamic coefficients, such as drag, lift, and the pitching moment, were compared with experimental data. Drag and lift coefficients of the UCAV were predicted well while the pitching moment coefficient was underpredicted at high angles of attack and influenced strongly by the selected turbulent models. After assessing the pressure distribution, skin friction lines and velocity field around UCAV configuration, it was found that the transition effect should be considered in the prediction of aerodynamic characteristics of vortical flow fields.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.7
• /
• pp.879-886
• /
• 2004

A numerical optimization has been carried out to determine the shape of the three-dimensional channel with oblique ribs attached on both walls to enhance turbulent heat transfer. The response surface based optimization is used as an optimization technique with Reynolds-averaged Navier-Stokes analysis of fluid flow and heat transfer. Shear stress transport (SST) turbulence model is used as a turbulence closure. Numerical results fur heat transfer rate show good agreements with experimental data. four dimensionless variables such as, rib pitch-to-rib height ratio, rib height-to-channel height ratio, streamwise rib distance on opposite wall to rib pitch ratio, and the attack angle of the rib are chosen as design variables. The objective function is defined as a linear combination of heat-transfer and friction-loss related coefficients with a weighting factor. D-optimal method is used to determine the training points as a means of design of experiment. Sensitivity of the objective parameters to each design variable has been analyzed. And, optimal values of the design variables have been obtained in a range of the weighting factor.

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

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.05a
• /
• pp.190-193
• /
• 2010

The turbulent mixing and combustion of a shear coaxial injector under supercritical pressures have been theoretically/numerically investigated. Turbulent numerical model is based on large eddy simulation with real-fluid transport and thermodynamics over the entire pressure; Soave modification of Redlich-Kwong equation of state, Chung's model for viscosity/conductivity, and Fuller's theorem for diffusivity to take account Takahashi's compressible effect. The results are compared with previous researcher's. The large-scale vortices shedding from the outer rim into the recirculation region to react with gaseous oxygen was investigated.

• The Bulletin of The Korean Astronomical Society
• /
• v.36 no.1
• /
• pp.34.1-34.1
• /
• 2011

Combining a detailed non-grey radiative transfer computation with the three dimensional hydrodynamics, we investigate a reliable numerical scheme for turbulent convection in the solar surface. The solar photosphere is the extremely turbulent region composed of partly ionized compressible gases in high temperature. Especially, the super adiabatic layer (SAL) near the solar photosphere is the shallow transition region where the energy transport varies steeply from convection to radiation. In order to describe physical processes accurately, a detailed treatment of radiative transfer should be considered as well as the high resolution computation of fluid dynamics. For a direct computation of radiation fields, the Accelerated Lambda Iteration (ALI) methods have been applied to hydrodynamical medium, incorporating the Opacity Distribution Function (ODF) as a realistic schemes for non-grey problems. Computational domain is the rectangular box of dimensions $42{\times}3Mn$ with the resolution of $1202{\times}190$ meshed grids, which covers several granules horizontally and 8 ~ 9 pressure scale heights vertically. During several convective turn-over times, the 3-D snapshots have been compiled with a second order accuracy. In addition, our radiation-hydrodynamical computation has been compared with the classical approximations such as grey atmospheres and Eddington approximation.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.04a
• /
• pp.103-108
• /
• 2011

The turbulent mixing of a kerosene/liquid oxygen coaxial swirl injector under supercritical pressures have been numerically investigated. Kerosene surrogate models are proposed for the kerosene thermodynamic properties. Turbulent numerical model is based on LES(Large Eddy Simulation) with real-fluid transport and thermodynamics over the entire pressure range; Soave modification of Redlich-Kwong equation of state, Chung's model for viscosity/conductivity, and Fuller's theorem for diffusivity to take account Takahashi's compressible effect. The effect of operating pressure on thermodynamic properties and mixing dynamics inside an injector and a combustion chamber are investigated. Power spectral densities of pressure fluctuations in the injector under various chamber pressure are analyzed.

• Journal of the Korean Society of Visualization
• /
• v.22 no.2
• /
• pp.63-73
• /
• 2024

The injection of polymer can significantly reduce drag, particularly in the turbulent flow region where the mutual interaction between the polymer and turbulent vortices occurs. In this study, Taylor-Couette flow of PEO-in-water solutions with a rotating inner cylinder was analyzed. Despite the shear-thinning behaviour of PEO-in-water solutions being well-documented, for a given range of shear rates their viscosity remains nearly constant. By varying the polymer concentration, we analyzed the torque evolution of different solutions followed by the viscoelasticity effects of the polymer on the interphase transition points. The torque was analyzed using a dimensionless torque scaling method, which allows for the assessment of the fluid's momentum transport capabilities. It was observed that for low concentrations of PEO, the flow behaviour exhibited only minor differences in comparison to that of water, the Newtonian fluid. However, once the PEO concentration exceeded the polymer overlap concentration, the flow behaviour was significantly altered.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.2
• /
• pp.169-178
• /
• 2011

The dispersion of a pollutant in a turbulent boundary layer has been described in this study by using a two-dimensional lattice Boltzmann method (LBM) and the Smagorinsky sub-grid-scale (SGS) model. The scalar transport equation corresponding to the pollutant concentration is adopted; the pollutant is considered to be in a continuous phase. The pollutant source is classified as ground-level source (GLS) and elevated-point source (ES). Air velocity and particle concentration profile for the pollutant are compared with the respective results and profiles obtained in the experiments of Fackrell and Robins (1982) and Raupach and Legg (1983). The numerical results obtained in this study, i.e., the simulation and the experimental data for the mean flow velocity profiles and the pollutant concentration profiles, are in good agreement with each other.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.3
• /
• pp.529-561
• /
• 2014

The two-phase turbulent flow past an interface-piercing circular cylinder is studied using a high-fidelity orthogonal curvilinear grid solver with a Lagrangian dynamic subgrid-scale model for large-eddy simulation and a coupled level set and volume of fluid method for air-water interface tracking. The simulations cover the sub-critical and critical and post critical regimes of the Reynolds and sub and super-critical Froude numbers in order to investigate the effect of both dimensionless parameters on the flow. Significant changes in flow features near the air-water interface were observed as the Reynolds number was increased from the sub-critical to the critical regime. The interface makes the separation point near the interface much delayed for all Reynolds numbers. The separation region at intermediate depths is remarkably reduced for the critical Reynolds number regime. The deep flow resembles the single-phase turbulent flow past a circular cylinder, but includes the effect of the free-surface and the limited span length for sub-critical Reynolds numbers. At different Froude numbers, the air-water interface exhibits significantly changed structures, including breaking bow waves with splashes and bubbles at high Froude numbers. Instantaneous and mean flow features such as interface structures, vortex shedding, Reynolds stresses, and vorticity transport are also analyzed. The results are compared with reference experimental data available in the literature. The deep flow is also compared with the single-phase turbulent flow past a circular cylinder in the similar ranges of Reynolds numbers. Discussion is provided concerning the limitations of the current simulations and available experimental data along with future research.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.1
• /
• pp.23-33
• /
• 2020

In this study, the effect of surface roughness distribution and its influence on the inflight icing code was investigated. Previous numerical studies focused on the magnitude of surface roughness, and the effects were only addressed in terms of changes in thermal boundary layers with fully turbulent assumption. In addition, the empirical formula was used to take account the turbulent transition due to surface roughness, which was regarded as reducing the accuracy of ice shape prediction. Therefore, in this study, the turbulent transition model based on the two-equation turbulence model was applied to consider the effects of surface roughness. In order to consider the effect of surface roughness, the transport equation for roughness amplification parameter was applied, and the surface roughness distribution model was implemented to consider the physical properties. For validation, the surface roughness, convective heat transfer coefficient, and ice shape were compared with experimental results and other numerical methodology. As a result, it was confirmed that the excessive prediction of the heat transfer coefficient at the leading edge and the ice horn shape at the bottom of the airfoil were improved accordingly.