1 Εισαγωγή

1.2 Επίλυση προβλημάτων συνοριακών τιμών και αρχικών και συνοριακών τιμών

Σε αυτή την παράγραφο θα δούμε ορισμένες ιδιότητες των διαφορικών εξισώσεων που θα συναντήσουμε, καθώς και των λύσεων τους.

1.2.1 Πρόβλημα δύο σημείων

Θα ασχοληθούμε πρώτα με μια απλή διαφορική εξίσωση σε μία διάσταση, το πρόβλημα δύο σημείων. Έστω u τέτοια, ώστε

$$ -u^{\prime\prime}(x)+q(x)u(x)=f(x),\quad\text{για }x\in[a,b].$$ (1.16)

Είναι γνωστό, βλ. π.χ. (Logan, (2014)), ότι η γενική λύση u της (1.16) προκύπτει ως γραμμικός συνδυασμός δύο γραμμικώς ανεξαρτήτων λύσων της αντίστοιχης ομογενούς, u1 και u2, και μιας λύσης up της μη ομογενούς (1.16),

u=c1u1+c2u2+up,

όπου c1,c2 πραγματικές σταθερές. Αν θεωρήσουμε τώρα ότι η λύση u του (1.16) ικανοποιεί, επιπλέον, τις συνθήκες

u(a)=a0,u(a)=a1, (1.17)

όπου a0,a1 δοσμένες σταθερές, τότε έχουμε ένα πρόβλημα αρχικών τιμών, το οποίο έχει μία ακριβώς λύση, σύμφωνα με το ακόλουθο θέωρημα, [βλ. π.χ. (Logan, (2014), Κεφάλαιο 3.4)].

Θεώρημα 1.1.

Αν οι συναρτήσεις f,q είναι συνεχείς στο διάστημα [a,b], τότε το πρόβλημα αρχικών τιμών (1.16)–(1.17) έχει μοναδική λύση στο [a,b].

Αν τώρα, αντί για τις αρχικές συνθήκες (1.17), θεωρήσουμε ότι η u ικανοποιεί τις ομογενείς συνοριακές συνθήκες Dirichlet,

u(a)=u(b)=0, (1.18)

τότε δεν είναι σίγουρο ότι υπάρχει μοναδική λύση u του προβλήματος (1.16) και (1.18), για κάθε συνεχή συνάρτηση f και q. Όμως, μπορούμε να δείξουμε το ακόλουθο θεώρημα, βλ. π.χ. (Brezis, (1997), Κεφάλαιο VIII).

Θεώρημα 1.2.

Αν οι συναρτήσεις f,q είναι συνεχείς στο διάστημα [a,b] και q0, τότε έχει μοναδική λύση ακόλουθο πρόβλημα: Ζητείται u:[a,b], τέτοια ώστε

$$-u^{\prime\prime}+qu=f,\quad\text{ στο }[a,b],\ \text{ με }u(a)=u(b)=0,$$ (1.19)

Η υπόθεση ότι q0 είναι σημαντική για την ύπαρξη και μοναδικότητα της λύσης του (1.19). Αν θεωρήσουμε, για παράδειγμα, την

$$-u^{\prime\prime}-\pi u=0,\quad\text{ στο }[0,1],\ \text{ με }u(0)=u(1)=0,$$ (1.20)

μπορούμε να δούμε ότι οι συναρτήσεις csin(πx), c αποτελούν λύσεις του (1.20) και επομένως, το πρόβλημα (1.20) έχει άπειρες λύσεις. Στη συνέχεια, αν τροποποιήσουμε το (1.20) και θεωρήσουμε το πρόβλημα

$$-u^{\prime\prime}-\pi u=\pi,\quad\text{ στο }[0,1],\ \text{ με }u(0)=u(1)=0$$ (1.21)

από τη θεωρία των Σ.Δ.Ε., βλ. (Logan, (2014)), έχουμε ότι η γενική λύση αυτής της Δ.Ε. είναι η c1sin(πx)+c2cos(πx)1, με c1,c2. Όμως για καμία επιλογή των c1,c2 δεν προκύπτει λύση του (1.21). Συνεπώς, το (1.21) δεν έχει καμμία λύση.

