• Elastoviscoplastic flow in porous media(1802.01993)

April 12, 2018 physics.flu-dyn
We investigate the elastoviscoplastic flow through porous media by numerical simulations. We solve the Navier-Stokes equations combined with the elastoviscoplastic model proposed by Saramito for the stress tensor evolution. In this model, the material behaves as a viscoelastic solid when unyielded, and as a viscoelastic Oldroyd-B fluid for stresses higher than the yield stress. The porous media is made of a symmetric array of cylinders, and we solve the flow in one periodic cell. We find that the solution is time-dependent even at low Reynolds numbers as we observe oscillations in time of the unyielded region especially at high Bingham numbers. The volume of the unyielded region slightly decreases with the Reynolds number and strongly increases with the Bingham number; up to 70% of the total volume is unyielded for the highest Bingham numbers considered here. The flow is mainly shear dominated in the yielded region, while shear and elongational flow are equally distributed in the unyielded region. We compute the relation between the pressure drop and the flow rate in the porous medium and present an empirical closure as function of the Bingham and Reynolds numbers. The apparent permeability, normalized with the case of Newtonian fluids, is shown to be greater than 1 at low Bingham numbers, corresponding to lower pressure drops due to the flow elasticity, and smaller than 1 for high Bingham numbers, indicating larger dissipation in the flow owing to the presence of the yielded regions. Finally we investigate the effect of the Weissenberg number on the distribution of the unyielded regions and on the pressure gradient.
• Inertial migration in dilute and semi-dilute suspensions of rigid particles in laminar square duct flow(1705.08734)

May 24, 2017 physics.flu-dyn
We study the inertial migration of finite-size neutrally buoyant spherical particles in dilute and semi-dilute suspensions in laminar square duct flow. We perform several direct numerical simulations using an immersed boundary method to investigate the effects of the bulk Reynolds number $Re_b$, particle Reynolds number $Re_p$ and duct to particle size ratio $h/a$ at different solid volume fractions $\phi$, from very dilute conditions to $20\%$. We show that the bulk Reynolds number $Re_b$ is the key parameter in inertial migration of particles in dilute suspensions. At low solid volume fraction ($\phi=0.4\%$) and low bulk Reynolds number $(Re_b= 144)$, particles accumulate at the center of the duct walls. As $Re_b$ is increased, the focusing position moves progressively towards the corners of the duct. At higher volume fractions, $\phi=5, 10$ and $20\%$, and in wider ducts with $Re_b=550$, particles are found to migrate away from the duct core towards the walls. In particular, for $\phi=5$ and $10\%$, particles accumulate preferentially at the corners. At the highest volume fraction considered, $\phi=20\%$, particles sample all the volume of the duct, with a lower concentration at the duct core. The presence of particles induces secondary cross-stream motions in the duct cross-section, for all $\phi$. The intensity of these secondary flows depends strongly on particle rotation rate, on the maximum concentration of particles in focusing positions, and on the solid volume fraction. We find that the secondary flow intensity increases with the volume fraction up to $\phi=5\%$. However, beyond $\phi=5\%$ excluded volume effects lead to a strong reduction of cross-stream velocities.
• Drag reduction in turbulent channel flow laden with finite-size oblate spheroids(1607.00679)

Jan. 24, 2017 physics.flu-dyn
We study suspensions of oblate rigid particles in a viscous fluid for different values of the particle volume fractions. Direct numerical simulations have been performed using a direct-forcing immersed boundary method to account for the dispersed phase, combined with a soft-sphere collision model and lubrication corrections for short-range particle-particle and particle-wall interactions. With respect to the single phase flow, we show that in flows laden with oblate spheroids the drag is reduced and the turbulent fluctuations attenuated. In particular, the turbulence activity decreases to lower values than those obtained by only accounting for the effective suspension viscosity. To explain the observed drag reduction we consider the particle dynamics and the interactions of the particles with the turbulent velocity field and show that the particle wall layer, previously observed and found to be responsible for the increased dissipation in suspensions of spheres, disappears in the case of oblate particles. These rotate significantly slower than spheres near the wall and tend to stay with their major axes parallel to the wall, which leads to a decrease of the Reynolds stresses and turbulence production and so to the overall drag reduction.
• Sedimentation of elongated non-motile prolate spheroids in homogenous isotropic turbulence(1611.10284)

