Your browser does not seem to support JavaScript. As a result, your viewing experience will be diminished, and you have been placed in read-only mode.
Please download a browser that supports JavaScript, or enable it if it's disabled (i.e. NoScript).
我们在做振荡流的数值模拟的时候，面法向压力梯度等于速度对时间的导数，但是如果速度对时间的导数存在不可导的点时，这种情况该怎么办呢？
我看你速度是sin(t)函数，这个一直可导，为什么会不可导？
@李东岳 老师，你好，就是我想表达得意思就是，如果我把um换成sin（wt）的绝对值，那他就会存在不可导的点，这种情况下该怎么办呢？
dp/dx=du/dt怎么得出来的
@李东岳 刚刚写的有点问题，dp/dx=-du/dt，
@李东岳 老师，我还想再请教一个问题，我如何确定物面法向梯度指向什么地方呢？是指向流域内还是流域外呢？跟我定义边界的时候定义有关么？
你上面贴图里面 ，对x的导数都不是0
面法向梯度指的是外面，$\bfn\cdot\nabla p =0$。边界条件是你人为给定的，不是从方程推出来的。一个出口，给固定值零梯度都可以，是用户自定义参数。但是边界条件要符合物理。给的不好会引起收敛的问题。
@李东岳 在 压力梯度中不可导的点如何处理呢？ 中说：
面法向梯度指的是外面
老师，你这句话的意思是我所有边界的面法向梯度都是指向流域外么？
是的，边界都指向外面。
@李东岳 好的，谢谢老师！ 老师，我还有一个小问题，我在定义的边界面的时候用了右手螺旋定则去写点的顺序，跟那个有关吗？还是说无论我怎么去定义边界面，他的面法向都指向流域外呢？
无论我怎么去定义边界面，他的面法向都指向流域外呢？
是的。边界的面方向都朝外。
@李东岳 明白了，谢谢老师
你好，我也在算振荡流的问题，进口的振荡流基本满足正弦形式，之前进口一直使用的零梯度压力边界条件，发现计算得到的压力值明显不对。最近才了解到进口处压力边界需要实时满足动量方程的形式，类似你图片里写的，请问你最后在openfoam里是怎样实现这样的压力进口边界的？
@yms00yms 这个边界条件需要自己写，类似于你写的振荡流入口的边界条件。
@Joker 感谢解答！我目前想了另一种技术路线，如果是用在场里面加体积力的方式来实现振荡流，然后进出口使用cyclic边界，请问你感觉这样做是否能够实现一样的效果呢
@yms00yms 不好意思哈，这我真不好说，那个cyclic边界条件我没用过。
