Mixing of miscible liquids: Dimensionless scaling for intermediate-to-large density differences in a stirred tank

本研究は、混合液を有する攪拌槽の数値シミュレーションを用いて、混合時間はリチャードソン数と正の相関を示す一方で、動力数、フルード数、およびリチャードソン数に基づいて導出された指数関数的スケーリング則が、中程度から大きな密度差を有するすべてのデータを単一のマスター曲線上に集約することを示す。

要約

巨大なタンクに、底部に重く粘り気のあるシロップ、上部に軽くて薄いジュースという2種類の液体を混ぜようとしている状況を想像してください。これらを混ぜ合わせるために、巨大な回転式パドル（インペラ）を投入します。

現実世界では、医薬品から廃水処理に至るまで、あらゆる分野の工場においてこれは一般的な作業です。しかし、ここには一つの難題があります。液体の密度（重さ）が異なるため、重い液体は底部に留まり、軽い液体は上部に浮かびたがります。これにより、混合しようとする回転パドルと、分離させようとする重力との間で「格闘」が生じます。

本論文は、科学者たちが強力なコンピュータシミュレーションを用いて、これら2種類の液体を完全に混合するのに正確にどれだけの時間がかかるかを突き止め、毎回物理的なタンクを建設して高価な試験を行うことなくその時間を予測する方法を解明する探偵物語のようなものです。

設定：デジタル実験キッチン

研究者たちは、標準的な産業用混合タンクの仮想モデルを構築しました。

• タンク: 壁と、液体が怠惰な川のように単に円を描いて回転するのを防ぐための4本の垂直なフィン（バッフル）を備えた大きな円筒です。
• パドル: 中央にある回転するブレードです。
• 液体: 重い液体と軽い液体を50対50で混合したものをシミュレーションしました。実際の化学物質を使用したのではなく、同じ粘度（粘性）を持つ「重い」流体と「軽い」流体として扱いました。
• 手法: 標準的な数学方程式の代わりに、格子ボルツマン法と呼ばれる巧妙なトリックを使用しました。これは、液体を連続した塊としてではなく、衝突し合う無数の目に見えないビリヤードの玉としてシミュレーションするものです。これにより、乱流（カオス的な渦）がどのように振る舞うかを正確に把握できました。

核心となる問い：混合はどれほど速く行えるか？

主な目的は、混合時間を予測する「魔法の公式」を見つけることでした。

• 変数: 彼らは2つの主要な要素を変更しました。
  1. パドルの回転速度（レイノルズ数）：回転が速いほど乱流が増え、混合が速くなります。
  2. 重量の差異（リチャードソン数）：液体の重量がほぼ同じであれば混合は容易です。一方、一方が他方よりも著しく重い場合、重力が混合に抵抗し、崩壊しにくい層を形成します。

発見：「重力対回転」の戦い

研究者たちはいくつかの興味深いパターンを発見しました。

1. 重力が関係しない場合（重量が同じ）:
   2種類の液体の重量が完全に同じであれば、混合時間は驚くほど一定です。パドルをどの程度速く回転させても（一定の範囲内であれば）、「無次元混合時間」（つまり「パドルが何回転する必要があるか」という意味の洗練された表現）は約20回転で一定に保たれます。これは自然の法則のようなものです。液体が十分に攪拌されれば、パドルの回転数を増やしても、回転数換算での混合速度はそれ以上速くならないのです。

2. 重力が抵抗する場合（重量が異なる）:
   液体の重量が異なると、重い液体は底部に留まりたがります。重量の差が大きいほど、混合は困難になります。

   • 傾向: 重量の差が大きいほど、混合にかかる時間は長くなります。
   • 意外な展開: 「重量の差」を一定に保ったままパドルを速く回転させても、混合時間が常に短くなるわけではありません。場合によっては、速く回転させることで、特定の混合状態に達するまでの時間がむしろ長くなることがあります。
   • なぜか？ 重い液体を分厚い毛布だと想像してください。パドルを速すぎると回転させると、多くのエネルギーが生まれますが、重い液体は安定した「蓋」を形成し、軽い液体が浸透できなくなります。エネルギーは上部の層を渦巻かせることに浪費され、底部の層は封じ込められたままになります。これは、重い野菜が底部に固まった塊として沈んでいる鍋のスープを混ぜようとするようなものです。スプーンを速く回しても、上部のブロッチを跳ねるだけで、野菜の塊を崩すことはできません。

解決策：新しい「マスターカーブ」

チームの最大の成果は、これらすべての要素を組み合わせた単一のシンプルな公式を作成したことでした。彼らは、混合時間をパワー、フルード数、リチャードソン数という3つの特定の数値の観点から見ると、すべての散らばったデータ点が1つの滑らかな指数曲線に収束することに気づきました。

次のように考えてみてください。以前は、エンジニアが新しい液体がどのように混合するかを判断するには、推測するか、数百回もの試験を行う必要がありました。しかし現在、彼らには「レシピ」があります。重量の差と回転速度を伝えれば、この公式は混合時間を高い精度で予測します。

結論

本論文は、これらの特定の産業用タンクについて以下の点を結論付けています。

• 乱流が鍵である: 液体が完全に攪拌されれば、混合の挙動は予測可能です。
• 重力が支配者である: 液体の密度が異なれば、重力は混合に抵抗する「層化」を生み出します。
• 速いことが常に良いわけではない: 密度差が大きいシステムでは、単にモーターの速度を上げても混合が速くなるとは限りません。場合によっては、より安定した分離を生み出すだけです。

著者らは、高価なプロトタイプを最初に構築する必要なく、エンジニアがより良い混合プロセスを設計できるよう、この新しい公式を提供しています。彼らは将来、この公式を異なるタンク形状やパドルの種類でテストする計画ですが、現時点ではシミュレーションした標準的なタンクに対して完璧に機能しています。

技術概要

技術的概要：撹拌槽における中～大密度差に対する無次元スケーリング

問題提起
互溶性液体の混合は、製薬生産から廃水処理に至るまで、工業プロセスにおける重要な単位操作である。単相・均質系における乱流混合はよく特徴付けられているが、流体間の密度差が存在すると、浮力駆動の成層が生じる。この現象は、重力が対流輸送と拮抗する垂直方向の層状領域を形成し、垂直混合を抑制して均質化に要する時間を増大させる可能性がある。単相トレーサー混合や特定の成層ケースに対する経験相関は存在するが、異なる流れ領域や幾何形状にわたって重力、慣性、粘性、および浮力の効果を一貫して考慮する、一般的に有効な枠組みは未だ見出されていない。本研究は、中～大の密度差を持つ二つの互溶性液体を含む仕切り付き撹拌槽における混合時間の予測という課題に取り組む。

手法
著者らは、格子ボルツマン法（LBM）を用いた計算流体力学（CFD）により、仕切り付き撹拌槽内の乱流混合をシミュレーションした。システムは、動粘度が同一だが密度（$\rho_l$ および $\rho_h$）が異なる二つの互溶性液体を 50/50 の体積比で構成し、軽い流体が重い流体の上に初期成層された状態で構成された。

• 数値アプローチ: シミュレーションでは、D3Q19 格子スキームを用いて非圧縮性ナビエ - ストークス方程式を解いた。高いレイノルズ数（$Re > 10^5$）を扱うため、サブグリッド乱流をモデル化するためにスマゴリンスキー - リリー型大渦シミュレーション（S-L LES）モデルを実装した。
• スカラー輸送: 軽い液体の濃度を追跡するために、質量流束ベースのスカラー輸送モデルを採用した。このモデルは、濃度依存性の重力体積力を通じてスカラー場を運動量方程式と結合させ、局所的な密度変動に起因する浮力効果をシミュレート可能にした。
• シミュレーション空間: インペラ回転速度（$N$）と軽い相の密度（$\rho_l$）を変化させる体系的なパラメータ研究を実施した。これにより、レイノルズ数（$Re$）が 50,000 から 125,000、リチャードソン数（$Ri$）が 0 から 14 という広範な無次元数の範囲が得られた。
• 検証: 数値精度は、グリッド収束指数（GCI）分析により検証され、生産グリッドにおける数値的不確実性は約 2.39% と推定された。さらに、シミュレーション結果は文献からの実験的トレーサーデータと比較され、プローブ応答および混合時間の定義において良好な一致を示した。

主要な結果
本研究では、90% 均質化に要する時間として定義される混合時間（$t_m$）を定量化し、無次元形式（$t^*_m = N \cdot t_m$）で表現した。

1. *均質ベースライン（$Ri = 0$）:** 密度差がない場合、無次元混合時間は調査されたレイノルズ数全体にわたって一定値$t^_m \approx 20$ に収束した。これは、流れが完全に乱流であり、慣性混合が支配的であることを示しており、仕切り付きタンクに関する既存の文献と一致する。
2. 成層の影響（$Ri > 0$）: リチャードソン数が増加し、密度差が大きくなるにつれて、無次元混合時間は系統的に増加した。データはすべてのレイノルズ数において明確な指数関数的傾向に従った。
3. 提案されるスケーリング則: 著者らは、すべてのシミュレーションデータを単一の指数スケーリング関係に収束させるマスター曲線を導出した：
   $$t^*_m = 20 \cdot e^{0.58 \cdot N_p \cdot Fr \cdot Ri}$$
   ここで、$N_p$ は動力数、$Fr$はフルード数、$Ri$ はリチャードソン数である。係数 20 は慣性ベースラインを表し、指数項は浮力の安定化効果を捉える。
4. 非単調な挙動: 一定の $Ri$における次元混合時間（$t_m$）を分析したところ、高いリチャードソン数（$Ri \geq 8$）において、$Re$に対する非単調な依存性が観察された。$Re$（すなわちインペラ速度）の増加は初期には乱流を強化して混合時間を短縮するが、さらに増加すると（シミュレーション設定におけるパラメータの結合により）より強い密度成層が生じ、垂直輸送を抑制する。この相互作用により、中間的なレイノルズ数で最小の混合時間が生じる。
5. 流れ構造: 速度場解析により、高い $Ri$ において軽い相がインペラからの動力入力から「遮断」され、垂直構造が主に高密度の液体に現れることが明らかになった。パワースペクトル密度（PSD）分析は、浮力が流れの空間的組織に影響を与える一方で、インペラ駆動のダイナミクスが依然としてエネルギーカスケードを支配していることを示した。

意義と主張
本論文は、完全に乱流領域にある初期に密度成層された系における乱流混合時間のためのコンパクトで予測的な相関を提供することを主張する。動力数、フルード数、およびリチャードソン数の複合的な影響に基づいて導出されたこのスケーリング則は、技術者および応用専門家にとって単純な解析式を提供する。

著者らは、完全に乱流領域における成層系に対する同様の相関が以前に報告されたことはないと述べている。この研究の意義は、高価な実世界テストへの依存を減らすことで混合プロセス設計を合理化する能力にある。提案されたモデルは、浮力効果が顕著な系における混合時間の推定を可能にし、強化された乱流と強化された成層との間の複雑な競合を捉える。著者らは控えめに、将来の研究では、提案されたスケーリングの堅牢性を評価するために、傾斜撹拌機を含む異なるタンク幾何形状およびインペラタイプでこの相関をテストすることを提案している。