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)開集合かいしゅうごう DR3a,b,FC1(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 は各初期点かくしょきてんちかくで局所的きょくしょてき一意いちいかいつ。F0 ならすべての特性上とくせいじょうz、すなわち u保存ほぞんされる。一般いっぱんには z も ODE にしたがって変化へんかする。

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

区間くかん I じょうx0,t0,hC1(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 である。したがって、hC1(I)げた曲線きょくせんD ないにあるかぎり、標準的ひょうじゅんてきt=0初期線しょきせん非特性的ひとくせいてきであり、局所古典解きょくしょこてんかい構成こうせいするためのデータをける。

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

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

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

gC1(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定数ていすうgC1(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変係数へんけいすうれい

gC1(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 方程式ほうていしき特性交差とくせいこうさ

非線形ひせんけいれいとして gC1(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
タブを全て閉じる