sonicFoamでのコアダンプについて

[kwck]

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

現在sonicFoamを使って層流計算をしているのですが、
「浮動小数点例外」が表示され計算が走らず困って折ります。

過去の投稿を参考にさせていただいたとろ、境界条件の設定に問題があることが多いとのことでした。
しかし自身が設定した境界条件を確認しましたが、得に問題はなさそうに思います。

この場合、どのように対処すればよいでしょうか。
下にログを添付します。
どのような情報でも構いませんので、どなたかご教授いただけますでしょうか。

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

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type laminar
Creating field kinetic energy K


PIMPLE: no residual control data found. Calculations will employ 2 corrector loops


Starting time loop

Time = 1

Courant Number mean: 0 max: 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 9.67356e-06, No Iterations 46
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 9.45992e-06, No Iterations 61
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 9.89304e-06, No Iterations 45
smoothSolver:  Solving for e, Initial residual = 0.0252458, Final residual = 0.0215646, No Iterations 1000
smoothSolver:  Solving for p, Initial residual = 1, Final residual = 0.0102958, No Iterations 1000
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = -1.05372, global = 1.05243, cumulative = 1.05243
PIMPLE: iteration 2
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::symGaussSeidelSmoother::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) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/sonicFoam"
#7 
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/sonicFoam"
#8 
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/sonicFoam"
#9 
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/sonicFoam"
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11 
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/sonicFoam"
浮動小数点例外 (コアダンプ)

[Youhei Takagi]

高木と申します。

sonicFoamは使ったことがないので的外れかもしれませんが、
出力から判断すると、エネルギー式の反復計算で最大反復の
1000回に達しても残差がほとんど減っておらず、そのため
計算が発散していると思われます。チュートリアルの設定から
変更した箇所を比較して非物理的な設定になっていないかチェック
することをおすすめします。

2016年1月12日 13:58 kwck <kwc...@gmail.com>:

--
このメールは Google グループのグループ「OpenFOAM」に登録しているユーザーに送られています。
このグループから退会し、グループからのメールの配信を停止するには openfoam+u...@googlegroups.com にメールを送信してください。
このグループに投稿するには open...@googlegroups.com にメールを送信してください。
https://groups.google.com/group/openfoam からこのグループにアクセスしてください。
その他のオプションについては https://groups.google.com/d/optout にアクセスしてください。

[kwck]

高木様

ご返信ありがとうございます。

境界条件の設定に物理的におかしい箇所があったため発散したようです。

再度境界条件を見直し設定すると無事計算が走りました。

ご教授頂き誠にありがとうございました。