类型转化

[hongjiewang]

我现在需要表示出dpc/da（volScalarField），一种方法是在界面处定义为无穷大的数（可以实现）；另一种方法是我想转化为和x的关系，但是涉及到dx，就会变成volVector（当然可以用mag把vector转化为scalar，再相除是可以得到scalar类型的，但是好像物理意义不对了,在运行算例的时候会报错~），想问一下应该怎么处理~
volScalarField dpcdx = mag(fvc::grad(pc_));
volScalarField dSedx = mag(fvc::grad(Se_));
volScalarField dpcdS_= dpcdx/dSedx;

[李东岳]

就是你有一个变量p，一个变量a，想求dp/da? 你要vol场还是surface场

[hongjiewang]

@李东岳 老师，我想要vol场~

[李东岳]

@hongjiewang $\frac{\mathrm{\partial }p}{\mathrm{\partial }\alpha }$可以写成$\frac{\mathrm{\partial }p}{\mathrm{\partial }x}/\frac{\mathrm{\partial }\alpha }{\mathrm{\partial }x}$，其中$\frac{\mathrm{\partial }p}{\mathrm{\partial }x}$就是$\mathrm{\nabla }\cdot p$的x方向分量，这样理论上可以计算出来。但感觉很迂回。

如果要surface场的话会简单一些，如果$\frac{\mathrm{\partial }p}{\mathrm{\partial }\alpha }$定义在面上，那就可以考虑这个面相邻的own与nei，把它离散就变成了$\frac{p_{nei}-p_{own}}{{\alpha }_{nei}-{\alpha }_{own}}$，这样就可以求出面的$\frac{\mathrm{\partial }p}{\mathrm{\partial }\alpha }$。如果进一步需要体的$\frac{\mathrm{\partial }p}{\mathrm{\partial }\alpha }$，可以reconstruct一下。

我感觉后者看起来更好一些。你要不要俩种方法都试试

[hongjiewang]

@李东岳 嗯嗯好的~所以我上面的采用的其实是老师提到的第一种方法，但是需要把x分量表示出来，这样就不是vector量了~
第二种方法我也尝试一下 谢谢老师~

[李东岳]

@hongjiewang 你试试 后续反馈 你这个还挺有意思