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

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

  • 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] Maryam Basiri, Frithjof Lutscher, Abbas Moameni . Traveling waves in a free boundary problem for the spread of ecosystem engineers. Mathematical Biosciences and Engineering, 2025, 22(1): 152-184. doi: 10.3934/mbe.2025008
    [2] Tracy L. Stepien, Erica M. Rutter, Yang Kuang . A data-motivated density-dependent diffusion model of in vitro glioblastoma growth. Mathematical Biosciences and Engineering, 2015, 12(6): 1157-1172. doi: 10.3934/mbe.2015.12.1157
    [3] Li-Jun Du, Wan-Tong Li, Jia-Bing Wang . Invasion entire solutions in a time periodic Lotka-Volterra competition system with diffusion. Mathematical Biosciences and Engineering, 2017, 14(5&6): 1187-1213. doi: 10.3934/mbe.2017061
    [4] M. B. A. Mansour . Computation of traveling wave fronts for a nonlinear diffusion-advection model. Mathematical Biosciences and Engineering, 2009, 6(1): 83-91. doi: 10.3934/mbe.2009.6.83
    [5] Youshan Tao, J. Ignacio Tello . Nonlinear stability of a heterogeneous state in a PDE-ODE model for acid-mediated tumor invasion. Mathematical Biosciences and Engineering, 2016, 13(1): 193-207. doi: 10.3934/mbe.2016.13.193
    [6] Guo Lin, Shuxia Pan, Xiang-Ping Yan . Spreading speeds of epidemic models with nonlocal delays. Mathematical Biosciences and Engineering, 2019, 16(6): 7562-7588. doi: 10.3934/mbe.2019380
    [7] Zhiyin Gao, Sen Liu, Weide Li . Biological control for predation invasion based on pair approximation. Mathematical Biosciences and Engineering, 2022, 19(10): 10252-10274. doi: 10.3934/mbe.2022480
    [8] Peter Hinow, Philip Gerlee, Lisa J. McCawley, Vito Quaranta, Madalina Ciobanu, Shizhen Wang, Jason M. Graham, Bruce P. Ayati, Jonathan Claridge, Kristin R. Swanson, Mary Loveless, Alexander R. A. Anderson . A spatial model of tumor-host interaction: Application of chemotherapy. Mathematical Biosciences and Engineering, 2009, 6(3): 521-546. doi: 10.3934/mbe.2009.6.521
    [9] Ardak Kashkynbayev, Yerlan Amanbek, Bibinur Shupeyeva, Yang Kuang . Existence of traveling wave solutions to data-driven glioblastoma multiforme growth models with density-dependent diffusion. Mathematical Biosciences and Engineering, 2020, 17(6): 7234-7247. doi: 10.3934/mbe.2020371
    [10] Max-Olivier Hongler, Roger Filliger, Olivier Gallay . Local versus nonlocal barycentric interactions in 1D agent dynamics. Mathematical Biosciences and Engineering, 2014, 11(2): 303-315. doi: 10.3934/mbe.2014.11.303
  • 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.



    The systematic analysis of reaction-diffusion models started with the seminal work by Kolmogorov, Petrovskii and Piskunov (see [1]), to predict the kinetic of gases in a model of combustion and by Fisher (see [2]) to understand the mechansims of genes interactions in bioloy. In both cases, the authors proposed a parabolic operator based on the classical order two Laplacian. This operator comes from the application of the Fick law, used to model the global particles motion in a domain. In parallel, the authors established a new conception in the search of solutions which has been of great interest in the last decades. We refer to Travelling Waves solutions. This was a remarkable idea inspired by the fact that any evolution, under the scope of a parabolic operator, shall propagate with a certain velocity (as main driver of the involved dynamics) which may define the shape of the solutions.

    Since its very early formulation, the Fisher-KPP model has attracted the interest of applied scientists, as well as theoreticians. In particular, some other principles in connection with diffusion have been considered in [3,4,5]. Further, the Fisher- KPP model has been analyzed with more complex forms of operators, see for example the heterogeneous diffusion of a higher-order [6], the density and gradient dependant p-Laplacian [7] and the non-local diffusion given by a fractional operator [8].

    Typically, the modelling of diffusion has been focused on the definition of the energy of motion. In the last few decades, the mathematical theories concerning nonlinear operators has developed strongly. This fact has permitted to model diffusion based on the mathematical properties introduced by such operators, along with an identification of the physical observations. In this regard, we can mention the references [9] and [10] where perturbation terms were introduced in the form of nonlinear diffusive operators. Further, in [11] and [12] the oscillatory behavior of solutions was shown to hold for the perturbed Fisher-Kolmogorov equation. Furthermore, other approaches to diffusion can be observed in the literature; in [13] fractional-in-space operators were considered to model anomalous diffusion. The authors started from the standard Brownian motion and they introduced fractional order operator in the sense of the Riesz derivative. In addition, the authors in [14] introduced a new Laplace residual power series algorithm to solved a nonlocal diffusion based in the Caputo derivative and it was applied to a biological problem in connection with bacteria growth systems. Some recent techniques based on homotopy and multigrid method have been considered recently to solve difficult diffusion-convection problems (e.g., [15,16,17]) as well as techniques based on a wavelet multiscale method [18]. In addition, in [19], the authors provided new conceptions based on Bayesian analysis to validate the selection of a model for biological applications.

    In the presented analysis, we focus our attention on the particular interaction between two biological species, namely invasive-invaded. To this end, we shall consider the definition by the Convention of Biological Diversity (p.1, Ch.1 [20]) in what regards with invasive species: "...those alien species which threaten ecosystems, habitats or species".

    The invasive-invaded dynamics has been an active area of mathematical research. In [21], the authors analyzed the haptotactic cell motion, based on an invasive dynamics to model the formation of melanoma. Further, the study in [22] examined the stability of travelling wave solutions based on spectral techniques to predict the cancer formation by haptotaxis processes. In this study, the authors made use of the Evans functions, formulated for a previously linearized diffusion. We should mention that, in the last two references, the diffusion was based on the classical order two Laplacian.

    The process of invasion has been typically modelled with the involvement of only two species. Given an invasive species concentration w that arrives at a new territory previously occupied by another species of concentration v, the following model concretes some of the further general works by Bramson about Fisher-KPP models (for further discussions, the reader is referred to the techniques exposed in [23] by the same author):

    vt=DvΔv+(rvαvvvαvww)vwt=DwΔw+(rwwαwvv)w, (1.1)

    where Dv and Dw are the specific diffusion rates of each species, rv and rw are the specific rates of increase and the α coefficients are defined the intraspecific and interspecific competition balances.

    To introduce a justification of the newly proposed model, we must start from the last model in (1.1) and we considered the following additional principles:

    ● An invasion may occur in a predominant direction. For instance, in the invasion of a cell, any virus can penetrate through a damage in the cell membrane. Similarly, this can be observed in the invasion of a country which may start by the zones where the supplying goods to the main cities are hosted.

    With the intention of modeling the existence of a preferred direction of invasion, we introduce an advection given by the terms cv and cw in the previously mentioned model studied by Bramson. The advection vector is considered to be the same for both species, as there exists a common front where the species fight.

    ● Once the invasive species arrives at a the domain to conquer, the invasion process related with diffusion shall be carefully assessed for finding accurate forms of motion. The selection of an adequate diffusion is supported by the mathematical properties of the associated parabolic operator and its possibilities to describe the observations in a real situation. Generally, we can think that the involved diffusion shall be inherently nonlinear and dependant on the population density. Indeed, once the invasive species arrives at the new domain, there exists a propagating dynamics that may not be fully reproduced by the Gaussian order two diffusion (we recall that this operator induces positive solutions everywhere). This observation makes the order of two diffusions imprecise as the positivity property is not fully compatible with the existence of a propagating front in the invasive species, wherein the diffusivity is further heterogeneous. As a consequence, a study has been conducted to select the most appropriate diffusive operator. Namely and for a general and sufficiently regular function v, we considered the following operators: Porous Medium Equation ((vm1v),m>1), Higher-order ((Δ)pv,p2), Thin film ((|v|nΔv),n>0) and p-Laplacian ((|v|p2v),p>1).

    In addition, we have carried out a review on the existing literature dealing with biological problems formulated with nonlinear diffusion operators. The celebrated model by Keller and Segel (refer to [24] for dedicated insights) was considered as a starting point. Herein, the diffusion is given by a density dependant diffusivity. Another related analyses are provided in [25,26,27,28] where additional dynamics are observed in other scenarios.

    Based on the mentioned ideas, the following problem is proposed:

    vt=DvΔvn+cv+(rvαvvvαvww)v,wt=DwΔwm+cw+(rwwαwvv)w,n,m>0,cRd,d>1,v0(x),w0(x)0L1(Rd)L(Rd), (1.2)

    Note that the second equation in (1.2) includes the reproducing term rw>0 that represents the renewal of forces in the invasive species to make the invasion prosper. It includes, as well, a positive intraspecific term w which includes all forms of generating new invasive agents to conquer the territory. In addition, the interspecific term αwv is the mortality of invasive species because of the action battle with the invaded one. The first equation in (1.2) includes a reproducing term rv in the invaded species. The intraspecific αvv and interspecific αvw coefficients represent the invaded species death rates because of the action of the invasion. The term αvv can be associated to the deaths due to possible disease and hunger conditions and αvw represents the mortality rate in the battle against the invasive species. All the mentioned coefficients shall be correctly rescaled in a particular application so that the carrying capacity in the medium does not impact the postulated equation.

    We should mention that the domain is given in the whole space Rd,d>1, no boundary conditions are considered. This is done to study the problem dynamics openly with no additional restrictions at any border. Our attention is then focused on the characterization of the system interaction and the species balance rather than geometrical restrictions. Precisely speaking, there is a restriction in what regards with the behavior at infinity, as we will discuss the problem with the support of a test function that shall be requested to be null, or small in the tail at |x|.

    The paper layout is as follows: First, the analysis introduces topics related with the regularity and boundedness of solutions together with their existence and uniqueness. The fact of modeling with a nonlinear diffusion (particularly with a density dependant diffusivity) makes us to introduce the analysis supported by a weak formulation. To this end, we will introduce a test function ΥC(Rd). Afterward, the profiles of solutions are explored within a travelling wave formulation. This permits to understand the involved dynamics given in the species during the invasion.

    Definition 1. (Parabolic operator). Consider a function C2+ν,1+ν/2(R3×(0,T)) with pC2+ν(Rd×(0,T)) for any ν>0 and such that p may be either n or m (refer to (1.2)). The following operator L is given by:

    L=tεΔpc. (2.1)

    Note that ε may be Dv or Dw according to the problem (1.2).

    The nonlinear diffusion introduced in (1.2) is density dependant with the associated diffusivity terms

    D(v)=vn1D(w)=wm1, (2.2)

    where n>0 and m>0. According to [29], in case of n, m or both are in the interval (0,1), the problem is driven by a fast diffusion while for any value higher than the unity, the diffusion is referred as slow. This last condition is particularly relevant as, according to [29], the slow diffusion induces the existence of a diffusive front for compactly supported functions that propagates in the domain with finite speed. The front may be regarded as a trace, such that in the inner region the invasive species exists while in the outer region the invasive species is null. This is, the invasive one does not arrive at that territory yet. The existence of such a front means that globally the invasive species reach the territory and expand following such a diffusive front.

    We should mention that both species are considered to be positive or slightly positive in the propagating front and locally in a time interval (0,T] for a T sufficiently small to preserve the parabolicity of the operator L. Furthermore, the physical interpretation of the problem leads us to assume that both solutions (referring to a volumetric or mass concentrations) are below one and positive at the moment of the interaction. Based on these principles, we can make use of the arguments about uniqueness and regularity of solutions provided in [30,31,32,33]. In our case, there exists a remarkable difference compared to such references given by the particular invasive-invaded dynamics considered in the independent terms of (1.2).

    Based on the discussed local parabolicity of the operator L for positive or slightly positive solutions, the following holds:

    Definition 2. (Maximum and minimum solutions). Assume a set of two pairs of solutions: Maximum solutions, (ˆv,ˆw) C2+ν,1+ν/2(Rd×(0,T)), and minimum solutions, (˜v,˜w) C2+ν,1+ν/2(Rd×(0,T)). The set (ˆv,ˆw) of maximum solutions satisfies:

    Lˆvrvˆvαvv˜v2αvw˜w˜v,Lˆwrwˆw2αwv˜v˜w. (2.3)

    The pair formed by (˜v,˜w) is a minimum solution, if the following condition holds:

    L˜vrv˜vαvvˆv2αvwˆwˆv,L˜wrw˜w2αwvˆvˆw. (2.4)

    Consider now:

    ˜v(x,0)v(x,0)ˆv(x,0),˜w(x,0)w(x,0)ˆw(x,0), (2.5)

    where:

    ˆv(x,0)=v(x,0)+δ0,ˆw(x,0)=w(x,0)+δ0,˜v(x,0)=v(x,0)δ0,˜w(x,0)=w(x,0)δ0. (2.6)

    and δ0 shall comply with:

    0<δ0<minxRd{v(x,0),w(x,0)}, (2.7)

    Given the inner parabolicity of the operator L for positive solutions and the continuity of the forcing terms, it is possible to state on the following preliminary regularity condition: (ˆv,ˆw) and (˜v,˜w) L1(Rd×(0,T)) (for additional insights refer to [34] and Ch. 3 in [35]).

    Now, consider a sequence of solutions given by (v(i),w(i))i=0,1,2,... and defined as per the following set of equations:

    v(i)t=DvΔ(v(i))n+cv(i)+(rvαvvv(i1)αvww(i1))v(i1),w(i)t=DwΔ(w(i))m+cw(i)+(rww(i1)αwvv(i1))w(i1)v(0)(x,0)=v0(x),w(0)(x,0)=w0(x),i=1,2,3... (2.8)

    The study of existence and uniqueness of each element of the sequence (v(i),w(i))i=0,1,2,... to (2.8) follows from standard analysis based on the mixed monotone properties of the independent terms and the inner parabolicity of the operator L (refer to the Ch. 3 of [35] and the Theorems 2.1 and 2.2, Ch. 7 in [36]). It should be noted that (2.8) apply to the maximum and minimum solutions in Definition 2 and for any initial condition as stated in (2.6). As a consequence, each element of the sequences satisfies:

    ˜v˜v(i+1)ˆv(i+1)ˆv,˜w˜w(i+1)ˆw(i+1)ˆw. (2.9)

    In addition, following well known techniques, it is straightforward to show that the sequence is ordered (see Ch. 12 of [36]). Then:

    ˆv(i)ˆv(i+1);˜v(i)˜v(i+1);ˆw(i)ˆw(i+1);˜w(i)˜w(i+1). (2.10)

    Hence, the sequence of maximum solutions (ˆvi),ˆw(i)) is nonincreasing while the sequence of minimum solutions (˜v(i),˜w(i)) is nondecreasing.

    Definition 3. (Weak Formulation). Given a test function ΥC(Rd×[0,t]), with 0<ϑ<t<T, the following weak formulation is defined in [ϑ,t]:

    Rdv(t)Υ(t)=Rdv(ϑ)Υ(ϑ)+tϑRd[vΥt+DvvnΔΥcΥv(rvαvvvαvww)vΥ]ds, (2.11)
    Rdw(t)Υ(t)=Rdw(ϑ)Υ(ϑ)+tϑRd[wΥt+DvwmΔΥcΥw(rwwαwvv)wΥ]ds. (2.12)

    Given any spatial ball with radium R, the test function is selected such that any spatial derivative is asymptotically null:

    βΥ(x=x0,t)xβ0+, (2.13)

    for any β=1,2,3... and being x0R.

    Consider the locating point x0R, and the interval Jx=(x0,). The interval Jx is taken for the definition of the following problem in the test function tail, this is in Jx×[0,t]:

    vΥt+DvvnΔΥcΥv(rvαvvvαvww)vΥ=0,wΥt+DvwmΔΥcΥw(rwwαwvv)wΥ=0, (2.14)

    with the following initial conditions:

    v0(x),w0(x),Υ0(x)L(Jx)L1loc(Jx). (2.15)

    Consider now the following truncation for any 0<δ1<1:

    hbδ1={hb,|h|<δ1δb1,|h|δ1},hδ1={h,|h|<δ1δ1,|h|δ1}, (2.16)

    where h may be v or w and b = {n,m}. Based on the defined truncation, the following problem (PΥδ1) shall be considered:

    vδ1Υt+Dvvnδ1ΔΥcΥvδ1(rvαvvvδ1αvwwδ1)vδ1Υ=0,wδ1Υt+Dvwmδ1ΔΥcΥwδ1(rwwδ1αwvvδ1)wδ1Υ=0, (2.17)

    together with (2.15).

    It should be noticed that the functions vδ1 and wδ1 are bounded as per the defined truncation. As a consequence, and based on the initial regularity in Υ, the inner parabolicity of L and the monotone behavior of the forcing terms; we can show the existence and uniqueness of Υ by well known techniques (refer to [36]) during the evolution in Jx=(x0,).

    The coming intention is to show that any weak solution, as stated in the Definition 3, is bounded given the invasive-invaded interaction expressed in the independent terms. This results is particularly interesting because the quadratic growth rate in the wequation leads to blow-up without competition between species. Even in the case of a high increase in the invasive species (due to the arrival of new troops to fight or new viruses to infect), the invaded species can make the invasion to slow down impeding the potential formation of global or local in time blow-up phenomena. To support this, we show that the solutions are sufficiently smooth under the scope of the weak formulation. That smoothness indicated that the problem can be studied asymptotically in the test function tail without risk of losing regularity.

    Theorem 1. Given a pair of weak solutions (v,w), in accordance with (2.11) and (2.12), such a pair is bounded for (x,t)Jx×[ϑ,t].

    Proof. Given an arbitrary positive γR, let us define a cut off function as follows:

    0Ωγ1,ΩγC0(Jx),Ωγ0,|Ωγ|0+and slowly forx,Ωγ=1,x(x0+γ,x1),x1x0,Ωγ=0,x(0,x0). (2.18)

    In addition, it holds that:

    |Ωγ|c1γ,|ΔΩγ|c1γ2. (2.19)

    Now, take the weak formulations as per (2.11) and (2.12) for (x,t)Jx×[ϑ,t] together with a test function verifying Υ(ϑ)=Υ(t)=0. Then, the following problem is given:

    tϑJxvΥtΩγ+tϑJxDvvnΔΥΩγtϑJxcΥvΩγtϑJx(rvαvvvαvww)vΥΩγ=0,tϑJxwΥtΩγ+tϑJxDwwmΔΥΩγtϑJxcΥΩγwtϑJx(rwwαwvv)wΥΩγ=0. (2.20)

    To continue, we shall consider that

    tϑJxubΔΥΩγ=tϑJxubΥΩγ, (2.21)

    where b = {n,m} and u can take the values of v and w, then the following holds:

    tϑJxvΥtΩγtϑJx(rvαvvvαvww)vΥΩγ=tϑJxDvvnΥΩγ+tϑJxcΥvΩγ,tϑJxwΥtΩγtϑJx(rwwαwvv)wΥΩγ=tϑJxDwwmΥΩγ+tϑJxcΥwΩγ. (2.22)

    In order to show the boundedness of v and w in Jx×[ϑ,t], we require the boundedness of the integrals in the right hand side of the last equality. To this end, considering the integration for any spatial variable RJx, it holds that (see [29]):

    tϑubc2(ϑ)R2bb1,tϑuc2(ϑ)R2b1. (2.23)

    Then, the following applies:

    tϑJxubΥΩγJxc2(ϑ)R2bb1|Υ|c1γc2(ϑ)c1JxR2bb11|Υ|. (2.24)

    Assume now the following test function satisfying the conditions: Υ(ϑ)=Υ(t)=0:

    Υ(x,τ)=(1+R2)η(τϑ)(τT), (2.25)

    where ϑτT. It should be noted that η is selected to ensure the convergence for RJx of the integral as per (2.24). Then, the following holds

    c2(ϑ)c1JxR2bb112η(τϑ)(τT)R2η1. (2.26)

    It suffices to consider η=1b1 for a finite integral.

    Based on the same form of a test function, the following assessment holds for the advection term:

    tϑJxcΥΩγuMscJx2ηR2η1u, (2.27)

    where Ms=maxτ[ϑ,T]{Υ(,τ)}. In addition, we shall consider that the solutions represent a mass or volumetric concentration and can take the values 0<v,w<1. The last integral converges for any value η>0.

    Considering now, the assessment done in (2.26), it suffices to take a value for η of the form:

    ηmax{1b1,0}. (2.28)

    In return to the expression (2.22), it applies that:

    tϑJxvΥtΩγ+tϑJx(rv+αvvv+αvww)vΥΩγ(c2(ϑ)c1+c)2ηMs, (2.29)
    tϑJxwΥtΩγ+tϑJx(rww+αwvv)wΥΩγ(c2(ϑ)c1+c)2ηMs. (2.30)

    The right hand terms are formed of finite constants, hence the previous two inequalities allow us to state on the global boundedness of the solutions (v,w) in Jx×[ϑ,t].

    Now, considering the Definition 2, and supported by the shown boundedness of solutions as stated in the last result, we show the convergence to a limit of the sequences of maximum and minimum solutions.

    Theorem 2. Consider the positive initial distributions (v0,w0)>(0,0) for xJx. The sequences of solutions given in Definition 2, (ˆv(i),ˆw(i)) and (˜v(i),˜w(i)) C2+ν,1+ν/2(Jx×(0,T)) converge with a monotone behavior.

    Proof. First, we consider the weak formulation in the Definition 3 for each element of the maximum and minimum sequences. For this purpose, assume the test function given in (2.25) and such that ΥC(Jx) with 0<ϑ<t<T. The following applies for (x,t)Jx×[ϑ,t]:

    Jxˆv(i)(t)Υ(t)=Jx(v0+δ1(i))Υ(0)+tϑJx[ˆv(i)Υt+Dvˆv(i)nΔΥcΥˆv(i)+(rvˆv(i)+αvv˜v(i)2+αvw˜w(i)˜v(i))Υ]ds, (2.31)
    Jx˜v(i)(t)Υ(t)=Jx(v0δ1(i))Υ(0)+tϑJx[˜v(i)Υt+Dv˜v(i)nΔΥcΥ˜v(i)+(rv˜v(i)+αvvˆv(i)2+αvwˆw(i)ˆv(i))Υ]ds, (2.32)

    such that limiδ1(i)=0. Taking the limit and considering the use of the dominated converge theorem, it holds that:

    Jxlimiˆv(i)(t)Υ(t)=Jx(v0+limiδ1(i))Υ(0)+tϑJx[limiˆv(i)Υt+limiDvˆv(i)nΔΥlimicΥˆv(i)+limi(rvˆv(i)+αvv˜v(i)2+αvw˜w(i)˜v(i))Υ]ds. (2.33)

    Now, for ˜v(i):

    Jxlimi˜v(i)(t)Υ(t)=Jx(v0limiδ1(i))Υ(0)+tϑJx[limi˜v(i)Υt+limiDv˜v(i)nΔΥlimicΥ˜v(i)+limi(rv˜v(i)+αvvˆv(i)2+αvwˆw(i)ˆv(i))Υ]ds. (2.34)

    The substraction of the last two equations leads to:

    Jxlimi(ˆv(i)˜v(i))(t)Υ(t)=tϑJx[limi(ˆv(i)˜v(i))Υt+limi(ˆv(i)n˜v(i)n)ΔΥlimi(cΥ(ˆv(i)˜v(i)))+limi(rv(ˆv(i)˜v(i))+αvv(˜v(i)+ˆv(i))(˜v(i)ˆv(i))+αvw(˜w(i)˜v(i)ˆw(i)ˆv(i)))Υ]ds. (2.35)

    Assume now that we depart from sufficiently close maximum and minimum solutions. It suffices to consider this assumption for one of the two solutions, say ˆw(i)˜w(i), then the equality (2.35) holds provided that ˆv(i)˜v(i).

    A similar procedure can be followed for the w-equation, so that we end into a similar expression to that in (2.35), but for the solution w. Then, we shall consider that the maximum and minimum solutions for v are sufficiently close, this is ˆv(i)˜v(i), then it will hold that ˆw(i)˜w(i).

    Further, the limit function preserves the same properties of each element of the converging sequence. As previously pointed, the inner regularity and positivity of the operator L, that is applied to each element of the sequence, allow us to state the following regularity condition for the limit functions:

    v(x,t),w(x,t)C2+ν,1+ν/2(Jx×(ϑ,T))L1((ϑ,T);W2,1(Jx))W1,1((ϑ,T);L1(Jx)). (2.36)

    The coming theorem has the intention of showing the uniqueness of the weak solutions as per the Definition 3.

    Theorem 3. Consider the following pair of minimum and positive solutions to (1.2), (v,w)>(0,0) in Jx×[ϑ,T). Such a pair of solutions coincides locally with the pair of maximum and positive solutions, (ˆv,ˆw)>(0,0), stating then, the local uniqueness of solutions.

    Proof. Consider the following pair of maximum solutions (ˆv,ˆw) in Jx×[ϑ,T) taking on from the initial distribution given by:

    (ˆv(x,0),ˆw(x,0))=(v0(x)+δ1,w0(x)+δ1), (2.37)

    being δ1>0 small. In addition, the pair (ˆv,ˆw) is defined as:

    ˆvt=DvΔˆvn+cˆv+rvˆvαvvv2αvwwv,ˆwt=DwΔˆwm+cˆw+rwˆw2αwvvw. (2.38)

    The pair formed of minimum solutions, and taking on from the initial data (N(x,0),w(x,0))=(v0(x),w0(x)), is given as:

    vt=DvΔvn+cv+rvvαvvˆv2αvwˆwˆv,wt=DwΔwm+cw+rww2αwvˆvˆw. (2.39)

    Assume again a test function satisfying ΥC(Jx) and the associated weak formulations for the last set of equations. Making the difference between the maximum and minimum weak solutions, the following applies:

    0Jx(ˆvv)(t)Υ(t)=tϑJx[(ˆvv)Υt+Dv(ˆvnvn)ΔΥcΥ(ˆvv)+(rv(ˆvv)+αvv(v+ˆv)(vˆv)+αvw(wvˆwˆv))Υ]ds, (2.40)
    0Jx(ˆww)(t)Υ(t)=tϑJx[(ˆww)Υt+Dw(ˆwmwm)ΔΥcΥ(ˆww)+(rw(ˆw+w)(ˆww)+αwv(vwˆvˆw))Υ]ds. (2.41)

    Consider now the definition of the continuous functions given in the following expressions:

    N(x,τ)={ˆvnvnˆvv,ˆvvnvn1,ˆvv},M(x,τ)=[ˆwmwmˆww,ˆwwmwm1,ˆww]. (2.42)

    Given a fixed value for x and for τ=ϑ<T, it holds that:

    0N(x,τ)C1( (2.43)

    Now, consider the following structure for a smooth test function defined in 0 < \vartheta < t < T and satisfying the conditions \Upsilon(\vartheta) = \Upsilon(t) = 0 :

    \begin{equation} \Upsilon(\vert x \vert , t) = e^{J_0(T-t)} (t-\vartheta)\,(t-T) \, (1 + \vert x \vert^2)^{-\kappa_0}, \end{equation} (2.44)

    for finite values of J_0 and \kappa_0 .

    It should be noticed that:

    \begin{equation} \Upsilon_t \leq - J_0 \, M_i \Upsilon(x,t), \, \, \, \, \lvert{\nabla_{\vert x \vert} \Upsilon}\rvert \leq C_3(\kappa_0) \Upsilon(x,t), \, \, \, \, \Delta_{\vert x \vert} \Upsilon\leq C_4(\kappa_0) \Upsilon(x,t), \end{equation} (2.45)

    where M_i = \inf_{t \in (\vartheta, T)} \, \lbrace D_t [(t-\vartheta) (t-T)]\rbrace . For the sake of simplicity in the coming expressions, let us define:

    \begin{equation} F_v (v, w) = -r_v (\hat{v} - v) + \alpha_{vv} (v + \hat{v})(v-\hat{v}) + \alpha_{vw} (wv - \hat{w} \hat{v}), \end{equation} (2.46)
    \begin{equation} F_w (v,w) = -r_w (\hat{w} + w) (\hat{w} - w) + \alpha_{wv} (vw - \hat{v} \hat{w}). \end{equation} (2.47)

    Based on the assessments presented up to now and in return to (2.40) and (2.41), it applies that:

    \begin{equation} \begin{split} & (\hat{v}-v)\,\Upsilon_t+ D_v (\hat{v}^{n}-v^{n})\, \Delta\Upsilon- c \cdot \nabla \Upsilon(\hat{v} - v)+ F_v(u,w) \, \Upsilon\\& \leq [-J_0 M_i + D_u n C_4 + c C_3 ] \Upsilon(\tilde{v} - v)+ F_v (v,w) \, \Upsilon, \end{split} \end{equation} (2.48)
    \begin{equation} \begin{split} & (\hat{w}-w)\,\Upsilon_t+ D_w (\hat{w}^{m}- w^{m})\, \Delta\Upsilon - c \cdot \nabla \Upsilon(\hat{w} - w)+ F_w(v,w) \, \Upsilon\\& \leq [-J_0 M_i + D_v m C_4 + c C_3 ] \Upsilon(\tilde{w} - w) +F_w(v,w) \, \Upsilon. \end{split} \end{equation} (2.49)

    The finite J_0 is selected to comply with the following inequality:

    \begin{equation} -J_0 M_i + D_v n C_4 + c C_3 \leq 0, \, \, \, \, -J_0 M_i + D_w m C_4 + c C_3 \leq 0. \end{equation} (2.50)

    Then, it suffices to take:

    \begin{equation} J_0 M_i \geq \max \lbrace D_v n C_4 + c C_3 , D_w m C_4 + c C_3 \rbrace, \end{equation} (2.51)

    so that:

    \begin{equation} [-J_0 M_i + D_v n C_4 + c C_3 ] \Upsilon(\tilde{v} - v)+F_v(v,w) \, \Upsilon \leq F_v(v,w) \, \Upsilon \end{equation} (2.52)

    and

    \begin{equation} [-J_0 M_i + D_w m C_4 + c C_3 ] \Upsilon(\tilde{w} - w) + F_w (v,w) \, \Upsilon \leq F_w (v,w) \, \Upsilon. \end{equation} (2.53)

    The expressions (2.40) and (2.41) are rewritten as follows:

    \begin{equation} 0 \leq \int_{J_x} (\hat{v}-v)(t)\,\Upsilon(t) \leq \int_{\vartheta}^{t} \int_{J_x} F_v (v,w) \, \Upsilon\, ds, \end{equation} (2.54)
    \begin{equation} 0 \leq \int_{J_x} (\hat{w}-w)(t)\,\Upsilon(t) \leq \int_{\vartheta}^{t} \int_{J_x} F_w (v,w) \, \Upsilon\, ds. \end{equation} (2.55)

    The uniqueness of the solutions is proved locally in time. For this purpose, assume that t \in (\vartheta, \varepsilon_0) for 0 < \varepsilon_0 \ll 1 . Furthermore, consider that in the mentioned t-interval, we can introduce, by simple translation, a \varepsilon_0 -contact between the minimum and maximum solutions for v , this is 0 < (\hat{v}-v)(\vartheta) \leq \varepsilon_0 . Now, let us make the first time derivative in (2.55), so that the following applies:

    \begin{equation} 0 \leq \frac{d}{dt} \int_{J_x} (\hat{w}-w)\,\Upsilon\leq \int_{J_x} (r_w \alpha_{w} 2 M_w + \varepsilon_0 \alpha_{wv}) ( \hat{w} - w ) \, \Upsilon\, , \end{equation} (2.56)

    where M_w = \max_{t \in (\vartheta, T)} \, \lbrace \hat{w} \rbrace . Note that the difference between the maximal and the minimal solution is certainly positive, but the difference between them may be decreasing over time to ensure the uniqueness of solutions, leading to \frac{d(\hat{w} - w)}{dt} < 0 . Assume that A_w = r_w \alpha_{w} 2 M_w + \varepsilon_0 \alpha_{wv} , then, and making a standard resolution, the following holds:

    \begin{equation} \int_{J_x} (\hat{w}-w)\,\Upsilon\leq (\hat{w}_0 - w_0) \, e^{-A_w \, t}, \end{equation} (2.57)

    where the minus sign in the exponential term is introduced to ensure that \frac{d(\hat{w} - w)}{dt} < 0 . Consider now the expression given in (2.37), \hat{w}(x, 0) - w(x, 0) = \delta_1 , then:

    \begin{equation} \int_{J_x} (\hat{w}-w)\,\Upsilon\leq \delta_1 \, e^{-A_w \, t}. \end{equation} (2.58)

    The local uniqueness is shown for \delta_1 \rightarrow 0^+ , then it holds that:

    \begin{equation} 0 \leq \int_{J_x} (\hat{w}-w)\,\Upsilon\leq 0^+. \end{equation} (2.59)

    The only possible condition in the last expression is given by \hat{w} = w , showing then the local uniqueness of the solutions. Now, return to (2.54) and define the constant A_v = r_v + 2 \alpha_{vv} M_v + \alpha_{vw} M_w where M_v = \max_{t \in (\vartheta, T)} \, \lbrace \hat{v} \rbrace . It holds that:

    \begin{equation} 0 \leq \int_{J_x} (\hat{v}-v)\,\Upsilon\leq \int_{\vartheta}^{t} \int_{J_x} A_v (\hat{v}-v) \Upsilon\, ds. \end{equation} (2.60)

    Performing the first time derivative and solving the resulting inequality by standard techniques, the following holds:

    \begin{equation} \int_{J_x} (\hat{v}-v)\,\Upsilon\leq \delta_1 \, e^{-A_v \, t}. \end{equation} (2.61)

    Making again \delta_1 \rightarrow 0^+ :

    \begin{equation} 0 \leq \int_{J_x} (\hat{v}-v)\,\Upsilon\leq 0^+, \end{equation} (2.62)

    where the only compatible result is to consider \hat{v} = v , locally in time.

    The analysis provided allows us to conclude on the local in time uniqueness of solutions as initially postulated.

    The intention now is to provide analytical solutions to describe the propagation of the invasive and invaded species in the domain while they are fighting. As mentioned, the nonlinear diffusion induces the existence of a support in the solutions that propagates with finite speed. The determination of a precise analytical expression to describe the evolution of the front will be of help to further understand the involved dynamics between the species.

    It should be noticed that the evolution of the support requires to study the system (1.2) for small values in v and w leading to the following alternative problem:

    \begin{equation} \begin{array}{c} v_t = D_v \Delta v^n + c \cdot \nabla v + r_v v , \\ \\ w_t = D_w \Delta w^m + c \cdot \nabla w + \varepsilon_3 , \\ \\ v_0(x), w_0(x) \geq 0 \in L^1( \mathbb{R}^d) \cap L^{\infty}( \mathbb{R}^d), \end{array} \end{equation} (3.1)

    where \vert \varepsilon_3 \vert \ll 1 and compiles the order two terms in the equation for w .

    The coming theorem aims at providing the spatio-temporal evolution of the support. As it will be shown, there are several propagating modes for each species.

    Theorem 4. As a result of the interaction between species and for \vert v \vert \ll 1, \, \, \vert w \vert \ll 1 , there exist two propagating modes for each species in the (x, t) variables, that have the following structures:

    For the concentration v :

    \begin{equation} \lvert{x}\rvert_+ = \vert c \vert t^2 + t \sqrt{ \frac{2 n}{1-n} \left( A_v - \frac{(v_f)^{n-1}}{t} \right)}, \end{equation} (3.2)
    \begin{equation} \lvert{x}\rvert_- = \vert c \vert t^2 - t \sqrt{ \frac{2 n}{1-n} \left( A_v - \frac{(v_f)^{n-1}}{t} \right)}. \end{equation} (3.3)

    For the concentration w :

    \begin{equation} \lvert{x}\rvert_+ = \vert c \vert t^2 + t \sqrt{ \frac{2 m}{1-m} \left( A_w - \frac{(w_f)^{m-1}}{t} \right)}, \end{equation} (3.4)
    \begin{equation} \lvert{x}\rvert_- = \vert c \vert t^2 - t \sqrt{ \frac{2 m}{1-m} \left( A_w - \frac{(w_f)^{m-1}}{t} \right)}, \end{equation} (3.5)

    where

    \begin{equation} A_v (\lVert{v_0}\rVert_\infty) > 0 \, , \, \, A_w (\lVert{w_0}\rVert_\infty) > 0, \end{equation} (3.6)

    and v_f , w_f refer to the levels of species concentration in the propagating front. This is mainly related with the level of individuals of each species effectively fighting.

    Proof. First, consider the equation:

    \begin{equation} v_t = D_v \Delta v^n + c \cdot \nabla v + r_v v. \end{equation} (3.7)

    The determination of an analytical expression for the support requires to solve the above equation, and this is done with the help of a selfsimilar solution. To this end, assume that the advection coefficient is small, \vert c \vert \ll 1 . This is a sufficiently accurate approximation in the proximity of the support, where the dynamics is focused on the interaction between species rather than in the transport motility induced by the advection. Nonetheless, the advection term will be kept in the equation during the resolution process as it will have an influence in the diffusive front, but we shall keep in mind the requirement of \vert c \vert \ll 1 . Then:

    \begin{equation} v(x,t) = t^{-\theta_1} F(\xi) \, \, , \, \, \, \, \, \xi = \lvert{x}\rvert t^{\gamma_1}. \end{equation} (3.8)

    Hence, in (3.7)

    \begin{equation} -\theta_1 t^{-\theta_1-1}F + \gamma_1 \xi t^{-\theta_1-1} F^{\prime} = M(n, D_v, F^{\prime \prime}, F^\prime, F) t^{-\theta_1 n + 2 \gamma_1} + c t^{-\theta_1 + \gamma_1} F^\prime + r_v t^{-\theta_1} F(\xi) . \end{equation} (3.9)

    Balancing the leading terms in the previous equations and considering the exponents in t , the following holds:

    \begin{equation} -\theta_1 - 1 = -\theta_1 n + 2 \gamma_1, \, \, \, \, \, \, \, \, -\theta_1 n + 2 \gamma_1 = -\theta_1 + \gamma_1. \end{equation} (3.10)

    The following values for \theta_1 and \gamma_1 are obtained:

    \begin{equation} \theta_1 = \frac{1}{1-n}, \, \, \, \, \, \gamma_1 = -1. \end{equation} (3.11)

    Consider again the leading terms in (3.9), then the selfsimilar profile F satisfies the following elliptic equation (where r_v = 1 for simplicity):

    \begin{equation} (\gamma_1 \xi - c) F^{\prime} = (F^{n})'' + \frac{2}{\xi} (F^{n})' + F. \end{equation} (3.12)

    The family of solutions to the postulated elliptic equation is (see [35] for additional insights):

    \begin{equation} F(\xi) = (A_v-B_v \xi^2)^ {\frac{1}{n-1}}. \end{equation} (3.13)

    Replacing this last expression in (3.12), and after performing standard assessments, the following values for A_v and B_v are obtained:

    \begin{equation} A_v (\Vert v_0 \Vert) > 0, \, \, \, \, \, \, B_v = \frac{n-1}{2n \gamma_1}. \end{equation} (3.14)

    Now, we make similar operations for the w -equation in (3.1) and considering a small \varepsilon_3 in comparison with the other involved terms (diffusion and advection):

    \begin{equation} \Theta(x,t) = t^{-\theta_2} G(\xi) \, \, , \, \, \, \, \, \xi = \lvert{x}\rvert t^{\gamma_2}, \, \, \, \, \, \theta_2 = \frac{1}{1-m}, \, \, \, \, \, \gamma_2 = -1, \end{equation} (3.15)

    where:

    \begin{equation} G(\xi) = (A_w-B_w \xi^2)^ {\frac{1}{m-1}}, \, \, \, \, \, \, \, \, A_w (\Vert \Theta_0 \Vert) > 0, \, \, \, \, \, \, B_w = \frac{m-1}{2m\gamma_2}. \end{equation} (3.16)

    The intention is to introduce an analytical expression for the propagating front. For this, we fix two concentration levels for each of the species involved. Let us consider that such fixed values are given by w = w_f (concentration level of the invasive species in the propagating front) and v = v_f (concentration level in the front for the invaded species). Coming again to the selfsimilar expression:

    \begin{equation} v_f = t^{-\theta_1} F(\xi) \, \rightarrow \, F(\xi) = t^{\frac{1}{1-n}} \, v_f, \end{equation} (3.17)

    where F(\xi) is as per (3.13) and \xi = \lvert{x}\rvert t^{-1} . After making the necessary standard assessments, the following equations describe the evolution of the support:

    \begin{equation} \lvert{x}\rvert_+ = \vert c \vert t^2 + t \sqrt{ \frac{2 n}{1-n} \left( A_v - \frac{(v_f)^{n-1}}{t} \right)}, \end{equation} (3.18)

    and

    \begin{equation} \lvert{x}\rvert_- = \vert c \vert t^2 - t \sqrt{ \frac{2 n}{1-n} \left( A_v - \frac{(v_f)^{n-1}}{t} \right)}, \end{equation} (3.19)

    And similarly for the equation in w :

    \begin{equation} \lvert{x}\rvert_+ = \vert c \vert t^2 + t \sqrt{ \frac{2 m}{1-m} \left( A_w - \frac{(w_f)^{m-1}}{t} \right)}, \end{equation} (3.20)
    \begin{equation} \lvert{x}\rvert_- = \vert c \vert t^2 - t \sqrt{ \frac{2 m}{1-m} \left( A_w - \frac{(w_f)^{m-1}}{t} \right)}. \end{equation} (3.21)

    It should be noticed that the constants involved in the equations for the front (n, m, c, A_v, A_w, v_f, w_f) shall be obtained supported by experimental observations about the particular dynamics. The constants A_v and A_w are related with the boundedness condition in the initial distributions and can simply take the values of the maximum level of initial invaded and invasive species. Some specific values for the levels v_f and w_f can be obtained from the numerical simulations (see Section 4) if we consider the interacting front where both species are fighting.

    The profiles of travelling waves are expressed as v(x, t) = f(\omega_1), \, \omega_1 = x \cdot n_d - \lambda_1 t and w(x, t) = g(\omega_2), \, \omega_2 = x \cdot n_d - \lambda_2 t \, , where n_d \in \mathbb{R}^d is a unitary vector that provides the travelling wave propagation direction and \lambda_1, \lambda_2 are defined the propagation speed for each travelling profile that comply with: f, g: \mathbb{R} \rightarrow (0, \infty) and are bounded as it will be shown afterward in the graphical representations. It should be noted that two travelling wave profiles are equivalent under a translation \omega \rightarrow \omega + \omega_0 and a symmetry \omega \rightarrow -\omega . For the sake of convenience, assume that the vector n_d = (1, 0, ..., 0) , then v(x, t) = f(\omega_1), \, \, \, \, \omega_1 = x - \lambda_1 t \, \in \mathbb{R} and w(x, t) = g(\omega_2), \, \, \, \, \omega_2 = x - \lambda_2 t \, \in \mathbb{R} .

    Consider now v(x, t) = f(\omega_1) and w(x, t) = g(\omega_2) , then the problem (1.2) is expressed under the travelling wave profiles formulation:

    \begin{equation} \begin{array}{c} -(\lambda_1 + c) f^\prime - D_v n(n-1)f^{n-1} (f^\prime)^2 + (\alpha_{vv} f + \alpha_{vw} g - r_v) f = n f^{n-1} f^{\prime \prime}, \\ \\ -(\lambda_2 + c) g^\prime - D_w m(m-1)g^{m-1} (g^\prime)^2 + (\alpha_{wv} f - r_w g) g = m g^{m-1} g^{\prime \prime}, \\ \\ f, g \in L^{\infty}( \mathbb{R}), \, \, \, \, f(-\infty) = 1, \, \, \, \, \, g(\infty) = 1. \end{array} \end{equation} (4.1)

    It is our aim to introduce a numerical analysis to depict the involved dynamics and to represent solutions to the problem (4.1). All the constants given in the model (1.2) are assumed to have particular values:

    \begin{equation} c = 1; \, n = 2; \,m = 2; \, r_v = 2; \, \alpha_{vv} = 1 \, \alpha_{vw} = 1; \, r_w = 1; \, \alpha_{wv} = 1. \end{equation} (4.2)

    We should mention that the solutions have been tested for other combinations in the given constants. In this regard, the behavior of the solutions follow similar profiles to that for the values in (4.2). This exercise was done during the ordinary trials in the numerical procedure. We do not reproduce such information here to account for a concise representation. Hence, the parameters, to be modified in the numerical assessments, are the travelling wave speeds in each of the profile, this is \lambda_1 and \lambda_2 .

    The numerical simulation is performed with the Matlab function bvp4c. For further details about the intrinsic numerical assessments executed by the function bvp4c, the reader can refer to [37].

    The domain size of integration has been assumed as R = 1000 (although the presented figures will be provided with a zoom in the interacting area) to minimize the impact of the collocation method at the borders. The number of nodes considered was 100,000 with an absolute error of 10^{-5} .

    The results are compiled in the Figures 16 for different combinations in the travelling wave speeds. It shall be noted that the invasive species comes from the right and interact with the invaded species. The interacting front moves from the right to the left (recall that the travelling wave profiles are invariant under a translation) so that the interacting front occupies the whole region until the invaded species cannot support the invasion. In addition, we should mention that the introduction of the nonlinear diffusion may induce oscillations that lead solutions to be outside of the interval [0, 1] . In a real application, it is not possible to get values outside of the mentioned interval. Then, for negative values in one of the species, we shall consider that the species concentration vanishes and for values higher than one, the species concentrations reach the maximum allowed in the region.

    Figure 1.  Representation of travelling waves solutions. It is interesting to observe the oscillations in the invaded species. Certainly this is related with the pressure density diffusive term which makes the invaded species to move forward and to resist the invasion. The interacting zone moves from the right to the left until the invasive species conquers the whole territory and leaves the invaded species with low capacity to support the fight. Note the less steep zone in the invasive species which is related with a rapid movement in the terrain until the confrontation with the invaded species occurs in the front.
    Figure 2.  Representation of travelling waves solutions. In this occasion, we observe certain instabilities in the system. There exists a sudden drop in the invasive species before reaching the interacting front. In this case, the dynamics in the front is quicker compared to the Figure 1 and the invaded species is more sensitive to the presence of the invasive, that introduces the conquering forces faster into the domain.
    Figure 3.  Representation of travelling waves solutions. Again, increasing the travelling wave speeds leads to a fast dynamics in the interactive front (the front is confined in the left margin of the picture).
    Figure 4.  Representation of travelling waves solutions for different values in the travelling wave speeds. Note that \omega_i refers to either w_1 or w_2 . It is possible to check that the invasive species dynamics is less steep compared with the previous figures.
    Figure 5.  Representation of travelling waves solutions for different values in the travelling wave speeds. Note that \omega_i refers to either w_1 or w_2 . It is important to remark that for different wave speeds the dynamics gets unstable leading to the formation of oscillatory profiles.
    Figure 6.  Representation of travelling wave solutions for significantly different values in the travelling wave speeds. Note that \omega_i refers to either w_1 or w_2 . It should be noted that the dynamics gets further unstable compared to that in the Figure 5. The fast dynamics in the invaded species makes the minimum invasive to react. Even in the region \omega \in [-4.82, -3.01] , the invaded species gets above the invasive one because of its quick dynamics, lately the invasive species enters into the front located on the left to continue the invasion.

    The initial objectives outlined have been accomplished based on the analytical and numerical assessments provided. The new set of equations in (1.2) were shown to comply with the conditions of regular and unique weak solutions. We obtained exact analytical expressions in the (x, t) variables for the propagating fronts, that modeled the interaction of the species while they were fighting. Eventually, the problem was analyzed within the scope of the travelling waves formulation, and a numerical assessment provided to further understand the interaction between the species. We observed here that a significant difference in the travelling waves speeds, between the involved species, provided an unstable behavior in the solutions.

    This work has been financed by the CEU San Pablo University, grant number MGI22AOQ.

    The author 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.
  • 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(1099) PDF downloads(41) Cited by(11)

Figures and Tables

Figures(16)  /  Tables(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog