Microscale Synthetic Schlieren

We develop the axisymmetric Synthetic Schlieren technique to study the wake of a microscale sphere settling through a density stratification. A video-microscope was used to magnify and image apparent displacements of a micron-sized random-dot pattern. Due to the nature of the wake, density gradient perturbations in the horizontal greatly exceed those in the vertical, requiring modification of previously developed axisymmetric techniques. We present results for 780 and 383 μm spheres, and describe the limiting role of noise in the system for a 157 μm sphere. This technique can be instrumental in understanding a range of ecological and environmental oceanic processes on the microscale.

---

A random synthetic jet array driven turbulence tank

We measure the flow above an array of randomly driven, upward-facing synthetic jets used to generate turbulence beneath a free surface. Compared to grid stirred tanks (GSTs), this system offers smaller mean flows at equivalent turbulent Reynolds numbers with fewer moving parts.

---

Two- and three-dimensional flows in nearly rectangular cavities driven by collinear motion of two facing walls

Two- and three-dimensional flows in nearly cuboidal cavities are investigated experimentally. A tight cavity is formed in the gap between two long and parallel cylinders of large radii by adding rigid top, bottom, and end walls. The cross-section perpendicular to the axes of the cylinders is nearly rectangular with aspect ratio Γ. The axial aspect ratio Λ > 10 is large to suppress end-wall effects. The fluid motion is driven by independent and steady rotation of the cylinders about their axes which defines two Reynolds numbers Re _1,2. Stability boundaries of the nearly two-dimensional steady flow have been determined as functions of Re _1,2 for Γ = 0.76 and Γ = 1. Up to six different three-dimensional supercritical modes have been identified. The critical thresholds for the onset of most of the three-dimensional modes, three of which have been observed for the first time, agree well with corresponding linear-stability calculations. Particular attention is paid to the flow for Γ = 1 under symmetric and parallel wall motion. In that case the basic flow consists of two mirror symmetric counter-rotating parallel vortices. They become modulated in span-wise direction as the driving increases. Detailed LDV measurements of the supercritical three-dimensional velocity field and the bifurcation show an excellent agreement with numerical simulations.

---

Quantitative planar Raman imaging through a spectrograph: visualisation of a supersonic wedge flow

Planar Raman imaging through a spectrograph is demonstrated as a diagnostic tool for quantitative flow visualisation of internal supersonic wedge flow. A dedicated Bayesian deconvolution filter is used to remove the spectral structure that is introduced by the spectrograph. The 2D density field is determined with ca. 10% precision using average images over 6,000 laser pulses, down to 0.5 mm from the surface of the wedge. Direct interpretations of Raman intensities provide more precise density data than indirect interpretations based on shock geometry in 2D inviscid flow.

---

Turbulence measurements in a swirling pipe flow

This paper reports laser-Doppler measurements of the mean flow and turbulence stresses in a swirling pipe flow. Experiments were carried out under well-controlled laboratory conditions in a refractive index-matched pipe flow facility. The results show pronounced asymmetry in mean and fluctuating quantities during the downstream decay of the swirl. Experimental data reveal that the swirl significantly modifies the anisotropy of turbulence and that it can induce explosive growth of the turbulent kinetic energy during its decay. Anisotropy invariant mapping of the turbulent stresses shows that the additional flow deformation imposed by initially strong swirling motion forces turbulence in the core region to tend towards the isotropic two-component state. When turbulence reaches this limiting state it induces rapid production of turbulent kinetic energy during the swirl decay.

---

Rosetta Stone Experiments

