Investigation of the Boundary Value Problem for an Extended System of Stationary Nernst–Planck–Poisson Equations in the Diffusion Layer

This article investigates the boundary value problem for an extended stationary system of Nernst–Planck–Poisson equations, corresponding to a mathematical model of the influence of changes in the equilibrium coefficient on the transport of ions of a binary salt in the diffusion layer. Dimensionless variables were introduced using characteristic parameter values. As a result, a dimensionless boundary value problem was obtained, which is singularly perturbed, containing a small parameter in the derivative of the Poisson equation and, additionally, another regular small parameter. A similarity theory was developed: trivial and non-trivial similarity criteria and their physical meaning were determined, which allowed for the identification of general properties of the solutions. A numerical investigation of the boundary value problem was conducted using the finite element method. With an increase in the initial solution concentration, the value of the small parameter entering singularly decreases, reaching values on the order of 10−12 and below, leading to computational difficulties that prevent a comprehensive analysis of the influence of changes in the equilibrium coefficient on salt ion transport. In this regard, an analytical solution to the problem was constructed, based on dividing the solution domain into several subdomains (regions of electroneutrality, extended space charge region, quasi-equilibrium region, recombination region, intermediate layer), in each of which the problem is solved differently, followed by matching these solutions. Verification of the analytical solution was carried out by comparing it with the numerical solution. The advantage of the obtained analytical solution is the possibility of a comprehensive analysis of the influence of the dissociation/recombination reaction of water molecules on salt ion transport over a wide range of real changes in the concentration and composition of the electrolyte solution and other input parameters. This boundary value problem serves as a benchmark for constructing asymptotic solutions for other singularly perturbed boundary value problems in membrane electrochemistry.