BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Luca Brandt

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

ResearcherUnclaimed source record
physics.flu-dynBiological Physicscond-mat.softphysics.comp-phComputational Engineering, Finance, and Sciencecond-mat.mtrl-sciphysics.med-ph

Catalog footprint

What is connected

44works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

44 published item(s)

preprint2023arXiv

Intermittency in turbulent emulsions

We investigate the statistics of turbulence in emulsions of two-immiscible fluids of same density. We compute for the first time velocity increments between points conditioned to be located in the same phase or in different phases and examine their probability density functions (PDF) and the associated structure functions (SF). This enables us to demonstrate that the the presence of the interface reduces the skewness of the PDF at scales below the Kolmogorov-Hinze scale and therefore the magnitude of the energy flux towards the dissipative scales, which is quantified by the third-order SF. The analysis of the higher order SFs shows that multiphase turbulence is more intermittent than single-phase turbulence. In particular, the local scaling exponents of the SFs display a saturation about the Kolmogorov-Hinze scale and below, which indicates the presence of large velocity gradients across the interface. Interestingly, the statistics approach of classic homogeneous isotropic turbulence when significantly increasing the viscosity of the dispersed phase.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

FluTAS: A GPU-accelerated finite difference code for multiphase flows

We present the Fluid Transport Accelerated Solver, FluTAS, a scalable GPU code for multiphase flows with thermal effects. The code solves the incompressible Navier-Stokes equation for two-fluid systems, with a direct FFT-based Poisson solver for the pressure equation. The interface between the two fluids is represented with the Volume of Fluid (VoF) method, which is mass conserving and well suited for complex flows thanks to its capacity of handling topological changes. The energy equation is explicitly solved and coupled with the momentum equation through the Boussinesq approximation. The code is conceived in a modular fashion so that different numerical methods can be used independently, the existing routines can be modified, and new ones can be included in a straightforward and sustainable manner. FluTAS is written in modern Fortran and parallelized using hybrid MPI/OpenMP in the CPU-only version and accelerated with OpenACC directives in the GPU implementation. We present different benchmarks to validate the code, and two large-scale simulations of fundamental interest in turbulent multiphase flows: isothermal emulsions in HIT and two-layer Rayleigh-Bénard convection. FluTAS is distributed through a MIT license and arises from a collaborative effort of several scientists, aiming to become a flexible tool to study complex multiphase flows.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The interaction of droplet dynamics and turbulence cascade

The dynamics of droplet fragmentation in turbulence is described in the Kolmogorov-Hinze framework. Yet, a quantitative theory is lacking at higher concentrations when strong interactions between the phases and coalescence become relevant, which is common in most flows. Here, we address this issue through a fully-coupled numerical study of the droplet dynamics in a turbulent flow at high Reynolds number. By means of time-space spectral statistics, not currently accessible to experiments, we demonstrate that the characteristic scale of the process, the Hinze scale, can be precisely identified as the scale at which the net energy exchange due to capillarity is zero. Droplets larger than this scale preferentially break up absorbing energy from the flow; smaller droplets, instead, undergo rapid oscillations and tend to coalesce releasing energy to the flow. Further, we link the droplet-size-distribution with the probability distribution of the turbulent dissipation. This shows that key in the fragmentation process is the local flux of energy which dominates the process at large scales, vindicating its locality.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Modulation of homogeneous and isotropic turbulence in emulsions

We present a numerical study of emulsions in homogeneous and isotropic turbulence at $Re_λ=137$. The problem is addressed via Direct Numerical Simulations (DNS), where the Volume of Fluid (VOF) is used to represent the complex features of the liquid-liquid interface. We consider a mixture of two iso-density fluids, where fluid properties are varied with the goal of understanding their role in turbulence modulation, in particular the volume fraction ($0.03<α<0.5$), viscosity ratio ($0.01<μ_d/μ_c<100$) and large scale Weber number ($10.6<We_\mathcal{L}<106.5$). The analysis, performed by studying integral quantities and spectral scale-by-scale analysis, reveals that energy is consistently transported from large to small scales by the interface, and no inverse cascade is observed. Furthermore, the total surface is found to be directly proportional to the amount of energy transported, while viscosity and surface tension alter the dynamic that regulates energy transport. We also observe the $-10/3$ and $-3/2$ scaling on droplet size distributions, suggesting that the dimensional arguments which led to their derivation are verified in HIT conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A volume-of-fluid method for interface-resolved simulations of phase-changing two-fluid flows

We present a numerical method for interface-resolved simulations of evaporating two-fluid flows based on the volume-of-fluid (VoF) method. The method has been implemented in an efficient FFT-based two-fluid Navier-Stokes solver, using an algebraic VoF method for the interface representation, and extended with the transport equations of thermal energy and vaporized liquid mass for the single-component evaporating liquid in an inert gas. The conservation of vaporizing liquid and computation of the interfacial mass flux are performed with the aid of a reconstructed signed-distance field, which enables the use of well-established methods for phase change solvers based on level-set methods. The interface velocity is computed with a novel approach that ensures accurate mass conservation, by constructing a divergence-free extension of the liquid velocity field onto the entire domain. The resulting approach does not depend on the type of interface reconstruction (i.e.\ can be employed in both algebraic and geometrical VoF methods). We extensively verified and validated the overall method against several benchmark cases, and demonstrated its excellent mass conservation and good overall performance for simulating evaporating two-fluid flows in two and three dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Aerodynamically-driven rupture of a liquid film by turbulent shear flow

The rupture of a liquid film due to co-flowing turbulent shear flows in the gas phase is studied using a volume-of-fluid method. To simulate this multiphase problem, we use a simplified numerical setup where the liquid film is 'sandwiched' between two fully developed boundary layers from a turbulent channel simulation. The film deforms and eventually ruptures within the shear zone created by the co-flows. This efficient setup allows systematic variation of physical parameters to gauge their role in the aerodynamically-driven deformation and rupture of a liquid film under fully developed sheared turbulence. The present work presents a detailed study of the developing pressure field over the deforming film and related aerodynamic effects, as previously suggested by other authors, in particular the role of the inviscid lift and drag forces. A cumulative lift force is introduced to capture the effect of the alternating pressure minima and maxima forming over the film which amplify and eventually rupture the film. A velocity scale derived from the lift-induced drag force reflects the state of the turbulent boundary layer over the film and collapses the temporal development of this cumulative lift force as well as the amplitude of film deformation with some success for the different film thicknesses and Reynolds numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An interface capturing method for liquid-gas flows at low-Mach number

Multiphase, compressible and viscous flows are of crucial importance in a wide range of scientific and engineering problems. Despite the large effort paid in the last decades to develop accurate and efficient numerical techniques to address this kind of problems, current models need to be further improved to address realistic applications. In this context, we propose a numerical approach to the simulation of multiphase, viscous flows where a compressible and an incompressible phase interact in the low-Mach number regime. In this frame, acoustics is neglected but large density variations of the compressible phase can be accounted for as well as heat transfer, convection and diffusion processes. The problem is addressed in a fully Eulerian framework exploiting a low-Mach number asymptotic expansion of the Navier-Stokes equations. A Volume of Fluid approach (VOF) is used to capture the liquid-gas interface, built on top of a massive parallel solver, second order accurate both in time and space. The second-order-pressure term is treated implicitly and the resulting pressure equation is solved with the eigenexpansion method employing a robust and novel formulation. We provide a detailed and complete description of the theoretical approach together with information about the numerical technique and implementation details. Results of benchmarking tests are provided for five different test cases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effects of surface nanostructure and wettability on pool boiling: A molecular dynamics study

We study the role of surface topology, surface chemistry, and wall superheat temperature on the onset of boiling, bubble nucleation and growth, and the possible formation of an insulating vapour film by means of large-scale MD simulations. In the numerical experiments, we control the system pressure by imposing a constant force on a moving piston. The simulations reveal that the presence of a nanostructure triggers the bubble formation, determines the nucleation site and facilitates the energy transfer from the hot substrate to the water. The surface chemistry, on the other hand, governs the shape of the formed bubble. A hydrophilic surface chemistry accelerates the bubble nucleation, however, decelerates the bubble expansion, thus postponing the formation of the film of vapour. Therefore, a hydrophilic surface provides better energy transfer from the hot wall to the water. By analysing the system energy, we show that irrespective of wall topology and chemistry, there is a wall temperature for which the amount of transferred energy is maximum.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GPU acceleration of CaNS for massively-parallel direct numerical simulations of canonical fluid flows

This work presents the GPU acceleration of the open-source code CaNS for very fast massively-parallel simulations of canonical fluid flows. The distinct feature of the many-CPU Navier-Stokes solver in CaNS is its fast direct solver for the second-order finite-difference Poisson equation, based on the method of eigenfunction expansions. The solver implements all the boundary conditions valid for this type of problems in a unified framework. Here, we extend the solver for GPU-accelerated clusters using CUDA Fortran. The porting makes extensive use of CUF kernels and has been greatly simplified by the unified memory feature of CUDA Fortran, which handles the data migration between host (CPU) and device (GPU) without defining new arrays in the source code. The overall implementation has been validated against benchmark data for turbulent channel flow and its performance assessed on a NVIDIA DGX-2 system (16 Tesla V100 32Gb, connected with NVLink via NVSwitch). The wall-clock time per time step of the GPU-accelerated implementation is impressively small when compared to its CPU implementation on state-of-the-art many-CPU clusters, as long as the domain partitioning is sufficiently small that the data resides mostly on the GPUs. The implementation has been made freely available and open-source under the terms of an MIT license.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Implementation note on a minimal hybrid lubrication/granular dynamics model for dense suspensions

We describe and summarize a class of minimal numerical models emerged from recent development of simulation methods for dense particle suspensions in overdamped linear flows. The main ingredients include (i) a frame-invariant, short-range lubrication model for spherical particles, and (ii) a soft-core, stick/slide frictional contact model activated when particles overlap. We implement a version of the model using a modified velocity-Verlet algorithm that explicitly solves the $N$-body dynamical system in $\mathcal{O}(cN)$ operations, where $c$ is a kernel constant depending on the cutoff of particle interactions. The implementation is validated against literature results on jamming transition and shear thickening suspensions from 40% to 64% volume fractions. Potential strategies to extend the present methodology to non-spherical particles are also suggested for very concentrated suspensions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Low Reynolds number turbulent flows over elastic walls

We study the laminar and turbulent channel flow over a viscous hyper-elastic wall and show that it is possible to sustain an unsteady chaotic turbulent-like flow at any Reynolds number by properly choosing the wall elastic modulus. We propose a physical explanation for this effect by evaluating the shear stress and the turbulent kinetic energy budget in the fluid and elastic layer. We vary the bulk Reynolds number from 2800 to 10 and identify two distinct mechanisms for turbulence production. At moderate and high Reynolds numbers, turbulent fluctuations activate the wall oscillations, which, in turn, amplify the turbulent Reynolds stresses in the fluid. At very low Reynolds number, the only production term is due to the energy input from the elastic wall, which increases with the wall elasticity. This mechanism may be exploited to passively enhance mixing in microfluidic devices.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical simulations of a sphere settling in simple shear flows of yield stress fluids

We perform $3$D numerical simulations to investigate the sedimentation of a single sphere in the absence and presence of a simple cross shear flow in a yield stress fluid with weak inertia. In our simulations, the settling flow is considered to be the primary flow, whereas the linear cross shear flow is a secondary flow with amplitude $10\%$ of the primary flow. To study the effects of elasticity and plasticity of the carrying fluid on the sphere drag as well as the flow dynamics, the fluid is modeled using the elastovisco-plastic (EVP) constitutive laws proposed by \cite{saramito2009new}. The extra non-Newtonian stress tensor is fully coupled with the flow equation and the solid particle is represented by an immersed boundary (IB) method. Our results show that the fore-aft asymmetry in the velocity is less pronounced and the negative wake disappears when a linear cross shear flow is applied. We find that the drag on a sphere settling in a sheared yield stress fluid is reduced significantly as compared to an otherwise quiescent fluid. More importantly, the sphere drag in the presence of a secondary cross shear flow cannot be derived from the pure sedimentation drag law owing to the non-linear coupling between the simple shear flow and the uniform flow. Finally, we show that the drag on the sphere settling in a sheared yield-stress fluid is reduced at higher material elasticity mainly due to the form and viscous drag reduction.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical Simulations of Vorticity Banding of Emulsions in Shear Flows

Multiphase shear flows often show banded structures that affect the global behavior of complex fluids e.g. in microdevices. Here we investigate numerically the banding of emulsions, i.e. the formation of regions of high and low volume fraction, alternated in the vorticity direction and aligned with the flow (shear bands). These bands are associated with a decrease of the effective viscosity of the system. To understand the mechanism of banding experimentally observed we have performed interface resolved simulations of the two-fluid system. The experiments were perfomed starting with a random distribution of droplets which, under the applied shear, evolves in time resulting in a phase separation. To numerically reproduce this process, the banded structures are initialized in a narrow channel confined by two walls moving in opposite direction. We find that the initial banded distribution is stable when droplets are free to merge and unstable when coalescence is prevented. In this case, additionally, the effective viscosity of the system increases, resembling the rheological behavior of suspensions of deformable particles. Droplets coalescence, on the other hand, allows emulsions to reduce the total surface of the system and hence the energy dissipation associated to the deformation, which in turn reduces the effective viscosity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Regimes of heat transfer in particle suspensions

We present results of interface-resolved simulations of heat transfer in suspensions of finite-size neutrally-buoyant spherical particles for solid volume fractions up to 35% and bulk Reynolds numbers from 500 to 5600. An Immersed Boundary-Volume of Fluid method is used to solve the energy equation in the fluid and solid phase. We relate the heat transfer to the regimes of particle motion previously identified, i.e. a viscous regime at low volume fractions and low Reynolds number, particle-laden turbulence at high Reynolds and moderate volume fraction and particulate regime at high volume fractions. We show that in the viscous dominated regime, the heat transfer is mainly due to thermal diffusion with enhancement due to the particle-induced fluctuations. In the turbulent-like regime, we observe the largest enhancement of the global heat transfer, dominated by the turbulent heat flux. In the particulate shear-thickening regime, however, the heat transfer enhancement decreases as mixing is quenched by the particle migration towards the channel core. As a result, a compact loosely-packed core region forms and the contribution of thermal diffusion to the total heat transfer becomes significant once again. The global heat transfer becomes, in these flows at volume fractions larger than 25%, lower than in single-phase turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Weakly adhesive suspension shows rate-dependence in oscillatory but not steady shear flows

We report rheological measurements of a noncolloidal particle suspension in a Newtonian solvent at 40% solid volume fraction. An anomalous, frequency-dependent complex viscosity is found under oscillatory shear (OS) flow, whereas a constant dynamic viscosity is found under the same shear rates in steady shear (SS) flow. We show that this contradiction arises from the underlying microstructural difference between OS and SS, mediated by weak interparticle forces. Discrete element simulations of proxy particle suspensions confirm this hypothesis and reveal an adhesion-induced, shear thinning mechanism with a -1/5 slope, only in OS, in agreement with experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Finite-size spherical particles in a square duct flow of an elastoviscoplastic fluid: an experimental study

The present experimental study addresses the flow of a Yield Stress Fluid with some elasticity (Carbopol gel) in a square duct. The behaviour of two fluids with lower and higher yield stress is investigated at multiple Reynolds numbers $Re^* \in$ (1, 200) and Bingham numbers $Bi \in$ (0.01, 0.35). A secondary flow consisting of eight vortices is observed to recirculate the fluid from the corners to the core. Their extent and intensity grows with increasing $Re^*$. The second objective of this study is to explore the change in flow in the presence of particles. Almost neutrally-buoyant finite-size spherical particles with duct height, $2H$, to particle diameter, $d_p$, ratio of 12 are used at two volume fractions $ϕ$ = 5 and 10\%. Particle Tracking Velocimetry (PTV) is used to measure the velocity of these refractive-index-matched spheres, and PIV to extract the fluid velocity. Simple shadowgraphy is also used for qualitatively visualising the development of the particle distribution along the streamwise direction. The particle distribution pattern changes from being concentrated at the four corners, at low flow rates, to being focussed along a diffused ring between the center and the corners, at high flow rates. The presence of particles induces streamwise and wall-normal velocity fluctuations in the fluid phase; however, the primary Reynolds shear stress is still very small compared to turbulent flows. The size of the plug in the particle-laden cases appears to be smaller than the corresponding single phase cases. Similar to Newtonian fluids, the friction factor increases due to the presence of particles, almost independently of the suspending fluid matrix. Interestingly, predictions based on an increased effective suspension viscosity agrees quite well with the experimental friction factor for the concentrations used in this study.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Effects of the finite particle size in turbulent wall-bounded flows of dense suspensions