Dynamic failure events such as armor penetration and explosive fragmentation are too complex to be treated by classical single-crack continuum fracture mechanics. In such cases deformation and fracture result from multiple cracks, voids, and shear bands acting simultaneously and influencing one another’s evolution. An alternative “meso” fracture mechanics is needed that treats microfailure activity while permitting fast and inexpensive predictive computations. This paper discusses the approach and experiments that elucidate and quantify failure physics on the micron level. “Rosetta Stone” experiments that isolate a damage mode, produce statistical distributions of damage features, and “freeze in” damage at various stages of development are described and illustrated. The observations and data lead to equations describing nucleation and growth of cracks, voids, and shear bands. The resulting mesomechanical material failure models link the microworld with the macroworld and can be used in continuum hydrocodes for fast, efficient simulations of dynamic fracture scenarios.

An application of Gappy POD

An iterative procedure, based on the proper orthogonal decomposition (POD), first proposed by Everson and Sirovich (J Opt Soc Am A 12(8):1657–1664, 1995) is applied to marred particle image velocimetry (PIV) data of shallow rectangular cavity flow at Mach 0.19, 0.28, 0.38, and 0.55. The procedure estimates the POD modes while simultaneously estimating the missing vectors in the PIV data. The results demonstrate that the absolute difference between the repaired vectors and the original PIV data approaches the experimental uncertainty as the number of included POD modes is increased. The estimation of the dominant POD modes is also shown to converge by examining the subspace spanned by the POD eigenfunctions.

---

Double-pulse planar-LIF investigations using fluorescence motion analysis for mixture formation investigation

A concept for dynamic mixture formation investigations of fuel/air mixtures is presented which can equally be applied to several other laser induced fluorescence (LIF) applications. Double-pulse LIF imaging was used to gain insight into dynamic mixture formation processes. The setup consists of a modified standard PIV setup. The "fuel/air ratio measurement by laser induced fluorescence (FARLIF)" approach is used for a quantification of the LIF images in order to obtain pairs of 2D fuel/air ratio maps. Two different evaluation concepts for LIF double pulse images are discussed. The first is based on the calculation of the temporal derivative field of the fuel/air ratio distribution. The result gives insight into the dynamic mixing process, showing where and how the mixture is changing locally. The second concept uses optical flow methods in order to estimate the motion of fluorescence (i.e., mixture) structures to gain insight into the dynamics, showing the distortion and the motion of the inhomogeneous mixture field. For this "fluorescence motion analysis" (FMA) two different evaluation approaches—the "variational gradient based approach" and the "variational cross correlation based approach"—are presented. For the validation of both, synthetic LIF image pairs with predefined motion fields were generated. Both methods were applied and the results compared with the known original motion field. This validation shows that FMA yields reliable results even for image pairs with low signal/noise ratio. Here, the "variational gradient based approach" turned out to be the better choice so far. Finally, the experimental combination of double-pulse FARLIF imaging with FMA and simultaneous PIV measurement is demonstrated. The comparison of the FMA motion field and the flow velocity field captured by PIV shows that both results basically reflect complementary information of the flow field. It is shown that the motion field of the fluorescence structures does not (necessarily) need to represent the actual flow velocity and that the flow velocity field alone can not illustrate the structure motion in any case. Therefore, the simultaneous measurement of both gives the deepest insight into the dynamic mixture formation process. The examined concepts and evaluation approaches of this paper can easily be adapted to various other planar LIF methods (with the LIF signal representing, e.g., species concentration, temperature, density etc.) broadening the insight for a wide range of different dynamic processes.

---

A pragmatic approach for deriving constitutive equations endowed with pom–pom attributes

The most important rheological and mathematical features of the pom–pom model are presently used to compare and improve other constitutive models such as the Giesekus and Phan-Thien–Tanner models. A pragmatic methodology is selected that allows derivation of simple constitutive equations, which are suited to possible software implementation. Alterations to the double convected pom–pom, Phan-Thien–Tanner and Giesekus models are proposed and assessed in rheometric flows by comparing model predictions to experimental data.

Measurement of the gradient field of a turbulent free surface

