phương trình đạo hàm

Cấp cao nhất của đạo hàm riêng của , có mặt trong phương trình được gọiuChẳng hạn, phương trình cấp một của hàm hai biến có dạng:F x,y,u,∂u∂y = 0.. 1.3Ví dụ: ∆u = 0 phương trình Laplace,

BỘ GIÁO DỤC VÀ ĐÀO TẠO

TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI

Mục lục

1 Đại cương về phương trình vật lý toán 4

1.1 Đại cương về phương trình vật lý toán 4

1.1.1 Các định nghĩa 4

1.1.2 Một số ví dụ giải phương trình đạo hàm riêng 5

1.1.3 Thiết lập một số phương trình đạo hàm riêng 6

1.2 Phân loại phương trình đạo hàm riêng cấp hai 8

1.2.1 Trường hợp hai biến số 8

1.2.2 Dạng chính tắc 8

1.2.3 Đưa phương trình về dạng chính tắc 9

1.2.4 Phân loại phương trình 10

1.3 Trường hợp nhiều biến số 13

1.3.1 Nhắc lại kết quả về đại số 13

1.3.2 Phân loại phương trình tuyến tính cấp 2 14

1.4 Khái niệm về mặt đặc trưng 15

1.4.1 Khái niệm về mặt đặc trưng 15

1.4.2 Bài toán Cauchy và bài toán Cauchy với dữ kiện cho trên các mặt đặc trưng 16

1.5 Ba loại phương trình vật lý toán cơ bản Bài toán đặt chỉnh 17

1.5.1 Ba loại phương trình cơ bản 17

1.5.2 Các dạng bài toán có điều kiện biên, điều kiện ban đầu 19

1.5.3 Bài toán đặt chỉnh 19

1.6 Phương pháp đặc trưng 20

2 Phương trình truyền sóng 21 2.1 Phương trình dao động của dây, của màng 21

2.2 Bài toán biên ban đầu (Bài toán hỗn hợp) 23

2.2.1 Định luật bảo toàn năng lượng 23

2.2.2 Tính duy nhất nghiệm 24

2.2.3 Đánh giá tiên nghiệm Sự phụ thuộc liên tục của nghiệm vào điều kiện ban đầu 25

1

2.3 Bài toán Cauchy Công thức Kirchoff 26

2.4 Công thức Poisson, công thức Dalembert 28

2.4.1 Công thức Poisson 28

2.4.2 Công thức Dalembert 28

2.5 Phương pháp tách biến (Fourier) 32

2.5.1 Dao động tự do của dây 32

2.5.2 Phương trình không thuần nhất 34

3 Phương trình truyền nhiệt 37 3.1 Thiết lập phương trình truyền nhiệt 37

3.2 Bài toán biên ban đầu (Bài toán hỗn hợp) 38

3.2.1 Nguyên lý cực đại 39

3.2.2 Tính duy nhất nghiệm 41

3.2.3 Sự phụ thuộc liên tục của nghiệm vào điều kiện biên và ban đầu 42

3.2.4 Phương pháp tách biến 42

3.2.5 Bài toán biên loại hai 47

3.3 Sự truyền nhiệt trong thanh vô hạn Bài toán Cauchy 48

3.3.1 Nguyên lý cực đại Tính duy nhất nghiệm 49

3.3.2 Biến đổi Fourier giải bài toán Cauchy 49

3.3.3 Giải bài toán Cauchy 51

3.3.4 Sự truyền nhiệt trong nửa thanh vô hạn 51

3.3.5 Bài toán Cauchy trong không gian hai chiều và ba chiều 53

4 Phương trình Laplace và phương trình Poisson 55 4.1 Phương trình Laplace Nghiệm cơ bản 55

4.2 Các bài toán biên: Bài toán biên Dirichlet, bài toán biên Neumann 58

4.2.1 Bài toán biên Dirichlet 58

4.2.2 Bài toán biên Neumann 58

4.2.3 Bài toán biên Robin 59

4.2.4 Bài toán Cauchy 59

4.3 Phương pháp hàm Green để giải bài toán Dirichlet Công thức Poisson 59

