Shallow Flows Over a Permeable Medium: The Hydrodynamics of Submerged Aquatic Canopies

Aquatic flow over a submerged vegetation canopy is a ubiquitous example of flow adjacent to a permeable medium. Aquatic canopy flows, however, have two important distinguishing features. Firstly, submerged vegetation typically grows in shallow regions. Consequently, the roughness sublayer, the region where the drag length scale of the canopy is dynamically important, can often encompass the entire flow depth. In such shallow flows, vortices generated by the inflectional velocity profile are the dominant mixing mechanism. Vertical transport across the canopy–water interface occurs over a narrow frequency range centered around f _v (the frequency of vortex passage), with the vortices responsible for more than three-quarters of the interfacial flux. Secondly, submerged canopies are typically flexible, coupling the motion of the fluid and canopy. Importantly, flexible canopies can exhibit a coherent waving (the monami) in response to vortex passage. This waving reduces canopy drag, allowing greater in-canopy velocities and turbulent stresses. As a result, the waving of an experimental canopy reduces the canopy residence time by a factor of four. Finally, the length required for the set-up and full development of mixing-layer-type canopy flow is investigated. This distance, which scales upon the drag length scale, can be of the same order as the length of the canopy. In several flows adjacent to permeable media (such as urban canopies and reef systems), patchiness of the medium is common such that the fully developed condition may not be representative of the flow as a whole. This paper has previously been published in Transp Porous Med (2009) 78: 309–326; DOI .     

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     

PIV investigation of flow behind surface mounted permeable ribs

The flow behind surface mounted permeable rib geometries, i.e. solid, slit, split-slit and inclined split-slit ribs have been studied using flow visualization and PIV (2-C and 3-C) technique in streamwise and cross-stream measurement planes. The objective behind this study is to understand the flow structures responsible for heat transfer/mixing enhancement with simultaneous pressure penalty reduction by permeable rib geometries. The Reynolds number based on the rib height has been set equal to 5,538 and the open area ratio of permeable ribs is equal to 20%. The permeable rib geometries have shorter reattachment length in comparison to the solid rib. The maximum 41% reduction in reattachment length is observed for the inclined split-slit rib. The splitter mounted inside the slit leads to two corner vortices behind it. The corner vortices drag the flow from the primary recirculation bubble region towards the rib resulting in drop of the reattachment length. Two horseshoe vortices are present in the flow through the slit at both sides of the splitter due to the upstream flow separation. The slit inclination moves these horseshoe vortices closer to the bottom wall. A film like flow through the slit is present near the downstream corner of the inclined split-slit rib. The spanwise velocity gradient due to the splitter leads to vorticity and turbulence enhancement by vortex stretching. The inclination of the slit and the use of a splitter inside the slit are two important design parameters responsible in generation of near-wall longitudinal vortices. The flow field behind permeable ribs is dominated by vortical structures with definable critical flow patterns, i.e. node, saddle and foci. These predominant swirling flow motions contribute to the mixing enhancement behind permeable rib geometries. On leave from Mechanical Engineering Department, IIT Kanpur, U.P. 208016, India     

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.     

Flow structures and force characteristics for flat plate in oscillatory flows withKc number from 2 to 40 and in combined flows

The evolution of wake structures and variation of the forces on a flat plate in harmonic oscillatory and in-line combined flows are obtained numerically by improved discrete vortex method. For the oscillatory oncoming flow cases, wyenKc number varies from 2 to 40, the vortex pattern changes from a “harmonic wave” shaped (in a range of smallKc number) to a slight inclined “harmonic wave” shaped (in a range of moderateKc numbers), then to inclined vortex clusters with an angle of 50° to the oncoming flow direction (atKc=20), at last, asKc number becomes large, the vortex pattern is like a normal Karman vortex street. The well predicted drag and inertia force coefficients are obtained, which are more close to the results of Keulegan & Carpenter's experiment as compared with previous vortex simulation by other authors. The existence of minimum point of inertia force coefficientC _m nearKc=20 is also well predicted and this phenomenon can be interpreted according to the vortex structure. For steady-oscillatory in-line combined flow cases, the vortex modes behave like a vortex street, exhibit a “longitudinal wave” structure, and a vortex cluster shape corresponding to the ratios ofU _m toU ₀ which are ofO (10⁻¹)O(1) andO(10), respectively. The effect on the prediction of forces on the flat plate from the disturbance component in a combined flow has been demonstrated qualitatively. In addition to this, the lock in phenomenon of vortex shedding has been checked. The project supported by National Natural Science Foundation of China & LNM, Institute of Mechanics, CAS     

