Polymer chains in turbulent flows are generally modelled as dumbbells, i.e. two beads joined by a nonlinear spring. The dumbbell only maps a single spatial configuration, described by the polymer end-to-end vector, thus a multi-bead FENE (finitely extensible nonlinear elastic) chain seems a natural improvement for a more accurate characterisation of the polymer spatial conformation. At a large Weissenberg number, a comparison with the more accurate Kuhn chain reveals that the multi-bead FENE chain drastically overestimates the probability of folded configurations. Surprisingly, the dumbbell turns out to be the only meaningful bead-spring model to coarse-grain a polymer macromolecule in turbulent pipe flows.

Two-way momentum-coupled direct numerical simulations of a particle-laden turbulent channel flow are addressed to investigate the effect of the particle Stokes number and of the particle-to-fluid density ratio on the turbulence modification. The exact regularised point-particle method is used to model the interphase momentum exchange in presence of solid boundaries, allowing the exploration of an extensive region of the parameter space. Results show that the particles increase the friction drag in the parameter space region considered, namely the Stokes number $St_+ \in [2,80]$, and the particle-to-fluid density ratio $\rho _p/\rho _f \in [90,5760]$ at a fixed mass loading $\phi =0.4$. It is noteworthy that the highest drag occurs for small Stokes number particles. A measurable drag increase occurs for all particle-to-fluid density ratios, the effect being reduced significantly only at the highest value of $\rho _p/\rho _f$. The modified stress budget and turbulent kinetic energy equation provide the rationale behind the observed behaviour. The particles’ extra stress causes an additional momentum flux towards the wall that modifies the structure of the buffer and of the viscous sublayer where the streamwise and wall-normal velocity fluctuations are increased. Indeed, in the viscous sublayer, additional turbulent kinetic energy is produced by the particles’ back-reaction, resulting in a strong augmentation of the spatial energy flux towards the wall where the energy is ultimately dissipated. This behaviour explains the increase of friction drag in particle-laden wall-bounded flows.

The exact regularised point particle (ERPP) method is extended to treat the inter-phase momentum coupling between particles and fluid in the presence of walls by accounting for vorticity generation due to particles close to solid boundaries. The ERPP method overcomes the limitations of other methods by allowing the simulation of an extensive parameter space (Stokes number, mass loading, particle-to-fluid density ratio and Reynolds number) and of particle spatial distributions that are uneven (few particles per computational cell). The enhanced ERPP method is explained in detail and validated by considering the global impulse balance. In conditions when particles are located close to the wall, a common scenario in wall-bounded turbulent flows, the main contribution to the total impulse arises from the particle-induced vorticity at the solid boundary. The method is applied to direct numerical simulations of particle-laden turbulent pipe flow in the two-way coupling regime to address turbulence modulation. The effects of the mass loading, the Stokes number and the particle-to-fluid density ratio are investigated. The drag is either unaltered or increased by the particles with respect to the uncoupled case. No drag reduction is found in the parameter space considered. The momentum stress budget, which includes an extra stress contribution by the particles, provides the rationale behind the drag behaviour. The extra stress produces a momentum flux towards the wall that strongly modifies the viscous stress, the culprit of drag at solid boundaries.

Selected samples of asbestos-containing material (ACM) with different Ca/Si ratios have been treated thermally at 1200°C for 15 min to obtain an ‘end of waste geo-inspired material’. Before and after treatment, micro-Raman spectroscopy allowed the investigation of both powdered and massive samples by directing the laser beam onto crystals with elongated morphology, thin fibres and the matrix. In the raw samples, chrysotile and/or crocidolite were detected. After the thermal treatment, no asbestos phases were identified in the Raman spectra collected on fibrous or fibre-like morphologies. The scanning electron microscopy/energy dispersive spectroscopy investigations confirmed the onset of a pseudomorphic process during annealing, leading to the complete transformation of asbestos minerals into non-hazardous magnesium or calcium magnesium silicates such as forsterite, monticellite, a˚kermanite and merwinite.

The identification of such mineral assemblages was inspired by the close inspection of a natural counterpart, the high-temperature contact metamorphic imprint due to the intrusion of a sill into carbonate rocks. The process turned out to occur largely at the solid state and involved substantial mobilization of Ca and Mg to form a spinel phase (namely MgFe2O4) which was recognized in the matrix and within, or close to elongated morphologies.

The generalised Kolmogorov equation is used to describe the scale-by-scale turbulence dynamics in the shear layer and in the separation bubble generated by a bulge at one of the walls in a turbulent channel flow. The second-order structure function, which is the basis of such an equation, is used as a proxy to define a scale-energy content, that is an interpretation of the energy associated with a given scale. Production and dissipation regions and the flux interchange between them, in both physical and separation space, are identified. Results show how the generalised Kolmogorov equation, a five-dimensional equation in our anisotropic and strongly inhomogeneous flow, can describe the turbulent flow behaviour and related energy mechanisms. Such complex statistical observables are linked to a visual inspection of instantaneous turbulent structures detected by means of the Q-criterion. Part of these turbulent structures are trapped in the recirculation where they undergo a pseudo-cyclic process of disruption and reformation. The rest are convected downstream, grow and tend to larger streamwise scales in an inverse cascade. The classical picture of homogeneous isotropic turbulence in which energy is fed at large scales and transferred to dissipate at small scales does not simply apply to this flow where the energy dynamics strongly depends on position, orientation and length scale.

