第153章 普适性积冰生长模型（1 / 2）

“多体耗散粒子动力学……”

祝兰重复了一遍这个名词，同时抬手理了一下头发。

这是她在思考时会做的一个无意识动作。

作为一个表面物理学教授，祝兰自然很早就开始接触分子动力学相关的计算模拟了。

其中自然也包括耗散粒子动力学（dpd）——一种在最近几年才被提出的、应对介观粒子运动的数值模拟方法。

但这个多体耗散粒子动力学又是什么？

dpd中的基本粒子本来就不是分子动力学中的单个原子或者分子，而是代表着一个个分子团的粗粒化粒子。

由于省去了分子动力学中对所有分子进行描述的麻烦，系统的自由度大大减少，计算量也随之呈指数级减少，非常适合介观尺度或者具有一定规模的宏观尺度的研究。

不过这个“一定规模”，在分子动力学领域指的也就是100nm，最多不会超过1μm这个数量级。

如果多体的意思是继续扩大研究的尺度，那么就意味着又要增加需要研究的分子团数量，和建立dpd方法的初衷相悖……

无法解开困惑的祝兰最终决定“不耻下问”——她和林国范都不是航空工业系统内部的人，对于常浩南这个名字自然不会有什么印象。

所以在他们看来，面前这个人应该只是182厂主要负责运8除冰装置改进的年轻工程师而已。

至少目前还是这样。

“好吧，你说的多体耗散粒子动力学和我所熟悉的耗散粒子动力学之间有什么关系么，区别又在哪里？”

实际上这个时候，常浩南也刚刚把瞬间涌入自己脑子里面的知识给整理完，关于mdpd的基本原理，他也就只比坐在旁边的祝兰早了大概两分钟知道。

“最大的不同在于，mdpd对状态方程进行了改造，基于此得到了新的控制方程。”

“dpd中的保守力是一个纯排斥作用的力，其大小随着距离的增加而减小，由这样的力控制的粒子会不断地彼此远离并最终充斥整个计算空间，形成一团密度、温度等各项热力学性质均匀的物质。”

说完，常浩南从旁边拿过纸笔，写下了dpd的状态方程。

p=kt+a^2

之所以这么做并不只是为了让林国范和祝兰更容易看懂，也是为了加深常浩南自己的理解，以及顺便把mdpd的理论原理给整理下来。

“这个状态方程里的流体密度的最高次项为二次项，而用于描述液体内部压力的状态方程需要含有流体密度的三次项，故这种形式的状态方程从本质上无法数值模拟带有自由液气界面的流体系统。”