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

Aerodynamic characteristics over Micro Aerial Vehicle(MAV) in low Reynolds number regime are numerically studied using 3-D unsteady, incompressible Navier-Stokes flow solver with single partitioning method for multi-block grid. For more efficient computation of unsteady flows, this flow solver is parallel-implemented with MPl(Message Passing Interface) programming method. Firstly, MAV wing with not complex geometry is considered and then, we analyze aerodynamic characteristics over full MAV configuration varying the angle of attack. Present computational results show a better agreement with the experimental data by MACDL(Micro Aerodynamic Control and Design Lab.), Seoul National University. We can also find the conceptually designed MAV by MACDL has the static stability.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.5
• /
• pp.731-740
• /
• 2001

A numerical study on the quantitative analogy between the fully developed turbulent flow in a straight square duct rotating about an axis perpendicular to that of the duct and that in a stationary curved duct of square cross-section is carried out. In order to clarify the similarity of the two flows, dimensionless parameters K(sub)TR=Re(sup)1/4/√Ro and Rossby number, Ro, in a rotating straight duct flow were used as a set corresponding to K(sub)TC=Re(sup)1/4/√λ and curvature ratio, λ, in a stationary curved duct flow so that they have the same dynamical meaning as those of the fully developed laminar flows. For the large values of Ro or λ, it is shown that the flow field satisfies the asymptotic invariance property, that is, there are strong quantitative similarities between the two flows such as flow patterns and friction factors for the same values of K(sub)TR and K(sub)TC.

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

• Transactions of the Korean Society of Automotive Engineers
• /
• v.2 no.1
• /
• pp.26-37
• /
• 1994

In this study, flow of a model intake port/valve system is analyzed by using low Reynolds number $k-{\varepsilon}$ model. Discharge coefficient was obtained from computational results for the various cases of valve lifts. Discharge coefficient becomes maximum when the valve lift is 20mm, and does not increase or decrease in proportional to valve lift. Most of pressure drop and production of turbulent kinetic energy occur at the edge points of the valve and the valve seat Thus, in order to improve discharge coefficient, rounding of edge points in valve and valve seat is recommended. As valve lift is increased, the velocity of the intake jet in the valve passage decreases, and the direction of the jet is more inclined toward the valve seat.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.8
• /
• pp.2650-2669
• /
• 1996

A new low Reynolds number nonlinear second moment turbulence closure was introduced to analyze a square sectioned 180.deg. bend flow. Inclusion of nonlinear return to isotropy term and cubic mean pressure strain term has brought out a marked improvement in the level of agreement with measured velocity profiles. Optimization of present closure was performed by comparison of computed velocity profiles with the experimental ones with variation of nonlinear return to isotropy term and quadratic and cubic pressure-strain model. Progressive vortex breakdown due to the interaction of primary and secondary flows was well captured by using the optimized second moment turbulence closure.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1971-1978
• /
• 1992

A turbulent flow in a conical diffuser with total divergence angle of 8.deg. was numerically studied. The low Reynolds number k-.epsilon. model(Launder-Sharma model) was adopted to simulate the turbulence. The continuity and time averaged Navier-Stokes equations in a nonorthogonal coordinate system were solved by a finite volume method based on the fully elliptic formulation. The low Reynolds number k-.epsilon. model reasonably simulates the pressure recovery and the mean velocity components. However, there are also considerable discrepancies between predicted and measured shear stress distribution on the wall and turbulent kinetic energy distributions. It is necessary to investigate the flow structure at the entry of the diffuser, numerically as well as experimentally.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.63-72
• /
• 1998

The numerical simulations of flowfield and pollutant dispersion over two-dimensional hills of various shapes are described. The Reynolds-averaged Wavier-Stokes equations and concentration diffusion equation based on the gradient diffusion theory have been applied to the atmospheric shear flow over the bell-shaped hills which are basic components of the complex terrain. The flow characteristics such as velocity profiles of the geophysical boundary layer, speed-up phenomena, mean pollutant concentration profiles are compared with experimental data to validate the present numerical procedure and it has been found that the present numerical results agree well with experiments and other numerical data. It has been also found that the distributions of ground level concentration are strongly influenced by the source location and height.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.37-42
• /
• 1995

In the present study, an improved version of 4-equation low-Reynolds-number 4-equation model is proposed. The equations of the temperature variance ($k_{\theta}$) and its dissipation rate(${\varepsilon}_{\theta}$) are solved, in concert with the equations of the turbulent kinetic energy(k) and its dissipation rate(${\varepsilon}$). In the present model, the near-wall effect and the non-equilibrium effect are fully taken into consideration. The validation of the model is then applied to the turbulent flow behind a backward-facing step and the flow over a blunt body. The predicted results of the present model are compared and evaluated with the relevant experiments.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.6
• /
• pp.733-739
• /
• 1999

A modified low-Re $k-\varepsilon$ model is used for the calculation of drag-reducing turbulent flow by polymer injection in a pipe. With the viscoelastic model, molecular viscosity in the definition of turbulent viscosity is related to elongations viscosity of the solution to account for the effects of drag reduction. Finite volume method is used for the discretization, and power-law scheme is used as a numerical scheme. Computed dimensionless velocity profiles are in good agreements with the experimental data in case of low drag reductions. However, in case of high drag reductions, they deviate largely from the measurements in the central zone of the flow field.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.2
• /
• pp.1-9
• /
• 1999

Applicability of the pressure gradient method which is formulated based on pressure gradient is verified against turbulent flow analysis. In the pressure gradient method, pressure gradient instead of pressure itself is obtained using continuity constraint. Since correct pressure gradient is found only when mass conservation is satisfied, pressure gradient method can reflect physics of flow field properly The pressure gradient method is formulated with semi-staggered grid system which locates each primitive variables on the same grid point but evaluates pressure gradient in-between. This grid system ensures easy programming and reflection of correct physics in analysis. For verifying applicability of this method, the pressure gradient method is applied to turbulent flow analysis with low Reynolds number $\kappa$-$\varepsilon$ model. Turbulent flows include fully developed channel flow, backward-facing step flow, and conical diffuser flow. Prediction results show that the pressure gradient method can be applied to turbulent flow analysis. However, the pressure gradient method requires somewhat long computation time. Proper way to find optimum under-relaxation factor, $\gamma$, is also need to be developed.