energy method の入口いりぐち

date2026-07-15document_iddoc_43a6bec521a85d16ea805c0223b35e4fdescriptionenergy method を、PDE の解の大きさや保存量を積分量で制御する方法として導入する。prerequisitesheat・wave・Laplace 方程式 / 多重積分と変数変換type講義content_typelecturestatusactive

mathpartial-differential-equationsenergy-methodlecture

1導入どうにゅう

この講義こうぎでは、PDE の解かいを点てんごとに追跡ついせきするのではなく、領域全体りょういきぜんたいの積分量せきぶんりょうによって制御せいぎょする方法ほうほうを説明せつめいする。前まえの maximum principle は解かいの最大値さいだいちを点てんごとに制御せいぎょしたが、energy method は二乗積分にじょうせきぶんを通とおして解全体かいぜんたいの大おおきさを評価ひょうかする。

2用語ようごと定義ていぎ

energy methodEnergy method は、解かいや微分びぶんの二乗積分にじょうせきぶんをエネルギーとして定義ていぎし、その時間変化じかんへんかを評価ひょうかする方法ほうほうである。

3方針ほうしん

波動方程式はどうほうていしきでは運動うんどうと歪ひずみのエネルギーが保存ほぞんされる。熱方程式ねつほうていしきでは二乗積分にじょうせきぶんが減衰げんすいする。積分せきぶんによる評価ひょうかは、明示解めいじかいがなくても有効ゆうこうである。

4接続せつぞく

energy method では部分積分ぶぶんせきぶんが反復はんぷくして現あらわれる。多次元たじげんでは、Gauss の発散定理はっさんていりを用もちいて領域積分りょういきせきぶんを境界積分きょうかいせきぶんへ変換へんかんする Green 公式こうしきが部分積分ぶぶんせきぶんに対応たいおうする。境界条件きょうかいじょうけんが境界項きょうかいこうを消去しょうきょする役割やくわりを持もつ。

5計算けいさんの前提ぜんてい

以下いかでは $c > 0$ 、 $κ > 0$ 、 $0 < x < L$ 、 $0 [PARSE ERROR: Undefined("Command(\"le\")")] t [PARSE ERROR: Undefined("Command(\"le\")")] T$ とする。また、積分せきぶんと時間微分じかんびぶんの順序交換じゅんじょこうかんおよび部分積分ぶぶんせきぶんができるだけ十分じゅうぶんに滑なめらかな古典解こてんかいを考かんがえる。ここで古典解こてんかいとは、方程式ほうていしきと条件じょうけんを各点かくてんで満みたすだけの微分可能性びぶんかのうせいを持もつ解かいである。弱解じゃくかいへ拡張かくちょうする場合ばあいは、まず滑なめらかな近似解きんじかいで同おなじ評価ひょうかを得えてから極限きょくげんを取とる、という追加ついかの議論ぎろんが必要ひつようである。

6波動方程式はどうほうていしきの例れい

$u_{t t} = c^{2} u_{x x}$ 、 $u (0, t) = u (L, t) = 0$ とする。エネルギーを

E (t) = \frac{1}{2} \int_{0}^{L} (u_{t^{2}} + c^{2} u_{x^{2}}) d x

と定義ていぎする。 $E^{'} (t)$ を計算けいさんし、部分積分ぶぶんせきぶんと境界条件きょうかいじょうけんを使用しようすると、 $E^{'} (t) = 0$ となる。これは明示解めいじかいを使用しようせずに保存量ほぞんりょうを得える方法ほうほうである。

計算けいさんは次つぎの通とおりである。

