类型转化
-
我现在需要表示出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场
-
@李东岳 老师，我想要vol场~
-
@hongjiewang $\frac{\p p}{\p \alpha}$可以写成$\frac{\p p}{\p x}/\frac{\p \alpha}{\p x}$，其中$\frac{\p p}{\p x}$就是$\nabla\cdot p$的x方向分量，这样理论上可以计算出来。但感觉很迂回。
如果要surface场的话会简单一些，如果$\frac{\p p}{\p \alpha}$定义在面上，那就可以考虑这个面相邻的own与nei，把它离散就变成了$\frac{p_{nei}-p_{own}}{\alpha_{nei}-\alpha_{own}}$，这样就可以求出面的$\frac{\p p}{\p \alpha}$。如果进一步需要体的$\frac{\p p}{\p \alpha}$，可以reconstruct一下。
我感觉后者看起来更好一些。你要不要俩种方法都试试
-
@李东岳 嗯嗯好的~所以我上面的采用的其实是老师提到的第一种方法，但是需要把x分量表示出来，这样就不是vector量了~
第二种方法我也尝试一下 谢谢老师~ -
@hongjiewang 你试试 后续反馈 你这个还挺有意思
