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 の解かいそのものを点源応答てんげんおうとうの積分せきぶんとして構成こうせいする。

Introduction to Energy Methods

1Introduction

This lecture explains how integral quantities over an entire domain can control a PDE solution without tracking it pointwise. The preceding maximum principle controls pointwise extrema, whereas an energy method evaluates the overall magnitude of a solution through quadratic integrals.

2Terminology and Definition

An energy methodEnergy method defines an integral involving squares of a solution and its derivatives as an energy and evaluates its evolution in time.

3Strategy

For the wave equation, kinetic and strain energy are conserved. For the heat equation, a quadratic integral decays. These integral estimates remain useful even when an explicit solution is unavailable.

4Connection to Earlier Material

Energy methods repeatedly use integration by parts. In several dimensions, Green's identities are the corresponding formulas obtained from Gauss's divergence theorem by converting domain integrals into boundary integrals. Boundary conditions determine whether the resulting boundary terms vanish and therefore form part of the estimate rather than an incidental assumption.

5Assumptions for the Calculations

Throughout this lecture, let $c > 0$ , $κ > 0$ , $0 < x < L$ , and $0 [PARSE ERROR: Undefined("Command(\"le\")")] t [PARSE ERROR: Undefined("Command(\"le\")")] T$ . We consider classical solutions sufficiently smooth to interchange integration and time differentiation and to apply integration by parts. A classical solution has enough differentiability to satisfy the equation and its conditions pointwise. Extending the estimates to weak solutions requires an additional approximation-and-limit argument.

6Example: The Wave Equation

Consider

u_{t t} = c^{2} u_{x x}, u (0, t) = u (L, t) = 0.

Define the energy by

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

Differentiation, integration by parts, and the boundary conditions give $E^{'} (t) = 0$ . The calculation begins with

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

Substituting $u_{t t} = c^{2} u_{x x}$ yields

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

Integration by parts in the first term gives

\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 .

The two remaining integrals cancel, so

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

This boundary term represents energy flux through the endpoints. Differentiating the fixed-end conditions in time gives $u_{t} = 0$ at both endpoints, so the boundary term vanishes and $E^{'} (t) = 0$ . Conservation therefore depends both on the equation and on the absence of energy flux through the boundary.

7Uniqueness for the Wave Equation

Suppose that two solutions $u$ and $v$ have the same initial data and fixed-end boundary conditions. Their difference $w = u - v$ satisfies the homogeneous wave equation with zero initial and boundary data. Its initial energy is zero, and conservation implies

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

for every $t$ . Under the classical-solution assumptions, the integrand is continuous and nonnegative; hence $w_{t} = w_{x} = 0$ at every point. Thus $w$ is constant in both time and space, and the fixed-end conditions force $w = 0$ . Therefore $u = v$ .

8Example: The Heat Equation

Consider

u_{t} = κ u_{x x}, u (0, t) = u (L, t) = 0.

Differentiating the squared norm gives

\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 .

Integration by parts yields

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 .

The Dirichlet conditions eliminate the boundary term, so

\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.

Thus the quadratic integral decays for heat flow, whereas the wave energy is conserved.

This estimate alone proves only that the quadratic integral is nonincreasing. To obtain a decay rate, use the Poincare inequality for functions that vanish at both endpoints:

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

The inequality controls the magnitude of a function by the magnitude of its derivative. Its constant follows from the sine eigenfunctions introduced in separation of variables. Expanding

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

and using orthogonality gives

\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.}

Since $n^{2} [PARSE ERROR: Undefined("Command(\"ge\")")] 1$ , it follows that $\int_{0}^{L} u_{x^{2}} d x [PARSE ERROR: Undefined("Command(\"ge\")")] (π / L)^{2} \int_{0}^{L} u^{2} d x$ , which is the stated Poincare inequality. This derivation uses the sine expansion of a sufficiently smooth function that vanishes at both endpoints.

Substituting this inequality into the decay estimate gives

\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,

and therefore

\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 .

Hence heat flow with homogeneous boundary data decays exponentially in the quadratic-integral sense. The factor $π / L$ is the wave number of the first sine eigenfunction in separation of variables.

9The Case with External Forcing

For comparison, consider

u_{t t} = c^{2} u_{x x} + F (x, t), u (0, t) = u (L, t) = 0.

The boundary flux still vanishes, but now

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

The right-hand side is the rate at which the external force performs work on the system. Define

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

The Cauchy--Schwarz inequality and $∥ u_{t} (\cdot, t) ∥_{2^{2}} [PARSE ERROR: Undefined("Command(\"le\")")] 2 E (t)$ imply

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

On intervals where $E (t) > 0$ , differentiating $\sqrt{E (t)}$ gives, provided $∥ F (\cdot, t) ∥_{2}$ is integrable on $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 .

If $E$ vanishes at some time, the same result follows by estimating $\sqrt{E + ε}$ and then letting $ε ↓ 0$ . This stability estimate shows that small forcing produces a controlled energy response. Energy methods therefore describe not only conservation but also growth and decay through inequalities.

10Scope of Validity

For nonlinear PDEs, the energy being estimated may not control every term on the right-hand side of the resulting inequality. One must determine how dissipation, conservation, forcing, and nonlinear terms interact before closing an estimate.

Energy methods can establish uniqueness, decay, and stability without constructing an explicit solution. The next lecture instead constructs solutions of linear PDEs as integrals of point-source responses using Green functions.