Re: エラーコードが出て計算が止まる件について

[kominami]

maさん

kominamiです。

エラーコートが出るタイミングから、経験上では境界条件の設定が誤っているような気がします。例えば流量収支が合っていないとか。

Time=1 の Solving for e, と Solving for p_rgh, で、 No Iterations 1000 となっているので、内部スイープ回数の上限になっています。

blockMeshDictの変更前でも、1000回になっていたでしょうか？

＞blockMeshを少しいじって今までと同じ方法で

＞計算を実行しようとしたところ、以下のエラーコードが出て

＞Time=2で計算が止まってしいます。

＞blockMesh自体は間違っていないと考えられます。

疑って悪いのですが、本当にメッシュは正しいのでしょうか？

2022年4月15日金曜日 11:09:58 UTC+9 ma:

お世話になっております。

maと申します。

blockMeshを少しいじって今までと同じ方法で

計算を実行しようとしたところ、以下のエラーコードが出て

Time=2で計算が止まってしいます。

blockMesh自体は間違っていないと考えられます。

計算の条件を変えてみたりしたのですが、

一向に計算が進みません。

どのように条件を変えたらうまく行くのかわかる方が

いらっしゃいましたらご教授いただけると幸いです。

ma

Time = 1

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 8.87402e-08, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 2.26712e-07, No Iterations 2
smoothSolver:  Solving for e, Initial residual = 1, Final residual = 0.000305756, No Iterations 1000
GAMG:  Solving for p_rgh, Initial residual = 0.996406, Final residual = 1.19162e+17, No Iterations 1000
time step continuity errors : sum local = 2.72563e+19, global = -3.00769e+18, cumulative = -3.00769e+18
ExecutionTime = 21.96 s  ClockTime = 22 s

Time = 2

smoothSolver:  Solving for Ux, Initial residual = 0.000230559, Final residual = 2.80173e-08, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.000350761, Final residual = 1.96329e-08, No Iterations 1
[2] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[2] #1  Foam::sigFpe::sigHandler(int) at ??:?
[2] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #3  Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
[2] #4  Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
[2] #5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[2] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[2] #7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[2] #8  Foam::fvMatrix<double>::solve() in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[2] #9  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[2] #10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #11  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[kam-Endeavor-MR4800E:05039] *** Process received signal ***
[kam-Endeavor-MR4800E:05039] Signal: Floating point exception (8)
[kam-Endeavor-MR4800E:05039] Signal code:  (-6)
[kam-Endeavor-MR4800E:05039] Failing at address: 0x3e8000013af
[kam-Endeavor-MR4800E:05039] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef10)[0x7fe54c210f10]
[kam-Endeavor-MR4800E:05039] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7fe54c210e87]
[kam-Endeavor-MR4800E:05039] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef10)[0x7fe54c210f10]
[kam-Endeavor-MR4800E:05039] [ 3] /opt/openfoam9/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam19GaussSeidelSmoother6smoothERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5_RKNS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEEhi+0x425)[0x7fe54d5da6d5]
[kam-Endeavor-MR4800E:05039] [ 4] /opt/openfoam9/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam19GaussSeidelSmoother6smoothERNS_5FieldIdEERKS2_hi+0x28)[0x7fe54d5daaa8]
[kam-Endeavor-MR4800E:05039] [ 5] /opt/openfoam9/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam12smoothSolver5solveERNS_5FieldIdEERKS2_h+0x416)[0x7fe54d5d2286]
[kam-Endeavor-MR4800E:05039] [ 6] /opt/openfoam9/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x17d)[0x7fe55123eddd]
[kam-Endeavor-MR4800E:05039] [ 7] buoyantSimpleFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x15b)[0x55eb009caaab]
[kam-Endeavor-MR4800E:05039] [ 8] buoyantSimpleFoam(_ZN4Foam8fvMatrixIdE5solveEv+0x12f)[0x55eb009cae4f]
[kam-Endeavor-MR4800E:05039] [ 9] buoyantSimpleFoam(+0x36158)[0x55eb00962158]
[kam-Endeavor-MR4800E:05039] [10] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7fe54c1f3c87]
[kam-Endeavor-MR4800E:05039] [11] buoyantSimpleFoam(+0x3a67a)[0x55eb0096667a]
[kam-Endeavor-MR4800E:05039] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 2 with PID 0 on node kam-Endeavor-MR4800E exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

[haruka tsubota @xsim.info]

> GAMG:  Solving for p_rgh, Initial residual = 0.996406, Final residual = 1.19162e+17, No Iterations 1000
> time step continuity errors : sum local = 2.72563e+19, global = -3.00769e+18, cumulative = -3.00769e+18

とあるので、1サイクル目で圧力が発散しています

blockMeshをどう変更したか書かれていないので原因は回答しようがないですが、もし「blockMesh自体は間違っていない」と確信があるのであれば、メッシュの粗密が変わったことによって計算が不安定化しているのでしょう。

とりあえず緩和係数（fvSolution内のrelaxationFactors）を小さくして （1/10など） 計算が安定化しないかを確認するといいと思います。

[ma]

kominamiさん

Solving for e, と Solving for p_rgh, を確認したところ、 

blockMeshの変更前もNo Iterations 1000となっていました。

blockMeshと境界条件をもう一度確認し直してみます。

ありがとうございます。

ma

2022年4月17日日曜日 1:51:56 UTC+9 kominami:

[ma]

harukaさん

ありがとうございます。

緩和係数を小さくして計算をしたところ、

Time=3までは計算が進んだのですが

結局同じようなエラーコードが出てきました。

blockMeshをもう一度確認してみます。

ma

2022年4月19日火曜日 8:48:37 UTC+9 haruka tsubota @xsim.info: