Air-water flows down stepped chutes: turbulence and flow structure observations

Interactions between turbulent waters and atmosphere may lead to strong air-water mixing. This experimental study is focused on the flow down a staircase channel characterised by very strong flow aeration and turbulence. Interfacial aeration is characterised by strong air-water mixing extending down to the invert. The size of entrained bubbles and droplets extends over several orders of magnitude, and a significant number of bubble/droplet clusters was observed. Velocity and turbulence intensity measurements suggest high levels of turbulence across the entire air-water flow. The increase in turbulence levels, compared to single-phase flow situations, is proportional to the number of entrained particles. (C) 2002 Elsevier Science Ltd. All rights reserved.

Dynamic fluid-structure interactions using finite volume unstructured mesh procedures

A three-dimensional finite volume, unstructured mesh (FV-UM) method for dynamic fluid–structure interaction (DFSI) is described. Fluid structure interaction, as applied to flexible structures, has wide application in diverse areas such as flutter in aircraft, wind response of buildings, flows in elastic pipes and blood vessels. It involves the coupling of fluid flow and structural mechanics, two fields that are conventionally modelled using two dissimilar methods, thus a single comprehensive computational model of both phenomena is a considerable challenge. Until recently work in this area focused on one phenomenon and represented the behaviour of the other more simply. More recently, strategies for solving the full coupling between the fluid and solid mechanics behaviour have been developed. A key contribution has been made by Farhat et al. [Int. J. Numer. Meth. Fluids 21 (1995) 807] employing FV-UM methods for solving the Euler flow equations and a conventional finite element method for the elastic solid mechanics and the spring based mesh procedure of Batina [AIAA paper 0115, 1989] for mesh movement. In this paper, we describe an approach which broadly exploits the three field strategy described by Farhat for fluid flow, structural dynamics and mesh movement but, in the context of DFSI, contains a number of novel features: a single mesh covering the entire domain, a Navier–Stokes flow, a single FV-UM discretisation approach for both the flow and solid mechanics procedures, an implicit predictor–corrector version of the Newmark algorithm, a single code embedding the whole strategy.

Bubbly Flow Structure in Hydraulic Jump

In an open channel, a hydraulic jump is the rapid transition from super- to sub-critical flow associated with strong turbulence and air bubble entrainment in the mixing layer. New experiments were performed at relatively large Reynolds numbers using phase-detection probes. Some new signal analysis provided characteristic air-water time and length scales of the vortical structures advecting the air bubbles in the developing shear flow. An analysis of the longitudinal air-water flow structure suggested little bubble clustering in the mixing layer, although an interparticle arrival time analysis showed some preferential bubble clustering for small bubbles with chord times below 3 ms. Correlation analyses yielded longitudinal air-water time scales Txx*V1/d1 of about 0.8 in average. The transverse integral length scale Z/d1 of the eddies advecting entrained bubbles was typically between 0.25 and 0.4, irrespective of the inflow conditions within the range of the investigations. Overall the findings highlighted the complicated nature of the air-water flow

Modeling the inelastic behavior of heat exchangers accounting for fluid-structure interactions

Magdeburg, Univ., Fak. für Maschinenbau, Diss., 2013

High-resolution numerical modelling of flow-vegetation interactions

In this paper, we present and apply a new three-dimensional model for the prediction of canopy-flow and turbulence dynamics in open-channel flow. The approach uses a dynamic immersed boundary technique that is coupled in a sequentially staggered manner to a large eddy simulation. Two different biomechanical models are developed depending on whether the vegetation is dominated by bending or tensile forces. For bending plants, a model structured on the Euler-Bernoulli beam equation has been developed, whilst for tensile plants, an N-pendula model has been developed. Validation against flume data shows good agreement and demonstrates that for a given stem density, the models are able to simulate the extraction of energy from the mean flow at the stem-scale which leads to the drag discontinuity and associated mixing layer.

An experimental study of discharge partitioning and flow structure at symmetrical bifurcations

Recent research has examined the factors controlling the geometrical configuration of bifurcations, determined the range of stability conditions for a number of bifurcation types and assessed the impact of perturbations on bifurcation evolution. However, the flow division process and the parameters that influence flow and sediment partitioning are still poorly characterized. To identify and isolate these parameters, three-dimensional velocities were measured at 11 cross-sections in a fixed-walled experimental bifurcation. Water surface gradients were controlled, and systematically varied, using a weir in each distributary. As may be expected, the steepest distributary conveyed the most discharge ( was dominant) while the mildest distributary conveyed the least discharge ( was subordinate). A zone of water surface super-elevation was co-located with the bifurcation in symmetric cases or displaced into the subordinate branch in asymmetric cases. Downstream of a relatively acute-angled bifurcation, primary velocity cores were near to the water surface and against the inner banks, with near-bed zones of lower primary velocity at the outer banks. Downstream of an obtuse-angled bifurcation, velocity cores were initially at the outer banks, with near-bed zones of lower velocities at the inner banks, but patterns soon reverted to match the acute-angled case. A single secondary flow cell was generated in each distributary, with water flowing inwards at the water surface and outwards at the bed. Circulation was relatively enhanced within the subordinate branch, which may help explain why subordinate distributaries remain open, may play a role in determining the size of commonly-observed topographic features, and may thus exert some control on the stability of asymmetric bifurcations. Further, because larger values of circulation result from larger gradient disadvantages, the length of confluence-diffluence units in braided rivers or between diffluences within delta distributary networks may vary depending upon flow structures inherited from upstream and whether, and how, they are fed by dominant or subordinate distributaries. Copyright (C) 2011 John Wiley & Sons, Ltd.

