tecktom.com

特異摂動の共通構造

2026年10月9日7分で読めますAI生成教科書『虎の巻』第7章 スケール分離と特異摂動：Born–Oppenheimer から境界層まで

AI 生成・木下が未勉強

近似の数理　第7章 §7.9　★★★ 境界層・WKB・断熱定理・BO を一つの表で並べ、物理の直観の位置と破綻の場所を特定する

この章で扱った近似はすべて、次の3段の構造をもつ。

  1. 小パラメータが方程式の最高階の項（あるいは最も速い過程）に掛かっていて、ε=0\varepsilon =0 で方程式の型が変わる。
  2. その結果、解には二つの尺度が現れる。ε=0\varepsilon =0 の方程式で記述できる「外部」（遅い・広い）と、伸長変数で記述する「内部」（速い・薄い）である。
  3. 外部と内部はそれぞれの方程式で解き、両者が重なる中間領域で接合して一つの近似にまとめる。
問題小パラメータ外部（遅い・広い）内部（速い・薄い）破綻する場所
モデル問題ε\varepsilon（y′′y'' の係数）1階方程式 y′+y=0y'+y=0幅 ε\varepsilon の層、Y′′+Y′=0Y''+Y'=0y′y' の係数が符号を変える点
多重尺度・平均化ε\varepsilon（減衰・非線形）振幅方程式（遅い時間 εt\varepsilon t）速い振動 e±it\mathrm{e}^{\pm\mathrm{i} t}時間 ≫1/ε\gg1/\varepsilon、共鳴
WKBℏ2\hbar^2（ψ′′\psi'' の係数）古典軌道と位相 ∫p dx/ℏ\int p\,\mathrm{d} x/\hbar転回点の層（幅 ℏ2/3\hbar^{2/3}）、Airy転回点・低い量子数
断熱定理ε\varepsilon（∂t\partial_t の係数）瞬間ごとの固有状態動的位相の速い回転ギャップ →0\to0（Landau–Zener）
Born–Oppenheimerκ4=me/M\kappa^4=m_e/M（ΔR\Delta_{\bm{R}} の係数）核の曲面上の運動電子の固有値問題円錐交差（g≲κg\lesssim\kappa）
Prandtl 境界層Re−1\mathrm{Re}^{-1}（Δu\Delta\bm{u} の係数）Euler 流幅 Re−1/2\mathrm{Re}^{-1/2} の層剥離・非単調な速度分布
Stokes／OseenRe\mathrm{Re}Oseen 流（r≳a/Rer\gtrsim a/\mathrm{Re}）Stokes 流（物体近傍）2次元で ln⁡\ln 型の展開
潤滑ε=H/L\varepsilon =H/L圧力の Reynolds 方程式隙間方向の放物線分布隙間の端・急な形状変化

表の「外部」と「内部」のどちらが物理屋の言う「遅いもの」「速いもの」に当たるかは問題によって入れ替わるが、共通しているのは、速い（薄い）方を先に解いて、その結果を遅い（広い）方の方程式の係数やポテンシャルとして使うという手続きである。BO では電子の固有値が核のポテンシャルに、潤滑では隙間方向の流量が圧力の方程式に、平均化では速い位相の平均が振幅の方程式に、断熱定理では瞬間の固有値が動的位相になる。物理の直観「速いものと遅いものを分ける」は、数学的にはこの手続きと、それが接合条件で矛盾なく閉じることに対応する。

直観が明示しないのは、分離が成り立つための一様性の条件である。どの例でも、外部と内部の尺度の比が全域で大きいことが必要で、それが局所的に破れる場所で近似は破綻する。断熱定理と BO ではギャップが閉じる点、WKB では運動量が 0 になる転回点、境界層では壁の摩擦が 0 になる剥離点がそれに当たる。これらはいずれも「二つの尺度が同じ大きさになる」場所であり、そこでは新たな内部問題（Airy 関数、Landau–Zener 模型、剥離の局所解析）を立てて解き直すか、もはや二尺度の分離を使えないかのどちらかになる。さらに、破綻の影響が冪でなく exp⁡(−c/ε)\exp(-c/\varepsilon ) の形で現れる場合（Landau–Zener 遷移、トンネル現象、Stirling 級数の最適打ち切り誤差）は、漸近展開のどの次数にも現れないので、展開を高次まで進めても検出できない。

例題 7.9　両端に層ができる反応拡散方程式 ★★

オリジナル0<ε≪10<\varepsilon \ll1 として、境界値問題

εy′′−y=−1,y(0)=y(1)=0\varepsilon y'' - y = -1,\qquad y(0)=y(1)=0

