Effects of Reynolds number on physiological‐type pulsatile flows in a pipe with ring‐type constrictions

The effects of Reynolds number on the physiological‐type of laminar pulsatile flow fields within the vicinity of mechanical ring‐type constriction in small pipes were studied numerically. The parameters considered are: the Reynolds number (Re) in the range of 50–1500; Strouhal number (St) in the range of 0.00156–3.98; Womersley number (Nw) from 0.0 to 50.0. The pulsatile flows considered were physiological‐type of simulated flows. Within a pulsating cycle, detailed flow characteristics were studied through the pulsating contours of streamline (ψ), vorticity (Ω), shear stress (τ) and isobar. The relations between the instantaneous flow rate (Q) and instantaneous pressure gradients (dp/dz) are observed to be elliptic. The relations between the instantaneous flow rate (Q) and pressure loss (P_loss) are quadratic. Linear relations were observed between the instantaneous flow rate (Q) and the maximum velocity, maximum vorticity and maximum shear stress. The Reynolds number of the flow in a pulsating cycle was found to have significant effects on the recirculation length and the pressure gradient within the pulsatile flow regime. Copyright © 1999 John Wiley & Sons, Ltd.     

NUMERICAL STUDIES OF TRANSITIONAL TURBULENT PULSATILE FLOW IN PIPES WITH RING–TYPE CONSTRICTIONS

Pulsatile flows in the vicinity of mechanical ring-type constrictions in pipes were studied for transitional turbulent flow with a Reynolds number (Re) of the order of 10⁴. The Womersley number (Nw) is in the range 30–50, with a corresponding Strouhal number (St) range of 0·0143–0·0398. The pulsatile flows considered are a pure sinusoidal flow, a physiological flow and an experimental pulsatile flow profile for mechanical aortic valve flow simulations. Transitional laminar and turbulent flow characteristics in an alternating manner within the pulsatile flow fields were studied numerically. It was observed that fluid accelerations tend to suppress the development of flow disturbances. All the instantaneous maximum values of turbulent kinetic energy, turbulent viscosity and turbulent shear stress are smaller during the acceleration phase than during the deceleration period. Various parametric equations have been formulated through numerical experimentation to better describe the relationships between the instantaneous flow rate (Q), the pressure loss (ΔP), the maximum velocity (V_max), the maximum vorticity (ζ_max), the maximum wall vorticity (ζ_w,max), the maximum shear stress (τ_max) and the maximum wall shear stress (τ_w,max) for turbulent pulsatile flow in the vicinity of constrictions in the vascular tube. An elliptic relationship has been found to exist between the instantaneous flow rate and the instantaneous pressure gradient. Other linear and quadratic relations between various flow parameters were also obtained.     

Numerical study of effects of pulsatile amplitude for transitional turbulent pulsatile flow in pipes with ring‐type constrictions

The effects of pulsatile amplitude on sinusoidal transitional turbulent flows through a rigid pipe in the vicinity of a sharp‐edged mechanical ring‐type constriction have been studied numerically. Pulsatile flows were studied for transitional turbulent flow with Reynolds number (Re) of the order of 10⁴, Womersley number (Nw) of the order of 50 with a corresponding Strouhal number (St) of the order of 0.04. The pulsatile flow considered is a sinusoidal flow with dimensionless amplitudes varying from 0.0 to 1.0. Transitional laminar and turbulent flow characteristics in an alternative manner within the pulsatile flow fields were observed and studied numerically. The flow characteristics were studied through the pulsatile contours of streamlines, vorticity, shear stress and isobars. It was observed that fluid accelerations tend to suppress the development of flow disturbances. All the instantaneous maximum values of turbulent kinetic energy, turbulent viscosity, turbulent shear stress are smaller during the acceleration phase when compared with those during deceleration period. Various parametric equations within a pulsatile cycle have also been formulated through numerical experimentations with different pulsatile amplitudes. In the vicinity of constrictions, the empirical relationships were obtained for the instantaneous flow rate (Q), the pressure gradient (dp/dz), the pressure loss (P_loss), the maximum velocity (V_max), the maximum vorticity (ζ_max), the maximum wall vorticity (ζ_w,max), the maximum shear stress (τ_max) and the maximum wall shear stress (τ_w,max). Elliptic relation was observed between flow rate and pressure gradient. Quadratic relations were observed between flow rate and the pressure loss, the maximum values of shear stress, wall shear stress, turbulent kinematic energy and the turbulent viscosity. Linear relationships exist between the instantaneous flow rate and the maximum values of vorticity, wall vorticity and velocity. The time‐average axial pressure gradient and the time average pressure loss across the constriction were observed to increase linearly with the pulsatile amplitude. Copyright © 1999 John Wiley & Sons, Ltd.     

