MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2014/180764 180764 Research Article On the Solutions of a Stefan Problem with Variable Latent Heat http://orcid.org/0000-0003-0189-9436 Lazhar Bougoffa Zenkour Ashraf M. 1 Department of Mathematics Faculty of Science Al Imam Mohammad Ibn Saud Islamic University (IMSIU) P.O. Box 90950, Riyadh 11623 Saudi Arabia imamu.edu.sa 2014 1372014 2014 17 04 2014 03 07 2014 13 7 2014 2014 Copyright © 2014 Bougoffa Lazhar. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper is concerned with a reliable treatment of the Stephan problem with variable latent heat. The power series solutions have been presented. The method is examined for computational efficiency and accuracy. Also, an analytical solution based on a similarity variable is presented in the case when the Dirichlet boundary condition at the water-ice interface depends on time.

1. Introduction

The solution of the Stefan problem with variable latent heat consists of finding u and the moving melt interface s such that (1)ut=2ux2,0xs(t), which is the governing equation, subject to the following Neumann boundary condition at the left end of the domain describing the inlet heat flux: (2)-ux(t,0)=h(t),t>0, the Dirichlet boundary condition at the water-ice interface: (3)u(t,s(t))=u*,t>0, the Stefan condition: (4)sdsdt=-ux(t,s(t)),s(0)=b,t>0, and the initial temperature distribution: (5)u(0,x)=φ(x),0<x<s(0)=b. On comparison of this condition with the one-phase Stefan condition, we observe that the latent heat term s is not a constant but, rather, a linear function of position . Problems such as this have been treated under the name of Stefan problems with variable latent heat .

An analytical solution for this Stefan problem is presented in , by introducing the similarity variables ξ=x/t and seeking the solution in the form u(t,x)=F(ξ) with F being an unknown function. Accordingly, it is natural that the interface location s(t) should be proportional to t; that is, s(t)=At, where A is a constant.

The purpose of this paper is to apply the Adomian decomposition method  to find the solution of (1), (2), and (5), that is, the temperature distribution of the water u(t,x), and then use a direct method to determine the position of the ice-water interface as a function of time. Also, an analytical solution based on a similarity variable is presented in the case when the Dirichlet boundary condition at the water-ice interface depends on time of the form u*(t)=αt, where α is a constant.

Based on the Adomian decomposition method, we write (1) in Adomian’s operator-theoretic notation as (6)Lxxu=ut, where (7)Lxx=2x2. We conveniently define the inverse linear operator as (8)Lxx-1(·)=0x(·)dxdx. Applying the inverse linear operator Lxx-1 to (6) and taking into account that (u/x)(t,0)=-h(x), we obtain (9)u(t,x)=A(t)-h(t)x+0xutdxdx, where the unknown boundary condition A(t)=u(t,0) will be determined.

Define the solution u(t,x) by an infinite series of components in the form (10)u(t,x)=n=0un(t,x). Consequently, the components un can be elegantly determined by setting the recursion scheme (11)u0=A(t)-h(t)x,un+1=0xuntdxdx,n0, for the complete determination of these components. In view of (11), the components u0(t,x),u1(t,x),u2(t,x), are immediately determined as (12)u0=A(t)-h(t)x,u1=A(t)x22!-h(t)x33!,u2=A′′(t)x44!-h′′(t)x55!,. Consequently, the solution is readily found to be (13)u(t,x)=n=0A(n)(t)x2n(2n)!-n=0h(n)(t)x2n+1(2n+1)!, obtained by substituting (12) into (10).

We remark here that the unknown boundary condition A(t) can be easily determined by using the initial condition equation (5). Substituting t=0 into (13) and using the Taylor expansion of φ(x) lead to (14)n=0A(n)(0)x2n(2n)!-n=0h(n)(0)x2n+1(2n+1)!=n=0φ(n)(0)xnn!. Equating the coefficients of like power of x in both sides of (14) and taking into account that the compatibility conditions (15)h(n)(0)=-φ(2n+1)(0),n0, yield (16)A(n)(0)=φ(2n)(0),n0. Thus (17)A(t)=n=0φ(2n)(0)tn(n)!. Accordingly, the solution equation (13) is completely determined by defining the function A(t).

Once the function u(t,x) is obtained, we can rewrite the Stefan condition equation (4) in terms of the known function u(t,x) including the Dirichlet and Neumann boundary conditions. For this, integrating (1) with respect to x from 0 to x and taking into account that -(u/x)(t,0)=h(t), we obtain (18)ux=-h(t)+0xutdx. Thus (4) can be replaced by (19)sdsdt=h(t)-0s(t)utdx,s(0)=b,t>0. Using the following Leibniz rule for differentiation under the integral sign: (20)ddt0s(t)u(t,x)dx=0s(t)utdx+u(t,s(t))dsdt, and taking into account that u(t,s(t))=u*, we obtain (21)0s(t)utdx=ddt0s(t)u(t,x)dx-u*dsdt. Substituting (21) into (19), we obtain (22)sdsdt-u*dsdt+ddt0s(t)u(t,x)dx=h(t). Since (ds2/dt)=2s(ds/dt), we have (23)ddt[s2-2u*s+20s(t)u(t,x)dx]=2h(t). Applying the inverse linear operator Lt-1(·)=0t(·)dt to (23), we obtain (24)s2-2u*s+20s(t)u(t,x)dx=b2-2bu*+20th(t)dt+20bu(0,x)dx. It follows that (25)s2-2u*s+20s(t)u(t,x)dx=b2-2bu*+20th(t)dt+20bφ(x)dx. Substituting the solution equation (13) of the heat equation into (25), we obtain (26)s2-2u*s+20s(t)[n=0A(n)(t)x2n(2n)!-n=0h(n)(t)x2n+1(2n+1)!]dx=g(t), where (27)g(t)=b2-2bu*+20th(t)dt+20bφ(x)dx. Thus (28)s2-2u*s+F(t,s(t))=g(t), where (29)F(t,s(t))=2[n=0A(n)(t)s2n+1(2n+1)!-n=0h(n)(t)s2n+2(2n+2)!]. We now can determine the shoreline s(t) with time by solving the nonlinear equation (28). In order to demonstrate the feasibility and efficiency of this method, we consider the following case: (30)h(t)=e-t. If we choose φ(x)=ex, then a simple calculation leads to A(t)=et, and (28) becomes (31)s2-2u*s+etn=0snn!=g(t), where F(t,s(t))=et+s.

