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 that is equal to along the boundary of the computational domain . Using a symmetrization of the differential operator based on the Maxwellian corresponding to , which vanishes along , we remove the unbounded drift coefficient at the expense of introducing a degeneracy, through , 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 norm on . 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.
@article{M2AN_2009__43_3_445_0,
author = {Knezevic, David J. and S\"uli, Endre},
title = {Spectral Galerkin approximation of Fokker-Planck equations with unbounded drift},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Mod\'elisation Math\'ematique et Analyse Num\'erique},
volume = {43},
year = {2009},
pages = {445-485},
doi = {10.1051/m2an:2008051},
mrnumber = {2536245},
zbl = {pre05574327},
language = {en},
url = {http://dml.mathdoc.fr/item/M2AN_2009__43_3_445_0}
}
Knezevic, David J.; Süli, Endre. Spectral Galerkin approximation of Fokker-Planck equations with unbounded drift. ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Tome 43 (2009) pp. 445-485. doi : 10.1051/m2an:2008051. http://gdmltest.u-ga.fr/item/M2AN_2009__43_3_445_0/
[1] , , and , 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.
[2] , , and , 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.
[3] and , Unified Poincaré and Hardy inequalities with sharp constants for convex domains. ZAMM Z. Angew. Math. Mech. 87 (2007) 632-642. | MR 2354734 | Zbl 1145.26005
[4] and , Existence of global weak solutions to kinetic models of dilute polymers. Multiscale Model. Simul. 6 (2007) 506-546. | MR 2338493 | Zbl pre05286335
[5] and , Existence of global weak solutions to dumbbell models for dilute polymers with microscopic cut-off. Math. Mod. Meth. Appl. Sci. 18 (2008) 935-971. | MR 2419205 | Zbl 1158.35070
[6] and , Numerical approximation of corotational dumbbell models for dilute polymers. IMA J. Numer. Anal. (2008) online. Available at http://imajna.oxfordjournals.org/cgi/content/abstract/drn022. | MR 2557051 | Zbl pre05642232
[7] , and , Existence of global weak solutions for some polymeric flow models. Math. Mod. Meth. Appl. Sci. 15 (2005) 939-983. | MR 2149930 | Zbl 1161.76453
[8] , Optimal finite-element interpolation on curved domains. SIAM J. Numer. Anal. 26 (1989) 1212-1240. | MR 1014883 | Zbl 0678.65003
[9] and , Spectral methods, in Handbook of Numerical Analysis V, P. Ciarlet and J. Lions Eds., Elsevier (1997). | MR 1470226
[10] and , The density of smooth functions in weight spaces. Czechoslova. Math. J. 18 (1968) 178-188. | MR 223877 | Zbl 0193.41502
[11] , and , Certain properties of weight classes. Dokl. Akad. Nauk SSSR 171 (1966) 514-516. | MR 204842 | Zbl 0168.11103
[12] , , and , Dynamics of Polymeric Liquids, Vol. 1, Fluid Mechanics. Second edition, John Wiley and Sons (1987).
[13] , , and , Dynamics of Polymeric Liquids, Vol. 2, Kinetic Theory. Second edition, John Wiley and Sons (1987).
[14] and , From Brunn-Minkowski to Brascamp-Lieb and to logarithmic Sobolev inequalities. Geom. Funct. Anal. 10 (2000) 1028-1052. | MR 1800062 | Zbl 0969.26019
[15] , , and , Spectral Methods: Fundamentals in Single Domains. Springer (2006). | MR 2223552 | Zbl 1093.76002
[16] , Second Order PDE's in Finite and Infinite Dimensions, A Probabilistic Approach, Lecture Notes in Mathematics 1762. Springer (2001). | MR 1840644 | Zbl 0983.60004
[17] and , Simulation of complex viscoelastic flows using Fokker-Planck equation: 3D FENE model. J. Non-Newtonian Fluid Mech. 122 (2004) 201-214. | Zbl 1131.76307
[18] and , Simulation of dilute polymer solutions using a Fokker-Planck equation. Comput. Fluids 33 (2004) 687-696. | Zbl 1100.76549
[19] and , 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. | MR 2148888 | Zbl 1162.35345
[20] , and , FENE dumbbell model and its several linear and nonlinear closure approximations. Multiscale Model. Simul. 4 (2005) 709-731. | MR 2203938 | Zbl 1108.76006
[21] , and , Spectral collocation methods and polar coordinate singularities. J. Comput. Phys. 96 (1991) 241-257. | MR 1128222 | Zbl 0731.65095
[22] and , Stable Mappings and Their Singularities. Springer (1973). | MR 341518 | Zbl 0294.58004
[23] , and , Numerical analysis of micro-macro simulations of polymeric fluid flows: A simple case. Math. Mod. Meth. Appl. Sci. 12 (2002) 1205-1243. | MR 1927023 | Zbl 1041.76003
[24] , Über die analytischen Methoden in der Wahrscheinlichkeitsrechnung. Math. Ann. 104 (1931). | Zbl 0001.14902
[25] , Weighted Sobolev Spaces, Teubner-Texte zur Mathematik. Teubner (1980). | MR 664599 | Zbl 0455.46034
[26] and , Mathematical analysis of multi-scale models of complex fluids. Commun. Math. Sci. 5 (2007) 1-51. | MR 2310632 | Zbl 1129.76006
[27] and , A fast solver for Fokker-Planck equation applied to viscoelastic flows calculation: 2D FENE model. J. Comput. Phys. 189 (2003) 607-625. | MR 1996059 | Zbl 1060.82525
[28] , and , On the best constant for Hardy’s inequality in . Trans. Amer. Math. Soc. 350 (1998) 3237-3255. | MR 1458330 | Zbl 0917.26016
[29] and , A spectral method for polar coordinates. J. Comput. Phys. 120 (1995) 365-374. | MR 1349468 | Zbl 0842.65051
[30] , Stochastic Processes in Polymeric Fluids. Springer (1996). | MR 1383323 | Zbl 0995.60098
[31] , Efficient spectral Galerkin methods III: Polar and cylindrical geometries. SIAM J. Sci. Comput. 18 (1997) 1583-1604. | MR 1480626 | Zbl 0890.65117
[32] , Navier-Stokes Equations: Theory and Numerical Analysis. Third edition, North-Holland, Amsterdam (1984). | MR 769654 | Zbl 0568.35002
[33] , Interpolation Theory, Function Spaces, Differential Operators. Second edition, Johan Ambrosius Barth, Heidelberg (1995). | MR 1328645 | Zbl 0830.46028
[34] , A spectral model for two-dimensional incompressible fluid flow in a circular basin I. Mathematical formulation. J. Comput. Phys. 136 (1997) 100-114. | MR 1468626 | Zbl 0889.76071