Lattice evolution solution for the nonlinear Poisson–Boltzmann equation in confined domains

The lattice evolution method for solving the nonlinear Poisson–Boltzmann equation in confined domain is developed by introducing the second-order accurate Dirichlet and Neumann boundary implements, which are consistent with the non-slip model in lattice Boltzmann method for fluid flows. The lattice evolution method is validated by comparing with various analytical solutions and shows superior to the classical numerical solvers of the nonlinear Poisson equations with Neumann boundary conditions. The accuracy and stability of the method are discussed. This lattice evolution nonlinear Poisson–Boltzmann solver is suitable for efficient parallel computing, complex geometry conditions, and easy extension to three-dimensional cases.