最后找到了解决办法,是边界条件的原因。本来上下圆柱连接的上下两个边界用的symmetry边界条件,改成cylic边界之后就能算了。。。具体原因还没研究出来
Fu
@Fu
Best posts made by Fu
Latest posts made by Fu
-
RE: 三维sixDoFRigidBodyMotion算例发散
-
RE: 用pimpleFOAM做圆柱的涡激振动遇到收敛性问题
@赵海盛 看起来没有问题的。。如果是自激振荡就是caculated,如果是强迫就要oscillatingdisplacement了。。。看边界条件是一个三维的圆形流场吗?。。
-
RE: 200万网格管道流qDNS并行算力测试
CPU型号: Intel(R) Xeon(R) Gold 6248R 系统:ubuntu 18.04 版本:openfoam-8 48 86.08 32 93.05 16 125.13 8 153.34 4 287.14 2 524.41 1 1013.26
正常了。。。
-
RE: 200万网格管道流qDNS并行算力测试
CPU型号: Intel(R) Xeon(R) Gold 6248R 系统:ubuntu 18.04 版本:openfoam-8 48 98.12 32 96.14 16 128.01 8 155.18 4 304.95 2 542.82 1 1068.49
没关超线程,物理核48 感觉不合理啊哈哈 明天再弄一次
-
三维sixDoFRigidBodyMotion算例发散
尝试使用sixDoFRigidBodyMotion计算一个弹簧连接的三维圆柱在来流下的振动情况,但是结果发散。
在算到发散的前几步之前,库朗数开始急剧增加。
Courant Number mean: 0.097781 max: 12712.2
然后后处理发现本来应该随着刚体一起运动的近壁网格发生了变形,图如下(红色是库朗数大的地方,同样是网格变形的地方)
二维圆柱动的算例之前已经跑通过了,所以我先以为是constraints约束不够,后来加了几个约束后还是发散。 然后考虑到是不是icem没有关联好,cylinder边界重新关联了一下还是没有解决。
现在不知道问题出在哪里:有可能是感觉是innerDistance 0.1没有起到作用,但是使用slice切片观察Co发现圆柱中部的库朗数没有边界处这么大。所以不知道是不是约束没有设置合理。
[0_1607610864136_pointDisplacement](正在上传 100%)
附件中是相关的文件,不知道有没有人能帮忙解答一下。。网格太大了传不上来 ,,,
dynamicDictandpointdisplacement.zip -
求助:Openfoam如何输出湍流耗散率
目前在用icoFoam算圆柱绕流,算例和http://www.dyfluid.com/icoFoam.html中设置基本相同。
现在想生成无量纲化平均应变张量$S_{ij}=\dfrac12\dfrac{k}{\epsilon}\left(\dfrac{\partial \bar{u_i}}{\partial \bar{x_j}}+\dfrac{\partial \bar{u_j}}{\partial \bar{x_i}}\right)$
OpenFOAM利用fieldAverage可以直接生成平均速度与雷诺应力$R_{ij}$。湍动能$k$即$tr(R_{ij})$
平均速度梯度可以后处理“postProcess -func “grad(UMean)””得到。
那耗散率在openfoam里面可以怎么求呢?$\epsilon=\nu \overline{\dfrac{\partial u_{i}^{\prime}}{\partial x_{k}} \dfrac{\partial u_{i}^{\prime}}{\partial x_{k}}}$ 看定义是可以通过改编fieldAverage求。但是试了一下没有成功。。请问一下能通过什么方法求得呢。