閉領域のinterFoam追加計算でエラー

[art.sugawara]

みなさまにはいつもお世話になっております。

主題のエラーなのですが原因・対処わからずアドバイスを頂けたら幸いです。
いままでは解析境界を開放状態で計算していたものを閉領域に変えたところ
エラー発生しました。

モデルはpatch条件をwall条件にしたのと、fvOptionsにpRefPointとpRefValue
を追加しただけです。

 

計算は流れるのですが、途中で止めて、最終時刻から追加計算させると
エラーとなります。途中で止めるのは材料特性を途中から変えたいためですが、
なにも修正しなくても追加計算ではエラーになります。途中で止めなければ
不思議ですが計算は続行します。

エラーメッセージを以下に添付します。（time=0.0001まで計算後の再実行）

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.2                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.2.2-9739c53ec43f
Exec   : interFoam
Date   : Dec 08 2014
Time   : 12:17:50
Host   : "dexcs2012"
PID    : 10446
Case   : /opt/OpenFOAM/me8000166-2.2.2/run/GC_shape3b
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0.0001
PIMPLE: Operating solver in PISO mode
Reading field p_rgh
Reading field U
Reading/calculating face flux field phi
Reading transportProperties
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
Reading g
Calculating field g.h
No finite volume options present
time step continuity errors : sum local = 1.52136e-07, global = -1.52039e-07, cumulative = -1.52039e-07
--> FOAM FATAL ERROR:
Continuity error cannot be removed by adjusting the outflow.
Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow.
Total flux              : 2.22806e-06
Specified mass inflow   : 3.28797e-09
Specified mass outflow  : 5.04532e-11
Adjustable mass outflow : 0
    From function adjustPhi(surfaceScalarField& phi, const volVectorField& U,const volScalarField& p
    in file cfdTools/general/adjustPhi/adjustPhi.C at line 118.
FOAM exiting

メッセージは「outflowの調整で連続性誤差を除去できなかった」と読めます。意味は図りかねますが。
閉領域と関係があるようにも感じますが。。。
よろしくお願いします。

 

[ohbuchi]

境界条件の設定に問題がありそうです。設定ファイルを提示頂けないと詳細はわかりません。

2014年12月9日火曜日 11時15分57秒 UTC+9 art.sugawara:

[art.sugawara]

ohbuchi 様

 

いつもお世話になっております。
解析モデルを添付いたしました。よろしくお願いします。

 

解析領域を閉領域にした目的は、開放境界での解析で
alpha1の積分値が減少する挙動に気が付き、閉領域に
すれば質量保存が保たれるのではないかと考えました。
しかし実際に解析してみると、やはりalpha1の変動は
起こってしまい、不適切な設定が他にあるのだと
想像しています。
これについてもアドバイスいただければと思います。

2014年12月9日火曜日 14時03分29秒 UTC+9 ohbuchi:

[ohbuchi]

p_rghの境界条件が間違っています。

zeroGradient -> fixedFluxPressure

壁面で静水圧を差し引いた圧力の勾配が０となる必要があります。
初期スタートで問題にならないのは、phiがゼロだからです。

エラーを発生しているのは、CorrectPhi.HのadjustPhi関数で

これは、計算初期にしか呼ばれません。

ご参考まで。

2014年12月9日火曜日 15時49分53秒 UTC+9 art.sugawara:

[art.sugawara]

ohbuchi 様

 

ご指摘ありがとうございます。

境界条件を修正したところ、エラーとならずに追加計算実行できました。

2014年12月9日火曜日 21時07分46秒 UTC+9 ohbuchi: