# Spectral Galerkin approximation of Fokker-Planck equations with unbounded drift

ESAIM: Mathematical Modelling and Numerical Analysis (2008)

- Volume: 43, Issue: 3, page 445-485
- ISSN: 0764-583X

## Access Full Article

top## Abstract

top## How to cite

topKnezevic, David J., and Süli, Endre. "Spectral Galerkin approximation of Fokker-Planck equations with unbounded drift." ESAIM: Mathematical Modelling and Numerical Analysis 43.3 (2008): 445-485. <http://eudml.org/doc/194458>.

@article{Knezevic2008,

abstract = {
This paper is concerned with the analysis and implementation of spectral Galerkin
methods for a class of Fokker-Planck equations that arises
from the kinetic theory of dilute polymers. A relevant feature of the class of equations
under consideration from the viewpoint of mathematical analysis and numerical approximation is
the presence of an unbounded drift coefficient, involving a smooth convex potential U that is equal to +∞ along
the boundary ∂D of the computational domain D.
Using a symmetrization of the differential operator based on the Maxwellian M corresponding to U,
which vanishes along ∂D, we remove the unbounded drift coefficient at the expense
of introducing a degeneracy, through M, in the principal part of the operator.
The general class of admissible potentials considered
includes the FENE (finitely extendible nonlinear elastic) model.
We show the existence of weak solutions to the initial-boundary-value problem, and develop a fully-discrete
spectral Galerkin method for such degenerate Fokker-Planck equations
that exhibits optimal-order convergence in the Maxwellian-weighted H1 norm on D.
In the case of the FENE model, we also discuss variants of these analytical results when the Fokker-Planck equation is subjected to an alternative class of transformations proposed by Chauvière and Lozinski; these map the original Fokker-Planck operator with an unbounded drift coefficient into Fokker-Planck operators with unbounded drift and reaction coefficients, that have improved coercivity properties in comparison with the
original operator. The analytical results are illustrated by numerical experiments for the FENE model in two space dimensions.
},

author = {Knezevic, David J., Süli, Endre},

journal = {ESAIM: Mathematical Modelling and Numerical Analysis},

keywords = {Spectral methods; Fokker-Planck equations; transport-diffusion problems; FENE.; spectral methods; finitely extendible nonlinear elastic model (FENE); unbounded coefficients; polymer chains},

language = {eng},

month = {12},

number = {3},

pages = {445-485},

publisher = {EDP Sciences},

title = {Spectral Galerkin approximation of Fokker-Planck equations with unbounded drift},

url = {http://eudml.org/doc/194458},

volume = {43},

year = {2008},

}

TY - JOUR

AU - Knezevic, David J.

AU - Süli, Endre

TI - Spectral Galerkin approximation of Fokker-Planck equations with unbounded drift

JO - ESAIM: Mathematical Modelling and Numerical Analysis

DA - 2008/12//

PB - EDP Sciences

VL - 43

IS - 3

SP - 445

EP - 485

AB -
This paper is concerned with the analysis and implementation of spectral Galerkin
methods for a class of Fokker-Planck equations that arises
from the kinetic theory of dilute polymers. A relevant feature of the class of equations
under consideration from the viewpoint of mathematical analysis and numerical approximation is
the presence of an unbounded drift coefficient, involving a smooth convex potential U that is equal to +∞ along
the boundary ∂D of the computational domain D.
Using a symmetrization of the differential operator based on the Maxwellian M corresponding to U,
which vanishes along ∂D, we remove the unbounded drift coefficient at the expense
of introducing a degeneracy, through M, in the principal part of the operator.
The general class of admissible potentials considered
includes the FENE (finitely extendible nonlinear elastic) model.
We show the existence of weak solutions to the initial-boundary-value problem, and develop a fully-discrete
spectral Galerkin method for such degenerate Fokker-Planck equations
that exhibits optimal-order convergence in the Maxwellian-weighted H1 norm on D.
In the case of the FENE model, we also discuss variants of these analytical results when the Fokker-Planck equation is subjected to an alternative class of transformations proposed by Chauvière and Lozinski; these map the original Fokker-Planck operator with an unbounded drift coefficient into Fokker-Planck operators with unbounded drift and reaction coefficients, that have improved coercivity properties in comparison with the
original operator. The analytical results are illustrated by numerical experiments for the FENE model in two space dimensions.

LA - eng

KW - Spectral methods; Fokker-Planck equations; transport-diffusion problems; FENE.; spectral methods; finitely extendible nonlinear elastic model (FENE); unbounded coefficients; polymer chains

UR - http://eudml.org/doc/194458

ER -

## References

top- A. Ammar, B. Mokdad, F. Chinesta and R. Keunings, A new family of solvers for some classes of multidimensional partial differential equations encountered in kinetic theory modeling of complex fluids. J. Non-Newtonian Fluid Mech.139 (2006) 153–176.
- A. Ammar, B. Mokdad, F. Chinesta and R. Keunings, A new family of solvers for some classes of multidimensional partial differential equations encountered in kinetic theory modelling of complex fluids. Part II: Transient simulation using space-time separated representations. J. Non-Newtonian Fluid Mech.144 (2007) 98–121.
- F.G. Avkhadiev and K.-J. Wirths, Unified Poincaré and Hardy inequalities with sharp constants for convex domains. ZAMM Z. Angew. Math. Mech.87 (2007) 632–642.
- J.W. Barrett and E. Süli, Existence of global weak solutions to kinetic models of dilute polymers. Multiscale Model. Simul.6 (2007) 506–546.
- J.W. Barrett and E. Süli, Existence of global weak solutions to dumbbell models for dilute polymers with microscopic cut-off. Math. Mod. Meth. Appl. Sci.18 (2008) 935–971.
- J.W. Barrett and E. Süli, Numerical approximation of corotational dumbbell models for dilute polymers. IMA J. Numer. Anal. (2008) online. Available at . URIhttp://imajna.oxfordjournals.org/cgi/content/abstract/drn022
- J.W. Barrett, C. Schwab and E. Süli, Existence of global weak solutions for some polymeric flow models. Math. Mod. Meth. Appl. Sci.15 (2005) 939–983.
- C. Bernardi, Optimal finite-element interpolation on curved domains. SIAM J. Numer. Anal.26 (1989) 1212–1240.
- C. Bernardi and Y. Maday, Spectral methods, in Handbook of Numerical AnalysisV, P. Ciarlet and J. Lions Eds., Elsevier (1997).
- O.V. Besov and A. Kufner, The density of smooth functions in weight spaces. Czechoslova. Math. J.18 (1968) 178–188.
- O.V. Besov, J. Kadlec and A. Kufner, Certain properties of weight classes. Dokl. Akad. Nauk SSSR171 (1966) 514–516.
- R.B. Bird, C.F. Curtiss, R.C. Armstrong and O. Hassager, Dynamics of Polymeric Liquids, Vol.1, Fluid Mechanics. Second edition, John Wiley and Sons (1987).
- R.B. Bird, C.F. Curtiss, R.C. Armstrong and O. Hassager, Dynamics of Polymeric Liquids, Vol.2, Kinetic Theory. Second edition, John Wiley and Sons (1987).
- S. Bobkov and M. Ledoux, From Brunn-Minkowski to Brascamp-Lieb and to logarithmic Sobolev inequalities. Geom. Funct. Anal.10 (2000) 1028–1052.
- C. Canuto, A. Quarteroni, M.Y. Hussaini and T.A. Zang, Spectral Methods: Fundamentals in Single Domains. Springer (2006).
- S. Cerrai, Second Order PDE's in Finite and Infinite Dimensions, A Probabilistic Approach, Lecture Notes in Mathematics1762. Springer (2001).
- C. Chauvière and A. Lozinski, Simulation of complex viscoelastic flows using Fokker-Planck equation: 3D FENE model. J. Non-Newtonian Fluid Mech.122 (2004) 201–214.
- C. Chauvière and A. Lozinski, Simulation of dilute polymer solutions using a Fokker-Planck equation. Comput. Fluids33 (2004) 687–696.
- G. Da Prato and A. Lunardi, On a class of elliptic operators with unbounded coefficients in convex domains. Atti Accad. Naz. Lincei Cl. Sci. Fis. Mat. Natur. Rend. Lincei (9) Mat. Appl.15 (2004) 315–326.
- Q. Du, C. Liu and P. Yu, FENE dumbbell model and its several linear and nonlinear closure approximations. Multiscale Model. Simul.4 (2005) 709–731.
- H. Eisen, W. Heinrichs and K. Witsch, Spectral collocation methods and polar coordinate singularities. J. Comput. Phys.96 (1991) 241–257.
- M. Golubitsky and V. Guillemin, Stable Mappings and Their Singularities. Springer (1973).
- B. Jourdain, T. Lelièvre and C. Le Bris, Numerical analysis of micro-macro simulations of polymeric fluid flows: A simple case. Math. Mod. Meth. Appl. Sci.12 (2002) 1205–1243.
- A.N. Kolmogorov, Über die analytischen Methoden in der Wahrscheinlichkeitsrechnung. Math. Ann.104 (1931).
- A. Kufner, Weighted Sobolev Spaces, Teubner-Texte zur Mathematik. Teubner (1980).
- T. Li and P.-W. Zhang, Mathematical analysis of multi-scale models of complex fluids. Commun. Math. Sci.5 (2007) 1–51.
- A. Lozinski and C. Chauvière, A fast solver for Fokker-Planck equation applied to viscoelastic flows calculation: 2D FENE model. J. Comput. Phys.189 (2003) 607–625.
- M. Marcus, V.J. Mizel and Y. Pinchover, On the best constant for Hardy's inequality in Rn. Trans. Amer. Math. Soc.350 (1998) 3237–3255.
- T. Matsushima and P.S. Marcus, A spectral method for polar coordinates. J. Comput. Phys.120 (1995) 365–374.
- H.C. Öttinger, Stochastic Processes in Polymeric Fluids. Springer (1996).
- J. Shen, Efficient spectral Galerkin methods III: Polar and cylindrical geometries. SIAM J. Sci. Comput.18 (1997) 1583–1604.
- R. Temam, Navier-Stokes Equations: Theory and Numerical Analysis. Third edition, North-Holland, Amsterdam (1984).
- H. Triebel, Interpolation Theory, Function Spaces, Differential Operators. Second edition, Johan Ambrosius Barth, Heidelberg (1995).
- W.T.M. Verkley, A spectral model for two-dimensional incompressible fluid flow in a circular basin I. Mathematical formulation. J. Comput. Phys.136 (1997) 100–114.

## Citations in EuDML Documents

top- John W. Barrett, Endre Süli, Finite element approximation of kinetic dilute polymer models with microscopic cut-off
- John W. Barrett, Endre Süli, Finite element approximation of kinetic dilute polymer models with microscopic cut-off
- John W. Barrett, Endre Süli, Finite element approximation of finitely extensible nonlinear elastic dumbbell models for dilute polymers
- David J. Knezevic, Endre Süli, A heterogeneous alternating-direction method for a micro-macro dilute polymeric fluid model
- John W. Barrett, Endre Süli, Finite element approximation of finitely extensible nonlinear elastic dumbbell models for dilute polymers

## NotesEmbed ?

topTo embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.