題目
問題 10.
(a) 考慮熱傳導方程在 J = ( ? ∞ , ∞ ) J = (-\infty, \infty) J=(?∞,∞) 上,證明“能量”
E ( t ) = ∫ J u 2 ( x , t ) d x E(t) = \int_{J} u^{2}(x,t) dx E(t)=∫J?u2(x,t)dx
(8)
不增加;進一步證明,除非 u ( x , t ) = 常數 u(x,t) = \text{常數} u(x,t)=常數,否則它確實減少。
(b) 考慮熱傳導方程在 J = ( 0 , l ) J = (0,l) J=(0,l) 上,帶有 Dirichlet 或 Neumann 邊界條件,證明 E ( t ) E(t) E(t) 不增加;進一步證明,除非 u ( x , t ) = 常數 u(x,t) = \text{常數} u(x,t)=常數,否則它確實減少。
? 考慮熱傳導方程在 J = ( 0 , l ) J = (0,l) J=(0,l) 上,帶有 Robin 邊界條件:
u x ( 0 , t ) ? a 0 u ( 0 , t ) = 0 , u_{x}(0,t) - a_{0}u(0,t) = 0, ux?(0,t)?a0?u(0,t)=0,
u x ( L , t ) + a L u ( L , t ) = 0. u_{x}(L,t) + a_{L}u(L,t) = 0. ux?(L,t)+aL?u(L,t)=0.
(9)
如果 a 0 > 0 a_{0} > 0 a0?>0 和 a L > 0 a_{L} > 0 aL?>0(注:原文中 a l a_{l} al? 應為 a L a_{L} aL?,已修正),證明端點對 E ( t ) = ∫ 0 L u 2 ( x , t ) d x E(t) = \int_{0}^{L} u^{2}(x,t) dx E(t)=∫0L?u2(x,t)dx 的減少有貢獻。這被解釋為部分能量在邊界處損失,因此我們稱邊界條件為輻射或耗散的。
提示. 為了證明 E ( t ) E(t) E(t) 的減少,考慮其對 t t t 的導數,用 k u x x k u_{xx} kuxx? 替換 u t u_{t} ut?,并進行分部積分。
備注 3.P.1. 在熱傳導(或擴散)方程的情況下,由 (8) 給出的能量更多是數學上的構造。
解決題目
熱傳導方程的標準形式為:
u t = k u x x , u_t = k u_{xx}, ut?=kuxx?,
其中 k > 0 k > 0 k>0 是熱擴散系數(常數)。能量泛函定義為:
E ( t ) = ∫ J u 2 ( x , t ) d x , E(t) = \int_{J} u^2(x,t) dx, E(t)=∫J?u2(x,t)dx,
其中 J J J 是定義域。證明的核心是計算 d E d t \frac{dE}{dt} dtdE? 并利用熱傳導方程和邊界條件。
(a) 無限域 J = ( ? ∞ , ∞ ) J = (-\infty, \infty) J=(?∞,∞)
證明:
計算 E ( t ) E(t) E(t) 的時間導數:
d E d t = d d t ∫ ? ∞ ∞ u 2 ( x , t ) d x . \frac{dE}{dt} = \frac{d}{dt} \int_{-\infty}^{\infty} u^2(x,t) dx. dtdE?=dtd?∫?∞∞?u2(x,t)dx.
假設 u u u 足夠光滑,且當 ∣ x ∣ → ∞ |x| \to \infty ∣x∣→∞ 時, u u u 及其一階導數趨于零(物理上合理的衰減條件),因此可以交換積分與導數順序:
d E d t = ∫ ? ∞ ∞ ? ? t ( u 2 ) d x = ∫ ? ∞ ∞ 2 u u t d x . \frac{dE}{dt} = \int_{-\infty}^{\infty} \frac{\partial}{\partial t} (u^2) dx = \int_{-\infty}^{\infty} 2u u_t dx. dtdE?=∫?∞∞??t??(u2)dx=∫?∞∞?2uut?dx.
由熱傳導方程 u t = k u x x u_t = k u_{xx} ut?=kuxx?,代入得:
d E d t = 2 k ∫ ? ∞ ∞ u u x x d x . \frac{dE}{dt} = 2k \int_{-\infty}^{\infty} u u_{xx} dx. dtdE?=2k∫?∞∞?uuxx?dx.
對積分進行分部積分:令 f = u f = u f=u, g ′ = u x x g' = u_{xx} g′=uxx?,則 f ′ = u x f' = u_x f′=ux?, g = u x g = u_x g=ux?,
∫ ? ∞ ∞ u u x x d x = [ u u x ] ? ∞ ∞ ? ∫ ? ∞ ∞ u x ? u x d x . \int_{-\infty}^{\infty} u u_{xx} dx = \left[ u u_x \right]_{-\infty}^{\infty} - \int_{-\infty}^{\infty} u_x \cdot u_x dx. ∫?∞∞?uuxx?