Accessible Requires Authentication Published by De Gruyter December 31, 2019

On convergent schemes for a two-phase Oldroyd-B type model with variable polymer density

Oliver Sieber

Abstract

The paper is concerned with a diffuse-interface model that describes two-phase flow of dilute polymeric solutions with a variable particle density. The additional stresses, which arise by elongations of the polymers caused by deformations of the fluid, are described by Kramers stress tensor. The evolution of Kramers stress tensor is modeled by an Oldroyd-B type equation that is coupled to a Navier–Stokes type equation, a Cahn–Hilliard type equation, and a parabolic equation for the particle density. We present a regularized finite element approximation of this model, prove that our scheme is energy stable and that there exist discrete solutions to it. Furthermore, in the case of equal mass densities and two space dimensions, we are able to pass to the limit rigorously as the regularization parameters and the spatial and temporal discretization parameters tend towards zero and prove that a subsequence of discrete solutions converges to a global-in-time weak solution to the unregularized coupled system. To the best of our knowledge, this is the first existence result for a two-phase flow model of viscoelastic fluids with an Oldroyd-B type equation. Additionally, we show that our finite element scheme is fully practical and we present numerical simulations.

JEL Classification: 35Q35; 76A05; 76D03; 76T99; 35Q30; 65M12

Acknowledgment

This work has been supported by the RTG 2339 ‘Interfaces, Complex Structures, and Singular Limits’ of the German Science Foundation (DFG). The author is grateful to Günther Grün for many valuable discussions.

References

[1] H. Abels, D. Depner, and H. Garcke, Existence of weak solutions for a diffuse interface model for two-phase flows of incompressible fluids with different densities, J. Math. Fluid Mech., 15 (2013), No. 3, 453–480. Search in Google Scholar

[2] H. Abels, H. Garcke, and G. Grün, Thermodynamically consistent, frame indifferent diffusive-interface models for incompressible two-phase flows with different densities, Math. Models Methods Appl. Sci., 22 (2012), No. 3, 1150013. Search in Google Scholar

[3] F. Armero and J.C. Simo, Formulation of a new class of fractional-step methods for the incompressible MHD equations that retains the long-term dissipativity of the continuum dynamical system, Fields Inst. Commun., 10 (1996), 1–24. Search in Google Scholar

[4] J. W. Barrett and S. Boyaval, Existence and approximation of a (regularized) Oldroyd-B model, Math. Models Methods Appl. Sci., 21 (2011), No. 9, 1783–1837. Search in Google Scholar

[5] J. W. Barrett and S. Boyaval, Finite element approximation of the FENE-P model, IMA J. Numer. Anal., 38 (2018), No. 4, 1599–1660. Search in Google Scholar

[6] J. W. Barrett, Y. Lu, and E. Süli, Existence of large-data finite-energy global weak solutions to a compressible Oldroyd-B model, ArXiv:1608.04229, 2016. Search in Google Scholar

[7] J. W. Barrett and E. Süli, Existence of global weak solutions to some regularized kinetic models for dilute polymers, Multiscale Model. Simul., 6 (2007), 505–546. Search in Google Scholar

[8] J. W. Barrett and E. Süli, Finite element approximation of kinetic dilute polymer models with microscopic cut-off, ESAIM: Math. Modelling Numer. Analysis, 45 (2011), No. 1, 39–89. Search in Google Scholar

[9] S. C. Brenner and L. R. Scott, The Mathematical Theory of Finite Element Methods, Texts in Applied Mathematics, 2007. Search in Google Scholar

[10] Q. Du and X. Wang, Convergence of numerical approximations to a phase field bending elasticity model of membrane deformations, Int. J. Numer. Anal. Mod., 4 (2007), No. 3-4, 441–459. Search in Google Scholar

[11] H. Garcke, M. Hinze, and C. Kahle, A stable and linear time discretization for a thermodynamically consistent model for two-phase incompressible flow, Appl. Numer. Math., 99 (2016), 151–171. Search in Google Scholar

[12] V. Girault and P. A. Raviart, Finite Element Methods for NavierStokes Equations: Theory and Algorithms, Springer Series in Computational Mathematics, Springer-Verlag, 1986. Search in Google Scholar

[13] G. Grün, On convergent schemes for diffuse interface models for two-phase flow of incompressible fluids with general mass densities, SIAM J. Numer. Analysis, 51 (2013), No. 6, 3036–3061. Search in Google Scholar

