An extended volume of fluid method and its application to single bubbles rising in a viscoelastic liquid

Matthias Niethammer, Günter Brenn, Holger Marschall, Dieter Bothe

Publikation: Beitrag in einer FachzeitschriftArtikelForschungBegutachtung

Abstract

An extended volume of fluid method is developed for two-phase direct numerical simulations of systems with one viscoelastic and one Newtonian phase. A complete set of governing equations is derived by conditional volume averaging the local instantaneous bulk equations and interface jump conditions. The homogeneous mixture model is applied for the closure of the volume-averaged equations. An additional interfacial stress term arises in this volume-averaged formulation which requires special treatment in the finite-volume discretization on a general unstructured mesh. A novel numerical scheme is proposed for the second-order accurate finite-volume discretization of the interface stress term. We demonstrate that this scheme allows for a consistent treatment of the interface stress and the surface tension force in the pressure equation of the segregated solution approach. Because of the high Weissenberg number problem, an appropriate stabilization approach is applied to the constitutive equation of the viscoelastic phase to increase the robustness of the method at higher fluid elasticity. Direct numerical simulations of the transient motion of a bubble rising in a quiescent viscoelastic fluid are performed for the purpose of experimental code validation. The well-known jump discontinuity in the terminal bubble rise velocity when the bubble volume exceeds a critical value is captured by the method. The formulation of the interfacial stress together with the novel scheme for its discretization is found crucial for the quantitatively correct prediction of the jump discontinuity in the terminal bubble rise velocity.

Originalspracheenglisch
Seiten (von - bis)326-355
Seitenumfang30
FachzeitschriftJournal of computational physics
Jahrgang387
DOIs
PublikationsstatusVeröffentlicht - 2019

Fingerprint

Bubbles (in fluids)
bubbles
Fluids
fluids
Direct numerical simulation
Liquids
liquids
direct numerical simulation
Constitutive equations
discontinuity
Surface tension
Elasticity
Stabilization
formulations
constitutive equations
closures
mesh
interfacial tension
elastic properties
stabilization

