@李东岳 在 Simple算法中通量phi的计算 中说:
也就是说在边界处phi是不正确的。不知道是我没有找到phi边界修正,
边界处的
phi
不需要计算,是用户给的。加入速度1m/s,乘以网格边界面积,就是phi
啊,并且不变。
如果是pressureInletVelocity 边界呢? 代码里可以看出pressureInletVelocity 需要用边界上的phi来更新出新的边界速度。如果是这种边界边界上的phi又应该如何确认呢?
@李东岳 在 Simple算法中通量phi的计算 中说:
也就是说在边界处phi是不正确的。不知道是我没有找到phi边界修正,
边界处的
phi
不需要计算,是用户给的。加入速度1m/s,乘以网格边界面积,就是phi
啊,并且不变。
如果是pressureInletVelocity 边界呢? 代码里可以看出pressureInletVelocity 需要用边界上的phi来更新出新的边界速度。如果是这种边界边界上的phi又应该如何确认呢?
@李东岳 这里我看先创建了湍流再创建了MRF区域所以说构造的时候还没有变成相对流量?但是如果phi是存在内存里的话,构造的时候传递了地址,实际上再动量方程中对phi有修正的话,湍流计算的时候也会相应的修正的吧。 我想把phi这个值给输出出来查看一下。李老师知道该怎么做吗?
如图所示的代码中,openfoam构造turbulence的时候用到了phi,想请问这个phi是不是和动量方程的phi是同一个?
cited:http://openfoamwiki.net/index.php/See_the_MRF_development
图为在旋转坐标系内解绝对速度的控制方程,其中 \vec \omega \times \vec U_i 这项称为科氏力,旋转坐标系内独有的力。他的效果如下面这个视频所示https://haokan.baidu.com/v?vid=4227630901810560602&pd=bjh&fr=bjhauthor&type=video
我的理解是可以将上面的式子理解为绝对速度产生组成的动量的输运方程,其中的通量是相对通量(\vec U_R \vec U_I),同时由于动量是一个矢量,因此在旋转坐标系中还受到科氏力的作用。等式坐标是压力的贡献,和由于运输的物理量U_I的扩散贡献。
因此,我在思考是不是诸如湍动能,湍流耗散率这些标量,也可以认为在旋转坐标系内,运输物理量为湍动能,湍流耗散率这些标量,通量为相对速度的通量,但是由于这些量都是标量,因此没有多的项需要添加。另外比如SA模型中的发生项需要计算到voricity的大小,此处的voricity是不是也应该用相对速度?
能不能将MRF理解为有一个新的速度场,该速度场是由绝对速度的速度场减去参考系的牵连速度得到的。该新的速度场运输了单位质量的流体,运输的物理量包括绝对速度的动量、温度、湍动能等。所以如果应用了MRF方法,界面的通量应该为相对速度的通量无论对于那个方程来说(动量,能量,质量,湍流物理量,sclar transport equation )。 请问我理解的对吗
@李东岳 在 SIMPLE(C)解一次動量方程後多次求解壓力方程,壓力方程中計算HbyA的疑問 中说:
pisoFoam
因爲我看到動量預測不是必須的,那我想我能不能某幾步求解一下壓力方程再做一次動量預測,因爲有時候動量方程的增量其實很小,解一次壓力方程對應一次動量方程,是不是有點浪費。
引用自李老師的simple解析。
如果現在每1部做一次動量預測,然後求解20次壓力方程,20次壓力方程的每一步是不是都需要構建新的hybA? 構建HbyA 需要用到An,同時壓力方程的系數需要用到Ap,是否每一步都需要去更新一次Ap與An?
您可以看看苏州舜云工程软件有限公司哦,主要业务就是做基于粒子法的流体多物理场数值仿真软件研发的,就在苏州市相城区,可以关注公众号:舜云多物理场仿真,了解企业的更多信息,简历投递邮箱是lijuan.cai@shoncloud.com