• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.8-15
• /
• 2005

Aerodynamic sounds generated by a uniform flow around a two-dimensional circular cylinder at Re=150 are simulated by applying the finite difference lattice Boltzmann method. Thethird-order-accurate up-wind scheme (UTOPIA) is used for the spatial derivatives, and the second-order-accurate Runge-Kutta scheme is applied for the time marching. We have succeed in capturing very small pressure fluctuations with the same frequency of the Karman vortex street compared with the pressure fluctuation around a circular cylinder. The propagation velocity of the acoustic waves shows that the points of peak pressure are biased upstream due to the Doppler effect in the uniform flow. For the downstream, on the other hand, it is faster. It is also apparent that the amplitude of sound pressure is proportional to r /sup -1/2/,r being the distance from the center of the circular cylinder. To investigate the effect of the lattice dependence, furthermore, 2D computations of the tone noises radiated by a square cylinder and NACA0012 with a blunt trailing edge at high incidence and low Reynolds number are also investigate.

• Structural Engineering and Mechanics
• /
• v.25 no.2
• /
• pp.161-180
• /
• 2007

A first endeavor is made to exploit the differential quadrature method (DQM) as a simple, accurate, and computationally efficient numerical tool for the large deformation analysis of thin laminated composite skew plates, which has very strong singularity at the obtuse vertex. The geometrical nonlinearity is modeled by using Green's strain and von Karman assumption. A recently developed DQ methodology is used to exactly implement the multiple boundary conditions at the edges of skew plates, which is a major draw back of conventional DQM. Using oblique coordinate system and the DQ methodology, a mapping-DQ discretization rule is developed to simultaneously transform and discretize the equilibrium equations and the related boundary conditions. The effects of skew angle, aspect ratio and different types of boundary conditions on the convergence and accuracy of the presented method are studied. Comparing the results with the available results from other numerical or analytical methods, it is shown that accurate results are obtained even when using only small number of grid points. Finally, numerical results for large deflection behavior of antisymmetric cross ply skew plates with different geometrical parameters and boundary conditions are presented.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.5
• /
• pp.1310-1321
• /
• 1994

Turbulent coherent structures in the intermediate wake of a stationary and rotating cylinder, spin rate S=0.7, situated in a uniform were experimentally investigated using a conditionalphase average technique. Measurements were carried out at a section of 8.5 diameters downstream form the center of cylinder and a Reynolds number of $Re=6.5{\times}10^{3}.$/TEX> The phase averaged velocity and velocity vector fields, contours of vorticity, turbulent intermittency function and velocity fluctuation energy are presented and discussed in relation to the large scale coherent structures by Karman vortices that shed periodically from the cylinder. Coherent wake structures of the rotating cylinder is almost identical with stationary cylinder, but the lateral displacement and shrinkage of turbulent wake region is occured by rotation. Rotation of the cylinder result in that the deflection of wake center to deceleration region(Y/D${\simeq}-0.3)$ and the decrease of mean velocity defect(10%), vorticity strength of large scale structures(19%), total velocity fluctuation energy(12%).

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1801-1806
• /
• 2004

Stochastic nature of subgrid-scale stress causes the predictability problem in large eddy simulation (LES) by which the LES solution field decorrelates with field from filtered directnumerical simulation (DNS). In order to evaluate the predictability limit in a priori sense, the information on the interplay between resolved scale and subgrid-scale (SGS) is required. In this study, the analysis on the inter-scale interaction is performed by applying tophat and cutoff filters to DNS database of flow over a circular cylinder at Reynolds number of 3900. The effect of filter shape is investigated on the interpretation of correlation between scales. A critique is given on the use of tophat filter for SGS analysis using DNS database. It is shown that correlations between Karman vortex and SGS kinetic energy drastically decrease when the cutoff filter is used, which implies that the small scale universality holds even in the presence of the large scale coherent structure.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.21 no.2
• /
• pp.9-17
• /
• 1984

