icoFoam运算出错

zzq704

各位老师和同学，本人是打算用icoFoam运算一个从Icem里面转化来的网格，在使用icoFoam试运行的时候，出现了如下错误，第一步就出错了。请问是什么问题呢？Icem里面的网格包含Interface,在openFoam里面设置了cyclicAMI代替，我猜测是不是这一步出了问题，因为interface相邻面面积不一致，我看cyclicAMI介绍好像有说相邻面大小要一样。
openFoam小白，还望各位赐教。

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  9
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 9-710dba47639b
Exec   : icoFoam
Date   : Mar 22 2022
Time   : 20:18:23
Host   : "nm"
PID    : 58789
I/O    : uncollated
Case   : /home/nm/OpenFOAM/nm/cavityceshi2/cavity
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi

AMI: Creating addressing and weights between 151 source faces and 134 target faces
AMI: Patch source sum(weights) min/max/average = 0, 1, 0.841191
AMI: Patch target sum(weights) min/max/average = 0, 1, 0.813281

Starting time loop

Time = 0.0005

Courant Number mean: 0 max: 0
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 9.17761e-06, No Iterations 77
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 9.92611e-06, No Iterations 97
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 9.44364e-06, No Iterations 71
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
#4  ? in "/home/nm/OpenFOAM/OpenFOAM-9/platforms/linux64GccDPInt32Opt/bin/icoFoam"
#5  ? in "/home/nm/OpenFOAM/OpenFOAM-9/platforms/linux64GccDPInt32Opt/bin/icoFoam"
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  ? in "/home/nm/OpenFOAM/OpenFOAM-9/platforms/linux64GccDPInt32Opt/bin/icoFoam"
浮点数例外 (核心已转储)
nm@nm:~/OpenFOAM/nm/cavityceshi2/cavity$

zzq704

@李东岳
大佬救命，搞不懂嘎

李东岳

@zzq704 感觉像是除以0了，你看看你边界条件都对么，以及网格什么样的

zzq704

@李东岳
感谢老师回复
顺便问一下，cyclicAMI两个面，其面积是可以不等的吗，我在纠结是不是这里有问题导致的出错。例如是下图两个面。

2022-03-25 10-52-21屏幕截图.png

2022-03-25 10-52-47屏幕截图.png

李东岳

@zzq704 不可以要相同

zzq704

@李东岳
那么如果我在icem里面画的网格，interface的两个面大小不相等，要导入到openFoam里面，有什么解决办法吗？

李东岳

@zzq704 可以把你第一个图的面，分割成3个面，中间的面不是cyclicAMI，上下的为cyclicAMI

zzq704

@李东岳好的，谢谢李老师