markdown
PDE の初期値問題と境界値問題md 9405466
lecture/math/partial-differential-equations/initial-and-boundary-value-problems.lecture.n.md
Download PDF

PDE の初期値問題しょきちもんだい境界値問題きょうかいちもんだい

date2026-07-16document_iddoc_e802fb3fa56b74151b415f16011b8a3cdescriptionPDE の初期値問題・境界値問題・初期境界値問題を区別し、境界条件、適切性、適合条件まで含めて問題設定を判定する。prerequisitesPDE とは何か / 勾配・発散・回転 / Green・Gauss・Stokes の定理type講義content_typelecturestatusactiverelateddata/lecture/math/partial-differential-equations/introduction-to-pdes.lecture.n.md / data/lecture/math/partial-differential-equations/method-of-characteristics.lecture.n.md / data/lecture/math/partial-differential-equations/classification-of-second-order-linear-pdes.lecture.n.md / data/lecture/math/partial-differential-equations/heat-wave-and-laplace-equations.lecture.n.md / data/lecture/math/partial-differential-equations/introduction-to-energy-methods.lecture.n.md / data/lecture/math/vector-calculus/green-gauss-and-stokes-theorems.lecture.n.md
mathpartial-differential-equationsboundary-value-problemslecture

1導入どうにゅう

PDE の問題もんだいは、方程式ほうていしきだけではさだまらない。未知関数みちかんすうかんがえる領域りょういきと、どの場所ばしょにどの条件じょうけんあたえるかをわせて、はじめてひとつの問題もんだいになる。

この講義こうぎでは、初期値問題しょきちもんだい境界値問題きょうかいちもんだい初期境界値問題しょききょうかいちもんだい区別くべつし、条件じょうけんおおいかすくないかだけでなく、存在そんざい一意性いちいせい安定性あんていせい適合条件てきごうじょうけんまで確認かくにんする。

2問題もんだい集合しゅうごう

ΩRn連結れんけつ空間領域くうかんりょういきΩ をその境界きょうかいT>0終端時刻しゅうたんじこくとする。この講義こうぎでは、Ω じょう外向そとむ単位法線たんいほうせん発散定理はっさんていり使用しようできる程度ていど境界きょうかいなめらかであると仮定かていする。時間発展じかんはってんあつかうときは、時空領域じくうりょういき

QT=Ω×(0,T)

u(x,t)かんがえる。t=0断面だんめん Ω×{0}初期面しょきめんΩ×(0,T)側面境界そくめんきょうかいである。初期条件しょきじょうけん境界条件きょうかいじょうけんは、ことなる集合しゅうごうかれる。

3みっつの問題形式もんだいけいしき

初期値問題しょきちもんだいInitial value problem は、初期面しょきめんでデータをあたえ、その時間発展じかんはってんもとめる問題もんだいである。Cauchy 問題もんだいはより一般いっぱんには初期超曲面しょきちょうきょくめんにデータを問題もんだいであり、このトラックでは全空間ぜんくうかん Rnt=0 にデータを標準例ひょうじゅんれいあつかう。

境界値問題きょうかいちもんだいBoundary value problem は、境界きょうかい領域りょういき境界きょうかいデータをあたえてかいもとめる問題もんだいである。このページでは、時間じかんふくまない Laplace 方程式ほうていしきや Poisson 方程式ほうていしき代表例だいひょうれいとしてあつかう。

初期境界値問題しょききょうかいちもんだいInitial-boundary value problem は、境界きょうかい空間領域くうかんりょういき時間発展じかんはってんあつかい、初期面しょきめん側面境界そくめんきょうかい両方りょうほうにデータをあたえる問題もんだいである。このページではおも有界領域ゆうかいりょういきかんがえ、有限区間上ゆうげんくかんじょう熱方程式ねつほうていしき波動方程式はどうほうていしき標準例ひょうじゅんれいとする。

形式けいしき方程式ほうていしき場所ばしょデータをあたえる場所ばしょ代表例だいひょうれい
初期値問題しょきちもんだいRn×(0,T)t=0全空間ぜんくうかんの transport・heat・wave
境界値問題きょうかいちもんだいΩΩLaplace・Poisson
初期境界値問題しょききょうかいちもんだいΩ×(0,T)t=0Ω×(0,T)有界領域ゆうかいりょういきの heat・wave

