
We investigate the elastoviscoplastic flow through porous media by numerical
simulations. We solve the NavierStokes 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 OldroydB 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 timedependent 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.

We study the inertial migration of finitesize neutrally buoyant spherical
particles in dilute and semidilute 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 crossstream motions in the duct crosssection,
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 crossstream
velocities.

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 directforcing immersed boundary method to account
for the dispersed phase, combined with a softsphere collision model and
lubrication corrections for shortrange particleparticle and particlewall
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.

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 encounterbased processes in the plankton including
preypredator 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
downwelling 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.

We report the large polar magnetooptical Kerr effect in L10MnxGa 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
spinorbit 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
compositiondependence that varies with the photon energy. These L10MnxGa
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.

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 neutrallybuoyant spheres at relatively high volume fractions in a plane
channel. Direct Numerical Simulation are performed in the range of volume
fractions Phi=00.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 loglaw. 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 particleinduced 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.

Shear flow significantly affects the transport of swimming algae in
suspension. For example, viscous and gravitational torques bias bottomheavy
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). Timeresolved
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
generalizedTaylordispersion. 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 nonzero drift. These results should significantly impact
photobioreactor design.

When used with coherent light, optical imaging systems, even
diffractionlimited, are inherently unable to reproduce both the amplitude and
the phase of a twodimensional field distribution because their impulse
response function varies slowly from point to point (a property known as
nonisoplanatism). 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 aberrationfree thin
lens. The effects predicted by our theorical analysis are also confirmed by
experimental observations.