Abstract

To evaluate the probable effects of the COVID-19 outbreak on Saudi Arabia, a novel dynamical model is developed. Using the most recent instances of COVID-19 infection that have been reported in Saudi Arabia, we examine the roles of quarantine and hospitalization. The model’s mathematical outcomes are displayed. The model’s infection-free equilibrium is displayed, and asymptotically, it is determined to be both locally and globally stable. We demonstrate that the model is locally asymptotically stable (LAS) for the basic reproduction . The model is globally asymptotically stable (GAS) when . To estimate the model parameters, recent COVID-19 instances in KSA that began between May 1, 2022, and August 4, 2022, are taken into account. We achieve the needed data fitting considering the approach of nonlinear least square, and we demonstrate that the predicted basic reproduction number is . Graphical representations of the calculated parameters and their effects on disease eradication are provided. The findings show that the most effective way to reduce the number of new instances of infection is to limit the contact of exposed, asymptomatic, symptomatic, and hospitalized people with vulnerable. The percentage of exposed people who are quarantined also plays a big part in lowering the number of infected cases.

1. Introduction

It is believed in epidemiology that mathematical models are essential for comprehending the dynamics and predicting the long-term behaviors of the sickness. Mathematical models are regarded as the main tools for studying the forecasts and eradication of disease outbreaks with the introduction of new infectious illnesses. Due to coronavirus infection, a huge number of death and infection cases were reported worldwide, while many nations in the world experienced financial crisis [1]. With COVID-19, mainly the low-income countries are badly suffered [2]. The coronavirus effect on the Saudi Arabia’s gross domestic product (GDP) can be seen in [3]. Records indicate that a significant number of infected cases in Saudi Arabia ended in death. There have been 9271 fatalities and 812093 total cases documented in the Kingdom to date. The recommendations of the World Health Organization (WHO) were successfully implemented by the Saudi government and as a result, the infected cases were found decreasing significantly, which is less prevalent than the earlier layers or waves of illness [4]. The outbreak of the coronavirus infection has been recorded in many countries of the world, where many of them experienced with more consecutive layers. In Saudi Arabia, there have been previously three waves of infection cases observed, and now, the fourth wave also occurred with a number of infections and deaths. The previous layers recorded more cases of infection and deaths compared to the recent layer.

A variety of mathematical models have been presented in science and engineering to understand the complicated nonlinear process of physical and biological problems [57]. A particular application of the modeling to engineering problems can be seen here [810]. For the application of mathematical modeling to smart phone and image processing, we refer the readers to visit [1113]. To learn more about mathematical models used to research infectious illnesses, and in particular to comprehend the COVID-19’s dynamics, visit [1423]. For instance, the work in [14] has included the lockdown in modeling coronavirus infection. The authors defined new fractal-fractional operators, and their application to coronavirus has been discussed. The researchers suggested from their results that how individuals can be saved with lockdown when there is no vaccination. In [15], a continuous Markov chain model of coronavirus infection is demonstrated. The authors studied the results for the herd immunity and presented results for the disease forecast. The Omicron variant and its mathematical modeling have been suggested in [16]. The authors used the homotopy perturbation method (HPM) and obtained the analytical and approximate solution for the COVID-19 model. The coronavirus infection and its mathematical modeling using the reported cases have been discussed in [17]. The variable order fractional differential equation and its application to COVID-19 have been used and obtained the theoretical results. Using various control measures to curtail the coronavirus infection has been considered in [18]. The authors considered the real data of eight countries and presented the results. In [19], the mathematical analysis to show the infected cases in Ethiopia has been described. They determined the model parameters and computed the threshold quantity. Also, they provided graphical results for possible elimination of the disease in Ethiopia. The mathematical model is designed to incorporate the stress and tension assumptions in the study given in [20]. In [24], the authors considered a fractional model to study COVID-19 with the impact of vaccination. The concept of quarantine together with the asymptomatic individuals in the formulation of COVID-19 has been studied in [25]. The fractional order using the concept of Mittag–Leffler kernel has been considered to investigate COVID-19 with real data in [26]. The COVID-19 disease with the age study has been formulated and discussed in [27]. They analyzed the stability of the model, obtained sensitivity analysis, and provided some simulation results regarding the disease elimination. A mathematical model to understand the infection cases’ trends in KSA has been suggested in [21]. Controlling the outbreak of coronavirus infection, the authors in [28] formulated a fractional model with vaccination. They discussed the role of vaccination in the control of the disease in the society. A numerical study has been organized to solve numerically the COVID-19 system using the Hermite wavelets [29]. The formulation of the coronavirus model in delay differential equations using the fractional order derivative has been discussed in [30]. A study has been organized on coronavirus using the face masks application in [31]. The authors used the real data and studied the impact of the face masks’ efficacy and their uses to control the infection. It has been suggested that with the face masks’ efficacy and their continuous use during the pandemic, the number of potential cases will decrease effectively. The authors used the compartmental model to estimate the parameters using the KSA data and presented the results. The authors explored the coronavirus infection under self-isolation study in [22]. They computed the threshold quantity for the proposed model and presented results for the infection elimination. The researchers in [23] constructed a mathematical model under the early reported cases of coronavirus in KSA. They obtained the numerical value of the threshold and found that the infection can be decreased faster when following the recommendation of the WHO. They studied the local and global dynamics of their proposed model and provided the most sensitive parameters that can increase or decrease the threshold quantity. The anxiety causes during the COVID-19 outbreak to people has been studied in [32], the bifurcation analysis involved in the fractional system in [33], and medical imaging in [34]. Some more applications can be seen in [35, 36].

