tecktom.com

何がどこまで証明されているか

2026年10月9日8分で読めますAI生成教科書『虎の巻』第8章 Boltzmann 方程式から流体方程式へ
目次
  1. Reynolds 数・Mach 数・Knudsen 数の関係
  2. 拡散スケーリングと Golse–Saint-Raymond の定理
  3. 粒子系から流体まで
  4. 木下の層

AI 生成・木下が未勉強

近似の数理　第8章 §8.7　★★ 非圧縮 Navier–Stokes が極限定理として得られるスケーリングと、証明の現状を正確に押さえる

Reynolds 数・Mach 数・Knudsen 数の関係

流れの性格は、Reynolds 数 Re=ρUL/μ\mathrm{Re}=\rho UL/\mu（慣性と粘性の比）、Mach 数 Ma=U/cs\mathrm{Ma}=U/c_s（流速と音速の比）、Knudsen 数 Kn=λ/L\mathrm{Kn}=\lambda/L で表される。運動論からは μ∼ρcˉλ\mu\sim\rho\bar c\lambda（cˉ\bar c は分子の平均の速さ）であり、cˉ\bar c と音速 csc_s は同じ桁なので、

Re∼MaKn\mathrm{Re}\sim\frac{\mathrm{Ma}}{\mathrm{Kn}}

が成り立つ（von Kármán の関係）。初等運動論の μ=12ρcˉλ\mu=\frac12\rho\bar c\lambda を使うと係数まで決まり、単原子気体（γ=5/3\gamma=5/3）で Re=πγ/2 Ma/Kn≈1.62 Ma/Kn\mathrm{Re}=\sqrt{\pi\gamma/2}\,\mathrm{Ma}/\mathrm{Kn}\approx1.62\,\mathrm{Ma}/\mathrm{Kn} である（検算スクリプト）。

この関係から、Kn→0\mathrm{Kn}\to0 の極限の取り方によって異なる方程式が出ることが分かる。

  • Ma\mathrm{Ma} を固定して Kn→0\mathrm{Kn}\to0 とすると、Re→∞\mathrm{Re}\to\infty となり、粘性は消えて圧縮性 Euler 方程式が出る（Hilbert 展開）。
  • Ma∼Kn=ε→0\mathrm{Ma}\sim\mathrm{Kn}=\varepsilon \to0 とすると、Re=O(1)\mathrm{Re}=\mathcal{O}(1) に保たれ、粘性が極限で残る。同時に Ma→0\mathrm{Ma}\to0 なので圧縮性は消え、非圧縮性 Navier–Stokes 方程式が出る。

拡散スケーリングと Golse–Saint-Raymond の定理

Ma∼Kn=ε\mathrm{Ma}\sim\mathrm{Kn}=\varepsilon に対応して、Boltzmann 方程式を

ε ∂tfε+v⋅∇xfε=1εQ(fε,fε),fε=M(1+εgε)\varepsilon \,\partial_tf_{\varepsilon }+\bm{v}\cdot\nabla_{\bm{x}}f_{\varepsilon }=\frac{1}{\varepsilon }Q(f_{\varepsilon },f_{\varepsilon }),\qquad f_{\varepsilon }=M(1+\varepsilon g_{\varepsilon })

とスケールする。MM は静止した一様な Maxwell 分布（u=0\bm{u}=\bm{0}、n=T=1n=T=1）、gεg_{\varepsilon } はそこからのずれで、ずれの大きさ ε\varepsilon が Mach 数に当たる。時間は ε−1\varepsilon ^{-1} 倍に引き延ばされている（流体の時間尺度まで見る）ので、拡散スケーリングと呼ばれる。

定理 8.11（非圧縮 Navier–Stokes 極限（主張） ★★）

