A Volume of Fluid Immersed Boundary Method for Industrial Polymer Mixing

本論文は、OpenFOAM に実装された新規ブロック結合型 VOF 浸没境界 (BC-VOF-IB) ソルバを提示するものであり、このソルバは高粘度比に起因する数値的不安定性を克服し、部分的に充填された産業用押出機における自由表面を有するポリマー混合を正確にシミュレートするものである。

要約

巨大なタンクの中で、空気と混ざり合う濃厚で粘り気のある蜂蜜を、回転する機械内で混ぜ合わせようとしていると想像してください。これは、ピレリのような企業がタイヤ、医療機器、自動車部品などを製造するために、溶融したプラスチックに添加剤を混合する工業的なポリマー混合において、本質的に起こっていることです。目標は、最終製品が強力かつ均一になるよう、すべてを完璧に混合することです。

しかし、このプロセスをコンピュータ上でシミュレーションすることは、数学者やエンジニアにとって悪夢のようなものです。その理由と、この論文がどのように単純なアナロジーを用いてそれを解決したかを以下に示します。

問題：「濃厚な蜂蜜対薄い空気」の格闘

これらの機械内には、非常に異なる二つの流体が存在します。

1. ポリマー溶融物：極めて濃厚で、粘り気があり、動きが遅い（冷えた蜂蜜のようなもの）。
2. 空気：非常に薄く、速く動く。

回転するスクリューを備えた機械内で、これら二つがどのように相互作用するかをシミュレーションしようとすると、標準的なコンピュータプログラムは混乱します。それは、同じトラック上でカタツムリとレーシングカーの動きを、同じルールセットを使って計算しようとするようなものです。コンピュータは「カタツムリ」（濃厚なプラスチック）が動きすぎないようにするために、極めて微小なステップを刻もうとします。その結果、シミュレーションは信じられないほど遅くなり、実際の混合時間の数秒分を完了するのに数日かかることもあります。

さらに、これらの機械には、固定された容器内で動く複雑な回転部品（スクリュー）があります。伝統的にこれをシミュレーションするには、回転するスクリューに完璧に巻き付くように、デジタルメッシュ（小さな箱のグリッド）を構築する必要があります。スクリューが回転するにつれて、このグリッドは常に再形成されなければなりません。それは、マラソンを走っている人が着ているセーターを編み直すようなものです。それは乱雑で困難であり、エラーを起こしやすいものです。

解決策：新しい「スマートグリッド」と「チームアプローチ」

この論文の著者らは、OpenFOAM というソフトウェアを使用して、これらのシミュレーションを実行する新しい方法を開発しました。彼らは、二つの強力な技術を組み合わせました。

1. 埋め込み境界法（「ゴーストウォール」のトリック）
回転するスクリューに合わせてグリッドを再形成する代わりに、彼らはグリッドを固定された剛体（氷の塊のようなもの）のままに保ちました。そして、コンピュータに「ねえ、この氷の塊の中に回転するスクリューがあるよ」と指示しました。

• アナロジー：底に固定されたタイルのグリッドを持つプールを想像してください。泳ぎ手に合わせてタイルを動かすのではなく、単に水に「泳ぎ手を通過するな」と指示するだけです。コンピュータは数学を用いてスクリューの周りに「ゴーストウォール」を作成し、流体がその周りを流れるように強制します。これにより、グリッドを再構築する必要なく、複雑で移動する形状を処理することが容易になります。

2. 流体体積法（VOF）（「追跡ペイント」のトリック）
濃厚なプラスチックと空気の境界がどこにあるかを把握するために、彼らはセルを埋める「ペイント」を使用します。

• アナロジー：コンピュータのグリッドを 3 次元のチェッカーボードだと想像してください。いくつかのマス目は 100% プラスチック、いくつかは 100% 空気、そしていくつかは混合状態です。コンピュータは、各マスに含まれる「プラスチックペイント」の量を追跡することで、液体の表面を把握します。

3. ブロック結合スキーム（「チームの集まり」）
これが最も重要な画期的な進歩です。標準的なシミュレーションでは、コンピュータは流体の速度を X、Y、Z 方向について、三人が順番に話すように、一つずつ順に解きます。流体が非常に濃厚な場合（ポリマーのように）、この「順番に話す」アプローチは、濃厚な流体がすべての方向を密に結合しているため、シミュレーションがクラッシュしたり、極端に遅くなったりします。

著者らはこれをブロック結合アプローチに変更しました。

