Fractional diffusion equation and diffusive stresses

  • Published on

  • View

  • Download

Embed Size (px)


<ul><li><p>Fractional diffusion equation and diffusive stressesYuriy Povstenko1,1 Institute of Mathematics and Computer Science, Jan Dugosz University of Czestochowa, Armii Krajowej 13/15,</p><p>42-200 Czestochowa, Poland.</p><p>A quasi-static uncoupled theory of diffusive stresses based on the time-fractional diffusion equation is considered. The Caputofractional derivative of order is used. In particular, the proposed theory interpolates the classical theory of diffusive stressesand that without energy dissipation. The fundamental solution to the second Cauchy problem for the fractional diffusionequation in a plane as well as the associated diffusive stresses are obtained in the case = 3/2.</p><p>1 IntroductionThe theory of diffusive stresses is an integration of the theory of diffusion and the theory of elasticity. In this paper we restrictthe discussion to the effect of diffusion on the deformation of a linear isotropic elastic solid. The inverse effect of deformationon diffusion will not be considered.</p><p>The conventional theory of diffusive stresses is formulated on the principles of the classical theory of diffusion, specificallyon the classical Fick law, which relates the matter flux vector to the concentration gradient. It is well known that frommathematical wievpoint the Fick law and the Fourier law as well as and the diffusion equation and the heat conductionequation are indentical. During the last three decades, nonclassical theories in which the Fick law and the Fourier law and,consequently, the diffusion and the heat conduction equations are replaced by more general equations, have been formulated.For an extensive bibliography on this subject see, for example, [1-5]. The time-nonlocal dependence between the flux vectorsand corresponding gradients with long-tale power kernel can be interpreted in terms of fractional integrals and derivativesand yields the fractional diffusion (or heat conduction) equation with time derivative of fractional order .</p><p>A quasi-static uncoupled theory of diffusive (or thermal) stresses based on this equation was proposed in [6, 7]. Becausethe time-fractional diffusion (or heat conduction) equation in the case 1 2 interpolates the standard heat conductionequation ( = 1) and the wave equation ( = 2), the proposed theory interpolates the classical thermoelasticity and thethermoelasticity without energy dissipation introduced by Green and Naghdi [8].</p><p>2 Statement of the problemThe stressed-strained state of a solid is governed by the equilibrium equation in terms of displacements</p><p>u + ( + )grad divu = Kgrad c, (1)the stress-strain-concentration relation</p><p> = 2e + ( tr e Kc)I, (2)the time-fractional diffusion equation</p><p>c</p><p>t= ac + Q(x, t), (3)</p><p>where u is the displacement vector, the stress tensor, e the linear strain tensor, c the concentration, a the diffusivity co-efficient, Q(x, t) the mass source, and are Lame constants, K = + 2/3, is the diffusive coefficient of volumetricexpansion, I denotes the unit tensor, </p><p>t is the Caputo fractional derivative [9,10]</p><p>dfdt =</p><p>1</p><p>(n ) t0</p><p>(t )n1 dnf()</p><p>dn d, n 1 &lt; &lt; n, (4)</p><p>with the following Laplace transform rule</p><p>L</p><p>{df(t)</p><p>dt</p><p>}= sL {f(t)} </p><p>n1k=0</p><p>f (k)(0+)s1k, n 1 &lt; &lt; n. (5)</p><p> Corresponding author E-mail:</p><p>PAMM Proc. Appl. Math. Mech. 7, 20400072040008 (2007) / DOI 10.1002/pamm.200700134</p><p> 2007 WILEY-VCH Verlag GmbH &amp; Co. KGaA, Weinheim</p><p> 2007 WILEY-VCH Verlag GmbH &amp; Co. KGaA, Weinheim</p></li><li><p>Equation (3) should be subject to initial conditions:t = 0 : c = U(x) for 0 &lt; 2 (6)</p><p>and</p><p>t = 0 :c</p><p>t= W (x) for 1 &lt; 2. (7)</p><p>3 The fundamental solution to the second Cauchy problemIn two-dimensional case the diffusive stresses exerted by a source of diffusion Q(r, t) = q2r +(r) +(t) were considered in[6], the stresses associated with the solution to the first Cauchy problem with U(r) = u2r +(r) were obtained in [7]. In thispaper, we consider the second Cauchy problem with W (r) = w2r +(r). The solution is obtained using the Laplace integraltransform with respect to time t and the Hankel transform with respect to the radial coordinate r and has the following form</p><p>c =w</p><p>2</p><p>s2</p><p>s + a2, (8)</p><p>where asterisk denotes the transforms, s is the Laplace transform variable, is the Hankel transform variable.For = 3/2 we obtain</p><p>c =w</p><p>2</p><p>0</p><p>H3/2(, t)J0(r) d. (9)</p><p>As in a case of classical thermoelasticity or diffusional elasticity we use the representation of stresses in terms of displace-ment potential and get</p><p>rr + = mw</p><p>0</p><p>H3/2(, t)J0(r) d, (10)</p><p>rr = mw</p><p>0</p><p>H3/2(, t)J2(r) d, (11)</p><p>where</p><p>H3/2(, t) = L1</p><p>{1</p><p>s(s3/2 + a2)</p><p>}=</p><p>t</p><p>3</p><p>{eerfc 2e/2 cos</p><p>(3</p><p>2 +</p><p>3</p><p>)</p><p>+4</p><p> 10</p><p>exp</p><p>[1</p><p>2(1 v2)</p><p>]cos</p><p>[3</p><p>2(1 v2) </p><p>3</p><p>]dv}</p><p>, = a2/34/3t, (12)</p><p>Jn(r) is the Bessel function of the first kind of the order n.</p><p>References[1] D. S. Chandrasekharaiah, Appl. Mech. Rev. 39, 355 (1986).[2] D. S. Chandrasekharaiah, Appl. Mech. Rev. 51, 705 (1998).[3] K. K. Tamma and X. Zhou, J. Thermal Stresses 21, 405 (1998).[4] R. B. Hetnarski and J. Ignaczak, J. Thermal Stresses 22, 451 (1999).[5] R. Metzler and J. Klafter, Phys. Rep. 339, 1 (2000).[6] Y. Z. Povstenko, Int. J. Engng Sci. 43, 977 (2005).[7] Y. Z. Povstenko, J. Thermal Stresses 28, 83 (2005).[8] A. E. Green and P. M. Naghdi, J. Elast. 31, 189 (1993).[9] R. Gorenflo and F. Mainardi, in: Fractals and Fractional Calculus in Continuum Mechanics, edited by A. Carpinteri and F. Mainardi</p><p>(Springer, New York, 1997), p. 223.[10] A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations (Elsevier, Amsterdam,</p><p>2006), p. 90.</p><p> 2007 WILEY-VCH Verlag GmbH &amp; Co. KGaA, Weinheim</p><p>ICIAM07 Contributed Papers 2040008</p></li></ul>