A numerical study of pulsatile laminar flows in a pipe with a ring-type constriction

Numerical simulations have been carried out to study pulsatile laminar flows in a pipe with an axisymmetric ringtype constriction. Three types of pulsatile flows were investigated, namely a physiological flow, a pure sinusoidal flow and a non-zero mean velocity sinusoidal flow. The laminar flow governing equations were solved by the SIMPLE algorithm on a non-staggered grid and a modified Crank-Nicolson approximation was used to discretrize the momentum equations with respect to time. The maximum flow Reynolds numer (Re) is 100. The Womersley number (N_w) ranges from 0 to 50, with the corresponding Strouhal number (St) ranging from 0 to 3·98. The constriction opening ratio (d/D) and thickness ratio (h/D) are fixed at 0·5 and 0·1 respectively. Within the time period investigated, all these pulsatile flows include both forward and backward flows. The unsteady recirculation region and the recirculation points change in size and location with time. For N_w ≤ 1 and St≤ 1·56 x 10^?3 the three pulsatile flows have the same simple relation between the instantaneous flow rate and pressure loss (Δp) across the constriction and the pressure gradient in the axial direction (dp/dz) in the fully developed flow region. The phase angles between the flow rate and pressure loss and the pressure gradient are equal to zero. With increasing N_w and St, the phase angle between the flow rate and the dp/dz becomes larger and has its maximum value of 90° at N_w = 50 and St = 3·98. The three pulsatile flows also show different relations between the flow rate and the pressure gradient. The pure sinusoidal flow has the largest maximum pressure gradient and the non-zero mean velocity sinusoidal flow has the smallest. For larger N_w and St the fully developed velocity profiles in the fully developed flow region have a smaller velocity gradient along the radial direction in the central region. The maximum recirculation length increases for N_w ranging from 0 to 4·2, while this length becomes very small at N_w = 50 and St = 3·98. The deceleration tends to enlarge the recirculation region and this effect appears for N_w ≥ 3 and St ≥ 1·43×10^?2. Linear relations exist between the flow rate and the instantaneous maximum values of velocity, vorticity and shear stress.     

Numerical simulation of shear effect on vortex shedding behind a square cylinder

This paper is concerned with the numerical simulation of the flow structure around a square cylinder in a uniform shear flow. The calculations were conducted by solving the unsteady 2D Navier–Stokes equations with a finite difference method. The effect of the shear parameter K of the approaching flow on the vortex-shedding Strouhal number and the force coefficients acting on the square cylinder is investigated in the range K=0·0–0·25 at various Reynolds numbers from 500 to 1500. The computational results are compared with some existing experimental data and previous studies. The effect of shear rate on the Strouhal number and the force acting on the cylinder has a tendency to reduce the oscillation. The Strouhal number, mean drag and amplitude of the fluctuating force tend to decrease as the shear rate increases, but show no significant change at low shear rate. Increasing the Reynolds number decreases the Strouhal number and increases the force acting on the cylinder. At high shear rate the shedding frequencies of the fluctuating drag and lift coefficients are identical. © 1997 John Wiley & Sons, Ltd.     

Turbulence in a skewed three‐dimensional wall‐bounded shear flow: effect of mean vorticity on structure modification