4境界条件きょうかいじょうけんみっつのかた

nΩ じょう外向そとむ単位法線たんいほうせんとし、法線微分ほうせんびぶん

nu=u·n

さだめる。代表的だいひょうてき境界条件きょうかいじょうけんつぎみっつである。

  • Dirichlet 条件じょうけん: u=g境界きょうかい未知関数みちかんすうあたい指定していする。
  • Neumann 条件じょうけん: nu=h境界きょうかい外向そとむ法線方向ほうせんほうこう微分びぶん指定していする。
  • Robin 条件じょうけん: αu+βnu=rあたい法線微分ほうせんびぶん線形結合せんけいけつごう指定していする。α,β同時どうじに 0 ではないとする。

Robin 条件じょうけんくだけで適切性てきせつせい自動的じどうてきしたがうわけではない。一意性いちいせい安定性あんていせいには、係数けいすう符号ふごう境界きょうかい分割ぶんかつ追加仮定ついかかてい必要ひつようである。

境界全体きょうかいぜんたいおなかた使つか必要ひつようはない。ΓDΓNかさならない境界部分きょうかいぶぶんとし、Ω=ΓDΓNけ、ΓD で Dirichlet 条件じょうけんΓN で Neumann 条件じょうけんすものを混合境界条件こんごうきょうかいじょうけんという。これはひとつの場所ばしょあたい法線微分ほうせんびぶん結合けつごうする Robin 条件じょうけんとはことなる。

Neumann 条件じょうけん流束りゅうそくとしてむときは符号ふごう注意ちゅういする。k>0熱伝導率ねつでんどうりつとし、熱流束ねつりゅうそくq=-kuさだめれば、外向そとむ流束りゅうそく

q·n=-knu

である。したがって「nu>0」と「外向そとむきの熱流束ねつりゅうそくせい」はおな意味いみではない。どちらをデータとしているかをしきから確認かくにんする必要ひつようがある。

5時間階数じかんかいすう初期条件しょきじょうけん

標準的ひょうじゅんてき発展方程式はってんほうていしきでは、時間じかんについての階数かいすう必要ひつよう初期しょきデータのかずしめす。

熱方程式ねつほうていしき

ut-κΔu=0

時間じかんについて一階いっかいなので、初期分布しょきぶんぷ u(x,0)=f(x)あたえる。波動方程式はどうほうていしき

utt-c2Δu=0

時間じかんについて二階にかいなので、初期変位しょきへんい初期速度しょきそくど

u(x,0)=f(x),ut(x,0)=v(x)

あたえる。ただし「階数かいすう同数どうすうけばかならずよい」という規則きそくではない。方程式ほうていしき拘束条件こうそくじょうけんがある場合ばあいや、データをめん方程式ほうていしき特性的とくせいてきめんである場合ばあいは、独立どくりつ指定していできるデータがわる。一階いっかい PDE でデータをける曲線きょくせん条件じょうけんは、つぎ特性曲線法とくせいきょくせんほうあつかう。

6transport 方程式ほうていしき: 条件じょうけん境界きょうかいえら

c>0 として、有界区間ゆうかいくかん 0<x<L の transport 方程式ほうていしき

ut+cux=0

かんがえる。情報じょうほう速度そくど cみぎすすむので、x=0情報じょうほう領域りょういきはい流入境界りゅうにゅうきょうかいx=L情報じょうほう流出境界りゅうしゅつきょうかいである。初期値しょきち u(x,0)=f(x)くわえて、流入境界りゅうにゅうきょうかいあたい u(0,t)=a(t)指定していするのが基本きほんである。流出境界りゅうしゅつきょうかい u(L,t)内部ないぶからはこばれてまるため、そこにも任意にんいあたい指定していすると一般いっぱんには過剰条件かじょうじょうけんになる。

