Mechanical Visualization of FLRW-Type Expansion Kinematics in Hele-Shaw Compression Flow: Tracer Trajectories, Comoving Scalar Modes, and Anisotropic Cells

技術要約：Hele-Shaw圧縮流におけるFLRW型膨張運動学の機械的可視化

問題提起
本論文は、フリードマン・ルメートル・ロバートソン・ウォーカー（FLRW）宇宙論に見られる計量膨張の運動学を、正確に可視化するマクロなベンチスケールの機械的システムの構築という課題に取り組んでいる。これまで、有効音響計量を通じて相対論的地平線や粒子生成をモデル化するために流体システムが用いられてきたが、本研究は、膨張の具体的な運動学的特徴であるスケール因子 $a(t)$、膨張率（ハッブルパラメータ）、共動分離、および重畳されたスカラー場の輸送を再現することを目指している。核心となる問題は、圧縮されたHele-Shaw流（スクイーズフロー）における物理的観測量が、質量保存則から導出される深さ平均の運動学的法則に厳密に従うものと、横方向の速度勾配の影響を受ける物理的なトレーサー粒子の挙動とを、どのように区別するかである。

手法
本研究では、解析的導出、潤滑方程式の厳密解、および高精度な三次元計算流体力学（CFD）を組み合わせて用いている。

1. 解析的枠組み： システムは、時間依存するギャップ $h(t)$ を持つ水平方向の軸対称（または楕円形）セル内における非圧縮ニュートン流体としてモデル化される。薄いギャップ（$\epsilon = h/R \ll 1$）およびクリーピング流（$Re_{squeeze} \ll 1$）の極限において、半径方向の速度プロファイルはギャップ方向に放物線状として導出される。
2. トレーサーの動力学： 物理的トレーサーの軌跡は、半径方向と垂直方向の運動方程式を連立させて解くことで解析される。これには、高さ依存の膨張率と中間面の安定性の決定が含まれる。
3. スカラー輸送： 受動的スカラー（温度や染料など）の輸送は、共動座標系で解析される。研究では、移流拡散方程式を、計量によって抑制された拡散率を持つ純粋な拡散へと簡約化し、波長の伸長と振幅の減衰を計算する。
4. 異方的アナログ： 幾何学的構造を楕円形セルへと拡張し、Bianchi Type-Iの異方性膨張をモデル化する。圧力場および速度場の厳密解を導出し、CFDシミュレーションによって検証する。
5. 数値的検証： 三次元CFDシミュレーション（OpenFOAMを使用）を行い、異方的スケーリング則を検証し、クリーピング流近似の限界をテストし、有限レイノルズ数による補正を定量化する。

主な貢献と結果

• 深さ平均場とトレーサーの区別：
  本論文は、厳密な共動スケーリング則 $r(t) = a(t)r_0$（ただし $a(t) = \sqrt{h_0/h(t)}$）が、物理的トレーサー粒子ではなく、深さ平均速度場に厳密に適用されることを確立した。

  • トレーサーの挙動： 物理的トレーサーは、高さ依存の膨張率 $H(\zeta) = 6\zeta(1-\zeta)\bar{H}$ を経験する。異なる高さにあるトレーサーは異なる速度で遠ざかり、$a=7.4$ への圧縮において半径で3.3倍の差が生じる。
  • 中間面アトラクター： 中間面（$\zeta = 0.5$）は漸近的に安定なアトラクターである。面からの変位は $a^{-1}$ として減衰する。中間面上のトレーサーは、修正された法則 $r(t) = a(t)^{3/2}r_0$ に従う。
  • 測定への示唆： 深さ分解能を持たない標準的な粒子追跡は、誤ったスケール因子をもたらす。中間面に焦点を当てたライトシート流速計（µPIV）のみが、クリーンな冪乗則の対応（$h_0/h$ に対して指数 3/4）を回収できる。

