就是你说的,计算得到的结果是一样的。
xiezhuoyu
帖子
-
-
你想在openfoam里面直接计算?有困难,虽然网格是结构化的形式,但是非结构网格求解器仍然不方便找到你想要的信息。
可能还是praview+宏,批量处理n条线,把数据提取出来会简单些。 -
我感觉比较简单的应该是,paraview用plotOverLine,将数据另存为文件,之后在对这条线上的量计算。
-
@李东岳 在 fluent用LES算水翼边界层选择中心差分引起的震荡如何消除 中说:
库朗数减小一些,比如0.3或者0.2左右来尝试纯粹的linear。
这个倒是找到一个很有意思的格式("Blends two schemes based on the local face-based Courant number"),
div(phi,U) Gauss CoBlended 0.01 // Co below which scheme1 is used linear // scheme1 0.05 // Co above which scheme2 is used LUST grad(U); // scheme2
这个格式可以根据Co数自动调整格式,可以实现按区域混合各种格式(虽然不是按权重混合)。
-
-
@心里的小漩涡 在 fluent用LES算水翼边界层选择中心差分引起的震荡如何消除 中说:
有试过二阶中心差分格式中耦合小部分一阶迎风格式吗,比如给个0.02%的一阶迎风格式,之前用fluent算3900的圆柱绕流,用纯二阶中心出现过数值振荡容易发散的情况,改成混合之后收敛性好很多,精度也还行,没有很明显的下降
@李东岳 李老师,这个在OpenFOAM里面能直接设置吗?
-
空的是指真空还是空气?
空气就两相流(入口是water,计算域本来是air),真空就知识盲区了。 -
假设问题中流体不可压缩,左侧有流体流进+1,右侧当然有流体流出-1,瞬时完成。
-
我的意思是
1.我只试过SST k-omega湍流模型,这个没问题。 ()
2.大涡模拟里面没有k方程,也就没这回事了。 -
SST k-omega湍流模型应该可以这么实现
-
按我的理解,如果你滤波器的尺寸跟网格尺寸一样,可以不用显式地使用滤波器,因为有限体积法用的体心值取的就是单元平均(?),相当于做了滤波操作。
如果要加滤波器,滤波器的尺寸应该是根据问题要捕捉的涡的大小来确定的吧。 -
@chapaofan
这里面有矛盾的设置。
每条边上的单元数改为10试试(下图),再慢慢调整。
-
@chapaofan 在 blockmesh生成网格总报错没有相邻面 中说:
hex (14 11 12 13 29 30 31 32) (10 10 1) simpleGrading (1 1 1) //3
上面的顺序错了。
可能是下面这个:
hex (14 11 12 13 32 29 30 31) (10 10 1) simpleGrading (1 1 1) //3 -
请教,symmetry边界,直观上理解是没有界面法向的速度变化吧,没有速度变化,为什么有切应力?
-
浮体为不规则形状,已经计算过浪向角为0°的算例。
现在想要计算其他浪向角,如60°工况,能否直接通过旋转浮体-60°来达到目的,可是旋转浮体,对应的转动惯量也要做相应的更改,不知道OF里面有没有什么设置可以方便地实现旋转浮体的功能并自动修改转动惯量? -
补充一下学长 @臧伟 的内容
2.关于时间步长的问题,考虑到波浪信息在单个时间步长内传播距离不宜太大,一般取T/500左右。以及上文提到的网格尺寸问题,可参考:
查晶晶. 基于OpenFOAM的数值造波与消波模型及其应用[D]. 上海交通大学, 2011.
如果单纯造波消波的话,时间项采用二阶离散格式也有利于抑制波幅衰减的问题。
3.如果是阻尼消波的话,阻尼区域建议为波长1~2倍。 -
@西极若木
邮箱:mechanics_xzy@163.com
微信:xzy11235813
多交流多交流 共同进步 -
@西极若木
我也很新手啊
我说的不一定对,你当作参考吧,我们一起讨论讨论
1、UMean,我以前用这个计算的时候,paraview看速度场,除了自由液面附近,发现大部分区域都是UMean这个速度,所以我以为它是流的速度(风浪流的“流”),不知道对不对;
2、scale和cross scale可以看看源代码里面的注释,应该是你理解的类似松弛区阻尼系数的概念。((4 1) (12 0))里面的4和12应该是之前定义的direction (-1 0 0)即-x方向上的位置(这里我有点迷,给的-x方向?),用来确定阻尼区的范围。
3、OF7造波不太确定了,消波应该是阻尼消波了吧。
OF7的这个波浪算例的话,除了源代码,相关资料看有没有前辈或者同辈或者后辈提供提供 我是只知道这么多了,还不一定对。
你是做振荡浮子的对吧,也是波物相互作用,我们很像啊(我研二),有机会多交流!!我之前用过一段时间的OF7,就是你这个算例,但是很头疼UMean的问题,我给UMean=0,就会出错,现在也不知道为什么,请老铁指点;然后转了OF1912,OF1912也自带造波消波,但是它消波针对浅水问题,有限水深效果也还行,深水无法保证;最后用了现在的waves2Foam+overInterDyMFoam,速度入口造波+阻尼消波+overset。总之,每一个都有不同的困难,希望你加油吧! -
1、UMean好像是流o_o;
2、位置“4”处开始消波,系数是1;位置“12”处消波末端,系数是0;cross scale忘了。。。
3、等大佬来分析
@西极若木 OF的波浪问题可以多交流啊 -
@CFDngu 谢谢建议!!
-
虽然我算例的浮体和你不一样,但是问题很类似:在波浪还没有传播到浮体附近之前,浮体就受力振荡,p_rgh的迭代次数不是递减的趋势,感觉我就是在踩你踩过的坑 @tower
overset: oversetMesh.zip 这个算例用的是OF1906自带的造波消波模块吧,不需要waves2Foam造波消波,@东岳 老师能不能给我们一点建议o_0,看看应该从什么角度着手尝试解决这个问题?谢谢!
也希望有其他前辈能够就这个问题指点一下,缩短小白踩坑的时间,感谢!! -
@tower 老哥,我也遇到了同样的问题
用overset,浮体z方向受力正负交替变化,然后越来越大。
我看结果文件是因为浮体底部的压力变化剧烈,忽正忽负。所以想着加密一下浮体附近的网格,试试能不能解决问题。
可是!我加密网格后,浮体壁面的Y+小于30了(浮体顶部(接触空气的部分)尤其小),其实我的网格也不是很密,网格尺寸大概就是浮体吃水的1/8,而且如果可能的话,我想用壁面函数啊,请问,老哥你遇到这种问题一般如何解决?谢谢
也希望有其他前辈能够指教!谢谢 -
OF-7这边,有interfaceHeight,可以很方便提取某处的液面高度。
OF-v1912这边的话,提取液面高度确实不是很方便,你可能得用你说的方法3。 -
@tower 老哥,不知道你找到问题原因没有? 求分享!我是模拟波物相互作用问题的,对你遇到的这个问题很好奇o_o
摩擦阻力
边界积分
边界积分
边界积分
fluent用LES算水翼边界层选择中心差分引起的震荡如何消除
fluent用LES算水翼边界层选择中心差分引起的震荡如何消除
fluent用LES算水翼边界层选择中心差分引起的震荡如何消除
为什么流体尚未到达的地方有速度?
为什么流体尚未到达的地方有速度?
不同湍流强度下的来流
不同湍流强度下的来流
如何确定大涡模拟非结构化网格滤波器尺寸
blockmesh生成网格总报错没有相邻面
blockmesh生成网格总报错没有相邻面
Stokes流中的圆球绕流的阻力系数,怎么设置才能算准。
OF如何设置并计算不同浪向角情况下的波物相互作用问题
关于waves2foam造波衰减问题
关于OpenFOAM里自带的waves算例的疑惑
关于OpenFOAM里自带的waves算例的疑惑
关于OpenFOAM里自带的waves算例的疑惑
动网格和overset物体速度快速震荡导致崩溃的问题
动网格和overset物体速度快速震荡导致崩溃的问题
动网格和overset物体速度快速震荡导致崩溃的问题
请教如何只提取某处的液面高度时程图,就像实验放置在水槽中的wave height probe那样
动网格和overset物体速度快速震荡导致崩溃的问题