Large aerodynamic forces on a sweeping wing at low Reynolds number

The aerodynamic forces and flow structure of a model insect wing is studied by solving the Navier-Stokes equations numerically. After an initial start from rest, the wing is made to execute an azimuthal rotation (sweeping) at a large angle of attack and constant angular velocity. The Reynolds number (Re) considered in the present note is 480 (Re is based on the mean chord length of the wing and the speed at 60% wing length from the wing root). During the constant-speed sweeping motion, the stall is absent and large and approximately constant lift and drag coefficients can be maintained. The mechanism for the absence of the stall or the maintenance of large aerodynamic force coefficients is as follows. Soon after the initial start, a vortex ring, which consists of the leading-edge vortex (LEV), the starting vortex, and the two wing-tip vortices, is formed in the wake of the wing. During the subsequent motion of the wing, a base-to-tip spanwise flow converts the vorticity in the LEV to the wing tip and the LEV keeps an approximately constant strength. This prevents the LEV from shedding. As a result, the size of the vortex ring increases approximately linearly with time, resulting in an approximately constant time rate of the first moment of vorticity, or approximately constant lift and drag coefficients. The variation of the relative velocity along the wing span causes a pressure gradient along the wingspan. The base-to-tip spanwise flow is mainly maintained by the pressure-gradient force. The project supported by the National Natural Science Foundation of China (10232010)     

Experimental study on the effect of pulsating inflow to an enclosure for improved mixing

Optimal control of inlet jet flows is of broad interest for enhanced mixing in ventilated rooms. The general approach in mechanical ventilation is forced convection by means of a constant flow rate supply. However, this type of ventilation may cause several problems such as draught and appearance of stagnation zones, which reduces the ventilation efficiency. A potential way to improve the ventilation quality is to apply a pulsating inflow, which has been hypothesised to reduce the stagnation zones due to enhanced mixing. The present study aims at testing this hypothesis, experimentally, in a small-scale two-dimensional water model using Particle Image Velocimetry with an in-house vortex detection program. We are able to show that for an increase in pulsation frequency or alternatively in the flow rate the stagnation zones are reduced in size and the distribution of vortices becomes more homogeneous over the considered domain. The number of vortices (all scales) increases by a factor of four and the swirl-strength by about 50% simply by turning on the inflow pulsation. Furthermore, the vortices are well balanced in terms of their rotational direction, which is validated by the symmetric Probability Density Functions of vortex circulation (Γ) around Γ = 0. There are two dominating vortex length scales in the flow, namely 0.6 and 0.8 inlet diameters and the spectrum of vortex diameters become broader by turning on the inflow pulsation. We conclude that the positive effect for enhanced mixing by increasing the flow rate can equally be accomplished by applying a pulsating inflow.     

Effects of the corner radius on the near wake of a square prism

The near wake of square cylinders with different corner radii was experimentally studied based on particle imaging velocimetry (PIV), laser doppler anemometry (LDA) and hotwire measurements. Four bluff bodies, i.e., r/d=0 (square cylinder), 0.157, 0.236, 0.5 (circular cylinder), where r is corner radius and d is the characteristic dimension of the bluff bodies, were examined. A conditional sampling technique was developed to obtain the phase-averaged PIV data in order to characterize quantitatively the effect of corner radii on the near-wake flow structure. The results show that, as r/d increases from 0 to 0.5, the maximum strength of shed vortices attenuates, the circulation associated with the vortices decreases progressively by 50%, the Strouhal number, St, increases by about 60%, the convection velocity of the vortices increases along with the widening of the wake width by about 25%, the vortex formation length and the wake closure length almost double in size. Meanwhile, both the vortex wavelength, λ _x , and the lateral spacing, λ _y , decrease as r/d increases, but the ratio of λ _y to λ _x is approximately 0.29, irrespective of r/d, which is close to the theoretical value of 0.281 for a stable Karman vortex street. The decrease in wavelength is probably responsible for the change in the flow structure from the approximately circular-shaped vortex at r/d=0 to the laterally stretched vortex at r/d=0.5. The leading edge corner radius is more important than the trailing one in influencing the near wake structure since it determines to a great extent the behavior of the streamlines, the separation angle and the base pressure. It is further found that the ratio of the mean drag coefficient to the total shed circulation, C _d/Γ₀, approaches a constant, about 0.25 for different bluff bodies in the subcritical flow regime. The streamwise evolution of vortices and the streamwise fluctuating velocity along the centerline for rounded cylinders are also discussed.     

