K. M. Furati

K. D. Olumoyin, A. Q. M. Khaliq, K. M. Furati

Mutating variants of COVID-19 have been reported across many US states since 2021. In the fight against COVID-19, it has become imperative to study the heterogeneity in the time-varying transmission rates for each variant in the presence of pharmaceutical and non-pharmaceutical mitigation measures. We develop a Susceptible-Exposed-Infected-Recovered mathematical model to highlight the differences in the transmission of the B.1.617.2 delta variant and the original SARS-CoV-2. Theoretical results for the well-posedness of the model are discussed. A Deep neural network is utilized and a deep learning algorithm is developed to learn the time-varying heterogeneous transmission rates for each variant. The accuracy of the algorithm for the model is shown using error metrics in the data-driven simulation for COVID-19 variants in the US states of Florida, Alabama, Tennessee, and Missouri. Short-term forecasting of daily cases is demonstrated using long short term memory neural network and an adaptive neuro-fuzzy inference system.

K. Mustapha, B. Abdallah, K. M. Furati et al.

We propose a piecewise-linear, time-stepping discontinuous Galerkin method to solve numerically a time fractional diffusion equation involving Caputo derivative of order $μ\in (0,1)$ with variable coefficients. For the spatial discretization, we apply the standard piecewise linear continuous Galerkin method. Well-posedness of the fully discrete scheme and error analysis will be shown. For a time interval~$(0,T)$ and a spatial domain~$Ω$, our analysis suggest that the error in $L^2\bigr((0,T),L^2(Ω)\bigr)$-norm is of order $O(k^{2-\fracμ{2}}+h^2)$ (that is, short by order $\fracμ{2}$ from being optimal in time) where $k$ denotes the maximum time step, and $h$ is the maximum diameter of the elements of the (quasi-uniform) spatial mesh. However, our numerical experiments indicate optimal $O(k^{2}+h^2)$ error bound in the stronger $L^\infty\bigr((0,T),L^2(Ω)\bigr)$-norm. Variable time steps are used to compensate the singularity of the continuous solution near $t=0$.