Kinetic derivation of thermal viscous models for nematic liquid crystal dynamics

技術要約：ネマチック液晶ダイナミクスにおける熱粘性モデルの動力学的導出

問題提起
近年の研究により、動力学理論から導出された、圧縮性かつ非粘性のネマチック液晶用レイスリー・エリクセン方程式が確立された。しかし、その手法は、すでに配向した状態の周囲でのモーメント閉鎖に依存しており、配向がダイナミクスから創発することを許容するのではなく、集団的な配向挙動をあらかじめ課していた。さらに、得られたモデルは厳密に非粘性であり、実際の液晶の流れに内在する散逸、熱、およびスピン拡散の効果を捉えることができなかった。本論文は、配向の創発を考慮し、熱および粘性効果を含む、圧縮性粘性版のレイスリー・エリクセン方程式を導出することによって、これらの限界に対処するものである。

手法
著者らは、秩序化された流体の動力学理論から出発して、マクロな熱力学理論を展開する。その手法は以下の段階を経て進行する。

1. 動力学的枠組み： 出発点は、Vlasovポテンシャルを介した平均場配向相互作用を組み込んだ、[7]で提案された秩序化された流体の動力学理論である。系は、配向緩和が並進運動量緩和よりも速いという時間スケールの分離を仮定している。分布関数 $f(x, p, \nu, \varsigma, t)$ は、位置、線形運動量、微視的構造構成（配向 $\nu$）、および共役運動量 $\varsigma$ に依存する。
2. 衝突演算子： 解析を扱いやすくするため、完全なボルツマン衝突演算子の代わりに、単一緩和BGK（Bhatnagar–Gross–Krook）演算子を用いる。この演算子は、完全な動力学理論の保存則とエントロピー生成構造を保持する。
3. バランス則： 質量、線形成分、角運動量、エネルギー、およびエントロピーに関するマクロなバランス方程式を、動力学方程式のモーメントを取ることによって導出する。重要な特徴は、角運動量のバランスにおける平均場トルク項の特定であり、これが配向の創発を駆動する。
4. 熱力学的構造： ヘルムホルツ自由エネルギーを熱力学的ポテンシャルとして特定する。著者らは、平衡分布（マクスウェル分布）がこの自由エネルギーを最小化することを証明し、関連するギブス分布を配向プロファイルに対して導出する。特定の相互作用ポテンシャル（Kuramoto型およびMaier–Saupe型）の下で、温度の低下に伴い系が等方相からネマチック相へと相転移することを実証し、配向した平衡状態においてVlasovトルクが消失することを示す。
5. チャップマン・エンスコッグ展開： 局所的なマクスウェル平衡の周囲で、一次のチャップマン・エンスコッグ展開を行う。標準的な導出とは異なり、展開における時間微分は、標準的な圧縮性オイラー方程式ではなく、ゼロ次バランス則（エリクセン応力およびカップル応力項を含む）を用いて代入される。これにより、配向弾性が一次補正へと正しく伝播することが保証される。
6. 構成則の閉鎖： 著者らは、RajagopalとSrinivasaによる制約付き最大化原理を用いる。動力学展開から得られたエントロピー生成率と、ヘルムホルツ自由エネルギーから推論される熱力学的構造を組み合わせる。バランス則の条件下でエントロピー生成を最大化することにより、コーシー応力、カップル応力、および熱流束の構成関係を導出する。

主要な貢献と結果

• 創発的配向： 配向を仮定していた従来の非粘性モデルとは異なり、本導出では、Vlasovトルクを介して等方・ネマチック転移が自然に創発することを許容する。解析により、平均場トルクが配向したギブス平衡において消失することが確認され、マクロなバランス則との整合性が確保されている。
• 粘性と熱効果： 圧縮性、粘性、熱伝導、およびスピン拡散特性を持つネマチック流体の完全な構成方程式のセットを導出する。導出された応力テンソルには以下が含まれる：
  • ニュートン粘性剪断粘性。
  • 並進・回転エネルギーの不均衡に由来する等方的な体積粘性。
  • エリクセンテンソルおよび速度勾配に結合した異方性粘性応力。
  • スピン拡散（Eringen型のマイクロポーラー項）およびディレクター勾配への結合を含むカップル応力テンソル。
• 輸送係数： 輸送係数（剪断粘性、熱伝導率、体積粘性、異方性粘性、およびスピン粘性）の明示的な表現が、BGK緩和時間 $\tau$ およびフランク弾性係数 $K$ を用いて導出される。熱伝導率は二原子分子ガスにおけるEuckenの関係を満たす。
• 非圧縮極限： 最大化の手順中にラグランジュ乗数を用いて非発散制約を課すことにより、本フレームワークは、別途の導出を必要とせずに、自然に非圧縮型の対応物をもたらす。
• 構造的新規性： 導出された粘性応力テンソルは、単一の粘性パラメータ（等方的なBGK近似による）によって支配されており、完全な5つの係数を持つレイスリー粘性テンソルを回復しない。著者らは、完全なレイスリーテンソルを回復するには、複数の緩和時間スケールを持つ異方的なBGK演算子が必要であることを指摘している。これは継続中の課題である。

意義と主張
本論文は、Farrell, Russo, および Zerbinati [1] による最近の非粘性動力学理論を、粘性、熱、およびスピン拡散効果を考慮するように一般化することを主張している。主な意義は、秩序化された流体における微視的な粒子相互作用とマクロな連続体力学との間の溝を埋め、これら散逸効果を動力学理論から熱力学的に整合した形で直接導出したことにある。

著者らは、配向をあらかじめ課すのではなく、自由エネルギーの最小化と平衡におけるVlas_ovトルクの消失を通じて配向が創発するというアプローチによって、従来モデルの限界を克服していることを強調している。得られた系（式 8.21）は、レイスリー・エリクセン方程式の圧縮性、粘性、熱伝導、およびスピン拡散版を代表するものである。非粘性極限（$\tau \to 0$）において、このモデルは [1] で導出された非粘性系へと還元される。

本研究は、より完全な液晶の動力学理論へのステップとして提示されており、著者らは、5つの独立した粘性を持つ完全なレイスリー粘性テンソルの回復には、より複雑な異方性衝突演算器が必要であるという未解決の問題についても謙虚に言及している。