E^{'} (t) = \int_{0}^{L} (u_{t} u_{t t} + c^{2} u_{x} u_{x t}) d x

方程式ほうていしき $u_{t t} = c^{2} u_{x x}$ を代入だいにゅうすると、

E^{'} (t) = c^{2} \int_{0}^{L} (u_{t} u_{x x} + u_{x} u_{x t}) d x

である。第一項だいいっこうを部分積分ぶぶんせきぶんすると、

\int_{0}^{L} u_{t} u_{x x} d x = {[u_{t} u_{x}]}_{0^{L}} - \int_{0}^{L} u_{t x} u_{x} d x

となる。残のこる二ふたつの積分せきぶんが打うち消けし合あうので、まず

E^{'} (t) = c^{2} {[u_{t} u_{x}]}_{0^{L}}

を得える。この境界項きょうかいこうは、端点たんてんを通とおるエネルギーの流出入りゅうしゅつにゅうを表あらわす。固定端条件こていたんじょうけん $u (0, t) = u (L, t) = 0$ を $t$ で微分びぶんすると端点たんてんで $u_{t} = 0$ となるため、境界項きょうかいこうが消きえて $E^{'} (t) = 0$ である。保存ほぞんは方程式ほうていしきだけの性質せいしつではなく、境界きょうかいからエネルギーが出入でいりしないことにも依存いぞんする。

7波動方程式はどうほうていしきの一意性いちいせい

同おなじ初期条件しょきじょうけんと固定端条件こていたんじょうけんを満みたす二ふたつの解かい $u, v$ があるとする。 $w = u - v$ と置おけば、 $w$ は斉次波動方程式せいじはどうほうていしきと零れいの初期条件しょきじょうけん・固定端条件こていたんじょうけんを満みたす。したがって $w$ の初期しょきエネルギーは 0 であり、エネルギー保存ほぞんからすべての $t$ で

\int_{0}^{L} (w_{t^{2}} + c^{2} w_{x^{2}}) d x = 0

となる。古典解こてんかいの仮定かていにより被積分関数ひせきぶんかんすうは連続れんぞくかつ非負ひふなので、積分せきぶんが 0 ならすべての点てんで $w_{t} = w_{x} = 0$ である。よって $w$ は時間じかんにも空間くうかんにも依存いぞんしない定数ていすうであり、固定端条件こていたんじょうけんから $w = 0$ 、すなわち $u = v$ である。

8熱方程式ねつほうていしきの例れい

$u_{t} = κ u_{x x}$ 、 $u (0, t) = u (L, t) = 0$ とする。 $\int_{0}^{L} u^{2} d x$ を時間微分じかんびぶんすると、

\frac{d}{d t} \int_{0}^{L} u^{2} d x = - 2 κ \int_{0}^{L} u_{x^{2}} d x [PARSE ERROR: Undefined("Command(\"le\")")] 0

を得える。熱ねつでは二乗積分にじょうせきぶんが減衰げんすいし、波なみではエネルギーが保存ほぞんされる。この差さが方程式ほうていしきの性質せいしつを示しめす。

導出どうしゅつは

\frac{d}{d t} \int_{0}^{L} u^{2} d x = 2 \int_{0}^{L} u u_{t} d x = 2 κ \int_{0}^{L} u u_{x x} d x

から始はじまる。部分積分ぶぶんせきぶんにより

2 κ \int_{0}^{L} u u_{x x} d x = 2 κ {[u u_{x}]}_{0^{L}} - 2 κ \int_{0}^{L} u_{x^{2}} d x

である。Dirichlet 条件じょうけんでは境界項きょうかいこうが 0 になるため、減衰評価げんすいひょうかが得えられる。

この評価ひょうかだけで二乗積分にじょうせきぶんが増加ぞうかしないことは分わかる。さらに、両端りょうたんで 0 になる関数かんすうに対たいする Poincare 不等式ふとうしき

\int_{0}^{L} u^{2} d x [PARSE ERROR: Undefined("Command(\"le\")")] {(\frac{L}{π})}^{2} \int_{0}^{L} u_{x^{2}} d x

をここでは用もちいると、減衰げんすいの速はやさまで得えられる。Poincare 不等式ふとうしきは、関数かんすうの大おおきさを、その変化率へんかりつの大おおきさで制御せいぎょする不等式ふとうしきである。

この定数ていすうは、変数分離法へんすうぶんりほうで導入どうにゅうした sine 固有関数こゆうかんすうから確認かくにんできる。 $u$ を

u (x) = \sum_{n = 1}^{\infty} b_{n} \sin (\frac{n π x}{L})

と展開てんかいすると、直交性ちょっこうせいにより

\int_{0}^{L} u^{2} d x = \frac{L}{2} \sum_{n = 1}^{\infty} | b_{n} |^{2}, \int_{0}^{L} u_{x^{2}} d x = \frac{L}{2} \sum_{n = 1}^{\infty} {(\frac{n π}{L})}^{2} | b_{n} |^{2}

である。 $n^{2} [PARSE ERROR: Undefined("Command(\"ge\")")] 1$ なので $\int_{0}^{L} u_{x^{2}} d x [PARSE ERROR: Undefined("Command(\"ge\")")] (π / L)^{2} \int_{0}^{L} u^{2} d x$ となり、上うえの Poincare 不等式ふとうしきが従したがう。この導出どうしゅつは、十分じゅうぶんに滑なめらかで両端りょうたんで 0 になる $u$ の sine 展開てんかいを使用しようしている。

この不等式ふとうしきを減衰評価げんすいひょうかへ代入だいにゅうすれば

\frac{d}{d t} \int_{0}^{L} u^{2} d x [PARSE ERROR: Undefined("Command(\"le\")")] - 2 κ {(\frac{π}{L})}^{2} \int_{0}^{L} u^{2} d x

となる。したがって

\int_{0}^{L} u (x, t)^{2} d x [PARSE ERROR: Undefined("Command(\"le\")")] e^{- 2 κ (π / L)^{2} t} \int_{0}^{L} u (x, 0)^{2} d x

であり、零境界条件れいきょうかいじょうけんの熱ねつは二乗積分にじょうせきぶんの意味いみで指数的しすうてきに減衰げんすいする。ここに現あらわれる $π / L$ は、変数分離法へんすうぶんりほうで最初さいしょの sine 固有関数こゆうかんすうに対応たいおうした波数はすうである。

9外力がいりょくがある場合ばあい

比較例ひかくれいとして、同おなじ固定端条件こていたんじょうけん $u (0, t) = u (L, t) = 0$ のもとで $u_{t t} = c^{2} u_{x x} + F (x, t)$ を考かんがえる。この場合ばあい、境界流束きょうかいりゅうそくが消きえるので

E^{'} (t) = \int_{0}^{L} F (x, t) u_{t} (x, t) d x

となり、エネルギーは保存ほぞんされない。右辺みぎへんは外力がいりょくが系けいへ注入ちゅうにゅうする仕事率しごとりつを表あらわす。ここで

∥ F (\cdot, t) ∥_{2} = {(\int_{0}^{L} | F (x, t) |^{2} d x)}^{1 / 2}

と定義ていぎする。Cauchy--Schwarz 不等式ふとうしきと $∥ u_{t} (\cdot, t) ∥_{2^{2}} [PARSE ERROR: Undefined("Command(\"le\")")] 2 E (t)$ から

E^{'} (t) [PARSE ERROR: Undefined("Command(\"le\")")] ∥ F (\cdot, t) ∥_{2} \sqrt{2 E (t)}

を得える。 $E (t) > 0$ の区間くかんで $\sqrt{E (t)}$ を微分びぶんし、 $∥ F (\cdot, t) ∥_{2}$ が $0 [PARSE ERROR: Undefined("Command(\"le\")")] t [PARSE ERROR: Undefined("Command(\"le\")")] T$ で積分可能せきぶんかのうであるとすれば、

\sqrt{E (t)} [PARSE ERROR: Undefined("Command(\"le\")")] \sqrt{E (0)} + \frac{1}{\sqrt{2}} \int_{0}^{t} ∥ F (\cdot, s) ∥_{2} d s

となる。 $E (t) = 0$ となる時刻じこくを含ふくむ場合ばあいも、 $\sqrt{E + ε}$ で評価ひょうかして $ε ↓ 0$ とすれば同おなじ結論けつろんを得える。これは外力がいりょくが小ちいさければ解かいのエネルギーも制御せいぎょされるという安定性あんていせいの評価ひょうかである。energy method は保存ほぞんだけでなく、増加ぞうかや減衰げんすいを不等式ふとうしきで制御せいぎょする枠組わくぐみでもある。

10どこまで成立せいりつするか

非線形ひせんけい PDE では、評価ひょうかしたいエネルギーだけで不等式ふとうしきの右辺うへんを制御せいぎょできるとは限かぎらない。散逸さんいつ・保存ほぞん・外力がいりょくの有無うむを確認かくにんする必要ひつようがある。

energy method は解かいの式しきを構成こうせいしなくても、一意性いちいせい・減衰げんすい・安定性あんていせいを示しめせる。次つぎの Green 関数かんすうでは、線形せんけい PDE の解かいそのものを点源応答てんげんおうとうの積分せきぶんとして構成こうせいする。