5 Εξίσωση θερμότητας: Πεπερασμένες διαφορές

5.6 Ασκήσεις

  1. 5.1.

    Θεωρούμε τη μέθοδο Uij+1=λUi+1j+(1-2λ-βk)Uij+λUi-1j για την επίλυση της εξίσωσης της θερμότητας (5.1).

    1. (a)

      Είναι η μέθοδος άμεση ή πεπλεγμένη;

    2. (b)

      Είναι η μέθοδος ευσταθής;

    3. (c)

      Για ποιες τιμές του β, αν υπάρχουν, είναι η μέθοδος συνεπής;

  2. 5.2.

    Θεωρούμε την εξίσωση της θερμότητας (5.1) όπου έχουμε τροποποιήσει τις συνοριακές συνθήκες ως εξής, ux(0,t)=α, και u(1,t)=0.

    1. (a)

      Γράψτε ένα πεπλεγμένο αριθμητικό σχήμα όπου το αντίστοιχο με την (5.28) σφάλμα διακριτοποίησης ηij ικανοποιεί μια σχέση όπως αυτή της (5.29), δηλαδή

      |ηij|C(k+h2),

      με C ανεξάρτητη των k και h. Στη συνέχεια, εκφράστε τη μέθοδο σε μορφή γραμμικού συστήματος, όπως η (5.33).

    2. (b)

      Γράψτε ένα πεπλεγμένο αριθμητικό σχήμα όπου το σφάλμα διακριτοποίησης ηij ικανοποιεί μια αντίστοιχη σχέση όπως αυτή της (5.44), δηλαδή

      |ηij|C(k2+h2),

      με C ανεξάρτητη των k και h. Επίσης, εκφράστε τη μέθοδο σε μορφή γραμμικού συστήματος.

  3. 5.3.

    Έστω N,M φυσικοί αριθμοί, h=L/(N+1), xi=ih, i=0,,N, μια διαμέριση του [0,L], καθώς και k=T/M, tj=jk, j=0,,M, μια διαμέριση του [0,T], με L,T>0. Θεωρούμε ένα αριθμητικό σχήμα που χρησιμοποιεί τις τιμές στα (xi-1,tj+1), (xi,tj+1), (xi+1,tj+1), (xi-1,tj) και (xi+1,tj) για την προσέγγιση της λύσης της εξίσωσης της θερμότητας (5.1) στο σημείο (xi,tj).

    1. (a)

      Δείξτε την ύπαρξη ένος τέτοιου αριθμητικού σχήματος, καθώς και την αντίστοιχη εκτίμηση του σφάλματος διακριτοποίησης ηij.

    2. (b)

      Είναι η μέθοδος συνεπής;

    3. (c)

      Είναι ευσταθής;

  4. 5.4.

    Έστω N,M φυσικοί αριθμοί, h=L/(N+1), xi=ih, i=0,,N, μια διαμέριση του [0,L], καθώς και k=T/M, tj=jk, j=0,,M, μια διαμέριση του [0,T], με L,T>0. Θεωρούμε ένα αριθμητικό σχήμα που χρησιμοποιεί τις τιμές στα (xi-1,tj+1), (xi+1,tj+1), (xi-1,tj), (xi-1,tj-1) και (xi+1,tj-1) για την προσέγγιση της λύσης της εξίσωσης της θερμότητας (5.1) στο σημείο (xi,tj).

    1. (a)

      Δείξτε την ύπαρξη ένος τέτοιου αριθμητικού σχήματος καθώς και την αντίστοιχη εκτίμηση του σφάλματος διακριτοποίησης ηij.

    2. (b)

      Είναι η μέθοδος συνεπής;

    3. (c)

      Είναι ευσταθής;

  5. 5.5.

    Θεωρούμε το πρόβλημα

    ut(x,t) =Duxx(x,t)-βu(x,t),(x,t)[0,L]×[0,T],
    u(x,0) =g(x),x[0,L],
    u(0,t) =u(L,t)=0,t[0,T],

    όπου D,β είναι θετικές σταθερές.

    1. (a)

      Γράψτε ένα πεπλεγμένο αριθμητικό σχήμα με σφάλμα διακριτοποίησης ηij το οποίο να ικανοποιεί μια σχέση της μορφής

      |ηij|C(k+h2),

      με C ανεξάρτητη των k και h. Στη συνέχεια, εκφράστε τη μέθοδο σε μορφή γραμμικού συστήματος.

    2. (b)

      Είναι αυτή η μέθοδος ευσταθής; Είναι ευσταθής υπό συνθήκες;

  6. 5.6.

    Θεωρούμε το πρόβλημα

    ut(x,t) =Duxx(x,t)-αux(x,t),(x,t)[0,L]×[0,T],
    u(x,0) =g(x),x[0,L],
    u(0,t) =u(L,t)=0,t[0,T]

    όπου D,α είναι σταθερές και D>0.

    1. (a)

      Γράψτε ένα πεπλεγμένο αριθμητικό σχήμα με σφάλμα διακριτοποίησης ηij το οποίο να ικανοποιεί μια σχέση της μορφής

      |ηij|C(k+h2),

      με C ανεξάρτητη των k και h. Στη συνέχεια, εκφράστε τη μέθοδο σε μορφή γραμμικού συστήματος.

    2. (b)

      Είναι αυτή η μέθοδος ευσταθής; Είναι ευσταθής υπό συνθήκες;

  7. 5.7.

    Θεωρούμε το πρόβλημα

    ut(x,t) =Duxx(x,t)+μ(x)ux(x,t),(x,t)[0,L]×[0,T],
    u(x,0) =g(x),x[0,L],
    u(0,t) =u(L,t)=0,t[0,T],

    όπου D,μ είναι θετικές σταθερές.

    1. (a)

      Γράψτε ένα πεπλεγμένο αριθμητικό σχήμα με σφάλμα διακριτοποίησης ηij το οποίο να ικανοποιεί μια σχέση της μορφής

      |ηij|C(k+h2),

      με C ανεξάρτητη των k και h. Στη συνέχεια, εκφράστε τη μέθοδο σε μορφή γραμμικού συστήματος.

    2. (b)

      Είναι αυτή η μέθοδος ευσταθής; Είναι ευσταθής υπό συνθήκες;

  8. 5.8.

    Θεωρούμε το πρόβλημα

    ut(x,t) =Duxx(x,t)-βu(x,t),(x,t)[0,1]×[0,T],
    u(x,0) =g(x),x[0,1],
    u(0,t) =h(t),u(1,t)=0,t[0,T],

    όπου D,β είναι θετικές σταθερές.

    1. (a)

      Γράψτε ένα πεπλεγμένο αριθμητικό σχήμα με σφάλμα διακριτοποίησης ηij το οποίο να ικανοποιεί μια σχέση της μορφής

      |ηij|C(k2+h2),

      με C ανεξάρτητη των k και h. Στη συνέχεια, εκφράστε τη μέθοδο σε μορφή γραμμικού συστήματος.

    2. (b)

      Είναι αυτή η μέθοδος ευσταθής; Είναι ευσταθής υπό συνθήκες;