We study the free surface above a turbulent channel flow. We describe a laser scanning technique that can be used to measure the space–time turbulent surface gradient field along a line. A harmonically swiveling laser beam is focused on the surface and its angle of refraction is measured using a position sensing device. The registered signals can be converted easily to the desired gradient field, and spectra and correlations can be measured. Examples of measured spectra and correlation functions of the surface above a turbulent channel flow (Reynolds number R _λ ≈ 250) demonstrate the viability of the technique. We further assess the validity of Taylor’s frozen turbulence hypothesis that implies that time-dependent signals measured along a line that is oriented perpendicularly to the mean channel velocity can be interpreted as 2D measurements of the surface slope. While Taylor’s hypothesis works for a turbulent velocity field, it does not work for its free surface.

Objective flow classification parameters and their use in general steady flows

Three quantitative flow classification parameters have been studied in the context of Tanner and Huilgol’s suggestion of strong and weak flows. Seen in this context, the different types of streamlines possible for general 3-D flows furnish no indication with respect to the flow strength. This is in total contrast to 2-D flows, where the type of the streamline and the strength of the flow go hand in hand. Astarita’s [J Non-Newton Fluid Mech, 6:69–76, 1979] flow classification parameter takes care of this fact and, if properly generalized, can be applied to more general flows: Two other flow classification parameters also have their basis in homogeneous 2-D flows, but their generalization leads, for general flows, to nonuniqueness and other unacceptable results. For 3-D flows, none of the parameters can quantitatively be used in general, and additional parameters, with their basis outside the 2-D flow regime, seem to be called for.

A model of a subcritical Joule–Thomson cryocooler with condensation inside the recuperator

To develop a tool for predicting of heat and mass transfer in Joule–Thomson cryocoolers working at subcritical pressures, we study a counter flow heat exchanger with condensation by employing the integral method. The effects of inlet pressure and working fluid are predicted. We also show that there is an optimal value of the enthalpy difference along the heat exchanger for which its length is minimal.

---

Enhanced boiling heat transfer using radial fins

A numerical bifurcation analysis is carried out in order to determine the solution structure of radial fins subjected to multi-boiling heat transfer mode. One-dimensional conduction is employed throughout the thermal analysis. The fluid heat transfer coefficient is temperature dependent on the three regimes of phase-change of the fluid. Six fin profiles, defined in the text, are considered. Multiplicity structure is obtained to determine different types of bifurcation diagrams, which describe the dependence of a state variable of the system like the temperature or the heat dissipation on the fin design parameters, conduction–convection parameter (CCP) or base temperature difference (ΔT). Specifically, the effects of ΔT, CCP and Biot number are analyzed. The results are presented graphically, showing the significant behavioral features of the heat rejection mechanism.

---

Variational optical flow estimation for particle image velocimetry

We introduce a novel class of algorithms for evaluating PIV image pairs. The mathematical basis is a continuous variational formulation for globally estimating the optical flow vector fields over the whole image. This class of approaches has been known in the field of image processing and computer vision for more than two decades but apparently has not been applied to PIV image pairs so far. We pay particular attention to a multi-scale representation of the image data so as to cope with the quite specific signal structure of particle image pairs. The experimental evaluation shows that a prototypical variational approach competes in noisy real-world scenarios with three alternative approaches especially designed for PIV-sequence evaluation. We outline the potential of the variational method for further developments.The publications of the CVGPR Group are listed under .

---

XPIV–Multi-plane stereoscopic particle image velocimetry

We introduce the three-dimensional measurement technique (XPIV) based on a Particle Image Velocimetry (PIV) system. The technique provides three-dimensional and statistically significant velocity data. The main principle of the technique lies in the combination of defocus, stereoscopic and multi-plane illumination concepts. Preliminary results of the turbulent boundary layer in a flume are presented. The quality of the velocity data is evaluated by using the velocity profiles and relative turbulent intensity of the boundary layer. The analysis indicates that the XPIV is a reliable experimental tool for three-dimensional fluid velocity measurements.More information at:

---