4.3.1 Hàm Green đối với hình cầu và công thức Poisson 60

4.4 Giải bài toán Dirichlet trong mặt tròn 62

4.5 Các tính chất của hàm điều hòa 65

4.6 Lý thuyết thế vị 67

4.7 Nghiệm yếu Phương pháp biến phân 67

4.7.1 Bài toán biên Dirichlet 67

4.7.2 Không gian Sobolev H1(Ω) 68

4.7.3 Toán tử Elliptic 68

4.7.4 Bài toán yếu 69

Chương 1

Đại cương về phương trình vật lý toán

1.1 Đại cương về phương trình vật lý toán

Ví dụ: Phương trình chuyển dịch tuyến tính: ut+ ni=1

biux i= 0, phương trình Burger (sửa đổi): ut+u.ux= 0.Phương trình cấp hai của hàm hai biến số có dạng

Hệ phương trình đạo hàm riêng là hệ gồm các phương trình đạo hàm riêng của một hay nhiều ẩn hàm vàcác đạo hàm của chúng

Phương trình đạo hàm riêng được gọi là tuyến tính, nếu như nó tuyến tính đối với ẩn hàm và tất cả cácđạo hàm của nó với các hệ số phụ thuộc vào các biến độc lậpx1,x , 2 Trong trường hợp tổng quát ta có thểviết nó dưới dạng

Ví dụ:

∆u = f phương trình Poisson,

ut− ∆u = f phương trình truyền nhiệt,

utt− ∆u = f phương trình truyền sóng,ihψt= −h22m∆ψ + Vψ phương trình Schrodinger,

∞

X

n=0

1n!

d un

dxn= g(x)

Phương trình đạo hàm riêng cấpmđược gọi là nửa tuyến tính (semi-linear) nếu nó tuyến tính với các đạohàm cấp m với các hệ số phụ thuộc vàox1,x , 2 , được gọi là tựa tuyến tính nếu nó tuyến tính với các đạohàm cấpmvới các hệ số phụ thuộc vàox1,x , 2 và các đạo hàm cấp bé hơn Phương trình đạo hàm riêngmkhông phải là phương trình tuyến tính được gọi là phương trình phi tuyến Ví dụ:

|Du| = 1 phương trình Eikonal,

ut+ cuux+uxxx= 0 phương trình Korteweg-de Vries,

det(D2u) = f phương trình Monge-Ampere, được đưa ra từ năm 1775

Nghiệm của phương trình đạo hàm riêng là một hệ bất kỳ các hàm sao cho khi thay vào phương trình (1.1)trở thành đồng nhất thức trong miềnΩnào đó của các biến số độc lập Để đơn giản, ta giả sử rằngx1,x , 2

là thực và các đạo hàm của trong phương trình (1.1) liên tục theo các biếnu x1,x , 2 trong miền ΩChú ý 1 Số chiều không gian của các biến độc lậpx1,x , 2 có thể là vô hạn Ngoài ra, cấp của đạo hàm

trongphươngtrình(1.1) cóthểlàvôhạn.Khiđó,tacóphươngtrìnhđạohàmriêngcấpvôhạn

1.1.2 Một số ví dụ giải phương trình đạo hàm riêng

Ví dụ 1: Giải phương trìnhuxy= 0, x,ylà các biến số độc lập

Zx+y

0 c(t) dt+ g(x − y) = h(x + y) + g(x − y)

Ví dụ 4: Phương trình tuyến tính cấp 1 hệ số hằng Xét phương trình

aux+buy= 0, (1.5)trong đóa,blà các hằng số không đồng thời bằng không

Phương pháp hình học: Biểu thứcau bux+ ylà đạo hàm theo hướng của hàm theo véctơu →−V = (a,b) = a→−i +b→j −

Do đó, từ phương trình (1.5), hàmu(x,y)là hằng số theo hướng→−V Đường thẳng song song với véctơ→−V cóphương trình làbx−ay= cvà dọc theo đường thẳng này, hàm là hằng số Như vậy,u u(x,y) = f(c) = f( − )bx ayvới c bất kỳ Vậy nghiệm là

u(x,y) = f(bx−ay) với mọi x,y

