• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.570-577
• /
• 2001

Three-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for Reynolds number 280 and 300. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme and the pressure and viscosity terms are applied by the 4th order central difference. And in spanwise, Navier-Stokes equation is distributed using of Spectral Method. At Reynolds number 259 the two-dimensional wake becomes linearly unstable to a second branch of modes with wavelength about 1.0 diameters at onset (B-mode). Present results of three-dimensional effects of in wake of a circular cylinder is represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Wind and Structures
• /
• v.17 no.3
• /
• pp.263-274
• /
• 2013

A numerical technique for fluid-structure interaction, which is based on the finite element method (FEM) and computational fluid dynamics (CFD), was developed for application to an industrial chimney equipped with a pendulum tuned mass damper (TMD). In order to solve the structural problem, a one-dimensional beam model (Navier-Bernoulli) was considered and, for the dynamical problem, the standard second-order Newmark method was used. Navier-Stokes equations for incompressible flow are solved in several horizontal planes to determine the pressure in the boundary of the corresponding cross-section of the chimney. Forces per unit length were obtained by integrating the pressure and are introduced in the structure using standard FEM interpolation techniques. For the fluid problem, a fractional step scheme based on a second order pressure splitting has been used. In each fluid plane, the displacements have been taken into account considering an Arbitrary Lagrangian Eulerian approach. The stabilization of convection and diffusion terms is achieved by means of quasi-static orthogonal subscales. For each period of time, the fluid problem was solved and the geometry of the mesh of each fluid plane is updated according to the structure displacements. Using this technique, along-wind and across-wind effects have been properly explained. The method was applied to an industrial chimney in three scenarios (with or without TMD and for different damping values) and for two wind speeds, showing different responses.

• International Journal of Aeronautical and Space Sciences
• /
• v.15 no.4
• /
• pp.356-365
• /
• 2014

This paper is the second in a series and aims to build a high-fidelity mathematical model for a propeller-driven airplane using the propeller's aerodynamics and inertial models, as developed in the first paper. It focuses on aerodynamic models for the fuselage, the main wing, and the stabilizers under the influence of the wake trailed from the propeller. For this, application of the vortex lattice method is proposed to reflect the propeller's wake effect on those aerodynamic surfaces. By considering the maneuvering flight states and the flow field generated by the propeller wake, the induced velocity at any point on the aerodynamic surfaces can be computed for general flight conditions. Thus, strip theory is well suited to predict the distribution of air loads over wing components and the viscous flow effect can be duly considered using the 2D aerodynamic coefficients for the airfoils used in each wing. These approaches are implemented in building a high-fidelity mathematical model for a propeller-driven airplane. Flight dynamic analysis modules for the trim, linearization, and simulation analyses were developed using the proposed techniques. The flight test results for a series of maneuvering flights with a scaled model were used for comparison with those obtained using the flight dynamics analysis modules to validate the usefulness of the present approaches. The resulting good correlations between the two data sets demonstrate that the flight characteristics of the propeller-driven airplane can be analyzed effectively through the integrated framework with the propeller and airframe aerodynamic models proposed in this study.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.3
• /
• pp.215-226
• /
• 2017

