Free-energy-dissipative schemes for the Oldroyd-B model
CERMICS, École Nationale des Ponts et Chaussées (ParisTech/Université Paris-Est), 6 & 8 avenue Blaise Pascal,
Cité Descartes, 77455 Marne-la-Vallée Cedex 2, France. firstname.lastname@example.org; email@example.com;
2 MICMAC team-project, INRIA, Domaine de Voluceau, BP. 105, Rocquencourt, 78153 Le Chesnay Cedex, France.
3 Institute of Mathematics, The Hebrew University, Jerusalem 91904, Israel.
Revised: 16 September 2008
In this article, we analyze the stability of various numerical schemes for differential models of viscoelastic fluids. More precisely, we consider the prototypical Oldroyd-B model, for which a free energy dissipation holds, and we show under which assumptions such a dissipation is also satisfied for the numerical scheme. Among the numerical schemes we analyze, we consider some discretizations based on the log-formulation of the Oldroyd-B system proposed by Fattal and Kupferman in [J. Non-Newtonian Fluid Mech. 123 (2004) 281–285], for which solutions in some benchmark problems have been obtained beyond the limiting Weissenberg numbers for the standard scheme (see [Hulsen et al. J. Non-Newtonian Fluid Mech. 127 (2005) 27–39]). Our analysis gives some tracks to understand these numerical observations.
Mathematics Subject Classification: 65M12 / 76M10 / 35B45 / 76A10 / 35B35
Key words: Viscoelastic fluids / Weissenberg number / stability / entropy / finite elements methods / discontinuous Galerkin method / characteristic method.
© EDP Sciences, SMAI, 2009