Source author record

Federico Toschi

Federico Toschi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

57works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

57 published item(s)

preprint2023arXiv

Stochastic fluctuations of diluted pedestrian dynamics along curved paths

As we walk towards our destinations, our trajectories are constantly influenced by the presence of obstacles and infrastructural elements: even in absence of crowding our paths are often curved. Over the last two decades pedestrian dynamics have been extensively studied aiming at quantitative models with both fundamental and technological relevance. Walking kinematics along straight paths have been experimentally investigated and quantitatively modeled in the diluted limit (i.e. in absence of pedestrian-pedestrian interactions). It is natural to expect that models for straight paths may be an accurate approximations of the dynamics even for paths with curvature radii much larger than the size of a single person. Conversely, as paths curvature increase one may expect larger and larger deviations. As no clear experimental consensus has been reached yet in the literature, here we accurately and systematically investigate the effect of paths curvature on diluted pedestrian dynamics. Thanks to a extensive and highly accurate set of real-life measurements campaign, we derive and validate via a Langevin-like social-force model capable of quantitatively describing both averages and fluctuations. Leveraging on the differential geometric notion of covariant derivative, we generalize previous work by some of the authors, effectively casting a Langevin social-force model for the straight walking dynamics in a curved geometric setting. We deem this the necessary first step to understand and model the more general and ubiquitous case of pedestrians following curved paths in the presence of crowd traffic.

preprint2022arXiv

Law of the wall in an unstably stratified turbulent channel flow

We perform direct numerical simulations of an unstably stratified turbulent channel flow to address the effects of buoyancy on the boundary layer dynamics and mean field quantities. We systematically span a range of parameters in the space of friction Reynolds number ($Re_τ$) and Rayleigh number ($Ra$). Our focus is on deviations from the logarithmic law of the wall due to buoyant motion. The effects of convection in the relevant ranges are discussed providing measurements of mean profiles of velocity, temperature and Reynolds stresses as well as of the friction coefficient. A phenomenological model is proposed and shown to capture the observed deviations of the velocity profile in the log-law region from the non-convective case.

preprint2022arXiv

Spatial population genetics with fluid flow

The growth and evolution of microbial populations is often subjected to advection by fluid flows in spatially extended environments, with immediate consequences for questions of spatial population genetics in marine ecology, planktonic diversity and origin of life scenarios. Here, we review recent progress made in understanding this rich problem in the simplified setting of two competing genetic microbial strains subjected to fluid flows. As a pedagogical example we focus on antagonsim, i.e., two killer microorganism strains, each secreting toxins that impede the growth of their competitors (competitive exclusion), in the presence of stationary fluid flows. By solving two coupled reaction-diffusion equations that include advection by simple steady cellular flows composed of characteristic flow motifs in two dimensions (2d), we show how local flow shear and compressibility effects can interact with selective advantage to have a dramatic influence on genetic competition and fixation in spatially distributed populations. We analyze several 1d and 2d flow geometries including sources, sinks, vortices and saddles, and show how simple analytical models of the dynamics of the genetic interface can be used to shed light on the nucleation, coexistence and flow-driven instabilities of genetic drops. By exploiting an analogy with phase separation with nonconserved order parameters, we uncover how these genetic drops harness fluid flows for novel evolutionary strategies, even in the presence of number fluctuations, as confirmed by agent-based simulations as well.

preprint2021arXiv

A quantitative comparison of physical accuracy and numerical stability of Lattice Boltzmann color gradient and pseudopotential multicomponent models for microfluidic applications

The performances of the Color-Gradient (CG) and of the Shan-Chen (SC) multicomponent Lattice Boltzmann models are quantitatively compared side-by-side on multiple physical flow problems where breakup, coalescence and contraction of fluid ligaments are important. The flow problems are relevant to microfluidic applications, jetting of microdroplets as seen in inkjet printing, as well as emulsion dynamics. A significantly wider range of parameters is shown to be accessible for CG in terms of density-ratio, viscosity-ratio and surface tension values. Numerical stability for a high density ratio $\mathcal{O}(1000)$ is required for simulating the drop formation process during inkjet printing which we show here to be achievable using the CG model but not using the SC model. In terms of physical accuracy, the CG model shows good agreement with analytical solutions for droplet oscillation and ligament contraction test-cases. The SC model is effective in accurately simulating ligament contraction, but less accurate in handling droplet oscillations. Rayleigh-Plateau instability simulations show that both the CG and SC models give physically realistic results when breakup occurs, although smaller satellite droplets quickly vanish in the SC model due to droplet mass evaporating into the ambient phase. Our results show that the CG model is a suitable choice for challenging simulations of droplet formation, due to a combination of both numerical stability and physical accuracy. We also present a novel approach to incorporate repulsion forces between interfaces for CG, with possible applications to the study of stabilized emulsions. Specifically, we show that the CG model can produce similar results to a known multirange potentials extension of the SC model for modelling a disjoining pressure, opening up its use for the study of dense stabilized emulsions.

preprint2021arXiv

Benchmarking high-fidelity pedestrian tracking systems for research, real-time monitoring and crowd control

High-fidelity pedestrian tracking in real-life conditions has been an important tool in fundamental crowd dynamics research allowing to quantify statistics of relevant observables including walking velocities, mutual distances and body orientations. As this technology advances, it is becoming increasingly useful also in society. In fact, continued urbanization is overwhelming existing pedestrian infrastructures such as transportation hubs and stations, generating an urgent need for real-time highly-accurate usage data, aiming both at flow monitoring and dynamics understanding. To successfully employ pedestrian tracking techniques in research and technology, it is crucial to validate and benchmark them for accuracy. This is not only necessary to guarantee data quality, but also to identify systematic errors. In this contribution, we present and discuss a benchmark suite, towards an open standard in the community, for privacy-respectful pedestrian tracking techniques. The suite is technology-independent and is applicable to academic and commercial pedestrian tracking systems, operating both in lab environments and real-life conditions. The benchmark suite consists of 5 tests addressing specific aspects of pedestrian tracking quality, including accurate crowd flux estimation, density estimation, position detection and trajectory accuracy. The output of the tests are quality factors expressed as single numbers. We provide the benchmark results for two tracking systems, both operating in real-life, one commercial, and the other based on overhead depth-maps developed at TU Eindhoven. We discuss the results on the basis of the quality factors and report on the typical sensor and algorithmic performance. This enables us to highlight the current state-of-the-art, its limitations and provide installation recommendations, with specific attention to multi-sensor setups and data stitching.

