scholarly journals Interplay of inertia and deformability on rheological properties of a suspension of capsules

2014 ◽  
Vol 751 ◽  
pp. 725-745 ◽  
Author(s):  
Timm Krüger ◽  
Badr Kaoui ◽  
Jens Harting
Keyword(s):  
Three Dimensional ◽  
Reynolds Numbers ◽  
Particle Collision ◽  
Immersed Boundary ◽  
Lateral Migration ◽  
Particle Deformation ◽  
Substantial Impact ◽  
Inertia Effects ◽  
Focusing Effect ◽  
Soft Particles

AbstractThe interplay of inertia and deformability has a substantial impact on the transport of soft particles suspended in a fluid. However, to date a thorough understanding of these systems is still missing, and only a limited number of experimental and theoretical studies are available. We combine the finite-element, immersed-boundary and lattice-Boltzmann methods to simulate three-dimensional suspensions of soft particles subjected to planar Poiseuille flow at finite Reynolds numbers. Our findings confirm that the particle deformation and inclination increase when inertia is present. We observe that the Segré–Silberberg effect is suppressed with respect to the particle deformability. Depending on the deformability and strength of inertial effects, inward or outward lateral migration of the particles takes place. In particular, for increasing Reynolds numbers and strongly deformable particles, a hitherto unreported distinct flow focusing effect emerges, which is accompanied by a non-monotonic behaviour of the apparent suspension viscosity and thickness of the particle-free layer close to the channel walls. This effect can be explained by the behaviour of a single particle and the change of the particle collision mechanism when both deformability and inertia effects are relevant.

A numerical study of the motion of drops in Poiseuille flow. Part 1. Lateral migration of one drop

2000 ◽  
Vol 411 ◽  
pp. 325-350 ◽  
Author(s):  
SAEED MORTAZAVI ◽  
GRÉTAR TRYGGVASON
Keyword(s):  
Reynolds Number ◽  
Steady State ◽  
Poiseuille Flow ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Lateral Position ◽  
Physical Parameters ◽  
Computational Domain ◽  
Two Dimensional ◽  
Lateral Migration

The cross-stream migration of a deformable drop in two-dimensional Hagen–Poiseuille flow at finite Reynolds numbers is studied numerically. In the limit of a small Reynolds number (< 1), the motion of the drop depends strongly on the ratio of the viscosity of the drop fluid to the viscosity of the suspending fluid. For viscosity ratio 0.125 a drop moves toward the centre of the channel, while for ratio 1.0 it moves away from the centre until halted by wall repulsion. The rate of migration increases with the deformability of the drop. At higher Reynolds numbers (5–50), the drop either moves to an equilibrium lateral position about halfway between the centreline and the wall – according to the so-called Segre–Silberberg effect or it undergoes oscillatory motion. The steady-state position depends only weakly on the various physical parameters of the flow, but the length of the transient oscillations increases as the Reynolds number is raised, or the density of the drop is increased, or the viscosity of the drop is decreased. Once the Reynolds number is high enough, the oscillations appear to persist forever and no steady state is observed. The numerical results are in good agreement with experimental observations, especially for drops that reach a steady-state lateral position. Most of the simulations assume that the flow is two-dimensional. A few simulations of three-dimensional flows for a modest Reynolds number (Re = 10), and a small computational domain, confirm the behaviour seen in two dimensions. The equilibrium position of the three-dimensional drop is close to that predicted in the simulations of two-dimensional flow.

scholarly journals Simulation of Fluid and Structure Interface with Immersed Boundary–Lattice Boltzmann Method Involving Turbulence Models

2018 ◽  
Vol 2018 ◽  
pp. 1-12
Author(s):  
Zhikai Wang ◽  
Xiongliang Yao ◽  
Nana Yang ◽  
Zhenhuan Xu
Keyword(s):  
Coordinate System ◽  
Turbulent Flows ◽  
Lattice Boltzmann ◽  
Three Dimensional ◽  
Turbulence Models ◽  
Slip Boundary Condition ◽  
Reynolds Numbers ◽  
Immersed Boundary ◽  
Square Cylinder ◽  
Slip Boundary