6Phương pháp đổi biến: Đổi biến số

x′= ax+ by, y′= bx− ay

Bài 1: Giải thích các đẳng thức vi phân sau có là phương trình đạo hàm riêng không?

1 cos(ux+ uy) − cosuxcosuy+ sin uxsin uy= 0

2 u2

xx+ u2

yy− (uxx− uyy)2= 0

3.sin (2uxx+ uyy) + cos2(uxx+ uyy) − u = 1

Bài 2: Xác định cấp của các phương trình sau

1.log|u uxx yy|−log|uxx|−log|uyy| + ux+ uy= 0

1.1.3 Thiết lập một số phương trình đạo hàm riêng

Phương trình vận tải đơn giản

Xét một luồng chất lỏng, chẳng hạn là nước, chảy với tốc độ -hằng số, dọc theo ống nằm ngang theo hướngcdương x Một lượng chất bẩn nổi trên nước Gọiu(x,t)là mật độ của nó (g/cm) ở thời điểm Ta biết rằng,tlượng chất bẩn trong khoảng[0,b]ở thời điểmtlà M =Rbu(x,t) dx(g) Ở thời điểm t + h, lượng chất đó đã

di chuyển sang phải một đoạnc.h(cm) Do đó,

M =

Zbu(x,t) dx=

Zb+ch

ch u(x,t + )h dx

Lấy đạo hàm theo , ta đượcb

u(b,t) = u(b +ch,t+ h )Lấy đạo hàm theo , và choh h = 0, ta được

0 = cux(b,t) + ut(b,t)

Do b là bất kỳ, nên phương trình mô tả bài toán vận chuyển đơn giản có dạng

ut+cux= 0 (1.6)

7Phương trình dao động

Phương trình khuếch tán (truyền nhiệt)

Bây giờ ta thiết lập phương trình truyền nhiệt hay phương trình khuếch tán Xét một vật rắn truyền nhiệtđẳng hướng Ký hiệuu(x,y,z,t)là nhiệt độ của một vật rắn tại điểm(x,y,z)ở thời điểm Chúng ta biếttrằng nhiệt sẽ truyền từ nơi có nhiệt độ cao sang nơi có nhiệt độ thấp hơn Sự truyền nhiệt này tuân theo địnhluật sau: nhiệt lượng∆Qđi qua một mảnh mặt khá bé∆Schứa điểm (x,y,z) trong một thời gian∆ttỷ lệvới ∆S, ∆t và đạo hàm theo pháp tuyến∂u , tức là

∆Q = −k(x,y,z t S.)∆∆ ∂u ∂n, (1.7)trong đó k(x,y,z) > 0 là hệ số truyền nhiệt trong (k không phụ thuộc vào hướng của pháp tuyến với∆Svì

sự truyền nhiệt là đẳng hướng), −→nlà vectơ pháp tuyến của∆Shướng theo chiều giảm của nhiệt độ.Gọi q là dòng nhiệt, tức là nhiệt lượng đi qua một đơn vị diện tích trong một đơn vị thời gian Từ (1.7) ta suyra

q = −k∂u

∂n.Xét một thể tích tùy ýVcủa vật rắn giới hạn bởi một mặt kín, trơnSvà xét sự biến thiên của nhiệt lượngtrong thể tích đó trong khoảng thời gian từt1đến t2 Từ (1.7) ta suy ra rằng nhiệt lượng qua mặtSvàotrong, từ thời điểm t1đến thời điểmt2bằng

Q1= −

t 2

Z

t 1 dtZZ

S

k(x,y,z)∂u

∂ndS,trong đó −→nlà vectơ pháp tuyến hướng vào trong của mặt Áp dụng công thức Ostrogradsky, chuyển tíchSphân mặt sang tích phân ba lớp, ta được

Q3=

Z Z Z

V

