tecktom.com

階層の全体像と Hilbert の第6問題

2026年10月9日8分で読めますAI生成教科書『虎の巻』第8章 Boltzmann 方程式から流体方程式へ
目次
  1. 三つの記述と二つの極限
  2. Hilbert の第6問題
  3. 木下の層

AI 生成・木下が未勉強

近似の数理　第8章 §8.1　★★★ 章全体の地図。どの矢印がどの極限で、どれが証明済みかを最初に押さえる

木下の問いに出てくる「流体力学での近似」には二つの層がある。表面の層は、Navier–Stokes 方程式を解くときの近似（境界層・乱流モデル・圧縮性の無視など）である。その下にもう一つ層があり、Navier–Stokes 方程式自体が「分子の集まりを連続体とみなす」という近似の産物である。この章は下の層を扱う。分子の運動（Newton 力学）から Boltzmann 方程式を経て流体方程式に至る道筋を、各段で何を小さいとみなし、どういう極限を取り、その極限が数学的にどこまで正当化されているかという観点で追う。

結論を先に書いておく。流体方程式は小さいパラメータ（Knudsen 数）の極限で現れる有効方程式であり、第1章の分類では典型的な B 型の近似である。ただしその途中には、Boltzmann 方程式を導く段（分子的混沌の仮定）と、流体方程式を閉じる段（応力と熱流束を勾配で書く段）の二か所に「閉じる」操作が入っている。前者は極限定理（Lanford の定理）によって短時間だけ正当化され、後者は Hilbert 展開・Chapman–Enskog 展開によって小さい Knudsen 数での漸近展開として正当化される。さらに、導かれた Navier–Stokes 方程式そのものについて、3 次元で滑らかな解が時間大域的に存在するかどうかが分かっていない。近似の正当化の連鎖は、最後の環で数学の未解決問題につながっている。

三つの記述と二つの極限

気体を記述する方法は、細かいほうから順に三つある。

  1. 粒子の記述。NN 個の分子の位置 x1,…,xN\bm{x}_1,\dots,\bm{x}_N と速度 v1,…,vN\bm{v}_1,\dots,\bm{v}_N を Newton の運動方程式で追う。初期条件に確率分布を入れると、6N6N 次元の位相空間上の確率密度 FN(t,x1,v1,…,xN,vN)F_N(t,\bm{x}_1,\bm{v}_1,\dots,\bm{x}_N,\bm{v}_N) の時間発展は Liouville 方程式で書ける。これは厳密だが、N∼1023N\sim 10^{23} では扱えない。
  2. 運動論的記述。1 粒子の位置と速度の分布 f(t,x,v)f(t,\bm{x},\bm{v})（6 次元の位相空間上の関数）だけを追う。これが Boltzmann 方程式の未知関数である。
  3. 流体の記述。各点での密度 ρ(t,x)\rho(t,\bm{x})、流速 u(t,x)\bm{u}(t,\bm{x})、温度 T(t,x)T(t,\bm{x}) だけを追う。3 次元空間上の 5 個の関数であり、Euler 方程式・Navier–Stokes 方程式の未知関数である。

1 から 2 へ移る段は、Liouville 方程式から 1 粒子分布の方程式を作る段であり、第6章で見た BBGKY 階層（kk 粒子分布の方程式に k+1k+1 粒子分布が現れる連鎖）を途中で切る操作になる。ここで使う極限は Boltzmann–Grad 極限と呼ばれ、粒子数 N→∞N\to\infty、粒子の直径 d→0d\to 0 を、Nd2Nd^2 が一定になるように取る（容器の体積は固定）。Nd2Nd^2 は 1 個の粒子が単位長さ進む間に出会う衝突の回数に比例するので、この極限では平均自由行程が一定に保たれる。一方で粒子の体積の総和 Nd3∼(Nd2) dNd^3\sim (Nd^2)\,d は 00 に行くので、気体は限りなく希薄になる。

2 から 3 へ移る段で使う小パラメータは Knudsen 数

Kn=λL\mathrm{Kn}=\frac{\lambda}{L}

である。λ\lambda は平均自由行程（分子が衝突から次の衝突までに進む平均距離）、LL は流れの巨視的な長さ（管の径、物体の大きさ、密度や温度が変化する距離）である。Kn→0\mathrm{Kn}\to 0 は、巨視的な変化が起きる距離の中で分子が非常に多く衝突し、各点で局所的な熱平衡に達している状況を表す。

図

図 8.1　気体の記述の階層。実線は本章で扱う極限、破線は流体方程式の中での低 Mach 数極限。

Hilbert の第6問題

1900 年に Hilbert が挙げた 23 の問題（パリの国際数学者会議ではそのうち 10 問を講演し、全体は論文として発表した）のうち、第6問題は「物理学の公理を数学的に扱うこと」である。その説明の中で Hilbert は、Boltzmann の気体論の仕事に触れ、原子論的な見方から連続体の運動法則へ至る極限過程を数学的に厳密に展開することを課題として挙げた。図 8.1 の上から下までの矢印を、すべて定理として結ぶことがこの課題の一つの読み方である。Hilbert 自身は 1912 年に Boltzmann 方程式を論じ、後で述べる Hilbert 展開を導入した。

