All publications

Journal article ·

Parametrization and Cartesian representation techniques for robust resolution of chemical equilibria

Maxime Jonval, Ibtihel Ben Gharbia, Clément Cancès, Thibault Faney, Quang-Huy Tran

Journal of Computational Physics

Abstract

Chemical equilibria computations, especially those with vanishing species in the aqueous phase, lead to nonlinear systems that are difficult to solve due to gradient blow up. Instead of the commonly used ad hoc treatments, we propose two reformulations of the single-phase chemical equilibrium problem which are in line with the spirit of preconditioning but whose actual aims are to guarantee a better stability of Newton's method. The first reformulation is a parametrization of the graph linking species mole fractions to their chemical potentials. The second is based on an augmented system where this relationship is relaxed for the iterates by means of a Cartesian representation. We theoretically prove the local quadratic convergence of Newton's method for both reformulations. From a numerical point of view, we demonstrate that the two techniques are accurate, allowing to compute equilibria with chemical species having very low concentrations. Moreover, the robustness of our methods combined with a globalization strategy is superior to that of the literature.

Keywords

  • Cartesian representation
  • Chemical equilibria
  • Newton's method
  • Parametrization
BibTeX citation
@article{jonval_parametrization_2025,
 abstract = {Chemical equilibria computations, especially those with vanishing species in the aqueous phase, lead to nonlinear systems that are difficult to solve due to gradient blow up. Instead of the commonly used ad hoc treatments, we propose two reformulations of the single-phase chemical equilibrium problem which are in line with the spirit of preconditioning but whose actual aims are to guarantee a better stability of Newton's method. The first reformulation is a parametrization of the graph linking species mole fractions to their chemical potentials. The second is based on an augmented system where this relationship is relaxed for the iterates by means of a Cartesian representation. We theoretically prove the local quadratic convergence of Newton's method for both reformulations. From a numerical point of view, we demonstrate that the two techniques are accurate, allowing to compute equilibria with chemical species having very low concentrations. Moreover, the robustness of our methods combined with a globalization strategy is superior to that of the literature.},
 author = {Jonval, Maxime and Ben Gharbia, Ibtihel and Cancès, Clément and Faney, Thibault and Tran, Quang-Huy},
 doi = {10.1016/j.jcp.2024.113596},
 file = {ScienceDirect Snapshot:C\:\\Users\\Tibo\\Zotero\\storage\\47E7V3EM\\S0021999124008441.html:text/html;Submitted Version:C\:\\Users\\Tibo\\Zotero\\storage\\9WPWG8QZ\\Jonval et al. - 2025 - Parametrization and Cartesian representation techniques for robust resolution of chemical equilibria.pdf:application/pdf},
 issn = {0021-9991},
 journal = {Journal of Computational Physics},
 keywords = {Cartesian representation, Chemical equilibria, Newton's method, Parametrization},
 month = {February},
 pages = {113596},
 title = {Parametrization and Cartesian representation techniques for robust resolution of chemical equilibria},
 url = {https://www.sciencedirect.com/science/article/pii/S0021999124008441},
 urldate = {2025-04-20},
 volume = {522},
 year = {2025}
}
Download citation