DEXCS2016-OF(Ver4.X)で、kOmegaSSTが計算できません。

小南秀彰

unread,

Dec 26, 2016, 5:18:34 AM12/26/16

to OpenFOAM

みなさま　　小南です。

私は、DEXCS2016-OpenFOAMを使っています（OpenFOAM Ver4.X ）

乱流モデルをkOmegaSSTにするとエラーが出ます。laminar,kEpsilon,RNGkEpsilon,realizableKE,kOmegaの全ては問題ありません。並列計算の有無に関わらず同じエラーが出ます。

ログではwallDict(壁関数？)が見つからないというエラーのように思えますが、kOmegaのほうは問題ないため不思議です。

対処方法を御存じの方は、御教示いただきたく。宜しくお願いします。

================================

kOmegaSST の場合の端末画面のログ

================================

/*---------------------------------------------------------------------------*\

| ========= | |

| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |

| \\ / O peration | Version: 4.x |

| \\ / A nd | Web: www.OpenFOAM.org |

| \\/ M anipulation | |

\*---------------------------------------------------------------------------*/

Build : 4.x-be7fba6cff9b

Exec : simpleFoam -parallel

Date : Dec 26 2016

Time : 18:57:47

Host : "dexcs-VirtualBox"

PID : 20640

Case : /home/dexcs/Desktop/tut/tut01_RANS

nProcs : 2

Slaves : 1("dexcs-VirtualBox.20641")

Pstream initialized with:

floatTransfer : 0

nProcsSimpleSum : 0

commsType : nonBlocking

polling iterations : 0

sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

fileModificationChecking : Monitoring run-time modified files using timeStampMaster

allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Create time

Create mesh for time = 0

SIMPLE: convergence criteria

field p tolerance 0.001

field U tolerance 0.001

field "(k|epsilon)" tolerance 0.001

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian

Selecting turbulence model type RAS

Selecting RAS turbulence model kOmegaSST

[1]

[1] --> FOAM FATAL IO ERROR:

[1] keyword wallDist is undefined in dictionary "/home/dexcs/Desktop/tut/tut01_RANS/processor1/system/fvSchemes"

[1]

[1] file: /home/dexcs/Desktop/tut/tut01_RANS/processor1/system/fvSchemes at line 0.

[1]

[1] From function const Foam::dictionary& Foam::dictionary::subDict(const Foam::word&) const

[1] in file db/dictionary/dictionary.C at line 642.

[1]

FOAM parallel run exiting

[1]

--------------------------------------------------------------------------

MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD

with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.

You may or may not see output from other processes, depending on

exactly when Open MPI kills them.

--------------------------------------------------------------------------

Killing PID 20632

PyFoam WARNING on line 271 of file /usr/local/lib/python2.7/dist-packages/PyFoam-0.6.6-py2.7.egg/PyFoam/Execution/FoamThread.py : Process 20632 was already dead

Getting LinuxMem: [Errno 2] No such file or directory: '/proc/20632/status'

================================

kOmegaSST の場合の端末画面のログ(END)

================================

==============================

kOmega の場合の端末画面のログ

==============================

/*---------------------------------------------------------------------------*\

| ========= | |

| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |

| \\ / O peration | Version: 4.x |

| \\ / A nd | Web: www.OpenFOAM.org |

| \\/ M anipulation | |

\*---------------------------------------------------------------------------*/

Build : 4.x-be7fba6cff9b

Exec : simpleFoam -parallel

Date : Dec 26 2016

Time : 18:55:50

Host : "dexcs-VirtualBox"

PID : 20094

Case : /home/dexcs/Desktop/tut/tut01_RANS

nProcs : 2

Slaves : 1("dexcs-VirtualBox.20095")

Pstream initialized with:

floatTransfer : 0

nProcsSimpleSum : 0

commsType : nonBlocking

polling iterations : 0

sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

fileModificationChecking : Monitoring run-time modified files using timeStampMaster

allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Create time

Create mesh for time = 0

SIMPLE: convergence criteria

field p tolerance 0.001

field U tolerance 0.001

field "(k|epsilon)" tolerance 0.001

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian

Selecting turbulence model type RAS

Selecting RAS turbulence model kOmega

kOmegaCoeffs

{

betaStar 0.09;

beta 0.072;

gamma 0.52;

alphaK 0.5;

alphaOmega 0.5;

}

No MRF models present

No finite volume options present

Starting time loop

===================================

kOmega の場合の端末画面のログ(END) このあと正常に　time loopに入ります。

==================================

nakagawa

unread,

Dec 26, 2016, 5:49:14 AM12/26/16

to OpenFOAM

小南様

version 4 から，壁面からの距離の算出方法を指定するようになったようです。
そのため，wallDistを設定する必要があります。
下記例題などが参考になります。

https://github.com/OpenFOAM/OpenFOAM-4.x/blob/master/tutorials/incompressible/simpleFoam/pitzDaily/system/fvSchemes

なかがわ

2016年12月26日月曜日 19時18分34秒 UTC+9 小南秀彰:

小南秀彰

unread,

Dec 27, 2016, 2:05:46 AM12/27/16

to OpenFOAM

nakagawa様

動作しました。ありがとうございました。

DEXCS2016のランチャーで出来るファイルがv2.0のままなのですね。

DEXCS2016のHDDの中のtutorialを調べると、pitzDailyだけでなく、windAroundBuildings、rotorDick、motorBikeなども同じようwallDistの行が増えていました。

2016年12月26日月曜日 19時49分14秒 UTC+9 nakagawa:

Fumiya Nozaki

unread,

Dec 28, 2016, 12:48:09 AM12/28/16

to open...@googlegroups.com

壁面からの距離を使用する乱流モデル（kOmegaSST、SpalartAllmaras、LienCubicKE など）では、

wallDist の設定が必要です。

対して、使用しない乱流モデル（kEpsilon など）では、設定不要です。

それぞれのモデルのソースコードにおいて、次の行の有無に対応しています。

y_(wallDist::New(this->mesh_).y())

ご参考まで。

野崎

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

Reply all

Reply to author

Forward