Parametrization and Cartesian representation techniques for robust resolution of chemical equilibria - INRIA 2
Article Dans Une Revue Journal of Computational Physics Année : 2025

Parametrization and Cartesian representation techniques for robust resolution of chemical equilibria

Résumé

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.
Fichier principal
Vignette du fichier
paper (1).pdf (928.28 Ko) Télécharger le fichier
Origine Fichiers produits par l'(les) auteur(s)
licence

Dates et versions

hal-04225504 , version 1 (02-10-2023)
hal-04225504 , version 2 (30-04-2024)

Licence

Identifiants

Citer

Maxime Jonval, Ibtihel Ben Gharbia, Clément Cancès, Thibault Faney, Quang Huy Tran. Parametrization and Cartesian representation techniques for robust resolution of chemical equilibria. Journal of Computational Physics, 2025, 522, pp.113596. ⟨10.1016/j.jcp.2024.113596⟩. ⟨hal-04225504v2⟩
210 Consultations
185 Téléchargements

Altmetric

Partager

More