Multiscale solution decomposition of nonlocal-in-time problems with application in numerical computation

Abstract: This work develops a multiscale solution decomposition (MSD) method for nonlocal-in-time problems to separate a series of known terms with multiscale singularity from the original singular solution such that the remaining unknown part becomes smoother. We demonstrate that the MSD provides a scenario where the smoothness assumption for solutions of weakly singular nonlocal-in-time problems, a commonly encountered assumption in numerous literature of numerical methods that is in general not true for original solutions, becomes appropriate such that abundant numerical analysis results therein become applicable. From computational aspect, instead of handling solution singularity, the MSD significantly reduces the numerical difficulties by separating and thus circumventing the solution singularity. We consider typical problems, including the fractional relaxation equation, Volterra integral equation, subdiffusion, integrodifferential equation and diffusion-wave equation, to demonstrate the universality of MSD and its effectiveness in improving the numerical accuracy or stability in comparison with classical methods.