Solution Second – Order Partial Differential Equation by Using Arbitrary Condition

DOI : 10.17577/IJERTV3IS030326

Download Full-Text PDF Cite this Publication

Text Only Version

Solution Second – Order Partial Differential Equation by Using Arbitrary Condition

Hayder Jabbar Abood (1) , Mushtaq Shakir Al-Shehani (2) , Oday Hatem Jalel (3)

(1) Department of Mathematics , College of Education for Pure Sciences , Babylon University , Babylon , Iraq .

(2)Vice Chancellor For Administrative Affairs , University of Babylon, Babylon , Iraq .

(3) Department of Mathematics , College of Education for Pure Sciences , Babylon University , Babylon , Iraq .

Abstract : We study the development of Burgers' equation by using the function depending on the time (t), and we prove the existence and uniqueness solution , also given some applications .

  1. INTRODUCTION

    In [1] consider the problem of stabilization of the inviscid Burgers partial differential equation (PDE) using boundary actuation ,and a solution to the problem using a Lyapunov approach and proved that the inviscid Burgers equation is stabilizable around a constant uniform state under an appropriate boundary control . In [2], In this paper we study the generalized Burgers equation ut +(u2/2)x = f (t)uxx , where f (t)>0 for t > 0 . The existence and uniqueness of classical solutions to the initial value problem of the generalized Burgers equation with rough initial data belonging to L (R), as well it is obtained the decay rates of u in Lp norm are algebra order for p [1,]. Burgers equation with uncertain initial and boundary conditions is approximated using a polynomial chaos expansion approach where the solution is represented as a series of stochastic, orthogonal polynomials, even though the analytical solution is smooth, a number of discontinuities emerge in the truncated system , in [3] . In [4] , a method for the solution of Burgers' equation is described ,the marker method relies on the definition of a convective field associated with the underlying partial differential equation ; the information about the approximate solution is associated with the response of an ensemble of markers to this convective field . An analysis of dispersive / dissipative features of the difference schemes used for simulations of the non-linear Burgers equation is developed based on the travelling wave asymptotic solutions of its differential approximation , It is showed that these particular solutions describe well deviations in the shock profile even outside the formal applicability of the asymptotic expansions, namely for shocks of moderate amplitudes , in [5] .In [6] , concerned with an interesting numerical anomaly associated with steady state solutions for the viscous Burgers' equation , also we considered Burgers' equation on the interval (0, 1) with Neumann boundary conditions , they showed that even for moderate values of the viscosity and for certain initial conditions, numerical solutions approach non – constant

    shock type stationary solutions .In [7] , consider Burgers equation with a time delay , by using the Liapunov function method , they showed that the delayed Burgers equation is exponentially stable if the delay parameter is sufficiently small , also give an explicit estimate of the delay parameter in terms of the viscosity and initial conditions, which indicates that the delay parameter tends to zero if the initial states tend to infinity or the viscosity tends to zero . The 1D Burgers equation is used as a toy model to mimick the resulting behaviour of numerical schemes when replacing a conservation law by a form which is equivalent for smooth solutions, such as the total energy by the internal energy balance in the Euler equations , if the initial Burgers equation is replaced by a balance equation for one of its entropies (the square of the unknown) and discredited by a standard scheme , the numerical solution converges, as expected, to a function which is not a weak solution to the initial problem , in [8] . In [9] , they proved that the temperature distribution in the limit one dimensional rod with time averaged sources of heat is the uniform asymptotic approximation of the temperature distribution in the initial problem in an arbitrary sub domain of the plane rod and in an arbitrary time interval , which are located at a positive distance from the ends of the rod and the initial time instance , respectively ; of course ,the temperature in the one dimensional rod , which is a function of the longitudinal coordinate x and the time t , is identified with the function of ( x , y ,t) , which is independent of the transversal coordinate y of the plane rod .In [10] , they are obtained an asymptotic expansion, containing regular boundary corner functions in the small parameter, for the solution of a second order partial differential equation, they are constructed the asymptotic expansion un(x, t, ) for the modified problem and prove it is the unique solution, also they have proved that the solution is valid uniformly in the domain, and the asymptotic approximation is within O(n+1)

    . In [11], they are studied the development wave equation with some conditions and proving the existence and uniqueness solution by using the reflection method .In [12], they are studied an modification of an initial boundary value problem in the critical case for the heat conduction equation in a thin domain , they are justify asymptotic

    expansions of the solution of the problems with respect to a small parameter >0 , also proved that the solution is uniform in the domain and the asymptotic approximation is within O(n+1). In [13], constructed asymptotic first order solution of a partial differential equation with small parameter , and have proven the solution is unique and

    uniform in the domain and, further , the asymptotic approximation is within. In this work , we proving the existence and uniqueness of second order partial differential function by using arbitrary conditions .

  2. STATEMENT OF THE PROBLEM :

    Consider the second order partial differential equation definitions of the following form :

    , + () , = 2 ,

    (1.2)

    2

    , > 0 , , 0 = ,

    where g(t)and f(x) are continuous functions define in , ={(, ),

    , > 0} also we suppose the following condition is valid , for given

    1 < 2

    =

    , = 1,2 ( 2.2)

    0

    , ~0 1 , > , ~0 2 ,

    < , 0 .

    We can change the equation by use the new variables

    , = 2 (,)

    ,

    (3.2)

    we have

    ,

    =

    2

    2

    = ()

    2

    + ()

    2

    2

    2

    ,

    =

    2

    2

    +

    2 2 2 2 2

    4

    =

    2

    +

    2

    2

    + 2

    2

    2

    2

    2

    2

    3

    =

    2

    + 3

    2

    2

    2

    2

    2

    3

    , = 2 = 2 + 2

    2

    2

    Compensate derivatives of x and t in the equation (1.2) , we get :

    2

    +

    2

    2

    2

    ()

    2

    2 =

    2

    2

    + 3

    2

    2

    2

    2

    2

    3

    2

    2

    +

    2

    2 + 2

    2

    2

    3 3

    3

    3 +

    2

    2 + 3 = 0

    ()

    2

    3

    2

    3

    2

    + ()

    2

    2

    2

    + 2

    2

    3

    3

    2 2 3 = 0

    2

    3 + + 2 1

    = 0

    3 + + 2 1

    = 0

    +

    = 0

    [ ] +

    = 0

    If we suppose u =u(x , t) , then we solving the equation (1.2), is :

    = 2 ,

    ,

    = 2 1

    () = 2 1

    = 1

    0

    1

    ()

    2

    1

    1 =

    2

    1 = 2 0

    2

    0

    and

    1 , 0 = = () = 1

    ,

    1 1 1 ( )2

    , =

    2 0

    2

    ,

    where

    = 1

    2

    1 (, ,)

    2

    ,

    0

    ( )2

    and

    , ,

    = + ,

    1 ,,

    2

    , =

    () (4.2)

    which agrees with what was found by the first approach for solving the Development Burgers equation. from (4.2)

    we get

    0

    , ~ =

    Can rewrite the solution to be asymptotically

    = 0

    = + 0 = 0 =

    Since = then we get :

    = 0 (5.2)

    The equation (5.2) is exactly the solution of the problem

    , + () , = 0 (6.2)

    , > 0 , , 0 = , .

    we found the equation (6.2) by the method of characteristics, the solution of equation (5.2) is smooth for t, if the function u0 (x) is smooth function and differentiate for x , then we have

    0

    = 1 ,

    where

    =

    ()

    0

    0

    1 + ,

    0

    1 + 0 ,

    If suppose u' (x) < 0 for every x , then ux

    = , if = 1

    ( )

    0

    is the first instant TO

    when ux

    = , known as gradient

    catastrophe, corresponds to a So where u'0(x) has a minimum

    = 1

    , " = 0

    0

    0 (0) 0 0

  3. WE CAN STUDY SOME APPLICATION ABOUT EQUATION (1.2):

    1. Find the instant of gradient catastrophe for the problem

      , + , =

      , > 0 , , 0 = 2 + 1 ,

      Solution: The solution of the problem

      For the function 0 = 2 + 1

      , = 2 , + 1

      0

      = 1 = 0

      Because 2 + 1 = 0 + 1 = 1

      Note that

      The graphs of the functions 2 + 1

      = (1)

      0

      0

      = (min 1) = 1

      lim( 2 + 1) = 1 + ()

      0

      = 0 < 0

      1 > 0

      Figure (1.3)

      Plot3D[1 + 2, {, 3.4,5.2}, {, 3.5,4.2}

    2. Solve the initial value problem foe u( x , t) at time t > 0in terms of t and a characteristic variable

      , + , = with

      u(x,0) =u (x)= 1 1 , do characteristic cross for any t >0 and if so where and when ?

      0 2

      Solution:

      , + , = 1 , = =

      , = 1

      1 1

      = + = , = = ()

      = = +

      = + = + = 0 =

      1

      =

      = 0

      1

      = 1 2 tanh

      1

      = 1 2 tanh + = [1 2 tanh ] +

      1

      0 = = 1 2 2

      2

      that is when = 2

      = + 2tanh

      2

      Figure (2.3)

      Plot3D[1 tanhx , {tanhx, 8,8}, {, 2.5,3.4}

      2

  4. PROCEDURE OF SOLVING THE PROBLEM : Consider the following problem :

    (,) + (,) = (1.4)

    , > 0 , , 0 = = 0 ,

    Which is a limit case of Development Burgers' equation as 0 .

    Definition 1.4 :

    Assume 0 1 , , 2 × 0, 1.2

    + +

    0 , 0 = 0

    0

    0

    1 × 0, .

    Now we need the following Proposition:

    Proposition 1.4 : [1]

    Let u (x, t) be a smooth solution of the problem , + , = 0, , > 0 .

    , 0 = 0 , . Then u (x ,t ) is a weak solution of the problem.

    Proof:

    Obviously , 1 × 0, . , 1 ( × 0, )

    0

    and , × 0, Multiplying + = 0 integrating in , × 0, and using

    ±, = (, ) = 0 we obtain :

    0 = , + ,

    0

    = (() + () +

    0

    = () +

    0 0

    = ( 0 , 0 + +

    Theorem 1.4

    0

    For a given function 2,1 the regularizations tend to u in 2,1 k for every compact

    , i ,e 2,1 0 0 . [1]

    For the theorem (1.4) and proposition (1.4) the following theorem is holes : Theorem 2.4

    Let 1 × 0, be a smooth solution of the equation + = 0 and a weak solution of the

    problem (1.2). If 0 is continuous at a point x0 , then u(x0 , 0) = u0(x) .

    Proof :

    0

    Let , 1( × 0, ) . As in proposition (1.4) we are led to

    , 0 0 , 0 = 0

    Suppose u(x0,0) > u0(x0) . By continuity there exist a neighborhood U such that

    , 0 > 0 ,

    0

    Take , 1 × 0, such that

    , 0 = , , > 0 , , ,

    = = 0

    Then

    , 0 0 , 0

    = , 0 0 , 0 > 0

    Which is a contradiction , similarly u(x0 , 0) < u0(x0) is impossible . then u(x0 , 0) = u0(x0)

    Now let us consider the initial data

    > 0

    0 = < 0

    where ul and u are constant .

    (2.4)

    The two cases > < are quite different with respect to the solvability of problem(1.4) . It can be proved that if > then the weak solution is unique, while if < then there exist infinitely many solutions.

  5. STUDYING THE FOLLOWING CASES FOR PROBLEM (1.4) :

    Case I : > then we have the following problem

    , + () , = ,

    , > 0 , , 0 = = 0() ,

    If > we are in a situation to apply the theorem (1.4) , Let x > 0 be fixed and

    = +

    2

    The instant of the theorem (2.4)and by using the condition (2.2) is determined by the slope k of the straight line

    through the points (x,0) and (0,m)

    then

    1

    = =

    ,

    = ,

    and

    >

    , ~ <

    0

    The unique solution of problem (1.4) is known as a shock wave , while = ( + )/ is a shock speed the speed at which the discontinuity of the solution travels .

    Theorem 1.5 :

    >

    Let the function , = &lt

    (1.5)

    2

    Is a weak solution of the problem (1.4) with initial data (2.4) , where > , = +

    Proof :

    0

    Let , 1 × 0, . Denote for simplicity

    = +

    0

    We have

    = , 0 , 0

    = + + +

    0 0

    1 = +

    0

    = ( ) + ( )

    0

    by

    0

    , =

    , , ,

    , = ( , ) , ) ,

    0

    0

    0

    (

    , ) = , , 0

    = , ,

    then

    0

    0

    , = ( , 0) , .

    And by

    0

    0

    then

    , = , , = (, ) ,

    0

    1 = ( , 0 + , ) + () , .

    Similarly

    0 0

    2 = + ,

    0

    = ( ) + ( ) ,

    0

    0

    = ( , 0 , ) () , .

    Because

    0 0 0

    , =

    , + , ,

    and

    then A = A1 + A2

    , = , , = (, ) ,

    0

    ( , 0 + , ) + () ,

    0 0

    ( , 0 , ) + () ,

    0 0 0

    ( , 0 + , ) ( , 0

    0 0 0

    , + () , + () ,

    0 0 0

    , 0 , , 0

    0 0 0

    + , + , + ,

    0 0 0

    , 0 , 0

    ,

    + , .

    0 0 0 0

    , 0 + , ( )

    0 0

    + ,

    0

    + ,

    = , 0 + , [ + ] = 0 .

    0 0

    On the other hand

    = , 0 , 0 ,

    0

    = , 0 , 0

    , 0 , 0 ,

    0

    0

    = , 0 , 0 = 0 .

    Then we obtain

    0

    A = B

    Case II : If ul < u then we have the following problem :

    In this case there exist more than one weak solution, we show the following theorem is valid .

    Theorem 2.5

    <

    Let the function , =

    >

    is a weak solution of the problem (1.4)with initial data (2.4) .

    Proof : Let , × 0, for simplicity we take ul = -1 and u = 1 and denote

    = +

    0

    0

    = , 0 , 0 = , 0 , 0 .

    Where ul = -1 and u = 1

    0

    the function satisfies the equation (,) + (,) = 0

    we have

    = (, ) + (, )

    0

    + ( (, ) + (, ))

    0

    + (, ) + (, ) = 1 + 2 + 3 ,

    where

    0

    1 = (, ) + (, ) ,

    0

    = (

    0

    ( , ) , (1)) ,

    by

    0 0

    (

    , ) = , , 0 = , 0

    and

    0

    , = ,

    0

    by

    0

    then

    , = , ( , ) = , ,

    , + , ,

    0

    0

    = , 0 , , . .

    0 0

    2 = ( (, ) + (, )) ,

    Y2 has a singularity at 0

    where

    0

    0

    2 = lim 2, ,

    2 = ( (, ) + (, )) ,

    = ( , ) + ( (, ))

    we have

    ( ) =

    () + , + , ,

    and

    () = , + , ,

    so

    ( (, ))

    = , + , ,

    we take the limit as 0

    then

    lim

    0

    , = 0

    2 = , + , + , + , .

    0 0

    Similarly

    3 = , + ,

    0

    = , 0 , , .

    Finally

    0 0 0

    Y = Y1 + Y2 + Y3

    0

    = , 0 , ,

    0 0

    + , + , + , + ,

    0 0

    0

    , 0 , , = , 0 , 0 =

    0 0 0

    The proof is complete .

    0

  6. WE HAVE SOME APPLICATIONS ABOUT THE CASE II :

1.6 : we solve the initial value problem for u(x,t) , t >0 in terms of t and a characteristic variable

, + , = 0 with

2 , 1

, 0 = 0 = 22 , 1 1

2 , 1

= ()

, =

=

= 0 =

() 0

= () = () +

hence in terms of t and L , u = B(L) with = () + = 0 = ,

and u = B(L) = 0

= 2 1 , 2 1 ,

2

2

so = 2 1 , 2 1 ,

2

with = + × 2 1 ,

2 1 ,

we have ,

0 = = 4 + ×

0 > 1 ,

4 0 < < 1 ,

4 1 < < 0 .

That is when at t = 0 1 < < 0 = + () = 0

Figure (1.6)

Plot 3D[22, {, 20.5,2.5}, {, 5.4,12}]

REFERENCES :

    1. Ioannis P. Stavroukis and stepan A. Tersian , Partial Differential Equations , world scientific, pp135 , 2003 .

    2. Jinghua Wang , Hui Zhang , Existence and decay rates of solutions to the generalized Burgers equation , Journal of mathematical analysis and applications , pp 213235 , 2003.

    3. Per Pettersson , Qaisar Abbasb, Gianluca Iaccarino and Jan Nordstr¨omb , Efficiency of shock capturing schemes for Burgers equation with boundary uncertainty , Seventh South African Conference on Computational and Applied Mechanics , 2010 .

    4. Jerome L. V. Lewandowski , Solution of Burgers, equation using the marker method , International journal of numerical analysis and modeling , Volume 3 , Number 1 , pp 80-93 , 2006

      .

    5. A.V. Porubov , D. Bouche and G.Bonnaud , Description of numerical shock profiles of non-linear Burgers equation by asymptotic solution of its differential approximations , 2003 .

    6. J. Burnsy , A. Balogh , D.S. Gilliamz and V.I. Shubov , Numerical Stationary Solutions for a Viscous Burgers' Equation

      , Journal of Mathematical Systems , Estimation, and Control , Volume 8 , number 2 , pp 1-16 1998 .

    7. Weijiu Liu , Asymptotic behavior of solutions of time delayed Burgers' equation , Discrete and continuous dynamical system , Volume 2 , Number 1, 2002 .

    8. T. Gallou¨et , R. Herbin , J.-C. Latch´e and T.T. Nguyen , Playing with Burgerss equation , 2011 .

    9. <>V. B. Levenshtam and H. D. Abood , Asymptotic Integration of The Problem on Heat Distribution in A Thin Rod Wirth Rapidly Varying Sources of Heat , Journal of Mathematical Science . Volume 129 , Number 1 , 2005 .

    10. Hayder Jabber Abood and Ali Abbas Jabir , Asymptotic Approximation of the Second Order Partial Differential Equation by Using Many Functions , British Journal of Science

      , Volume 5 , 2012.

    11. Hayder Jabber Abood and Ahmed Hadi Hussain , The Existence Solution to the Development Wave Equation With Arbitrary Conditions , Mathematical Theory and Modeling , Volume 3 , Number 4 , 2013.

    12. Hayder Jabber Abood, On The Formation of an Asymptotic Solution of the Heat Equation With Small Parameter, European Journal of Scientific Research , Volume 56, Number 4 ,pp471- 481 , 2011 .

    13. Hayder Jabber Abood, Estimation of Partial Differential Equations Depending on a Small Parameter , American Journal of Scientific Research , pp32-39 , 2011

Leave a Reply