@五好青年 非常非常感谢您!我的邮箱是dydsj_168@163.com
DY大世界
帖子
-
-
@wwzhao 老师您好,目前我也在用重叠网格模拟固体入水,但是我的水面是静止的没有波浪。我也出现了类似的问题,就是固体还没有接触到水面的时候,水面就已经发生了巨大的扰动,不知道为什么固体还没有接触到水面的时候就会对水体产生这么大的非正常影响。
-
@五好青年 您好,我最近也在学习利用重叠网格模拟自由入水,可以看看您的算例吗?
-
@迟铁 好的,非常感谢大佬的指点!
-
有大佬帮我看看这个问题吗
-
@李东岳 啊,这。。。
-
@队长别开枪 好的,谢谢队长!盼望师姐的答疑解惑
-
利用OF内的overInterDyMFoam,基于重叠网格模拟二维方形结构物以恒定的速度垂直入水。
数值模型包含了两套网格,一是固定不动的背景网格,二是做恒速垂直下落运动的组件网格。
目前的问题是,组件网格在下落的时候,整个组件网格里面的air、water一直都与初始时刻的相同,无任何变化,如下图所示。感觉是背景网格里的air、 water与组件网格里的air、water是相互独立的两部分。
计算四五个时间步后就发散了。恳请老师们能予以解答,感激不尽!
-
@kennyadapter 您好,打扰您了。请问您的问题解决了吗,我现在也需要模拟结构入水的案例,请问您的dynamicMeshDict是怎么写的呢,能否看一下您的dynamicMeshDict字典?
-
例如下图所示的波浪-圆柱体耦合计算
根据理论波浪力方程,圆柱所受到的波浪力(水平方向、垂直方向)为
Fy表达式中的最后一项ρgV0为圆柱体的静水浮力(initial buoyancy in still water);
请问:OF输出的Fy(垂直方向上圆柱受到的Force)里面包含ρgV0这部分吗?想要获得圆柱受到的垂直方向上的波浪力,是不是需要从OF获得的计算结果中手动减去ρgV0这部分?
-
@yuno 好的,感谢回复指导!
-
@东岳 哈哈,明白,谢谢李老师~
-
@东岳 好的,感谢李老师的回复!之前老是想不太明白同一个算例,只是用了加速运算,老想着应该是仅效率提高,结果应该相同才对。这样看结果就是会略有差异的。
-
波浪-固体耦合方面的算例,开了8核的并行计算,计算效率很高,用时大概是大概是单核的1/3。但是计算的结果与单核的始终不能完全吻合,总存在一点差异,无论是使用scotch还hierarchical,计算结果都与单核的有差异。两者之间是必然会有误差吗?可以通过什么方法消除?请大神们赐教!!
图1 固体受到的水平波浪作用力。
-
@闻久STU 收到~感谢大神指点!
-
@xpqiu 非常感谢大佬的回答,那不用extrudeMesh命令所得到的确实是不正确的二维网格。而且我对比了使用extrudeMesh与不使用extrudeMesh的计算案例,得到的结果也是不完全吻合。
-
@东岳 感谢李老师回复!
例如这篇文章没有提到extrudeMesh,说的是将三维转换成二维只需要将边界指定为empty。
现在我的案例是关于波浪-固体耦合的,如果不用extrudeMesh还好,一使用计算就发散,不知道使用extrudeMesh后还有什么要求,比如setFields什么的。
-
@xpqiu 你好大神,我想再问您一下,这两天我看了一些文献,有很多文章都是在使用snappyHexMesh后,没有再用extrudeMesh命令,而是直接把前后面的type设置成empty进行计算。但是仿真结果与试验的对比也挺好。
我想问一下这种使用snappyHexMesh后通过直接设置front和back的type为empty的方法,能否表示二维模型,这种网格是算二维的还是三维的。因为他们的计算结果精度挺好,是否表明这种方法对于二维模型是可行的?
-
@史浩 好的,很感谢大佬的回复,我对比看看结果!
-
@xpqiu 好的,谢谢老师的指点回复!
-
老师不好意思又打扰您了,按照您的指点现在划分的网格在empty方向上处处都有只有一层网格了;
但是不知道为什么现在网格是二维的了,但计算确发散了。利用extrudeMesh将三维网格转换为二维网格需要注意什么吗,现在不管怎么调试都是出现下图所示。
-
@xpqiu 收到!感谢老师的指导!巨谢!
-
@东岳
谢谢李老师的回复!!这两种方法都是在数值计算前把front和back两个面(sides)的边界条件设置为empty后才进行计算的,边界设置成empty就是表明模型是二维的了。
在网格划分的过程中,我在blockMesh中将这两个面的边界设置为empty,
两种方法的区别就是:
方法一是:在使用snappyHexMesh前先将这两个面的边界改为wall,这样在snappyHexMesh过程中就没有图2的警示了,而且网格划分消耗的时间很少。在使用完snappyHexMesh后再把边界由wall改为empty进行数值计算。
方法二是:在使用snappyHexMesh时不修改front和back的边界条件,两个面的type依然是empty,这样就会出现图二的警示,而且计算时间很长。
但是这两种方法得到的cells和faces是相同的,就是points稍有不同,不过数值计算结果也是略有差异。所以实际最好采用哪种方法来划分网格进行二维数值计算呢。
方法一中我把front和back的边界改为wall后再使用snappyHexMesh,是考虑snappyHexMesh只能划分三维网格,所以边界的type如果是empty应该就不合适了,所以系统也会出现警示。但是仿真的是二维模型,所以在网格划分完成后再将这两个面的type改为empty。
-
利用二维水槽来研究波浪-圆柱固体耦合,水槽的前后两面(sides)的边界条件为empty,以下这两种网格划分方式有什么区别?:
方法一:
(1)先使用blockMesh划分背景网格;
(2)在生成的‘’polyMesh‘’文件夹下的‘’boundary''文件中,将sides的type由empty修改为wall;
(3)再执行snappyHexMesh生成带有圆柱固体的网格;
(4)再将‘’polyMesh‘’文件夹下的‘’boundary''文件中 sides的type由wall修改为empty,从而开展后续的二维数值计算。这种网格生成的方法只需13.33s,如下图所示
方法二:
(1)先使用blockMesh划分背景网格;
(2)直接执行snappyHexMesh生成带有圆柱固体的网格,以开展后续的二维数值计算。这种网格生成的方法在生成过程中会有这样的提示
网格生成时间需要417秒
这两种方法所需的时间差别很大,生成网格的cells和faces数相同,points数略有不同,但是计算结果也不完全相同,有差别。
请老师们指点,两种网格划分的区别,二维模型的流固耦合计算,应该选择哪种方法?
-
@东岳 感谢李老师的解答
-
@东岳 想看一下加速度云图应该怎么设置
-
在controldict文件中添加了grad(U)字典了,想看加速度场的云图,这里grad(U)有8个magnitude,这是表示什么意思,想看加速度应该选择哪一个?
-
现在很多研究采用的是waves2foam造波工具箱集成到OF中。我用的是OpenFoam1906版本,这套版本中有自带的造波工具箱即原来的ihFoam,所以我就不想再用waves2foam了。但是这个版本中造波后如何提取任意一点处的波高随时间的变化呢,有没有用过的大神。
目前我在controldic文件中添加了一段代码能追踪自由波面,但是提取的文件也是vtk的格式无法打开,恳请会的大佬不吝赐教!感激不尽!
-
@cccrrryyy 哦哦,好的。感谢大神的解答!
-
关于二维数值波浪水槽的流固耦合问题,如下图所示的波浪与半浸没的二维圆柱体相互作用,圆柱体的半径为a.
二维数值水槽的圆柱体长度L1=0.04m,对应的一个三维水池试验的圆柱的长度L2=3m,(数值和试验除了圆柱长度这一参数不同,其余均相同)。我想用三维的试验结果验证二维的数值结果,请问采用如下的方法是否是可行?
数值计算和水池试验除了圆柱的长度不同,其余所有参数均相同,请问一下如果数值计算和水池试验最终得到的“无量纲化后的波浪载荷振幅”基本一致,是否可能证明数值仿真的正确性。
-
@东岳 感谢李老师回复!圆柱周围的网格细化我明白,但是在整个水槽的自由液面处网格细化也是用snappyHexMesh做的吗,而且这里的自由液面区域和圆柱周围区域的衔接也很好,看上去像是一个整体加密区域。但是我总感觉自由液面处的加密像是通过blockMesh实现的,我在这里还是有点困惑,李老师您能详细的指点我一下吗。
-
关于二维数值波浪水槽,波浪-圆柱耦合问题,二维流固耦合简图如下所示:
我想请教各位大神,下图的这种网格是怎么样得到的?下图的网格是在整个水槽的自由液面处和圆柱体的周围进行网格加密,而且也能使加密区域和稀疏区域的网格都是正方形,从加密区域到稀疏区域网格的过渡也比较好,这种网格是如何处理才能得到,请不吝赐教!感谢感谢!
-
@cccrrryyy 好的好的,十分感谢老师的解答!
-
@cccrrryyy 感谢老师的解答,很受教!
另外我还有两个疑问;
一个是在这种波槽内,关于雷诺数里面的特征速度应该如何取值,是取波浪的波速吗,感觉只有波速是已知的,但有些人都说采用平均速度,我不太明白波浪运动里的平均速度是指的什么。
二是,我看文献很多学者在研究这种波浪绕流的问题都选择层流模型而没有用湍流,按理说像这种问题应该会发生波浪破碎,搞不明白为什么使用层流模型,选择层流是否合适,感觉这样就不能捕捉到波浪绕流时发生的扰动了。以下是列举的两篇使用层流研究此类问题的文章。谢谢老师!
-
请教一下:二维数值波槽里水流动时的雷诺数应该怎么求,流速U是指波浪水质点的流速吗?特征尺寸应该为多少?
看了如下一篇文章有一点不太明白,就是为什么作者把雷诺数公式里的特征尺寸用平板的高度D代替,雷诺数和水槽里的固体有关系吗?关于二维下的数值水槽,水槽内有无固体对雷诺数有影响吗?并且特征尺寸和特征速度应该取什么?
请大家赐教!感谢!
-
@Cp_Zhao 明白了,感谢指点!!
-
请问一下图里面黄色矩形框里面说的这个方法,banana method找所有可用选项,这个方法是什么意思?
-
@队长别开枪 队长您好,您可否帮我解答一下这个帖子里问的问题,就是为什么snappyHexMesh画网格始终无法移除固体,感谢队长!
http://www.cfd-china.com/topic/3408/snappyhexmesh为什么总是无法移除固体 -
画网格过程中出现这个提示是什么意思
-
@东岳 李老师您好!几何模型应该是封闭的,因为solidworks画图草图不封闭的话也不能拉伸成实体,请问李老师还有没有可能是哪里出问题?谢谢李老师!
-
@东岳 谢谢李老师回复,请问这个怎么检查是否封闭或者使之封闭,我这是用solidworks软件画的外部几何模型,导出来以后就是一系列点坐标,不知道应该怎么改。
-
想用snappyHexMesh划分一个二维的流域,流域中有一个矩形薄板,如下图所示(别人文献中的图)
但是为什么我画网格的时候始终无法把这个固体部分从网格中移除,如下图所示,红色区域部分本该有一个被删除的固体,但是现在全部都被画上了网格。
并且locationInMesh也是在固体外部并且不在网格面上,请问大神们问题出在了哪里呢? -
@东岳 李老师,我试了一下,把出口和入口速度设置的速度一样,在开始的时候水池的水就不稳定,在入口处水面上升,在出口处水面下降。得好久才能接近一种平衡,而且导致液面晃荡。
-
@东岳 好的,感谢李老师解答,我再试试。
-
做一个二维两项圆柱绕流,想模仿实验室的造流机,入口以一个恒定的速度1m/s造流,出口处水的高度不变。
但是现在要么就水涨船高,水池里的水越来越多。
要么水就都从出口处流走了。
alpha.warer文件,p_rgh文件,U文件处的出口边界,应该如何设置,才能模拟实验室的造流效果呢?
-
@东岳 感谢李老师!已解决~
-
@队长别开枪 十分感谢队长,按照您的方法问题已解决!
-
编译安装paraview5.6.3时,键入 ./makeParaView一直提示bad variable name
请教一下大家,出错是因为什么原因呢?
OpenFOAM中动网格(或者重叠网格)计算后,怎么提取浮体的位移呢?
OpenFOAM重叠网格
OpenFOAM中动网格(或者重叠网格)计算后,怎么提取浮体的位移呢?
基于OF重叠网格模拟结构物入水的问题
基于OF重叠网格模拟结构物入水的问题
基于OF重叠网格模拟结构物入水的问题
基于OF重叠网格模拟结构物入水的问题
基于OF重叠网格模拟结构物入水的问题
基于OF重叠网格模拟结构物入水的问题
OpenFOAM重叠网格
OF中关于波浪-固体耦合的计算,输出的Force包含哪些力?
使用OpenFoam1906,软件自带的造波工具包造波,如何提取波面曲线(surfaceElevation)
并行计算与单核计算的结果不能完全吻合吗?
并行计算与单核计算的结果不能完全吻合吗?
并行计算与单核计算的结果不能完全吻合吗?
如何看流场中的加速度云图?
关于二维模型的计算,这两种snappyHexMesh方法有什么不同,应该采用哪一个?
关于二维模型的计算,这两种snappyHexMesh方法有什么不同,应该采用哪一个?
关于二维模型的计算,这两种snappyHexMesh方法有什么不同,应该采用哪一个?
关于二维模型的计算,这两种snappyHexMesh方法有什么不同,应该采用哪一个?
关于二维模型的计算,这两种snappyHexMesh方法有什么不同,应该采用哪一个?
关于二维模型的计算,这两种snappyHexMesh方法有什么不同,应该采用哪一个?
关于二维模型的计算,这两种snappyHexMesh方法有什么不同,应该采用哪一个?
关于二维模型的计算,这两种snappyHexMesh方法有什么不同,应该采用哪一个?
关于二维模型的计算,这两种snappyHexMesh方法有什么不同,应该采用哪一个?
关于二维模型的计算,这两种snappyHexMesh方法有什么不同,应该采用哪一个?
如何看流场中的加速度云图?
如何看流场中的加速度云图?
如何看流场中的加速度云图?
OpenFoam1906版本自带的造波工具箱如何提取任意点波高
用三维水池试验验证二维数值试验的问题
用三维水池试验验证二维数值试验的问题
请教如何划分才能得到这种网格?
请教如何划分才能得到这种网格?
二维数值波槽的雷诺数如何求
二维数值波槽的雷诺数如何求
二维数值波槽的雷诺数如何求
banana method是什么意思
banana method是什么意思
安装paraview时,./makeParaView出错
SnappyHexMesh为什么总是无法移除固体
SnappyHexMesh为什么总是无法移除固体
SnappyHexMesh为什么总是无法移除固体
SnappyHexMesh为什么总是无法移除固体
模仿实验室的造流水槽,出口边界怎么设置
模仿实验室的造流水槽,出口边界怎么设置
模仿实验室的造流水槽,出口边界怎么设置
安装paraview时,./makeParaView出错
安装paraview时,./makeParaView出错
安装paraview时,./makeParaView出错