An experimental investigation of the effects of nozzle ellipticity on the flow structure of co-flow jet diffusion flames

The flow structure of cold and ignited jets issuing into a co-flowing air stream was experimentally studied using a laser Doppler velocimeter. Methane was employed as the jet fluid discharging from circular and elliptic nozzles with aspect ratios varying from 1.29 to 1.60. The diameter of the circular nozzle was 4.6 mm and the elliptic nozzles had approximately the same exit area as that of the circular nozzle. These non-circular nozzles were employed in order to increase the stability of attached jet diffusion flames. The time-averaged velocity and r.m.s. value of the velocity fluctuation in the streamwise and transverse directions were measured over the range of co-flowing stream velocities corresponding to different modes of flame blowout that are identified as either lifted or attached flames. On the basis of these measurements, attempts were made to explain the existence of an apparent optimum aspect ratio for the blowout of attached flames observed at higher values of co-flowing stream velocities. The insensitivity of the blowout limits of lifted flames to nozzle geometry observed in our previous work at low co-flowing stream velocities was also explained. Measurements of the fuel concentration at the jet centerline indicated that the mixing process was enhanced with the 1.38 aspect ratio jet compared with the 1.60 aspect ratio jet. On the basis of the obtained experimental data, it was suggested that the higher blowout limits of attached flames for an elliptic jet of 1.38 aspect ratio was due to higher entrainment rates.

Flow structure and variability in the subtropical Indian Ocean: Instability of the South Indian Ocean Countercurrent

The origin of the eddy variability around the 25°S band in the Indian Ocean is investigated. We have found that the surface circulation east of Madagascar shows an anticyclonic subgyre bounded to the south by eastward flow from southwest Madagascar, and to the north by the westward flowing South Equatorial Current (SEC) between 15° and 20°S. The shallow, eastward flowing South Indian Ocean Countercurrent (SICC) extends above the deep reaching, westward flowing SEC to 95°E around the latitude of the high variability band. Applying a two-layer model reveals that regions of large vertical shear along the SICC-SEC system are baroclinically unstable. Estimates of the frequencies (3.5–6 times/year) and wavelengths (290–470 km) of the unstable modes are close to observations of the mesoscale variability derived from altimetry data. It is likely then that Rossby wave variability locally generated in the subtropical South Indian Ocean by baroclinic instability is the origin of the eddy variability around 25°S as seen, for example, in satellite altimetry.

Flow structure and near-field dispersion in arrays of building-like obstacles

Dispersion in the near-field region of localised releases in urban areas is difficult to predict because of the strong influence of individual buildings. Effects include upstream dispersion, trapping of material into building wakes and enhanced concentration fluctuations. As a result, concentration patterns are highly variable in time and mean profiles in the near field are strongly non-Gaussian. These aspects of near-field dispersion are documented by analysing data from direct numerical simulations in arrays of building-like obstacles and are related to the underlying flow structure. The mean flow structure around the buildings is found to exert a strong influence over the dispersion of material in the near field. Diverging streamlines around buildings enhance lateral dispersion. Entrainment of material into building wakes in the very near field gives rise to secondary sources, which then affect the subsequent dispersion pattern. High levels of concentration fluctuations are also found in this very near field; the fluctuation intensity is of order 2 to 5.

Stability and structure: interactions in economic growth

Includes bibliography

Deep chlorophyll maximum and upper ocean structure interactions: Case of the Guinea Thermal Dome

Deep Chlorophyll Maximum (DCM) modifies the upper ocean heat capture distribution and thus impacts water column temperature and stratification, as well as biogeochemical processes. This energetical role of the DCM is assessed using a 1 m-resolution 1D physical-biogeochemical model of the upper ocean, using climatological forcing conditions of the Guinea Dome (GD). This zone has been chosen among others because a strong and shallow DCM is present all year round. The results show that the DCM warms the seasonal thermocline by +2 degrees C in September/October and causes an increase of heat transfer from below into the mixed layer (ML) by vertical diffusion and entrainment, leading to a ML warming of about 0.3 degrees C in October. In the permanent thermocline, temperature decreases by up to 2 degrees C. The result is a stratification increase of the water column by 0.3 degrees C m(-1) which improves the thermocline realism when compared with observations. At the same time, the heating associated with the DCM is responsible for an increase of nitrate (+300%, 0.024 mu M), chlorophyll (+50%, 0.02 mu g l(-1)) and primary production (+45%: 10 mg C m(-2) day(-1)) in the ML during the entrainment period of October. The considered concentrations are small but this mechanism could be potentially important to give a better explanation of why there is a significant amount of nitrate in the ML. The mechanisms associated with the DCM presence, no matter which temperature or biogeochemical tracers are concerned, are likely to occur in a wide range of tropical or subpolar regions; in these zones a pronounced DCM is present at least episodically at shallow or moderate depths. These results can be generalized to other thermal dome regions where relatively similar physical and biogeochemical structures are encountered. After testing different vertical resolutions (10 m, 5 m, 2.5 m, 1 m and 0.5 m), we show that using at least a 1 to vertical resolution model is mandatory to assess the energetical importance of the DCM.