Laminar and turbulent boundary layers on a rotating sector and a helical blade are calculated by differential method. The estimation of three dimensional viscous flows provide quite useful informations for the design of propellers and turbo-machinery. A general method of calculation is presented in this paper. Calculated laminar boundary layer on a sector shows smooth development of flows from Blasius' solution at the leading edge to von Karman's solution of a rotating disk at the down-stream. Eddy viscosity model is adopted for the calculation of turbulent flows. Turbulent flows on a rotating blade show similar characters as laminar flows. But cross-flow angle of turbulent flows are reduced in comparison with laminar boundary layers. Effects of rotation make flow structures significantly different from two-dimensional flows. In the range of Reynolds number of model scale propellers, large portion of the blade are still in the transition region from laminar to turbulent flows. Therefore viscous flow pattern might be quite different on the blade of model propeller. The present method of calculation is to be useful for the research of scale effects, cavitation, and roughness effects of propeller blades.

• Wind and Structures
• /
• v.19 no.4
• /
• pp.389-404
• /
• 2014

The flow characteristics of a circular cylinder surrounded by an outer permeable cylinder were experimentally investigated using Particle Image Velocimetry Technique in deep water flow. In order to consider the effects of diameter and porosity of the outer cylinder on flow structures of the inner cylinder, five different outer cylinder diameters (D=37.5, 52.5, 60, 75 and 90 mm) and eight different porosities (${\beta}$=0.4, 0.5, 0.6, 0.65, 0.7, 0.75, 0.8 and 0.85) were selected. During the experiments, the diameter of inner cylinder was kept constant as d=30 mm. The depth-averaged free-stream velocity was adjusted as U=0.156 m/s, which corresponds to the Reynolds number of Re=5000 based on the inner cylinder diameter. It has been concluded that both the outer permeable cylinder diameter and the porosity have important influences on the attenuation of vortex shedding in the wake region. The presence of outer permeable cylinder decreases the magnitude of Reynolds shear stress and turbulent kinetic energy compared to the bare cylinder case. Moreover, the spectral analysis of vortex shedding frequency has revealed that the dominant frequency of vortex shedding downstream of the cylinder arrangement also reduces substantially due to the weakened Karman shear layer instability.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.1
• /
• pp.26-32
• /
• 2019

The vortex flow behind a bluff body has been a subject of interest for a very long time because of its engineering applicability such as to vortex induced vibration. In the near wake of a bluff body, vortices are periodically shed in two shear layers, which originate in the trailing edges. The far wake is made up of the classical Karman vortices, which are connected together by streamwise and spanwise vortices. These vortex formations have been studied in many experimental and numerical ways. However, most of the studies considered non-cavitating flow. In this study, we investigated cavitating flow in the wake of a two-dimensional wedge. Experiments were conducted in a cavitation tunnel of Chungnam National University. Using a particle image velocimetry (PIV), we measured the velocity fields under two different flow conditions: non-cavitating and cavitating regimes. We also investigated the vortex shedding frequencies using an absolute pressure transducer mounted on the top of the test window. Throughout the experiments, it was found that the shedding frequency of the vortex was strongly affected by cavitation, and the Strouhal number could exceed its value in the non-cavitating regime.

• Journal of Mechanical Science and Technology
• /
• v.17 no.9
• /
• pp.1388-1396
• /
• 2003