かど (0,0) では f(0)=a(0)適合条件てきごうじょうけんになる。また、c<0 なら流入境界りゅうにゅうきょうかいx=Lわる。このれいは、境界条件きょうかいじょうけん場所ばしょ個数こすう幾何学的きかがくてき端点たんてんかずではなく、特性とくせいがどちらから情報じょうほうはこぶかでまることをしめす。

7適切性てきせつせい

適切てきせつ問題もんだいWell-posed problem とは、指定していしたかいのクラスでつぎ三条件さんじょうけんたす問題もんだいである。

  1. 存在そんざい: 条件じょうけんたすかいがある。
  2. 一意性いちいせい: 条件じょうけんたすかいひとつにさだまる。
  3. 安定性あんていせい: データをすこえたとき、かいえらんだ尺度しゃくどすこしだけわる。

三条件目さんじょうけんめは Hadamard の意味いみでの連続依存性れんぞくいぞんせいである。「公式こうしきける」ことだけでは適切性てきせつせい証明しょうめいにならない。どの関数空間かんすうくうかんでデータとかいおおきさをはかるかも問題設定もんだいせってい一部いちぶである。このページでは古典解こてんかい念頭ねんとう条件じょうけん配置はいち整理せいりし、energy estimate による一意性いちいせい連続依存性れんぞくいぞんせいのちの energy method であつかう。

たとえば熱方程式ねつほうていしき時間じかん逆向ぎゃくむきにいて過去かこ温度おんど復元ふくげんする問題もんだいでは、空間的くうかんてきこまかく振動しんどうする高周波成分こうしゅうはせいぶんちいさな誤差ごさ指数的しすうてき増幅ぞうふくされる。形式的けいしきてきかい逆算ぎゃくさんできても、データへの連続依存性れんぞくいぞんせいうしなえば適切てきせつ問題もんだいではない。存在そんざい一意性いちいせい安定性あんていせいけて確認かくにんする理由りゆうがここにある。

data/lecture/math/partial-differential-equations/introduction-to-energy-methods.lecture.n.md

8初期面しょきめん境界きょうかい適合条件てきごうじょうけん

初期面しょきめん側面境界そくめんきょうかいΩ×{0}まじわる。このかどふたつのデータが矛盾むじゅんしないための条件じょうけん適合条件てきごうじょうけんという。

たとえば熱方程式ねつほうていしき

u(x,0)=f(x),u(x,t)=g(x,t)(xΩ)

という Dirichlet 条件じょうけんし、閉領域へいりょういきまで連続れんぞく古典解こてんかいもとめるなら、すくなくとも

f(x)=g(x,0)(xΩ)

必要ひつようである。Neumann 条件じょうけん nu=h なら、十分じゅうぶんなめらかな初期値しょきちたいして nf=h(·,0)必要ひつようになる。さらにたか正則性せいそくせい要求ようきゅうすれば、方程式ほうていしき使つかってられる高階こうかい適合条件てきごうじょうけん必要ひつようになる。

固定端こていたん一次元波動方程式いちじげんはどうほうていしきu(0,t)=u(L,t)=0場合ばあいなめらかな古典解こてんかいには

f(0)=f(L)=0

必要ひつようであり、境界条件きょうかいじょうけん時間微分じかんびぶんできる正則性せいそくせいまで要求ようきゅうするなら

v(0)=v(L)=0

必要ひつようである。適合条件てきごうじょうけんやぶれていても、よわ意味いみt>0かぎればかいあつかえる場合ばあいがあるため、「かいまったくない」と一括ひとくくりにはしない。どの解概念かいがいねん境界きょうかいでの連続性れんぞくせい要求ようきゅうするかをさき固定こていする。

9熱方程式ねつほうていしき: 境界きょうかい総量そうりょうえる

ut=κΔuなめらかな有界領域ゆうかいりょういきかんがえる。発散定理はっさんていりにより

ddtΩudx=κΩnudS

である。したがって斉次せいじ Neumann 条件じょうけん nu=0 では総量そうりょう Ωudx保存ほぞんされる。一方いっぽう、Dirichlet 条件じょうけんでは境界きょうかいつうじた流出入りゅうしゅつにゅう一般いっぱんには 0 でなく、総量そうりょう保存ほぞんされない。このちがいは、境界条件きょうかいじょうけんたんなる計算上けいさんじょう付属物ふぞくぶつではなく、モデルの外部がいぶとの接続せつぞくあらわすことをしめす。