• 深さ平均法則の物理的実現：
  深さ平均FLRWの対応関係は、セルが膨張するよりも速くギャップを横断して混合する観測量によってのみ、物理的に実現される。これには、条件 $\Pi = \bar{H}h^2/\kappa \ll 1$（ここで $\kappa$ は拡散係数）が必要である。

  • 粒子トレーサーはこの条件を満たさない。
  • ミリメートル単位のギャップにおける分子染料も、3桁分条件を満たしていない。
  • 薄いギャップ（例：$h=100\,\mu\text{m}$）における熱的スカラーは、この条件を満たしており、深さ平均場に対する正しい物理的アナログとなる。

• 共動スカラーモードと減衰：
  よく混合されたスカラーについて、共動座標系における輸送方程式は、時間依存する拡散率 $\kappa/a^2(t)$ を持つ純粋な拡散へと簡約化される。

  • 波長の伸長： 波長は $\lambda(t) = a(t)\lambda_0$ として伸長する。
  • シル型減衰： 振幅は $A(t) = A_0 \exp\left(-\kappa k_0^2 \int a^{-2} dt\right)$ に従って減衰する。
  • 凍結基準： 残留振幅が生き残るためには、スケール因子が漸近的に $t^{1/2}$ よりも速く成長する必要がある。放射型の場合（$a \propto t^{1/2}$）は限界的であり、完全な消去をもたらす。

• 異方的（Bianchi Type-I）アナログ：
  楕円形セル（半軸 $a_e$ および $b_e$）において、流れは厳密に均質な異方的ハッブル法則を生成する：
  $$\frac{H_x}{H_y} = \frac{b_e^2}{a_e^2}$$
  この比率は、質量保存則のみからではなく、潤滑解の真の予測である。

  • CFD検証： 三次元シミュレーションにより、クリーピング流の極限において、この比率が0.1%以内の精度で一致することが確認された。
  • 有限レイノルズ効果： この比率は、単位スクイーズ・レイノルズ数 $Re_{squeeze}$ あたり約 $+15\%$ の補正を受けることが定量化されており、経験的なフィットによって記述される。この異方性は、ダイナミックなBianchi-I宇宙論における剪断（シア）の進化とは異なり、境界幾何学によって「凍結」されている。

• 設計ウィンドウ：
  本論文は、すべての制約（薄いギャップ、クリーピング流、横断混合、およびテイラー・アリス分散の限界）が同時に満たされる、狭いが空ではない設計空間を定義している。実行可能な構成として、半径 $50\,\text{mm}$ のセルにおける $100\,\mu\text{m}$ ギャップ内の熱的スカラーが特定された。

意義と範囲
本論文は、正しい観測量（深さ平均場またはよく混合されたスカラーであり、弾道的なトレーサーではない）に対して適用される限り、定量的に正確な機械的可視化を提供するものであると主張している。

その意義は、膨張の履歴 $h(t)$ を機械的に「プログラム」し、その結果としての運動学的効果（後退、波長の伸長、減衰、および異方性）をベンチトップで直接観察できる能力にある。このシステムは、スケール因子、膨張率、および共動輸送を再現する運動学的アナログとして機能する。

明示的な限界：
著者は、このシステムがモデル化しない事項を明示している：

• これは純粋に運動学的であり、有効計量、波動方程式、または因果的地平線は存在しない。
• アインシュタイン方程式やフリードマン方程式を再現しない（膨張は規定されるものであり、応力エネルギーによって決定されるものではない）。
• 真のシル減衰（音響振動を伴うもの）ではなく、単なるモードの拡散的消去をモデル化している。
• 異方性は幾何学的に凍結されており、動的に進化するものではない。
• 実験的な実現はまだ行われておらず、設計ウィンドウは理論的な提案の段階である。

要約すると、本研究は流体力学と宇宙論的運動学を架け橋し、どのような条件下でスクイーズされた粘性流が宇宙の膨張を模倣するかを厳密に定義することで、将来の教育的または実験的なデモンストレーションのための精密な枠組みを提供している。