[ (u x,y,z,t2) − u x,y,z,t( 1)] (c x,y,z ρ x,y,z) dxdy dz,) (

trong đóc(x,y,z)là nhiệt dung, ρ(x,y,z)là mật độ của vật Ta có

Q2=

t 2

Zdt

Z Z Z cρ∂u

∂tdxdy dz.

∂t−div(k.gradu) − F x,y,z,t( ) dxdy dz = 0.

Vì khoảng thời gian(t1,t2)và thể tíchVđược chọn tùy ý, nên biểu thức dưới dấu tích phân bằng không,

cρ∂u∂t=div(k.gradu) + F x,y,z,t( )hay

cρlà hệ số truyền nhiệt độ, f(x,y,z,t) =F (x,y,z,t)cρ

1.2 Phân loại phương trình đạo hàm riêng cấp hai

1.2.1 Trường hợp hai biến số

2uxx+ 7uyy= 0là elip trên toàn mặt phẳng

uxx− 5uyy+ ux− 2uy+ u = 0là hypecbôn trên toàn mặt phẳng

uxx− 2uy+ 3u = 0là parabôn trên toàn mặt phẳng

1.2.2 Dạng chính tắc

a) Phương trình uxx−uyy+F (x,y,u,u ,ux y) = 0gọi là dạng chính tắc của phương trình loại hypecbôn Người

ta cũng gọi phương trìnhuxy+ F (x,y,u,u ,ux y) = 0là dạng chính tắc của phương trình loại hypecbôn.Chú ý rằng từ dạng đầu có thể suy ra dạng thứ hai bằng cách đổi biến sốξ = x+ y, η= − yx b) Phương trìnhuxx+ uyy+ F (x,y,u,u ,ux y) = 0là dạng chính tắc của phương trình loại elip

c) Phương trình uxx+ F (x,y,u,u ,ux y) = 0hay uyy+ G(x,y,u,u ,ux y) = 0gọi là dạng chính tắc củaphương trình parabôn

a dy( )2− bdxdy+ ( )2 c dx2= 0 (1.12)Chứngminh.a) Giả sử φ(x,y) là một nghiệm riêng của phương trình (1.11) và giả sửφy = 0 Khi đó từ (1.11)

φy.Vậy từ (1.13) suy ra a( )y′ 2− 2by′+ c = 0haya dy( )2− 2bdydx+ ( )c dx2= 0

Trường hợpφy= 0, do giả thiết suy raφx = 0 Khi đó ta cũng lý luận tương tự như trường hợpφy = 0.b) Giả sửφ(x,y) = Clà tích phân tổng quát của phương trình (1.12) Ta sẽ chứng minhz = φ(x,y)là nghiệm

10riêng của phương trình (1.11) trong miền Ω

Thật vậy, giả sửφ(x0,y0) = C0với (x0,y0) ∈ Ω Xét quan hệφ(x,y) =C0; giả sửφy = 0 Khi đó quan hệ nàyxác định một hàm ẩny = y(x)khả vi liên tục vày′(x) = −φx

và do vậy tại(x0,y0)ta có đồng nhất thức (1.11) Trường hợp tại(x0,y0)mà φy= 0thì từ giả thiết suy ra

φx = 0 Ta lại lập luận tương tự như đối vớiφy = 0

1.2.4 Phân loại phương trình

Trên cơ sở bổ đề 1 ta sẽ tìm cách đổi biến số sao cho một trong các hệ sốa1,b ,c1 1của phương trình (1.10)

bị triệt tiêu Để làm điều đó ta sẽ nghiên cứu phương trình (1.12) mà cơ sở của nó là biệt thức∆ = b2−ac.Trường hợp∆ > 0:

Khi đó phương trình (1.8) là loại hypecbôn (hyperbolic) Lúc này phương trình

a(y′)2−2by′+ c = 0 (1.14)

có hai nghiệm phân biệt

a) Nếua = 0: Từ phương trình (1.14) suy ra

là các tích phân tổng quát của phương trình (1.12)

Khi đó ta dùng phép đổi biến ξ = φ(x,y ,η) = ψ(x,y)với φ,ψ ∈ C2(R2)và

D(ξ,η)D(x,y)=

φx φy

ψx ψy = 0(theo lý thuyết phương trình vi phân), vì y1′(x) = −φx

φy = y2(x) = −′ ψxψ

y.Với phép đổi biến này thì các hệ số a1và c1của phương trình (1.10) bằng 0 Khi đó phương trình (1.10) trởthành