10Laplace・Poisson 方程式ほうていしき: じゅん Neumann 問題もんだい例外れいがい

Dirichlet 問題もんだいでは、適切てきせつ仮定かていのもとで境界値きょうかいち調和関数ちょうわかんすう一意いちいさだめる。これにたいじゅん Neumann 問題もんだいでは、uかいなら u+Cおな法線微分ほうせんびぶんつため、かい定数分ていすうぶんだけしかさだまらない。

さらに

-Δu=FinΩ,nu=honΩ

積分せきぶんすると、発散定理はっさんていりから必要条件ひつようじょうけん

ΩFdx=-ΩhdS

る。Laplace 方程式ほうていしき F=0 なら、ΩhdS=0必要ひつようである。この条件じょうけんたさないデータにはかいがない。たす場合ばあい一意性いちいせいるには、たとえば Ωudx=0 のような正規化せいきか追加ついかする。

11条件じょうけん過不足かぶそくなく

  • 方程式ほうていしき空間領域くうかんりょういき時間区間じかんくかん解概念かいがいねん最初さいしょ固定こていする。
  • 初期値問題しょきちもんだい境界値問題きょうかいちもんだい初期境界値問題しょききょうかいちもんだいのどれかを区別くべつする。
  • 境界きょうかい外向そとむ法線ほうせん流束りゅうそく符号規約ふごうきやく確認かくにんする。
  • 時間階数じかんかいすうだけで機械的きかいてき条件数じょうけんすうめず、特性とくせい拘束条件こうそくじょうけん確認かくにんする。
  • 初期面しょきめん側面境界そくめんきょうかいまじわりで適合条件てきごうじょうけん確認かくにんする。
  • 存在そんざい一意性いちいせい安定性あんていせい別々べつべつ主張しゅちょうとして確認かくにんする。
  • じゅん Neumann 問題もんだいでは積分条件せきぶんじょうけん定数ていすう不定性ふていせい確認かくにんする。

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

このページはなめらかな領域りょういき標準的ひょうじゅんてき線形せんけい PDE を中心ちゅうしんに、条件じょうけん配置はいち説明せつめいした。境界きょうかいなめらかでない場合ばあい係数けいすう不連続ふれんぞく場合ばあい非線形方程式ひせんけいほうていしき衝撃波しょうげきはしょうじる場合ばあいには、古典解こてんかいではなく弱解じゃくかいやトレースの理論りろん必要ひつようになる。また、適切性てきせつせい証明しょうめいには maximum principle、energy estimate、固有関数展開こゆうかんすうてんかいなど、方程式ほうていしきかたおうじた道具どうぐ必要ひつようである。

13つぎむページ

一階いっかい PDE では、初期しょきデータをどの曲線きょくせんけるかを特性曲線法とくせいきょくせんほう判定はんていする。二階にかい PDE では主部しゅぶ分類ぶんるいし、適切てきせつ追加仮定ついかかていのもとで自然しぜん問題設定もんだいせってい情報伝播じょうほうでんぱ考察こうさつする。その、heat・wave・Laplace 方程式ほうていしき条件じょうけんかい性質せいしつをどうえるかを比較ひかくする。

data/lecture/math/partial-differential-equations/method-of-characteristics.lecture.n.md data/lecture/math/partial-differential-equations/classification-of-second-order-linear-pdes.lecture.n.md data/lecture/math/partial-differential-equations/heat-wave-and-laplace-equations.lecture.n.md data/lecture/math/vector-calculus/green-gauss-and-stokes-theorems.lecture.n.md

Initial-Value and Boundary-Value Problems for PDEs

1Introduction

A PDE problem is not determined by the equation alone. It is specified only after combining the domain of the unknown with the locations and kinds of prescribed data.

This lecture distinguishes initial-value, boundary-value, and initial-boundary-value problems and explains how to assess existence, uniqueness, stability, and compatibility rather than merely counting conditions.

2Sets on which the problem is posed