Αν τροποποιήσουμε τώρα τις συνοριακές συνθήκες Dirichlet στο πρόβλημα (1.19) και θεωρήσουμε το αντίστοιχο πρόβλημα με ομογενείς συνοριακές συνθήκες Neumann, δηλαδή: Ζητείται, u τέτοια ώστε

$$-u^{\prime\prime}+qu=f,\quad\text{ στο }[a,b],\ \text{ με }u^{\prime}(a)=u^{\prime}(b)=0,$$ (1.22)

τότε, όπως και στο Θεώρημα 1.2 μπορούμε να δείξουμε μοναδικότητα της λύσης αν q,f είναι συνεχείς συναρτήσεις στο [a,b] και q>0. Διαφορετικά, αν π.χ. q=f=0, τότε το (1.22) γίνεται

$$-u^{\prime\prime}=0,\quad\text{ στο }[a,b],\ \text{ με }u^{\prime}(a)=u^{% \prime}(b)=0.$$ (1.23)

Είναι προφανές ότι όλες οι συναρτήσεις της μορφής u(x)=c, με c αποτελούν λύση και, άρα, έχουμε άπειρες λύσεις του προβλήματος.

1.2.2 Εξίσωση της θερμότητας

Θα θεωρήσουμε τώρα μια απλή παραβολική διαφορική εξίσωση και συγκεκριμένα την εξίσωση της θερμότητας (1.4), όπου η χωρική μεταβλητή x ανήκει στο διάστημα [0,L], με L>0 και θα θεωρήσουμε, επιπλέον, αρχικές και ομογενείς συνοριακές συνθήκες Dirichlet. Έτσι, έχουμε το ακόλουθο πρόβλημα: Ζητείται u:[0,L]×[0,+), τέτοια ώστε

ut =uxx, στο [0,L]×[0,+) (1.24)
u(0,t) =u(L,t)=0, για t0, (1.25)
u(x,0) =g(x), στο [0,L], (1.26)

με g μια συνεχή συνάρτηση.

Ένας τρόπος για να λύσουμε την (1.24) είναι να υποθέσουμε ότι u(x,t)=X(x)T(t), δηλαδή να θεωρήσουμε τη μέθοδο των χωριζομένων μεταβλητών. Σύμφωνα με αυτή τη μέθοδο, βλ. π.χ. (Logan, (2014)), οι συναρτήσεις un, με un=XnTn, n=1,2,, όπου

un=XnTn, με Xn(x)=csin(λnx), και Tn(t)=de-λn2t με λn=nπL,c,d (1.27)

ικανοποιούν την εξίσωση (1.24) και τις συνοριακές συνθήκες (1.25). Λόγω της αρχής της υπέρθεσης, κάθε γραμμικός συνδυασμός λύσεων της μορφής (1.27) ικανοποιεί τις (1.24)-(1.25). Επομένως, αν θεωρήσουμε μια άπειρη σειρά συναρτήσεων της μορφής

u(x,t)=n=1cne-λn2tsin(λnx),με cn, (1.28)

η οποία συγκλίνει για κάθε σημείο x[0,L] και t0, τότε αποτελεί και λύση του προβλήματος (1.24)-(1.25). Μια λύση u της μορφής (1.28) θα ικανοποιεί και την αρχική συνθήκη (1.26), μόνο αν

u(x,0)=n=1cnsin(λnx)=g(x),με cn. (1.29)

Ένα ανάπτυγμα της μορφής (1.29) καλείται σειρά Fourier ημιτόνων της g, βλ. π.χ. την Παράγραφο 1.3.

Είναι γνωστό ότι το προβλήμα (1.24)-(1.26) έχει μοναδική λύση, βλ. π.χ. (Logan, (2014)).

Θεώρημα 1.3.

Αν η συνάρτηση g είναι συνεχής στο [0,L], τότε υπάρχει μοναδική λύση u του προβλήματος (1.24)–(1.26).

Έστω τώρα ότι γνωρίζουμε το ανάπτυγμα στη μορφή (1.29) της συνάρτησης g(x), της αρχικής συνθήκης (1.26), δηλαδή είναι γνωστές οι σταθερές cn, n=1,2,, στην (1.29). Τότε, η λύση u του προβλήματος (1.24)–(1.26) δίνεται από την (1.28), με σταθερές cn αυτές της (1.29).

