Analysis of a finite element method for the Stokes--Poisson--Boltzmann equations

1. Problem Statement

The paper addresses the mathematical modeling and numerical simulation of electro-osmotic flows in micro-channels and nanopores. These flows occur when an electrically charged fluid is subjected to an external electric field, causing motion due to the interaction with the electric double layer.

The governing system couples three physical phenomena:

1. Fluid Dynamics: Described by the stationary Stokes equations (incompressible flow).
2. Electrostatics: Described by a regularized Poisson–Boltzmann equation, which relates the electric potential to the charge density of the electrolyte.
3. Coupling Mechanism:
   • The electric potential influences the fluid momentum via an electric body force (drag).
   • The fluid velocity influences the potential distribution via an advection term (transport of ions).

The specific system of partial differential equations (PDEs) considered on a domain $\Omega$ involves:

• Momentum: $-\mu\Delta \mathbf{u} + \nabla p = \mathbf{f} - \varepsilon \Delta \psi \mathbf{E}$
• Mass Conservation: $\nabla \cdot \mathbf{u} = 0$
• Potential: $\kappa(\psi) + \mathbf{u} \cdot \nabla \psi - \varepsilon \Delta \psi = g$
  • Where $\kappa(\psi) = k_0 \sinh(k_1 \psi)$ represents the nonlinear charge density.
  • $\mathbf{E}$ is an externally applied electric field.
  • $\psi$ is the electrostatic double layer potential.

2. Methodology

A. Novel Formulation

The authors propose a specific reformulation of the coupling term in the momentum equation. Instead of treating the electric force term $-\varepsilon \Delta \psi \mathbf{E}$ directly (which would require integration by parts of the Laplacian of $\psi$), they use the Poisson–Boltzmann equation to rewrite this term as a weighted advection term.

• This results in a trilinear form $A_{\hat{\psi}}(\mathbf{u}, \mathbf{v}) = \int_\Omega [\mathbf{u} \cdot \nabla \hat{\psi}] (\mathbf{E} \cdot \mathbf{v})$.
• Advantage: This structure avoids the need for higher-order regularity on the potential in the momentum equation and facilitates the fixed-point analysis.

B. Functional Analysis and Well-Posedness

The paper establishes the existence and uniqueness of the weak solution using a fixed-point argument combined with classical functional analysis tools:

1. Spaces: Standard Sobolev spaces $H^1(\Omega)$ for velocity and potential, and $L^2(\Omega)$ for pressure. The potential is restricted to a bounded set $\Phi$ to handle the nonlinearity of the $\sinh$ function.
2. Decoupling Strategy:
   • Define an operator $S_{flow}$ that maps a fixed potential $\hat{\psi}$ to the velocity/pressure solution $(\mathbf{u}, p)$ of the Stokes problem.
   • Define an operator $S_{elec}$ that maps a fixed velocity $\hat{\mathbf{u}}$ to the potential solution $\psi$ of the nonlinear transport equation.
3. Theorems Used:
   • Babuška–Brezzi Theory: To ensure the well-posedness of the Stokes sub-problem (inf-sup condition).
   • Minty–Browder Theorem: To prove the existence and uniqueness of the solution to the nonlinear potential equation (using monotonicity of the $\sinh$ operator).
   • Banach's Contraction Principle: To prove that the composite operator $T = S_{flow} \circ S_{elec}$ has a unique fixed point, provided the data (forces, source terms) are sufficiently small (small-data assumption).

C. Finite Element Discretization

The authors propose a Galerkin Finite Element Method using generalized Taylor–Hood elements:

• Velocity ($\mathbf{u}_h$): Continuous piecewise polynomials of degree $k+1$.
• Pressure ($p_h$): Continuous piecewise polynomials of degree $k$.
• Potential ($\psi_h$): Continuous piecewise polynomials of degree $k+1$.
• Discrete Well-Posedness: The existence and uniqueness of the discrete solution are proven using techniques analogous to the continuous case.

D. Error Analysis

• Céa-type Estimates: The paper derives a priori error bounds showing that the error in the discrete solution is bounded by the best approximation error in the finite element spaces.
• Convergence Rates: Under standard regularity assumptions, the method achieves optimal convergence rates of order $O(h^s)$ in the $H^1$ norm for velocity and potential, and $O(h^s)$ in the $L^2$ norm for pressure, where $s$ depends on the polynomial degree.

3. Key Contributions

1. Novel Coupling Formulation: The rewriting of the electric body force as a weighted advection term simplifies the mathematical structure and avoids technical difficulties associated with the Laplacian of the potential in the momentum balance.
2. Rigorous Mathematical Analysis: The paper provides a complete proof of well-posedness (existence and uniqueness) for the coupled nonlinear system using a fixed-point approach, which is non-trivial due to the bidirectional coupling.
3. Discrete Stability and Convergence: The authors prove that the proposed finite element scheme is well-posed and derive optimal error estimates.
4. Numerical Validation: The theoretical results are validated through numerical experiments demonstrating optimal convergence rates and physical applicability.

4. Numerical Results

The authors implemented the scheme using the Gridap library and a Newton–Raphson solver. Three test cases were presented:

1. Convergence Test (Unit Square):
   • Used manufactured solutions to verify theoretical convergence rates.
   • Result: Optimal convergence rates were observed for polynomial degrees $k=1$ and $k=2$. For $k=1$, rates approached 2.0 in $H^1$; for $k=2$, rates approached 3.0.
2. Electro-osmotic Flow in a Micro-annulus:
   • Simulated flow in an eccentric micro-tube with different boundary potentials.
   • Result: The simulation captured higher fluid mobility near the narrow part of the channel and potential distributions consistent with previous literature, validating the model's ability to handle complex geometries.
3. Nanopore Sensor with Obstacles:
   • Simulated a 3D flow in a nanopore sensor with obstacles, incorporating both advection and convective nonlinearity.
   • Result: The model successfully simulated recirculation patterns and potential distributions induced by a slightly angled electric field, demonstrating the method's robustness in complex, realistic scenarios.

5. Significance

This work is significant for the computational physics and engineering communities studying microfluidics and lab-on-a-chip devices.

• Theoretical Rigor: It fills a gap in the literature by providing a rigorous mathematical foundation for the coupled Stokes–Poisson–Boltzmann system, which is often treated heuristically in engineering applications.
• Practical Applicability: The proposed numerical scheme is shown to be efficient and accurate for simulating electro-osmotic flows, which are critical for designing biomedical devices, water purification systems, and nanopore sensors.
• Methodological Innovation: The specific reformulation of the coupling term offers a new perspective that could be adapted for other coupled fluid-electrostatic problems, potentially simplifying implementation and analysis in future research.