空間を R3\mathbb{R}^3 または周期的な箱とし、剛体球を含む一群の衝突演算子を考える。fεf_{\varepsilon } を Boltzmann 方程式の再正規化解（DiPerna–Lions による時間大域的な弱い意味の解）で、初期値の相対エントロピーが O(ε2)\mathcal{O}(\varepsilon ^2) に抑えられるものとする。このとき gεg_{\varepsilon } は（部分列で）弱い意味で

g=ρ+u⋅v+θ ∣v∣2−52g=\rho+\bm{u}\cdot\bm{v}+\theta\,\frac{\lvert \bm{v}\rvert ^2-5}{2}

の形の関数に収束し、(u,θ)(\bm{u},\theta) は非圧縮性 Navier–Stokes–Fourier 方程式

∇⋅u=0,∂tu+(u⋅∇)u+∇P=νΔu,∂tθ+u⋅∇θ=κ′Δθ,∇(ρ+θ)=0\nabla\cdot\bm{u}=0,\quad \partial_t\bm{u}+(\bm{u}\cdot\nabla)\bm{u}+\nabla P=\nu\Delta\bm{u},\quad \partial_t\theta+\bm{u}\cdot\nabla\theta=\kappa'\Delta\theta,\quad \nabla(\rho+\theta)=\bm{0}

の Leray の意味の弱解（8.8 節）になる。ν,κ′\nu,\kappa' は線型化衝突演算子 L\mathcal{L} の逆から決まる定数である。

この主張は、Bardos・Golse・Levermore が 1990 年代初頭に提示したプログラム（弱解から弱解への極限として非圧縮流体方程式を導く方針）を、Golse と Saint-Raymond が 2004 年に完成させたものである。ここで「弱い意味」とは、関数そのものが各点で収束するのでなく、滑らかな試験関数 ψ\psi との積分 ∫gεψ\int g_{\varepsilon }\psi が収束することを言う（8.9 節で詳しく述べる）。最後の式 ∇(ρ+θ)=0\nabla(\rho+\theta)=\bm{0} は Boussinesq の関係と呼ばれ、極限で圧力の揺らぎ ρ+θ\rho+\theta が一様になること（音波が消えること）を表す。

この定理の位置

図 8.1 の「Boltzmann →\to 非圧縮性 NS」の矢印は、時間大域的に、大きなデータ（初期値の大きさに小ささの仮定がない）で、弱い意味で証明されている。代わりに得られる流体方程式の解も弱解であり、その一意性や滑らかさは分からない。これは次節の Navier–Stokes 方程式自体の未解決問題を引き継いでいる。

粒子系から流体まで

図 8.1 を上から下まで一本の定理で結ぶには、Boltzmann 方程式の導出（Lanford の定理）が流体の時間尺度まで延びている必要がある。拡散スケーリングでは平均自由時間の ε−2\varepsilon ^{-2} 倍の時間を見るので、平均自由時間の数分の 1 しか保証しない Lanford の定理では届かない。8.3 節の注意で述べたとおり、2025 年に Deng・Hani・Ma のプレプリントが、Boltzmann 方程式の長時間導出を経て流体方程式まで到達したと主張している。これが確定すれば Hilbert の第6問題のこの部分に一つの解答を与えることになるが、査読の状況は本書では確定していない。

例題 8.8　空気の流れで Reynolds 数・Mach 数・Knudsen 数を見積もる ★★

