Numerisk studie av partikelsuspensioner i Newtonska och icke-Newtonska vätskor

Solid or deformable particles suspended in a viscous fluid are of scientific and technological interest in a broad range of applications. Pyroclastic flows from volcanoes, sedimentation flows in river bed, food industries, oil-well drilling, as well as blood flow in the human body and the motion of suspended micro-organisms in water (like plankton) are among the possible examples. Often, in these particulate flows, the carrier fluid might exhibit an inelastic or a visco-elastic non-Newtonian behavior. Understanding the behavior of these suspensions is a very difficult task. Indeed, the complexities and challenges of multiphase flows are mainly due to the large number of governing parameters such as the physical properties of the particles (e.g., shape, size, stiffness, density difference with suspended fluid, solid volume fraction), the large set of interactions among particles and the properties of the carrier fluid (Newtonian or non-Newtonian); variations of each of these parameters may provide substantial quantitative and qualitative changes in the behavior of the suspension and affect the overall dynamics in several and sometimes surprising ways. The aim of this work is therefore to provide a deeper understanding of the behavior of particle suspensions in laminar Newtonian and non-Newtonian (inelastic and/or visco-elastic) fluid flows for a wide range of different parameters. To this purpose, particle-resolved direct numerical simulations of spherical particles are performed, using an efficient and accurate numerical tool. The code is based on the Immersed Boundary Method (IBM) for the fluid-solid interactions with lubrication, friction and collision models for the close range particle-particle (particle-wall) interactions. Both inelastic (Carreau and power-law), and visco-elastic models (Oldroyd-B and Giesekus) are employed to investigate separately the shear-thinning, shear-thickening, viscoelastic and combined shear-thinning visco-elastic features of the most commonly encountered non-Newtonian fluids. Moreover, a fully Eulerian numerical algorithm based on the one-continuum formulation is used to examine the case of an hyper-elastic neo-Hookean deformable particle suspended in a Newtonian flows. Firstly, we have investigated suspensions of solid spheres in Newtonian, shear thinning and shear thickening fluids in the simple shear flow created by two walls moving in opposite directions, considering various solid volume fractions and particle Reynolds numbers, thus including inertial effects. The results show that that the non-dimensional relative viscosity of of the suspension and the mean value of the local shear-rate can be well predicted by homogenization theory, more accurately for lower particle concentrations. Moreover, we show that in the presence of inertia, the effective viscosity of these suspensions deviates from that of Stokesian suspensions. We also examine the role of fluid elasticity, shear-thinning and combined shear-thinning visco-elastic effects on the simple linear Couette shear flow of neutrally-buoyant rigid spherical particles. It is found that the effective viscosity grows monotonically with the solid volume fraction and that all the Non-Newtonian cases exhibit a lower effective viscosity than the Newtonian ones; in addition, we show that elastic effects dominate at low elasticity whereas shear thinning is predominant at high applied shear rates. These variations in the effective viscosity are mainly due to changes in the particle-induced shear stress component. We then study the settling of spherical particles in quiescent wall-bounded Newtonian and shear-thinning fluids at three different solid volume fractions. We find that the mean settling velocities decrease with the particle concentration as a consequence of the hindering effect and thatthe mean settling speed is always larger in the shear thinning fluid than in the Newtonian one, due to the reduction of the local fluid viscosity around the particles which leads to a lower drag force acting on the particles. Finally, the inertial migration of hyper-elastic deformable particle in laminar pipe flows is also investigated. We consider different flow rates and various levels of particle elasticity. We observe that the particle deforms and experiences a lateral movement while traveling downstream through the pipe, always finding a stable position at the pipe centerline.Suspensioner av solida eller deformerbara partiklar iviskösa vätskor är av vetenskapligt och teknologiskt intresse för ett stortspann av applikationer. Några typiska exempel inkluderar pyroklastiskaflöden från vulkaner, sedimenterande flöden i flodbäddar,livsmedelsindustrin, oljebrunnsborrning, blodflödet i människokroppen samtrörelsen hos mikroorganismer (till exempel plankton) i havet. I dessapartikelflöden kan den bärande vätskan ha ett icke-elastiskt ellerviskoelastiskt icke-Newtonskt beteende. Att förstå beteendet hos dessasuspensioner är en mycket svår uppgift. Komplexiteten hos, och utmaningenmed, multifasflöden beror till största delen på det stora antal styrandeparametrar. Dessa inkluderar de fysikaliska partikelegenskaperna (tillexempel form, storlek, styvhet, densitetsskillnad mot det bärande medietsamt volymfraktion), den stora mängden interaktioner mellan partiklarnasamt egenskaperna hos den bärande fluiden (Newtonsk eller icke-Newtonsk).Variationer i vardera av dessa parametrar kan leda till stora kvantitativaoch kvalitativa förändringar i suspensionens beteende och kan påverka denövergripande dynamiken på många, ibland överraskande, sätt. Målet meddenna avhandling är därför att ge en djupare förståelse avpartikelsuspensioner i laminära, Newtonska och icke-Newtonska(icke-elastiska och/eller visko-elasiska), flöden för ett stort spann avparametrar. För detta används ett effektivt och precist simuleringsverktygsom tillåter partikelupplösta, numeriska simuleringar av sfäriskapartiklar. Koden är baserad på Immersed boundary-metodiken (IBM) förfluid-strukturinteraktion med lubrikations-, friktions- ochkollisionsmodeller för partikel-partikel- och partikel-vägginteraktioner.Både icke-elastiska (Carreau och power-law) och viskoelastiska (Oldroyd-Boch Giesekus) modeller betraktades för att, i isolering, undersökaeffekterna av skjuvförtunnande, skjuvförtjockande, viskoelasticitet samtkombinationen av skjuvförtunning och viskoelastik, vilka vanligen förekommerhos icke-Newtonska fluider. Därutöver användes en Eulerisk numeriskalgoritm baserad på en en-kontinuumformulering för att undersöka fallet meden hyperelastisk, neo-Hookisk och deformerbar partikel i en Newtonsk vätska. Till att börja med undersöks suspensioner av solida sfärer i Netwonska,skjuvförtunnande samt skjuvförtjockande fluider i ett skjuvflöde genereratmellan två väggar som rör sig i motsatt riktning. Varierandevolymfraktioner (av partiklar) och partikel-Reynoldstal (dvs inkluderandeav fluidtröghet) betraktas. Resultaten visar att den dimensionslösarelativa viskositeten hos suspensionen och medelvärdet av den lokalaskjuvhastigheten kan väl förutsägas av homogeniseringsteori, speciellttillförlitligt vid låga partikelkoncentrationer. Därutöver visas att deneffektiva viskositeten hos dessa suspensioner avviker från suspensioner iStokesflöde när flödeströghet inkluderas. Därutöver undersöktes rollen hos elasticitet, skjuvförtunnande samtkombinerad skjuvförtunnande och viskoelasticitet i det bärande mediet påett linjärt Couetteflöde med densitetsmatchade, rigida och sfäriskapartiklar. Den effektiva viskositeten växer monotont medpartikelvolymfraktionen och alla icke-Newtonska fall uppvisar en lägreeffektiv viskositet än de motsvarande Newtonska fallen. Det visas även attelastiska effekter dominerar vid låg elasticitet medan skjuvförtunnandeeffekter dominerar vid höga skjuvhastigheter. Dessa variationer i effektivviskositet beror främst på förändringar i den partikelinduceradekomponenten av skjuvspänningen. Efter detta studeras sedimentering av sfäriska partiklar i ettstillastående flöde mellan två väggar. Både Newtonska och skjuvförtunnandevätskor betraktas vid tre olika partikelvolymfraktioner. Det visas attmedelvärdet av sedimenteringshastigheten minskar med partikelkoncentrationpå grund av den hindrande effekten av omgivande partiklar. Därutöver ärmedelsedimentationshastigheten alltid större i en skjuvförtunnande än enNewtonsk vätska på grund av reduktionen i lokal fluidviskositet runtpartiklarna, vilket leder till en lägre motståndskraft. Slutligen undersöks även tröghetsinducerad migration av hyperelastiska ochdeformerbara partiklar i ett laminärt rörflöde. Olika flöden och nivåer avelasticitet hos partikeln betraktas. Partikeldeformation och lateralrörelse observeras för partiklarna när de rör sig nedströms längs röret,vilket leder till att partiklarna alltid finner en stabil position vidrörets centerlinje.QC 20191114