key: cord-0697079-4pttfjxd authors: Mohammad, Mutaz; Trounev, Alexander title: On the Dynamical Modeling of Covid-19 involving Atangana-Baleanu fractional derivative and based on Daubechies framelet simulations date: 2020-07-28 journal: Chaos Solitons Fractals DOI: 10.1016/j.chaos.2020.110171 sha: 2c1b8e54333c1b3a88bbfbfd3ab9144c435f3c6e doc_id: 697079 cord_uid: 4pttfjxd In this paper, we present a novel fractional order COVID-19 mathematical model by involving fractional order with specific parameters. The new fractional model is based on the well-known Atangana-Baleanu fractional derivative with non-singular kernel. The proposed system is developed using eight fractional-order nonlinear differential equations. The Daubechies framelet system of the model is used to simulate the nonlinear differential equations presented in this paper. The framelet system is generated based on the quasi-affine setting. In order to validate the numerical scheme, we provide numerical simulations of all variables given in the model. The novel corona-virus is a new strain of coronavirus which may cause illness, fever, dry cough where these symptoms are usually mild and begin gradually. The world health organization has declared this virus as a pandemic in early March of 2020 where many countries have taken serious actions and implemented curfew, quarantine and lock-down measures as a plan to control the rapid spread of COVID-19. The first case of COVID-19 was detected in Wuhan city in China at the end of the year of 2019 where it is suggested that the COVID-19 virus might be originated from bats and it's transmission might related to a seafood market exposure. Many researchers worldwide started to work on developing mathematical models that best describe the dynamics of this pandemic. It is known in biological system with memory it would be suitable to use fractional derivatives to describe evolution of the system [1] [2] [3] [4] [5] [6] [7] [8] . Furthermore, Atangana-Baleanu fractional derivative (ABFD) has been one of the most useful operators for modeling non-local behaviors by fractional differential equations. The advantage of using such derivative lies on its properties such as the non-locality and non-singularity of its kernel, and the crossover behavior in the model can only be best described using this derivative. Additionally, it allows traditional and various types of initial conditions to be consider in the creation of the dynamical model. Many scientists proposed new models to best describe the dynamics of all possible parameters responsible for the daily cases reported including deaths, control the fatality rate, and prediction of COVID-19 behavior in future within a specific region. It is known that several models can describe the same system, which is a challenging step. In this paper we intend to formulate a new mathematical model of Corona virus based on the model presented in [9] based on ABFD. The numerical method simulation is conducted via the framelet system generated using Daubechies scaling functions. Daubechies wavelets have been proven as a useful tool in a variety of various applications such as filter banks constructions in image painting. This is largely due to the fact that wavelets have the right structure to capture the sparsity in physical images, perfect mathematical properties such as its multiscale structure, sparsity, smoothness, compactly supported, and high vanishing moments properties. It has many applications in fractional integral and differential equations (see for example [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] . Framelets have been used extensively in the context of both pure and numerical methods in several applications, due to their well prevailing and recognized theory and its natural properties such as sparsity and stability which lead to a well-conditioned scheme. In this paper, an effective and accurate technique based on Daubechies wavelets is presented for solving the transmission model of COVID-19 based on Caputo fractional derivative. The advantage of using such wavelets, lies on its simple structure of the reduced systems and in the powerfulness of obtaining approximated solutions for such equations that have weakly singular kernels. The proposed method shows a good performance and high accuracy orders. where a[k] ∈ 2 (Z) is finitely supported sequence and is called the refinement mask of φ. The corresponding wavelet function is defined by is finitely supported sequence and is called the high pass filter of ψ. For a function f ∈ L 1 (R) (which can be naturaly extended to L 2 (R)), we use the following Fourier transform defined byf The Fourier series of the sequence a is defined bŷ 2 Daubechies framelets using the unitary extension principle (UEP) If g is a wavelet function that has q vanishing moments such that t m g (t) dt = 0, m = 0, 1, . . . , q − 1. Suppose that the function g generates an orthonormal basis of L 2 (R), then the constructed wavelet will be compactly supported within the domain [0, 2q − 1]. Daubechies wavelets do not have explicit form but defined recursively as follows One of the important features of this wavelet is its smoothness as it increases for a higher q. We present the graphs of φ and its corresponding wavelet ψ when q = 1, 2 and 3, 4 in Figures 1 and 2 respectively. The constants r, R are called frame bounds [28] . A frame is called tight if we have r = R as frame bounds, and it is Parseval frame if r = R = 1. The idea is to construct framelet system based on Daubechies scaling function φ and its corresponding wavelet function ψ. is a finitely supported sequence. Define the wavelet system Theorem 2.2 (UEP [28] ). Assume that φ ∈ L 2 (R) be a compactly supported scaling function. Let be a set of finitely supported sequences, then According to Theorem 2.2, for any constructed framelet system we have the following representation given by This system can be truncated by U n f as follows Here we provide some examples on the construction of framelet systems basd on several orders of Daubechies scaling functions of different orders. Then based on the UEP, we can find two finitely supported sequences b 1 [k], b 2 [k] such that the following two functions generate a framelet system of Note that, according to the UEP we need to solve the following system of equation written in MATLAB software to be able to get the required sequences end; and for q = 2 we have the following where x(k) is the nonzero value of the compactly supported sequences of both b 1 and b 2 . Note that, when q = 2, we have the following low pass filter The graphs of Daubechies scaling functions of order one and two along with their corresponding framelets are depicted in Figures 3. Then based on the UEP, we can find two finitely supported sequences b 1 [k], b 2 [k] such that two functions ψ 1 , ψ 2 can generate a framelet system of L 2 (R). Note that, according to the UEP we need to solve the following system of equation written in MATLAB software to be able to get the required (1)) 2 + (x(2)) 2 + (x(3)) 2 + (x(4)) 2 + (x(5)) 2 + (x(6)) 2 + (x(7)) 2 + (x(8)) 2 + (x(9)) 2 + (x(10)) 2 + (x(11)) 2 + (x(12)) 2 − .5; x(8) * x(9) + x(9) * x(10) + x(10) * x(11) + x(11) * x(12) + (75/256); (1)) 2 + (x(2)) 2 + (x(3)) 2 + (x(4)) 2 + (x(5)) 2 + (x(6)) 2 + (x(7)) 2 + (x(8)) 2 + (x(9)) 2 + (x(10)) 2 + (x(11)) 2 + (x(12)) 2 + (x(13)) 2 + (x(14)) 2 + (x(15)) 2 + (x(16) (x(1)) 2 − (x(2)) 2 + (x(3)) 2 − (x(4)) 2 + (x(5)) 2 − (x(6)) 2 + (x(7)) 2 − (x(8)) 2 + (x(9)) 2 − (x(10)) 2 + (x(11)) 2 − (x(12)) 2 + (x(13)) 2 − (x(14)) 2 + (x(15)) 2 − (x(16)) 2 + (−26/2023); end; We present the graphs of Daubechies scaling functions of order three and four along with their corresponding framelets in Figures 4. Given the construction in the first part above and to simulate the resulting equations, now we are ready to introduce the new COVID-19 fractional model of nonlinear differential equations by applying the Atangana-Baleanu derivative. The advantage of using such framelet lies in its properties such as the highest number of vanishing moments, redundancy and its applications in solving a broad range of problems such as fractal problems and function discontinuities, see e.g., [12] Herein, we consider the model presented in [9] using ABFD. The model has eight nonlinear DEs. To simulate the system and for simplicity, we consider Daubechies framelet system of order one. Hence, the new modified model that obtained by changing the left hand side of the system presented in [9] by involving ABFD. Before presenting the new model in fractional sense, let us provide the definition of ABFD and its associated integral. The advantage of using such framelets lies Definition 3.1. For a real function u(t) where t, α > 0 and n ∈ N, we have the following fractional operators of order α, namely: where B(α) is a normalization function such that B(0) = B(1) = 0 and M α is the MittagLeffler function. • The integral operator corresponding to this definition is given by We refer the reader to [4, 5] for more details and properties of the fractional derivative. Therefore, the new model can be written as follows ABC a D α t R(t) = H(t)γ r + γ i (P (t) + I(t)); with the initial conditions where the model parameters and its values are given in Table 1 for which the reproduction number We provide a numerical scheme based on the collocation technique by discretizing the domain function across the Daubechies framelet system being used to solve the proposed COVID-19 model. Therefore, by truncating each unknown variable using the truncated partial sum given in Equation The susceptible cases -E(t) The exposed cases - Symptomatic and infectious class - Super-spreaders class - Infectious but asymptomatic class - Recovery class - Fatality class β Transmission coefficient from infected individuals 2.55 ι Relative transmissibility of hospitalized patients 1.56 β Transmission coefficient due to super-spreaders 7.65 κ Rate at which exposed become infectious 0.25 ρ 1 Rate at which exposed people become infected I 0.580 ρ 2 Rate at which exposed people become super-spreaders 0.001 γ a Rate of being hospitalized 0.94 γ i Recovery rate without being hospitalized 0.27 γ r Recovery rate of hospitalized patients 0.50 δ i Disease induced death rate due to infected class 3.5 δ p Disease induced death rate due to super-spreaders 1.00 δ h Disease induced death rate due to hospitalized class 0.30 (2.6) generated using Daubechies framelet, our new model will take the following structure such that the derivative of each variable takes the following approximation and the coefficient C S , C E , . . . , C F to be determined. Applying the algorithm proposed in [29] yields the following In the present paper, we presented a COVID-19 model with new fractional operator using ABFD. This mathematical and dynamical model is more suitable to describe the biological phenomena with memory than the integer order model. To test the behavior of all variables of the model, we simulated the resulting nonlinear fractional differential equations model by involving ABFD based on Daubechies framelet systems and obtained various graphical illustrations. It turns out that, increasing of the fractional value of the parameters resulting a decrease in the infection rates. A new application of fractional AtanganaBaleanu derivatives: Designing ABC-fractional masks in image processing A fractional order optimal 4D chaotic financial model with Mittag-Leffler law Fractional differential and integral operators with non-singular and non-local kernel with application to nonlinear dynamical systems New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model Chaos in a simple nonlinear system with AtanganaBaleanu derivatives with fractional order On the new fractional derivative and application to nonlinear Fishers reactiondiffusion equation Decolonisation of fractional calculus rules: breaking commutativity and associativity to capture more natural phenomena, the Modelling the spread of COVID-19 with new fractal-fractional operators: Can the lockdown save mankind before vaccination? Mathematical modeling of Covid-19, 205, transmission dynamics with a case study of Wuhan Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative Gibbs phenomenon in tight framelet expansions Gibbs Effects Using Daubechies and Coiflet Tight Framelet Systems, Contemporary Mathematics Special B-spline Tight Framelet and It's Applications On the Gibbs Effect Based on the Quasi-Affine Dual Tight Framelets System Generated Using the Mixed Oblique Extension Principle Wavelets Based Simulation and Visualization Approach for Unmixing of Hyperspectral Data Biorthogonal-Wavelet-Based Method for Numerical Solution of Volterra Integral Equations M Mohammad A Numerical Solution of Fredholm Integral Equations of the Second Kind Based on Tight Framelets Generated by the Oblique Extension Principle A collocation method via the quasi-affine biorthogonal systems for solving weakly singular type of Volterra-Fredholm integral equations Applications of bi-framelet systems for solving fractional order differential equations The dynamics of COVID-19 in the UAE based on fractional derivative modeling using Riesz wavelets simulation New results on nondensely characterized integrodifferential equations with fractional order Haci Mehmet Baskonus, Cancer treatment model with the Caputo-Fabrizio fractional derivative Haci Mehmet Baskonus, Solutions of partial differential equations using the fractional operator involving Mittag-Leffler kernel Haci Mehmet Baskonus Analytic study for a fractional model of HIV infection of CD4+TCD4+T lymphocyte cells Mathematics in Natural Science Haci Mehmet Baskonus, Existence results of Hilfer integro-differential equations with fractional order Framelets and wavelets: Algorithms, analysis, and applications. Applied and Numerical Harmonic Analysis New numerical approximation of fractional derivative with non-local and non-singular kernel: application to chaotic models Based on a specific division, we create collocation points as followsand by substituting them to the model we have following simplified equations given byWe approximate the integrals in the above model using the composite trapezoidal rule. Therefore,By simulating the above equations and as an illustration of the proposed numerical algorithm, we present some graphical illustrations for all variables of the new COVID-19 model in Figures 5, 6, 7 , and 8 . ☒ The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.☐The authors declare the following financial interests/personal relationships which may be considered as potential competing interests:CRediT author statement Mutaz Mohammad: Conceptualization, Methodology, Visualization, Software and coding, Investigation, Supervision, Validation, Writing-Reviewing and Editing. Alex Trounev: Software and coding.