markdown
特性曲線法md e4d4769
lecture/math/partial-differential-equations/method-of-characteristics.lecture.n.md
Download PDF

特性曲線法とくせいきょくせんほう

date2026-07-16document_iddoc_254189ad3c7267d2aa811d17bba947b0description一階 PDE の特性 ODE、非特性的な初期曲線、局所解の定義域、特性交差による古典解の限界を整理する。prerequisitesPDE とは何か / PDE の初期値問題と境界値問題 / 接平面・連鎖律・Jacobian / 一階微分方程式の解法診断type講義content_typelecturestatusactiverelateddata/lecture/math/partial-differential-equations/introduction-to-pdes.lecture.n.md / data/lecture/math/partial-differential-equations/initial-and-boundary-value-problems.lecture.n.md / data/lecture/math/partial-differential-equations/classification-of-second-order-linear-pdes.lecture.n.md / data/lecture/math/partial-differential-equations/transport-equations-and-conservation-laws.lecture.n.md / data/lecture/math/multivariable-calculus/tangent-planes-chain-rule-and-jacobian.lecture.n.md / data/lecture/math/differential-equations/first-order-ode-solution-diagnostics.lecture.n.md
mathpartial-differential-equationscharacteristicslecture

1導入どうにゅう

この講義こうぎでは、一階いっかい PDE を特性曲線とくせいきょくせんに沿そって ODE へ還元かんげんし、その ODE から局所的きょくしょてきな古典解こてんかいを再構成さいこうせいする方法ほうほうを説明せつめいする。ただし、データを置おく曲線きょくせんが特性方向とくせいほうこうと平行へいこうな場合ばあいや、異ことなる特性とくせいが交差こうさする場合ばあいには、存在そんざい・一意性いちいせい・滑なめらかさが失うしなわれうる。

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

特性曲線とくせいきょくせんCharacteristic curve は、PDE の微分方向びぶんほうこうに沿そって進すすむ曲線きょくせんであり、その上うえで偏微分方程式へんびぶんほうていしきが ODE に還元かんげんされる。

3一般形いっぱんけいと特性とくせい ODE

(x,t,z) の開集合かいしゅうごう D⊂R3 で a,b,F∈C1(D) とする。準線形じゅんせんけいQuasilinearとは、最高階さいこうかいの微分びぶんが線形せんけいに現あらわれ、その係数けいすうが独立変数どくりつへんすうと未知関数みちかんすうに依存いぞんしてよい形かたちである。このページでは準線形じゅんせんけいの一階いっかい PDE

a(x,t,u)ux+b(x,t,u)ut=F(x,t,u)

に対たいし、z(s)=u(x(s),t(s)) と置おく。連鎖律れんさりつから

dzds=uxdxds+utdtds

なので、特性方程式とくせいほうていしき

dxds=a(x,t,z),dtds=b(x,t,z),dzds=F(x,t,z)

を得える。C1 係数けいすうなら、この ODE は各初期点かくしょきてんの近ちかくで局所的きょくしょてきに一意いちいな解かいを持もつ。F≡0 ならすべての特性上とくせいじょうで z、すなわち u が保存ほぞんされる。一般いっぱんには z も ODE に従したがって変化へんかする。

4データを置おく曲線きょくせん

区間くかん I 上じょうで x0,t0,h∈C1(I) とし、(x0(r),t0(r),h(r))∈D を満みたす初期曲線しょききょくせんとデータを

γ(r)=(x0(r),t0(r)),u(x0(r),t0(r))=h(r)

とする。曲線上きょくせんじょうの特性方向とくせいほうこうは (a,b) である。初期曲線しょききょくせんと特性方向とくせいほうこうが横断おうだんする、すなわち

x0′(r)b(x0,t0,h)-t0′(r)a(x0,t0,h)≠0

であるとき、この初期曲線しょききょくせんを非特性的ひとくせいてきNon-characteristicという。この行列式ぎょうれつしきは、(r,s) から (x,t) への特性座標写像とくせいざひょうしゃぞうが初期曲線しょききょくせんの近ちかくで局所的きょくしょてきに逆写像ぎゃくしゃぞうを持もつための条件じょうけんである。したがって特性とくせい ODE で得えた z(r,s) を u(x,t) として戻もどせる。