オリジナル初等運動論の粘性率 μ=12ρcˉλ\mu=\frac12\rho\bar c\lambda、平均の速さ cˉ=8kBT/(πm)\bar c=\sqrt{8 k_{\mathrm B} T/(\pi m)}、音速 cs=γkBT/mc_s=\sqrt{\gamma k_{\mathrm B} T/m} を使う。20 ∘C20\,{}^\circ\mathrm{C}、11 気圧の空気では γ=1.4\gamma=1.4、cs≈343 m/sc_s\approx343\,\mathrm{m/s}、動粘性率 ν=μ/ρ≈1.50×10−5 m2/s\nu=\mu/\rho\approx1.50\times10^{-5}\,\mathrm{m^2/s}、平均自由行程（μ\mu から逆算した値）λ≈65 nm\lambda\approx65\,\mathrm{nm} とする。

  1. Re=ρUL/μ=πγ/2  Ma/Kn\mathrm{Re}=\rho UL/\mu=\sqrt{\pi\gamma/2}\;\mathrm{Ma}/\mathrm{Kn} を示せ。
  2. 長さ L=1 cmL=1\,\mathrm{cm} で Re=100\mathrm{Re}=100 となる流速 UU と、そのときの Ma\mathrm{Ma}、Kn\mathrm{Kn} を求めよ。
  3. 長さ L=1 mL=1\,\mathrm{m} で Ma=0.3\mathrm{Ma}=0.3 の流れの Re\mathrm{Re} を求めよ。
  4. (2)(3) はそれぞれ本節の二つの極限のどちらに近いか。
解答

(1) Re=ρUL12ρcˉλ=2 Ucs Lλ cscˉ=2 MaKnγkBT/m8kBT/(πm)=πγ2 MaKn\mathrm{Re}=\dfrac{\rho UL}{\frac12\rho\bar c\lambda}=2\,\dfrac{U}{c_s}\,\dfrac{L}{\lambda}\,\dfrac{c_s}{\bar c}=2\,\dfrac{\mathrm{Ma}}{\mathrm{Kn}}\sqrt{\dfrac{\gamma k_{\mathrm B} T/m}{8 k_{\mathrm B} T/(\pi m)}}=\sqrt{\dfrac{\pi\gamma}{2}}\,\dfrac{\mathrm{Ma}}{\mathrm{Kn}}。空気では係数は 0.7π≈1.48\sqrt{0.7\pi}\approx1.48 である。

(2) U=Re ν/L=100×1.50×10−5/0.01≈0.15 m/sU=\mathrm{Re}\,\nu/L=100\times1.50\times10^{-5}/0.01\approx0.15\,\mathrm{m/s}。Ma=0.15/343≈4.4×10−4\mathrm{Ma}=0.15/343\approx4.4\times10^{-4}、Kn=65×10−9/0.01≈6.5×10−6\mathrm{Kn}=65\times10^{-9}/0.01\approx6.5\times10^{-6}。比は Ma/Kn≈67\mathrm{Ma}/\mathrm{Kn}\approx67 で、1.48×67≈1001.48\times67\approx100 と (1) に合う。

(3) U=0.3×343≈103 m/sU=0.3\times343\approx103\,\mathrm{m/s}、Kn=6.5×10−8\mathrm{Kn}=6.5\times10^{-8} なので Re≈1.48×0.3/(6.5×10−8)≈6.9×106\mathrm{Re}\approx1.48\times0.3/(6.5\times10^{-8})\approx6.9\times10^{6}。

(4) (2) は Ma\mathrm{Ma} と Kn\mathrm{Kn} がともに小さく、比が O(102)\mathcal{O}(10^2) にとどまって Re\mathrm{Re} が有限な流れであり、Ma∼Kn→0\mathrm{Ma}\sim\mathrm{Kn}\to0 の拡散スケーリング（非圧縮 Navier–Stokes）の側にある。(3) は Ma\mathrm{Ma} を有限に保ったまま Kn\mathrm{Kn} が非常に小さく、Re\mathrm{Re} が 10710^7 近い。これは Kn→0\mathrm{Kn}\to0 で粘性が消える側（Euler）であり、粘性は壁の近くの薄い境界層にだけ残る。空気の λ\lambda が短いために、日常の大きさの流れはほとんどが後者に入る。粘性が主役の流れを得るには、流れを遅く小さくして Ma\mathrm{Ma} を Kn\mathrm{Kn} と同じ程度まで下げる必要がある（検算スクリプト）。

木下の勉強ステータス

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

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

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

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

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