Heat Conduction Basic Research Part 6 potx

114 Heat Conduction – Basic Research and it follows from (2.7) and (2.8), that m  G '  U    i i    G  i 1  i1 i  G'   G'        G  G i 1  ,   i2 i1 i   G'   G'   G'   (i  1)     (2i  1)    i (   )    m G G G U    i i    i 1 i 2 i 1    (2 i  1)  G '    ( i  1)  G '         G G   (2.10) and so on Here, the prime denotes the derivative with respective to  To determine u explicitly, we take the following four steps: Step Determine the integer m by substituting Eq (2.7) along with Eq (2.8) into Eq (2.5) or (2.6), and balancing the highest-order nonlinear term(s) and the highest-order partial derivative Step Substitute Eq (2.7) with the value of m determined in Step 1, along with Eq (2.8) into  G'  Eq (2.5) or (2.6) and collect all terms with the same order of   together; the left-hand G  G'  side of Eq (2.5) or (2.6) is converted into a polynomial in   Then set each coefficient of G this polynomial to zero to derive a set of algebraic equations for k , c ,  and  i , for i  1, 2, , m Step Solve the system of algebraic equations obtained in Step 2, for k , c ,  and  i , for i  1, 2, , m , by use of Maple Step Use the results obtained in the above steps to derive a series of fundamental  G'  solutions u( ) of Eq (2.5) or (2.6) depending on   ; since the solutions of Eq (2.8) have G been well known for us, we can obtain exact solutions of Eqs (2.1) and (2.2) 2.2 The Exp-function method According to the classic Exp-function method, it is assumed that the solution of ODEs (2.5) or (2.6) can be written as g  u( )  n  f q  m  p an exp(n )  bm exp(m ) a f exp( f  )    a g exp(  g ) bp exp( p )    b q exp( q ) , (2.11) where f , g , p and q are positive integers which are unknown, to be further determined, and an and bm are unknown constants Exact Travelling Wave Solutions for Generalized Forms of the Nonlinear Heat Conduction Equation 115 A generalized form of the nonlinear heat conduction equation 3.1 Application of the (G'/G)-expansion method Introducing a complex variable  defined as Eq (2.3), Eq (1.1) becomes an ordinary differential equation, which can be written as  kcU   ak (U n )  U  U n  0, a0 (3.1) or, equivalently,  kcU   ak n(n  1)U n  2U 2  ak nU n  1U   U  U n  0, (3.2) To get a closed-form analytic solution, we use the transformation (Kabir & Khajeh, 2009; Wazwaz, 2005) 1 U ( )  V n  ( ), (3.3) kc(n  1)V V  ak n(1  2n)V 2  ak n(n  1)VV   (n  1)2 V  (n  1)2 V  0, (3.4) which will convert Eq (3.2) into According to Step 1, considering the homogeneous balance between VV  and V V in Eq (3.4) gives m   3m  1, (3.5) m  (3.6) so that  G'  Suppose that the solutions of (3.4) can be expressed by a polynomial in   as follows: G  G'  V ( )       , G   (3.7) where  and  , are constants which are unknown, to be determined later Substituting Eq (3.7) along with Eq (2.8) into Eq (3.4) and collecting all terms with the same  G'  power of   together, the left-hand side of Eq (3.4) is converted into a polynomial in G  G'    Equating each coefficient of this polynomial to zero yields a set of simultaneous G algebraic equations for  ,  , k , c ,  and  Solving the system of algebraic equations with the aid of Maple 12, we obtain the following 116 Heat Conduction – Basic Research Case A: When     Case A-1 0   n1  , 1  , k , c a 2 2   4 n a   4   4 (3.8) where  and  are arbitrary constants By using Eq (3.8), expression (3.7) can be written as V ( )    G'     , 2   4   4  G  (3.9) Substituting the general solution of (2.9) into Eq (3.9), we get the generalized travelling wave solution as follows:      4     4      C sinh     C cosh        2  1    , V ( )     2    4     4         C sinh      C cosh       2          (3.10) where  n1 x  at n a   4   inserting Eq (3.10) into Eq (3.3), it yields the following exact solution of Eq (1.1): 1     n1   n1    n1 x  at   C cosh   x  at        C sinh    2n a   2n a    u( x , t )      2    n1   n1    x  at   C sinh   x  at        C cosh    2n a   2n a                 (3.11) in which C and C are arbitrary parameters that can be determined by the related initial and boundary conditions Now, to obtain some special cases of the above general solution, we set C  ; then (3.11) leads to the formal solitary wave solution to (1.1) as follows: 1 1  n1   n1 u( x , t )       x  at    ,  2n a   2    (3.12) Exact Travelling Wave Solutions for Generalized Forms of the Nonlinear Heat Conduction Equation 117 and, when C  , the general solution (3.11) reduces to 1 1  n1   n1 u( x , t )     coth   x  at    ,  2n a   2    (3.13) Comparing the particular cases of our general solution, Eqs (3.12) and (3.13), with Wazwaz’s results (2005), Eqs (73) and (74), it can be seen that the results are exactly the same Case A-2 0  n1  1  , 1  , k , c a 2 2   4 n a   4   4 (3.14) Inserting Eq (3.14) into (3.7) yields V ( )    G'     , 2   4   4  G  (3.15) Substituting the general solution of (2.9) into Eq (3.15), we obtain      4     4      C sinh     C cosh        2  1    , V ( )     2    4     4         C sinh      C cosh       2          where    (3.16) n1 x  at n a   4   Substituting Eq (3.16) into the transformation (3.3) leads to the generalized solitary wave solution of Eq (1.1) as follows: 1     n1   n1    n1 x  at   C cosh   x  at        C sinh    2n a   2n a    u( x , t )      2    n1   n1    x  at   C sinh   x  at        C cosh    2n a   2n a                 (3.17) Similarly, to derive some special cases of the above general solution (3.17), we choose C  ; then (3.17) leads to 1 1  n1   n1 u( x , t )       x  at    ,  2n a   2    (3.18) 118 Heat Conduction – Basic Research and, when C  , the general solution (3.17) reduces to 1 1  n1   n1 u( x , t )     coth   x  at    ,  2n a   2    (3.19) Validating our results, Eqs (3.18) and (3.19), with Wazwaz’s solutions (2005), Eqs (71) and (72), we can conclude that the results are exactly the same Case B: When     Case B-1 0  i  , 1  2 4   i 4   , k n1 i , c a n a 4   (3.20) Inserting Eq (3.20) into (3.7) results V ( )  i i  G'  ,    G  2 4   4     (3.21) Substituting the general solution of (2.9) for     into Eq (3.21), we get    4     4        C sin     C cos         2  1    , V ( )  1 i 2  4     4           C sin      C cos       2          (3.22) where  n1 i x  at n a 4     Using the following transformation,   i ,  4     4    sinh     i sin   ,     2      4     4    cosh     cos        2     (3.23) in Eq (3.22) and substituting the result into (3.3), we obtain the following exact solution of Eq (1.1): Exact Travelling Wave Solutions for Generalized Forms of the Nonlinear Heat Conduction Equation 119 1     n1   n1    n1 x  at   C i cosh   x  at        C sinh    2n a   2n a    u( x , t )      2    n1   n1       C cosh   n a x  at   C i sinh   n a x  at                        (3.24) We note that if we set C  and C  in the general solution (3.24), we can recover the solutions (3.12) and (3.13), respectively Case B-2 0  i  , 1  2 4   i 4   , k n1 i , c a n a 4   (3.25) Inserting Eq (3.25) into (3.7) leads to V ( )  i i  G'     , 2 4   4    G  (3.26) Substituting the general solution of (2.9) for     into Eq (3.26), we have    4     4        C sin     C cos         2  1    , V ( )  1i 2  4     4           C sin      C cos       2          in which    (3.27) n1 i x  at n a 4     Using the transformation (3.23) into Eq (3.27), and substituting the result into (3.3) yields the following exact solution: 1     n1   n1    n1 x  at   C i cosh   x  at        C sinh   1  2n a   2n a    u( x , t )     2    n1   n1    x  at   C i sinh   x  at        C cosh    2n a   2n a                 (3.28) Similarly, if we set C  and C  in the general solution (3.28), we arrive at the same solutions (3.18) and (3.19), respectively 3.2 Application of the Exp-function method In order to determine values of f and p , we balance the term v with vv in Eq (3.4); we have v3  c1 exp(3 f  )   , c exp(3 p )   (3.29) 120 Heat Conduction – Basic Research vv  c exp([2 f  p ] )   , c exp(5 p )   (3.30) where ci are determined coefficients only for simplicity Balancing the highest order of the Exp-function in Eqs (3.29) and (3.30), we have f  p  f  3p, (3.31) p f, (3.32) which leads to the result Similarly, to determine values of g and q , we have   d1 exp( 3 g ) ,   d2 exp( 3q ) (3.33)   d3 exp( [2 g  3q ] ) ,   d4 exp( 5 p ) (3.34) v3  vv  where di are determined coefficients for simplicity Balancing the lowest order of the Expfunction in Eqs (3.33) and (3.34), we have g  q  g  3q , (3.35) q  g (3.36) which leads to the result Case A: p  f  1, q  g  We can freely choose the values of p and q For simplicity, we set p  f  and q  g  , so Eq (2.11) reduces to v( )  a1 exp( )  a0  a1 exp(  ) , exp( )  b0  b1 exp(  ) (3.37) Substituting Eq (3.37) into Eq (3.4), and making use of Maple, we arrive at [c exp(4 )  c exp(3 )  c exp(2 )  c1 exp( )  c0  c 1 exp(  ) A  c 2 exp( 2 )  c 3 exp( 3 )  c 4 exp( 4 )]  0, (3.38) A  [exp( )  b0  b1 exp(  )]4 , (3.39) in which Exact Travelling Wave Solutions for Generalized Forms of the Nonlinear Heat Conduction Equation 121 And the cn are coefficients of exp( n ) Equating to zero the coefficients of all powers of exp( n ) yields a set of algebraic equations for a0 , b0 , a1 , a1 , b1 , k , and c Solving the system of algebraic equations with the aid of Maple 12, we obtain: Case a0  0, b0  0, a1  0, a1  b1 , b1  b1 , k   n1 , c a 2n a (3.40) Substituting Eq (3.40) into (3.37) and inserting the result into the transformation (3.3), we get the generalized solitary wave solution of Eq (1.1) as follows: 1   n1 b1 exp(  ) u( x , t )   ,  exp( )  b1 exp(  )   where    n1 x  at 2n a   (3.41) and b1 is an arbitrary parameter which can be determined by the initial and boundary conditions If we set b1  and b1  1 in (3.41), the solutions (3.18) and (3.19) can be recovered, respectively Case a0  0, b0  0, a1  1, a1  0, b1  b1 , k   n1 , c a 2n a (3.42) By the same procedure as illustrated above, we obtain 1   n1 exp( ) , u( x , t )     exp( )  b1 exp(  )  in which    (3.43) n1 x  at and b1 is a free parameter 2n a   If we set b1  and b1  1 in (3.43), then it can be easily converted to the same solutions (3.12) and (3.13), respectively Case a1  0, b1  0, a0  a0 , b0  b0 , a1  a0 b0 , k   and consequently we get n1 , c  n a n a (3.44) 122 Heat Conduction – Basic Research 1  a  a b exp(  )  n  u( x , t )   0 ,   exp( )  b0  where    n1 x  n at n a   (3.45) and a0 , b0 , are arbitrary parameters; for example, if we put b0  , solution (3.45) reduces to 1 u( x , t )   a0  cosh  sinh   n  ,   (3.46) a1  0, a0  a0 , b0  b0 , b1   a0 ( a0  b0 ), a1   a0 ( a0  b0 ), n1 , c a k n a (3.47) Case and 1   n1 a0  a0 ( a0  b0 )exp(  ) u( x , t )   ,   exp( )  b0  a0 ( a0  b0 )exp(  )  where    n1 x  at n a   (3.48) and a0 , b0 are free parameters; for example, if we set a0  1, b0  in Eq (3.48), it can be easily converted to 1 1  n1 u( x , t )   (1  coth  csc h ) , 2  (3.49) Case a1  1, a0  0, b0  b0 , b1  0, a1  0, k   n1 , c a n a (3.50) and finally we obtain 1  exp( )  n  u( x , t )     exp( )  b0  in which    n1 n a  x  at  and b Case B: p  f  2, q  g  is a free parameter (3.51) Exact Travelling Wave Solutions for Generalized Forms of the Nonlinear Heat Conduction Equation 123 Since the values of g and f can be freely chosen, we can put p  f  and q  g  , the trial function, Eq (2.11) becomes v( )  a2 exp(2 )  a1 exp( )  a0  a1 exp(  ) , exp(2 )  b1 exp( )  b0  b1 exp(  ) (3.52) By the same manipulation as illustrated above, we have the following sets of solutions: Case a1  0, a0  a0 , b0  0, b1  0, a1  0, a2  0, b1  0, k   n1 , c  n a 2n a (3.53) Substituting Eq (3.53) into (3.52), we have v( )  a0 exp( 2 ), (3.54) Substituting Eq (3.54) into Eq (3.3), we get the generalized solitary wave solution of Eq (1.1) as 1 u( x , t )  [ a0 exp( 2 )]n  , where    n1 2n a (3.55) ( x  n at ) and a0 is an arbitrary parameter Using the transformation exp( )  cosh  sinh , Eq (3.55) yields the same solution (3.46)  exp(  )  cosh  sinh Case a1  0, a0  b0 , b0  b0 , b1  0, a1  0, a2  0, b1  0, k   n1 , c a 2n a (3.56) Substituting Eq (3.56) into (3.52), we have v( )  b0 , exp(2 )  b0 (3.57) Inserting Eq (3.57) into (3.3), it admits to the generalized solitary wave solution of Eq (1.1) as follows: 1   n1 b0 u( x , t )   ,   exp(2 )  b0  where    n1 2n a ( x  at ) and b0 is a free parameter (3.58) 124 Heat Conduction – Basic Research We note that if we set a0  b0 in Eq (3.48), we can recover the solution (3.58) Case a1  0, a0  0, b0  0, b1  b1 , a1  b1 , a2  0, b1  0, k   n1 3n a , c a (3.59) Substituting Eq (3.59) into (3.52) we obtain v( )  b1 exp(  ) , exp(2 )  b1 exp(  ) (3.60) and by inserting Eq (3.60) into (3.3), we get the generalized solitary wave solution of (1.1) as 1   n1 b1 exp(  ) u( x , t )   ,   exp(2 )  b1 exp(  )  in which    (3.61) n1 ( x  at ) and b1 is a free parameter that can be determined by the 3n a initial and boundary conditions The generalized nonlinear heat conduction equation in two dimensions 4.1 Application of the (G'/G)-expansion method Using the wave variable (2.4) transforms Eq (1.2) to the ODE  kcU   ak (U n )  U  U n  0, a0 (4.1) or, equivalently,  kcU   ak n(n  1)U n  2U 2  ak nU n  1U   U  U n  0, (4.2) Then we use the transformation (3.3), which will convert Eq (4.2) into kc(n  1)V V  ak n(1  n)V 2  ak n(n  1)VV   ( n  1)2 V  ( n  1)2 V  0, (4.3) By the same manipulation as illustrated in Section 3.1, we obtain the following sets of solutions Case A: When     Case A-1 0  n1  1  , 1  , k , c   2a 2 2   4 n 2a   4   4 (4.4) Exact Travelling Wave Solutions for Generalized Forms of the Nonlinear Heat Conduction Equation 125 By the same procedure as illustrated in Case A-1 of Section 3.1, Eqs (3.9) and (3.10), we can finally find the generalized solitary wave solution of Eq (1.2) as     n1   n1 x  y  at   C cosh   x  y  at    C sinh    2n a   2n a 1 1   u( x , y , t )  2    n1   n1    C cosh   n a x  y  at   C sinh   n a x  y  at        1    n1       (4.5)                   in which C and C are arbitrary parameters that can be determined by the related initial and boundary conditions Now, to obtain some special cases of the above general solution, we set C  ; then (4.5) leads to 1 1  n1   n1 u( x , y , t )       x  y  at    ,  2n a   2    (4.6) and, when C  , the exact solution (4.5) reduces to 1 1  n1   n1 u( x , y , t )     coth   x  y  at    , 2  2n a      (4.7) Comparing the particular cases of our general solution, Eqs (4.6) and (4.7), with Wazwaz’s results (2005), Eqs (87) and (88), it can be seen that the results are exactly the same Case A-2 0   1 n1  , 1  , k , c   2a 2   4 n 2a   4   4 (4.8) By the similar process as illustrated in Case A-2 of Section 3.1, Eqs (3.15) and (3.16), we can easily gain the following exact solution of Eq (1.2):     n1   n1 x  y  at   C cosh   x  y  at    C sinh   1  2n a   2n a u( x , y , t )     2    n1   n1 x  y  at   C sinh   x  y  at    C cosh   2n a    2n a     1    n1       (4.9)                   Similarly, to derive some special cases of the above general solution, we choose C  ; then (4.9) leads to the formal solitary wave solution as follows: 1  n1 u( x , y , t )       x  y  at 2  2n a   1    n1 ,     (4.10) 126 Heat Conduction – Basic Research and, when C  , the general solution (4.9) reduces to 1 1  n1   n1 u( x , y , t )     coth   x  y  at    ,  2n a   2    (4.11) Validating our results, Eqs (4.10) and (4.11), with Wazwaz’s solutions (2005), Eqs (85) and (86), it can be seen that the results are exactly the same Case B: When     Case B-1 0  i  , 1  2 4   i 4   , k n1 i , c   2a n a 4   (4.12) By the same manipulation as illustrated in Case B-1 of Section 3.1, Eqs (3.21)-(3.23), we can finally obtain the following exact solution: 1     n1   n1 x  y  at   C i cosh   x  y  at    C sinh   1  2n a   2n a  1  u( x , y , t )  2    n1   n1 x  y  at   C i sinh   x  y  at    C cosh    2n a   2n a        n1       (4.13)                   We note that, if we set C  and C  in the general solution (4.13), we can recover the solutions (4.6) and (4.7), respectively Case B-2 0  i  , 1  2 4   i 4   , k n1 i , c   2a n 2a 4   (4.14) Similar to Case B-2 of Section 3.1, we can find the following result:     n1   n1 x  y  at   C i cosh   x  y  at    C sinh    2n a   2n a 1 1   u( x , y , t )  2    n1   n1    C cosh   2n a x  y  at   C i sinh   n a x  y  at        1    n1       (4.15)                  In particular, if we take C  and C  in the general solution (4.15), we arrive at the same solutions (4.10) and (4.11), respectively 4.2 Application of the Exp-function method By the same manipulation as illustrated in Section 3.2, we obtain the following sets of solutions Exact Travelling Wave Solutions for Generalized Forms of the Nonlinear Heat Conduction Equation 127 Case a1  a1 , a1  0, a0  0, b0  0, b1  a1 , k   n1 , c   2a 2n a (4.16) Substituting Eq (4.16) into (3.37) and inserting the result into the transformation (3.3), we get the generalized solitary wave solution of Eq (1.2) as follows: 1   n1 a1 exp(  ) u( x , y , t )   ,   exp( )  a1 exp(  )  where    n1 2n a x  y  at  and a1 (4.17) is an arbitrary parameter which can be determined by the initial and boundary conditions If we set a1  and a1  1 in (4.17), the solutions (4.10) and (4.11) can be recovered, respectively Case a0  0, b0  0, a1  1, a1  0, b1  b1 , k   n1 2n a , c   2a (4.18) By the same process as illustrated in the previous case, we obtain 1   n1 exp( ) u( x , y , t )   ,   exp( )  b1 exp(  )  in which    n1 2n a x  y  (4.19)  at and b1 is a free parameter If we set b1  and b1  1 in (4.19), then it can be easily converted to the same solutions (4.6) and (4.7), respectively Case a1  a1 , a1  0, a0  0, b1  0, b0  0, k   n1 2n a , c  n a (4.20) and consequently we get 1 1 u( x , y , t )   a1 exp( 2 ) n    a1  cosh 2  sinh 2   n  ,   where    n1 x  y  n at and a1 is an arbitrary parameter 2n a   (4.21) 128 Heat Conduction – Basic Research Case a1  1, a0  a0 , a1  0, b1  b1 , b0  b1  a0 n1 , k , c   2a a0 n 2a (4.22) and 1   n1   exp( )  a0  u( x , y , t )   ,   b1  a0  b1 exp(  )   exp( )  a0   where    (4.23) n1 x  y  at and a0 , b1 are free parameters n 2a   Case a1  0, a1  a1 , a0  a0 , b1  a1 , b0  a1  a0 n1 , k , c   2a a0 n 2a (4.24) and finally we obtain 1   n1   a0  a1 exp(  )  u( x , y , t )     a1  a0  a1 exp(  )   exp( )  a0   (4.25) n1 ( x  y  at ) and a0 , a1 are free parameters n 2a Remark We have verified all the obtained solutions by putting them back into the original equations (1.1) and (1.2) with the aid of Maple 12 Remark The solutions (3.12), (3.13), (3.18), (3.19), (4.6), (4.7), (4.10), (4.11) have been obtained by the method (Wazwaz, 2005); the other solutions are new and more general solutions for the generalized forms of the nonlinear heat conduction equation in which    Conclusions To sum up, the purpose of the study is to show that exact solutions of two generalized forms of the nonlinear heat conduction equation can be obtained by the (G'/G)-expansion and the Exp-function methods The final results from the proposed methods have been compared and verified with those obtained by the method New exact solutions, not obtained by the previously available methods, are also found It can be seen that the Exp-function method yields more general solutions in comparison with the other method Overall, the results reveal that the (G'/G)-expansion and the Exp-function methods are powerful mathematical tools to solve the nonlinear partial differential equations (NPDEs) in the terms Exact Travelling Wave Solutions for Generalized Forms of the Nonlinear Heat Conduction Equation 129 of accuracy and efficiency This is important, since systems of NPDEs have many applications in engineering References Abbasbandy, S (2010) Homotopy analysis method for the Kawahara equation Nonlinear Analysis: Real World Applications, 11, 1, 307-312 Bekir, A., Cevikel, C (2009) New exact travelling wave solutions of nonlinear physical models Chaos, Solitons and Fractals, 41, 1733–1739 Borhanifar, A., Kabir, MM (2009) New periodic and soliton solutions by application of Expfunction method for nonlinear evolution equations Journal of Computational & Applied Mathematics, 229, 158-167 Borhanifar, A., Kabir, MM., Vahdat Lasemi, M (2009) New periodic and soliton wave solutions for the generalized Zakharov system and (2+1)-dimensional Nizhnik– Novikov–Veselov system Chaos, Solitons & Fractals, 42, 1646–1654 Borhanifar, A., Kabir, MM (2010) Soliton and Periodic solutions for (3+1)-dimensional nonlinear evolution equations by Exp-function method Applications and Applied Mathematics: International Journal (AAM), 5, 1, 59-69 Fan, E (2002) Traveling wave solutions for nonlinear equations using symbolic computation Comput Math Appl., 43, 671–680 He, JH (1998) Approximate analytical solution for seepage flow with fractional derivatives in porous media Comput Methods Appl Mech Eng., 167, 57-68 He, JH (2000) A coupling method of a homotopy technique and a perturbation technique for non-linear problems Int J Non-Linear Mechanics, 35, 37-43 He, JH (2006) New interpretation of homotopy perturbation method Int J Mod Phys B, 20, 18, 2561–2568 He, JH., Wu, XH (2006) Exp-function method for nonlinear wave equations Chaos, Solitons & Fractals, 30, 3, 700–708 He, JH., Zhang, LN (2008) Generalized solitary solution and compacton-like solution of the Jaulent–Miodek equations using the Exp-function method Physics Letters A, 372, 1044–1047 He, JH (2009) Nonlinear science as a fluctuating research frontier Chaos, Solitons and Fractals, 41, 2533–2537 He, JH., Wu, GC., Austin, F (2010) The variational iteration method which should be followed Nonlinear Science Letters A, 1, 1, 1–30 Kabir, MM., Khajeh, A (2009) New explicit solutions for the Vakhnenko and a generalized form of the nonlinear heat conduction equations via Exp-function method International Journal of Nonlinear Sciences & Numerical Simulation, 10, 10, 1307-1318 Kabir, MM., Khajeh, A., Abdi Aghdam, E., Yousefi Koma, A (2011) Modified Kudryashov method for finding exact solitary wave solutions of higher-order nonlinear equations Mathematical Methods in the Applied Sciences, 34, 213-219 Kabir, MM., Borhanifar, A., Abazari, R (2011) Application of (G'/G)-expansion method to Regularized Long Wave (RLW) equation Computers and Mathematics with Applications, 61, 8, 2044–2047 130 Heat Conduction – Basic Research Wang, M., Li, X., Zhang, J (2008) The (G'/G)-expansion method and traveling wave solutions of nonlinear evolution equations in mathematical physics Phys Lett A, 372, 417–423 Wazwaz, AM (2004) A sine-cosine method for handling nonlinear wave equations Math Comput Model, 40, 499–508 Wazwaz, AM (2005) The method for generalized forms of nonlinear heat conduction and Burgers–Fisher equations Applied Mathematics and Computation, 169, 321–338 Wazwaz, AM (2006) New solitary wave solutions to the Kuramoto–Sivashinsky and the Kawahara equations Appl Math Comput., 182, 1642-1650 Zedan, HA (2010) New classes of solutions for a system of partial differential equations by G'/G)-expansion method Nonlinear Science Letters A, 1, 3, 219–238 Zhang, S., Wang, W., Tong, J (2009) A generalized (G'/G)-expansion method and its application to the (2+1)-dimensional Broer-Kaup equations Appl Math Comput., 209, 399-404 6 Heat Conduction Problems of Thermosensitive Solids under Complex Heat Exchange Roman M Kushnir and Vasyl S Popovych Pidstryhach Institute for Applied Problems of Mechanics and Mathematics, Ukrainian National Academy of Sciences Ukraine Introduction To provide efficient investigations for engineering problems related to heating/cooling process in solids, the effect of thermosensitivity (the material characteristics depend on the temperature) should be taken into consideration when solving the heat conductivity problems (Carslaw & Jaeger, 1959; Noda, 1986; Nowinski, 1962; Podstrihach & Kolyano, 1972) It is important to construct the solutions to the aforementioned heat conduction problems in analytical form This requirement is motivated, for instance, by the need to solve the thermoelasticity problems for thermosensitive bodies, for which the determined temperature is a kind of input data, and thus, is desired in analytical form In general, the model of a thermosensitive body leads to a nonlinear heat conductivity problem It is mentioned in (Carslaw & Jaeger, 1959) that the exact solutions of such problems can be determined when the temperature or heat flux is given on the surface by assuming the material to be “simply nonlinear” (thermal conductivity t and volumetric volumetric heat capacity c v depend on the temperature, but the relation, called thermal diffusivity a  t c v , is assumed to be constant) For construction of the solution in this case, it is sufficient to use the Kirchhoff’s transformation to obtain the corresponding linear problem for the Kirchhoff’s variable This problem can be solved (Ditkin & Prudnikov, 1975; Galitsyn & Zhukovskii, 1976; Sneddon, 1951) by application of classical methods (separation of variables, integral transformations, etc.) The solutions to the heat conductivity problems for crystal bodies, whose thermal characteristics are proportional to the third power of the absolute temperature, can be constructed in a similar manner for the case of radiation heat exchange with environment In the case of complex heat exchange, the Kirchhoff transform makes the heat conductivity problem to be linear only in part In the heat conductivity problem for the Kirchhoff’s variable, the heat conduction equation is nonlinear due to dependence of the thermal diffusivity on the Kirchhoff’s variable The boundary condition of the complex heat exchange is also nonlinear due to a nonlinear expression of the temperature on the surface Herein we discuss several approaches, developed by the authors for determining temperature distribution in thermosensitive bodies of classical shape under complex (convective, radiation or convective-radiation) heat exchange on the surface (Kushnir & Popovych, 2006, 2007, 2009; Kushnir & Protsiuk, 2009; Kushnir et al., 2001, 2008; Popovych, 132 Heat Conduction – Basic Research 1993a, 1993b; Popovych & Harmatiy, 1996, 1998; Popovych & Sulym, 2004; Popovych et al 2006) Note that the necessity of these investigations is emphasized in (Carslaw & Jaeger, 1959) The step-by-step linearization method for solving the one-dimensional transient heat conductivity problems with simple thermal non-linearity Let us consider the step-by-step method for determining one-dimensional transient temperature field t( x , ) , which can be found from the following non-linear heat conduction equation:   m t  t  x t (t )   c v (t )  W , m x  x   x (1) where t (t ) is the thermal conductivity; c v (t ) is the volumetric heat capacity; m  0; 1; corresponds to Cartesian, cylindrical and spherical coordinate systems, respectively; a  x  b , a  0, a  b   The thermosensitive body of consideration is made of a material with simple nonlinearity The density of heat sources W is a function of coordinate x and time  Let the surface x  a , for instance, is exposed to convective-radiation heat exchange t   0 t (t ) x   a (t )(t  t a )   a (t )(t  t a )  xa (2) with the environment of constant temperature ta , where  a (t ) is the temperature dependent coefficient of heat exchange between the surface and the environment;  a (t ) is the temperature dependent emittance;  is the Stefan-Boltzmann constant The surface x  b is heated with constant temperature tb or constant heat flux qb : t x  b  tb or t (t ) t x x b  qb (3) At the initial moment of time, the temperature is uniformly distributed within the body: t  0  (4) The key point of the solution method for the formulated non-linear heat conductivity problem (1)–(4), which is presented below, consists in the step-by-step linearization involving the Kirchhoff transformation along with linearization of the nonlinear term in the boundary conditions by means of the spline approximation By introducing the dimensionless coordinates x  x l0 , temperature T  t t0 , and time Fo  a l0 (the Fourier number), we can present the functional parameters t (t ) , c v (t ) ,  a (t ) , and  a (t ) in the form  (t )     (T ) , where  is a reference value and   (T ) stands for the dimensionless function; t0 is a reference temperature and l0 is a characteristic dimension The density of heat sources can be presented as W  q0 q( x ,Fo) , where q0 is the Heat Conduction Problems of Thermosensitive Solids under Complex Heat Exchange 133 dimensional constants, q( x ,Fo) is the dimensionless function describing the time variation of the heat sources As a result, the problem (1)–(4) takes the form   m  T   T  Po q( x ,Fo) ,  x t (T )   c v (T ) m x  x  Fo x (5) T     4  t (T ) x  Bi a a (T )(T  Ta )  Sk a a (T )(T  T a ) x  a  ,  (6) T x  b  Tb or t (T ) T x x b  Ki b , T Fo   Tp (8) (0) Here P0  q0l0 /(t0t ) (the Pomerantsev number), Bi a   a r0 t (0) Sk a   a l0t0 t (7) (the Biot number), Ki b  qb l0 /(t0 t ) (the Kirpichev number), (the Starc number), Tb  tb / t0 , Tp  t p t0 Let us apply the Kirchhoff’s integral transformation (Carslaw & Jaeger, 1959; Noda, 1986; Podstrihach & Kolyano, 1972)  T  t (t )dt * (9) Tp to the problem (5)–(8) By taking into account the feature of simple nonlinearity  T  T T ( t (T )  c  (T ) ) and expressions , the equation  t (T )  t (T )  c  (T ) , v v x x Fo Fo Fo   m     Po q( x ,Fo) x  x  Fo x m x  (10) follows from the nonlinear heat conductivity equation (5) The boundary condition of convective-radiation heat exchange (6) can be partially linearized and presented as    0,  x  Qa T ( )   x a (11)   where Qa T ( )  Bi a a T ( ) T ( )  Ta   Sk a a T ( ) (T ( )4  Ta4 ) The boundary condi- tions (7) and initial condition (8) yield  x b  b or  Fo   x x b 0,  Ki b , (12) (13) 134 where  b  Heat Conduction – Basic Research Tb   t (T )dT , T ( ) denotes the temperature expressed through the Kirchhoff’s Tp variable and determined for certain t (T ) by means of the integral equation (9) Application of the Kirchhoff’s variable allows us to linearize the nonlinear heat conductivity equation (5) and the second boundary condition (7) completely, whereas the convective-radiation heat exchange condition is linearized in a part Due to the nonlinear expression Qa T ( ) , it is impossible to apply any classical method to solve the boundary problem (10)–(13) Therefore, it is necessary to linearize the boundary condition (11) In (Nedoseka, 1988; Podstrihach & Kolyano, 1972), the convective heat exchange condition has been considered Therefore, the nonlinear expression T ( ) is simply replaced by  As a result, the nonlinear convective heat exchange condition on  becomes linear However, it has been shown in (Kushnir & Popovych, 2009; Popovych, 1993b; Popovych & Harmatiy, 1996) that this unsubstantiated linearization leads to the numerically or physically incorrect results In our case, when we take into account the radiation constituent (which is nonlinear even for a non-thermosensitive material) and dependence of the heat transfer coefficient and emittance on the temperature, the considered substitution does not provide the complete linearization of the condition (11) Instead, the boundary condition (11) can be linearized by means of interpolation of the nonlinear expression Qa T ( ) by special splines with order or For x  a , the expression Qa T ( ) is a function of Fo only Let us select a finite set of points Foi ( i  1, n ;  Fo0  Fo1  Fo    Fon ) , which divides the region of time variation into n  intervals Let us construct the spline S(0) (Fo) with order 0, whose values coincide a with the values of expression Qa (Fo)  Qa T ( )  at Fo  Foi and x a n1 ( a (0) Sa (Fo)  Q1a )   (Qi( )  Qi( a ) ) S (Fo  Foi ) ; (14) i 1    Qi( a )  Bi a a (Ti( a ) )(Ti( a )  Ta )  Sk a a (Ti( a ) ) (Ti( a ) )4  Ta4 on the every interval of interpolation Here Ti( a )  (15) (i  , n ) are the values of temperature T ( x ,Fo) , which are to be found on the surface x  a at the moments of time Foi (the unknown parameters of spline approximation), S () denotes the asymmetric unit Heaviside function (H Korn & T Korn, 1977) Having presented the nonlinear expression Qa T ( ) by spline (14), the boundary xa condition (11) becomes linear  x xa  S(0) (Fo)  a (16) Similarly, the first-order spline S(0) (Fo) , whose values coincide with values of expression a Qa (Fo) at the points Foi and on every segment of decomposition approximates Qa (Fo) by Heat Conduction Problems of Thermosensitive Solids under Complex Heat Exchange 135 the linear polynom Pi( a ) (Fo)  ki( a )Fo  bi( a ) , can be constructed by the abovementioned decomposition This spline can be written as n1   ( ) (1) Sa (Fo)  P1 a ) (Fo)   Pi(a1 (Fo)  Pi( a ) (Fo) S (Fo  Foi ) i 1 (17) Here the coefficients ki( a ) , bi( a ) of polynom Pi( a ) (Fo) are calculated by formulae ki( a )  a Qi( a )  Qi( ) a , bi( a )  Qi( )  ki( a )Fo i  , Foi  Fo i  (18) where Qi( a ) is expressed through Ti( a ) by means of formula (15) If Qa T ( ) x  a is expressed as the first-order spline (17), then boundary condition (11) becomes linear  x xa  S(1) (Fo)  a (19) Having solved the obtained linear problem (10), (12), (13), (16) or (10), (12), (13), (19) by means of the classical methods, the Kirchhoff’s variable is found as a function of x and Fo Besides the input data of the problem, this variable contains Foi and unknown values Ti( a )  T ( a ,Foi ) : (    ( x , Fo ,Fo1 ,  ,Fon , T1( a ) , , Tn a ) ) (20) By substitution  into the expression for T ( ) (for specific dependence t (T ) ), the formula for determination of the temperature ( T  f ( x , Fo ,Fo1 ,  ,Fon , T1( a ) , , Tn a ) (21) can be obtained at arbitrary point x and arbitrary moment of time Fo For determination of unknown values Ti( a ) in the expressions for temperature (21), the collocation method is used Assuming Fo  Foi (i  1, n ) in (21), the system of equation for determination Ti( a )  ( a) T1  f ( a , Fo1 , T1 ) ,  ( a) ( a) T2  f ( a , Fo1 , Fo , T1 , T2 ),   T  f ( a ,Fo , ,Fo , T ( a ) , , T ( a ) ) n n  n (22) is obtained The structure of system (22) makes it possible to determine all unknown values ( Ti( a ) , starting from T1 a ) Substitution of values, determined from (22), into the formula (21) completes the solution procedure The temperature at given point x and moment of time can be calculated in accordance to the following scheme: 136 a Heat Conduction – Basic Research to divide the time axis by Foi and then to determine the approximation parameters Ti( a ) ( from the system (22); as a result, the value of temperature (21) Tn a ) is obtained; b to divide every interval in two; to compute the values of parameters Ti( a ) for this new (a time-segmentation and then to obtain the values of temperature Tn )1 ; c (a ( to calculate the difference Tn )1  Tn a ) If Tn   Tn   , where  is the accuracy, then the calculation is over Otherwise, we shall return to the stage b The temperature can be computed with any given accuracy  for arbitrary segmentation of the time axis However, the increasing of number of time-segments decreases the convergence of the proposed scheme An appropriate choice of the initial moment of time can be done by means of the estimated ‘a priory’ time-dependence of the temperature on the surface x  a We can also use the solution of corresponding boundary value problem for the body of the same shape with constant characteristics Then the initial choice for values Foi can be used as the appropriate one for the thermosensitive body The method of step-by-step linearization is applicable for determination of the temperature fields in thermosensitive plates, half-space, solid and hollow cylinders or spheres, space with cylindrical or spherical cavities, on the surfaces of which, the conditions of convective, radiation or convective-radiation heat exchange may be given This method has been efficiently used for solving the two-dimensional steady problem in thermosensitive body Method of linearizing parameters The method of step-by-step linearization makes it possible to determine the solutions to the two-dimensional heat conductivity problems in thermosensitive bodies with simple nonlinearity, when the nonlinear term in the condition of complex heat exchange for the Kirchhoff’s variable depends on one (spatial or time) variable only In this section, we consider an efficient method for solving the steady-state and transient heat conductivity problems of arbitrary dimension those describe the propagation of heat in thermosensitive bodies with simple nonlinearity under the convective heat exchange with environment Let the body occupies region D with surface S The surface (whole or a part) is subjected to the convective heat exchange with the environment of temperature t p From the moment of time   , the heat sources W ( x , y , z , ) are acting in the body The temperature in the body shall be determined from the following heat conduction equation: div  t (t )grad t   c v (t ) t W  (23) and the boundary t   t (t ) n   (t  tc )   s (24) and initial t  0  (25) Heat Conduction Problems of Thermosensitive Solids under Complex Heat Exchange 137 conditions, where  is the constant heat transfer coefficient; n is the external normal to surface S By making use of the above-introduced presentation for the material characteristics, heat sources, and dimensionless variables, the boundary value problem (23)–(25) can be reduced to the dimensionless form After application of the Kirchhoff’s transformation, the following boundary value problem for variable  divgrad    Po q( X , Y , Z ,Fo) , Fo     n  Bi T ( )  Tc    ,  s  Fo  0 (26) (27) (28) X  x l0 , Y  y l0 , Z  z l0 are dimensionless coordinates; is obtained, where n  n l0 , q ( X , Y , Z ,Fo) is the dimensionless function of heat sources As a result, the initial problem is partially linearized, meanwhile the condition (27) remains nonlinear The latter conditions have been obtained from the conditions of convective heat exchange due to nonlinear expression T ( ) on the surface S For solving the problem (26)–(28) by using an analytical method, it is necessary to linearize this condition Let us prove the possibility of such linearization Consider the simplest case of linear dependence of heat conductivity coefficient on the temperature: t (t )  to t (T )  to 1  k(T  Tp ) ,   (29) where k is a constant From the equation (9), the formula k   (T  Tp )  (T  Tp )2 (30) T ( )  k 1 (  k  1)  Tp (31) follows, where From the physical standpoint, the square root is chosen to be positive After substitution of the equation (31) into the boundary condition (27), the last one takes the form     k    Bi   Tp  Tc       k  n  s   (32) Be decomposing the square root in (32) into the series and restricting this series with two terms, the boundary condition     n  Bi   (Tc  Tp )    s   (33) 138 Heat Conduction – Basic Research is obtained The solution of equation (26) with boundary conditions (28), (33) is an approximate solution to the boundary value problem (26), (28), (32) To determine the exact solution, the equation (26) is to be solved under initial condition (28) and the following linear boundary condition     n  Bi (1   )  (Tc  Tp )    s   (34) instead of the nonlinear condition (32), where  is an unknown constant (linearized parameter) Note that the boundary condition (34) coincides at   with the condition (33) Since the problem (26), (28), (34) is linear, the appropriate classical analytical method can be used for its solution In addition to the original parameters of the problem ( Po,Bi, Tc , Tp , dimensions of the body, coordinates and time), the solution involves the unknown linearized parameter  :    ( X , Y , Z ,Fo,  ) (35) For an arbitrary value of  , the solution (35) meets the equation (26) and the initial condition (28) In order the solution (35) to satisfy the nonlinear conditions (32) and (34), the parameter  is to be the solution of the equation   k    (1   )    k  s   After some transformations, this equation can be given as  s  2 k(1   )2 (36) This equation holds for every moment of time Fo After the paramenter  is found, we substitute it into (35) In such manner, the expression for Kirchhoff’s variable is obtained The temperature in the body is then calculated by means of the relation (31) Note that the boundary condition (34) can be represented as       n  Bi (  Tc )  ,  s (37) where Bi  Bi(1   ); Tc  (Tc  Tp ) (1   ) This condition can be interpreted as a condition of convective heat exchange with certain parameters (the Biot number Bi and the temperature Tc of external environment) depending on the unknown parameter  The equation (36) is nonlinear It provides analytical solutions only for some cases of steadystate problems with substantional use of the numerical methods Therefore, these solutions can be regarded as analytico-numerical solutions Let us consider the non-linear dependence of the heat conductivity coefficient on the temperature For linearization of the boundary condition (27), we shall find the Kirhoff’s variable for the case when the surface temperature of the thermosensitive body is equal to ... convective-radiation) heat exchange on the surface (Kushnir & Popovych, 20 06, 2007, 2009; Kushnir & Protsiuk, 2009; Kushnir et al., 2001, 2008; Popovych, 132 Heat Conduction – Basic Research 1993a,...  (33) 138 Heat Conduction – Basic Research is obtained The solution of equation ( 26) with boundary conditions (28), (33) is an approximate solution to the boundary value problem ( 26) , (28), (32)... Regularized Long Wave (RLW) equation Computers and Mathematics with Applications, 61 , 8, 2044–2047 130 Heat Conduction – Basic Research Wang, M., Li, X., Zhang, J (2008) The (G''/G)-expansion method and