preprint2021arXiv

Build up of yield stress fluids via chaotic emulsification

Stabilized dense emulsions display a rich phenomenology connecting microstructure and rheology. In this work we study how an emulsion with a finite yield stress can be built via large-scale stirring. By gradually increasing the volume fraction of the dispersed minority phase, under the constant action of a stirring force, we are able to achieve volume fractions close to $80\%$. Despite the fact that our system is highly concentrated and not yet turbulent we observe a droplet size distribution consistent with the $-10/3$ scaling, often associated to inertial range droplets breakup. We report that the polydispersity of droplet sizes correlates with the dynamics of the emulsion formation process. Additionally we quantify the visco-elastic properties of the dense emulsion finally obtained and we demonstrate the presence of a finite yield stress. The approach reported can pave the way to a quantitative understanding of the complex interplay between the dynamics of mesoscale constituents and the large scale flow properties yield-stress fluids.

preprint2021arXiv

On the drag and lift coefficients of ellipsoidal particles under rarefied flow conditions

The capability to simulate a two-way coupled interaction between a rarefied gas and an arbitrary-shaped colloidal particle is important for many practical applications, such as aerospace engineering, lung drug deliver and semiconductor manufacturing. By means of numerical simulations based on the Direct Simulation Monte Carlo (DSMC) method, we investigate the influence of the orientation of the particle and rarefaction on the drag and lift coefficients, in the case of prolate and oblate ellipsoidal particles immersed in a uniform ambient flow. This is done by modelling the solid particles using a cut-cell algorithm embedded within our DSMC solver. In this approach, the surface of the particle is described by its analytical expression and the microscopic gas-solid interactions are computed exactly using a ray-tracing technique. The measured drag and lift coefficients are used to extend the correlations available in the continuum regime to the rarefied regime, focusing on the transitional and free-molecular regimes. The functional forms for the correlations for the ellipsoidal particles are chosen as a generalisation from the spherical case. We show that the fits over the data from numerical simulations can be extended to regimes outside the simulated range of $Kn$ by testing the obtained predictive model on values of $Kn$ that where not included in the fitting process, allowing to achieve an higher precision when compared with existing predictive models from literature. Finally, we underline the importance of this work in providing new correlations for non-spherical particles that can be used for point-particle Euler-Lagrangian simulations to address the problem of contamination from finite-size particles in high-tech mechanical systems.

preprint2021arXiv

Self Similar Properties of Avalanche Statistics in a Simple Turbulent Model

In this paper, we consider a simplified model of turbulence for large Reynolds numbers driven by a constant power energy input on large scales. In the statistical stationary regime, the behaviour of the kinetic energy is characterised by two well defined phases: a laminar phase where the kinetic energy grows linearly for a (random) time $t_w$ followed by abrupt avalanche-like energy drops of sizes $S$ due to strong intermittent fluctuations of energy dissipation. We study the probability distribution $P[t_w]$ and $P[S]$ which both exhibit a quite well defined scaling behaviour. Although $t_w$ and $S$ are not statistically correlated, we suggest and numerically checked that their scaling properties are related based on a simple, but non trivial, scaling argument. We propose that the same approach can be used for other systems showing avalanche-like behaviour such as amorphous solids and seismic events.

preprint2020arXiv

Global and local statistics in turbulent emulsions

Turbulent emulsions are complex physical systems characterized by a strong and dynamical coupling between small-scale droplets and large-scale rheology. By using a specifically designed Taylor-Couette (TC) shear flow system, we are able to characterize the statistical properties of a turbulent emulsion made of oil droplets dispersed in an ethanol-water continuous solution, at the oil volume fraction up to 40%. We find that the dependence of the droplet size on the Reynolds number of the flow at the volume fraction of 1% can be well described by Hinze's criterion. The distribution of droplet sizes is found to follow a log-normal distribution, hinting at a fragmentation process as the possible mechanism dominating droplet formation. Additionally, the effective viscosity of the turbulent emulsion increases with the volume fraction of the dispersed oil phase, and decreases when the shear strength is increased. We find that the dependence of the effective viscosity on the shear rate can be described by the Herschel-Bulkley model, with a flow index monotonically decreasing with increasing the oil volume fraction. This finding indicates that the degree of shear thinning systematically increases with the volume fraction of the dispersed phase. The current findings have important implications for bridging the knowledge on turbulence and low-Reynolds-number emulsion flows to turbulent emulsion flows.

preprint2020arXiv

Lattice Boltzmann simulations of droplet breakup in confined and time-dependent flows

We study droplet dynamics and breakup in generic time-dependent flows via a multicomponent Lattice Boltzmann algorithm, with emphasis on flow start up conditions. We first study droplet breakup in a confined oscillatory shear flow via two different protocols. In one set up, we start from an initially spherical droplet and turn on the flow abruptly ("shock method"); in the other protocol, we start from an initially spherical droplet as well, but we progressively increase the amplitude of the flow, by allowing the droplet to relax to the steady state for each increase in amplitude, before increasing the flow amplitude again ("relaxation method"). The two protocols are shown to produce substantially different breakup scenarios. The mismatch between these two protocols is also studied for variations in the flow topology, the degree of confinement and the inertia of the fluid. All results point to the fact that under extreme conditions of confinement the relaxation protocols can drive the droplets into metastable states, which break only for very intense flow amplitudes, but their stability is prone to external perturbations, such as an oscillatory driving force.

preprint2020arXiv

Pedestrian orientation dynamics from high-fidelity measurements

We investigate in real-life conditions and with very high accuracy the dynamics of body rotation, or yawing, of walking pedestrians - an highly complex task due to the wide variety in shapes, postures and walking gestures. We propose a novel measurement method based on a deep neural architecture that we train on the basis of generic physical properties of the motion of pedestrians. Specifically, we leverage on the strong statistical correlation between individual velocity and body orientation: the velocity direction is typically orthogonal with respect to the shoulder line. We make the reasonable assumption that this approximation, although instantaneously slightly imperfect, is correct on average. This enables us to use velocity data as training labels for a highly-accurate point-estimator of individual orientation, that we can train with no dedicated annotation labor. We discuss the measurement accuracy and show the error scaling, both on synthetic and real-life data: we show that our method is capable of estimating orientation with an error as low as 7.5 degrees. This tool opens up new possibilities in the studies of human crowd dynamics where orientation is key. By analyzing the dynamics of body rotation in real-life conditions, we show that the instantaneous velocity direction can be described by the combination of orientation and a random delay, where randomness is provided by an Ornstein-Uhlenbeck process centered on an average delay of 100ms. Quantifying these dynamics could have only been possible thanks to a tool as precise as that proposed.