Large Eddy Simulation of Turbulent Flow Through Submerged Vegetation

Large Eddy Simulations (LES) are performed for an open channel flow through idealized submerged vegetation with a water depth (h) to plant height (h _p) ratio of h/h _p = 1.5 according to the experimental configuration of Liu et al. (J Geophys Res Earth Sci, 2008). They used a 1D laser Doppler velocimeter (LDV) to measure longitudinal and vertical velocities as well as turbulence intensities along several verticals in the flow and the data are used for the validation of the present simulations. The code MGLET is used to solve the filtered Navier–Stokes equations on a Cartesian non-uniform grid. In order to represent solid objects in the flow, the immersed boundary method is employed. The computational domain is idealized with a box containing 16 submerged circular cylinders and periodic boundary conditions are applied in both longitudinal and transverse directions. The predicted streamwise as well as vertical mean velocities are in good agreement with the LDV measurements. Furthermore, fairly good agreement is found between calculated and measured streamwise and vertical turbulence intensities. Large-scale flow structures of different shapes are present in the form of vortex rolls above the vegetation tops as well as locally generated trailing and von- Karman-type vortices due to flow separation at the free end and the sides of the cylinders. In this paper, the flow field is analyzed statistically and evidence is provided for the existence of these structures based on the LES.     

A study of the aerodynamics of a generic container freight wagon using Large-Eddy Simulation

In this work simulations using the Large Eddy Simulation technique have been made of the flow around a generic container freight wagon model. The model consists of one 11.8 m standard length container placed on a wagon. Details of the undercarriage such as wheels are included, but the container is generic and smoothed in comparison to a real freight wagon. The Reynolds number of the flow is 10⁵ based on the container width of 2.354 m. Two cases have been considered in the study, one case where the wagon is standing alone and one case where it is submerged into a train set with wagons ahead and behind the wagon. The latter case is simulated using periodic boundary condition. Both the time-averaged and the instantaneous flow around the wagon for the two cases are described. For the single wagon case, it is found that the separation bubble formed on the roof of the container oscillates back and forth in the streamwise direction and that this oscillation is in phase with oscillations found in the upper shear layer of the ring vortex in the wake. The mechanism that is causing the synchronization of the oscillations of the separation bubble at the front and the upper shear layers in the wake is found to be waves of vorticity being shed from the separation bubble. The time-averaged ring vortex in the near wake of the single wagon is found to be inclined due to the disturbance of the undercarriage details on flow in the lower shear layer. The lower center of the ring vortex is located closer to the base face than the upper center. The drag coefficient of the wagon in the periodic case was found to be only 10% of that of the single wagon case. This is due to two symmetrical counter-rotating vortices found in the gaps which make the train set appear as a single body to the oncoming flow and shielding the wagon from any direct impingement of the flow. The counter-rotating vortices in the gap are found to inhibit periodic oscillations in the lateral direction. These oscillations cause vortical structures to form by the air that is pushed out from the gap and these flow structures cause a dominating oscillation of non-dimensional frequency St=0.12 in the side force signal.     

Wall proximity effects on the effectiveness of upstream control rod

Two dimensional flow over a circular cylinder with an upstream control rod of same diameter is simulated in unbound condition and in wall bounded conditions. The cylinders are placed at various heights from the wall and the inter-distance between cylinders is also varied. The control rod is subjected to different rotation rates. It is found that, in unbound condition, rotating the control rod decreases the critical pitch length (S/D_cr) and increases the drag and Strouhal number of the main cylinder. In presence of plane wall, the shielding provided by the separated shear layers from the control rod in cavity regime is deteriorated due to deflection of shear layers which results in higher drag and large fluctuation of lift coefficient. However, in wake impingement regime, the binary vortices from the control rod are weakened due to diffusion of vorticity and hence, the main cylinder experiences a lower drag and small lift fluctuations than that of unbound condition. The critical height of vortex suppression (H/D_cr) is higher in cavity regime than that of wake impingement regime due to the single extended-bluff body like configuration. The rotation of control rod energizes the wall boundary layer and increases the critical height of vortex suppression. Increasing the rotational rate of control rod decreases the drag force and reduces the amplitude of lift fluctuation. Analysis of the wall shear stress distribution reveals that it suffers a sudden drop at moderate height where the normal Karman vortex shedding changes to irregular shedding consisting of single row of negative vortices. Modal structures obtained from dynamic mode decomposition (DMD) reveal that the flow structures behind the main cylinder are suppressed due to wall and the flow is dominated by the wake of control rod.     

