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

Heat transfer analysis: convective-radiative moving exponential porous fins with internal heat generation


  • The efficiency, temperature distribution, and temperature at the tip of straight rectangular, growing and decaying moving exponential fins are investigated in this article. The influence of internal heat generation, surface and surrounding temperatures, convection-conduction, Peclet number and radiation-conduction is studied numerically on the efficiency, temperature profile, and temperature at the tip of the fin. Differential transform method is used to investigate the problem. It is revealed that thermal and thermo-geometric characteristics have a significant impact on the performance, temperature distribution, and temperature of the fin's tip.The results show that the temperature distribution of decaying exponential and rectangular fins is approximately 15 and 7% higher than growing exponential and rectangular fins respectively. It is estimated that the temperature distribution of the fin increases by approximately 6% when the porosity parameter is increased from 0.1 to 0.5. It is also observed that the decay exponential fin has better efficiency compared to growing exponential fin which offers significant advantages in mechanical engineering.

    Citation: Zia Ud Din, Amir Ali, Zareen A. Khan, Gul Zaman. Heat transfer analysis: convective-radiative moving exponential porous fins with internal heat generation[J]. Mathematical Biosciences and Engineering, 2022, 19(11): 11491-11511. doi: 10.3934/mbe.2022535

    Related Papers:

    [1] Minlong Lin, Ke Tang . Selective further learning of hybrid ensemble for class imbalanced increment learning. Big Data and Information Analytics, 2017, 2(1): 1-21. doi: 10.3934/bdia.2017005
    [2] Subrata Dasgupta . Disentangling data, information and knowledge. Big Data and Information Analytics, 2016, 1(4): 377-390. doi: 10.3934/bdia.2016016
    [3] Qinglei Zhang, Wenying Feng . Detecting Coalition Attacks in Online Advertising: A hybrid data mining approach. Big Data and Information Analytics, 2016, 1(2): 227-245. doi: 10.3934/bdia.2016006
    [4] Tieliang Gong, Qian Zhao, Deyu Meng, Zongben Xu . Why Curriculum Learning & Self-paced Learning Work in Big/Noisy Data: A Theoretical Perspective. Big Data and Information Analytics, 2016, 1(1): 111-127. doi: 10.3934/bdia.2016.1.111
    [5] Xin Yun, Myung Hwan Chun . The impact of personalized recommendation on purchase intention under the background of big data. Big Data and Information Analytics, 2024, 8(0): 80-108. doi: 10.3934/bdia.2024005
    [6] Pankaj Sharma, David Baglee, Jaime Campos, Erkki Jantunen . Big data collection and analysis for manufacturing organisations. Big Data and Information Analytics, 2017, 2(2): 127-139. doi: 10.3934/bdia.2017002
    [7] Zhen Mei . Manifold Data Mining Helps Businesses Grow More Effectively. Big Data and Information Analytics, 2016, 1(2): 275-276. doi: 10.3934/bdia.2016009
    [8] Ricky Fok, Agnieszka Lasek, Jiye Li, Aijun An . Modeling daily guest count prediction. Big Data and Information Analytics, 2016, 1(4): 299-308. doi: 10.3934/bdia.2016012
    [9] M Supriya, AJ Deepa . Machine learning approach on healthcare big data: a review. Big Data and Information Analytics, 2020, 5(1): 58-75. doi: 10.3934/bdia.2020005
    [10] Sunmoo Yoon, Maria Patrao, Debbie Schauer, Jose Gutierrez . Prediction Models for Burden of Caregivers Applying Data Mining Techniques. Big Data and Information Analytics, 2017, 2(3): 209-217. doi: 10.3934/bdia.2017014
  • The efficiency, temperature distribution, and temperature at the tip of straight rectangular, growing and decaying moving exponential fins are investigated in this article. The influence of internal heat generation, surface and surrounding temperatures, convection-conduction, Peclet number and radiation-conduction is studied numerically on the efficiency, temperature profile, and temperature at the tip of the fin. Differential transform method is used to investigate the problem. It is revealed that thermal and thermo-geometric characteristics have a significant impact on the performance, temperature distribution, and temperature of the fin's tip.The results show that the temperature distribution of decaying exponential and rectangular fins is approximately 15 and 7% higher than growing exponential and rectangular fins respectively. It is estimated that the temperature distribution of the fin increases by approximately 6% when the porosity parameter is increased from 0.1 to 0.5. It is also observed that the decay exponential fin has better efficiency compared to growing exponential fin which offers significant advantages in mechanical engineering.



    For a continuous risk outcome 0<y<1, a model with a random effect has potentially a wide application in portfolio risk management, especially, for stress testing [1,2,7,16,19], capital allocation, conditional expected shortfall estimation [3,11,17].

    Given fixed effects x=(x1,x2,,xk), two widely used regression models to estimate the expected value E(y|x) are: the fraction response model [10] and Beta regression model [4,6,8]. There are cases, however, where tail behaviours or severity levels of the risk outcome are relevant. In those cases, a regression model may no longer fit in for the requirements. In addition, a fraction response model of the form E(y|x)=Φ(a0+a1x1++akxk) may not be adequate when data exhibits significant heteroscedasticity, where Φ is a map from R1 to the open interval (0,1).

    In this paper, we assume that the risk outcome y is driven by a model:

    y=Φ(a0+a1x1++akxk+bs), (1.1)

    where s is a random continuous variable following a known distribution, independent of fixed effects (x1,x2,,xk). Parameters a0,a1,,ak are constant, while parameter b can be chosen to be dependent on (x1,x2,,xk) when required, for example, for addressing data heteroscedasticity.

    Given random effect model (1.1), the expected value E(y|x) can be deduced accordingly. It is given by the integral ΩΦ(a0+a1x1++akxk+bs)f(s)ds over the domain Ω of s, where f is the probability density of s. Given the routine QUAD implemented in SAS and Python, this integral can be evaluated as quickly as other function calls. Relative error tolerance for QUAD is 1.49e-8 in Python and is 1e-7 in SAS. But one can rescale the default tolerance to a desired level when necessary. This leads to an alternative regression tool to the fraction response model and Beta regression model.

    We introduce a family of interval distributions based on variable transformations. Probability densities for these distributions are provided (Proposition 2.1). Parameters of model (1.1) can then be estimated by maximum likelihood approaches assuming an interval distribution. In some cases, these parameters get an analytical solution without the needs for a model fitting (Proposition 4.1). We call a model with a random effect, where parameters are estimated by maximum likelihood assuming an interval distribution, an interval distribution model.

    In its simplest form, the interval distribution model y=Φ(a+bs), where a and b, are constant, can be used to model the loss rate as a random distribution for a homogeneous portfolio. Let yα and sα denote the α -quantiles for y and s at level α, 0<α<1. Then yα=Φ(a+bsα). The conditional expected shortfall for loss rate y, at level α, can then be estimated as the integral 11α[sα,+)Φ(a+bs)f(s)ds, where f is the density of s. Meanwhile, a stress testing loss estimate, derived from a model on a specific scenario, can be compared in loss rate to severity yα(=Φ(a+bsα)), to position its level of severity. A loss estimate may not have reached the desired, for example, 99% level yet, if it is far below y0.99, and far below the maximum historical loss rate. In which case, further recalibrations for the model may be required.

    The paper is organized as follows: in section 2, we introduce a family of interval distributions. A measure for tail fatness is defined. In section 3, we show examples of interval distributions and investigate their tail behaviours. We propose in section 4 an algorithm for estimating the parameters in model (1.1).

    Interval distributions introduced in this section are defined for a risk outcome over a finite open interval (c0,c1), where c0< c1 are finite numbers. These interval distributions can potentially be used for modeling a risk outcome over an arbitrary finite interval, including interval (0, 1), by maximum likelihood approaches.

    Let D=(d0,d1), d0<d1, be an open interval, where d0 can be finite or and d1 can be finite or +.

    Let

    Φ:D(c0,c1) (2.1)

    be a transformation with continuous and positive derivatives Φ(x)=ϕ(x). A special example is (c0,c1)=(0,1), and Φ:D(0,1) is the cumulative distribution function (CDF) of a random variable with a continuous and positive density.

    Given a continuous random variable s, let f and F be respectively its density and CDF. For constants a and b>0, let

    y=Φ(a+bs), (2.2)

    where we assume that the range of variable (a+bs) is in the domain D of Φ. Let g(y,a,b) and G(y,a,b) denote respectively the density and CDF of y in (2.2).

    Proposition 2.1. Given Φ1(y), functions g(y,a,b) and G(y,a,b) are given as:

    g(y,a,b)=U1/(bU2) (2.3)
    G(y,a,b)=F[Φ1(y)ab]. (2.4)

    where

    U1=f{[Φ1(y)a]/b},U2=ϕ[Φ1(y)] (2.5)

    Proof. A proof for the case when (c0,c1)=(0,1) can be found in [18]. The proof here is similar. Since G(y,a,b) is the CDF of y, it follows:

    G(y,a,b)=P[Φ(a+bs)y]
    =P{s[Φ1(y)a]/b}
    =F{[Φ1(y)a]/b}.

    By chain rule and the relationship Φ[Φ1(y)]=y, the derivative of Φ1(y) with respect to y is

    Φ1(y)y=1ϕ[Φ1(y)]. (2.6)

    Taking the derivative of G(y,a,b) with respect to y, we have

    G(y,a,b)y=f{[Φ1(y)a]/b}bϕ[Φ1(y)]=U1bU2.

    One can explore into these interval distributions for their shapes, including skewness and modality. For stress testing purposes, we are more interested in tail risk behaviours for these distributions.

    Recall that, for a variable X over (− ,+), we say that the distribution of X has a fat right tail if there is a positive exponent α>0, called tailed index, such that P(X>x)xα. The relation refers to the asymptotic equivalence of functions, meaning that their ratio tends to a positive constant. Note that, when the density is a continuous function, it tends to 0 when x+. Hence, by L’Hospital’s rule, the existence of tailed index is equivalent to saying that the density decays like a power law, whenever the density is a continuous function.

    For a risk outcome over a finite interval (c0,c1), c0,<c1, however, its density can be + when approaching boundaries c0 and c1. Let y0 be the largest lower bound for all values of y under (2.2), and y1 the smallest upper bound. We assume y0=c0 and y1=c1.

    We say that an interval distribution has a fat right tail if the limit limyy1  g(y,a,b)=+, and a fat left tail if limyy+0  g(y,a,b)=+, where yy+0 and yy1 denote respectively y approaching y0 from the right-hand-side, and y1 from the left-hand-side. For simplicity, we write yy0 for yy+0, and yy1 for yy1.

    Given α>0, we say that an interval distribution has a fat right tail with tailed index α if limyy1  g(y,a,b)(y1y)β=+ whenever 0<β<α, and limyy1  g(y,a,b)(y1y)β=0 for β>α. Similarly, an interval distribution has a fat left tail with tailed index α if limyy0  g(y,a,b)(yy0)β=+ whenever 0<β<α, and limyy0  g(y,a,b)(yy0)β=0 for β>α. Here the status at β=α is left open. There are examples (Remark 3.4), where an interval distribution has a fat right tail with tailed index α, but the limit limyy1  g(y,a,b)(y1y)α can either be + or 0. Under this definition, a tailed index of an interval distribution with a continuous density is always larger than 0 and less or equal to 1, if it exists.

    Recall that, for a Beta distribution with parameters α>0 and β>0, its density is given by f(x)=xα1(1x)β1B(α,β), where B(α,β) is the Beta function . Under the above definition, Beta distribution has a fat right tail with tailed index (1β) when 0<β<1, and a fat left tail with tailed index (1α) when 0<α<1.

    Next, because the derivative of Φ is assumed to be continuous and positive, it is strictly monotonic. Hence Φ1(y) is defined. Let

    z=Φ1(y) (2.7)

    Then limyy0z exists (can be ), and the same for limyy1z (can be +). Let limyy0  z=z0, and limyy1  z=z1. Rewrite g(y,a,b) as g(Φ(z),a,b) by (2.7). Let [g(Φ(z),a,b)]1β/z denote the derivative of [g(Φ(z),a,b)]1/β with respect to z.

    Lemma 2.2. Given β>0, the following statements hold:

    (ⅰ) limyy0  g(y,a,b)(yy0)β=limzz0  g(Φ(z),a,b)(Φ(z)y0)β and limyy1  g(y,a,b)(y1y)β=limzz1  g(Φ(z),a,b)(y1Φ(z))β.

    (ⅱ) If limyy0  g(y,a,b)=+ and limzz0{[g(Φ(z),a,b)]1β/z}/ϕ(z) is 0 (resp. +), then limyy0  g(y,a,b)(yy0)β=+ (resp. 0).

    (ⅲ) If limyy1  g(y,a,b)=+ and limzz1{[g(Φ(z),a,b)]1β/z}/ϕ(z)) is 0 (resp. +), then limyy1  g(y,a,b)(y1y)β=+ (resp. 0).

    Proof. The first statement follows from the relationship y=Φ(z). For statements (ⅱ) and (ⅲ), we show only (ⅲ). The proof for (ⅱ) is similar. Notice that

    [g(y,a,b)(y1y)β]1/β=[g(y,a,b)]1/βy1y=[g(Φ(z),a,b)]1/βy1Φ(z). (2.8)

    By L’Hospital’s rule and taking the derivatives of the numerator and the denominator of (2.8) with respect to z, we have limyy1[g(y,a,b)(y1y)β]1/β=0 (resp. +) if limzz0{[g(Φ(z),a,b)]1/β/z}/ϕ(z) is 0 (resp. +). Hence limyy1  g(y,a,b)(y1y)β=+ (resp. 0).

    For tail convexity, we say that the right tail of an interval distribution is convex if g(y,a,b) is convex for y1є<y<y1 for sufficiently small є>0. Similarly, the left tail is convex if g(y,a,b) is convex for y0<y<y0+є for sufficiently small є>0. One sufficient condition for convexity for the right (resp. left) tail is gyy(y,a,b)0 when y is sufficiently close to y1 (resp. y0).

    Again, write g(y,a,b)=g(Φ(z),a,b). Let

    h(z,a,b)=log[g(Φ(z),a,b)], (2.9)

    where log(x) denotes the natural logarithmic function. Then

    g(y,a,b)=exp[h(z,a,b)]. (2.10)

    By (2.9), (2.10), using (2.6) and the relationship z=Φ1(y), we have

    gy=[hz(z)/ϕ(z)]exp[h(Φ1(y),a,b)],gyy=[hzz(z)ϕ2(z)hz(z)ϕz(z)ϕ3(z)+hz(z)hz(z)ϕ2(z)]exp[h(Φ1(y),a,b)]. (2.11)

    The following lemma is useful for checking tail convexity, it follows from (2.11).

    Lemma 2.3. Suppose ϕ(z)>0, and derivatives hz(z),hz(z), and ϕz(z), with respect to z, all exist. If hzz(z)0 and hz(z)ϕz(z)0, then gyy(y,a,b)0.

    In this section, we focus on the case where (c0,c1)=(0,1), and Φ:D(0,1) in (2.2) is the CDF of a continuous distribution . This includes, for example, the CDFs for standard normal and standard logistic distributions.

    One can explore into a wide list of densities with different choices for Φ and s under (2.2). We consider here only the following four interval distributions:

    A. sN(0,1) and Φ is the CDF for the standard normal distribution.

    B. s follows the standard logistic distribution and Φ is the CDF for the standard normal distribution.

    C. s follows the standard logistic distribution and Φ is its CDF.

    D.D. sN(0,1) and Φ is the CDF for standard logistic distribution.

    Densities for cases A, B, C, and D are given respectively in (3.3) (section 3.1), (A.1), (A.3), and (A5) (Appendix A). Tail behaviour study is summarized in Propositions 3.3, 3.5, and Remark 3.6. Sketches of density plots are provided in Appendix B for distributions A, B, and C.

    Using the notations of section 2, we have ϕ=f and Φ=F. We claim that y=Φ(a+bs) under (2.2) follows the Vasicek distribution [13,14].

    By (2.5), we have

    log(U1U2)=z2+2aza2+b2z22b2 (3.1)
    =(1b2)(za1b2)2+b21b2a22b2. (3.2)

    Therefore, we have

    g(y,a,b)=1bexp{(1b2)(za1b2)2+b21b2a22b2}. (3.3)

    Again, using the notations of section 2, we have y0=0 and y1=1. With z=Φ1(y), we have limy0  z= and limy1  z=+. Recall that a variable 0<y<1 follows a Vasicek distribution [13,14] if its density has the form:

    g(y,p,ρ)=1ρρexp{12ρ[1ρΦ1(y)Φ1(p)]2+12[Φ1(y)]2} (3.4)

    where p is the mean of y , and ρ is a parameter called asset correlation.

    Proposition 3.1. Density (3.3) is equivalent to (3.4) under the relationships:

    a=Φ1(p)1ρ  and  b=ρ1ρ. (3.5)

    Proof. A similar proof can be found in [19]. By (3.4), we have

    g(y,p,ρ)=1ρρexp{1ρ2ρ[Φ1(y)Φ1(p)/1ρ]2+12[Φ1(y)]2}
    =1bexp{12[Φ1(y)ab]2}exp{12[Φ1(y)]2}
    =U1/(bU2)=g(y,a,b).

    The following relationships are implied by (3.5):

    ρ=b21+b2, (3.6)
    a=Φ1(p)1+b2. (3.7)

    Remark 3.2. The mode of g(y,p,ρ) in (3.4) is given in [14] as Φ(1ρ12ρΦ1(p)). We claim this is the same as Φ(a1b2). By (3.6), 12ρ=1b21+b2 and 1ρ=11+b2. Therefore, we have

    1ρ12ρΦ1(p)=1+b21b2Φ1(p)=a1b2.

    This means Φ(1ρ12ρΦ1(p))=Φ(a1b2).

    Proposition 3.3. The following statements hold for g(y,a,b) given in (3.3):

    (ⅰ) g(y,a,b) is unimodal if 0<b<1 with mode given by Φ(a1b2), and is in U-shape if b>1.

    (ⅱ) If b>1,then  g(y,a,b) has a fat left tail and a fat right tail with tailed index (11/b2).

    (ⅲ) If b>1, both tails of g(y,a,b) are convex , and is globally convex if in addition a=0.

    Proof. For statement (ⅰ), we have (1b2)<0 when 0<b<1. Therefore by (3.2) function log(U1U2) reaches its unique maximum at z=a1b2, resulting in a value for the mode at Φ(a1b2). If b>1, then (1b2)>0, thus by (3.2), g(y,a,b) is first decreasing and then increasing when y varying from 0 to 1. This means (y,a,b) is in U-shape.

    Consider statement (ⅱ). First by (3.3), if b>1, then limy1  g(y,a,b)=+ and limy0  g(y,a,b)=+. Thus g(y,a,b) has a fat right and a fat left tail. Next for tailed index, we use Lemma 2.2 (ⅱ) and (ⅲ). By (3.1),

    [g(Φ(z),a,b)]1/β=b1/βexp((b21)z2+2aza22βb2) (3.8)

    By taking the derivative of (3.8) with respect to z and noting that ϕ(z)=12πexp(z22), we have

    {[g(Φ(z),a,b)]1β/z}/ϕ(z)=2πb1β(b21)z+aβb2exp((b21)z2+2aza22βb2+z22). (3.9)

    Thus limz+{[g(Φ(z),a,b)]1β/z}/ϕ(z) is 0 if b21βb2>1, and is + if b21βb2<1. Hence by Lemma 2.2 (ⅲ), g(y,a,b) has a fat right tail with tailed index (11/b2). Similarly, for the left tail, we have by (3.9)

    {[g(Φ(z),a,b)]1β/z}/ϕ(z)=2πb1β(b21)z+aβb2exp((b21)z2+2aza22βb2+z22). (3.10)

    Thus limz{[g(Φ(z),a,b)]1β/z}/ϕ(z) is 0 if b21βb2>1, and is + if b21βb2<1. Hence g(y,a,b) has a fat left tail with tailed index (11/b2) by Lemma 2.2 (ⅱ).

    For statement (ⅲ), we use Lemma 2.3. By (2.9) and using (3.2), we have

    h(z,a,b)=log(U1bU2)=(1b2)(za1b2)2+b21b2a22b2log(b).

    When b>1, it is not difficult to check out that hzz(z)0 and hz(z)ϕz(z)0 when z± or when a=0.

    Remark 3.4. Assume β=(11/b2) and b>1. By (3.9), we see

    limz+{[g(Φ(z),a,b)]1β/z}/ϕ(z)

    is + for a=0, and is 0 for a>0. Hence for this β, the limit limy1  g(y,a,b)(1y)β can be either 0 or +, depending on the value of a.

    For these distributions, we again focus on their tail behaviours. A proof for the next proposition can be found in Appendix A.

    Proposition 3.5. The following statements hold:

    (a) Density g(y,a,b) has a fat left tail and a fat right tail for case B for all b>0, and for case C if b>1. For case D, it does not have a fat right tail nor a fat left tail for any b>0.

    (b) The tailed index of g(y,a,b) for both right and left tails is 1 for case B for all b>0, and is (11b) for case C for B for b>1.

    Remark 3.6. Among distributions A, B, C, and Beta distribution, distribution B gets the highest tailed index of 1, independent of the choices of b>0.

    In this section, we assume that Φ in (2.2) is a function from R1 to (0,1) with positive continuous derivatives. We focus on parameter estimation algorithms for model (1.1).

    First, we consider a simple case, where risk outcome y is driven by a model:

    y=Φ(v+bs), (4.1)

    where b>0 is a constant, v=a0+a1x1++akxk, and sN(0,1), independent of fixed effects x=(x1,x2,,xk). The function Φ does not have to be the standard normal CDF. But when Φ is the standard normal CDF, the expected value E(y|x) can be evaluated by the formula ES[Φ(a+bs)]=Φ(a1+b2) [12].

    Given a sample {(x1i,x2i,,xki,yi)}ni=1, where (x1i,x2i,,xki,yi) denotes the ith data point of the sample, let zi=Φ1(yi). and vi=a0+a1x1i++akxki. By (2.3), the log-likelihood function for model (4.1) is:

    LL=ni=1{logf(zivib)logϕ(zi)logb} (4.2)

    where f is the density of s. The part of ni=1logϕ(zi) is constant, which can be dropped off from the maximization.

    Recall that the least squares estimators of a0,a1,,ak, as a row vector, that minimize the sum squares

    SS=ni=1(zivi)2 (4.3)

    has a closed form solution given by the transpose of (XTX)1XTZ [5,9] whenever the design matrix X has a rank of k, where

    X=1x11xk11x12xk21x1nxkn,Z=z1z2zn.

    The next proposition shows there exists an analytical solution for the parameters of model (4.1).

    Proposition 4.1. Given a sample {(x1i,x2i,,xki,yi)}ni=1, assume that the design matrix has a rank of k. If sN(0,1), then the maximum likelihood estimates of parameters (a0,a1,,a), as a row vector, and parameter b are respectively given by the transpose of (XTX)1XTZ, and b2=1nni=1(zivi)2. In absence of fixed effects {x1,x2,,xk}, parameters a0 and b2 degenerate respectively to the sample mean and variance of z1, z2,,zn.

    Proof. Dropping off the constant term from (4.2) and noting f(z)=12πexp(z22), we have

    LL=12b2ni=1(zivi)2nlogb, (4.4)

    Hence the maximum likelihood estimates (a0,a1,,ak) are the same as least squares estimators of (4.3), which are given by the transpose of (XTX)1XTZ. By taking the derivative of (4.4) with respect to b and setting it to zero, we have b2=1nni=1(zivi)2.

    Next, we consider the general case of model (1.1), where the risk outcome y is driven by a model:

    y=Φ[v+ws], (4.5)

    where parameter w is formulated as w=exp(u), and u=b0+b1x1++bkxk. We focus on the following two cases:

    (a) sN(0,1),

    (b) s is standard logistic.

    Given a sample {(x1i,x2i,,xki,yi)}ni=1, let wi=exp(b0+b1x1i++bkxki) and ui=b0+b1x1i++bkxki. The log-likelihood functions for model (4.5), dropping off the constant part log(U2), for cases (a) and (b) are given respectively by (4.6) and (4.7):

    LL=ni=112[(zivi)2/w2iui], (4.6)
    LL=ni=1{(zivi)/wi2log[1+exp[(zivi)/wi]ui}, (4.7)

    Recall that a function is log-concave if its logarithm is concave. If a function is concave, a local maximum is a global maximum, and the function is unimodal. This property is useful for searching maximum likelihood estimates.

    Proposition 4.2. The functions (4.6) and (4.7) are concave as a function of (a0,a1,,ak). As a function of (b0,b1,,bk), (4.6) is concave.

    Proof. It is well-known that, if f(x) is log-concave, then so is f(Az+b), where Az+b : RmR1 is any affine transformation from the m-dimensional Euclidean space to the 1-dimensional Euclidean space. For (4.6), the function f(x)=(zv)2exp(2u) is concave as a function of v, thus function (4.6) is concave as a function of (a0,a1,,ak). Similarly, this function f(x) is concave as a function of u, so (4.6) is concave as a function of (b0,b1,,bk).

    For (4.7), the linear part (zivi)exp(ui), as a function of (a0,a1,,ak), in (4.7) is ignored. For the second part in (4.7), we know log{1+exp[(zv)/exp(u)]}, as a function of v, is the logarithm of the CDF of a logistic distribution. It is well-known that the CDF for a logistic distribution is log-concave. Thus (4.7) is concave with respect to (a0,a1,,ak).

    In general, parameters (a0,a1,,ak) and (b0,b1,,bk) in model (4.5) can be estimated by the algorithm below.

    Algorithm 4.3. Follow the steps below to estimate parameters of model (4.5):

    (a) Given (b0,b1,,bk), estimate (a0,a1,,ak) by maximizing the log-likelihood function;

    (b) Given (a0,a1,,ak), estimate (b0,b1,,bk) by maximizing the log-likelihood function;

    (c) Iterate (a) and (b) until a convergence is reached.

    With the interval distributions introduced in this paper, models with a random effect can be fitted for a continuous risk outcome by maximum likelihood approaches assuming an interval distribution. These models provide an alternative regression tool to the Beta regression model and fraction response model, and a tool for tail risk assessment as well.

    Authors are very grateful to the third reviewer for many constructive comments. The first author is grateful to Biao Wu for many valuable conversations. Thanks also go to Clovis Sukam for his critical reading for the manuscript.

    We would like to thank you for following the instructions above very closely in advance. It will definitely save us lot of time and expedite the process of your paper's publication.

    The views expressed in this article are not necessarily those of Royal Bank of Canada and Scotiabank or any of their affiliates. Please direct any comments to Bill Huajian Yang at h_y02@yahoo.ca.



    [1] M. Hatami, D. D. Ganji, M. Gorji-Bandpy, Experimental and numerical analysis of the optimized finned-tube heat exchanger for OM314 diesel exhaust exergy recovery, Energy Convers. Manage., 97 (2014), 26–41. https://doi.org/10.1016/j.enconman.2015.03.032 doi: 10.1016/j.enconman.2015.03.032
    [2] M. Ghazikhani, M. Hatami, B. Safari, The effect of alcoholic fuel additives on exergy parameters and emissions in a two-stroke gasoline engine, Arab. J. Sci. Eng., 39 (2014), 2117–2125. https://doi.org/10.1007/s13369-013-0738-3 doi: 10.1007/s13369-013-0738-3
    [3] M. Hatami, D. D. Ganji, Thermal performance of circular convective-radiative porous fins with different section shapes and materials, Energy Convers. Manage., 76 (2013), 185–193. https://doi.org/10.1016/j.enconman.2013.07.040 doi: 10.1016/j.enconman.2013.07.040
    [4] M. Turkyilmazoglu, Heat transfer from moving exponential fins exposed to heat generation, Int. J. Heat Mass Transfer, 116 (2018), 346–351. https://doi.org/10.1016/j.ijheatmasstransfer.2017.08.091 doi: 10.1016/j.ijheatmasstransfer.2017.08.091
    [5] E. Cuce, P. M. Cuce, Homotopy perturbation method for temperature distribution, fin efficiency and fin effectiveness of convective straight fins with temperature-dependent thermal conductivity, J. Mech. Eng. Sci., 227 (2013), 1754–1760. https://doi.org/10.1177/0954406212469579 doi: 10.1177/0954406212469579
    [6] A. Y. Cengel, Introduction to Thermodynamics and Heat Transfer, Second Edition, McGraw-Hill Companies, 2008.
    [7] S. A. Atouei, K. Hosseinzadeh, M. Hatamic, S. E. Ghasemid, S. A. R. Sahebi, D. D. Ganji, Heat transfer study on convective-radiative semi-spherical fins with temperature-dependent properties and heat generation using efficient computational methods, Appl. Therm. Eng., 89 (2015), 299–305. https://doi.org/10.1016/j.applthermaleng.2015.05.084 doi: 10.1016/j.applthermaleng.2015.05.084
    [8] K. Hosseinzadeh, E. Montazer, M. B. Shafii, A. R. D. Ganji, Solidification enhancement in triplex thermal energy storage system via triplets fins configuration and hybrid nanoparticles, J. Energy Storage, 34 (2021), 102177. https://doi.org/10.1016/j.est.2020.102177 doi: 10.1016/j.est.2020.102177
    [9] M. Hatami, D. D. Ganji, Optimization of the longitudinal fins with different geometries for increasing the heat transfer, in ISER 10th International Conference, Kuala Lumpur, Malaysia, 2015.
    [10] M. Turkyilmazoglu, Heat transfer from moving exponential fins exposed to heat generation, Int. J. Heat Mass Transfer, 116 (2018), 346–351. https://doi.org/10.1016/j.ijheatmasstransfer.2017.08.091 doi: 10.1016/j.ijheatmasstransfer.2017.08.091
    [11] B. Kundu, D. Bhanja, K. S. Lee, A model on the basis of analytics for computing maximum heat transfer in porous fins, Int. J. Heat Mass Transfer, 55 (2012), 7611–7622. https://doi.org/10.1016/j.ijheatmasstransfer.2012.07.069 doi: 10.1016/j.ijheatmasstransfer.2012.07.069
    [12] Z. Din, A. Ali, S. Ullah, G. Zaman, Investigation of heat transfer from convective and radiative stretching/shrinking rectangular fins, Math. Probl. Eng., 2022 (2022). https://doi.org/10.1155/2022/1026698
    [13] W. Ahmad, K. S. Syed, M. Ishaq, A. Hassan, Z. Iqbal, Numerical study of conjugate heat transfer in a double-pipe with exponential fins using DGFEM, Appl. Therm. Eng., 111 (2017), 1184–1201. https://doi.org/10.1016/j.applthermaleng.2016.09.171 doi: 10.1016/j.applthermaleng.2016.09.171
    [14] M. M. Rashidi, T. Hayat, T. Keimanesh, H. Yousefian, A study on heat transfer in a second-grade fluid through a porous medium with the modified differential transform method, Heat Transfer Asian Res., 42 (2013), 31–45. https://doi.org/10.1002/htj.21030 doi: 10.1002/htj.21030
    [15] E. Erfani, M. M. Rashidi, A. B. Parsa. The modified differential transform method for solving off-centered stagnation flow toward a rotating disc, Int. J. Comput. Methods, 7 (2010), 655–670. https://doi.org/10.1142/S0219876210002404
    [16] Y. Huang, X. Li, Exact and approximate solutions of convective-radiative fins with temperature-dependent thermal conductivity using integral equation method, Int. J. Heat Mass Transfer, 150 (2020), 119303. https://doi.org/10.1016/j.ijheatmasstransfer.2019.119303 doi: 10.1016/j.ijheatmasstransfer.2019.119303
    [17] M. M. Rashidi, E. Erfani, New analytical method for solving Burgers and nonlinear heat transfer equations and comparison with HAM, Comput. Phys. Commun., 180 (2009), 1539–1544.
    [18] S. Panda, A. Bhowmik, R. Das, R. Repaka, S. C. Martha, Application of Homotopy analysis method and inverse solution of a rectangular wet fin, Energy Convers. Manage., 80 (2014), 303–318. https://doi.org/10.1016/j.cpc.2009.04.009 doi: 10.1016/j.cpc.2009.04.009
    [19] R. K. Singla, R. Das, Application of decomposition method and inverse prediction of parameters in a moving fin, Energy Convers. Manage., 84 (2014), 268–281. https://doi.org/10.1016/j.enconman.2014.04.045 doi: 10.1016/j.enconman.2014.04.045
    [20] C. Y. Zhang, X. F. Li, Temperature distribution of conductive-convective-radiative fins with temperature-dependent thermal conductivity, Int. Comm. Heat Mass Transfer, 117 (2020), 104799. https://doi.org/10.1016/j.icheatmasstransfer.2020.104799 doi: 10.1016/j.icheatmasstransfer.2020.104799
    [21] S. W. Sun, X. F. Li, Exact solution of the nonlinear fin problem with exponentially temperature-dependent thermal conductivity and heat transfer coefficient, Pramana J. Phys., 94 (2020), 1–10. https://doi.org/10.1007/s12043-020-01971-4 doi: 10.1007/s12043-020-01971-4
    [22] A. K. Asl, S. Hossainpour, M. M. Rashidi, M. A. Sheremet, Z. Yang, Comprehensive investigation of solid and porous fins influence on natural convection in an inclined rectangular enclosure, Int. J. Heat Mass Transfer, 133 (2019), 729–744. https://doi.org/10.1016/j.ijheatmasstransfer.2018.12.156 doi: 10.1016/j.ijheatmasstransfer.2018.12.156
    [23] S. Maalej, A. Zayoud, I. Abdelaziz, I. Saad, M. C. Zaghdoudi, Thermal performance of finned heat pipe system for Central Processing Unit cooling, Energy Convers. Manage., 218 (2020), 112977. https://doi.org/10.1016/j.enconman.2020.112977 doi: 10.1016/j.enconman.2020.112977
    [24] A. A. Joneidi, D. D. Ganji, M. Babaelahi, Differential Transformation Method to determine fin efficiency of convective straight fins with temperature dependent thermal conductivity, Int. Commun. Heat Mass Transfer, 36 (2009), 757–762. https://doi.org/10.1016/j.icheatmasstransfer.2009.03.020 doi: 10.1016/j.icheatmasstransfer.2009.03.020
    [25] C. H. Chiu, C. K. Chen, Applications of Adomian decomposition procedure to the analysis of convective radiative fins, J. Heat Transfer, 125 (2003), 312–316. https://doi.org/10.1115/1.1532012 doi: 10.1115/1.1532012
    [26] D. Lesnic, P. J. Heggs, A decomposition method for power-law fin-type problems, Int. Commun. Heat Mass Transfer, 31 (2004), 673–682. https://doi.org/10.1016/S0735-1933(04)00054-5 doi: 10.1016/S0735-1933(04)00054-5
    [27] R. Das, B. Kundu, Prediction of heat-generation and electromagnetic parameters from temperature response in porous fins, J. Thermophys. Heat Transfer, 35 (2021), 761–769. https://doi.org/10.2514/1.T6224 doi: 10.2514/1.T6224
    [28] R. Das, B. Kundu, Simultaneous estimation of heat generation and magnetic field in a radial porous fin from surface temperature information, Int. Commun. Heat Mass Transfer, 127 (2021), 105497. https://doi.org/10.1016/j.icheatmasstransfer.2021.105497 doi: 10.1016/j.icheatmasstransfer.2021.105497
    [29] D. Bhanja, B. Kundu, Thermal analysis of a constructal T-shaped porous fin with radiation effects, Int. J. Refrig., 31 (2011), 337–352. https://doi.org/10.1016/j.ijrefrig.2011.04.003 doi: 10.1016/j.ijrefrig.2011.04.003
    [30] B. Kundu, D. Bhanja, An analytical prediction for performance and optimum design analysis of porous fins, Int. J. Refrig., 31 (2011), 1483–1496. https://doi.org/10.1016/j.ijrefrig.2010.06.011 doi: 10.1016/j.ijrefrig.2010.06.011
    [31] R. Das, B. Kundu, Prediction of heat generation in a porous fin from surface temperature, J. Thermophys. Heat Transfer, 31 (2017), 781–790. https://doi.org/10.2514/1.T5098 doi: 10.2514/1.T5098
    [32] R. Das, Forward and inverse solutions of a conductive, convective and radiative cylindrical porous fin, Energy Convers. Manage., 87 (2014), 96–106. https://doi.org/10.1016/j.enconman.2014.06.096 doi: 10.1016/j.enconman.2014.06.096
    [33] R. Das, D. K. Prasad. Prediction of porosity and thermal diffusivity in a porous fin using differential evolution algorithm, Swarm Evol. Comput., 23 (2015), 27–39. https://doi.org/10.1016/j.swevo.2015.03.001
    [34] B. Kundu, S. J. Yook, An accurate approach for thermal analysis of porous longitudinal, spine and radial fins with all nonlinearity effects-analytical and unified assessment, Appl. Math. Comput., 402 (2021), 126124. https://doi.org/10.1016/j.amc.2021.126124 doi: 10.1016/j.amc.2021.126124
    [35] G. A. Oguntala, R. A. Abd-Alhameed, G. M. Sobamowo, N. Eya, Effects of particles deposition on thermal performance of a convective-radiative heat sink porous fin of an electronic component, Therm. Sci. Eng. Prog., 6 (2018), 177–185. https://doi.org/10.1016/j.tsep.2017.10.019
    [36] M. A. Vatanparast, S. Hossainpour, A. Keyhani-Asl, S. Forouzi, Numerical investigation of total entropy generation in a rectangular channel with staggered semi-porous fins, Int. Commun. Heat Mass Transfer, 111 (2020), 104446. https://doi.org/10.1016/j.icheatmasstransfer.2019.104446 doi: 10.1016/j.icheatmasstransfer.2019.104446
    [37] M. Turkyilmazoglu, Exact solutions to heat transfer in straight fins of varying exponential shape having temperature dependent properties, Int. J. Therm. Sci., 55 (2012), 69–79. https://doi.org/10.1016/j.ijthermalsci.2011.12.019 doi: 10.1016/j.ijthermalsci.2011.12.019
    [38] Z. U. Din, A. Ali, G. Zaman, Entropy generation in moving exponential porous fins with natural convection, radiation and internal heat generation, Arch. Appl. Mech., 92 (2022), 933–944. https://doi.org/10.1007/s00419-021-02081-2 doi: 10.1007/s00419-021-02081-2
    [39] Z. U. Din, A. Ali, M. D. la Sen, G. Zaman, Entropy generation from convective-radiative moving exponential porous fins with variable thermal conductivity and internal heat generations, Sci. Rep., 12 (2022), 1791. https://doi.org/10.1038/s41598-022-05507-1
    [40] M. Hatami, A. Hasanpour, D. D. Ganji, Heat transfer study through porous fins (Si3N4 and AL) with temperature-dependent heat generation, Energy Convers. Manage., 74 (2013), 9–16. https://doi.org/10.1016/j.enconman.2013.04.034 doi: 10.1016/j.enconman.2013.04.034
    [41] M. Hatami, D. D. Ganji, Thermal and flow analysis of microchannel heat sink (MCHS) cooled by Cu-water nanofluid using porous media approach and least square method, Energy Convers. Manage., 78 (2014), 347–358. https://doi.org/10.1016/j.enconman.2013.10.063 doi: 10.1016/j.enconman.2013.10.063
    [42] M. Turkyilmazoglu, Thermal performance of optimum exponential fin profiles subjected to a temperature jump, Int. J. Numer. Methods Heat Fluid Flow, 32 (2021), 1002–1011. https://doi.org/10.1108/HFF-02-2021-0132 doi: 10.1108/HFF-02-2021-0132
    [43] A. R. A. Khaled, Thermal characterizations of exponential fin systems, Math. Probl. Eng., (2010), 765729. https://doi.org/10.1155/2010/765729
    [44] M. F. Najafabadi, H. T. Rostami, K. Hosseinzadeh, D. D. Ganji, Thermal analysis of a moving fin using the radial basis function approximation, Heat Transfer, 50 (2021), 7553–7567. https://doi.org/10.1002/htj.22242 doi: 10.1002/htj.22242
    [45] S. Hosseinzadeh, K. Hosseinzadeh, A. Hasibi, D. D. Ganji, Thermal analysis of moving porous fin wetted by hybrid nanofluid with trapezoidal, concave parabolic and convex cross sections, Case Stud. Therm. Eng., 30 (2022), 101757. https://doi.org/10.1016/j.csite.2022.101757 doi: 10.1016/j.csite.2022.101757
    [46] M. A. E. Moghaddam, M. R. H. S. Abandani, K. Hosseinzadeh, M. B. Shafii, D. D. Ganji, Metal foam and fin implementation into a triple concentric tube heat exchanger over melting evolution, Theor. Appl. Mech. Lett., (2022), 100332. https://doi.org/10.1016/j.taml.2022.100332
    [47] B. Jalili, N. Aghaee, P. Jalili, D. D. Ganji, Novel usage of the curved rectangular fin on the heat transfer of a double-pipe heat exchanger with a nanofluid, Case Stud. Therm., (2022), 102086. https://doi.org/10.1016/j.csite.2022.102086
    [48] B. Jalili, S. Sadighi, P. Jalili, D. D. Ganji, Characteristics of ferrofluid flow over a stretching sheet with suction and injection, Case Stud. Therm., 14 (2019), 100470. https://doi.org/10.1016/j.csite.2022.102086 doi: 10.1016/j.csite.2022.102086
    [49] B. Jalili, S. Sadighi, P. Jalili, D. D. Ganji, Effect of magnetic and boundary parameters on flow characteristics analysis of micropolar ferrofluid through the shrinking sheet with effective thermal conductivity, Chin. J. Phys., 71 (2021), 136–150. https://doi.org/10.1016/j.cjph.2020.02.034 doi: 10.1016/j.cjph.2020.02.034
    [50] P. Jalili, D. D. Ganji, B. Jalili, D. D. Ganji, Evaluation of electro-osmotic flow in a nanochannel via semi-analytical method, Therm. Sci., 16 (2012), 1297–1302. http://DOI:10.2298/TSCI1205297J doi: 10.2298/TSCI1205297J
    [51] M. Turkyilmazoglu, Efficiency of heat and mass transfer in fully wet porous fins: exponential fins versus straight fins, Int. J. Refrig., 46 (2014), 158–164. https://doi.org/10.1016/j.ijrefrig.2014.04.011 doi: 10.1016/j.ijrefrig.2014.04.011
    [52] B. Kundu, K. S. Lee, Analytic solution for heat transfer of wet fins on account of all nonlinearity effects, Energy, 41 (2012), 354–367. https://doi.org/10.1016/j.energy.2012.03.004 doi: 10.1016/j.energy.2012.03.004
    [53] R. das, K. T. Ooi, Predicting multiple combination of parameters for designing a porous fin subjected to a given temperature requirement, Energy Convers. Manage., 66 (2013), 211–219. https://doi.org/10.1016/j.enconman.2012.10.019 doi: 10.1016/j.enconman.2012.10.019
    [54] M. Torabi A. Aziz, K. Zhang, A comparative study of longitudinal fins of rectangular, trapezoidal and concave parabolic profiles with multiple nonlinearities, Energies, 51 (2013), 243–256. https://doi.org/10.1016/j.energy.2012.11.052 doi: 10.1016/j.energy.2012.11.052
  • Reader Comments
  • © 2022 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(1900) PDF downloads(63) Cited by(9)

Figures and Tables

Figures(11)  /  Tables(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog