关于在浓度方程对流项中添加U+sudu(sudu添加无效果）的问题
-
东岳老师及各位大神们早上好！
本人的课题遇到了一个新的瓶颈，在此发表希望大伙给与指导与点拨。
我想在浓度方程对流项里引入一个速度，并且通过简单的学习将原速度U和新速度Unew引入到phi中，将phi定义为phi1（编译成功），但是我将phi1放入浓度方程的对流项（动量方程与压力泊松方程的对流项仍保持为phi）中时，发现浓度的扩散并没有Unew的效果。
在这里我想问的是，我这样做是不是也需要考虑将动量方程和压力泊松方程的phi也改成phi1?并且对浓度方程里新的phi1应该怎样进行调整（是否需要在动量方程和压力泊松方程里也同样调整phi1)？才能使我加入的phi1能够对浓度的扩散有一个流动的效果？
下面我将改动的求解器代码发到这里，希望大神们能帮我搞清楚这个问题ps:一直对压力方程理解的不是很好。希望大神们给与指导，这些文件我是基于buoyantBoussinesqSimpleFoam的基础上加入浓度方程进行修改的。
-
@东岳 我做的测吹方式要比刘仲秋老师（东岳先生您发给我的文章）的模型简单，我是用浓度方程去代替含气率方程（他们在方程形式上有很高的类似度）而且以前有前人用混合单项流模型去模拟的E-E模型，现在我已经做好了滑移速度在模型中的分布文件Unew，和导师探讨后认为Unew文件应写入到我之前写的浓度方程的对流项中去。然后我发现pEqn.H中关于phi 进行了一系列眼花缭乱的操作，我进行了简单的复制后发现wmake报错（我的想法是按照pEqn.H简单的复写一个phi1,给phi1变成U+Unew的形式），这种操作对于我来说十分吃力，我弄不清楚如果添加这个新的phi1应该对除了createFields和浓度方程之外还有哪些文件进行改动（比如说或者pEqn.H)说实话CFD理论理解十分浅显，根本看不懂Weller的操作，希望东岳大神帮我指点一下迷津
ps:我将混合单项流模型的论文里含气率守恒公式和我做的简单的pEqn.H粘贴在下方。
这个是湍流nu
pEqn.tar 这个是我做简单改动后的pEqn.H
标蓝处我做的简单添加，在文件的第48行
我用的OF版本是2.4.0的，基于buoyantBousssinesqSimpleFoam求解器进行的改动
中去的，而且第二个phi想达到一个Unew乘以网格插值出来一个新的面通量，如果压力方程里做phi = 2*phi1的计算的话，会不会影响连续性方程和动量方程呢？我并不想去影响这两个方程
这是我写的phi1声明。
这是我写的浓度方程文件，在括号里加入了phi1 = phi +phi;