We use interface-resolved simulations to study finite-size effects in turbulent channel flow of neutrally-buoyant spheres. Two cases with particle sizes differing by a factor of 2, at the same solid volume fraction of 20% and bulk Reynolds number are considered. These are complemented with two reference single-phase flows: the unladen case, and the flow of a Newtonian fluid with the effective suspension viscosity of the same mixture in the laminar regime. As recently highlighted in Costa et al. (PRL 117, 134501), a particle-wall layer is responsible for deviations of the statistics from what is observed in the continuum limit where the suspension is modeled as a Newtonian fluid with an effective viscosity. Here we investigate the fluid and particle dynamics in this layer and in the bulk. In the particle-wall layer, the near wall inhomogeneity has an influence on the suspension micro-structure over a distance proportional to the particle size. In this layer, particles have a significant (apparent) slip velocity that is reflected in the distribution of wall shear stresses. This is characterized by extreme events (both much higher and much lower than the mean). Based on these observations we provide a scaling for the particle-to-fluid apparent slip velocity as a function of the flow parameters. We also extend the flow scaling laws in to second-order Eulerian statistics in the homogeneous suspension region away from the wall. Finite-size effects in the bulk of the channel become important for larger particles, while negligible for lower-order statistics and smaller particles. Finally, we study the particle dynamics along the wall-normal direction. Our results suggest that 1-point dispersion is dominated by particle-turbulence (and not particle-particle) interactions, while differences in 2-point dispersion and collisional dynamics are consistent with a picture of shear-driven interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Experimental investigation of turbulent suspensions of spherical particles in a square duct

We report experimental observations of turbulent flow with spherical particles in a square duct. Three particle sizes namely: $2H/d_{p}$ = 40, 16 and 9 ($2H$ being the duct full height and $d_{p}$ being the particle diameter) are investigated. The particles are nearly neutrally-buoyant with respect to the suspending fluid. Refractive Index Matched - Particle Image Velocimetry is used for fluid velocity measurement even at the highest particle volume fraction (20%) and Particle Tracking Velocimetry for the particle velocity statistics for the flows seeded with particles of the two largest sizes, whereas only pressure measurements are reported for the smallest particles. Settling effects are seen at the lowest bulk Reynolds number $Re_{2H}\approx$ 10000 whereas, at the highest $Re_{2H}\approx$ 27000, particles are in almost full suspension. The friction factor of the suspensions is found to be significantly larger than that of single-phase duct flow at the lower $Re_{2H}$ investigated; however, the difference decreases when increasing the flow rate and the total drag approaches the values of the single phase flow at the higher Reynolds number considered. The pressure drop is found to decrease with the particle diameter for volume fractions lower than $ϕ$ = 10% for nearly all $Re_{2H}$. However, at the highest volume fraction $ϕ$ = 20%, we report a peculiar non-monotonic behavior: the pressure drop first decreases and then increases with increasing particle size. The decrease of the turbulent drag with particle size at the lowest volume fractions is related to an attenuation of the turbulence. The drag increase for the two largest particle sizes at $ϕ$ = 20%, however, occurs despite this large reduction of the turbulent stresses, and it is therefore due to significant particle-induced stresses.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Turbulent channel flow of a dense binary mixture of rigid particles

We study turbulent channel flow of a binary mixture of finite-size neutrally-buoyant rigid particles by means of interface-resolved direct numerical simulations. We fix the bulk Reynolds number and total solid volume fraction, $Re_b = 5600$ and $Φ=20\%$, and vary the relative fraction of small and large particles. The binary mixture consists of particles of two different sizes, $2h/d_l=20$ and $2h/d_s=30$ where $h$ is the half channel height and $d_l$ and $d_s$ the diameter of the large and small particles. While the particulate flow statistics exhibit a significant alteration of the mean velocity profile and turbulent fluctuations with respect to the unladen flow, the differences between the mono-disperse and bi-disperse cases are small. However, we observe a clear segregation of small particles at the wall in binary mixtures, which affects the dynamics of the near wall region and thus the overall drag. This results in a higher drag in suspensions with a larger amount of large particles. As regards bi-disperse effects on the particle dynamics, a non-monotonic variation of the particle dispersion in the spanwise (homogeneous) direction is observed when increasing the percentage of small/large particles. Finally, we note that particles of the same size tend to cluster more at contact whereas the dynamics of the large particles gives highest collision kernels due to a higher approaching speed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Aspect ratio effect on particle transport in turbulent duct flows

The dynamics of dilute micron-sized spherical inertial particles in turbulent duct flows is studied by means of direct numerical simulations of the carrier phase turbulence with one-way coupled Lagrangian particles. The geometries are a square and a rectangular duct with width-to-height aspect ratio $AR$ of 3 operating at $Re_{τ,c}=360$ . The present study is designed to determine the effect of turbulence-driven secondary motion on the particle dynamics. Our results show that a weak cross-flow secondary motion significantly changes the cross-sectional map of the particle concentration, mean velocity and fluctuations. As the geometry of the duct is widened from $AR=1$ to 3, the secondary vortex on the horizontal wall significantly expands in the spanwise direction, and although the kinetic energy of the secondary flow increases close to the corner, it decays towards the duct centreplane in the $AR=3$ case so as the turbulent carrier phase approaches the behavior in spanwise-periodic channel flows, a fact that significantly affects the particle statistics. In the square duct the particle concentration in the viscous sublayer is maximum at the duct centreplane, whereas the maximum is found closer to the corner, at a distance of $|z/h|\approx1.25$ from the centreplane, in the $AR=3$ case. Interestingly the centreplane concentration in the rectangular duct is around 3 times lower than that in the square duct. Moreover, a second peak in the accumulation distribution is found right at the corners for both ducts. At this location the concentration increases with particle inertia. The secondary motion changes also the cross-stream map of the particle velocities significantly in comparison to the fluid flow statistics. These directly affect the particle velocity fluctuations such that multiple peaks appear near the duct walls for the particle streamwise and wall-normal velocity fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Heat transfer in laminar Couette flow laden with rigid spherical particles