Παρόμοια αποτελέσματα ισχύουν αν αντί για τις συνοριακές συνθήκες Dirichlet (1.25), θεωρήσουμε ομογενείς συνοριακές συνθήκες Neumann, δηλαδή αν αναζητήσουμε u, τέτοια ώστε

ut =uxx, στο [0,L]×[0,+) (1.30)
ux(0,t) =ux(L,t)=0, για t0, (1.31)
u(x,0) =g(x), στο [0,L], (1.32)

τότε οι συναρτήσεις

un(x,t)=e-λn2tcos(λnx),λn=nπL,n=0,1,2,,

ικανοποιούν τις (1.30)–(1.31). Έστω τώρα ότι υπάρχουν σταθερές cn, τέτοιες ώστε

u(x,0)=g(x)=n=0cncos(λnx), (1.33)

όπου η σειρά συγκλίνει για κάθε x[0,L]. Αν επιπλέον για κάθε x[0,L] και t0 η σειρά n=1cne-λn2tcos(λnx) συγκλίνει, τότε η

u(x,t)=n=1cne-λn2tcos(λnx),με cn, (1.34)

είναι η μοναδική λύση του προβλήματος αρχικών και συνοριακών τιμών (1.30)–(1.32). Από τη μορφή της λύσης (1.28) ή (1.34), μπορούμε να παρατηρήσουμε ότι η u φθίνει κατά απόλυτη τιμή, καθώς το t τείνει στο άπειρο.

Παράδειγμα 1.1.

Αν u(x,t)=e-4π2tsin(2πx), τότε εύκολα βλέπουμε ότι η u είναι λύση του προβλήματος (1.24)–(1.26), με L=1 και g(x)=sin(2πx). Στο Σχήμα 1.1 απεικονίζεται η u στο [0,1]×[0,T], με T=0.05.

Σχήμα 1.1: Η ακριβής λύση u του Παραδείγματος 1.1.

1.2.3 Εξίσωση μεταφοράς

Θα θεωρήσουμε τώρα μια απλή μερική διαφορική εξίσωση πρώτης τάξης, την εξίσωση μεταφοράς και το αντίστοιχο πρόβλημα αρχικών τιμών: Ζητείται μια συνάρτηση u C1(×[0,T]), T>0, τέτοια ώστε

ut(x,t)+αux(x,t) =0, x,t[0,T], (1.35)
u(x,0) =g(x), x, (1.36)

όπου α και g μια δοσμένη συνάρτηση. Αν θεωρήσουμε την ακόλουθη αλλαγή μεταβλητών

x=αr+s,t=r,

τότε εύκολα μπορούμε να δούμε ότι αν η u είναι λύση του (1.35), τότε

ur=xrux+trut=αux+ut=0.

Επομένως, αν θεωρήσουμε ότι

z(s,r):=u(x-αt,t),

έχουμε zr=0, δηλαδή η z δεν εξαρτάται από το r και είναι σταθερή ως προς r. Συνεπώς, η u πάνω σε όλες τις ευθείες x-αt=c, με c σταθερά, δεν μεταβάλλει την τιμή της και, έτσι, u(x,t)=z(s,r)=z(s,0)=u(x-αt,0)=g(x-αt). Λόγω αυτού του γεγονότος, οι ευθείες x-αt=c καλούνται χαρακτηριστικές για τη διαφορική εξίσωση (1.35), βλ. Σχήμα 1.2.

Σχήμα 1.2: Οι χαρακτηριστικές ευθείες x-αt=c πάνω στις οποίες είναι σταθερή η λύση u της διαφορικής εξίσωσης (1.35).

Παράδειγμα 1.2.

Ας θεωρήσουμε το πρόβλημα (1.35) και ότι η λύση u για t=0 δίνεται από τη συνάρτηση

