• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.1
• /
• pp.72-78
• /
• 2008

In this study, advanced (fluid-structure interaction (FSI)) analysis system has been developed in order to predict turbine cascade performance with blade deformation effect due to aerodynamic loads. Intereference effects due to the relative movement of the rotor cascade with respect to the stator cascade are also considered. Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation k-ω SST turbulence models are solved to accurately predict fluid dynamic loads considering flow separation effects. A fully implicit time marching scheme based on the (coupled Newmark time-integration method) with high artificial damping is efficiently used to compute the complex fluid-structure interaction problem. Predicted aerodynamic performance considering structural deformation effect of the blade shows somewhat different results compared to the case of rigid blade model. Cascade performance evaluations for different elastic axis positions are importantly presented and its aeroelastic effects are investigated.

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

In the present study, the Reynolds-averaged Navier-Stokes equations, together with the equations of the $k-{\varepsilon}$ model of turbulence, were solved numerically in a general body-fitted coordinate system for three-dimensional turbulent flows around the six basic shapes of the magnetically levitated train(MAGLEV). The numerical computations were conducted on the MAGLEV model configurations to provide information on shapes of this type very near the elevated track at a constant Reynolds number of $1.48{\times}10^{6}$ based on the body length. The coordinate system was generated by numerically solving a set of Poisson equations. The convective transport equations were discretized using the finite-analytic scheme which employed analytic solutions of the locally-linearized equations. A time marching algorithm was employed to enable future extensions to be made to handle unsteady and fully-elliptic problems. The pressure-velocity coupling was treated with the SIMPLER-algorithm. Of particular interests were wall effect by the elevated track on the aerodynamic forces and flow characteristics of the six models calculated. The results indicated that the half-circle configuration with extended sides and with smooth curvature of sides was desirable because of the low aerodynamic forces and pitching moment. And it was found that the separation bubble was occured at wake region in near the elevated track.

• Wind and Structures
• /
• v.15 no.6
• /
• pp.481-494
• /
• 2012

The greenhouse type metal structures are increasingly used in modern construction of livestock farms because they are less laborious to construct and they provide a more favorable microclimate for the growth of animals compared to conventional livestock structures. A key stress factor for metal structures is the wind. The external pressure coefficient ($c_{pe}$) is used for the calculation of the wind effect on the structures. A high pressure coefficient value leads to an increase of the construction weight and subsequently to an increase in the construction cost. The EC1 in conjunction with EN 13031-1:2001, which is specialized for greenhouses, gives values for this coefficient. This value must satisfy two requirements: the safety of the structure and a reduced construction cost. In this paper, the Navier - Stokes and continuity equations are solved numerically with the finite element method (Galerkin Method) in order to simulate the two dimensional, incompressible, viscous air flow over the vaulted roofs of single span and twin-span with eaves livestock greenhouses' structures, with a height of 4.5 meters and with length of span of 9.6 and 14 m. The simulation was carried out in a wind tunnel. The numerical results of pressure coefficients, as well as, the distribution of them are presented and compared with data from Eurocodes for wind actions (EC1, EN 13031-1:2001). The results of the numerical experiment were close to the values given by the Eurocodes mainly on the leeward area of the roof while on the windward area a further segmentation is suggested.

• Journal of Mechanical Science and Technology
• /
• v.19 no.11
• /
• pp.2032-2039
• /
• 2005

The lattice Boltzman method (LBM) and the finite difference-based lattice Boltzmann method (FDLBM) are quite recent approaches for simulating fluid flow, which have been proven as valid and efficient tools in a variety of complex flow problems. They are considered attractive alternatives to conventional finite-difference schemes because they recover the Navier-Stokes equations and are computationally more stable, and easily parallelizable. However, most models of the LBM or FDLBM are for incompressible fluids because of the simplicity of the structure of the model. Although some models for compressible thermal fluids have been introduced, these models are for monatomic gases, and suffer from the instability in calculations. A lattice BGK model based on a finite difference scheme with an internal degree of freedom is employed and it is shown that a diatomic gas such as air is successfully simulated. In this research we present a 2-dimensional edge tone to predict the frequency characteristics of discrete oscillations of a jet-edge feedback cycle by the FDLBM in which any specific heat ratio $\gamma$ can be chosen freely. The jet is chosen long enough in order to guarantee the parabolic velocity profile of a jet at the outlet, and the edge is of an angle of $\alpha$=23$^{o}$. At a stand-off distance w, the edge is inserted along the centerline of the jet, and a sinuous instability wave with real frequency is assumed to be created in the vicinity of the nozzle exit and to propagate towards the downstream. We have succeeded in capturing very small pressure fluctuations resulting from periodic oscillation of the jet around the edge.

• Wind and Structures
• /
• v.17 no.6
• /
• pp.647-670
• /
• 2013

The paper presents a numerical approach to study of fluid flow characteristics and to predict performance of wind turbines. The numerical model is based on Finite-volume method (FVM) discretization of unsteady Reynolds-averaged Navier-Stokes (URANS) equations. The movement of turbine blades is modeled using moving mesh technique. The turbulence is modeled using commonly used turbulence models: Renormalization Group (RNG) k-${\varepsilon}$ turbulence model and the standard k-${\varepsilon}$ and k-${\omega}$ turbulence models. The model is validated with the experimental data over a large range of tip-speed to wind ratio (TSR) and blade pitch angles. In order to demonstrate the use of numerical method as a tool for designing wind turbines, two dimensional (2-D) and three-dimensional (3-D) simulations are carried out to study the flow through a small scale Darrieus type H-rotor Vertical Axis Wind Turbine (VAWT). The flows predictions are used to determine the performance of the turbine. The turbine consists of 3-symmetrical NACA0022 blades. A number of simulations are performed for a range of approaching angles and wind speeds. This numerical study highlights the concerns with the self-starting capabilities of the present VAWT turbine. However results also indicate that self-starting capabilities of the turbine can be increased when the mounted angle of attack of the blades is increased. The 2-D simulations using the presented model can successfully be used at preliminary stage of turbine design to compare performance of the turbine for different design and operating parameters, whereas 3-D studies are preferred for the final design.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.12 no.1
• /
• pp.9-14
• /
• 2009

A numerical investigation was performed based on the Reynolds-Averaged Navier-Stokes(RANS) equations for the two-dimensional unsteady flow around a vertical axis turbine(VAT) with three or four blades. VAT is one of the promising devices for tidal current energy conversion. The geometry of the turbine blade was $NACA65_3$-018 airfoil, for which CFD analysis using Fluent was carried out at several angles of attack and the results were compared with the corresponding experimental data for validation and calibration. Then CFD simulations were carried out for the whole vertical axis turbine with a two-dimensional setup. The CFD simulation demonstrated the usefulness of the method to study the typical unsteady flows around VATs and the results showed that the optimum turbine efficiency could be achieved for carefully selected combinations of the number of blade and Tip-Speed Ratio(TSR).

• Journal of Energy Engineering
• /
• v.21 no.3
• /
• pp.237-242
• /
• 2012

This study presents numerical solutions of the two-dimensional Navier-Stokes equations for supersonic unsteady flow in a convergent-divergent nozzle with a isolator. The TVD scheme in generalized coordinates is employed in order to calculate the moving shock waves caused by thermal choking. We discuss on transient characteristics, unstart phenomena, fluctuations of specific thrust caused by thermal choking and effects of isolator. The adverse pressure gradient caused by heat addition brings about separation of the wall boundary layers and formation of the oblique shock wave that proceed to upstream. The proceeding speed of the oblique shock wave to upstream direction for the convergent-divergent nozzle with isolator is lower than that for the nozzle without isolator.

• Tribology and Lubricants
• /
• v.33 no.5
• /
• pp.207-213
• /
• 2017

In sliding bearings, viscous friction due to high shear acting on the bearing surface raises the oil temperature. One of the mechanisms responsible for generating the load-carrying capacity in parallel surfaces is known as the viscosity wedge effect. In this paper, we investigate the effect of film-temperature boundary conditions on the thermohydrodynamic (THD) lubrication of parallel slider bearings. For this purpose, the continuity equation, Navier-Stokes equation, and the energy equation with temperature-viscosity-density relations are numerically analyzed using the commercial computational fluid dynamics (CFD) code FLUENT. Two different film-temperature boundary conditions are adopted to investigate the pressure generation mechanism. The temperature and viscosity distributions in the film thickness and flow directions were obtained, and the factors related to the pressure generation in the equation of motion were examined in detail. It was confirmed that the temperature gradients in the film and flow directions contribute heavily to the thermal wedge effect, due to which parallel slider bearing can not only support a considerable load but also reduce the frictional force, and its effect is significantly changed with the film-temperature boundary conditions. The present results can be used as basic data for THD analysis of surface-textured sliding bearings; however, further studies on various film-temperature boundary conditions are required.

• International Journal of Aeronautical and Space Sciences
• /
• v.17 no.2
• /
• pp.260-267
• /
• 2016

The design process of blended inlet body (BIB) for the preliminary design of a near-space high supersonic unmanned aerial vehicle (HSUAV) is presented. The mass flow rate and cowl area of inlet at a design point are obtained according to the cruise condition of the HSUAV. A mixed-compression axisymmetric supersonic inlet section with a fixed geometry reasonably matching the high supersonic cruise state is created by using the inviscid theory of aerodynamics. The inlet section is optimized and used as a baseline section for the BIB design. Three BIB concepts for the HSUAV are proposed, and their internal aerodynamic characteristics of inlet are evaluated using Euler computational fluid dynamics (Euler CFD) solver. The preferred concept is identified, in which the straight leading edge of the baseline HSUAV configuration is modified into the convex leading edge to accommodate the inlet and meet the requirements of the cowl area to capture the sufficient air flow. The total recovery of inlet for the preferred BIB concept and the aerodynamic characteristics of the modified HSUAV configuration are verified using Navier-Stokes computational fluid dynamics (NS CFD) solver. The validation indicates that the preferred BIB concept can meet both the requirements of the inlet and aerodynamic performance of the HSUAV.