# A heterogeneous alternating-direction method for a micro-macro dilute polymeric fluid model

ESAIM: Mathematical Modelling and Numerical Analysis (2009)

- Volume: 43, Issue: 6, page 1117-1156
- ISSN: 0764-583X

## Access Full Article

top## Abstract

top## How to cite

topKnezevic, David J., and Süli, Endre. "A heterogeneous alternating-direction method for a micro-macro dilute polymeric fluid model." ESAIM: Mathematical Modelling and Numerical Analysis 43.6 (2009): 1117-1156. <http://eudml.org/doc/250613>.

@article{Knezevic2009,

abstract = {
We examine a heterogeneous alternating-direction method for the approximate solution of the FENE Fokker–Planck equation from polymer fluid dynamics and we use this method to solve a coupled (macro-micro) Navier–Stokes–Fokker–Planck system for dilute polymeric fluids. In this context the Fokker–Planck equation is posed on a high-dimensional domain and is therefore challenging from a computational point of view. The heterogeneous alternating-direction scheme combines a spectral Galerkin method for the Fokker–Planck equation in configuration space with a finite element method in physical space to obtain a scheme for the high-dimensional Fokker–Planck equation. Alternating-direction methods have been considered previously in the literature for this problem (e.g. in the work of Lozinski, Chauvière and collaborators [J. Non-Newtonian Fluid Mech.122 (2004) 201–214; Comput. Fluids33 (2004) 687–696; CRM Proc. Lect. Notes41 (2007) 73–89; Ph.D. Thesis (2003); J. Computat. Phys.189 (2003) 607–625]), but this approach has not previously been subject to rigorous numerical analysis. The numerical methods we develop are fully-practical, and we present a range of numerical results demonstrating their accuracy and efficiency. We also examine an advantageous superconvergence property related to the polymeric extra-stress tensor. The heterogeneous alternating-direction method is well suited to implementation on a parallel computer, and we exploit this fact to make large-scale computations feasible.
},

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

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

keywords = {Multiscale modelling; kinetic models; dilute polymers; alternating-direction methods; spectral methods; finite element methods; high-performance computing.; multiscale modelling; high-performance computing},

language = {eng},

month = {8},

number = {6},

pages = {1117-1156},

publisher = {EDP Sciences},

title = {A heterogeneous alternating-direction method for a micro-macro dilute polymeric fluid model},

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

volume = {43},

year = {2009},

}

TY - JOUR

AU - Knezevic, David J.

AU - Süli, Endre

TI - A heterogeneous alternating-direction method for a micro-macro dilute polymeric fluid model

JO - ESAIM: Mathematical Modelling and Numerical Analysis

DA - 2009/8//

PB - EDP Sciences

VL - 43

IS - 6

SP - 1117

EP - 1156

AB -
We examine a heterogeneous alternating-direction method for the approximate solution of the FENE Fokker–Planck equation from polymer fluid dynamics and we use this method to solve a coupled (macro-micro) Navier–Stokes–Fokker–Planck system for dilute polymeric fluids. In this context the Fokker–Planck equation is posed on a high-dimensional domain and is therefore challenging from a computational point of view. The heterogeneous alternating-direction scheme combines a spectral Galerkin method for the Fokker–Planck equation in configuration space with a finite element method in physical space to obtain a scheme for the high-dimensional Fokker–Planck equation. Alternating-direction methods have been considered previously in the literature for this problem (e.g. in the work of Lozinski, Chauvière and collaborators [J. Non-Newtonian Fluid Mech.122 (2004) 201–214; Comput. Fluids33 (2004) 687–696; CRM Proc. Lect. Notes41 (2007) 73–89; Ph.D. Thesis (2003); J. Computat. Phys.189 (2003) 607–625]), but this approach has not previously been subject to rigorous numerical analysis. The numerical methods we develop are fully-practical, and we present a range of numerical results demonstrating their accuracy and efficiency. We also examine an advantageous superconvergence property related to the polymeric extra-stress tensor. The heterogeneous alternating-direction method is well suited to implementation on a parallel computer, and we exploit this fact to make large-scale computations feasible.

LA - eng