From the coronavirus literature, we highlight here further research output. The scientists took into account a mathematical model for the instances that were recorded in India and produced the best control strategy for any potential disease eradication [37]. The use of face masks is one of the best controls for COVID-19 infection, according to the research published in [38]. Infection with COVID-19 and cholera has both been documented in [39]. Using a mathematical model from [40], specific applications are applied to Turkish data. Coronavirus with decline in immunity has been considered in [41]. The authors in [42] proposed a mathematical model for coronavirus infection and obtained results regarding the infection elimination. A clinical study-based model about the coinfection of dengue and COVID-19 has been explored in [43]. The diffusion process involved in the coronavirus model has been modeled, and the link between the two has been demonstrated in [44]. A delay differential equations model is suggested in [45] to investigate the COVID-19 disease. Some more work that recently studied the coronavirus and other epidemic diseases through nonlinear differential equations are as follows: a mathematical model to study COVID-19 with the crowding effect has been analyzed in [46]. The authors in [47] used the concept of diffusion to analyze the coronavirus infection. Mathematical modeling of the Nipah virus and analysis using the new numerical scheme has been shown in [48]. The impact of treatment on gonorrhea has been discussed in [49]. Mathematical modeling of Lassa fever in terms of stochastic modeling has been shown in [50]. The dynamics of cancer disease through virotherapy under a new numerical scheme has been considered in [51]. The COVID-19 infection and their control in India have been given in [52], where the authors suggested control strategies for the minimization of the infection. The efficacy of preventive and optimal control analysis for the COVID-19 infection in India has been considered by the authors in [53].

The primary goal of this research is to examine the COVID-19 disease’s dynamics in light of a recent report on cases in Saudi Arabia between May 1 and August 4, 2022. The contact of healthy people with exposed, asymptomatic, symptomatic, and hospitalized people is the four transmission channels that are responsible for the disease spread of coronavirus infection in the population. Such force of infection has not yet been considered in the literature to analyze the COVID-19 infection dynamics with real data. The potential transmission paths for COVID-19 are addressed and pertinent references are given. We developed the model and acquired the requisite mathematical and simulation findings based on these transmissions.

The manuscript will be divided into the following sections: a thorough literature review on COVID-19 and other outcomes of the mathematical modeling have been presented in Section 1. The remaining findings are presented in the following sections. The formulations of the problem are described in Section 2. The model’s mathematical analysis and the stability findings are presented in Sections 3 and 4, respectively. The estimate of the parameters and their numerical outcomes are presented in Section 5. The final findings are presented in Section 6.

2. Model Formulation

