@东岳 成功解决了，在头文件里加了些东西，好有成就感 。准备复盘，写个教程
开心

今年九月份博士入学，研究生做了2年半的OpenFOAM，做的课题不难，在算法上对求解器也没有变动。博士阶段，自身非常想往CFD这方面靠一靠。感觉自己现在就属于那种在OpenFOAM门口站的那种人（不知道有没有人跟我有一样的感觉，就是OF的水平没多高，就是会调调参数，在控制方程里面做点小改动什么的，都能做到，但是再就不行了），十分想想听听老师前辈们的故事 目的就是能打破我目前的状态，正在看那本 The Finite Volume Method in Computational Fluid Dynamics,那本书反正已经拿起来放下好几次了，这回心里立了flag了，我就是啥也没看懂，我也一页一页把那本书翻一遍。

@Calf-Z-DNS 之前硕士的时候在做示踪剂在水中传输过程中，密度随浓度变化的问题，所处理的问题和你遇到的差不多，我的压力也是收敛的不太好，大概在0.01，然后硕士毕业时间紧迫问题就放下了，看见有人也在做这个问题，突然感觉我以前做的工作还挺有意义的（以前觉得没什么意义 ），现在突然想捡起来做一做了

解决了，把1换成1.0就好了，感谢 @麦迪文 师兄的帮助。

在使用OpenFOAM编译公式的时候，发现了个问题：pow(N, 1/2)与pow(N,0.5)输出的结果不一样，OpenFOAM编译公式时幂函数是使用pow()函数吗？我翻看了自带的一些湍流模型中发现了pow4()或者pow3()这种写法，难道是OpenFOAM的幂函数编写规则和C++不一样？下面上图，图中分别是pow(N, 1/2)与pow(N,0.5)的结果，明显量纲不一样了。

@东岳 光看流场的矢量图可能确实看着还凑合对称，但是假如我要在出口得到RTD曲线（停留时间分布曲线），就会得到对称出口明显不同的RTD曲线？

对于一些网格划分较复杂的对称结构来说，左右网格不对称是在所难免的吧，比如在网格划分的时候移动一些线。但是这样会严重影响数值结果的对称性？这种情况怎么处理呢？只画对称结构的一半进行计算？

在reactingtwophaseEulerFoam的默认算例中气相与液相的湍流模型均设置为mixtureKEpsilon，在东岳流体的验证算例讲解中气相的湍流模型设置为层流，液相为kEpsilon湍流模型，我自己算了个钢包底吹的算例，气相与液相都设置为kEpsilon湍流模型算例一步也算不下去。我看了一下文献感觉大家说明的都不是很详细啊（可能我看的文献偏应用），另外附图按照气相为层流模型，液相为湍流模型计算结果如下：



本来以为就是画个模型套用下算例就行，真正做起来还挺费劲，气相分数那个图为啥细细的一根而且上部的气相分数的分布还断断续续的呀，气相速度下部靠近入口处的分布还不是很规则

@东岳 ok

使用reactingTwoPhaseEulerFoam进行计算时，在设置算例时，液相直径是啥意思呀？与气泡直径有什么关系吗？不同工况下我该凭借什么原则进行设置呢？

按照twoPhaseEulerFoam中的EEqns.H写出能量方程如下（算例文件中能量形式是sensibleInternalEnergy）大神们先看看这么写对吗？

fvScalarMatrix E1Eqn
(
    fvm::ddt(alpha1, rho1, he1) + fvm::div(alphaRhoPhi1, he1)
  - fvm::Sp(contErr1, he1)

  + fvc::ddt(alpha1, rho1, K1) + fvc::div(alphaRhoPhi1, K1)
  - contErr1*K1
  + (
        he1.name() == thermo1.phasePropertyName("e")
      ? fvc::div(fvc::absolute(alphaPhi1, alpha1, U1), p)
      + p*fvc::ddt(alpha1)
      : -alpha1*dpdt
    )

  - fvm::laplacian
    (
        fvc::interpolate(alpha1)
       *fvc::interpolate(thermo1.alphaEff(phase1.turbulence().mut())),
        he1
    )
);

E1Eqn.relax();

E1Eqn -=
(
    Kh*(thermo2.T() - thermo1.T())
  + Kh*he1/Cpv1
  - fvm::Sp(Kh/Cpv1, he1)
  + alpha1*rho1*(U1&g)
  + fvOptions(alpha1, rho1, he1)
);

那如果不考虑传热的话，只把导热系数设为0和温度设为一样，这样做对吗？
我为了不考虑传热，是直接把twoPhaseEulerFoam.C中的EEqns.H给注释掉了，我这种做法对吗？

@win 稳态流场的求解的话，为啥不用buoyantBoussinesqSimpleFoam呢？人家重力都给你植入好了，你想前面的系数是多少就是多少。

@麦迪文 可以可以，小师兄进度很快，我们老师给我博士课题分的要用这个求解器，我最近也刚刚开始做这个，找到小伙伴了