Let ΩRn be a connected spatial domain, let Ω be its boundary, and let T>0. Throughout this lecture, assume that the boundary is sufficiently smooth for an outward unit normal and the divergence theorem to be available. For time evolution, consider u(x,t) on the spacetime domain

QT=Ω×(0,T).

The section Ω×{0} is the initial surface, while Ω×(0,T) is the lateral boundary. Initial and boundary conditions are prescribed on different sets.

3Three problem formulations

An {initial-value problem} prescribes data on an initial surface and asks for the subsequent evolution. More generally, a Cauchy problem prescribes data on an initial hypersurface; this track uses the standard example of data on t=0 in the whole space Rn.

A {boundary-value problem} prescribes data on the boundary of a domain. The principal static examples here are Laplace's and Poisson's equations.

An {initial-boundary-value problem} concerns time evolution in a spatial domain with boundary and prescribes data on both the initial surface and the lateral boundary. Standard examples are the heat and wave equations on a bounded interval.

FormulationWhere the equation is imposedWhere data are prescribedRepresentative examples
Initial-value problemRn×(0,T)t=0whole-space transport, heat, and wave equations
Boundary-value problemΩΩLaplace and Poisson equations
Initial-boundary-value problemΩ×(0,T)t=0 and Ω×(0,T)heat and wave equations on bounded domains

4Three principal boundary conditions

Let n be the outward unit normal on Ω and define

nu=u·n.

The principal boundary conditions are:

  • Dirichlet: u=g, prescribing the value of the unknown on the boundary;
  • Neumann: nu=h, prescribing its outward normal derivative;
  • Robin: αu+βnu=r, prescribing a linear combination of the value and normal derivative, with α,β not simultaneously zero.

Writing a Robin condition does not by itself imply well-posedness. Uniqueness and stability may require sign assumptions on the coefficients and conditions on the boundary decomposition.

The same condition need not be used on the entire boundary. If disjoint boundary portions satisfy Ω=ΓDΓN, imposing Dirichlet data on ΓD and Neumann data on ΓN gives a mixed boundary condition. This differs from a Robin condition, which combines a value and a normal derivative at the same location.

The sign convention matters when Neumann data are interpreted as flux. If k>0 is the thermal conductivity and q=-ku is the heat flux, then the outward flux is

q·n=-knu.

Thus nu>0 does not mean that the outward heat flux is positive. The equation must be inspected to determine which quantity is prescribed.

5Time order and initial conditions

For standard evolution equations, the order in time indicates the number of initial data. The heat equation

ut-κΔu=0

is first order in time and is supplied with an initial distribution u(x,0)=f(x). The wave equation

utt-c2Δu=0

is second order in time and is supplied with initial displacement and velocity,

u(x,0)=f(x),ut(x,0)=v(x).

This is not a universal rule that exactly as many arbitrary conditions as the time order may always be imposed. Constraints in the equation or a characteristic data surface can change which data are independent. The next lecture treats the admissible data curves for first-order PDEs.

6Transport equation: selecting the boundary carrying data

Let c>0 and consider on 0<x<L

ut+cux=0.

Information travels to the right with speed c. Hence x=0 is the inflow boundary and x=L is the outflow boundary. In addition to u(x,0)=f(x), one ordinarily prescribes the inflow value u(0,t)=a(t). The value u(L,t) is transported from the interior, so prescribing an arbitrary value there as well generally overdetermines the problem.

At the corner (0,0), compatibility requires f(0)=a(0). If c<0, the inflow boundary changes to x=L. Thus the location and number of boundary conditions are determined by the direction in which characteristics carry information, not merely by the number of geometric endpoints.

7Well-posedness

A {well-posed problem} satisfies the following three conditions in the specified solution class:

  1. existence: at least one solution satisfies the conditions;
  2. uniqueness: at most one solution satisfies them;
  3. stability: small changes in the data produce small changes in the solution in the selected norms.

The third condition is continuous dependence in Hadamard's sense. The existence of an explicit formula alone does not prove well-posedness. The function spaces and norms used to measure data and solutions are part of the problem formulation. This lecture organizes the placement of conditions with classical solutions in mind; a later lecture uses energy estimates for uniqueness and continuous dependence.

