An efficient Gauss-Newton-like method for the numerical solution of the Ornstein-Zernike integral equation for a class of fluid models

Journal Article · Tue Nov 01 04:00:00 EST 1994 · Journal of Computational Physics

DOI:https://doi.org/10.1006/jcph.1994.1174· OSTI ID:105468

Labik, S; Pospisil, R; Malijevsky, A ^[1]

---

1. Institute of Chemical Technology, Prague (Czechoslovakia); and others

A numerical algorithm for solving the Ornstein-Zernike (OZ) intergal equation of statistical mechanics is described for the class of fluids composed of molecules with axially symmetric interactions. Since the OZ equation is a nonlinear second-kind Fredholm equation whose key feature for the class of problems of interest is the highly computationally intensive nature of the kernel, the general approach employed in this paper is thus potentially useful for similar problems with this characteristic. The algorithm achieves a high degree of computational efficiency by combining iterative linearization of the most complex portion of the kernel with a combination of Newton-Raphson and Picard iteration methods for the resulting approximate equation. This approach makes the algorithm analogous to the approach of the classical Gauss-Newton method for nonlinear regression, and we call our method the GN algorithm. An example calculation is given illustration the use of the algorithm for the hard prolate ellipsoid fluid and its results are compared directly with those of the Picard iteration method. The GN algorithm is four to ten times as fast as the Picard method, and we present evidence that it is the most efficient general method currently available.

• OSTI ID: 105468
• Journal Information: Journal of Computational Physics, Journal Name: Journal of Computational Physics Journal Issue: 1 Vol. 115; ISSN 0021-9991; ISSN JCTPAH
• Country of Publication: United States
• Language: English

Related Subjects

42 ENGINEERING
66 PHYSICS
99 GENERAL AND MISCELLANEOUS
ALGORITHMS
FLOW MODELS
INTEGRAL EQUATIONS
NEWTON METHOD