Loading [MathJax]/jax/output/SVG/jax.js
Research article

A difference scheme for a triangular system of conservation laws with discontinuous flux modeling three-phase flows

  • Received: 29 August 2022 Revised: 05 November 2022 Accepted: 10 November 2022 Published: 28 November 2022
  • A triangular system of conservation laws with discontinuous flux that models the one-dimensional flow of two disperse phases through a continuous one is formulated. The triangularity arises from the distinction between a primary and a secondary disperse phase, where the movement of the primary disperse phase does not depend on the local volume fraction of the secondary one. A particular application is the movement of aggregate bubbles and solid particles in flotation columns under feed and discharge operations. This model is formulated under the assumption of a variable cross-sectional area. A monotone numerical scheme to approximate solutions to this model is presented. The scheme is supported by three partial theoretical arguments. Firstly, it is proved that it satisfies an invariant-region property, i.e., the approximate volume fractions of the three phases, and their sum, stay between zero and one. Secondly, under the assumption of flow in a column with constant cross-sectional area it is shown that the scheme for the primary disperse phase converges to a suitably defined entropy solution. Thirdly, under the additional assumption of absence of flux discontinuities it is further demonstrated, by invoking arguments of compensated compactness, that the scheme for the secondary disperse phase converges to a weak solution of the corresponding conservation law. Numerical examples along with estimations of numerical error and convergence rates are presented for counter-current and co-current flows of the two disperse phases.

    Citation: Raimund Bürger, Stefan Diehl, M. Carmen Martí, Yolanda Vásquez. A difference scheme for a triangular system of conservation laws with discontinuous flux modeling three-phase flows[J]. Networks and Heterogeneous Media, 2023, 18(1): 140-190. doi: 10.3934/nhm.2023006

    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
  • A triangular system of conservation laws with discontinuous flux that models the one-dimensional flow of two disperse phases through a continuous one is formulated. The triangularity arises from the distinction between a primary and a secondary disperse phase, where the movement of the primary disperse phase does not depend on the local volume fraction of the secondary one. A particular application is the movement of aggregate bubbles and solid particles in flotation columns under feed and discharge operations. This model is formulated under the assumption of a variable cross-sectional area. A monotone numerical scheme to approximate solutions to this model is presented. The scheme is supported by three partial theoretical arguments. Firstly, it is proved that it satisfies an invariant-region property, i.e., the approximate volume fractions of the three phases, and their sum, stay between zero and one. Secondly, under the assumption of flow in a column with constant cross-sectional area it is shown that the scheme for the primary disperse phase converges to a suitably defined entropy solution. Thirdly, under the additional assumption of absence of flux discontinuities it is further demonstrated, by invoking arguments of compensated compactness, that the scheme for the secondary disperse phase converges to a weak solution of the corresponding conservation law. Numerical examples along with estimations of numerical error and convergence rates are presented for counter-current and co-current flows of the two disperse phases.



    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] Adimurthi, G. D. Veerappa Gowda, S. Mishra, Optimal entropy solutions for conservation laws with discontinuous flux, J. Hyperbolic Differ. Equ., 2 (2005), 1–56. https://doi.org/10.1142/S0219891605000361 doi: 10.1142/S0219891605000361
    [2] B. Andreianov, K. H. Karlsen, N. H. Risebro, A theory of L1-dissipative solvers for scalar conservation laws with discontinuous flux, Arch. Ration. Mech. Anal., 201, (2011), 1–60. https://doi.org/10.1007/s00205-010-0389-4 doi: 10.1007/s00205-010-0389-4
    [3] B. Andreianov, C. Donadello, S. S. Ghoshal, U. Razafison, On the attainable set for a class of triangular systems of conservation laws, J. Evol. Equ., 15 (2015), 503–532. https://doi.org/10.1007/s00028-014-0267-x doi: 10.1007/s00028-014-0267-x
    [4] B. Andreianov, D. Mitrović, Entropy conditions for scalar conservation laws with discontinuous flux revisited, Ann. Inst. H. Poincaré C Anal. Non Linéaire, 32 (2015), 1307–1335. https://doi.org/10.1016/j.anihpc.2014.08.002 doi: 10.1016/j.anihpc.2014.08.002
    [5] A. Bressan, G. Guerra, W. Shen, Vanishing viscosity solutions for conservation laws with regulated flux, J. Differ. Equ., 266 (2018), 312–351. https://doi.org/10.1016/j.jde.2018.07.044 doi: 10.1016/j.jde.2018.07.044
    [6] R. Bürger, A. García, K. H. Karlsen, J.D. Towers, A family of numerical schemes for kinematic flows with discontinuous flux, J. Eng. Math., 60 (2008), 387–425. https://doi.org/10.1007/s10665-007-9148-4 doi: 10.1007/s10665-007-9148-4
    [7] S. Boscarino, R. Bürger, P. Mulet, G. Russo, L. M. Villada, Linearly implicit IMEX Runge-Kutta methods for a class of degenerate convection-diffusion problems, SIAM J. Sci. Comput., 37 (2015), B305–B331. https://doi.org/10.1137/140967544 doi: 10.1137/140967544
    [8] R. Bürger, S. Diehl, M. C. Martí, A conservation law with multiply discontinuous flux modelling a flotation column, Netw. Heterog. Media, 13 (2018), 339–371. https://doi.org/10.3934/nhm.2018015 doi: 10.3934/nhm.2018015
    [9] R. Bürger, S. Diehl, M. C. Martí, A system of conservation laws with discontinuous flux modelling flotation with sedimentation, IMA J. Appl. Math., 84 (2019), 930–973. https://doi.org/10.1093/imamat/hxz021 doi: 10.1093/imamat/hxz021
    [10] R. Bürger, S. Diehl, M. C. Martí, Y. Vásquez, Flotation with sedimentation: Steady states and numerical simulation of transient operation, Minerals Eng., 157 (2020), 106419. https://doi.org/10.1016/j.mineng.2020.106419 doi: 10.1016/j.mineng.2020.106419
    [11] R. Bürger, S. Diehl, M.C. Martí, Y. Vásquez, Simulation and control of dissolved air flotation and column froth flotation with simultaneous sedimentation, Water Sci. Tech., 81 (2020), 1723–1732. https://doi.org/10.2166/wst.2020.258 doi: 10.2166/wst.2020.258
    [12] R. Bürger, S. Diehl, M. C. Martí, Y. Vásquez, A degenerating convection-diffusion system modelling froth flotation with drainage, arXiv: 2202.04679, [Preprint], (2022), [cited 2022 Nov 22 ]. Available from: https://doi.org/10.48550/arXiv.2202.04679
    [13] R. Bürger, K. H. Karlsen, N. H. Risebro, J. D. Towers, Well-posedness in BVt and convergence of a difference scheme for continuous sedimentation in ideal clarifier-thickener units, Numer. Math., 97 (2004), 25–65. https://doi.org/10.1007/s00211-003-0503-8 doi: 10.1007/s00211-003-0503-8
    [14] R. Bürger, K. H. Karlsen, H. Torres, J. D. Towers, Second-order schemes for conservation laws with discontinuous flux modelling clarifier-thickener units, Numer. Math., 116 (2010), 579–617. https://doi.org/10.1007/s00211-010-0325-4 doi: 10.1007/s00211-010-0325-4
    [15] R. Bürger, K. H. Karlsen, J. D. Towers, An Engquist-Osher-type scheme for conservation laws with discontinuous flux adapted to flux connections, SIAM J. Numer. Anal., 47 (2009), 1684–1712. https://doi.org/10.1007/s00211-010-0325-4 doi: 10.1007/s00211-010-0325-4
    [16] R. Bürger, K. H. Karlsen, J. D. Towers, On some difference schemes and entropy conditions for a class of multi-species kinematic flow models with discontinuous flux, Netw. Heterog. Media, 5 (2010), 461–485. https://doi.org/10.3934/nhm.2010.5.461 doi: 10.3934/nhm.2010.5.461
    [17] G. M. Coclite, K. H. Karlsen, S. Mishra, N. H. Risebro, Convergence of vanishing viscosity approximations of 2×2 triangular systems of multi-dimensional conservation laws, Boll. Unione Mat. Ital., 2 (2009), 275–284.
    [18] G. M. Coclite, S. Mishra, N. H. Risebro, Convergence of an Engquist-Osher scheme for a multi-dimensional triangular system of conservation laws, Math. Comp., 79 (2010), 71–94. https://doi.org/10.1090/S0025-5718-09-02251-0 doi: 10.1090/S0025-5718-09-02251-0
    [19] M. G. Crandall, L. Tartar, Some relations between nonexpansive and order preserving mappings, Proc. Amer. Math. Soc., 78 (1980), 385–390. https://doi.org/10.1090/S0002-9939-1980-0553381-X doi: 10.1090/S0002-9939-1980-0553381-X
    [20] V. G. Danilov, D. Mitrovic, Delta shock wave formation in the case of triangular hyperbolic system of conservation laws. J. Differ. Equ., 245 (2008), 3704–3734. https://doi.org/10.1090/S0025-5718-09-02251-0 doi: 10.1090/S0025-5718-09-02251-0
    [21] J. E. Dickinson, K. P. Galvin, Fluidized bed desliming in fine particle flotation, Part Ⅰ, Chem. Eng. Sci., 108 (2014), 283–298. https://doi.org/10.1016/j.ces.2013.11.006 doi: 10.1016/j.ces.2013.11.006
    [22] S. Diehl, On scalar conservation laws with point source and discontinuous flux function, SIAM J. Math. Anal., 26 (1995), 1425–1451. https://doi.org/10.1137/S0036141093242533 doi: 10.1137/S0036141093242533
    [23] S. Diehl, A conservation law with point source and discontinuous flux function modelling continuous sedimentation, SIAM J. Appl. Math., 56 (1996), 388–419. https://doi.org/10.1137/S0036139994242425 doi: 10.1137/S0036139994242425
    [24] S. Diehl, Operating charts for continuous sedimentation Ⅰ: control of steady states, J. Eng. Math., 41 (2001), 117–144. https://doi.org/10.1023/A:1011959425670 doi: 10.1023/A:1011959425670
    [25] S. Diehl, The solids-flux theory—confirmation and extension by using partial differential equations, Water Res. 42 (2008), 4976–4988. https://doi.org/10.1016/j.watres.2008.09.005 doi: 10.1016/j.watres.2008.09.005
    [26] B. Engquist, S. Osher, One-sided difference approximations for nonlinear conservation laws, Math. Comp., 36 (1981), 321–351. https://doi.org/10.1090/S0025-5718-1981-0606500-X doi: 10.1090/S0025-5718-1981-0606500-X
    [27] J. A. Finch, G. S. Dobby, Column Flotation, London: Pergamon Press, 1990.
    [28] K. P. Galvin, J. E. Dickinson, Fluidized bed desliming in fine particle flotation Part Ⅱ: Flotation of a model feed, Chem. Eng. Sci., 108 (2014), 299–309. https://doi.org/10.1016/j.ces.2013.11.027 doi: 10.1016/j.ces.2013.11.027
    [29] T. Gimse, N. H. Risebro, Riemann problems with a discontinuous flux function. Third International Conference on Hyperbolic Problems, Theory, Numerical Methods and Applications, 1 (1991), 488–502.
    [30] H. Holden, N. H. Risebro, Front Tracking for Hyperbolic Conservation Laws, 2nd edition, Berlin: Springer Verlag, 2015.
    [31] E. L. Isaacson, J. B. Temple, Analysis of a singular hyperbolic system of conservation laws, J. Differ. Equ., 65 (1986), 250–268.
    [32] K. H. Karlsen, S. Mishra, N. H. Risebro, Semi-Godunov schemes for general triangular systems of conservation laws, J. Eng. Math., 60 (2008), 337–349. https://doi.org/10.1007/s10665-007-9163-5 doi: 10.1007/s10665-007-9163-5
    [33] K. H. Karlsen, S. Mishra, N. H. Risebro, Convergence of finite volume schemes for triangular systems of conservation laws, Numer. Math., 111 (2009), 559–589. https://doi.org/10.1007/s00211-008-0199-x doi: 10.1007/s00211-008-0199-x
    [34] K. H. Karlsen, N. H. Risebro, Convergence of finite difference schemes for viscous and inviscid conservation laws with rough coefficients, ESAIM: Math. Model. Numer. Anal., 35 (2001), 239–269. https://doi.org/10.1051/m2an:2001114 doi: 10.1051/m2an:2001114
    [35] K. H. Karlsen, N. H. Risebro, J. D. Towers, L1 stability for entropy solutions of nonlinear degenerate parabolic convection-diffusion equations with discontinuous coefficients, Skr. K. Nor. Vidensk. Selsk., 3 (2003), 1–49.
    [36] K. H. Karlsen, J. D. Towers, Convergence of a Godunov scheme for conservation laws with a discontinuous flux lacking the crossing condition, J. Hyperbolic Differ. Equ., 14 (2017), 671–701. https://doi.org/10.1142/S0219891617500229 doi: 10.1142/S0219891617500229
    [37] S. N. Kružkov, First order quasi-linear equations in several independent variables, Math. USSR Sb., 10 (1970), 217–243.
    [38] G. J. Kynch, A theory of sedimentation, Trans. Faraday Soc., 48 (1952), 166–176. https://doi.org/10.1039/tf9524800166 doi: 10.1039/tf9524800166
    [39] D. Mitrovic, New entropy conditions for scalar conservation laws with discontinuous flux, Discrete Contin. Dyn. Syst., 30 (2011), 1191–1210. https://doi.org/10.3934/dcds.2011.30.1191 doi: 10.3934/dcds.2011.30.1191
    [40] D. Mitrovic, V. Bojkovic, V. G. Danilov, Linearization of the Riemann problem for a triangular system of conservation laws and delta shock wave formation process, Math. Meth. Appl. Sci., 33 (2010), 904–921. https://doi.org/10.1002/mma.1226 doi: 10.1002/mma.1226
    [41] R. Pal, J. H. Masliyah, Flow characterization of a flotation column, Canad. J. Chem. Eng., 67 (1989), 916–923. https://doi.org/10.1002/cjce.5450670608 doi: 10.1002/cjce.5450670608
    [42] R. Pal, J. H. Masliyah, Oil recovery from oil in water emulsions using a flotation column, Canad. J. Chem. Eng., 68 (1990), 959–967. https://doi.org/10.1002/cjce.5450680611 doi: 10.1002/cjce.5450680611
    [43] E. Yu. Panov, Existence and strong pre-compactness properties for entropy solutions of a first-order quasilinear equation with discontinuous flux, Arch. Rational Mech. Anal., 195 (2010), 643–673. https://doi.org/10.1007/s00205-009-0217-x doi: 10.1007/s00205-009-0217-x
    [44] E. Yu. Panov, Erratum to: Existence and strong pre-compactness properties for entropy solutions of a first-order quasilinear equation with discontinuous flux, Arch. Rational Mech. Anal., 196 (2010), 1077–1078. https://doi.org/10.1007/s00205-009-0217-x doi: 10.1007/s00205-009-0217-x
    [45] P. Quintanilla, S. J. Neethling, P. R. Brito-Parada, Modelling for froth flotation control: A review, Min. Eng., 162 (2021), 106718. https://doi.org/10.1016/j.mineng.2020.106718 doi: 10.1016/j.mineng.2020.106718
    [46] J. F. Richardson, W. N. Zaki, Sedimentation and fluidisation: Part Ⅰ, Trans. Instn. Chem. Engrs. (London), 32 (1954), 35–53.
    [47] M. D. Rosini, Systems of conservation laws with discontinuous fluxes and applications to traffic, Ann. Univ. Mariae Curie-Skłodowska Sect. A, 73 (2019), 135–173.
    [48] C. O. R. Sarrico, A distributional product approach to δ-shock wave solutions for a generalized pressureless gas dynamic system, Int. J. Math., 25, (2014), 145007.
    [49] W. Shen, On the Cauchy problems for polymer flooding with gravitation, J. Differ. Equ., 261 (2016), 627–653. https://doi.org/10.1016/j.jde.2016.03.020 doi: 10.1016/j.jde.2016.03.020
    [50] P. Stevenson, P. S. Fennell, K. P. Galvin, On the drift-flux analysis of flotation and foam fractionation processes, Canad. J. Chem. Eng., 86 (2008), 635–642. https://doi.org/10.1002/cjce.20076 doi: 10.1002/cjce.20076
    [51] J. Vandenberghe, J. Chung, Z. Xu, J. Masliyah, Drift flux modelling for a two-phase system in a flotation column, Canad. J. Chem. Eng., 83 (2005), 169–176.
    [52] Y. Vásquez, Conservation Laws with Discontinuous Flux Modeling Flotation Columns, (Spanish), Doctoral Thesis of Universidad de Concepción, Concepción, 2022.
  • 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
  • © 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(2487) PDF downloads(120) Cited by(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog