• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.2
• /
• pp.43-55
• /
• 1995

The turbulent flow around the strut mounted on the plate is studied numerically. The main objective of this paper is to validate the numerical scheme by the comparison of the computed results with the measured one, especially, to investigate the applicability of the Baldwin-Lomax(B-L) model to the juncture flow. Computations are made by solving Reynolds-averaged wavier-Stokes equation with MAC method. The computed results are compared with experimental data of Dickinson, collected in the wind tunnel at DTRC. Comparisons show good agreements generally except at the region of wake and very near the juncture. Reynolds stress model seems to be required to improve the accuracy applicable to the juncture flow in spite of the many simplification of the turbulence modelling in B-L model.

• Journal of Mechanical Science and Technology
• /
• v.20 no.8
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split with van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge- Kutta method. The aspect ratios of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two- dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the first cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

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

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.22 no.5
• /
• pp.576-582
• /
• 2016

To cope with international regulations, such as Energy Efficiency Design Index (EEDI), Ship Energy Efficiency Management Plan (SEEMP) and so forth, and to enforce limitations on $CO_2$ emissions, green-ship technology to lower fuel consumption has been actively researched, and the development of an energy-saving device (ESD) is being pursued. In order to design an ESD for small and medium-sized domestic vessels, an analysis on flow characteristics has been performed in the present study. Through a model test and numerical analyses, the characteristics of flow around the stern bilge and bulb have been compared to improve wake quality and resistance performance. As a result of these comparisons and analyses, a vertical plate has been adopted,, as a new ESD. Design criteria for the proposed ESD are also suggested. By applying this new ESD, it is expected that the total resistance and average nominal wake can be reduced by 3.04 % and 18.8 %, respectively.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.2
• /
• pp.1-10
• /
• 2002

Hull-propeller-rudder interactions are studied by the iterative computational procedures. Hull effects on the propeller are reflected through the effective velocities computed by the vortex ring method which used the measured nominal wake as input data. A potential based panel method has been developed to solve the propeller-rudder interactions using the obtained effective velocities. Steady flow characteristics around the rudder surface can be obtained by computing the induced velocities on the rudder by the propeller and vice versa are computed by the iterative manner until the converged solutions are obtained. Flow characteristics around the propeller and the rudder are measured by Laser Doppler Velocimetry(L.D.V.) in large cavitation tunnel at Samsung Heavy industries. The gap flow model is adopted to solve the characteristics of the horn rudder. Numerical results are compared with the experimental values and the computed velocity fields and pressure distributions with rudder angle on the horn rudder surface show good agreement with measured ones in large cavitation tunnel.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.157-164
• /
• 2015

The present work performs three-dimensional flow calculations based on Reynolds Averaged Navier-Stokes (RANS) and Delayed Detached Eddy Simulation (DDES) to investigate the flow field of a transonic rotor (NASA Rotor 37) at near-stall condition. It is found that the DES approach is likely to predict well the complex flow characteristics such as secondary vortex or turbulent flow phenomenon than RANS approach, which is useful to describe the flow mechanism of a transonic compressor. Especially, the DES results show improvement of predicting the flow field in the wake region and the model captures reasonably well separated regions compared to the RANS model. Besides, it is discovered that the three-dimensional vortical flows after the vortex breakdown from the rotor tip region are widely distributed and its vortex structures are clearly present. Near the rotor leading edge, a part of the tip leakage flows in DES solution spill over into next passage of the blade owing to the separation vortex flow and the backflow is clearly seen around the trailing edge of rotor tip. Furthermore, the DES solution shows strong turbulent eddies especially in the rotor hub, rotor tip section and the downstream of rotor trailing edge compared to the RANS solution.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.235-238
• /
• 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 - w turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in cavities. Numerical method is upwind TVD scheme based on the flux vector split with the Van Leer limiters, and time accuracy is used explicit 4th stage Runge-Kutta scheme. Cavity flows are Comparison of two- and three-dimensional. The cavity has a L/D ratio of 3 for two-dimensional case. and same L/D and W/D ratio is 1 for three-dimensional case. The Mach and Reynolds numbers are held constant at 1.5 and 450000 respectively. For the three-dimensional case, the flow field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follow Rossiter's formula. On the other hand, the self-sustained oscillating flow transitions to a 'wake mode' for the two-dimensional simulation, with more violent fluctuations inside the cavity.

• Journal of Korean Society of Occupational and Environmental Hygiene
• /
• v.10 no.2
• /
• pp.140-149
• /
• 2000

In industrial field, there are several operations where a horizontal unidirectional airflow is used to control airborne contaminants. When a worker is immersed in a uniform freestream, a recirculating airflow can be created downstream of the worker by the phenomenon of boundary layer seperation. If the contaminant source and the breathing zone are within this near-wake region, high exposure may occur. The investigation for the effect of contaminant source location on worker exposure was performed by using CFD(Computational Fluid Dynamics). The airflow field was numerically calculated by assuming a steady flow and using the standard $k-{\varepsilon}$ turbulence model. As the results were compared with experimental data, the applicability of CFD was successfully verified. Subsequently, the breathing zone concentrations of the worker were predicted and compared with experimental data. The effects of contaminant density and turbulence intensity of freestream on worker exposure were evaluated.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.407-413
• /
• 2011

VIC (Vortex-In-Cell) method for viscous incompressible flow is presented to simulate the wake behind a modified NACA16 foil. With uniform rectangular grid, the velocity in field is calculated using streamfunction from vorticity field by solving the Poisson equation in which FFT(Fast Fourier Transform) is combined with 2nd order finite difference scheme. Here, LES(Large Eddy Simulation) with Smagorinsky model is applied for turbulence calculation. Effective viscosity is formulated using magnitude of strain tensor(or vorticity). Then the turbulent diffusion as well as viscous diffusion becomes particle strength exchange(PSE) with averaged eddy viscosity. The well-established panel method is combined to obtain the irrotational velocity and to apply the no-penetration boundary condition on the body panel. And wall diffusion is used for no-slip condition numerical results of turbulent stresses are compared with experimental results (Bourgoyne, 2003). Before comparing process, LES(Large Eddy Simulation) SGS(Subgrid scale) stress is transformed Reynolds averaged stress (Winckelmans, 2001).

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

Aerodynamic loads for a horizontal axis wind turbine of the National Renewable Energy Laboratory (NREL) Phase VI rotor in yawed condition were predicted by using the blade element momentum theorem. The classical blade element momentum theorem was complemented by several aerodynamic corrections and models including the Pitt and Peters' yaw correction, Buhl's wake correction, Prandtl's tip loss model, Du and Selig's three-dimensional (3-D) stall delay model, etc. Changes of the aerodynamic loads according to the azimuth angle acting on the span-wise location of the NREL Phase VI blade were compared with the experimental data with various yaw angles and inflow speeds. The computational flow chart for the classical blade element momentum theorem was adequately modified to accurately calculate the combined functions of additional corrections and models stated above. A successive under-relaxation technique was developed and applied to prevent possible failure during the iteration process. Changes of the angle of attack according to the azimuth angle at the specified radial location of the blade were also obtained. The proposed numerical procedure was verified, and the predicted data of aerodynamic loads for the NREL Phase VI rotor bears an extremely close resemblance to those of the experimental data.