We study heat transfer in plane Couette flow laden with rigid spherical particles by means of direct numerical simulations using a direct-forcing immersed boundary method to account for the dispersed phase. A volume of fluid approach is employed to solve the temperature field inside and outside of the particles. We focus on the variation of the heat transfer with the particle Reynolds number, total volume fraction (number of particles) and the ratio between the particle and fluid thermal diffusivity, quantified in terms of an effective suspension diffusivity. We show that, when inertia at the particle scale is negligible, the heat transfer increases with respect to the unladen case following an empirical correlation recently proposed. In addition, an average composite diffusivity can be used to predict the effective diffusivity of the suspension the inertialess regime when varying the molecular diffusion in the two phases. At finite particle inertia, however, the heat transfer increase is significantly larger, saturating at higher volume fractions. By phase-ensemble averaging we identify the different mechanisms contributing to the total heat transfer and show that the increase of the effective conductivity observed at finite inertia is due to the increase of the transport associated to fluid and particle velocity. We also show that the heat conduction in the solid phase reduces when increasing the particle Reynolds number and so that particles of low thermal diffusivity weakly alter the total heat flux in the suspension at finite particle Reynolds numbers. On the other hand, a higher particle thermal diffusivity significantly increase the total heat transfer.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Motion of an elastic capsule in a constricted microchannel

We study the motion of an elastic capsule through a microchannel characterized by a localized constriction. We consider a capsule with a stress-free spherical shape and impose its steady state configuration in an infinitely long straight channel as the initial condition for our calculations. We report how the capsule deformation, velocity, retention time, and maximum stress of the membrane are affected by the capillary number, Ca, and the constriction shape. We estimate the deformation by measuring the variation of the three-dimensional surface area and a series of alternative quantities easier to extract from experiments. These are the Taylor parameter, the perimeter and the area of the capsule in the spanwise plane. We find that the perimeter is the quantity that reproduces the behavior of the three-dimensional surface area the best. We observe that, in general, area-deformation correlated quantities grow linearly with Ca, while velocity-correlated quantities saturate for large Ca but display a steeper increase for small Ca. The velocity of the capsule divided by the velocity of the flow displays, surprisingly, two different qualitative behaviors for small and large capillary numbers. Finally, we report that longer constrictions and spanwise wall bounded (versus spanwise periodic) domains cause larger deformations and velocities. If the deformation and velocity in the spanwise wall bounded domains are rescaled by the initial equilibrium deformation and velocity, their behavior is undistinguishable from that in a periodic domain. In contrast, a remarkably different behavior is reported in sinusoidally shaped and smoothed rectangular constrictions indicating that the capsule dynamics is particularly sensitive to abrupt changes in the cross section. In a smoothed rectangular constriction larger deformations and velocities occur over a larger distance.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical Study of the Sedimentation of Spheroidal Particles

The gravity-driven motion of rigid particles in a viscous fluid is relevant in many natural and industrial processes, yet this has mainly been investigated for spherical particles. We therefore consider the sedimentation of non-spherical (spheroidal) isolated and particle pairs in a viscous fluid via numerical simulations using the Immersed Boundary Method. The simulations performed here show that the critical Galileo number for the onset of secondary motions decreases as the spheroid aspect ratio departs from 1. Above this critical threshold, oblate particles perform a zigzagging motion whereas prolate particles rotate around the vertical axis while having their broad side facing the falling direction. Instabilities of the vortices in the wake follow when farther increasing the Galileo number. We also study the drafting-kissing-tumbling associated with the settling of particle pairs. We find that the interaction time increases significantly for non-spherical particles and, more interestingly, spheroidal particles are attracted from larger lateral displacements. This has important implications for the estimation of collision kernels and can result in increasing clustering in suspensions of sedimenting spheroids.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Particle transport in turbulent curved pipe flow

Direct numerical simulations (DNS) of particle-laden turbulent flow in straight, mildly curved and strongly bent pipes are performed in which the solid phase is modelled as small heavy spherical particles. A total of seven populations of dilute particles with different Stokes numbers, one-way coupled with their carrier phase, are simulated. The objective is to examine the effect of the curvature on micro-particle transport and accumulation. It is shown that even a slight non-zero curvature in the flow configuration strongly impact the particle concentration map such that the concentration of inertial particles with bulk Stokes number 0.45 (based on bulk velocity and pipe radius) at the inner-bend wall of mildly curved pipe becomes 12.8 times larger than that in the viscous sublayer of the straight pipe. Near-wall helicoidal particle streaks are observed in the curved configurations with their inclination varying with the strength of the secondary motion of the carrier phase. A reflection layer, as previously observed in particle laden turbulent S-shaped channels, is also apparent in the strongly curved pipe with heavy particles. In addition, depending on the curvature, the central regions of the mean Dean vortices appear to be completely depleted of particles, as observed also in the partially re-laminarised region at the inner bend. The turbophoretic drift of the particles is shown to be affected by weak and strong secondary motions of the carrier phase and geometry-induced centrifugal forces. The first and second-order moments of the velocity and acceleration of the particulate phase in the same configurations are addressed in a companion paper by the same authors. The current data-set will be useful for modelling particles advected in wall-bounded turbulent flows where the effects of the curvature are not negligible.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Reduced particle settling speed in turbulence

We study the settling of finite-size rigid spheres in sustained homogeneous isotropic turbulence (HIT) by direct numerical simulations using an immersed boundary method to account for the dispersed solid phase. We study semi-dilute suspensions at different Galileo numbers, Ga. The Galileo number is the ratio between buoyancy and viscous forces, and is here varied via the solid-to-fluid density ratio. The focus is on particles that are slightly heavier than the fluid. We find that in HIT, the mean settling speed is less than that in quiescent fluid; in particular, it reduces by 6%-60% with respect to the terminal velocity of an isolated sphere in quiescent fluid as the ratio between the latter and the turbulent velocity fluctuations is decreased. Analysing the fluid-particle relative motion, we find that the mean settling speed is progressively reduced while reducing the density ratio due to the increase of the vertical drag induced by the particle cross-flow velocity. Unsteady effects contribute to the mean overall drag by about 6%-10%. The probability density functions of particle velocities and accelerations reveal that these are closely related to the features of the turbulent flow. The particle mean-square displacement in the settling direction is found to be similar for all Ga if time is scaled by (2a)/u' (where 2a is the particle diameter and u' is the turbulence velocity root mean square).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effect of particle density in turbulent channel flow laden with finite size particles in semi-dilute conditions