Mean‐flow three‐dimensionalities affect both the turbulence level and the coherent flow structures in wall‐bounded shear flows. A tailor‐made flow configuration was designed to enable a thorough investigation of moderately and severely skewed channel flows. A unidirectional shear‐driven plane Couette flow was skewed by means of an imposed spanwise pressure gradient. Three different cases with 8°, 34°and 52°skewing were simulated numerically and the results compared with data from a purely two‐dimensional plane Couette flow. The resulting three‐dimensional flow field became statistically stationary and homogeneous in the streamwise and spanwise directions while the mean velocity vector V and the mean vorticity vector Ω remained parallel with the walls. Mean flow profiles were presented together with all components of the Reynolds stress tensor. The mean shear rate in the core region gradually increased with increasing skewing whereas the velocity fluctuations were enhanced in the spanwise direction and reduced in the streamwise direction. The Reynolds shear stress is known to be closely related to the coherent flow structures in the near‐wall region. The instantaneous and ensemble‐averaged flow structures were turned by the skewed mean flow. We demonstrated for the medium‐skewed case that the coherent structures should be examined in a coordinate system aligned with V to enable a sound interpretation of 3D effects. The conventional symmetry between Case 1 and Case 2 vortices was broken and Case 1 vortices turned out to be stronger than Case 2. This observation is in conflict with the common understanding on the basis of the spanwise (secondary) mean shear rate. A refined model was proposed to interpret the structure modifications in three‐dimensional wall‐flows. What matters is the orientation of the mean vorticity vector Ω relative to the vortex vorticity vector ω _v, that is, the sign of Ω · ω _v. In the present situation, Ω · ω _v > 0 for the Case 1 vortices causing a strengthening relative to the Case 2 vortices. Copyright © 2011 John Wiley & Sons, Ltd.     

THE STRUCTURE OF TURBULENT SHEAR FLOW AROUND A TWO-DIMENSIONAL POROUS FENCE HAVING A BOTTOM GAP

The flow field behind porous fences of geometric porosity ε=38·5% with various bottom gaps (G) has been investigated using a hybrid PTV velocity field measurement technique. Four gap ratiosG /H=0·0, 0·1, 0·2 and 0·3 were tested in this study. One thousand instantaneous velocity vector fields in the x–y plane were consecutively measured for each gap ratio. The free-stream velocity was fixed at 10cm/s and the corresponding Reynolds number based on the fence height (H) was Re=2985. The results show that the gap ratio G/H=0·1 gives the best shelter effect among the four gap ratios tested in this study, having a small shelter parameter ψ in a large area behind the fence. As the gap ratio increases, the region of mean velocity reduction decreases and the lower shear layer developed from the bottom gap expands upward. From the spatial distributions of turbulence statistics including turbulence intensities, Reynolds shear stress and turbulent kinetic energy, the wake characteristics can be divided into two categories depending on the gap ratio. When the gap ratio is aboveG /H=0·2, the turbulence statistics have large values in the lower shear layer. For the gap ratio G/H≤0·1, however, the lower shear layer displays small turbulence-statistics values and approach those of the no-gap case (G/H=0) with increasing distance downstream. In the upper shear layer separated from the fence top, the turbulence statistics are nearly independent of the gap ratio.     

Numerical modelling of confined flow past a cylinder of square cross-section at various orientations

Results are presented for the unsteady, two-dimensional flow and heat transfer due to a square obstruction of diameter d located asymmetrically between the parallel sliding walls of a channel with length-to-height ratio W/H = 6·44. Analysis is based on the numerical solution of spatially and temporally second-order accurate finite difference approximations of the transport equations expressed in curvilinear co-ordinates. Laminar, constant property flow is assumed for obstruction configurations in which the blockage ratio is d/H = 0·192, the nearest-wall distances are g/d = 0·2, 0·5 and 1, the orientation angles are α=0°, 10° and 20° and the Reynolds numbers are Re=100, 500, and 1000. Preparatory testing of the numerical procedure was performed for a variety of documented flows to verify its physiconumerical accuracy and obtain estimates of the residual grid-dependent uncertainties in the variables calculated. Heat transfer, drag and lift coefficients and Strouhal numbers for the present flow were finally calculated to within 4%–7% of their grid-dependent values using non-uniformly spaced grids consisting of (x=99, y=55) nodes. Above a critical value of the Reynolds number, which depends on the geometrical parameters, the flow is characterized by alternate vortex shedding from the obstruction top and bottom surfaces. Streamline, vorticity and particle streakline plots provide qualitative impressions of the unsteady vortical flow. Especially noteworthy are the extremes in the lift coefficient which ranges from large positive values for an obstruction with g/d=0·2 and α=10° to negative values for one with g/d=0·5 and α=0°. Both the drag and lift coefficients as well as the Strouhal number exhibit non-monotonic variations with respect to the parameters explored. Asymmetries in the obstruction location and orientation account for relatively large vortex-induced periodic variations in heat transfer, especially along the wall nearest the obstruction. Notable differences are also predicted for the heat transfer coefficients of the individual obstruction surfaces as a function of the orientation angle.     