preprint2019arXiv

Modelling sea ice and melt ponds evolution: sensitivity to microscale heat transfer mechanisms

We present a mathematical model describing the evolution of sea ice and meltwater during summer. The system is described by two coupled partial differential equations for the ice thickness $h$ and pond depth $w$ fields. We test the sensitivity of the model to variations of parameters controlling fluid-dynamic processes at the pond level, namely the variation of turbulent heat flux with pond depth and the lateral melting of ice enclosing a pond. We observe that different heat flux scalings determine different rates of total surface ablations, while the system is relatively robust in terms of probability distributions of pond surface areas. Finally, we study pond morphology in terms of fractal dimensions, showing that the role of lateral melting is minor, whereas there is evidence of an impact from the initial sea ice topography.

preprint2019arXiv

The statistical properties of turbulence in the presence of a smart small-scale control

By means of high-resolution numerical simulations, we compare the statistical properties of homogeneous and isotropic turbulence to those of the Navier-Stokes equation where small-scale vortex filaments are strongly depleted, thanks to a non-linear extra viscosity acting preferentially on high vorticity regions. We show that the presence of such smart small-scale drag can strongly reduce intermittency and non-Gaussian fluctuations. Our results pave the way towards a deeper understanding on the fundamental role of degrees of freedom in turbulence as well as on the impact of (pseudo)coherent structures on the statistical small-scale properties. Our work can be seen as a first attempt to develop smart-Lagrangian forcing/drag mechanisms to control turbulence.

preprint2019arXiv

Topological structure and dynamics of three dimensional active nematics

Point-like motile topological defects control the universal dynamics of diverse two-dimensional active nematics ranging from shaken granular rods to cellular monolayers. A comparable understanding in higher dimensions has yet to emerge. We report the creation of three-dimensional active nematics by dispersing extensile microtubule bundles in a passive colloidal liquid crystal. Light-sheet microscopy reveals the millimeter-scale structure of active nematics with a single bundle resolution and the temporal evolution of the associated nematic director field. The dominant excitations of three-dimensional active nematics are extended charge-neutral disclination loops that undergo complex dynamics and recombination events. These studies introduce a new class of non-equilibrium systems whose turbulent-like dynamics arises from the interplay between internally generated active stresses, the chaotic flows and the topological structure of the constituent defects.

preprint2019arXiv

Unified theoretical and experimental view on transient shear banding

Dense emulsions, colloidal gels, microgels, and foams all display a solid-like behavior at rest characterized by a yield stress, above which the material flows like a liquid. Such a fluidization transition often consists of long-lasting transient flows that involve shear-banded velocity profiles. The characteristic time for full fluidization, $τ_\text{f}$, has been reported to decay as a power-law of the shear rate $\dot γ$ and of the shear stress $σ$ with respective exponents $α$ and $β$. Strikingly, the ratio of these exponents was empirically observed to coincide with the exponent of the Herschel-Bulkley law that describes the steady-state flow behavior of these complex fluids. Here we introduce a continuum model, based on the minimization of a "free energy", that captures quantitatively all the salient features associated with such \textit{transient} shear-banding. More generally, our results provide a unified theoretical framework for describing the yielding transition and the steady-state flow properties of yield stress fluids.

preprint2018arXiv

Accurate pedestrian localization in overhead depth images via Height-Augmented HOG

We tackle the challenge of reliably and automatically localizing pedestrians in real-life conditions through overhead depth imaging at unprecedented high-density conditions. Leveraging upon a combination of Histogram of Oriented Gradients-like feature descriptors, neural networks, data augmentation and custom data annotation strategies, this work contributes a robust and scalable machine learning-based localization algorithm, which delivers near-human localization performance in real-time, even with local pedestrian density of about 3 ped/m2, a case in which most state-of-the art algorithms degrade significantly in performance.

preprint2016arXiv

A Lagrangian model of copepod dynamics: Clustering by escape jumps in turbulence

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 per- turbations, 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 behaviour 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 behaviour 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 tur- bulence. 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

preprint2016arXiv

Coherent structures and extreme events in rotating multiphase turbulent flows

By using direct numerical simulations (DNS) at unprecedented resolution we study turbulence under rotation in the presence of simultaneous direct and inverse cascades. The accumulation of energy at large scale leads to the formation of vertical coherent regions with high vorticity oriented along the rotation axis. By seeding the flow with millions of inertial particles, we quantify -for the first time- the effects of those coherent vertical structures on the preferential concentration of light and heavy particles. Furthermore, we quantitatively show that extreme fluctuations, leading to deviations from a normal-distributed statistics, result from the entangled interaction of the vertical structures with the turbulent background. Finally, we present the first-ever measurement of the relative importance between Stokes drag, Coriolis force and centripetal forces along the trajectories of inertial particles. We discover that vortical coherent structures lead to unexpected diffusion properties for heavy and light particles in the directions parallel and perpendicular to the rotation axis.

preprint2016arXiv

Earthquake statistics inferred from plastic events in soft-glassy materials

We propose a new approach for generating synthetic earthquake catalogues based on the physics of soft glasses. The continuum approach produces yield-stress materials based on Lattice-Boltzmann simulations. We show that, if the material is stimulated below yield stress, plastic events occur, which have strong similarities with seismic events. Based on a suitable definition of displacement in the continuum, we show that the plastic events obey a Gutenberg-Richter law with exponents similar to those for real earthquakes. We further find that average acceleration, energy release, stress drop and recurrence times scale with the same exponent. The approach is fully self-consistent and all quantities can be calculated at all scales without the need of ad hoc friction or statistical laws. We therefore suggest that our approach may lead to new insight into understanding of the physics connecting the micro and macro scale of earthquakes.

preprint2016arXiv

Finite-size effects on bacterial population expansion under controlled fow conditions

