[论文解读] 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 是一个基于 C++ 的面向对象程序，利用高精度的分子间势能曲线（PECs）计算低密度二元气体混合物的输运性质（粘度、导热系数、扩散系数和热扩散系数）。它通过数值求解碰撞积分和玻尔兹曼方程的 Chapman–Enskog 解法至任意阶次，实现仅受计算资源限制的高精度结果，并可通过动态共享库实现可扩展性。

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 计算输运性质的问题。
• 提供一种高精度、面向对象且可动态嵌入的碰撞积分与 Chapman–Enskog 输运性质计算解决方案。
• 克服传统代码的局限性，包括硬编码精度、单体式设计和可维护性差等问题。

提出的方法

• 该程序使用精确的 PECs，在一系列温度范围内数值计算碰撞积分 $\Omega^{(\ell,s)}$。
• 通过 Chapman–Enskog 方法求解玻尔兹曼方程至任意阶次，避免硬编码精度限制。
• 软件使用 C++ 实现，采用面向对象设计原则，以提高模块化、清晰度和可扩展性。
• 支持动态共享库编译，便于集成到其他应用程序中。
• 该框架允许用户在不重构整体架构的前提下，通过添加新的势能函数或积分方法来扩展代码。
• 程序在统计力学框架内处理张量代数和偏导数，确保数学一致性。

实验结果

研究问题

• RQ1 如何从从头计算得到的分子间 PECs 高精度计算低密度二元气体混合物的输运性质？
• RQ2 哪些设计原则能够使输运性质程序在现代高性能计算环境中实现高度可扩展性和可维护性？
• RQ3 面向对象的 C++ 代码在多大程度上能够提升传统统计力学软件的可维护性和可重用性？
• RQ4 动态加载库接口能否有效用于将输运性质计算嵌入到更大的模拟工作流中？
• RQ5 去除硬编码精度限制对输运性质计算的精度和可扩展性有何影响？

主要发现

• Peng 允许在 Chapman–Enskog 展开中任意阶次计算输运性质，其精度仅受可用内存和计算时间的限制。
• 该程序的面向对象设计实现了组件间的清晰分离，使其比传统代码更易于修改和扩展。
• 该软件被编译为动态共享库，便于集成到外部应用程序和工作流中。
• 它克服了旧程序的常见局限，如单体式结构、依赖全局变量以及过时的编程风格。
• 该框架支持对碰撞积分进行高精度数值积分，确保高质量 PECs 下的计算准确性。
• 该程序与现代系统兼容，可在 GNU/Linux 系统上开箱即用，解决了先前工具的可移植性问题。