This paper presents numerical results of the performance of a marin propeller in cavitating and non-cavitating flow conditions. The geometry and experimental validation data of the propeller are provided in Potsdam Propeller Test Case(PPTC) in the framework of the second International Symposium on Marine Propulsors 2011(SMP'11) workshop. The PPTC includes open water tests, velocity field measurements and cavitation tests. The present numerical analysis was carried out by using the Reynolds averaged Navier-Stokes(RANS) method on a wall-resolved grid ensuring a y+=1, where the SST k-${\omega}$ model was mainly used for turbulence closure. The influence of the turbulence model was investigated in the prediction of the wake field under a non-cavitating flow condition. The propeller tip vortex flows in both cavitating and non-cavitating conditions were captured through adaptation of additional grids. For the cavitation flows at three operation points, Schnerr-Sauer's cavitation model was used with a Volume-Of Fluid(VOF) approach to capture the two-phase flows. The present numerical results for the propeller wake and cavitation predictions including the open water performance showed a qualitatively reasonable agreement with the model test results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.2 s.245
• /
• pp.126-135
• /
• 2006

The pressure drop characteristics in a rotating two-pass duct with rib turbulators are investigated in the present study. The square duct has a hydraulic diameter $(D_h)$ of 26.7 mm, and $1.5mm{\times}1.5mm$ square $90^{\circ}-rib$ turbulators are attached on the leading and trailing walls. The pitch-to-rib height ratio (p/e) is 10. The distance between the tip of the divider and the outer wall of the duct is $1.0D_h$ and the width of divider wall is 6.0mm or $0.225D_h$. The Reynolds number (Re) based on the hydraulic diameter is kept constant at 10,000 to exclude the Reynolds effect, and the rotation number (Ro) is varied from 0.0 to 0.20. The pressure drop distribution, the friction factor and thermal performance are presented for the leading, trailing and the outer surfaces. It is found that the curvature of the $180^{\circ}$-turn produces Dean vortices that cause high pressure drop in the turn. The channel rotation results in pressure drop discrepancy between leading and trailing surfaces so that non-dimensional pressure drops are higher on the trailing surface in the first-pass and on the leading and side surfaces in the second-pass. In the turning region, Dean vortices shown in the stationary case transform into one large asymmetric vortex cell, and subsequent pressure drop characteristics also change. As the rotation number increases, the pressure drop discrepancy enlarges.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.8 s.227
• /
• pp.910-920
• /
• 2004

The heat/mass transfer characteristics in a rotating two-pass duct with and without rib turbulators are investigated in the present study. The square duct has a hydraulic diameter ($D_h$) of 26.7 mm, and $1.5\;mm{\times}1.5\;mm$ square $90^{\circ}$-rib turbulators are attached on the leading and trailing walls. The pitch-to-rib height ratio (p/e) is 10. The Reynolds number based on the hydraulic diameter is kept constant at 10,000 to exclude the Reynolds effect, and the rotation number is varied from 0.0 to 0.20. In the smooth duct, the curvature of the $180^{\circ}$-turn produces Dean vortices that enhance heat/mass transfer in the post-turn region. When rib turbulators are installed, heat/mass transfer is augmented 2.5 times higher than that of the smooth duct since the main flow is turbulated by reattaching and separating in the vicinity of the duct surfaces. The duct rotation results in heat/mass transfer discrepancy so that Sherwood number ratios are higher on the trailing surface in the first-pass and on the leading surface in the second-pass. In the turning region, Dean vortices shown in the stationary case transform into one large asymmetric vortex cell, and subsequent heat/mass transfer characteristics also change. As the rotation number increases, the heat/mass transfer discrepancy enlarges.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.1
• /
• pp.20-29
• /
• 2010

The complicated flow characteristics of upper propeller wake influenced by hull wake are investigated in detail in the present study. A two-frame PIV (particle image velocimetry) technique was employed to visualize the upper propeller wake region. As the upper hull wake affects strongly propeller inflow, upper propeller wake shows much unstable vortical behavior, especially in the tip vortices. Velocity field measurements were conducted in a cavitation tunnel with a simulated hull wake. Generally, the hull wake generated by the hull of a marine ship may cause different loading distributions on the propeller blade in both upper and lower propeller planes. The unstable upper propeller wake caused by the ship's hull is expressed in terms of turbulent kinetic energy (TKE) and is identified by using the proper orthogonal decomposition (POD) method to characterize the coherent flow structure in it. Instabilities appeared in the eigen functions higher than the second one, giving unsteadiness to the downstream flow characteristics. The first eigen mode would be useful to find out the tip vortex positions immersed in the unstable downstream region.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.8
• /
• pp.829-837
• /
• 2012

Direct numerical simulations were carried out for turbulent channel flows with $Re_{\tau}$ = 180, 395 and 590 to investigate the turbulent flow structure related to the Reynolds shear stress. By examining the probability density function, the second quadrant (Q2) events with the largest contribution to the mean Reynolds shear stress were identified. The change in the inclination angle of Q2 events varies with wall units in $y^+<50$ and with the channel half height in y/h > 0.5. Conditionally averaged flow fields for the Q2 event show that the flow structures associated with Reynolds shear stress are a quasi-streamwise vortex in the buffer layer and a hairpin-shaped vortex in the outer layer. Three-dimensional visualization of the distribution of high Reynolds shear stress reveals that the organization of hairpin vortices in the outer layer having a size of 1.5~3 h is associated with large-scale motions with high Reynolds shear stress in the outer layer.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.4
• /
• pp.439-448
• /
• 2009

Three inducers were designed to avoid cavitation instabilities. This was accomplished by avoiding the interaction of tip cavity with the leading edge of the next blade. The first one was designed with extremely larger leading edge sweep, the second and third ones were designed with smaller incidence angle by reducing the inlet blade angle or increasing the design flow rate, respectively. The inducer with larger design flow rate has larger outlet blade angle to obtain sufficient pressure rise. The inducer with larger sweep could suppress the cavitation instabilities in higher flow rates more than 95% of design flow coefficient, owing to weaker tip leakage vortex cavity with stronger disturbance by backflow vortices. The inducer with larger outlet blade angle could avoid the cavitation instabilities at higher flow rates, owing to the extension of the tip cavity along the suction surface of the blade. The inducer with smaller inlet blade angle could avoid the cavitation instabilities at higher flow rates, owing to the occurrence of the cavity first in the blade passage and its extension upstream. The cavity shape and suction performance were reasonably simulated by three dimensional CFD computations under the steady cavitating condition, except for the backflow vortex cavity. The difference in the growth of cavity for each inducer is explained from the difference of the pressure distribution on the suction side of the blades.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.4
• /
• pp.363-374
• /
• 2009

The 2D flow around 13 similar stay-vane profiles with different trailing edge geometries is investigated to determinate the main characteristics of the excitation forces for each one of them and their respective dynamic behaviors when modeled as a free-oscillating system. The main goal is avoid problems with cracks of hydraulic turbines components. A stay vane profile with a history of cracks was selected as the basis for this work. The commercial finite-volume code $FLUENT^{(R)}$ was employed in the simulations of the stationary profiles and, then, modified to take into account the transversal motion of elastically mounted profiles with equivalent structural stiffness and damping. The k-$\omega$ SST turbulence model is employed in all simulations and a deforming mesh technique used for models with profile motion. The static-model simulations were carried out for each one of the 13 geometries using a constant far field flow velocity value in order to determine the lift force oscillating frequency and amplitude as a function of the geometry. The free-oscillating stay-vane simulations were run with a low mass-damping parameter ($m^*{\xi}=0.0072$) and a single mean flow velocity value (5m/s). The structural bending stiffness of the stay-vane is defined by the Reduced Velocity parameter (Vr). The dynamic analyses were divided into two sets. The first set of simulations was carried out only for one profile with $2{\leq}Vr{\leq}12$. The second set of simulations focused on determining the behavior of each one of the 13 profiles in resonance.