Simulation of finite-size fibers in turbulent channel flows

Assessment of Analytical Orientation Prediction Models for Suspensions Containing Fibers and Spheres

Analytical orientation models like the Folgar Tucker (FT) model are widely applied to predict the orientation of suspended non-spherical particles. The accuracy of these models depends on empirical model parameters. In this work, we assess how well analytical orientation models can predict the orientation of suspensions not only consisting of fibers but also of an additional second particle type in the shape of disks, which are varied in size and filling fraction. We mainly focus on the FT model, and we also compare its accuracy to more complex models like Reduced-Strain Closure model (RSC), Moldflow Rotational Diffusion model (MRD), and Anisotropic Rotary Diffusion model (ARD). In our work, we address the following questions. First, can the FT model predict the orientation of suspensions despite the additional particle phase affecting the rotation of the fibers? Second, is it possible to formulate an expression for the sole Folgar Tucker model parameter that is based on the suspension composition? Third, is there an advantage to choose more complex orientation prediction models that require the adjustment of additional model parameters?

Influence of numerical resolution on the dynamics of finite-size particles with the lattice Boltzmann method

We investigate and compare the accuracy and efficiency of different numerical approaches to model the dynamics of finite-size particles using the lattice Boltzmann method (LBM). This includes the standard bounce-back (BB) and the equilibrium interpolation (EI) schemes. To accurately compare the different implementations, we first introduce a boundary condition to approximate the flow properties of an unbounded fluid in a finite simulation domain, taking into account the perturbation induced by a moving particle. We show that this boundary treatment is efficient in suppressing detrimental effects on the dynamics of spherical and ellipsoidal particles arising from the finite size of the simulation domain. We then investigate the performances of the BB and EI schemes in modeling the dynamics of a spherical particle settling under Stokes conditions, which can now be reproduced with great accuracy thanks to the treatment of the exterior boundary. We find that the EI scheme outperforms the BB scheme in providing a better accuracy scaling with respect to the resolution of the settling particle, while suppressing finite-size effects due to the particle discretization on the lattice grid. Additionally, in order to further increase the capability of the algorithm in modeling particles of sizes comparable to the lattice spacing, we propose an improvement to the EI scheme, the complete equilibrium interpolation (CEI). This approach allows us to accurately capture the boundaries of the particle also when located between two fluid nodes. We evaluate the CEI performance in solving the dynamics of an under-resolved particle under analogous Stokes conditions and also for the case of a rotating ellipsoid in a shear flow. Finally, we show that EI and CEI are able to recover the correct flow solutions also at small, but finite, Reynolds number. Adopting the CEI scheme it is not only possible to detect particles with zero lattice occupation, but also to increase up to one order of magnitude the accuracy of the dynamics of particles with a size comparable to the lattice spacing with respect to the BB and the EI schemes.

Fiber Orientation Predictions—A Review of Existing Models

Fiber reinforced polymers are key materials across different industries. The manufacturing processes of those materials have typically strong impact on their final microstructure, which at the same time controls the mechanical performance of the part. A reliable virtual engineering design of fiber-reinforced polymers requires therefore considering the simulation of the process-induced microstructure. One relevant microstructure descriptor in fiber-reinforced polymers is the fiber orientation. This work focuses on the modeling of the fiber orientation phenomenon and presents a historical review of the different modelling approaches. In this context, the article describes different macroscopic fiber orientation models such as the Folgar-Tucker, nematic, reduced strain closure (RSC), retarding principal rate (RPR), anisotropic rotary diffusion (ARD), principal anisotropic rotary diffusion (pARD), and Moldflow rotary diffusion (MRD) model. We discuss briefly about closure approximations, which are a common mathematical element of those macroscopic fiber orientation models. In the last section, we introduce some micro-scale numerical methods for simulating the fiber orientation phenomenon, such as the discrete element method (DEM), the smoothed particle hydrodynamics (SPH) method and the moving particle semi-implicit (MPS) method.

Effect of particle shape on fluid statistics and particle dynamics in turbulent pipe flow

Anisotropic particles are present in many natural and industrial flows. Here we perform direct numerical simulation (DNS) of turbulent pipe flows with dispersed finite-size prolate spheroids simulated by means of the lattice Boltzmann method (LBM). We consider three different particle shapes: spheroidal (aspect ratio 2 and 3) and spherical. These three simulations are complemented with a reference simulation of a single-phase flow. For the sake of comparison, all simulations, laden or unladen have the same energy input. The flow geometry used is a straight pipe with length eight times its radius where the fluid is randomly seeded with 256 finite-size particles. The volume fraction of particles in the flow has been kept fixed at 0.48% by varying the major and minor axis of each particle such that their volume remains the same. We studied the effect of different particle shapes on particle dynamics and orientation, as well as on the flow modulation. We show that the local accumulation of spheres close to the wall decreases for spheroids with increasing aspect ratio. These spheroidal particles rotate slower than spheres near to the wall and tend to stay with their major axes aligned to the flow streamwise direction. Despite the lower rotation rates, a higher intermittency in the rotational rates was observed for spheroids and this increase at increasing the aspect ratio. The drag reduction observed for particles with higher aspect ratio have also been investigated using the one-dimensional energy and dissipation spectra. These results point to the relevance of particle shapes on their dynamics and their influence on the turbulent flow.