The hydrodynamic instability of the three-dimensional boundary layer on a rotating disk introduces a periodic modulation of the mean flow in the form of stationary cross flow vortices. The instability labeled Type II by Faller occurs first at lower Reynolds number than that of well known Type I instability. Detailed numerical values of the amplification rates, neutral curves and other characteristics of the two instabilities have been calculated over a wide range of parameters. Presented are the neutral stability results concerning the two instability modes by solving the appropriate linear stability equations reformulated not only by considering whole convective terms but also by correcting some errors in the previous stability equations. The present stability results agree with the previously known ones within reasonable limit. Consequently, the flow is found to be always stable for a disturbance whose dimensionless wave number is greater than 0.75. Some spatial amplification contours have been computed for the stationary disturbance wave, whose azimuth angle $\varepsilon$= 11.29$^{\circ}$ to 15$^{\circ}$ and for the moving disturbance wave, whose azimuth angle $\varepsilon$ = 12.5$^{\circ}$ to 15$^{\circ}$. Also, some temporal amplification contours have been computed for the stationary disturbance wave, whose azimuth angle $\varepsilon$= 11.29$^{\circ}$ to 15$^{\circ}$ and for the moving disturbance wave, whose azimuth angle $\varepsilon$= 12$^{\circ}$ to 15$^{\circ}$. The flow instability was observed by using a white titanium tetrachloride gas over rotating disk system. When the numerical results are compared to the present experimental data, the numerical results agree quantitatively, indicating the existence of the selective frequency mechanism.

• Structural Engineering and Mechanics
• /
• v.45 no.5
• /
• pp.677-691
• /
• 2013

A semi-analytical finite strip method is developed for analyzing the post-buckling behavior of rectangular composite laminated plates of arbitrary lay-up subjected to progressive end-shortening in their plane and to normal pressure loading. In this method, all the displacements are postulated by the appropriate harmonic shape functions in the longitudinal direction and polynomial interpolation functions in the transverse direction. Thin or thick plates are assumed and correspondingly the Classical Plate Theory (CPT) or Higher Order Plate Theory (HOPT) is applied. The in-plane transverse deflection is allowed at the loaded ends of the plate, whilst the same deflection at the unloaded edges is either allowed to occur or completely restrained. Geometric non-linearity is introduced in the strain-displacement equations in the manner of the von-Karman assumptions. The formulations of the finite strip methods are based on the concept of the principle of the minimum potential energy. The Newton-Raphson method is used to solve the non-linear equilibrium equations. A number of applications involving isotropic plates, symmetric and unsymmetric cross-ply laminates are described to investigate the through-thickness shearing effects as well as the effect of pressure loading, end-shortening and boundary conditions. The study of the results has revealed that the response of the composite laminated plates is particularly influenced by the application of the Higher Order Plate Theory (HOPT) and normal pressure loading. In the relatively thick plates, the HOPT results have more accuracy than CPT.

• Steel and Composite Structures
• /
• v.25 no.1
• /
• pp.105-116
• /
• 2017

This paper presents the dynamic stability study of laminated composite plates with different force combinations and aspect ratios. Optimum non-diverging stacking is obtained for certain loading combination and aspect ratio. In addition, the stability force is maximized for a definite operating frequency. A dynamic version of the principle of virtual work for laminated composites is used to obtain force-frequency relation. Since dynamic stiffness governs the divergence or flutter, an efficient optimization method is necessary for the response functional and the relevant constraints. In this way, a model based on the ant colony optimization (ACO) algorithm is proposed to search for the proper stacking. The ACO algorithm is used since it treats with large number of dynamic stability parameters. Governing equations are formulated using classic laminate theory (CLT) and von-Karman plate technique. Load-frequency relations are explicitly obtained for fundamental and secondary flutter modes of simply supported composite plate with arbitrary aspect ratio, stacking and boundary load, which are used in optimization process. Obtained results are compared with the finite element method results for validity and accuracy convince. Results revealed that the optimum stacking with stable dynamic response and maximum critical load is in angle-ply mode with almost near-unidirectional fiber orientations for fundamental flutter mode. In addition, short plates behave better than long plates in combined axial-shear load case regarding stable oscillation. The interaction of uniaxial and shear forces intensifies the instability in long plates than short ones which needs low-angle layup orientations to provide required dynamic stiffness. However, a combination of angle-ply and cross-ply stacking with a near-square aspect ratio is appropriate for the composite plate regarding secondary flutter mode.