[論文レビュー] Peng: A program for transport properties of low-density binary gas mixtures

Yu Zhai, You Li+2 more|arXiv (Cornell University)|Aug 6, 2022

Phase Equilibria and Thermodynamics被引用数 10

ひとこと要約

Pengは、高精度な分子間ポテンシャルエネルギー曲線（PECs）を用いて、低密度二成分気体混合物の輸送特性（粘性率、熱伝導率、拡散率、熱拡散率）を計算するC++ベースのオブジェクト指向プログラムである。衝突積分とボルツマン方程式のチャップマン–エンスコグ解法を任意の次数で数値的に解くことで、計算リソースの制限のみを除き高精度な結果が得られ、動的共有ライブラリによる拡張性も備えている。

ABSTRACT

The fundamental properties of molecules bridge experiment and theory. Transport properties (diffusion, thermal diffusion, thermal conductivity and viscosity) of binary mixtures are measurable in experiments, and well-defined in theory, but difficult to compute with high accuracy. In addition to high-accuracy inter-molecular potential energy curves (PECs), a reliable and high-order solution program that compute the properties based on the PECs is required. In this work, we present a computer program called Peng that performs the collision integration numerically, and solves the Boltzmann equation in Chapman-Enskog fashion. The program has been devised to perform both parts of the solution procedure to arbitrary order, so that no hard-coded limitation will prevent a user from computing at higher precision, except the amount of RAM and the required computational time. Peng is well-designed in an Object-Oriented Programming (OOP) fashion, which make the program clear and easy to modify. In addition to the end-user oriented program, Peng is also compiled as a dynamic shared library that may readily be extended and embedded in users' programs.

研究の動機と目的

• 理論的分子間ポテンシャルと低密度二成分気体混合物の測定可能な輸送特性の間の橋渡しを図ること。
• PECsから輸送特性を計算するための公開可能で保守可能で拡張可能なソフトウェアの不足に対処すること。
• 衝突積分とチャップマン–エンスコグ輸送特性計算のための高精度でオブジェクト指向的かつ動的に埋め込めるソリューションを提供すること。
• 従来のコードが抱える制限、例えばハードコードされた精度、モノリシックな設計、保守性の低さを克服すること。

提案手法

• プログラムは、正確なPECsを用いて、温度範囲にわたる衝突積分 $\Omega^{(\ell,s)}$ を数値的に計算する。
• ボルツマン方程式をチャップマン–エンスコグ法で任意の次数にまで解き、ハードコードされた精度制限を回避する。
• モジュール性、明確さ、拡張可能性を高めるために、オブジェクト指向設計原則に基づいてC++で実装されている。
• 動的共有ライブラリのコンパイルをサポートしており、他のアプリケーションへの統合が可能である。
• ユーザーは、コードのアーキテクチャを再設計せずに、新しいポテンシャル関数や統合手法を追加できる。
• 運動論的枠組みにおいてテンソル代数と偏微分を処理し、数学的整合性を保証している。

実験結果

リサーチクエスチョン

• RQ1 どのようにして、ab initio の分子間PECsから低密度二成分気体混合物の輸送特性を高精度に計算できるか？
• RQ2 どのような設計原則が、現代のHPC環境における高機能な拡張性と保守可能性を両立させる輸送特性プログラムを可能にするか？
• RQ3 オブジェクト指向C++コードは、従来の運動論的理論ソフトウェアの保守性と再利用性をどの程度向上できるか？
• RQ4 動的ロード可能なライブラリインターフェースは、大規模なシミュレーションワークフローへの輸送特性計算の埋め込みに実際に有効に使えるか？
• RQ5 ハードコードされた精度制限が欠落していることで、輸送特性計算の精度とスケーラビリティにどのような影響が生じるか？

主な発見

• Pengは、利用可能なRAMと計算時間の制限のみを除き、チャップマン–エンスコグ展開の任意の次数での輸送特性計算を可能にする。
• プログラムのオブジェクト指向設計により、コンponentsの明確な分離が実現され、従来のコードよりも変更・拡張が容易である。
• ソフトウェアは動的共有ライブラリとしてコンパイルされ、外部アプリケーションやワークフローへの統合が容易である。
• モノリシック構造、グローバル変数への依存、古くなったプログラミングスタイルといった、一般的な旧来型プログラムの制限を克服している。
• 高精度な数値統合により衝突積分を高精度に計算でき、高品質なPECsに対して正確な結果を保証している。
• 現代のシステムと互換性があり、GNU/Linux上で即座にコンパイル可能であり、以前のツールのポータビリティ問題を解決している。