关于在浓度方程对流项中添加U+sudu(sudu添加无效果）的问题

[LiuZhuang]

@东岳 果然是我的phi1写的方式不准确，我在createFields.H中声明了一个phi1，如下图
我的想法是创立一个新的phi1让phi1=2倍的phi，不知道这种写法错误在哪，做出来的phi和phi1两个文件值相同。
如何创立一个新的phi1呢？除了在createFields里填写还需要在哪注意呢？ 我那个phi+phi的写法应该错误很大吧，应该怎么去实现phi1=phi+phi的感觉呢？

[李东岳]

我觉得你应该在压力方程之后每一步都进行phi = 2*phi1的计算，你那个只是声明

[LiuZhuang]

@东岳 这个新的phi1我是想放在浓度方程 中去的，而且第二个phi想达到一个Unew乘以网格插值出来一个新的面通量，如果压力方程里做phi = 2*phi1的计算的话，会不会影响连续性方程和动量方程呢？我并不想去影响这两个方程

[LiuZhuang]

@东岳 或者说我想在浓度方程的对流项中将phi改为phi1，phi1的表达式为phi+phi(第二个phi想改动成自己写的速度场文件(即Unew)&mesh.Sf)

[麦迪文]

小藤， 没变化是因为得重新赋值 ，在每个C 方程 前面 来一句 这个 phi1 = phi + sudu

[麦迪文]

不影响其他方程 跟其他方程 没耦合

[LiuZhuang]

@东岳 确实在方程中添加了phi1=phi+phi,也影响了浓度场的分布，现在我遇到的问题是phi 和 phi1文件里的值还是一摸一样的，正常phi1的值也应该和phi 有所区别，我如何才能让phi和phi1的值有区别呢？问题是否出在声明文件中呢？ 这是我写的phi1声明。
这是我写的浓度方程文件，在括号里加入了phi1 = phi +phi;

[李东岳]

把你的代码+算例处理干净（代码wclean，算例网格5千以下）上传一下我试试

[LiuZhuang]

@东岳 @麦迪文
感谢东岳先生和大佬，问题已经解决了，phi1=phi+phi也实现了，浓度场也变化了；
现在想问一下关于声明里有一段代码理解的依然不是很好

这个是我仿照U的声明做的一个Uslip，目的是为了读取0文件夹下的Uslip。

这个是我新生命的phi1场，其目的是为了让我的Uslip从体心插值到面上，这个思路是否正确呢？

这句代码的含义是将速度从体心插值到面上并且初始化（没理解错的话）
ps:我的网格文件实在是有点大，而且处理成5000以下计算效果很不好，就不麻烦东岳先生帮我验证了。

感谢东岳先生百忙之中帮我研究我的问题，这些问题很基础，但对于小白来说也很棘手，感激东岳先生，拜谢！

[李东岳]

目前来看你的思路是正确的