压力梯度中不可导的点如何处理呢？

[李东岳]

你上面贴图里面 ，对x的导数都不是0

面法向梯度指的是外面，$\bfn\cdot\nabla p =0$。边界条件是你人为给定的，不是从方程推出来的。一个出口，给固定值零梯度都可以，是用户自定义参数。但是边界条件要符合物理。给的不好会引起收敛的问题。

[Joker]

@李东岳 在 压力梯度中不可导的点如何处理呢？ 中说：

面法向梯度指的是外面

老师，你这句话的意思是我所有边界的面法向梯度都是指向流域外么？

[李东岳]

是的，边界都指向外面。

[Joker]

@李东岳 好的，谢谢老师！
老师，我还有一个小问题，我在定义的边界面的时候用了右手螺旋定则去写点的顺序，跟那个有关吗？还是说无论我怎么去定义边界面，他的面法向都指向流域外呢？

[李东岳]

无论我怎么去定义边界面，他的面法向都指向流域外呢？

是的。边界的面方向都朝外。

[Joker]

@李东岳 明白了，谢谢老师

[yms00yms]

你好，我也在算振荡流的问题，进口的振荡流基本满足正弦形式，之前进口一直使用的零梯度压力边界条件，发现计算得到的压力值明显不对。最近才了解到进口处压力边界需要实时满足动量方程的形式，类似你图片里写的，请问你最后在openfoam里是怎样实现这样的压力进口边界的？

[Joker]

@yms00yms 这个边界条件需要自己写，类似于你写的振荡流入口的边界条件。

[yms00yms]

@Joker 感谢解答！我目前想了另一种技术路线，如果是用在场里面加体积力的方式来实现振荡流，然后进出口使用cyclic边界，请问你感觉这样做是否能够实现一样的效果呢

[Joker]

@yms00yms 不好意思哈，这我真不好说，那个cyclic边界条件我没用过。