Turbulent flow separation induced by a protuberance on one of the walls of an otherwise planar channel is investigated using direct numerical simulations. Different bulge geometries and Reynolds numbers – with the highest friction Reynolds number simulation reaching a peak of $Re_{\unicode[STIX]{x1D70F}}=900$ – are addressed to understand the effect of the wall curvature and of the Reynolds number on the dynamics of the recirculating bubble behind the bump. Global quantities reveal that most of the drag is due to the form contribution, whilst the friction contribution does not change appreciably with respect to an equivalent planar channel flow. The size and position of the separation bubble strongly depends on the bump shape and the Reynolds number. The most bluff geometry has a larger recirculation region, whilst the Reynolds number increase results in a smaller recirculation bubble and a shear layer more attached to the bump. The position of the reattachment point only depends on the Reynolds number, in agreement with experimental data available in the literature. Both the mean and the turbulent kinetic energy equations are addressed in such non-homogeneous conditions revealing a non-trivial behaviour of the energy fluxes. The energy introduced by the pressure drop follows two routes: part of it is transferred towards the walls to be dissipated and part feeds the turbulent production hence the velocity fluctuations in the separating shear layer. Spatial energy fluxes transfer the kinetic energy into the recirculation bubble and downstream near the wall where it is ultimately dissipated. Consistently, anisotropy concentrates at small scales near the walls irrespective of the value of the Reynolds number. In the bulk flow and in the recirculation bubble, isotropy is restored at small scales and the isotropy recovery rate is controlled by the Reynolds number. Anisotropy invariant maps are presented, showing the difficulty in developing suitable turbulence models to predict separated turbulent flow dynamics. Results shed light on the processes of production, transfer and dissipation of energy in this relatively complex turbulent flow where non-homogeneous effects overwhelm the classical picture of wall-bounded turbulent flows which typically exploits streamwise homogeneity.

The hydrodynamics of a flagellated micro-organism is investigated when swimming close to a planar free-slip surface by means of numerical solutions of the Stokes equations obtained via a boundary element method. Depending on the initial conditions, the swimmer can either escape from the free-slip surface or collide with the boundary. Interestingly, the micro-organism does not exhibit a stable orbit. Independently of escape or attraction to the interface, close to a free-slip surface, the swimmer follows a counter-clockwise trajectory, in agreement with experimental findings (Di Leonardo et al., Phys. Rev. Lett., vol. 106 (3), 2011, 038101). The hydrodynamics is indeed modified by the free surface. In fact, when the same swimmer moves close to a no-slip wall, a set of initial conditions exists which result in stable orbits. Moreover, when moving close to a free-slip or a no-slip boundary, the swimmer assumes a different orientation with respect to its trajectory. Taken together, these results contribute to shed light on the hydrodynamical behaviour of micro-organisms close to liquid–air interfaces which are relevant for the formation of interfacial biofilms of aerobic bacteria.

Particulate flows have mainly been studied under the simplifying assumption of a one-way coupling regime where the disperse phase does not modify the carrier fluid. A more complete view of multiphase flows can be gained calling into play two-way coupling effects, i.e. by accounting for the inter-phase momentum exchange, which is certainly relevant at increasing mass loading. In this paper we present a new methodology rigorously designed to capture the inter-phase momentum exchange for particles smaller than the smallest hydrodynamical scale, e.g. the Kolmogorov scale in a turbulent flow. The momentum coupling mechanism exploits the unsteady Stokes flow around a small rigid sphere, where the transient disturbance produced by each particle is evaluated in a closed form. The particles are described as lumped point masses, which would lead to the appearance of singularities. A rigorous regularization procedure is conceived to extract the physically relevant interactions between the particles and the fluid which avoids any ‘ad hoc’ assumption. The approach is suited for high-efficiency implementation on massively parallel machines since the transient disturbance produced by the particles is strongly localized in space. We will show that hundreds of thousands of particles can be handled at an affordable computational cost, as demonstrated by a preliminary application to a particle-laden turbulent shear flow.

