Authors: Wegener, Katharina
Kuzmin, Dmitri
Turek, Stefan
Title: Efficient numerical solution of the Fokker-Planck equation using physics-conforming finite element methods
Language (ISO): en
Abstract: We consider the Fokker-Planck equation (FPE) for the orientation proba­bility density of fiber suspensions. Using the continuous Galerkin method, we express the numerical solution in terms of Lagrange basis functions that are associated with N nodes of a computational mesh for a domain in the 3D physical space and M nodes of a mesh for the surface of a unit sphere representing the configuration space. The NM time-dependent unknowns of our finite element approximations are probabilities corresponding to dis­crete space locations and orientation angles. The framework of alternating­direction methods enables us to update the numerical solution in parallel by solving N evolution equations on the sphere and M three-dimensional advection equations in each (pseudo-)time step. To ensure positivity preservation as well as the normalization property of the probability density, we perform algebraic flux correction for each equation and synchronize the correction factors corresponding to different orientation angles. The velocity field for the spatial advection step is obtained using a Schur complement method to solve a generalized system of the incompressible Navier-Stokes equations (NSE). Fiber-induced subgrid-scale effects are taken into account using an effective stress tensor that depends on the second- and fourth-order mo­ments of the orientation density function. Numerical studies are performed for individual subproblems and for the coupled FPE-NSE system.
Subject Headings: fiber suspensions
flux limiting
positivity preservation
alternating-direction methods
finite elements
Fokker-Planck equation
Issue Date: 2023-01
Appears in Collections:Ergebnisberichte des Instituts für Angewandte Mathematik

Files in This Item:
File Description SizeFormat 
Ergebnisbericht Nr. 658.pdfDNB17.94 MBAdobe PDFView/Open

This item is protected by original copyright

Items in Eldorado are protected by copyright, with all rights reserved, unless otherwise indicated.