The multiple-relaxation-time (MRT) version of the immersed boundary–lattice Boltzmann (IB-LB) method is developed to simulate fluid-structure interfaces. The innovations include the implicit velocity correction to ensure no-slip boundary conditions and the incorporated Smagorinsky’s algebraic eddy viscosity for simulating turbulent flows. Both straight and curved interfaces are investigated. The streamlines penetration can be well prevented, which means the no-slip boundary condition can be guaranteed. Due to the existence of two coordinate systems: the Lagrangian coordinate system and the Eulerian coordinate system, the velocity and force properties on the structure can be easily calculated. Several benchmark simulation cases are carried out to verify the correctness of the model, including flow around circular cylinder at Re = 20, 150, and 3900 and flow around square cylinder at Re = 150 and 1000. The results agree well with previous studies, especially in the events of lower Reynolds numbers. Due to the three-dimensional turbulence vortex effects, the discrepancy increases are associated with higher Reynolds numbers. In addition, the effect of rotating velocity on the interaction process of the square cylinder in flows is researched, coupled with the capability of dealing with the moving boundaries.

3D simulation of emulsion flow using the boundary element method on the heterogeneous computational systems

2012 ◽  
Vol 9 (1) ◽  
pp. 142-146
Author(s):  
O.A. Solnyshkina
Keyword(s):  
Boundary Element Method ◽  
Boundary Element ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Problem Size ◽  
Low Reynolds Numbers ◽  
Infinite Domains ◽  
The Matrix ◽  
Computational Systems ◽  
Element Method

In this work the 3D dynamics of two immiscible liquids in unbounded domain at low Reynolds numbers is considered. The numerical method is based on the boundary element method, which is very efficient for simulation of the three-dimensional problems in infinite domains. To accelerate calculations and increase the problem size, a heterogeneous approach to parallelization of the computations on the central (CPU) and graphics (GPU) processors is applied. To accelerate the iterative solver (GMRES) and overcome the limitations associated with the size of the memory of the computation system, the software component of the matrix-vector product

scholarly journals Analogue tuning of particle focusing in elasto-inertial flow

Meccanica ◽  
2021 ◽  
Author(s):  
I. Banerjee ◽  
M. E. Rosti ◽  
T. Kumar ◽  
L. Brandt ◽  
A. Russom
Keyword(s):  
Numerical Simulations ◽  
Cross Section ◽  
Immersed Boundary Method ◽  
Circular Cross Section ◽  
Finite Size ◽  
Reynolds Numbers ◽  
Immersed Boundary ◽  
Particle Focusing ◽  
Inertial Flow ◽  
Particle Behaviour

AbstractWe report a unique tuneable analogue trend in particle focusing in the laminar and weak viscoelastic regime of elasto-inertial flows. We observe experimentally that particles in circular cross-section microchannels can be tuned to any focusing bandwidths that lie between the “Segre-Silberberg annulus” and the centre of a circular microcapillary. We use direct numerical simulations to investigate this phenomenon and to understand how minute amounts of elasticity affect the focussing of particles at increasing flow rates. An Immersed Boundary Method is used to account for the presence of the particles and a FENE-P model is used to simulate the presence of polymers in a Non-Newtonian fluid. The numerical simulations study the dynamics and stability of finite size particles and are further used to analyse the particle behaviour at Reynolds numbers higher than what is allowed by the experimental setup. In particular, we are able to report the entire migration trajectories of the particles as they reach their final focussing positions and extend our predictions to other geometries such as the square cross section. We believe complex effects originate due to a combination of inertia and elasticity in the weakly viscoelastic regime, where neither inertia nor elasticity are able to mask each other’s effect completely, leading to a number of intermediate focusing positions. The present study provides a fundamental new understanding of particle focusing in weakly elastic and strongly inertial flows, whose findings can be exploited for potentially multiple microfluidics-based biological sorting applications.

A three-dimensional immersed boundary method based on an algebraic forcing-point-searching scheme for water impact problems

2021 ◽  
Vol 233 ◽  
pp. 109189
Author(s):  
Bin Yan ◽  
Wei Bai ◽  
Sheng-Chao Jiang ◽  
Peiwen Cong ◽  
Dezhi Ning ◽  
...  
Keyword(s):  
Immersed Boundary Method ◽  
Three Dimensional ◽  
Immersed Boundary ◽  
Water Impact ◽  
Boundary Method ◽  
Impact Problems

Turbulent three-dimensional dielectric electrohydrodynamic convection between two plates

2012 ◽  
Vol 696 ◽  
pp. 228-262 ◽  
Author(s):  
A. Kourmatzis ◽  
J. S. Shrimpton
Keyword(s):  
Spatial Data ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Two Dimensions ◽  
Three Dimensions ◽  
Energy Budgets ◽  
Turbulent Regime ◽  
Scalar Flux ◽  
Wide Range ◽  
Scalar Variance

AbstractThe fundamental mechanisms responsible for the creation of electrohydrodynamically driven roll structures in free electroconvection between two plates are analysed with reference to traditional Rayleigh–Bénard convection (RBC). Previously available knowledge limited to two dimensions is extended to three-dimensions, and a wide range of electric Reynolds numbers is analysed, extending into a fully inherently three-dimensional turbulent regime. Results reveal that structures appearing in three-dimensional electrohydrodynamics (EHD) are similar to those observed for RBC, and while two-dimensional EHD results bear some similarities with the three-dimensional results there are distinct differences. Analysis of two-point correlations and integral length scales show that full three-dimensional electroconvection is more chaotic than in two dimensions and this is also noted by qualitatively observing the roll structures that arise for both low (${\mathit{Re}}_{E} = 1$) and high electric Reynolds numbers (up to ${\mathit{Re}}_{E} = 120$). Furthermore, calculations of mean profiles and second-order moments along with energy budgets and spectra have examined the validity of neglecting the fluctuating electric field ${ E}_{i}^{\ensuremath{\prime} } $ in the Reynolds-averaged EHD equations and provide insight into the generation and transport mechanisms of turbulent EHD. Spectral and spatial data clearly indicate how fluctuating energy is transferred from electrical to hydrodynamic forms, on moving through the domain away from the charging electrode. It is shown that ${ E}_{i}^{\ensuremath{\prime} } $ is not negligible close to the walls and terms acting as sources and sinks in the turbulent kinetic energy, turbulent scalar flux and turbulent scalar variance equations are examined. Profiles of hydrodynamic terms in the budgets resemble those in the literature for RBC; however there are terms specific to EHD that are significant, indicating that the transfer of energy in EHD is also attributed to further electrodynamic terms and a strong coupling exists between the charge flux and variance, due to the ionic drift term.

scholarly journals Steady flow separation patterns in a 45 degree junction

2000 ◽  
Vol 411 ◽  
pp. 1-38 ◽  
Author(s):  
C. ROSS ETHIER ◽  
SUJATA PRAKASH ◽  
DAVID A. STEINMAN ◽  
RICHARD L. LEASK ◽  
GREGORY G. COUCH ◽  
...  
Keyword(s):  
Flow Separation ◽  
Stokes Equations ◽  
Bypass Graft ◽  
Three Dimensional ◽  
Symmetry Plane ◽  
Axial Flow ◽  
Reynolds Numbers ◽  
Secondary Flows ◽  
Measured Velocity ◽  
Curved Tube

