@coolhhh 谢谢您!我尝试做一下。
strive_DUT
帖子
-
-
@coolhhh 我的速度分布是从文件中读取的,射流入口速度是按照下图给出:
-
@coolhhh 您提到的这个
turbulentDFSEMInlet
应该跟InflowGenerator差不多,我看是在V版本中的。但是我的速度分布跟时间相关,我想在这个基础上人为的添加随机扰动。 -
如题,我想实现一个依时类边界与人工湍流的耦合,比如以下两种。请问各位大佬该如何实现呢?给个提示也行
依时类边界:type timeVaryingUniformFixedValue; fileName "$FOAM_CASE/150MPa-U"; outOfBounds clamp;
人工湍流:
value uniform (600 0 0); type decayingTurbulenceInflowGenerator; direction 1; LField uniform 0.00011;//integral scale = width refField uniform (600 0 0); RField uniform (0.0729 0 0 0.0729 0 0.0729); Umean (600 0 0);
-
@soulx7 按照模版正常使用就行,或者你搜索一下这个边界,有很多相关说明。
-
@soulx7 就是OpenFOAM自带的无反射边界条件waveTransmissive。但后来我好像把上边界定义为了壁面计算的。
-
@李东岳 是的李老师,都是涡量输运方程,但不是一个实验室哈哈哈
-
@尚善若水 后处理可以直接导出涡度,但是我的计算中涡度是个中间量,所以我写出了计算公式,再把它加入到求解器中了
-
@李东岳 哈哈是的,李老师。我自己提前写了一下,担心不准确,来请教一下各位。
(1)体积膨胀效应中的涡度omega是不是应该定义为fvc::curl(U)
呀?
(2)斜压效应中的原式是▽ρ×▽p
(密度梯度和压力梯度),fvc::grad(rho) ^ fvc::grad(p)
中间是叉乘吧?
(3)粘性效应中的粘性应力张量tau是不是应该定义为tau=-muEff()*dev(twoSymm(fvc::grad(U_)))
呢?
关于涡量和粘性应力张量的定义是要在求解器中重新声明吧?关于上述三个效应的量化,我不确定能否在计算结果中直接导出呢?还请李老师指点。 -
如下图红框所示,我想在OpenFOAM中输出这三项所代表的场数据,请问该如何实现呢?跪求
-
如题,请教各位老师同学,我在LES模拟后,对数据进行后处理时,发现一个问题:定压比热的场值经过平均后被缩小了,但是其他场的平均值没有这个问题,请问各位这是怎么回事呢?平均后的峰值(8000)明显低于真实值(120000),如图所示:
-
求问各位,以往的OpenFOAM Workshop 的论文集在哪可以找到?在这个网站里没有呢?https://openfoamworkshop.org/previous-workshops/
我想找一个OpenFOAM Workshop 2016年的一些报告,请问有没有知晓如何获取的? -
@cccrrryyy 感谢感谢!给我指了一个好的修改方向!我做过RANS模型的对比,确实影响不大,但我看文献里有的是Smagorinsky算的很好,比如文献1,有的就是算的误差比较大,比如文献2,一头雾水 还是得自己测试呀!
-
@cccrrryyy 您说的这个方法我在文献也看过,不过跟我这个方法完全不同,就没太关注。
我也按照您的方式调整一下特征长度,看看效果怎么样。特征长度我现在给的是管径的10%。您调整后与实验对的怎么样?
我试了of自带的turbulentInlet,效果不好,别的我也没再找过。我看我这个方向的文献大多基于inflowGenerator。我现在觉得问题有两个:一个是湍流强度,就是雷诺应力影响比较大;再一个就是LES模型,我想试试动态Smagorinsky呢,无奈没找到extend版本的,也没再做修改。 -
@李东岳 李老师,这个是超临界压力下的射流,工质是氮气。具体工况如下:
-
@cccrrryyy 您好!我算的是可压缩射流,物性变化比较剧烈,具有高密度梯度。射流与环境存在传热过程。
蓝色的线,射流的液核长度比较短,跟实验差别比较大,其他区域还算理想。我正在测试一个Prt=1.0的Case,看看结果。
您说的对!是给的雷诺应力,雷诺应力计算的时候是(UI)^2,这个I应该指的湍流强度吧? -
最近做了一些算例测试,更新一下结果,欢迎各位提出宝贵意见!
入口脉动通过InflowGenerator产生,主要修改了入口的湍流强度和湍流普朗特数。有关残差问题,我尝试缩小至1e-9,结果并没有变化。
另外,图中的红色曲线是两篇不同文献的结果,同样使用Smag模型,看起来差距很大。这里附上参考文献2的一些模拟结果,感觉是Smagorinsky模型不太合适呢。
240万网格模拟结果:
330万网格模拟结果:
参考文献2的模拟结果(SM:Smagorinsky;WL:Wale;VR:Verman;SMD:Dynamic Smagorinsky):
-
东岳老师,各位同行,我近期也一直在修改参数重新计算,我抽空整理一下,放到论坛,还望各位大佬多指点。
-
关于OpenFOAM中在采用K-Omega SST模型时,使用/不使用壁面函数的Y+及边界设置问题,可以参考这个链接链接文本,解释的很详细。
-
@李东岳 谢谢李老师,我把残差设置的小一些试试。
-
@李东岳 我设置的残差是1e-6,没太关注每个时间步是否收敛呢,李老师,需要把PIMPLE外循环次数增多吗?
-
@李东岳 李老师,LES结果是我取的时间平均,不过取的间隔时间比较少,曲线有些震荡。您清楚可能哪个环节出问题了吗
-
@李东岳 李老师,那造成我的LES跟实验数据差别的问题应该不是它的原因。求个大佬能给分析一下问题所在
-
@李东岳 李老师,是液氮。
另外,我看另外一个相关帖子中解释说:有效热扩散率alphaEff=alphat+alpha,其中alphat=nut/prt,这里面的nut在RANS中是湍流粘度,可是在LES中是亚格子粘度,会非常小,这就会导致大涡模拟计算alphat非常小。我看了一下自己的RANS和LES中的热扩散率,RANS的alphat的确比LES的alphaSgs高了几个量级,而且跟网格密度也有关。这个问题应该怎么解决呢?难道LES模拟都会遇到这个问题吗? -
如题,我正在算一个自由射流,LES模型采用的Smagorinsky,同时计算了RANS,几何模型就是一个三维圆柱。LES的计算结果显示,在过渡区,射流轴向中心线的密度值总是比实验值高,如下图所示。麻烦各位大佬,可以帮忙分析一下可能的原因吗?
-
@李东岳 多谢李老师提醒,我测试一下这样组合会不会好一些!
-
2D几何模型如图所示,对于图中的对称轴,在OpenFOAM中应该设置成什么边界呢?我尝试了以下方案:
(1)对称轴设置为symmetryPlane,前后面为empty,但是计算出的结果有问题,应该是物理场受到了symmetryPlane这个边界的影响。
(2)对称轴设置为empty,前后面设置为wedge,但是当我的喷嘴处(Adiabatic/Isothermal)上边界设置为等温边界时(研究壁面传热对射流影响,给的固定温度),计算结果仍然存在问题。
请问有么有清楚地大神指导一下这个喷嘴的上下壁面具体应该怎么设置?谢谢啦!
-
@李东岳 李老师,信赖域算法也是一种搜索最优解的方法。目前我用的迭代方式如下,其中h_fix是固定值,h_metlab和Cp需要根据更新的温度进行重新计算:
F = (h_fix - h_metlab) / h_fix; // Objective function to be minimized scalar eps = Foam::mag(F*F); // Residuum if (eps <= 1e-3) { converged = true; TEQ = T; break; } //如果残差不被满足,更新temp,重新计算H_mix else { // Update temperature scalar alpha = 0.1; scalar delta_h = (h_metlab - h_fix); scalar dT = - delta_h / Cp; scalar C = 1.0 / (1 + Foam::mag(dT * alpha)); T = T + dT * C; }
-
如题,最近遇到一个通过更新焓求解温度的迭代问题,使用了牛顿迭代等方式,因为热物性变化比较剧烈,导致牛顿迭代非常不稳定,总是在某个工况下不收敛,看文献中使用了Intel的MKL库中的信赖域算法,没有出现稳定性问题,请问有没有了解如何在OpenFOAM中如何应用这个算法的,求指点!
-
@bestucan 非常感谢老师您的耐心解答!受教了!再次感谢!
-
@bestucan 老师,我好像明白您的意思了。通过将我的求解器求出的压力和温度传递给函数
FXc(float temperature, pressure)
,然后在函数中可以得到具体的Xc或者Yc值,接下来进行对比,再按照流程就可以了。因为求解过程中每个网格单元的温度和压力都是在不断变化的,那我的拟合函数应该是一个曲面的数据吧? -
@bestucan 感谢老师的耐心解答!第一帖的两张图的确是通过MATLAB计算出的散点数据然后进行拟合得到的。我可以通过MATLAB计算出这些原始的散点数据,按照您说的方式将其拟合为函数表达式。您说的“把求解器里的Xc出现的位置换成拟合出的函数”这个地方我不太明白,拟合好的函数应该怎么用OpenFOAM识别呢?
-
@bestucan 老师,您好,求解器是已经编译好的了,只是需要加入这个热力学状态判定过程。您分析的是对的,您说的将散点数据输入到OpenFOAM进行插值,可以具体一点吗?
我的具体问题是:在数值模拟计算过程中,根据压力从混合物临界温度脉谱中读取混合物临界温度;根据压力和温度在饱和液相正庚烷质量分数脉谱中读取液相中正庚烷质量分数;根据压力和温度在饱和气相正庚烷质量分数脉谱中读取气相中正庚烷质量分数,然后在OpenFoam 中实现正庚烷-氮气混合物热力学状态的判定过程。最后,根据判定结果可获得混合物中气相、液相和超临界相的质量分数,用于进一步混合物热力参数的计算。
判定流程如下:
-
@bestucan 感谢老师的解答。是这样的,每执行一个时间步,我都需要调用这个脉谱数据来对比计算域的每个单元的数据,然后对于需要部分单元的数据采用两相参数进行重新求解。不知道您说的这种方式可不可以实现,您可以再具体给一些提示吗?
-
最近看一篇博士论文,提到“利用 C++程序语言将正庚烷-氮气混合物气液平衡数据库转化为 OpenFoam 可以调用和读取的脉谱数据”,请问各位,这个应该怎么实现呢,给个提示也行。谢谢啦~
脉谱数据如下:
-
@流体中二少年 上面回复中提到了,你可以取变量p的最大值,然后进行比较,即
scalar pmax = max(p2).value(); if ( pmax < scalar(5.0) )
或者你将volScalarField无量纲化,应该也可以比较。
-
@李东岳 感谢李老师!这样写是更有深度吗?
-
请问各位大佬,最近看一篇论文中提到“without sub-grid scale modelling for turbulence or heat/species diffusion”,这应该怎么理解呢,是指层流模型吗?
-
@chien 解决了就好
-
@chien 还可将 volScalarField无量纲化
p2.dimensions().reset(dimless)
-
@chien 你好!抱歉回复晚了。可以用最值替换来解决,例如p2 < 0,即体标量场p2的最大值小于0,这样可以编译通过,但是是否具体可行我没有做测试呢。
-
@李东岳 谢李老师指导,我来拜读一下您推荐的大作
-
@stardust 感谢分享!
-
@zhouxu 据我了解,隐式大涡模拟没有附加的亚格子模型,更多的依赖于数值格式。你说的是不是Implicit filtering LES?
-
如题,在OpenFOAM中存在隐式大涡模拟相关的代码或者模型吗?刚开始接触,不太熟悉,还望各位同行指点或者推荐一下相关内容
-
@tens 再次感谢您!重新尝试了一下,可以输出了!
-
@李东岳 是的老师,我用的是foam-extend版本,数据放到probe文件中了。可以实现对应时间的数据输出!
-
@tens 感谢指点!计算过程中添加probe导出的数据存放在了输出的时间步文件夹中,能否实现直接导出随时间的变化呢?probe可以添加希望输出的时间点吗?我没有找到呢
各位CFDer,请问速度的依时类边界条件与人工湍流如何同时施加呢?
各位CFDer,请问速度的依时类边界条件与人工湍流如何同时施加呢?
各位CFDer,请问速度的依时类边界条件与人工湍流如何同时施加呢?
各位CFDer,请问速度的依时类边界条件与人工湍流如何同时施加呢?
各位CFDer,请问速度的依时类边界条件与人工湍流如何同时施加呢?
求助大家:图中所示的边界viscous wall和farField在of中应该怎么设置?
求助大家:图中所示的边界viscous wall和farField在of中应该怎么设置?
请教各位大佬,下图公式中的三项怎么在OpenFOAM或者其他后处理软件中导出呢?
请教各位大佬,下图公式中的三项怎么在OpenFOAM或者其他后处理软件中导出呢?
请教各位大佬,下图公式中的三项怎么在OpenFOAM或者其他后处理软件中导出呢?
请教各位大佬,下图公式中的三项怎么在OpenFOAM或者其他后处理软件中导出呢?
请问LES模拟后平均值与真实值不对应是什么原因导致?
请问以往的OpenFOAM Workshop 的论文集在哪可以找到?
LES模拟求助
LES模拟求助
LES模拟求助
LES模拟求助
LES模拟求助
LES模拟求助
关于OpenFOAM中K-Omega SST模型使用壁面函数时的Y+问题
LES模拟求助
LES模拟求助
LES模拟求助
LES模拟求助
LES模拟求助
LES模拟求助
请教各位一个关于轴对称边界的设置问题
请教各位一个关于轴对称边界的设置问题
OpenFOAM在并行运算结束后 能否并行reconstructPar?或者有没有其他方法能够加快reconstructPar?
请教各位关于Intel的MKL库中的函数在OpenFOAM中的使用问题
请教各位关于Intel的MKL库中的函数在OpenFOAM中的使用问题
求问:OpenFOAM如何读取脉谱数据呢?
求问:OpenFOAM如何读取脉谱数据呢?
求问:OpenFOAM如何读取脉谱数据呢?
求问:OpenFOAM如何读取脉谱数据呢?
求问:OpenFOAM如何读取脉谱数据呢?
求问:OpenFOAM如何读取脉谱数据呢?
编译新的湍流模型遇到no match for ‘operator<’
求助“without sub-grid scale modelling for turbulence...”怎么理解?
求助“without sub-grid scale modelling for turbulence...”怎么理解?
编译新的湍流模型遇到no match for ‘operator<’
编译新的湍流模型遇到no match for ‘operator<’
编译新的湍流模型遇到no match for ‘operator<’
请教大家一个问题:OpenFOAM中有隐式大涡模拟(ILES)相关的内容吗?
请教大家一个问题:OpenFOAM中有隐式大涡模拟(ILES)相关的内容吗?
请教大家一个问题:OpenFOAM中有隐式大涡模拟(ILES)相关的内容吗?
请教大家一个问题:OpenFOAM中有隐式大涡模拟(ILES)相关的内容吗?
OpenFOAM后处理怎么提取某点随时间变化的数据
OpenFOAM后处理怎么提取某点随时间变化的数据
OpenFOAM后处理怎么提取某点随时间变化的数据