关于在浓度方程对流项中添加U+sudu(sudu添加无效果）的问题
-
@东岳 我做的测吹方式要比刘仲秋老师（东岳先生您发给我的文章）的模型简单，我是用浓度方程去代替含气率方程（他们在方程形式上有很高的类似度）而且以前有前人用混合单项流模型去模拟的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求解器进行的改动 -
@东岳 @麦迪文
感谢东岳先生和大佬，问题已经解决了，phi1=phi+phi也实现了，浓度场也变化了；
现在想问一下关于声明里有一段代码理解的依然不是很好
这个是我仿照U的声明做的一个Uslip，目的是为了读取0文件夹下的Uslip。
这个是我新生命的phi1场，其目的是为了让我的Uslip从体心插值到面上，这个思路是否正确呢？
这句代码的含义是将速度从体心插值到面上并且初始化（没理解错的话）
ps:我的网格文件实在是有点大，而且处理成5000以下计算效果很不好，就不麻烦东岳先生帮我验证了。感谢东岳先生百忙之中帮我研究我的问题，这些问题很基础，但对于小白来说也很棘手，感激东岳先生，拜谢！
中去的，而且第二个phi想达到一个Unew乘以网格插值出来一个新的面通量，如果压力方程里做phi = 2*phi1的计算的话，会不会影响连续性方程和动量方程呢？我并不想去影响这两个方程
这是我写的phi1声明。
这是我写的浓度方程文件，在括号里加入了phi1 = phi +phi;