Research article

Boundedness of fractional integrals on grand weighted Herz spaces with variable exponent

  • Received: 21 August 2022 Revised: 22 September 2022 Accepted: 22 September 2022 Published: 12 October 2022
  • MSC : Primary 46E30, Secondary 47B38

  • In this paper, we introduce grand weighted Herz spaces with variable exponent and prove the boundedness of fractional integrals on these spaces.

    Citation: Babar Sultan, Mehvish Sultan, Mazhar Mehmood, Fatima Azmi, Maryam Ali Alghafli, Nabil Mlaiki. Boundedness of fractional integrals on grand weighted Herz spaces with variable exponent[J]. AIMS Mathematics, 2023, 8(1): 752-764. doi: 10.3934/math.2023036

    Related Papers:

    [1] Robert Carlson . Myopic models of population dynamics on infinite networks. Networks and Heterogeneous Media, 2014, 9(3): 477-499. doi: 10.3934/nhm.2014.9.477
    [2] Shi Jin, Min Tang, Houde Han . A uniformly second order numerical method for the one-dimensional discrete-ordinate transport equation and its diffusion limit with interface. Networks and Heterogeneous Media, 2009, 4(1): 35-65. doi: 10.3934/nhm.2009.4.35
    [3] Sabrina Bonandin, Mattia Zanella . Effects of heterogeneous opinion interactions in many-agent systems for epidemic dynamics. Networks and Heterogeneous Media, 2024, 19(1): 235-261. doi: 10.3934/nhm.2024011
    [4] Nastassia Pouradier Duteil . Mean-field limit of collective dynamics with time-varying weights. Networks and Heterogeneous Media, 2022, 17(2): 129-161. doi: 10.3934/nhm.2022001
    [5] Xu Yang, François Golse, Zhongyi Huang, Shi Jin . Numerical study of a domain decomposition method for a two-scale linear transport equation. Networks and Heterogeneous Media, 2006, 1(1): 143-166. doi: 10.3934/nhm.2006.1.143
    [6] Christian Budde, Marjeta Kramar Fijavž . Bi-Continuous semigroups for flows on infinite networks. Networks and Heterogeneous Media, 2021, 16(4): 553-567. doi: 10.3934/nhm.2021017
    [7] Rong Huang, Zhifeng Weng . A numerical method based on barycentric interpolation collocation for nonlinear convection-diffusion optimal control problems. Networks and Heterogeneous Media, 2023, 18(2): 562-580. doi: 10.3934/nhm.2023024
    [8] Raul Borsche, Axel Klar, T. N. Ha Pham . Nonlinear flux-limited models for chemotaxis on networks. Networks and Heterogeneous Media, 2017, 12(3): 381-401. doi: 10.3934/nhm.2017017
    [9] Paulo Amorim, Alessandro Margheri, Carlota Rebelo . Modeling disease awareness and variable susceptibility with a structured epidemic model. Networks and Heterogeneous Media, 2024, 19(1): 262-290. doi: 10.3934/nhm.2024012
    [10] Elisabeth Logak, Isabelle Passat . An epidemic model with nonlocal diffusion on networks. Networks and Heterogeneous Media, 2016, 11(4): 693-719. doi: 10.3934/nhm.2016014
  • In this paper, we introduce grand weighted Herz spaces with variable exponent and prove the boundedness of fractional integrals on these spaces.



    Mathematical modeling in epidemiology has certainly experienced an impressive increase in recent times driven by the overwhelming effects of the COVID-19 pandemic. A large part of the research has been directed towards the construction of models capable of describing specific characteristics associated with the pandemic, in particular the presence of asymptomatic individuals, which were largely underestimated, especially in the early stages of the spread of the disease [14,23,38,36,41]. Another line of research was aimed at the construction of models capable of describing the spatial characteristics of the epidemic, in order to be able to properly assess the impact of containment measures, in particular with regard to the mobility on the territory and restrictions in areas with higher risk of infection [13,43,42,16,1,39]. Other approaches took into account the network structure of connections [10], social heterogeneity aspects [21], and the multiscale nature of the pandemic [6]. See also [37,32] and the survey [2] for some recent developments in epidemiological modelling based on the use of kinetic equations.

    Regardless of the characteristics of the model, a common aspect concerns the uncertainty of the data, which, especially in the first phase of the pandemic, largely underestimated the number of infected individuals. More precisely, the difficulty in correctly identifying all infected individuals in the early stages of the pandemic, due to structural limitations in performing large-scale screening and the inability to track the number of contacts, required the introduction of stochastic parameters into the models and the construction of related techniques for quantifying uncertainty [3,4,9,7].

    Among the various techniques of uncertainty quantification, the approaches based on stochastic strategies that do not necessarily require a-priori knowledge of the probability distribution of the uncertain parameters, as needed in the case of methods based on generalized polynomial chaos [46], are particularly interesting in view of a comparison with experimental data. On the other hand, the low convergence rate of Monte Carlo type sampling techniques poses serious limitations to their practical use.

    In this context, multi-fidelity methods have shown to be able to efficiently alleviate such limitations through control variate techniques based on an appropriate use of low-fidelity surrogate models able to accelerate the convergence of stochastic sampling [48,35,30,18,19]. Specifically, general multi-fidelity approaches for kinetic equations have been developed in [18,19], while bi-fidelity techniques with greedy sample selection in [30,31]. We refer also to the recent survey in [17].

    The present paper is devoted to extend the bi-fidelity method for transport equations in the diffusive limit developed in [30] to the case of compartmental systems of multiscale equations designed to model mobility dynamics in an epidemic setting with uncertainty [13,7]. For this purpose, the corresponding hyperbolic system recently introduced in [10,9] will be used as a reduced low-fidelity model. The two models allow to correctly describe the hyperbolic dynamics of the movement of individuals over long distances together with the small-scale, high-density, diffusive nature typical of urban areas [13]. In addition, both models share the same diffusive limit in which it is possible to recover classical models of diffusive type [26]. From a numerical viewpoint, a space-time asymptotic-preserving discretization, which work uniformly in all regimes, has been adopted in combination with the bi-fidelity approach. This permits to obtain efficient stochastic asymptotic-preserving methods.

    The rest of the manuscript is organized as follows. In Section 2 we introduce the epidemic transport model with uncertainty in the simplified SIR compartmental setting, together with its diffusive limit. The corresponding reduced order two-velocity model used as low-fidelity surrogate is also discussed. Next, we extend the previous modeling to more realistic compartmental settings based on the introduction of the exposed and asymptomatic compartments in Section 3. The details of the bi-fidelity method and the asymptotic-preserving IMEX Finite Volume scheme are then given in Section 4. Section 5 contains several numerical experiments that illustrate the performance of the bi-fidelity approach. Some conclusions are reported at the end of the manuscript.

    For simplicity, we first illustrate the modeling in the case of a classic SIR compartmental dynamic and subsequently we will extend our arguments to a more realistic SEIAR model, designed to take into account specific features of the COVID-19 pandemic, in Section 3.

    Let us consider a random vector z=(z1,,zd)TRd characterizing possible sources of uncertainty due to the independent stochastic parameters z1,,zd, which may affect variables of the mathematical model, as well as parameters or initial conditions. Individuals at position xΩR at time t moving with velocity v[1,1] are denote by fS=fS(x,v,t,z), fI=fI(x,v,t,z) and fR=fR(x,v,t,z), which are the respective kinetic densities of susceptible S (individuals who may be infected by the disease), infectious I (individuals who may transmit the disease) and removed R (individuals healed or died due to the disease). We assume to have a population with subjects having no prior immunity and neglect the vital dynamics represented by births and deaths due to the time scale considered.

    The kinetic distribution is then given by

    f(x,v,t,z)=fS(x,v,t,z)+fI(x,v,t,z)+fR(x,v,t,z),

    and we recover their total density by integration over the velocity space

    ρ(x,t,z)=11f(x,v,t,z)dv.

    As a consequence,

    S(x,t,z)=11fS(x,v,t,z)dvI(x,t,z)=11fI(x,v,t,z)dvR(x,t,z)=11fR(x,v,t,z)dv, (1)

    with ρ(x,t,z)=S(x,t,z)+I(x,t,z)+R(x,t,z), denote the density fractions of the population at position x and time t>0 that are susceptible, infected and recovered respectively. In this setting, the kinetic densities satisfy the epidemic transport equations [13]

    fSt+vSfSx=F(fS,I)+1τS(S2fS)fIt+vIfIx=F(fS,I)γfI+1τI(I2fI)fRt+vRfRx=γfI+1τR(R2fR), (2)

    where vS=λS(x)v, vI=λI(x)v, vR=λR(x)v, λS,λI,λR0 take into account the heterogeneities of geographical areas, and are thus chosen dependent on the spatial location. Similarly, also the relaxation times τS=τS(x), τI=τI(x) and τR=τR(x). The quantity γ=γ(x,z) is the recovery rate of infected, while the transmission of the infection is governed by an incidence function F(,I) modeling the transmission of the disease [25]

    F(g,I)=βgIp1+κI, (3)

    with the classic bi-linear case corresponding to p=1, κ=0, even though it has been observed that an incidence rate that increases more than linearly with respect to the number of infected I can occur under certain circumstances [15,5,29]. The parameter β=β(x,t,z) characterizes the average number of contacts per person per time, multiplied by the probability of disease transmission in a contact between a susceptible and an infectious subject; whereas κ=κ(x,t,z)>0 acts as an incidence damping parameter based on the self-protective behavior of the individual that arises from awareness of the epidemic risk as the disease progresses [44,22]. Notice that both β and κ may vary in time as a consequence of governmental control actions, such as mandatory wearing of masks or full lockdowns, and the increasing awareness of the epidemic risks among the population.

    The standard threshold of epidemic models is the well-known reproduction number R0, which defines the average number of secondary infections produced when one infected individual is introduced into a host population in which everyone is susceptible [25]. This number determines when an infection can invade and persist in a new host population. For many deterministic infectious disease models, an infection begins in a fully susceptible population if and only if R0>1. Assuming no inflow/outflow boundary conditions in Ω, integrating over velocity/space and summing up the second equation in (2) we have

    tΩI(x,t,z)dx=ΩF(S,I)dxΩγ(x,z)I(x,t,z)dx0

    when

    R0(t,z)=ΩF(S,I)dxΩγ(x,z)I(x,t,z)dx1. (4)

    The above quantity, therefore, defines the stochastic reproduction number for system (2) describing the space averaged instantaneous variation of the number of infectious individuals at time t>0. This definition naturally extends locally by integrating over any subset of the computational domain Ω if one ignores the boundary flows.

    Let us introduce the flux functions

    JS=λS11vfS(x,v,t,z)dv,JI=λI11vfI(x,v,t,z)dv,JR=λR11vfR(x,v,t,z)dv. (5)

    Then, integrating the system (2) against v, it is straightforward to get the following set of equations for the macroscopic densities of commuters

    St+JSx=F(S,I)It+JIx=F(S,I)γIRt+JRx=γI (6)

    whereas the flux functions satisfy

    JSt+λS211v2fSxdv=F(JS,I)JSτSJIt+λI211v2fIxdv=λIλSF(JS,I)γJIJIτIJRt+λR211v2fRxdv=λRλIγJIJRτR. (7)

    By introducing the space dependent diffusion coefficients

    DS=13λ2SτS,DI=13λ2IτI,DR=13λ2RτR, (8)

    which characterize the diffusive transport mechanism of susceptible, infectious and removed, respectively, and keeping the above quantities fixed while letting the relaxation times τS,I,R to zero, we get from the r.h.s. in (2)

    fS=S/2,fI=I/2,fR=R/2,

    and consequently, from (7), we obtain a proportionality relation between the fluxes and the spatial derivatives (Fick's law):

    JS=DSSx,JI=DIIx,JR=DRRx. (9)

    Thus, substituting (9) into (6) we get the diffusion system [34,40,45,26]

    St=F(S,I)+x(DSSx)It=F(S,I)γI+x(DIIx)Rt=γI+x(DRRx). (10)

    We remark that the model's capability to account for different regimes, hyperbolic or parabolic, accordingly to the space dependent values τS, τI, τR, makes it suitable for describing the dynamics of populations composed of human beings. Indeed, it is clear that the daily routine is a complex mixing of individuals moving at the scale of a city and individuals moving among different urban centers. In this situation, it seems reasonable to avoid, due to the lack of microscopic information and the high complexity, the description of the details of movements within an urban area and to describe this aspect through a diffusion operator. On the other hand, individuals when moving from one city to another follow well established connections for which a hyperbolic setting is certainly more appropriate.

    The low-fidelity model, is based on considering individuals moving in two opposite directions (indicated by signs "+" and "-"), with velocities ±λS for susceptible, ±λI for infectious and ±λR for removed, we can describe the spatio-temporal dynamics of the population through the following two-velocity epidemic model [10]

    S±t+λSS±x=F(S±,I)12τS(S+S)I±t+λII±x=F(S±,I)γI±12τI(I+I)R±t+λRR±x=γI±12τR(R+R). (11)

    In the above system, individuals S(x,t,z), I(x,t,z) and R(x,t,z) are defined as

    S=S++S,I=I++I,R=R++R.

    The transmission of the infection is governed by the same incidence function as in the high-fidelity model, defined in (3). Also the definition of the basic reproduction number R0 results the same previously introduced in (4).

    If we now introduce the fluxes, defined by

    JS=λS(S+S),JI=λI(I+I),JR=λR(R+R), (12)

    we obtain a hyperbolic model equivalent to (11), but presenting a macroscopic description of the propagation of the epidemic at finite speeds, for which the densities follow system (6) and equations of fluxes read

    JSt+λ2SSx=F(JS,I)JSτSJIt+λ2IIx=λIλSF(JS,I)γJIJIτIJRt+λ2RRx=λRλIγJIJRτR. (13)

    Let us now consider the behavior of the low-fidelity model in diffusive regimes. To this aim, we introduce the diffusion coefficients

    DS=λ2SτS,DI=λ2IτI,DR=λ2RτR. (14)

    As for the previous model, the diffusion limit of the system is formally recovered letting the relaxation times τS,I,R0, while keeping the diffusion coefficients (14) finite. Under this scaling, from (13) we get equations (9), which inserted into (6) lead again to the parabolic reaction-diffusion system (10).

    Remark 1. An important aspect in the bi-fidelity approach here proposed, is that the high-fidelity model and the low-fidelity one exactly coincide in the diffusive limit. The only difference lays in the definition of the diffusion coefficients (see eqs. (8) and eqs. (14)). As a consequence, in such a regime the bi-fidelity method achieves the maximum accuracy.

    To account for more realistic models to analyze the evolution of the ongoing COVID-19 pandemic, we consider extending the simple SIR compartmentalization by taking into account two additional population compartments, E and A, resulting in a SEIAR model [9,7]. Subjects in the E compartment are the exposed, hence infected but not yet infectious, being in the latent period. Moreover, among the infectious subjects, we distinguish the population between a group of individuals I who will develop severe symptoms and a group of individuals A who will never develop symptoms or, if they do, these will be very mild. This feature turns out to be essential to correctly analyze the evolution of COVID-19. In fact, it has been shown that individuals belonging to the A group are very difficult to detect and isolate, contributing more strongly to the spread of the virus than the more easily detectable I individuals [23,36,41].

    Let us now consider a population at position xΩ moving with velocity directions v[1,1], still taking into account possible uncertainties related to the random vector z. Defining the kinetic densities of susceptible fS=fS(x,v,t,z), exposed fE=fE(x,v,t,z), severe symptomatic infected fI=fI(x,v,t,z), mildly symptomatic or asymptomatic infected fA=fA(x,v,t,z) and removed (healed or deceased) fR=fR(x,v,t,z), the kinetic distribution of the population results

    f(x,v,t,z)=fS(x,v,t,z)+fE(x,v,t,z)+fI(x,v,t,z)+fA(x,v,t,z)+fR(x,v,t,z).

    In addition to (1), we have that

    E(x,t,z)=11fE(x,v,t,z)dv,A(x,t,z)=11fA(x,v,t,z)dv.

    In this setting, the kinetic densities satisfy the following transport equations [7]

    fSt+vSfSx=F(fS,I)FA(fS,A)+1τS(S2fS)fEt+vEfEx=F(fS,I)+FA(fS,A)afE+1τE(E2fE)fIt+vIfIx=aσfEγIfI+1τI(I2fI)fAt+vAfAx=a(1σ)fEγAfA+1τA(A2fA)fRt+vRfRx=γIfI+γAfA+1τR(R2fR) (15)

    The quantities γI=γI(x,z) and γA=γA(x,z) are the recovery rates of symptomatic and asymptomatic infected (inverse of the infectious periods), respectively, while a=a(x,z) represents the inverse of the latency period and σ=σ(x,z) is the probability rate of developing severe symptoms [41,23,14]. In this model, the transmission of the infection is governed by two different incidence functions, F(,I) and FA(,A), simply to distinguish between the behavior of I and A individuals. Analogously to (3),

    FA(g,A)=βAgAp1+κAA, (16)

    where a different contact rate, βA=βA(x,t,z), and coefficient κA=κA(x,t,z) are taken into account for mildly/no symptomatic people. For the derivation of the reproduction number R0 of this SEIAR kinetic model, which results

    R0(t,z)=ΩFI(S,I)dxΩγI(x,z)I(x,t,z)dxΩa(x,z)σ(x,z)E(x,t,z)dxΩa(x,z)E(x,t,z)dx+ΩFA(S,A)dxΩγA(x,z)A(x,t,z)dxΩa(x,z)(1σ(x,z))E(x,t,z)dxΩa(x,z)E(x,t,z)dx, (17)

    the reader is invited to refer to [9,7].

    When introducing the same definition (5) of flux for the additional compartments, JE and JA, integrating system (2) in v, we get the following set of equations for the macroscopic densities

    St+JSx=F(S,I)FA(S,A)Et+JEx=F(S,I)+FA(S,A)aEIt+JIx=aσEγIIAt+JAx=a(1σ)EγAARt+JRx=γII+γAA (18)

    and for the fluxes

    JSt+λS211v2fSxdv=F(JS,I)FA(JS,A)JSτSJEt+λE211v2fExdv=λEλS(F(JS,I)+FA(JS,A))aJEJEτEJIt+λI211v2fIxdv=λIλEaσJEγIJIJIτIJAt+λA211v2fAxdv=λAλEa(1σ)JEγAJAJAτAJRt+λR211v2fRxdv=λRλIγIJI+λRλAγAJAJRτR. (19)

    Moreover, defining also DE=13λ2EτE and DA=13λ2AτA and following the same procedure discussed in Section 2.1.1, we recover this SEIAR system in the diffusive regime:

    St=F(S,I)FA(S,A)+x(DSSx)Et=F(S,I)+FA(S,A)aE+x(DEEx)It=aσEγII+x(DIIx)At=a(1σ)EγAA+x(DAAx)Rt=γII+γAA+x(DRRx). (20)

    The low-fidelity model is obtained, also with the more complex SEIAR compartmentalization, considering the discrete-velocity case with only two opposite velocities ±λ [7]:

    S±t±λSS±x=F(S±,I)FA(S±,A)+12τS(SS±)E±t±λEE±x=F(S±,I)+FA(S±,A)aE±+12τE(EE±)I±t±λII±x=aσE±γII±+12τI(II±)A±t±λAA±x=a(1σ)E±γAA±+12τA(AA±)R±t±λRR±x=γII±+γAA±+12τR(RR±). (21)

    When defining fluxes JE and JA as in (12), the equivalent hyperbolic model underlying the macroscopic formulation of the spatial propagation of an epidemic is obtained. The evolution of the densities follows (18), while for the fluxes we have

    JSt+λ2SSx=F(JS,I)FA(JS,A)JSτSJEt+λ2EEx=λEλS(F(JS,I)+FA(JS,A))aJEJEτEJIt+λ2IIx=λIλEaσJEγIJIJIτIJAt+λ2AAx=λAλEa(1σ)JEγAJAJAτAJRt+λ2RRx=λRλIγIJI+λRλAγAJAJRτR. (22)

    As for the previous cases, the diffusion limit of the system is formally recovered letting the relaxation times τS,E,I,A,R0, while keeping the diffusion coefficients, Di=λ2iτi,i{S,E,I,A,R}, finite. Thus, the same procedure presented for the SIR compartmentalization in Section 2.2.1 lead to the SEIAR parabolic system (20).

    In this Section, we present the details of the numerical method used to solve the stochastic problem following the bi-fidelity asymptotic-preserving scheme proposed in [35,48]. For the high-fidelity model, the numerical scheme is structured in agreement with a discrete ordinate method in velocity with the even and odd parity formulation [20,28]. Both the high-fidelity and low-fidelity solvers use Finite Volume Method (FVM) in space and achieve asymptotic preservation in time using suitable IMEX Runge-Kutta schemes [11,12]. This permits to obtain a numerical scheme able to deal with the diffusion limit of the mathematical models without loosing consistency, and for which the time step size of the temporal discretization is not subject to excessive restrictions related to the smallness of the scaling parameters τi. For simplicity, we illustrate the numerical method in the case of the simpler SIR model. The application of the same numerical scheme results straightforward for the case of the SEIAR compartmentalization.

    In this Section, we present the AP-IMEX Finite Volume scheme adopted to solve the SIR model at each stochastic collocation point selected for the bi-fidelity approximation.

    The asymptotic-preserving IMEX method and the corresponding even and odd parities formulation was introduced in [13] for an SIR kinetic transport model in 2D domains. According to [28,24], for v>0, we can define the even and odd parities for the high-fidelity SIR kinetic transport model (2) as follows:

    rS(v)=12(fS(v)+fS(v)),jS(v)=λS2(fS(v)fS(v)),rI(v)=12(fI(v)+fI(v)),jI(v)=λI2(fI(v)fI(v)),rR(v)=12(fR(v)+fR(v)),jR(v)=λR2(fR(v)fR(v)).

    An equivalent formulation of (2) can be written as

    rSt+vjSx=F(rS,I)+1τS(12SrS)rIt+vjIx=F(rS,I)γrI+1τI(12IrI)rRt+vjRx=γrI+1τR(12RrR)jSt+λ2SvrSx=F(jS,I)1τSjSjIt+λ2IvrIx=F(jS,I)γjI1τIjIjRt+λ2RvrRx=γjI1τRjR, (23)

    where

    S=210rSdv,I=210rIdv,R=210rRdv. (24)

    The above densities can be approximated by a Gauss-Legendre quadrature rule. This leads to a discrete velocity setting, usually referred to as the discrete ordinate method, where we approximate

    SSM=NGi=1wirS(ζi)IIM=NGi=1wirI(ζi)RRM=NGi=1wirR(ζi)

    where wi and ζi are the NG standard Gauss-Legendre quadrature weights and points in [1,1], and NG=Nv, number of chosen discrete velocities.

    Assuming for simplicity of notation that τS,I,R=τ, we can write the above system in the following compact form

    rt+vjx=E(r)1τ(rR2)jt+Λ2vrx=E(j)1τj, (25)

    with

    r=(rS,rI,rR)T,j=(jS,jI,jR)T,R=(S,I,R)T,Λ=diag{λS,λI,λR},
    E(r)=(F(rS,I),F(rS,I)γrI,γrI)T,E(j)(F(jS,I),F(jS,I)γjI,γjI)T.

    Following [12], the Implicit-Explicit (IMEX) Runge-Kutta discretization that we consider for system (25) consists in computing the internal stages

    r(k)=rnΔtkj=1akj(vj(j)x+1τ(r(j)R(j)2))+Δtk1j=1˜akjE(r(j))j(k)=jnΔtk1j=1˜akj(Λ2vr(j)xE(j(j)))Δtkj=1akj1τj(j), (26)

    and evaluating the final numerical solution

    rn+1=rnΔtsk=1bk(vj(k)x+1τ(r(k)R(k)2))+Δtsk=1˜bkE(r(k))jn+1=jnΔtsk=1˜bk(Λ2vr(k)xE(j(k)))Δtsk=1bk1τj(k). (27)

    To properly compute the implicit terms R(j)2 at each Runge-Kutta internal step explicitly, we refer to [27]. Matrices ˜A=(˜akj), with ˜akj=0 for jk, and A=(akj), with akj=0 for j>k are s×s matrices, with s number of Runge-Kutta stages, defining respectively the explicit and the implicit part of the scheme, and vectors ˜b=(˜b1,...,˜bs)T and b=(b1,...,bs)T are the quadrature weights. Furthermore, referring to [11,12], if the following relations hold,

    akj=bj,j=1,,s,˜akj=˜bj,j=1,,s1,

    the method is said to be globally stiffly accurate (GSA). It is worth to notice that this definition states also that the numerical solution of a GSA IMEX Runge-Kutta scheme coincides exactly with the last internal stage of the scheme. Since the GSA property is fundamental to preserve the correct diffusion limit and to achieve asymptotic-preservation stability in stiff regimes [13,10], in the sequel, the GSA BPR(4, 4, 2) scheme presented in [12] is chosen, characterized by s=4 stages for the implicit part, 4 stages for the explicit part and 2nd order of accuracy.

    At each internal stage of the IMEX scheme (26), we apply a Total-Variation-Diminishing (TVD) Finite Volume discretization to evaluate the numerical fluxes [10,8]. To achieve second order accuracy also in space, while avoiding the occurrence of spurious oscillations, a classical minmod slope limiter has been adopted.

    The same AP-IMEX Finite Volume scheme is adopted also to solve the low-fidelity SIR model (6)-(13), as fully presented in [10]. The reader is invited to refer to [10,9,7,13] for further details on the properties of the chosen numerical scheme applied to epidemic models.

    In this Section, we briefly review the bi-fidelity method developed in [35,48]. Bi-fidelity methods make use of low-fidelity models to effectively inform the selection of representative points in the parameter space and then employ this information to construct accurate approximations to high-fidelity solutions. To facilitate future discussion, we denote the expensive high-fidelity solution uH(z) and the cheap low-fidelity solution uL(z) for any given random parameter zIzRd, where Iz is the domain of the random parameter z.

    The basic idea of the bi-fidelity approximation is to construct an inexpensive surrogate uB(z) of the high-fidelity solution in the following non-intrusive manner:

    uB(z)=nk=1ck(z)uH(zk), (28)

    where n is the number of the selected parameter points in the parameter space. If n is small and the coefficient ck(z) can be efficiently and accurately approximated, an efficient bi-fidelity approximation can be constructed. Once the surrogate is constructed, the statistics of high-fidelity solutions can be quickly approximated by evaluating the bi-fidelity surrogate via Monte Carlo methods or quadrature rules.

    There are two major key questions for the performance of the above bi-fidelity algorithms: (a) how to select the representative points zk effectively? (b) how to efficiently construct the bi-fidelity approximation for any new given z but avoiding requiring a high-fidelity simulation?

    Subset Selection. Existing predefined or structured nodes (e.g., sparse grids, cubature rules, etc.) often grows fast in high dimensions. Therefore, these options cannot easily accommodate the current situation where we would like the size n of γn to be small and also arbitrary. Alternatively, adaptive approaches based on high-fidelity models to explore the parameter space suit our needs. However, this typically requires a large number of high-fidelity samples, which might not be computationally affordable. In contrast, the low-fidelity model is inexpensive to evaluate and it mimics the variations of high-fidelity solutions in the parameter space. This motivates us to employ the inexpensive low-fidelity model to learn and explore the behaviors of the high-fidelity model in the parameter space.

    We shall identify important points iteratively by a greedy approach [35,48]. Specifically, we denote a candidate sample set ΓN={z1,z2,,zN}, which is assumed to be large enough to cover the parameter space Iz. Initially, denote γ0={} and assume we have the first k important points γk={zi1,zi2,,zik} available at the k-th iteration. Denote the snapshot matrix uL(γk)={u(z)|zγk} and the corresponding spanned approximation space UL(γk)=span{uL(z)|zγk}. The corresponding high-fidelity approximation space can be defined similarly, UH(γn)=span{uH(z)|zγn}. Then we pick the point (from the candidate set ΓN) so that the corresponding low-fidelity solution is farthest away from the existing spanned low-fidelity approximation space UL(γk), to be the next sampling point:

    zik+1=argmaxzΓNdL(uL(z),UL(γk)),γk+1=γkzik+1, (29)

    where dL(v,W) is the distance between a function vuL(ΓN) and the space WuL(γk). We then repeat this step to select all n important points. The whole procedure can be efficiently implemented by apply the pivoted Cholesky decomposition on uL(ΓN) [35,48].

    Bi-fidelity approximation. For any given z, to efficiently compute bi-fidelity approximation in (28), it is desirable to find a cheap yet reasonably accurate approximation of ck(z). The bi-fidelity approach developed in [35,48] learns these coefficients from the inexpensive low-fidelity model uL(z). Specifically, for any given z, we shall compute the low-fidelity solution uL(z) and then construct its best approximation in the spanned low-fidelity approximation space UL(γn) by orthogonal projection:

    uL(z)PUL(γn)uL(z)=nk=1cLk(z)uL(zk),zkγn, (30)

    where the low-fidelity projection coefficients can be computed as follows:

    GLcL=f,fL=(fLk)1kn,fLk=uL(z),uL(zk), (31)

    and GL is the Gramian matrix of uL(γn),

    (GL)ij=uL(zi),uL(zj)L,1i,jn, (32)

    where ,L is the standard inner product associated with UL(γn).

    Under certain conditions [35,48], cLk(z) can be reasonably good approximation for the high-fidelity coefficients in (28). Consequently, we can construct the bi-fidelity approximation of the high-fidelity approximation solution uH(z) as follows:

    uB(z)=nk=1cLk(z)uH(zk). (33)

    Algorithm 1:A bi-fidelity approximation for a high-fidelity solution at given z
    1Given a candidate sample set ΓN={z1,z2,,zN}Iz, run the low-fidelity model uL(zj) for each zjΓN.
    2Select n "important" points γN from ΓN, where γN={zi1,zin}ΓN and the low-fidelity approximation space by UL(γn).
    3Run high-fidelity simulation only at the point in the selected sample set γn.
    4For any given z, run the low-fidelity solver to get the low-fidelity solution uL(z) and compute its low-fidelity coefficients in (30):
    uL(z)PUL(γn)uL(z)=nk=1cLk(z)uL(zk),zkγn.
    5Construct the bi-fidelity approximation by applying the same approximation rule as in low-fidelity model with (33):
    uB(z)=nk=1cLk(z)uH(zk).

    To put things together, we outline the major steps for the bi-fidelity approximation of the high-fidelity sample for a given z in Algorithm 1, reported in the following.

    Remark 2. As we mentioned above, the bi-fidelity method relies on the assumption that the low-fidelity coefficients are similar to the high-fidelity coefficients in the parameter space under certain conditions, as stated in [35,48]. This may not hold for some problems. In this case, a correction mapping between low-fidelity and high-fidelity coefficient can be constructed by leveraging approximation power of neural network to further improve the accuracy, if the additional high-fidelity data are available. We refer readers to [33] for details on this approach.

    To construct the bi-fidelity approximation of the high-fidelity mean, with the bi-fidelity surrogate uB(z), we can employ the Monte Carlo or other quadrature-based methods to compute the statistical moments quickly by evaluating bi-fidelity surrogates. To further reduce the number of bi-fidelity surrogate evaluations, a more efficient extension to this approach is developed in [47]. The general procedure is as follows:

    ● Compute the low-fidelity sample mean (via the Monte Carlo or quadrature rules):

    μL=Mi=1wiuL(zi), (34)

    where the high-fidelity mean μH is defined similarly.

    ● Compute its best approximation on the low-fidelity approximation space UL(γn) by orthogonal projection,

    μLPUL(γn)μL=ni=1cLiuL(zi), (35)

    where the expansion coefficients cL are computed by solving the following linear system:

    GLcL=gL,gL=μL,uL(zj)L,1jn,zjγn. (36)

    ● Using this coefficient cLk as the surrogate of the high-fidelity coefficient of the high-fidelity mean μH, the bi-fidelity approximation of the high-fidelity mean can be constructed as follows:

    μB=nk=1cLkuH(zk),zkγn. (37)

    Note that in this way, only one bi-fidelity surrogate evaluation is required. The standard deviation can be computed similarly. We refer readers to [47] for additional details.

    To examine the performance of the proposed methodology, two benchmark tests are considered: the first concerning the kinetic transport model with the SIR compartmentalization discussed in Section 2, and the second one regarding the extension to the SEIAR modeling presented in Section 3. It is worth to highlight that, in these tests, we collect all the quantities of interests together in a single vector to choose the zk points of the bi-fidelity algorithm, for either the 3 compartments S,I,R in Test 1 or the five compartments S, E, I, A, R in Test 2. From the experiments, we also find that the results are similar if we choose the zk points separately for the different epidemic compartments.

    We first consider the initial distributions of the high-fidelity kinetic SIR model (2) as follows:

    fi(x,v,0)=ci(x,0)ev22,i{S,I,R} (38)

    where c=12NGi=1wieζ2i2 is a re-normalization constant, with NG number of Gauss-Legendre quadrature points as defined in Section 4.1, and

    S(x,0)=1I(x,0),I(x,0)=0.01e(x10)2,R(x,0)=0,

    in the physical domain L=[0,20]. The initial fluxes JS(x,0),JI(x,0) and JR(x,0) are null and we consider periodic boundary conditions. The same initial conditions for S,I,R and JS,JI,JR are imposed in the low-fidelity SIR model (6)-(13). In the following tests, we consider a 2-dimensional random vector z=(z1,z2)T, with independent random parameters z1 and z2 that follow a uniform distribution, zjU(1,1), j=1,2.

    We analyze the behavior of the proposed methodology concerning spatially heterogeneous environments, considering a spatially variable contact rate [10,44]

    β(x,z)=β0(z)(1+0.05sin(13πx20)),

    where

    β0(z)=11(1+0.6z1),

    and a recovery rate

    γ(z)=10(1+0.4z2).

    In the incidence function, we set κ=0 and p=1. Therefore, we simulate an infectious disease characterized by an R0 perturbed around the value 1 depending on the random fluctuations. We discretize the spatial domain, both in the high-fidelity and low-fidelity models, with Nx=150 cells and consider Nv=8 velocities for the high-fidelity model. To compute the reference solutions for the mean and standard deviation, we use a 3-rd level sparse grid quadrature based on Clenshaw-Curtis rules for the choice of the stochastic collocation nodes, with a total of 29 points, for both the high-fidelity and low-fidelity models.

    Test 1 (a): In this case, a parabolic configuration of speeds and relaxation parameters is considered, setting λ2i=105, i{S,I,R}, and τi=105 in the low-fidelity model and τi=3×105 in the high-fidelity model, to maintain consistency of the two simulations (see Remark 1). The time step size results Δt=0.89×102 in both models, nevertheless the low-fidelity code is almost 5 times faster than the high-fidelity one (tHFCPU70.0s, tLFCPU14.3s to simulate until t=5 in a machine with processor Intel(R) Core(TM) i7-9750H CPU @2.60 GHz – 6 cores and 12 logical processors).

    In the first row of Figure 1, the expectation and standard deviation of the solution of compartment I for the high-fidelity model, the low-fidelity model and the bi-fidelity approximation at time t=5 are shown. As expected, since high-fidelity and low-fidelity models share the same diffusive limit, a perfect agreement of the solutions can be observed. To confirm this, in the second row of the same figure, we also plot the L2 errors of the mean and standard deviation between the bi-fidelity and high-fidelity solutions at T=5 with respect to the number of selected "important" points n of the bi-fidelity algorithm. A fast error decay is clearly observed. With only n=8 hi-fidelity sample points, bi-fidelity approximation can achieved a relative error of O(106) for both the mean and standard deviation.

    Figure 1. 

    Test 1 (a): SIR model in diffusive regime. First row: expectation (left) and standard deviation (right) obtained at t=5 for the variable I with the three methodologies, by using n=8 points for the bi-fidelity approximation. Second row: relative L2 errors of the bi-fidelity approximation for the mean (left) and standard deviation (right) of density I with respect to the number of "important" points n used in the bi-fidelity algorithm, compared with low-fidelity errors

    .

    Test 1 (b): In this second case, we consider the hyperbolic regime by letting λi=1, i{S,I,R}, while τi=1 in the low-fidelity model and τi=3 in the high-fidelity model, for consistency. The time step size results Δt=0.12 in both models, but the low-fidelity simulation is again almost 5 times faster than the high-fidelity one (tHFCPU5.6s, tLFCPU1.3s to simulate until t=5 in the same machine previously mentioned).

    In Figure 2, the results of Test 1(b) are reported for the infectious compartment I at time t=5. The first row shows expectation and standard deviation of the high-fidelity, the low-fidelity and the bi-fidelity solutions. While low-fidelity model fails to capture details of the hi-fidelity solutions around peaks due to the different velocity setting considered in the two models, bi-fidelity approximation are in almost perfect agreement with the high-fidelity ones, confirming the validity of the proposed methodology even with hyperbolic configurations of the scaling parameters. In addition, we plot the relative L2 error convergence of the mean and standard deviation of the bi-fidelity and high-fidelity solutions at t=5 with the number of high-fidelity samples, in the second row of Figure 2. Again, we observe a fast error decay for mean and standard deviation with O(10) high-fidelity samples.

    Figure 2. 

    Test 1 (b): SIR model in hyperbolic regime. First row: expectation (left) and standard deviation (right) obtained at t=5 for the variable I with the three methodologies, by using n=14 points for the bi-fidelity approximation. Second row: relative L2 errors of the bi-fidelity approximation for the mean (left) and standard deviation (right) of density I with respect to the number of "important" points n used in the bi-fidelity algorithm, compared with low-fidelity errors

    .

    Next, we analyze the effectiveness of the proposed methodology also with the extended SEIAR compartmentalization examining a more realistic epidemic scenario. Let us now consider a 2-dimensional random vector z, this time with zjU(0,1), j=1,2. We design an initial condition for the low-fidelity SEIAR model (18)-(22) that simulates the presence of 3 cities aligned in the spatial domain L=[0,20] with a different number of exposed (infected but not yet infectious) individuals, subjected to randomness, distributed following a Gaussian function,

    E(x,0,z)=α1(z)e(xx1)2+α2(z)e(xx2)2+α3(z)e(xx3)2,

    where x1=10/3, x2=10, x3=50/3 are the coordinates of the city centers and α1=0.01(1+z1), α2=0.001(1+z1), α3=0.004(1+z1) are the stochastic amplitudes. In this case, we set the large stochastic amplitude to model large uncertainty present on the initial amount of exposed people when an epidemic starts spreading. Indeed, it is worth to underline that the initial number of exposed individuals is certainly one of the variables most affected by uncertainty during epidemic outbreaks [23]. We consider that there are no infectious people at t=0, being

    S(x,0,z)=1E(x,0,z),I(x,0)=0,A(x,0)=0,R(x,0)=0,

    and JS(x,0)=JE(x,0)=JI(x,0)=JA(x,0)=JR(x,0)=0, with periodic boundary conditions, to allow a connection also between cities 1 and 3. The initial distributions of the high-fidelity kinetic SEIAR model (15) then read as defined in (38), for i{S,E,I,A,R}.

    Concerning epidemic parameters, to simulate the more challenging scenario in which the incidence function presents sinusoidal oscillations in space as well as being greater in the most populated areas, we consider the following distribution of the contact rate related to asymptomatic individuals:

    βA(x,z)=β0A(z)(1+12e(xx1)2+14e(xx2)2+12e(xx3)2)+0.05sin(2πx),

    where β0A is affected by the following random fluctuations:

    β0A(z)=0.5(1+0.5z2).

    Assuming that highly infectious subjects are mostly detected in the most optimistic scenario, being subsequently quarantined or hospitalized, we set βI(x,z)=0.03βA(x,z). Then, we fix γI=1/14, γA=1/7, a=1/3, σ=1/12.5, considering these clinical parameters deterministic according to values adopted in [9,7,14,23] to simulate the COVID-19 spread. Finally, in the incidence functions, we set κI=κA=0, hence assuming that initially individuals are not aware of the epidemic outbreak, and p=1 to work with the standard bi-linear case. The resulting setting permits to simulate and epidemic characterized by a baseline reproduction number R0=3. We consider Nx=150 cells for both high-fidelity and low-fidelity simulations, and use Nv=8 velocities for the high-fidelity solution.

    Test 2 (a): In the first case, an intermediate regime between parabolic and hyperbolic is considered, setting λ2i=10, i{S,E,A,R}, λI=0 and τi=0.25, i{S,E,I,A,R}, in the low-fidelity model and τi=0.75 in the high-fidelity model. The characteristic speed of compartment I is fixed to zero because we assumed that infectious people with severe symptoms are generally detected and isolated [9,7].

    The baseline space-time evolution of the compartments is shown in Figure 3, where the fast propagation of the epidemic can be observed, especially starting from the first city on the left, as expected. In Figure 4, the expectation and standard deviation of infected individuals, divided by each compartment, are presented for each methodology adopted, highlighting the validity of the bi-fidelity approach. In Figure 5, L2 error decay of the bi-fidelity approximations, with respect to the number of selected "important" points n, is compared with the error that would be obtained simply using the low-fidelity model.

    Figure 3. 

    Test 2 (a): SEIAR model in intermediate regime. The baseline temporal and spatial evolution of compartments S (first row, left), E (first row, right), I (second row, left) and A (second row, right) in the high-fidelity model

    .
    Figure 4. 

    Test 2 (a): SEIAR model in intermediate regime. Expectation (left) and standard deviation (right) of densities E (first row), I (second row) and A (third row) at time t=5, obtained with the three methodologies, using n=6 for the bi-fidelity solution

    .
    Figure 5. 

    Test 2 (a): SEIAR model in intermediate regime. Relative L2 error decay of the bi-fidelity approximation of expectation (left) and standard deviation (right) for the density A with respect to the number of selected "important" points n, compared with low-fidelity errors

    .

    Test 2 (b): In the second case, we consider a fully hyperbolic regime choosing λ2i=1, i{S,E,A,R}, λI=0, and τi=10, i{S,E,I,A,R}, and τi=30 respectively in the low-fidelity and high-fidelity model.

    The baseline space-time evolution of the compartments of major interest is presented in Figure 6. In this figure, the spatial heterogeneity of the epidemic spread related to the prescribed contact rate can be clearly appreciated. Mean and standard deviation obtained with the 3 approaches are shown for compartments E, I and A in Figure 7. Again, although the low-fidelity and high-fidelity solutions show clear discordant trends, the effectiveness of the bi-fidelity method is plainly confirmed. In Figure 8, the decay of the L2 relative error norm of the bi-fidelity approximation with respect to the number of points n shows that a proper accuracy can be achieved for both mean and standard deviation even in this challenging test case.

    Figure 6. 

    Test 2 (b): SEIAR model in hyperbolic regime. Baseline temporal and spatial evolution of compartments S (first row, left), E (first row, right), I (second row, left) and A (second row, right) in the high-fidelity model

    .
    Figure 7. 

    Test 2 (b): SEIAR model in hyperbolic regime. Expectation (left) and standard deviation (right) of densities E (first row), I (second row) and A (third row) at time t=5, obtained with the three methodologies, using n=7 for the bi-fidelity solution

    .
    Figure 8. 

    Test 2 (b): SEIAR model in hyperbolic regime. Relative L2 error decay of the bi-fidelity approximation of expectation (left) and standard deviation (right) for the density A with respect to the number of selected "important" points n, compared with low-fidelity errors

    .

    In this work we introduced a bi-fidelity method for the quantification of uncertainty in epidemiological transport models based on an asymptotic-preserving space-time discretization. In detail, after presenting the high-fidelity epidemiological model, we considered the corresponding bi-fidelity model. Both models share the same diffusive limit and permits to recover classical epidemic models based on diffusion equations. The numerical scheme used allows to have an efficient quantification of the uncertainty in different regimes, using few simulations of the high-fidelity model and several runs of the low-fidelity model for points selection in the random space. Results for one-dimensional transport problems based on realistic compartmentalization in relation to the recent COVID-19 pandemic, which also include asymptomatic individuals, show the validity of the presented approach. Further research will be directed toward extending the present approach to realistic contexts such as those studied in [7,9].

    G.B. and L.P. acknowledge the support of MIUR PRIN 2017, Project No. 2017KKJP4X Innovative numerical methods for evolutionary partial differential equations and applications. G.B. holds a Research Fellowship from the Italian National Institute of High Mathematics, INdAM (GNCS). L.L. holds the Direct Grant for Research supported by Chinese University of Hong Kong and Early Career Scheme 2021/22, No. 24301021, from Research Grants Council of Hong Kong. X.Z. was supported by the Simons Foundation (504054).



    [1] A. Almeida, D. Drihem, Maximal, potential and singular type operators on Herz spaces with variable exponents, J. Math. Anal. Appl., 394 (2012), 965–974. https://doi.org/10.1016/j.jmaa.2012.04.043 doi: 10.1016/j.jmaa.2012.04.043
    [2] D. Cruz-Uribe, L. Diening, P. Hästö, The maximal operator on weighted variable Lebesgue spaces, Fract. Calc. Appl. Anal., 14 (2011), 361–374. https://doi.org/10.2478/s13540-011-0023-7 doi: 10.2478/s13540-011-0023-7
    [3] D. Cruz-Uribe, A. Fiorenza, J. M. Martell, C. Pérez, The boundedness of classical operators on variable Lp spaces, Ann. Acad. Sci. Fenn. Math., 31 (2006), 239–264. https://doi.org/10.1146/annurev.energy.31.020105.100323 doi: 10.1146/annurev.energy.31.020105.100323
    [4] D. Cruz-Uribe, A. Fiorenza, Variable Lebesgue spaces: Foundations and harmonic analysis, Applied and Numerical Harmonic Analysis, Heidelberg: Birkhäuser Basel, 2013. https://doi.org/10.1007/978-3-0348-0548-3
    [5] D. Cruz-Uribe, A Fiorenza, C. J. Neugebauer, Weighted norm inequalities for the maximal operator on variable Lebesgue spaces, J. Math. Anal. Appl., 394 (2012), 744–760. https://doi.org/10.1016/j.jmaa.2012.04.044 doi: 10.1016/j.jmaa.2012.04.044
    [6] D. Cruz-Uribe, L. A. Wang, Extrapolation and weighted norm inequalities in the variable Lebesgue spaces, Trans. Amer. Math. Soc., 369 (2017), 1205–1235. https://doi.org/10.1090/tran/6730 doi: 10.1090/tran/6730
    [7] L. Diening, P. Hästö, A. Nekvinda, Open problems in variable exponent Lebesgue and Sobolev spaces, Fsdona04 Proc., 2004, 38–58.
    [8] L. Diening, P. Harjulehto, P. Hästö, M. Ruzicka, Lebesgue and Sobolev spaces with variable exponents, Springer, 2011. https://doi.org/10.1007/978-3-642-18363-8
    [9] H. G. Feichtinger, F. Weisz, Herz spaces and summability of Fourier transforms, Math. Nachr., 281 (2008), 309–324. https://doi.org/10.1002/mana.200510604 doi: 10.1002/mana.200510604
    [10] L. Grafakos, X. Li, D. Yang, Bilinear operators on Herz-type Hardy spaces, Trans. Amer. Math. Soc., 350 (1998), 1249–1275. https://doi.org/10.1090/S0002-9947-98-01878-9 doi: 10.1090/S0002-9947-98-01878-9
    [11] E. Hernández, D. Yang, Interpolation of Herz spaces and applications, Math. Nachr., 205 (1999), 69–87. https://doi.org/10.1002/mana.3212050104 doi: 10.1002/mana.3212050104
    [12] C. S. Herz, Lipschitz spaces and Bernstein's theorem on absolutely convergent Fourier transforms, J. Math. Mech., 18 (1968), 283–323. https://doi.org/10.1512/iumj.1969.18.18024 doi: 10.1512/iumj.1969.18.18024
    [13] M. Izuki, Boundedness of sublinear operators on Herz spaces with variable exponent and application to wavelet characterization, Anal. Math., 36 (2010), 33–50. https://doi.org/10.1007/s10476-010-0102-8 doi: 10.1007/s10476-010-0102-8
    [14] M. Izuki, Boundedness of vector-valued sublinear operators on Herz-Morrey spaces with variable exponents, Math. Sci. Res. J., 13 (2009), 243–253.
    [15] M. Izuki, T. Noi, An intrinsic square function on weighted Herz spaces with variable exponents, arXiv, 2016. https://doi.org/10.48550/arXiv.1606.01019
    [16] M. Izuki, T. Noi, Boundedness of fractional integrals on weighted Herz spaces with variable exponents, J. Inequal. Appl., 2016 (2016), 199.
    [17] R. Johnson, Lipschitz spaces, Littlewood-Paley spaces, and convoluteurs, Proc. Lond. Math. Soc., 3 (1974), 127–141. https://doi.org/10.1112/plms/s3-29.1.127
    [18] R. Johnson, Temperatures, Riesz potentials, and the Lipschitz spaces of Herz, Proc. Lond. Math. Soc., 3 (1973), 290–316. https://doi.org/10.1112/plms/s3-27.2.290 doi: 10.1112/plms/s3-27.2.290
    [19] V. Kokilashvili, S. Samko, On Sobolev theorem for Riesz-type potentials in Lebesgue spaces with variable exponents, Z. Anal. Anwend., 22 (2003), 899–910. https://doi.org/10.4171/zaa/1178 doi: 10.4171/zaa/1178
    [20] V. M. Kokilashvili, A. Meskhi, H. Rafeiro, S. Samko, Integral operators in non-standard function spaces, Vol. 1, Variable Exponent Lebesgue and Amalgam Spaces Birkhäuser Cham, 2016. https://doi.org/10.1007/978-3-319-21015-5
    [21] Y. Komori, Notes on singular integrals on some inhomogeneous Herz spaces, Taiwan. J. Math., 8 (2004), 547–556. https://doi.org/10.11650/twjm/1500407672 doi: 10.11650/twjm/1500407672
    [22] A. Meskhi, H. Rafeiro, M. A. Zaighum, On the boundedness of Marcinkiewicz integrals on continual variable exponent Herz spaces, Georgian Math. J., 26 (2019), 105–116. https://doi.org/10.1515/gmj-2017-0050 doi: 10.1515/gmj-2017-0050
    [23] C. B. Morrey, On the solutions of quasi-linear elliptic partial differential equations, Trans. Amer. Math. Soc., 43 (1938), 126–166. https://doi.org/10.2307/1989904 doi: 10.2307/1989904
    [24] B. Muckenhoupt, Weighted norm inequalities for the Hardy maximal function, Trans. Amer. Math. Soc., 165 (1972), 207–226.
    [25] H. Nafis, H. Rafeiro, M. Zaighum, A note on the boundedness of sublinear operators on grand variable Herz spaces, J. Inequal. Appl., 2020 (2020), 1. https://doi.org/10.1186/s13660-019-2265-6 doi: 10.1186/s13660-019-2265-6
    [26] H. Nafis, H. Rafeiro, M. A. Zaighum, Boundedness of the Marcinkiewicz integral on grand Herz spaces, J. Math. Inequal., 15 (2021), 739–753. http://dx.doi.org/10.7153/jmi-2021-15-52 doi: 10.7153/jmi-2021-15-52
    [27] H. Rafeiro, S. Samko, Riesz potential operator in continual variable exponents Herz spaces, Math. Nachr., 288 (2015), 465–475. https://doi.org/10.1002/mana.201300270 doi: 10.1002/mana.201300270
    [28] S. Samko, Variable exponent Herz spaces, Mediterr. J. Math., 10 (2013), 2007–2025. https://doi.org/10.1007/s00009-013-0285-x doi: 10.1007/s00009-013-0285-x
    [29] S. Umarkhadzhiev, The boundedness of the Riesz potential operator from generalized grand Lebesgue spaces to generalized grand Morrey spaces, In: M. Bastos, A. Lebre, S. Samko, I. Spitkovsky, Operator theory, operator algebras and applications, Vol. 242, Birkhäuser, Basel, 2014,363–373. https://doi.org/10.1007/978-3-0348-0816-3_22
    [30] J. L. Wu, W. J. Zhao, Boundedness for higher order commutators of fractional integrals on variable exponent Herz–Morrey spaces, Mediterr. J. Math., 14 (2017), 198. https://doi.org/10.1007/s00009-017-1002-y doi: 10.1007/s00009-017-1002-y
    [31] S. Polidoro, M. A Ragusa, Harnack inequality for hypoelliptic ultraparabolic equations with a singular lower order term, Rev. Mat. Iberoamericana, 24 (2008), 1011–1046.
    [32] M. A. Ragusa, Homogeneous Herz spaces and regularity results, Nonlinear Anal.: Theory Methods Appl., 71 (2009), e1909–e1914. https://doi.org/10.1016/j.na.2009.02.075 doi: 10.1016/j.na.2009.02.075
    [33] A. Scapellato, Regularity of solutions to elliptic equations on Herz spaces with variable exponents, Bound. Value Probl., 2019 (2019), 2. https://doi.org/10.1186/s13661-018-1116-6 doi: 10.1186/s13661-018-1116-6
    [34] A. Scapellato, Homogeneous Herz spaces with variable exponents and regularity results, Electron. J. Qual. Theory Differ. Equ., 2018 (2018), 1–11. https://doi.org/10.14232/ejqtde.2018.1.82 doi: 10.14232/ejqtde.2018.1.82
    [35] A. Abdalmonem, A. Scapellato, Fractional operators with homogeneous kernels in weighted Herz spaces with variable exponent, J. Appl. Anal., 101 (2022), 1953–1962. https://doi.org/10.1080/00036811.2020.1789602 doi: 10.1080/00036811.2020.1789602
  • This article has been cited by:

    1. Giulia Bertaglia, Chuan Lu, Lorenzo Pareschi, Xueyu Zhu, Asymptotic-Preserving Neural Networks for multiscale hyperbolic models of epidemic spread, 2022, 32, 0218-2025, 1949, 10.1142/S0218202522500452
    2. Rinaldo M. Colombo, Mauro Garavello, Infectious Disease Spreading Fought by Multiple Vaccines Having a Prescribed Time Effect, 2023, 71, 0001-5342, 10.1007/s10441-022-09452-4
    3. Giulia Bertaglia, Andrea Bondesan, Diletta Burini, Raluca Eftimie, Lorenzo Pareschi, Giuseppe Toscani, New trends on the systems approach to modeling SARS-CoV-2 pandemics in a globally connected planet, 2024, 34, 0218-2025, 1995, 10.1142/S0218202524500301
    4. Giulia Bertaglia, 2023, Chapter 2, 978-3-031-29874-5, 23, 10.1007/978-3-031-29875-2_2
    5. Giulia Bertaglia, Lorenzo Pareschi, Multiscale Constitutive Framework of One-Dimensional Blood Flow Modeling: Asymptotic Limits and Numerical Methods, 2023, 21, 1540-3459, 1237, 10.1137/23M1554230
    6. Xiaojing Xu, Wenjun Sun, Qi Li, The asymptotic preserving unified gas kinetic scheme for the multi-scale kinetic SIR epidemic model, 2024, 174, 08981221, 298, 10.1016/j.camwa.2024.09.021
    7. Liu Liu, 2023, Chapter 7, 978-3-031-29874-5, 139, 10.1007/978-3-031-29875-2_7
    8. Giulia Bertaglia, Lorenzo Pareschi, Giuseppe Toscani, Modelling contagious viral dynamics: a kinetic approach based on mutual utility, 2024, 21, 1551-0018, 4241, 10.3934/mbe.2024187
    9. Giulia Bertaglia, Lorenzo Pareschi, Russel E. Caflisch, Gradient-Based Monte Carlo Methods for Relaxation Approximations of Hyperbolic Conservation Laws, 2024, 100, 0885-7474, 10.1007/s10915-024-02614-1
  • Reader Comments
  • © 2023 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(2105) PDF downloads(138) Cited by(16)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog