Research article Special Issues

Global weak solutions and asymptotic limits of a Cahn–Hilliard–Darcy system modelling tumour growth

  • Received: 30 August 2016 Accepted: 17 October 2016 Published: 28 October 2016
  • We study the existence of weak solutions to a Cahn-Hilliard-Darcy system coupled witha convection-reaction-diffusion equation through the fluxes, through the source terms and in Darcy’slaw. The system of equations arises from a mixture model for tumour growth accounting for transportmechanisms such as chemotaxis and active transport. We prove, via a Galerkin approximation, theexistence of global weak solutions in two and three dimensions, along with new regularity results forthe velocity field and for the pressure. Due to the coupling with the Darcy system, the time derivativeshave lower regularity compared to systems without Darcy flow, but in the two dimensional case weemploy a new regularity result for the velocity to obtain better integrability and temporal regularityfor the time derivatives. Then, we deduce the global existence of weak solutions for two variantsof the model; one where the velocity is zero and another where the chemotaxis and active transportmechanisms are absent.

    Citation: Harald Garcke, Kei Fong Lam. Global weak solutions and asymptotic limits of a Cahn–Hilliard–Darcy system modelling tumour growth[J]. AIMS Mathematics, 2016, 1(3): 318-360. doi: 10.3934/Math.2016.3.318

    Related Papers:

    [1] Imtiaz Ahmad, Muhammad Nawaz Khan, Rashid Jan, Normy Norfiza Abdul Razak . Efficient numerical method for pricing multi-asset options with the time-fractional Black-Scholes model: focus on American and digital options. Mathematical Modelling and Control, 2025, 5(2): 147-163. doi: 10.3934/mmc.2025011
    [2] Mrutyunjaya Sahoo, Dhabaleswar Mohapatra, S. Chakraverty . Wave solution for time fractional geophysical KdV equation in uncertain environment. Mathematical Modelling and Control, 2025, 5(1): 61-72. doi: 10.3934/mmc.2025005
    [3] Biresh Kumar Dakua, Bibhuti Bhusan Pati . A frequency domain-based loop shaping procedure for the parameter estimation of the fractional-order tilt integral derivative controller. Mathematical Modelling and Control, 2024, 4(4): 374-389. doi: 10.3934/mmc.2024030
    [4] Kexin Ouyang, Xinmin Qu, Huiqin Lu . Sign-changing and signed solutions for fractional Laplacian equations with critical or supercritical nonlinearity. Mathematical Modelling and Control, 2025, 5(1): 1-14. doi: 10.3934/mmc.2025001
    [5] Lusong Ding, Weiwei Sun . Neuro-adaptive finite-time control of fractional-order nonlinear systems with multiple objective constraints. Mathematical Modelling and Control, 2023, 3(4): 355-369. doi: 10.3934/mmc.2023029
    [6] Ihteram Ali, Imtiaz Ahmad . Applications of the nonlinear Klein/Sinh-Gordon equations in modern physics: a numerical study. Mathematical Modelling and Control, 2024, 4(3): 361-373. doi: 10.3934/mmc.2024029
    [7] Abduljawad Anwar, Shayma Adil Murad . On the Ulam stability and existence of Lp-solutions for fractional differential and integro-differential equations with Caputo-Hadamard derivative. Mathematical Modelling and Control, 2024, 4(4): 439-458. doi: 10.3934/mmc.2024035
    [8] Bengisen Pekmen, Ummuhan Yirmili . Numerical and statistical approach on chemotaxis-haptotaxis model for cancer cell invasion of tissue. Mathematical Modelling and Control, 2024, 4(2): 195-207. doi: 10.3934/mmc.2024017
    [9] Hui Li, Nana Jin, Yu Zhang . Existence of nonoscillatory solutions for higher order nonlinear mixed neutral differential equations. Mathematical Modelling and Control, 2024, 4(4): 417-423. doi: 10.3934/mmc.2024033
    [10] Xueling Fan, Ying Li, Wenxv Ding, Jianli Zhao . H-representation method for solving reduced biquaternion matrix equation. Mathematical Modelling and Control, 2022, 2(2): 65-74. doi: 10.3934/mmc.2022008
  • We study the existence of weak solutions to a Cahn-Hilliard-Darcy system coupled witha convection-reaction-diffusion equation through the fluxes, through the source terms and in Darcy’slaw. The system of equations arises from a mixture model for tumour growth accounting for transportmechanisms such as chemotaxis and active transport. We prove, via a Galerkin approximation, theexistence of global weak solutions in two and three dimensions, along with new regularity results forthe velocity field and for the pressure. Due to the coupling with the Darcy system, the time derivativeshave lower regularity compared to systems without Darcy flow, but in the two dimensional case weemploy a new regularity result for the velocity to obtain better integrability and temporal regularityfor the time derivatives. Then, we deduce the global existence of weak solutions for two variantsof the model; one where the velocity is zero and another where the chemotaxis and active transportmechanisms are absent.


    Record statistics are interest and importance in many areas of real-life applications including data relating to meteorology and weather, economics, sport, athletic events, oil, mining surveys and life testing studies and so on. In the midst of practical life there are many situations that depend entirely on the record statistics (record values) involving: Guinness World Records where only record values are observed (lower or upper record values). Shortest ever tennis matches both in terms of number of games (number of runs) and duration in terms of time (lower record values). News items like fastest time taken to recite the periodic table of the elements (lower record values). Weightlifting competitions where the highest weight is lifted in the fewest number of attempts (upper record values). Swimming competitions where the swimmer travels a certain distance in the least time of swimming (lower record values). Mountaineering, fastest marathon, longest time to hop, etc are of great importance to sports federations and athletes.

    From the above, we note that to create a record, several attempts are made, and when the attempt is successful, the record is created, and accordingly we do not get the data for all the attempts that were made to break records. The data that we have are records. Many scientists specially the statisticians have become interested in record values and discussed the statistical inferences for different models based on lower and upper record values including Chandler [1] studied the distribution and frequency of lower record values, Resnick [2] compared the behavior of the sequence of record values, Shorrock [3] proposed asymptotic results on record values and inter-record times, Glick [4] discussed the breaking records and breaking boards, Nevzorov [5] and Nagaraja [6] listed some theoretical characteristics for record values for various univariate model, Kumar [7] developed some recurrence relations for generalized logistic distribution based on lower record values, Mahmoud et al. [8] studied the estimation problem for the behavior of the coefficient of variation of Lomax distribution based on upper record values, Kumar and Saran [9] investigated the ratio and inverse moments of record values from Marshall-Olkin log-logistic distribution, and EL-Sagheer [10] studied the estimation of the parameters for generalized logistic distribution based on the lower record values. Arnold et al. [11] summarized the study of upper and lower record values and presented many of the basic properties of them as follows: Assume X1,X2,X3,... a sequence of independent and identically distributed random variables with probability density function (PDF) f(x) and cumulative distribution function (CDF) F(x). Setting Yj=max(X1,X2,X3,...,Xj), for j1, we say Xj an upper record and denoted by XU(j) if Yj > Yj1 , j>1. Also, let Yj=min(X1,X2,X3,...,Xj), for j1 we say that Xj is a lower record and denoted by XL(j) if Yj < Yj1 , j>1. Thus the joint density function of the first n upper record values is given by

    f1,2,3,...,n(xU(1),xU(2),xU(3),...,xU(n))=f(xU(n))n1i=1f(xU(i))1F(xU(i)), (1)

    and the joint density function of the first n lower record values is given by

    f1,2,3,...,n(xL(1),xL(2),xL(3),...,xL(n))=f(xL(n))n1i=1f(xL(i))F(xL(i)). (2)

    Dagum [12] introduced a statistical distribution closely fitting income distributions which accommodating the heavy tails, called the Dagum distribution (DD). In the actuarial literature, is also called the inverse Burr. Furthermore, DD admits a mixture representation in terms of generalized Gamma and inverse Weibull distributions. The DD can also be obtained as a compound generalized Gamma distribution whose scale parameter follows an inverse Weibull distribution with identical shape parameters. Since Dagum introduced his model as income distribution, many authors have discussed its properties in economics and financial fields, reliability, and survival data, etc. See for instance, Kleiber and Kotz [13], Kleiber [14], Domma et al. [15], Domma and Condino [16], among others. Recently, Alotaib et al. [17] discussed parameter estimation for the DD based on progressive type-Ⅰ interval censored data. A random variable X has the Dagum distribution and denoted by XD(λ,θ,β), if it has PDF, CDF, reliability function r(x) and hazard rate function h(x), given by

    f(x)=λθβx(β+1)(1+λxβ)(θ+1),   x>0, (3)
    F(x)=(1+λxβ)θ,   x>0, (4)
    r(x)=1(1+λxβ)θ,  x>0, (5)

    and

    h(x)=λθβx(β+1)(1+λxβ)(θ+1)1(1+λxβ)θ,  x>0, (6)

    where β, θ>0 are the shape parameters and λ>0 is the scale parameter.

    In this article, we investigate the Bayes and maximum likelihood estimates for the unknown quantities of the DD using lower record values. The Lindley's approximation and the Metropolis-Hasting (M-H) algorithm within Gibbs sampler are proposed to construct the Bayes estimates. To this end, the results of the Bayes estimates are obtained under both squared error (SE) and linear exponential (LINEX) loss functions. Also, the corresponding credible intervals (CRIs) are constructed under MCMC technique.

    The rest of this article is organized as follows. Section 2 deals with the maximum likelihood estimates and asymptotic confidence intervals. Bayesian estimates using Lindley's approximation and MCMC technique are provided in Section 3. In Section 4 a simulation study is provided to examine the performance of the proposed estimation methods. Section 5 discusses an application to real life data. Finally, a brief conclusion is given in Section 6.

    In this section, we discuss the maximum likelihood estimates (MLEs) and approximate confidence intervals (ACIs) of the parameters as well as reliability and hazard rate functions of DD given in (3) when the available data are lower record values. Suppose that x_=xL(1),xL(2),xL(3),...,xL(n) be the lower record values of size n from D(λ,θ,β), then the likelihood function for observed record x_ is given by

    L(λ,θ,β|x_)=λnβnθn(1+λxβL(n))θni=1x(β+1)L(i)(1+λxβL(i))1. (7)

    The logarithm of the likelihood function may then be written as

    (λ,θ,β|x_)=nlogλ+nlogβ+nlogθθlog(1+λxβL(n))(β+1)ni=1logxL(i)ni=1log(1+λxβL(i)). (8)

    The estimators ˆλ, ˆθ and ˆβ of the parameters λ, θ, and β, respectively, can be then obtained as the solution likelihood equations

    (λ,θ,β|x_)λ=nλθxβL(n)1+λxβL(n)ni=1xβL(i)1+λxβL(i)=0, (9)
    (λ,θ,β|x_)β=nβ+θλxβL(n)lnxL(n)1+λxβL(n)ni=1logxL(i)+ni=1xβL(i)lnxL(i)1+λxβL(i)=0, (10)

    and

    (λ,θ,β|x_)θ=nθlog(1+λxβL(n))=0. (11)

    Since (9)–(11) cannot be solved analytically for ˆλ, ˆθ and ˆβ, some numerical methods such as Newton-Raphson method must be employed to obtain the desired MLEs in such situations. Once we get desired MLEs of λ, θ and β using the invariant property of MLEs, for given mission time t, the MLEs of r(t) and h(t) can be obtained after replacing λ, θ and β by ˆλ, ˆθ and ˆβ as

    ˆr(t)=1(1+ˆλtˆβ)ˆθ, ˆh(t)=ˆλˆθˆβt(ˆβ+1)(1+ˆλtˆβ)(ˆθ+1)1(1+ˆλtˆβ)ˆθ, t>0. (12)

    From the log-likelihood function (λ,θ,β|x_) in (8), we have

    2(λ,θ,β|x_)λ2=nλ2+θ(xβL(n))2(1+λxβL(n))2+ni=1(xβL(i))2(1+λxβL(i))2, (13)
    2(λ,θ,β|x_)λθ=2(λ,θ,β|x_)θλ=xβL(n)1+λxβL(n), (14)
    2(λ,θ,β|x_)λβ=2(λ,θ,β|x_)βλ=θxβL(n)lnxL(n)[(1+λxβL(n))λxβL(n)](1+λxβL(n))2+ni=1xβL(i)lnxL(i)[(1+λxβL(i))λxβL(i)](1+λxβL(i))2, (15)
    2(λ,θ,β|x_)β2=nβ2+θλ(xβL(n))2(lnxL(n))2[λxβL(n)(1+λxβL(n))](1+λxβL(n))2+ni=1xβL(i)(lnxL(i))2[λxβL(i)(1+λxβL(i))](1+λxβL(i))2, (16)
    2(λ,θ,β|x_)βθ=2(λ,θ,β|x_)θβ=λxβL(n)lnxL(n)1+λxβL(n), 2(λ,θ,β|x_)θ2=nθ2. (17)

    The Fisher information matrix I(λ,θ,β) is then obtained by taking expectation of negative Eqs (13)–(17). Based on some mild regularity conditions, (ˆλ,ˆθ,ˆβ) is approximately bivariate normal with mean (λ,θ,β) and covariance matrix I1(λ,θ,β). Practically, we usually estimate I1(λ,θ,β) by I1(ˆλ,ˆθ,ˆβ). Moreover, a simpler and equally valid procedure is to use the approximation

    (ˆλ,ˆθ,ˆβ)N((λ,θ,β),I1(ˆλ,ˆθ,ˆβ)),

    where

    I(ˆλ,ˆθ,ˆβ)=(2(λ,θ,β|x_)λ22(λ,θ,β|x_)λθ2(λ,θ,β|x_)λβ2(λ,θ,β|x_)θλ2(λ,θ,β|x_)θ22(λ,θ,β|x_)θβ2(λ,θ,β|x_)βλ2(λ,θ,β|x_)βθ2(λ,θ,β|x_)β2)(ˆλ,ˆθ,ˆβ).

    ACIs for λ, θ and β can be obtained by to be bivariate normal distributed with mean (λ,θ,β) and variance-covariance matrix I1(λ,θ,β). Therefore, the 100(1γ)% ACIs for λ, θ and β, respectively, are

    ˆλzγ/2Λ11,  ˆθzγ/2Λ22,  ˆβzγ/2Λ33, (18)

    where Λ11, Λ22 and Λ33 are the elements on the main diagonal of the variance-covariance matrix I1(λ,θ,β), and zγ/2 is the percentile of the standard normal distribution with right-tail probability γ/2.

    In order to find the approximate estimates of the variance of r(t) and h(t), we use the delta method, see Agresti [18]. Assume that G1=(r(t)λ,r(t)θ,r(t)β) and G2=(h(t)λ,h(t)θ,h(t)β), where r(t)λ, r(t)θ, r(t)β, h(t)λ, h(t)θ, and h(t)β are the first partial derivatives of r(t) and h(t) with respect to the parameters λ, θ and β, respectively. The approximate asymptotic variances of ˆr(t) and ˆh(t) are given by

    Var(ˆr(t))[GT1I1(ˆλ,ˆθ,ˆβ)G1](ˆλ,ˆθ,ˆβ) and Var(ˆh(t))[GT2I1(ˆλ,ˆθ,ˆβ)G2](ˆλ,ˆθ,ˆβ),

    where and GTi is the transpose of Gi, i=1,2. Thus, the (1γ)100% ACIs for ˆr(t) and ˆh(t) can be written as

    ˆr(t)zγ/2Var(ˆr(t))  and  ˆh(t) zγ/2Var(ˆh(t)). (19)

    Besides being an accurate analysis and not easy to handle mathematically, the use of the Bayesian method allows the incorporation of previous knowledge of the parameters through informative prior densities. In the Bayesian approach, the information referring to the model parameters is obtained through posterior marginal distributions. In this section, two Bayesian inference procedures (Lindley's approximation and MCMC technique) are proposed to estimate the parameters λ, θ and β as well as the reliability r(t) and hazard rate h(t) functions. For estimating these quantities, we assume mainly both SE and LINEX loss functions. However, any other loss function can be easily implemented. From the perspective of frequentist, a natural choice for the prior distributions of λ, θ and β would be to assume that the three quantities are independent gamma distributions G(a1,b1), G(a2,b2) and G(a3,b3) respectively, where ai and bi, i=1,2,3 reflect the knowledge of prior about (λ,θ,β) and they are assumed to be known and nonnegative hyper-parameters. Further, the joint prior distribution can be written as

    g(λ,θ,β)λa11θa21βa31exp{(b1λ+b2θ+b3β)},λ>0,θ>0,β>0. (20)

    Subsequently, via Bayes' theorem across combining the likelihood function given (7) with the joint prior given in (20), the joint posterior density function of λ, θ and β can be written as follows

    g(α,β,λ|x_)λn+a11θn+a21βn+a31exp{b1λb3βθ[b2+log(1+λxβL(n))]}×ni=1x(β+1)L(i)(1+λxβL(i))1. (21)

    Thus, the conditional posterior densities of λ, θ and β can be computed and written, respectively, as

    g1(λ|θ,β,x_)λn+a11exp{b1λθlog(1+λxβL(n))}ni=1(1+λxβL(i))1, (22)
    g2(θ|λ,β,x_)θn+a21exp{θ[b2+log(1+λxβL(n))]}, (23)

    and

    g3(β|θ,λ,x_)βn+a31exp{b3βθlog(1+λxβL(n))}ni=1x(β+1)L(i)(1+λxβL(i))1. (24)

    Therefore, for any function, say u(λ,θ,β), the Bayes estimates under SE and LINEX loss functions can be written as

    E[u(λ,θ,β)|x_]=000u(λ,θ,β)L(λ,θ,β|x_)g(λ,θ,β)dλdθdβ000L(λ,θ,β|x_)g(λ,θ,β)dλdθdβ,E[ecu(λ,θ,β)|x_]=000ecu(λ,θ,β)L(λ,θ,β|x_)g(λ,θ,β)dλdθdβ000L(λ,θ,β|x_)g(λ,θ,β)dλdθdβ.}. (25)

    There is a difficulty to compute (25) analytically. Therefore, we propose to approximate it by using both Lindley's approximation and MCMC technique with Metropolis-Hasting algorithm within Gibbs sampler, to obtain the Bayes estimates of λ, θ, β, r(t) and h(t).

    There are various methods to approximate the ratio of integrals of the above form, may be the simplest one is Lindley approximation method, which was introduced by Lindley [19] to approximate the Bayes estimates into a form containing no integrals. This method has been used with various lifetime distributions by several authors for instance, Singh et al. [20], Ying and Gui [21], Ramin and Kohansal [22]. Briefly, it can be described as follows: If the sample size is sufficiently large, according to Lindley [19], any ratio of the integral of the form

    I(x)=E[u(λ,θ,β)|x_]=(λ,θ,β)u(λ,θ,β)e(λ,θ,β)+ρ(λ,θ,β)d(λ,θ,β)(λ,θ,β)e(λ,θ,β)+ρ(λ,θ,β)d(λ,θ,β), (26)

    where u(λ,θ,β) is the function of λ, θ and β only, (λ,θ,β) is the log-likelihood function and ρ(λ,θ,β)=logg(λ,θ,β). Hence, I(x) can be estimated as

    I(x)=u(ˆλ,ˆθ,ˆβ)+(ˆuλa1+ˆuθa2+ˆuβa3+a4+a5)+12[A(ˆuλˆσλλ+ˆuθˆσλθ+ˆuβˆσλβ)+B(ˆuλˆσθλ+ˆuθˆσθθ+ˆuβˆσθβ)+C(ˆuλˆσβλ+ˆuθˆσβθ+ˆuβˆσββ)], (27)

    where ˆλ, ˆθ,and ˆβ are the MLEs of λ, θ and β respectively, and

    ai=ˆρλˆσiλ+ˆρθˆσiθ+ˆρβˆσiβ for i=1,2,3,a4=ˆuλθˆσλθ+ˆuλβˆσλβ+ˆuθβˆσθβ ,          a5=12(ˆuλλˆσλλ+ˆuθθˆσθθ+ˆuββˆσββ),      }, (28)
    A=ˆσλλˆλλλ+2ˆσλθˆλθλ+2ˆσλβˆλβλ+2ˆσθβˆθβλ+ˆσθθˆθθλ+ˆσββˆββλB=ˆσλλˆλλθ+2ˆσλθˆλθθ+2ˆσλβˆλβθ+2ˆσθβˆθβθ+ˆσθθˆθθθ+ˆσββˆββθ  C=ˆσλλˆλλβ+2ˆσλθˆλθβ+2ˆσλβˆλββ+2ˆσθβˆθββ+ˆσθθˆθθβ+ˆσββˆβββ }, (29)
    ρi=ρϕi,ui=u(ϕ1,ϕ2,ϕ3)ϕi,uij=u(ϕ1,ϕ2,ϕ3)ϕiϕj,ij=(ϕ1,ϕ2,ϕ3)ϕiϕj,ijl=(ϕ1,ϕ2,ϕ3)ϕiϕjϕl, (30)

    where ϕ1=λ, ϕ2=θ, ϕ3=β, i,j,l=1,2,3 and σij is (i,j)th elements of the matrix ˆI1(ˆλ,ˆθ,ˆβ) in Subsection (2.1).

    ˆσij={1ˆIij,i=j0,ij,for i,j=1,2,3.

    From (20), ρ(λ,θ,β) can be written as

    ρ(λ,θ,β)=logg(λ,θ,β)=(a11)logλ+(a21)logθ+(a31)logβ(b1λ+b2θ+b3β).

    Hence, we get

    ˆρλ=logg(λ,θ,β)λ|(λ,θ,β)=(ˆλ,ˆθ,ˆβ)=(a11)ˆλb1ˆρθ=logg(λ,θ,β)θ|(λ,θ,β)=(ˆλ,ˆθ,ˆβ)=(a21)ˆθb2ˆρβ=logg(λ,θ,β)β|(λ,θ,β)=(ˆλ,ˆθ,ˆβ)=(a31)ˆβb3}. (31)

    Now, under SE and LINEX loss function, the approximate Bayes estimators λ are given by

    ˆλBS=E(λ|x_)=ˆλ+(ˆuλa1+ˆuθa2+ˆuβa3+a4+a5)+12[A(ˆuλˆσλλ+ˆuθˆσλθ+ˆuβˆσλβ)+B(ˆuλˆσθλ+ˆuθˆσθθ+ˆuβˆσθβ)+C(ˆuλˆσβλ+ˆuθˆσβθ+ˆuβˆσββ)], (32)

    and

    ˆλBL=1clogE(ecλ|x_)=1clog{ecˆλ+(ˆuλa1+ˆuθa2+ˆuβa3+a4+a5)+12[A(ˆuλˆσλλ+ˆuθˆσλθ+ˆuβˆσλβ)+B(ˆuλˆσθλ+ˆuθˆσθθ+ˆuβˆσθβ)+C(ˆuλˆσβλ+ˆuθˆσβθ+ˆuβˆσββ)]}. (33)

    Similarly, we can obtain the approximate Bayes estimators of θ, β, r(t) and h(t).

    It is known that there are several procedures of MCMC technique available in which samples are generated from the conditional posterior densities. One of the simplest MCMC procedures is the Gibbs sampling procedure. Another procedure is considered the Metropolis-Hastings (M-H) algorithm and later extended by Hastings [23]. A more general procedure of MCMC procedures which we will use is considered the M-H within Gibbs sampling. Gibbs sampler is required to decompose the joint posterior distribution into full conditional distributions for each parameter and then sample from them. Since, g2(θ|λ,β,x_) follows gamma [n+a2,b2+log(1+λxβL(n))], it is quite simple to generate samples of θ from g2(θ|λ,β,x_) by implementing any gamma generating routine. In addition, the conditional posterior densities g1(λ|θ,β,x_) and g3(β|θ,λ,x_) can't be reduced analytically to well-known distributions. So, according to Tierney [24] M-H algorithm within Gibbs sampling with normal proposal distribution is used to conduct the MCMC methodology. The hybrid M-H algorithm and Gibbs sampler works as follows:

    (1) Use the MLEs as the initial value, denoted by ˆλ(0),ˆθ(0)and ˆβ(0).

    (2) Set j=1.

    (3) Generate θ(j) from Gamma(n+a2,b2+log(1+λxβL(n))).

    (4) Using M-H algorithm, generate λ(j) and β(j) from g1(λ(j1)|θ(j),β(j1),x_), and g3(β(j1)|θ(j),λ(j1),x_), respectively, with the normal distributions N(λ(j1),Λ11) and N(β(j1),Λ33).

    (a) Generate λ from N(λ(j1),Λ11) and β from N(β(j1),Λ33).

    (b) Evaluate the acceptance probabilities

    Ωλ=min[1,g1(λ|θ(j),β(j1),x_)g1(λ(j1)|θ(j),β(j1),x_)],Ωβ=min[1,g3(β|θ(j),λ(j),x_)g3(β(j1)|θ(j),λ(j),,x_)],

    (c) Generate a ρ1 and ρ2 from a Uniform [0,1] distribution.

    (d) If ρ1<Ωλ accept the proposal and set λ=λ(j), else set λ(j)=λ(j1).

    (e) If ρ2<Ωβ accept the proposal and set β=β(j), else set β(j)=β(j1).

    (5) Compute r(t) and h(t) as

    r(j)(t)=1(1+λ(j)tβ(j))θ(j),t>0,
    h(j)(t)=λ(j)θ(j)β(j)t(β(j)+1)(1+λ(j)tβ(j))(θ(j)+1)1(1+λ(j)tβ(j))θ(j),  t>0.

    (6) Set j=j+1.

    (7) Repeat Steps (3)(6) M times. The first M0 simulated varieties are ignored to remove the affection of the selection of initial value and to guarantee the convergence. Then the selected sample are λ(j),θ(j),β(j),r(j)(t) and h(j)(t), j=M0+1,...,M, for sufficiently large M, forms an approximate posterior sample which can be used to develop the Bayesian inference.

    (8) Under SE and LINEX Loss functions the approximate Bayes estimate of η (where η= λ, θ, β, r(t) and h(t)) can be obtained by

    ˆηBS=1MM0Mj=M0+1η(j),ˆηBL=1clog(1MM0Mj=M0+1ecη(j)),c0,

    where M0 is the burn-in period and η(j)= λ(j), θ(j), β(j), r(j)(t) and h(j)(t).

    (9) To compute the CRIs of η, order {ηM0+1,ηM0+2,,ηM} as {η[1],η[2],,η[MM0]}. Then, the (1γ)100% CRIs of η can be given by

    [η(MM0)(γ2), η(MM0)(1γ2)].

    In this section, through a simulation study using R language, we compare the performances of the MLEs and Bayes estimates of the unknown parameters as well as reliability and hazard rate functions of the DD under SE and LINEX loss functions proposed in the previous sections. Here, we consider the simulation results in the case of (λ,θ,β)=(3.2,5.9,1.75). Then, the true values of r(t) and h(t) at time t=1.3 are computed to be 0,963895 and 0,234171, respectively. The performance of estimators is evaluated in terms of mean squared error (MSE) which computed as

    MSE=1NNk=1(ˆϖ(k)iϖi)2,N=1000,i=1,2,...,5,

    where ϖ1=λ, ϖ2=θ, ϖ3=β, ϖ4=r(t) and ϖ5=h(t) for the point estimates (MLE, Lindley, MCMC), also average lengths (ALs) and coverage probability (CPs), which computed as the number of CIs that covered the true values divided by 1000, for interval estimates (ACIs and CRIs). Bayes estimates and the highest posterior density CRIs are computed based on M=12000 MCMC samples and discard the first values M0=2000 as burn-in. In addition, we assume the informative gamma priors for λ, θ and β that is, when the hyper-parameters are ai=1.3 and bi=2.2, i=1,2,3. Moreover, 95% CRIs were computed for each simulated sample. In our study, we consider different sample sizes n=5, 7, 9, 12, 15, 18, 20, 23 and 25. The gamma distribution was chosen as prior distributions of the parameters because it is the most appropriate one that matches the maximum likelihood function. Moreover, they are from the same family. The evidence for this is that two of the full conditional posterior distributions of the parameters resulted in a gamma distribution, which proves the validity of the chose. In addition, choosing another prior distribution or dependent prior will increase the complexity and difficulty of mathematical equations. Gamma distribution is one of the rich distributions, as when changing its parameters (hyper-parameters), we get new data with new information, so it is the focus of attention of most statisticians. Gamma distribution is flexible model with any values of hyper-parameters. A special case: When all hyper-parameters of gamma distribution are zero, we obtain Jaffrey prior. The results of means, MSEs, ALs and CPs of estimates are displayed in Tables 15.

    Table 1.  Means estimates of λ (first row) with their MSEs (second row) and ALs (first row) for 95% ACIs and CRIs of λ with their CPs (second row).
    n MLE Lindley MCMC
    ALs SE LINEX SE LINEX ALs
    c=2 c=2 c=2 c=2
    5 3.3114 4.2564 3.2687 3.1475 3.2456 3.1457 3.0147 3.2001 3.9452
    0.0754 0.938 0.0715 0.0726 0.0699 0.0684 0.0689 0.0655 0.941
    7 3.1547 4.0452 3.3124 3.1479 2.9988 2.9945 3.0145 2.9997 3.7658
    0.0721 0.939 0.0694 0.0696 0.0673 0.0668 0.0670 0.0622 0.947
    9 3.2568 3.8779 2.9745 3.0145 2.9968 3.1245 3.0897 3.0368 3.6877
    0.0667 0.939 0.0645 0.0639 0.0617 0.0610 0.0599 0.0574 0.952
    12 3.3110 3.5674 3.0124 3.1478 3.2010 2.9997 3.1457 3.1278 3.3947
    0.0635 0.945 0.0628 0.0629 0.0609 0.0593 0.0595 0.0551 0.958
    15 3.4789 3.2458 3.2478 2.9965 3.0124 3.1457 3.1211 3.2001 3.1547
    0.0576 0.951 0.0556 0.0555 0.0539 0.0528 0.0521 0.0503 0.962
    18 2.9999 2.9847 3.1245 2.9984 3.2658 3.2547 3.3142 3.2874 2.8845
    0.0514 0.949 0.0497 0.0492 0.0465 0.0457 0.0449 0.0425 0.959
    20 3.1978 2.8759 3.1847 3.0988 3.1598 2.9987 3.0567 3.0696 2.5546
    0.0425 0.952 0.0419 0.0407 0.0388 0.0374 0.0368 0.0347 0.961
    23 3.2457 2.5588 3.2581 3.3146 2.9478 3.0126 3.2145 3.1897 2.2665
    0.0339 0.949 0.0305 0.0298 0.0267 0.0258 0.0249 0.0229 0.965
    25 3.1457 2.3691 3.1879 3.0457 2.9968 3.0458 3.1254 3.2045 1.9974
    0.0215 0.955 0.0197 0.0191 0.0176 0.0168 0.0157 0.0126 0.967

     | Show Table
    DownLoad: CSV
    Table 2.  Means estimates of θ (first row) with their MSEs (second row) and ALs (first row) for 95% ACIs and CRIs of θ with their CPs (second row).
    n MLE Lindley MCMC
    ALs SE LINEX SE LINEX ALs
    c=2 c=2 c=2 c=2
    5 6.1454 8.3335 6.1245 6.2014 6.1689 6.1445 6.2134 5.9987 7.4568
    0.1247 0.945 0.1186 0.1175 0.1055 0.1023 0.1014 0.0998 0.949
    7 6.2145 7.9451 6.2457 6.2479 5.8479 5.9974 6.1102 6.1058 7.1639
    0.1064 0.951 0.1032 0.1021 0.0974 0.0968 0.0961 0.0943 0.956
    9 6.1892 6.8654 6.1003 6.1450 5.9984 6.2354 6.1020 5.9989 6.3547
    0.0989 0.948 0.0975 0.0969 0.0948 0.0931 0.0926 0.0911 0.961
    12 5.9999 6.2457 6.1002 6.1109 5.9476 6.1120 6.2001 6.1036 5.8411
    0.0812 0.954 0.0806 0.0795 0.0776 0.0769 0.0758 0.0724 0.958
    15 6.1369 5.6774 6.2457 6.1582 6.1894 5.8997 6.1457 6.0236 4.9876
    0.0786 0.939 0.0775 0.0763 0.0751 0.0748 0.0739 0.0699 0.951
    18 6.0258 5.1462 6.1457 6.0458 6.1243 6.0587 6.1245 6.0039 4.4698
    0.0743 0.948 0.0725 0.0714 0.0697 0.0691 0.0687 0.0663 0.963
    20 6.1001 4.7695 6.0457 6.1212 6.0124 6.2100 6.1239 6.0358 3.9974
    0.0672 0.951 0.0653 0.0647 0.0624 0.0620 0.0618 0.0586 0.957
    23 6.0989 4.3588 6.1258 6.0847 5.9978 6.1245 6.0845 5.9778 3.7656
    0.0553 0.953 0.0534 0.0527 0.0499 0.0492 0.0487 0.0452 0.959
    25 6.0002 4.1235 6.1024 6.0647 6.0006 6.1354 6.0874 6.0110 3.4657
    0.0422 0.958 0.0413 0.0405 0.0387 0.0381 0.0375 0.0329 0.962

     | Show Table
    DownLoad: CSV
    Table 3.  Means estimates of β (first row) with their MSEs (second row) and ALs (first row) for 95% ACIs and CRIs of β with their CPs (second row).
    n MLE Lindley MCMC
    ALs SE LINEX SE LINEX ALs
    c=2 c=2 c=2 c=2
    5 1.9124 2.4314 1.8475 1.9562 1.8989 1.8974 1.9974 1.6987 2.1099
    0.0312 0.947 0.0297 0.0278 0.0256 0.0249 0.0241 0.0222 0.954
    7 2.0132 1.9947 1.9564 1.8794 1.7658 1.8122 1.9445 1.7458 1.9547
    0.0286 0.949 0.0269 0.0254 0.0236 0.0231 0.0233 0.0198 0.961
    9 1.8456 1.9648 1.7452 1.9325 1.6452 1.9745 1.6999 1.7421 1.8945
    0.0259 0.945 0.0238 0.0229 0.0197 0.0195 0.0189 0.0175 0.958
    12 1.7564 1.9044 1.6952 1.8475 1.6994 1.7569 1.8002 1.7236 1.8263
    0.0226 0.952 0.0215 0.0208 0.0179 0.0172 0.0169 0.0149 0.952
    15 1.8452 1.8362 1.8023 1.7459 1.7256 1.7320 1.8001 1.7136 1.7956
    0.0198 0.951 0.0185 0.0181 0.0165 0.0160 0.0152 0.0128 0.963
    18 1.7938 1.7639 1.6945 1.7365 1.7214 1.6998 1.7412 1.7365 1.5547
    0.0166 0.949 0.0158 0.0147 0.0125 0.0119 0.0120 0.0108 0.961
    20 1.7422 1.5345 1.7695 1.7469 1.6947 1.6879 1.7189 1.7511 1.2636
    0.0125 0.953 0.0113 0.0105 0.0099 0.0097 0.0096 0.0091 0.959
    23 1.8001 1.3458 1.7598 1.7423 1.7124 1.6978 1.8011 1.7433 1.1009
    0.0105 0.958 0.0098 0.0093 0.0085 0.0083 0.0081 0.0077 0.962
    25 1.7542 1.1258 1.8451 1.8001 1.7645 1.7423 1.8022 1.7449 1.0658
    0.0092 0.961 0.0087 0.0085 0.0081 0.0079 0.0078 0.0072 0.964

     | Show Table
    DownLoad: CSV
    Table 4.  Means estimates of r(t) (first row) with their MSEs (second row) and ALs (first row) for 95% ACIs and CRIs of r(t) with their CPs (second row).
    n MLE Lindley MCMC
    ALs SE LINEX SE LINEX ALs
    c=2 c=2 c=2 c=2
    5 0.9865 0.4325 0.9745 0.9877 0.9687 0.9547 0.96359 0.9584 0.4155
    0.0065 0.939 0.0061 0.0061 0.0057 0.0056 0.0055 0.0052 0.954
    7 0.9963 0.4136 0.9876 0.9911 0.9689 0.9584 0.9741 0.9711 0.3847
    0.0062 0.948 0.0059 0.0058 0.0056 0.0054 0.0053 0.0049 0.949
    9 0.9687 0.3965 0.9587 0.9522 0.9713 0.9573 0.9505 0.9619 0.3684
    0.0058 0.949 0.0055 0.0054 0.0051 0.0049 0.0049 0.0044 0.956
    12 0.9764 0.3577 0.9745 0.9612 0.9589 0.9645 0.9781 0.9587 0.3365
    0.0054 0.951 0.0051 0.0050 0.0048 0.0047 0.0046 0.0041 0.961
    15 0.9846 0.3248 0.9589 0.9611 0.9745 0.9639 0.9745 0.9599 0.3157
    0.0049 0.948 0.0047 0.0046 0.0044 0.0043 0.0042 0.0038 0.959
    18 0.9645 0.2984 0.9587 0.9642 0.9639 0.9587 0.9647 0.9711 0.2769
    0.0043 0.956 0.0041 0.0041 0.0038 0.0037 0.0037 0.0034 0.962
    20 0.9746 0.2569 0.9589 0.9645 0.9678 0.9784 0.9642 0.9599 0.2168
    0.0039 0.952 0.0037 0.0038 0.0036 0.0035 0.0034 0.0029 0.958
    23 0.9599 0.2268 0.9641 0.9522 0.9517 0.9625 0.9589 0.9602 0.1987
    0.0036 0.954 0.0034 0.0035 0.0033 0.0032 0.0031 0.0026 0.963
    25 0.9615 0.1899 0.9658 0.9647 0.9578 0.9678 0.9565 0.9584 0.1634
    0.0032 0.959 0.0031 0.0029 0.0027 0.0026 0.0025 0.0021 0.965

     | Show Table
    DownLoad: CSV
    Table 5.  Means estimates of h(t) (first row) with their MSEs (second row) and ALs (first row) for 95% ACIs and CRIs of h(t) with their CPs (second row).
    n MLE Lindley MCMC
    ALs SE LINEX SE LINEX ALs
    c=2 c=2 c=2 c=2
    5 0.2547 0.3145 0.2478 0.2412 0.2369 0.2456 0.2514 0.2399 0.2834
    0.0027 0.951 0.0025 0.0023 0.0021 0.0019 0.0019 0.0017 0.959
    7 0.2345 0.2965 0.2547 0.2311 0.2245 0.2199 0.2236 0.2147 0.2563
    0.0025 0.947 0.0023 0.0022 0.0019 0.0017 0.0016 0.0014 0.954
    9 0.2456 0.2463 0.2235 0.2344 0.2469 0.2399 0.2147 0.2187 0.2134
    0.0022 0.945 0.0021 0.0021 0.0017 0.0016 0.0015 0.0012 0.956
    12 0.2364 0.2184 0.2145 0.2287 0.2345 0.2289 0.2198 0.2301 0.1966
    0.0020 0.955 0.0018 0.0017 0.0015 0.0014 0.0014 0.0011 0.961
    15 0.2258 0.1875 0.2391 0.2457 0.2589 0.2469 0.2231 0.2189 0.1644
    0.0017 0.951 0.0016 0.0015 0.0013 0.0012 0.0013 0.0010 0.964
    18 0.2541 0.1563 0.2236 0.2647 0.2345 0.2239 0.2199 0.2311 0.1367
    0.0014 0.948 0.0013 0.0013 0.0011 0.0010 0.0012 0.0009 0.959
    20 0.2475 0.1359 0.2345 0.2265 0.2274 0.2189 0.2314 0.2188 0.1199
    0.0012 0.951 0.0011 0.0010 0.0009 0.0008 0.0008 0.0007 0.958
    23 0.2451 0.1187 0.2514 0.2244 0.2314 0.2174 0.2314 0.2296 0.1098
    0.0011 0.957 0.0010 0.0009 0.0008 0.0007 0.0006 0.0005 0.963
    25 0.2399 0.1059 0.2246 0.2187 0.2236 0.2378 0.2198 0.2298 0.0978
    0.0009 0.959 0.0008 0.0008 0.0007 0.0006 0.005 0.0004 0.966

     | Show Table
    DownLoad: CSV

    From the results in Tables 15, we observe the following:

    (1) It is clear that from all Tables, as sample size n increases, the MSEs and average interval lengths decrease, also Bayes estimates are performed better than the MLEs of λ, θ, β, r(t) and h(t) in terms of MSEs and average interval lengths, as expected.

    (2) Bayes estimates under MCMC technique perform better than Bayes estimates under Lindley approximation in the sense of having smaller MSEs.

    (3) Bayes estimates under Lindley approximation perform better than MLEs in the sense of having smaller MSEs.

    (4) Bayes estimates under LINEX loss function perform better than Bayes estimates under SE loss function in the sense of having smaller MSEs.

    (5) Bayes estimate under LINEX with c=2 provides better estimates for λ, θ, β, r(t) and h(t) because of having the smallest MSEs.

    (6) Bayes estimates under LINEX for the choice c=2 are performed better than their estimates for the choice c=2 in the sense of having smaller MSEs.

    (7) Both MLE and Bayesian methods have very close estimates and their ACIs have quite high CPs (around 0.95). Also, the Bayesian CRIs have the highest CPs.

    (8) In general, we can conclude that the best estimation method is the Bayes method under MCMC technique with LINEX loss function, especially if the prior information about the problem under study is available.

    To illustrate the inferential procedures developed in the preceding sections. we use the data set from Almongy [25], see Table 6. A complete sample from the data represents a COVID-19 data belonging to the Netherlands of 30 days, which recorded from 31 March to 30 April 2020. This data formed of rough mortality rate. The data are as follows

    Table 6.  The mortality rate of COVID-19.
    14.918 10.656 12.274 10.289 10.832 7.099 5.928 13.211 7.968 7.584
    5.555 6.027 4.097 3.611 4.960 7.498 6.940 5.307 5.048 2.857
    2.254 5.431 4.462 3.883 3.461 3.647 1.974 1.273 1.416 4.235

     | Show Table
    DownLoad: CSV

    For the goodness of fit test, we compute the Kolmogorov-Smirnov (K-S) distance. It was found that the K-S is 0.072956 and the associated p-value is 0.9936. Therefore, according to the p-value we can say that the DD fits quite well to the above data. Empirical, QQ and PP plots are shown in Figure 1, which clear that the DD fits the data very well. Table 7 lists the lower record values from the observed data.

    Figure 1.  Graphical fitting of the DD.
    Table 7.  The lower record values for the real data.
    XL(1) XL(2) XL(3) XL(4) XL(5) XL(6)
    14.918 10.656 10.289 7.099 5.928 5.555
    XL(7) XL(8) XL(9) XL(10) XL(11) XL(12)
    4.097 3.611 2.857 2.254 1.974 1.273

     | Show Table
    DownLoad: CSV

    Based on these 12 lower record values, the MLEs and ACIs for λ, θ, β, r(10) and h(10) are computed to be as in Tables 6 and 7. Moreover, to compute the Bayesian estimates under Lindley's approximation and MCMC technique, the prior distributions of the parameters λ, θ and β are needed to specify. Since we do not have any prior information available, we assume that the non-informative gamma priors for λ, θ and β that is, when the hyperparameters are ai=0.0001 and bi=0.0001, i=1,2,3. Under MCMC technique, the posterior analysis was done across combining M-H algorithm within Gibbs sampler. To conduct the MCMC algorithm, which was described in Subsection 3.2, the initial values for the parameters λ, θ and β were taken to be their MLEs. In addition, we run the chain for 12000 times and discard the first 2000 values as 'burn-in' to avoid the effect of the initial values (starting point). The Bayesian estimates as well as 95% CRIs for λ, θ, β, r(10) and h(10) are presented in Tables 8 and 9.

    Table 8.  Point estimations of λ, θ, β, r(t) and h(t).
    Parameter MLE Lindley MCMC
    SE LINEX SE LINEX
    c=2 c=2 c=2 c=2
    λ 229.025 223.147 224.532 219.967 218.365 217.542 214.658
    θ 2.4484 2.4786 2.5123 2.3954 2.3457 2.3369 2.2987
    β 2.3128 2.2567 2.3014 2.2111 2.2014 2.2155 2.1996
    r(t) 0.8402 0.8246 0.8257 0.8365 0.8154 0.8022 0.8001
    h(t) 0.0568 0.05432 0.0554 0.05311 0.05278 0.05266 0.0499

     | Show Table
    DownLoad: CSV
    Table 9.  95% confidence intervals of λ, θ, β, r(t) and h(t).
    Parameter ACIs CRIs
    Interval Length Interval Length
    λ (140.557,317.493) 176,937 (137.647,298.346) 160.699
    θ (2.2190, 2.6779) 0.4589 (2.1012, 2.5113) 0.4101
    β (1.6068, 3.0187) 1.4119 (1.9231, 3.1631) 1.2400
    r(t) (0.4962, 1.1842) 0.6880 (0.5359, 0.9998) 0.4639
    h(t) (0.0000, 0.1672) 0.1672 (0.0073, 0.1378) 0.1305

     | Show Table
    DownLoad: CSV

    In this paper, we have estimated the reliability and hazard rate functions of the Dagum distribution based record statistics. To get the best estimators for this purpose, the maximum likelihood and Bayes estimations methods have been utilized. The MLEs and the asymptotic confidence intervals based on the observed Fisher information matrix have been discussed. For Bayesian estimation approach, we assume the gamma priors for the unknown parameters and provide the Bayes estimators under the assumptions of SE and LINEX loss functions. Since the Bayes estimates cannot be obtained in explicit form, we have applied both Lindley's approximation and MCMC technique with Metropolis–Hasting algorithm within Gibbs sampler to compute the approximate Bayes estimates and constructed the credible intervals. The performance of the proposed techniques has been compared via a Monte Carlo simulation study with different sample sizes n. It was found that if we have no prior information on the unknown parameters, then it is always better to use the MLEs rather than the Bayes estimators, because the Bayes estimators are computationally more expensive. Finally, for illustrative purpose, the details have been explained using a real-life data set of COVID-19 mortality rates.

    The authors declare no conflict of interest.

    [1] R.A. Adams and J.J.F. Fournier, Sobolev spaces. Pure and applied mathematics, Vol. 140. 2 Ed., Elsevier/Academic Press, Amsterdam, 2003.
    [2] H.W. Alt, Linear Functional Analysis. An Application-Oriented Introduction. Translated from the German edition by Robert Nürnberg. Universitext. Springer Berlin London, 2016.
    [3] S. Bosia, M. Conti, and M. Grasselli, On the Cahn-Hilliard-Brinkman system, Commun. Math. Sci., 13 (2015), no. 6, 1541-1567.
    [4] E.A. Coddington and N. Levinson, Theory of Ordinary Di erential Equations. International series in pure and applied mathematics. Tata McGraw-Hill, New York, 1955.
    [5] P. Colli, G. Gilardi, and D. Hilhorst, On a Cahn–Hilliard type phase field model related to tumor growth. Discrete Contin. Dyn. Syst., 35 (2015), no. 6, 2423-2442.
    [6] P. Colli, G. Gilardi, E. Rocca, and J. Sprekels, Vanishing viscosities and error estimate for a Cahn-Hilliard type phase field system related to tumor growth, Nonlinear Anal. Real World Appl., 26 (2015), 93-108.
    [7] P. Colli, G. Gilardi, E. Rocca, and J. Sprekels, Asymptotic analyses and error estimates for a Cahn-Hilliard type phase field system modelling tumor growth, Discrete Contin. Dyn. Syst. Ser. S., in press (2016).
    [8] V. Cristini and J. Lowengru, Multiscaled modeling of cancer. An Integrated Experiemental and Mathematical Modeling Approach. Cambridge University Press, 2010.
    [9] M. Dai, E. Feireisl, E. Rocca, G. Schimperna, and M. Schonbek, Analysis of a di use interface model for multispecies tumor growth, Preprint arXiv:1507.07683 (2015).
    [10] F. Della Porta and M. Grasselli, On the nonlocal Cahn-Hilliard-Brinkman and Cahn-Hilliard-Hele-Shaw systems, Commun. Pure Appl. Anal., 15 (2016), 299-317.
    [11] E. DiBenedetto, Degenerate Parabolic Equations. Universitext. Springer–Verlag New York, 1993.
    [12] L.C. Evans, Partial Differential Equations. Graduate Studies in Mathematics, Volume 19. AMS,Providence, Rhode Island, 2002.
    [13] A. Fasano, A. Bertuzzi, and A. Gandolfi, Mathematical modeling of tumour growth and treatment.Complex Systems in Biomedicine. Springer Milan, 2006.
    [14] X. Feng and S.Wise, Analysis of a Darcy-Cahn-Hilliard diffuse interface model for the Hele-Shawflow and its fully discrete finite element approximation, SIAM J. Numer. Anal., 50 (2012), no. 3,1320-1343.
    [15] A. Friedman, Partial Differential Equations. Holt, Rinehart and Winston, New York, 1969.
    [16] S. Frigeri, M. Grasselli, and E. Rocca, On a diffuse interface model of tumor growth, European J.Appl. Math., 26 (2015), 215-243.
    [17] H. Garcke and K.F. Lam, Well-posedness of a Cahn-Hilliard system modelling tumour growth withchemotaxis and active transport, European J. Appl. Math., in press (2016). DOI http://dx.doi.org/10.1017/S0956792516000292.
    [18] H. Garcke and K.F. Lam, Analysis of a Cahn-Hilliard system with non zero Dirichlet conditionsmodelling tumour growth with chemotaxis, Preprint arXiv:1604.00287 (2016).
    [19] H. Garcke, K.F. Lam, and E. Rocca, Optimal control of treatment time in a diffuse interface modelfor tumor growth, Preprint arXiv:1608.00488 (2016).
    [20] H. Garcke, K.F. Lam, E. Sitka, and V. Styles, A Cahn-Hilliard-Darcy model for tumour growth withchemotaxis and active transport, Math. Models Methods Appl. Sci., 26 (2016), no. 6, 1095-1148.
    [21] A. Hawkins-Daarud, K.G. van der Zee, and J.T. Oden, Numerical simulation of a thermodynamicallyconsistent four-species tumor growth model, Int. J. Numer. Methods Biomed. Eng., 28 (2012),3-24.
    [22] J. Jiang, H.Wu, and S. Zheng, Well-posedness and long-time behavior of a non-autonomous Cahn-Hilliard-Darcy system with mass source modeling tumor growth, J. Differential Equ., 259 (2015),no. 7, 3032-3077.
    [23] H. Lee, J. Lowengrub, and J. Goodman, Modeling pinchoff and reconnection in a Hele–Shaw cell.I. The models and their calibration, Phys. Fluids, 14 (2002), no. 2, 492-513.
    [24] H. Lee, J. Lowengrub, and J. Goodman, Modeling pinchoff and reconnection in a Hele-Shaw cell.II. Analysis and simulation in the nonlinear regime, Phys. Fluids, 14 (2002), no. 2, 514-545.
    [25] J.S. Lowengrub, E. Titi, and K. Zhao, Analysis of a mixture model of tumor growth, European J.Appl. Math., 24 (2013), 691-734.
    [26] M. Renardy and R.C. Rogers, An Introduction to Partial Differential Equations. Texts in AppliedMathematics. 2 Eds., Springer-Verlag New York, 2004.
    [27] H.L. Royden and P. Fitzpatrick, Real Analysis. Featured Titles for Real Analysis Series. 4 Eds.,Pearson Prentice Hall, Boston, 2010.
    [28] J. Simon, Compact sets in space Lp(0; T; B), Ann. Mat. Pura Appl., 146 (1986), no. 1, 65-96.
    [29] X. Wang and H. Wu, Long-time behavior for the Hele-Shaw-Cahn-Hilliard system, Asymptot.Anal., 78 (2012), no. 4, 217-245.
    [30] X. Wang and Z. Zhang, Well-posedness of the Hele-Shaw-Cahn-Hilliard system, Ann. Inst. H. Poincar´e Anal. Non Lin´eaire, 30 (2013), no. 3, 367-384.
  • This article has been cited by:

    1. Rashid Jan, Normy Norfiza Abdul Razak, Sania Qureshi, Imtiaz Ahmad, Salma Bahramand, Modeling Rift Valley fever transmission: insights from fractal-fractional dynamics with the Caputo derivative, 2024, 4, 2767-8946, 163, 10.3934/mmc.2024015
    2. Abdulrahman Obaid Alshammari, Imtiaz Ahmad, Rashid Jan, Sahar Ahmed Idris, Fractional-calculus analysis of the dynamics of CD4+
    T cells and human immunodeficiency viruses, 2024, 1951-6355, 10.1140/epjs/s11734-024-01192-5
    3. Zakieh Avazzadeh, Hossein Hassani, Mohammad Javad Ebadi, Ali Bayati Eshkaftaki, Roghayeh Katani, Alireza Rezvani, Generalization of Bernoulli polynomials to find optimal solution of fractional hematopoietic stem cells model, 2024, 99, 0031-8949, 085015, 10.1088/1402-4896/ad5e44
    4. Sehrish Ramzan, Syeda Alishwa Zanib, Muzamil Abbas Shah, Nadeem Abbas, Wasfi Shatanawi, On bifurcation analysis and numerical investigation of alcohol consumption dynamics, 2024, 0217-9849, 10.1142/S021798492550040X
    5. Alamgir Khan, Ikram Ullah, Javed Iqbal, Kamal Shah, Muhammad Bilal, An innovative method for solving the nonlinear fractional diffusion reaction equation with quadratic nonlinearity analysis, 2025, 100, 0031-8949, 015209, 10.1088/1402-4896/ad952b
    6. Imtiaz Ahmad, Muhammad Nawaz Khan, Rashid Jan, Normy Norfiza Abdul Razak, Efficient numerical method for pricing multi-asset options with the time-fractional Black-Scholes model: focus on American and digital options, 2025, 5, 2767-8946, 147, 10.3934/mmc.2025011
    7. S. Hajiollow, F. Zabihi, Using the spectral meshless radial basis functions method for solving time fractional Burgers’ equation, 2025, 2025, 1687-2770, 10.1186/s13661-025-02075-x
  • Reader Comments
  • © 2016 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(7744) PDF downloads(1367) Cited by(47)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog