We present a particle-scale simulation method for non-Brownian suspension flows at low Reynolds number, based on the Fictitious Domain Method and supplemented by sub-grid lubrication corrections. In their usual form, sub-grid corrections involve the underlying linear flow. In the present work, the conditions required to avoid considering this ambient flow while keeping frame indifference are determined, and a sub-grid correction matrix is built for the particle-particle and particle-wall hydrodynamic interactions such that the mentioned conditions inherently hold. This procedure extends the correct use of such sub-grid corrections to nonlinear flows. The method is validated against various particle-flow configurations, involving a few particles in linear and nonlinear flows with various boundary conditions, or concentrated suspensions in simple shear flow. Then, we study the role of adhesive forces in frictional suspensions, by varying both the intensity of adhesion and the volume fraction: after presenting the difficulties of simulating adhesive suspensions at low volume fractions due to particle depletion at the wall and shear-banding, we show that the relative viscosity of the suspension is a function of both volume fraction and shear stress. The variation of the viscosity with these two parameters may be understood in the usual frame of suspension jamming, provided that the jamming volume fraction now depends on the dimensionless suspension stress that involves the adhesive force between particles. The variation of the jamming volume fraction with stress may be interpreted as the variation of the yield stress with volume fraction. This curve separates the plane (volume fraction, shear stress) into two regions: one in which the suspension flows and the other in which it jams. We continue by investigating the shear-induced particle migration in a pressure-driven channel suspension flow. In such a system, particles are driven toward the channel center, resulting in a volume fraction gradient across the channel and a plug region where the jamming volume fraction as measured in shear flow may be exceeded. The flow profile is thus highly nonlinear, both due to the pressure gradient and the time-varying concentration gradient, which justifies the use of the present numerical method. We compare the results to a modified version of the Suspension Balance Model, confirming the well-known limits of such a model in the plug region. Finally, we show that the computed stresses obey the usual momentum balance in the channel: the tangential stress is indeed driven by the pressure gradient as predicted by theory and the second normal stress (parallel to the velocity gradient) does not vary over the channel width, although contact and hydrodynamic contributions significantly do over both time and space. This confirms the good behavior of the present numerical method when nonlinear flows are tackled.
Simulation of concentrated non-Brownian frictional and adhesive suspensions in linear and nonlinear flows / Orsi, Michel. - (2022).
Simulation of concentrated non-Brownian frictional and adhesive suspensions in linear and nonlinear flows
Michel Orsi
2022
Abstract
We present a particle-scale simulation method for non-Brownian suspension flows at low Reynolds number, based on the Fictitious Domain Method and supplemented by sub-grid lubrication corrections. In their usual form, sub-grid corrections involve the underlying linear flow. In the present work, the conditions required to avoid considering this ambient flow while keeping frame indifference are determined, and a sub-grid correction matrix is built for the particle-particle and particle-wall hydrodynamic interactions such that the mentioned conditions inherently hold. This procedure extends the correct use of such sub-grid corrections to nonlinear flows. The method is validated against various particle-flow configurations, involving a few particles in linear and nonlinear flows with various boundary conditions, or concentrated suspensions in simple shear flow. Then, we study the role of adhesive forces in frictional suspensions, by varying both the intensity of adhesion and the volume fraction: after presenting the difficulties of simulating adhesive suspensions at low volume fractions due to particle depletion at the wall and shear-banding, we show that the relative viscosity of the suspension is a function of both volume fraction and shear stress. The variation of the viscosity with these two parameters may be understood in the usual frame of suspension jamming, provided that the jamming volume fraction now depends on the dimensionless suspension stress that involves the adhesive force between particles. The variation of the jamming volume fraction with stress may be interpreted as the variation of the yield stress with volume fraction. This curve separates the plane (volume fraction, shear stress) into two regions: one in which the suspension flows and the other in which it jams. We continue by investigating the shear-induced particle migration in a pressure-driven channel suspension flow. In such a system, particles are driven toward the channel center, resulting in a volume fraction gradient across the channel and a plug region where the jamming volume fraction as measured in shear flow may be exceeded. The flow profile is thus highly nonlinear, both due to the pressure gradient and the time-varying concentration gradient, which justifies the use of the present numerical method. We compare the results to a modified version of the Suspension Balance Model, confirming the well-known limits of such a model in the plug region. Finally, we show that the computed stresses obey the usual momentum balance in the channel: the tangential stress is indeed driven by the pressure gradient as predicted by theory and the second normal stress (parallel to the velocity gradient) does not vary over the channel width, although contact and hydrodynamic contributions significantly do over both time and space. This confirms the good behavior of the present numerical method when nonlinear flows are tackled.File | Dimensione | Formato | |
---|---|---|---|
Thesis.pdf
accesso aperto
Descrizione: Ph.D. Thesis
Tipologia:
Tesi di dottorato
Licenza:
Pubblico - Tutti i diritti riservati
Dimensione
11.63 MB
Formato
Adobe PDF
|
11.63 MB | Adobe PDF | Visualizza/Apri |
Pubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.
https://hdl.handle.net/11583/2996642