1
|
Yang W, Peng L, Zhu Y, Hong L. When machine learning meets multiscale modeling in chemical reactions. J Chem Phys 2020; 153:094117. [DOI: 10.1063/5.0015779] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Wuyue Yang
- Yau Mathematical Sciences Center, Tsinghua University, Beijing 100084, People’s Republic of China
| | - Liangrong Peng
- College of Mathematics and Data Science, Minjiang University, Fuzhou 350108, People’s Republic of China
| | - Yi Zhu
- Yau Mathematical Sciences Center, Tsinghua University, Beijing 100084, People’s Republic of China
| | - Liu Hong
- School of Mathematics, Sun Yat-sen University, Guangzhou 510275, People’s Republic of China
| |
Collapse
|
2
|
Leal AMM, Kulik DA, Smith WR, Saar MO. An overview of computational methods for chemical equilibrium and kinetic calculations for geochemical and reactive transport modeling. PURE APPL CHEM 2017. [DOI: 10.1515/pac-2016-1107] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
AbstractWe present an overview of novel numerical methods for chemical equilibrium and kinetic calculations for complex non-ideal multiphase systems. The methods we present for equilibrium calculations are based either on Gibbs energy minimization (GEM) calculations or on solving the system of extended law of mass-action (xLMA) equations. In both methods, no a posteriori phase stability tests, and thus no tentative addition or removal of phases during or at the end of the calculations, are necessary. All potentially stable phases are considered from the beginning of the calculation, and stability indices are immediately available at the end of the computation to determine which phases are actually stable at equilibrium. Both GEM and xLMA equilibrium methods are tailored for computationally demanding applications that require many rapid local equilibrium calculations, such as reactive transport modeling. The numerical method for chemical kinetic calculations we present supports both closed and open systems, and it considers a partial equilibrium simplification for fast reactions. The method employs an implicit integration scheme that improves stability and speed when solving the often stiff differential equations in kinetic calculations. As such, it requires compositional derivatives of the reaction rates to assemble the Jacobian matrix of the resultant implicit algebraic equations that are solved at every time step. We present a detailed procedure to calculate these derivatives, and we show how the partial equilibrium assumption affects their computation. These numerical methods have been implemented in Reaktoro (reaktoro.org), an open-source software for modeling chemically reactive systems. We finish with a discussion on the comparison of these methods with others in the literature.
Collapse
Affiliation(s)
- Allan M. M. Leal
- 1Geothermal Energy and Geofluids Group, Department of Earth Sciences, ETH Zürich, Switzerland
| | - Dmitrii A. Kulik
- 2Laboratory for Waste Management, Nuclear Energy and Safety Research Department, Paul Scherrer Institut, 5232 Villigen PSI, Switzerland
| | - William R. Smith
- 3Department of Mathematics and Statistics, University of Guelph, Guelph, Ontario NIG2 WI, Canada
| | - Martin O. Saar
- 1Geothermal Energy and Geofluids Group, Department of Earth Sciences, ETH Zürich, Switzerland
- 4Department of Earth Sciences, University of Minnesota, Minneapolis, USA
| |
Collapse
|