1549301742-The_Theory_of_Difference_Schemes__Samarskii

(jair2018) #1
416 Stability Theory of Difference Schemes


  1. Stability of the weighted scheme. As an example of applying the
    theorems just established the weighted scheme comes first
    ( 46) Yt + A ( O"Y + ( 1 - O") y) = 'P, y(O) =Yo.
    Recall that in Section 1 scheme ( 46) has been already reduced to the
    canonical fonn
    ( 47) ( E + () T A) Yt + Ay = 'P J y(O) =Yo ·
    By comparing (46) with (1) we see that B
    exists an inverse operator A-^1. Applying A-^1
    canonical form of the weighted scheme:


E + O"TA. Suppose there
to ( 4 7) reveals the second




    • By 1 + Ay = <p, y(O) = Yo ,
      (48)
      B=A-^1 +0"rE,
      In the sequel the form ( 47) is more convenient for the case of self-adjoint
      operators A and the form ( 48) - for the case of non-self-adjoint positive
      definite operators A = A( t). We will elaborate on this later.
      The first analysis is connected with the case when A is a constant
      self-adjoint positive operator A= A* > 0. As we have shown in Section 2,
      a necessary and sufficient condition for the stability of the weighted scheme
      (47) with respect to the initial data is
      1 1
      O"o = - - ---
      2 rll All.
      Under this con di ti on estimate ( 15) holds true for a solution of problem
      (47). In particular, for an explicit scheme (for O" = 0) the condition O" > 0" 0
      implies T < 2/llAll, that is, an explicit scheme is stable in the space HA
      for T < 2/llAll· A sche1ne with T >~is unconditionally stable, that is,
      for any T. In Section 2.4 a model exan1ple with Ay = -Ay = -Yxx for
      0
      y E Dh = Hh has been considered in full details, in which II A II< 4/h^2 and
      the appropriate explicit scheme is stable for T < ~h^2.
      For the heat conduction equation with a variable coefficient k( x) we
      might have




o ( ou)
Lu = ox k( x) ox '

A y = (a( x) Yx) x ,


and 0" 0 < ~ - h^2 /(4c 2 r), while an explicit scheme is stable for T < ~h^2 /c 2.

Free download pdf