Scholarlink

<p>We demonstrate the use of computer algebra for facilitating the derivation of thin-film reduced-order models. We focus on the weighted-residual integral boundary layer (WRIBL) method, which has proven to be a very effective technique for developing reduced-order models by averaging the Navier-Stokes equations over the thin-gap direction. In particular, SymPy (the symbolic computing library in Python) is used to derive the core-annular WRIBL model of Dietze and Ruyer-Quil in (J Fluid Mech 762:60, 2015); the derivation is especially involved due to the inclusion of second-order terms, the presence of two hydrodynamically active phases, the enforcement of interfacial boundary conditions, and the cylindrical geometry. We show, using excerpts of code, how each step of the derivation can be broken down into substeps that are amenable to symbolic computation. To illustrate the application of the derived model, we solve it numerically using scientific computing libraries in Python, and briefly explore the dynamics of the Rayleigh-Plateau instability. The use of open-source computer algebra, in the manner described here, greatly eases the derivation of averaged models, thereby facilitating their use for the study of multiscale flows, as well as for computationally efficient prediction and optimization.</p>

错误：搜索内容不能为空，请输入英文关键词
错误：关键词超出字数限制，请精简
高级检索

Deriving thin-film averaged equations using computer algebra

  • 期刊论文

  • 发表日期：

    2025年9月10日

查看全文
  • Swarnaditya Hazra,
  • Jason R. Picardo

摘要

We demonstrate the use of computer algebra for facilitating the derivation of thin-film reduced-order models. We focus on the weighted-residual integral boundary layer (WRIBL) method, which has proven to be a very effective technique for developing reduced-order models by averaging the Navier-Stokes equations over the thin-gap direction. In particular, SymPy (the symbolic computing library in Python) is used to derive the core-annular WRIBL model of Dietze and Ruyer-Quil in (J Fluid Mech 762:60, 2015); the derivation is especially involved due to the inclusion of second-order terms, the presence of two hydrodynamically active phases, the enforcement of interfacial boundary conditions, and the cylindrical geometry. We show, using excerpts of code, how each step of the derivation can be broken down into substeps that are amenable to symbolic computation. To illustrate the application of the derived model, we solve it numerically using scientific computing libraries in Python, and briefly explore the dynamics of the Rayleigh-Plateau instability. The use of open-source computer algebra, in the manner described here, greatly eases the derivation of averaged models, thereby facilitating their use for the study of multiscale flows, as well as for computationally efficient prediction and optimization.