We study the effect of varying the mass and volume fraction of a suspension of rigid spheres dispersed in a turbulent channel flow. We performed several Direct Numerical Simulations using an Immersed Boundary Method for finite-size particles changing the solid to fluid density ratio, the mass fraction and the volume fraction. We find that varying the density ratio between 1 and 10 at constant volume fraction does not alter the flow statistics as much as when varying the volume fraction at constant and at constant mass fraction. Interestingly, the increase in overall drag found when varying the volume fraction is considerably higher than that obtained for increasing density ratios at same volume fraction. The main effect at density ratios of the order of 10 is a strong shear-induced migration towards the centerline of the channel. When the density ratio is further increased up to 1000 the particle dynamics decouple from that of the fluid. The solid phase behaves as a dense gas and the fluid and solid phase statistics drastically change. In this regime, the collision rate is high and dominated by the normal relative velocity among particles.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Universal Scaling Laws for Dense Particle Suspensions in Turbulent Wall-Bounded Flows

The macroscopic behavior of dense suspensions of neutrally-buoyant spheres in turbulent plane channel flow is examined. We show that particles larger than the smallest turbulence scales cause the suspension to deviate from the continuum limit in which its dynamics is well described by an effective suspension viscosity. This deviation is caused by the formation of a particle layer close to the wall with significant slip velocity. By assuming two distinct transport mechanisms in the near-wall layer and the turbulence in the bulk, we define an effective wall location such that the flow in the bulk can still be accurately described by an effective suspension viscosity. We thus propose scaling laws for the mean velocity profile of the suspension flow, together with a master equation able to predict the increase in drag as function of the particle size and volume fraction.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Channel flow of rigid sphere suspensions: particle dynamics in the inertial regime

We consider suspensions of neutrally-buoyant finite-size rigid spherical particles in channel flow and investigate the relation between the particle dynamics and the mean bulk behavior of the mixture for Reynolds numbers $500 \le Re \le 5000$ and particle volume fraction $0\le Φ\le 0.3$, via fully resolved numerical simulations. Analysis of the momentum balance reveals the existence of three different regimes: laminar, turbulent and inertial shear-thickening depending on which of the stress terms, viscous, Reynolds or particle stress, is the major responsible for the momentum transfer across the channel. We show that both Reynolds and particle stress dominated flows fall into the Bagnoldian inertial regime and that the Bagnold number can predict the bulk behavior although this is due to two distinct physical mechanisms. A turbulent flow is characterized by larger particle dispersion and a more uniform particle distribution, whereas the particulate-dominated flows is associated with a significant particle migration towards the channel center where the flow is smooth laminar-like and dispersion low.Interestingly, the collision kernel shows similar values in the different regimes, although the relative particle velocity and clustering clearly vary with inertia and particle concentration.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Continuous growth of droplet size variance due to condensation in turbulent clouds

We use a stochastic model and direct numerical simulation to study the impact of turbulence on cloud droplet growth by condensation. We show that the variance of the droplet size distribution increases in time as t^{1/2}, with growth rate proportional to the large-to-small turbulent scale separation and to the turbulence integral scales but independent of the mean turbulent dissipation. Direct numerical simulations confirm this result and produce realistically broad droplet size spectra over time intervals of 20 minutes, comparable with the time of rain formation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Rheology of confined non-Brownian suspensions

We study the rheology of confined suspensions of neutrally buoyant rigid monodisperse spheres in plane-Couette flow using Direct Numerical Simulations. We find that if the width of the channel is a (small) integer multiple of the sphere's diameter, the spheres self-organize into two-dimensional layers that slide on each other and the suspension's effective viscosity is significantly reduced. Each two-dimensional layer is found to be structurally liquid-like but their dynamics is frozen in time.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sedimentation of finite-size spheres in quiescent and turbulent environments

Sedimentation of a dispersed solid phase is widely encountered in applications and environmental flows, yet little is known about the behavior of finite-size particles in homogeneous isotropic turbulence. To fill this gap, we perform Direct Numerical Simulations of sedimentation in quiescent and turbulent environments using an Immersed Boundary Method to account for the dispersed rigid spherical particles. The solid volume fractions considered are 0.5-1%, while the solid to fluid density ratio 1.02. The particle radius is chosen to be approximately 6 Komlogorov lengthscales. The results show that the mean settling velocity is lower in an already turbulent flow than in a quiescent fluid. The reduction with respect to a single particle in quiescent fluid is about 12\% and 14\% for the two volume fractions investigated. The probability density function of the particle velocity is almost Gaussian in a turbulent flow, whereas it displays large positive tails in quiescent fluid. These tails are associated to the intermittent fast sedimentation of particle pairs in drafting-kissing-tumbling motions. The particle lateral dispersion is higher in a turbulent flow, whereas the vertical one is, surprisingly, of comparable magnitude as a consequence of the highly intermittent behavior observed in the quiescent fluid. Using the concept of mean relative velocity we estimate the mean drag coefficient from empirical formulas and show that non stationary effects, related to vortex shedding, explain the increased reduction in mean settling velocity in a turbulent environment.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The dynamics of a capsule in a wall-bounded oscillating shear flow