Numerical and experimental techniques were used to study the physics of flow separation for steady internal flow in a 45° junction geometry, such as that observed between two pipes or between the downstream end of a bypass graft and an artery. The three-dimensional Navier–Stokes equations were solved using a validated finite element code, and complementary experiments were performed using the photochromic dye tracer technique. Inlet Reynolds numbers in the range 250 to 1650 were considered. An adaptive mesh refinement approach was adopted to ensure grid-independent solutions. Good agreement was observed between the numerical results and the experimentally measured velocity fields; however, the wall shear stress agreement was less satisfactory. Just distal to the ‘toe’ of the junction, axial flow separation was observed for all Reynolds numbers greater than 250. Further downstream (approximately 1.3 diameters from the toe), the axial flow again separated for Re [ges ] 450. The location and structure of axial flow separation in this geometry is controlled by secondary flows, which at sufficiently high Re create free stagnation points on the model symmetry plane. In fact, separation in this flow is best explained by a secondary flow boundary layer collision model, analogous to that proposed for flow in the entry region of a curved tube. Novel features of this flow include axial flow separation at modest Re (as compared to flow in a curved tube, where separation occurs only at much higher Re), and the existence and interaction of two distinct three-dimensional separation zones.

scholarly journals Fish larvae exploit edge vortices along their dorsal and ventral fin folds to propel themselves

2016 ◽  
Vol 13 (116) ◽  
pp. 20160068 ◽  
Author(s):  
Gen Li ◽  
Ulrike K. Müller ◽  
Johan L. van Leeuwen ◽  
Hao Liu
Keyword(s):  
Larval Fish ◽  
Fish Larvae ◽  
Fluid Dynamic ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Bony Fish ◽  
The Body ◽  
Functional Explanation ◽  
Image Velocimetry ◽  
Median Fins

Larvae of bony fish swim in the intermediate Reynolds number ( Re ) regime, using body- and caudal-fin undulation to propel themselves. They share a median fin fold that transforms into separate median fins as they grow into juveniles. The fin fold was suggested to be an adaption for locomotion in the intermediate Reynolds regime, but its fluid-dynamic role is still enigmatic. Using three-dimensional fluid-dynamic computations, we quantified the swimming trajectory from body-shape changes during cyclic swimming of larval fish. We predicted unsteady vortices around the upper and lower edges of the fin fold, and identified similar vortices around real larvae with particle image velocimetry. We show that thrust contributions on the body peak adjacent to the upper and lower edges of the fin fold where large left–right pressure differences occur in concert with the periodical generation and shedding of edge vortices. The fin fold enhances effective flow separation and drag-based thrust. Along the body, net thrust is generated in multiple zones posterior to the centre of mass. Counterfactual simulations exploring the effect of having a fin fold across a range of Reynolds numbers show that the fin fold helps larvae achieve high swimming speeds, yet requires high power. We conclude that propulsion in larval fish partly relies on unsteady high-intensity vortices along the upper and lower edges of the fin fold, providing a functional explanation for the omnipresence of the fin fold in bony-fish larvae.

Turbulent flow between counter-rotating concentric cylinders: a direct numerical simulation study

2008 ◽  
Vol 615 ◽  
pp. 371-399 ◽  
Author(s):  
S. DONG
Keyword(s):  
Turbulent Flow ◽  
Large Scale ◽  
Outer Cylinder ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Taylor Vortex ◽  
Small Scale ◽  
Mean Flow ◽  
Concentric Cylinders ◽  
High Reynolds

We report three-dimensional direct numerical simulations of the turbulent flow between counter-rotating concentric cylinders with a radius ratio 0.5. The inner- and outer-cylinder Reynolds numbers have the same magnitude, which ranges from 500 to 4000 in the simulations. We show that with the increase of Reynolds number, the prevailing structures in the flow are azimuthal vortices with scales much smaller than the cylinder gap. At high Reynolds numbers, while the instantaneous small-scale vortices permeate the entire domain, the large-scale Taylor vortex motions manifested by the time-averaged field do not penetrate a layer of fluid near the outer cylinder. Comparisons between the standard Taylor–Couette system (rotating inner cylinder, fixed outer cylinder) and the counter-rotating system demonstrate the profound effects of the Coriolis force on the mean flow and other statistical quantities. The dynamical and statistical features of the flow have been investigated in detail.

Sign in / Sign up

Export Citation Format

Share Document