
To limit the huge damage caused by plastic pollution, major changes need to be made in the food and beverage packaging sector. In this context, a new packaging system for dry-cured ham slices, containing natural antioxidants, was developed; it consists of a reusable polymer tray sealed with disposable polymer film. The life cycle of the packaging was assessed to compare its environmental impacts with a reference disposable packaging system already in use. The life cycle assessment was performed in accordance with the ISO 14040-14044 series; the system was model using the Gabi software and the ILCD PEF method was used to evaluate the impacts. The functional unit chosen was to pack 1000 batches of 4 slices of dry-cured ham in France. Three packaging scenarios were compared: a reference disposable packaging system, incinerated at end of life; the partially reusable packaging, recycled at end of life and the same partially reusable packaging, incinerated at end of life. The study of the relative impacts of each scenario revealed that for the reference packaging, the production of the tray was the highest-impact stage. With the reusable packaging, the highest-impact stages were the thermosealing process and the production of the trays and films. A significant reduction in all impacts was observed with the use of the reusable packaging. Sensitivity analysis was carried out to study the influence of the breakage rate of the tray during reuse and the number of reuse cycles of the tray. Except for freshwater resource depletion, the reusable packaging had lower environmental impacts even with a high tray breakage rate or a low number of reuses. This paper demonstrates the interest of this reusable and recyclable food contact packaging to lower the environmental footprint of packaging; the reuse and recycling stages now need to be tested in real situations for the packaging prototype to confirm the feasibility of the reuse process.
Citation: Joana Beigbeder, Ahmed Allal, Nathalie Robert. Ex-ante life cycle assessment of a partially reusable packaging system for dry-cured ham slices[J]. Clean Technologies and Recycling, 2022, 2(3): 119-135. doi: 10.3934/ctr.2022007
[1] | J. O. Akanni, S. Ajao, S. F. Abimbade, Fatmawati . Dynamical analysis of COVID-19 and tuberculosis co-infection using mathematical modelling approach. Mathematical Modelling and Control, 2024, 4(2): 208-229. doi: 10.3934/mmc.2024018 |
[2] | Yaxin Ren, Yakui Xue . Modeling and optimal control of COVID-19 and malaria co-infection based on vaccination. Mathematical Modelling and Control, 2024, 4(3): 316-335. doi: 10.3934/mmc.2024026 |
[3] | Ihtisham Ul Haq, Nigar Ali, Shabir Ahmad . A fractional mathematical model for COVID-19 outbreak transmission dynamics with the impact of isolation and social distancing. Mathematical Modelling and Control, 2022, 2(4): 228-242. doi: 10.3934/mmc.2022022 |
[4] | Paride O. Lolika, Mlyashimbi Helikumi . Global stability analysis of a COVID-19 epidemic model with incubation delay. Mathematical Modelling and Control, 2023, 3(1): 23-38. doi: 10.3934/mmc.2023003 |
[5] | Ankur Jyoti Kashyap, Arnab Jyoti Bordoloi, Fanitsha Mohan, Anuradha Devi . Dynamical analysis of an anthrax disease model in animals with nonlinear transmission rate. Mathematical Modelling and Control, 2023, 3(4): 370-386. doi: 10.3934/mmc.2023030 |
[6] | Anita T. Kurniawati, Fatmawati, Chidozie W. Chukwu, Windarto, Faishal F. Herdicho . Optimal control of dengue fever model with a logistically growing human population. Mathematical Modelling and Control, 2025, 5(1): 48-60. doi: 10.3934/mmc.2025004 |
[7] | Salma Al-Tuwairqi, Asma Badrah . A qualitative analysis of a model on alpha-synuclein transport and aggregation in neurons. Mathematical Modelling and Control, 2023, 3(2): 104-115. doi: 10.3934/mmc.2023010 |
[8] | Monica Veronica Crankson, Olusegun Olotu, Ayodeji Sunday Afolabi, Afeez Abidemi . Modeling the vaccination control of bacterial meningitis transmission dynamics: a case study. Mathematical Modelling and Control, 2023, 3(4): 416-434. doi: 10.3934/mmc.2023033 |
[9] | Abdul-Fatawu O. Ayembillah, Baba Seidu, C. S. Bornaa . Mathematical modeling of the dynamics of maize streak virus disease (MSVD). Mathematical Modelling and Control, 2022, 2(4): 153-164. doi: 10.3934/mmc.2022016 |
[10] | Eminugroho Ratna Sari, Nikken Prima Puspita, R. N. Farah . Prevention of dengue virus transmission: insights from host-vector mathematical model. Mathematical Modelling and Control, 2025, 5(2): 131-146. doi: 10.3934/mmc.2025010 |
To limit the huge damage caused by plastic pollution, major changes need to be made in the food and beverage packaging sector. In this context, a new packaging system for dry-cured ham slices, containing natural antioxidants, was developed; it consists of a reusable polymer tray sealed with disposable polymer film. The life cycle of the packaging was assessed to compare its environmental impacts with a reference disposable packaging system already in use. The life cycle assessment was performed in accordance with the ISO 14040-14044 series; the system was model using the Gabi software and the ILCD PEF method was used to evaluate the impacts. The functional unit chosen was to pack 1000 batches of 4 slices of dry-cured ham in France. Three packaging scenarios were compared: a reference disposable packaging system, incinerated at end of life; the partially reusable packaging, recycled at end of life and the same partially reusable packaging, incinerated at end of life. The study of the relative impacts of each scenario revealed that for the reference packaging, the production of the tray was the highest-impact stage. With the reusable packaging, the highest-impact stages were the thermosealing process and the production of the trays and films. A significant reduction in all impacts was observed with the use of the reusable packaging. Sensitivity analysis was carried out to study the influence of the breakage rate of the tray during reuse and the number of reuse cycles of the tray. Except for freshwater resource depletion, the reusable packaging had lower environmental impacts even with a high tray breakage rate or a low number of reuses. This paper demonstrates the interest of this reusable and recyclable food contact packaging to lower the environmental footprint of packaging; the reuse and recycling stages now need to be tested in real situations for the packaging prototype to confirm the feasibility of the reuse process.
Coronavirus Disease (COVID-19) is a highly contagious and viral disease that spreads easily from person to person through contact [1]. It can also be contracted through respiratory droplets released when an infected person coughs, sneezes, breathes, sings or talks [2]. The common symptoms of COVID-19 include fever, cough, chills, headache, muscles aches, vomiting and diarrhea [3]. Other symptoms include breathing difficulties, loss of speech and in severe cases, pneumonia, stroke and blood clots have been the most common problems in COVID-19 patients [4]. The disease affects most people with high blood pressure, cancer, heart failure, overweight, obesity, liver disease and weakened immune systems [5]. Evidence from literature shows that old age has high risk of serious illness from COVID-19 and the risk increases with age [6].
The World Health Organization (WHO) declared COVID-19 a public health problem on 30th January 2020 [7] and a pandemic on 11th March 2020 [7,8]. The first case report of COVID-19 in Tanzania was declared in Arusha before spreading to other parts of the country. Then, on 29th April 2020, the Ministry of Health, Community Development, Gender, Elderly and Children reported a total of 480 confirmed cases and 16 deaths from the disease [9].
Public health prevention measures, including banning all large gatherings, limiting the number of people attending burials, physical distancing and wearing masks, were implemented all over the world [10]. Unlike other East African countries, Tanzania did not enforce lock-down, and people were allowed to continue with income-generating activities normally [9]. However, health education campaigns on prevention measures, such as wearing masks, hygiene practices, avoiding public gatherings and keeping physical distance, were intensively encouraged through mass media and social media [11]. Despite these prevention measures, vaccination remains the most powerful tool in preventing the spread of COVID-19 in the population [12]. Globally, WHO reported that more than 116,135,492 and 2,581,976 confirmed cases and deaths, respectively, and 249,160,837 vaccine doses have been administrated worldwide, including Tanzania [9,13]. However, the disease still persists in the population and there are several cases of resurgences globally. Understanding the dynamics of COVID-19 and its control strategies is a substantially important in order to minimize the spread of the disease in the community [14,15].
Mathematical models of disease transmission have been widely studied since the work of Kermack and McKendrick [16], Greenwood and Yule [17], Ross [18], Bernoulli [19], Brownlee [20], Soper [21], Greenwood [22] and also for the details of history in disease modeling (see, [23,24]) and the references therein. Mathematical models with fractional-order differential equations have particulary received greater attention (see, [25,26,27,28,29,30,31]) and are widely used in disease modeling compared to the integer-order derivatives [32,33]. Fractional-order operators are more accurate in modeling dynamical systems than classical operators as fractional-order models allow more degree of freedom [34,35]. It is worth mentioning here that fractional-order operators have more advantages than classical-order operators as fractional-order models adequately capture hereditary properties, long-range interactions and memory effects that exist in many biological systems [32,33,34,35]. In contrast, it has been documented in the literature that models that utilize integer-order derivatives do not adequately capture hereditary properties, and memory effects that exist in many biological systems [33]. Besides, when compared to the classical integer-order models, fractional-order models provide a higher level of precision and give a better fit to the actual data [36].
Recently, Bushnaq et al. [37], Owusu et al. [38], Ahmed et al. [39], Baba et al. [40], Omame et al. [41] and Aslam et al. [42] utilized the fractional-order derivatives to investigate the effect of memory on COVID-19 disease transmission. Mathematical studies of fractional order differential equations in disease modeling are also found in [43,44,45,46,47,48] and the references therein. For instance, Singh et al. [44] proposed and studied a fractional order model using Atangana-Baleanu Caputo sense to investigate the effect of quarantine on the spread of COVID-19 disease in the population. In conclusion the authors mentioned that quarantine of the infected individuals is effective to minimize the spread of disease in the population. Rehman et al. [43] formulated a fractional order model using the Caputo derivative to investigate the dynamics of COVID-19 and dengue co-infections in the population. The authors simulated the graphs of both COVID-19 and dengue co-infection to compare the results in the sense of Caputo, Caputo-Fabrizio and Atangana-Baleanu. From numerical simulations, their results demonstrated that Caputo sense had better results in the form of stability compared to other operators. Anggriani and Beay [47], formulated a mathematical model of COVID-19 to study the impact of self-isolation and hospitalization. The authors performed a global sensitivity analysis of the model using Latin Hypercube sampling and partial correlation coefficient methods. Their results revealed that parameters representing self-isolation and hospitalization have negative relations. Lolika and Mlyashimbi [48] proposed and studied a COVID-19 epidemic model with incubation delay. The authors computed the basic reproduction number and used to establish the conditions for global stability of equilibrium points. Furthermore, the authors concluded that quarantine of asymptomatic and symptomatic individuals have an impact on minimizing the spread of COVID-19 in the population.
Fractional derivatives have many definitions. In this study, we have chosen to utilize the renowned Caputo fractional operator due to the fact that, the ability to use classical initial conditions in the model formulation is the key benefits of the Caputo fractional derivatives compared to other fractional operators [49,50,51,52,53]. It is worthwhile to mention that the Caputo fractional derivative of any constant is zero. Furthermore, the Caputo fractional operator has a singular non-kernel, which is missing in other operators [54,55,56]. Therefore, we proposed and studied a Caputo fractional-order model of COVID-19 disease transmission to assess the effects of vaccination and quarantine in order to minimize the spread of disease in the population.
The rest of the paper is organized as follows: In Section 2, the proposed model and its analytical results are presented. Results discussion are provided in Section 3. Finally, the concluding remarks are presented in Section 4.
In this section, the Caputo fractional-order derivative has been used to define the model differential equations for COVID-19 transmission. The compartments proposed in this study are used to represent the epidemiological status of each human population. The proposed model consists of six sub-divided compartments: susceptible S(t), exposed E(t), vaccinated V(t), infectious I(t), quarantine Q(t) and recovered R(t) human populations. Thus, the total human population is denoted by N(t) which is, N(t)=S(t)+E(t)+V(t)+Iv(t)+Q(t)+R(t). Throughout the article, variables and parameters are assumed to be none-negative and are defined as follows: Λ and ν1 represent the rate of new recruitment and transition rate from susceptible to vaccinated classes respectively. Thus; ϕ represents the efficacy of vaccination for vaccinated susceptible individuals; 1α denotes the average time humans spend in incubation period; ω represents the rate of quarantine in the population; following successful treatment, I(t) infected humans recover from disease after 1κ days; μ and γ represent the human natural mortality and the transition rate of quarantined people to recovered classes, respectively; d represents the death rate of infected humans. Additionally, it is assumed that once the people become aware of COVID-19 transmission, they change their behavior and take precautions, such as hand-washing, wearing masks, keeping social distances and even quarantining themselves. Thus, the parameter ϵ represents the reduction rate of COVID-19 transmission of susceptible individuals due to health education campaigns. Furthermore, it was assumed that β represents the probability of disease transmission following the successful contact rate σ between infected and susceptible individuals.
Our assumptions on the dynamics of COVID-19 in this study are illustrated in figure 1 and the corresponding model differential equations are presented in model (2.1):
{cbDθtS(t)=Λθ−(1−ϵθ)σθβθS(t)I(t)−(νθ1+μθ)S(t),cbDθtV(t)=νθ1S(t)−(1−ϕq)σθβθI(t)V(t)−μθV(t),cbDθtE(t)=(1−ϵθ)σθβθI(t)S(t)+(1−ϕθ)σθβθI(t)V(t)−(αθ+μθ)E(t),cbDθtQ(t)=ωθαθE(t)−(γθ+μθ)Q(t),cbDθtI(t)=(1−ωθ)αθE(t)−(κθ+μθ+dθ)I(t),cbDθtR(t)=γθQ(t)+κθI(t)−μθR(t). | (2.1) |
We begin by introducing the definition of Caputo fractional derivative and state the related theorems (see, [56,57,58,59]) that we will utilize to derive important results in this work.
Definition 2.1. Suppose that θ>0,t>b,θ,b,t∈R, the Caputo fractional derivative is given by:
cbDθtf(t)1Γ(n−θ)∫tbfn(ξ)(t−ξ)θ+1−ndξ,n−1<θ,n∈N. | (2.2) |
Definition 2.2. (Linearity property [56]). Let f(t),g(t):[b,d]→R be such that cbDθtf(t) and cbDθtg(t) exist almost everywhere and let c1,c2∈R. Then, cbDθt(c1f(t))+cbDθt(c2g(t)) exists everywhere, and
cbDθt(c1f(t)+c2g(t))=c1 cbDqtf(t)+c2 cbDθtg(t). | (2.3) |
Definition 2.3. (Caputo derivative of a constant [59]). The fractional derivative for a constant function f(t)=c is zero, that is:
cbDθtc=0. | (2.4) |
Let us consider the following general type of fractional differential equations involving the Caputo derivative:
cbDθtx(t)=f(t,x(t)),θ∈(0,1), | (2.5) |
with initial condition x0=x(t0).
Definition 2.4. (see [56]). The constant x∗ is an equilibrium point of the Caputo fractional dynamic system (2.5) if and only if f(t,x∗)=0.
In what follows, we present an extension of the Lyapunov direct method for Caputo type fractional order for nonlinear systems [56,60].
Theorem 2.1. (Uniform Asymptotic Stability [56,60]). Let x∗ be an equilibrium point for the non-autonomous fractional order system (2.5) and Ω⊂Rn be a domain containing x∗. Let L:[0,∞)×Ω→R be a continuously differentiable function such that:
M1(x)≤N(t,x(t))≤M2(x), |
and:
cbDθtN(t,x(t))≤M3(x), |
for all q∈(0,1) and all x∈Ω, where M1(x), M2(x) and M3(x) are continuous positive definite functions on Ω. Then, the equilibrium point of system (2.5) is uniformly asymptotically stable.
The following theorem summarizes a lemma proved in [56], where a Volterra-type Lyapunov function is obtained for fractional-order epidemic systems.
Lemma 2.1. (see [56]. Let x(⋅) be a continuous and differentiable function with x(t)∈R+. Then, for any time instant t≥b, one has:
cbDθt(x(t)−x∗−x∗lnx(t)x∗)≤(1−x∗x(t)) cbDθtx(t),x∗∈R+,∀θ∈(0,1). |
Theorem 2.2. For the model (2.1), there exists a unique solution in (0,∞); however, the solution is always positive for all values of t≥0 and remains in R6+.
Proof. From the model (2.1), we first show that
R6+= {N(t) ∈R6+:N(t)≥0} is a positive invariant set. Then, we have to demonstrate that each hyper-plane bounding the positive orthant and the vector field points to R6+. Now consider the following: let us assume that there exists a t∗>t0 such that N(t∗)=0, and N(t)<0 for t∈(t∗,t1), where t1 is sufficiently close to t∗, if N(t∗)=0, then we have that,
cbDθtN(t∗)−Λθ>0. This implies that cbDθtN(t)>0 for all t∈[t∗,t1]. The above discussion shows that the three hyper-plane bounding the orthants that is the vector field points to R6+. This shows that all the solutions of the model (2.1) remains positive for all t≥0.
Theorem 2.3. Let Φ(t)=N(t) be the unique solution of the model (2.1) for all t≥0. Then, the solution Φ(t) is bounded above, that is, Φ(t)∈Ω, where Ω which is the feasible region is defined as,
Ω={N(t)∈R6+0≤N(t)≤CN.} |
and its interior denoted by int(Ω) is given by,
int(Ω)={N(t)∈R6+0≤N(t)≤CN}. |
Proof. Here, we prove that the solutions of model (2.1) are bounded for all t≥0. Biologically, the lowest possible value of each state of model (2.1) is zero. Next, we determine the upper-bound of states. Based on this discussion, it is easy to show that the following condition holds for biological relevance of species. 0≤N(t)≤CN. From this condition one gets:
cbDθtN(t)≤Λθ−μθN(t). |
From the Laplace transformation condition one gets:
SθL[N(t)]−Sθ−1N(0)≤ΛθS−μθL[N(t)]. |
Collecting the likely terms we have:
L[N(t)]≤ΛθS−1Sθ+μθ+N(0)Sq−1Sθ+μθ=ΛθSθ−(1+θ)Sθ+μθ+N(0)Sθ−1Sq+μθ. |
Using the inverse Laplace transform we have:
N(t)≤L−1{pθΛθSθ−(1+θ)Sθ+μθ}−N(0)L−1{Sθ−1Sθ+μθ}≤ΛθtθEq,θ+1(−μθ)tθ+N(0)Eθ,1(−μθ)tθ≤ΛθμθtθEθ,θ+1(−μθ)tθ+N(0)Eθ,1(−μθ)tθ≤Max{Λθμθ,N(0)}(tθEθ,θ+1(−μθ)tθ+Eθ,1(−μθ)tθ)=CΓ(1)=CN. |
Where, CN=Max{Λθμθ,N(0)}. Therefore, N(t) is bounded above and this completes the proof.
Since R(t) does not appear in all the equations in model (2.1), it is sufficient to analyze the solutions of model (2.6) for the behavior of model differential equations (2.1).
{cbDθtS(t)=Λθ−(1−ϵθ)σθβθS(t)I(t)−(νθ1+μθ)S(t),cbDθtV(t)=νθ1S(t)−(1−ϕθ)σθβθI(t)V(t)−μθV(t),cbDθtE(t)=(1−ϵθ)σθβθI(t)S(t)+(1−ϕq)σθβθI(t)V(t)−(αθ+μθ)E(t),cbDθtQ(t)=ωθαθE(t)−(γθ+μθ)Q(t),cbDθtI(t)=(1−ωθ)αθE(t)−(κθ+μθ+dθ)I(t). | (2.6) |
In what follows, we compute the threshold quantity, R0 which determines the power of the disease to spread in the population. The model (2.6) always has a disease-free equilibrium, E0 given by:
E0:(S0,V0,E0,Q0,I0,R0)=(Λθνθ1+μθ,νθ1Λθμθ(νθ1+μθ),0,0,0). |
Following the next generation matrix approach as used in [31,61], the non-negative matrix F that denotes the generation of new infections and the non-singular matrix V that denotes the disease transfer among compartments evaluated at E0 are defined as follows:
F=(00(1−ϵθ)σθβθS0+(1−ϕ)σθβθV0000000), | (2.7) |
V=(αθ+μθ00−ωθαθγθ+μθ0−(1−ωθ)αθ0κθ+μθ+dθ). | (2.8) |
Therefore, from (2.7) and (2.8) it can easily be verified that the basic reproduction number R0 of model (2.1) is:
R0=Λθνθ1+μθ(1−ωθ)αθαθ+μθ((1−ϵθ)σθβθ(κθ+μθ+dθ)+1−ϕθκθ+μθ+dθν1μθ). |
The basic reproduction number R0 is defined as the expected number of secondary cases of humans produced in a completely susceptible population by one infected individual during its lifetime as infectious. The terms Λθνθ1+μθ, νθ1μθ1 and (1−ωθ)αθαθ+μθ represent the total life span of humans and the average life span of vaccinated and quarantined individuals respectively.
Our goal in this section is to investigate the global stability of the disease-free equilibrium and the endemic equilibrium of the model (2.6).
Theorem 2.4. If R0<1, the disease free-equilibrium point of the model (2.1) is locally asymptotically stable and unstable if R0>1.
Proof. To prove theorem (2.4), we evaluate the Jacobian matrix of the model (2.6) at the disease-free equilibrium and investigate the behavior of eigenvalues. In what follows, the Jacobian matrix of the model (2.6) evaluated at the disease free-equilibrium is given by:
JDFE=(−(νθ1+μθ)000−(1−ϵθ)σθβθS0νθ1−μq00−(1−ϵθ)σθβθV000−(αθ+μθ)0(1−ϵθ)σθβθS0+(1−ϕθ)σθβθV000ωθαθ−(γθ+μθ)000(1−ωθ)αθ0−(κθ+μθ+dθ)). | (2.9) |
The first three eigenvalues of matrix (2.9) are λ1=−(νθ1+μθ),λθ2=−μθ, and λ3=−(γθ+μθ) which are non-positive. The remaining two eigenvalues are obtained in the following matrix;
M=(−(αθ+μθ)(1−ϵθ)σθβθS0+(1−ϕθ)σθβθV0(1−ωθ)αθκθ+μθ+dθ). | (2.10) |
It follows that we find the characteristic polynomial of the matrix (2.10) which is given as follows:
λ2+(αθ+2μθ+κθ+dθ)λ+(1−R0)=0. | (2.11) |
Since the coefficients of characteristic polynomial (2.11) are all non-negative for R0<1, we conclude that the disease-free equilibrium E0 of the model (2.2) is locally asymptotically stable and this completes the proof.
Theorem 2.5. The disease-free equilibrium E0 is globally asymptotically stable if R0≤1, otherwise is unstable.
Proof. To prove the theorem (2.5), we first evaluate the model (2.6) at the point E0 and this leads to the following model;
{cbDθtS(t)=S(t)(Λθ(1S(t)−1S0)−(1−ϵθ)βθI(t)),cbDθtV(t)=V(t)(σθνθ1(S(t)V(t)−S0V0)−(1−ν1)σθI(t)),cbDθtE(t)=(1−ϵθ)(S0+(S(t)−S0))+(1−ν1)βθI(t)(V0+(V(t)−V0)−(αθ+μθ)E(t)),cbDθtQ(t)=ωθαθE(t)−(γθ+μθ)Q(t),cbDθtI(t)=(1−ωθ)αθE(t)−(κθ+μθ+dθ)I(t). | (2.12) |
In the followings, we consider the following Lyapunov function:
L0(t)={S(t)−S0−S0lnS(t)S0}+{V(t)−V0−V0lnV(t)V0}+Q(t)+E(t)+(1−ωθ)αθ+μθ(1−ωθ)αθI(t). |
Taking the derivative of L0(t) along the model (2.12) and making simplifications, one gets:
cbDθtL0(t)≤Λθ{2−S(t)S0−S0S}−νθ1S0{1+S(t)S0−S(t)S0V0V(t)−V(t)V0}−{γq+μθ}Q+{{1−ωθ}αθ{κθ+μθ+dθ}}{1−ωθ}αq{R0−1}I(t). |
Since all the parameters and variables in system (2.13) are non-negative, it follows that cbDθtL0(t)<0 holds if R0<1. Moreover, cbDθtL0(t)=0 if and only if S(t)=0, V(t)=0, E(t)=0, Q(t)=0, I(t)=0, for all t≥0. Thus, L0(t) is Lyapunov function on Ω. Using Lasalle Invariance principle [62] it implies that every solution of the system (2.6) approaches the disease-free equilibrium E0 as t→∞. Therefore, we conclude that the disease-free equilibrium of system (2.6) is globally asymptotically stable whenever R0≤1. This completes the proof.
Theorem 2.6. The Model (2.6) has endemic equilibrium E∗ point which is globally asymptotically stable for R0>1.
Proof. To prove the theorem (2.6), we consider the following Lyapunov functional:
L1(t)=A1{S(t)−S∗−S∗lnS(t)S∗}+A2{V(t)−V∗−V∗lnV(t)V0}+A3{E(t)−E∗−E∗lnE(t)E∗}+A4{Q(t)−Q∗−Q∗lnQ(t)Q∗}+A5{I(t)−I∗−I∗lnI(t)I∗}. |
Differentiating L1(t) one gets the following:
cbDθtL1(t)≤A1(1−S∗S)cbDθtS(t)+A2(1−V∗V(t))cbDθtV(t)+A3(1−E∗E(t))cbDθtE(t)+A4(1−Q∗Q(t))cbDθtQ(t)+A5(1−I∗I(t))cbDθtI(t). | (2.13) |
In what follows, we substitute (2.1) in (2.13) and get the following:
cbDθtL1(t)≤A1(1−S∗S(t))(Λθ−(1−ϵθ)σθβθS(t)I(t)−(ηθρθ+ϵθ+μθ)S(t))+A2(1−V∗V(t))(νθ1S(t)−(1−ϕθ)σθβθI(t)V(t)−μθV(t))+A3(1−E∗E(t))((1−ϕθ)σθβθI(t)S(t)+(1−ϕθ)σθβθI(t)V(t)−(αθ+μθ)E(t))+A4(1−Q∗Q(t))(ωθαθE(t)−(γθ+μθ)Q(t))+A5(1−I∗I(t))((1−ωθ)αθE(t)−(κθ+μθ+dθ)I(t)). | (2.14) |
Setting the model (2.6) at the endemic equilibrium point,
{νθ1+μq=ΛθS∗−(1−ϵθ)σθβθI∗μθ=νθ1S∗V∗−(1−ν1)θθβθI∗(αθ+μθ)=(1−ϵθ)σθβθS∗I∗E∗+(1−ϕθ)σθI∗V∗E∗,(γθ+μθ)=ωθαθE∗Q∗,(κθ+μθ+dθ)=(1−ωθ)αθE∗I∗. | (2.15) |
By substituting (2.15) in (2.14) and solving the constants Ai for i=1,2,...5, one gets the following after simplifications:
cbDθtL1(t)≤νθ1S∗(2−S(t)S∗−S∗S(t))+(1−ωθ)αθE∗(3−V(t)V∗−S(t)S∗−S(t)S∗V∗V)+(1−ϵθ)σθβθ(3−S(t)S∗−E(t)E∗I∗I(t)−S(t)S∗I(t)I∗E∗E(t))+(1−ϕθ)σθI∗V∗(4−S∗S(t)−S∗S(t)V∗V(t)−E(t)E∗I(t)I∗−I(t)I∗V(t)V∗E∗E(t)). | (2.16) |
Since the arithmetic mean is greater than or equal to the geometrical mean, it follows that, from (2.16) we have the following:
(2−S(t)S∗−S∗S(t))≤0. | (2.17) |
Furthermore, let Φ(z)=1−z−ln(z) for z>0. One can note that Φ(z)≤0 if and only if z=1. Using the aforementioned properties of Φ(z), from (2.16) one can note that:
(3−V(t)V∗−S(t)S∗−S(t)S∗V∗V)=Φ(S(t)S∗V(t)V∗)−V(t)V∗−S(t)S∗≤ln(V(t)V∗)−V(t)V∗+ln(S(t)S∗)−S(t)S∗≤0. | (2.18) |
(3−S(t)S∗−E(t)E∗I∗I−S(t)S∗II∗E∗E(t))=Φ(S(t)S∗E(t)E∗I∗I)−S(t)S∗−E(t)E∗I∗I(t)≤ln(S(t)S∗)−S(t)S∗+ln(E(t)E∗I∗I(t))−E(t)E∗I∗I(t)≤0. | (2.19) |
(4−S∗S(t)−S∗S(t)V∗V(t)−E(t)E∗II∗−I(t)I∗V(t)V∗E∗E(t))=Φ(I(t)I∗V(t)V∗E∗E(t))+Φ(E(t)E∗II∗)−S∗S(t)−S∗S(t)V∗V(t)≤ln(S(t)S∗)−S(t)S∗+ln(S∗S(t)V∗V(t))−S∗S(t)V∗V(t)≤0. | (2.20) |
From (2.17), (2.19), (2.18), and (2.20), one can note that cbDθtL1(t)≤0 whenever R0>1. Therefore, using Lasalle Invariance principle [62], the model (2.1) has a global asymptotically stable equilibrium point for all R0≥1 and this completes the proof.
In this section, we perform the numerical simulations of the model (2.1) to justify the analytical results. Most of the parameter values that are not available in the literature have been estimated. Additionally, for the simulation initial conditions of the model (2.1), was assumed to be S(0)=900,V(0)=0,E(0)=200,Q(0)=0 and I(0)=2. Using the similar concept in [63], the fractional Adam-Bashforth-Moulton scheme for the model (2.1) has the following form:
{S(tn+1)=S0+hθΓ(θ+2)fS(tn+1,Sp(tn+1),Vp(tn+1),Ep(tn+1),Ip(tn+1),Qp(tn+1),Rp(tn+1))+hθΓ(θ+2)n∑m=0am,n+1fS(tm,S(tm),V(tm),E(tm),I(tm),Q(tm),R(tm)),V(tn+1)=V0+hθΓ(θ+2)fV(tn+1,Sp(tn+1),Vp(tn+1),Ep(tn+1),Ip(tn+1),Qp(tn+1),Rp(tn+1))+hqΓ(θ+2)n∑m=0am,n+1fV(tm,S(tm),V(tm),E(tm),I(tm),Q(tm),R(tm)),E(tn+1)=E0+hθΓ(θ+2)fE(tn+1,Sp(tn+1),Vp(tn+1),Eθ(tn+1),Ip(tn+1),Qp(tn+1),Rp(tn+1))+hθΓ(q+2)n∑m=0am,n+1fE(tm,S(tm),V(tm),E(tm),I(tm),Q(tm)),R(tm)I(tn+1)=I0+hθΓ(θ+2)fI(tn+1,Sp(tn+1),Vp(tn+1),Ep(tn+1),Ip(tn+1),Qp(tn+1),Rp(tn+1))+hθΓ(θ+2)n∑m=0am,n+1fI(tm,S(tm),V(tm),E(tm),I(tm),Q(tm),R(tm)),Q(tn+1)=Q0+hθΓ(θ+2)fQ(tn+1,Sp(tn+1),Vp(tn+1),Ep(tn+1),Ip(tn+1),Qp(tn+1),Rp(tn+1))+hθΓ(θ+2)n∑m=0am,n+1fQ(tm,S(tm),V(tm),E(tm),I(tm),Q(tm),R(tm)),R(tn+1)=R0+hθΓ(θ+2)fR(tn+1,Sp(tn+1),Vp(tn+1),Ep(tn+1),Ip(tn+1),Qp(tn+1),Rp(tn+1))+hθΓ(θ+2)n∑m=0am,n+1fR(tm,S(tm),V(tm),E(tm),I(tm),Q(tm),R(tm)). | (3.1) |
Where:
{Sp(tn+1)=S0+1Γθn∑m=0bm,n+1fS(tm,S(tm),V(tm),E(tm),I(tm),Q(tm)),R(tm)Vp(tn+1)=V0+1Γθn∑m=0bm,n+1fV(tm,S(tm),V(tm),E(tm),I(tm),Q(tm)),R(tm)Ep(tn+1)=E0+1Γθn∑m=0bm,n+1fE(tm,S(tm),V(tm),E(tm),I(tm),Q(tm)),R(tm)Ip(tn+1)=I0+1Γθn∑m=0bm,n+1fI(tm,S(tm),V(tm),E(tm),I(tm),Q(tm)),R(tm)Qp(tn+1)=Q0+1Γθn∑m=0bm,n+1fQ(tm,S(tm),V(tm),E(tm),I(tm),Q(tm)),R(tm)Rp(tn+1)=R0+1Γθn∑m=0bm,n+1fR(tm,S(tm),V(tm),E(tm),I(tm),Q(tm)),R(tm). | (3.2) |
In what follows we have:
{fS(tm,S(tm),V(tm),E(tm),I(tm),Q(tm),R(tm))=cbDαtS(t),fV(tm,S(tm),V(tm),E(tm),I(tm),Q(tm),R(tm))=cbDαtV(t),fE(tm,S(tm),V(tm),E(tm),I(tm),Q(tm),R(tm))=cbDαtE(t),fI(tm,S(tm),V(tm),E(tm),I(tm),Q(tm),R(tm))=cbDαtI(t),fQ(tm,S(tm),V(tm),E(tm),I(tm),Q(tm),R(tm))=cbDαtQ(t),fR(tm,S(tm),V(tm),E(tm),I(tm),Q(tm),R(tm))=cbDαtR(t). | (3.3) |
Additionally, the quantities:
{fS(tn+1,Sp(tn+1),Vp(tn+1),Ep(tn+1),Ip(tn+1),Qp(tn+1),Rp(tn+1)),fV(tn+1,Sp(tn+1),Vp(tn+1),Ep(tn+1),Ip(tn+1),Qp(tn+1),Rp(tn+1)),fE(tn+1,Sp(tn+1),Vp(tn+1),Ep(tn+1),Ip(tn+1),Qp(tn+1),Rp(tn+1)),fI(tn+1,Sp(tn+1),Vp(tn+1),Ep(tn+1),Ip(tn+1),Qp(tn+1),Rp(tn+1)),fQ(tn+1,Sp(tn+1),Vp(tn+1),Ep(tn+1),Ip(tn+1),Qp(tn+1),Rp(tn+1)),fR(tn+1,Sp(tn+1),Vp(tn+1),Ep(tn+1),Ip(tn+1),Qp(tn+1),Rp(tn+1)). | (3.4) |
are the derivatives from (3.3) at the point tn+1,n=1,2,3,...m..
Symbol | Description | Value | Units | Source |
Λ | Per capita human recruitment rate | 11826 | day −1 | [64] |
d | Disease induced death rate | day −1 | 0.0413 | [9] |
μ | Natural death rate | day −1 | 0.2 | [9] |
β | Force of infections | day −1 | Fitted | |
α | Incubation period | day −1 | 0.5171 | [65] |
σ | Average per capita contact rate | day −1 | Fitted | |
ϕ | Efficacy rate of vaccination | day −1 | 0:854 | [64] |
ν1 | Vaccination rate | day −1 | 0.0313 | [65] |
γ | Recovery rate | day−1 | 0.362 | [64] |
ϵ | Human awareness rate | day−1 | Vary | Assumed |
κ | Treatment rate of infected individuals | unit-less | Vary | Assumed |
ω | Quarantine rate of suspected humans | unit-less | Vary | Assumed |
In this section, the sensitivity analysis of the model (2.1) has been performed to demonstrate the influence of each parameter on the magnitude of the threshold quantity R0. Most of the parameters used in this study have been drawn from literature and some are estimated using reasonable ranges for the purpose of simulation.
Definition 3.1. (See, [66]) The normalized sensitivity index of R0, which depends on differentiability of parameter, ω is defined as follows:
ΨR0ω=∂R0∂ω×ωR0 | (3.5) |
Figure 2 demonstrates the relationship between the basic reproduction number R0 and the model parameters of the model (2.1). Overall, one can note that the model parameters β, γ, σ and Λ have a positive influence on the R0, that is whenever they are increased, the size of R0 increases. On the other-hand, parameters with negative index values have a negative influence on R0, that is, whenever they increased, the value of R0 decreases.
In this section, we use the daily cases of COVID-19 from Wuhan in China as reported in [67] and estimate the parameters (β,σ) that minimize the deviation of real data from prediction of model system (2.1). The main advantage of fractional-order differential equation is that the order of fractional can be any real positive number, so one can choose the one that has best fit of real data to the model and predict the future evolution of the disease in the population. Therefore, in this study, we use both the least squares and Nelder mead algorithm methods as presented in [71] to fit and estimate the parameters (d,β,θ) of the model (2.1). The real data used in this study are daily reported cases as shown in table (2), and the commutative new infections predicted by the model (2.1) is obtained using the equation (3.6)
cbDθtC(t)=(1−ϵ)σθβθI(t)S(t)+(1−ϕ)σθβθI(t)V(t) | (3.6) |
Day | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
Cases | 6 | 12 | 19 | 25 | 31 | 38 | 44 | 60 |
Day | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 |
Cases | 80 | 131 | 259 | 3839 | 469 | 688 | 776 | 1776 |
Day | 17 | 18 | 19 | 20 | 21 | 22 | 23 | |
Cases | 1460 | 1739 | 1984 | 2101 | 2590 | 2827 | 3233 |
We use the following function to compute the best fitting:
F:R2(β,σ)→R(β,σ) | (3.7) |
where β,σ are parameters such that:
(1) For a given (β,σ), we numerically solve the system (2.1) to get a solution ˆYi(t)=(ˆS,ˆV,ˆE,ˆQ,ˆI,ˆR which is an estimation of daily reported cases Y(t) of COVID-19 from Wuhan in China.
(2) Set t0=1 (the fitting process starts in day 1) and for t=2,3,...,23, corresponding to daily in where data are available, evaluate the computed numerical solution for I(t); that is., ˆI(1), ˆI(2), ˆI(3),.....,ˆI(23).
(3) Compute the root mean square (RMSE) of the difference between ˆI(1),ˆIh(2),....,ˆI(23) and real data. This function F returns the root-mean-square error (RMSE) where:
RMSE=√1n23∑k=1(I(k)−ˆI(k))2, | (3.8) |
(4) Determine a global minimum for the RMSE using Nelder-Mead algorithm. The function F takes values in R2 and returns a positive real number.
Using the formula (3.8), we computed the RMSE that measures the closeness of the model prediction to the real data and was found to be 0.1186. This shows that the system (2.1) has a good fit to the daily reported cases of COVID-19. On performing the fitting process we set the following initial conditions S(0)=9999, V(0)=7990, E(0)=10, I(0)=5, Q(0)=5 R(0)=3 and the model parameters are in Table (1). Note that since the fractional-order is any positive real number θ∈(0,1], one can choose the one that better fits the model to the real data. Based on this assertion, the fractional-order θ were assumed to be 0.59,0.6,0.61,0.62 and 0.63 that had a better fit of model to the real data reported in [67].
The outbreak of the COVID-19 depends on the close contact between infected and susceptible individuals in the community. Vaccination and quarantine are proposed as the most effective control strategies that minimize the spread of the COVID-19 disease in the population. In this study, a fractional-order model for COVID-19 has been proposed and studied to assess the effects of the aforementioned control strategies. To analyze the model, the important threshold parameter, R0, has been computed to investigate the stability analysis of the steady states of the model. The results from the analysis demonstrated that both the disease-free equilibrium and endemic equilibrium points are globally stable whenever the basic reproduction number is less and greater than unity, respectively. A Normalized sensitivity index of the basic reproduction number was been performed to establish the relationship between the threshold and model parameters. Overall, one can note that model parameters with a positive index increased the disease persistence in the population. Furthermore, to support the analytical results, matlab software was used to simulate the proposed model and the results demonstrated that both vaccination and quarantine have the potential to minimize the spread of disease in the population. In particular, the disease can die out from the population by vaccinating 70 susceptible people and quarantine 60 of suspected individuals. Besides, real data of COVID-19 disease from Wuhan city in China has been used to fit and validate the proposed model. From the numerical findings, it can be deduced70 that the model fits well with reported cases of COVID-19 in China. Additional simulation of the model to assess the effect of memory on the spread of COVID-19 disease has been performed and the results demonstrated that memory effect has an influence on the spread of COVID-19 in the population. In future, the proposed model presented in this study will be improved by incorporating time delay and assessing its effect on the spread of COVID-19 disease in the population.
All authors are grateful to their respective institutions for the support during preparation of the manuscript. Paride O. Lolika acknowledges the support from the University of Juba, South Sudan.
The authors declare that they have no conflicts of interest in this paper.
All authors have equal contributions and they read and approved the final version of the paper.
[1] | Horton AA (2022) Plastic pollution: When do we know enough? J Hazard Mater 422: 126885. |
[2] |
Iroegbu AOC, Ray SS, Mbarane V, et al. (2021) Plastic Pollution: A Perspective on Matters Arising: Challenges and Opportunities. ACS Omega 6: 19343-19355. https://doi.org/10.1016/j.jhazmat.2021.126885 doi: 10.1016/j.jhazmat.2021.126885
![]() |
[3] | Ellen MacArthur Foundation, The New Plastics Economy—Rethinking the future of plastics, 2016. Available form: https://ellenmacarthurfoundation.org/the-new-plastics-economy-rethinking-the-future-of-plastics-and-catalysing. https://doi.org/10.1021/acsomega.1c02760. |
[4] |
Phelan AA, Meissner K, Humphrey J, et al. (2022) Plastic pollution and packaging: Corporate commitments and actions from the food and beverage sector. J Clean Prod 331: 129827. https://doi.org/10.1016/j.jclepro.2021.129827 doi: 10.1016/j.jclepro.2021.129827
![]() |
[5] |
Persson L, Carney Almroth BM, Collins CD, et al. (2022) Outside the safe operating space of the planetary boundary for novel entities. Environ Sci Technol 56: 1510-1521. https://doi.org/10.1021/acs.est.1c04158 doi: 10.1021/acs.est.1c04158
![]() |
[6] | Karlsson T, Brosché S, Alidoust M, et al. (2021) Plastic pellets found on beaches all over the world contain toxic chemicals. Available form: https://ipen.org/documents/plastic-pellets-found-beaches-all-over-world-contain-toxic-chemicals. |
[7] |
Wiesinger H, Wang Z, Hellweg S (2021) Deep dive into plastic monomers, additives, and processing aids. Environ Sci Technol 55: 9339-9351. https://doi.org/10.1021/acs.est.1c00976 doi: 10.1021/acs.est.1c00976
![]() |
[8] |
Senko JF, Nelms SE, Reavis JL, et al. (2020) Understanding individual and population-level effects of plastic pollution on marine megafauna. Endanger Species Res 43: 234-252. https://doi.org/10.3354/esr01064 doi: 10.3354/esr01064
![]() |
[9] |
Gall SC, Thompson RC (2015) The impact of debris on marine life. Mar Pollut Bull 92: 170-179. https://doi.org/10.1016/j.marpolbul.2014.12.041 doi: 10.1016/j.marpolbul.2014.12.041
![]() |
[10] |
Kasavan S, Yusoff S, Rahmat Fakri MF, et al. (2021) Plastic pollution in water ecosystems: A bibliometric analysis from 2000 to 2020. J Clean Prod 313: 127946. https://doi.org/10.1016/j.jclepro.2021.127946 doi: 10.1016/j.jclepro.2021.127946
![]() |
[11] |
Kurniawan SB, Said NSM, Imron MF, et al. (2021) Microplastic pollution in the environment: Insights into emerging sources and potential threats. Environ Technol Innov 23: 101790. https://doi.org/10.1016/j.eti.2021.101790 doi: 10.1016/j.eti.2021.101790
![]() |
[12] |
Jaiswal KK, Dutta S, Banerjee I, et al. (2022) Impact of aquatic microplastics and nanoplastics pollution on ecological systems and sustainable remediation strategies of biodegradation and photodegradation. Sci Total Environ 806: 151358. https://doi.org/10.1016/j.scitotenv.2021.151358 doi: 10.1016/j.scitotenv.2021.151358
![]() |
[13] |
Fusi A, Guidetti R, Benedetto G (2014) Delving into the environmental aspect of a Sardinian white wine: From partial to total life cycle assessment. Sci Total Environ 472: 989-1000. https://doi.org/10.1016/j.scitotenv.2013.11.148 doi: 10.1016/j.scitotenv.2013.11.148
![]() |
[14] |
Gazulla C, Raugei M, Fullana-i-Palmer P (2010) Taking a life cycle look at crianza wine production in Spain: Where are the bottlenecks? Int J Life Cycle Assess 15: 330-337. https://doi.org/10.1007/s11367-010-0173-6 doi: 10.1007/s11367-010-0173-6
![]() |
[15] |
Manfredi M, Vignali G (2014) Life cycle assessment of a packaged tomato puree: A comparison of environmental impacts produced by different life cycle phases. J Clean Prod 73: 275-284. https://doi.org/10.1016/j.jclepro.2013.10.010 doi: 10.1016/j.jclepro.2013.10.010
![]() |
[16] |
Pasqualino J, Meneses M, Castells F (2011) The carbon footprint and energy consumption of beverage packaging selection and disposal. J Food Eng 103: 357-365. https://doi.org/10.1016/j.jfoodeng.2010.11.005 doi: 10.1016/j.jfoodeng.2010.11.005
![]() |
[17] |
Louis D, Lombart C, Durif F (2021) Packaging-free products: A lever of proximity and loyalty between consumers and grocery stores. J Retail Consum Serv 60: 102499. https://doi.org/10.1016/j.jretconser.2021.102499 doi: 10.1016/j.jretconser.2021.102499
![]() |
[18] |
Flury M, Narayan R (2021) Biodegradable plastic as an integral part of the solution to plastic waste pollution of the environment. Curr Opin Green Sustain Chem 30: 100490. https://doi.org/10.1016/j.cogsc.2021.100490 doi: 10.1016/j.cogsc.2021.100490
![]() |
[19] |
Silva DAL, Renó GWS, Sevegnani G, et al. (2013) Comparison of disposable and returnable packaging: A case study of reverse logistics in Brazil. J Clean Prod 47: 377-387. https://doi.org/10.1016/j.jclepro.2012.07.057 doi: 10.1016/j.jclepro.2012.07.057
![]() |
[20] |
Coelho PM, Corona B, ten Klooster R, et al. (2020) Sustainability of reusable packaging-Current situation and trends. Resour Conserv Recycl X 6: 100037. https://doi.org/10.1016/j.rcrx.2020.100037 doi: 10.1016/j.rcrx.2020.100037
![]() |
[21] | European Union, Directive 2006/12/EC of the European parliament and of the council of 5 April 2006 on waste. OJEU, 2006. Available form: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32006L0012&from=RO. |
[22] | The European Parliament and the Council of the European Union, Directive (EU) 2018/852 of the European Parliament and of the Council of 30 May 2018 amending Directive 94/62/EC on packaging and packaging waste, 2018. Available form: http://data.europa.eu/eli/dir/2018/852/oj. |
[23] |
Fogt Jacobsen L, Pedersen S, Thøgersen J (2022) Drivers of and barriers to consumers' plastic packaging waste avoidance and recycling—A systematic literature review. Waste Manage 141: 63-78. https://doi.org/10.1016/j.wasman.2022.01.021 doi: 10.1016/j.wasman.2022.01.021
![]() |
[24] |
Hesser F (2015) Environmental advantage by choice: Ex-ante LCA for a new Kraft pulp fibre reinforced polypropylene composite in comparison to reference materials. Compos Part B-Eng 79: 197-203. https://doi.org/10.1016/j.compositesb.2015.04.038 doi: 10.1016/j.compositesb.2015.04.038
![]() |
[25] | European Commission, Circular Plastics Alliance, 2022. Available from: https://ec.europa.eu/growth/industry/strategy/industrial-alliances/circular-plastics-alliance_en. |
[26] | ISO 14040: 2006(en): Environmental management—life cycle assessment—Principles and framework. The International Organization for Standardization, 2006. Available from: https://www.iso.org/obp/ui#iso:std:iso:14040:ed-2:v1:en. |
[27] | ISO 14044: 2006: Environmental management—Life cycle assessment—Requirements and guidelines. The International Organization for Standardization, 2006. Available from: https://www.iso.org/obp/ui/#iso:std:iso:14044:ed-1:v1:en. |
[28] | Camps-Posino L, Batlle-Bayer L, Bala A, et al. (2021) Potential climate benefits of reusable packaging in food delivery services. A Chinese case study. Sci Total Environ 794: 148570. https://doi.org/10.1016/j.scitotenv.2021.148570 |
[29] |
Postacchini L, Mazzuto G, Paciarotti C, et al. (2018) Reuse of honey jars for healthier bees: Developing a sustainable honey jars supply chain through the use of LCA. J Clean Prod 177: 573-588. https://doi.org/10.1016/j.jclepro.2017.12.240 doi: 10.1016/j.jclepro.2017.12.240
![]() |
[30] |
Accorsi R, Cascini A, Cholette S, et al. (2014) Economic and environmental assessment of reusable plastic containers: A food catering supply chain case study. Int J Prod Econ 152: 88-101. https://doi.org/10.1016/j.ijpe.2013.12.014 doi: 10.1016/j.ijpe.2013.12.014
![]() |
[31] |
Arunan I, Crawford RH (2021) Greenhouse gas emissions associated with food packaging for online food delivery services in Australia. Resour Conserv Recy 168: 105299. https://doi.org/10.1016/j.resconrec.2020.105299 doi: 10.1016/j.resconrec.2020.105299
![]() |
[32] |
Tamburini E, Costa S, Summa D, et al. (2021) Plastic (PET) vs bioplastic (PLA) or refillable aluminium bottles—What is the most sustainable choice for drinking water? A life-cycle (LCA) analysis. Environ Res 196: 110974. https://doi.org/10.1016/j.envres.2021.110974 doi: 10.1016/j.envres.2021.110974
![]() |
[33] |
Sazdovski I, Bala A, Fullana-i-Palmer P (2021) Linking LCA literature with circular economy value creation: A review on beverage packaging. Sci Total Environ 771: 145322. https://doi.org/10.1016/j.scitotenv.2021.145322 doi: 10.1016/j.scitotenv.2021.145322
![]() |
[34] | Ligthart TN, Ansems AMM (2007) Single use cups or reusable (coffee) drinking systems: an environmental comparison. Available form: http://www.bekerrecycling.nl/nl/file/20110104172024/3/tno-eng.pdf. |
[35] |
Changwichan K, Gheewala SH (2020) Choice of materials for takeaway beverage cups towards a circular economy. Sustain Prod Consum 22: 34-44. https://doi.org/10.1016/j.spc.2020.02.004 doi: 10.1016/j.spc.2020.02.004
![]() |
[36] |
Gallego-Schmid A, Mendoza JMF, Azapagic A (2019) Environmental impacts of takeaway food containers. J Clean Prod 211: 417-427. https://doi.org/10.1016/j.jclepro.2018.11.220 doi: 10.1016/j.jclepro.2018.11.220
![]() |
[37] |
Cottafava D, Costamagna M, Baricco M, et al. (2021) Assessment of the environmental break-even point for deposit return systems through an LCA analysis of single-use and reusable cups. Sustain Prod Consum 27: 228-241. https://doi.org/10.1016/j.spc.2020.11.002 doi: 10.1016/j.spc.2020.11.002
![]() |
[38] |
Buyle M, Audenaert A, Billen P, et al. (2019) The future of ex-ante LCA? Lessons learned and practical recommendations. Sustain 11: 5456. https://doi.org/10.3390/su11195456 doi: 10.3390/su11195456
![]() |
[39] | Issart A (2019) Potential of natural antioxidants for the stabilization of polymers for food packaging and the development of methods for evaluating their migration[PhD's thesis]. Université de Pau et des Pays de l'Adour, France. |
[40] | Allal AM, Issart A (2021) Patent material for food packaging and method for the preparation thereof. European Patent Office, WO/2020/260829. Available form: https://patentscope.wipo.int/search/en/detail.jsf?docId=WO2020260829. |
[41] | Shonfield P (2008) LCA of management options for mixed waste plastics. Available form: https://silo.tips/download/final-report-lca-of-management-options-for-mixed-waste-plastics. |
[42] |
Krehula LK, Siročić AP, Dukić M, et al. (2012) Cleaning efficiency of poly (ethylene terephthalate) washing procedure in recycling process. J Elastom Plast 45: 429-444. https://doi.org/10.1177/0095244312457798 doi: 10.1177/0095244312457798
![]() |
[43] | Doka G (2003) Life cycle inventories of waste treatment services. Available form: https://www.doka.ch/13_I_WasteTreatmentGeneral.pdf. |
[44] | European Commission, Product environmental footprint category rules guidance—version 6.3, 2018. Available form: https://ec.europa.eu/environment/eussd/smgp/pdf/PEFCR_guidance_v6.3.pdf. |
1. | Ahmed Refaie Ali, Daniyal Ur Rehman, Najeeb Alam Khan, Muhammad Ayaz, Asmat Ara, M. Ijaz Khan, Integrating fractional-order SEI1I2I3QCR model with awareness and non-pharmaceutical interventions for optimal COVID-19 pandemic, 2025, 25, 1471-2288, 10.1186/s12874-024-02452-7 | |
2. | Mohamad Chakroun, Samir Haddad, Wafaa M.R. Shakir, Jinane Sayah, Chadi Kallab, Salah Falou, 2024, Recent Monte Carlo-Based Simulations for Modeling Disease Spread, 979-8-3503-6756-0, 1, 10.1109/ICCA62237.2024.10927993 |
Symbol | Description | Value | Units | Source |
Λ | Per capita human recruitment rate | 11826 | day −1 | [64] |
d | Disease induced death rate | day −1 | 0.0413 | [9] |
μ | Natural death rate | day −1 | 0.2 | [9] |
β | Force of infections | day −1 | Fitted | |
α | Incubation period | day −1 | 0.5171 | [65] |
σ | Average per capita contact rate | day −1 | Fitted | |
ϕ | Efficacy rate of vaccination | day −1 | 0:854 | [64] |
ν1 | Vaccination rate | day −1 | 0.0313 | [65] |
γ | Recovery rate | day−1 | 0.362 | [64] |
ϵ | Human awareness rate | day−1 | Vary | Assumed |
κ | Treatment rate of infected individuals | unit-less | Vary | Assumed |
ω | Quarantine rate of suspected humans | unit-less | Vary | Assumed |
Day | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
Cases | 6 | 12 | 19 | 25 | 31 | 38 | 44 | 60 |
Day | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 |
Cases | 80 | 131 | 259 | 3839 | 469 | 688 | 776 | 1776 |
Day | 17 | 18 | 19 | 20 | 21 | 22 | 23 | |
Cases | 1460 | 1739 | 1984 | 2101 | 2590 | 2827 | 3233 |
Symbol | Description | Value | Units | Source |
Λ | Per capita human recruitment rate | 11826 | day −1 | [64] |
d | Disease induced death rate | day −1 | 0.0413 | [9] |
μ | Natural death rate | day −1 | 0.2 | [9] |
β | Force of infections | day −1 | Fitted | |
α | Incubation period | day −1 | 0.5171 | [65] |
σ | Average per capita contact rate | day −1 | Fitted | |
ϕ | Efficacy rate of vaccination | day −1 | 0:854 | [64] |
ν1 | Vaccination rate | day −1 | 0.0313 | [65] |
γ | Recovery rate | day−1 | 0.362 | [64] |
ϵ | Human awareness rate | day−1 | Vary | Assumed |
κ | Treatment rate of infected individuals | unit-less | Vary | Assumed |
ω | Quarantine rate of suspected humans | unit-less | Vary | Assumed |
Day | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
Cases | 6 | 12 | 19 | 25 | 31 | 38 | 44 | 60 |
Day | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 |
Cases | 80 | 131 | 259 | 3839 | 469 | 688 | 776 | 1776 |
Day | 17 | 18 | 19 | 20 | 21 | 22 | 23 | |
Cases | 1460 | 1739 | 1984 | 2101 | 2590 | 2827 | 3233 |