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.

---

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 .

---

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.

---

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.

---

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.

---

Comment on the Clauser chart method for determining the friction velocity

A known difficulty with using the Clauser chart method to determine the friction velocity in wall bounded flows is that it assumes, a priori, a logarithmic law for the mean velocity profile. Using both experimental and DNS data in the literature, this note explicitly shows how friction velocities obtained using the Clauser chart method can potentially mask subtle Reynolds-number-dependent behavior.

---

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.

---

An alternative approach to the heat equation

The technique of conformal transformations of the metric, widely used in general relativity and in cosmology, is applied to the analysis of heat conduction in an anisotropic medium, in which the thermal conductivity is described by a tensor instead of a scalar. The anisotropic medium is reduced to an effective equivalent one, which is isotropic. The simplification is achieved for a particular conformal factor of the transformation, uniquely determined under physically reasonable assumptions on the thermal conductivity, density, and specific heat of the medium. Another application consists in the formal elimination of source or sink terms from the heat equation by using a suitable conformal transformation.

---

A new interpretation for the dynamic behaviour of complex fluids at the sol–gel transition using the fractional calculus

We propose to analyse power law shear stress relaxation modulus observed at the sol–gel transition (SGT) in many gelling systems in terms of fractional calculus. We show that the critical gel (gel at SGT) can be associated to a single fractional element and the gel in the post-SGT state to a fractional Kelvin–Voigt model. In this case, it is possible to give a physical interpretation to the fractional derivative order. It is associated to the power law exponent of the shear modulus related to the fractal dimension of the critical gel. A preliminary experimental application to silica alkoxide-based systems is given.

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.

---

Analysis of velocity interpolation schemes for image deformation methods in PIV

The spatial resolution of PIV can be increased significantly by using an image deformation method (IDM) and very small grid distance (i.e. the final distance between vectors), therefore, also increasing the processing time. By using an interpolation scheme with a good spectral response, in the dense predictor step of the algorithm, it is possible to increase the grid distance without decreasing the spatial resolution therefore decreasing the total processing time.

---

Linear rheology of comb polymers with star-like backbones: melts and solutions

We investigate the linear rheology of model star-comb homopolymers consisting of star-like backbone chains with grafted branches. We show that the tube-based theory in the framework of full dynamic dilution, appropriately modified to account for the effects of the fluctuations of the free segments (the segments of the star arms between the outer branching points and the arms’ free end-monomers) of the backbone star arms, and the polydispersity accurately describes the linear viscoelastic spectrum. For these branched polymers, the relaxation is found to proceed hierarchically, similar to combs with linear backbones. However, in contrast to the latter, here, there is no reptation. The higher functionality star combs with sparse branching are particularly interesting because they resemble a Cayley tree structure and are treated in the context of asymmetric star polymers. We also employ solutions of star combs and test our model in static dilution conditions. In the region where despite the presence of both static (solvent) and dynamic dilution the polymers remain entangled, our model is particularly successful without adjusting any parameters.

Time-resolved reconstruction of the full velocity field around a dynamically-scaled flapping wing

The understanding of the physics of flapping flight has long been limited due to the obvious experimental difficulties in studying the flow field around real insects. In this study the time-dependent three-dimensional velocity field around a flapping wing was measured quantitatively for the first time. This was done using a dynamically-scaled wing moving in mineral oil in a pattern based on the kinematics obtained from real insects. The periodic flow is very reproducible, due to the relatively low Reynolds number and precise control of the wing. This repeatability was used to reconstruct the full evolving flow field around the wing from separate stereoscopic particle image velocimetry measurements for a number of spanwise planes and time steps. Typical results for two cases (an impulsive start and a simplified flapping pattern) are reported. Visualizations of the obtained data confirm the general picture of the leading-edge vortex that has been reported in recent publications, but allow a refinement of the detailed structure: rather than a single strand of vorticity, we find a stable pair of counter-rotating structures. We show that the data can also be used for quantitative studies, such as lift and drag prediction.

---

Recent experimental probes of shear banding

Recent experimental techniques used to investigate shear banding are reviewed. After recalling the rheological signature of shear-banded flows, we summarize the various tools for measuring locally the microstructure and the velocity field under shear. Local velocity measurements using dynamic light scattering and ultrasound are emphasized. A few results are extracted from current works to illustrate open questions and directions for future research.

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.

---

Density stratified environments: the double-tank method

We present an alternative method of producing density stratifications in the laboratory based on the ‘double-tank’ method proposed by Oster (Sci Am 213:70–76, 1965). We refer to Oster’s method as the ‘forced-drain’ approach, as the volume flow rates between connecting tanks are controlled by mechanical pumps. We first determine the range of density profiles that may be established with the forced-drain approach other than the linear stratification predicted by Oster. The dimensionless density stratification is expressed analytically as a function of three ratios: the volume flow rate ratio n, the ratio of the initial liquid volumes λ and the ratio of the initial densities ψ. We then propose a method which does not require pumps to control the volume flow rates but instead allows the connecting tanks to drain freely under gravity. This is referred to as the ‘free-drain’ approach. We derive an expression for the density stratification produced and compare our predictions with saline stratifications established in the laboratory using the ‘free-drain’ extension of Oster’s method. To assist in the practical application of our results we plot the region of parameter space that yield concave/convex or linear density profiles for both forced-drain and free-drain approaches. The free-drain approach allows the experimentalist to produce a broad range of density profiles by varying the initial liquid depths, cross-sectional and drain opening areas of the tanks. One advantage over the original Oster approach is that density profiles with an inflexion point can now be established.

---

Influence of longitudinal conduction in the matrix on effectiveness of rotary heat regenerator used in air-conditioning

The paper presents an efficient analytical method of solving the problem of a rotary heat regenerator taking into account longitudinal heat conduction in the matrix. The small parameter method, Laplace transform as well as one of the spline functions have been applied for approximation of an initial condition in the reversion time. In the application part, the solution for a model in analysis of an influence of longitudinal conduction in the matrix on effectiveness of rotary heat regenerator in a wide range of dimensionless parameters as well as for the particular matrix applied in air-conditioning was used.

---

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.