TY - JOUR

T1 - Numerical simulation of the time-fractional Fokker–Planck equation and applications to polymeric fluids

AU - Beddrich, Jonas

AU - Süli, Endre

AU - Wohlmuth, Barbara

N1 - Publisher Copyright:
© 2023 Elsevier Inc.

PY - 2024/1/15

Y1 - 2024/1/15

N2 - We introduce a new approach to the numerical approximation of the time-fractional Navier–Stokes–Fokker–Planck (TFNSFP) system, which involves the coupling of the incompressible Navier–Stokes equations to a time-fractional Fokker–Planck equation. The model arises in the context of dilute polymeric fluids, and it enhances the standard integer-derivative version of the model by including memory effects. The key challenge associated with the numerical solution of the TFNSFP system is that, in addition to it being nonlocal in time, it is, even in its simplest form, posed on a spatial domain that is the Cartesian product of two d-dimensional domains, for d∈{2,3}, the d-dimensional flow domain and the d-dimensional configuration space domain. We combine a kernel compression technique based on rational approximation of the integral kernel of the time-fractional derivative with a space-splitting method. By doing so, we transform the time-fractional partial differential equation into a fixed number of integer-order in time partial differential equations. The Fokker–Planck equation posed on the 2d-dimensional domain is decoupled into two d-dimensional problems, and a standard combination of the Hermite spectral method on the configuration space domain and a finite element method on the flow domain is then applied. Finally, we combine our numerical scheme for the time-fractional Fokker–Planck equation (TFFPE) with a standard Navier–Stokes solver. We propose an efficient implementation strategy based on the algebraic structure of the discretized time-fractional Fokker–Planck equation, which significantly reduces the computational cost.

AB - We introduce a new approach to the numerical approximation of the time-fractional Navier–Stokes–Fokker–Planck (TFNSFP) system, which involves the coupling of the incompressible Navier–Stokes equations to a time-fractional Fokker–Planck equation. The model arises in the context of dilute polymeric fluids, and it enhances the standard integer-derivative version of the model by including memory effects. The key challenge associated with the numerical solution of the TFNSFP system is that, in addition to it being nonlocal in time, it is, even in its simplest form, posed on a spatial domain that is the Cartesian product of two d-dimensional domains, for d∈{2,3}, the d-dimensional flow domain and the d-dimensional configuration space domain. We combine a kernel compression technique based on rational approximation of the integral kernel of the time-fractional derivative with a space-splitting method. By doing so, we transform the time-fractional partial differential equation into a fixed number of integer-order in time partial differential equations. The Fokker–Planck equation posed on the 2d-dimensional domain is decoupled into two d-dimensional problems, and a standard combination of the Hermite spectral method on the configuration space domain and a finite element method on the flow domain is then applied. Finally, we combine our numerical scheme for the time-fractional Fokker–Planck equation (TFFPE) with a standard Navier–Stokes solver. We propose an efficient implementation strategy based on the algebraic structure of the discretized time-fractional Fokker–Planck equation, which significantly reduces the computational cost.

KW - Dilute polymeric fluids

KW - Hermite spectral method

KW - Hookean bead-spring-model

KW - Navier–Stokes–Fokker–Planck system

KW - Numerical simulation

KW - Time-fractional PDE

UR - http://www.scopus.com/inward/record.url?scp=85176275363&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2023.112598

DO - 10.1016/j.jcp.2023.112598

M3 - Article

AN - SCOPUS:85176275363

SN - 0021-9991

VL - 497

JO - Journal of Computational Physics

JF - Journal of Computational Physics

M1 - 112598

ER -