tecktom.com

三つの型の中で C 型が残る理由

2026年10月9日6分で読めますAI生成教科書『虎の巻』第9章 均質化と「閉じない」近似

AI 生成・木下が未勉強

近似の数理 第9章 §9.8 ★★★ 近似が閉じるための構造の有無で、正当化できるものとできないものを分ける

この章で扱った近似を、第1章の3型の地図に置き直す。

表 9.2 この章の近似の型と、それを閉じさせている構造
近似型閉じさせる構造数学が保証するもの
均質化・有効係数B周期性またはエルゴード性と ε/L→0\varepsilon /L\to0極限方程式、収束の意味(強/弱)、速さ、上下界
Darcy 則B同上(細孔の周期構造)透過率の公式と正定値性
Lorenz の 3 モードA の形臨界点の近くでのモードの分離臨界点の近くの分岐。長時間・強非線型では保証なし
渦粘性・kk–ε\varepsilonCなし(スケール分離なし)実現可能性などの必要条件のみ
4/5 則厳密一様性・等方性・異常散逸の仮定仮定の下での厳密な関係
Kohn–Sham の枠組み厳密定義による書き換えExcE_{\mathrm{xc}} が厳密なら厳密
LDA・GGA・混成Cなし(LDA はゆっくり変わる密度の極限で厳密)厳密条件の充足

B 型の近似が閉じるのは、小さいパラメータ(ε/L\varepsilon /L、1/N1/N、Kn\mathrm{Kn}、質量比)があり、しかもその極限で消える自由度(微細構造の振動、粒子間相関、非平衡の分布)を既知の構造が決めているからである。均質化では周期性がセル問題を与え、平均場極限では対称性とカオスの伝播が相関を消し、流体力学極限では局所 Maxwell 分布が高次モーメントを決めた。極限で消える自由度の行き先が分かっているので、残った自由度だけで方程式が閉じる。

C 型が残るのは、(i) 小さいパラメータがない(乱流では渦の大きさと平均流の尺度が同程度、分子の電子では交換相関の穴の大きさと密度変化の尺度が同程度)か、(ii) 小さいパラメータはあっても極限で消える自由度の行き先を決める構造がない(Re→∞\mathrm{Re}\to\infty の乱流)、のどちらかの場合である。このとき「未知の項を既知の量で表す」関係はどこからも出てこないので、物理的直観・次元解析・実験への当てはめで作るしかない。

C 型の近似を評価する方法は、それでも空ではない。この章で見た方法をまとめると次の三つである。

  1. 厳密に解ける極限での一致:LDA の一様電子気体、渦粘性モデルの壁近くの対数領域。C 型の近似の多くは、ある B 型の極限で正しい式を外挿したものである。
  2. 厳密条件:スケーリング則、Lieb–Oxford の下界、Reynolds 応力の実現可能性、4/5 則。近似の正しさは保証しないが、満たさない近似は排除できる。
  3. ベンチマーク:信頼できる高精度の数値解(直接数値計算、高精度の量子化学計算)や実験値の集合との系統的な比較。

物理の近似の多くが「直観に基づく」と感じられるのは、C 型の近似が少なからず含まれているからである。B 型・A 型の近似についてはその直観を定理で裏付けられ、C 型の近似についてはその直観がどの極限で正しくどこから外挿なのかを特定できる。この区別が、微分方程式の理論が物理の近似に与える意味付けである。

例題 9.9 消す自由度の大きさで小さいパラメータを測る ★★

オリジナル直径 D=0.10 mD=0.10\,\mathrm{m} の円管を、空気(動粘性率 ν=1.5×10−5 m2/s\nu=1.5\times10^{-5}\,\mathrm{m^2/s}、平均自由行程 λ=6.8×10−8 m\lambda=6.8\times10^{-8}\,\mathrm{m})がReynolds 数 Re=UD/ν=1.0×105\mathrm{Re}=UD/\nu=1.0\times10^5 で流れている。

  1. 平均流速 UU と Knudsen 数 Kn=λ/D\mathrm{Kn}=\lambda/D を求めよ。分子の自由度を消して流体方程式に移る近似(第8章)は、どの小さいパラメータで閉じているか。
  2. Kolmogorov の長さ η∼D Re−3/4\eta\sim D\,\mathrm{Re}^{-3/4} を見積もれ。η/D\eta/D も小さいのに、Reynolds 平均の方程式がこの比で閉じないのはなぜか。
  3. 壁からの距離 yy の所で、Reynolds 応力を運ぶ渦の大きさを混合距離 ℓm=κy\ell_m=\kappa y(κ≈0.41\kappa\approx0.41)とする。平均流速は対数則 dU/dy∝1/y\mathrm{d} U/\mathrm{d} y\propto1/y に従い、尺度 yy で変わる。消す自由度の大きさと残す自由度の変化の尺度の比を求め、Re→∞\mathrm{Re}\to\infty でこの比がどうなるかを述べよ。
解答

(1) U=Re ν/D=1.0×105×1.5×10−5/0.10=15 m/sU=\mathrm{Re}\,\nu/D=1.0\times10^5\times1.5\times10^{-5}/0.10=15\,\mathrm{m/s}。Kn=6.8×10−8/0.10=6.8×10−7\mathrm{Kn}=6.8\times10^{-8}/0.10=6.8\times10^{-7}。流体方程式への移行で消すのは分子の速度分布の非平衡部分で、それが緩和する長さは平均自由行程 λ\lambda、残す密度・流速・温度が変わる長さは DD である。比 Kn\mathrm{Kn} が 10−610^{-6} 程度と非常に小さいので、Chapman–Enskog 展開は Kn\mathrm{Kn} の冪で整理でき、Navier–Stokes 方程式は B 型の近似として閉じている。

(2) η∼0.10×(105)−3/4≈0.10×1.78×10−4≈1.8×10−5 m\eta\sim0.10\times(10^5)^{-3/4}\approx0.10\times1.78\times10^{-4}\approx1.8\times10^{-5}\,\mathrm{m}、η/D≈1.8×10−4\eta/D\approx1.8\times10^{-4}(検算スクリプト)。η\eta は粘性で散逸する最も小さい渦の大きさである。Reynolds 平均で消すのは揺らぎ ui′u_i' のすべてで、その中には η\eta 程度の小さい渦だけでなく、DD に近い大きさの渦まで含まれる。Reynolds 応力 RijR_{ij} を主に担うのはエネルギーを多くもつ大きい渦の方なので、消す自由度の大きさは η\eta でなく大きい渦の大きさで測らなければならない。η/D\eta/D が小さいことは、散逸の起きる尺度が分離していることを意味するだけで、平均流の方程式を閉じさせる比にはならない。

(3) 比は ℓm/y=κ≈0.41\ell_m/y=\kappa\approx0.41 である。これは yy にも Re\mathrm{Re} にもよらない定数で、Re→∞\mathrm{Re}\to\infty でも 00 にならない。均質化の ε/L\varepsilon /L や流体力学極限の Kn\mathrm{Kn} は、消す自由度が残す自由度に比べて限りなく細かくなる極限をとれたが、壁乱流の Reynolds 応力では、消す渦と平均流の変化がどこでも同じ尺度にある。展開のパラメータになる小さい比が無く、本節の (i) の場合に当たる。混合距離モデルや渦粘性モデルが定理でなく C 型の閉じ方にとどまるのは、このためである。

木下の勉強ステータス

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

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

← AI生成教科書『虎の巻』の記事一覧へ