In the backward heat problem, small errors in rapidly oscillating spatial components are amplified exponentially when reconstructing an earlier temperature. Even if the equation can be inverted formally, loss of continuous dependence makes the problem ill-posed. This illustrates why existence, uniqueness, and stability must be assessed separately.

data/lecture/math/partial-differential-equations/introduction-to-energy-methods.lecture.n.md

8Compatibility of the initial surface and boundary

The initial surface and lateral boundary meet on Ω×{0}. Conditions preventing the two data from contradicting each other there are called compatibility conditions.

For the heat equation with

u(x,0)=f(x),u(x,t)=g(x,t)(xΩ),

a classical solution continuous up to the closed domain must satisfy at least

f(x)=g(x,0)(xΩ).

For Neumann data nu=h and a sufficiently smooth initial value, one needs nf=h(·,0). Higher regularity requires higher-order compatibility conditions obtained using the equation.

For the one-dimensional wave equation with fixed endpoints u(0,t)=u(L,t)=0, a smooth classical solution requires

f(0)=f(L)=0.

If the boundary condition may be differentiated in time at the required regularity, then one also needs

v(0)=v(L)=0.

Failure of compatibility does not always imply that no solution exists: weak solutions or solutions restricted to t>0 may still be meaningful. The solution concept and required boundary continuity must be fixed first.

9Heat equation: the boundary changes the total quantity

Let ut=κΔu on a smooth bounded domain. The divergence theorem gives

ddtΩudx=κΩnudS.

Under homogeneous Neumann data nu=0, the total quantity Ωudx is conserved. Under Dirichlet data, the boundary flux is generally nonzero, so the total quantity need not be conserved. Boundary conditions therefore describe interaction with the exterior rather than being merely computational accessories.

10Laplace and Poisson equations: the pure Neumann exception

Under suitable assumptions, boundary values in a Dirichlet problem determine a harmonic function uniquely. In a pure Neumann problem, however, if u is a solution, then u+C has the same normal derivative, so the solution is determined only up to a constant.

Integrating

-Δu=FinΩ,nu=honΩ

and applying the divergence theorem gives the necessary condition

ΩFdx=-ΩhdS.

For Laplace's equation, F=0, so ΩhdS=0 is necessary. If it fails, no solution exists. Even when it holds, uniqueness requires a normalization such as Ωudx=0.

11Reading conditions without excess or deficiency

  • Fix the equation, spatial domain, time interval, and solution concept first.
  • Distinguish initial-value, boundary-value, and initial-boundary-value problems.
  • Check the outward-normal and flux sign conventions.
  • Do not determine the number of conditions mechanically from the time order; inspect characteristics and constraints.
  • Check compatibility where the initial surface meets the lateral boundary.
  • Treat existence, uniqueness, and stability as separate claims.
  • For a pure Neumann problem, check the integral condition and indeterminacy by constants.

12Scope of validity

This lecture focused on smooth domains and standard linear PDEs. Nonsmooth boundaries, discontinuous coefficients, and shocks in nonlinear equations require weak solutions and trace theory rather than only classical solutions. Proving well-posedness also requires tools appropriate to the equation's type, such as maximum principles, energy estimates, and eigenfunction expansions.

13Subsequent lectures

For first-order PDEs, the method of characteristics determines which curves can carry initial data. For second-order PDEs, classify the principal part and then examine natural problem formulations and information propagation under appropriate additional assumptions. The model comparison then examines how conditions affect heat, wave, and Laplace equations.

data/lecture/math/partial-differential-equations/method-of-characteristics.lecture.n.md data/lecture/math/partial-differential-equations/classification-of-second-order-linear-pdes.lecture.n.md data/lecture/math/partial-differential-equations/heat-wave-and-laplace-equations.lecture.n.md data/lecture/math/vector-calculus/green-gauss-and-stokes-theorems.lecture.n.md
raw .n.md をコピー
loc をコピー (filepath:line ~ line)
copy share link
copy encoded share link
path をコピー
copy share link
copy encoded share link
copy share link
copy encoded share link
タブを全て閉じる