The motion of an initially spherical capsule in a wall-bounded oscillating shear flow is investigated via an accelerated boundary integral implementation. The neo-Hookean model is used as the constitutive law of the capsule membrane. The maximum wall-normal migration is observed when the oscillation period of the imposed shear is of the order of the relaxation time of the elastic membrane; hence, the optimal capillary number scales with the inverse of the oscillation frequency and the ratio agrees well with the theoretical prediction in the limit of high-frequency oscillation. The migration velocity decreases monotonically with the frequency of the applied shear and the capsule-wall distance. We report a significant correlation between the capsule lateral migration and the normal stress difference induced in the flow. The periodic variation of the capsule deformation is roughly in phase with that of the migration velocity and normal stress difference, with twice the frequency of the imposed shear. The maximum deformation increases linearly with the membrane elasticity before reaching a plateau at higher capillary numbers when the deformation is limited by the time over which shear is applied in the same direction and not by the membrane deformability. The maximum membrane deformation scales as the distance to the wall to the power 1/3 as observed for capsules and droplets in near-wall steady shear flows.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The motion of a deforming capsule through a corner

A three-dimensional deformable capsule convected through a square duct with a corner is studied via numerical simulations. We develop an accelerated boundary integral implementation adapted to general geometries and boundary conditions. A global spectral method is adopted to resolve the dynamics of the capsule membrane developing elastic tension according to the neo-Hookean constitutive law and bending moments in an inertialess flow. The simulations show that the trajectory of the capsule closely follows the underlying streamlines independently of the capillary number. The membrane deformability, on the other hand, significantly influences the relative area variations, the advection velocity and the principal tensions observed during the capsule motion. The evolution of the capsule velocity displays a loss of the time-reversal symmetry of Stokes flow due to the elasticity of the membrane. The velocity decreases while the capsule is approaching the corner as the background flow does, reaches a minimum at the corner and displays an overshoot past the corner due to the streamwise elongation induced by the flow acceleration in the downstream branch. This velocity overshoot increases with confinement while the maxima of the major principal tension increase linearly with the inverse of the duct width. Finally, the deformation and tension of the capsule are shown to decrease in a curved corner.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Transition and self-sustained turbulence in dilute suspensions of finite-size particles

We study the transition to turbulence of channel flow of finite-size particle suspensions at low volume fraction, i.e. $Φ\approx 0.001$. The critical Reynolds number above which turbulence is sustained reduces to $Re \approx 1675$, in the presence of few particles, independently of the initial condition, a value lower than that of the corresponding single-phase flow, i.e. $Re\approx1775$. In the dilute suspension, the initial arrangement of the particles is important to trigger the transition at a fixed Reynolds number and particle volume fraction. As in single phase flows, streamwise elongated disturbances are initially induced in the flow. If particles can induce oblique disturbances with high enough energy within a certain time, the streaks breakdown, flow experiences the transition to turbulence and the particle trajectories become chaotic. Otherwise, the streaks decay in time and the particles immigrate towards the channel core in a laminar flow.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A microfluidic device to sort capsules by deformability: A numerical study

Guided by extensive numerical simulations, we propose a microfluidic device that can sort elastic capsules by their deformability. The device consists of a duct embedded with a semi-cylindrical obstacle, and a diffuser which further enhances the sorting capability. We demonstrate that the device can operate reasonably well under changes in the initial position of the the capsule. The efficiency of the device remains essentially unaltered under small changes of the obstacle shape (from semi-circular to semi-elliptic cross-section). Confinement along the direction perpendicular to the plane of the device increases its efficiency. This work is the first numerical study of cell sorting by a realistic microfluidic device.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Laminar, turbulent and inertial shear-thickening regimes in channel flow of neutrally buoyant particle suspensions

The aim of this Letter is to characterize the flow regimes of suspensions of finite-size rigid particles in a viscous fluid at finite inertia. We explore the system behavior as function of the particle volume fraction and the Reynolds number (the ratio of flow and particle inertia to viscous forces). Unlike single phase flows where a clear distinction exists between the laminar and the turbulent states, three different regimes can be identified in the presence of a particulate phase, with smooth transitions between them. At low volume fractions, the flow becomes turbulent when increasing the Reynolds number, transitioning from the laminar regime dominated by viscous forces to the turbulent regime characterized by enhanced momentum transport by turbulent eddies. At larger volume fractions, we identify a new regime characterized by an even larger increase of the wall friction. The wall friction increases with the Reynolds number (inertial effects) while the turbulent transport is weakly affected, as in a state of intense inertial shear-thickening. This state may prevent the transition to a fully turbulent regime at arbitrary high speed of the flow.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The effect of the Basset history force on particle clustering in Homogeneous and Isotropic Turbulence

We study the effect of the Basset history force on the dynamics of small particles transported in homogeneous and isotropic turbulence and show that this term, often neglected in previous numerical studies, reduces the small-scale clustering typical of inertial particles. The contribution of this force to the total particle acceleration is, on average, responsible for about 10% of the total acceleration and particularly relevant during rare strong events. At moderate density ratios, i.e. sand or metal powder in water, its presence alters the balance of forces determining the particle acceleration.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The lift-up effect: the linear mechanism behind transition and turbulence in shear flows

The formation and amplification of streamwise velocity perturbations induced by cross-stream disturbances is ubiquitous in shear flows. This disturbance growth mechanism, so neatly identified by Ellingsen and Palm in 1975, is a key process in transition to turbulence and self-sustained turbulence. In this review, we first present the original derivation and early studies and then discuss the non-modal growth of streaks, the result of the lift-up process, in transitional and turbulent shear flows. In the second part, the effects on the lift-up process of additives in the fluid and of a second phase are discussed and new results presented with emphasis on particle-laden shear flows. For all cases considered, we see the lift-up process to be a very robust process, always present as a first step in subcritical transition.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Accumulation of motile elongated micro-organisms in turbulence