実際じっさい、逆写像ぎゃくしゃぞうを (r,s)=(r(x,t),s(x,t)) とし、u(x,t)=z(r(x,t),s(x,t)) と定義ていぎする。特性とくせい ODE と連鎖律れんさりつから

zs=uxxs+utts=a(x,t,u)ux+b(x,t,u)ut=F(x,t,u)

が成立せいりつする。また s=0 では (x,t,z)=(x0(r),t0(r),h(r)) なので、u(x0(r),t0(r))=h(r) も成立せいりつする。したがって、ODE から再構成さいこうせいした u は PDE と初期条件しょきじょうけんの双方そうほうを満みたす。

正規形せいきけい

ut+A(x,t,u)ux=C(x,t,u)

で x0(r)=r,t0(r)=0 とすれば、上うえの行列式ぎょうれつしきは 1 である。したがって、h∈C1(I) で持もち上あげた曲線きょくせんが D 内ないにある限かぎり、標準的ひょうじゅんてきな t=0 の初期線しょきせんは非特性的ひとくせいてきであり、局所古典解きょくしょこてんかいを構成こうせいするためのデータを置おける。

この構成こうせいが有効ゆうこうなのは、持もち上あげた特性とくせい (x,t,z) が D の中なかにあり、ODE 解かいが存在そんざいし、特性座標写像とくせいざひょうしゃぞうの Jacobian が 0 でない範囲はんいである。解かいの定義域ていぎいきはこの写像しゃぞうの像ぞうであり、自動的じどうてきに全空間ぜんくうかん・全時刻ぜんじこくへ広ひろがるわけではない。大域解たいいきかいには、ODE の大域存在たいいきそんざいと写像しゃぞうの大域的たいいきてきな一対一性いったいいちせい・到達性とうたつせいを別べつに示しめす必要ひつようがある。

行列式ぎょうれつしきが 0 の特性的とくせいてきな曲線きょくせんに任意にんいのデータを置おくと、同おなじ特性上とくせいじょうでデータが矛盾むじゅんして解かいが存在そんざいしない場合ばあいや、特性とくせいの外側そとがわへ値あたいを決きめられず一意性いちいせいを失うしなう場合ばあいがある。条件じょうけんは曲線きょくせんの本数ほんすうだけでなく、PDE の微分方向びぶんほうこうとの位置関係いちかんけいで判断はんだんする。

5定係数ていけいすう transport の具体例ぐたいれい

g∈C1(R) とし、全空間ぜんくうかんの Cauchy 問題もんだい

ut+2ux=0,u(x,0)=g(x)

を考かんがえる。特性曲線とくせいきょくせんは x-2t=ξ であり、その上うえでは du/dt=ut+2ux=0 である。したがって

u(x,t)=g(x-2t)

となる。初期形状しょきけいじょうが速度そくど 2 で右方向みぎほうこうへ移動いどうする。

6生成せいせい・消滅項しょうめつこうがある場合ばあい

c,λ∈R を定数ていすう、g∈C1(R) とし、

ut+cux=λu,u(x,0)=g(x)

を考かんがえる。特性とくせいは x=ξ+ct であり、その上うえでは

ddtu(ξ+ct,t)=λu(ξ+ct,t)

となる。したがって

u(x,t)=eλtg(x-ct)

である。形かたちは速度そくど c で移動いどうするが、値あたいは λ に応おうじて増幅ぞうふくまたは減衰げんすいする。「特性とくせいに沿そって一定いってい」は右辺うへんが 0 の場合ばあいだけであり、特性曲線法とくせいきょくせんほうそのものは値あたいが変かわる ODE にも適用てきようできる。

7変係数へんけいすうの例れい

g∈C1(R) とし、ut+xux=0 を考かんがえる。特性曲線とくせいきょくせんは x′=x,t′=1 を満みたす。したがって x(t)=Cet であり、xe-t が保存ほぞんされる。初期条件しょきじょうけん u(x,0)=g(x) なら、u(x,t)=g(xe-t) となる。

ここでは初期点しょきてん ξ と時刻じこく t から x=ξet を得える。逆写像ぎゃくしゃぞう ξ=xe-t が存在そんざいするため、特性上とくせいじょうの値あたいを u(x,t) として戻もどせる。一般いっぱんの問題もんだいでも、特性とくせい ODE を解とくだけでなく、特性とくせいのラベルと時刻じこくから空間点くうかんてんへの写像しゃぞうが一対一いったいいちである範囲はんいを確認かくにんする必要ひつようがある。その範囲はんいが古典解こてんかいの定義域ていぎいきになる。