The expansion of biological species in natural environments is usually described as the combined effect of individual spatial dispersal and growth. In the case of aquatic ecosystems flow transport can also be extremely relevant as an extra, advection induced, dispersal factor. There is a lack of reproducible experimental studies on biological fronts of living organisms in controlled streaming habitats. It is thus not clear if, and to which extent, the current theoretical and experimental knowledge on advective-reactive-diffusive fronts for chemical reactions can also apply to the expansion of biological populations. We designed and assembled a dedicated microfluidic device to control and quantify the expansion of populations of $E.coli$ bacteria under both co-flowing and counter-flowing conditions, measuring the front speed at varying intensity of the imposed flow. At variance with respect to the case of autocatalytic reactions, we measure that almost irrespective of the counter-flow velocity, the front speed remains finite at a constant positive value. A simple model incorporating growth, dispersion and drift on finite-size hard beads allows to explain this finding as due to a finite volume effect of the bacteria. This indicates that models based on the Fisher-Kolmogorov-Petrovsky-Piscounov equation (FKPP) that ignore the finite size of organisms may be inaccurate to describe the physics of spatial growth dynamics of bacteria.

preprint2016arXiv

Preferential rotation of chiral dipoles in isotropic turbulence

Particles in the shape of chiral dipoles show a preferential rotation in three dimensional homogeneous isotropic turbulence. A chiral dipole consists of a rod with two helices of opposite handedness, one at each end. We can use 3d printing to fabricate these particles with length in the inertial range and track their rotations in a turbulent flow between oscillating grids. High aspect ratio chiral dipoles will align with the extensional eigenvectors of the strain rate tensor and the helical ends will respond to the strain field by spinning around its long axis. The mean of the measured spinning rate is non-zero and reflects the average stretching the particles experience. We use Stokesian dynamics simulations of chiral dipoles in pure strain flow to quantify the dependence of spinning on particle shape. Based on the known response to pure strain, we build a model that gives the spinning rate of small chiral dipoles using Lagrangian velocity gradients from high resolution direct numerical simulations. The statistics of chiral dipole spinning determined with this model show surprisingly good agreement with the measured spinning of much larger chiral dipoles in the experiments.

preprint2016arXiv

Turbulence on a Fractal Fourier set

A novel investigation of the nature of intermittency in incompressible, homogeneous and isotropic turbulence is performed by a numerical study of the Navier-Stokes equations constrained on a fractal Fourier set. The robustness of the energy transfer and of the vortex stretching mechanisms is tested by changing the fractal dimension, D, from the original three dimensional case to a strongly decimated system with D=2.5, where only about $3\%$ of the Fourier modes interact. This is a unique methodology to probe the statistical properties of the turbulent energy cascade, without breaking any of the original symmetries of the equations. While the direct energy cascade persists, deviations from the Kolmogorov scaling are observed in the kinetic energy spectra. A model in terms of a correction with a linear dependency on the co-dimension of the fractal set, $E(k) \sim k^{-5/3 + 3-D}$, explains the results. At small scales, the intermittency of the vorticity field is observed to be quasi-singular as a function of the fractal mode reduction, leading to an almost Gaussian statistics already at $D \sim 2.98$. These effects must be connected to a genuine modification in the triad-to-triad nonlinear energy transfer mechanism.

preprint2015arXiv

Asymmetric pedestrian dynamics on a staircase landing from continuous measurements

We investigate via extensive experimental data the dynamics of pedestrians walking in a corridor-shaped landing in a building at Eindhoven University of Technology. With year-long automatic measurements employing a Microsoft KinectTM 3D-range sensor and ad hoc tracking techniques, we acquired few hundreds of thousands pedestrian trajectories in real-life conditions. Here we discuss the asymmetric features of the dynamics in the two walking directions with respect to the flights of stairs (i.e. ascending or descending). We provide a detailed analysis of position and speed fields for the cases of pedestrians walking alone undisturbed and for couple of pedestrians in counter-flow. Then, we show average walking velocities exploring all the observed combinations in terms of numbers of pedestrians and walking directions.

preprint2015arXiv

Axisymmetric multiphase Lattice Boltzmann method for generic equations of state

We present an axisymmetric lattice Boltzmann model based on the Kupershtokh et al. multiphase model that is capable of solving liquid-gas density ratios up to $10^3$. Appropriate source terms are added to the lattice Boltzmann evolution equation to fully recover the axisymmetric multiphase conservation equations. We validate the model by showing that a stationary droplet obeys the Young-Laplace law, comparing the second oscillation mode of a droplet with respect to an analytical solution and showing correct mass conservation of a propagating density wave.

preprint2015arXiv

Clustering of vertically constrained passive particles in homogeneous, isotropic turbulence

We analyze the dynamics of small particles vertically confined, by means of a linear restoring force, to move within a horizontal fluid slab in a three-dimensional (3D) homogeneous isotropic turbulent velocity field. The model that we introduce and study is possibly the simplest description for the dynamics of small aquatic organisms that, due to swimming, active regulation of their buoyancy, or any other mechanism, maintain themselves in a shallow horizontal layer below the free surface of oceans or lakes. By varying the strength of the restoring force, we are able to control the thickness of the fluid slab in which the particles can move. This allows us to analyze the statistical features of the system over a wide range of conditions going from a fully 3D incompressible flow (corresponding to the case of no confinement) to the extremely confined case corresponding to a two-dimensional slice. The background 3D turbulent velocity field is evolved by means of fully resolved direct numerical simulations. Whenever some level of vertical confinement is present, the particle trajectories deviate from that of fluid tracers and the particles experience an effectively compressible velocity field. Here, we have quantified the compressibility, the preferential concentration of the particles, and the correlation dimension by changing the strength of the restoring force. The main result is that there exists a particular value of the force constant, corresponding to a mean slab depth approximately equal to a few times the Kolmogorov length scale, that maximizes the clustering of the particles.

preprint2015arXiv

On the Peterlin approximation for turbulent flows of polymer solutions

We study the impact of the Peterlin approximation on the statistics of the end-to-end separation of poly- mers in a turbulent flow. The FENE and FENE-P models are numerically integrated along a large number of Lagrangian trajectories resulting from a direct numerical simulation of three-dimensional homogeneous isotropic turbulence. Although the FENE-P model yields results in qualitative agreement with those of the FENE model, quantitative differences emerge. The steady-state probability of large extensions is overesti- mated by the FENE-P model. The alignment of polymers with the eigenvectors of the rate-of-strain tensor and with the direction of vorticity is weaker when the Peterlin approximation is used. At large Weissenberg numbers, both the correlation times of the extension and of the orientation of polymers are underestimated by the FENE-P model.

preprint2014arXiv

A new assessment of the second order moment of Lagrangian velocity increments in turbulence

