Progress In Electromagnetics Research
ISSN: 1070-4698, E-ISSN: 1559-8985
Home | Search | Notification | Authors | Submission | PIERS Home | EM Academy
Home > Vol. 57 > pp. 237-252


By K. Chatterjee and J. Poggie

Full Article PDF (5,942 KB)

This paper presents a new three-dimensional floating random-walk (FRW) algorithm for the solution of the Nonlinear Poisson-Boltzmann (NPB) equation. The FRW method has not been previously used in the numerical solution of the NPB equation (and other nonlinear equations) because of the non-availability of analytical expressions for volumetric Green's functions. In the past, numerical studies using the FRW method have examined only the linearized Poisson-Boltzmann equation, producing solutions that are only accurate for small values of the potential. No such linearization is required for this algorithm. An approximate expression for a volumetric Green's functions has been calculated with the help of a novel use of perturbation theory, and the resultant integral form has been incorporated within the FRW framework. The algorithm requires no discretization of either the volume or the surface of the problem domains, and hence the memory requirements are expected to be lower than approaches based on spatial discretization, such as finite-difference methods. Another advantage of this algorithm is that each random walk is independent, so that the computational procedure is inherently parallelizable and an almost linear increase in computational speed is expected with increase in the number of processors. We have recently published the preliminary results for benchmark problems in one and two dimensions. In this work, we present our results for benchmark problems in three dimensions and demonstrate excellent agreement between the FRW- and finite-difference based algorithms. We also present the results of parallelization of the newly developed FRW algorithm. The solution of the NPB equation has applications in diverse branches of science and engineering including (but not limited to) the modeling of plasma discharges, semiconductor device modeling and the modeling of biomolecular structures and dynamics.

Citation: (See works that cites this article)
K. Chatterjee and J. Poggie, "A parallelized 3D floating random-walk algorithm for the solution of the nonlinear poisson-boltzmann equation," Progress In Electromagnetics Research, Vol. 57, 237-252, 2006.

1. Mitchner, M. and C. H. Kruger, Partial ly Ionized Gases, John Wiley & Sons, New York, 1973.

2. Sze, S. M., Physics of Semiconductor Devices, 366-368, 2nd ed., 366-368, John Wiley & Sons, New York, 1999.

3. Davis, M. E. and J. A. McCammon, "Electrostatics in bio-molecular structure and dynamics," Chemical Reviews, Vol. 90, 509-521, 1990.

4. Brown, G. M., Modern Mathematics for Engineers, E. F. Beckenbach (ed.), McGraw-Hill, New York, 1956.

5. Le Coz, Y. L., H. J. Greub, and R. B. Iverson, "Performance of random walk capacitance extractors for IC interconnects: a numerical study," Solid-State Electronics, Vol. 42, 581-588, 1998.

6. Le Coz, Y. L., R. B. Iverson, T. L. Sham, H. F. Tiersten, and M. S. Shepard, "Theory of a floating random walk algorithm for solving the steady-state heat equation in complex materially inhomogeneous rectilinear domains," Numerical Heat Transfer, Vol. 26, 353-366, 1994.

7. Haberman, R., Elementary Applied Partial Differential Equations, 3rd ed., Prentice-Hall, Englewood Cliffs, NJ, 1998.

8. Hwang, C.-O. and M. Mascagni, "Efficient modified 'walk on spheres' algorithm for the linearized Poisson-Boltzmann equation," App. Phys. Lett., Vol. 78, No. 6, 787-789, 2001.

9. Mascagni, M. and N. A. Simonov, "Monte carlo methods for calculating the electrostatic energy of a molecule," Proceedings of the 2003 International Conference on Computational Science (ICCS 2003)..

10. Chatterjee, K. and J. Poggie, "A meshless stochastic algorithm for the solution of the nonlinear Poisson-Boltzmann equation in the context of plasma discharge modeling: 1D analytical benchmark," Proceedings of the 17th AIAA Computational Fluid Dynamics Conference, 6-9.

11. Chatterjee, K. and J. Poggie, "A two-dimensional floating random-walk algorithm for the solution of the nonlinear Poisson-Boltzmann equation: application to the modeling of plasma sheaths," Proceedings of the 3rd MIT Conference on Computational Fluid and Solid Mechanics, 14-17, 2005.

12. Sobol, I. M., A Primer for the Monte Carlo Method, CRC Press, Boca Raton, 1994.

13. Curtiss, J. H., "Monte carlo methods for the iteration of linear operators," J. Math. and Phys., Vol. 32, 209-232, 1954.

14. Sabelfeld, K. K. and N. A. Simonov, Random Walks on Boundary for Solving PDEs, 119-120, 119-120, VSP, Utrecht, The Netherlands, 1994.

15. Chatterjee, K., "Development of a floating random walk algorithm for solving Maxwell's equations in complex IC-interconnect structures," Rensselaer Polytechnic Institute, 48106-1346.

16. Chatterjee, K. and Y. L. Le Coz, "A floating random-walk algorithm based on iterative perturbation theory: solution of the 2D vector-potential maxwell-helmholtz equation," Applied Computational Electromagnetics Society Journal, Vol. 18, No. 1, 48-57, 2003.

17. Hammersley, J. M. and D. C. Handscomb, Monte Carlo Methods, 50-54, 50-54, John Wiley & Sons, New York, NY, 1964.

© Copyright 2014 EMW Publishing. All Rights Reserved