In this paper we propose a “Discrete Duality Finite Volume” method (DDFV for short) for the diffuse interface modelling of incompressible two-phase flows. This numerical method is, conservative, robust and is able to handle general geometries and meshes. The model we study couples the Cahn−Hilliard equation and the unsteady Stokes equation and is endowed with particular nonlinear boundary conditions called dynamic boundary conditions. To implement the scheme for this model we have to derive new discrete consistent DDFV operators that allows an energy stable coupling between both discrete equations. We are thus able to obtain the existence of a family of solutions satisfying a suitable energy inequality, even in the case where a first order time-splitting method between the two subsystems is used. We illustrate various properties of such a model with some numerical results.
Mots clés : Cahn–Hilliard/Stokes model, dynamic boundary conditions, contact angle dynamics, finite volume method
@article{M2AN_2017__51_5_1691_0, author = {Boyer, Franck and Nabet, Flore}, title = {A {DDFV} method for a {Cahn\ensuremath{-}Hilliard/Stokes} phase field model with dynamic boundary conditions}, journal = {ESAIM: Mathematical Modelling and Numerical Analysis }, pages = {1691--1731}, publisher = {EDP-Sciences}, volume = {51}, number = {5}, year = {2017}, doi = {10.1051/m2an/2016073}, mrnumber = {3731546}, zbl = {1391.35196}, language = {en}, url = {http://www.numdam.org/articles/10.1051/m2an/2016073/} }
TY - JOUR AU - Boyer, Franck AU - Nabet, Flore TI - A DDFV method for a Cahn−Hilliard/Stokes phase field model with dynamic boundary conditions JO - ESAIM: Mathematical Modelling and Numerical Analysis PY - 2017 SP - 1691 EP - 1731 VL - 51 IS - 5 PB - EDP-Sciences UR - http://www.numdam.org/articles/10.1051/m2an/2016073/ DO - 10.1051/m2an/2016073 LA - en ID - M2AN_2017__51_5_1691_0 ER -
%0 Journal Article %A Boyer, Franck %A Nabet, Flore %T A DDFV method for a Cahn−Hilliard/Stokes phase field model with dynamic boundary conditions %J ESAIM: Mathematical Modelling and Numerical Analysis %D 2017 %P 1691-1731 %V 51 %N 5 %I EDP-Sciences %U http://www.numdam.org/articles/10.1051/m2an/2016073/ %R 10.1051/m2an/2016073 %G en %F M2AN_2017__51_5_1691_0
Boyer, Franck; Nabet, Flore. A DDFV method for a Cahn−Hilliard/Stokes phase field model with dynamic boundary conditions. ESAIM: Mathematical Modelling and Numerical Analysis , Tome 51 (2017) no. 5, pp. 1691-1731. doi : 10.1051/m2an/2016073. http://www.numdam.org/articles/10.1051/m2an/2016073/
On a diffuse interface model for two-phase flows of viscous, incompressible fluids with matched densities. Arch. Ration. Mech. Anal. 194 (2009) 463–506. | DOI | MR | Zbl
,On a diffuse interface model for a two-phase flow of compressible viscous fluids. Indiana Univ. Math. J. 57 (2008) 659–698. | DOI | MR | Zbl
and ,Thermodynamically consistent, frame indifferent diffuse interface models for incompressible two-phase flows with different densities. Math. Models Methods Appl. Sci. 22 (2012) 1150013, 40. | DOI | MR | Zbl
, and ,On 3D DDFV discretization of gradient and divergence operators. I. meshing, operators and discrete duality. IMA J. Numer. Anal. 32 (2012) 1574–1603. | DOI | MR | Zbl
, , and ,Discrete duality finite volume schemes for Leray-Lions-type elliptic problems on general 2D meshes. Numer. Methods Partial Differ. Equ. 23 (2007) 145–195. | DOI | MR | Zbl
, and ,On discrete functional inequalities for some finite volume schemes. IMA J. Numer. Anal. 35 (2015) 1125–1149. | DOI | MR | Zbl
, and .A theoretical and numerical model for the study of incompressible mixture flows. Comput. Fluids 31 (2002) 41–68. | DOI | Zbl
,Nonoverlapping Schwarz algorithm for solving two-dimensional m-DDFV schemes. IMA J. Numer. Anal. 30 (2010) 1062–1100. | DOI | MR | Zbl
, and ,Inf-sup stability of the discrete duality finite volume method for the 2D stokes problem. Math. Comput. 84 (2015) 2705–2742. | DOI | MR | Zbl
, and ,Cahn−Hilliard/Navier−Stokes model for the simulation of three-phase flows. Trans. Porous Media 82 (2010) 463–483. | DOI | MR
, , , and ,Numerical schemes for a three component Cahn−Hilliard model. ESAIM: M2AN 45 (2011) 697–738. | DOI | Numdam | MR | Zbl
and ,Dissipation in rapid dynamic wetting. J. Fluid Mech. 682 (2011) 213–240. | DOI | Zbl
, and ,A numerical analysis of the Cahn−Hilliard equation with dynamic boundary conditions. Discrete Contin. Dyn. Syst. 27 (2010) 1511–1533. | DOI | MR | Zbl
, and ,Convergence to steady state of solutions of the Cahn−Hilliard and Caginalp equations with dynamic boundary conditions. Math. Nachr. 279 (2006) 1448–1462. | DOI | MR | Zbl
, ,A finite volume method for the Laplace equation on almost arbitrary two-dimensional grids. M2AN: M2AN 39 (2005) 1203–1249. | DOI | Numdam | MR | Zbl
and .On imposing dynamic contact-angle boundary conditions for wall-bounded liquid-gas flows. Comput. Methods Appl. Mech. Engrg. 247/248 (2012) 179–200. | DOI | MR | Zbl
,An outflow boundary condition and algorithm for incompressible two-phase flows with phase field approach. J. Comput. Phys. 266 (2014) 47–73. | DOI | MR | Zbl
,R. Eymard, T. Gallouët and R. Herbin, Finite volume methods. In Handbook of numerical analysis, Vol. VII. Edited by Ph. Ciarlet and J.L. Lions. North-Holland, Amsterdam (2000) 715–1022. | MR | Zbl
R. Eymard, G. Henry, R. Herbin, F. Hubert, R. Klofkorn and G. Manzini, 3d benchmark on discretization schemes for anisotropic diffusion problems on general grids. In Proc. of Finite Volumes for Complex Applications, Vol. VI. Springer (2011) 895–930. | Zbl
Novel surface modes in spinodal decomposition. Phys. Rev. Lett. 79 (1997) 893–896. | DOI
, and ,Diverging time and length scales of spinodal decomposition modes in thin films. EPL Europhys. Lett. 42 (1998) 49–54. | DOI
, and ,T. Goudon and S. Krell, A DDFV Scheme for Incompressible Navier−Stokes equations with variable density. In Proc. of Finite Volumes for Complex Applications VII, edited by J. Fuhrmann, M. Ohlberger and C. Rohde. In Vol. 77 and 78. Springer Proceedings in Mathematics and Statistics. Springer, Berlin, Allemagne (2014). | MR
R. Herbin and F. Hubert, Benchmark on discretization schemes for anisotropic diffusion problems on general grids. In Proc. of Finite Volumes for Complex Applications V, Aussois, France. Edited by R. Eymard and J.M. Herard. Hermès (2008). | MR
Quantitative benchmark computations of two-dimensional bubble dynamics. Int. J. Numer. Methods Fluids 60 (2009) 1259–1288. | DOI | MR | Zbl
, , , , , and ,Calculation of two-phase Navier−Stokes flows using phase-field modeling. J. Comput. Phys. 155 (1999) 96–127. | DOI | MR | Zbl
,Contact-line dynamics of a diffuse fluid interface. J. Fluid Mech. 402 (2000) 57–88. | DOI | MR | Zbl
,Finite element approximation of a Cahn–Hilliard–Navier−Stokes system. Interfaces Free Bound. 10 (2008) 15–43. | DOI | MR | Zbl
, and ,Phase separation in confined geometries: Solving the Cahn–Hilliard equation with generic boundary conditions. Comput. Phys. Commun. 133 (2001) 139–157. | DOI | MR | Zbl
, , , , , and ,Stabilized DDFV schemes for Stokes problem with variable viscosity on general 2D meshes. Numer. Methods Partial Differ. Equ. 27 (2011) 1666–1706. | DOI | MR | Zbl
,The discrete duality finite volume method for stokes equations on three-dimensional polyhedral meshes. SIAM J. Numer. Anal. 50 (2012) 808–837. | DOI | MR | Zbl
and ,A phase field model for the mixture of two incompressible fluids and its approximation by a Fourier-spectral method. Physica D 179 (2003) 211–228. | DOI | MR | Zbl
and ,On numerical schemes for phase-field models for electrowetting with electrolyte solutions. Proc. Appl. Math. Mech. 15 (2015) 715–718. | DOI
,An adaptive pressure correction method without spurious velocities for diffuse-interface models of incompressible flows. J. Comput. Phys. 236 (2013) 143–156. | DOI | MR | Zbl
,An unconditionally stable uncoupled scheme for a triphasic Cahn−Hilliard/Navier−Stokes model. Numer. Methods Partial Differ. Equ. 29 (2013) 584–618. | DOI | MR | Zbl
,Exponential attractors for the Cahn–Hilliard equation with dynamic boundary conditions. Math. Methods Appl. Sci. 28 (2005) 709–735. | DOI | MR | Zbl
and ,F. Nabet, Schémas volumes finis pour des problèmes multiphasiques. Ph.D. thesis, Aix-Marseille université (2014).
Convergence of a finite-volume scheme for the cahn–hilliard equation with dynamic boundary conditions. IMA J. Numer. Anal. 36 (2015) 1898–1942. | DOI | MR | Zbl
,F. Nabet, An error estimate for a finite-volume scheme for the Cahn−Hilliard equation with dynamic boundary conditions (2016). Preprint (2016). | HAL
Maximal regularity and asymptotic behavior of solutions for the Cahn−Hilliard equation with dynamic boundary conditions. Ann. Mat. Pura Appl. 185 (2006) 627–648. | DOI | MR | Zbl
, and ,A variational approach to moving contact line hydrodynamics. J. Fluid Mech. 564 (2006) 333. | DOI | MR | Zbl
, and ,The Cahn−Hilliard equation with dynamic boundary conditions. Adv. Differ. Equ. 8 (2003) 83–110. | MR | Zbl
and ,A diffuse interface fractional time-stepping technique for incompressible two-phase flows with moving contact lines. ESAIM: M2AN 47 (2013) 743–769, 5. | DOI | Numdam | MR | Zbl
,Efficient energy stable numerical schemes for a phase field moving contact line model. J. Comput. Phys. 284 (2015) 617–630. | DOI | MR | Zbl
, and ,X-.Moving contact line on chemically patterned surfaces. J. Fluid Mech. 605 (2008) 59–78. | DOI | MR | Zbl
, and ,Convergence to equilibrium for the Cahn−Hilliard equation with dynamic boundary conditions. J. Differ. Equ. 204 (2004) 511–531. | DOI | MR | Zbl
and ,Cité par Sources :