Numerical studies of slow viscous rotating flow past a sphere. III

The flow of steady incompressible viscous fluid rotating about the z-axis with angular velocity ω and moving with velocity u past a sphere of radius a which is kept fixed at the origin is investigated by means of a numerical method for small values of the Reynolds number Re_ω. The Navier–Stokes equations governing the axisymmetric flow can be written as three coupled non-linear partial differential equations for the streamfunction, vorticity and rotational velocity component. Central differences are applied to the partial differential equations for solution by the Peaceman–Rachford ADI method, and the resulting algebraic equations are solved by the ‘method of sweeps’. The results obtained by solving the non-linear partial differential equations are compared with the results obtained by linearizing the equations for very small values of Re_ω. Streamlines are plotted for Ψ = 0·05, 0·2, 0·5 for both linear and non-linear cases. The magnitude of the vorticity vector near the body, i.e. at z = 0·2, is plotted for Re_ω = 0·05, 0·24, 0·5. The correction to the Stokes drag as a result of rotation of the fluid is calculated.     

Grid‐free surface vorticity method applied to flow induced vibration of flexible cylinders

In order to study cross flow induced vibration of heat exchanger tube bundles, a new fluid–structure interaction model based on surface vorticity method is proposed. With this model, the vibration of a flexible cylinder is simulated at Re=2.67 × 10⁴, the computational results of the cylinder response, the fluid force, the vibration frequency, and the vorticity map are presented. The numerical results reproduce the amplitude‐limiting and non‐linear (lock‐in) characteristics of flow‐induced vibration. The maximum vibration amplitude as well as its corresponding lock‐in frequency is in good agreement with experimental results. The amplitude of vibration can be as high as 0.88D for the case investigated. As vibration amplitude increases, the amplitude of the lift force also increases. With enhancement of vibration amplitude, the vortex pattern in the near wake changes significantly. This fluid–structure interaction model is further applied to simulate flow‐induced vibration of two tandem cylinders and two side‐by‐side cylinders at similar Reynolds number. Promising and reasonable results and predictions are obtained. It is hopeful that with this relatively simple and computer time saving method, flow induced vibration of a large number of flexible tube bundles can be successfully simulated. Copyright © 2004 John Wiley & Sons, Ltd.     

Interference effect of two equal‐sized square cylinders in tandem arrangement: With planar shear flow

Numerical simulations have been performed for flow past two equal‐sized square cylinders in tandem arrangement subjected to incoming planar shear flow. Effect of L/d ratio and the shear parameter has been studied. The range of L/d ratio (ratio of center‐to‐center distance (L) to cylinder width (d)) is varied from 2 to 7 and the non‐dimensional shear parameter (K) is varied from 0.0 to 0.4 in steps of 0.1. For all the cases the Reynolds number (Re) based on centerline velocity and cylinder width is fixed at 100. The results are compared with that of isolated square cylinder with uniform flow. Strouhal number decreases with increasing shear parameter. There are more than one shedding frequency at high shear parameters and L/d ratios. The mean drag coefficient is decreased with shear parameter and lesser than that of the single cylinder. The root mean square (RMS) value of both lift and drag coefficients is higher for the downstream cylinder for all values of shear parameter. With increasing L/d ratio, for both lift and drag, the RMS value increases and then decreases for upstream cylinder, whereas it continuously increases for the downstream cylinder. The stagnation point is moved towards the top leading edge with increasing shear. The critical L/d ratio, which is defined as the distance between two cylinders, beyond which the vortex shedding from the upstream cylinder occurs, decreases with increasing shear parameter. Copyright © 2007 John Wiley & Sons, Ltd.     