Let us write s(t) and g(t) in series forms s(t)=n=0sntn and g(t)=n=0gntn. Thus (32)s(t)=n=0sntn=n=01n![sntn]n. To compute [sntn]n we need the following theorem .

Theorem 1.

If u=n=0antn and F(u)=n=0bntn are convergent, then F(u)=n=0βntn is convergent, where (33)β0=b0+n=0bnA0n(a0),β1=n=0bnA1n(a0,a1),β2=n=0bnA2n(a0,a1,a2),, and Amn(a0,a1,,am) are the Adomian polynomials.

Using this theorem with the given formula bn=1/n!, we see that (34)s(t)=n=0sntn=n=0βntn,(35)β0=n=01n!s0n,β1=n=11n!(ns1s0n-1),. Consequently, we obtain the recurrence relations for the coefficients (36)cn-2u*sn+γn=gn,n0, where cn=k=0nsksn-k and γn=k=0n(1/k!)βn-k.

This solution, of course, is (37)s(t)=n=0sntn.

3. Similarity Solution

We begin our approach by introducing the similarity variable (38)ξ=xt and look for solutions of (1) in the form (39)u(t,x)=tpμ(ξ(t,x)), where the number p and the function μ are to be determined. Substituting (39) into (1) we find that (40)tp-1(pμ-d2μdξ2-ξ2dμdξ)=0 and so (41)d2μdξ2+ξ2dμdξ=pμ. This is also difficult to solve for arbitrary values of p but for special values we can do something. Define F(ξ)=μeξ4/4, and thus (41) can be transformed into (42)d2Fdξ2-ξ2dFdξ=2p+12F. This has the general solution (43)F(ξ)=Cξeξ2/4+D[-πξe-ξ2/4erf(ξ2)-2], where C and D are constants, provided p=1/2. Hence (44)μ(ξ)=Cξ+D[-πξe-ξ2/2erf(ξ2)-2e-ξ2/4]. This gives a full solution for u(t,x) in the form (45)u(t,x)=Cx+D[-πxe-x2/2terf(x2t)-2te-x2/4t]. Taking the spatial derivative of the solution u given by (45) (46)ux(t,x)=C+D[πe-x2/2terf(x2t)(x2t-1)+xt(e-x2/4t-e-3x2/4t)(x2t-1)]. From the Neumann boundary condition at x=0, we obtain (47)ux(t,0)=C. It follows that h(t)=q must be a constant. Hence (48)C=-q, and from the Dirichlet condition we get (49)-qs+D[-πse-s2/2terf(s2t)-2te-s2/4t]=u*(t). Let λ=s/2t, and thus (49) becomes (50)-2qλt+D[-2πλe-2λ2erf(λ)-2e-λ2]t=u*(t). The analytical solution based on a similarity variable  requires that u*=0. This means that u* does not depend on time. In this treatment we will assume that u* is a function of t, that is, u*(t)=αt, where α is a constant, which needs to be specified to realize an explicit solution. Thus (51)-2qλ+D[-2πλe-2λ2erf(λ)-2e-λ2]=α. Since D in (51) is a constant, it follows that λ must also be constant. Thus (52)D=α+2qλ-2πλe-2λ2erf(λ)-2e-λ2,s(t)=2λt. The Stefan condition at the free boundary x=s(t) is (53)sdsdt=-ux(t,s), where (54)sdsdt=2λ2,(55)ux(t,s)=-q+D[πe-2λ2erf(λ)(4λ2-1)+2λ(e-λ2-e-3λ2)]. By putting (53) and (54) in the Stefan condition, we get the following transcendental equation: (56)2λ2+D[πe-2λ2erf(λ)(4λ2-1)+2λ(e-λ2-e-3λ2)]=q. Thus we have proved the following theorem.

Theorem 2.

For the Stefan problem with variable latent heat (57)ut=2ux2,0xs(t),ux(t,0)=-q,t>0,u(t,s(t))=αt,t>0,sdsdt=-ux(t,s(t)),t>0. The solution is given by (58)u(t,x)=-qx+D[-πxe-x2/2terf(x2t)-2te-x2/4t],(59)s(t)=2λt, where (60)D=α+2qλ-2πλe-2λ2erf(λ)-2e-λ2,(61)q=2λ2+D[(e-λ2-e-3λ2)πe-2λ2erf(λ)(4λ2-1)+2λ(e-λ2-e-3λ2)].

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.