Analysis of super compact finite difference method and application to simulation of vortex-shock interaction

Turbulence and aeroacoustic noise high-order accurate schemes are required, and preferred, for solving complex flow fields with multi-scale structures. In this paper a super compact finite difference method (SCFDM) is presented, the accuracy is analysed and the method is compared with a sixth-order traditional and compact finite difference approximation. The comparison shows that the sixth-order accurate super compact method has higher resolving efficiency. The sixth-order super compact method, with a three-stage Runge-Kutta method for approximation of the compressible Navier-Stokes equations, is used to solve the complex flow structures induced by vortex-shock interactions. The basic nature of the near-field sound generated by interaction is studied.

Flow-Induced Vibrations Of A Cylinder With Two Degrees Of Freedom Near Rigid Plane Boundary

Based on similarity analyses, the flow-induced vibrations of a near-wall cylinder with 2 degrees of freedom are investigated experimentally by employing a hydroelastic apparatus in conjunction with a flume. The cylinder's vibration amplitude, vibration frequency and vortex shedding frequency were measured and analyzed. The effects of gap-to-diameter ratio (e,ID) upon the vibration responses are further investigated. The experimental results indicate that, when the reduced velocity (Vr) is small (e.g. Vr = 1.2 similar to 2.6), only streamwise vibration occurs, and its frequency is quite close to its natural frequency in still water. When increasing Vr (e.g. Vr > 3.4), both streamwise and transverse vibrations of the near-wall cylinder may occur. In the examined range of gap-to-diameter ratio (0.42 < e(0)/D < 2.68), 2 vibration stages (in terms of Vr) of streamwise vibrations usually exist: First Streamwise Vibration (FSV) and Second Streamwise Vibration (SSV). In the SSV stage, the vortex shedding frequency may either undergo a jump to that of the streamwise vibration, or stay consistent with that of the transverse vibration. The amplitudes of transverse vibration are usually much larger than those of streamwise vibration for the same value of e(0)/D. The maximum amplitudes of both streamwise and transverse vibration get larger with the increase of e(0)/D (0.42 < e(0)/D < 2.68).

Control Of Vortex Shedding From A Square Cylinder

Small circular, square, and thin-strip cross-sectional elements are used to suppress vortex shedding from a square cylinder at Reynolds numbers in the range of 1.12 x 10(4)-1.02 x 10(5). The axes of the element and cylinder are parallel. The element's size, position, and angle of attack are varied. Measurements of the fluctuating surface pressures and wake velocities, together with smoke flow visualization, show that vortex shedding from both sides of the cylinder is suppressed and the mean drag and fluctuating lift on the cylinder is reduced if the element is installed in an effective zone downstream of the cylinder. The effective zone of the circular element is shown to be much smaller than those of the other elements. The effects of Reynolds number and blockage ratio are investigated. A phenomenon of monoside vortex shedding is observed. The role of the element's bluffness is investigated and the suppression mechanism is discussed.

A Note On The Numerical Simulations Of Flow Past A Wavy Square-Section Cylinder

The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Karman vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.

VISUALIZATION RESULTS OF AXISYMMETRICAL WAKE FLOW WITH/WITHOUT ACOUSTIC EXCITATION

Visualization results demonstrate the evolution of Kelvin-Helmholtz unstable waves into vortex pairing in a separated shear layer of a blunf circular. The results with acoustic excitation are quite different from that without acoustic excitation, and the phenomenon with excitation in a separated shear layer follows the rule of Devil s staircase, which always occurs in a non-linear dynamical system of two coupling vibrators.

MILU-CG method and the numerical study on the flow around a rotating circular cylinder

A hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient method (MILU-CG), and a high order implicit difference algorithm. The flow around a rotating circular cylinder at Reynolds number R-e = 1000, 200 and the angular to rectilinear speed ratio alpha is an element of (0.5, 3.25) is studied numerically. The long-time full developed features about the variations of the vortex patterns in the wake, and drag, lift forces on the cylinder are given. The calculated streamline contours agreed well with the experimental visualized flow pictures. The existence of critical states and the vortex patterns at the states are given for the first time. The maximum lift to drag force ratio can be obtained nearby the critical states.

Attractor crisis and bursting in a fluid flow with two no-slip directions

Characteristic burtsing behavior is observed in a driven, two-dimensional viscous flow, confined to a square domain and subject to no-slip boundaries. Passing a critical parameter value, an existing chaotic attractor undergoes a crisis, after which the flow initially enters a transient bursting regime. Bursting is caused by ejections from and return to a limited subdomain of the phase space, whereas the precrisis chaotic set forms the asymptotic attractor of the flow. For increasing values of the control parameter the length of the bursting regime increases progressively. Passing another critical parameter value, a second crisis leads to the appearance of a secondary type of bursting, of very large dynamical range. Within the bursting regime the flow then switches in irregular intervals from the primary to the secondary type of bursting. Peak enstrophy levels for both types of bursting are associated to the collapse of a primary vortex into a quadrupolar state.

Periodic vortex breakdown in wide spherical gaps

The flow field with vortex breakdown in wide spherical gaps was studied numerically by a finite difference method under the axisymmetric condition. The result shows that the flow bifurcates to periodic motion as the Reynolds number or the eccentricity of the spheres increases. (C) 1997 American Institute of Physics.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.

Numerical study of bluff body flow structures

Our recent progress in numerical studies of bluff body flow structures and a new method for the numerical analysis of near wake flow field for high Reynolds number flow are introduced. The paper consists of three parts. In part one, the evolution of wake vortex structure and variation of forces on a flat plate in harmonic oscillatory flows and in in-line steady-harmonic combined flows are presented by an improved discrete vortex method, as the Keulegan-Carpenter number (KC) varies from 2 to 40 and ratios of U-m to U-0 are of O(10(-1)), O(10) and O(10), respectively. In part 2, a domain decomposition hybrid method, combining the finite-difference and vortex methods for numerical simulation of unsteady viscous separated flow around a bluff body, is introduced. By the new method, some high resolution numerical visualization on near wake evolution behind a circular cylinder at Re = 10(2), 10(3) and 3 x 10(3) are shown. In part 3, the mechanism and the dynamic process for the three-dimensional evolution of the Karman vortex and vortex filaments in braid regions as well as the early features of turbulent structure in the wake behind a circular cylinder are presented numerically by the vortex dynamics method.

Two-time scales inner solutions and motion of a geostrophic vortex

Based on the authors' previous work, in this paper the systematical analyses on the motion and the inner solutions of a geostrophic vortex have been presented by means of thematched asymptotic expansion method with multiple time scales (S/gh001/2 and α S/gh001/2) and space scales. It has been shown that the leading inner solutions to the core structure in two-time scales analyses are identified with the results in normal one-time scale analyses. The time averages of the first-order solutions on short time variable τ are the same as the first-order solutions obtained in one normal time scale analyses. The geostrophic vortex induces an oscillatory motion in addition to moving with the background flow. The period, amplitude andthe deviation from the mean trajectory depend on the core structure and the initial conditions. The velocity of the motion of vortex center varies periodically and the time average of the velocity on short time variable τ is equal to the value of the local mean velocity.

The motion and the core structure of a geostrophic vortex - one-time scale inner solution and the motion of the vortex

By means of the matched asymptotic expansion method with one-time scale analysis we have shown that the inviscid geostrophic vortex solution represents our leading solution away from the vortex. Near the vortex there is a viscous core structure, with the length scale O(a). In the core the viscous stresses (or turbulent stresses) are important, the variations of the velocity and the equivalent height are finite and dependent of time. It also has been shown that the leading inner solutions of the core structure are the same for two different time scales of S/(ghoo)1/2 and S/a (ghoo)1/2. Within the accuracy of O(a) the velocity of a geostrophic vortex center is equal to the velocity of the local background flow, where the vortex is located, in the absence of the vortex. Some numerical examples demonstrate the contributions of these results.

Characteristics of Wake Flow in Long Bluff Bodies

Vortex shedding is the main characteristics of bluff bodies, which will cause the bluff bodies to vibrate and sometimes result in the structures failure. In this paper the wake flow characteristics of 21 bluff bodies with rectangular, rounded and angular profiles and the length-to-width ratio in the range of 4~12 were deeply studied by Micro ADV. Two parameters, namely the relative intensity of the load due to Karman vortices and the large scale vortex intensity, were introduced to measure the wake flow intensity. Generally, the values of these parameters for different bluff bodies are consistent with each other. The experiment results showed that the key factor affecting the wake flow characteristics is the bluff edge, especially the leading edge geometry. The wake flow in bluff bodies with rounded edge profiles has more regular vortices and becomes more periodic than that in bluff bodies with rectangular ones. A bluff body with angular edged profile was witnessed to have not only small wake loading but small hydraulic resistance also.

Experimental Study on Flow Induced Vibration of a Cylinder with Two Degrees of Freedom Near a Rigid Wall

The flow-induced vibration of a cylinder with two degrees of freedom near a rigid wall under the action of steady flow is investigated experimentally. The vibration amplitude and frequency of the cylinder and the vortex shedding frequency at the wake flow region of the cylinder are measured. The influence of gap-to-diameter ratio upon the amplitude response is analyzed. The experimental results indicate that when the reduced velocity (Vr) is in the range of 1.2 < Vr < 2.6, only streamwise vibration with small amplitude occurs, whose frequency is quite close to its natural frequency in the still water. When the reduced velocity Vr > 3.4, both the streamwise and transverse vibrations of the cylinder occur. In this range, the amplitudes of transverse vibration are much larger than those of streamwise vibrations, and the amplitudes of the streamwise vibration also get larger than those at the range of 1.2 < Vr < 2.6. At the range of Vr > 3.4, the frequency of streamwise vibration undergoes a jump at certain values of Vr, at which the streamwise vibrating frequency is twice as much as the transverse one. However, when the streamwise vibration does not experience a jump, its frequency is the same as that of the transverse vibration. The maximum values of second streamwise and transverse vibration amplitudes increase with increasing gap-to-diameter ratios.

Numerical study on heavy rigid particle motion in a plane wake flow by spectral element method

Experimental particle dispersion patterns in a plane wake flow at a high Reynolds number have been predicted numerically by discrete vortex method (Phys. Fluids A 1992; 4:2244-2251; Int. J. Multiphase Flow 2000; 26:1583-1607). To address the particle motion at a moderate Reynolds number, spectral element method is employed to provide an instantaneous wake flow field for particle dynamics equations, which are solved to make a detail classification of the patterns in relation to the Stokes and Froude numbers. It is found that particle motion features only depend on the Stokes number at a high Froude number and depend on both numbers at a low Froude number. A ratio of the Stokes number to squared Froude number is introduced and threshold values of this parameter are evaluated that delineate the different regions of particle behavior. The parameter describes approximately the gravitational settling velocity divided by the characteristic velocity of wake flow. In order to present effects of particle density but preserve rigid sphere, hollow sphere particle dynamics in the plane wake flow is investigated. The evolution of hollow particle motion patterns for the increase of equivalent particle density corresponds to that of solid particle motion patterns for the decrease of particle size. Although the thresholds change a little, the parameter can still make a good qualitative classification of particle motion patterns as the inner diameter changes.