浮動小数点例外が出てしまう件

maho sawada

unread,

Dec 3, 2023, 5:24:07 PM12/3/23

to OpenFOAM

初めまして。

OpenFfoam初心者の沢田と申します。

私は現在大学でOpenFoamを用いて混相流の解析を行っております。

内容としましては、空気で充満された室内にアンモニアを流入させ、拡散の様子を見るものです。

Xsimにて設定を行った後にOpenFoamで実行をしようとしたところ、浮動小数点例外が出てしまい、困難に陥っています。自分では解決できなかったため、皆さまのお力をお借りしたく思いご連絡させていただきました。

エラー内容は以下のようになってます。

また添付ファイルが実行したものとなっております。

他に必要な情報ございましたらお知らせ願います。

No MRF models present

No fvModels present
No fvConstraints present
#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::adjustPhi(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) at ??:?
#4 ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam"
#5 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam"
Floating point exception

room10+2.zip

H‐Tsubota＠xsim.info

unread,

Dec 4, 2023, 1:47:53 AM12/4/23

to OpenFOAM

手元で計算して確認してみました。計算が発散しています。
おそらくAssembly1にAssemblyからNH3が流入してAssembly0から流出するという計算だと思いますが、以下の点が気になりました。

1. x、y 、ｚ方向の長さがそれぞれ8km、6km、10kmあるがこの縮尺で間違いないか。
2. Assembly1の空間にメッシュが切られていない。XSimの「メッシュ」の「計算領域」設定の座標が境界面上にあるためだと思います。設定を見直すべきです。

とりあえず上記「2」の「計算領域」設定を(0.0, 0.0, 5000.0)から(0.0 , 0.0, 5100.0)に変更すれば計算は流れるようです。

2023年12月4日月曜日 7:24:07 UTC+9 maho sawada:

fune sawada

unread,

Dec 4, 2023, 8:42:59 PM12/4/23

to OpenFOAM

ご返信いただきましてありがとうございます。

ご指導いただきました通り、

1→形状の編集にてkmをmにしておりませんでしたので、修正いたしました。

2→計算領域の設定を(0 0 5100)にいたしました。計算領域の設定等を理解しておりませんでした。

以上の修正をしたところ、最後まで実行され結果を得ることができました。

そのような設定が必要であることを自分だけでは気づくことができませんでした。

迅速かつ的確なアドバイスをいただきましたことを、心より感謝申し上げます。

今後も質問させていただくことがございましたら、その際は何卒よろしくお願いいたします。

最後になりますが重ねてお礼申し上げます。誠にありがとうございました。