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

[李东岳]

dp/dx=du/dt怎么得出来的

[Joker]

@李东岳 刚刚写的有点问题，dp/dx=-du/dt，

[Joker]

@李东岳 老师，我还想再请教一个问题，我如何确定物面法向梯度指向什么地方呢？是指向流域内还是流域外呢？跟我定义边界的时候定义有关么？

[李东岳]

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

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

[Joker]

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

面法向梯度指的是外面

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

[李东岳]

是的，边界都指向外面。

[Joker]

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

[李东岳]

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

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

[Joker]

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

[yms00yms]

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

[Joker]

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