• アナロジー：三人が順番に話すのではなく、全員が集まって同時に問題を解決します。すべての方向の動きを、一つの巨大で相互接続されたチームとして扱うことで、コンピュータは濃厚なプラスチックと薄い空気の巨大な差に立ち向かい、行き詰まることなく処理できます。

結果：数時間から数分へ

チームは、二つのシナリオで新しい方法をテストしました。

1. ドッグボーン形状のチャネル：濃厚なプラスチックが狭くねじれたチャネルに注入されるテストケース。

   • 旧来の方法：標準的なコンピュータプログラムはクラッシュするか、微小なステップを強いられたため、数秒のシミュレーションに7 時間を要しました。
   • 新しい方法：彼らの新しい「チームの集まり」法は、同じ作業をわずか16 分で完了し、プラスチックが極めて濃厚になった場合でもクラッシュしませんでした。

2. 実際の工業用機械：彼らは、プラスチックペレット製造に使用される機械である、実際の単スクリューおよび二重スクリュー押出機をシミュレーションしました。

   • 彼らは、プラスチックが機械をどのように充填し、圧力がどのように上昇し、空気がどのように押し出されるかを成功裏に示しました。
   • 彼らは、彼らの「ゴーストウォール」法が、グリッドを再形成する古く困難な方法と同じくらいよく機能し、はるかに高速で設定が容易であることを証明しました。

次は何が起きるか？

この論文は、これが業界にとって大きな前進であると結論付けています。これは、学術的な数学と実際の工場ニーズの間のギャップを埋めます。ただし、著者らは、現在のモデルは温度が一定（等温）であると仮定していることに注意を促しています。実際には、プラスチックの混合は熱を発生させ、それがプラスチックの粘度を変化させます。温度効果やより複雑な「伸縮性」のあるプラスチックの挙動を追加することが、今後の研究の次のステップです。

要約すると：彼らは、回転する機械内で濃厚なプラスチックが空気と混合する様子をコンピュータ上で監視するための、より高速で安定した方法を開発しました。これにより、かつて数時間を要していたプロセスが数分で完了するようになり、スクリューが回転するたびにデジタル世界を再構築する必要がなくなりました。

技術概要

技術概要：産業用ポリマー混合のための体積流体法と浸没境界法

問題提起
ポリマー混合プロセスの数値シミュレーションは、非ニュートンレオロジー、多相自由表面ダイナミクス、および移動する複雑な幾何形状の複雑な相互作用により、重大な課題に直面しています。単軸スクリュー（SSE）や二軸スクリューエクストルーダー（TSE）などの産業用ミキサーは、しばしば高粘性のポリマー溶融体が空気と共存する部分充填条件下で運転されます。これらのシナリオをシミュレートするには、回転する固体境界を処理しながら、相間の界面を捉える必要があります。

標準的な計算流体力学（CFD）アプローチは、この文脈において 2 つの主要な障壁に直面します：

1. 幾何学的複雑性：複雑な幾何形状を持つ回転スクリュー用の体適合（コンフォーミング）メッシュを生成することは、計算コストが高く、しばしばメッシュ歪みを引き起こす可能性のある動的メッシュ技術を必要とします。
2. 数値的不安定性：ポリマー溶融体と空気との間の極端な粘性比（しばしば $10^5$ から $10^8$ を超える）は、標準的な分離型体積流体（VOF）ソルバーにおいて深刻な数値的不安定性を引き起こします。具体的には、粘性拡散項における転置速度勾配の陽的处理は、厳しい時間ステップ制約（$\Delta t \propto h^2$）を課し、高粘性流れのシミュレーションを計算的に不可能なものにしたり、不安定にしたりします。

手法
本研究は、OpenFOAM 有限体積ライブラリ内で実装された、界面捕捉法である体積流体（VOF）アプローチと非適合浸没境界（IB）法を統合した新しい数値フレームワークを提案します。

