• 대기
• /
• 제27권4호
• /
• pp.445-453
• /
• 2017

This study investigates turbulent flow around a square cylinder mounted on a flat surface at high Reynolds number using a large-eddy simulation (LES) model, particularly focusing on vortex streets behind the square cylinder. Total 9 simulation cases with different inflow wind directions, inflow wind speeds, and cylinder widths in the x- and y-directions are considered to examine the effects of inflow wind direction, speed, and cylinder widths on turbulent flow and vortex streets. In the control case, the inflow wind parallel to the x-direction has a maximum speed of $5m\;s^{-1}$ and the width and height of the cylinder are 50 m and 200 m, respectively. In all cases, down-drafts in front of the cylinder and updrafts, wakes, and vortex streets behind the cylinder appear. Low-speed flow below the cylinder height and high-speed flow above it are mixed behind the cylinder, resulting in strong negative vertical turbulent momentum flux at the boundary. Accordingly, the magnitude of the vertical turbulent momentum flux is the largest near the cylinder top. In the case of an inflow wind direction of $45^{\circ}$, the height of the boundary is lower than in other cases. As the inflow wind speed increases, the magnitude of the peak in the vertical profile of mean turbulent momentum flux increases due to the increase in speed difference between the low-speed and high-speed flows. As the cylinder width in the y-direction increases, the height of the boundary increases due to the enhanced updrafts near the top of the cylinder. In addition, the magnitude of the peak of the mean turbulent momentum flux increases because the low-speed flow region expands. Spectral analysis shows that the non-dimensional vortex generation frequency in the control case is 0.2 and that the cylinder width in the y-direction and the inflow wind direction affect the non-dimensional vortex generation frequency. The non-dimensional vortex generation frequency increases as the projected width of the cylinder normal to the inflow direction increases.

• 한국전산유체공학회지
• /
• 제21권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.

• 대한기계학회논문집A
• /
• 제29권5호
• /
• pp.707-715
• /
• 2005

In this work, Scalar Passive code in Lattice Boltzmann Method is employed to simulate two-phase flow of low Reynolds number in a micro-channel. The mixing characteristics in a micro-channel is a function of Peclet number. The mixing length increases with the Peclet number. It is found that with the inclusion of static elements at the channel, rapid mixing of two liquids can be achieved, as shown by the results of computer simulations. The enhancement in mixing performance is thought to be caused by the generation of eddies and by lateral velocity component when the mixture flows past static elements. The results indicate that the size of static element has more effect on the mixing than the number of static element.

• 대한기계학회논문집B
• /
• 제23권3호
• /
• pp.374-387
• /
• 1999

Two-dimensional turbulent fluid flow and solidification were investigated in a continuous casting process of a steel slab with electromagnetic field. The electromagnetic field was described by the Maxwell equations. The enthalpy-porosity relation was employed to suppress the velocity within a mushy region. A revised low-Reynolds number $k-{\varepsilon}$ turbulence model was used to consider the turbulent effects. It is shown that the temperature gradient in the casting direction in the case with EMBR becomes very weak compared to that of the case without EMBR. The results also show that the velocity profiles of the case with solidification are quite different from those of the case without solidification.

• International Journal of Aeronautical and Space Sciences
• /
• 제16권2호
• /
• pp.165-176
• /
• 2015

A numerical analysis was conducted to investigate the inlet buzz and combustion oscillation in an axisymmetric ramjet engine with wedge-type flame holders. The physical model of concern includes the entire engine flow path, extending from the leading edge of the inlet center-body through the exhaust nozzle. The theoretical formulation is based on the Farve-averaged conservation equations of mass, momentum, energy, and species concentration, and accommodates finite-rate chemical kinetics and variable thermo-physical properties. Turbulence closure is achieved using a combined scheme comprising of a low-Reynolds number k-${\varepsilon}$ two-equation model and Sarkar's compressible turbulence model. Detailed flow phenomena such as inlet flow aerodynamics, flame evolution, and acoustic excitation as well as their interactions, are investigated. Mechanisms responsible for driving the inlet buzz are identified and quantified for the engine operating at subcritical conditions.

• 대한기계학회논문집B
• /
• 제22권5호
• /
• pp.573-583
• /
• 1998

Effects of the oscillating incoming flow on vortex shedding and lock-on behind a square cylinder are investigated using numerical simulations at a Reynolds number of 100. Vortex shedding occurred at low forcing frequencies of the incoming flow similar to the natural vortex shedding. As the forcing frequency further increases, the shedding frequency decreases to the half of the forcing freqnency. For a sufficiently large frequency, vortex shedding returns to the natural vortex shedding irrespective of the forcing amplitude. Also, the lock-on region becomes wider with higher forcing amplitudes. The phase diagram between the drag and lift shows a simple periodic behavior in the lock-on region, while a complicated periodic phase relation is observed when there is no lock-on.

• Nuclear Engineering and Technology
• /
• 제28권4호
• /
• pp.405-414
• /
• 1996

In this paper, the unsteady 2-dimensional turbulent flow model for thermal stratification in a pressurizer surge line of PWR plant is proposed to numerically investigate the heat transfer and flow characteristics. The turbulence model is adapted to the low Reynolds number K-$\varepsilon$ model (Davidson model). The dimensionless governing equations are solved by using the SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm. The results are compared with simulated experimental results of TEMR Test. The time-dependent temperature profiles in the fluid and pipe nil are shown with the thermal stratification occurring in the horizontal section of the pipe. The corresponding thermal stresses are also presented. The numerical result for thermal stratification by the outsurge during heatup operation of PWR shows that the maximum dimensionless temperature difference is about 0.83 between hot and cold sections of pipe well and the maximum thermal stress is calculated about 322MPa at the dimensionless time 28.5 under given conditions.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.32-36
• /
• 2010

This paper dealt with the flow simulation for the optimum designed propeller for Micro Aerial Vehicle, using a commercial CFD program(FLUENT). The propeller was modeled by the Multiple Reference Frame(MRF) method. For the validation of the computational method, the flow field analysis results for the propeller were compared with the flow analysis results, which are using Xfoil, for the optimum design, and with the wind tunnel data of a similar propeller model. By these validation processes, the reliability of MRF method was confirmed.

• 대한기계학회논문집
• /
• 제18권4호
• /
• pp.967-976
• /
• 1994

A tensor invariant model equation for the turbulent energy dissipation rate is proposed in the present study, which is able to simulate secondary straining effects such as curvature effects without the introduction of additional empirical input. The source term in this model has a combined form of the generation term due to the mean vorticity with the conventional one due to the mean strain rate. An extended low-Reynolds-number $k-\epsilon$ turbulence model involving this new model equation is tested for a turbulent Coutte flow between coaxial cylinders with inner cylinder rotated, which is a well defined example of curved flows. The predicted results indicate that the present model works much better for this flow, compared with previous models.

• 한국분무공학회지
• /
• 제19권4호
• /
• pp.167-173
• /
• 2014

The spray characteristics of Gas-Centered Swirl Coaxial Injector was investigated that there were different characteristics with or without gas flow. As gas flow was accelerated, the momentum of gas was transferred to the momentum of liquid in the low liquid Reynolds number. Therefore, the axial velocity of liquid was increased and the measured value was smaller than without gas flow. However, in the high momentum flux ratio, the momentum transfer hardly occurred and the results had constant values. As the recess length was increased, the mixing area of gas and liquid also was increased, the results were decreased.