We present a study of Eulerian and Lagrangian statistics from a high-resolution numerical simulation of isotropic and homogeneous turbulence using the FLASH code, with an estimated Taylor microscale ...Reynolds number of around 600. Statistics are evaluated over a data set with 18563 spatial grid points and with 2563 = 16.8 million particles, followed for about one large-scale eddy turnover time. We present data for the Eulerian and Lagrangian structure functions up to the tenth order. We analyze the local scaling properties in the inertial range. The Eulerian velocity field results show good agreement with previous data and confirm the puzzling differences previously found between the scaling of the transverse and the longitudinal structure functions. On the other hand, accurate measurements of sixth-and-higher-order Lagrangian structure functions allow us to highlight some discrepancies from earlier experimental and numerical results. We interpret this result in terms of a possible contamination from the viscous scale, which may have affected estimates of the scaling properties in previous studies. We show that a simple bridge relation based on a multifractal theory is able to connect scaling properties of both Eulerian and Lagrangian observables, provided that the small differences between intermittency of transverse and longitudinal Eulerian structure functions are properly considered.
Human blood flow is a multiscale problem: in first approximation, blood is a dense suspension of plasma and deformable red cells. Physiological vessel diameters range from about one to thousands of ...cell radii. Current computational models either involve a homogeneous fluid and cannot track particulate effects or describe a relatively small number of cells with high resolution but are incapable to reach relevant time and length scales. Our approach is to simplify much further than existing particulate models. We combine well-established methods from other areas of physics in order to find the essential ingredients for a minimalist description that still recovers hemorheology. These ingredients are a lattice Boltzmann method describing rigid particle suspensions to account for hydrodynamic long-range interactions and-in order to describe the more complex short-range behavior of cells-anisotropic model potentials known from molecular-dynamics simulations. Paying detailedness, we achieve an efficient and scalable implementation which is crucial for our ultimate goal: establishing a link between the collective behavior of millions of cells and the macroscopic properties of blood in realistic flow situations. In this paper we present our model and demonstrate its applicability to conditions typical for the microvasculature.
The relative dispersion of pairs of inertial point particles in incompressible, homogeneous and isotropic three-dimensional turbulence is studied by means of direct numerical simulations at two ...values of the Taylor-scale Reynolds number Reλ ~ 200 and Reλ ~ 400, corresponding to resolutions of 5123 and 20483 grid points, respectively. The evolution of both heavy and light particle pairs is analysed by varying the particle Stokes number and the fluid-to-particle density ratio. For particles much heavier than the fluid, the range of available Stokes numbers is St ∈ 0.1 : 70, while for light particles the Stokes numbers span the range St ∈ 0.1 : 3 and the density ratio is varied up to the limit of vanishing particle density. For heavy particles, it is found that turbulent dispersion is schematically governed by two temporal regimes. The first is dominated by the presence, at large Stokes numbers, of small-scale caustics in the particle velocity statistics, and it lasts until heavy particle velocities have relaxed towards the underlying flow velocities. At such large scales, a second regime starts where heavy particles separate as tracers' particles would do. As a consequence, at increasing inertia, a larger transient stage is observed, and the Richardson diffusion of simple tracers is recovered only at large times and large scales. These features also arise from a statistical closure of the equation of motion for heavy particle separation that is proposed and is supported by the numerical results. In the case of light particles with high density ratio, strong small-scale clustering leads to a considerable fraction of pairs that do not separate at all, although the mean separation increases with time. This effect strongly alters the shape of the probability density function of light particle separations.
In this work, we show how to control the morphology of femtosecond laser melted gold nanosphere aggregates. A careful choice of both laser fluence and wavelength makes it possible to selectively ...excite different aggregate substructures to produce larger spherical nanoparticles, nanorods, and nanoprisms or necklace-like 1D nanostructures in which the nanoparticles are interlinked by bridges. Finite integral technique calculations have been performed on the near-field concentration of light in the nanostructures which confirm the wavelength dependence of the light concentration and suggest that the resulting localized high intensities lead to nonthermal melting. We show that by tuning the wavelength of the melting light it is possible to choose the spatial extension of the ensembles of NPs heated thus allowing us to exhibit control over the morphology of the nanostructures formed by the melting process. By a proper combination of this method with self-assembly of chemically synthesized nanoparticles, one can envisage the development of an innovative high-throughput high-resolution nanofabrication technique.
The dynamics of particles in turbulence when the particle size is larger than the dissipative scale of the carrier flow are studied. Recent experiments have highlighted signatures of particles' ...finiteness on their statistical properties, namely a decrease of their acceleration variance, an increase of correlation times (at increasing the particles size) and an independence of the probability density function of the acceleration once normalized to their variance. These effects are not captured by point-particle models. By means of a detailed comparison between numerical simulations and experimental data, we show that a more accurate description is obtained once Faxén corrections are included.
Planktonic copepods are small crustaceans that have the ability to swim by quick powerful jumps. Such an aptness is used to escape from high shear regions, which may be caused either by flow ...perturbations, produced by a large predator (i.e., fish larvae), or by the inherent highly turbulent dynamics of the ocean. Through a combined experimental and numerical study, we investigate the impact of jumping behavior on the small-scale patchiness of copepods in a turbulent environment. Recorded velocity tracks of copepods displaying escape response jumps in still water are here used to define and tune a Lagrangian copepod (LC) model. The model is further employed to simulate the behavior of thousands of copepods in a fully developed hydrodynamic turbulent flow obtained by direct numerical simulation of the Navier-Stokes equations. First, we show that the LC velocity statistics is in qualitative agreement with available experimental observations of copepods in turbulence. Second, we quantify the clustering of LC, via the fractal dimension D_{2}. We show that D_{2} can be as low as ∼2.3 and that it critically depends on the shear-rate sensitivity of the proposed LC model, in particular it exhibits a minimum in a narrow range of shear-rate values. We further investigate the effect of jump intensity, jump orientation, and geometrical aspect ratio of the copepods on the small-scale spatial distribution. At last, possible ecological implications of the observed clustering on encounter rates and mating success are discussed.
.
Particle-laden turbulent flows occur in a variety of industrial applications as well as in naturally occurring flows. While the numerical simulation of such flows has seen significant advances in ...recent years, it still remains a challenging problem. Many studies investigated the rheology of dense suspensions in laminar flows as well as the dynamics of point-particles in turbulence. Here we employ a fully-resolved numerical simulation based on a lattice Boltzmann scheme, to investigate turbulent flow with large neutrally buoyant particles in a pipe flow at low Reynolds number and in dilute regimes. The energy input is kept fixed resulting in a Reynolds number based on the friction velocity around 250. Two different particle radii were used giving a particle-pipe diameter ratio of 0.05 and 0.075. The number of particles is kept constant resulting in a volume fraction of 0.54% and 1.83%, respectively. We investigated Eulerian and Lagrangian statistics along with the stresslet exerted by the fluid on the spherical particles. It was observed that the high particle-to-fluid slip velocity close to the wall corresponds locally to events of high energy dissipation, which are not present in the single-phase flow. The migration of particles from the inner to the outer region of the pipe, the dependence of the stresslet on the particle radial positions and a proxy for the fragmentation rate of the particles computed using the stresslet have been investigated.
Graphical abstract
Hybrid particle-continuum computational frameworks permit the simulation of gas flows by locally adjusting the resolution to the degree of non-equilibrium displayed by the flow in different regions ...of space and time. In this work, we present a new scheme that couples the direct simulation Monte Carlo (DSMC) with the lattice Boltzmann (LB) method in the limit of isothermal flows. The former handles strong non-equilibrium effects, as they typically occur in the vicinity of solid boundaries, whereas the latter is in charge of the bulk flow, where non-equilibrium can be dealt with perturbatively, i.e. according to Navier-Stokes hydrodynamics. The proposed concurrent multiscale method is applied to the dilute gas Couette flow, showing major computational gains when compared with the full DSMC scenarios. In addition, it is shown that the coupling with LB in the bulk flow can speed up the DSMC treatment of the Knudsen layer with respect to the full DSMC case. In other words, LB acts as a DSMC accelerator.
This article is part of the themed issue ‘Multiscale modelling at the physics–chemistry–biology interface’.
A key ecological parameter for planktonic copepod studies is their encounter rates within the same population as well as with other species. The encounter rate is partly determined by copepod’s ...swimming behaviour and is strongly influenced by turbulence of the surrounding environment. A distinctive feature of copepods’ motility is their ability to perform quick displacements, often dubbed jumps, by means of powerful swimming strokes. Such a reaction has been associated to an escape behaviour from flow disturbances due to predators or other external signals. In the present study, we investigate the encounter rates of copepods from the same species in a developed turbulent flow with intensities comparable to those encountered in their natural habitat. This is done by means of a Lagrangian copepod (LC) model that mimics the jump escape reaction behaviour from localisedhigh-deformation rate fluctuations in the turbulent flows. Our analysis shows that the encounter rate for copepods of typical perception radius of ∼ η, where η is the dissipative scale of turbulence, can be increased by a factor up to ∼ 102 compared to the one experienced by passively transported fluid tracers of the same size. Furthermore, we address the effect of a minimal waiting time between consecutive jumps. It is shown that any encounter-rate enhancement is lost if such time goes beyond the dissipative time-scale of turbulence, τη. Because typically in the ocean η ∼ 1 mm and τη ∼ 1 s,this provides stringent constraints on the turbulent-driven enhancement of contact-rate due to apurely mechanical induced escape reaction. The implications of our results in the context of ecologyof copepods were discussed.