The behavior of the second-order Lagrangian structure functions on state-of-the-art numerical data both in two and three dimensions is studied. On the basis of a phenomenological connection between Eulerian space-fluctuations and the Lagrangian time-fluctuations, it is possible to rephrase the Kolmogorov $4/5$-law into a relation predicting the linear (in time) scaling for the second order Lagrangian structure function. When such a function is directly observed on current experimental or numerical data, it does not clearly display a scaling regime. A parameterization of the Lagrangian structure functions based on Batchelor model is introduced and tested on data for $3d$ turbulence, and for $2d$ turbulence in the inverse cascade regime. Such parameterization supports the idea, previously suggested, that both Eulerian and Lagrangian data are consistent with a linear scaling plus finite-Reynolds number effects affecting the small- and large-time scales. When large-time saturation effects are properly accounted for, compensated plots show a detectable plateau already at the available Reynolds number. Furthermore, this parameterization allows us to make quantitative predictions on the Reynolds number value for which Lagrangian structure functions are expected to display a scaling region. Finally, we show that this is also sufficient to predict the anomalous dependency of the normalized root mean squared acceleration as a function of the Reynolds number, without fitting parameters.

preprint2013arXiv

Accurate lubrication corrections for spherical and non-spherical particles in discretized fluid simulations

Discretized fluid solvers coupled to a Newtonian dynamics method are a popular tool to study suspension flow. As any simulation technique with finite resolution, the lattice Boltzmann method, when coupled to discrete particles using the momentum exchange method, resolves the diverging lubrication interactions between surfaces near contact only insufficiently. For spheres, it is common practice to account for surface-normal lubrication forces by means of an explicit correction term. A method that additionally covers all further singular interactions for spheres is present in the literature as well as a link-based approach that allows for more general shapes but does not capture non-normal interactions correctly. In this paper, lattice-independent lubrication corrections for aspherical particles are outlined, taking into account all leading divergent interaction terms. An efficient implementation for arbitrary spheroids is presented and compared to purely normal and link-based models. Good consistency with Stokesian dynamics simulations of spheres is found. The non-normal interactions affect the viscosity of suspensions of spheres at volume fractions Φ>= 0.3 but already at Φ>= 0.2 for spheroids. Regarding shear-induced diffusion of spheres, a distinct effect is found at 0.1 <= Φ<= 0.5 and even increasing the resolution of the radius to 8 lattice units is no substitute for an accurate modeling of non-normal interactions.

preprint2013arXiv

Axisymmetric multiphase lattice Boltzmann method

A lattice Boltzmann method for axisymmetric multiphase flows is presented and validated. The method is capable of accurately modeling flows with variable density. We develop the classic Shan-Chen multiphase model [ Phys. Rev. E 47 1815 (1993)] for axisymmetric flows. The model can be used to efficiently simulate single and multiphase flows. The convergence to the axisymmetric Navier-Stokes equations is demonstrated analytically by means of a Chapmann-Enskog expansion and numerically through several test cases. In particular, the model is benchmarked for its accuracy in reproducing the dynamics of the oscillations of an axially symmetric droplet and on the capillary breakup of a viscous liquid thread. Very good quantitative agreement between the numerical solutions and the analytical results is observed.

preprint2013arXiv

Heat flux scaling in turbulent Rayleigh-Bénard convection with an imposed longitudinal wind

We present a numerical study of Rayleigh-Bénard convection disturbed by a longitudinal wind. Our results show that under the action of the wind, the vertical heat flux through the cell initially decreases, due to the mechanism of plumes-sweeping, and then increases again when turbulent forced convection dominates over the buoyancy. As a result, the Nusselt number is a non-monotonic function of the shear Reynolds number. We provide a simple model that captures with good accuracy all the dynamical regimes observed. We expect that our findings can lead the way to a more fundamental understanding of the of the complex interplay between mean-wind and plumes ejection in the Rayleigh-Bénard phenomenology.

preprint2013arXiv

Lagrangian single particle turbulent statistics through the Hilbert-Huang Transform

The Hilbert-Huang transform is applied to analyze single particle Lagrangian velocity data from numerical simulations of hydrodynamic turbulence. The velocity trajectory is described in terms of a set of intrinsic mode functions, C_{i}(t), and of their instantaneous frequency, ω_{i}(t). On the basis of this decomposition we define the ω-conditioned statistical moments of the C_{i} modes, named q-order Hilbert Spectra (HS). We show that such new quantities have enhanced scaling properties as compared to traditional Fourier transform- or correlation-based (Structure Functions) statistical indicators, thus providing better insights into the turbulent energy transfer process. We present a clear empirical evidence that the energy-like quantity, i.e. the second-order HS, displays a linear scaling in time in the inertial range, as expected from dimensional analysis and never observed before. We also measure high order moment scaling exponents in a direct way, without resorting the Extended Self Similarity (ESS) procedure. This leads to a new estimate of the Lagrangian structure functions exponents which are consistent with the multifractal prediction in the Lagrangian frame as proposed in [Biferale et al., Phys. Rev. Lett. vol. 93, 064502 (2004)].

preprint2013arXiv

Lattice Boltzmann method to study the contraction of a viscous ligament

We employ a recently formulated axisymmetric version of the multiphase Shan-Chen (SC) lattice Boltzmann method (LBM) [Srivastava et al, in preparation (2013)] to simulate the contraction of a liquid ligament. We compare the axisymmetric LBM simulation against the slender jet (SJ) approximation model [T. Driessen and R. Jeurissen, IJCFD {\bf 25}, 333 (2011)]. We compare the retraction dynamics of the tail-end of the liquid ligament from the LBM simulation, the SJ model, Flow3D simulations and a simple model based on the force balance (FB). We find good agreement between the theoretical prediction (FB), the SJ model, and the LBM simulations.

preprint2013arXiv

Spinodal decomposition in homogeneous and isotropic turbulence

We study the competition between domain coarsening in a symmetric binary mixtures below the critical temperature and turbulent fluctuations. We find that the coarsening process is arrested in presence of turbulence. The physics of the process shares remarkable similarities with the behaviour of diluted turbulent emulsions and the arrest length scale can be estimated with an argument similar to the one proposed by Kolmogorov and Hinze for the maximal stability diameter of droplets in turbulence. Although in the absence of flow the microscopic diffusion constant is negative, turbulence does effectively arrest the inverse cascade of concentration fluctuations by making the low wavelength diffusion constant positive for scales above the Hinze length.

preprint2013arXiv

Stability of viscous long liquid filaments

