chương 1 báo cáo bài toán biên của phương trình truyền nhiệt. Phương trình truyền nhiệt dasjlaksjgapghawkljdwiofhogewhaaw;lkgjawoiehtawklngas,.gnwaoiwutawegahsldkgnalsdgnaodgjwaiwa;gknamsdna;vbahnowiewgawgas;lkdgasgnw,and;ghươawiouioaghadlkgnawoegawgwal;keghawlekghawoieuwaetoiwahgds;lkgha;dhgoiwaeghadklghadlkghaskjghwiaugdklhhhhhhhhaa;doghs;g
I)Giới thiệu !!"# !$%&'(!)&*+,&-(,./012 !!"#!$3456(2!&-7'& (#!$%89: 3006 ;<$;=>?@(3>%6AB30060 C ;<$D/E0(/FE9&0G$ %,. /&-HD)IJ"!!!!"8 K L(0M&()I=@!N!;O/E()ICP )+!D,Q&-RC0 M,./(,&-/EBJ 7!!"#!$0 9:)& *+!:2(& 7!!"#!$ "-%2+-%8(S%2+-%2+-%? -3)+)&?6%)+%,8J-%9-% -%,*/9$,8%!)%&HA(,8%-JT0 )I!:U2(D&9:(V)I J!"8(W!$ 0/!X7 F0 W0Y+8 ;<$;=>"9+!0MQJ (9&)+D2R'!"T0Z;9[( 0\R)+9[#)+'R!$T"/& H2QR2!$TJD)&!$T$UH0 MQJJ9[)+,9[!"T,X9 D,D0M&(2!"T9&29D#]() ],%H)IQJR!"T"/UH0Y+ /J-(0GD;<$C 7@ %)I^ 7_J!9[(HH$ -U)+0)` 7R)+9[U<% a9]/RK)+&3A( ,8)6()+9[0 .2)Sự truyền nhiệt trên thanh đồng chất Mật độ năng lương nhiệt%;<$>N2H98;a< ,&ObHU`&>3#>cd>ce6 0F00&QJ V?D8(2;T;(-H(/2QJ 0 e(x,t) = thermal energy density f)T=QJ2?D8(,XOg(2 0)J+D&(&(98;aR 0Y+!:2(&>(7U_J!,XJ(H 2Sg/2QJO`&!=0 Z0F0hU)+(&( 2<A Năng lượng nhiệt0>N2<A #?8>?8 >ij> 0F00M/2QJ (,XO &(;2%OJ ;a;=8R/2Q( ,J0+%/2J(O0 j>(RA `3>%6HDk&(,XO &(9&H Heat energy = e(x,t)A∆x, &H,JRa<(bj> Sự bảo toàn năng lượng. MQJ&&'>>ij> O`&_^9&)+`&; Ra<3>(j>6(Q J'&; &a<0lQJ ;a; (,XO 9&H)T;a; J0mJX;K ! Tốc độ thay đổi nhiệt lượng theo thời gian = nhiệt lượng chạy qua biên mỗi đơn vị thời gian + nhiệt lượng được tạo ra bên trong trên một đơn vị thời gian G"J-(sự bảo toàn nhiệt lượng. GUa<A%2 OJ`&_;= &H!"J9]>(?0 Thông lượng nhiệt0M J& -H(XJ Φ(x,t) = thông lượng nhiệt(Nhiệt lượng mỗi đơn vị thời gian chảy trên một đơn vị thể tích) Mn3>%6od*(J;?0MJ; 2 ?_Ra<(n3>%6bpn3>i%6b%XJ;= J 2?D8"U98a<0Mn3>%6qd( n3>i%6qd 0F0%J 2?_'>( QJ a<%&,J'>i0 Nhiệt năng.-OJ; & Q(x,t) = nhiệt lượng chảy trên một đơn vị thể tích trong một đơn vị thời gian m3>%6,XO<RA(OJ)&2? _(2,&;=m3>%6b0 Sự bảo toàn nhiệt lượng( trong lát cắt mỏng).2; Q(9&QU(S2;2 0F0(,X8>;K?; &, ??,XO<RA0U,8'&(%C)I 8@EK;&'&(%0M/% &j>US)H8U',J \98a<;R;A0CX&=,((8 >3,XH2VA(&6(9&Hk&(k9&X7r30F06 ;=9c&X730F0F60\8U'(%%(?04&H# ?*R!"#!$ Sự bảo toàn nhiệt lượng( chính xác)0Y+;='&(R! ;&&(J,X^>U_J!<A0>`>N 2)@'&' 3#>c>c;60C)I`7 )+;&&(J&&'(0OJ;=%(H^ O9&J'; 3>c(>c;6(QJ) &(9&H3),>H;A=)b6 la,./>&X730F0s6,^!:2(&( ,X!(&>0 ,(;(=)3(` :60G(C,&8!"'&( _J@,>?%(9&HHJ;=2!"#!$0 1-U'&X730F0s6(28!"X_H 4&H 8!"(;=dU-(;]Eg/9UR_&!Kd0 G(C^,!&98!";=d04&HH 30F0s6%`&?/;&&(8!"%(9';R9' !"&X730F0t6030F0t6J!&_J! X_,H;( :0 l89:%qdUu>u;%XJn(2(QR>0 f?RH';U'3U,;q604&HJ;K @&'>c(>c;%R9#&!30F0t6 M2(9 0X/;=2RH7,& !;=/2QJ0ZB!";,2 iG"(2?E;R80 iZ$;=7R,(,X:D0f)T=v3>6( :(cdU%]E0$7v3>6cdU->0HD7 ;=!7S'2D)&&v3Mv3d(v3>6( : S'@D"/)&&v3>6]90-(;&,& ((9&Hwd,v3>6(0G("[Ucd(9& HH(,XDUv3wd0l(&@@Qxdd(/EBH ;?T8>&!N!&)J,9]DQ JRF/,#2! &0G(yA)+_ R,9 3y-(X)6 c9 3J$D!&2?& JH6 O%#+3(?*69 R2/!: 2(&20l89:%JDQ2?,J# )I ,)&UJDQ#&]20l& H9 !:2(&2R/(2;(&!7'!0 3;(/!0F0)IDL(60X_U,&2;? '9 ,X!:2(&20 %T '&U@/,$?J, D(H 0l>"9+2X78>&- &H(!$RHDO#?8(?8 ,%9 )I!2(&>%c3>60&R S%&(_J!)I&9(=)0 +%$ $J/!&)(3z) ,6]2?=)>!>^0 Nhiệt năng0MQ&2a<A;=`3>%6bj>0 % HzJk&(JQJ$DQ2R[/# 2+3>%60\9 ,X!:2(&2% QJ 2?,J;=3>63>%60$U, /2,!3>6 !3>6c/2,3,J 2?D86 &?>O%HD(9&S(,XS9'0O ,JRa<A(!bj>0OQ V<A;= 3>63>%6!bj>%9&H `3>%6bj>c3>63>%6!bj>0 {=(CB8)+ ,;@Q(2 `3>%6c3>6!3>63>%630F0|6 Q 2?,J;=Q 2?,J (2?2_ 2/2,J3 , J 2?D86\/2Q;?&';A;=)T9:X7 30F0|6%?/;&SQJ%30F0W630F0t6%K( Định luật fourier 30F0x6(2!HF;%2 3>%6(XJ3J 2?98 2? _6n3>%60')&((&DJ}MH, $>"9+2;D7;D9~)+!:2R9yJ(& _20G$ CH<2(8R9y (CB; 1. M2(=)&2(&H,XH9y 0 2. MS')+ 2J#]H2 &>]H2! 3. Y+ (U3&]2/6H( 4. Y+)IOU@/,%,D,CH ])+ 20 •&`3x|€•€Wd6B#/s8 (H<';=X7 3HJ_+6 GJ;/v&`)+9[0K"()+; R 2g(2OR23(2(R>&_ ?6g30F0€6&XJ‚/U)+ 23 2?29(60M2>z% (JC#30F0€60 >?)70Z)(&,Q9[R/ (J-()9[0+&=/, ,Q9[RCz,g!:2(&#&' /0+'(%9yURJ,K]2, 20M@/H)!9[,N0lU2J (#&'/)I((R>0Z@+&=, Q9[R$/,,K2,0 U9 )+!:2(&2_,X -&04&H&!$()9[^!: 2(&>0+_)T9:H(=)0 Phương trình nhiệt0`&?Ev&`30F0€6%(&! ;&&(J30F0x6%J!!") _*=SmBJ&U(!$,&;^( 23>%60l()%!%!:2(&/((R>0& _J!a;%%!% là hằng số, phương trình vi phân từng phần trở thành: f)TS;=d),&=)!!!"_ ( &H,(=)J-(29[(;= 30F0d6_J-(!%&HS ;=d(8(,XO0MJ;$/!(&2 V%30F0d6)IX7J&(%/8( J-()+,0M'J/E,&()+ &(zJ&7+0l/30F0d6zJ ;U !ƒ&,0l89:S2H3>%638 9:U&"X~6AB!,30F0€6 &2)_J!0 Điều kiện ban đầu0!"#!$X)+, RJ%30F0„630F0d6%('&(`&_0\! !"H'&('2D%?$ ;&S !!"U,;$0G/M`ƒ&)+O?8 R>R2''2!!";/F%c+0MH '&(!F0l?;$;&S! !"#!$UF,;$%?8;$>(/;$0# @,H%!!"U,;$%HD& JD2R'`&U`&>0(JHU !!"R%&J2&0#! '&(2;%)I&H2,;$3…63_'cd6% 2;$0f?(,X!(=)%H!:2(&>04& HB()2R?2;$0 3>%d6cv3>6 M(BRD9+;&2)2,&_}; =)+!"!2;$()+O2`&! 30F0„630F0d60 %$;?'F_D>cd(>ce0M ,X;JH(,XD&J0,($!:2 (&'&(;/J &30F0„6(30F0d60)I 7 ,((&!$)0 1.3 Điều kiện biên \!%30F0„630F0d6%2,; 3{6($ 0G,8J!!:2(&/E'_D0_ ,?!:2(&/; &(; &(0GD &-!7'!%)#=B;X_; &(%(,XO,D0 Quy định về nhiệt độ0&2)_J!2R%8 9:'>cd%&D>!>^2? 3d%6c3630W06 &H36(2R/7 UV0 fU'0&_J!,%HD?9y 7,X!2 \n36&U0G(U2,&'& ($ %%'>cd0G29J&'>cd030W0F6,XD 8!"'>U^;29'2?R>0l89: &?,; (&(&(3X,C;A †)+&(&‡60&_J!(,XH)+'; 0M '>cd G?/('R`ƒ&0\2!>CUA 203`0ˆ%,X,86%)H,X/(&9yy8J!0 l89:BKJ2!>CU,X,8'0M)I& ,A((H,X,80\X,8)IJ0Y+ `&(J-()+04&H,X,8)IH,K$ 0Z@"(2!7'!g2,X,8)IO`& ,&)&U08&%',8J!% J&(^/U)+ 2@(2 ; &(0G,; (J-(?/`ƒ&('0MH AB'>cd &H=)ZJ-()3-()60G ,; )+,J!8R0MH; ‰3d%6q%(J_&('>cd04&H) ; (&_J!(("0GH(89&')&&9" &X730W0s63UZqd60\?3d%6o)IJ@ ,/0,D;9&!30W0s6(? 3d%6q0'>cd2; !)IH(,Š-=2; !)I!:Q04&H!9'>cd030W0s6,X" [UE/(0&;(/!0W0;'!X7%&]2 ,%&?/('R`ƒ&'D>ce( &H(2; &('>ce0C/!7@9 )+,9@; ; 30W0s6(; ; !30W0t60 [...]... T(x,y,z,t) ở mọi chỗ trên biên, khi đó T là một hàm đã biết theo biến t tại mỗi điểm o biên Điều đó có thể là sự truyền dọc theo biên được quy định Thông thường, ta có biên ( hay một phần của biên) được cách nhiệt nó có nghĩa là không có dòng nhiệt tại phần này của biên Khi vecto thông lượng nhiệt là -∇u, nhiệt lượng truyền ra ngoài là đơn vị hướng ngoại của vecto thành phần của vecto truyền nhiệt, -∇u · , trong... suy ra vecto thông lượng nhiệt Φ tỉ lệ với độ dốc của nhiệt độ (∇u ): Biết như Định luật Fourier về truyền nhiệt, trong đó là tính dẫn nhiệt Do đó trong 3 chiều ∇u thay cho Phương trình nhiệt khi vecto thông lượng nhiệt ,(1.5.7), thay thế trong phương trình bảo toàn nhiệt năng,(1.5.6), vi phân từng phần của nhiệt độ bằng Trong trường hợp không có nguồn nhiệt (Q = 0) và hệ số nhiệt là hằng số, (1.5.8)... điều kiện biên xảy ra ở mỗi biên Điều này là không quan trọng nếu cả hai điều kiện đáp ứng được những điều kiện biên giống nhau Ví dụ, điều này hợp lí khi x = 0 tồn tại sự dao động của nhiệt độ u(0,t) = 100 – 25cost và với x = L, trở thành cách nhiệt 1.4 Trạng thái cân bằng nhiệt 1.4.1 Quy định về nhiệt độ Bây giờ chúng ta xây dựng một phương trình truyền nhiệt đơn giản,nhưng tiêu biểu Nếu hệ số nhiệt. .. thái cân bằng thỏa mãn Được gọi là phương trình Poisson Nếu,trong phép cộng, không có nguồn nhiệt ( Q=0) thì Phép biến đổi laplace của phân phối nhiệt độ bằng 0 Phương trình (1.5.14) được gọi là phương trình laplace Nó còn được gọi là phương trình thế, khi lực hút và lực tĩnh điện thế thỏa mãn (1.5.14) nếu không có nguồn nhiệt ta sẽ làm sang tỏ một số vấn đề về Phuong trình la place ở chương sau Vấn... nhiệt độ tăng nhanh nhất theo chiều của gradient Khi sự truyền nhiệt theo chiều giảm của nhiệt độ, vecto truyền nhiệt sẽ ngươc chiều với gradient nhiệt. Theo đó Khi | | bằng độ lớn của tỷ lệ thay đổi của u ( chiều của gradient) Hay nó còn gọi là định luật của Fourier về dẫn nhiệt Như vậy, trong ba chiều, gradient thay cho Một tính chất cơ bản khác của gradient là pháp tuyến ( vuông góc ) với mặt phẳng... lượng nhiệt và vecto bình thường Ta cần một sự diễn tả cho sự truyền nhiệt Trong vấn đề một chiều thông lượng nhiệt Φ hướng sang bên phải (Φ < 0 nghĩa là hướng sang bên trái ) Trong vấn đề 3D nhiệt truyền theo nhiều hướng và do đó thông lượng nhiệt là một vecto Φ Độ lớn của Φ số lượng nhiệt lượng truyền ra một đơn vị diện tích trong một đơn vị thời gian Tuy nhiên khi xét sự bảo tồn nhiệt lượng chỉ có nhiệt. .. thái ổn định) nhiệt độ và sự truyền nhiệt tại x = 0 1.4.6 Hai đầu thanh có độ dai L được cách nhiệt Được mắc vào một nguồn nhiệt ≠ 0 không đổi có nhiệt độ ban đầu u(x,0) = f(x) (a) Bằng toán học hãy chỉ ra không tồn tại sự phân phối nhiệt độ tại trạng thái cân bằng Giải thích ngắn gọn bằng vật lý (b) Tính tổng nhiệt năng trong thanh 1.4.7 Với vần đề truyền nhiệt, xác định sự phân phối nhiệt độ tại trạng... được gọi là độ khuếch tán nhiệt Từ định nghĩa , ta tính sự chênh lệnh độ dốc của u Biểu thức được định nghĩa là toán tử laplace của u Do đó trong trường hợp này Phương trình (1.5.11) thường được gọi là phương trình nhiệt hay phương trình khuếch tán trong không gian ba chiều Ký hiệu thường được dùng để nhấn mạnh ∇: Lưu ý rằng ∇u là ∇ hoạt động tại u, khi là một vecto chấm sản phẩm của ∇ với A Hơn nữa là... toàn nhiệt năng Tại mỗi điểm số lượng nhiệt năng truyền ra ngoài khoảng R trên một đơn vị diện tích trong một đơn vị thời gian là thành phần bên ngoài bình thường của vecto thong lượng nhiệt Tại điểm B trong hình 1.5.2 thành phần hướng ngoại của vecto thông lượng nhiệt là |Φ|cosθ = Φ· = Φ · Nếu vecto thông lượng hướng vào trong thì Φ · < 0 và nhiệt lượng truyền ra ngoài là âm Để tính tổng nhiệt lượng truyền. .. điều kiện biên và bất kì nguồn nhiệt nào phụ thuộc vào thời gian, điều này đúng rằng tồn tại cách giải trạng thái cân bằng cho phương trình nhiệt thỏa mãn điều kiện biên ổn định đã cho: Lưu ý rằng sự phân phối nhiệt độ ở trạng thái cân bằng u(x,y,z) thỏa mãn một phương trình vi phân từng phần khi nhiều hơn khoảng một chiều tham gia Trong trường hợp tính chất nhiệt không đổi, sự phân phối nhiệt độ ở . độ thay đổi nhiệt lượng theo thời gian = nhiệt lượng chạy qua biên mỗi đơn vị thời gian + nhiệt lượng được tạo ra bên trong trên một đơn vị thời gian G"J-(sự bảo toàn nhiệt lượng 7R)+9[U<% a9]/RK)+&3A( ,8)6()+9[0 .2)Sự truyền nhiệt trên thanh đồng chất Mật độ năng lương nhiệt% ;<$>N2H98;a< ,&ObHU`&>3#>cd>ce6. GUa<A%2 OJ`&_;= &H!"J9]>(?0 Thông lượng nhiệt0 M J& -H(XJ Φ(x,t) = thông lượng nhiệt( Nhiệt lượng mỗi đơn vị thời gian chảy trên một