Streamwise evolution of an inclined cylinder wake

The streamwise evolution of an inclined circular cylinder wake was investigated by measuring all three velocity and vorticity components using an eight-hotwire vorticity probe in a wind tunnel at a Reynolds number Re_d of 7,200 based on free stream velocity (U _∞) and cylinder diameter (d). The measurements were conducted at four different inclination angles (α), namely 0°, 15°, 30°, and 45° and at three downstream locations, i.e., x/d = 10, 20, and 40 from the cylinder. At x/d = 10, the effects of α on the three coherent vorticity components are negligibly small for α ≤ 15°. When α increases further to 45°, the maximum of coherent spanwise vorticity reduces by about 50%, while that of the streamwise vorticity increases by about 70%. Similar results are found at x/d = 20, indicating the impaired spanwise vortices and the enhancement of the three-dimensionality of the wake with increasing α. The streamwise decay rate of the coherent spanwise vorticity is smaller for a larger α. This is because the streamwise spacing between the spanwise vortices is bigger for a larger α, resulting in a weak interaction between the vortices and hence slower decaying rate in the streamwise direction. For all tested α, the coherent contribution to [`(v<sup>2</sup>)] \overline{{v^{2}}} is remarkable at x/d = 10 and 20 and significantly larger than that to [`(u²)] \overline{{u^{2}}} and [`(w<sup>2</sup>)]. \overline{{w^{2}}}. This contribution to all three Reynolds normal stresses becomes negligibly small at x/d = 40. The coherent contribution to [`(u²)] \overline{{u^{2}}} and [`(v²)] \overline{{v^{2}}} decays slower as moving downstream for a larger α, consistent with the slow decay of the coherent spanwise vorticity for a larger α.     

Laminar fluid convection between concentric and eccentric heated horizontal rotating cylinders for low-Prandtl-number fluids

Numerical experiments are performed to study rotational effects on the mixed convection of low-Prandtlnumber fluids enclosed between the annuli of concentric and eccentric horizontal cylinders. The inner cylinder is assumed to be heated and rotating. The rotational Reynolds number considered is in the range where the effect of Taylor vortices is negligible. The Prandtl number of the fluid considered is in the range 0·01–1·0. The Rayleigh number considered is up to 10⁶. A non-uniform mesh transformation technique coupled with the introduction of ‘false transient’ parameters to the vorticity and streamfunction-vorticity expressions was used to solve the governing set of equations. Results show that when the inner cylinder is made to rotate, the multicellular flow patterns observed in stationary cylindrical annuli subside in a manner depending on the Prandtl number of the fluids. Eventually the flow tends toward a uniform flow similar to that of a solid body rotation. For a fixed Rayleigh number and with a Prandtl number of the order of 1·0, when the inner cylinder is made to rotate, the mean Nusselt number is observed to decrease throughout the flow. For lower Prandtl number of the order 0·1–0·01 the mean Nusselt number remained fairly constant when the inner cylinder was made to rotate. The mean Nusselt numbers obtained were also compared with available data from other investigators.     

Study on the effect of steady,simple pulsatile and physiological pulsatile flows through a stenosed artery

In the present paper, the comparison of steady, simple pulsatile flow and physiological pulsatile flow on flow reversal zone and hemodynamic wall parameters [wall shear stress (WSS) and oscillatory shear index (OSI)] for the progression of the disease, atherosclerosis has been investigated numerically. The governing equations have been solved by finite volume method. For the numerical analysis, Womersley number, Reynolds number and percentage of restriction are taken as 10, 100 and 50 % respectively. From this study it is revealed that the separated flow from the stenosis strongly depends on inlet flow situations, the maximum chance of deposition of plaque material due to streamline contour is higher at time step t* = 0.75 for simple pulsatile flow and at time step t* = 0 for physiological pulsatile flow and this chance is lower in case of steady flow. The effect of WSS on the disease is higher in physiological pulsatile flow compared to steady as well as simple pulsatile flow. The maximum possibility of initiation as well as progression for atherosclerosis in arterial wall due to high WSS takes place at t* = 0.25 for physiological pulsatile flow. OSI indicates same length of separation for two cases of transient flow, but the rate of cyclic departure of WSS is higher in case of physiological pulsatile flow.     

Numerical investigations into the asymmetric effects on the aerodynamic response of a pitching airfoil

The effects of asymmetric sinusoidal motion on pitching airfoil aerodynamics were studied by numerical simulations for 2-D flow around a NACA0012 airfoil at Re=1.35×10⁵. Various unsteady parameters (amplitude of oscillation, d; reduced frequency, k) were applied to investigate the effect of asymmetry parameter S on the instantaneous force coefficients and flow patterns. The results reveal that S has a noticeable effect on the aerodynamic performance, as it affects the instantaneous force coefficient, maximum lift and drag coefficient, hysteresis loops and the flow structures.     

Three-dimensional wake structure measurement using a modified PIV technique

Three-dimensional vortical structures have been measured in a circular-cylinder wake using particle imaging velocimetry (PIV) for the Reynolds number range of 2×10³ to 1×10⁴. The PIV was modified, compared with the conventional one, in terms of its light sheet arrangement to capture reliably streamwise vortices. While in agreement with previous reports, the presently measured spanwise structures complement the data in the literature in the streamwise evolution of the near-wake spanwise vortex in size, strength, streamwise and lateral convection velocities, shedding new light upon vigorous interactions between oppositely signed spanwise structures. The longitudinal vortices display mushroom patterns in the (x, z)-plane in the immediate proximity to the cylinder. Their most likely inclination in the (x, y)-plane is inferred from the measurements in different (x, z)-planes. The longitudinal vortices in the (y, z)-plane show alternate change in sign, though not discernible at x/d > 15. They decay in the maximum vorticity and circulation rapidly from x/d = 5 to 10 and slowly for x/d > 10, and are further compared with the spanwise vortices in size, strength and rate of decay.     

FLUID–STRUCTURE RESONANCE PRODUCED BY ONCOMING ALTERNATING VORTICES

The present investigation examines a simple fluid–structure interaction problem, which is represented by the unsteady response of an airfoil/blade to a Karman vortex street in an inviscid uniform flow. Two different cases were examined; one with a rigid airfoil/blade, where the structural stiffness is infinite, another with an elastic blade. In both cases, the flow remains attached to the airfoil/blade surface. A time-marching technique solving the Euler equations and a two-degree-of-freedom structural dynamic model is used to examine the interactions between the fluid and the structure. The interactions between the convected vortices and the structure modify the shed wake whose energy, in turn, feeds into the forces and moments acting on the structure. For a rigid airfoil/blade, it is found that the amplitude of the aerodynamic response is not proportional to the density of the oncoming vortex street, but depends on c/d , the ratio of the chord length (c) to the axial spacing (d) of the convected vortices. When the number of vortices per unit length is increased, the amplitudes of the aerodynamic response increase and then decrease even though the density of the vorticity keeps increasing and so is the energy of the excitation wake. Maxima are observed at c/d=0·5, 1·5 and 2·5. This behaviour is analogous to the structural resonance phenomenon and is labeled “aerodynamic resonance”. The existence of such an “aerodynamic resonance” is important to turbomachinery applications where the blade is elastic, the flow is unsteady and the shed vortices from the previous row are convected downstream by the mean flow. Thus, “aerodynamic resonance” alone or in conjunction with structural resonance could impact negatively on the fatigue life of turbine blades and their combined effects should be accounted for in blade design. A preliminary attempt to assess this impact has been carried out. It is found that the relative fatigue life of a blade could be reduced by four orders of magnitude as a result.     

Landslide generated impulse waves. 2. Hydrodynamic impact craters

Landslide generated impulse waves were investigated in a two-dimensional physical laboratory model based on the generalized Froude similarity. Digital particle image velocimetry (PIV) was applied to the landslide impact and wave generation. Areas of interest up to 0.8 m by 0.8 m were investigated. PIV provided instantaneous velocity vector fields in a large area of interest and gave insight into the kinematics of the wave generation process. Differential estimates such as vorticity, divergence, and elongational and shear strain were extracted from the velocity vector fields. At high impact velocities flow separation occurred on the slide shoulder resulting in a hydrodynamic impact crater, whereas at low impact velocities no flow detachment was observed. The hydrodynamic impact craters may be distinguished into outward and backward collapsing impact craters. The maximum crater volume, which corresponds to the water displacement volume, exceeded the landslide volume by up to an order of magnitude. The water displacement caused by the landslide generated the first wave crest and the collapse of the air cavity followed by a run-up along the slide ramp issued the second wave crest. The extracted water displacement curves may replace the complex wave generation process in numerical models. The water displacement and displacement rate were described by multiple regressions of the following three dimensionless quantities: the slide Froude number, the relative slide volume, and the relative slide thickness. The slide Froude number was identified as the dominant parameter.List of symbols a wave amplitude (L) - b slide width (L) - c wave celerity (LT^–1) - d _g granulate grain diameter (L) - d _p seeding particle diameter (L) - F slide Froude number - g gravitational acceleration (LT^–2) - h stillwater depth (L) - H wave height (L) - l _s slide length (L) - L wave length (L) - M magnification - m _s slide mass (M) - n _por slide porosity - Q _d water displacement rate (L³) - Q _D maximum water displacement rate (L³) - Q _s maximum slide displacement rate - s slide thickness (L) - S relative slide thickness - t time after impact (T) - t _D time of maximum water displacement volume (L³) - t _qD time of maximum water displacement rate (L³) - t _si slide impact duration (T) - t _sd duration of subaqueous slide motion (T) - T wave period (T) - v velocity (LT^–1) - v _p particle velocity (LT^–1) - v _px streamwise horizontal component of particle velocity (LT^–1) - v _pz vertical component of particle velocity (LT^–1) - v _s slide centroid velocity at impact (LT^–1) - V dimensionless slide volume - V _d water displacement volume (L³) - V _D maximum water displacement volume (L³) - V _s slide volume (L³) - x streamwise coordinate (L) - z vertical coordinate (L) - slide impact angle (°) - bed friction angle (°) - x mean particle image x-displacement in interrogation window (L) - _x random displacement x error (L) - _tot total random velocity v error (LT^–1) - _xx streamwise horizontal elongational strain component (1/T) - _xz shear strain component (1/T) - _zx shear strain component (1/T) - _zz vertical elongational strain component (1/T) - water surface displacement (L) - density (ML^–3) - _g granulate density (ML^–3) - _p particle density (ML^–3) - _s mean slide density (ML^–3) - _w water density (ML^–3) - granulate internal friction angle (°) - _y vorticity vector component (out-of-plane) (1/T)     

用流体双折射方法研究人工心脏瓣膜的脉动流场

本文选用1.5%的MillingYellow溶液,在对二维心瓣低频脉动流场的分析中,仍采用稳态流场的线性光-力学关系。通过摄象机记录下了一个脉动周期内脉动流量Q从36升/小时至55升/小时流场的变化情况,直观地看到主动脉瓣流场中最大切变率的瞬态分布。记录了同一模型在相同条件下稳态流场流量Q为50升/小时的流场,对两种流场下几个典型截面处的最大切变率的数值做了比较。通过对图线的分析,表明在稳态流量与脉动流量相同的情况下,脉动流场的最大切变率数值高于稳态流场,这就说明脉动使流场的剪切作用增强。     

Experimental investigation of the interaction between a plane wall jet and a parallel offset jet

The dual-jet flow generated by a plane wall jet and a parallel offset jet at an offset ratio of d/w = 1.0 has been investigated using Particle Image Velocimetry (PIV). The particle images are captured, processed, and subsequently used to characterize the flow in terms of the 2D velocity and vorticity distributions. Statistical characteristics of the flow are obtained through ensemble averaging of 360 instantaneous velocity fields. Also presented is a time series of instantaneous flow fields to illustrate the dynamic interaction between the two jets. Results reveal that the near field of the flow is characterized by a periodic large-scale Karman-like vortex shedding similar to what would be expected in the wake of a bluff body. The existence of the Karman-like vortices results in periodic interactions between the two jets; in addition, these vortices produce noticeable impact on the jet outer layers, i.e., the free shear layer of the offset jet and the wall boundary layer of the wall jet. A schematic of vortex/shear layer interaction is proposed to illustrate the flow pattern.