2b1uξη+ F1(ξ,η,u,u ,uξ η) = 0 (1.15)với b1= aξxηx+ (b ξxηy+ξyηx) +cξyηy Lúc này b2

1− a1c1= (b2− ac)(ξxηy− ξyηx)2> 0 Vậyb1 = 0nênphương trình (1.15) có thể đưa về dạng

uξη= F2(ξ,η,u,u ,uξ η) (1.16)b) Nếua = 0:

+) Khi c = 0, do phương trình (1.8) là phương trình cấp 2 nên suy rab = 0 Khi đó phương trình (1.8) là có

11dạng chính tắc.

+) Khi c = 0 thì phương trình (1.12) tương đương với phương trình c(x′(y))2− bx2 ′(y) = 0 Khi đó ta làmgiống như trường hợpa = 0

Chú ý 2 Nếudùngphépđổibiếnξ = α+ β,η = − β thìphươngtrìnhα (4.42)sẽtrởthành

uαα−uββ= F3(α,β,u,uα,u β)

Trường hợp ∆ = 0:

Khi đó phương trình (1.8) thuộc dạng parabôn (parabolic)

a) Nếub = 0, do ∆ = b2− ac= 0suy raac= 0 Khi đó hoặc a = 0 hoặc c = 0 Do đó phương trình (1.8) códạng chính tắc cuyy+ F = 0hoặcauxx+ F = 0

b) Nếub = 0, do∆ = b2− ac= 0suy raavà c khác 0 Khi đó do phương trình (1.14) có nghiệm kép là

y′=b

a⇒ y =

R b

adx+ C1suy ra (x,y) = y −φ R b a

dx= C1là tích phân tổng quát của phương trình (1.12)

Ta chọn phép đổi biến sốξ = φ(x,y)còn η = ψ(x,y) ∈ C2(R2)tùy ý sao choD(ξ,η)

D(x,y)nên b1= 0 Ta sẽ chỉ rac1 = 0 Thậtvậy, nếu c1= 0thì ta cóc1= (ηa x)2+ 2bηxηy+ (ηc y)2= 0 Do giả thiết∆ = b2− ac= 0nênb2= ac>0

Giả sửb>0, suy ra b =√ac Vậy

c1= (a ψx)2+ 2√acψxψy+ (c ψy)2= 0

aψx+ cψy)2= 0 nếu a>0,c > 0(√−aψx+ −cψy) = 0 nếu a< ,c < 0 0Tương tự ta có (√

aφx+ cφy)2= 0 nếu a>0,c > 0(√−aφx+ −cφy) = 0 nếu a< ,c < 0 0

Từ đó ta có hệ (√

aψx+ cψy)2= 0(√aφx+ cφy) = 0 nếu a>0,c > 0

−aψx+ −cψy)2= 0(√−aφx+ −cφy)2= 0 nếu a<0,c < 0

Từ đó suy raD(φ,ψ)D(x,y)= J = 0(mâu thuẫn) Khib<0ta lập luận tương tự

Vậy c1 = 0 Do đó dạng chính tắc là c u1 ηη= −F (ξ,η,u,u ,uξ η)hay

Do đó φ(x,y) = y−Rb +i√ac− b 2

a dx= C1hay φ(x,y) = α(x,y)+ (x,y) = Ciβ 1và ψ(x,y) = y−Rb − i√ac − b 2 a dx=

C2hay ψ(x,y) = α x,y) − iβ(x,y) = C( 2là các tích phân tổng quát của phương trình (1.12), trong đóα(x,y ,β x,y) ( )là các hàm thực và (φ′

aφx= −bαy+ −∆βy, (1.17)

aβx= −bβy− −∆αy (1.18)Sau khi nhân βyvào hai vế của (1.17) vàαyvào hai vế của (1.18) và trừ vế với vế ta đượcaJ=√−∆(α2 y+β2)

13Vậy ta nhận được dạng chính tắc là

4uξη+ 8uξ= 0 (1.19)Chú ý 4 Ta có thể giải phương trình(1.19) bằng cách đặtv = uξ, ta có phương trìnhvη+ 2v = 0, suy ra

ln |v| = −2η + C1.Từđósuyrav = e−2ηC(ξ) vàdođó

uξ= e−2η φ ,

ZC(ξ) dξ+ (η)ởđóC(ξ) làhàmtùyýchỉphụthuộcvàoξ cònφ(η) làhàmtúyýchỉphụthuộcvào η

1.3 Trường hợp nhiều biến số

1.3.1 Nhắc lại kết quả về đại số

Cụ thể là số các phần tử dương, số các phần tử âm, số các phần tử bằng 0 trên đường chéo củaBlần lượt bằng

số các nghiệm dương, số các nghiệm âm, số các nghiệm bằng 0 của phương trình đặc trưng det(A − λE) = 0,

a∗

ilà nghiệm của phương trình det(A −λE) = 0

Nếu dùng phép đổi biến thích hợp ta có thể đưa dạng (1.23) về dạng

aij=aji, x= (x1,x2, ,xn) ∈ Ω ⊂ Rn.Lấy điểm x0= (x0 ,x0, ,x0 )cố định trong Ω Ma trận A = (aij(x0))là ma trận hằng đối xứng cấp nXét phương trình đặc trưng

det(A − λE) = 0, với E là ma trận đơn vị (1.26)Định nghĩa 3 a) Phương trình(1.25)gọi là thuộc loại elliptic tại điểmx0nếu như tất cả nghiệm củan

phương trình đặc trưng(1.26)đều khác không và cùng dấu (Khi đó dạng toàn phươngPaij( )x0titjlà xác

địnhdươnghayâm)

b)Phươngtrình(1.25)gọilàthuộcloạihyperbolictạiđiểmx0nếunhưtấtcả nghiệmcủaphươngtrìnhđặcn

trưng(1.26) đềukháckhôngvàcó(n − 1)nghiệmcùngdấu,cònnghiệmcuốicùngcònlạikhácdấu

c) Phương trình(1.25)gọi là thuộc loại parabolic tại điểmx0nếu như trong nghiệm của phương trình đặcn

trưng(1.26) cómộtnghiệmbằng0,còn(n − 1)nghiệmcònlạikháckhôngvàcócùngmộtdấu

Tương tự như đối với phương trình tuyến tính cấp 2 hai biến ta cũng có khái niệm phương trình thuộc loạielliptic, hyperbolic, parabolic trong miền Liên hệ với phương trình tuyến tính cấp 2 hai biếnΩ

a(x,y)uxx+ 2b(x,y)uxy+ c(x,y)uyy+ F (x,y,u,u ,ux y) = 0

Xét phương trình đặc trưngdet(A −λE) = 0

A = a(x0,y0) b x(0,y0)

b(x0,y0) c(x0,y0) , (x0,y0) ∈ Ω

Phương trình trên tương đương với λ2− (a + c)λ + ac− b2= 0, có ∆ = (a+ c)2− 4ac+4b2= ( − c)a 2+4b2≥ 0.Chú ý: Khia = c,b = 0(khi đó b2− ac<0) thì phương trình thuộc loại elliptic Do vậy ta chỉ xét với∆ > 0.+) Khib2− ac<0thì phương trình đặc trưng có hai nghiệm cùng dấu nên phương trình thuộc loại elliptic.+) Khi b2− ac>0thì phương trình đặc trưng có hai nghiệm trái dấu nên phương trình thuộc loại hyperbolic.+) Khi b2− ac= 0thì phương trình đặc trưng có một nghiệm bằng 0 và một nghiệm khác không, vì khôngxảy ra (a + c) = 0 Do đó phương trình thuộc loại parabolic

1.4 Khái niệm về mặt đặc trưng

1.4.1 Khái niệm về mặt đặc trưng

n

X

i,j=1

aij(x)ω ωx i x j= 0 (1.28)Phương trình (1.28) được gọi là phương trình các mặt đặc trưng (khin = 2được gọi là phương trình cácđường đặc trưng) của phương trình (1.27)

Mặt S nào đó có phương trìnhω(x1,x2, ,xn) = 0được gọi là mặt đặc trưng của phương trình (1.27) nếunhư trên mặt S hàm ω(x1,x2, ,xn)là nghiệm của phương trình (1.28) vàPn i=1ω2

i = 0

Rõ ràng rằngω(x1,x2, ,xn)là một nghiệm của phương trình (1.28) thì họ các mặt cong xác định bởi phươngtrình ω = C, C là hằng số tùy ý, là họ các đường (n = 2), họ các mặt (n ≥ 3) đặc trưng của phương trình(1.27)

Vectơ ξ = (ξ1,ξ2, ,ξn) = 0được gọi là vectơ có hướng đặc trưng của phương trình (1.28) nếu ta có

Pni,j=1aij(x)ξiξj= 0

Do vậy mặt congS(không có điểm kỳ dị, nghĩa là tại mọi điểmx ∈ Snó đều có vectơ pháp tuyếnn(x) = 0)

là mặt đặc trưng của phương trình (1.27) nếu như tại mọi điểmx ∈ S, vectơ pháp tuyếnn(x)có hướng đặctrưng

Bây giờ ta mô tả một số mặt đặc trưng của một vài phương trình quen thuộc Xét phương trình (1.27)trong một miềnΩnào đó

a) Giả sử phương trình (1.27) thuộc loại elliptic trong miền Khi đó, dạng toàn phươngΩ Paij(x t)itjlà xácđịnh dương (hoặc âm), nênPaij(x)ti jt = 0chỉ xảy ra khi t1=t2= =tn= 0 Vậy phương trình (1.28)không xác định một mặt đặc trưng nào Ta nói rằng phương trình thuộc loại elliptic không có mặt đặc trưngthực

b) Xét phương trình (1.27) (n = 2), thuộc loại hyperbolic Khi đó, phương trình đặc trưng có dạnga(ωx) +2

2bωxωy+ c(ωy) = 0

Giả sử z = φ(x,y) là hàm nghiệm, thì các đường cong φ(x,y) = C là các đường đặc trưng (Do phương trình

at2+ 2bt+ c = 0có hai nghiệm phân biệt)

Chẳng hạn, phương trình có dạngutt=a2uxx+ F (x,t,u,u ,ux y) = 0 a>0, ( ) Phương trình đặc trưng của nó

là (ωt)2−a ω2( x)2= 0 ⇔ωt= ±aωx, suy ra phương trình này xác định các mặt đặc trưng dạngt = ±Ax+ C,

t ±1p(x − x0)2+ (y − y0)2= C, Clà hằng số,

là các mặt đặc trưng LấyC = t0, ta viết phương trình này dưới dạng

(x − x0) + (y − y0) − a(t − t0) = 0

Đây là phương trình mặt nón tròn xoay có đỉnh(x0,y ,t0 0)và trục song song vớiOt

d) Xét phương trình (1.27), (khin = 2), thuộc loại parabolic Xét phương trình truyền nhiệtut−a u2

Trong lân cận của mặt , tìm một nghiệm của phương trình (1.29) thỏa mãn điều kiệnS

u|S= (x),φ ∂u

∂r|S= ψ(x), (1.30)trong đó φ( )x, ψ(x) là những hàm tùy ý cho trên mặt Ta giả thiếtS ψ(x)là hàm liên tục cònφ(x)là mộthàm khả vi liên tục trên S

Các hàm φ( )x, ψ(x) gọi là các dữ kiện Cauchy, cònSgọi là mặt mang dữ kiện Cauchy Ta có hai kết quảquan trọng sau:

φ1(x)có thể cho tùy ý

Ví dụ 2: Bài toán Cauchy đối với phương trình vi phân thường

y′′= f(x,y,y′)y(x0) = y0

y′(x0) = y1.Các dữ kiện Cauchyy0,y1cho tùy ý

17Bài toán Cauchy đối với phương trình cấp 2

1.5 Ba loại phương trình vật lý toán cơ bản Bài toán đặt chỉnh

1.5.1 Ba loại phương trình cơ bản

Fđi qua ∂Ω bị triệt tiêu, nghĩa là ta có∂ΩR

div Du( ) = ∆u = 0

Phương trình Laplace được Laplace đưa ra vào khoảng năm 1780

Phương trình Poisson

−∆u = f