Loading [MathJax]/jax/element/mml/optable/BasicLatin.js
Research article Special Issues

Performance of concrete paving materials incorporating biomass olive oil waste ash and nano-silica

  • Received: 27 July 2024 Revised: 23 September 2024 Accepted: 18 October 2024 Published: 30 October 2024
  • This study evaluates the utilization of biomass olive oil waste ash (OA) as a concrete paving material. Concrete pavement was produced by replacing a portion of Portland cement with OA at different percentages up to 15%. An additional set of concrete pavement was prepared by incorporating OA with nano-silica (NS) at various contents up to 1.5%. The optimal replacement contents of OA or OA and NS were investigated in terms of workability, compressive strength, strength development rate, and durability. The results showed that the optimal replacement level of OA content was 7.5%. The incorporation of NS with OA increased the optimal replacement level to 15%. The incorporation of NS with OA improved the strength, durability, and workability of all mixes. The utilization of OA with NS at optimal levels can produce concrete pavements. Using NS and OA, approximately 10 % cost savings could be achieved, together with a sustainable, environmentally friendly disposal method of olive oil waste.

    Citation: Hashem Al-Mattarneh, Musab Abuaddous, Rabah Ismail, Ahmad B. Malkawi, Yaser Jaradat, Hamsa Nimer, Mohanad Khodier. Performance of concrete paving materials incorporating biomass olive oil waste ash and nano-silica[J]. AIMS Materials Science, 2024, 11(5): 1035-1055. doi: 10.3934/matersci.2024049

    Related Papers:

    [1] Yan Zhang, Shujing Gao, Shihua Chen . Modelling and analysis of a stochastic nonautonomous predator-prey model with impulsive effects and nonlinear functional response. Mathematical Biosciences and Engineering, 2021, 18(2): 1485-1512. doi: 10.3934/mbe.2021077
    [2] Chunmei Zhang, Suli Liu, Jianhua Huang, Weiming Wang . Stability and Hopf bifurcation in an eco-epidemiological system with the cost of anti-predator behaviors. Mathematical Biosciences and Engineering, 2023, 20(5): 8146-8161. doi: 10.3934/mbe.2023354
    [3] Manoj Kumar Singh, Anjali., Brajesh K. Singh, Carlo Cattani . Impact of general incidence function on three-strain SEIAR model. Mathematical Biosciences and Engineering, 2023, 20(11): 19710-19731. doi: 10.3934/mbe.2023873
    [4] Wenjie Qin, Yue Xia, Yi Yang . An eco-epidemic model for assessing the application of integrated pest management strategies. Mathematical Biosciences and Engineering, 2023, 20(9): 16506-16527. doi: 10.3934/mbe.2023736
    [5] Raimund Bürger, Gerardo Chowell, Elvis Gavilán, Pep Mulet, Luis M. Villada . Numerical solution of a spatio-temporal predator-prey model with infected prey. Mathematical Biosciences and Engineering, 2019, 16(1): 438-473. doi: 10.3934/mbe.2019021
    [6] Xue Liu, Xin You Meng . Dynamics of Bacterial white spot disease spreads in Litopenaeus Vannamei with time-varying delay. Mathematical Biosciences and Engineering, 2023, 20(12): 20748-20769. doi: 10.3934/mbe.2023918
    [7] Yun Li, Hongyong Zhao, Kai Wang . Dynamics of an impulsive reaction-diffusion mosquitoes model with multiple control measures. Mathematical Biosciences and Engineering, 2023, 20(1): 775-806. doi: 10.3934/mbe.2023036
    [8] Bruno Buonomo, Rossella Della Marca, Manalebish Debalike Asfaw . Modelling human response to information in voluntary vaccination behaviour using epidemic data. Mathematical Biosciences and Engineering, 2025, 22(5): 1185-1205. doi: 10.3934/mbe.2025043
    [9] Zixiao Xiong, Xining Li, Ming Ye, Qimin Zhang . Finite-time stability and optimal control of an impulsive stochastic reaction-diffusion vegetation-water system driven by Lˊevy process with time-varying delay. Mathematical Biosciences and Engineering, 2021, 18(6): 8462-8498. doi: 10.3934/mbe.2021419
    [10] Luyao Li, Licheng Fang, Huan Liang, Tengda Wei . Observer-based event-triggered impulsive control of delayed reaction-diffusion neural networks. Mathematical Biosciences and Engineering, 2025, 22(7): 1634-1652. doi: 10.3934/mbe.2025060
  • This study evaluates the utilization of biomass olive oil waste ash (OA) as a concrete paving material. Concrete pavement was produced by replacing a portion of Portland cement with OA at different percentages up to 15%. An additional set of concrete pavement was prepared by incorporating OA with nano-silica (NS) at various contents up to 1.5%. The optimal replacement contents of OA or OA and NS were investigated in terms of workability, compressive strength, strength development rate, and durability. The results showed that the optimal replacement level of OA content was 7.5%. The incorporation of NS with OA increased the optimal replacement level to 15%. The incorporation of NS with OA improved the strength, durability, and workability of all mixes. The utilization of OA with NS at optimal levels can produce concrete pavements. Using NS and OA, approximately 10 % cost savings could be achieved, together with a sustainable, environmentally friendly disposal method of olive oil waste.



    Recently, Whitman and Jayaprakash [1] published a study of a simple, stochastic, agent-based model of influenza infection, Infectious diseases have been and continue to be a major public-health problem [2], disrupting people's quality of life and reducing their chances of survival. The appearance of new diseases and the repetition of mutant strains have added to their massive unhelpful blow. The control of the influenza virus, a foremost international health condition poses a scientific challenge at many levels. Influenza is caused by a virus. According to [3], it is characterized by a severe cytopathogenic respiratory disease that is infectious in nature. Based on matrix protein and nucleoprotein differences, it can be subdivide into A, B and C [4].

    Humans and animals are both infected with Types A and B. The most virulent type, i.e., the Type A virus which is now known as one of the most problematic viruses to tackle [5] is further classified based on the hemagglutinin and neuraminidase proteins located on the virus's surface. Hemagglutinin H1 to H16 and neuraminidase N1 to N9 are the two groups of proteins, the combination of which classifies the influenza subtypes. Influenza A's nomenclature is determined by a mix of hemagglutinin and neuraminidase. The H1N1 virus, for example, is influenza A with both H1 and N1 proteins; its spread is usually not epidemic, and that is why it is difficult to distinguish it from a regular cold [6].

    Another way to characterize them is by the influenza A and B strains [7]. New influenza strains emerge as a result of genetic changes or drift [8]. Antigenic drift is caused by progressive changes in the virus over time. Antibodies have a difficult time recognizing new strains. Antigenic shift, on the other hand, is characterized by fast changes in the virus that result in the emergence of a completely new strain. Influenza A can go through either of the modifications, whereas Influenza B only goes through antigenic drift [9,10].

    Vaccination is incapable of protecting against the new influenza strain. The 2009 H1N1 influenza pandemic showed this. Antivirals are thus required to halt the spread of the influenza outbreak [11]. Recently, incidences of influenza virus resistance have been discovered. The H3N2 virus resistance to aminoadamantanes and the H1N1 virus resistance to oseltamivir are two examples. Resistance is lethal and has the potential to create numerous pandemics in the future [12].

    When compared to the original strain, the transmission rate of a new strain is thought to be quite low. The fact that mutation reduces viral strength is linked to this phenomenon [13]. In the instance of H1N1 influenza virus resistance to oseltamivir, however, it was discovered that these alterations do not always affect virus transmission [14,15].

    For at least the last several centuries, the influenza virus has been responsible for periodic outbreaks of acute febrile illness every 1 to 4 years. The first influenza like sickness outbreak was documented in 1173 and 1174 [16], while the first true epidemic occurred in 1694. [17]. Between 1918 and 1919, the globe was struck by the worst epidemic in recorded history, with an estimated 21 million victims [18]. As indicated by records, it was one of the most terrible events in mankind's set of experiences. Three additional pandemics occurred in the 20th century, namely the 1957 H2N2 pandemic, the 1968 H3N2 pandemic and the 2009 flu A (H1N1) infection (pH1N1) pandemic.

    In the most famous model, a flu strain was found in Mexico and later in the United States of America with a blend of various qualities not recently seen in pig or human flu infection strains [19]. The pandemic was announced in August 2010 after the involvement of different zones [20].

    Young people have the highest rates of influenza infection, while older adults have the highest fatality rates. Mortality and morbidity are particularly high for those with certain high risk medical conditions, such as extreme aging, cardiovascular disorders and metabolic diseases such as diabetes mellitus. During the 2009 pandemic, there was an elevated risk of influenza morbidity and mortality in pregnant women [21]. Furthermore, evidence from prior pandemic and seasonal influenza outbreaks suggests that the risk of influenza complications is higher in the second half of pregnancy than in the first.

    To investigate the influence of environmental conditions on influenza transmission and make the results more realistic, we first developed a stochastic mathematical influenza model. Following that, the factors required for extinction and persistence were analyzed. The threshold of the suggested stochastic influenza model has also been established. When there is tiny or large noise, it plays a critical role in the mathematical models as a backbone [22,23]. Finally, we visualized the numerical simulations using MATLAB.

    In this section, we provide our new stochastic influenza model in the form of differential equations.

    ● The total inhabitant (t) is distributed in four compartment, i.e., St, INt, IRt and Rt, which represent the susceptible and infected people with resistance, infected peoples with non-resistance and recovered people respectively.

    ● The variables and parameters of the proposed stochastic model are non-negative.

    ● We deliberate that the variability of μ and γ is subject to stochastic white noise disturbance, i.e., μμ+σ1B1(t) and γγ+σ2B2; where B1(t)andB2(t) represent the Brownian motion with the property B1(0)=0=B2(0) and the intensities σ21andσ22 are positive.

    Remark 2.1. The deterministic general epidemic study estimates that if Rr<1, a small outburst will arise, and if Rr>1, a large outbreak will occur, infecting a large chunk of the populace. The results are based on the assumption that the community is homogeneous and that individuals mingle evenly. However, if the hypothesis of an evenly mixed society is accepted, this model may not be appropriate in particular situations. When contemplating a tiny population, such as an epidemic outbreak in a daycare center or school, it appears logical to presume that the eventual number of infected will be unpredictable or random. Also, even if Rr>1 and the society is huge, if the outbreak is started by only one (or a few) early infectives, the epidemic may never take off by accident. The formulation of a related stochastic epidemic model is motivated by these two aspects. It allows parameter estimations from disease outbreak data to include standard errors, and the subject of disease extinction is better suited for stochastic models for researching epidemic diseases.

    In light of the above speculations, we established the following new stochastic influenza model;

    dS=(bdSSINαSβIRkIR+1)dt,dIN=(αSIN(μ+d)IN)dtσ1INdB1(t),dIR=(βSIRkIR+1(γ+d)IR)dtσ2IRdB2(t),dR=(μIN+γIRdR)dt+σ1INdB1(t)+σ2IRdB2(t). (2.1)
    Table 1.  Parametric description of the model.
    Symbol Description Value
    α The rate of infection by a non-resistant strain 0.15
    β The rate of infection by a resistant strain 0.10
    γ The rate at which resistant individuals are eliminated from the inhabitants 0.75
    μ Represents the pace at which non-resistant strains are evicted from the population 0.80
    b Denotes recruitment into the susceptible group 1.00
    d The death rate 0.20
    k The effect of mutation on the resistant strain 0.10

     | Show Table
    DownLoad: CSV

    Also, we have the compartment table below:

    Table 2.  Compartments and description.
    Symbol Description Value
    S Susceptible 20
    IN Infected peoples with resistance 2
    IR Infected peoples with non-resistance 2
    R Recovered peoples 1

     | Show Table
    DownLoad: CSV

    The authors of [24] developed the following deterministic model:

    dSdt=bdSαSINSIRβkIR+1,dINdt=αSIN(d+μ)IN,dIRdt=βSIRkIR+1(d+γ)IR,dRdt=μIN+γIRdR, (2.2)

    and

    d(t)=bdt (2.3)

    where t=St+INt+IRt+Rt indicates all the constant residents for bμ and 0=S0+IN0+R0+IR0. Equation (2.3) has the exact solution

    t=edt[0+bdedt] (2.4)

    Also we have 0IN0,0S0,IR00,R00INt0,St0,IRt0andRt0, so the result has a positivity property. If Rr<1 then the model (2.2) will be locally stable and unstable otherwise. Similarly for b=0 the model (2.2) will be globally asymptotically stable.

    In this article, we will establish the computational and analytical aspects of the stochastic influenza model, and we will mathematically correlate our results with the deterministic model for future forecasts by using various parametric variables. As a result, this research can help the local community become more aware of the disease's spread.

    Here, we made the following speculations:

    Rd+={ςRd:0<ςi,d>1}.

    ● Suppose a complete probability space (Ω,,{t}t0,P) with filtration {t}t0 satisfies the usual condition.

    We reflect a common four-dimensional stochastic differential equation for the existence of the solution of our model which is described by (2.1):

    dς(t)=Θ(ς(t),t)dt)+Θ(ς(t),t)dB(t)),fort=t0 (3.1)

    with the initial condition ς(t0)=ς0Rd. By defining the differential operator Ł using (3.1), we get

    Ł=t+4i=1Θi(ς,t)ςi+125i,j=1[ΘT(ς,t)ΘT(ς,t)]ij(2ςiςj). (3.2)

    If Ł acts on the function V=(RdטR+;˜R+), then

    ŁV(ς,t)=Vt(ς,t)+Vς(ς,t)Θ(ς,t)+12trace[ΘT(ς,t)Vςς(ς,t)Θ(ς,t)]. (3.3)

    In this section, our discussion will be on the solution of the stochastic influenza model (2.1).

    Theorem 4.1. There is a unique positive solution (INt,St,IRt,Rt) of the system (2.1) for t0 with (IN0,S0,IR0,R0)R4+, and the solution will be left in R4+, with a probability equal to one.

    Proof. Because (2.1) satisfies the local Lipschitz condition [25], formally for (INt,St,IRt,Rt)R4+, we do have (INt,St,IRt,Rt)R4+ a distinctive local solution on t[0,τe), where τe is the flare-up time. Next, our aim is to show τe= for the global solution of (2.1). Assume 00 is very large so that (IN0,S0,IR0,R0) lies in [10,0]. For 0, define

    τ=inf{t[0,τe):1min{(INt,St,IRt,Rt)} or max{(INt,St,IRt,Rt)}}.

    Let inf (because typically \emptyset is the empty set). Since \tau_\ell is increasing for \ell\rightarrow \infty , let \tau_\infty = lim_{\ell\rightarrow \infty}\tau_\ell ; then, we have \tau_\infty\leq\tau_e almost surely. Next, we need to confirm \tau_\infty = \infty a.s. If this assertion is incorrect, then there exist a constant \mathcal{T} > 0 and \in (0, 1) such that \mathbb{P}\{\tau_\infty\leq \mathcal{T}\} > \in . As a consequence, we have \ell_1\geq \ell_0 such that

    \begin{eqnarray} \mathbb{P}\{\tau_{\ell}\leq \mathcal{T}\}\geq\in , \forall k\geq \ell_1 \quad \text{for}\quad t\leq\tau_\ell. \end{eqnarray} (4.1)

    We outline a \mathbb{C}^2 -function U:\mathbb{R}^4_+\rightarrow \mathbb{R}_+ by using the resulting formulation

    \begin{eqnarray} U(\mathcal{I}_{N_t}, \mathcal{S}_t, \mathcal{I}_{R_t}, \mathcal{R}_t)& = &(\mathcal{S}_t-\underline{c}\tilde{}-\underline{c}\tilde{}ln\frac{\mathcal{S}_t}{\underline{c}\tilde{}})+(\mathcal{I}_{N_t}-(\frac{3}{4}+\frac{1}{4})-ln\mathcal{I}_{N_t})\\&&+(\mathcal{I}_{R_t}-(\frac{3}{4}+\frac{1}{4})-ln\mathcal{I}_{R_t})\\ &+&(\mathcal{R}_t-(\frac{3}{4}+\frac{1}{4})-ln\mathcal{R}_t). \end{eqnarray} (4.2)

    Obviously the function U is non-negative which can follow from {z}-(\frac{3}{4}+\frac{1}{4})-\log z \geq 0\quad \forall \quad z > 0 . Suppose \ell\geq \ell_0 and \mathcal{T}\geq0 are arbitrary. Applying Ito's formula to (4.2) we get

    \begin{eqnarray} dU(\mathcal{I}_{N_t}, \mathcal{S}_t, \mathcal{I}_{R_t}, \mathcal{R}_t) & = &(1-\frac{1}{\mathcal{I}_{N_t}})d\mathcal{I}_{N_t}+\frac{1}{{2}\mathcal{I}_{N_t}^2}(d\mathcal{I}_{N_t})^2+(1-\frac{1}{\mathcal{I}_{R_t}})d\mathcal{I}_{R_t}+\frac{1}{{2}\mathcal{I}_{R_t}^2}(d\mathcal{I}_{R_t})^2\\ &+&(1-\frac{1}{\mathcal{R}_t})d\mathcal{R}_t+\frac{1}{{2}\mathcal{R}_t^2}(d\mathcal{R}_t)^2+(1-\frac{\underline{c}\tilde{}}{\mathcal{S}_t})d\mathcal{S}_t \end{eqnarray} (4.3)
    \begin{eqnarray} & = &{\rm{Ł}}^\star U(\mathcal{I}_{N_t}, \mathcal{S}_t, \mathcal{I}_{R_t}, \mathcal{R}_t)+\sigma_1(\mathcal{I}_{N_t}-\mathcal{S}_t)dB_1(t)+\sigma_2(\mathcal{I}_{R_t}-\mathcal{S}_t)dB_2(t), \end{eqnarray} (4.4)

    where {\rm{Ł}}^\star U:\mathbb{R}^4_+\rightarrow \mathbb{R}_+ is defined by

    \begin{eqnarray} {\rm{Ł}}^\star U(\mathcal{I}_{N_t}, \mathcal{S}_t, \mathcal{I}_{R_t}, \mathcal{R}_t)& = &(1-\frac{\underline{c}\tilde{}}{\mathcal{S}_t})(b-d\mathcal{S}_t-\alpha\mathcal{S}_t\mathcal{I}_{N_t}-\frac{\mathcal{S}_t\beta\mathcal{I}_R}{1+k\mathcal{I}_{R_t}})\\ &+&(1-\frac{1}{\mathcal{I}_{N_t}})(\alpha\mathcal{S}_t\mathcal{I}_{N_t}-(d+\mu)\mathcal{I}_{N_t})+\frac{1}{2}\sigma^2_1\\ &+&(1-\frac{1}{\mathcal{I}_{R_t}})(\frac{\beta\mathcal{S}_t\mathcal{I}_{R_t}}{1+k\mathcal{I}_{R_t}}-(d+r)\mathcal{I}_{R_t})+\frac{1}{2}\sigma^2_2\\ &+&(1-\frac{1}{\mathcal{R}_t})(\mu\mathcal{I}_{N_t}+\gamma\mathcal{I}_{R_t}-d\mathcal{R}_t)+\frac{1}{2}\sigma^2_1+\frac{1}{2}\sigma^2_2\\ & = &b-d\mathcal{S}_t-\alpha\mathcal{S}_t\mathcal{I}_{N_t}-\frac{\beta\mathcal{S}_t\mathcal{I}_{R_t}}{1+k\mathcal{I}_{R_t}}-\frac{\underline{c}\tilde{} b}{\mathcal{S}_t}+\underline{c}\tilde{} d\\ &+&\underline{c}\tilde{}\alpha\mathcal{I}_{N_t}+\frac{\underline{c}\tilde{}\beta\mathcal{I}_{R_t}}{1+k\mathcal{I}_{R_t}}+\alpha\mathcal{S}_t\mathcal{I}_{N_t}+(d+\gamma)\mathcal{I}_{N_t}\\ &-&\alpha\mathcal{S}_t+(d+\mu)+\frac{\mathcal{S}_t\beta\mathcal{I}_{R_t}}{1+k\mathcal{I}_{R_t}}-(d+\gamma)\mathcal{I}_{R_t}+(d+\gamma)\\ &-&\frac{\beta\mathcal{S}_t}{1+k\mathcal{I}_{R_t}}+\mu\mathcal{I}_{N_t}+\gamma\mathcal{I}_{R_t}-d\mathcal{R}_t-\gamma +d-\frac{\mu\mathcal{I}_{N_t}}{\mathcal{R}_t}\\ &+&\sigma^2_1+\sigma^2_2\\ &\leq& b-d\mathcal{S}_t-\alpha\mathcal{S}_t\mathcal{I}_{N_t}+\underline{c}\tilde{} d+\underline{c}\tilde{}\alpha\mathcal{I}_{N_t}+\alpha\mathcal{S}_t\mathcal{I}_{N_t}\\ &+&(d+\gamma)\mathcal{I}_{N_t}-\alpha\mathcal{S}_t+(2d+\mu+\gamma)-d\mathcal{I}_{R_t}\gamma\mathcal{I}_{R_t}\\ &+&\mu\mathcal{I}_{N_t}+\gamma\mathcal{I}_{R_t}-d\mathcal{R}_t-\gamma +d+\sigma^2_1+\sigma^2_2\\ &\leq&b+\underline{c}\tilde{} d+\underline{c}\tilde{}\alpha\mathcal{I}_{N_t}+\alpha\mathcal{S}_t\mathcal{I}_{N_t}+(d+\gamma)\mathcal{I}_{N_t}+(d+\gamma)\\ &+&\gamma\mathcal{I}_{R_t}+\mu\mathcal{I}_{N_t}+d+\sigma^2_1+\sigma^2_2\triangleq B^*. \end{eqnarray} (4.5)

    Thus, we have

    \begin{eqnarray} \bar{E}^*[U(\mathcal{S}_t(\tau_{\jmath}\wedge \rightthreetimes), (\mathcal{I}_{N_t}(\tau_{\jmath}\wedge \rightthreetimes), (\mathcal{I}_{R_t}(\tau_{\jmath}\wedge \rightthreetimes), (\mathcal{R}(\tau_{\jmath}\wedge \rightthreetimes)]&\leq& U(\mathcal{S}_0, \mathcal{I}_{N_0}, \mathcal{I}_{R_0}, \mathcal{R}_0\\ &+&\bar{E}^*[\int_0^{\tau_{\jmath}\wedge \rightthreetimes}\texttt{K}dt]\\ &\leq& U(\mathcal{S}_0, \mathcal{I}_{N_0}, \mathcal{I}_{R_0}, \mathcal{R}_0)+\texttt{K}\rightthreetimes . \end{eqnarray} (4.6)

    We consider that \Omega_\jmath = \{\sigma_\jmath\leq \rightthreetimes\} for all \jmath\geq \jmath_1 and by (2.3), P(\Omega_\jmath)\geq\epsilon . We comment that for every \omega\in\Omega_\jmath there exist at least \mathcal{S}_t(\tau_\jmath, \omega), \mathcal{I}_{N_t}(\tau_\jmath, \omega), \mathcal{I}_{R_t}(\tau_\jmath, \omega)\; and\; \mathcal{R}(\tau_\jmath, \omega) , equaling the value \jmath or \frac{1}{\jmath} then we get that U(\mathcal{S}_t(\tau_\jmath), \mathcal{I}_{N_t}(\tau_\jmath), \mathcal{I}_{R_t}(\tau_\jmath), \mathcal{R}(\tau_\jmath)) is not less than \jmath-1-\log\jmath or (\frac{1}{\jmath})-1+\log\jmath. Accordingly,

    \begin{eqnarray} U(\mathcal{S}_t(\tau_\jmath), \mathcal{I}_{N_t}(\tau_\jmath), \mathcal{I}_{R_t}(\tau_\jmath), \mathcal{R}(\tau_\jmath))\geq\bar{E}^*(\jmath-1-\log\jmath)\wedge (\frac{1}{\jmath})-1+\log\jmath). \end{eqnarray} (4.7)

    From (4.1) and (4.6), we get the following relation

    \begin{eqnarray} U(\mathcal{S}_0, \mathcal{I}_{N_0}, \mathcal{I}_{R_0}, \mathcal{R}_0)+\texttt{K}\rightthreetimes\geq\bar{E}^*[1_{\Omega_\jmath}U(\mathcal{S}_t(\tau_\jmath), \mathcal{I}_{N_t}(\tau_\jmath), \mathcal{I}_{R_t}(\tau_\jmath), \mathcal{R}_t(\tau_\jmath))]\\ \geq\xi[(\jmath-1-\log\jmath)\wedge((\frac{1}{\jmath})-1+\log\jmath)], \end{eqnarray} (4.8)

    where 1_{\Omega_\jmath} denotes the indicator function. We observe that k\rightarrow \infty leads to the ambiguity \infty > U(\mathcal{S}_0, \mathcal{I}_{N_0}, \mathcal{I}_{R_0}, \mathcal{R}_0)+\mathcal{M^\prime}\rightthreetimes = \infty , which implies that \tau_\infty = \infty a.s.

    In this part, we determine when the sickness will be cured and when it will be revived. As a result, the system's (2.1) vital reproduction is demonstrated. Based on the proof in [26], we can deduce the subsequent lemmas:

    Lemma 5.1. Let (\mathcal{S}_t, \mathcal{I}_{N_t}, \mathcal{I}_{R_t}, \mathcal{R}_t) be the solution of the model (2.1) with the initial values given by (\mathcal{S}_0, \mathcal{I}_{N_0}, \mathcal{I}_{R_0}, \mathcal{R}_0)\in\mathbb{R}^4_+ , ; then, \lim_{t\rightarrow \infty}\frac{\mathcal{I}_{N_t}+\mathcal{S}_t+\mathcal{I}_{R_t}+\mathcal{R}_t}{t} = 0 is almost certain.

    In fact, together with the positivity of the solution and the system (2.1), we have that \lim_{t\rightarrow \infty}\frac{\mathcal{S}_t}{t} = 0 \lim_{t\rightarrow \infty}\frac{\mathcal{I}_{N_t}}{t} = 0 , \lim_{t\rightarrow \infty}\frac{\mathcal{I}_{R_t}}{t} = 0 and \lim_{t\rightarrow \infty}\frac{\mathcal{R}_t}{t} = 0 a.s.

    Lemma 6.1. Suppose d > (\frac{1}{2}\sigma^2_1\vee\frac{1}{2}\sigma^2_2) . Assume (\mathcal{S}_t, \mathcal{I}_{N_t}, \mathcal{I}_{R_t}, \mathcal{R}_t) is the solution of the model (2.1) with initial values given by (\mathcal{S}_0, \mathcal{I}_{N_0}, \mathcal{I}_{R_0}, \mathcal{R}_0)\in\mathbb{R}^4_+ ; then,

    \begin{eqnarray} \lim\limits_{t\rightarrow \infty}\frac{\int_0^t \mathcal{I}_R(s)dB_2(s)}{t} = 0, \quad \lim\limits_{t\rightarrow \infty}\frac{\int_0^t \mathcal{I}_N(s)dB_1(s)}{t} = 0. \end{eqnarray} (6.1)

    Let

    \begin{eqnarray} \mathcal{R}_r^\star = \frac{\alpha b}{d(d+\mu+\frac{1}{2}\sigma^2_1)} = \mathcal{R}_r-\frac{\alpha b}{2d(d+\mu)(d+\mu+\frac{1}{2}\sigma^2_1)}\sigma^2_1, \end{eqnarray}

    where

    \begin{eqnarray} \mathcal{R}_r = \frac{\alpha b}{d(d+\mu)} \end{eqnarray} (6.2)

    is the deterministic model's fundamental reproduction number.

    In this segment, we scrutinize the condition for the disappearance of the influenza model (2.1), we lead with the following representation and definition. Let \langle \jmath(t)\rangle = \frac{1}{t}\int_0^t \jmath(r)dr , then the following are the outcomes for the disease's termination.

    Theorem 7.1. Let (\mathcal{S}_t, \mathcal{I}_{N_t}, \mathcal{R}, \mathcal{I}_{R_t}) be the solution of the stochastic influenza model (2.1), with the initial values given by (\mathcal{S}_0, \mathcal{I}_{N_0}, \mathcal{R}_0, \mathcal{I}_{R_0})\in \Omega^\odot . If R_r < 1 then \lim_{t\rightarrow \infty}(\frac{\log\mathcal{I}_{N_t}}{t}) < 0 \quad and \quad \lim_{t\rightarrow \infty}(\frac{\log\mathcal{I}_{R_t}}{t}) < 0 ; almost surely, \mathcal{I}_{N_t}\rightarrow 0 and \mathcal{I}_{R_t}\rightarrow 0 exponentially a.s which means that the disease terminates with a probability of one). Also \lim_{t\rightarrow \infty}\int_0^t\mathcal{S}_t = (\frac{b}{d}) , \quad \lim_{t\rightarrow \infty}\mathcal{I}_{N_t}(t)(t) = 0 , \quad \lim_{t\rightarrow \infty}\mathcal{I}_{R_t}(t) = 0 and \quad \lim_{t\rightarrow \infty}\int_0^t\mathcal{R}_t(t) = 0 .

    Proof. After integrating (2.1), we can get the following system of equations

    \begin{eqnarray} \frac{\mathcal{S}_t-\mathcal{S}_0}{t}& = &b-d\langle\mathcal{S}_t\rangle-\alpha\langle\mathcal{S}_t\mathcal{I}_{N_t}\rangle-\frac{\beta\langle\mathcal{S}_t\mathcal{I}_{R_t}\rangle}{1+k\langle\mathcal{I}_{R_t}\rangle}, \\ \frac{\mathcal{I}_{N_t}-\mathcal{I}_{N_0}}{t}& = &\alpha\langle\mathcal{S}_t\mathcal{I}_{N_t}\rangle-(d+\mu)\langle\mathcal{I}_{N_t}\rangle-\frac{1}{t}\sigma_1\int_0^t\mathcal{I}_{N_t}(r)dB_1(r), \\ \frac{\mathcal{I}_{R_t}-\mathcal{I}_{R_0}}{t}& = &\frac{\beta\langle\mathcal{S}_t\mathcal{I}_{R_t}\rangle}{1+k\langle\mathcal{I}_{R_t}\rangle}-(d+\gamma)\langle\mathcal{I}_{R_t}\rangle-\frac{1}{t}\sigma_2\int_0^t\mathcal{I}_{R_t}(r)dB_2(r), \\ \frac{\mathcal{R}_t-\mathcal{R}_0}{t}& = &\mu\langle\mathcal{I}_{N_t}\rangle+\gamma\langle\mathcal{I}_{R_t}\rangle-d\langle\mathcal{R}\rangle+\frac{1}{t}\sigma_1\int_0^t\mathcal{I}_{N_t}(r)dB_1(r)\\ &+&\frac{1}{t}\sigma_2\int_0^t\mathcal{I}_{R_t}(r)dB_2(r), \end{eqnarray} (7.1)
    \begin{eqnarray} \frac{\mathcal{S}_t-\mathcal{S}_0}{t}+\frac{\mathcal{I}_{N_t}-\mathcal{I}_{N_0}}{t}+\frac{\mathcal{I}_{R_t}-\mathcal{I}_{R_0}}{t}+\frac{\mathcal{R}_t-\mathcal{R}_0}{t}& = &b-d\langle\mathcal{S}_t\rangle-d\langle\mathcal{I}_{N_t}\rangle\\ &-&d\langle\mathcal{I}_{R_t}\rangle-d\langle\mathcal{R}\rangle, \end{eqnarray} (7.2)
    \begin{eqnarray} \langle\mathcal{S}_t\rangle = \frac{b}{d}-\langle\mathcal{I}_{N_t}\rangle-\langle\mathcal{I}_{R_t}\rangle-\langle\mathcal{R}_t\rangle+\Phi(t), \end{eqnarray} (7.3)

    where

    \begin{eqnarray} \Phi(t)& = &-\frac{1}{d}\bigg[\frac{\mathcal{S}_t-\mathcal{S}_0}{t}+\frac{\mathcal{I}_{N_t}-\mathcal{I}_{N_0}}{t}+\frac{\mathcal{I}_{R_t}-\mathcal{I}_{R_0}}{t}+\frac{\mathcal{R}_t-\mathcal{R}_0}{t}\bigg] \end{eqnarray} (7.4)

    Obviously \Phi(t)\rightarrow 0 as t\rightarrow \infty . Applying Ito's formula to the second equation of (2.1) gives

    \begin{eqnarray} d\log\mathcal{I}_{N_t}& = &(\alpha\mathcal{S}_t-(d+\mu)+\frac{1}{2}\sigma^2_1)dt-\sigma_1dB_1(t) \end{eqnarray} (7.5)

    If we Integrate (7.5) from 0 to t and divide by t , we get

    \begin{eqnarray} \frac{\log \mathcal{I}_{N_t}-\log \mathcal{I}_{N_0}}{t}& = &\alpha\langle\mathcal{S}_t\rangle-(d+\mu)+\frac{1}{2}\sigma^2_1-\frac{1}{t}\sigma_1\int_0^tdB_1(r) \end{eqnarray} (7.6)

    Substituting (7.3) in (7.6), we have

    \begin{eqnarray} \frac{\log \mathcal{I}_{N_t}-\log \mathcal{I}_{N_0}}{t}& = &\alpha\bigg[\frac{b}{d}-\langle\mathcal{I}_{N_t}\rangle-\langle\mathcal{I}_{R_t}\rangle-\langle\mathcal{R}_t\rangle+\Phi(t)\bigg]\\ &-&(d+\mu)+\frac{1}{2}\sigma^2_1-\frac{1}{t}\sigma_1\int_0^tdB_1(r)\\ &&\leq\frac{\alpha b}{d}-(d+\mu)+\frac{1}{2}\sigma^2_1-\alpha\langle\mathcal{I}_{N_t}\rangle-\alpha\langle\mathcal{I}_{R_t}\rangle-\\ &&\alpha\langle\mathcal{R}_t\rangle-\frac{1}{t}\sigma_1\int_0^tdB_1(r)+\Phi(t)\\ & = & -(d+\mu-\frac{1}{2}\sigma^2_1)(1-\mathcal{R}^*_r)-\alpha\langle\mathcal{I}_{N_t}\rangle-\alpha\langle\mathcal{I}_{R_t}\rangle-\alpha\langle\mathcal{R}_t\rangle+\daleth, \end{eqnarray} (7.7)

    where

    \begin{eqnarray} \daleth& = &\Phi(t)-\frac{1}{t}\sigma_1\int_0^tdB_1(r). \end{eqnarray} (7.8)

    For \daleth = 0 and t\rightarrow \infty , we have

    \begin{eqnarray} \lim\limits_{t\rightarrow \infty}\sup\frac{\log \mathcal{I}_{N_t}}{t}\leq-(d+\mu-\frac{1}{2}\sigma^2_1)(1-\mathcal{R}^\star_r)-\alpha\langle\mathcal{I}_{N_t}\rangle-\alpha\langle\mathcal{I}_{R_t}\rangle-\alpha\langle\mathcal{R}_t\rangle. \end{eqnarray} (7.9)

    Equation (7.9) implies that

    \begin{eqnarray} \lim\limits_{t\rightarrow \infty}\mathcal{I}_{N_t}& = &0. \end{eqnarray} (7.10)

    Similarly, it may also be proved that

    \begin{eqnarray} \lim\limits_{t\rightarrow \infty}\mathcal{I}_{R_t}& = &0. \end{eqnarray} (7.11)

    Now for \mathcal{S}_t , we have the following from the first equation of the model (2.1)

    \begin{eqnarray} \frac{\mathcal{S}_t-\mathcal{S}_0}{t}& = &b-d\int_0^t\mathcal{S}_tdt-\alpha\int_0^t\mathcal{S}_t\mathcal{I}_{N_t}dt-\int_0^t\frac{\beta\mathcal{S}_t\mathcal{I}_{N_t}}{1+k\mathcal{I}_{R_t}}, \\ d\int_0^t\mathcal{S}_tdt& = &b-\alpha\int_0^t\mathcal{S}_t\mathcal{I}_{N_t}dt-\int_0^t\frac{\beta\mathcal{S}_t\mathcal{I}_{N_t}}{1+k\mathcal{I}_{R_t}}-\frac{\mathcal{S}_t-\mathcal{S}_0}{t}, \\ \int_0^t\mathcal{S}_tdt& = &\frac{b}{d}-\frac{\alpha}{d}\int_0^t\mathcal{S}_t\mathcal{I}_{N_t}dt-\int_0^t\frac{\frac{\beta}{d}\mathcal{S}_t\mathcal{I}_{N_t}}{1+k\mathcal{I}_{R_t}}-\frac{1}{b}(\frac{\mathcal{S}_t-\mathcal{S}_0}{t}). \end{eqnarray} (7.12)

    This implies that \lim_{t\rightarrow \infty}\int_0^t\mathcal{S}_t = \frac{b}{d} . Now from the fourth equation of the system (2.1), it follows that

    \begin{eqnarray} \mathcal{R}_t = e^{-dt}\bigg[\mathcal{R}_0+\int_0^t\mu\mathcal{I}_N(r)e^{dt}+\int_0^t\gamma\mathcal{I}_R(r)e^{dt}\bigg]. \end{eqnarray} (7.13)

    By applying the L'Hospital's rule to the above result, we get

    \begin{eqnarray} \lim\limits_{t\rightarrow \infty}\int_0^t\mathcal{R}dt& = &0, \end{eqnarray} (7.14)

    which completes the proof.

    In this section, we will investigate the necessary conditions for the persistence of the disease.

    Theorem 8.1. Assume d > (\frac{1}{2}\sigma^2_1\vee\frac{1}{2}\sigma^2_2) . Let (\mathcal{S}_t, \mathcal{I}_{N_t}, \mathcal{I}_{R_t}, \mathcal{R}_t) be the solution of the stochastic influenza model (2.1), with the initial values given by (\mathcal{S}_0, \mathcal{I}_{N_0}, \mathcal{I}_{R_0}, \mathcal{R}_0)\in\mathbb{R}^4_+ . If R^\star_r > 1 , then

    \begin{eqnarray} &&\lim\limits_{t\rightarrow \infty}\int_0^t\mathcal{S}(s)ds = \frac{b}{d\mathcal{R}^\star_r} \quad a.s, \\ &&\lim\limits_{t\rightarrow \infty}\int_0^t\mathcal{I}_N(s)ds = \frac{d(d+\mu+\frac{1}{2}\sigma^2_1)}{\alpha(d+\mu)}(\mathcal{R}^\star_r-1)\quad a.s., \\ &&\lim\limits_{t\rightarrow \infty}\int_0^t \mathcal{R}(s)ds = \frac{\mu(d+\mu+\frac{1}{2}\sigma^2_1)}{\alpha(d+\mu)}(\mathcal{R}^\star_r-1)\quad a.s., \\ &&\lim\limits_{t\rightarrow \infty}\int_0^t \mathcal{I}_R(s)ds = \bigg[\frac{d(d+\mu+\frac{1}{2}\sigma^2_1)((d+\mu+\frac{1}{2}\sigma^2_1)-(d+\gamma+\frac{1}{2}\sigma^2_1)\alpha d))}{\alpha d k(d+\gamma(b\alpha-d(d+\mu+\frac{1}{2}\sigma^2_1)))}\bigg](\mathcal{R}^\star_r-1)\quad a.s. \end{eqnarray}

    Proof. If \mathcal{R}^\star_r > 1 , then by (7.9) and Lemmas 5.1 and 5.2 in [27], we have

    \begin{eqnarray} \lim\limits_{t\rightarrow \infty}\int_0^t\mathcal{I}_N(s)ds & = &\frac{\frac{b\alpha}{d}-({d+\mu+\frac{1}{2}\sigma^2_1})}{\frac{\alpha(d+\mu)}{d}} = \frac{d(d+\mu+\frac{1}{2}\sigma^2_1)}{d+\mu}(\mathcal{R}^\star_r-1). \end{eqnarray} (8.1)
    \begin{eqnarray} \lim\limits_{t\rightarrow \infty}\int_0^t\mathcal{S}(s)ds& = &\frac{b}{d}-\frac{d+\mu+\frac{1}{2}\sigma^2_1}{\alpha}(\mathcal{R}^\star_r-1) = \frac{b}{d\mathcal{R}^\star_r}. \end{eqnarray}

    Now from the fourth equation of model (2.1), we obtain

    \begin{eqnarray} \frac{\mathcal{R}_t-\mathcal{R}_0}{t}& = &\frac{\mu}{t}\int_0^t\mathcal{I}_N(s)ds+\frac{\gamma}{t}\int_0^t\mathcal{I}_R(s)ds-\frac{d}{t}\int_0^t\mathcal{R}(s)ds\\ &+&\frac{\sigma_1}{t}\int_0^t\mathcal{I}_N(s)dB_1(s)+\frac{\sigma_2}{t}\int_0^t\mathcal{I}_R(s)dB_1(s), \\ &&\frac{1}{t}\int_0^t\mathcal{R}(s)ds = \frac{\mu}{t}\int_0^t\mathcal{I}_N(s)ds+\eth, \end{eqnarray} (8.2)

    where

    \begin{eqnarray} \eth& = &\frac{\gamma}{t}\int_0^t\mathcal{I}_R(s)ds-\frac{d}{t}\int_0^t\mathcal{R}(s)ds+\frac{\sigma_1}{t}\int_0^t\mathcal{I}_N(s)dB_1(s)\\ &+&\frac{\sigma_2}{t}\int_0^t\mathcal{I}_R(s)dB_1(s)-\frac{\mathcal{R}_t-\mathcal{R}_0}{t}, \end{eqnarray}

    \eth(t) has the property that

    \begin{eqnarray} &&\lim\limits_{t\rightarrow \infty}\eth = 0. \end{eqnarray} (8.3)

    By substituting (8.1) in (8.2), we have

    \begin{eqnarray} \lim\limits_{t\rightarrow \infty}\int_0^t \mathcal{R}(s)ds& = &\frac{\mu(d+\mu+\frac{1}{2}\sigma^2_1)}{\alpha(d+\mu)}(\mathcal{R}^\star_r-1). \end{eqnarray}

    Now from the third equation of the model (2.1), using Ito's formula yields

    \begin{eqnarray} d(ln\mathcal{I}_{R_t}-k\mathcal{I}_{R_t})& = &\bigg[\beta\mathcal{S}_t-(d+\gamma)-k(d+\gamma)\mathcal{I}_{R_t}-\frac{1}{2}\sigma^2_2\bigg] dt+\sigma_2 dB_2(s). \end{eqnarray} (8.4)

    Integrating (8.4) from 0 to t , we have

    \begin{eqnarray} \frac{\ln\mathcal{I}_{R_t}-\mathcal{I}_{R_0}}{t}+k(\frac{\mathcal{I}_{R_t}-\mathcal{I}_{R_0}}{t})& = &\frac{\beta}{t}\int_0^t\mathcal{S}(s)ds-(d+\gamma+\frac{1}{2}\sigma^2_2)\\ &-&\frac{k(d+\gamma)}{t}\int_0^t\mathcal{I}_R(s)ds+\frac{1}{t}\int_0^t\sigma_2dB_2(s)\\ & = &\frac{b\beta}{d\mathcal{R}^\star_r}-\frac{k(d+\gamma)}{t}\int_0^t\mathcal{I}_R(s)ds+\frac{1}{t}\int_0^t\sigma_2dB_2(s)\\ &-&(d+\gamma)-\frac{1}{2}\sigma^2_2\\ \frac{k(d+\gamma)}{t}\int_0^t\mathcal{I}_R(s)ds& = &\frac{b\beta}{d\mathcal{R}^\star_r}+\frac{1}{t}\int_0^t\sigma_2dB_2(s)-(d+\gamma)\\ &-&\frac{1}{2}\sigma^2_2-\bigg[\frac{\ln\mathcal{I}_{R_t}-\mathcal{I}_{R_0}}{t}+k(\frac{\mathcal{I}_{R_t}-\mathcal{I}_{R_0}}{t})\bigg], \end{eqnarray}
    \begin{eqnarray} \frac{1}{t}\int_0^t\mathcal{I}_R(s)ds& = &\frac{1}{k(d+\gamma)}\bigg[\frac{b\beta}{d\mathcal{R}^\star_r}-\frac{1}{2}\sigma^2_2-(d+\gamma)+\curlywedge(s)\bigg], \end{eqnarray} (8.5)

    where \curlywedge(s) = \frac{1}{t}\int_0^t\sigma_2dB_2(s)-(\frac{\ln\mathcal{I}_{R_t}-\mathcal{I}_{R_0}}{t}+k(\frac{\mathcal{I}_{R_t}-\mathcal{I}_{R_0}}{t}) ; \eth(t) has the property that

    \begin{eqnarray} &&\lim\limits_{t\rightarrow \infty}\curlywedge(s) = 0. \end{eqnarray} (8.6)

    Taking the limit of (8.5) and incorporating the value \mathcal{R}^\star_r we have

    \begin{eqnarray} \lim\limits_{t\rightarrow \infty}\int_0^t \mathcal{I}_R(s)ds = \bigg[\frac{d(d+\mu+\frac{1}{2}\sigma^2_1)((d+\mu+\frac{1}{2}\sigma^2_1)-(d+\gamma+\frac{1}{2}\sigma^2_1)\alpha d))}{\alpha d k(d+\gamma(b\alpha-d(d+\mu+\frac{1}{2}\sigma^2_1)))}\bigg](\mathcal{R}^\star_r-1)\quad a.s. \end{eqnarray}

    This completes the proof.

    We have accomplished our analysis of disease extinction and persistence. We will now perform some numerical simulations of (2.1) to illustrate the applicability of our findings. The Milstein technique [28] was used to generate the numerical simulations. Consider the model's discretization equation:

    \begin{eqnarray} \label{scheme} \mathcal{S}_{k+1}& = &\mathcal{S}_k+(b-d\mathcal{S}_k-\alpha\mathcal{S}_k\mathcal{I}_{N_k}-\frac{\beta\mathcal{S}_k\mathcal{I}_{R_k}}{1+k\mathcal{I}_{R_k}})\Delta t, \\ \mathcal{I}_{N_{k+1}}& = &\mathcal{I}_{N_k}+(\alpha\mathcal{S}_k\mathcal{I}_{N_k}-(d+\mu)\mathcal{I}_{N_k})\Delta t-\sigma_1 \mathcal{I}_{N_k} \sqrt{\Delta t}\tau_k-\frac{\sigma_1^2}{2}\mathcal{I}_{N_k}(\tau_k^2-1)\Delta t, \\ \mathcal{I}_{R_{k+1}} & = &\mathcal{I}_{R_k}+(\frac{\beta\mathcal{S}_k\mathcal{I}_{R_k}}{1+k\mathcal{I}_{R_k}}-(d+\gamma)\mathcal{I}_{R_k})\Delta t -\sigma_2 \mathcal{I}_{R_k} \sqrt{\Delta t}\tau_k-\frac{\sigma_2^2}{2}\mathcal{I}_{R_k}(\tau_k^2-1)\Delta t, \\ \mathcal{R}_{k+1}& = &\mathcal{R}_k+(\mu\mathcal{I}_{N_k}+\gamma\mathcal{I}_{R_k}-d\mathcal{R}_k)\Delta t+\sigma_1 \mathcal{I}_{N_k} \sqrt{\Delta t}\tau_k+\frac{\sigma_1^2}{2}\mathcal{I}_{N_k}(\tau_k^2-1)\Delta t, \\ &+&\sigma_2 \mathcal{I}_{R_k} \sqrt{\Delta t}\tau_k+\frac{\sigma_2^2}{2}\mathcal{I}_{R_k}(\tau_k^2-1)\Delta t. \end{eqnarray}

    Here, we shall discuss the graphical description of the model (2.1). In Figure 1, we have illustrated a numerical solution of the model (2.1) and those obtained in comparative studies of all the classes for the white noise values \sigma_1 = \sigma_2 = 0.0, \, 0.05, \, 0.10, \, 0.12 and S(0) = 20, \, I_N(1) = 2, \, I_R(1) = 2, \, R(1) = 1 , b = 1.00, \, d = 0.20, \, \alpha = 0.15, \, \mu = 0.80, \, \beta = 0.10, \, \gamma = 0.75, \ and\; k = 0.10 . In Figure 2, we have illustrated a numerical solution of the model (2.1) and those obtained in comparative studies of all the classes for the white noise values \sigma_1 = \sigma_2 = 0.0, \, 0.05, \, 0.10, \, 0.12 and S(0) = 20, \, I_N(1) = 2, \, I_R(1) = 2, \, R(1) = 1 , b = 1.00, \, d = 0.20, \, \alpha = 0.60, \, \mu = 0.20, \, \beta = 0.10, \, \gamma = 0.75, \ and\; k = 0.10 . We have observed that there is an important role in the dynamics of the values for \alpha and \mu . As \alpha was increased from 0.15 to 0.60 and \mu was decreased from 0.80 to 0.20 , we observed a rapid fall in the susceptible class and I_N increased for the cases with the white noise and without white noise. In this case the recovery rate also increased. The role of white noise has presented a change in the dynamics more accurately. Figure 3 represents the joint solutions for the model (2.1) at zero noise for different values of \mu and \alpha . The left side corresponds to \mu = 0.80 and \alpha = 0.15 while the right side graphs show the joint solution at \mu = 0.15 and \alpha = 0.60 . This change shows a clear difference in the dynamics.

    Figure 1.  Numerical results for the model (2.1) at \sigma_1 = \sigma_2 = 0.0, \, 0.05, \, 0.10, \, 0.12 .
    Figure 2.  Numerical results for the model (2.1) at \sigma_1 = \sigma_2 = 0.0, \, 0.05, \, 0.10, \, 0.12 .
    Figure 3.  Joint solution of the model (2.1) at zero noise.

    In this work, we have explored the dynamic behavior of an \mathcal{S}\mathcal{I}_N\mathcal{I}_R\mathcal{R} influenza stochastic model that ponders the effects of information interference and environmental noise. Information interventions and white noise have been discovered to have a significant impacts on the condition. Hereafter, we present our primary findings.

    We have measured the effects of environmental white noise on the disease. We have shown that {{\underset{\scriptscriptstyle\centerdot}{\text{R}}}}{_r} is a threshold of the model (2.1) for the disease to die out or persist, and that noise strength can change the value of the stochastic reproduction number {{\underset{\scriptscriptstyle\centerdot}{\text{R}}}}{_r} . If {{\underset{\scriptscriptstyle\centerdot}{\text{R}}}}{_r} < 1 , the disease will die out with a probability of one. On the other hand, if {{\underset{\scriptscriptstyle\centerdot}{\text{R}}}}{_r} > 1 , there is a stationary distribution for the model (2.1), which means that the disease will prevail. The discretization approach was used to construct a numerical scheme for the model simulations. The results of the simulations are presented throughout the article in the form of graphs that are divided into three sections. In Figure 2, we have shown a numerical solution of the model (2.1) as well as those obtained in comparative studies of all classes for the white noise values \sigma_1 = \sigma_2 = 0.0, \, 0.05, \, 0.10, \, 0.12 and S(0) = 20, \, I_N(0) = 2, \, I_R(0) = 2, \, R(0) = 1 , b = 1.00, \, d = 0.20, \, \alpha = 0.60, \, \mu = 0.20, \ and\; \beta = 0.10 . We have seen that the dynamics of the values for \alpha and \mu play an essential role. We noticed a rapid fall in the susceptible class as the \alpha value was increased from 0.15 to 0.60 and \mu dropped from 0.80 to 0.20 ; additionally, I_N was increased for both instances with and without white noise. The healing rate was also boosted in this instance. White noise has played a larger role in influencing the dynamics.

    The joint solutions for the model (2.1) at zero noise are shown in figure 3 for different values of \mu and \alpha . The left side graphs correspond to \mu = 0.80 and \alpha = 0.15 , whereas the right side graphs correspond to \mu = 0.15 and \alpha = 0.60 for the combined solution. This shift reveals a significant shift in the dynamics.

    J. Alzabut is thankful to Prince Sultan University and OSTİM Technical University for their endless support. G. Alobaidi was supported by Faculty research grant from the American University of Sharjah (Project number FRG21-S-S05).

    The authors declare there is no conflict of interest.



    [1] Jaradat Y, Matalkah F (2021) Olive biomass ash-based geopolymer composite: Development and characterisation. Adv Appl Ceram 120: 1–9. https://doi.org/10.1080/17436753.2020.1839848 doi: 10.1080/17436753.2020.1839848
    [2] Ribeiro TB, Campos D, Oliveira A, et al. (2021) Study of olive pomace antioxidant dietary fibre powder throughout gastrointestinal tract as multisource of phenolics, fatty acids and dietary fibre. Food Res Int 142: 110032. https://doi.org/10.1016/j.foodres.2020.110032 doi: 10.1016/j.foodres.2020.110032
    [3] Rueda MP, Domínguez-Vidal A, Llorent-Martínez EJ, et al. (2024) Monitoring organic matter transformation of olive oil production residues in a full-scale composting plant by fluorescence spectroscopy. Environ Technol Inno 35: 103695. https://doi.org/10.1016/j.eti.2024.103695 doi: 10.1016/j.eti.2024.103695
    [4] Tsoupras A, Panagopoulou E, Kyzas GZ (2024) Olive pomace bioactives for functional foods and cosmetics. AIMS Agric Food 9: 743–766. https://doi.org/10.3934/agrfood.2024040 doi: 10.3934/agrfood.2024040
    [5] Pachaiappan R, Cornejo-Ponce L, Rajendran R, et al. (2022) A review on biofiltration techniques: Recent advancements in the removal of volatile organic compounds and heavy metals in the treatment of polluted water. Bioengineered 13: 8432–8477. https://doi.org/10.1080/21655979.2022.2050538 doi: 10.1080/21655979.2022.2050538
    [6] Arvanitoyannis IS, Kassaveti A (2007) Current and potential uses of composted olive oil waste. Int J Food Sci Tech 42: 281–295. https://doi.org/10.1111/j.1365-2621.2006.01211.x doi: 10.1111/j.1365-2621.2006.01211.x
    [7] Keskes MA, Zouari A, Houssin R, et al. (2022) An overview on olive oil waste valorization scenarios: Life cycle approach. IFAC-Pap 55: 1477–1482. https://doi.org/10.1016/j.ifacol.2022.09.599 doi: 10.1016/j.ifacol.2022.09.599
    [8] Capablo J, Jensen PA, Pedersen KH, et al. (2009) Ash properties of alternative biomass. Energy Fuels 23: 1965–1976. https://doi.org/10.1021/ef8008426 doi: 10.1021/ef8008426
    [9] Alyami M, Hakeem IY, Amin M, et al. (2023) Effect of agricultural olive, rice husk and sugarcane leaf waste ashes on sustainable ultra-high-performance concrete. J Build Eng 72: 106689. https://doi.org/10.1016/j.jobe.2023.106689 doi: 10.1016/j.jobe.2023.106689
    [10] Ismail R, Dahim M, Jaradat A, et al. (2021) Field dielectric sensor for soil pollution application. IOP Conf Ser Earth Environ Sci 801: 012003. https://doi.org/10.1088/1755-1315/801/1/012003 doi: 10.1088/1755-1315/801/1/012003
    [11] Al-Mattarneh H, Dahim M (2021) Comparison of nondestructive testing method for strength prediction of asphalt concrete material. Civ Eng J 7: 165–178. https://doi.org/10.28991/cej-2021-03091645 doi: 10.28991/cej-2021-03091645
    [12] Rosario RD, Cruz ADL, Guzman MPD (2024) A review of biomineralization as solution for roads and infrastructures concrete sustainability. Civ Eng J 10: 2745–2760. http://dx.doi.org/10.28991/CEJ-2024-010-08-020 doi: 10.28991/CEJ-2024-010-08-020
    [13] Tayeh BA, Hadzima-Nyarko M, Zeyad AM, et al. (2021) Properties and durability of concrete with olive waste ash as a partial cement replacement. Adv Concr Constr 11: 59–71. https://doi.org/10.12989/acc.2021.11.1.059 doi: 10.12989/acc.2021.11.1.059
    [14] Malkawi AB, Nuruddin MF, Fauzi A. et al. (2017) Effect of plasticizers and water on properties of HCFA geopolymers. Key Eng Mater 733: 76–79. https://doi.org/10.4028/www.scientific.net/KEM.733.76 doi: 10.4028/www.scientific.net/KEM.733.76
    [15] Nuruddin MF, Malkawi AB, Fauzi A, et al. (2016) Effects of alkaline solution on the microstructure of HCFA geopolymers, In: Wan Abdullah Zawawi NA, Engineering Challenges for Sustainable Future, London: CRC Press.
    [16] Alkheder S, Obaidat YT, Taamneh M (2016) Effect of olive waste (husk) on behavior of cement paste. Case Stud Constr Mat 5: 19–25. https://doi.org/10.1016/j.cscm.2016.05.001 doi: 10.1016/j.cscm.2016.05.001
    [17] Aras U, Kalaycıoğlu H, Yel H, et al. (2022) Utilization of olive mill solid waste in the manufacturing of cement-bonded particleboard. J Build Eng 49: 104055. https://doi.org/10.1016/j.jobe.2022.104055 doi: 10.1016/j.jobe.2022.104055
    [18] Mohamed AM, Tayeh BA, Aisheh YIA, et al. (2023) Utilising olive-stone biomass ash and examining its effect on green concrete: A review paper. J Mater Res Technol 24: 7091–7107. https://doi.org/10.1016/j.jmrt.2023.05.039 doi: 10.1016/j.jmrt.2023.05.039
    [19] Thomas BS, Yang J, Mo KH, et al. (2021) Biomass ashes from agricultural wastes as supplementary cementitious materials or aggregate replacement in cement/geopolymer concrete: A comprehensive review. J Build Eng 40: 102332. https://doi.org/10.1016/j.jobe.2021.102332 doi: 10.1016/j.jobe.2021.102332
    [20] Hakeem IY, Agwa IS, Tayeh BA, et al. (2022) Effect of using a combination of rice husk and olive waste ashes on high-strength concrete properties. Case Stud Constr Mat 17: 01486. https://doi.org/10.1016/j.cscm.2022.e01486 doi: 10.1016/j.cscm.2022.e01486
    [21] Wang C, Zhang P, Guo J (2023) Durability and microstructure of cementitious composites under the complex environment: Synergistic effects of nano-SiO2 and polyvinyl alcohol fiber. Constr Build Mater 400: 132621. https://doi.org/10.1016/j.conbuildmat.2023.132621 doi: 10.1016/j.conbuildmat.2023.132621
    [22] Zhang XH, Ahmad J, Jebur YM, et al. (2024) A review on partial substitution of nanosilica in concrete. Rev Adv Mater Sci 63: 20230157. https://doi.org/10.1515/rams-2023-0157 10.1515/rams-2023-0157 doi: 10.1515/rams-2023-015710.1515/rams-2023-0157
    [23] Althoey F, Zaid O, Martínez-García R, et al. (2023) Impact of nano-silica on the hydration, strength, durability, and microstructural properties of concrete: A state-of-the-art review. Case Stud Constr Mat 18: e01997. https://doi.org/10.1016/j.cscm.e01997
    [24] Sun Y, Zhang P, Guo J (2024) Rheological properties and workability of PVA fiber and nano-SiO2 modified cement-based materials. Dev Built Environ 18: 100396. https://doi.org/10.1016/j.dibe.2024.100396 doi: 10.1016/j.dibe.2024.100396
    [25] Zhang P, SunY, Wu J, et al. (2023) Mechanical properties and microstructure of nano-modified geopolymer concrete containing hybrid fibers after exposure to elevated temperature. Constr Build Mater 409: 134044. https://doi.org/10.1016/j.conbuildmat.2023.134044 doi: 10.1016/j.conbuildmat.2023.134044
    [26] Zhang X, Zhang P, Yuan W, et al. (2023) Durability prediction of geopolymer mortar reinforced with nanoparticles and PVA fiber using particle swarm optimized BP neural network. Nanotechnol Rev 13: 1. https://doi.org/10.1515/ntrev-2023-0214 doi: 10.1515/ntrev-2023-0214
    [27] Chen JJ, Ng PL, Xu L, et al. (2024) Use of nano-silica sol in concrete: Performance and influence mechanisms. Constr Build Mater 411: 134582. https://doi.org/10.1016/j.conbuildmat.2023.134582 doi: 10.1016/j.conbuildmat.2023.134582
    [28] Dahim M, Abuaddous M, Ismail R, et al. (2021) Using a dielectric capacitance cell to determine the dielectric properties of pure sand artificially contaminated with Pb, Cd, Fe, and Zn. Appl Environ Soil Sci 2021: 8838054. https://doi.org/10.1155/2020/8838054 doi: 10.1155/2020/8838054
    [29] Zain MFM, Karim MR, Islam MN, et al. (2015) Prediction of strength and slump of silica fume incorporated high-performance concrete. Asian J Sci Res 8: 264–277. https://doi.org/10.3923/ajsr.2015.264.277 doi: 10.3923/ajsr.2015.264.277
    [30] Mohammed BS, Nuruddin MF, Aswin M, et al. (2016) Structural behavior of reinforced self-compacted engineered cementitious composite beams. Adv Mater Sci Eng 2016: 1–12. https://doi.org/10.1155/2016/5615124 doi: 10.1155/2016/5615124
    [31] Abdullahi M, Al-Mattarneh H, Mohammed B (2009) Equations for mix design of structural lightweight concrete. Eur J Sci Res 31: 132–141.
    [32] Abdullahi M, Al-Mattarneh H, Mohammed B (2009) Statistical modeling of lightweight concrete mixtures. Eur J Sci Res.
    [33] Malkawi AB, Habib M, Alzubi Y, et al. (2020) Engineering properties of lightweight geopolymer concrete using palm oil clinker aggregate. Int J Geomate 18: 132–139. https://geomatejournal.com/geomate/article/view/414
    [34] Hanjitsuwan S, Chindaprasirt P, Pimraksa K (2011) Electrical conductivity and dielectric property of fly ash geopolymer pastes. Int J Miner Metall Mater 18: 94–99. https://doi.org/10.1007/s12613-011-0406-0 doi: 10.1007/s12613-011-0406-0
    [35] Nuruddin MF, Malkawi AB, Fauzi A, et al. (2016) Geopolymer concrete for structural use: Recent findings and limitations. IOP Conf Ser Mater Sci Eng 133: 012021. https://doi.org/10.1088/1757-899X/133/1/012021 doi: 10.1088/1757-899X/133/1/012021
    [36] Cabrera M, Martinez-Echevarria MJ, López-Alonso M, et al. (2021) Self-compacting recycled concrete using biomass bottom ash. Materials 14: 6084. https://doi.org/10.3390/ma14206084 doi: 10.3390/ma14206084
    [37] Mahameid AA, Yasin AA, Malkawi AB (2024) Prediction of the compressive and tensile strengths of geopolymer concrete using artificial neural networks. Civ Eng Archit 12: 3193–3206. https://doi.org/10.13189/cea.2024.120506 doi: 10.13189/cea.2024.120506
    [38] Rojo-López G, Nunes S, González-Fonteboa B, et al. (2020) Fernando Martínez-Abella, Quaternary blends of portland cement, metakaolin, biomass ash and granite powder for production of self-compacting concrete. J Clean Prod 266: 121666. https://doi.org/10.1016/j.jclepro.2020.121666 doi: 10.1016/j.jclepro.2020.121666
    [39] Abuaddous M, Dahim M, Ismail R, et al. (2021) Sustainable asphalt concrete for road construction and building material. IOP Conf Ser Earth Environ Sci 801: 012023. https://doi.org/10.1088/1755-1315/801/1/012023 doi: 10.1088/1755-1315/801/1/012023
    [40] Woszuk A, Bandura L, Franus W (2019) Fly ash as low cost and environmentally friendly filler and its effect on the properties of mix asphalt. J Clean Prod 235: 493–502. https://doi.org/10.1016/j.jclepro.2019.06.353 doi: 10.1016/j.jclepro.2019.06.353
    [41] Dahim M, Abuaddous M, Al-Mattarneh H, et al. (2021) Enhancement of road pavement material using conventional and nano-crude oil fly ash. Appl Nanosci 11: 2517–2524. https://doi.org/10.1007/s13204-021-02103-z doi: 10.1007/s13204-021-02103-z
    [42] El Boukhari M, Merroun O, Maalouf C, et al. (2023) Mechanical performance of cement mortar with olive pomace aggregates and olive mill wastewater: An experimental investigation. Cogent Eng 10: 2212522. https://doi.org/10.1080/23311916.2023.2212522 doi: 10.1080/23311916.2023.2212522
    [43] Lila K, Belaadi S, Solimando R, et al. (2020) Valorisation of organic waste: Use of olive kernels and pomace for cement manufacture. J Clean Prod 277: 123703. https://doi.org/10.1016/j.jclepro.2020.123703 doi: 10.1016/j.jclepro.2020.123703
    [44] Antoun M, Issa CA, Aouad G, et al. (2021) Sustainable masonry blocks: Olive wood waste as substitute for fine aggregates. Case Stud Constr Mat 15: e00590. https://doi.org/10.1016/j.cscm.2021.e00590
    [45] ASTM International (2012) Standard Specification for Portland Cement, West Conshohocken, PA, ASTM C150-07.
    [46] ASTM International (2016) Standard Specification for Concrete Aggregates, West Conshohocken, PA, ASTM C33/C33M-16e1.
    [47] British Standards Institution (2000) Testing fresh concrete—Part 2: Slump test.
    [48] British Standards Institution (2019) Testing hardened concrete—Part 3: Compressive strength of test specimens.
    [49] British Standards Institution (2020) Testing concrete—Method for determination of water absorption.
  • This article has been cited by:

    1. Wenxuan Li, Suli Liu, Dynamic analysis of a stochastic epidemic model incorporating the double epidemic hypothesis and Crowley-Martin incidence term, 2023, 31, 2688-1594, 6134, 10.3934/era.2023312
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1159) PDF downloads(44) Cited by(11)

Figures and Tables

Figures(16)  /  Tables(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog