A decoupled Crank-Nicolson leap-frog scheme for the unsteady bioconvection flows problem with concentration dependent viscosity
Chenyang Li
TL;DR
This paper addresses unsteady bioconvection with concentration-dependent viscosity by formulating a coupled velocity–pressure–concentration system. It develops a decoupled Crank-Nicolson Leap-Frog (CNLF) finite element scheme using mini elements for $(\\mathbf{u},p)$ and a linear element for $c$, and proves unconditional stability via energy methods. Under Assumptions $A1$–$A2$ and a projection-based error framework, it establishes an $L^2$-error bound $\\max_{0\\le i \\le N-1} (\\|\\mathbf{u}^i-\\mathbf{u}_h^i\\|_{L^2}^2+\\|c^i-c_h^i\\|_{L^2}^2) \\le C(\\tau^4+h^4)$. Numerical experiments confirm second-order convergence in $L^2$ for $\\mathbf{u}$ and $c$, first-order for $p$, and robustness across viscosity models, highlighting improved efficiency over CN-based schemes.
Abstract
A fully discrete Crank--Nicolson Leap--Frog (CNLF) scheme is proposed and analyzed for the unsteady bioconvection flow problem with concentration-dependent viscosity. Spatial discretization is handled via the Galerkin finite element method (FEM), while temporal discretization employs the CNLF method for the linear terms and a semi-implicit approach for the nonlinear terms. The scheme is proven to be unconditionally stable, i.e., the time step is not subject to a restrictive upper bound. Using the energy method, $L^2$-optimal error estimates are derived for the velocity and concentration . Finally, numerical experiments are presented to validate the theoretical results.
