Lagrangian stochastic model for the orientation of inertialess spheroidal particles in turbulent flows: an efficient numerical method for CFD approach

Abstract

In this work, we propose a model for the orientation of inertialess spheroidal particles suspended in turbulent flows. This model consists in a stochastic version of the Jeffery equation that can be included in a statistical Lagrangian description of particles suspended in a flow. It is compatible and coherent with turbulence models that are widely used in CFD codes for the simulation of the flow field in practical large-scale applications. In this context, we propose and analyze a numerical scheme based on a splitting scheme algorithm that decouples the orientation dynamics into its main contributions: stretching and rotation. We detail its implementation in an open-source CFD software. We analyze the weak and strong convergence of both the global scheme and of each sub-part. Subsequently, the splitting technique yields to a highly efficient hybrid algorithm coupling pure probabilistic and deterministic numerical schemes. Various numerical experiments were implemented and the results were compared with analytical predictions of the model to assess the algorithm efficiency and accuracy.

References in corpus (11)

• An Immersed Boundary Method with Direct Forcing for the Simulation of Particulate Flows
• Rotation Rate of Rods in Turbulent Fluid Flow
• Shape-dependence of particle rotation in isotropic turbulence
• Alignment of vorticity and rods with Lagrangian fluid stretching in turbulence
• Tumbling of small axisymmetric particles in random and turbulent flows
• Orientation dynamics of small, triaxial-ellipsoidal particles in isotropic turbulence
• A closure for Lagrangian velocity gradient evolution in turbulence using recent deformation mapping of initially Gaussian fields
• Stretching and buckling of small elastic fibers in turbulence
• A note on the exact simulation of spherical Brownian motion
• Brownian Motion on a Sphere: Distribution of Solid Angles
• On the weak convergence rate of an exponential Euler scheme for SDEs governed by coefficients with superlinear growth