8Burgers 方程式ほうていしきと特性交差とくせいこうさ

非線形ひせんけいの例れいとして g∈C1(R) を仮定かていし、

ut+uux=0,u(x,0)=g(x)

を考かんがえる。初期点しょきてん ξ から出でる特性上とくせいじょうでは u=g(ξ) が一定いっていなので、

x=ξ+tg(ξ),u=g(ξ)

である。ξ から x への写像しゃぞうの微分びぶんは

∂x∂ξ=1+tg′(ξ)

となる。この量りょうが 0 でない間あいだは、ξ を (x,t) から局所的きょくしょてきに復元ふくげんでき、古典解こてんかいを構成こうせいできる。ある ξ で 1+tg′(ξ)=0 になると特性座標写像とくせいざひょうしゃぞうが局所可逆性きょくしょかぎゃくせいを失うしない、一価いっかで滑なめらかな古典解こてんかいをこの構成こうせいのまま続つづけられない。圧縮あっしゅくが進すすむ場合ばあいは、その後ごに異ことなる特性とくせいが同おなじ点てんへ到達とうたつして交差こうさする。

さらに陰関数いんかんすうを微分びぶんすると、古典解こてんかいの範囲はんいで

ux(x,t)=g′(ξ)1+tg′(ξ)

となる。したがって座標写像ざひょうしゃぞうの退化たいかと同時どうじに、分母ぶんぼが 0 へ近ちかづいて勾配こうばいが発散はっさんしうる。

特とくに

m=infξg′(ξ)

が -∞<m<0 なら、最初さいしょの座標退化ざひょうたいか・勾配破綻時刻こうばいはたんじこくの候補こうほは

T*=-1m

である。ただしこの式しきを厳密げんみつな破綻時刻はたんじこくとして使つかうには、g の正則性せいそくせいや下限かげんの達成たっせいなどを確認かくにんする必要ひつようがある。g′(ξ)[PARSE ERROR: Undefined("Command(\"ge\")")]0 が全すべての ξ で成立せいりつするなら、正せいの時刻じこくにはこの機構きこうによる特性交差とくせいこうさは起おこらない。交差後こうさごは保存則ほぞんそくの積分形せきぶんけいに基もとづく弱解じゃくかいと、解かいを選えらぶ entropy 条件じょうけんが必要ひつようになる。これらは後のちの transport と保存則ほぞんそくの講義こうぎで定義ていぎする。

9不適用ふてきようの例れい

二階にかい PDE では、一階いっかいの方向微分ほうこうびぶんを ODE へ変換へんかんする構造こうぞうが直接ちょくせつには成立せいりつしない。たとえば熱方程式ねつほうていしき ut=κuxx は拡散かくさんを含ふくむため、単一たんいつの曲線きょくせんに沿そって値あたいが保存ほぞんされる問題もんだいではない。

二階にかい PDE にも特性方向とくせいほうこうという概念がいねんはあるが、ここで扱あつかった三本さんぼんの特性とくせい ODE による直接ちょくせつの再構成さいこうせいとは別べつの理論りろんである。二階主部にかいしゅぶの特性とくせいと型かたは次つぎの分類講義ぶんるいこうぎで扱あつかう。

data/lecture/math/partial-differential-equations/classification-of-second-order-linear-pdes.lecture.n.md

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

特性曲線法とくせいきょくせんほうで古典解こてんかいを構成こうせいするには、特性とくせい ODE の局所存在きょくしょそんざいと一意性いちいせい、初期曲線しょききょくせんの非特性条件ひとくせいじょうけん、座標写像ざひょうしゃぞうの局所可逆性きょくしょかぎゃくせいを確認かくにんする必要ひつようがある。係数けいすうが滑なめらかでない場合ばあい、ODE 自体じたいの一意性いちいせいが失うしなわれる場合ばあい、特性とくせいが交差こうさする場合ばあいには、この手順てじゅんをそのまま続つづけられない。

11関連かんれんリンク

data/lecture/math/differential-equations/first-order-ode-solution-diagnostics.lecture.n.md data/lecture/math/partial-differential-equations/initial-and-boundary-value-problems.lecture.n.md data/lecture/math/partial-differential-equations/transport-equations-and-conservation-laws.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
タブを全て閉じる