を考える。

  • (1) 上の表の1行（小パラメータ・外部・内部・破綻する場所）をこの問題について埋めよ。境界層の厚さは主要釣り合いで決めよ。
  • (2) 外部解・内部解を求め、合成解を作れ。
  • (3) 厳密解と比べて合成解の誤差を求め、モデル問題の誤差 O(ε)\mathcal{O}(\varepsilon ) と違う理由を述べよ。
解答

(1) ε=0\varepsilon =0 と置くと y=1y=1 という代数方程式になり、境界条件を一つも課せない。端の近くで X=x/δX=x/\delta と伸長すると(ε/δ2)Y′′−Y=−1(\varepsilon /\delta^2)Y''-Y=-1 で、Y′′Y'' の項が他と釣り合うのは δ=ε\delta=\sqrt\varepsilon のときである。1階微分の項が無いので移流の向きが無く、層は両端にできる。表の行は、小パラメータが ε\varepsilon（y′′y'' の係数）、外部が y=1y=1、内部が幅 ε\sqrt\varepsilon の層での Y′′−Y=−1Y''-Y=-1、破綻する場所は無し、となる。yy の係数が区間の内部で 0 になる場合（εy′′=q(x)y−1\varepsilon y''=q(x)y-1 で qq が符号を変える場合）は、WKB の転回点と同じ型の内部層が現れて別扱いになる。

(2) 外部解は y0=1y_0=1 である。左端では Y′′−Y=−1Y''-Y=-1、Y(0)=0Y(0)=0 で有界な解は Y=1−e−XY=1-\mathrm{e}^{-X}（X=x/εX=x/\sqrt\varepsilon）で、X→∞X\to\infty で外部解の値 1 に自動的に接合する。右端も同様に ξ=(1−x)/ε\xi=(1-x)/\sqrt\varepsilon で 1−e−ξ1-\mathrm{e}^{-\xi} である。外部解と二つの内部解を足し、二重に数えた共通部分 1 を2回引いて

yc(x)=1−e−x/ε−e−(1−x)/ε.y_c(x) = 1 - \mathrm{e}^{-x/\sqrt\varepsilon } - \mathrm{e}^{-(1-x)/\sqrt\varepsilon }.

(3) 厳密解は y=1−cosh⁡((x−12)/ε)/cosh⁡(1/(2ε))y = 1-\cosh\bigl((x-\tfrac12)/\sqrt\varepsilon \bigr)/\cosh\bigl(1/(2\sqrt\varepsilon )\bigr) である。分子と分母にe−1/(2ε)\mathrm{e}^{-1/(2\sqrt\varepsilon )} を掛けると、r=e−1/εr=\mathrm{e}^{-1/\sqrt\varepsilon } として

y=1−e−x/ε+e−(1−x)/ε1+r,y−yc=(e−x/ε+e−(1−x)/ε)r1+ry = 1-\frac{\mathrm{e}^{-x/\sqrt\varepsilon }+\mathrm{e}^{-(1-x)/\sqrt\varepsilon }}{1+r},\qquad y-y_c = \bigl(\mathrm{e}^{-x/\sqrt\varepsilon }+\mathrm{e}^{-(1-x)/\sqrt\varepsilon }\bigr)\frac{r}{1+r}

となる。括弧の最大は端での 1+r1+r なので sup⁡∣y−yc∣=r=e−1/ε\sup\lvert y-y_c\rvert = r = \mathrm{e}^{-1/\sqrt\varepsilon } であり、ε=0.01\varepsilon =0.01 で 4.5×10−54.5\times10^{-5} である（検算スクリプト）。誤差は ε\varepsilon のどの冪よりも小さい。

モデル問題では外部解 e1−x\mathrm{e}^{1-x} に ε\varepsilon の補正（特性根 λ1=−1−ε+⋯\lambda_1=-1-\varepsilon +\cdots のずれ）があり、それが O(ε)\mathcal{O}(\varepsilon ) の誤差になった。この問題では外部解 y=1y=1 が外部方程式をすべての次数で厳密に満たし、内部方程式も省いた項が無いので厳密である。残る誤差は、片方の層の裾がもう片方の端の境界条件をわずかに乱すことだけで、その大きさは層の裾の値 e−1/ε\mathrm{e}^{-1/\sqrt\varepsilon } になる。表の「破綻する場所」が無い問題でも、誤差が冪で出るか指数で出るかは外部解に補正が要るかどうかで決まる。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI（Claude）が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先：support@tecktom.com

PDF第7章の PDF近似の数理の PDF（全章）

この記事が含まれている連載

連載AI生成教科書『虎の巻』授業とは別に、分野をまたいで気になったことを AI（Claude）が教材にし、木下が勉強しながら公開しています。上の段は AI が書いた本文、下の段は木下が勉強したあとに書いた文です。1787本の記事・木下 智貴
← AI生成教科書『虎の巻』の記事一覧へ