@yxg2018 谢谢,我将网格按照单元类型分离并将边界层网格也设定为刚体运动区域,计算依旧会报错,提示Divergence detected in AMG solver: x-momentum,我将迭代次数增加并改小了动量的松弛因子,计算还是发散。
ycx1042341217
帖子
-
-
@yxg2018 我现在动网格是采用弹簧更新和网格重构,将圆柱边界设置为刚体运动
-
各位老师,我在做亚临界雷诺数下二维圆柱的涡激振动(y向),目前模拟的工况是Re=5000,圆柱直径0.0381m,速度0.132m/s,约化速度为8,质量比2.6,阻尼因子0.0036,振动方程通过编写四阶龙格库塔的UDF导入到fluent中,振动方程进行了无量纲化。
模型如下,上游10D,下游20D,宽度上下各15D,圆柱附近采用四边形网格,其余为三角形网格。第一层网格高度0.0023D,时间步长0.0014s,湍流模型k-wSST。
在低雷诺数时基于编好的UDF进行过计算(此时的模型不一样,直径0.0016m),与实验数据接近,但在亚临界时,计算开始一段时间后圆柱位移一直增大,直到出现了负体积,此时的位移已接近10D。自己更改过时间步长,但是问题也是一样,想问下有老师做过类似的模拟吗?这个问题排查了很久还是没能解决,谢谢大家!
编写的UDF如下。
-
@random_ran 好的谢谢老师
-
@random_ran 老师可以帮忙看看吗?我后面也换了k-omega SST模型,结果也相差不多,看了一下也并没有漩涡脱落产生。
-
谢谢东岳老师
-
大家好,我最近在模拟Re=10000时,二维圆柱的升力系数和阻力系数。采用的湍流模型为k-omega,网格划分采用非结构化网格,边界层采用四边形网格,圆柱壁面y+值在1附近。换过不同的网格,时间步长等,阻力系数为1.0-1.2左右,与文献可以对应上,但是升力系数模拟出来在1.25左右,文献数据为0.4-0.6附近。大家能帮着看看问题出在哪里嘛?fluent里再选择k-omega模型后还可以选择低雷诺数修正或者自由剪切修正,我在计算的时候并没有勾选,大家知道这个一般怎么用吗?
-
各位老师好,我在画圆柱扰流网格的时候想在圆柱周围加边界层,再在一定范围内加密,该怎么实现呢?我现在试了两种方法。
1.分成了两个区域,但是由于后面圆柱需要运动,当运动范围超过紫色的部分时,就会报错。两个区域中间使用了interior,静止时没影响,但网格运动应该是全局网格,如果有interior其实是不正确的。
2.采用了密度盒,但是该方法要求网格划分采用Patch Independent,这样就不会遵守我设置的全局网格参数,而且设置的圆柱边界层网格也生成不了,直接全部都是三角形网格。
想请问各位老师,应该如何实现呢? -
或者有老师知道如何将ICEM的模型输出为.stp格式吗?
-
在ICEM中建立模型如下图
通过两个圆面和圆柱侧面形成一个名为soild的part。ICEM生成.stl文件后导入到Design Modeler中。
但是在DM中,圆面和圆柱侧面均变成了属性为soild的part。
请问在ICEM中如何设置part可以使得圆柱的三个面在DM中以一个soild形式生成? -
@东岳 嗯嗯,确实是尺度问题,谢谢东岳老师!
-
@ibelief 确实是scale的问题,谢谢老师!
-
各位老师,我在用fluent模拟圆柱二维扰流时,采用瞬态求解器,雷诺数200,但是无法看到漩涡脱落现象,加大雷诺数后也没有变化,只能在圆柱周围看到少量的涡量线。请问老师们问题可能出现在什么地方呢?网格和相关设置都是跟着书里算例走的,所以不清楚问题在哪里,谢谢!
-
大家好,我用OpenFOAM计算三维风力机的转矩和转矩系数时,结果和文献相差大约10倍,风速10m/s,转速31.4rad/s,风机叶片高度3m,旋转轴高4m,风轮半径1m。计算转矩时对整个风机求解。大家可以帮忙看看问题出在哪里吗?谢谢! (用gambit划分网格,圆柱区域是滑移区。) (fvschemes文件) (controldict文件)
-
@dzw05 湍流模型中有一个运动黏度nu,正常情况下在后面赋值即可。我当时在网上查询如何提取升力系数和阻力系数时有篇帖子提示要在nu后面先再输入一个nu字符再赋值,由于之后的计算我都用的是这个错误的湍流模型文件,导致阻力系数求解一直出现问题。后来一个一个文件排查的时候发现是这里的问题。
-
@dzw05 老师我排查出来是我的湍流模型文件里面的参数设置出了问题,所以导致阻力系数一直存在问题,现在已经解决了。十分感谢老师这几天的帮助!
-
@dzw05 老师,我修改了边界条件,对网格也进行了一些修改,但是重新计算后阻力系数还是基本没变。会不会是采用的SA模型本身对于这种情况模拟就不是很好?
-
@东岳 老师您说的分辨率指的是gambit中划分网格时这个方向的successive ratio吗?
-
@dzw05 老师,我这两天重新划分网格的时候发现之前计算过程中郎库数非常大,现在我根据郎库数为1重新设定了时间步长,但阻力系数依旧不正确。我在gambit网格质量检查的时候,远场质量为0,翼型附近质量为0.1-0.2。不知道这样的网格质量可以吗,或者还有什么其他的可能问题吗(比如说我的nut,nutlida值会不会设定有问题),因为我现在修改了很多次网格,但是阻力系数总是维持在10以上。谢谢老师!
-
@dzw05 老师您好,我的第一层网格高度大约为0.18mm,昨天跟您说错了。刚才重新划分网格后,checkMesh显示正常,但是进行计算后升力系数和阻力系数都与实验数据有偏差。老师我把我的网格文件,controlDict文件和初始条件传上来,您能帮我看看问题出在哪里吗?谢谢老师!4_1523608055301_U.pdf 3_1523608055301_p.pdf 2_1523608055301_nuTilda.pdf 1_1523608055301_nut.pdf 0_1523608055300_controlDict.pdf [网格](链接: https://pan.baidu.com/s/129VWe7pswJfTMWmG4oHxrg 密码: x7g4)
-
@dzw05 老师,我采用的是NACA0018翼型,弦长1米,速度1.463m/s,雷诺数100000,采用的是SA模型。不同攻角下升力系数较为吻合。阻力系数基本维持在14.多。之前划分网格的时候是跟着网上的一个教程划分的,没有考虑到y+。刚才我重新划分了网格(根据我的数据,按照y+等于1计算,我的壁面第一层网格高度大约为0.08mm。),但是在OpenFOAM里面checkMesh后,提示我High aspect ratio cells found, Max aspect ratio: 10536.4, number of cells 1020。老师我是不是需要重新划分网格?那为了避免这种情况,网格的纵横比最好不要超过多少呢?老师您看看还需要我的哪些数据吗?谢谢老师!
-
@dzw05 明白了,十分感谢!不知道您对于我不同攻角下算出的Cd系数基本不变这种情况有了解吗?-
-
@dzw05 您好,想请问一下这个方向矢量在写的时候有什么规范吗?我在给定liftDir的方向时第一回给定(0 1 0),第二回给定(0 2 0),从矢量角度看是相同的方向,结果第二回计算出的Cl是第一回的二倍。想请问这种情况应该怎么处理,应该如何确定准确的方向矢量数值?谢谢!
-
您好,想请问一下这个方向矢量在写的时候有什么规范吗?我在给定liftDir的方向时第一回给定(0 1 0),第二回给定(0 2 0),从矢量角度看是相同的方向,结果第二回计算出的Cl是第一回的二倍。想请问这种情况应该怎么处理,应该如何确定准确的方向矢量数值?谢谢!
-
大家好,我现在用simpleFoam求解NACA二维翼型升力系数和阻力系数,不同攻角下的升力系数值和实验数据基本吻合,但是不同攻角下的阻力系数值却基本不变且和实验数据相差较大,我采用的是C型网格,计算域长度是弦长的33倍左右,网格数量6840。大家有知道是什么原因造成的吗?谢谢大家!
还有一个问题是在提取升力系数和阻力系数时,liftDir和DragDir是要随着攻角不同而变化还是一直为(0 1 0)和(1 0 0)? -
@程迪 老师您好,想请问一下您说的二维翼型的参考长度和面积怎么取值呢?我在算的时候参考长度取的是弦长,参考面积取的是翼型的截面积,但是Cl和Cd都差很多。
-
学生在用simpleFoam 求解NACA0012翼型的升力系数和阻力系数时,所得结果Cd100多,Cl 20多。学生采用速度1.6m/s,攻角8度,弦长1m,SA模型。U文件中入口和出口类型均为freestream,翼型为fixedValue;P文件中入口和出口类型均为 freestreamPressure;nut和nuTlida文件中参数为0.14。controlDict文件密度输入1.225,lref1,aref0.08。学生是刚学习openfoam,很多东西还不是很清楚。请各位老师帮忙看一下问题出在哪里,不胜感谢!
亚临界雷诺数下二维圆柱涡激振动出现负体积
亚临界雷诺数下二维圆柱涡激振动出现负体积
亚临界雷诺数下二维圆柱涡激振动出现负体积
k-omega湍流模型模拟二维圆柱升力系数问题
k-omega湍流模型模拟二维圆柱升力系数问题
k-omega湍流模型模拟二维圆柱升力系数问题
k-omega湍流模型模拟二维圆柱升力系数问题
ICEM网格加密
ICEM模型向ANSYS DesignModeler转换
ICEM模型向ANSYS DesignModeler转换
fluent新手模拟二维圆柱扰流无漩涡脱落现象
fluent新手模拟二维圆柱扰流无漩涡脱落现象
fluent新手模拟二维圆柱扰流无漩涡脱落现象
三维风力机转矩系数Cm值较大
simpleFoam求解二维翼型Cd不变
simpleFoam求解二维翼型Cd不变
simpleFoam求解二维翼型Cd不变
simpleFoam求解二维翼型Cd不变
simpleFoam求解二维翼型Cd不变
simpleFoam求解二维翼型Cd不变
simpleFoam求解二维翼型Cd不变
simpleFoam求解二维翼型Cd不变
simpleFoam求解二维翼型Cd不变
simpleFoam求解二维翼型Cd不变
simpleFoam求解二维翼型Cd不变
openfoam SA模型计算空翼出现Cl和Cd算不准的问题?
simpleFoam求解NACA0012翼型