Laboratory investigation of flow structure and resistance of high- and low-angle dunes

A prominent control on the flow over subaqueous dunes is the slope of the downstream leeside. While previous work has focused on steep (~30°), asymmetric dunes with permanent flow separation, little is known about dunes with lower lee-slope angles for which flow separation is absent or intermittent. Here, we present a laboratory investigation where we systematically varied the dune lee-slope, holding other geometric parameters and flow hydraulics constant, to explore effects on the turbulent flow field and flow resistance. Three sets of fixed dunes (lee-slopes of 10°, 20° and 30°) were separately installed in a 15 m long and 1 m wide flume and subjected to 0.20 m deep flow. Measurements consisted of high-frequency, vertical profiles collected with a Laser Doppler Velocimeter (LDV). We show that the temporal and spatial occurrence of flow separation decreases with dune lee-slope. Velocity gradients in the dune leeside depict a free shear layer downstream of the 30° dunes and a weaker shear layer closer to the bed for the 20° and 10° dunes. The decrease in velocity gradients leads to lower magnitude of turbulence production for gentle lee-slopes. Aperiodic, strong ejection events dominate the shear layer, but decrease in strength and frequency for low-angle dunes. Flow resistance of dunes decreases with lee-slope; the transition being non-linear. Over the 10°, 20° and 30° dunes, shear stress is 8%, 33% and 90 % greater than a flat bed, respectively. Our results demonstrate that dune lee-slope plays an important, but often ignored role in flow resistance.

Tranversal motion and flow structure of fully nonlinear streaks in a laminar boundary layer

Typical streak computations present in the literature correspond to linear streaks or to small amplitude nonlinear streaks computed using DNS or nonlinear PSE. We use the Reduced Navier-Stokes (RNS) equations to compute the streamwise evolution of fully non-linear streaks with high amplitude in a laminar flat plate boundary layer. The RNS formulation provides Reynolds number independent solutions that are asymptotically exact in the limit $Re \gg 1$, it requires much less computational effort than DNS, and it does not have the consistency and convergence problems of the PSE. We present various streak computations to show that the flow configuration changes substantially when the amplitude of the streaks grows and the nonlinear effects come into play. The transversal motion (in the wall normal-streamwise plane) becomes more important and strongly distorts the streamwise velocity profiles, that end up being quite different from those of the linear case. We analyze in detail the resulting flow patterns for the nonlinearly saturated streaks and compare them with available experimental results.

On the wake-induced vibration of tandem circular cylinders: the vortex interaction excitation mechanism

The mechanism of wake-induced vibrations (WIV) of a pair of cylinders in a tandem arrangement is investigated by experiments. A typical WIV response is characterized by a build-up of amplitude persisting to high reduced velocities; this is different from a typical vortex-induced vibration (VIV) response, which occurs in a limited resonance range. We suggest that WIV of the downstream cylinder is excited by the unsteady vortex-structure interactions between the body and the upstream wake. Coherent vortices interfering with the downstream cylinder induce fluctuations in the fluid force that are not synchronized with the motion. A favourable phase lag between the displacement and the fluid force guarantees that a positive energy transfer from the flow to the structure sustains the oscillations. If the unsteady vortices are removed from the wake of the upstream body then WIV will not be excited. An experiment performed in a steady shear flow turned out to be central to the understanding of the origin of the fluid forces acting on the downstream cylinder.

The role of wake stiffness on the wake-induced vibration of the downstream cylinder of a tandem pair

When a pair of tandem cylinders is immersed in a flow the downstream cylinder can be excited into wake-induced vibrations (WIV) due to the interaction with vortices coming from the upstream cylinder. Assi, Bearman & Meneghini ( J. Fluid Mech. , vol. 661, 2010, pp. 365–401) concluded that the WIV excitation mechanism has its origin in the unsteady vortex–structure interaction encountered by the cylinder as it oscillates across the wake. In the present paper we investigate how the cylinder responds to that excitation, characterising the amplitude and frequency of response and its dependency on other parameters of the system. We introduce the concept of wake stiffness , a fluid dynamic effect that can be associated, to a first approximation, with a linear spring with stiffness proportional to Re and to the steady lift force occurring for staggered cylinders. By a series of experiments with a cylinder mounted on a base without springs we verify that such wake stiffness is not only strong enough to sustain oscillatory motion, but can also dominate over the structural stiffness of the system. We conclude that while unsteady vortex–structure interactions provide the energy input to sustain the vibrations, it is the wake stiffness phenomenon that defines the character of the WIV response