This part investigates the development of a mathematical model for the dynamics of the SARS-CoV-2 illness under isolation and quarantine. We consider the model of total human population denoted by and divide it further into susceptible (individuals that can contract the disease), exposed (individuals that attracted the disease and undergoes the disease incubation period), asymptomatic (individuals that do not show disease symptoms of COVID-19), symptomatic (individuals that show disease symptoms of COVID-19), quarantined (individuals that are quarantined during the exposed period), hospitalized individuals (individuals that show clinical symptoms and quarantined individuals which are infected and being hospitalized), and recovered individuals (individuals in different compartments that are recovered from COVID-19 infection). So, we write the total population as , and the force of the infection term is shown bywhere the healthy individuals after getting in close contact with the exposed, asymptomatic, symptomatic, and hospitalized individuals are infected with the contact rates given by , , , and , respectively. The contact between healthy and those hospitalized is shown by , as in many countries of the world with fewer medical facilities [54]. The healthy individuals and their contact with the exposed individuals are possible because the virus can easily transfer to healthy individuals through coughing, shaking hands, flu, etc., and this fact has been documented in the literature, see [5558]. It is well known that individuals in the asymptomatic class have strong immunity and they do not show the disease symptoms until it is tested to be identified or hard to screen for and can easily transmit the infection to healthy people. It has been recorded that most of the COVID-19 cases are from asymptomatic individuals [59, 60]. The interaction of symptomatic people that do not follow the guidelines of COVID-19 and their interaction with vulnerable people increase the number of infections. The hospitalized individuals while facing susceptible individuals in hospitals and other healthcare centers transmit the infections to healthcare workers. It is recommended that healthcare workers should follow the necessary guidelines, such as using the kits, to reduce the disease’s further spread. The growth rate of the susceptible population is given by the parameter ; as opposed to this, provides the natural death rate for each compartment in the model. The vulnerable people after getting an infection undergo an incubation period and then split the individual progress into asymptomatic and symptomatic classes, respectively, with the parameters and . The exposed individuals during their exposed period are quarantined at a rate . Individuals who are asymptomatic, symptomatic, quarantined, and hospitalized are shown to recover in turn by the rate , and . The symptomatic and quarantined people are hospitalized by the rate and , respectively. The parameters and , respectively, represent the disease mortality rate of symptomatic and hospitalized individuals. In our model, we do not give a transition from asymptomatic to symptomatic class as there are no visible facts that asymptomatic develops symptoms or it may be later, so we follow the suggestion of [61]. The abovementioned details of the parameters and their compartments are shown in details in Figure 1. The abovementioned brief discussion leads to the development of the evolutionary nonlinear ordinary differential equations given bywith the non-negative initial conditions (ICs)

The following biological feasible region which is positively invariant is considered for the model (2):and any system trajectory under an IC shall remain in for each time . Since the region is positively invariant, may be used to study the dynamical outcomes.

3. Mathematical Analysis of the Model

This section studies the mathematical analysis of model (2). For a dynamical system, we first need to determine the possible equilibrium points of model (2). Obviously, two equilibrium points that often involve in the diseases model, such as the disease-free equilibrium (DFE) and the endemic equilibrium (EE). The DFE shall be represented by of system (2) and can be obtained through the following way:

3.1. Basic Reproduction Number

The basic reproduction number has an important role in disease epidemiology. It characterizes how the disease spread in the population and how it can be controlled. One can easily define the number of average COVID-19 infection cases introduced to the healthy population who have not yet infected and produced a number of new secondary cases. Usually, in the literature, it is commonly denoted by . For compartmental models of the disease, the approach of the next-generation matrix is often used to obtain the expression for the basic reproduction number [62]. For the given model (2), with the abovementioned approach, we obtain the required matrices as follows:

Furthermore, we need to compute the spectral radius of that shall gives the basic reproduction number of system (2) given by

Here, the expression given by accounts for the generation of secondary infection due to the interaction of the exposed individuals with healthy individuals. The expressions and account for the cases due to hospitalized individuals. The secondary infected cases due to asymptomatic contact are generated through the expression . The secondary cases generated by the symptomatic individuals by contacting the healthy people are shown by . It should be noted that .

