[Paper Review] Modeling electrochemical systems with weakly imposed Dirichlet boundary conditions

Sun-Gu Kim, Khanwale, Makrand A.+2 more|arXiv (Cornell University)|Oct 17, 2020

Nanopore and Nanochannel Transport Studies43 references4 citations

TL;DR

This paper presents a finite element framework that uses weak imposition of Dirichlet boundary conditions via a Dirichlet-to-Neumann transformation to accurately compute current fluxes in electrochemical systems governed by the coupled Navier-Stokes-Poisson-Nernst-Planck equations. The method enables accurate resolution of thin Debye layers with coarse meshes—demonstrating superior convergence and reduced computational cost compared to strong imposition, even in complex 3D and unstable electrokinetic flows.

ABSTRACT

Finite element modeling of charged species transport has enabled analysis, design, and optimization of a diverse array of electrochemical and electrokinetic devices. These systems are represented by the Poisson-Nernst-Planck equations coupled with the Navier-Stokes equation, with a key quantity of interest being the current at the system boundaries. Accurately computing the current flux is challenging due to the small critical dimension of the boundary layers (small Debye layer) that require fine mesh resolution at the boundaries. We resolve this challenge by using the Dirichlet-to-Neumanntransformation to weakly impose the Dirichlet conditions for the Poisson-Nernst-Planck equations. The results obtained with weakly imposed Dirichlet boundary conditions showed excellent agreement with those obtained when conventional boundary conditions with highly resolved mesh we reemployed. Furthermore, the calculated current flux showed faster mesh convergence using weakly imposed conditions compared to the conventionally imposed Dirichlet boundary conditions. We illustrate the approach on canonical 3D problems that otherwise would have been computationally intractable to solve accurately. This approach substantially reduces the computational cost of model-ing electrochemical systems.

Motivation & Objective

• To address the computational challenge of resolving thin Debye layers in electrochemical systems with high aspect-ratio boundary layers.
• To reduce the computational cost of simulating accurate current fluxes at boundaries in multi-scale electrokinetic systems.
• To develop a robust numerical framework that maintains accuracy with coarse meshes near boundaries, avoiding the need for extreme mesh refinement.

Proposed method

• Uses a weak imposition of Dirichlet boundary conditions via a Dirichlet-to-Neumann transformation for the Poisson-Nernst-Planck equations.
• Employs a block-iterative solution strategy to couple Navier-Stokes and PNP equations efficiently.
• Applies geometric mesh grading with a progression ratio of 1.021 to cluster elements near boundaries without requiring extremely fine resolution.
• Utilizes manufactured solutions and analytical benchmarks (e.g., electroosmotic flow) for validation.
• Applies the method to 1D, 2D, and 3D problems, including electrokinetic instabilities near perm-selective membranes.
• Solves the coupled system with a time step of Δt = 1×10⁻⁶ in non-dimensionalized form.

Experimental results

Research questions

• RQ1 Can weak imposition of Dirichlet conditions produce accurate current fluxes in electrochemical systems with thin Debye layers, even when mesh resolution is coarse?
• RQ2 How does the convergence rate of the weak imposition method compare to strong imposition under mesh refinement near boundaries?
• RQ3 Can the framework accurately simulate complex electrokinetic phenomena such as ion concentration polarization and electroconvective instabilities with minimal mesh refinement?
• RQ4 Does the method maintain accuracy across varying flow conditions and boundary layer thicknesses, including in 3D geometries?

Key findings

• The weak imposition method achieves excellent agreement with strong imposition results, even with only 1–2 elements across the Debye layer (Λ = 1×10⁻³), while strong imposition requires significantly finer meshes.
• The method converges at a faster rate with mesh refinement near boundaries compared to strong imposition, indicating superior numerical efficiency.
• Boundary fluxes computed via weak imposition show negligible error (residual error ~10⁻⁸) compared to analytical and manufactured solutions.
• The framework successfully captures electrokinetic instabilities near a cation-selective membrane with a coarse mesh (1280×180 quad elements), reproducing finger-like charge density structures seen in prior benchmarks.
• Simulations of 3D branching microchannels for desalination demonstrate the method’s scalability and robustness in complex geometries.
• The approach reduces computational cost substantially by avoiding the need for extremely fine boundary-layer meshes, making previously intractable 3D electrochemical simulations feasible.