Efficient algorithms for constraining orientation tensors in Galerkin methods for the Fokker-Planck equation. (English) Zbl 1443.65214
Summary: This paper deals with the problem of positivity preservation in numerical algorithms for simulating fiber suspension flows. In contrast to fiber orientation models based on the Advani-Tucker evolution equations for even-order orientation tensors, the probability distribution function of fiber orientation is approximated using the Galerkin discretization of the Fokker-Planck equation with Fourier basis functions or spherical harmonics. This procedure leads to a natural generalization of orientation tensor models replacing ad hoc closure approximations by Galerkin equations for the fine-scale components. After introducing an operator splitting approach to solving the discretized Fokker-Planck equation, we present conditions and correction techniques that guarantee physically correct distribution functions. As the reader will see, the derivation of these conditions is independent of the space dimension and their applicability is not limited to the simulation of fiber suspensions.
MSC:
| 65M60 | Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs |
| 82C31 | Stochastic methods (Fokker-Planck, Langevin, etc.) applied to problems in time-dependent statistical mechanics |
Keywords:
fiber suspension flows; Fokker-Planck equation; orientation tensors; Galerkin approximation; Fourier analysis; spherical harmonicsReferences:
| [1] | Advani, Suresh G.; Tucker, Charles L., The use of tensors to describe and predict fiber orientation in short fiber composites, J. Rheol., 31, 8, 751-784 (1987), (1978-present) |
| [2] | Donea, Jean; Huerta, Antonio, Finite Element Methods for Flow Problems (2003), John Wiley & Sons |
| [3] | Kuzmin, Dmitri; Hämäläinen, Jari, Finite Element Methods for Computational Fluid Dynamics: A Practical Guide (2014), SIAM · Zbl 1317.76001 |
| [4] | Löhner, Rainald, Applied CFD Techniques: An Introduction Based on Finite Element Methods (2008), John Wiley & Sons · Zbl 1151.76002 |
| [5] | Akbar, Sameer; Cengiz Altan, M., On the solution of fiber orientation in two-dimensional homogeneous flows, Polym. Eng. Sci., 32, 12, 810-822 (1992) |
| [6] | Fokker, Adriaan Daniël, Die mittlere energie rotierender elektrischer dipole im strahlungsfeld, Ann. Phys., 348, 5, 810-820 (1914) |
| [7] | Montgomery-Smith, Stephen; He, Wei; Jack, David A.; Smith, Douglas E., Exact tensor closures for the three-dimensional Jeffery’s equation, J. Fluid Mech., 680, 321-335 (2011) · Zbl 1241.76402 |
| [8] | Jeffery, George Barker, The motion of ellipsoidal particles immersed in a viscous fluid, Proc. R. Soc. Lond. Ser. A, 102, 161-179 (1922) · JFM 49.0748.02 |
| [9] | Folgar, Fransisco; Tucker, Charles L., Orientation behavior of fibers in concentrated suspensions, J. Reinf. Plast. Compos., 3, 2, 98-119 (1984) |
| [10] | Glenn Lipscomb, G.; Denn, Morton M.; Hur, D. U.; Boger, David V., The Flow of Fiber Suspensions in Complex Geometries, J. Non-Newton. Fluid Mech., 26, 3, 297-325 (1988) |
| [11] | Cintra, Joaquim S.; Tucker, Charles L., Orthotropic closure approximations for flow-induced fiber orientation, J. Rheol., 39, 6, 1095-1122 (1995) |
| [12] | Lohmann, Christoph, Galerkin-Spektralverfahren für die Fokker-Planck-Gleichung (Galerkin spectral methods for the Fokker-Planck equation) (2016), TU Dortmund University, in press, Springer Spektrum |
| [13] | Vincenzi, Dario, Orientation of non-spherical particles in an axisymmetric random flow, J. Fluid Mech., 719, 3, 465-487 (2013) · Zbl 1284.76385 |
| [14] | Albert, Abraham Adrian, An inductive proof of descartes’ rule of signs, Amer. Math. Monthly, 50, 3, 178-180 (1943) · Zbl 0060.05004 |
| [15] | Anderson, Bruce; Jackson, Jeffrey; Sitharam, Meera, Descartes’ rule of signs revisited, Amer. Math. Monthly, 105, 5, 447-451 (1998) · Zbl 0913.12001 |
| [16] | Kershaw, David S., Flux limiting nature’s own way—A new method for numerical solution of the transport equation. Technical Report UCRL-78378 (1976), Lawrence Livermore National Laboratory |
| [17] | Cengiz Altan, M.; Tang, Lan, Orientation tensors in simple flows of dilute suspensions of non-Brownian rigid ellipsoids, comparison of analytical and approximate solutions, Rheol. Acta, 32, 3, 227-244 (1993) |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
