Linear, decoupled, second-order and structure-preserving scheme for Carreau fluid equations coupled with steric Poisson-Nernst-Planck model

Abstract

In this paper, to study ionic steric effects, we present a linear, decoupled, second-order accurate in time and structure-preserving scheme with finite element approximations for Carreau fluid equations coupled with steric Poisson-Nernst-Planck (SPNP) model. The logarithmic transformation for the ion concentration is used to preserve positivity property. To deal with the nonlinear coupling terms in fluid equation, a nonlocal auxiliary variable with respect to the free energy of SPNP equations and its associated ordinary differential equation are introduced. The obtained system is equivalent to the original system. The fully discrete scheme is proved to be mass conservative, positivity-preserving for ion concentration and energy dissipative at discrete level. Some numerical simulations are provided to demonstrate its stability and accuracy. Moreover, the ionic steric effects are numerically investigated.