Nov. 30, 2016 physics.flu-dyn
Phytoplankton are the foundation of aquatic food webs. Through photosynthesis, phytoplankton draw down CO2 at magnitudes equivalent to forests and other terrestrial plants and convert it to organic material that is then consumed by other organisms of phytoplankton in higher trophic levels. Mechanisms that affect local concentrations and velocities are of primary significance to many encounter-based processes in the plankton including prey-predator interactions, fertilization and aggregate formation. We report results from simulations of sinking phytoplankton, considered as elongated spheroids, in homogenous isotropic turbulence to answer the question of whether trajectories and velocities of sinking phytoplankton are altered by turbulence. We show in particular that settling spheroids with physical characteristics similar to those of diatoms weakly cluster and preferentially sample regions of down-welling flow, corresponding to an increase of the mean settling speed with respect to the mean settling speed in quiescent fluid. We explain how different parameters can affect the settling speed and what underlying mechanisms might be involved. Interestingly, we observe that the increase in the aspect ratio of the prolate spheroids can affect the clustering and the average settling speed of particles by two mechanisms: first is the effect of aspect ratio on the rotation rate of the particles, which saturates faster than the second mechanism of increasing drag anisotropy.
• Composition-tuned magneto-optical Kerr effect in L10-MnxGa films with giant perpendicular anisotropy(1601.07939)

Nov. 12, 2016 cond-mat.mtrl-sci
We report the large polar magnetooptical Kerr effect in L10-MnxGa epitaxial films with giant perpendicular magnetic anisotropy in a wide composition range. The Kerr rotation was enhanced by a factor of up to 10 by decreasing Mn atomic concentration, which most likely arises from the variation of the effective spin-orbit coupling strength, compensation effect of magnetic moments at different Mn atom sites, and overall strain. The Kerr ellipticity and the magnitude of the complex Kerr angle is found to have more complex composition-dependence that varies with the photon energy. These L10-MnxGa films show large Kerr rotation of up to 0.10o, high reflectivity of 35%-55% in a wide wavelength range of 400~850 nm, and giant magnetic anisotropic field of up to 210 kOe, making them an interesting material system for emerging spintronics and terahertz modulator applications.
• Turbulent channel flow of dense suspensions of neutrally-buoyant spheres(1409.1996)

Oct. 14, 2015 physics.flu-dyn
Dense particle suspensions are widely encountered in many applications and in environmental flows. While many previous studies investigate their rheological properties in laminar flows, little is known on the behaviour of these suspensions in the turbulent/inertial regime. The present study aims to fill this gap by investigating the turbulent flow of a Newtonian fluid laden with solid neutrally-buoyant spheres at relatively high volume fractions in a plane channel. Direct Numerical Simulation are performed in the range of volume fractions Phi=0-0.2 with an Immersed Boundary Method used to account for the dispersed phase. The results show that the mean velocity profiles are significantly altered by the presence of a solid phase with a decrease of the von Karman constant in the log-law. The overall drag is found to increase with the volume fraction, more than one would expect just considering the increase of the system viscosity due to the presence of the particles. At the highest volume fraction here investigated, Phi=0.2, the velocity fluctuation intensities and the Reynolds shear stress are found to decrease. The analysis of the mean momentum balance shows that the particle-induced stresses govern the dynamics at high Phi and are the main responsible of the overall drag increase. In the dense limit, we therefore find a decrease of the turbulence activity and a growth of the particle induced stress, where the latter dominates for the Reynolds numbers considered here.
• List of contributions from the CTA Consortium presented at the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.
• Dispersion of swimming algae in laminar and turbulent channel flows: consequences for photobioreactors(1205.2840)

Shear flow significantly affects the transport of swimming algae in suspension. For example, viscous and gravitational torques bias bottom-heavy cells to swim towards regions of downwelling fluid (gyrotaxis). It is necessary to understand how such biases affect algal dispersion in natural and industrial flows, especially in view of growing interest in algal photobioreactors. Motivated by this, we here study the dispersion of gyrotactic algae in laminar and turbulent channel flows using direct numerical simulation (DNS) and the analytical swimming dispersion theory of Bees and Croze (2010). Time-resolved dispersion measures are evaluated as functions of the Peclet and Reynolds numbers in upwelling and downwelling flows. For laminar flows, DNS results are compared with theory using competing descriptions of biased swimming cells in shear flow. Excellent agreement is found for predictions that employ generalized-Taylor-dispersion. The results highlight peculiarities of gyrotactic swimmer dispersion relative to passive tracers. In laminar downwelling flow the cell distribution drifts in excess of the mean flow, increasing in magnitude with Peclet number. The cell effective axial diffusivity increases and decreases with Peclet number (for tracers it merely increases). In turbulent flows, gyrotactic effects are weaker, but discernable and manifested as non-zero drift. These results should significantly impact photobioreactor design.
• Coherent imaging of extended objects(0807.1718)

July 10, 2008 physics.optics
When used with coherent light, optical imaging systems, even diffraction-limited, are inherently unable to reproduce both the amplitude and the phase of a two-dimensional field distribution because their impulse response function varies slowly from point to point (a property known as non-isoplanatism). For sufficiently small objects, this usually results in a phase distortion and has no impact on the measured intensity. Here, we show that the intensity distribution can also be dramatically distorted when objects of large extension or of special shapes are imaged. We illustrate the problem using two simple examples: the pinhole camera and the aberration-free thin lens. The effects predicted by our theorical analysis are also confirmed by experimental observations.