Schlagwörter

    ASJC Scopus subject areas

    • !!Computer Science Applications
    • !!Physics and Astronomy (miscellaneous)

    Fields of Expertise

    • Sonstiges

    Dies zitieren

    An extended volume of fluid method and its application to single bubbles rising in a viscoelastic liquid. / Niethammer, Matthias; Brenn, Günter; Marschall, Holger; Bothe, Dieter.

    in: Journal of computational physics, Jahrgang 387, 2019, S. 326-355.

    Publikation: Beitrag in einer FachzeitschriftArtikelForschungBegutachtung

    @article{b9c1638fcd0049e08a0c943d38d8ee29,
    title = "An extended volume of fluid method and its application to single bubbles rising in a viscoelastic liquid",
    abstract = "An extended volume of fluid method is developed for two-phase direct numerical simulations of systems with one viscoelastic and one Newtonian phase. A complete set of governing equations is derived by conditional volume averaging the local instantaneous bulk equations and interface jump conditions. The homogeneous mixture model is applied for the closure of the volume-averaged equations. An additional interfacial stress term arises in this volume-averaged formulation which requires special treatment in the finite-volume discretization on a general unstructured mesh. A novel numerical scheme is proposed for the second-order accurate finite-volume discretization of the interface stress term. We demonstrate that this scheme allows for a consistent treatment of the interface stress and the surface tension force in the pressure equation of the segregated solution approach. Because of the high Weissenberg number problem, an appropriate stabilization approach is applied to the constitutive equation of the viscoelastic phase to increase the robustness of the method at higher fluid elasticity. Direct numerical simulations of the transient motion of a bubble rising in a quiescent viscoelastic fluid are performed for the purpose of experimental code validation. The well-known jump discontinuity in the terminal bubble rise velocity when the bubble volume exceeds a critical value is captured by the method. The formulation of the interfacial stress together with the novel scheme for its discretization is found crucial for the quantitatively correct prediction of the jump discontinuity in the terminal bubble rise velocity.",
    keywords = "Extended volume of fluid method, Negative wake, Rising bubble, Velocity jump discontinuity, Viscoelastic liquid",
    author = "Matthias Niethammer and G{\"u}nter Brenn and Holger Marschall and Dieter Bothe",
    year = "2019",
    doi = "10.1016/j.jcp.2019.02.021",
    language = "English",
    volume = "387",
    pages = "326--355",
    journal = "Journal of computational physics",
    issn = "0021-9991",
    publisher = "Academic Press",

    }

    TY - JOUR

    T1 - An extended volume of fluid method and its application to single bubbles rising in a viscoelastic liquid

    AU - Niethammer, Matthias

    AU - Brenn, Günter

    AU - Marschall, Holger

    AU - Bothe, Dieter

    PY - 2019

    Y1 - 2019

    N2 - An extended volume of fluid method is developed for two-phase direct numerical simulations of systems with one viscoelastic and one Newtonian phase. A complete set of governing equations is derived by conditional volume averaging the local instantaneous bulk equations and interface jump conditions. The homogeneous mixture model is applied for the closure of the volume-averaged equations. An additional interfacial stress term arises in this volume-averaged formulation which requires special treatment in the finite-volume discretization on a general unstructured mesh. A novel numerical scheme is proposed for the second-order accurate finite-volume discretization of the interface stress term. We demonstrate that this scheme allows for a consistent treatment of the interface stress and the surface tension force in the pressure equation of the segregated solution approach. Because of the high Weissenberg number problem, an appropriate stabilization approach is applied to the constitutive equation of the viscoelastic phase to increase the robustness of the method at higher fluid elasticity. Direct numerical simulations of the transient motion of a bubble rising in a quiescent viscoelastic fluid are performed for the purpose of experimental code validation. The well-known jump discontinuity in the terminal bubble rise velocity when the bubble volume exceeds a critical value is captured by the method. The formulation of the interfacial stress together with the novel scheme for its discretization is found crucial for the quantitatively correct prediction of the jump discontinuity in the terminal bubble rise velocity.

    AB - An extended volume of fluid method is developed for two-phase direct numerical simulations of systems with one viscoelastic and one Newtonian phase. A complete set of governing equations is derived by conditional volume averaging the local instantaneous bulk equations and interface jump conditions. The homogeneous mixture model is applied for the closure of the volume-averaged equations. An additional interfacial stress term arises in this volume-averaged formulation which requires special treatment in the finite-volume discretization on a general unstructured mesh. A novel numerical scheme is proposed for the second-order accurate finite-volume discretization of the interface stress term. We demonstrate that this scheme allows for a consistent treatment of the interface stress and the surface tension force in the pressure equation of the segregated solution approach. Because of the high Weissenberg number problem, an appropriate stabilization approach is applied to the constitutive equation of the viscoelastic phase to increase the robustness of the method at higher fluid elasticity. Direct numerical simulations of the transient motion of a bubble rising in a quiescent viscoelastic fluid are performed for the purpose of experimental code validation. The well-known jump discontinuity in the terminal bubble rise velocity when the bubble volume exceeds a critical value is captured by the method. The formulation of the interfacial stress together with the novel scheme for its discretization is found crucial for the quantitatively correct prediction of the jump discontinuity in the terminal bubble rise velocity.

    KW - Extended volume of fluid method

    KW - Negative wake

    KW - Rising bubble

    KW - Velocity jump discontinuity

    KW - Viscoelastic liquid

    UR - http://www.scopus.com/inward/record.url?scp=85063061986&partnerID=8YFLogxK

    U2 - 10.1016/j.jcp.2019.02.021

    DO - 10.1016/j.jcp.2019.02.021

    M3 - Article

    VL - 387

    SP - 326

    EP - 355

    JO - Journal of computational physics

    JF - Journal of computational physics

    SN - 0021-9991

    ER -