您好,@shirouchen
感谢您的回复,这里补充了显示单元后的图片,应该可以比较清楚地看到网格变形了。
首先,液面部分和顶部边界附近的网格已进行过调整,但并没有改善。此外,由于该案例中的柱体并非刚性体,所以并没有采用其他的如SolidBodyMotion这样的动网格方法。想请问您对这类问题有无其他可用的动网格方法的建议?谢谢。
您好,@shirouchen
感谢您的回复,这里补充了显示单元后的图片,应该可以比较清楚地看到网格变形了。
首先,液面部分和顶部边界附近的网格已进行过调整,但并没有改善。此外,由于该案例中的柱体并非刚性体,所以并没有采用其他的如SolidBodyMotion这样的动网格方法。想请问您对这类问题有无其他可用的动网格方法的建议?谢谢。
向各位大佬请教一个OpenFOAM动网格方面的问题。想要模拟的是一个多相流的流固耦合案例,固体域的结构物底部震动,使得结构物上部出现位移,进而震荡出水波。计算中所采用的动网格方法是displacementLaplacian,扩散方式采用的diffusivity quadratic inverseDistance。在结构物发生位移的过程中,顶部边界出现了如下面两幅图所示的波纹状的网格畸变,随后计算就崩溃了。
目前调整了很多次网格尺寸,但该问题依旧会出现,还请各位大佬多多分享经验,谢谢~
@队长别开枪
感谢前辈,采用几何重构方法后确实解决了这个问题~
@李东岳
感谢东岳老师的回复。目前我们尝试将上边界加高了一定距离,并采用interIsoFoam,上述的气泡和内部的线条就消失了。
在使用olaFlow和interFoam进行波浪模拟时,靠近自由面的区域出现了如图所示的气泡,内部还出现了一条条的线。目前推测可能是湍流模型和离散格式的问题,用k-wSST的话气泡很多,若改为层流后气泡几乎没有了。当时间离散格式改为C-N格式后,不论用哪种湍流模型,气泡都会减少很多,请问有没有大佬遇到过类似问题呢?
可以尝试用functionObject中的surfaces提取物体表面的压力
@李东岳
感谢东岳老师的分享,我先照着该教程试一试
各位大佬好,刚尝试在一台服务器上进行OpenFOAM跨节点的并行计算,测试案例是dambreak,所用版本为OpenFOAM8。我参照该链接中的做法,对/etc路径下的hosts以及/networking路径下的interfaces文件进行了修改:
hosts
interfaces
随后,在dambreak案例目录下创建了machines文件:
machines
执行'mpirun -np 4 --hostfile machines interFoam -parallel'后计算并未开始,而是出现如下报错信息:
请问各位大佬,这可能是哪个文件的设置出了问题呢?
另外,为何lscpu能看到很多NUMA节点,但在ifconfig时只有三个ip(包括下图中的enP1p130s0f1),其他的都是如下面这种情况,这是硬件配置的问题吗?
lscpu.log
问题已解决,在瞬态结构的设置里选中流固耦合界面,在左下方的“数据传输”选项中选择“所有系统耦合数据传递”即可。
各位前辈好,在使用ANSYS Fluent和Mechanical耦合计算时,设置系统耦合的时候,瞬态结构里明明有“增量位移”,但是选择时却提示Invalid option:增量位移。请问这是哪里的设置出了问题呢?
各位前辈好,关于外流场计算中的动压、总压等有些疑问,想和大家请教一下:
对于不可压缩的流体,总压=静压+动压。OF的多相流计算中,初始条件中需要给一个参考压力p_rgh,然后经计算得到总压p。
那么,以一个初始位置在水位以上的结构物(即没有静水压力)为例,当结构物受到波浪冲击时,结构物所受到的流体压力可以认为全部是动压吗?或者说,OF计算出的在自由液面处的压力都是动压?
在计算多相流动网格案例时,计算一段时间后库朗数突然增大并导致计算发散。起初以为是网格变形过大导致的发散,但查看日志文件和结果发现,计算并非是结构物受到波浪冲击后达到最大变形时终止的,而是在逐渐回到原位置时终止的,此时相分数alpha.water超过了1并在急剧增大。这时结构物附近有一定的越流(如下图所示,截图截得不太整齐请大家见谅),且计算停止时(44.2s)的气液交界面显然是有很大问题的。因此我认为应该是相方程计算不收敛进而导致的整个计算终止。
日志文件中相分数在44.076s开始大于1:
最后计算崩溃:
因此,想请问大家是否在多相流、动网格的计算中是否遇到过类似的问题,有什么办法可以解决这种相分数越界的问题吗?
想了一下,虽然网格有变形但整个过程中网格拓扑不变,是否有办法根据该单元的拓扑来提取这个单元的压力时程呢?
以结构物受波浪冲击导致结构物产生变形为例,想要得到距离结构物最近的一层流体网格上某一单元的压力变化情况。在以往的非动网格案例中,我是在system目录下设置一个sampleDict文件并在其中指定提取方式、提取点坐标,但是动网格案例中位置会变,这样的做法显然不行。请问各位是否有比较好用的提取方法?
@guohuiqun
目前,我们是通过渗透试验测出孔隙介质的压差、流速,拟合出I=AU+BU^2中的A和B,然后反算出来a和b的。至于c,OlaFlow作者在其论文中提到该参数的影响很小,故一般是取经验值0.34。
如果是计算时要实时计算边界条件,这个可能符合你的要求:
http://openfoamwiki.net/index.php/Contrib/groovyBC#Usage_of_the_boundary_condition
感谢回复,需要的确实是在计算中实时计算的边界条件,我先了解一下
目前在尝试用质量流率\体积流率边界flowRateInletVelocity来实现,但是发现官方给出的说明中好像只能施加一个常数的质量\体积流率边界,然后根据此流率边界得出此边界上的速度,无法根据边界上的流速得到体积流率(和我要实现的恰巧相反。。)
请问OF中有无根据边界上的流速来增加流体体积的方法呢?
各位好,想请教一个关于多相流的流固耦合问题。具体要模拟的情况就是一个孔隙介质挡板,左侧有水右侧无水,左侧的水进入挡板并从右侧渗出。
目前已在固体求解器中算出了固体边界上的变形和渗流速度并将二者加到了OpenFOAM流体域的相应边界上,请问怎样才能通过该渗流速度在边界上相应的位置补充上因渗流而增加的水呢?
再和各位前辈请教一个相关的问题,在用snappyHexMesh处理网格量比较大的案例时(网格数量200w左右,模型和上图类似,有一些拐角):
@李东岳
东岳前辈好,我也没找到那个帖子中的Use VTK Polyhedron(通过paraview和paraFoam打开都没有),但是有这个Decompose polyhedra,默认是勾选上的。如果将它取消勾选,显示的网格会有一些不同,但是扭曲的部分还是存在,如下面的对比。
另外,即使用了那个帖子中的第2、3点,显示的结果也是没什么区别。
各位前辈好,我在用snappyHexMesh划分外流场网格时遇到了如下图所示的问题,即边界层网格和内部网格之间出现了一块扭曲比较严重的区域。尽管调整了几次snappyHexMeshDict中Explicit feature edge refinement和refinementSurfaces的细化等级,这种问题还是存在。
目前我猜测该问题是同时启用snap和addLayers导致的,但是这种斜面不使用snap功能显然会使网格质量下降,请问有没有什么办法可以处理呢?
(snappyHexMeshDict文件已附上)
snappyHexMeshDict
@samuel-tu
前辈好,请问交界面网格畸形的问题您是否解决了?我也遇到了类似的问题
@number44
看起来确实是相关的,我再详细看下这个文件,谢谢前辈!
@number44
你好,您可能误解了我的意思,目前我已经这样改过并得到了二进制格式的结果文件,想跟您请教的是如何查看它的输出格式,以便从其他软件读取OpenFOAM的结果。
各位前辈好,请问在哪里可以找到OpenFOAM所输出的二进制结果文件的格式呢,想采用二进制来读取OpenFOAM的计算结果
请教各位老师和同学们:
一、由于计算的流体域形状比较复杂,所以尝试用Ansys划分流体域网格,在Ansys中对流体域进行分块,对复杂形状部分进行网格加密,并且在复杂部分采用四面体网格进行贴合。最终网格数量在八百万左右。通过转化后形成openfoam的网格,在checkMesh中出现网格报错信息。下面为checkMesh的输出信息。请教老师和同学们这种情况如何去调整网格,或者这种网格会影响计算的收敛性吗?
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create polyMesh for time = 0
Time = 0
Mesh stats
points: 4200703
internal points: 3796151
faces: 21129446
internal faces: 20649965
cells: 8878252
faces per cell: 4.70581
boundary patches: 10
point zones: 0
face zones: 0
cell zones: 0
Overall number of cells of each type:
hexahedra: 2937143
prisms: 28266
wedges: 0
pyramids: 363851
tet wedges: 0
tetrahedra: 5548992
polyhedra: 0
Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology
wall-__-2 102854 103255 ok (non-closed singly connected)
wall-__-6 60678 61164 ok (non-closed singly connected)
wall-__-3 12304 12645 ok (non-closed singly connected)
wall-__-10 46213 46721 ok (non-closed singly connected)
wall-__-7 23224 23830 ok (non-closed singly connected)
wall-__-8 24985 25355 ok (non-closed singly connected)
wall-__-4 10978 11907 ok (non-closed singly connected)
wall-__-5 11537 12317 ok (non-closed singly connected)
wall-__-__ 38084 38454 ok (non-closed singly connected)
wall-__-9 148624 74571 ok (non-closed singly connected)
Checking geometry...
Overall domain bounding box (-150 -5.7e-14 -5.82643e-12) (277.36 65.8 5.8)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (1.18923e-16 -6.74968e-16 -5.50711e-16) OK.
Max cell openness = 1.04585e-15 OK.
Max aspect ratio = 262.763 OK.
Minimum face area = 2.10595e-05. Maximum face area = 3.10956. Face area magnitudes OK.
Min volume = 1.0868e-07. Max volume = 1.19771. Total volume = 93583.6. Cell volumes OK.
Mesh non-orthogonality Max: 88.7853 average: 15.3095
*Number of severely non-orthogonal (> 70 degrees) faces: 25145.
Non-orthogonality check OK.
<<Writing 25145 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
***Max skewness = 5.26267, 9 highly skew faces detected which may impair the quality of the results
<<Writing 9 skew faces to set skewFaces
Coupled point location match (average 0) OK.
Failed 1 mesh checks.
End
二、在Ansys对流体域划分网格的过程中,由于需要局部加密,将流体域分为若干块,但整体网格共点。将该网格导入openfoam并转化后,发现识别出来的每个面为当时分块时每个块的外表面,这种情况无法进行边界条件的设置。请教各位老师和同学,如何在Ansys中,对整体流体域的各个外表面进行标注,可以在导入openfoam后找到标注的面进行边界条件设置。
下面为后续的报错信息,辛苦各位同学和老师,感谢!!!
[ubuntu:330341] *** Process received signal ***
at ??:?
[ubuntu:330453] *** Process received signal ***
[ubuntu:330341] Signal: Floating point exception (8)
[ubuntu:330341] Signal code: (-6)
[ubuntu:330341] Failing at address: 0x3e800050a65
at ??:?
[94] #7 at ??:?
[78] #8 __libc_start_main[ubuntu:330341] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7fcc84b7e4c0]
[ubuntu:330341] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7fcc84b7e438]
[ubuntu:330341] [ 2] [ubuntu:330453] Signal: Floating point exception (8)
[ubuntu:330453] Signal code: (-6)
[ubuntu:330453] Failing at address: 0x3e800050ad5
[ubuntu:330453] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f6677a5d4c0]
[ubuntu:330453] [ 1] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7fcc84b7e4c0]
[ubuntu:330341] [ 3] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f6677a5d438]
[ubuntu:330453] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f6677a5d4c0]
[ubuntu:330453] [ 3] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7fcc86035f1d]
[ubuntu:330341] [ 4] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7f6678f14f1d]
[ubuntu:330453] [ 4] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7fcc85e4234b]
[ubuntu:330341] [ 5] at ??:?
[ubuntu:330402] *** Process received signal ***
/opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7f6678d2134b]
[ubuntu:330453] [ 5] [ubuntu:330402] Signal: Floating point exception (8)
[ubuntu:330402] Signal code: (-6)
[ubuntu:330402] Failing at address: 0x3e800050aa2
[ubuntu:330402] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f923bdb74c0]
[ubuntu:330402] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f923bdb7438]
[ubuntu:330402] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f923bdb74c0]
[ubuntu:330402] [ 3] ?/opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7fcc8855fbab]
[ubuntu:330341] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330341] [ 7] olaFlow[0x44db4e]
[ubuntu:330341] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7fcc84b69840]
[ubuntu:330341] [ 9] olaFlow[0x456269]
[ubuntu:330341] *** End of error message ***
/opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7f667b43ebab]
[ubuntu:330453] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330453] [ 7] olaFlow[0x44db4e]
[ubuntu:330453] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f6677a48840]
[ubuntu:330453] [ 9] olaFlow[0x456269]
[ubuntu:330453] *** End of error message ***
/opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7f923d26ef1d]
[ubuntu:330402] [ 4] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7f923d07b34b]
[ubuntu:330402] [ 5] ? in "/lib/x86_64-linux-gnu/libc.so.6"
[91] #9 at ??:?
[75] #7 at ??:?
[95] #8 __libc_start_main/opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7f923f798bab]
[ubuntu:330402] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330402] [ 7] olaFlow[0x44db4e]
[ubuntu:330402] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f923bda2840]
[ubuntu:330402] [ 9] olaFlow[0x456269]
[ubuntu:330402] *** End of error message ***
at ??:?
[ubuntu:330384] *** Process received signal ***
[ubuntu:330384] Signal: Floating point exception (8)
[ubuntu:330384] Signal code: (-6)
[ubuntu:330384] Failing at address: 0x3e800050a90
[ubuntu:330384] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f861028e4c0]
[ubuntu:330384] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f861028e438]
[ubuntu:330384] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f861028e4c0]
[ubuntu:330384] [ 3] at ??:?
[ubuntu:330417] *** Process received signal ***
[ubuntu:330417] Signal: Floating point exception (8)
[ubuntu:330417] Signal code: (-6)
[ubuntu:330417] Failing at address: 0x3e800050ab1
[ubuntu:330417] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f72e021f4c0]
[ubuntu:330417] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f72e021f438]
[ubuntu:330417] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f72e021f4c0]
[ubuntu:330417] [ 3] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7f8611745f1d]
[ubuntu:330384] [ 4] at ??:?
[79] #8 __libc_start_main/opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7f861155234b]
[ubuntu:330384] [ 5] ?/opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7f72e16d6f1d]
[ubuntu:330417] [ 4] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7f72e14e334b]
[ubuntu:330417] [ 5] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7f8613c6fbab]
[ubuntu:330384] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330384] [ 7] olaFlow[0x44db4e]
[ubuntu:330384] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f8610279840]
[ubuntu:330384] [ 9] olaFlow[0x456269]
[ubuntu:330384] *** End of error message ***
in "/lib/x86_64-linux-gnu/libc.so.6"
[78] #9 /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7f72e3c00bab]
[ubuntu:330417] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330417] [ 7] olaFlow[0x44db4e]
[ubuntu:330417] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f72e020a840]
[ubuntu:330417] [ 9] olaFlow[0x456269]
[ubuntu:330417] *** End of error message ***
at ??:?
[ubuntu:330407] *** Process received signal ***
[ubuntu:330407] Signal: Floating point exception (8)
[ubuntu:330407] Signal code: (-6)
[ubuntu:330407] Failing at address: 0x3e800050aa7
at ??:?
[ubuntu:330316] *** Process received signal ***
[ubuntu:330316] Signal: Floating point exception (8)
[ubuntu:330316] Signal code: (-6)
[ubuntu:330316] Failing at address: 0x3e800050a4c
[ubuntu:330407] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f7ed4e984c0]
[ubuntu:330407] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f7ed4e98438]
[ubuntu:330407] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f7ed4e984c0]
[ubuntu:330407] [ 3] [ubuntu:330316] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f43261c64c0]
[ubuntu:330316] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f43261c6438]
[ubuntu:330316] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f43261c64c0]
[ubuntu:330316] [ 3] in "/lib/x86_64-linux-gnu/libc.so.6"
[95] #9 /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7f7ed634ff1d]
[ubuntu:330407] [ 4] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7f432767df1d]
[ubuntu:330316] [ 4] ?/opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7f7ed615c34b]
[ubuntu:330407] [ 5] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7f432748a34b]
[ubuntu:330316] [ 5] in "/lib/x86_64-linux-gnu/libc.so.6"
[79] #9 /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7f7ed8879bab]
[ubuntu:330407] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330407] [ 7] olaFlow[0x44db4e]
[ubuntu:330407] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f7ed4e83840]
[ubuntu:330407] [ 9] olaFlow[0x456269]
[ubuntu:330407] *** End of error message ***
/opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7f4329ba7bab]
[ubuntu:330316] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330316] [ 7] olaFlow[0x44db4e]
[ubuntu:330316] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f43261b1840]
[ubuntu:330316] [ 9] olaFlow[0x456269]
[ubuntu:330316] *** End of error message ***
?? at ??:?
[94] #8 __libc_start_main?? at ??:?
[ubuntu:330379] *** Process received signal ***
[ubuntu:330379] Signal: Floating point exception (8)
[ubuntu:330379] Signal code: (-6)
[ubuntu:330379] Failing at address: 0x3e800050a8b
[ubuntu:330379] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7fd834b6e4c0]
[ubuntu:330379] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7fd834b6e438]
[ubuntu:330379] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7fd834b6e4c0]
[ubuntu:330379] [ 3] in "/lib/x86_64-linux-gnu/libc.so.6"
[94] #9 /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7fd836025f1d]
[ubuntu:330379] [ 4] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7fd835e3234b]
[ubuntu:330379] [ 5] at ??:?
[ubuntu:330323] *** Process received signal ***
[ubuntu:330323] Signal: Floating point exception (8)
[ubuntu:330323] Signal code: (-6)
[ubuntu:330323] Failing at address: 0x3e800050a53
[ubuntu:330323] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7fdf7c7694c0]
[ubuntu:330323] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7fdf7c769438]
[ubuntu:330323] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7fdf7c7694c0]
[ubuntu:330323] [ 3] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7fdf7dc20f1d]
[ubuntu:330323] [ 4] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7fd83854fbab]
[ubuntu:330379] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330379] [ 7] olaFlow[0x44db4e]
[ubuntu:330379] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7fd834b59840]
[ubuntu:330379] [ 9] olaFlow[0x456269]
[ubuntu:330379] *** End of error message ***
/opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7fdf7da2d34b]
[ubuntu:330323] [ 5] at ??:?
[75] #8 __libc_start_main/opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7fdf8014abab]
[ubuntu:330323] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330323] [ 7] olaFlow[0x44db4e]
[ubuntu:330323] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7fdf7c754840]
[ubuntu:330323] [ 9] olaFlow[0x456269]
[ubuntu:330323] *** End of error message ***
? at ??:?
[ubuntu:330389] *** Process received signal ***
[ubuntu:330389] Signal: Floating point exception (8)
[ubuntu:330389] Signal code: (-6)
[ubuntu:330389] Failing at address: 0x3e800050a95
at ??:?
[ubuntu:330311] *** Process received signal ***
[ubuntu:330311] Signal: Floating point exception (8)
[ubuntu:330311] Signal code: (-6)
[ubuntu:330311] Failing at address: 0x3e800050a47
[ubuntu:330311] [ 0] [ubuntu:330389] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7fe3de54b4c0]
[ubuntu:330389] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7fe3de54b438]
[ubuntu:330389] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7fe3de54b4c0]
[ubuntu:330389] [ 3] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7fa2535b84c0]
[ubuntu:330311] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7fa2535b8438]
[ubuntu:330311] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7fa2535b84c0]
[ubuntu:330311] [ 3] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7fe3dfa02f1d]
[ubuntu:330389] [ 4] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7fa254a6ff1d]
[ubuntu:330311] [ 4] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7fa25487c34b]
[ubuntu:330311] [ 5] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7fe3df80f34b]
[ubuntu:330389] [ 5] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7fa256f99bab]
[ubuntu:330311] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330311] [ 7] olaFlow[0x44db4e]
[ubuntu:330311] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7fa2535a3840]
[ubuntu:330311] [ 9] olaFlow[0x456269]
[ubuntu:330311] *** End of error message ***
/opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7fe3e1f2cbab]
[ubuntu:330389] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330389] [ 7] olaFlow[0x44db4e]
[ubuntu:330389] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7fe3de536840]
[ubuntu:330389] [ 9] olaFlow[0x456269]
[ubuntu:330389] *** End of error message ***
in "/lib/x86_64-linux-gnu/libc.so.6"
[75] #9 [93] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[ubuntu:330400] *** Process received signal ***
[ubuntu:330400] Signal: Floating point exception (8)
[ubuntu:330400] Signal code: (-6)
[ubuntu:330400] Failing at address: 0x3e800050aa0
[ubuntu:330400] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7ff1e15af4c0]
[ubuntu:330400] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7ff1e15af438]
[ubuntu:330400] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7ff1e15af4c0]
[ubuntu:330400] [ 3] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7ff1e2a66f1d]
[ubuntu:330400] [ 4] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7ff1e287334b]
[ubuntu:330400] [ 5] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7ff1e4f90bab]
[ubuntu:330400] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330400] [ 7] olaFlow[0x44db4e]
[ubuntu:330400] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7ff1e159a840]
[ubuntu:330400] [ 9] olaFlow[0x456269]
[ubuntu:330400] *** End of error message ***
? at ??:?
[ubuntu:330305] *** Process received signal ***
[ubuntu:330305] Signal: Floating point exception (8)
[ubuntu:330305] Signal code: (-6)
[ubuntu:330305] Failing at address: 0x3e800050a41
[ubuntu:330305] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7fac1b5164c0]
[ubuntu:330305] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7fac1b516438]
[ubuntu:330305] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7fac1b5164c0]
[ubuntu:330305] [ 3] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7fac1c9cdf1d]
[ubuntu:330305] [ 4] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7fac1c7da34b]
[ubuntu:330305] [ 5] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7fac1eef7bab]
[ubuntu:330305] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330305] [ 7] olaFlow[0x44db4e]
[ubuntu:330305] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7fac1b501840]
[ubuntu:330305] [ 9] olaFlow[0x456269]
[ubuntu:330305] *** End of error message ***
[109] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[93] #1 Foam::sigFpe::sigHandler(int) at ??:?
[93] #2 ? at ??:?
[109] #1 Foam::sigFpe::sigHandler(int) in "/lib/x86_64-linux-gnu/libc.so.6"
[93] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[109] #2 ? at ??:?
[93] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/lib/x86_64-linux-gnu/libc.so.6"
[109] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[93] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[109] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[109] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[93] #6 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[93] #7 at ??:?
[109] #6 ?Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[93] #8 __libc_start_main at ??:?
[109] #7 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[93] #9 --------------------------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
at ??:?
[109] #8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[109] #9 ? at ??:?
[ubuntu:330431] *** Process received signal ***
? at ??:?
[ubuntu:330390] *** Process received signal ***
[ubuntu:330431] Signal: Floating point exception (8)
[ubuntu:330431] Signal code: (-6)
[ubuntu:330431] Failing at address: 0x3e800050abf
[ubuntu:330431] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f9bda55a4c0]
[ubuntu:330431] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f9bda55a438]
[ubuntu:330431] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f9bda55a4c0]
[ubuntu:330431] [ 3] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7f9bdba11f1d]
[ubuntu:330431] [ 4] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7f9bdb81e34b]
[ubuntu:330431] [ 5] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7f9bddf3bbab]
[ubuntu:330431] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330431] [ 7] olaFlow[0x44db4e]
[ubuntu:330431] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f9bda545840]
[ubuntu:330431] [ 9] olaFlow[0x456269]
[ubuntu:330431] *** End of error message ***
[ubuntu:330390] Signal: Floating point exception (8)
[ubuntu:330390] Signal code: (-6)
[ubuntu:330390] Failing at address: 0x3e800050a96
[ubuntu:330390] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f944a96c4c0]
[ubuntu:330390] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f944a96c438]
[ubuntu:330390] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354c0)[0x7f944a96c4c0]
[ubuntu:330390] [ 3] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEEdRKNS_5UListIT_EES5_+0x2d)[0x7f944be23f1d]
[ubuntu:330390] [ 4] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0xb0b)[0x7f944bc3034b]
[ubuntu:330390] [ 5] /opt/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7f944e34dbab]
[ubuntu:330390] [ 6] olaFlow(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x4db5f1]
[ubuntu:330390] [ 7] olaFlow[0x44db4e]
[ubuntu:330390] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f944a957840]
[ubuntu:330390] [ 9] olaFlow[0x456269]
[ubuntu:330390] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 95 with PID 0 on node ubuntu exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
请教各位老师和同学们:
在一个三维复杂流体域案例的计算中,部分流体域形状较为复杂,网格采用snappyhexmesh生成,网格数量六百五十万左右,checkMesh结果OK,边界设置检查无误。计算到30s左右时出现不收敛的情况。下面为计算过程输出信息中的最后的几个不收敛时间步以及最后的部分报错信息(因为字符限制,后续部分在评论中)。请教各位老师和同学们,报错信息表示了计算模型哪方面的问题以及如何去调整。
Courant Number mean: 0.0811704 max: 339375
Interface Courant Number mean: 0.0331014 max: 90043
deltaT = 3.36389e-15
Time = 33.5117
PIMPLE: iteration 1
alpha.water BC on patch inlet
MULES: Solving for alpha.water
alpha.water BC on patch inlet
Phase-1 volume fraction = 0.647158 Min(alpha.water) = -6.71024e-20 Max(alpha.water) = 1.00032
alpha.water BC on patch inlet
MULES: Solving for alpha.water
alpha.water BC on patch inlet
Phase-1 volume fraction = 0.647158 Min(alpha.water) = -3.42344e-19 Max(alpha.water) = 1.00032
Velocity BC on patch inlet
3D_3D Absorption BC on patch outlet
"Correction Levels" 25( 1.75485 1.85601 1.55984 1.59708 1.56351 1.52502 1.52423 1.52653 1.39188 1.24439 1.50892 1.18572 0.917495 0.860397 1.13003 1.02722 1.37729 1.65376 1.86368 2.10405 2.20481 1.7503 1.49513 1.47808 1.62577 )
3D_3D Absorption BC on patch frontAndBack
"Correction Levels" 25( -4.53749 -4.14273 -0.130794 1.00985 2.2047 3.62849 4.42181 3.58968 1.44737 -1.07649 -3.04945 -3.68013 -1.06624 0.682384 1.53162 2.34454 2.25524 0.924583 0.698287 1.05531 0.708991 0.277332 -0.286046 -0.200047 1.10636 )
DICPCG: Solving for p_rgh, Initial residual = 0.999205, Final residual = 0.0470571, No Iterations 149
DICPCG: Solving for p_rgh, Initial residual = 0.000705545, Final residual = 2.50405e-05, No Iterations 3
time step continuity errors : sum local = 2.35881e+12, global = 0.00174181, cumulative = 0.00137343
DICPCG: Solving for p_rgh, Initial residual = 0.000172039, Final residual = 7.00718e-06, No Iterations 8
DICPCG: Solving for p_rgh, Initial residual = 0.00146033, Final residual = 1.3504e-05, No Iterations 1
Velocity BC on patch inlet
3D_3D Absorption BC on patch outlet
"Correction Levels" 25( 1.75485 1.85601 1.55984 1.59708 1.56351 1.52502 1.52423 1.52653 1.39188 1.24439 1.50892 1.18572 0.917495 0.860397 1.13003 1.02722 1.37729 1.65376 1.86368 2.10405 2.20481 1.7503 1.49513 1.47808 1.62577 )
3D_3D Absorption BC on patch frontAndBack
"Correction Levels" 25( -4.53749 -4.14273 -0.130794 1.00985 2.2047 3.62849 4.42181 3.58968 1.44737 -1.07649 -3.04945 -3.68013 -1.06624 0.682384 1.53162 2.34454 2.25524 0.924583 0.698287 1.05531 0.708991 0.277332 -0.286046 -0.200047 1.10636 )
time step continuity errors : sum local = 1.56951e+12, global = 0.000252718, cumulative = 0.00162615
DICPCG: Solving for p_rgh, Initial residual = 0.000155355, Final residual = 6.95743e-06, No Iterations 8
DICPCG: Solving for p_rgh, Initial residual = 0.00108953, Final residual = 9.66053e-08, No Iterations 238
Velocity BC on patch inlet
3D_3D Absorption BC on patch outlet
"Correction Levels" 25( 1.75485 1.85601 1.55984 1.59708 1.56351 1.52502 1.52423 1.52653 1.39188 1.24439 1.50892 1.18572 0.917495 0.860397 1.13003 1.02722 1.37729 1.65376 1.86368 2.10405 2.20481 1.7503 1.49513 1.47808 1.62577 )
3D_3D Absorption BC on patch frontAndBack
"Correction Levels" 25( -4.53749 -4.14273 -0.130794 1.00985 2.2047 3.62849 4.42181 3.58968 1.44737 -1.07649 -3.04945 -3.68013 -1.06624 0.682384 1.53162 2.34454 2.25524 0.924583 0.698287 1.05531 0.708991 0.277332 -0.286046 -0.200047 1.10636 )
time step continuity errors : sum local = 1.34352e+10, global = -6.4846e+08, cumulative = -6.4846e+08
smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 9.98979e-09, No Iterations 675
bounding epsilon, min: -3.65743e+29 max: 3.75536e+54 average: 2.82799e+48
smoothSolver: Solving for k, Initial residual = 6.76611e-05, Final residual = 1.66848e-08, No Iterations 1000
bounding k, min: 3.91848e-41 max: 6.52827e+46 average: 3.19235e+40
ExecutionTime = 43544.6 s ClockTime = 43560 s
Courant Number mean: 1.67126e+13 max: 3.55514e+19
Interface Courant Number mean: 3.87388e+12 max: 3.55514e+19
deltaT = 4.25793e-35
Time = 33.5117
PIMPLE: iteration 1
alpha.water BC on patch inlet
MULES: Solving for alpha.water
alpha.water BC on patch inlet
Phase-1 volume fraction = 0.647158 Min(alpha.water) = -1.56077e-21 Max(alpha.water) = 1.00032
alpha.water BC on patch inlet
MULES: Solving for alpha.water
alpha.water BC on patch inlet
Phase-1 volume fraction = 0.647158 Min(alpha.water) = -1.56077e-21 Max(alpha.water) = 1.00032
Velocity BC on patch inlet
3D_3D Absorption BC on patch outlet
"Correction Levels" 25( 1.75485 1.85601 1.55984 1.59708 1.56351 1.52502 1.52423 1.52653 1.39188 1.24439 1.50892 1.18572 0.917495 0.860397 1.13003 1.02722 1.37729 1.65376 1.86368 2.10405 2.20481 1.7503 1.49513 1.47808 1.62577 )
3D_3D Absorption BC on patch frontAndBack
"Correction Levels" 25( -4.53749 -4.14273 -0.130794 1.00985 2.2047 3.62849 4.42181 3.58968 1.44737 -1.07649 -3.04945 -3.68013 -1.06624 0.682384 1.53162 2.34454 2.25524 0.924583 0.698287 1.05531 0.708991 0.277332 -0.286046 -0.200047 1.10636 )
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 47.0446, No Iterations 1001
DICPCG: Solving for p_rgh, Initial residual = 0.000216109, Final residual = 9.48986e-06, No Iterations 3
time step continuity errors : sum local = 4.08786e+21, global = 873386, cumulative = -6.47586e+08
DICPCG: Solving for p_rgh, Initial residual = 4.95933e-05, Final residual = 2.30123e-06, No Iterations 12
DICPCG: Solving for p_rgh, Initial residual = 0.000358529, Final residual = 1.65707e-05, No Iterations 5
Velocity BC on patch inlet
3D_3D Absorption BC on patch outlet
"Correction Levels" 25( 1.75485 1.85601 1.55984 1.59708 1.56351 1.52502 1.52423 1.52653 1.39188 1.24439 1.50892 1.18572 0.917495 0.860397 1.13003 1.02722 1.37729 1.65376 1.86368 2.10405 2.20481 1.7503 1.49513 1.47808 1.62577 )
3D_3D Absorption BC on patch frontAndBack
"Correction Levels" 25( -4.53749 -4.14273 -0.130794 1.00985 2.2047 3.62849 4.42181 3.58968 1.44737 -1.07649 -3.04945 -3.68013 -1.06624 0.682384 1.53162 2.34454 2.25524 0.924583 0.698287 1.05531 0.708991 0.277332 -0.286046 -0.200047 1.10636 )
time step continuity errors : sum local = 2.78772e+21, global = 1.35102e+07, cumulative = -6.34076e+08
DICPCG: Solving for p_rgh, Initial residual = 8.00534e-05, Final residual = 3.79662e-06, No Iterations 26
DICPCG: Solving for p_rgh, Initial residual = 0.000430357, Final residual = 9.90187e-08, No Iterations 152
Velocity BC on patch inlet
3D_3D Absorption BC on patch outlet
"Correction Levels" 25( 1.75485 1.85601 1.55984 1.59708 1.56351 1.52502 1.52423 1.52653 1.39188 1.24439 1.50892 1.18572 0.917495 0.860397 1.13003 1.02722 1.37729 1.65376 1.86368 2.10405 2.20481 1.7503 1.49513 1.47808 1.62577 )
3D_3D Absorption BC on patch frontAndBack
"Correction Levels" 25( -4.53749 -4.14273 -0.130794 1.00985 2.2047 3.62849 4.42181 3.58968 1.44737 -1.07649 -3.04945 -3.68013 -1.06624 0.682384 1.53162 2.34454 2.25524 0.924583 0.698287 1.05531 0.708991 0.277332 -0.286046 -0.200047 1.10636 )
time step continuity errors : sum local = 7.72207e+18, global = 1.04902e+18, cumulative = 1.04902e+18
smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 9.99364e-09, No Iterations 967
bounding epsilon, min: -1.30423e+57 max: 3.78086e+105 average: 1.05636e+99
smoothSolver: Solving for k, Initial residual = 0.999975, Final residual = 3.70498e-05, No Iterations 1000
bounding k, min: 1.10387e-72 max: 5.9365e+94 average: 3.93865e+88
ExecutionTime = 43557.5 s ClockTime = 43573 s
Courant Number mean: 6.08378e+22 max: 5.54233e+28
Interface Courant Number mean: 8.58371e+21 max: 4.62459e+28
deltaT = 3.45715e-64
Time = 33.5117
PIMPLE: iteration 1
alpha.water BC on patch inlet
MULES: Solving for alpha.water
alpha.water BC on patch inlet
Phase-1 volume fraction = 0.647158 Min(alpha.water) = -1.56077e-21 Max(alpha.water) = 1.00032
alpha.water BC on patch inlet
MULES: Solving for alpha.water
alpha.water BC on patch inlet
Phase-1 volume fraction = 0.647158 Min(alpha.water) = -1.56077e-21 Max(alpha.water) = 1.00032
Velocity BC on patch inlet
3D_3D Absorption BC on patch outlet
"Correction Levels" 25( 1.75485 1.85601 1.55984 1.59708 1.56351 1.52502 1.52423 1.52653 1.39188 1.24439 1.50892 1.18572 0.917495 0.860397 1.13003 1.02722 1.37729 1.65376 1.86368 2.10405 2.20481 1.7503 1.49513 1.47808 1.62577 )
3D_3D Absorption BC on patch frontAndBack
"Correction Levels" 25( -4.53749 -4.14273 -0.130794 1.00985 2.2047 3.62849 4.42181 3.58968 1.44737 -1.07649 -3.04945 -3.68013 -1.06624 0.682384 1.53162 2.34454 2.25524 0.924583 0.698287 1.05531 0.708991 0.277332 -0.286046 -0.200047 1.10636 )
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0479307, No Iterations 15
[75] #0 Foam::error::printStack(Foam::Ostream&)[78] #0 Foam::error::printStack(Foam::Ostream&)[79] #0 Foam::error::printStack(Foam::Ostream&)[82] #0 Foam::error::printStack(Foam::Ostream&)[83] #0 Foam::error::printStack(Foam::Ostream&)[90] #0 Foam::error::printStack(Foam::Ostream&)[91] #0 Foam::error::printStack(Foam::Ostream&)[94] #0 Foam::error::printStack(Foam::Ostream&)[95] #0 Foam::error::printStack(Foam::Ostream&)[96] #0 Foam::error::printStack(Foam::Ostream&)[97] #0 Foam::error::printStack(Foam::Ostream&)[100] #0 Foam::error::printStack(Foam::Ostream&)[119] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[95] #1 Foam::sigFpe::sigHandler(int) at ??:?
at ??:?
at ??:?
at ??:?
[83] #1 Foam::sigFpe::sigHandler(int)[90] #1 Foam::sigFpe::sigHandler(int)[91] #1 Foam::sigFpe::sigHandler(int) at ??:?
[97] #1 Foam::sigFpe::sigHandler(int) at ??:?
[78] #1 Foam::sigFpe::sigHandler(int) at ??:?
at ??:?
at ??:?
[82] #1 Foam::sigFpe::sigHandler(int)[100] #1 Foam::sigFpe::sigHandler(int)[96] #1 Foam::sigFpe::sigHandler(int)[119] #1 Foam::sigFpe::sigHandler(int) at ??:?
at ??:?
[94] #1 Foam::sigFpe::sigHandler(int) at ??:?
[75] #1 Foam::sigFpe::sigHandler(int)[79] #1 Foam::sigFpe::sigHandler(int) at ??:?
[83] #2 ? at ??:?
[91] #2 ? at ??:?
[96] #2 ? at ??:?
[119] #2 ? at ??:?
[90] #2 ? at ??:?
[97] #2 ? at ??:?
[100] #2 ? at ??:?
[82] #2 ? at ??:?
[78] #2 ? at ??:?
[95] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
in "/lib/x86_64-linux-gnu/libc.so.6"
[91] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&)[83] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
in "/lib/x86_64-linux-gnu/libc.so.6"
[97] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&)[119] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[79] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
in "/lib/x86_64-linux-gnu/libc.so.6"
[90] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
in "/lib/x86_64-linux-gnu/libc.so.6"
[100] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&)[82] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&)[96] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[75] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[95] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[94] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[78] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
in "/lib/x86_64-linux-gnu/libc.so.6"
[91] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const[79] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
at ??:?
at ??:?
in "/lib/x86_64-linux-gnu/libc.so.6"
[97] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const[119] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[83] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const[75] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&)[90] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
in "/lib/x86_64-linux-gnu/libc.so.6"
[82] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const[94] #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[100] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[96] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[95] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[78] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
at ??:?
at ??:?
[119] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&)[97] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&)[79] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
at ??:?
[90] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[82] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&)[96] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
at ??:?
at ??:?
[94] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[83] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&)[100] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&)[91] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[78] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[75] #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
at ??:?
[95] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&)[79] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[94] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[97] #6 at ??:?
[119] #6 at ??:?
[82] #6 at ??:?
[96] #6 at ??:?
[90] #6 at ??:?
[83] #6 at ??:?
[100] #6 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
Foam::fvMatrix<double>::solve(Foam::dictionary const&)[75] #5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&)Foam::fvMatrix<double>::solve(Foam::dictionary const&)Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[91] #6 Foam::fvMatrix<double>::solve(Foam::dictionary const&)Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[78] #6 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[97] #7 at ??:?
[119] #7 at ??:?
[82] #7 at ??:?
[95] #6 at ??:?
[79] #6 at ??:?
[90] #7 at ??:?
[100] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&)Foam::fvMatrix<double>::solve(Foam::dictionary const&)??? at ??:?
[83] #7 at ??:?
[96] #7 ?? at ??:?
[91] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[94] #6 Foam::fvMatrix<double>::solve(Foam::dictionary const&)?? at ??:?
[119] #8 __libc_start_main at ??:?
[82] #8 __libc_start_main at ??:?
[90] #8 __libc_start_main at ??:?
[97] #8 __libc_start_main at ??:?
[100] #8 __libc_start_main at ??:?
[78] #7 at ??:?
[75] #6 ? at ??:?
[95] #7 at ??:?
[79] #7 in "/lib/x86_64-linux-gnu/libc.so.6"
[82] #9 in "/lib/x86_64-linux-gnu/libc.so.6"
[119] #9 at ??:?
[83] #8 __libc_start_main at ??:?
[96] #8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[97] #9 in "/lib/x86_64-linux-gnu/libc.so.6"
[90] #9 in "/lib/x86_64-linux-gnu/libc.so.6"
[100] #9 ?Foam::fvMatrix<double>::solve(Foam::dictionary const&)Foam::fvMatrix<double>::solve(Foam::dictionary const&)?? at ??:?
[91] #8 __libc_start_main?? in "/lib/x86_64-linux-gnu/libc.so.6"
[83] #9 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[96] #9 ?? at ??:?
@李东岳 了解了,谢谢东岳老师
@myjimp
你好,请问snappyHexMesh并行剖切网格后如何将网格合起来呢?
在找资料的过程中发现了关于sampleDict中可用type的说明,包括uniform, face, midPoint, midPointAndFace, curve, cloud(之前在苏军伟前辈的博客中看到过相关介绍,但当时并无curve类型)。
可以看到是有针对曲线的类型(curve),但是没有给出具体的做法,后面会尝试用这个提取一下曲面上的压力看看。
各位前辈好,想请教一下有没有类似midPoint这种提一条线上各单元物理量的方法,来提取像下图这种曲面上的各个单元的物理量,例如压力p。
(虽然用cloud方法也能提出来,但需要一个个找曲面上的单元坐标,单元多了的话感觉会相当麻烦 )
@cfd_上弦月 已经解决了,点击左侧的选项''Zone Style'',确定内部流场对应的zone(以上面的图片为例,我这里是Step 1737 Incr 0)。选择''Surfaces''标签,然后在''Surfaces to Plot''处右键选择All surfaces,内部的流体域就能显示出来了。
@huoyanxiadebing 在 说实话,根本不知道选择什么湍流模型 中说:
看到你的帖子,发现你做的都是外流场,对于外部流动,沿表面位置分布的雷诺数大于500000时,通常认为流动状态为湍流;沿障碍物的雷诺数大于20000时认为流动状态为湍流,如果你的模型中雷诺数不是很大的话,你用层流模型计算即可,这就是有时候你用层流模型算的比较准确的原因;对于雷诺数的计算,最重要的一点就是你要选好特征长度,对于没有结构物时,你要以水槽的长度或高度来算雷诺数,加上了结构物,就要判断是否产生了湍流现象,以结构物的某个长度及其附近的速度来计算雷诺数。
前辈您好,请问这里提到的两个临界值是否有相应文献的支撑?
如果要计算防波堤(例如沉箱防波堤和水下防波堤)附近的雷诺数来选取相应的湍流模型(适用于低雷诺数或高雷诺数的),这个特征长度应该取防波堤的高度吗?对于二维和三维情况,特征长度的选取是否有不同呢?
另外,结构物附近的速度是否有较为通用的估算方式?对于防波堤附近的速度,能否用基于浅水波理论的公式Uh=cη进行估算呢(水深h,波速c,自由面高度η)?
@李东岳
东岳老师,您这里说的直接给点的位移是指在下一步求解前,对constant/polyMesh路径下的points文件里相应边界上的点坐标进行覆写吗?
然后在0文件夹下给出该边界的各个初始条件,同时在pointDisplacement指定相应的movingWall边界,内部再使用某种动网格求解器进行处理?
@李东岳
感谢东岳老师,起初我以为边界点和内部是动网格方法一起求解,看来之前的理解是有问题的。
请问这样的边界点移动在OpenFOAM中是如何实现的呢?
这是做的波浪、海床与结构物的相互作用,由于结构物下沉后海床变形(有限元部分),所以想用算出的位移来更新流体域网格从而更准确地计算下一步的流场。
preCICE is the leading open-source coupling library for partitioned multi-physics simulations, including, but not restricted to fluid-structure interaction and conjugate heat transfer simulations.
Partitioned means that preCICE couples existing programs/solvers capable of simulating a subpart of the complete physics involved in a simulation. This allows for the high flexibility that is needed to keep a decent time-to-solution for complex multi-physics scenarios.
The software offers convenient methods for transient equation coupling, communication, and data mapping.
visualisation of how preCICE couples different solvers
感谢前辈,我刚刚看了preCICE官网的演示视频,感觉值得一试。在介绍里看到这个团队已经开发了OpenFOAM的adapter,之后要做的应该就是参考这些现有的adapter并针对自己采用的有限元软件做一个adapter来实现和OpenFOAM之间的数据传递了。
@bestucan
感谢前辈的分享,不过现阶段可能还不会用到重叠网格。
目前想实现的是通过读取其他软件计算的位移来更新局部边界处各个单元的网格,最主要的可能还是一个接口的问题。此外,想要更新的边界处各单元的位移是不同的,不知道如附图这样的网格变形是否可以用displacementLaplacian实现?
看到pimpleDyMFoam和interDyMFoam的求解过程中调用了一个mesh.update()函数来更新网格,但是我在C++ Source Code Guide中没有找到这个函数具体是怎样实现的,还请各位前辈指点一下。
@fireztw 在 如何将6DOF和动边界一起使用? 中说:
楼主您好,动边界是施加的主动运动,浮体是在流体作用下的被动运动,这样理解对吗?请问这种流固耦合问题,使用xxDyMFoam求解器可以计算吗?不需要用专门的FsiFoam吗?谢谢!另外,如果只是浮体运动,参考哪个算例比较好啊
各位前辈好,最近在看与OpenFOAM动网格技术相关的资料,想将有限元软件计算得到的结构物位移施加到OpenFOAM的流体网格上来使流体网格变形,并计算下一步的压力场。想请问OpenFOAM自带的动网格求解器能否处理这种问题或有没有其他的方法可以解决呢?
各位前辈好,我最近也在学习fsiFoam,对于其中的fluidSolidInterface部分有些疑问,不知道是否方便提供邮箱等联系方式来请教一下
@wwzhao 明白了 谢谢前辈
@赵一铭
你好,请问对于并行计算的案例该如何用gnuplot实时监测残差呢?
@队长别开枪 好的 感谢!
@队长别开枪 了解了,谢谢!
@队长别开枪 感谢您的建议!
另外想请问您说的每个patch单独保存一个stl文件,这个用CAD软件可以实现吗 ?
(之前一直是用的CAD导出stl文件,应该是您说的整个区域为一个stl文件)