Implicit-explicit schemes for incompressible flow problems with variable viscosity

Barrenechea, Gabriel R. and Castillo, Ernesto and Pacheco, Douglas R.Q. (2024) Implicit-explicit schemes for incompressible flow problems with variable viscosity. SIAM Journal on Scientific Computing. ISSN 1064-8275 (In Press)

[thumbnail of Barrenechea-etal-SIAMJSC-2024-Implicit-explicit-schemes-for-incompressible-flow-problems] Text. Filename: Barrenechea-etal-SIAMJSC-2024-Implicit-explicit-schemes-for-incompressible-flow-problems.pdf
Accepted Author Manuscript
Restricted to Repository staff only until 1 January 2099.

Download (1MB) | Request a copy

Abstract

This article investigates different implicit-explicit (IMEX) methods for incompressible flows with variable viscosity. The viscosity field may depend on space and time alone, or, for example, on velocity gradients. Unlike most previous works on IMEX schemes, which focus on the convective term, we propose treating also parts of the diffusive term explicitly, which can reduce the coupling between the velocity components. We present different IMEX alternatives for the variable-viscosity Navier-Stokes system, analysing their theoretical and algorithmic properties. Temporal stability is proven for all the methods presented, including monolithic and fractional-step variants. These results are unconditional, except for one of the fractional-step discretisations, whose stability is shown for time-step sizes under an upper bound that depends solely on the problem data. The key finding of this work is a class of IMEX schemes whose steps decouple the velocity components and are fully linearised (even if the viscosity depends nonlinearly on the velocity), without requiring any CFL condition for stability. Moreover, in the presence of Neumann boundaries, some of our formulations lead naturally to conditions involving normal pseudo-tractions. This generalises to the variable viscosity case what happens for the standard Laplacian form with constant viscosity. Our analysis is supported by a series of numerical experiments.