[14] G. Grün, F. Guillén-González, and S. Metzger, On fully decoupled, convergent schemes for diffuse interface models for two-phase flow with general mass densities, Commun. Comp. Phys., 19 (2016), No. 5, 1473–1502. Search in Google Scholar

[15] G. Grün and F. Klingbeil, Two-phase flow with mass density contrast: Stable schemes for a thermodynamic consistent and frame-indifferent diffuse-interface model, J. Comp. Phys., 257 (2014), 708–725. Search in Google Scholar

[16] G. Grün and S. Metzger, On micro-macro-models for two-phase flow with dilute polymeric solutions — modeling and analysis, Math. Models Meth. Appl. Sci., 26 (2016), No. 5, 823–866. Search in Google Scholar

[17] F. Guillén-González and G. Tierra, Splitting schemes for a Navier–Stokes–Cahn–Hilliard model for two fluids with different densities, J. Comp. Math., 32 (2014), No. 6, 643–664 Search in Google Scholar

[18] D. Harvie, M. Davidson, and J. Cooper-White, Simulating the deformation of Newtonian and non-Newtonian drops through a micro-fluidic contraction, In: 15th Australasian Fluid Mechanics Conference, 2004. Search in Google Scholar

[19] D. Hu and T. Lelièvre, New entropy estimates for the oldroyd-b model and related models, Commun. Math. Sci. 5 (2007), No. 4, 909–916. Search in Google Scholar

[20] D. Izbassarov and M. Muradoglu, A front-tracking method for computational modeling of viscoelastic two-phase flow systems, J. Non-Newtonian Fluid Mech., 223 (2015), 122–140. Search in Google Scholar

[21] F. Klingbeil, On the numerics of diffuse-interface models for two-phase flow with species transport, Ph.D. Thesis, Friedrich-Alexander-Universität Erlangen-Nürnberg, 2014. Search in Google Scholar

[22] PL Lions and N. Masmoudi, Global solutions for some Oldroyd models of non-Newtonian flows, Chinese Annals Math., 21 (2000), No. 2, 131–146. Search in Google Scholar

[23] C. McIlroy, O. Harlen, and N. Morrison, Modelling the jetting of dilute polymer solutions in drop-on-demand inkjet printing, J. Non-Newtonian Fluid Mech., 201 (2013), 17–28. Search in Google Scholar

[24] S. Metzger, Diffuse Interface Models for Complex Flow Scenarios: Modeling, Analysis, and Simulation, Dr. Hut Verlag, 2017. Search in Google Scholar

[25] S. Metzger, On convergent schemes for two-phase flow of dilute polymeric solutions, ESAIM: Math. Modelling Numer. Analysis52 (2018), No. 6, 2357–2408. Search in Google Scholar

[26] S. Minjeaud, An unconditionally stable uncoupled scheme for a triphasic Cahn–Hilliard/Navier–Stokes model, Numer. Methods Partial Diff. Equ., 29 (2013), No. 2, 584–618. Search in Google Scholar

[27] S. B. Pillapakkam and P. Singh, A level-set method for computing solutions to viscoelastic two-phase flow, J. Comp. Phys., 174 (2001), No. 2, 552–578. Search in Google Scholar

[28] M. Ruzicka, Nichtlineare Funktionalanalysis: Eine Einführung, Springer-Lehrbuch Masterclass, Springer, Berlin–Heidelberg, 2004. Search in Google Scholar

[29] X. Yang, J. J. Feng, C. Liu, and J. Shen, Numerical simulations of jet pinching-off and drop formation using an energetic variational phase-field method, J. Comp. Phys., 218 (2006), No. 1, 417–428. Search in Google Scholar

[30] P. Yue, J. J. Feng, C. Liu, and J. Shen, A diffuse-interface method for simulating two-phase flows of complex fluids, J. Fluid Mech., 515 (2004), 293–317. Search in Google Scholar

[31] Y. Zhang, H. Wang, and T. Tang, Simulating two-phase viscoelastic flows using moving finite element methods, Commun. Comp. Phys., 7 (2010), 333–349. Search in Google Scholar

Received: 2019-01-31
Revised: 2019-10-14
Accepted: 2019-10-14
Published Online: 2019-12-31
Published in Print: 2020-06-25

© 2020 Walter de Gruyter GmbH, Berlin/Boston