Multidomain Chebyshev pseudo-spectral method applied to the Poisson–Boltzmann equation for two parallel plates

We address a boundary-value problem involving a Poisson–Boltzmann equation that models the electrostatic potential of a channel formed by parallel plates with an electrolyte solution confined between the plates. We show the existence and uniqueness of solution to the problem, with special (particular) solutions as bounds, namely, a Debye–Hückel type solution as lower bound and a Gouy–Chapman type solution as upper bound. Our results are based on the maximum principle for elliptic equations and are useful for characterizing the behavior of the solutions. Also, we introduce a numerical scheme based on the Chebyshev pseudo-spectral method to calculate approximate solutions. This method is applied in conjunction with a multidomain procedure that attempts to capture the dramatic exponential increase/decay of the solution near the plates.