@五好青年 我用cfd post ,直接提取就好
宝丁
帖子
-
-
@五好青年 你是用什么软件做的计算和后处理
-
@五好青年 我一般是在cfd post输入马赫数求解公式算的
-
@东岳 只求毕业啊,不说了,要去写论文了,现在在等一波计算结果出来,只求是最后一波了
-
@宝丁 P*=115600
-
@paopaoai11
首先,目前咱们共同的观点在于吸力面近前缘处局部加速主要是由该处流线曲率决定的,而不同的观点在于曲率导致前缘加速现象的机理是什么。
我的观点在于曲率的存在使得来流虚拟流道被压缩,因而流动流速增加。这里对虚拟流道加以解释:壁面的存在使得流场中产生了相应的扰动,不同壁面所对应的扰动传播的距离是不同的,因此将壁面至扰动传播的边界称之为虚拟流道,其是变化的,而非一成不变的。目前这个观点的问题在于我还没有定量的去描述这个虚拟流道,但就流动细节,这又是客观存在的。
但是凹凸曲面的虚拟边界相同吗?用凹凸曲面举例只是对为什么前缘加速只在吸力面一侧显著的复述罢了。注意,这并不是说压力面一侧没有。
其次,我不同意你的这一观点:
叶栅前缘吸力面局部加速跟流动滞止没太大关系,因为前缘滞止点压力最大也不会超过气流总压,所以由顺压梯度造成的加速最大也不会超过主流速度。
这个观点暗含了一个条件:流场中静压处处小于进口静压,但是这并不现实。在本问题中,叶栅前缘的最小静压可是要比进口静压小很多,因而并不一定能保证顺压梯度造成的加速最大也不会超过主流速度。
此外,关于前缘加速的原因我想通过圆柱扰流的算例加以说明:
下图给出了圆柱扰流瞬态静压与速度云图,流体首先在圆柱前缘滞止,静压升高,而在尾迹区,分离涡的存在使得该处为一个低压区,因此在圆柱上下表面存在一个较大的静压梯度,气流加速。
综上所述,曲率的存在使得来流虚拟流道被压缩,流体流速增加,其动力来源是静压差。 -
@paopaoai11 另外,关于前缘流动的压缩,这个是指的流体虚拟流道的面积被压缩,因而流体流速增加,压强降低,密度减小,(比体积增大 )
-
@paopaoai11
另外,关于前缘流动的压缩,这个是指的流体虚拟流道的面积被压缩,因而流体流速增加,压强降低,密度增大。 -
@paopaoai11
非常感谢您指出我的一些错误,受教了。我再结合您的回答说一下我的理解哈。
第一点是前缘滞止所造成的顺压梯度只是一个流动加速的有利因素。
第二点是流体在前缘加速,动压升高,静压降低,确实是对应一个膨胀过程,说压缩确实不合适。
第三点是曲率对流动的影响,我当时的理解是受到前缘壁面的扰动,局部存在一个收缩的过程,也就对应曲率的影响。 -
@东岳 等我写完大论文再好好整理一下
-
@东岳 流体流至叶片前缘,受到阻滞而形成驻点,速度滞止,静压升高,因而会有一个顺压梯度;叶栅吸力面是一个凸曲面,流体流至叶栅前缘,存在一个局部压缩的现象,流体速度增加。在以上两者的共同作用下,叶栅前缘会出现局部高速。
叶栅流道是一个扩张型通道,起到减速增压的效果,因而整体而言速度会降低,同时静压升高,通道内会存在一个逆压梯度。在吸力面一侧,尤其是叶栅中后部,当流体惯性力不足以克服当地逆压梯度所产生的“阻力”时,会存在流动分离现象,在速度云图上就会出现一个低速区。
-
这个算的感觉很好啊,东岳大神能不能把文件发到我邮箱,baoding@mail.nwpu.edu.cn
还没有毕业,准备写打毕业论文呢,就是换了一个题目而已
-
两个面设置完交界面后会生成一个内部面,设置成rans-les-interface就可以了
-
在用fluent中ELES/ZonalLES 发现其交界面的处理很是不方便,下图是fluent帮助文档中关于交界面设置的说明,很显然这个需要内部面(interior),也就是交界面节点一一对应,但实际中这很难做到,因此常用interface,但是这样就选不到rans-les-interface。如果用一般的交界面感觉算出来的结果有点怪,有人遇到过类似的情况吗?
-
@东岳 嗯嗯,好的
-
@东岳 好的 好的
13 six 36 seven 12 seven three two
微信和手机号同号
-
@东岳 嗯嗯 我这边会继续跟进的。
原先想着毕设做这个的,但是没想到卡在这里。现在因为毕业的压力,所以只能用Fluent或者CFX先往前做事情。
如果有新的进展,我也会发到这里跟您讨论的。
-
@东岳 是15楼的图,原来的叶栅试验结果和计算结果对不上,后来就放弃了。
从文献找了NACA65-810的试验和计算结果(下图),发现fluent和试验对的非常好,但是rhoPimpleFoam在趋势上就差很多,计算结果就没有放在图上。
-
@东岳 是这样,求解格式如果全部都是二阶迎风求解确实很稳定,但是和试验或者Fluent相差很大,所以想着改成中心差分,但是中心差分格式不稳定,第一个叶栅可以收敛,换一下就发散了,异常诡异。
-
@东岳 计我这里计算一开始没什么问题,但是一段时间以后就会发散啊
-
最后一次自己顶一下这个帖子
最近尝试了很多,从差分格式到求解和算法控制,一直都没有得到合理的结果,应该是求解器本身的问题了。但是自己时间限制,马上研三了,再做下去心里实在没底,准备去用自己ANSYS 那一系列的软件做事情去了。
但是还是有些不甘心,毕竟了尝试了那么久。我的求解文件和一些网格都在上面,如果有人知道我的错误在哪里,请一定要告诉我,不论时间过去多久,我都愿意拿自己现在一个月的助研费作为酬劳(400,不多哈)。
最后还是感谢东岳老师一直以来的帮助。
-
@东岳 首先非常感谢东岳老师
就我目前能找到的试验数据,Naca 65-810,fluent 计算和试验值相比较准一些(以下是该叶型计算值与试验值的对比),但是OpenFOAM 的计算结果在趋势上就差很多,所以个人倾向于fluent的计算结果。下面给出NACA65-810的计算网格以及fluent的计算结果与试验值的对比。
NACA65-810网格(y+~30-60) 链接:https://pan.baidu.com/s/1tZPl2NfspkOqjmoV4rSpvQ
提取码:jhg4 -
最近将离散格式全部改为中心差分格式(原先全部是二阶迎风),计算的结果还比较让人满意,但是这个差分格式还是不太让人满意,并不稳定,对于这个算例最终可以收敛,收敛的结果如下图所示。但是只要换一个叶片马上又发散了,现在真的是很疑惑,前缘的这类现象究竟要怎么调,或者需要动求解器内部的结构了吗?
divSchemes { default none; div(phi,U) Gauss limitedLinearV 1;//Gauss linearUpwind limited; div(phi,k) Gauss limitedLinear 1; div(phi,omega) Gauss limitedLinear 1; div(phi,e) Gauss limitedLinear 1; div(phi,h) Gauss limitedLinear 1; div(phi,K) Gauss limitedLinear 1; div(phi,Ekp) Gauss limitedLinear 1; div(phiv,p) Gauss limitedLinear 1;//Gauss linearUpwind limited; div(phid,p) Gauss limitedLinear 1;//Gauss linearUpwind limited; div((phi|interpolate(rho)),p) Gauss limitedLinear 1;//Gauss linearUpwind limited; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; }
-
我把我的求解文件放到这里,大家有兴趣的可以帮忙瞧一下。
链接:https://pan.baidu.com/s/15zFl76UsfD_XyWwUYEdyGQ
提取码:o3ae -
@Foamer24 用过层流、S_A和kOmegeSST,结果并没有本质的区别,前缘加速不明显,低速区域影响范围过大
-
@东岳 利用NACA0012,算了一把外流,速度和压力分布都很合理;
但是对与我要算的叶栅而言,不论如何调节差分格式,更改网格,计算结果都不让人满意,现在怀疑是算法的差异 -
@Foamer24 +1
-
@Foamer24 哈哈,既然粗网格可以插值,那相同的网格不是更没有问题吗
-
@Foamer24 有一个mapField,可以试一下,另一种也是我现在使用直接修改计算收敛值最后一个时刻的边界条件(比如总压),作为初始流场计算
-
@Foamer24 这个做法只是辅助性的,应该不会影响本质的算法不稳定
-
@Foamer24 温度超限或者参数错误
-
@东岳 这个离SCI还很遥远,应该是两个求解器之间的差别吧,实验和计算对起来也很麻烦,我再做两个算例
-
这是叶片表面等熵马赫数分布,可以明显看出叶片前缘那一块fluent和rhoPimpleFoam有很明显的差别 -
来流马赫数都是0.69
-
一直尝试用OF求解叶栅这类的跨音速内流流动,原先用rhoCentralFoam和sonicFoam一直不收敛,然后尝试利用rhoPimpleFoam 计算,本来不抱什么希望,但是计算竟然收敛的挺好,但是计算的结果并不理想,和Fluent 计算得到的结果相比,其叶片前缘的局部高速区域明显很弱,这是不太符合一般的气动规律。
原先认为是壁面yplus的影响,但是将yplus从(30,60)(标准壁面函数)调到了1(低雷诺数壁面函数),但是情况并没有太大改观。
感觉每走一步都会遇到很多困难,这应该是与我原先只是将OF 作为一个求解器使用的缘故吧。
-
是不是你的几何精度问题,壁面网格这么密,画网格的时候如果几何精度太差可能画出来的网格可能会有问题。这里的几何精度既包括建模软件里面的精度,也包括画网格软件转换的精度。如果直接在画网格的软件里面画几何精度应该会好一些,ICEM里面会有调精度的选项,Pointwise不清楚
-
这个链接里的算例是我现在迭代的最新的算例设置了,我再整理一下问题,希望浏览的人如果遇到过类似的情况的人能给些建议,我只想知道我错在哪里了。
算例可以稳定的收敛一段时间,但是出口背压不能固定,出口边界是 waveTransmissive,出口背压不断减小,进口静压也是减小,进口总压一定,进口速度也会不断增大,当进口Ma数大于1的时候整个计算就会崩掉。
如果强行换成fixedMean或者fixedValue出口背压就会不断增大,出口就像被堵住了一样。链接:https://pan.baidu.com/s/1hzOJcuYQw0GOYWyNK_TqOw
提取码:sr9c -
@东岳 好的,谢谢老师的关注,原先您能给的指导我已经很感动了,希望您的工作顺利
-
@东岳 老师,您能否给些建议,我现在应该从哪方面去想办法,现在我觉得已经试了很多,但是依然找不到方向。
-
感觉就差最后一点,但是就是做不出来
-
@东岳 还是用的sonicFoam, 就是出口背压一种下降,虽然可以收敛一会,但是紧接着进口Ma 数超音以后整个计算就崩掉了
-
目前来看这个问题是不是已经要修改OF的一些底层代码
不知道这个难度如何,应该要放弃这个算例了 -
这应该与我的出口边界条件有关
/*type fixedValue; //meanValue $pOut; value uniform 1e5;*/ type waveTransmissive; field p; psi thermo:psi; gamma 1.4; fieldInf 120000; lInf 0.2; value uniform 110000;
原先想用fixedMean这个边界条件,但是一用出口边界的压力值是固定了,但是出口附近的压力会异常的高,好像出口把流道给堵住了一样;
但是用 waveTransmissive这个边界条件,边界的压力值就会一直下降,进口的压力也下降,速度升高,Ma数达到0.9左右整个计算就会崩掉。 -
@宝丁 出口静压不能固定
-
计算虽然可以收敛,但是出口总压却一直不能固定,而且进口马赫数一旦达到0.9,整个计算就会崩掉,这难道是求解器的问题?
-
@东岳 我把库朗数降到了0.006,计算发现边界还是发散。rhoCentralFoam本身是一个跨声速求解器,为什么要满足高超音速的假定?
-
@东岳 好像理解错了,降低时间步长试一下
-
@东岳 我的边界条件远没有超音,进口马赫数约0.55, 但是rhoCentralFoam只能用于高超声速吗
-
入门不久,尝试着答一下。
我认为openFoam中的边界条件大致分两类:- 强制性边界条件,直接规定数值,如fixedValue,fixedMean,这类边界条件的信息会一直存在于求解所得的结果文件,算很多步以后查看文件就会发现边界条件并没有变化;
- 参考性边界条件,这类边界条件大多是第二类边界条件,如出口定义inletOutlet, 对于流出的气体而言,该边界其实是zeroGradient,这时定义的value只是初始迭代的时候给的一个初值,该边界的数值由流场内部的点插值得到
-
上一个帖子是关于sonicFoam使用总温总压入口时遇到了一些问题,问题似乎已经解决,再次感谢李东岳老师的帮助。
现在换一个求解器,用rhoCentralFoam又有些别的问题。
还是给出我要用到的边界条件
计算的问题在于进口的第二层网格会出现数值的异常,如下
我将文件放到这里,大家有兴趣的可以看一下:420_rhoCentralFoam.zip
这样的问题在CFD_online上面也是有人问到过,但是并没有解决,这是链接:
https://www.cfd-online.com/Forums/openfoam-solving/81423-rhocentralfoam-totalpressure-totaltemperature-inlet-subsonic-channel.html还是说rhoCentralFoam不支持总温总压的入口边界条件呢?
还有,我将总温总压的入口边界条件换成了速度入口,整个计算收敛的非常好。
有没有人同样遇到这样的问题呢?虽然以后用不到这个求解器,但是还是想把这个问题解决。
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
关于Zonal LES交界面设置的问题
关于Zonal LES交界面设置的问题
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
rhoSimpleFoam 求解可压缩流动时一个奇怪的现象
OpenFOAM中cyclicAMI周期边界在有边界层时出现问题
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
rhoCentralFoam使用总温总压入口时遇到的问题
rhoCentralFoam使用总温总压入口时遇到的问题
rhoCentralFoam使用总温总压入口时遇到的问题
给定初值时,边界条件如何发挥作用?
rhoCentralFoam使用总温总压入口时遇到的问题