We study the collapse of an axisymmetric liquid filament both analytically and by means of a numerical model. The liquid filament, also known as ligament, may either collapse stably into a single droplet or break up into multiple droplets. The dynamics of the filament are governed by the viscosity and the aspect ratio, and the initial perturbations of its surface. We find that the instability of long viscous filaments can be completely explained by the Rayleigh-Plateau instability, whereas a low viscous filament can also break up due to end pinching. We analytically derive the transition between stable collapse and breakup in the Ohnesorge number versus aspect ratio phase space. Our result is confirmed by numerical simulations based on the slender jet approximation and explains recent experimental findings by Castrejon-Pita et al., PRL 108, 074506 (2012).

preprint2013arXiv

Towards a continuum model for particle-induced velocity fluctuations in suspension flow through a stenosed geometry

Non-particulate continuum descriptions allow for computationally efficient modeling of suspension flows at scales that are inaccessible to more detailed particulate approaches. It is well known that the presence of particles influences the effective viscosity of a suspension and that this effect has thus to be accounted for in macroscopic continuum models. The present paper aims at developing a non-particulate model that reproduces not only the rheology but also the cell-induced velocity fluctuations, responsible for enhanced diffusivity. The results are obtained from a coarse-grained blood model based on the lattice Boltzmann method. The benchmark system comprises a flow between two parallel plates with one of them featuring a smooth obstacle imitating a stenosis. Appropriate boundary conditions are developed for the particulate model to generate equilibrated cell configurations mimicking an infinite channel in front of the stenosis. The averaged flow field in the bulk of the channel can be described well by a non-particulate simulation with a matched viscosity. We show that our proposed phenomenological model is capable to reproduce many features of the velocity fluctuations.

preprint2012arXiv

Analytical modeling for the heat transfer in sheared flows of nanofluids

We developed a model for the enhancement of the heat flux by spherical and elongated nano- particles in sheared laminar flows of nano-fluids. Besides the heat flux carried by the nanoparticles the model accounts for the contribution of their rotation to the heat flux inside and outside the particles. The rotation of the nanoparticles has a twofold effect, it induces a fluid advection around the particle and it strongly influences the statistical distribution of particle orientations. These dynamical effects, which were not included in existing thermal models, are responsible for changing the thermal properties of flowing fluids as compared to quiescent fluids. The proposed model is strongly supported by extensive numerical simulations, demonstrating a potential increase of the heat flux far beyond the Maxwell-Garnet limit for the spherical nanoparticles. The road ahead which should lead towards robust predictive models of heat flux enhancement is discussed.

preprint2012arXiv

Growth, competition and cooperation in spatial population genetics

We study an individual based model describing competition in space between two different alleles. Although the model is similar in spirit to classic models of spatial population genetics such as the stepping stone model, here however space is continuous and the total density of competing individuals fluctuates due to demographic stochasticity. By means of analytics and numerical simulations, we study the behavior of fixation probabilities, fixation times, and heterozygosity, in a neutral setting and in cases where the two species can compete or cooperate. By concluding with examples in which individuals are transported by fluid flows, we argue that this model is a natural choice to describe competition in marine environments.

preprint2012arXiv

How gravity and size affect the acceleration statistics of bubbles in turbulence

We report results from the first systematic Lagrangian experimental investigation in the regime of very light (air bubbles in water) and large particles (D/η >> 1) in turbulence. Using a traversing camera setup and particle tracking, we study the Lagrangian acceleration statistics of ~3 mm diameter (D) bubbles in a water tunnel with nearly homogeneous and isotropic turbulence generated by an active-grid. The Reynolds number (Re_λ) is varied from 145 to 230, resulting in size ratios, D/η in the range of 7.3--12.5, where η is the Kolmogorov length scale. The experiments reveal that gravity increases the acceleration variance and reduces the intermittency of the PDF in the vertical direction. Once the gravity offset is subtracted, the variances of both the horizontal and vertical acceleration components are about 5\pm2 times larger than the one measured in the same flow for fluid tracers. Moreover, for these light particles, the experimental acceleration PDF shows a substantial reduction in intermittency at growing size ratios, in contrast to neutrally buoyant or heavy particles. All these results are closely matched by numerical simulations of finite-size bubbles with the Faxén corrections.

preprint2012arXiv

Population dynamics in compressible flows

Organisms often grow, migrate and compete in liquid environments, as well as on solid surfaces. However, relatively little is known about what happens when competing species are mixed and compressed by fluid turbulence. In these lectures we review our recent work on population dynamics and population genetics in compressible velocity fields of one and two dimensions. We discuss why compressible turbulence is relevant for population dynamics in the ocean and we consider cases both where the velocity field is turbulent and when it is static. Furthermore, we investigate populations in terms of a continuos density field and when the populations are treated via discrete particles. In the last case we focus on the competition and fixation of one species compared to another

preprint2012arXiv

Rotation Rate of Rods in Turbulent Fluid Flow

The rotational dynamics of anisotropic particles advected in a turbulent fluid flow are important in many industrial and natural setting. Particle rotations are controlled by small scale properties of turbulence that are nearly universal, and so provide a rich system where experiments can be directly compared with theory and simulations. Here we report the first three-dimensional experimental measurements of the orientation dynamics of rod-like particles as they are advected in a turbulent fluid flow. We also present numerical simulations that show good agreement with the experiments and allow extension to a wide range of particle shapes. Anisotropic tracer particles preferentially sample the flow since their orientations become correlated with the velocity gradient tensor. The rotation rate is heavily influenced by this preferential alignment, and the alignment depends strongly on particle shape.

preprint2012arXiv

Unified framework for a side-by-side comparison of different multicomponent algorithms: lattice Boltzmann vs. phase field model

Lattice Boltzmann Models (LBM) and Phase Field Models (PFM) are two of the most widespread approaches for the numerical study of multicomponent fluid systems. Both methods have been successfully employed by several authors but, despite their popularity, still remains unclear how to properly compare them and how they perform on the same problem. Here we present a unified framework for the direct (one-to-one) comparison of the multicomponent LBM against the PFM. We provide analytical guidelines on how to compare the Shan-Chen (SC) lattice Boltzmann model for non-ideal multicomponent fluids with a corresponding free energy (FE) lattice Boltzmann model. Then, in order to properly compare the LBM vs. the PFM, we propose a new formulation for the free energy of the Cahn-Hilliard/Navier-Stokes equations. Finally, the LBM model is numerically compared with the corresponding phase field model solved by means of a pseudo-spectral algorithm. This work constitute a first attempt to set the basis for a quantitative comparison between different algorithms for multicomponent fluids. We limit our scope to the few of the most common variants of the two most widespread methodologies, namely the lattice Boltzmann model (SC and FE variants) and the Phase Field Model.

