TwoPhaseEulerFoam求解器算例设置

Jack_kerouac

用两相欧拉求解器模拟固液两相流时，phaseProperties文件里面连续相为什么还需要设置粒径？？？？

李东岳

想象一下：如果你的水变成了离散相，就需要设置粒径了。比如你的水变成了雾。

Jack_kerouac

@李东岳也就说如果固液两相流不涉及的相变，那么模拟计算过程中，这个连续相粒径参数其实是被忽略的，可以这样理解吗？东岳老师

Jack_kerouac

@李东岳东岳老师，两相欧拉自带流化床算例中，仅在thermophysicalproperties.particles里面提到密度设置，thermophysicalproperties.air里面并没有出现密度相关的设置项，这是为何？若固液两相流不涉及传热计算时，不置密度的话没法体现分层啊。

李东岳

一直在用twoPhaseEulerFoam求解器模拟固液两相流动（固体颗粒为水合物颗粒，入口体积分数为0.2，填充率为0.53）特性，针对求解器做了如下修改:
1.注释掉温度方程，因为所研究的固液流动不涉及传热
2.注释掉升力，虚拟质量力，壁面润滑力，仅保留曳力和湍流扩散力
由于想先从层流模拟入手，所以一直调试层流算例（对流化床层流算例做针对性的修改），调试算例近一月，仍无进展，一筹莫展，每当算例计算到中途得时候，终端提示计算曳力出问题，如下所示。请教一下东岳老师，能否帮我看看哪里出问题了。

这个是你给我发的？最近我在搬家，焦头烂额，等我安顿下来给你具体回复，实在抱歉。

PS divive错误一般是因为你除了0导致的，分母不能为0.

Jack_kerouac

@李东岳是的我邮件发给你的谢谢东岳老师。

李东岳

每当算例计算到中途得时候，终端提示计算曳力出问题

你用的什么曳力模型，试过其他曳力模型么？你的最小ahlph的值是多少？

Jack_kerouac

@李东岳层流计算中，曳力用的是GidaspowErgunWenYu模型，其余模型暂时没试过，alpha.particles设置中，入口设为0.2，内部场为0.1，出口和壁面为zeroGradient.
相关文件设置如图，依次为：颗粒速度文件，液体速度文件（颗粒和液体的密度已经在热物理特性文件中设置了），P_rgh以及颗粒相分数文件。两相的湍流特性文件中我均设置为laminar。

Jack_kerouac

@李东岳东岳老师，两相欧拉自带的流化床算例中没有在热物理特性文件中给定颗粒相的粘度为0，在固液两相计算过程中，给定颗粒相粘度为0会导致曳力计算出问题？

李东岳

曳力出错的根本问题在于在计算曳力的时候需要调用相分数、粒径等参数。如果某个参数在分母是0，曳力就算不了了。我怀疑你的alpha有问题，你测试你那个最小的ahlph值是多少？

李东岳

或许 @xpqiu 会给你一些关于颗粒相粘度的提示。

因为最近在做欧拉拉格朗日的研究，获取了一些气固模拟的知识：

首先是DPM模拟，直接跟踪颗粒，但是需要考虑颗粒之间的碰撞作用，通常选用硬球模型或者软球模型；
由于DPM计算要求较高（粒子数量较多），因此可以将气相和固相都看成流体，这就是TFM，也就是OpenFOAM里面的用于模拟气固的twoPhaseEulerFoam。在TFM中，颗粒的碰撞采用kinetic theory of granular flow来模化，因此在控制方程中出现了$\tau$这一项，也即solid stress。

李东岳

http://pubs.acs.org/doi/pdf/10.1021/acs.iecr.7b01862

这个文章前面的讨论很易懂，类似的基本理论介绍在CFD-DEM，DPM，MPPIC的文章中都有，可以看一下。

LiuZhuang

东岳老师您的邮箱多少？关于twoPhaseEulerFoam的问题我有几点想跟您请教

CFD中文网

TwoPhaseEulerFoam求解器算例设置