• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.5
• /
• pp.1300-1307
• /
• 1995

The effect of thermal stratification on the turbulent dispersion from a fine cylindrical heat source was experimentally examined in a wind tunnel with and without a strong temperature gradient. A 0.5 mm dia. nichrome wire was used as a line heat source. Turbulent intensities, r.m.s. value of temperature and convective heat fluxes were measured by using a hot-wire and cold-wire combination probe. The results show that the peack value and the spread of the vertical turbulent intensity for the stratified case are far lower than those in the neutral case, which indicates that the stable temperature gradient suppresses the vertical velocity component. All of the third order moments including heat fluxes measured in the stable condition have very small values than those of the neutral case. This nature suggests that the decrease of scalar fluctuations in the stably stratified flow is mainly due to the suppression ofthe turbulent diffusion processes by the stable stratification. A simple gradient model with a composite timescale which has a simple weighted algebraic mean between dynamic and thermal time scale yields reasonably good numerical values in comparison with the experimental data.

• Journal of computational fluids engineering
• /
• v.5 no.2
• /
• pp.38-46
• /
• 2000

It has been indicated that the rotor/stator interaction has distinct causes of unsteadiness, such as the viscous vortex shedding, wake/stator interaction and potential rotor/stator interaction. In this paper, the mechanism of unsteady potential interaction in one stage axial compressor is numerically investigated for blade row ratio 1:1 and 2:3 at design point and for blade row ratio 2:3 at off-design point in two-dimensional view point. The numerical technique used is the upwind scheme of Van-Leer's Flux Vector Splitting(FVS) and Cubic spline interpolation is applied on zonal interface. In this study the flow unsteadiness due to potential interaction are found to be larger in blade row ratio 2:3 than in 1:1. The total pressure rise in blade row ratio 2:3 is closer to the real value in design point than that in 1:1. The change of unsteady pressure amplitude according to the variation of stator exit pressure is very small.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.7
• /
• pp.495-505
• /
• 2009

A two-dimensional circular cylinder freely falling in a channel has been simulated by using immersed boundary - lattice Boltzmann method (IB-LBM) in order to analyze the characteristics of motion originated by the interaction between the fluid flow and the cylinder. The wide range of the solid/fluid density ratio has been considered to identify the effect of the solid/fluid density ratio on the motion characteristics such as the falling time, the transverse force and the trajectory in the streamwise and transverse directions. In addition, the effect of the gap between the cylinder and the wall on the motion of a two-dimensional freely falling circular cylinder has been revealed by taking into account a various range of the gap size. As the cylinder is close to the wall at the initial dropping position, vortex shedding in the wake occurs early since the shear flow formed in the spacing between the cylinder and the wall drives flow instabilities from the initial stage of freely falling. In order to consider the characteristics of transverse motion of the cylinder in the initial stage of freely falling, quantitative information about the cylinder motion variables such as the transverse force, trajectory and settling time has been investigate.

• Journal of the Society of Naval Architects of Korea
• /
• v.55 no.1
• /
• pp.1-8
• /
• 2018

A cavitator plays an important role to generate the supercavity. Most previous numerical and experimental studies have been focused on the presence of cavitators alone. However, the body behind the cavitator causes a change in the wake flow and hence it affects generation and growth of the supercavity. In this paper, we present a boundary elementary method based on a potential flow analysis, and calculate characteristics of the supercavity formation depending on the change of the body shape of two-dimensional submerged objects. Various parameters such as cone angle of the cavitator, length of the forehead and diameter of the body are considered. The results show that the longer the forepart length, the longer the cavity is created under the same conditions, and also the change in the diameter of the body is the most influential factor on the growth of the supercavity. As a result, we suggest that it is necessary to carefully consider the influence of the body shape during the initial design stage of the supercavitating underwater vehicle.

• Proceedings of the KIEE Conference
• /
• 2015.07a
• /
• pp.247-248
• /
• 2015

Wind turbines causes instabilities on the grid as their penetration increase. To mitigate harmful effects from wind turbines, transmission system operator(TSO) set up some requirements to obligate for wind generation operator for grid connection. So wind farm management system(WFMS) has important role to follow requirement from TSO, WFMS calculates available real power by considering wake effects, and dispatches real power order to each wind turbine in wind farm to optimize for decreasing load fatigue. To verify operation of WFMS, real-time simulator(RTS) is necessary. This paper deals with RTS configuration to verify WFMS operation. RTS includes wind farm model and power flow code. Normally, wind farm equivalent simple model makes wind turbines in wind farm to one wind turbine mode which cannot verify power flow in wind farm and WFMS operation. Thus, this paper makes wind farm model using simple wind turbine model with transfer function. Matlab is used for make power flow code and wind farm model to impose RTS and those model is certified by PSCAD/EMTDC.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.3
• /
• pp.267-279
• /
• 1998

