Structure of wake of a sharp-edged bluff body in a shallow channel flow

The flow field downstream of a bluff body in a typical open channel flow was explored by two-dimensional particle image velocimetry. Measurements are obtained in horizontal planes at the near-bed, mid-depth and near-surface locations downstream of the body up to a streamwise distance of 10D, where D is the width of the body. The dimensionless streamwise defect velocity profile of the wake flow matches well with the data of a previous investigation and does not reflect any dependency on the distance from the bed. However, the nature of development of the recirculation region is found to be different at the three vertical locations. The time-averaged streamline pattern indicates the existence of a unique nodal pattern close to the bed. The variation of the half-width is also found to be affected by the presence of the bed and the free surface. The bed friction arrests the transverse growth of the shear layer, and the free-surface helps to redistribute the turbulent kinetic energy in the streamwise and transverse directions. Swirling strength analysis is carried out to compare the behavior and statistics of the vortex population in the vertical direction. The prevailing magnitude of the swirling strength is found to be different at the three vertical locations. Bed friction assists to dissipate vorticity rapidly, and therefore reduces the probability of appearance of strong vortices close to the bed.     

Numerical simulation of motion and deformation of ring bubble along body surface

Numerical simulation for fluid flow over an attached rigid body with a deformable ring bubble is analyzed based on the velocity potential theory together with the boundary element method （BEM）. The analysis is focused on the axisymmetric case. The bubble surface is treated as a well defined air-liquid interface and is tracked by a mixed Eulerian-Lagrangian method. The points of intersection between the bubble and body are treated, specially in the numerical procedure. The auxiliary function method is adopted to calculate the pressure on the body surface and in the flow field. The convergence study is undertaken to assess the developed numerical method and the computation code. Some case studies are undertaken in which the interactions between the bubble/body and the incoming flow field are simulated. The effects of various physical parameters on the interactions are investigated.     

3-D scanning-particle-image-velocimetry: Technique and application to a spherical cap wake flow

A 3-D time-resolving and whole-volume Digital-Particle-Image-Velocimetry (DPIV) technique based on the concept of a scanning light-sheet is presented and applied here to the 3-D transient wake phenomena in the spherical cap wake flow. The technique uses a scanning light-sheet for rapid sampling of the flow in depth and a two-camera recording system for stereoscopic 3-D DPIV. Application of a correlation technique in combination with a calibration yields, aside from the correct in-plane displacement, also the out of plane component and thus the total velocity vectors within the planes of the scanning light-sheet. With a high scanning rate in comparison to the characteristic time scales the method provides the 3-D velocity field in space and time. Through the use of conventional video-techniques the temporal evolution of the complete velocity and vorticity field can be obtained quantitatively from experiments. This is demonstrated for the 3-D starting flow around a spherical cap at Re=300. During the starting process, the flow in the wake evolves into a spherical vortex ring where the velocity distribution is very close to the theoretical solution of the Hill-type vortex. Later on, the Hill-type vortex ring deforms and the flow changes from a rotational symmetric stage to a planar symmetric flow with a double-threaded vortical structure which consists of two counter-rotating streamwise vortices similar to the ones observed in sphere wake flow.Presented at the EUROMECH Coll. 335, Image Techniques and Analysis in Fluid Dynamics, 5–7 June 1995, Roma, Italy. A version of this paper has been published in Proc. 7th Int. Symp. Flow Visualization (ed. J.P. Crowder), Begell House Inc., New York, 1995, pp. 715–720.     

The wake flow structure of an open-slit V gutter

This work experimentally investigates the near-wake flow structure behind an open-slit V gutter at airflow speeds between 1 m/s and 20 m/s. With the aid of Schlieren photography and a Dantec three-beam, two-component laser-Doppler anemometry system, the phenomena of vortex shedding and flow recirculation behind the flameholder are well investigated. The results indicate that the interaction between the flow penetrating through the open slit and the shear layer results in an asymmetric wake flow structure. The lower shear layer develops more stably and smoothly than the upper shear layer. Besides, the vortex formation along the lower wing is delayed and at a location farther downstream. The size of the entire recirculation zone is enlarged, and its center shifts toward one of the two wings. Measurements of wake pressure distribution show that the open-slit V gutter generates higher back pressure and thus induces less drag force than the regular V gutter. Moreover, the maximum values of the pressure fluctuation of the wake flow behind the open-slit V gutter reduce 27% and 9%, for the upper side and the lower side, respectively, much lower than those of the regular V gutter. In general, the application of mass bleed from the open-slit V gutter favors both the flame-holding mechanism and the suppression of the flow-induced oscillation.     

无网格MLPG法求解轴对称弹性体扭转问题

应用无网格局部彼得洛夫-伽辽金法(MLPG)研究轴对称弹性体扭转问题,给出了矩阵形式的控制方程,发展了MLPG求解轴对称体弹性扭转问题的数值计算方法。算例分析表明:此方法对求解此类问题具有良好的适应性,数值解能达到理想的计算精度。     

Three-dimensional evolution of vortices and early features of coherent structure in the turbulent wake behind a 2-D circular cylinder

3-D evolution of Kármán vortex filaments and vortex filaments in braid regions in the turbulent wake of a 2-D circular cylinder is investigated numerically based on inviscid vortex dynamics by analyzing the response of the initially 2-D spanwise vortex filaments to periodic spanwise disturbance of varying magnitude, wavelength and initial phase angles. Our results reveal a kind of 3-D vortex system in the wake which consists of large scale horseshoe-shaped vortices and small scale γ-shaped vortex filaments as well as vortex loops. The mechanism and the dynamic process about the generation of streamwise vortical structure and the 3-D coherent structure are reported. currently published in the Chinese Edition of Acta Mechanica Sinica, Vol.25, No.3, 1993 The project supported by National Natural Science Foundation of China and the National Basic Research Project “Nonlinear Science”     

Stability of a supersonic flat-plate wake (Comparison of numerical and experimental results)

The stability of the wake downstream of a plat plate with a blunt trailing edge is numerically investigated at a supersonic freestream velocity. The calculated stability characteristics are compared with the experimental data.     

Behavior of the vorticity vector in supersonic axisymmetric swirl flows behind a detonation wave

The behavior of the vorticity vector on a discontinuity surface arising in a supersonic nonuniform combustible gas flow with the formation of a shock or detonation wave is studied. In the general case, it is a vortex flow with prescribed distributions of parameters. It is demonstrated that the ratio of the tangential component of vorticity to density remains continuous in passing through the discontinuity surface, while the quantities proper become discontinuous. Results calculated for flow vorticity behind a steady-state detonation wave in an axisymmetric supersonic flow of a combustible mixture of gases are presented. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 6, pp. 15–21, November–December, 2007     

Investigations of wake flows of a flat plate in steady,oscillatory and combined flows

Numerical study on near wake flows of a flat plate in three kinds of oncoming flows is made by using the discrete vortex model and improved vorticity creation method. For steady oncoming flow, both gross and detailed features of the wake flow are calculated and discussed. Then, in harmonic oscillatory oncoming flow two different wake flow patterns withK _c=2,4 and 10 are obtained respectively. Our results present a new wake flow pattern for lowKc numbers (Kc<5) describing vortex shedding, pairing and moving in a period of the oscillatory flow starting from rest. The calculated drag and inertia force coefficients are closer to experimental data from the U-tube than the previous results of vortex simulation. For in-line combined oncoming flow the vortex lock-in and dynamic characteristics are simulated. The results are shown to be in good agreement with experiments. The project supported by National Natural Science Fundation of China and LNM of Institute of Mechanics. CAS     

Dynamical mode decomposition of Gurney flap wake flow

The present work uses dynamic mode decomposition (DMD) to analyze wake flow of NACA0015 airfoil with Gurney flap. The physics of DMD is first introduced. Then the PIV-measured wake flow velocity field is decomposed into dynamical modes. The vortex shedding pattern behind the trailing edge and its high-order harmonics have been captured with abundant information such as frequency, wavelength and convection speed. It is observed that high-order dynamic modes convect faster than low-order modes; moreover the wavelength of the dynamic modes scales with the corresponding frequency in power law.     

Wavelet analysis of unsteady flows: Application on the determination of the Strouhal number of the transient wake behind a single cylinder

