IMA Journal of Numerical Analysis Advance Access published online on September 18, 2008
IMA Journal of Numerical Analysis, doi:10.1093/imanum/drn022
| ||||||||||||||||||||||||||||||||||||||||||||||||||
Numerical approximation of corotational dumbbell models for dilute polymers
Department of Mathematics, Imperial College London, London SW7 2AZ, UK

Oxford University Computing Laboratory, Numerical Analysis Group, University of Oxford, Parks Road, Oxford OX1 3QD, UK
Email: endre.suli{at}comlab.ox.ac.uk
Received on 28 July 2007. Revised on 4 January 2008.
| Abstract |
|---|
We construct a general family of Galerkin methods for the numerical approximation of weak solutions to a bead-spring model that arises from the kinetic theory of dilute solutions of polymeric liquids with noninteracting polymer chains. The model consists of the unsteady incompressible Navier–Stokes equations, for the velocity and the pressure of the fluid, with an elastic extra-stress tensor as right-hand side in the momentum equation. The extra-stress tensor stems from the random movement of the polymer chains and is defined through the associated probability density function satisfying a Fokker–Planck-type parabolic equation. We focus on finitely extensible nonlinear elastic-type dumbbell models. In the case of a corotational drag term, we perform a rigorous passage to the limit as the spatial and temporal discretization parameters tend to zero and show that a (sub)sequence of numerical solutions converges to a weak solution of this coupled Navier–Stokes–Fokker–Planck system.
Key Words: Galerkin methods; polymeric flow models; existence of weak solutions; Navier–Stokes equations; Fokker–Planck equations; FENE