A physical decomposition of the stress tensor for complex flows

Traditionally, the components of the stress with respect to a relevant coordinate system are used for the purpose of stress visualisation and interpretation. A case for using a flow dependent measure to interpret and visualise stress is made for two dimensional flow, together with a suggestion for extending the idea to three dimensions. The method is illustrated for Newtonian and Oldroyd B fluids in both the eccentrically rotating cylinder and flow past a cylinder benchmark problems. In the context of a generalised Newtonian fluid, the relation between the flow-dependent stress measure to other field variables under certain flow conditions, is examined and is indicative of its importance in complex flow.

Attenuation technique for measuring sediment displacement levels

A technique for obtaining accurate, high (spatial) resolution measurements of sediment redeposition levels is described. In certain regimes, the method may also be employed to provide measurements of sediment layer thickness as a function of time. The method uses a uniform light source placed beneath the layer, consisting of transparent particles, so that the intensity of light at a point on the surface of the layer can be related to the depth of particles at that point. A set of experiments, using the impact of a vortex ring with a glass ballotini particle layer as the resuspension mechanism, are described to test and illustrate the technique.

---

Turbulence statistics from optical whole-field measurements in particle-laden turbulence

Obtaining turbulence statistics in particle-laden flows using optical whole-field measurements is complicated due to the inevitable data loss. The effects of this data loss are first studied using synthetic data and it is shown that the interpolation of missing data leads to biased results for the turbulence spectrum and its derived quantities. It is also shown that the use of overlapping interrogation regions in images with a low image density can lead to biased results due to oversampling. The slotting method is introduced for the processing of particle image velocimetry (PIV) data fields with missing data. Next to this, it is extended to handle unstructured data. Using experimental data obtained by a dual-camera PIV/PTV (particle tracking velocimetry) system in particle-laden grid turbulence, the performance of the new approach is studied. Some preliminary two-phase results are presented to indicate the significant improvement in the statistics, as well as to demonstrate the unique capabilities of the system.An erratum to this article can be found at

---

The influence of wall interaction on dynamic particle modelling of magneto-rheological suspensions between shearing plates

We report on the modelling of a magneto-rheological (MR) suspension bound between shearing parallel plates using a particle-level numerical simulation. The simulation is similar to an approach used previously but includes particle hydrodynamic interaction using elements of the Stokesian-dynamic method. Observations of initially chain-like aggregations are reported, and the evolving morphology of suspension particle clusters is explored. Our early-strain observations concur with the prevailing ideas of experimentalists on the important role that the microstructure has on bulk viscosity. We then study in particular the effects of simulation size and strain on viscosity. While initial viscous response is similar to previously reported observations in the literature, when left to run for longer strains, suspensions evolved into markedly different microstructures from those observed experimentally, or in electro-rheological suspensions, or MR simulations with artificial wall interaction. Substantial qualitative and quantitative divergence was observed over long strains. We argue that this divergence is due to the lack of a particle–wall interaction model for MR fluids. While current theories in MR modelling do not justify the requirement for a particle–wall interaction, these results suggest that one is required in order to match experimental observations.

Subsonic jet aeroacoustics: associating experiment,modelling and simulation

An overview of jet noise research is presented wherein the principal movements in the field are traced since its beginnings. Particular attention is paid to the evolution of our understanding of what we call a “source mechanism” in free shear flows; to the theoretical, experimental and numerical studies which have nurtured this understanding; and to the currently unresolved conceptual difficulties which render analysis of experimental and numerical data so difficult. As it is clear that accelerated progress in this field of research can be made possible by a more effective synergy between the theoretical, experimental and numerical disciplines—one which draws in particular on the impressive recent progress in experimental and numerical techniques—we endeavour to elucidate the various “source” characteristics identified by these different means of study; the points on which the studies agree or disagree, and the significance of such accord or discord; and, the new analysis possibilities which can now be realised by effectively associating experiment, modelling and simulation.

---