@李东岳 是单塔式电脑
Fu
帖子
-
-
可以用newcommand
\newcommand\xx{\begin{equation}} \newcommand\xxx{\end{equation}}
\xx 1+1=2 \xxx
输出:
-
是这种吗
-
最近在研究圆柱绕流。 用的算例是李老师提供的算例:http://www.dyfluid.com/icofoam.html
现在我把圆柱从二维(Z方向一层网格)拉伸到三维计算圆柱绕流。 但是现在的forceCoeffs 和Forces function只能生成整个三维圆柱的升力和阻力。 如果我想生成沿管长方向的升力分布,例如均分成100个分段的升力。 Openfoam能够见到做到这点吗
隔壁有人提供了这种后处理方法,但是看起来需要每个时间步单独处理,而且三维的数据量太大了,,这样好像不能每个时间步自动处理。
These are my paraviewsteps in short:
1)Extracted Surface
2)Generate Surface Normals
3)Calculated:
cp = normals*(p2-p1); ->Pressure Drag
cf=tangentialsmuwallGradU_y -> Friction Drag (or calculated directly by foamCalc components wallGradU
4) Build Integral Drag=S(cp_X+cf_X)dA2/rhovvA -
最近也在研究这个函数。 看了半天没有明白。看了一下2.X的pisofoam中用的是turbulence->divDevReff(U). 这个比较好找。。 请问一下老师这两个可以当做是一个东西吗
-
@五好青年 tutorials里面的wingMotion
-
@dzw05 您好~ 那我可以理解成图上的差距就是SST模型模化的雷诺应力和真实的雷诺应力的差距吗?
-
@cccrrryyy 您好~
我将URANS的结果缩比成与Quasi-DNS相同范围则不能显示了。
关于收敛的问题。使用的都是非定常计算,fvSolution里面的值都是默认的。时间步长收敛性上,库朗数基本上在0.5左右。网格收敛性都经过验证了,而且用的网格都非常密了,这样一个流场用了15W个网格计算的。
所以这个结果是不合理的是吗。。。
-
算圆柱绕流,利用icofoam分别不用湍流模型(laminar)和SSTkomega湍流模型算出了两种工况
层流模型利用 fieldAverage生成雷诺应力
fieldAverage { type fieldAverage; functionObjectLibs ("libfieldFunctionObjects.so"); enabled true; outputControl outputTime; fields ( U { mean on; prime2Mean on; base time; } ); }
sstkomega模型利用
Postprocss -func R
生成雷诺应力但是两种方法生成的雷诺应力完全不同,quasiDNS的雷诺应力XX场呈现下图且随时间基本不变
SSTkomega的雷诺应力XX场为
且像涡街一样随着时间而变化。两者算出来的雷诺应力差别有这么大吗? 还是我哪里的计算出了问题?
-
gradU输出格式是
-
没错吧。。
《Fluid Mechanics》 Frank M.White P239
《工程流体力学》 李玉柱弹性力也是这么定义的
-
最后找到了解决办法,是边界条件的原因。本来上下圆柱连接的上下两个边界用的symmetry边界条件,改成cylic边界之后就能算了。。。具体原因还没研究出来
-
@赵海盛 看起来没有问题的。。如果是自激振荡就是caculated,如果是强迫就要oscillatingdisplacement了。。。看边界条件是一个三维的圆形流场吗?。。
-
二维还是三维呢。。边界条件问题吧一般是
-
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
正常了。。。
-
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计算一个弹簧连接的三维圆柱在来流下的振动情况,但是结果发散。
在算到发散的前几步之前,库朗数开始急剧增加。
Courant Number mean: 0.097781 max: 12712.2
然后后处理发现本来应该随着刚体一起运动的近壁网格发生了变形,图如下(红色是库朗数大的地方,同样是网格变形的地方)
二维圆柱动的算例之前已经跑通过了,所以我先以为是constraints约束不够,后来加了几个约束后还是发散。 然后考虑到是不是icem没有关联好,cylinder边界重新关联了一下还是没有解决。
现在不知道问题出在哪里:有可能是感觉是innerDistance 0.1没有起到作用,但是使用slice切片观察Co发现圆柱中部的库朗数没有边界处这么大。所以不知道是不是约束没有设置合理。
[0_1607610864136_pointDisplacement](正在上传 100%)
附件中是相关的文件,不知道有没有人能帮忙解答一下。。网格太大了传不上来 ,,,
dynamicDictandpointdisplacement.zip -
目前在用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求。但是试了一下没有成功。。请问一下能通过什么方法求得呢。
-
可以在BIOS里面关掉超线程,计算速度感觉快了一点,但是没有详细测试。。。
200万网格并行算力测试(OpenFOAM版本)
latex这个代码怎么省略
三角形面积计算公式
三维圆柱绕流升力沿管长分布
三维圆柱绕流升力沿管长分布
turbulence->divDevSigma(U),不可能完成的任务
OpenFOAM中开发浮体系泊(锚链)模块
使用SSTkomega和laminar算出来的雷诺应力问题
使用SSTkomega和laminar算出来的雷诺应力问题
使用SSTkomega和laminar算出来的雷诺应力问题
关于postProcess得到张量的排序
提取log中数据问题
《无痛苦N-S方程笔记》勘误
三维sixDoFRigidBodyMotion算例发散
sixDoFRigidBodyMotion设置参数
用pimpleFOAM做圆柱的涡激振动遇到收敛性问题
用pimpleFOAM做圆柱的涡激振动遇到收敛性问题
200万网格并行算力测试(OpenFOAM版本)
200万网格并行算力测试(OpenFOAM版本)
三维sixDoFRigidBodyMotion算例发散
求助:Openfoam中的LES/DNS如何输出湍流耗散率
超线程到底要不要开启