We study the effect of turbulence on marine life by performing numerical simulations of motile microorganisms, modelled as prolate spheroids, in isotropic homogeneous turbulence. We show that the clustering and patchiness observed in laminar flows, linear shear and vortex flows, are significantly reduced in a three-dimensional turbulent flow mainly because of the complex topology; elongated micro-orgamisms show some level of clustering in the case of swimmers without any preferential alignment whereas spherical swimmers remain uniformly distributed. Micro-organisms with one preferential swimming direction (e.g. gyrotaxis) still show significant clustering if spherical in shape, whereas prolate swimmers remain more uniformly distributed. Due to their large sensitivity to the local shear, these elongated swimmers react slower to the action of vorticity and gravity and therefore do not have time to accumulate in a turbulent flow. These results show how purely hydrodynamic effects can alter the ecology of microorganisms that can vary their shape and their preferential orientation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Low-Reynolds number swimming in a capillary tube

We use the boundary element method to study the low-Reynolds number locomotion of a spherical model microorganism in a circular tube. The swimmer propels itself by tangen- tial or normal surface motion in a tube whose radius is on the order of the swimmer size. Hydrodynamic interactions with the tube walls significantly affect the average swimming speed and power consumption of the model microorganism. In the case of swimming parallel to the tube axis, the locomotion speed is always reduced (resp. increased) for swimmers with tangential (resp. normal) deformation. In all cases, the rate of work nec- essary for swimming is increased by confinement. Swimmers with no force-dipoles in the far field generally follow helical trajectories, solely induced by hydrodynamic interactions with the tube walls, and in qualitative agreement with recent experimental observations for Paramecium. Swimmers of the puller type always display stable locomotion at a lo- cation which depends on the strength of their force dipoles: swimmers with weak dipoles (small α) swim in the centre of the tube while those with strong dipoles (large α) swim near the walls. In contrast, pusher swimmers and those employing normal deformation are unstable and end up crashing into the walls of the tube. Similar dynamics is observed for swimming into a curved tube. These results could be relevant for the future design of artificial microswimmers in confined geometries.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Micropropulsion and microrheology in complex fluids via symmetry breaking

Many biological fluids have polymeric microstructures and display non-Newtonian rheology. We take advantage of such nonlinear fluid behavior and combine it with geometrical symmetry-breaking to design a novel small-scale propeller able to move only in complex fluids. Its propulsion characteristics are explored numerically in an Oldroyd-B fluid for finite Deborah numbers while the small Deborah number limit is investigated analytically using a second-order fluid model. We then derive expressions relating the propulsion speed to the rheological properties of the complex fluid, allowing thus to infer the normal stress coefficients in the fluid from the locomotion of the propeller. Our simple mechanism can therefore be used either as a non-Newtonian micro-propeller or as a micro-rheometer.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Shear thickening in non-Brownian suspensions: an excluded volume effect

Shear thickening appears as an increase of the viscosity of a dense suspension with the shear rate, sometimes sudden and violent at high volume fraction. Its origin for noncolloidal suspension with non-negligible inertial effects is still debated. Here we consider a simple shear flow and demonstrate that fluid inertia causes a strong microstructure anisotropy that results in the formation of a shadow region with no relative flux of particles. We show that shear thickening at finite inertia can be explained as an increase of the effective volume fraction when considering the dynamically excluded volume due to these shadow regions

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Locomotion by tangential deformation in a polymeric fluid

In several biologically relevant situations, cell locomotion occurs in polymeric fluids with Weissenberg {number} larger than one. Here we present results of three-dimensional numerical simulations for the steady locomotion of a self-propelled body in a model polymeric (Giesekus) fluid at low Reynolds number. Locomotion is driven by steady tangential deformation at the surface of the body (so-called squirming motion). In the case of a spherical squirmer, we show that the swimming velocity is systematically less than that in a Newtonian fluid, with a minimum occurring for Weissenberg numbers of order one. The rate of work done by the swimmer always goes up compared to that occurring in the Newtonian solvent alone, but is always lower than the power necessary to swim in a Newtonian fluid with the same viscosity. The swimming efficiency, defined as the ratio between the rate of work necessary to pull the body at the swimming speed in the same fluid and the rate of work done by swimming, is found to always be increased in a polymeric fluid. Further analysis reveals that polymeric stresses break the Newtonian front-back symmetry in the flow profile around the body. In particular, a strong negative elastic wake is present behind the swimmer, which correlates with strong polymer stretching, and its intensity increases with Weissenberg number and viscosity contrasts. {The velocity induced by the squirmer is found to decay in space faster than} in a Newtonian flow, with a strong {dependence} on the polymer relaxation time and viscosity. Our computational results are also extended to prolate spheroidal swimmers and smaller polymer stretching are obtained for slender shapes compared to bluff swimmers. The swimmer with an aspect ratio of two is found to be the most hydrodynamically efficient.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Statistics of polymer extensions in turbulent channel flow

We present direct numerical simulations of turbulent channel flow with passive Lagrangian polymers. To understand the polymer behavior we investigate the behavior of infinitesimal line elements and calculate the probability distribution function (PDF) of finite-time Lyapunov exponents and from them the corresponding Cramer's function for the channel flow. We study the statistics of polymer elongation for both the Oldroyd-B model (for Weissenberg number $\Wi <1$) and the FENE model. We use the location of the minima of the Cramer's function to define the Weissenberg number precisely such that we observe coil-stretch transition at $\Wi\approx1$. We find agreement with earlier analytical predictions for PDF of polymer extensions made by Balkovsky, Fouxon and Lebedev [Phys. Rev. Lett., 84, 4765 (2000).] for linear polymers (Oldroyd-B model) with $\Wi<1$ and by Chertkov [Phys. Rev. Lett., 84, 4761 (2000).] for nonlinear FENE-P model of polymers. For $\Wi>1$ (FENE model) the polymer are significantly more stretched near the wall than at the center of the channel where the flow is closer to homogenous isotropic turbulence. Furthermore near the wall the polymers show a strong tendency to orient along the stream-wise direction of the flow but near the centerline the statistics of orientation of the polymers is consistent with analogous results obtained recently in homogeneous and isotropic flows.

0 citations0 reviewsNo code linkedOpen access