补充说明一下,涡脱频率是根据forceCoeffsIncompressible监测的升力时程得到的,虽然频率不对,但是其他阻力升力的统计特性都是能和实验对的上的
Mania
帖子
-
-
最近在OF里用LES模拟宽高比5:1的矩形,计算出来的力的统计特性都是对的上试验结果的,但是旋涡脱落的频率只有正常结果的四分之一。
用的CFD的设置也算过其他模型,结果也都是对的。
在这里向各位老师请教一下可能是哪里出了问题?
因为和大论文相关,非常需要各位大佬的帮助,希望大家不吝指教,谢谢大家了! -
@mania 发现是我搞错了,用成了以前的不是DFSEM的算例,搞了两三天头晕了,捂脸
问题还是没解决@李东岳
这个二楼能删了么李老师 -
问题已解决。
发现原因可能是points文件需要的是一条入口inlet面的边上的线。
比如下图,我本来用的是线CD上所有cell的坐标,然后报错如主楼所说,我稍微看了一下代码,发现 if (!fullMatched) 就会输出错误信息如上。然后我试了试把points中的数据换成AB线上的cell坐标就可以了。
具体原因本工程专业的代码小白现在先不深究了,等有空了再看看。
-
最近在学习使用DFSEM湍流入口,需要在入口设置一些参数U、R、L,用的方法是mappedFile,相应的U、R、L和points文件在constant/ boundaryData/ inlet中。
在运行时候报错:
【
FOAM FATAL ERROR:(openfoam-2106)
Did not find a corresponding surfacePoint for every face centre】
之所以确定是因为mappedFile报的错,是因为我把入口改为均匀入口,不使用这些mappedFile就没问题。因此我猜想,是不是因为points文件中写的坐标的问题。(我选了inlet面200个点的坐标写入了points文件里,然后再其他三个文件中写入对应的U、R、L值。)
问题如上所述,希望各位前辈能指点一二。
-
@马乔 前辈你好,谢谢您的回复,也祝您中秋快乐。
看了您的回复,我的理解是,这个DFSEM算法,是根据我输入的R、L文件的值来计算出脉动速度(也就是您说的波动速度),然后加上U文件的值来生成相应的风场。
我看了tutorials中的case后,现在还有一些细节的地方不知道怎么理解,麻烦您指导一下:
1、对于U、R、L和point四个文件中数据的个数,可以随便定义若干个点就行?比如我inlet面(y×z)有5000个点,我就挑y方向一列800个点的数据,来输入对应的四个文件就行?2、这个U文件应该是输入的平均风速吧?比如我想要生成的流场是平均速度为(1 0 0),波动速度为0.5,那么U文件中所有的数据都应该是(1 0 0)对么?
3、相应的R和L文件,雷诺应力R应该就是[ui'uj']的一个对称张量,湍流尺度L应该怎么定义呢,看case中是一个列向量,那假如一个三维各向异性的湍流场是否需要输入一个n行3列的数据,如
{ L11 L22 L33 L11' L22' L33' L11'' L22'' L33'' L11''' L22''' L33''' ...... }
再次感谢您的回复,希望您空闲的时候能再指点一二,谢谢!
-
关于OF中DFSEM的设置,有一些问题请教各位前辈;
Q :0文件中U的设置
看OF用户手册里给的一个case:channel395DFSEM tutorial,地址是链接文本,它的入口代码是inlet { type turbulentDFSEMInlet; delta 1; U { type mappedFile; mapMethod nearest; } R { type mappedFile; mapMethod nearest; } L { type mappedFile; mapMethod nearest; } d 1; nCellPerEddy 1; scale 1; value $internalField;
这个case中U、R、L用的是把初始插值数据放入到constant文件中然后在计算开始时候赋值给入口(不知道我理解的对不对,如有错处恳请指导一二)
这里我有两个问题:
1、入口inlet面有多少个cell就在U、R、L三个文件中赋多少个值么?
2、问题1的基础上,如果是用U、R、L文件赋值的话,又该怎么给定这些数值。比如入口有4个点ABCD,如果每个点赋值一个UX 1 1 1 1,为什么不是1.1 1.2 1.3 1.4呢? -
@小考拉 您好,请问您成功编译这个代码了么,我现在也要做展向平均,有位前辈给我推荐了这个程序,但是我不知道怎么编译😂您能指教一下么
-
@benqing 您好,非常感谢您的回复!
我的确是想要做展向平均,看到了你给的链接太激动了,本来想着装好这个程序后再拿着成果来这个帖子跟您回复,也和大家分享一下,结果整了一晚上没有编译成这个代码😂折腾了一晚上没弄好还是厚着脸皮跟您请教一下。
我按照他的方法,先是把代码文件拷到根目录里面,
然后 wmake ,但是提示 command not found,
我又搜了一圈,说可能是没有编译环境变量,我就 sudo vi /etc/profile ——export PATH=<你要加入的路径>:$PATH,但还是不行
不好意思打扰您了,谢谢。以下是他的安装方法
info:- tool is based on the OpenFOAM utility 'postChannel'
- to compile the utility: source code should be copied to
WM_PROJECT_DIR/applications/utilities/postProcessing/miscellaneous/
'wmake' to compile the source code -
@benqing
您好,谢谢你的意见,我试了一下还是不行,如图。
现在的问题似乎是,就算新的变量计算成功了,也不能正确的生成图像,现在是困在这一步,也找不到原因😂您建议的用matlab算,我也试过了,问题是我的case中有好几个钝体,在matlab、origin里面作图的话,在本应该是空白的钝体的位置处,会产生很多杂乱的图(我猜应该是自动插值生成地)
很感谢您的回复和建议!!谢谢
-
@bestucan 大意了,没考虑到你说的问题😂
我也试过各种范围,都不能成。
问题似乎是,虽然算出来了正确的结果,但是显示出来图中每个单元的值都是0,实在没办法就不算这个东西了。
谢谢您的回复 -
有没有遇到过类似问题的大佬啊,再求助一次
被卡在这一步真的很难受T T
最后顶一次贴希望能被知道怎么解决这个问题的大佬看到,谢谢大家! -
求教啊,大佬们来帮帮忙
-
openfoam算完后用tecplot后处理,要对几个slice求平均,但是算出来之后的新的物理量无法正常显示,如图所示:
1 是我算的新的变量,就是对三个切面求平均
2 是算完之后显示有正常的数据范围range -740~720,但是怎么调整都没办法正常显示
3 是这个新变量的信息,也挺正常的求大佬帮忙解答一下,谢谢!!
-
@东岳 李老师好,我这边做的网格也很少,300万,renumbermesh后感觉速度上加快了没多少😂没太留意这个
-
@drink-water 果然是这个问题,师兄🐂🍺
最后还是麻烦师兄才解决了😂
谢谢师兄 -
@炎冰灬灬 检查网格看看呢,你也可以参考下这个帖子的设置
https://www.cfd-china.com/topic/664/cd计算不准的问题-les-re3900-三维圆柱绕流-pisofoam -
@炎冰灬灬 我之前用的codedFixedValue,直接把风速式子打出来。现在用oscillatingFixedValue,报错说不存在这个边界😂,我用的是OF7.0版本
-
@drink-water 计算完后就是reconstructPar。
”之前的“指的是什么 -
@drink-water 传了😂😂😂就是太离谱了所以刚开始不太想传
-
前提:算例做的是正弦来流边界。入口速度u=U0+Um*sin(2πωT),压力P为zerogradient;出口速度U为zerogradient,压力P为fixedvalue 0;上下面是symmetry;前后面是cyclic。
过程:算例完成后,输出的单个点的速度、压力数据都是对的,但是用parafoam得到的各种云图都是乱七八糟的图。
个人猜测:,①因为其他均匀流算例处理结果都是正常的,正弦边界条件设置的有问题;②其他原因
希望能得到大佬们的解答,谢谢!
-
@炎冰灬灬 你好,请问最后怎么调整成功的呢,我最近也在做正弦来流,遇到了一些问题。
-
怎么看自己有多少C币
😂回复完就看到了
-
@都别动就我能动 就在0文件夹里的U文件啊
-
@东岳 谢谢东岳老师!您的说法太形象了一下就理解了!
我刚刚也找到cyclic的设置方法了,一会儿就试试。还有一个问题,关于残差residual的也请教一下您。
在controldict设置了residual输出的残差文件,和图片里命令运行时候终端在每个时间步中显示的final residual一样么?
我得到的residual文件中的残差都是10-3量级的,和图片里的查了3个量级。最后再感谢东岳老师!
-
我在用cfd模拟正弦流场的钝体绕流,我自己是把上(up)下(down)前(front)后(back)四个部分都设置的symmetry对称边界,但是看论文里面一般都是上下symmetry,前后两部分(front back)设置的periodic。
实验室没有做cfd的人,来这里请教一下大家。谢谢!
问题总结
1、周期性边界和对称边界有什么区别呢,我应该怎么选择?
2、我自己先在openfoam(版本6)里面改成periodic试了一下,但是报错提示说没有这种边界,用cyclic也不行,报错说没有相匹配的面。所以如果要用周期性边界,应该用什么呢?希望又大佬能指导一下,谢谢!或者提示一下我能在哪里找到这些问题的答案?
-
@东岳 感谢东岳老师,意思就是 在 P 文件中写成
inlet
type codedfixedvalue
value -ρ*(V-V'/δt)
这种格式就好么
东岳老师写的这个Δt是时间步长么,V-V(t-Δt)要怎么表述呢
感谢东岳老师,谢谢你的指导 -
求大大们有空指导一下呀
-
想设置一个正弦来流,在论坛里搜了一下发现速度边界可以用oscillatingFixedValue,
想着先看看代码是怎么样的,但是在电脑里搜索oscillating关键字搜不到相关的文件我现在就有两个问题,一个是哪里能找到这个代码能看看
二一个就是速度入口是正弦来流,压力入口就是∂P/∂x=-ρ*∂V/∂t(应该是吧,望指正),那要怎么设置呢,以前只设置过0梯度的情况。我意思是这个压力条件的type该怎么写呢,就写个 codedfixedvalue?然后把∂P/∂x的结果写下么?萌新求指导,谢谢大家
-
@bestucan 好的,谢谢你的指点~
-
@东岳 还是谢谢东岳老师,我准备用一个小点儿的case在这个台式机上算一算,小的case只有300万个网格,看能不能成
-
@bestucan 请问怎么更新库呢,针对你说的 用了旧的库 的问题
-
@bestucan 谢谢您的指点,我这个同一个case在我笔记本里就算的很正常,所以对比看的话会怀疑是电脑问题
-
@东岳 东岳老师好,我用的case在我笔记本上就能正常算,但是换到这个电脑上就出现这种情况不行了
我用的是ras k-w SST simplefoam 模型,东岳老师说的重新编译是指的哪些方面呢?谢谢东岳老师 -
笔记本和台式机都是8g内存,i7处理器,
笔记本就算得很顺利,大概8个小时能算300步。
要是只用笔记本的话,就得电脑不眠不休算10天 -
这是因为网格太多了么,大概1300万个网格
但是我用我的笔记本就能算,因为模型太多了,想快点算完就又找了一个台式机来提高效率,结果台式机里面刚开始算就自动停止了
想请问有什么解决方法 -
@东岳 好的,谢谢东岳老师~~~
-
以前没学过c+,这几天看了看琢磨一下,
最后我是这样设置的INLET { type codedFixedValue; value uniform (0 0 0); name IN; code #{ const vectorField& Cf = patch().Cf(); vectorField& field = *this; forAll(Cf,faceI) { if( Cf[faceI].z()<300) { field[faceI] = vector(20*pow((Cf[faceI].z()/10), 0.15), 0, 0); }else{ field[faceI] = vector(30, 0, 0); } } #};
-
网上没搜到什么信息,在这里向大家求教一下,谢谢了!!
-
@calrel 好的~~~谢谢你,我前两天电脑画网格卡死了,没能上网及时回复你,现在试试,非常感谢
-
发现用代码高亮反而让人看的不方便
dimensions [0 1 -1 0 0 0 0];internalField uniform (0 0 0);
boundaryField
{
INLET
{
type codedfixedValue;
value uniform (x 0 0);
name uuu;//我看有的帖子里说name是一定要有的,但是后面的可以随便写,是这样的么Q2;还有的写的是redirectType XXX;,所以也想请教一下有什么区别,还是说是不一样的情况Q3code #{ const fvPatch& boundaryPatch = patch();//这个意思是获取inlet这个边界面么Q4 scalar z = boundaryPatch.z();//这里我是想定义高度z,但是不知道怎么定义,就这样写的,感觉肯定不对Q5 forAll( patch()) if z<300 x==20*(z/25)^0.15; if 300=<z<1000 x=29.03 #}; }
}
-
想把入口设置成大气边界层,
对于高度z小于300,速度x=20*(z/25)^0.15,
高度300=<z<1000,x=29.03
(整个流场高度就是1000m,所以这里就有一个小疑问,需不需要写上<1000的条件,Q1)代码我是先再论坛里搜了一下,然后自己琢磨的,没有写过代码所以希望哪位大神能指点一下
dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { INLET { type codedfixedValue; value uniform (x 0 0); name uuu;//我看有的帖子里说name是一定要有的,但是后面的可以随便写,是这样的么Q2;还有的写的是redirectType XXX;,所以也想请教一下有什么区别,还是说是不一样的情况Q3 code #{ const fvPatch& boundaryPatch = patch();//这个意思是获取inlet这个边界面么Q4 scalar z = boundaryPatch.z();//这里我是想定义高度z,但是不知道怎么定义,就这样写的,感觉肯定不对Q5 forAll( patch()) if z<300 x==20*(z/25)^0.15; if 300=<z<1000 x=29.03 #}; } }
希望有大佬能指点一下,非常感谢
-
@东岳 好的,了解了!!谢谢您
-
@东岳 在 请教,LES模拟中需要设置的参数 亚格子粘性 怎么计算 中说:
亚格子尺度
东岳大佬好!谢谢您的回复指导。
我是看一篇一个已经毕业很久的师兄写的文章,里面有这个说法——亚格子粘性nuSGS,不知道它是什么定义 -
看王福军的《CFD软件原理及应用》里没看到这个的设定,所以来请教大家,希望能有大佬来指点一下呀
-
@wwzhao 太感谢了,新手上路,昨天遇到好多首字母大小写报错的问题,最后看得眼花了没看到这个标点错误
感谢指正 -
新手刚刚入门,按照一个算例来做,但是碰到了如下图这样的问题。
搜了一下,发现有人的建议是
You have at least 2 choices for fixing this:
1.You properly set-up the boundary conditions inside the field files.
2.Or you rename the folder "0" to something else, such as "0.org" (it's the standard used in OpenFOAM's tutorials for the original boundary conditions to be used in the case).
但还是不太理解该怎么处理,求解决方法
OF模拟结果涡脱频率是正常结果的四分之一
OF模拟结果涡脱频率是正常结果的四分之一
请教mappedFile的相关问题
请教mappedFile的相关问题
请教mappedFile的相关问题
请教各位前辈OF中湍流入口DFSEM的设置问题
请教各位前辈OF中湍流入口DFSEM的设置问题
paraView-求展向平均问题
tecplot算一个新变量后无法显示出来怎么解决?
tecplot算一个新变量后无法显示出来怎么解决?
tecplot算一个新变量后无法显示出来怎么解决?
tecplot算一个新变量后无法显示出来怎么解决?
tecplot算一个新变量后无法显示出来怎么解决?
tecplot算一个新变量后无法显示出来怎么解决?
正弦来流边界算例的后处理出现错误,求解决
正弦来流边界算例的后处理出现错误,求解决
振荡来流的流场倒流问题
振荡来流的流场倒流问题
正弦来流边界算例的后处理出现错误,求解决
正弦来流边界算例的后处理出现错误,求解决
正弦来流边界算例的后处理出现错误,求解决
振荡来流的流场倒流问题
C币运行机制测试
请教一个周期性边界的设置问题
请教一个周期性边界的设置问题
请教一个周期性边界的设置问题
求指导,请教些关于oscillatingFixedValue边界的问题
求指导,请教些关于oscillatingFixedValue边界的问题
求指导,请教些关于oscillatingFixedValue边界的问题
openfoam计算出现new cannot satisfy memory request的问题
openfoam计算出现new cannot satisfy memory request的问题
openfoam计算出现new cannot satisfy memory request的问题
openfoam计算出现new cannot satisfy memory request的问题
openfoam计算出现new cannot satisfy memory request的问题
openfoam计算出现new cannot satisfy memory request的问题
openfoam计算出现new cannot satisfy memory request的问题
请教OF里面自带的算例有没有关羽k-w SST的求解器设置呢
请教一个入口非均匀风速设置的问题
请教OF里面自带的算例有没有关羽k-w SST的求解器设置呢
请教一个入口非均匀风速设置的问题
请教一个入口非均匀风速设置的问题
请教一个入口非均匀风速设置的问题
请教,LES模拟中需要设置的参数 亚格子粘性 怎么计算
请教,LES模拟中需要设置的参数 亚格子粘性 怎么计算
请教,LES模拟中需要设置的参数 亚格子粘性 怎么计算
Can not find patchField entry for XXX 的问题该怎么解决呢
Can not find patchField entry for XXX 的问题该怎么解决呢