High order positivity-preserving numerical methods for a non-local photochemical model

In this paper we design high-order positivity-preserving approximation schemes for an integro-differential model describing photochemical reactions. Specifically, we introduce and analyze three classes of dynamically consistent methods, encompassing non-standard finite difference schemes, direct quadrature techniques and predictor- corrector approaches. The proposed discretizations guarantee the positivity, monotonicity and boundedness of the solution regardless of the temporal, spatial and frequency stepsizes. Comprehensive numerical experiments confirm the theoretical findings and demonstrate the efficacy of the proposed methods in simulating realistic photochemical phenomena.