A numerical study was made of flow behind a circular cylinder in a uniform flow, where the cylinder was rotationally oscillated in time. The temporal behavior of vortex formation was scrutinized over broad ranges of the two externally specified parameters, i.e., the dimensionless rotary oscillating frequency (.110.leq. $S_{f}$.leq..220) and the maximum angular amplitude of rotation (.theta.$_{max}$=15 deg., 30 deg. and 60 deg.). The Reynolds number (Re= $U_{{\inf}D}$.nu.) was fixed at Re=110. A fractional-step method was utilized to solve the Navier-Stokes equations with a generalized coordinate system. The main emphasis was placed on the initial vortex formations by varying $S_{f}$ and .theta.$_{max}$. Instantaneous streamlines and pressure distributions were displayed to show the vortex formation patterns. The vortex formation modes and relevant phase changes were characterized by measuring the lift coefficient ( $C_{L}$) and the time of negative maximum $C_{L}$( $t_{-C}$$_{Lmax}$) with variable forcing conditions.s.tions.s.s.s.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.3
• /
• pp.336-345
• /
• 1998

The substantial loss behind axial flow rotor was generated by wake, various vortices in the hub region and the leakage vortex in the tip region. Particularly, the leakage vortex formed near blade tip was one of the main causes of the reduction of performance, the generation of noise and the aerodynamic vibration in rotor downstream. In this study, the three-dimensional flowfields in an axial flow rotor for various tip clearances were calculated, and the numerical results were compared with the experimental ones. The numerical technique was based on SIMPLE algorithm using standard k-.epsilon. model (WFM). Through calculations, the effects of the tip clearance on the overall performance of rotor and the loss distributions, and the increase in the displacement, momentum, and blade-force-deficit thickness of the casing wall boundary layer were investigated. The mass-averaged flow variables behind rotor agreed well with the experimental results. The presence of the tip leakage vortex behind rotor was described well. Although the loci of leakage vortex by calculation showed some differences compared with the experimental results, its behavior for various tip clearances was clarified by examining the loci of vortex center.

• Journal of the Korean Society for Aviation and Aeronautics
• /
• v.21 no.1
• /
• pp.45-50
• /
• 2013

In this study, lumped-vortex element method and thin airfoil theory were used to analyze aerodynamic characteristics of airfoils with relative motion that had camber lines of NACA $44{\times}{\times}$ airfoil in 2-dimensional unsteady incompressible potential flow. Velocity disturbance due to airfoil was calculated by lumped-vortex element model and force distribution on airfoil by unsteady Bernoulli's equation. Variables in relative motion were considered the period p, the amplitude of flapping $A_f$ and pitching $A_p$, and the phase difference between flapping and pitching ${\phi}_p$ and the angle of attack ${\alpha}$. Due to movement of an airfoil, dag was induced in 2-dimensional unsteady incompressible potential flow. The numerical results show that the aerodynamic characteristics of the airfoil with flapping and pitching at the same time are illustrated. Especially the mean lift coefficient became smaller, but drag coefficient became larger.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.2
• /
• pp.899-913
• /
• 2019

In order to provide a complementary perspective to the effects of the maneuvering motions on the unsteady propeller performance, the numerical simulation of the flow field of the hull-rudder- propeller system is performed by Unsteady Reynolds-averaged Naiver-Stokes (URANS) method. Firstly, the flow fields around the submarine model without the presence of propeller in straight ahead motion and the steady diving maneuvers with submergence rudder deflections of 4°, 8° and 12° are predicted numerically. The non-uniformity characteristic of the nominal wake field is exacerbated with the increase submergence rudder angle. Then the flow field around the SUBOFF-G submarine fitted with the 4381 propeller is simulated. The axial, transverse and vertical unsteady propeller forces in different maneuvering conditions are compared. In general, as the submarine maneuvers more violently, the harmonic amplitudes of the unsteady force at the 2BPF and 3BPF increased more significantly than that at BPF.

• Journal of the Korean Solar Energy Society
• /
• v.33 no.1
• /
• pp.1-6
• /
• 2013

It is difficult to avoid measurement errors caused by the shading effect of the meteorological tower, which is used for wind resource assessment according to the IEC Standard. This paper presents a validation of the computational flow analysis results by comparing the results with the wind tunnel experiment conducted for Reynolds numbers in the $10^4$ to $10^5$ range, for the preparation of a database for use in an automatic method of correcting met-tower shading errors. A three-dimensional simulation employing the MP (Modified Production) $k-{\varepsilon}$ turbulence model predicted a wind speed deficit in the wake region according to minimum wind speed ratio, within an MAE (Mean Absolute Error) of 2.4%.