preprint2011arXiv

Axially-homogeneous Rayleigh-Benard convection in a cylindrical cell

Previous numerical studies have shown that the "ultimate regime of thermal convection" can be attained in a Rayleigh-Benard cell when the kinetic and thermal boundary layers are eliminated by replacing the walls with periodic boundary conditions (homogeneous Rayleigh-Benard convection). Then, the heat transfer scales like Nu ~ Ra^{1/2} and turbulence intensity as Re ~ Ra^{1/2}, where the Rayleigh number Ra indicates the strength of the driving force. However, experiments never operate in unbounded domains and it is important to understand how confinement might alter the approach to this ultimate regime. Here we consider homogeneous Rayleigh-Benard convection in a laterally confined geometry - a small aspect-ratio vertical cylindrical cell - and show evidence of the ultimate regime as Ra is increased: In spite of the confinement and the resulting kinetic boundary layers, we still find Nu ~ Re ~ Ra^{1/2}. The system supports exact solutions composed of modes of exponentially growing vertical velocity and temperature fields, with Ra as the critical parameter determining the properties of these modes. Counterintuitively, in the low Ra regime, or for very narrow cylinders, the numerical simulations are susceptible to these solutions which can dominate the dynamics and lead to very high and unsteady heat transfer. As Ra is increased, interaction between modes stabilizes the system, evidenced by the increasing homogeneity and reduced fluctuations in the r.m.s. velocity and temperature fields. We also test that physical results become independent of the periodicity length of the cylinder, a purely numerical parameter, as the aspect ratio is increased.

preprint2011arXiv

Cumulative compressibility effects on population dynamics in turbulent flows

Bacteria and plankton populations living in oceans and lakes reproduce and die under the in- fluence of turbulent currents. Turbulent transport interacts in a complex way with the dynamics of populations because the typical reproduction time of microorganism is within the inertial range of turbulent time scales. In the present manuscript we quantitatively investigate the effect of flow compressibility on the dynamics of populations. While a small compressibility can be induced by several physical mechanisms, like density mismatch or the finite size of microorganisms with respect to the fluid turbulence, its effect on the the carrying capacity of the ecosystem can be dramatic. We report, for the first time, how a small compressibility can produce a sizeable reduction in the carrying capacity, due to an integrated effect made possible by the long replication times of the organisms with respect to turbulent time scales. A full statistical quantification of the fluctuations of population concentration field leads to data collapse over a broad range in parameter space.

preprint2011arXiv

Droplet size distribution in homogeneous isotropic turbulence

We study the physics of droplet breakup in a statistically stationary homogeneous and isotropic turbulent flow by means of high resolution numerical investigations based on the multicomponent lattice Boltzmann method. We verified the validity of the criterion proposed by Hinze (1955) for droplet breakup and we measured the full probability distribution function (pdf) of droplets radii at different Reynolds numbers and for different volume fraction. By means of a Lagrangian tracking we could follow individual droplets along their trajectories, define a local Weber number based on the velocity gradients and study its cross-correlation with droplet deformation.

preprint2011arXiv

Inverse energy cascade in three-dimensional isotropic turbulence

In turbulent flows kinetic energy is spread by nonlinear interactions over a broad range of scales. Energy transfer may proceed either toward small scales or in the reverse direction. The latter case is peculiar of two-dimensional (2D) flows. Interestingly, a reversal of the energy flux is observed also in three-dimensional (3D) geophysical flows under rotation and/or confined in thin layers. The question is whether this phenomenon is enforced solely by external anisotropic mechanisms or it is intimately embedded in the Navier-Stokes (NS) equations. Here we show that an inverse energy cascade occurs also in 3D isotropic flow. The flow is obtained from a suitable surgery of the NS equations, keeping only triadic interactions between sign-defined helical modes, preserving homogeneity and isotropy and breaking reflection invariance. Our findings highlight the role played by helicity in the energy transfer process and show that both 2D and 3D properties naturally coexist in all flows in nature.

preprint2011arXiv

Multi-Time Multi-Scale Correlation Functions in Hydrodynamic Turbulence