表 8.1　階層の各矢印。詳細は本章の各節で述べる。
矢印小パラメータ何を捨てるか数学的な状況
Liouville →\to BBGKYなし何も捨てない（周辺化）厳密な恒等式
BBGKY →\to Boltzmannd→0d\to0（Nd2Nd^2 一定）衝突前の 2 粒子相関（分子的混沌）Lanford（1975）が平均自由時間の数分の 1 の時間まで証明。長時間は 2025 年のプレプリントが主張
Boltzmann →\to EulerKn→0\mathrm{Kn}\to0局所平衡からのずれ滑らかな Euler 解が存在する時間内で証明済み
Boltzmann →\to NS（圧縮性）Kn\mathrm{Kn} の 1 次までO(Kn2)\mathcal{O}(\mathrm{Kn}^2) の項漸近展開としての整合性。極限定理の形にはならない（後述）
Boltzmann →\to 非圧縮性 NSMa∼Kn→0\mathrm{Ma}\sim\mathrm{Kn}\to0圧縮性・非線形の高次弱解の意味で証明済み（Golse–Saint-Raymond 2004）

例題 8.1　Boltzmann–Grad 極限で何が小さくなるか ★★

オリジナル直径 dd の剛体球が数密度 nn で詰まった気体を考える。平均自由行程は λ=1/(2 nπd2)\lambda=1/(\sqrt2\,n\pi d^2)、球が占める体積の割合（体積分率）は ϕ=n⋅π6d3\phi=n\cdot\frac{\pi}{6}d^3 である。

  1. ϕ=d62 λ\phi=\dfrac{d}{6\sqrt2\,\lambda} を示し、Boltzmann–Grad 極限（λ\lambda を固定して d→0d\to0）で ϕ→0\phi\to0 となることを確かめよ。
  2. 20 ∘C20\,{}^\circ\mathrm{C}、11 気圧の空気（n≈2.50×1025 m−3n\approx2.50\times10^{25}\,\mathrm{m^{-3}}、d=0.37 nmd=0.37\,\mathrm{nm}）について ϕ\phi と λ/d\lambda/d を求めよ。
  3. 液体の水（数密度 n≈3.34×1028 m−3n\approx3.34\times10^{28}\,\mathrm{m^{-3}}、分子の有効直径 d≈0.28 nmd\approx0.28\,\mathrm{nm} とする）について同じ量を求め、Boltzmann 方程式が使えるかどうかを判断せよ。
解答

(1) λ\lambda の式から nπd2=1/(2λ)n\pi d^2=1/(\sqrt2\lambda) なので

ϕ=π6nd3=d6 nπd2=d62 λ.\phi=\frac{\pi}{6}nd^3=\frac{d}{6}\,n\pi d^2=\frac{d}{6\sqrt2\,\lambda}.

λ\lambda を固定すれば ϕ\phi は dd に比例して 00 に行く。平均自由行程（衝突の頻度）を保ったまま、粒子が占める体積を消していく極限であり、これが「希薄」の意味である。

(2) λ=1/(1.414×2.50×1025×π×(0.37×10−9)2)≈6.6×10−8 m\lambda=1/(1.414\times2.50\times10^{25}\times\pi\times(0.37\times10^{-9})^2)\approx6.6\times10^{-8}\,\mathrm{m} なので λ/d≈1.8×102\lambda/d\approx1.8\times10^2、ϕ=0.37/(6×1.414×66)≈6.6×10−4\phi=0.37/(6\times1.414\times66)\approx6.6\times10^{-4}。空間の 0.07%0.07\% ほどしか分子が占めておらず、分子は自分の直径の約 180180 倍の距離を進んでから次の分子に当たる。2 個の分子がちょうど衝突している瞬間に第 3 の分子が近くにいる確率は ϕ\phi 程度なので、3 体衝突は無視でき、衝突は 2 体の独立な出来事として扱える。

(3) λ=1/(1.414×3.34×1028×π×(0.28×10−9)2)≈8.6×10−11 m\lambda=1/(1.414\times3.34\times10^{28}\times\pi\times(0.28\times10^{-9})^2)\approx8.6\times10^{-11}\,\mathrm{m} で、λ/d≈0.31\lambda/d\approx0.31 と直径より短い。ϕ≈0.38\phi\approx0.38 で、空間の 4 割近くを分子が占める。分子は常に複数の隣の分子と接していて、「自由に飛んでまれに 2 体で衝突する」という Boltzmann 方程式の前提が成り立たない。表 8.1 の BBGKY から Boltzmann への矢印は液体には使えず、液体の流体方程式（粘性率など）は運動論とは別の方法で求める必要がある（検算スクリプト）。

木下の勉強ステータス

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

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

目次PDF第8章の PDF近似の数理の PDF（全章）

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

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