Journal Home Online First Current Issue Archive For Authors Journal Information 中文版

Frontiers of Chemical Science and Engineering >> 2024, Volume 18, Issue 3 doi: 10.1007/s11705-024-2390-8

A projected Newton algorithm based on chemically allowed interval for chemical equilibrium computations

Institute of Process Systems Engineering, College of Chemical Engineering, Qingdao University of Science and Technology, Qingdao 266042, China;Institute of Process Systems Engineering, College of Chemical Engineering, Qingdao University of Science and Technology, Qingdao 266042, China;Institute of Process Systems Engineering, College of Chemical Engineering, Qingdao University of Science and Technology, Qingdao 266042, China;Institute of Process Systems Engineering, College of Chemical Engineering, Qingdao University of Science and Technology, Qingdao 266042, China;Institute of Process Systems Engineering, College of Chemical Engineering, Qingdao University of Science and Technology, Qingdao 266042, China

Received: 2023-09-08 Available online: 2023-09-08

Next Previous

Abstract

The chemical equilibrium equations utilized in reactive transport modeling are complex and nonlinear, and are typically solved using the Newton-Raphson method. Although this algorithm is known for its quadratic convergence near the solution, it is less effective far from the solution, especially for ill-conditioned problems. In such cases, the algorithm may fail to converge or require excessive iterations. To address these limitations, a projected Newton method is introduced to incorporate the concept of projection. This method constrains the Newton step by utilizing a chemically allowed interval that generates feasible descending iterations. Moreover, we utilize the positive continuous fraction method as a preconditioning technique, providing reliable initial values for solving the algorithms. The numerical results are compared with those derived using the regular Newton-Raphson method, the Newton-Raphson method based on chemically allowed interval updating rules, and the bounded variable least squares method in six different test cases. The numerical results highlight the robustness and efficacy of the proposed algorithm.

Related Research