High Reynolds numbers Navier-Stokes equations are believed to break self-similarity concerning both spatial and temporal properties: correlation functions of different orders exhibit distinct decorrelation times and anomalous spatial scaling properties. Here, we present a systematic attempt to measure multi-time and multi-scale correlations functions, by using high Reynolds numbers numerical simulations of fully homogeneous and isotropic turbulent flow. The main idea is to set-up an ensemble of probing stations riding the flow, i.e. measuring correlations in a reference frame centered on the trajectory of distinct fluid particles (the quasi-Lagrangian reference frame introduced by Belinicher & L'vov, Sov. Phys. JETP 66, 303 (1987)). In this way we reduce the large-scale sweeping and measure the non-trivial temporal dynamics governing the turbulent energy transfer from large to small scales. We present evidences of the existence of dynamic multiscaling: multi-time correlation functions are characterized by an infinite set of characteristic times.

preprint2011arXiv

Simulations of Blood Flow in Plain Cylindrical and Constricted Vessels with Single Cell Resolution

Understanding the physics of blood is challenging due to its nature as a suspension of soft particles and the fact that typical problems involve different scales. This is valid also for numerical investigations. In fact, many computational studies either neglect the existence of discrete cells or resolve relatively few cells very accurately. The authors recently developed a simple and highly efficient yet still particulate model with the aim to bridge the gap between currently applied methods. The present work focuses on its applicability to confined flows in vessels of diameters up to 100 micrometres. For hematocrit values below 30 percent, a dependence of the apparent viscosity on the vessel diameter in agreement with experimental literature data is found.

preprint2010arXiv

A simplified particulate model for coarse-grained hemodynamics simulations

Human blood flow is a multi-scale 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.

preprint2010arXiv

Droplet breakup in homogeneous and isotropic turbulence

This fluid dynamics video shows the breakup of a droplet in a stationary homogeneous and isotropic turbulent flow. We consider droplets with the same density of the transporting fluid. The droplets and the fluid are numerically modelled by means of a multicompo- nent Lattice-Boltzmann method. The turbulent fluid is maintained through a large scale stirring force and the radius of stable droplets, for the parameters in our simulation, is larger than the Kolmogorov scale. Events of droplet deformation, break-up and aggregation are clearly visible from the movie. With the present database droplet evo- lution can be studied from both an Eulerian and Lagrangian point of view. The Kolmogorov-Hinze criteria for droplets break-up can be tested also by means of simulations with different viscosity contrast between the two components.

preprint2010arXiv

Impact of trailing wake drag on the statistical properties and dynamics of finite-sized particle in turbulence

We study by means of an Eulerian-Lagrangian model the statistical properties of velocity and acceleration of a neutrally-buoyant finite-sized particle in a turbulent flow statistically homogeneous and isotropic. The particle equation of motion, beside added mass and steady Stokes drag, keeps into account the unsteady Stokes drag force - known as Basset-Boussinesq history force - and the non-Stokesian drag based on Schiller-Naumann parametrization, together with the finite-size Faxen corrections. We focus on the case of flow at low Taylor-Reynolds number, Re_lambda ~ 31, for which fully resolved numerical data which can be taken as a reference are available (Homann & Bec 651 81-91 J. Fluid Mech. (2010)). Remarkably, we show that while drag forces have always minor effects on the acceleration statistics, their role is important on the velocity behavior. We propose also that the scaling relations for the particle velocity variance as a function of its size, which have been first detected in fully resolved simulations, does not originate from inertial-scale properties of the background turbulent flow but it is likely to arise from the non-Stokesian component of the drag produced by the wake behind the particle. Furthermore, by means of comparison with fully resolved simulations, we show that the Faxen correction to the added mass has a dominant role in the particle acceleration statistics even for particle with size in the inertial range.

preprint2010arXiv

Lattice Boltzmann Methods for thermal flows: continuum limit and applications to compressible Rayleigh-Taylor systems

We compute the continuum thermo-hydrodynamical limit of a new formulation of lattice kinetic equations for thermal compressible flows, recently proposed in [Sbragaglia et al., J. Fluid Mech. 628 299 (2009)]. We show that the hydrodynamical manifold is given by the correct compressible Fourier- Navier-Stokes equations for a perfect fluid. We validate the numerical algorithm by means of exact results for transition to convection in Rayleigh-Bénard compressible systems and against direct comparison with finite-difference schemes. The method is stable and reliable up to temperature jumps between top and bottom walls of the order of 50% the averaged bulk temperature. We use this method to study Rayleigh-Taylor instability for compressible stratified flows and we determine the growth of the mixing layer at changing Atwood numbers up to At ~ 0.4. We highlight the role played by the adiabatic gradient in stopping the mixing layer growth in presence of high stratification and we quantify the asymmetric growth rate for spikes and bubbles for two dimensional Rayleigh- Taylor systems with resolution up to Lx \times Lz = 1664 \times 4400 and with Rayleigh numbers up to Ra ~ 2 \times 10^10.

preprint2010arXiv

Population dynamics at high Reynolds number

We study the statistical properties of population dynamics evolving in a realistic two-dimensional compressible turbulent velocity field. We show that the interplay between turbulent dynamics and population growth and saturation leads to quasi-localization and a remarkable reduction in the carrying capacity. The statistical properties of the population density are investigated and quantified via multifractal scaling analysis. We also investigate numerically the singular limit of negligibly small growth rates and delocalization of population ridges triggered by uniform advection.

preprint2010arXiv

Rotational behavior of red blood cells in suspension---a mesoscale simulation study

The nature of blood as a suspension of red blood cells makes computational hemodynamics a demanding task. Our coarse-grained blood model, which builds on a lattice Boltzmann method for soft particle suspensions, enables the study of the collective behavior of the order of 10^6 cells in suspension. After demonstrating the viscosity measurement in Kolmogorov flow, we focus on the statistical analysis of the cell orientation and rotation in Couette flow. We quantify the average inclination with respect to the flow and the nematic order as a function of shear rate and hematocrit. We further record the distribution of rotation periods around the vorticity direction and find a pronounced peak in the vicinity of the theoretical value for free model cells even though cell-cell interactions manifest themselves in a substantial width of the distribution.

preprint2008arXiv

Acceleration statistics of finite-sized particles in turbulent flow: the role of Faxen forces

The dynamics of particles in turbulence when the particle-size is larger than the dissipative scale of the carrier flow is 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 model is obtained once Faxen corrections are included.

preprint2008arXiv

Dimensionality and morphology of particle and bubble clusters in turbulent flow

We conduct numerical experiments to investigate the spatial clustering of particles and bubbles in simulations of homogeneous and isotropic turbulence. Varying the Stokes parameter and the densities, striking differences in the clustering of the particles can be observed. To quantify these visual findings we use the Kaplan--Yorke dimension. This local scaling analysis shows a dimension of approximately 1.4 for the light bubble distribution, whereas the distribution of very heavy particles shows a dimension of approximately 2.4. However, clearly separate parameter combinations yield the same dimensions. To overcome this degeneracy and to further develop the understanding of clustering, we perform a morphological (geometrical and topological) analysis of the particle distribution. For such an analysis, Minkowski functionals have been successfully employed in cosmology, in order to quantify the global geometry and topology of the large-scale distribution of galaxies. In the context of dispersed multiphase flow, these Minkowski functionals -- being morphological order parameters -- allow us to discern the filamentary structure of the light particle distribution from the wall-like distribution of heavy particles around empty interconnected tunnels.

preprint2001arXiv

Evidences of Bolgiano scaling in 3D Rayleigh-Benard convection

We present new results from high-resolution high-statistics direct numerical simulations of a tri-dimensional convective cell. We test the fundamental physical picture of the presence of both a Bolgiano-like and a Kolmogorov-like regime. We find that the dimensional predictions for these two distinct regimes (characterized respectively by an active and passive role of the temperature field) are consistent with our measurements.

preprint1997arXiv

FFT for the APE Parallel Computer

We present a parallel FFT algorithm for SIMD systems following the `Transpose Algorithm' approach. The method is based on the assignment of the data field onto a 1-dimensional ring of systolic cells. The systolic array can be universally mapped onto any parallel system. In particular for systems with next-neighbour connectivity our method has the potential to improve the efficiency of matrix transposition by use of hyper-systolic communication. We have realized a scalable parallel FFT on the APE100/Quadrics massively parallel computer, where our implementation is part of a 2-dimensional hydrodynamics code for turbulence studies. A possible generalization to 4-dimensional FFT is presented, having in mind QCD applications.