KW - Multiscale modelling; kinetic models; dilute polymers; alternating-direction methods; spectral methods; finite element methods; high-performance computing.; multiscale modelling; high-performance computing

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

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.
- S. Balay, K. Buschelman, V. Eijkhout, W.D. Gropp, D. Kaushik, M.G. Knepley, L.C. McInnes, B.F. Smith and H. Zhang, PETSc users manual. Tech. Rep. ANL-95/11 – Revision 2.1.5, Argonne National Laboratory (2004).
- J.W. Barrett and E. Süli, Existence of global weak solutions to dumbbell models for dilute polymers with microscopic cut-off. Math. Models Methods Appl. Sci.18 (2008) 935–971.
- B. Bialecki and R. Fernandes, An orthogonal spline collocation alternating direction implicit Crank-Nicolson method for linear parabolic problems on rectangles. SIAM J. Numer. Anal.36 (1999) 1414–1434.
- P.B. Bochev, M.D. Gunzburger and J.N. Shadid, Stability of the SUPG finite element method for transient advection-diffusion problems. Comput. Methods Appl. Mech. Engrg.193 (2004) 2301–2323.
- S.C. Brenner and L.R. Scott, The Mathematical Theory of Finite Element Methods. Second Edn., Springer (2002).
- M. Celia and G. Pinder, An analysis of alternating-direction methods for parabolic equations. Numer. Methods Part. Differ. Equ.1 (1985) 57–70.
- M. Celia and G. Pinder, Generalized alternating-direction collocation methods for parabolic equations. I. Spatially varying coefficients. Numer. Methods Partial Differ. Equ.3 (1990) 193–214.
- 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.
- P. Clément, Approximation by finite element functions using local regularization. Rev. Française Automat. Informat. Recherche Opérationnelle Sér. RAIRO Anal. Numér.9 (1975) 77–84.
- P. Delaunay, A. Lozinski and R.G. Owens, Sparse tensor-product Fokker–Planck-based methods for nonlinear bead-spring chain models of dilute polymer solutions. CRM Proc. Lect. Notes41 (2007) 73–89.
- J. Douglas and T. Dupont, Alternating-direction Galerkin methods on rectangles. Numer. Solution Partial Differ. Equ. II (SYNSPADE 1970) (1971) 133–214.
- H. Eisen, W. Heinrichs and K. Witsch, Spectral collocation methods and polar coordinate singularities. J. Comput. Phys.96 (1991) 241–257.
- H. Elman, D. Silvester and A. Wathen, Finite elements and fast iterative solvers. Oxford Science Publications, UK (2005).
- C. Helzel and F. Otto, Multiscale simulations of suspensions of rod-like molecules. J. Comp. Phys.216 (2006) 52–75.
- W. Huang and B. Guo, Fully discrete Jacobi-spherical harmonic spectral method for Navier-Stokes equations. Appl. Math. Mech.29 (2008) 453–476 (English Ed.).
- B. Jourdain, T. Lelièvre and C. Le Bris, Existence of solution for a micro-macro model of polymeric fluid: the FENE model. J. Funct. Anal.209 (2004) 162–193.
- B.S. Kirk, J.W. Peterson, R.M. Stogner and G.F. Carey, libMesh: A C++ library for parallel adaptive mesh refinement/coarsening simulations. Eng. Comput.23 (2006) 237–254.
- D.J. Knezevic, Analysis and implementation of numerical methods for simulating dilute polymeric fluids. Ph.D. Thesis, University of Oxford, UK (2008), . URIhttp://www.comlab.ox.ac.uk/people/David.Knezevic
- D.J. Knezevic and E. Süli, Spectral Galerkin approximation of Fokker–Planck equations with unbounded drift. ESAIM: M2AN43 (2009) 445–485.
- A.N. Kolmogorov, Über die analytischen Methoden in der Wahrscheinlichkeitsrechnung. Math. Ann.104 (1931).
- T. Li and P. Zhang, Mathematical analysis of multi-scale models of complex fluids. Commun. Math. Sci.5 (2007) 1–51.
- C. Liu and H. Liu, Boundary conditions for the microscopic FENE models. SIAM J. Appl. Math.68 (2008) 1304–1315.
- A. Lozinski, Spectral methods for kinetic theory models of viscoelastic fluids. Ph.D. Thesis, École Polytechnique Fédérale de Lausanne, Suisse (2003).
- A. Lozinski and C. Chauvière, A fast solver for Fokker–Planck equation applied to viscoelastic flows calculation: 2D FENE model. J. Computat. Phys.189 (2003) 607–625.
- J.N. Lyness and D. Jespersen, Moderate degree symmetric quadrature rules for the triangle. J. Inst. Math. Appl.15 (1975) 19–32.
- 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).
- R.G. Owens and T.N. Phillips, Computational Rheology. Imperial College Press (2002).
- C. Schwab, E. Süli and R.A. Todor, Sparse finite element approximation of high-dimensional transport-dominated diffusion problems. ESAIM: M2AN42 (2008) 777–820.
- L.R. Scott and S. Zhang, Finite element interpolation of nonsmooth functions satisfying boundary conditions. Math. Comp.54 (1990) 483–493.
- 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.
- N.J. Walkington, Quadrature on simplices of arbitrary dimension. nw0z/publications/00-CNA-023/023abs/. URIhttp://www.math.cmu.edu/
- H.R. Warner, Kinetic theory and rheology of dilute suspensions of finitely extendible dumbbells. Ind. Eng. Chem. Fundamentals11 (1972) 379–387.
- H. Zhang and P. Zhang, Local existence for the FENE-dumbbell model of polymeric fluids. Arch. Ration. Mech. Anal.181 (2006) 373–400.

## 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
- 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.