g(x)={1 αν 0x10 διαφορετικά.

Τότε η λύση είναι

u(x,t)={1 αν 0x-αt10 διαφορετικά,

ή ισοδύναμμα

u(x,t)={1 αν αtx1+αt0 διαφορετικά,
Σχήμα 1.3: Η αρχική συνθήκη g και η ακριβής λύση u του Παραδείγματος 1.2 στην περίπτωση που α>0 ή α<0.

Παρατηρούμε ότι για να προσδιορίσουμε με μοναδικό τρόπο τη λύση της διαφορικής εξίσωσης (1.35), αρκεί να γνωρίζουμε την αρχική συνθήκη u(x,0)=g(x). Επίσης, αν α>0, μπορούμε να διατυπώσουμε το πρόβλημα (1.35) ως πρόβλημα αρχικών και συνοριακών τιμών ως εξής: Ζητείται μια συνάρτηση u C1([a,b]×[0,T]), a<b, T>0, τέτοια ώστε

ut(x,t)+αux(x,t)=0,x[a,b],t[0,T]u(x,0)=g(x),x[a,b]u(a,t)=ϕ0(t),t[a,b] (1.37)

όπου ϕ0(t)=g(a-αt). Συνεπώς, σε αντίθεση με τα προβλήματα που είδαμε στις Παραγράφους 1.2.1 και 1.2.2, αρκεί μια συνοριακή συνθήκη για το σημείο x=a για να προσδιορίσουμε μοναδικά τη λύση του προβλήματος (1.37). Όμως, στην περίπτωση τώρα που α<0, οι χαρακτηριστικές θα έχουν διαφορετική κατεύθυνση και διατυπώνουμε το πρόβλημα (1.35) ως πρόβλημα συνοριακών τιμών ως εξής: Ζητείται μια συνάρτηση u C1([a,b]×[0,T]), a<b, T>0, τέτοια ώστε

ut(x,t)+αux(x,t)=0,x[a,b],t[0,T]u(x,0)=g(x),x[a,b]u(b,t)=ϕ1(t),t[a,b] (1.38)

όπου ϕ1(t)=g(a-αt). Και σε αυτή την περίπτωση αρκεί μια συνοριακή συνθήκη για x=b για να προσδιορίσουμε μοναδικά τη λύση του προβλήματος (1.38).

Είναι προφανές ότι η λύση u του προβλήματος (1.35) διατηρεί τη μορφή της αρχικής συνάρτησης, καθώς ο χρόνος t αυξάνει και ιδιαίτερα η μέγιστη κατά απόλυτο τιμή παραμένει σταθερή. Στο Σχήμα 1.3 βλέπουμε τη λύση του Παραδείγματος 1.2, για t=3, αν η σταθερά α είναι θετική ή αρνητική.


1.2.4 Εξίσωση κύματος

Θα θεωρήσουμε τώρα μια απλή υπερβολική διαφορική εξίσωση και συγκεκριμένα την εξίσωση του κύματος (1.6), όπου η χωρική μεταβλητή x ανήκει στο διάστημα [0,L], με L>0 και θα θεωρήσουμε, επιπλέον, αρχικές και ομογενείς συνοριακές συνθήκες Dirichlet. Έτσι, έχουμε το ακόλουθο πρόβλημα: Ζητείται u:[0,L]×[0,+), τέτοια ώστε

utt =α2uxx, στο [0,L]×[0,+) (1.39)
u(0,t) =u(L,t)=0, για t0, (1.40)
u(x,0) =g(x), και ut(x,0)=g1(x), στο [0,L], (1.41)

με α>0 και g,g1 συνεχείς συναρτήσεις.

Χρησιμοποιώντας τη μέθοδο των χωριζομένων μεταβλητών, βλ. π.χ. (Logan, (2014)), οι συναρτήσεις un, με un=XnTn, n=1,2,, όπου

un =XnTn, (1.42)
Xn(x) =csin(λnx), και (1.43)
Tn(t) =d1sin(λnαt)+d2cos(λnαt) με λn=nπL,c,d1,d2, (1.44)

ικανοποιούν την εξίσωση (1.39) και τις συνοριακές συνθήκες (1.40). Λόγω της αρχής της υπέρθεσης, κάθε γραμμικός συνδυασμός λύσεων της μορφής (1.42) ικανοποιεί τις (1.39)-(1.40). Επομένως, αν θεωρήσουμε μια άπειρη σειρά συναρτήσεων της μορφής

u(x,t)=n=1sin(λnx)(cnsin(λnαt)+dncos(λnαt)),με cn,dn, (1.45)

η οποία συγκλίνει για κάθε x[0,L] και t0, τότε αυτή αποτελεί και λύση του προβλήματος (1.39)-(1.41). Μια λύση u της μορφής (1.45) θα ικανοποιεί και την αρχική συνθήκη (1.41), αν

u(x,0) =n=1dnsin(λnx)=g(x),με dn, (1.46)
ut(x,0) =n=1cnλnαsin(λnx)=g1(x),με cn. (1.47)

Ένας διαφορετικός τρόπος για την περιγραφή της λύσης του προβλήματος (1.39)– (1.41) είναι η αναπαράσταση d’Alembert της λύσης, η οποία οδηγεί στην

u(x,t)=Φ(x-αt)+Ψ(x+αt), (1.48)

όπου Φ,Ψ είναι συνεχείς συναρτήσεις που εξαρτώνται από τις αρχικές συνθήκες (1.41). Μάλιστα, εύκολα μπορούμε να δούμε ότι στην περίπτωση που στην αρχική συνθήκη (1.41) η συνάρτηση g1 μηδενίζεται στο [0,L], τότε η ακριβής λύση u του προβλήματος (1.39)–(1.41) δίνεται από

u(x,t)=12(g(x-αt)+g(x-αt)). (1.49)
Παράδειγμα 1.3.

Αν u(x,t)=sin(2πx)(sin(2πt)+cos(2πt)), τότε εύκολα βλέπουμε ότι η u είναι λύση του προβλήματος (1.24)–(1.26), με L=α=1, g(x)=sin(2πx) και g1(x)=2πsin(2πx). Στο Σχήμα 1.4, απεικονίζεται η u στο [0,1]×[0,T], με T=1.

Σχήμα 1.4: Η ακριβής λύση u του Παραδείγματος 1.3.

1.2.5 Ελλειπτική εξίσωση

Θεωρούμε το ακόλουθο πρόβλημα: Ζητείται μια συνάρτηση u:Ω=(0,1)×(0,1), τέτοια ώστε

-(uxx(x,y)+uyy(x,y))+q(x,y)u(x,y)=f(x,y),(x,y)Ωu(x,y)=g(x,y),(x,y)Ω (1.50)

όπου Ω είναι το σύνορο του Ω, q,fC(Ω), gC(Ω) και q(x,y)0, για κάθε (x,y)Ω. Στη συνέχεια, θα συμβολίζουμε με qmin=min(x,y)Ωq(x,y) και Ω¯=[0,1]×[0,1].

Η υπόθεση q0 είναι σημαντική για την ύπαρξη και μοναδικότητα της λύσης του προβλήματος (1.50). Παραδείγματος χάριν, αν θεωρήσουμε την

-(uxx(x,y)+uyy(x,y))-2π2u(x,y)=0,(x,y)Ωu(x,y)=g(x,y),(x,y)Ω (1.51)

τότε οι συναρτήσεις u(x,y)=csin(πx)sin(πy), c, αποτελούν λύσεις και, επομένως, το πρόβλημα (1.51) έχει άπειρες λύσεις.

Για τη διαφορική εξίσωση (1.50) ισχύει η αρχή του μεγίστου, η οποία δίνεται από το ακόλουθο θεώρημα, βλ. π.χ. (Larsson and V. Thomée, (2009)).

Θεώρημα 1.4.

Έστω ότι η λύση u του (1.50) είναι τέτοια, ώστε uC2(Ω¯). Τότε

  1. 1.

    Αν

    f(x,y)0,για κάθε(x,y)Ω,

    έχουμε ότι

    max(x,y)Ω¯u(x,y){max(x,y)Ωg(x,y),αν q0,max{max(x,y)Ωg(x,y),0},αν q0. (1.52)
  2. 2.

    Αν

    f(x,y)0,για κάθε(x,y)Ω,

    έχουμε ότι

    min(x,y)Ω¯u(x,y){min(x,y)Ωg(x,y),αν q0,min{min(x,y)Ωg(x,y),0},αν q0. (1.53)
Παράδειγμα 1.4.

Αν u(x,t)=sin(πx)sin(πy), τότε εύκολα βλέπουμε ότι η u είναι λύση του προβλήματος (1.50), με g=0, q=1 και f(x,y)=(2π2+1)sin(πx)sin(πy). Στο Σχήμα 1.5, απεικονίζεται η u στο [0,1]×[0,1].

Σχήμα 1.5: Η ακριβής λύση u του Παραδείγματος 1.4.