@Afr1yne 如果模型什么都没变,只是网格加密, u k e正常也是一样的。把fvSchemes换成最稳定的格式计算看看,再逐个改二阶格式试算,看哪个设置导致的发散
coolhhh
帖子
-
同一个case更换网格后计算发散 -
同一个case更换网格后计算发散@Afr1yne 如果普通的均匀网格都发散,有点难搞,得再仔细查一查什么原因。fvSchemes和fvSolution改为更稳定的格式试试
-
同一个case更换网格后计算发散基本不用修改,发散大概率是网格问题。看看
fvSchemes
和fvSolution
怎么设置的,有的格式在好点的网格能算,复杂的网格就发散 -
wallShearStress计算结果问题对于不可压缩湍流,of输出的wallShearStress还要乘以空气密度。另外是否可能实验的数据无量纲化处理了?找下类似的研究,别人的文献数量级是多少,判断哪个计算错了
-
在流场内部加一个速度场我没具体操作过,正常画网格时候就要像你那个图下面挖空一样,额外挖空一个小长方形
-
在流场内部加一个速度场@Afr1yne 设置边界要掏空一个区域,源项参考这个帖子:关于OpenFOAM中的fvOptions的作用
-
在流场内部加一个速度场@Afr1yne 我的理解是旋翼在这个地方本身就是有阻挡作用,只是很薄而已,而不是凭空有速度。因为加的是速度,不是源项,强行把这部分体网格速度改为目标速度,那么继续用pisoFoam之类的计算,就根据边界条件和初始条件计算更新整个U场,这部分体网格速度也就会改变,无法保持原来的速度。只有边界条件能一直保持原来的速度。
-
在流场内部加一个速度场能否在方框这里掏空一个很小的矩形,模拟旋翼,四条边有的设置为速度边界,没有产生速度的设置为wall
-
foamMeshToFluent如果是Fluent Meshing画的网格,要用
fluent3DMeshToFoam
-
有偿提问:snappyHexMesh多级加密,能否指定单个方向进行加密?@chenboyao 参考 这个帖子
-
在fluent里面画的网格植入到了of中,修改周期性边界发现导致映射的id是倒序的@yzxwh 没关注过cell的id是否改变,但感觉应该没变,你再确认一下
-
在fluent里面画的网格植入到了of中,修改周期性边界发现导致映射的id是倒序的@coolhhh 这是of2206写法。如果是org版本,记得of8关键字和平移向量写法记得有点不同,类似与下面,平移向量的值与of2206的值互为相反数
transformType translational; separation (0 -0.3 0);
-
在fluent里面画的网格植入到了of中,修改周期性边界发现导致映射的id是倒序的@yzxwh 增加
translational
关键字试一试,由于算例没有Allrun,没试算对不对。patches ( { name CYCLIC1; patchInfo { type cyclic; neighbourPatch CYCLIC2; matchTolerance 0.02; transform translational; separationVector (0.3 0 0); } constructFrom patches; patches (left); // 源面 // set f0; } { name CYCLIC2; patchInfo { type cyclic; neighbourPatch CYCLIC1; matchTolerance 0.02; transform translational; separationVector (-0.3 0 0); } constructFrom patches; patches (right); // set f0; } );
转换后boundary文件:
-
OpenFOAM RANS方程采用涡粘模型时算得的压力是不是包含2/3kOpenFOAM没有做处理,可以参考 这个帖子 中44#楼的两个链接,有大佬讨论过这个问题:
)
-
openfoam7用fieldAverage函数无法获得热释放率Qdot的平均值?Qdot在计算的时间文件夹中是否存在,不存在就无法用fieldAverage统计
-
fluent速度入口造波 波幅是0.07m 但跑出的报告 没有波峰可以达到0.07m,请问怎么解决呢@2682577721 没做过造波,但猜测会不会跟网格粗糙有关。加密网格,或者把波长变大,即波数减小,可能结果会好
-
中性大气环境湍流动能的自保持 | 附有算例下载@李东岳
李老师,我之前测试的Yang模型算全尺寸无法自保持,湍动能完全衰减没了(见46#楼结果)。HW07我没算过,但主要看图片应该是算的全尺寸(因为k的高度已经500m了),但k的数量级感觉偏小,这主要是平均速度取值偏小导致的,猜测HW07算全尺寸的目标值取得不太符合实际 -
中性大气环境湍流动能的自保持 | 附有算例下载@小李的幸福生活 我没模拟过大跨结构,平时主要用LES,RANS模拟少。看模拟高层建筑的更多使用 SST k-w模型,主要看类似的论文怎么设置,跟着设置基本问题不大。下面两个链接讨论了如何选择模型,可以看看:
选用合适的湍流模型
如何选择湍流 CFD 模型? -
中性大气环境湍流动能的自保持 | 附有算例下载@小李的幸福生活 最好就跟风洞试验一样设置缩尺模拟对比,风洞试验怎么设置,数值风洞就怎么设置。全尺寸模拟的话y+很难满足
-
中性大气环境湍流动能的自保持 | 附有算例下载@小李的幸福生活 由于我主要不是研究这个内容,后面没再尝试了。我的理解是全尺寸和缩尺转换,相关的参数应该都要转换。一个简单方法验证变换的对不对,可以把缩尺后的
k
剖面画出来,看下数量级是否合理的。比如现在算全尺寸,平均风速取$U_{avg}=40\ m/s$,湍流强度取$I_u=10\%$,近似取三个方向湍流强度一样大,那么估计的湍动能数量级应该为:$k=\frac{3}{2} \sigma_u^2=\frac{3}{2} (U_{avg}I_u)^2=24\ m^2/s^2$。比如7#楼中HW07方法用的目标值,感觉有点偏小,主要原因全尺寸的平均速度剖面取值也偏小只到12m/s这个方法简称HW07,参考的是
D.M. Hargreaves and N.G. Wright 2007
,目前这一套方法计算出来的结果可以把k保持住。