3.2. Endemic Equilibria

The endemic equilibrium denoted by can be determined for model (2) by making the time rate of change equal to zero and is given by

Inserting (8) into the force of the infection term givesand we get the following:where

We see that and can be positive if . It can be observed from (11) that a positive endemic equilibrium exists only if .

4. Stability Analysis

For the disease-free situation, both locally and globally, the stability analysis of model (2) shall be explored in the present section. The following theorems are given.

Theorem 1. The COVID-19 system is LAS if and the Routh–Hurtwiz conditions are met.

Proof. We shall compute the Jacobian matrix of system (2) and its evaluation based on the disease-free case , which is computed and is given bywhere . We have two roots for that contains negative real parts , while the rest of the five eigenvalues can be determined through the following equations:whereIt can be observed from the coefficients of equation (13) that all for are positive depending on the value of less than 1. Furthermore, one can establish the Routh–Hurtwiz criteria given byto ensure that equation (13) contains all eigenvalues with negative real parts. So, we can conclude that the infection-free equilibrium of model (2) is locally asymptotically stable provided that .
Next, we give the global stability of the infection-free equilibrium.

Theorem 2. If , then at the infection-free state , system (2) is GAS.

Proof. We obtain the global stability results by defining the Lyapunov function and is given byHere, are constants and positive, and their numerical value in later stepping of the proof shall be obtained. Now inserting system (2) into (16) and simplifying, we getIn the last step of (17), we used . Furthermore, considering the value of the constants , , , , and in (17), we getHere, if , then and if , then . So, system (2) at the infection-free state is globally asymptotically stable on provided that .

5. Parameters Estimations and Numerical Simulation

5.1. Parameters Estimations

Here, we focus on the most recent instances of COVID-19 recorded in Saudi Arabia between May 1 and August 4, 2022; for more information, visit [4]. The actual cases are obtained from Worldometer on daily bases. According to the data details given in [4], it is the fourth wave in the country, which has fewer cases as compared to the other three waves in the country. We used the data of the given period and apply the curve fitting technique and desired results are obtained in time units per day. In this curve fitting analysis, we consider the birth by estimating from the data while the natural death rate has been obtained from the literature. The total population of Saudi Arabia in 2022 has been considered to be . The following are the starting conditions used to acquire the data fitting: , , , , , , and . The infected cases (symptomatic case) , which is the first case of COVID-19 at the beginning of the fourth wave on May 1, 2022. The other initial conditions are subjected to the data fitting. Up until the required desired fitting to the data is attained, the least square curve fitting approach has been used to simulate the data fitting to the model. The results of the data fitting to the model have been shown graphically in Figure 2, whereas the fitted/estimated values of the parameters are given in Tables 1 and 2. Table 2 contains the model parameters fitted to the model with 95% confidence interval (CI). Considering the listed values of the parameters given in Tables 1 and 2, one can estimate the numerical value of the basic reproduction number to be . According to the nature of the data, the present fitting is considered to be good and its estimated parameters can be used to further study the behavior of the model and the possible disease elimination. Figure 3 shows the simulation of versus the model parameters. It can be seen that increase or decrease in these parameters as a function of will increase or decrease the basic reproduction number.

5.2. Numerical Simulation

The numerical values listed in Tables 1 and 2 and the initial conditions of the system variables used in the estimations of parameters have been utilized to obtain simulation results of model (2). We obtain the simulation of model (2) using the sensitive parameters, which are beneficial for disease elimination. In Figure 4, we consider the impact of the parameter on the infected population, and it can be observed that minimizing the contact of the healthy individuals with those who have disease symptoms can best decrease the infected cases further in the community. One of the importance of awareness among people is to be isolated and make social distance from the people who are considered to threat to the spread of the infection.

The impact of the parameter on the infected population is shown in Figure 5. It can be seen that protection from asymptomatic people is critical as it cannot be identified that the individuals are infected unless it is being tested. So, the general guidelines suggested by the World Health Organization (WHO) should be followed strictly to reduce future cases of infection, see Figure 5.

