• Journal of computational fluids engineering
• /
• v.10 no.4 s.31
• /
• pp.51-58
• /
• 2005

The unsteady supersonic flow over two- and three-dimensional cavities has been analyzed by the integration of unsteady Reynolds-Averaged Navier-Stokes(RANS) with the k-$\omega$ turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in the cavity. An explicit 4th order Runge-Kutta scheme and an upwind TVD scheme based on the flux vector split with the van Leer limiters are used for time and space discritizations, respectively. The cavity has a L/D ratio of 3 for two-dimensional case, and same L/D and W/D ratio of I for three-dimensional case. The Mach and Reynolds numbers are 1.5 and 450000 respectively. In the three-dimensional flow, the field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follows Rossiter's formula. In the two-dimensional simulation, the self-sustained oscillating flow has more violent fluctuation inside the cavity. The primary fluctuating frequencies of two- and three- dimensional flow agree very well with the 2nd mode of Rossiter's frequency. In the three-dimensional flow, the 1st mode of frequency could be seen.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.24 no.5
• /
• pp.219-232
• /
• 2020

The slip-weakening model developed by Ohnaka and Yamashita is extended over the breakdown zone by equating the scaling relationships for the breakdown zone and the whole rupture area. For the extension, the study uses the relationship between rupture velocity and radiation efficiency, which was derived in the theory of linear elastic fracture mechanics, and the definition of f_max given in the specific barrier model proposed by Papageorgiou and Aki. The results clearly show that the extended scaling relationship is governed by the ratio of rupture velocity to S wave velocity, and the velocity ratio can be determined by the ratio of characteristic frequencies of a Fourier amplitude spectrum, which are corner frequency, f_c, and source-controlled cut-off frequency, f_max, or vice versa. The derived relationship is tested by using the characteristic frequencies extracted from previous studies of more than 130 shallow crustal events (focal depth less than 25 km, M_W 3.0~7.5) that occurred in Japan. Under the assumption of a dynamic similarity, the rupture velocity estimated from f_max/f_c and the modified integral timescale give quite similar scale-dependence of the rupture area to that given by Kanamori and Anderson. Also, the results for large earthquakes show good agreement to the values from a kinematic inversion in previous studies. The test results also indicate the unavailability of the spectral self-similarity proposed by Aki because of the scale-dependent rupture velocity and the rupture velocity-dependent f_max/f_c; however, the results do support the local similarity asserted by Ohnaka. It is also remarkable that the relationship between the rupture velocity and f_max/f_c is quite similar to Kolmogorov's hypothesis on a similarity in the theory of isotropic turbulence.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.8 no.4
• /
• pp.561-575
• /
• 1996

The interaction of turbulent mixed convection and surface radiation in a three-dimensional channel with the heated blocks is analyzed numerically. Two blocks are maintained at high temperature and the other bottom and horizontal walls are insulated. S-4 method is employed to calculate the effect of the radiative heat transfer. The low Reynolds number k-$\varepsilon$ model proposed by Launder and Sharma is used to estimate the turbulent influence on the heat transfer enhancement. From above modeling, the effects of various channel specifications on the flow and heat transfer characteristics are investigated. The variables used for the present study are Reynolds number, block spacing, the channel height spacing for block and the emissivity. Average Nusselt numbers along the block surfaces are correlated and presented in terms of Reynolds number, emissivity and dimensionless geometric parameters. For the range of conditions in this study, average Nusselt numbers along the block surfaces are strongly influenced by the Reynolds numbers and channel height spacing for block but weakly influenced by the block spacing and the emissivity of the adiabatic walls.

• 유체기계공업학회:학술대회논문집
• /
• 2002.12a
• /
• pp.152-157
• /
• 2002

A cross-flow fan is widely used on many industrial fields: a blower for the general industry, mining industry, automobile and home appliances. The design point of the cross-flow fan is generally chosen by based on the region within low static pressure and high flow rate. It relatively makes high dynamic pressure at low speed because a working fluid passes through an impeller blade twice. However, it has low static pressure efficiency between $30\%$ and $40\%$ because of relative high impact loss. Recently, in the air-conditioning systems, the operating behaviors at the off-design points are highly regarded to broaden the application area for various air-cooling loads. Especially, at the low flow rate, there exists a rapid pressure head reduction, a noise increase and an irregular flow against a rearguider as a scroll of centrifugal fan. Numerical analyses are carried out for cross-flow fan including the impeller, the rearguider and the stabilizer. Numerical domains are discretized by hexahedral cells. Three-dimensional, unsteady governing equations are solved using FVM, SIMPLE algorithm, sliding grid system and standard k-$\epsilon$ turbulence model.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.195-200
• /
• 2010

Numerical investigation was conducted to study the effects of after-body configurations and nozzle lip on the PIFS(Plume Induced Flow Separation) and eat flux to the base face. Two dimensional and axi-symmetric non-equilibrium Navier-Stoke's solver with $k-{\omega}$ SST turbulence model was used to solve the launching vehicle type configuration with propulsive jet. The experimental result of Robert J. McGhee was compared with our computational results for code validation. Three types of the after-body configurations (Straight, Boat-tail, Flare type) were simulated for this study. And the nozzle lip effect was studies using the three types of base configurations same simulation conditions. As a result of numerical investigations, higher pressure ratio condition and boat-tail after-body configuration caused severe PIFS phenomenon but the flare type after-body configuration and low pressure ratio suppressed PIFS. Flare type after-body configuration and low pressure ratio case reduced heat flux to base face. The nozzle lip dispersed the heat flux widely along the base face and the nozzle lip.

