A decoupled energy-stable mixed finite element method for Poisson-Nernst-Planck-Navier-Stokes equations
This work provides a computationally efficient and structure-preserving numerical method for simulating complex electrokinetic flows, which is important for applications in microfluidics and biophysics.
The authors propose a linearized mixed finite element method for the Poisson-Nernst-Planck-Navier-Stokes system that eliminates nonlinear solvers via staggered time discretization and preserves physical properties like mass conservation and energy stability. Numerical experiments with mimetic spectral elements verify accuracy and effectiveness.
We propose a novel linearized mixed finite element method for the Poisson-Nernst-Planck-Navier-Stokes (PNPNS) system. Specifically, the method combines a staggered time discretization that eliminates the need for expensive nonlinear solvers by carefully treating nonlinear terms in a time-staggered manner, with a mimetic spatial discretization that preserves the exact structure of the discrete de Rham complex. Both semi-discrete scheme and its fully discrete counterpart are developed, which preserve key physical properties, including conservation of the total mass and energy stability. Under appropriate assumptions on the initial data, a rigorous theoretical analysis is carried out for the fully discrete scheme. Numerical experiments using mimetic spectral elements are presented to demonstrate the properties and to verify the accuracy and effectiveness of the proposed decoupled approach.