Stereoscopic PIV measurement of unsteady flows in an axial compressor stage

Stereoscopic particle-image velocimetry (SPIV) measurements are conducted in a Low-speed Large-scale Axial Compressor. During the experiment the two CCD cameras are placed at the different sides of the laser light sheet and it is proved that this configuration is more suitable for the investigation in multi-stage turbomachines. The measured results, including the overall performances of many typical flow structures near the rotor tip region and the phase locked unsteady flows inside the stator passage at both the design and near-stall conditions, are introduced. Some new features of the complicated flow structures, such as the breakdown of the tip leakage vortex, the formation of the compound corner vortex at the rotor suction tip corner, the interactions between the hub stall and the tip separation and the rotor wakes, and the evolutions of the tip corner anti-rotating streamwise vortices inside the stator passage, are revealed.     

Reynolds number effects on leading edge vortex development on a waving wing

The waving wing experiment is a fully three-dimensional simplification of the flapping wing motion observed in nature. The spanwise velocity gradient and wing starting and stopping acceleration that exist on an insect-like flapping wing are generated by rotational motion of a finite span wing. The flow development around a waving wing at Reynolds number between 10,000 and 60,000 has been studied using flow visualization and high-speed PIV to capture the unsteady velocity field. Lift and drag forces have been measured over a range of angles of attack, and the lift curve shape was similar in all cases. A transient high-lift peak approximately 1.5 times the quasi-steady value occurred in the first chord length of travel, caused by the formation of a strong attached leading edge vortex. This vortex appears to develop and shed more quickly at lower Reynolds numbers. The circulation of the leading edge vortex has been measured and agrees well with force data.     

Three dimensional flow at the junction between a turbine blade and end-wall

A visualization of the flow on the suction side and end-wall of a passage between two neighboring turbine blades is compared with mass (heat) transfer measurements on the same surfaces. Besides the horseshoe and passage vortices, there are several smaller vortices formed near the junction of blade and end-wall whose origins are discussed. The vortices detach from the end-wall and move up the blade's span. These vortices, sometimes in counter rotating pairs, are responsible for substantial local variations of heat transfer.

---

Drei-dimensionale Strömung an der Naht zwischen Turbinenschaufel und Endwand

Zusammenfassung Der Durchfluß auf der Saugseite und an der Endwand eines Kanals zwischen zwei benachbarten Turbinenschaufeln wird sichtbar gemacht und mit gemessenen Massen-(Wärme)strömen der selben Oberflächen verglichen. Neben den hufförmigen Wirbeln und Durchflußwirbeln werden mehrere kleinere Wirbel in der Nähe der Verbindungsstelle von Schaufel und Endwand gebildet und die Ursache ihrer Entstehung wird diskutiert. Die Wirbel lösen sich von der Endwand ab und bewegen sich über die Spannweite der Schaufel. Diese Wirbel, manchmal als gegenläufige rotierende Paare, sind für die wesentlichen lokalen Variationen des Wärmeflusses verantwortlich.

---

Nomenclature AT Attachment line between counter-rotating vortices - C Chord length of the test blade, c.f. Fig. 1 b;C=16.91 cm in the present study - DE Detachment line between pairs of counter-rotating vortices - DL Flow dividing line-line along which surface flow divides into neighboring blade passages - Re ₁ Reynolds number based on cascade inlet velocity and blade chord length,U ₁ C/ - S _s Suction side curvilinear distance from the stagnation line of the test blade (see Fig. 1 b). Note thatS _s /C=1.355 at the trailing edge of the blades in the present study - U ₁ inlet velocity to the cascade - V _h the horseshoe vortex - V _Lc the leading edge corner vortex - V _p the passage vortex - V _ph pressure side leg of the horseshoe vortex - V _pLc pressure side leg of the leading edge corner vortex - V _s1 leading edge stagnation region vortex 1 - V _s2 leading edge stagnation region vortex 2 - V _sc1 suction side corner vortex 1 - V _sc1s portion ofV _sc1 which climbs up the blade suction surface - V _sc2 suction side corner vortex 2 - V _sc3 suction side corner vortex 3 - V _sh suction side leg of the horseshoe vortex - V _sLc suction side leg of the leading edge corner vortex - kinematic viscosity Dedicated to Prof. Dr.-Ing. F. Mayinger's 60th birthday     

Time resolved investigations on flow field and quasi wall shear stress of an impingement configuration with pulsating jets by means of high speed PIV and a surface hot wire array

