CFD中文网

fvc::flux(phi, alpha1, alphaScheme) = fv::convectionScheme<Type>::New ( mesh(), phi, alphaScheme )().flux(phi, alpha1) = phi*interpolate(phi, alpha1) = phi*alpha1f

因此，代码fvc::flux(phi, alpha1, alphaScheme)即为$\nabla\cdot\left(\alpha \mathbf{U}\right)$

https://github.com/OpenFOAM/OpenFOAM-dev/blob/master/src/finiteVolume/finiteVolume/convectionSchemes/gaussConvectionScheme/gaussConvectionScheme.C

https://github.com/OpenFOAM/OpenFOAM-dev/blob/master/src/finiteVolume/finiteVolume/fvc/fvcFluxTemplates.C

@yhdthu 好的，多谢您的回答

sigma是表面张力

需要给真实的压力，可以从单位判别出来，如果是压力/密度的，就不需要，如果是压力的单位，建议给真实的压力。

@nanxuan 请问你解决这个问题了吗？能否分享一下？谢谢

@李东岳 那我该怎么提取速度梯度中的量呢？

gradU.component(tensor::ZX)

这样对吗？

@程迪
rhoE应该是$\frac{\partial \rho e}{\partial t} + \frac{\partial \rho K}{\partial t}$
phiEp应该是$\nabla \cdot (\rho \mathbf{U} e)+\nabla \cdot (\rho \mathbf{U} K)+\nabla\cdot(p\mathbf{U})$
sigmaDotU应该是$\nabla \cdot(\tau \cdot \mathbf{U})$
感觉拆不了。

另外，请问对定常的求解器也有时间步限制吗？如果没有的话，是不是定常的求解器中，比如buoyanySimpleFoam，就不用管隐式还是显式了？

@李东岳 就是考虑粘性耗散项，粘性耗散热也是热源

@李东岳 问题已经解决了，是我在录制脚本的时候进行了多余的不必要操作，再次感谢李东岳老师，您的建议大大的减少了我的工作量。

不写的话默认是yes

@小龙 openfoam后处理forcesIncompressible可以用吗？因为我没有找到这个的详细解释，可以给我说一下这个后处理计算指定patch上压力的原理是什么吗？

flowRatePatch自动适配，phi是质量通量就计算的质量通量，phi是体积通量就计算的体积通量

你好，我也遇到了你一样斗问题，请问你最后是怎么设置成功的呢，是给定压力进出口边界条件，压力固定值进口，固定值出口，速度全部零法向梯度吗

@XIN 请问您解决了吗，我也遇到这个问题了

0_1495693485604_upload-f68f59bc-d32a-4cba-93cb-7d5862ab8826 如图箭头所指的面，想用toposet设置为边界，该用什么命令实现呢？几何文件上传失败了。

@赵一铭 润滑理论的雷诺方程就是NS方程简化了，不考虑厚度方向的压力变化。因为在润滑中动静环之间的液膜很薄，可以认为在厚度方向没有压力变化，是雷诺提出的，就叫雷诺方程了。不知道和那个雷诺是不是一个雷诺。

@赵一铭
你好，请问对于并行计算的案例该如何用gnuplot实时监测残差呢？

您好！我遇到类似问题，准备采用上述的第二种方法即topoSet + createPatch，但不太清楚具体应该如何操作，希望能指点一下，谢谢！@赵一铭 @小龙

@李东岳 好的，谢谢您的帮助