Turbulent fluctuations induce the common phenomenon known as clustering in the spatial arrangement of small inertial particles transported by the fluid. Particles spread non-uniformly, and form clusters where their local concentration is much higher than in nearby rarefaction regions. The underlying physics has been exhaustively analysed in the so-called one-way coupling regime, i.e. negligible back-reaction of the particles on the fluid, where the mean flow anisotropy induces preferential orientation of the clusters. Turbulent transport in suspensions with significant mass in the disperse phase, i.e. particles back-reacting in the carrier phase (the two-way coupling regime), has instead been much less investigated and is still poorly understood. The issue is discussed here by addressing direct numerical simulations of particle-laden homogeneous shear flows in the two-way coupling regime. Consistent with previous findings, we observe an overall depletion of the turbulent fluctuations for particles with response time of the order of the Kolmogorov time scale. The depletion occurs in the energy-containing range, while augmentation is observed in the small-scale range down to the dissipative scales. Increasing the mass load results in substantial broadening of the energy cospectrum, thereby extending the range of scales driven by anisotropic production mechanisms. As discussed throughout the paper, this is due to the clusters which form the spatial support of the back-reaction field and give rise to a highly anisotropic forcing, active down to the smallest scales. A certain impact on two-phase flow turbulence modelling is expected from the above conclusions, since the frequently assumed small-scale isotropy is poorly recovered when the coupling between the phases becomes significant.

Recently, clustering of inertial particles in turbulence has been thoroughly analysed for statistically homogeneous isotropic flows. Phenomenologically, spatial homogeneity of particle configurations is broken by the advection of a range of eddies determined by the Stokes relaxation time of the particles. This in turn results in a multi-scale distribution of local particle concentration and voids. Much less is known concerning anisotropic flows. Here, by addressing direct numerical simulations (DNS) of a statistically steady particle-laden homogeneous shear flow, we provide evidence that the mean shear preferentially orients particle patterns. By imprinting anisotropy on large-scale velocity fluctuations, the shear indirectly affects the geometry of the clusters. Quantitative evaluation is provided by a purposely designed tool, the angular distribution function (ADF) of particle pairs, which allows to address the anisotropy content of particle aggregates on a scale-by-scale basis. The data provide evidence that, depending on the Stokes relaxation time of the particles, anisotropic clustering may occur even in the range of scales in which the carrier phase velocity field is already recovering isotropy. The strength of the singularity in the anisotropic component of the ADF quantifies the level of fine-scale anisotropy, which may even reach values of more than 30% direction-dependent variation in the probability to find two closeby particles at viscous-scale separation.

We discuss how large-eddy simulation (LES) can be properly employed to predict the statistics of the resolved velocity fluctuations in shear turbulence. To this purpose an a posteriori comparison of LES data against filtered direct numerical simulation (DNS) is used to establish the necessary conditions that the filter scale LF must satisfy to achieve the preservation of the statistical properties of the resolved field. In this context, by exploiting the physical role of the shear scale LS, the Kármán–Howarth equation allows for the assessment of LES data in terms of scale-by-scale energy production, energy transfer and subgrid energy fluxes. Even higher-order statistical properties of the resolved scales such as the probability density function of longitudinal velocity increments are well reproduced, provided the relative position of the filter scale with respect to the shear scale is properly selected. We consider here the homogeneous shear flow as the simplest non-trivial flow which fully retains the basic mechanism of turbulent kinetic energy production typical of any shear flow, with the advantage that spatial homogeneity implies a well-defined value of the shear scale while numerical difficulties related to resolution requirements in the near wall region are avoided.

Turbulent shear flows, such as those occurring in the wall region of turbulent boundary layers, show a substantial increase of intermittency in comparison with isotropic conditions. This suggests a close link between anisotropy and intermittency. However, a rigorous statistical description of anisotropic flows is, in most cases, hampered by the inhomogeneity of the field. This difficulty is absent for homogeneous shear flow. For this flow the scale-by-scale budget is discussed here by using the appropriate form of the Kámán–Howarth equation, to determine the range of scales where the shear is dominant. The resulting generalization of the four-fifths law is then used to extend to shear-dominated flows the Kolmogorov–Oboukhov theory of intermittency. The procedure leads naturally to the formulation of generalized structure functions, and the description of intermittency thus obtained reduces to the K62 theory for vanishing shear. The intermittency corrections to the scaling exponents are related to the moments of the coarse-grained energy dissipation field. Numerical experiments give indications that the dissipation field is statistically unaffected by the shear, supporting the conjecture that the intermittency corrections are universal. This observation together with the present reformulation of the theory gives a reason for the increased intermittency observed in the classical longitudinal velocity increments.

A case of double bilateral renal vessels in both twins of a MZ pair is reported for the first time. Even in non-twins, the anomalies reported so far involved the inferior polar arteries in agreement with the embryological development. The two male twins, examined at the age of 14 years, had simultaneously developed a marked hypertension at the age of 7 years. Zygosity was determined by blood group and HLA analysis and various clinical tests were carried out to diagnose the condition. It is suggested that the anomaly is the result of a genetically induced early block in embryonic development.

Raman scattering has been used to examine strained-layer GaSb/AlSb superlattices for a variety of growth imperfections. Analysis of zone-folded acoustic phonon spectra permit a determination of the superlattice period, the individual layer widths, some estimate of the interfacial widths, and evidence for variations in layer thicknesses.