The effects of jet pulsation on flow field and quasi wall shear stress of an impingement configuration were investigated experimentally. The excitation Strouhal number and amplitude were varied as the most influential parameters. A line-array with three submerged air jets, and a confining plate were used. The flow field analysis by means of time resolved particle image velocimetry shows that the controlled excitation can considerably affect the near-field flow of an impinging jet array. These effects are visualized as organization of the coherent flow structures. Augmentation of the Kelvin–Helmholtz vortices in the jet shear layer depends on the Strouhal number and pulsation magnitude and can be associated with pairing of small scale vortices in the jet. A total maximum of vortex strength was observed when exciting with Sr = 0.82 and coincident high amplitudes.Time resolved interaction between impinging vortices and impingement plate boundary layer due to jet excitation was verified by using an array of 5 μm surface hot wires. Corresponding to the global flow field modification due to periodic jet pulsation, the impact of the vortex rings on the wall boundary layer is highly influenced by the above mentioned excitation parameters and reaches a maximum at Sr = 0.82.     

LES and RANS for Turbulent Flow over Arrays of Wall-Mounted Obstacles

Large-eddy simulation (LES) has been applied to calculate the turbulent flow over staggered wall-mounted cubes and staggered random arrays of obstacles with area density 25%, at Reynolds numbers between 5 × 10³ and 5 10⁶, based on the free stream velocity and the obstacle height. Re = 5 × 10³ data were intensively validated against direct numerical simulation (DNS) results at the same Re and experimental data obtained in a boundary layer developing over an identical roughness and at a rather higher Re. The results collectively confirm that Reynolds number dependency is very weak, principally because the surface drag is predominantly form drag and the turbulence production process is at scales comparable to the roughness element sizes. LES is thus able to simulate turbulent flow over the urban-like obstacles at high Re with grids that would be far too coarse for adequate computation of corresponding smooth-wall flows. Comparison between LES and steady Reynolds-averaged Navier-Stokes (RANS) results are included, emphasising that the latter are inadequate, especially within the canopy region.     

Development and stability of swirling flows

The purpose of this research is to investigate steady axisymmetric swirling flows in channels and free vortices and also to establish the role of hydrodynamic instability in the formation of the sharp changes in flow structure associated with an increased rate of rotation. On the basis of numerical solutions of the complete Navier-Stokes equations obtained by a finite-difference method swirling flows in pipes with impermeable and permeable walls and in a free vortex are investigated. The stability of the swirling axisymmetric flows is considered on the assumption of local parallelism: the problem of the normal modes developing against the background of the axisymmetric flow determined by the velocity profiles in local cross sections of the flow is solved. Attention is mainly concentrated on free vortex flows with reverse current zones, their structure and stability.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 3–11, July–August, 1988.     

The lowest stability and bifurcation in supercritical Taylor vortices

Here, we numerically investigate the lowest stability and bifurcation boundary of supercritical Taylor vortices in flows with different wavenumbers and for various radius ratios; the radius ratios range from those corresponding to axisymmetrical Taylor vortex flow (TVF) to those corresponding to wavy vortex flow (WVF). The variation in the wavenumber of a supercritical TVF is found to affect the stability of the flow, because the wavenumber of the Taylor vortices remains constant only when the flow is quasi-static. The variation in the wavenumber is examined and found to be significant when the radius ratio is less than 0.7842. The results for TVF are compared with those for the flow during the quasi-static transition from TVF to WVF.     

Vortex shedding around a near-wall circular cylinder induced by a solitary wave

This study developed a two-dimensional generalized vortex method to analyze the shedding of vortices and the hydrodynamic forces resulting from a solitary wave passing over a submerged circular cylinder placed near a flat seabed. Numerical results for validation are compared with other numerical and experimental results, and satisfactory agreement is found. A series of simulations were performed to study the effects of gap-to-diameter ratio and incident wave height on vorticity pattern as well as the forces exerted on the cylinder. The range of the heights of incident waves is from 0.3h to 0.7h, where h is the still water depth. The range of the gap-to-diameter ratios is from 0.1 to 0.8. The results indicate that the flow pattern and the pressure distribution change significantly because of the close proximity of the seabed where the vorticity flux on the seabed-side surface of the cylinder is suppressed. Placing the cylinder nearer the seabed increases the drag and the positive lift on the cylinder. When the gap-to-diameter ratio increases, the pattern of vortices changes because of the interaction between the main recirculation zone and the shear layers separated from the gap. The maxima of drag, lift and total force increase linearly with the height of the incident wave.