Computational study of radial particle migration and stresslet distributions in particle-laden turbulent pipe flow

Particle-laden turbulent flows occur in a variety of industrial applications as well as in naturally occurring flows. While the numerical simulation of such flows has seen significant advances in recent years, it still remains a challenging problem. Many studies investigated the rheology of dense suspensions in laminar flows as well as the dynamics of point-particles in turbulence. Here we employ a fully-resolved numerical simulation based on a lattice Boltzmann scheme, to investigate turbulent flow with large neutrally buoyant particles in a pipe flow at low Reynolds number and in dilute regimes. The energy input is kept fixed resulting in a Reynolds number based on the friction velocity around 250. Two different particle radii were used giving a particle-pipe diameter ratio of 0.05 and 0.075. The number of particles is kept constant resulting in a volume fraction of 0.54% and 1.83%, respectively. We investigated Eulerian and Lagrangian statistics along with the stresslet exerted by the fluid on the spherical particles. It was observed that the high particle-to-fluid slip velocity close to the wall corresponds locally to events of high energy dissipation, which are not present in the single-phase flow. The migration of particles from the inner to the outer region of the pipe, the dependence of the stresslet on the particle radial positions and a proxy for the fragmentation rate of the particles computed using the stresslet have been investigated.

Effects of particle-fluid density ratio on the interactions between the turbulent channel flow and finite-size particles

A parallel direct-forcing fictitious domain method is employed to perform fully resolved numerical simulations of turbulent channel flow laden with finite-size particles. The effects of the particle-fluid density ratio on the turbulence modulation in the channel flow are investigated at the friction Reynolds number of 180, the particle volume fraction of 0.84%, and the particle-fluid density ratio ranging from 1 to 104.2. The results show that the variation of the flow drag with the particle-fluid density ratio is not monotonic, with a larger flow drag for the density ratio of 10.42, compared to those of unity and 104.2. A significant drag reduction by the particles is observed for large particle-fluid density ratios during the transient stage, but not at the statistically stationary stage. The intensity of particle velocity fluctuations generally decreases with increasing particle inertia, except that the particle streamwise root-mean-square velocity and streamwise-transverse velocity correlation in the near-wall region are largest at the density ratio of the order of 10. The averaged momentum equations are derived with the spatial averaging theorem and are used to analyze the mechanisms for the effects of the particles on the flow drag. The results indicate that the drag-reduction effect due to the decrease in the fluid Reynolds shear stress is counteracted by the drag-enhancement effect due to the increase in the total particle stress or the interphase drag force for the large particle-inertia case. The sum of the total Reynolds stress and particle inner stress contributions to the flow drag is largest at the density ratio of the order of 10, which is the reason for the largest flow drag at this density ratio. The interphase drag force obtained from the averaged momentum equation (the balance theory) is significantly smaller than (but agrees qualitatively with) that from the empirical drag formula based on the phase-averaged slip velocity for large density ratios. For the neutrally buoyant case, the balance theory predicts a positive interphase force on the particles arising from the negative gradient of the particle inner stress, which cannot be predicted by the drag formula based on the phase-averaged slip velocity. In addition, our results show that both particle collision and particle-turbulence interaction play roles in the formation of the inhomogeneous distribution of the particles at the density ratio of the order of 10.

Orientational dynamics of a triaxial ellipsoid in simple shear flow: Influence of inertia

The motion of a single ellipsoidal particle in simple shear flow can provide valuable insights toward understanding suspension flows with nonspherical particles. Previously, extensive studies have been performed on the ellipsoidal particle with rotational symmetry, a so-called spheroid. The nearly prolate ellipsoid (one major and two minor axes of almost equal size) is known to perform quasiperiodic or even chaotic orbits in the absence of inertia. With small particle inertia, the particle is also known to drift toward this irregular motion. However, it is not previously understood what effects from fluid inertia could be, which is of highest importance for particles close to neutral buoyancy. Here, we find that fluid inertia is acting strongly to suppress the chaotic motion and only very weak fluid inertia is sufficient to stabilize a rotation around the middle axis. The mechanism responsible for this transition is believed to be centrifugal forces acting on fluid, which is dragged along with the rotational motion of the particle. With moderate fluid inertia, it is found that nearly prolate triaxial particles behave similarly to the perfectly spheroidal particles. Finally, we also are able to provide predictions about the stable rotational states for the general triaxial ellipsoid in simple shear with weak inertia.