Energy-based operator splitting approach for the time discretization of coupled systems of partial and ordinary differential equations for fluid flows: The Stokes case

JOURNAL OF COMPUTATIONAL PHYSICS

卷:364

Article

Carichino, Lucia¹  Guidoboni, Giovanna²  Szopos, Marcela³ 

[1] Worcester Polytech Inst, Dept Math Sci, 100 Inst Rd, Worcester, MA 01609 USA

[2] Univ Missouri, Dept Elect Engn & Comp Sci, 201 Naka Hall, Columbia, MO 65211 USA

[3] Univ Strasbourg, CNRS, IRMA UMR 7501, F-67000 Strasbourg, France

关键词: Multiscale fluid flow;    Operator splitting;    Partial and ordinary differential equations;    Blood flow simulations;   

DOI  :  10.1016/j.jcp.2018.02.030

来源: Elsevier

【 摘 要 】

The goal of this work is to develop a novel splitting approach for the numerical solution of multiscale problems involving the coupling between Stokes equations and ODE systems, as often encountered in blood flow modeling applications. The proposed algorithm is based on a semi-discretization in time based on operator splitting, whose design is guided by the rationale of ensuring that the physical energy balance is maintained at the discrete level. As a result, unconditional stability with respect to the time step choice is ensured by the implicit treatment of interface conditions within the Stokes substeps, whereas the coupling between Stokes and ODE substeps is enforced via appropriate initial conditions for each substep. Notably, unconditional stability is attained without the need of subiterating between Stokes and ODE substeps. Stability and convergence properties of the proposed algorithm are tested on three specific examples for which analytical solutions are derived. (C) 2018 Elsevier Inc. All rights reserved.

【 授权许可】

Free