We present a systematic derivation of asymptotic expansions of nonequilibrium thermodynamics for chemical reaction networks (CRNs) based on singular perturbation theory. For a general reversible CRN with fast–slow kinetics, we obtain the first and second laws of thermodynamics for the asymptotic expansion model. We derive composite expansions of the enthalpy, entropy, entropy production rate, and relative entropy. The slow-varying outer parts of these thermodynamic quantities capture the long-time trend, while the fast-varying corrected inner parts decay to zero as the fast time variable tends to infinity. The convergence order of these quantities is determined by the local Lipschitz properties of the respective functions. The enthalpy retains the same convergence order as the kinetic variables, whereas the entropy, relative entropy, and entropy production rate involve logarithmic terms that cause their gradients to diverge when some concentrations approach zero, reducing their theoretical convergence order by one. The general theory is validated on the reversible Michaelis–Menten reaction, for which both leading-order and first-order matched asymptotic expansions are obtained analytically. Numerical simulations confirm the uniform accuracy of the composite thermodynamic approximations and further reveal that the entropy production rate converges with a higher order than theoretically predicted. The results demonstrate that the composite expansion provides a rigorous and physically consistent tool for analyzing energy and entropy balances in multiscale CRNs.
Peng et al. (Mon,) studied this question.