@李东岳
是的,我确实还需要多加学习!谢谢李老师
是嘿嘿牛啊
帖子
-
-
@李东岳
内部我就想模拟为一个装满水的容器,弧面我视作固体壁面,然后通过映射把力映射到结构求解器中,在结构求解器中给予材料属性来计算结构变形,然后把变形映射回OpenFOAM。底部边界我也看作壁面,就如同水槽中底部边界那样设置。我现在把边界条件设置成这样alpha.water p_rgh pointDisplacement U, 这样的边界条件无法计算,报错是这样的log.interFoam 。 -
@李东岳
李老师,如果我采用分区域耦合的话,整个模拟分成Fluid-outer, Solid, Fluid-inner三部分。那我的Fluid-inner的边界条件该如何给定?Fluid-inner有两个边界弧(interface)和底部边界。弧(interface)是需要用来耦合交换数据的。 -
@李东岳
我用interfoam计算,算例会很快发散。查看结果发现结构物附近的p_rgh值有异常的变化, 而且结构物上的液面有流速,位置正是液面锯齿处。但是我用wave2foam来计算的话,计算不会发散,可以计算完20s.我把两个算例放在了github(https://github.com/NUGUY2022/OpenFOAM_revise.git). 请李老师给点意见 -
@李东岳
之前计算的单一流体域的是计算的弹性板
您的意思是让我只计算流体域的上半部分看看吗?把圆弧结构内部填实吗? -
@李东岳
流体域只有一个的时候是可以很好的计算的,我之前一直计算的单一流体域。 -
@李东岳
算变形
我想知道我的流体求解器设置是不是有问题。 -
@李东岳
李老师,我这个网格有大大的问题啊。结构物上方的初始液面呈现锯齿状, 然后我单独计算Openfoam可以计算,但是流固耦合的时候会发散。这应该是网格导致的问题把, 李老师可以给我一些建议来解决这个问题吗?
我的算例放在了github(https://github.com/NUGUY2022/OpenFOAM_Calculix_membrane.git) -
@李东岳
感谢李老师的回复。
两个区域是分开的,我用其它的结构求解器来计算结构的变形,然后传递到这两个区域,但流体的两个区域不会互相影响。
我想做二维的模拟,于是我把前后面设置为empty,这时checkMesh出现下面的问题。
log.checkMesh
该如何解决这个问题? -
Fluent3DMeshToFoam将fluent网格转换成openfoam网格时出现这样一条内容“ <<Writing 158278 conflicting points to set nonManifoldPoints”,想请问一下大家nonManifoldPoints会对计算产生影响吗?整个网格checkMesh显示的是Mesh OK。网格如图所示
另外还想问一下,我中间结构将流体域分成了内外两部分,网格这样划分是否合适?
polyMesh和fluent.msh也都上传了,希望大家可以提供一些建议。fluent.rar polyMesh.rar -
@李东岳 看到这些质量不好的网格点,好想像hypermesh一样去手移
-
请教一下各位老师,利用snappyHexMesh贴合STL模型时,有什么办法可以让我STL模型的边界网格点(在不改变细化等级和边界网格总数的前提下)等距分布吗?如图所示。
-
@李东岳 正在学习,谢谢李老师!
-
@李东岳 明白了,谢谢李老师!
-
@李东岳 李老师,下面那一行代码alpha()*rho()*betaStar_*omegaInf_*kInf_和alpha()*rho()*betaStar_*omega_()*k有什么区别吗?
我看上面代码中有对omegaInf_和kInf_定义,kInf_ ( dimensioned<scalar>::lookupOrAddToDict ( "kInf", this->coeffDict_, k_.dimensions(), 0 ) ), omegaInf_ ( dimensioned<scalar>::lookupOrAddToDict ( "omegaInf", this->coeffDict_, omega_.dimensions(), 0 ) )
那是否可以理解成omegaInf_和kInf_就是0?不明白alpha()*rho()*betaStar_*omegaInf_*kInf_这串代码的意义,请教一下老师,谢谢!
-
@李东岳 公式没问题,只是左边红框内代码fvm::Sp(alpha()*rho()*epsilonByk(F1, tgradU()), k_)不知道是什么意思,是公式里面的哪一项?刚开始学of编程,求李老师指点一下,十分感谢!
-
@李东岳 李老师您好,我目前用的版本是v1812
-
编程小白请教各位老师两个不成熟的问题。。。
图一中湍流模型方程左右都乘了一个相分数,那如果在两相流模型中,都是空气时,相分数为0,并且相分数的时间导数和空间导数也是0,这方程怎么解?
图二中左边的代码和右边的公式感觉不匹配,求老师们指点一下,十分感谢!
-
最近在做风浪相互作用的模拟时遇到一个问题:在气液剪切过程中,湍流的研究是非常重要的,我在加入湍流模型之后,计算结果耗散特别严重,但是没有气体作用(nowind)时,SSTk-w湍流模型能够较准确的模拟,从下面的水位图中也可以看出来。并且湍动能k和湍流粘度nut在出口outlet逐渐向计算域内部扩散导致波高衰减(见下图),但是不知道具体原因。单纯的平底地形结果也没问题,但是遇到到这种斜坡地形就有衰减。
目前我已经尝试过of自带的几乎所有RANS类的湍流模型,计算结果均耗散严重(以SSTk-w为例),也尝试过对边界条件和壁面边界层的各种处理方式,还是效果不理想,想请教一下各位大佬,想解决这个问题,我还应该从哪些地方入手?
二维模型如图:
水位图:
不同时刻湍动能k分布:
不同时刻湍流粘度nut分布:
-
@李东岳 是的李老师,我理解成了这个边界用在入口时相当于固定值,由inletValue指定,用在出口相当于零梯度,应该是我理解错了。。。
-
@李东岳 好的,谢谢李老师!
-
@李东岳 李老师您好,这是SA模型里面nuTilda的边界设置,除了入口这里取值不同,其他设置一模一样,虽然取值比较极端,但是确实得到的液面结果有很大不同。
-
@李东岳 想请教一下老师,在v1812版本中,inletValue值一定的情况下,如果value取不同值,结果有一定差异,那value值有什么含义呢?
-
@李东岳 hyperMesh的ployMesh里面多4个Zones文件,然后faces、points里面的顺序也不同
-
@队长别开枪 提问中的图一是示意图,不是实际网格分布,实际网格如图,上图为blockMesh,下图为hyperMesh
-
@李东岳 是的很奇怪,我是在一次使用局部加密网格中发现这个问题的,当时利用两种方法生成了两种相似网格结构,计算结果有较大差异。这次为了让网格完全一样,使用了最简单的矩形计算域和均匀结构网格,结果竟然还有差距...
-
@李东岳 非稳态波浪场
-
一模一样的网格,只是生成工具不同,算出来的结果为什么会有差异?
-
分别用hyperMesh和blockMesh画一个长20米,高1米,宽(empty)0.006米的长方形二维计算域,长方向3333层网格,高方向167层网格,宽方向1层网格,算出来的结果为什么会有差异?请教一下各位老师。
-
@tens 明白了,十分感谢!
-
我在书中看到AlphaCo的定义,想请教一下各位老师,图中的pos指的什么意思?
-
@cccrrryyy过去的研究中,我还未发现过层流发散的情况,但是会出现湍流强的时候,层流模型计算效率低于湍流模型。湍流弱的时候,层流模型计算效率高于湍流模型。
-
@cccrrryyy 层流模型应该接近于直接数值模拟(DNS),在有限体积法中只是精度上不去,波浪传播中湍流较小,只有在破碎的时候有较强的湍动特征,可能湍流强度还没大到层流模型会发散的地步。
-
@tens 你好,暂时不能确定是不是壁面处理的问题,因为低风速下两个模型的结果基本一致。
-
最近做了高风速场下的波浪变形规律研究,发现标准k-e湍流模型和SSTk-w湍流模型算出来的结果相差甚远,经过验证已经确定SSTk-w的计算结果较为准确,想请教一下各位老师,标准k-e模型是不是不适合高流速场?还是说不适合我研究的流动情况?十分感谢!
图1为层流模型云图,图2为SSTk-w模型云图,图3为k-e模型云图,(图1、图2、图3为同时刻同位置)图4为SSTk-w和k-e模型算出来的水位时变对比图。
-
@zhouxu 横坐标时间T 纵坐标力F
-
这是垂向力
这是水平力,没有出现垂向力的阶梯
-
@李东岳 非常感谢东岳老师的回复,这个问题出现在均匀网格,有的时候甚至同样的位置,但是气域的速度不同,他最后采样的范围也不一样。。。
-
如图,我用line函数监测水位,检测范围给的0-1m,但是输出的数据只到0.6m,有人遇到过这种情况吗?
-
测力时出现很强烈的震荡 想请教下是什么原因导致的 黑色的是想要验证的图(没有震荡),红色的是我计算的。用的Openfoam自带的测力工具。
-
刚盖假定在OpenFOAM中要使用什那种边界条件来实现。
-
@是嘿嘿牛啊 按照我的理解,这样分的size不应该都是10吗?为什么size会变来变去呢?求大佬指点一下,拜托啦!
-
@zhouxu 感谢回答,我还是不太明白,您有参考文档吗,如何在openfoam中表现
-
@李东岳 流体是运动的 ,K-H不稳定也是有的,只是想在初始时刻给这个混合区指定一个从ρ1→ρ2密度分布。
-
如图上下两层流体密度不同,想要使得中间混合区的密度根据某一方程变化,应该如何实现
-
@bestucan 好的,我去按照这个思路研究研究,非常感谢!
-
上次有大佬建议我做库朗数的收敛性验证,我回去尝试发现,最大库朗数maxCo在设置很小的时候模型会发散,不是很明白,求各位老师解释一下,十分感谢!
如图,maxCo设置0.05时有结果,但是设置为0.01时,模型算了2秒多后就会发散。
-
@yingdong 明白了,谢谢老师。
nonManifoldPoints与多域网格问题
nonManifoldPoints与多域网格问题
nonManifoldPoints与多域网格问题
nonManifoldPoints与多域网格问题
nonManifoldPoints与多域网格问题
nonManifoldPoints与多域网格问题
nonManifoldPoints与多域网格问题
nonManifoldPoints与多域网格问题
nonManifoldPoints与多域网格问题
nonManifoldPoints与多域网格问题
关于snappyHexMesh网格划分问题
关于snappyHexMesh网格划分问题
湍流模型方程和代码问题
湍流模型方程和代码问题
湍流模型方程和代码问题
湍流模型方程和代码问题
湍流模型方程和代码问题
湍流模型方程和代码问题
气液剪切流的湍流耗散问题
inletOutlet边界条件
inletOutlet边界条件
inletOutlet边界条件
inletOutlet边界条件
hyperMesh与blockMesh的网格生成问题
hyperMesh与blockMesh的网格生成问题
hyperMesh与blockMesh的网格生成问题
hyperMesh与blockMesh的网格生成问题
hyperMesh与blockMesh的网格生成问题
hyperMesh与blockMesh的网格生成问题
AlphaCo问题
AlphaCo问题
高速度场下的标准k-e和SSTk-w湍流模型问题
高速度场下的标准k-e和SSTk-w湍流模型问题
高速度场下的标准k-e和SSTk-w湍流模型问题
高速度场下的标准k-e和SSTk-w湍流模型问题
求助:为什么openfoam测出的垂向力会出现这种阶梯状
求助:为什么openfoam测出的垂向力会出现这种阶梯状
interFoam水位监测
interFoam水位监测
测受力的震荡问题
OpenFOAM中刚盖假定
topoSet的faceSet问题
topoSet的faceSet问题
topoSet的faceSet问题
如何实现两层流体密度过渡区的随某一方程变化
如何实现两层流体密度过渡区的随某一方程变化
如何实现两层流体密度过渡区的随某一方程变化
两相流模型库朗数问题(二)
两相流模型库朗数问题(二)
关于两相流模型中的库朗数取值问题