• The KSFM Journal of Fluid Machinery
• /
• v.7 no.4 s.25
• /
• pp.7-15
• /
• 2004

The design point of the cross-flow fan is generally based on the region within low static pressure and high flow rate. It relatively makes high dynamic pressure at low speed because a working fluid passes through an impeller blade twice. However, it has low static pressure efficiency between $30\%$ and $40\%$ because of relative high impact loss. The purpose of this study is to research the reciprocal relation among each parameter. Experiments and numerical analyses are conducted on effects of a stabilizer and a rearguider on performance analysis of a cross-flow fan. Two-dimensional, unsteady governing equations are solved using FVM, PISO algorithm, sliding grid system and standard $k-{\epsilon}$ turbulence model. Experiments are also carried out to estimate the performance of the modeled cross-flow fan. It is clarified that the rearguider of Archimedes type has excellent results for the most part.

• Journal of the Korean Solar Energy Society
• /
• v.35 no.1
• /
• pp.69-80
• /
• 2015

A series of numerical analyses was conducted to predict the thermal performance of a solar air heater depending on the hole configuration and geometry in the absorber plate. The planar dimensions of the prototype were 1 m (W) by 1.6 m (H), and the maximum air flow considered was $187m^3/h$. It was considered that protruding holes with a triangular opening in the absorber plate would invoke turbulence in the air flow to enhance the convection heat transfer. Six different hole configurations were investigated and compared with each other, while the hole opening height was considered as a design variable. Three-dimensional transient analyses were performed with a commercial software package on the airflow and heat transfer in the model. The numerical results were analyzed and compared from the view point of the outlet air temperature and its time response to derive the optimal hole pattern and hole opening height.

• Journal of computational fluids engineering
• /
• v.21 no.2
• /
• pp.70-80
• /
• 2016

The comparison of two commercial codes(FLUENT and STAR-CCM+) and an open-source code(OpenFOAM) are carried out for the aerodynamic analysis of flight vehicles at low speeds. Tailless blended-wing-body UCAV, main wing and propeller of HALE UAV(EAV-3) are chosen as geometries for the investigation. Using the same mesh, incompressible flow simulations are carried out and the results from three different codes are compared. In the linear region, the maximum difference of lift and drag coefficients of UCAV are found to be less than 2% and 5 counts, respectively and shows good agreement with wind tunnel test data. In a stall region, however, the reliability of RANS simulation is found to become poor and the uncertainty according to code also increases. The effect of turbulence models and meshes generated from different tools are also examined. The transition model yields better results in terms of drag which are much closer to the test data. The pitching moment is confirmed to be sensitive to the existence and the location of transition. For the case of EAV-3 wing, the difference of results with ${\kappa}-{\omega}$ SST model is increased when Reynolds number becomes low. The results for the propeller show good agreement within 1% difference of thrust. The reliability and uncertainty of three codes is found to be reasonable for the purpose of engineering use. However, the physical validity and reliability of results seem to be carefully examined when ${\kappa}-{\omega}$ SST model is used for aerodynamic simulation at low speeds or low Reynolds number conditions.

• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.1
• /
• pp.30-45
• /
• 2013

The present work focuses on the unsteady aerodynamics and aeroelastic properties of a small-medium sized wind-turbine blade operating under ideal conditions. A tapered/twisted blade representative of commercial blades used in an experiment setup at the National Renewable Energy Laboratory is considered. The aerodynamic loads are computed using Computational Fluid Dynamics (CFD) techniques. For this purpose, FLUENT$^{(R)}$, a commercial finite-volume code that solves the Navier-Stokes and the Reynolds-Averaged Navier-Stokes (RANS) equations, is used. Turbulence effects in the 2D simulations are modeled using the Wilcox k-w model for validation of the CFD approach. For the 3D aerodynamic simulations, in a first approximation, and considering that the intent is to present a methodology and workflow philosophy more than highly accurate turbulent simulations, the unsteady laminar Navier-Stokes equations were used to determine the unsteady loads acting on the blades. Five different blade pitch angles were considered and their aerodynamic performance compared. The structural dynamics of the flexible wind-turbine blade undergoing significant elastic displacements has been described by a nonlinear flap-lag-torsion slender-beam differential model. The aerodynamic quasi-steady forcing terms needed for the aeroelastic governing equations have been predicted through a strip-theory based on a simple 2D model, and the pertinent aerodynamic coefficients and the distribution over the blade span of the induced velocity derived using CFD. The resulting unsteady hub loads are achieved by a first space integration of the aeroelastic equations by applying the Galerkin's approach and by a time integration using a harmonic balance scheme. Comparison among two- and three- dimensional computations for the unsteady aerodynamic load, the flap, lag and torsional deflections, forces and moments are presented in the paper. Results, discussions and pertinent conclusions are outlined.

• Journal of computational fluids engineering
• /
• v.16 no.1
• /
• pp.53-59
• /
• 2011

In this study, steady and unsteady aerodynamic analyses of a huge rocket configuration have been conducted in a transonic flow region. The launch vehicle structural response are coupled with the transonic flow state transitions at the nose of the payload fairing. Before performing the coupled fluid-structure transonic aeroealstic simulations transonic aerodynamic characteristics are investigated for the pitching motions of the rocket at finite angle-of-attack. An unsteady CFD analysis method with a moving grid technique based on the Reynolds-averaged Navier-Stokes equations with the k-w SST transition turbulence model is applied to accurately predict the transonic loads of the rocket at pitching motion. It is shown that the fluctuating amplitude of the lateral aerodynamic loads imposed on the rocket due to the pitching motion can be significantly increased in the transonic flow region.