The infected individuals after getting treatment in hospitals can also infect healthcare workers [6567] and a number of deaths have been reported in the aforementioned references. It is also noted from the literature that a number of medical workers has been infected during treating the infected people in hospitals and some have died. So, the proper guidelines suggested for the healthcare workers must be implemented to reduce the death of the healthcare workers and the future spread of the infection, see Figure 6. It can be observed from Figure 6 that proper guidelines followed by the healthcare workers may decrease well the future infection.

We give Figure 7 to see the impact of the parameter on the dynamics of the infected compartments of the model. The parameter is the fraction of individuals; while completing the incubation period, the individuals are separated based on the infection, asymptomatic or symptomatic. The literature suggests that asymptomatic people have great contributions in the disease progress; a high percentage of cases have been reported due to asymptomatic infection. Therefore, if possible, many individuals should be tested and those identified to be asymptomatic must be isolated to reduce future cases. The parameter value of is obtained from data fitting and is reasonable to the literature. By lowering this value, it becomes clear from Figure 7 that there are fewer asymptomatic instances while there are more symptomatic cases, which is the result of .

Similarly, decreasing the value of , the number of infected cases increases for exposed, quarantined, and hospitalized individuals while decreasing for the asymptomatic and symptomatic people, see Figure 8. The role of quarantine of the exposed individuals has been considered useful in many infectious diseases such as SARS and Ebola and also found useful for SARS-CoV-2 infection. Increasing the number of quarantine individuals will minimize the number of future cases.

In Figure 9, we plotted the role of the quarantine rate of exposed individuals shown by , which shows a significant decrease in the future infected cases in different compartments.

6. Conclusion

We formulated and analyzed a new mathematical model for understanding the COVID-19 dynamics with hospitalization and quarantine with the recent data of Saudi Arabia. The real data of the fourth wave which are the recent cases are considered.

We studied the model-related properties in detail. We obtained the equilibrium points of the system. We discussed that the infection-free equilibrium of the model is found to be locally asymptotically stable when the basic reproduction number . We studied the global asymptotical stability of the model at the infection-free state and found it globally asymptotically stable for . The global asymptotical stability of the model ensures the impossibilities of the backward bifurcation in the present COVID-19 infection model. Also, the endemic expressions indicate that there is no possible occurrence of the backward bifurcation for the proposed model. The possibility of the backward bifurcation would lead to that the infection will lost in the community and hence endemic.

In light of the COVID-19 cases that have been reported in Saudi Arabia, we used the least square curve fitting techniques to parameterize the model. We considered some of the parameters from the literature and have cited them in Table 1. The rest of the parameters are fitted to the infected cases in Saudi Arabia. The realistic parameter values obtained through the data fitting have been used to estimate the basic reproduction number . The parameters impact on the model has been shown graphically, and it can be seen that the infection can be minimized well if the number of individuals identified exposed to the virus could be quarantined early. This early quarantine of the exposed individuals will minimize the future cases of the virus in the country. We observed that quarantine and hospitalization are both effective for the eradication of the disease. The impact of the parameters and the model’s infected compartments has been simulated and we found that the variation in these contacts/parameters rarely reduce the disease future cases significantly. It has been concluded from the results that exposed, asymptomatic, hospitalized, and symptomatic individuals’ interaction with healthy individuals increased the future cases of infection. It is better to minimize such contact with healthy individuals to reduce the future risk and increase of the infection. The simulation results indicate that we may certainly decrease the number of infected cases by following the suggestions of the WHO well; as shown in a mathematical model, the contacts such as must be minimized. For COVID-19, quarantine and hospitalization are important in order to control the possible spread of the disease, so in this work, we used both quarantine and hospitalization and also the asymptomatic class, which can be used for any country data to get the result accurately. This model contains all the possible routes of the disease transmission that take place in COVID-19, so it is a perfect model.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The Deanship of Scientific Research (DSR) at King Abdulaziz University (KAU), Jeddah, Saudi Arabia, has funded this Project under grant no. (G: 154-130-1443). Open access funding was enabled and organized by SANLiC Gold.