• 浸没境界（IB）アプローチ：この手法は、固体境界（例：回転スクリュー）を三角化された STL 表面として表現する固定背景グリッドを利用します。重み付き最小二乗法（WLS）補間演算子を用いて、浸没境界によって交差するセルに対して境界条件（速度に対するディリクレ条件、体積分率に対するノイマン条件）を課します。これにより、適合メッシュ生成の必要性が排除され、複雑な回転幾何形状に沿った移動接触線の処理が可能になります。
• VOF 定式化：界面は、相体積分率（$\alpha$）の輸送方程式を解くことで追跡されます。この手法は、界面の鮮明さと有界性を維持するために、数値拡散を相殺する圧縮速度項を含む MULES（多次元汎用リミッターと陽的解法）アルゴリズムを採用しています。
• ブロック結合（BC）スキーム：高い粘性比に起因する安定性の問題に対処するため、著者らは運動量方程式に対してブロック結合スキームを導入します。速度成分を逐次的に扱う標準的な分離型ソルバーとは異なり、BC アプローチは完全な粘性拡散項 $\nabla \cdot (\mu(\nabla \mathbf{u} + \nabla^\top \mathbf{u}))$ を陰的に扱います。これには、単一の線形システム内で 3 つの速度成分を結合することが含まれます。接線勾配項の離散化には、弾性応用から適応された [11]、拡張されたステンシル（面隣接点だけでなく点隣接点を含む）上での WLS 補間を組み合わせた有限面積（FA）法が利用されます。
• 解法アルゴリズム：ソルバーは、圧力 - 速度結合を処理するために修正された PIMPLE アルゴリズム（SIMPLE と PISO の組み合わせ）を採用しています。IB 制約は代数システムに統合され、浸没表面上の境界条件を満たしつつ質量保存を維持するために、行列係数とソース項が修正されます。

主な貢献

1. BC-VOF-IB ソルバーの開発：本論文は、非適合 IB 法とブロック結合 VOF アプローチを組み合わせる UCInterIbFoam ソルバーを導入します。これは、二相流向けに OpenFOAM-10 フレームワーク内で実装された、そのような結合手法の最初の実装です。
2. ブロック結合による安定性向上：この研究は、粘性拡散の完全な陰的処理が時間ステップの安定性制約を大幅に緩和することを示しています。これにより、クーラン数違反によるシミュレーションのクラッシュなしに、粘性比が $10^8$ までの高粘性ポリマー流れのシミュレーションが可能になります。
3. 移動接触線の処理：このフレームワークは、二相界面が浸没境界と交差する移動接触線における幾何学的およびアルゴリズム的な困難を成功裡に解決しており、部分充填混合装置にとって重要な要件を満たしています。

数値結果
提案された手法は、2 つのカテゴリーの問題で検証およびテストされました。

• ベンチマークケース（ドッグボーン型射出成形）：

  • 高粘性パワー則流体の射出ケースにおいて、分離型（interFoam）、ブロック結合適合型（UCInterFoam）、およびブロック結合非適合型（UCInterIbFoam）のソルバー間での比較が行われました。
  • 標準的な分離型ソルバーは、厳しい時間ステップ制約（$\Delta t \approx 10^{-6}$ 秒）により、合理的な時間枠内で安定した結果を生成できませんでした。
  • ブロック結合ソルバーは、最大 $10^{-3}$ 秒までの時間ステップで安定した解を達成し、計算時間を 7 時間から 16 分（適合型）および 58 分（非適合型）に短縮しました。
  • 非適合 IB 解は、適合解と定性的な一致を示しましたが、幾何学的制限付近の粘性ピークにおいて、グリッドの非適合性に起因するわずかな不一致が見られました。

• 産業応用（SSE および TSE）：

  • 単純化された単軸スクリューエクストルーダー（SSE）と現実的な二軸スクリューエクストルーダー（TSE）のシミュレーションが、完全充填および部分充填（スターブフィーディング）条件の両方で行われました。
  • ソルバーは、自由表面の進化、充填率、および圧力場の捕捉に成功しました。
  • 結果は物理的な期待と一致しました：完全充填シナリオ（$Q_{in} > Q_0$）は入口から出口への負の圧力勾配を示し、部分充填シナリオ（$Q_{in} < Q_0$）は出口付近で圧力上昇を示しました。
  • IB 法は、動的メッシュ再生を必要とすることなく、TSE の複雑な噛み合い幾何形状を処理する上で効果的であることが証明されました。

意義と主張
本論文は、提案された BC-VOF-IB フレームワークが、学術的な CFD 研究と産業用ポリマー加工の要求との間のギャップを埋めるための有意義な一歩を表すと主張しています。高い粘性比に関連する数値的不安定性を克服し、回転幾何形状のための複雑な体適合メッシュ生成の必要性を排除することで、この手法は、産業的に関連する部分充填混合装置における速度場および圧力場の物理的に整合的な予測を可能にします。

著者らは、現在の仕事が等温プロセスおよび一般化ニュートンレオロジーに限定されていることを認めています。彼らは、産業用ポリマー混合の物理を完全に捉えるためには、将来の開発に熱効果（エネルギー方程式）および粘弾性構成則モデルを含める必要があると述べています。しかし、現在のフレームワークはこれらの拡張のための厳密な基盤を確立し、複雑な移動幾何形状における自由表面流れのシミュレーションのための堅牢なツールを提供しています。