This paper presents experimental results of the accelerating and decelerating flow in the wake of a cylinder obtained by means of hot wire anemometry measurements in a wind tunnel with high blockage ratio. The analysis was done in Fourier and wavelet spaces. The Strouhal number for Reynolds numbers up to 3 × 10⁴ was studied in a transient flow and compared with the results obtained from steady flows at several velocities uniformly distributed from Re = 1.7 × 10³ to 3 × 10⁴. Results show that the wavelet analysis is a valuable tool to deal with both transient and stationary random phenomena and that is able to capture the characteristics of the transient flow as well as the Fourier analysis can do with the steady state acquisitions.     

Computation of turbulent asymmetric wake

The development of asymmetric wake behind an aerofoil in turbulent incompressible flow has been computed using finite volume scheme for solving two-dimensional Navier-Stokes equations along with the k-ε model of turbulence. The results are compared with available experimental data. It is observed that the computed shift of the point of minimum velocity with distance is sensitive to the prescribed value of the normal component of velocity at the trailing edge of the aerofoil. Making the model constant C_u as a function of streamline curvature and changing the production term in the equation for ε, has only marginal influence on the results.     

Inertial motion of a body in an ideal fluid from the state at rest

The motion of a body in an ideal incompressible fluid flow without vortices in the absence of external forces is considered. It is demonstrated that the body can move inertially from the state at rest if its shape satisfies certain conditions. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 4, pp. 214–219, July–August, 2008.     

Quasi-Euler-Poinsot motion of a rigid body

The phase-plane method of nonlinear oscillation is used to discuss the influence of the small dissipation upon the Euler-Poinsot motion of a rigid body about a fixed point. The equations of phase coordinates are applied instead of Eulerian equations, and the global characteristics of the motion of rigid body are analysed according to the distribution and the type of the singular points. A Chaplygin's sphere on a rough plane, a rigid body in viscous medium and one with a cavity filled with viscous fluid are discussed as examples. It is shown that the motions of rigid bodies dissipated by various physical factors have a common qualitative character. The rigid body tends to make a permanent rotation about the principal axis of the largest moment of inertia. The transitive process can change from oscillatory to aperiodic with the decrease in dissipation.     

A new analytical solution to axisymmetric Blot's consolidation of a finite soil layer

A new analytical method is presented to study the axisymmetric Blot's consolidation of a finite soil layer. Starting from the governing equations of axisymmetric Blot's consolidation, and based on the property of Laplace transform, the relation of basic variables for a point of a finite soil layer is established between the ground surface (z= 0) and the depth z in the Laplace and Hankel transform domains. Combined with the boundary conditions of the finite soil layer, the analytical solution of any point in the transform domain can be obtained. The actual solution in the physical domain can be obtained by inverse Laplace and Hankel transforms. A numerical analysis for the axisymmetric consolidation of a finite soil layer is carried out.     

The oseen flow of finite clusters of spheres and the wake effect on the drag factor

In this paper, the wake effect on drag factor in the axisymmetric Oseen flow of the finite clusters of equally spaced spheres with same size is studied. Putting the Oseen lets on the centres of all the spheres, the series solution of the problem is obtained. By truncating the infinite series and applying the collocation method to solve a set of the linear algebraic equations, the approximate solution of the Oseen flow of finite clusters of spheres and the drag factor for each sphere are presented. The effect of the sphere number and spacing on the drag factor of each sphere under different Reynolds numbers are calculated and the wake effect as well as the shielding effect and the end effect are revealed. The influence of various parameters on the effects is considered and compared with the corresponding results of the Stokes flow. The convergence of the method is also studied numerically in this paper.     

Determining the limiting solutions of nonstationary axisymmetric Hele-Shaw problems

Numerical solution of the Hele-Shaw problem reduces to solution of three boundary-value problems of determining analytic functions of a complex variable in each time step: conformal mapping of the range of the parametric variable to the physical plane, the Dirichlet problems for determining the electric-field strength, and the Riemann-Hilbert problem for calculating partial time derivatives of the coordinates of points of the interelectrode space (the images of the points on the boundary of the parametric plane are fixed). Unlike in the two-dimensional problem, the electric-field strength is determined using integral transformations of an analytic function. Approximation by spline function is performed, and more accurate and steady (than the well-known ones) general solution algorithms for the nonstationary axisymmetric problems are described. Results of a numerical study of the formation of stationary and self-similar configurations are presented. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 87–99, July–August, 2009.