OpenFOAMによる要素数約1億の解析に関しまして
2,600 views
Skip to first unread message
Konishi
Dec 19, 2012, 4:11:40 AM12/19/12
to open...@googlegroups.com
Google Groups 「OpenFOAM」の皆様
お世話に成って居ります。konishiと申します。今回はOpenFOAMによる約1億要素の解析(100並列)に関しまして質問が御座います。
現在、東京大学様のサーバーをお借りして要素数約1億の解析を行おうとしております。解析対象は自動車で、抵抗係数に注目した解析を目指しております。所属研究室サーバーで行った同様の解析(要素数16000000へ削減し8並列)では参考値と抵抗係数が合致しましたが、同様の解析条件を用いても要素数約1億の解析が発散してしまいます。以下に各種解析条件等を挙げますので、お忙しい中大変申し訳御座いませんが御教示頂けると幸いで御座います。
また、商用メッシュジェネレーターからOpenFOAM形式に出力する際、自動車のサイズにズレがあり、境界条件設定で通常ではあり得ない設定を行っております。詳細に追記ましては下部※1に記載致しますので併せて御確認下さい。また、これにつきましてもOpenFOAM内で修正する方法を御存知でしたら御教示頂けると幸いで御座います。
1. 各種サイズ (下部※1に補足あり)
1.1 自動車のサイズ
全長4000m、全幅2000m、全高1000m
1.2 解析領域のサイズ
各方向に15倍～20倍を確保
2. 各種境界条件
2.1 epsilon
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -3 0 0 0 0];
internalField uniform 3.68e+5;
boundaryField
{
inlet
{
type fixedValue;
value uniform 3.68e+5;
}
outlet
{
type zeroGradient;
}
sidewall
{
type slip;
}
floor
{
type epsilonWallFunction;
value uniform 3.68e+5;
}
car
{
type epsilonWallFunction;
value uniform 3.68e+5;
}
ceiling
{
type slip;
}
}
// ************************************************************************* //
2.2 k
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 1.5e+6;
boundaryField
{
inlet
{
type fixedValue;
value uniform 1.5e+6;
}
outlet
{
type zeroGradient;
}
sidewall
{
type slip;
}
floor
{
type kqRWallFunction;
value uniform 1.5e+06;
}
car
{
type kqRWallFunction;
value uniform 1.5e+06;
}
ceiling
{
type slip;
}
}
// ************************************************************************* //
2.3 nut
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
sidewall
{
type calculated;
value uniform 0;
}
floor
{
type nutkWallFunction;
value uniform 0;
}
car
{
type nutkWallFunction;
value uniform 0;
}
ceiling
{
type calculated;
value uniform 0;
}
}
// ************************************************************************* //
2.4 nuTilda
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];
internalField uniform 2.0e-7;
boundaryField
{
inlet
{
type fixedValue;
value uniform 2.0e-7;
}
outlet
{
type zeroGradient;
}
sidewall
{
type zeroGradient;
}
floor
{
type zeroGradient;
}
car
{
type zeroGradient;
}
ceiling
{
type zeroGradient;
}
}
// ************************************************************************* //
2.5 omega
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 -1 0 0 0 0];
internalField uniform 8.444;
boundaryField
{
inlet
{
type fixedValue;
value uniform 8.444;
}
outlet
{
type zeroGradient;
}
car
{
type omegaWallFunction;
value uniform 8.444;
}
sidewall
{
type slip;
}
floor
{
type omegaWallFunction;
value uniform 8.444;
}
ceiling
{
type slip;
}
}
// ************************************************************************* //
2.6 p
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
sidewall
{
type slip;
}
floor
{
type zeroGradient;
}
car
{
type zeroGradient;
}
ceiling
{
type slip;
}
}
// ************************************************************************* //
2.7 R
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform (0 0 0 0 0 0);
boundaryField
{
inlet
{
type fixedValue;
value uniform (0 0 0 0 0 0);
}
outlet
{
type zeroGradient;
}
sidewall
{
type slip;
}
floor
{
type kqRWallFunction;
value uniform ( 0 0 0 0 0 0 );
}
car
{
type kqRWallFunction;
value uniform ( 0 0 0 0 0 0 );
}
ceiling
{
type slip;
}
}
// ************************************************************************* //
2.8 U
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 -20000);
boundaryField
{
inlet
{
type fixedValue;
value uniform (0 0 -20000);
}
outlet
{
type zeroGradient;
}
sidewall
{
type slip;
}
floor
{
type fixedValue;
value uniform (0 0 0);
}
car
{
type fixedValue;
value uniform (0 0 0);
}
ceiling
{
type slip;
}
}
// ************************************************************************* //
3. fvSolution
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver GAMG;
tolerance 1e-8;
relTol 0.0001;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
cacheAgglomeration on;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 10;
mergeLevels 1;
};
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
};
k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
};
omega
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
};
epsilon
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
};
}
SIMPLE
{
nNonOrthogonalCorrectors 1;
pRefCell 0;
pRefPoint 0;
pRefValue 0;
residualControl
{
p 1e-8;
U 1e-8;
"(k|epsilon|omega)" 1e-8;
}
}
potentialFlow
{
nNonOrthogonalCorrectors 10;
}
relaxationFactors
{
p 0.01;
U 0.3;
k 0.3;
epsilon 0.3;
omega 0.3;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
4. fvScheme
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default CrankNicholson 0.3;
}
gradSchemes
{
default leastSquares QUICK;
grad(p) leastSquares QUICK;
grad(U) leastSquares QUICK;
//Gauss linear
}
divSchemes
{
default none;
div(phi,U) Gauss QUICK;
div(phi,k) Gauss QUICK;
div(phi,epsilon) Gauss QUICK;
div(phi,R) Gauss QUICK;
div(R) Gauss QUICK;
div(phi,nuTilda) Gauss QUICK;
div((nuEff*dev(T(grad(U))))) Gauss linear;
div(phi,omega) Gauss QUICK;
}
laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear limited 0.5;
laplacian((1|A(U)),p) Gauss linear limited 0.5;
laplacian(DkEff,k) Gauss linear limited 0.5;
laplacian(DepsilonEff,epsilon) Gauss linear limited 0.5;
laplacian(DREff,R) Gauss linear limited 0.5;
laplacian(DnuTildaEff,nuTilda) Gauss linear limited 0.5;
laplacian(DomegaEff,omega) Gauss linear limited 0.5;
}
interpolationSchemes
{
default upwind phi;
interpolate(U) upwind phi;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p ;
}
// ************************************************************************* //
5. RASProperties
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
RASModel kOmegaSST;
turbulence on;
printCoeffs on;
// ************************************************************************* //
6. transportProperties
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
transportModel Newtonian;
nu nu [ 0 2 -1 0 0 0 0 ] 15.11;
// ************************************************************************* //
※1
本来、自動車のサイズは全長4m、全幅2m、全高1mでありますが、前述の通り出力時にサイズが×10e+3されてしまう不具合を解消出来ておりません。同様の問題は今迄行って来た複数の解析ケースの内1つで発生しましたが、出力されるサイズをそのまま用い、流速に×10e+3した値を用いて乱流エネルギー等を計算し設定した所、抵抗係数は参考値と合致しました。この経験から、今回も各種境界条件の設定が通常ではあり得ない値と成って居ります。
Konishi
Dec 19, 2012, 5:09:31 AM12/19/12
to open...@googlegroups.com
Google Groups 「OpenFOAM」の皆様
お世話に成って居ります。konishiと申します。先ほどの投稿に関しまして追記致します。
解析ログが残って居なかった為、再度simpleFoamを実行している最中で御座います。ログが取れましたら再度投稿致します。
宜しくお願い申し上げます。
satoshi
Dec 19, 2012, 8:35:01 AM12/19/12
to open...@googlegroups.com
konishi様
onayaと申します。
解析設定はあまりよく見てませんが、
メッシュの単位変換(サイズ変更)は、
2012年12月19日水曜日 19時09分31秒 UTC+9 Konishi:
onayaと申します。
解析設定はあまりよく見てませんが、
メッシュの単位変換(サイズ変更)は、
transformPoints -scale "(1e-3 1e-3 1e-3)"
でできます。
また、fvSchemeを見るとスキームをQUICKを使用していますが、
計算の安定のため、初めはupwindを用いてはどうでしょうか？
また、初期値をpotentialFoamなどをもちいて作成して、
計算の安定化を図ってみてはどうでしょうか？
以上、
2012年12月19日水曜日 19時09分31秒 UTC+9 Konishi:
Konishi
Dec 20, 2012, 12:31:09 AM12/20/12
to open...@googlegroups.com
onaya様
御教示有難う御座います。konishiで御座います。
早速transformPointsを実行された所、メッシュの単位変換に成功致しました。有難う御座います。
スキームにつきましては、当初upwindで実行して居りましたが発散してしまった為、様々なスキームを模索していたという経緯が御座います。potentialFoamに関しましても解析が実行されません。次の投稿で各種ログを提示致します。
宜しくお願い申し上げます。
konishi
Konishi
Dec 20, 2012, 12:43:32 AM12/20/12
to open...@googlegroups.com
Google Groups 「OpenFOAM」の皆様
お世話に成って居ります。konishiで御座います。
本スレッドにおける解析のログを得ましたので提示致します。一つ目はpotentialFoam、二つ目はupwindでのsimpleFoam、三つ目はQUICKでのsimpleFoamに関するログで御座います。三つの解析全てにおきまして解析が発散若しくは原因不明の停止をしてしまっておりまして、提示致しましたログは解析停止5ステップ前からの物で御座います。
本解析の問題解決に関しまして御教示頂けると幸いで御座います。宜しくお願い申し上げます。
1. potentialFoam
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
relaxation factors:
fields:
{
p 0.3;
}
equations:
{
p 0.3;
U 0.7;
k 0.7;
epsilon 0.7;
omega 0.7;
}
Reading field p
Reading field U
Lookup interpolationScheme for interpolate(U)
Calculating potential flow
[38] [[43] [37] [44] [39] ####[42] [40] #00 0 [41] 0 0 ###0 0 0 36] #Foam:Foam::eFoam::erFoam::errFoam::erro:error:rror::priror::prinor::printr::printSFoam::error::printFoam::error::printFoam::error::printStackntStack(FoamtStack(Foam:Stack(Foam::Ostack(Foam::OstStack(Foam::OstrStack(Foam::O:printStack((Foam::Ostream::Ostream&):Ostream&)tream&)ream&)eam&)stream&)Foam::Os&)tream&)0 Foam::error::printStack(Foam::Ostream&)[83] [[84] [82] ###[88[86] 0 [85] ] 0 [89] 0 [87] ##0## 0 # 0 0 0 81] FoaFoaFoamm::em::e::erFoam::erFoam::rror::prror::ror::prror::perror::pFoam::error::printStack(FoaprintStack(FointStack(Foam:rintStack(FoarintStack(FoarintStack(FoaFoam::error::printStack(Foam::error::printStam::Ostream&)am::Ostream&):Ostream&)m::Ostream&)m::Ostream::OstreamFoam::Ostck(Foam::Ostm&)&)ream&)ream&)#0 Foam::error::printStack(Foam::Ostream&)[[73] [79] [77] [75] ####[74] #[80] #[78] 00 0 0 [76] 0 0 #0# 70 2] #FoaFoam:Foam::Foam::errm::error:error:or::pFoam::error::prinFoam::error::printSFoam::error::prin::printStack(Fo:printStack(FoamrintStack(Foam:tStack(Foam::tack(Foam::OtStack(Ferror::printStack(Foam::OstrFoam::error::printStack(Foamam::Ostream&)::Ostream&):Ostream&)Ostream&)stream&)oam::Ostream&eam&)::Ostream)&)0 Foam::error::printStack(Foam::Ostream&)[48] [52] [50] [47] [51] #[49] #[46] #####0[53] #0 0 [0 0 0 0 0 45] FoamFoamFoam::#::error::p::error::error::Foam::error::printFoam::error::prinFoam::error::prinFoam::error::rintStack(FoprintStaprintStacStack(FoamtStack(Foam:tStack(Foam:printStackFoam::error::printStackam::Ostream&)ck(Foam::Osk(Foam::Ostre::Ostream&):Ostream&):Ostream&)(Foam::Os(Foam::Osttream&)am&)treaream&m&))0 Foam::error::printStack(Foam::Ostream&)[1] [2] [4] [6][7]## #0[3] 0 0 # [5] ##0 #[8] 0 0 0 # 0 [0] FoaFoamFoam:Foam::em::error::::error::p:error::rror::Foam::erroFoam::error::Foam::error::prprintStack(FoarintStack(FoamprintStack(FoaprintStack(Fr::printStacprintStack(intStack(FoamFoam::error::printStackm::Ostream&)::Ostream&)m::Ostream&)oam::Ostreak(Foam::OstrFoam::Ost::Ostream&(Foam::Ostrm&)eam&)ream&))eam&)#0 Foam::error::printStack(Foam::Ostream&)[[25] [19] [20] [21] #[23] #[22] ###0#0[24] #0 [26] #0 0 0 0 0 18] #0FoaFoamm::e::erFoam::errFoam::error::pFoam::error::pFoam::error::prFoam::error::priFoam::error::printStarror::printStack(Foaror::printStack(Foaor::printStarintStack(Foam:rintStack(Foam:intStack(Foam::ntStack(Fock(Foam::Ostm::Ostream&m::Ostream&)ck(Foam::Ostre:Ostream&):Ostream&)Ostream&)am::Ostream&)ream&))am&) Foam::error::printStack(Foam::Ostream&)[30] [34] [31] [[29] #[28] #[33] #[32] ###0#[35] #0 0 0 0 0 0 0 27] #FoaFoamFoam:Foam::Foam::erm::error::p::error::pr:error::prierror::prror::priFoam::error::printStFoam::error::printStaFoam::error::printStrintStack(Foam::intStack(Foam::OstntStack(Foam::OstreaintStack(Foam::OntStack(Foam::Ostrack(Foam::Ostreack(Foam::Ostreamack(Foam::OstreOstream&)ream&)m&)stream&)eam&)m&)&)am&)0 Foam::error::printStack(Foam::Ostream&)[99] #0 Foam::error::printStack(Foam::Ostream&)[64] [69] [68] [70] [67] #[65] #####[66] #[71] #00 0 0 0 0 0 0 [63] FoaFoamFoam:Foam::m::err::err:errorerror:Foam::error::prFoam::error::prFoam::error::prFoam::error::por::printStaor::printStac::printStack:printStack(intStack(FoamintStack(FointStack(Foam:ck(Foam::k(Foam::Ostr(Foam::OstreFoam::Ostream::Ostream&)am::Ostream:Ostream&)rintStack(Foam::OstrOstream&eam&)am&)&)&)eam&))#0 Foam::error::printStack(Foam::Ostream&)[12] [15] [16] [11] ####[14] [10] [13] [17] 0## 00 #00 #0 0 0 [9] FoaFoam::eFoam::eFoam::erroFoam::errorm::error::rror::prinrror::prir::printSt::printStaFoam::error::printStaFoam::error::printStFoam::error::printStaprintStack(FotStack(Foam::ntStack(Foamack(Foam::Ostreck(Foam::Ostreack(Foam::Ostreack(Foam::Ostreamck(Foam::Ostream&)am::Ostream&)Ostream&)::Ostream&)am&)m&)am&)&)#0 Foam::error::printStack(Foam::Ostream&)[96] [95] [93] [92] [91] ####[94] [98] [97] 0 0 #00 0 #0 ## 00 [90] FoaFoamFoam:m::er::err:erroFoam::errFoam::error:Foam::error::pFoam::error::priFoam::error::prinror::printStack(For::printStack(Foamr::printStack(Foamor::printStack(F:printStack(ForintStack(FontStack(FoamtStack(Foam::oam::Ostr::Ostream&)::Ostream&)oam::Ostreamam::Ostream&am::Ostream&::Ostream&)Ostream&)eam&)&)))#0 Foam::error::printStack(Foam::Ostream&)[55] [61] [[59] #[56] #[60] #[57] #[58] ###[62] #0 0 00 0 0 0 0 54] #FoaFoamFoam:Foam::em::error::p::error::pr:error::rror::pFoam::error::pFoam::error::pFoam::error::prinFoam::error::pririntStack(FoaintStack(FoamprintStackrintStack(rintStack(ForintStack(FoatStack(Foam::OsntStack(Foam::m::Ostream&)::Ostream&)(Foam::Ostream&)Foam::Ostreamam::Ostream&)m::Ostream&)tream&)Ostream&)&)0 Foam::error::printStack(Foam::Ostream&) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[99] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[99] #2 _IO_funlockfile in "/opt/FJSVxosmmm/lib64/libmpgpthread.so.1"
[99] #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>&) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libfiniteVolume.so"
[99] #4 main in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/bin/potentialFoam"
[99] #5 __libc_start_main in "/lib64/libc.so.6"
[99] #6 in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[31] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[65] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[26] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[52] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[54] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[79] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[42] #1 Foam::sigFpe::sigHandler(int) in in "/u"/ussrr/lo/lcalocal/Op/OpenFOAMenF/2.1.OAM0//2Ope.1.0nFO/OAM-2.1.penFOA0/pM-2lat.1.0/fplorms/S6atfo4FrmsXFcc/S64FDPOXFcpt/cDlib/libPOptOpenF/libOAM.s/libo"Ope
nFO[39]AM. #so"1
[40]Foam #::sigF1 peFo::sam:igH:sanigFdleper(::sintig)Handler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[12] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[84] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[85] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[96] #1 Foam::sigFpe::sigHandler(int)_start in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[6] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[25] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[74] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[55] #1 Foam::sigFpe::sigHandler(int) in in in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[13] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[15] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[93] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[80] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[56] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[59] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[44] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/ in in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[3] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[70] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[19] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[34] #1 Foam::sigFpe::sigHandler(int) in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[53] #1 Foam::sigFpe::sigHandler(int) in in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[9] #1 Foam::sigFpe::sigHandler(int)platforms/S64FXFccDPOpt/bin/potentialFoam"
[b03t11058:13410] *** Process received signal ***
[b03t11058:13410] Signal: Floating point exception (8)
[b03t11058:13410] Signal code: (-6)
[b03t11058:13410] Failing at address: 0x346200008c34
[b03t11058:13410] [ 0] /opt/FJSVxosmmm/lib64/libmpgpthread.so.1(_IO_funlockfile+0x5c) [0xffffffff0a4b1974]
[b03t11058:13410] [ 1] /opt/FJSVxosmmm/lib64/libmpgpthread.so.1(raise+0x34) [0xffffffff0a4b17cc]
[b03t11058:13410] [ 2] /usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6sigFpe10sigHandlerEi+0xfc) [0xffffffff086bf23c]
[b03t11058:13410] [ 3] /opt/FJSVxosmmm/lib64/libmpgpthread.so.1(_IO_funlockfile+0x5c) [0xffffffff0a4b1974]
[b03t11058:13410] [ 4] /usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libfiniteVolume.so(_ZN4Foam9adjustPhiERNS_14GeometricFieldIdNS_13fvsPatchFieldENS_11surfaceMeshEEERKNS0_INS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERNS0_IdS7_S8_EE+0x7c0) [0xffffffff05c3da40]
[b03t11058:13410] [ 5] potentialFoam(main+0x1734) [0x159e74]
[b03t11058:13410] [ 6] /lib64/libc.so.6(__libc_start_main+0x194) [0xffffffff0ae3381c]
[b03t11058:13410] [ 7] potentialFoam [0x1581ac]
[b03t11058:13410] *** End of error message ***
"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[97] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[95] #1 Foam::sigFpe::sigHandler(int) in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[5] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[61] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[62] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[20] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[18] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[71] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[33] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[35] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[27] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[32] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[30] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[41] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[38] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[51] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[89] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[83] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[86] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[78] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[77] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[90] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[64] #1 Foam::sigFpe::sigHandler(int) in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[48] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[47] #1 Foam::sigFpe::sigHandler(int) in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[73] #1 Foam::sigFpe::sigHandler(int) in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[68] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[45] #1 Foam::sigFpe::sigHandler(int) in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[2] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[50] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[49] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[46] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[66] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOA in in M/2.1.0/Op"/uenFsr/lOAMoc-2.1"/usr.0//loplcaal/OpenFatforms/Sl/OpOAM64FenFOAM/2.1.X/2.1.0/OpeFccD0/OpenPOFOAMnFOAM--2.2.1.01.0/p/platforpt/lib/llams/ibOtforS64FpenFms/XFcOAMS64FcDPO.so"XFcpt/cDPlib
[Opt/lib/l/libOpenFOA16] #1 ibOpenFOM.so"AM.sFoam::
[17o"
sigFpe::si] #gHa1 [10] #1ndler(int )Foam::Foam::sigFpe:sigFpe::s:sigHandigHlerand(iler(innt)t) in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[76] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[21] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[11] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[24] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[1] #1 Foam::sigFpe::sigHandler(int) in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[7] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[63] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[87] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[98] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[92] #1 Foam::sigFpe::sigHandler(int) in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[94] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[14] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[81] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[60] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[58] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFO inA M/2.1.0/OpenFOAM"/-usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
2[37] #.1 1Foa.m::sigFpe::sigHandler(int)0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[36] #1 Foam::sigFpe::sigHandler(int) in in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[88] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[43] #1 Foam::sigFpe::sigHandler(int) in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[8] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[22] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[23] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[57] #1 Foam::sigFpe::sigHandler(int) i in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[75] #1 Foam::sigFpe::sigHandler(int)n "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[72] #1 Foam::sigFpe::sigHandler(int) in in in in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[69] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[28] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[4] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[91] #1 Foam::sigFpe::sigHandler(int)"/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[67] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/platforms/S64FXFccDPOpt/lib/libOpenFOAM.so"
[29] #1 Foam::sigFpe::sigHandler(int)[ERR.] PLE 0017 plexec The process terminated with the signal.(rank=99)(nid=0x010c0011)(sig=8)
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
2. upwindによるsimpleFoam
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Time = 8.5
Lookup divScheme for div(phi,U)
Lookup laplacianScheme for laplacian(nuEff,U)
Lookup fluxRequired for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 93132
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 93132
Lookup divScheme for div((nuEff*dev(T(grad(U)))))
Find equation relaxation factor for U
Lookup equation relaxation factor for U
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 93102
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for U
DILUPBiCG: Solving for Ux, Initial residual = 0.323156, Final residual = 0.00068826, No Iterations 7
DILUPBiCG: Solving for Uy, Initial residual = 0.354448, Final residual = 0.000432722, No Iterations 7
DILUPBiCG: Solving for Uz, Initial residual = 0.371166, Final residual = 0.00050478, No Iterations 7
Lookup interpolationScheme for interpolate(HbyA)
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 93915
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.815898, Final residual = 0.0075377, No Iterations 13
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 94063
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.551099, Final residual = 0.00530596, No Iterations 11
Lookup fluxRequired for p
time step continuity errors : sum local = 1.02793, global = -7.59921e-15, cumulative = -7.82499e-16
Find variable relaxation factor for p
Lookup variable relaxation factor for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 94217
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 94247
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 93621
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 94286
Lookup ddtScheme for ddt(omega)
Lookup divScheme for div(phi,omega)
Lookup laplacianScheme for laplacian(DomegaEff,omega)
Lookup fluxRequired for omega
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 94286
Find equation relaxation factor for omega
Lookup equation relaxation factor for omega
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for omega
DILUPBiCG: Solving for omega, Initial residual = 0.538462, Final residual = 1.37308e-19, No Iterations 1
bounding omega, min: -2.85354e+19 max: 4.21131e+49 average: 6.48688e+41
Lookup ddtScheme for ddt(k)
Lookup divScheme for div(phi,k)
Lookup laplacianScheme for laplacian(DkEff,k)
Lookup fluxRequired for k
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 93621
Find equation relaxation factor for k
Lookup equation relaxation factor for k
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for k
DILUPBiCG: Solving for k, Initial residual = 5.75753e-09, Final residual = 5.75753e-09, No Iterations 0
ExecutionTime = 7466.42 s ClockTime = 7483 s
Find variable relaxation factor for nu
Find variable relaxation factor for y
Find variable relaxation factor for k
Find variable relaxation factor for k_0
Find variable relaxation factor for ddt0(omega)
Find variable relaxation factor for omega_0_0
Find variable relaxation factor for ddt0(k)
Find variable relaxation factor for omega_0
Find variable relaxation factor for p
Find variable relaxation factor for omega
Find variable relaxation factor for nut
Find variable relaxation factor for k_0_0
Find variable relaxation factor for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 94247
forceCoeffs output:
Cd = 1.1859e+15
Cl = -9.75653e+14
Cm = -1.48643e+14
Time = 8.6
Lookup divScheme for div(phi,U)
Lookup laplacianScheme for laplacian(nuEff,U)
Lookup fluxRequired for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 94247
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 94247
Lookup divScheme for div((nuEff*dev(T(grad(U)))))
Find equation relaxation factor for U
Lookup equation relaxation factor for U
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 94217
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for U
DILUPBiCG: Solving for Ux, Initial residual = 0.36185, Final residual = 0.00185434, No Iterations 5
DILUPBiCG: Solving for Uy, Initial residual = 0.578874, Final residual = 0.00403744, No Iterations 5
DILUPBiCG: Solving for Uz, Initial residual = 0.577615, Final residual = 0.00466314, No Iterations 5
Lookup interpolationScheme for interpolate(HbyA)
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 95030
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.778917, Final residual = 0.00769274, No Iterations 11
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 95178
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.488951, Final residual = 0.00472946, No Iterations 9
Lookup fluxRequired for p
time step continuity errors : sum local = 0.79329, global = 1.3102e-15, cumulative = 5.27703e-16
Find variable relaxation factor for p
Lookup variable relaxation factor for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 95332
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 95362
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 94736
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 95401
Lookup ddtScheme for ddt(omega)
Lookup divScheme for div(phi,omega)
Lookup laplacianScheme for laplacian(DomegaEff,omega)
Lookup fluxRequired for omega
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 95401
Find equation relaxation factor for omega
Lookup equation relaxation factor for omega
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for omega
DILUPBiCG: Solving for omega, Initial residual = 0.538462, Final residual = 7.30268e-17, No Iterations 1
bounding omega, min: -4.83629e+16 max: 1.26339e+49 average: 1.94606e+41
Lookup ddtScheme for ddt(k)
Lookup divScheme for div(phi,k)
Lookup laplacianScheme for laplacian(DkEff,k)
Lookup fluxRequired for k
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 94736
Find equation relaxation factor for k
Lookup equation relaxation factor for k
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for k
DILUPBiCG: Solving for k, Initial residual = 5.75753e-09, Final residual = 5.75753e-09, No Iterations 0
ExecutionTime = 7549.86 s ClockTime = 7567 s
Find variable relaxation factor for nu
Find variable relaxation factor for y
Find variable relaxation factor for k
Find variable relaxation factor for k_0
Find variable relaxation factor for ddt0(omega)
Find variable relaxation factor for omega_0_0
Find variable relaxation factor for ddt0(k)
Find variable relaxation factor for omega_0
Find variable relaxation factor for p
Find variable relaxation factor for omega
Find variable relaxation factor for nut
Find variable relaxation factor for k_0_0
Find variable relaxation factor for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 95362
forceCoeffs output:
Cd = 5.88604e+14
Cl = -4.71925e+14
Cm = -6.73357e+13
Time = 8.7
Lookup divScheme for div(phi,U)
Lookup laplacianScheme for laplacian(nuEff,U)
Lookup fluxRequired for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 95362
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 95362
Lookup divScheme for div((nuEff*dev(T(grad(U)))))
Find equation relaxation factor for U
Lookup equation relaxation factor for U
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 95332
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for U
DILUPBiCG: Solving for Ux, Initial residual = 0.42731, Final residual = 0.000906428, No Iterations 6
DILUPBiCG: Solving for Uy, Initial residual = 0.517154, Final residual = 0.00156311, No Iterations 6
DILUPBiCG: Solving for Uz, Initial residual = 0.51426, Final residual = 0.00385563, No Iterations 5
Lookup interpolationScheme for interpolate(HbyA)
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 96145
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.739858, Final residual = 0.00680917, No Iterations 11
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 96293
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.575671, Final residual = 0.00502396, No Iterations 6
Lookup fluxRequired for p
time step continuity errors : sum local = 1.01787, global = -1.92081e-15, cumulative = -1.39311e-15
Find variable relaxation factor for p
Lookup variable relaxation factor for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 96447
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 96477
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 95851
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 96516
Lookup ddtScheme for ddt(omega)
Lookup divScheme for div(phi,omega)
Lookup laplacianScheme for laplacian(DomegaEff,omega)
Lookup fluxRequired for omega
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 96516
Find equation relaxation factor for omega
Lookup equation relaxation factor for omega
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for omega
DILUPBiCG: Solving for omega, Initial residual = 0.538462, Final residual = 2.80877e-18, No Iterations 1
bounding omega, min: -1.5898e+18 max: 3.79018e+48 average: 5.83819e+40
Lookup ddtScheme for ddt(k)
Lookup divScheme for div(phi,k)
Lookup laplacianScheme for laplacian(DkEff,k)
Lookup fluxRequired for k
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 95851
Find equation relaxation factor for k
Lookup equation relaxation factor for k
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for k
DILUPBiCG: Solving for k, Initial residual = 5.75753e-09, Final residual = 5.75753e-09, No Iterations 0
ExecutionTime = 7628.41 s ClockTime = 7646 s
Find variable relaxation factor for nu
Find variable relaxation factor for y
Find variable relaxation factor for k
Find variable relaxation factor for k_0
Find variable relaxation factor for ddt0(omega)
Find variable relaxation factor for omega_0_0
Find variable relaxation factor for ddt0(k)
Find variable relaxation factor for omega_0
Find variable relaxation factor for p
Find variable relaxation factor for omega
Find variable relaxation factor for nut
Find variable relaxation factor for k_0_0
Find variable relaxation factor for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 96477
forceCoeffs output:
Cd = 2.98948e+14
Cl = -2.5483e+14
Cm = -4.09709e+13
Time = 8.8
Lookup divScheme for div(phi,U)
Lookup laplacianScheme for laplacian(nuEff,U)
Lookup fluxRequired for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 96477
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 96477
Lookup divScheme for div((nuEff*dev(T(grad(U)))))
Find equation relaxation factor for U
Lookup equation relaxation factor for U
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 96447
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for U
DILUPBiCG: Solving for Ux, Initial residual = 0.341987, Final residual = 0.00154238, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.446989, Final residual = 0.00312287, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.3851, Final residual = 0.00252988, No Iterations 4
Lookup interpolationScheme for interpolate(HbyA)
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 97260
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.81282, Final residual = 0.00769668, No Iterations 13
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 97408
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.574601, Final residual = 0.005717, No Iterations 5
Lookup fluxRequired for p
time step continuity errors : sum local = 1.92299, global = -7.70046e-15, cumulative = -9.09357e-15
Find variable relaxation factor for p
Lookup variable relaxation factor for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 97562
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 97592
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 96966
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 97631
Lookup ddtScheme for ddt(omega)
Lookup divScheme for div(phi,omega)
Lookup laplacianScheme for laplacian(DomegaEff,omega)
Lookup fluxRequired for omega
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 97631
Find equation relaxation factor for omega
Lookup equation relaxation factor for omega
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for omega
DILUPBiCG: Solving for omega, Initial residual = 0.538462, Final residual = 6.95815e-17, No Iterations 1
bounding omega, min: -6.1885e+16 max: 1.13705e+48 average: 1.75146e+40
Lookup ddtScheme for ddt(k)
Lookup divScheme for div(phi,k)
Lookup laplacianScheme for laplacian(DkEff,k)
Lookup fluxRequired for k
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 96966
Find equation relaxation factor for k
Lookup equation relaxation factor for k
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for k
DILUPBiCG: Solving for k, Initial residual = 5.75753e-09, Final residual = 5.75753e-09, No Iterations 0
ExecutionTime = 7704.49 s ClockTime = 7722 s
Find variable relaxation factor for nu
Find variable relaxation factor for y
Find variable relaxation factor for k
Find variable relaxation factor for k_0
Find variable relaxation factor for ddt0(omega)
Find variable relaxation factor for omega_0_0
Find variable relaxation factor for ddt0(k)
Find variable relaxation factor for omega_0
Find variable relaxation factor for p
Find variable relaxation factor for omega
Find variable relaxation factor for nut
Find variable relaxation factor for k_0_0
Find variable relaxation factor for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 97592
forceCoeffs output:
Cd = -4.8483e+14
Cl = 3.69579e+14
Cm = 4.73914e+13
Time = 8.9
Lookup divScheme for div(phi,U)
Lookup laplacianScheme for laplacian(nuEff,U)
Lookup fluxRequired for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 97592
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 97592
Lookup divScheme for div((nuEff*dev(T(grad(U)))))
Find equation relaxation factor for U
Lookup equation relaxation factor for U
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 97562
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for U
DILUPBiCG: Solving for Ux, Initial residual = 0.304593, Final residual = 0.00182562, No Iterations 5
DILUPBiCG: Solving for Uy, Initial residual = 0.417584, Final residual = 0.000345409, No Iterations 6
DILUPBiCG: Solving for Uz, Initial residual = 0.429483, Final residual = 0.00337662, No Iterations 5
Lookup interpolationScheme for interpolate(HbyA)
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 98375
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
3. QUICKによるsimpleFoam
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Time = 9.6
Lookup divScheme for div(phi,U)
Lookup gradScheme for grad(magSqr(U))
Cache: Calculating grad(magSqr(U)), magSqr(U) event No. 116193
Lookup laplacianScheme for laplacian(nuEff,U)
Lookup fluxRequired for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 115546
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 115546
Lookup divScheme for div((nuEff*dev(T(grad(U)))))
Find equation relaxation factor for U
Lookup equation relaxation factor for U
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 115516
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for U
DILUPBiCG: Solving for Ux, Initial residual = 0.194504, Final residual = 0.000164898, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.193035, Final residual = 0.000163908, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.19177, Final residual = 0.000161784, No Iterations 2
Lookup interpolationScheme for interpolate(HbyA)
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 116419
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.942304, Final residual = 6.29957e-05, No Iterations 48
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 116567
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.488158, Final residual = 4.05629e-05, No Iterations 21
Lookup fluxRequired for p
time step continuity errors : sum local = 680.526, global = 4.45855e-07, cumulative = 4.21319e-07
Find variable relaxation factor for p
Lookup variable relaxation factor for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 116721
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 116751
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 116095
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 116790
Lookup ddtScheme for ddt(omega)
Lookup divScheme for div(phi,omega)
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 116790
Lookup laplacianScheme for laplacian(DomegaEff,omega)
Lookup fluxRequired for omega
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 116790
Find equation relaxation factor for omega
Lookup equation relaxation factor for omega
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for omega
DILUPBiCG: Solving for omega, Initial residual = 0.176471, Final residual = 5.67719e-21, No Iterations 1
bounding omega, min: -1.68263e+41 max: 3.33112e+58 average: 3.63398e+50
Lookup ddtScheme for ddt(k)
Lookup divScheme for div(phi,k)
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 116095
Lookup laplacianScheme for laplacian(DkEff,k)
Lookup fluxRequired for k
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 116095
Find equation relaxation factor for k
Lookup equation relaxation factor for k
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for k
DILUPBiCG: Solving for k, Initial residual = 6.65244e-06, Final residual = 6.65244e-06, No Iterations 0
ExecutionTime = 18188.3 s ClockTime = 18212 s
Find variable relaxation factor for nu
Find variable relaxation factor for y
Find variable relaxation factor for k
Find variable relaxation factor for k_0
Find variable relaxation factor for ddt0(omega)
Find variable relaxation factor for omega_0_0
Find variable relaxation factor for ddt0(k)
Find variable relaxation factor for omega_0
Find variable relaxation factor for p
Find variable relaxation factor for omega
Find variable relaxation factor for nut
Find variable relaxation factor for k_0_0
Find variable relaxation factor for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 116751
forceCoeffs output:
Cd = -3.89745e+23
Cl = 4.42125e+23
Cm = -1.97373e+23
Time = 9.7
Lookup divScheme for div(phi,U)
Lookup gradScheme for grad(magSqr(U))
Cache: Calculating grad(magSqr(U)), magSqr(U) event No. 117398
Lookup laplacianScheme for laplacian(nuEff,U)
Lookup fluxRequired for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 116751
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 116751
Lookup divScheme for div((nuEff*dev(T(grad(U)))))
Find equation relaxation factor for U
Lookup equation relaxation factor for U
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 116721
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for U
DILUPBiCG: Solving for Ux, Initial residual = 0.195255, Final residual = 0.000167527, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.19375, Final residual = 0.000166544, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.192335, Final residual = 0.000164394, No Iterations 2
Lookup interpolationScheme for interpolate(HbyA)
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 117624
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.961468, Final residual = 6.44001e-05, No Iterations 31
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 117772
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.535971, Final residual = 3.97685e-05, No Iterations 23
Lookup fluxRequired for p
time step continuity errors : sum local = 720.531, global = 2.34729e-07, cumulative = 6.56048e-07
Find variable relaxation factor for p
Lookup variable relaxation factor for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 117926
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 117956
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 117300
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 117995
Lookup ddtScheme for ddt(omega)
Lookup divScheme for div(phi,omega)
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 117995
Lookup laplacianScheme for laplacian(DomegaEff,omega)
Lookup fluxRequired for omega
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 117995
Find equation relaxation factor for omega
Lookup equation relaxation factor for omega
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for omega
DILUPBiCG: Solving for omega, Initial residual = 0.176471, Final residual = 5.96788e-21, No Iterations 1
bounding omega, min: -1.82869e+46 max: 2.33178e+58 average: 2.54568e+50
Lookup ddtScheme for ddt(k)
Lookup divScheme for div(phi,k)
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 117300
Lookup laplacianScheme for laplacian(DkEff,k)
Lookup fluxRequired for k
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 117300
Find equation relaxation factor for k
Lookup equation relaxation factor for k
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for k
DILUPBiCG: Solving for k, Initial residual = 6.65244e-06, Final residual = 6.65244e-06, No Iterations 0
ExecutionTime = 18323.1 s ClockTime = 18347 s
Find variable relaxation factor for nu
Find variable relaxation factor for y
Find variable relaxation factor for k
Find variable relaxation factor for k_0
Find variable relaxation factor for ddt0(omega)
Find variable relaxation factor for omega_0_0
Find variable relaxation factor for ddt0(k)
Find variable relaxation factor for omega_0
Find variable relaxation factor for p
Find variable relaxation factor for omega
Find variable relaxation factor for nut
Find variable relaxation factor for k_0_0
Find variable relaxation factor for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 117956
forceCoeffs output:
Cd = -8.52386e+23
Cl = 9.66362e+23
Cm = -4.31793e+23
Time = 9.8
Lookup divScheme for div(phi,U)
Lookup gradScheme for grad(magSqr(U))
Cache: Calculating grad(magSqr(U)), magSqr(U) event No. 118603
Lookup laplacianScheme for laplacian(nuEff,U)
Lookup fluxRequired for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 117956
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 117956
Lookup divScheme for div((nuEff*dev(T(grad(U)))))
Find equation relaxation factor for U
Lookup equation relaxation factor for U
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 117926
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for U
DILUPBiCG: Solving for Ux, Initial residual = 0.194292, Final residual = 0.000164749, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.192853, Final residual = 0.000163788, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.19165, Final residual = 0.000161875, No Iterations 2
Lookup interpolationScheme for interpolate(HbyA)
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 118829
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.946739, Final residual = 4.81717e+15, No Iterations 1000
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 118977
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.803774, Final residual = 0.197134, No Iterations 1000
Lookup fluxRequired for p
time step continuity errors : sum local = 7.9708e+21, global = -2.51477e+08, cumulative = -2.51477e+08
Find variable relaxation factor for p
Lookup variable relaxation factor for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 119131
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 119161
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 118505
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 119200
Lookup ddtScheme for ddt(omega)
Lookup divScheme for div(phi,omega)
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 119200
Lookup laplacianScheme for laplacian(DomegaEff,omega)
Lookup fluxRequired for omega
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 119200
Find equation relaxation factor for omega
Lookup equation relaxation factor for omega
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for omega
DILUPBiCG: Solving for omega, Initial residual = 0.176471, Final residual = 2.38579e-17, No Iterations 1
bounding omega, min: -3.66234e+46 max: 1.63225e+58 average: 1.78203e+50
Lookup ddtScheme for ddt(k)
Lookup divScheme for div(phi,k)
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 118505
Lookup laplacianScheme for laplacian(DkEff,k)
Lookup fluxRequired for k
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 118505
Find equation relaxation factor for k
Lookup equation relaxation factor for k
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for k
DILUPBiCG: Solving for k, Initial residual = 6.65249e-06, Final residual = 6.65249e-06, No Iterations 0
ExecutionTime = 21958.7 s ClockTime = 21983 s
Find variable relaxation factor for nu
Find variable relaxation factor for y
Find variable relaxation factor for k
Find variable relaxation factor for k_0
Find variable relaxation factor for ddt0(omega)
Find variable relaxation factor for omega_0_0
Find variable relaxation factor for ddt0(k)
Find variable relaxation factor for omega_0
Find variable relaxation factor for p
Find variable relaxation factor for omega
Find variable relaxation factor for nut
Find variable relaxation factor for k_0_0
Find variable relaxation factor for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 119161
forceCoeffs output:
Cd = 2.31461e+28
Cl = 5.51246e+29
Cm = -5.08054e+28
Time = 9.9
Lookup divScheme for div(phi,U)
Lookup gradScheme for grad(magSqr(U))
Cache: Calculating grad(magSqr(U)), magSqr(U) event No. 119808
Lookup laplacianScheme for laplacian(nuEff,U)
Lookup fluxRequired for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 119161
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 119161
Lookup divScheme for div((nuEff*dev(T(grad(U)))))
Find equation relaxation factor for U
Lookup equation relaxation factor for U
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 119131
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for U
DILUPBiCG: Solving for Ux, Initial residual = 0.22276, Final residual = 0.0010759, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.224566, Final residual = 0.00147002, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.217666, Final residual = 0.000526327, No Iterations 2
Lookup interpolationScheme for interpolate(HbyA)
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 120034
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 1, Final residual = 8.38685e-05, No Iterations 28
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 120182
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.68524, Final residual = 5.24467e-05, No Iterations 18
Lookup fluxRequired for p
time step continuity errors : sum local = 1.8957e+16, global = -453.082, cumulative = -2.51478e+08
Find variable relaxation factor for p
Lookup variable relaxation factor for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 120336
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 120366
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 119710
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 120405
Lookup ddtScheme for ddt(omega)
Lookup divScheme for div(phi,omega)
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 120405
Lookup laplacianScheme for laplacian(DomegaEff,omega)
Lookup fluxRequired for omega
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 120405
Find equation relaxation factor for omega
Lookup equation relaxation factor for omega
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for omega
DILUPBiCG: Solving for omega, Initial residual = 0.176471, Final residual = 1.04307e-22, No Iterations 1
bounding omega, min: -7.7595e+50 max: 1.14257e+58 average: 1.24927e+50
Lookup ddtScheme for ddt(k)
Lookup divScheme for div(phi,k)
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 119710
Lookup laplacianScheme for laplacian(DkEff,k)
Lookup fluxRequired for k
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 119710
Find equation relaxation factor for k
Lookup equation relaxation factor for k
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for k
DILUPBiCG: Solving for k, Initial residual = 6.65249e-06, Final residual = 6.65249e-06, No Iterations 0
ExecutionTime = 22081 s ClockTime = 22106 s
Find variable relaxation factor for nu
Find variable relaxation factor for y
Find variable relaxation factor for k
Find variable relaxation factor for k_0
Find variable relaxation factor for ddt0(omega)
Find variable relaxation factor for omega_0_0
Find variable relaxation factor for ddt0(k)
Find variable relaxation factor for omega_0
Find variable relaxation factor for p
Find variable relaxation factor for omega
Find variable relaxation factor for nut
Find variable relaxation factor for k_0_0
Find variable relaxation factor for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 120366
forceCoeffs output:
Cd = 5.49112e+45
Cl = 4.49192e+42
Cm = 1.56355e+45
Time = 10
Lookup divScheme for div(phi,U)
Lookup gradScheme for grad(magSqr(U))
Cache: Calculating grad(magSqr(U)), magSqr(U) event No. 121013
Lookup laplacianScheme for laplacian(nuEff,U)
Lookup fluxRequired for U
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 120366
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 120366
Lookup divScheme for div((nuEff*dev(T(grad(U)))))
Find equation relaxation factor for U
Lookup equation relaxation factor for U
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 120336
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for U
DILUPBiCG: Solving for Ux, Initial residual = 0.297636, Final residual = 0.00043465, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.27013, Final residual = 0.000425807, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.281012, Final residual = 0.00043721, No Iterations 3
Lookup interpolationScheme for interpolate(HbyA)
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 121239
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.634057, Final residual = 5.57307e-05, No Iterations 24
Lookup laplacianScheme for laplacian((1|A(U)),p)
Lookup fluxRequired for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 121387
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for p
GAMG: Solving for p, Initial residual = 0.317457, Final residual = 2.68174e-05, No Iterations 23
Lookup fluxRequired for p
time step continuity errors : sum local = 4.16185e+15, global = -5249.27, cumulative = -2.51483e+08
Find variable relaxation factor for p
Lookup variable relaxation factor for p
Lookup gradScheme for grad(p)
Cache: Calculating grad(p), p event No. 121541
Lookup gradScheme for grad(U)
Cache: Calculating grad(U), U event No. 121571
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 120915
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 121610
Lookup ddtScheme for ddt(omega)
Lookup divScheme for div(phi,omega)
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 121610
Lookup laplacianScheme for laplacian(DomegaEff,omega)
Lookup fluxRequired for omega
Lookup gradScheme for grad(omega)
Cache: Calculating grad(omega), omega event No. 121610
Find equation relaxation factor for omega
Lookup equation relaxation factor for omega
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for omega
DILUPBiCG: Solving for omega, Initial residual = 0.176471, Final residual = 3.87951e-19, No Iterations 1
bounding omega, min: -5.33538e+42 max: 7.99801e+57 average: 8.74491e+49
Lookup ddtScheme for ddt(k)
Lookup divScheme for div(phi,k)
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 120915
Lookup laplacianScheme for laplacian(DkEff,k)
Lookup fluxRequired for k
Lookup gradScheme for grad(k)
Cache: Calculating grad(k), k event No. 120915
Find equation relaxation factor for k
Lookup equation relaxation factor for k
From function solution::solverDict(const word&)
in file matrices/solution/solution.C at line 353
Lookup solver for k
DILUPBiCG: Solving for k, Initial residual = 6.65249e-06, Final residual = 6.65249e-06, No Iterations 0
[79]
[79]
[79] --> FOAM FATAL IO ERROR:
[79] error in IOstream "/home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor79/10/omega" for operation Ostream& operator<<(Ostream&, const Scalar&)
[79]
[79] file: /home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor79/10/omega at line 749688.
[79]
[79] From function IOstream::check(const char*) const
[79] in file db/IOstreams/IOstreams/IOstream.C at line 96.
[79]
FOAM parallel run exiting
[79]
[34]
[34]
[34] --> FOAM FATAL IO ERROR:
[34] error in IOstream "/home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor34/10/omega" for operation Ostream& operator<<(Ostream&, const Scalar&)
[34]
[34] file: /home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor34/10/omega at line 588725.
[34]
[34] From function IOstream::check(const char*) const
[34] in file db/IOstreams/IOstreams/IOstream.C at line 96.
[34]
FOAM parallel run exiting
[34]
[17]
[17]
[17] --> FOAM FATAL IO ERROR:
[17] error in IOstream "/home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor17/10/omega" for operation Ostream& operator<<(Ostream&, const Scalar&)
[17]
[17] file: /home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor17/10/omega at line 396901.
[17]
[17] From function IOstream::check(const char*) const
[17] in file db/IOstreams/IOstreams/IOstream.C at line 96.
[17]
FOAM parallel run exiting
[17]
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 17 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.
--------------------------------------------------------------------------
[70]
[70]
[70] --> FOAM FATAL IO ERROR:
[70] error in IOstream "/home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor70/10/omega" for operation Ostream& operator<<(Ostream&, const Scalar&)
[70]
[70] file: /home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor70/10/omega at line 596554.
[70]
[70] From function IOstream::check(const char*) const
[70] in file db/IOstreams/IOstreams/IOstream.C at line 96.
[70]
FOAM parallel run exiting
[70]
[44]
[44]
[44] --> FOAM FATAL IO ERROR:
[44] error in IOstream "/home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor44/10/omega" for operation Ostream& operator<<(Ostream&, const Scalar&)
[44]
[44] file: /home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor44/10/omega at line 588401.
[44]
[44] From function IOstream::check(const char*) const
[44] in file db/IOstreams/IOstreams/IOstream.C at line 96.
[44]
FOAM parallel run exiting
[44]
[99]
[99]
[99] --> FOAM FATAL IO ERROR:
[99] error in IOstream "/home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor99/10/omega" for operation Ostream& operator<<(Ostream&, const Scalar&)
[99]
[99] file: /home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor99/10/omega at line 837904.
[99]
[99] From function IOstream::check(const char*) const
[99] in file db/IOstreams/IOstreams/IOstream.C at line 96.
[99]
FOAM parallel run exiting
[99]
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 99 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.
--------------------------------------------------------------------------
[88]
[88]
[88] --> FOAM FATAL IO ERROR:
[88] error in IOstream "/home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor88/10/omega" for operation Ostream& operator<<(Ostream&, const Scalar&)
[88]
[88] file: /home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor88/10/omega at line 404199.
[88]
[88] From function IOstream::check(const char*) const
[88] in file db/IOstreams/IOstreams/IOstream.C at line 96.
[88]
FOAM parallel run exiting
[88]
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 88 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.
--------------------------------------------------------------------------
[47]
[47]
[47] --> FOAM FATAL IO ERROR:
[47] error in IOstream "/home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor47/10/omega" for operation Ostream& operator<<(Ostream&, const Scalar&)
[47]
[47] file: /home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor47/10/omega at line 588252.
[47]
[47] From function IOstream::check(const char*) const
[47] in file db/IOstreams/IOstreams/IOstream.C at line 96.
[47]
FOAM parallel run exiting
[47]
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 47 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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 79 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.
--------------------------------------------------------------------------
[97]
[97]
[97] --> FOAM FATAL IO ERROR:
[97] error in IOstream "/home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor97/10/omega" for operation Ostream& operator<<(Ostream&, const Scalar&)
[97]
[97] file: /home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor97/10/omega at line 485142.
[97]
[97] From function IOstream::check(const char*) const
[97] in file db/IOstreams/IOstreams/IOstream.C at line 96.
[97]
FOAM parallel run exiting
[97]
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 97 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.
--------------------------------------------------------------------------
[1]
[1]
[1] --> FOAM FATAL IO ERROR:
[1] error in IOstream "/home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor1/10/omega" for operation Ostream& operator<<(Ostream&, const Scalar&)
[1]
[1] file: /home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor1/10/omega at line 576284.
[1]
[1] From function IOstream::check(const char*) const
[1] in file db/IOstreams/IOstreams/IOstream.C at line 96.
[1]
FOAM parallel run exiting
[1]
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 34 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.
--------------------------------------------------------------------------
[57]
[57]
[57] --> FOAM FATAL IO ERROR:
[57] error in IOstream "/home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor57/10/omega" for operation Ostream& operator<<(Ostream&, const Scalar&)
[57]
[57] file: /home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor57/10/omega at line 661484.
[57]
[57] From function IOstream::check(const char*) const
[57] in file db/IOstreams/IOstreams/IOstream.C at line 96.
[57]
FOAM parallel run exiting
[57]
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 57 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.
--------------------------------------------------------------------------
[25]
[25]
[25] --> FOAM FATAL IO ERROR:
[25] error in IOstream "/home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor25/10/omega" for operation Ostream& operator<<(Ostream&, const Scalar&)
[25]
[25] file: /home/c20055/OpenFOAM/c20055-2.1.0/run/Automobile/AudiTT001/processor25/10/omega at line 666852.
[25]
[25] From function IOstream::check(const char*) const
[25] in file db/IOstreams/IOstreams/IOstream.C at line 96.
[25]
FOAM parallel run exiting
[25]
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 25 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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 70 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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 44 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.
--------------------------------------------------------------------------
[ERR.] PLE 0019 plexec One of MPI processes was aborted.(rank=1)(nid=0x02040030)(CODE=1878,793702259720650752,256)
Masashi Imano
Dec 20, 2012, 1:12:44 AM12/20/12
to open...@googlegroups.com
今野です。
東京大学情報基盤センターのどのスパコンを使われているのでしょうか?
HA8000 (T2K東大)でしょうか?
FX10 (Oakleaf-FX)でしょうか?
貴研究室サーバーで行った同様の解析(要素数16000000へ削減し8並列)は、
情報基盤センターで動かすと同じ結果が出ているのでしょうか?
2012/12/20 Konishi <ryo.kon...@gmail.com>:
> --
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> このディスカッションをウェブ上で閲覧するには、https://groups.google.com/d/msg/openfoam/-/PBB6l2HsBAwJ
> にアクセスしてください。
>
> このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
> このグループから退会するには、openfoam+u...@googlegroups.com にメールを送信してください。
> 詳細については、http://groups.google.com/group/openfoam?hl=ja からこのグループにアクセスしてください。
--
IMANO Masashi, Ph.D.
Konishi
Dec 20, 2012, 2:48:34 AM12/20/12
to open...@googlegroups.com
今野先生
お世話に成っております。konishiで御座います。
現在お借りしているサーバーはOakleaf-FXで御座います。要素数16000000の解析は、若干のズレはあるものの抵抗係数は似た数字をしてしておりました。これに設定した各種条件を約1億のメッシュへ適用し解析致しましたが、発散してしまうのが現状で御座います。
konishi
Masashi Imano
Dec 20, 2012, 4:30:43 AM12/20/12
to open...@googlegroups.com
今野です。
* Oakleaf-FXと研究室のPCで挙動が異なる件について
> 現在お借りしているサーバーはOakleaf-FXで御座います。要素数16000000の解析は、若干のズレはあるものの抵抗係数は似た数字をしてしておりました。これに設定した各種条件を約1億のメッシュへ適用し解析致しましたが、発散してしまうのが現状で御座います。
若干のズレがあるということは、プログラムソースが異なるか、
コンパイラの挙動が異なるということなりますが、研究室のソースも
FXと同じ2.1.0でしょうか？
もし研究室のほうが2.1.xであるなど、ソースが違うのであれば、
ソースが異なることが原因かもしれないので、特に乱流モデルのソースの
変更をgithubで調べることをお勧めします。
もし、ソースが同じなら、コンパイラの挙動の違いによると思いますので、
情報基盤センターに問いあわせてみてください。
* potentialFoamが動かない件について
functionがCd値などを計算するがcontrolDictにあると思うのですが、これが
potentialFoamではエラーになるかもしれないので、もしfunctionを無効にして
いなければ、以下のようにして無効してみてください。
potentialFoam -noFunctionObjects
* simpleFoamが発散する件について
離散化スキームやソルバーの設定がかなり複雑になっていますので、
基本はmoterBikeのチュートリアルの設定に戻して、fvSchemesで
Uのdivスキームだけupwindに変更して、発散しないか試してみることを
お勧めします。
divSchemes
{
default none;
div(phi,U) Gauss upwind;
また、ddtSchemesがクランクニコルソンになっておりましたが、
simpleFoamでは時間項が無いので意味がありませんのでsteadyStateにてください。
非定常で解きたいのであれば、pimpleFoam等の非定常用のソルバーを使ってください。
また、格子数が約1億ということで、物体表面の格子がかなり薄くなっていると思い
ますが、y+がかなり小さくなっており、kOmegaSSTのように高レイノルズ型乱流モデル
では解けない格子になっている可能性があるので、発散する前の計算結果を用いて
yPlusRASユーティリティ で y+ を確認してください。
以上です。
2012年12月20日 16:48 Konishi <ryo.kon...@gmail.com>:
> --
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> このディスカッションをウェブ上で閲覧するには、https://groups.google.com/d/msg/openfoam/-/RqnYPT-CU-QJ
S.kita
Dec 20, 2012, 10:39:12 PM12/20/12
to OpenFOAM
Konishi様
s.kitaといいます。
私も１ヶ月前にKonishi様同様に東大FX-10を用いてOpenFOAMの大規模計算を
流そうとしましたが、予備解析として取り組んだ100万セル程度の単相乱流計算でも
パフォーマンスが出ませんでした。並列効率ではなく、純粋に1stepあたりに要する
演算時間が長すぎるという意味です。
それは事前にシステム側にインストールされているモジュール(ver2.1.0)を用いてやったのですが、
今回Konishi様がやられているのはソースから最新バージョンの富士通コンパイラを用いて
buildしたものを使って、計算されようとしているのでしょうか？また何かチューニング等は
されていらっしゃいますか？
ソースからのbuildでしたら、是非予備解析等のパフォーマンスに関する情報を
簡単で結構ですので、教えて頂けると幸いです。
宜しくお願い致します。
On 12月20日, 午後6:30, Masashi Imano <masashi.im...@gmail.com> wrote:
> 今野です。
>
> * Oakleaf-FXと研究室のPCで挙動が異なる件について
>
> > 現在お借りしているサーバーはOakleaf-FXで御座います。要素数16000000の解析は、若干のズレはあるものの抵抗係数は似た数字をしてしておりました。これに設定した各種条件を約1億のメッシュへ適用し解析致しましたが、発散してしまうのが現状で御座います。
>
> 若干のズレがあるということは、プログラムソースが異なるか、
> コンパイラの挙動が異なるということなりますが、研究室のソースも
> FXと同じ2.1.0でしょうか？
>
> もし研究室のほうが2.1.xであるなど、ソースが違うのであれば、
> ソースが異なることが原因かもしれないので、特に乱流モデルのソースの
> 変更をgithubで調べることをお勧めします。
>
> kOmegaSST.Cの例https://github.com/OpenFOAM/OpenFOAM-2.1.x/commits/master/src/turbule...
>
> もし、ソースが同じなら、コンパイラの挙動の違いによると思いますので、
> 情報基盤センターに問いあわせてみてください。
>
> * potentialFoamが動かない件について
>
> functionがCd値などを計算するがcontrolDictにあると思うのですが、これが
> potentialFoamではエラーになるかもしれないので、もしfunctionを無効にして
> いなければ、以下のようにして無効してみてください。
>
> potentialFoam -noFunctionObjects
>
> * simpleFoamが発散する件について
>
> 離散化スキームやソルバーの設定がかなり複雑になっていますので、
> 基本はmoterBikeのチュートリアルの設定に戻して、fvSchemesで
> Uのdivスキームだけupwindに変更して、発散しないか試してみることを
> お勧めします。
>
> divSchemes
> {
> default none;
> div(phi,U) Gauss upwind;
>
> また、ddtSchemesがクランクニコルソンになっておりましたが、
> simpleFoamでは時間項が無いので意味がありませんのでsteadyStateにてください。
> 非定常で解きたいのであれば、pimpleFoam等の非定常用のソルバーを使ってください。
>
> また、格子数が約1億ということで、物体表面の格子がかなり薄くなっていると思い
> ますが、y+がかなり小さくなっており、kOmegaSSTのように高レイノルズ型乱流モデル
> では解けない格子になっている可能性があるので、発散する前の計算結果を用いて
> yPlusRASユーティリティ で y+ を確認してください。
>
> 以上です。
>
> 2012年12月20日 16:48 Konishi <ryo.konishi....@gmail.com>:
>
>
>
>
>
> > 今野先生
>
> > お世話に成っております。konishiで御座います。
>
> > 現在お借りしているサーバーはOakleaf-FXで御座います。要素数16000000の解析は、若干のズレはあるものの抵抗係数は似た数字をしてしておりました。これに設定した各種条件を約1億のメッシュへ適用し解析致しましたが、発散してしまうのが現状で御座います。
>
> > konishi
>
> > --
> > このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> > このディスカッションをウェブ上で閲覧するには、https://groups.google.com/d/msg/openfoam/-/RqnYPT-CU-QJ
> > にアクセスしてください。
>
> > このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
> > このグループから退会するには、openfoam+u...@googlegroups.com にメールを送信してください。
> > 詳細については、http://groups.google.com/group/openfoam?hl=jaからこのグループにアクセスしてください。
>
> --
> IMANO Masashi, Ph.D.- 引用テキストを表示しない -
>
> - 引用テキストを表示 -
Konishi
Dec 22, 2012, 12:49:54 AM12/22/12
to open...@googlegroups.com
今野先生
お世話に成っております。konishiで御座います。
先ずソースに関しましてお答え致します。今回の解析に際しまして、要素数約16000000の解析を行ったバージョンは2.0.1で御座います。しかし本研究室では2.1.0も使用可能で御座いますので、こちらでも解析を行ってみます。これに関連しまして、要素数約1億のメッシュを100分割する際、2.0.1では無く2.1.0で分割するとOakleaf-FXで解析を行った際に改善する可能性はあるでしょうか。
次に、御教示頂いたpotentialFoam及びsimpleFoamに関する情報を適用し再度解析を実行致しました。potentialFoamに関しましては先に提示致しましたログと同様のエラーが吐き出され、解析が停止してします。fvSchemeに関しましては、解析に使用していないnut等も「定義しろ」と吐き出された為、これらも全てupwindとして再度実行しました。しかし、この設定を適用しても解析が発散してしまいました。
また、kOmegaSST.Cは2.0.1と2.1.0で大きく異なって居りました。
そこで一点の質問が御座います。現在抵抗係数を得ようとしている車両は全体が少々複雑に成っており、これ以上格子点を削減する事が出来ません。yPulsRASを実行した所、「min: 0.051948 max: 754.75 average: 33.5026」となり、最小部分は非常に小さい値と成って居ります。そこで現状のメッシュでも解析可能な乱流モデルは御座いますでしょうか。
お忙しい中大変申し訳御座いませんが、御教示頂けると幸いで御座います。
konishi
Konishi
Dec 22, 2012, 12:55:17 AM12/22/12
to open...@googlegroups.com
s.kita様
お忙しい中有難う御座います。konishiで御座います。
システムに関する知識に疎い為、申し訳御座いませんが正しいお答えが出来ない可能性が有る事を御容赦頂けると幸いで御座います。
予備計算は所属研究室で所有しているサーバーで、2.0.1を使用しております。また、Oakelaf-FXでの解析は特にチューニング等はせずに解析を行っております。
予備解析(要素数約16000000)では収束までに約100ステップ(6時間少々)だったと記憶しております。
このご回答でよろしいでしょうか。
konishi
Masashi Imano
Dec 23, 2012, 6:25:37 PM12/23/12
to open...@googlegroups.com
今野です。
今回の件では、発散の原因が東大FX10のコンパイラにある可能性もあるので、
研究室でも東大FX10と同じ2.1.0を用いて計算結果がほとんど
一致することを確認して、発散の原因の切り分けすることをお勧め致します
2.0.1と2.1.0で領域分割( decomposePar )が変更されている可能性はありますが、
それによって並列化効率が変わっても、発散の原因になることはないと思います。
potentialFoamが動かない原因は良くわかりませんが、ログを見るとsimpleFoam
とは別の緩和係数になっており、異なるfvSolutionを使っていると思いますので、
potentialFoamとsimpleFoamで実行ディレクトリを変えているのでしょうか?
もしそうであれば、potentialFoamを動かしている場所の0/p,Uや
system/{fv*, controlDict} を貼ってください。
y+に関して最小値がかなり小さいようですが、平均値や最大値は十分大きいので
逆に低レイノルズ数型モデルは使えないようです。
もし可能であれば、格子生成の段階でyPlusが極端に小さい領域の格子の厚みを
再調整してみてください。
また、FX10で発散しないように試行している段階では、乱流モデルは
kEpsilonを使うほうが良いのではないでしょうか?
以上です。
2012年12月22日 14:49 Konishi <ryo.kon...@gmail.com>:
> --
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> このディスカッションをウェブ上で閲覧するには、https://groups.google.com/d/msg/openfoam/-/Q0Vin_JMFcAJ
S.kita
Dec 24, 2012, 6:15:26 AM12/24/12
to OpenFOAM
konishi様
s.kitaです。
ご回答有難うございます。
non tuningならば富士通C++コンパイラがまだ途上段階にあることからFX-10よりも
Redhat OSでOpteronCPUのT2Kの方が素直にパフォーマンスが出るはずなんですが、
T2Kは現在公開利用をしていないのが非常に残念なところです。
大学所属の方なら、科研費絡めてひょっとすると利用申請できるのかもしれません。
FX-10でのOpenFOAMのパフォーマンスがもし振るわなかったら、T2Kや他大の
スパコン利用を考えてみてはいかがでしょう？
s.kitaです。
ご回答有難うございます。
non tuningならば富士通C++コンパイラがまだ途上段階にあることからFX-10よりも
Redhat OSでOpteronCPUのT2Kの方が素直にパフォーマンスが出るはずなんですが、
T2Kは現在公開利用をしていないのが非常に残念なところです。
大学所属の方なら、科研費絡めてひょっとすると利用申請できるのかもしれません。
FX-10でのOpenFOAMのパフォーマンスがもし振るわなかったら、T2Kや他大の
スパコン利用を考えてみてはいかがでしょう？
Konishi
Dec 25, 2012, 12:23:31 AM12/25/12
to open...@googlegroups.com
今野先生
お忙しい中有難う御座います。konishiで御座います。
FX10と所属研究室での解析ですが、同じ2.1.0を用いても出力される抵抗係数には差がありました。
potentialFoamに関しましては、全て同一ディレクトリで解析を行っております。御指定のファイルを以下に提示致しますので御確認頂けると幸いで御座います。
また、メッシュの能解度ですが、形状を妥協してメッシュを粗くしましたが、やはり同様に発散してしまいます。更にkEpsilonですと10ステップ目前後からe+15と成ってしまいます。乱流エネルギーや散逸率は以前Ohbuchi様に御教示頂きましたURL (http://www.cfd-online.com/Wiki/Turbulence_free-stream_boundary_conditions) を参考に算出しております。
何度もお手数をお掛けしてしまい大変申し訳御座いません。御返事を頂けると幸いで御座います。
konishi
1. p
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
sidewall
{
type slip;
}
floor
{
type zeroGradient;
}
car
{
type zeroGradient;
}
ceiling
{
type slip;
}
}
// ************************************************************************* //
2. U
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 -80);
boundaryField
{
inlet
{
type fixedValue;
value uniform (0 0 -80);
}
outlet
{
type zeroGradient;
}
sidewall
{
type slip;
}
floor
{
type fixedValue;
value uniform (0 0 0);
}
car
{
type fixedValue;
value uniform (0 0 0);
}
ceiling
{
type slip;
}
}
// ************************************************************************* //
3. controlDict
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application potentialFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 15000;
deltaT 1;
writeControl timeStep;
writeInterval 100;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression off;
timeFormat general;
timePrecision 6;
runTimeModifiable true;
// functions
// {
// #include "readFields"
// #include "streamLines"
// #include "cuttingPlane"
// #include "forceCoeffs"
// }
// ************************************************************************* //
4. fvScheme
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}
divSchemes
{
default none;
div(phi,U) Gauss linearUpwindV grad(U);
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss upwind;
div(phi,nuTilda) Gauss upwind;
div((nuEff*dev(T(grad(U))))) Gauss linear;
div(phi,omega) Gauss upwind;
}
laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DREff,R) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
laplacian(DomegaEff,omega) Gauss linear corrected;
}
interpolationSchemes
{
default linear;
interpolate(U) linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p ;
}
// ************************************************************************* //
5. fvSolution
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver smoothSolver;
tolerance 1e-7;
relTol 0.001;
smoother DICGaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
cacheAgglomeration on;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 10;
mergeLevels 1;
};
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
omega
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1;
}
}
SIMPLE
{
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefPoint 0;
pRefValue 0;
}
potentialFlow
{
nNonOrthogonalCorrectors 10;
pRefCell 0;
pRefPoint 0;
pRefValue 10;
}
relaxationFactors
{
p 0.3;
U 0.7;
k 0.7;
epsilon 0.7;
omega 0.7;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Konishi
Dec 25, 2012, 12:28:25 AM12/25/12
to open...@googlegroups.com
s.kita様
konishiで御座います。
以前T2Kの利用申請を行って居た筈ですので、T2Kで試す事も可能で御座います。OpenFOAMをインストールし、試行致します。
目標達成の為にも、FX10の発展を願って止みません。
お忙しい中有難う御座いました。
konishi
Masashi Imano
Dec 25, 2012, 2:26:14 AM12/25/12
to open...@googlegroups.com
今野です。
貼って頂いた設定を用いて tutorials/mesh/snappyHexMesh/motorBike のチュートリアル
で potentialFoam を実行すると
potentialFoam -noFunctionObjects -writep
--> FOAM FATAL IO ERROR:
keyword laplacian(1,p) is undefined in dictionary
"/Users/imano/OpenFOAM/imano-2.1.x/run/tutorials/mesh/snappyHexMesh/motorBike/system/fvSchemes::laplacianSchemes"
とエラーが出ました。
これからわかる通り system/fvSchmes の laplacianSchemes ブロックに
以下が必要です。
system/fvSchmes:
laplacian(1,p) Gauss linear corrected;
上記を追加したら一応動きますが、pの線形ソルバーが smoothSolver となっており、
ディフォルトの反復上限回数の1000回になっても指定の収束判定条件を満たしませんし、
線形ソルバーの他のパラメータ設定がGAMG用ですから、線形ソルバーはGAMGに
します。
また、GAMGでの最粗レベル格子数のnCellsInCoarsestLevelが10になっていますが、
一億規模の格子で10まで再帰的に下げるとかなり多くのレベルが必要となり、計算コスト
がかかると思いますので、もう少し大きな値(例えば1000)にしたほうが良いと思います。
system/fvSolver:
solvers
{
p
{
solver GAMG;
tolerance 1e-7;
略
nCellsInCoarsestLevel 1000;
また、些細な話ですが、圧力参照点の格子番号を指定するpRefCell と格子位置を指定
する pRefPoint を同時に指定するのは意味が無く、ここでは最初に指定しているpRefCell
のみが有効になっています。また、元々pRefPointでは(1 1 1)といった位置を指定する必要
があるので、ここでの指定方法は間違いです。
さらに、圧力基準値を指定するpRefValueの値が10になっていますが、
この圧力の値は相対値であり絶対値には意味を持たないので、特に理由が無ければ0にします。
以上より、以下のように指定したほうが良いでしょう。
system/fvSolver:
potentialFlow
{
nNonOrthogonalCorrectors 10;
pRefCell 0;
pRefValue 0;
}
あと、一番重要な事をお聞きするのを忘れていたのですが、
checkMesh -constant でFATALなエラーは出てないのですよね?
2012年12月25日 14:23 Konishi <ryo.kon...@gmail.com>:
貼って頂いた設定を用いて tutorials/mesh/snappyHexMesh/motorBike のチュートリアル
で potentialFoam を実行すると
potentialFoam -noFunctionObjects -writep
--> FOAM FATAL IO ERROR:
"/Users/imano/OpenFOAM/imano-2.1.x/run/tutorials/mesh/snappyHexMesh/motorBike/system/fvSchemes::laplacianSchemes"
とエラーが出ました。
これからわかる通り system/fvSchmes の laplacianSchemes ブロックに
以下が必要です。
system/fvSchmes:
laplacian(1,p) Gauss linear corrected;
上記を追加したら一応動きますが、pの線形ソルバーが smoothSolver となっており、
ディフォルトの反復上限回数の1000回になっても指定の収束判定条件を満たしませんし、
線形ソルバーの他のパラメータ設定がGAMG用ですから、線形ソルバーはGAMGに
します。
また、GAMGでの最粗レベル格子数のnCellsInCoarsestLevelが10になっていますが、
一億規模の格子で10まで再帰的に下げるとかなり多くのレベルが必要となり、計算コスト
がかかると思いますので、もう少し大きな値(例えば1000)にしたほうが良いと思います。
system/fvSolver:
solvers
{
p
{
solver GAMG;
tolerance 1e-7;
略
nCellsInCoarsestLevel 1000;
また、些細な話ですが、圧力参照点の格子番号を指定するpRefCell と格子位置を指定
する pRefPoint を同時に指定するのは意味が無く、ここでは最初に指定しているpRefCell
のみが有効になっています。また、元々pRefPointでは(1 1 1)といった位置を指定する必要
があるので、ここでの指定方法は間違いです。
さらに、圧力基準値を指定するpRefValueの値が10になっていますが、
この圧力の値は相対値であり絶対値には意味を持たないので、特に理由が無ければ0にします。
以上より、以下のように指定したほうが良いでしょう。
system/fvSolver:
potentialFlow
{
nNonOrthogonalCorrectors 10;
pRefCell 0;
}
あと、一番重要な事をお聞きするのを忘れていたのですが、
checkMesh -constant でFATALなエラーは出てないのですよね?
2012年12月25日 14:23 Konishi <ryo.kon...@gmail.com>:
> --
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> このディスカッションをウェブ上で閲覧するには、https://groups.google.com/d/msg/openfoam/-/eEXeZyzeKLoJ
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
Konishi
Dec 25, 2012, 1:36:00 PM12/25/12
to open...@googlegroups.com
今野先生
御教示有難う御座います。konishiで御座います。
御指定頂きました通りに再設定しpotentialFoamを回した所、以下の様なログが吐き出されました。
[64] --> FOAM FATAL IO ERROR:
[64] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[64]
[64] file: IOstream::interpolationSchemes
[64]
[64] From function dictionary::lookupEntry(const word&, bool, bool) const
[64] in file db/dictionary/dictionary.C at line 396.
[64]
FOAM parallel run exiting
これを元にinterpolationSchemesのinterpolate(U)をコメントアウト、interpolationSchemes自体をコメントアウトの二通りを行った所、次は以下の様なログが吐き出されました。
Pleamoved by adjusting td by adjusting thby adjustinsting thjustingse checkhe oute outflowg the outfe outflow. the outfl the velocity flow.
Plea.
Pleaselow.
Plea
Please cow.
Plboundary se check the check the velocse check the vheck the velocityease check the vconditions and/or ruvelocity boundity boundary coelocity bou boundary conditelocity bounn potentialFoam ary conditions nditions and/or ions and/dary cto initialand/or run porun potentialFoaor run potentindary conditions and/or ronditions and/orise the outflow.tentialFoamm to initialalFoam to inun potentia run poten
Total flux to initialise the outitialislFoam ttialFo ise theflow.
Total e the outflow.
To initialise theam to initialise : 1e-300
Spe outflow.
Total flflux otal flux outflow.
Tota the outflow.
Tcified mass inflow ux : 1e-300
l flux otal flux : 1.64447e+0 : 1e-300
SpSpecified : 1e-30 : 6
Specified ecified masmass inflow0
Specifi : 1e-3mass outfs inflow : ed mas1e-300
Specifi00
Speciflow : 0
Adj : 1.644471.64447e+06
Ss inflow ed mass infied massustable mae+06
Specifpecified mas: 1.6444low : 1 inflow ss outfloied mass s outflow7e+06
Spe.64447e+0 : 1w : 0
outflow : 0
Adjucified ma6
Specifi.64447e+0[56] : 0
Adjustable mss outfed mass6
Specif
ログが散らかっておりますが、恐らく出口側へ流れて行っていないという旨のログであると考えております。各種境界条件を再度確認致しましたが、出口側へ流出しない設定には成って居りません。試行錯誤を重ねておりますが、未だ解決出来ておりません。
また、checkMeshに関しましては以下の様にログが吐き出され、不適切なメッシュは存在して居ない様です。お伝えし忘れておりましたが、今野先生からご教示頂いた通り、今回のメッシュは車体壁面のメッシュを粗くした物で、要素数約70000000、並列数は70と成って居ります。併せてご報告致します。
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create polyMesh for time = 0
Time = 0
Mesh stats
points: 11989325
faces: 140543075
internal faces: 138825937
cells: 69842253
boundary patches: 6
point zones: 0
face zones: 0
cell zones: 0
Overall number of cells of each type:
hexahedra: 0
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 69842253
polyhedra: 0
Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface topology
car 186710 93365 ok (closed singly connected)
ceiling 623896 313216 ok (non-closed singly connected)
floor 698588 350563 ok (non-closed singly connected)
inlet 21227 10916 ok (non-closed singly connected)
outlet 26427 13516 ok (non-closed singly connected)
sidewall 160290 82219 ok (non-closed singly connected)
Checking geometry...
Overall domain bounding box (-200 -1.55 -1000) (198.757 50 500)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-6.64509e-16 -1.77658e-14 -1.63062e-16) OK.
Max cell openness = 2.40386e-15 OK.
Max aspect ratio = 122.076 OK.
Minumum face area = 1.03533e-06. Maximum face area = 4.91035. Face area magnitudes OK.
Min volume = 8.09346e-10. Max volume = 3.02402. Total volume = 3.08338e+07. Cell volumes OK.
Mesh non-orthogonality Max: 88.9434 average: 20.7102
*Number of severely non-orthogonal faces: 517.
Non-orthogonality check OK.
<<Writing 517 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 1.06114 OK.
Coupled point location match (average 0) OK.
Mesh OK.
お忙しい中いつも有難う御座います。宜しくお願い申し上げます。
konishi
Masashi Imano
Dec 25, 2012, 9:23:16 PM12/25/12
to open...@googlegroups.com
今野です。
> これを元にinterpolationSchemesのinterpolate(U)をコメントアウト、interpolationSchemes自体をコメントアウトの二通りを行った所、次は以下の様なログが吐き出されました。
このエラーの原因は良くわかりませんが、とりあえずpotentialFoam や simpleFoam が
FX10で動くまでは、前にも書きました通り、motorBikeチュートリアルの fvSchemes や
fvSolution をコピーして使ったほうが良いと思います。
さらに、fvSchemes において div(phi,U) は upwind にしたほうが安全でしょう。
divSchemes
{
default none;
div(phi,U) Gauss linearUpwindV grad(U);
> ログが散らかっておりますが、恐らく出口側へ流れて行っていないという旨のログであると考えております。各種境界条件を再度確認致しましたが、出口側へ流出しない設定には成って居りません。試行錯誤を重ねておりますが、未だ解決出来ておりません。
このように重複したメッセージがログに出るというのは、OpenFOAMのアプリケーションを
-parallel オプションを付けないので並列実行した時に現われる現象ですが、potentialFoam を
並列実行する際に -parallel オプション付けてますか?
2012/12/26 Konishi <ryo.kon...@gmail.com>:
> --
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> このディスカッションをウェブ上で閲覧するには、https://groups.google.com/d/msg/openfoam/-/Qb94Za-DA-cJ
Masashi Imano
Dec 25, 2012, 9:40:07 PM12/25/12
to open...@googlegroups.com
今野です。
> FX10と所属研究室での解析ですが、同じ2.1.0を用いても出力される抵抗係数には差がありました。
何%程度差があるのでしょうか?
最終桁が付近が異なる程度でしたら良いのですが、それ以上違うとなると
FX10でのOpenFOAMの不具合かもしれないので、早急に情報基盤センターに問い
あわせたがほう良いと思います。
あと、motorBike のチュートリアルについても、研究室とFX10で結果が同じになるか
チェックしたほうが良いと思います。
2012年12月26日 11:23 Masashi Imano <masash...@gmail.com>:
--
IMANO Masashi, Ph.D.
Masashi Imano
Dec 25, 2012, 10:05:40 PM12/25/12
to open...@googlegroups.com
今野です。
> あと、motorBike のチュートリアルについても、研究室とFX10で結果が同じになるか
> チェックしたほうが良いと思います。
これらのチェックですが、単一ノード並列と複数ノード並列とでチェックして
みてください。
ただ、ログに出力される forceCoeffs の値が異なるか確認できれば良いので、
FX10のトークンを無駄に使わないよう、simpleFoamは 数ステップ動かすだけ
で良いかもしれません。
forceCoeffs output:
Cd = 2.25098
Cl = 0.0207447
Cm = 0.885132
2012年12月26日 11:40 Masashi Imano <masash...@gmail.com>:
--
IMANO Masashi, Ph.D.
Konishi
Dec 26, 2012, 2:18:53 AM12/26/12
to open...@googlegroups.com
今野先生
お世話に成っております。konishiで御座います。
先ず、fvSchemesのdiv(phi.U)ですが、upwindにして再度解析を行います。またpotentialFoamですが、これは-parallelオプションを付けて実行しております。
続きましてFX10と所属研究室との解析の差に関しまして、正確な値を得る為にmotorBikeを含め幾つかのケースを回しご報告致します。
宜しくお願い申し上げます。
konishi
Konishi
Jan 1, 2013, 5:31:01 AM1/1/13
to open...@googlegroups.com
Google Groups 「OpenFOAM」の皆様
konishiで御座います。旧年は一方ならぬ御厚情を賜り、誠に有難う御座います。本年も変わらず宜しくお願い申し上げます。皆様の御健康と御多幸をお祈り申し上げます。
さて、今野先生より御推奨頂きました「ソルバーの挙動検証」に関しまして、データを揃える事が出来ましたので添付ファイルにて御報告致します。解析ケースはチュートリアルの「motorBike」と、以前解析を行っていた船舶で御座います。結論から申し上げますと、Oakleaf-FXと所属研究室との間に差は見られませんでした。
さて、今野先生より御推奨頂きました「ソルバーの挙動検証」に関しまして、データを揃える事が出来ましたので添付ファイルにて御報告致します。解析ケースはチュートリアルの「motorBike」と、以前解析を行っていた船舶で御座います。結論から申し上げますと、Oakleaf-FXと所属研究室との間に差は見られませんでした。
そこで、今回の問題の原因は乱流エネルギーや比散逸率の計算及び設定に問題があるのかと考えております。解析領域を、解析対象の車体に対して20倍以上の大きさで設定しておりますが、乱流エネルギー等の設定値は僅かな差(例えば0.5か0.6か)によって要素数約1億では発散してしまうものなのでしょうか。
新年早々申し訳ございませんが、御手隙に成りましたら御教示頂けると幸いで御座います。
konishi
Masashi Imano
Jan 1, 2013, 8:57:00 PM1/1/13
to open...@googlegroups.com
あけましておめでとうございます。
今野です。
> さて、今野先生より御推奨頂きました「ソルバーの挙動検証」に関しまして、データを揃える事が出来ましたので添付ファイルにて御報告致します。解析ケースはチュートリアルの「motorBike」と、以前解析を行っていた船舶で御座います。結論から申し上げますと、Oakleaf-FXと所属研究室との間に差は見られませんでした。
確かにmotorBikeや船舶のsingleはCdの解析結果に差が無いのですが、
添付の通り、船舶のparallelは反復過程を平均しても80%もの差があります。
これの差異が流れ場の解析結果自体の差によるものなのか、
それともCdを計算するfunctionがFXで正しく動いていないのか確かめたいので、
とりあえず、研究室とFXでの船舶parallel計算のログを全て添付してください。
2013年1月1日 19:31 Konishi <ryo.kon...@gmail.com>:
> --
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> このディスカッションをウェブ上で閲覧するには、https://groups.google.com/d/msg/openfoam/-/5t-HtjcNZJ8J
Masashi Imano
Jan 2, 2013, 5:12:57 AM1/2/13
to open...@googlegroups.com
今野です。
> 確かに、先生のファイルを拝見致しますと不自然に感じます。御指定のファイルを添付いたしますので御確認頂けると幸いで御座います。
今回Word形式で送って頂きましたが、このようなテキストファイルの場合、
テキスト形式のほうが比較検討がしやすいで、テキスト形式にしたものを添付します。
また、FXのログは以下のファイルでデバッグスイッチが1になっているものがある
ために、デバッグログが大量に含まれているので、これらは削除してあります。
/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/etc/controlDict:
DebugSwitches
{
fvSchemes 1;
solution 1;
}
二つのlogを比較すると、解析結果も異なるのですが、以下の緩和係数の
出力がFXのほうにしかありません。
relaxation factors:
fields:
{
p 0.3;
}
equations:
{
nuTilda 0.7;
}
これはLaboのOpenFOAMが2.1.0より古いためだと思いますが、以下のようにログの
先頭部分のBuild行を確認してください。
Build : 2.1.0-bd7367f93311
以上です。
2013年1月2日 17:33 Konishi <ryo.kon...@gmail.com>:
> 確かに、先生のファイルを拝見致しますと不自然に感じます。御指定のファイルを添付いたしますので御確認頂けると幸いで御座います。
今回Word形式で送って頂きましたが、このようなテキストファイルの場合、
テキスト形式のほうが比較検討がしやすいで、テキスト形式にしたものを添付します。
また、FXのログは以下のファイルでデバッグスイッチが1になっているものがある
ために、デバッグログが大量に含まれているので、これらは削除してあります。
/usr/local/OpenFOAM/2.1.0/OpenFOAM-2.1.0/etc/controlDict:
DebugSwitches
{
fvSchemes 1;
solution 1;
}
二つのlogを比較すると、解析結果も異なるのですが、以下の緩和係数の
出力がFXのほうにしかありません。
relaxation factors:
fields:
{
p 0.3;
}
equations:
{
p 0.3;
U 0.7;
k 0.7;
epsilon 0.7;
R 0.7;
U 0.7;
k 0.7;
epsilon 0.7;
nuTilda 0.7;
}
これはLaboのOpenFOAMが2.1.0より古いためだと思いますが、以下のようにログの
先頭部分のBuild行を確認してください。
Build : 2.1.0-bd7367f93311
以上です。
2013年1月2日 17:33 Konishi <ryo.kon...@gmail.com>:
> --
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> このディスカッションをウェブ上で閲覧するには、https://groups.google.com/d/msg/openfoam/-/93No0KFFJZEJ
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
Konishi
Jan 2, 2013, 10:05:49 AM1/2/13
to open...@googlegroups.com
今野先生
お世話に成って居ります。konishiで御座います。
二つのファイルを添付した際、容量オーバーの為かファイルを含んだ投稿が弾かれてしまいました。以後、御教示頂いた処置を行いテキストファイルにて添付致します。有難う御座います。
さて、先生から御指摘頂いた「所属研究室側のBuild行の確認」ですが、以下の様に成って居る事をご報告致します。
2.1.0-0bc225064152
また蛇足ではありますが、Oakleaf-FXで解析を実行する際、初回実行時はinterpolateSchemesのinterpolate(U)や(HbyA)を設定せず解析出来たにも関わらず二回目以降は両者の設定を求められ、設定した後も同様のエラーで解析が停止してしまう事が多々有ります。残存ストレージ容量を考えても余裕が有る為、この点に関しましては非常に不可解に感じております。
[21] keyword interpolate(HbyA) is undefined in dictionary "IOstream::interpolationSchemes"
[21]
[21] file: IOstream::interpolationSchemes
[21]
[21] From function dictionary::lookupEntry(const word&, bool, bool) const
[21] in file db/dictionary/dictionary.C at line 396.
[21]
FOAM parallel run exiting
konishi
Masashi Imano
Jan 2, 2013, 9:35:55 PM1/2/13
to open...@googlegroups.com
今野です。
> さて、先生から御指摘頂いた「所属研究室側のBuild行の確認」ですが、以下の様に成って居る事をご報告致します。
>
> 2.1.0-0bc225064152
では研究室でもやはり同じ2.1.0を使っているのですね。
こちらでも、緩和係数の値のデバッグログは2.1.0でも通常出力しないことを確認しました。
こちらの勘違いだったようです。すみません。
> また蛇足ではありますが、Oakleaf-FXで解析を実行する際、初回実行時はinterpolateSchemesのinterpolate(U)や(HbyA)を設定せず解析出来たにも関わらず二回目以降は両者の設定を求められ、設定した後も同様のエラーで解析が停止してしまう事が多々有ります。残存ストレージ容量を考えても余裕が有る為、この点に関しましては非常に不可解に感じております。
確かに不可解ですね。
バッチジョブスクリプトが原因とも考えられるので、良かったらバッチジョブ全体
を貼ってください。
2013年1月3日 0:05 Konishi <ryo.kon...@gmail.com>:
> --
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> このディスカッションをウェブ上で閲覧するには、https://groups.google.com/d/msg/openfoam/-/G2cucwiLcMgJ
Konishi
Jan 2, 2013, 11:30:14 PM1/2/13
to open...@googlegroups.com
今野先生
お世話に成って居ります。konishiで御座います。
「バッチジョブ」とは以下の2点の事でよろしいでしょうか。
1. simpleFoam
#!/bin/sh
#PJM -L "rscgrp=regular"
#PJM -j
#PJM -L "node=12"
#PJM --mpi "proc=70"
#PJM -L "elapse=48:00:00"
cd $PJM_JOBDIR
mpirun -np 70 simpleFoam -parallel
2. potentialFoam
#!/bin/sh
#PJM -L "rscgrp=regular"
#PJM -j
#PJM -L "node=12"
#PJM --mpi "proc=70"
#PJM -L "elapse=48:00:00"
cd $PJM_JOBDIR
mpirun -np 70 potentialFoam -noFunctionObjects -parallel
konishi
Masashi Imano
Jan 3, 2013, 12:40:04 AM1/3/13
to open...@googlegroups.com
今野です。
> 「バッチジョブ」とは以下の2点の事でよろしいでしょうか。
はいそうです。
> 1. simpleFoam
>
> #!/bin/sh
>
> #PJM -L "rscgrp=regular"
> #PJM -j
>
> #PJM -L "node=12"
> #PJM --mpi "proc=70"
> #PJM -L "elapse=48:00:00"
> cd $PJM_JOBDIR
> mpirun -np 70 simpleFoam -parallel
OakleafのFX-10は1ノード12コアなので、12ノード使って70コア使用というのは
間違っていないかもしれませんが、勿体ないです^^;
16の倍数のコアを使う設定にはできないでしょうか?
また、--mpiオプションでコア数を指定しているので
mpirunには-npオプションが不要だったと思います。
さらに、以下のOpenFOAM設定行を入れるほうが良いと思います。
module load OpenFOAM
source $WM_PROJECT_DIR/etc/bashrc
simpleFoam用バッチジョブスクリプト (ノード数やコア数は適当です)
---
#PJM --mpi "proc=80"
source $WM_PROJECT_DIR/etc/bashrc
cd $PJM_JOBDIR
mpirun simpleFoam -parallel
---
只今電車内でOakleafにログインできないので試せないので、間違っているかもしれませんが、
とりあえず以上試してみてください。
また、上記のpotentialFoamが動かない件ですが、流出側の圧力がディレクレ境界条件でなく、
勾配ゼロの場合、元々potentialFoamが動かないようです。
流出側の圧力の境界条件をmotorBikeのようにディレクレ型(fixedValue)にして、
fvSolutionのpRef*をコメントアウトすればpotentialFoamが動くと思いますが、
現状の境界条件で粗い格子の計算結果が既にあるので、そこまでする必要はないかも
しれません。
また、本来はpotentialFoamを用いるよりも粗い格子の結果をmapFieldsを用いて
マッピングするほうが収束は早いと思います。
とりあえず以上です。
2013年1月3日 13:30 Konishi <ryo.kon...@gmail.com>:
> 「バッチジョブ」とは以下の2点の事でよろしいでしょうか。
はいそうです。
> 1. simpleFoam
>
> #!/bin/sh
>
> #PJM -L "rscgrp=regular"
> #PJM -j
>
> #PJM -L "node=12"
> #PJM --mpi "proc=70"
> #PJM -L "elapse=48:00:00"
> cd $PJM_JOBDIR
> mpirun -np 70 simpleFoam -parallel
間違っていないかもしれませんが、勿体ないです^^;
16の倍数のコアを使う設定にはできないでしょうか?
また、--mpiオプションでコア数を指定しているので
mpirunには-npオプションが不要だったと思います。
さらに、以下のOpenFOAM設定行を入れるほうが良いと思います。
module load OpenFOAM
source $WM_PROJECT_DIR/etc/bashrc
simpleFoam用バッチジョブスクリプト (ノード数やコア数は適当です)
---
#!/bin/sh
#PJM -L "rscgrp=regular"
#PJM -j
#PJM -L "node=5"
#PJM -L "rscgrp=regular"
#PJM -j
#PJM --mpi "proc=80"
#PJM -L "elapse=48:00:00"
module load OpenFOAM
source $WM_PROJECT_DIR/etc/bashrc
cd $PJM_JOBDIR
mpirun simpleFoam -parallel
---
只今電車内でOakleafにログインできないので試せないので、間違っているかもしれませんが、
とりあえず以上試してみてください。
勾配ゼロの場合、元々potentialFoamが動かないようです。
流出側の圧力の境界条件をmotorBikeのようにディレクレ型(fixedValue)にして、
fvSolutionのpRef*をコメントアウトすればpotentialFoamが動くと思いますが、
現状の境界条件で粗い格子の計算結果が既にあるので、そこまでする必要はないかも
しれません。
また、本来はpotentialFoamを用いるよりも粗い格子の結果をmapFieldsを用いて
マッピングするほうが収束は早いと思います。
とりあえず以上です。
2013年1月3日 13:30 Konishi <ryo.kon...@gmail.com>:
> --
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> このディスカッションをウェブ上で閲覧するには、https://groups.google.com/d/msg/openfoam/-/CS1AJDCH4-4J
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
Konishi
Jan 4, 2013, 12:22:10 AM1/4/13
to open...@googlegroups.com
今野先生
お世話に成って居ります。konishiで御座います。解析の様子見及び投稿時のエラーに見舞われ、ご報告が滞ってしまい申し訳御座いません。
御教示頂きました通りに[module load OpenFOAM source $WM_PROJECT_DIR/etc/bashrc]を追加し、fvSchemesとfvSolutionを以下の様に設定しました所、potentialFoam及びsimpleFoamが回り、今迄で最も良い抵抗係数を得られております。下記1, 2の設定では解析を実行出来ましたが、3のfvSolutionを用いると必ずinterpolateSchemesのエラーが表示され解析が停止する所迄分かりました事をご報告致します。また、4は抵抗係数の推移で御座います。添付ファイルとしての投稿を試みましたが、投稿時にエラーが発生する為この様な形式でのご報告とさせて頂きます。併せて御確認頂けると幸いで御座います。
1. fvSchemes
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,omega) Gauss upwind;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p;
}
// ************************************************************************* //
2. fvSolution
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver GAMG;
tolerance 1e-7;
relTol 0.01;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
cacheAgglomeration on;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 10;
mergeLevels 1;
};
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
omega
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
epsilon
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
R
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
nuTilda
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
}
SIMPLE
{
nNonOrthogonalCorrectors 0;
// pRefCell 0;
// pRefPoint 0;
// pRefValue 0;
}
potentialFlow
{
nNonOrthogonalCorrectors 0;
// pRefCell 0;
// pRefPoint 0;
// pRefValue 0;
}
relaxationFactors
{
p 0.05;
U 0.15;
k 0.15;
epsilon 0.15;
omega 0.15;
}
// ************************************************************************* //
3. fvSolution
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver GAMG;
tolerance 1e-7;
relTol 0.1;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
cacheAgglomeration on;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 1000;
mergeLevels 1;
}
U
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1;
}
k
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1;
}
omega
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1;
}
}
SIMPLE
{
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefPoint 0;
pRefValue 0;
}
potentialFlow
{
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefPoint 0;
pRefValue 0;
}
relaxationFactors
{
fields
{
p 0.05;
}
equations
{
U 0.15;
k 0.15;
omega 0.15;
}
}
cache
{
grad(U);
}
// ************************************************************************* //
4. Cd number
| 0.1 | -0.49417 |
| 0.2 | -1.36719 |
| 0.3 | -2.51861 |
| 0.4 | -3.38012 |
| 0.5 | -3.48472 |
| 0.6 | -3.50734 |
| 0.7 | -3.40804 |
| 0.8 | -3.20634 |
| 0.9 | -2.94749 |
| 1 | -2.68481 |
| 1.1 | -1.57977 |
| 1.2 | -0.48696 |
| 1.3 | 0.573262 |
| 1.4 | 1.62382 |
| 1.5 | 2.6425 |
| 1.6 | 3.33302 |
| 1.7 | 3.57137 |
| 1.8 | 3.57731 |
| 1.9 | 3.45156 |
| 2 | 3.20908 |
| 2.1 | 2.83924 |
| 2.2 | 2.33181 |
| 2.3 | 1.69263 |
| 2.4 | 0.954182 |
| 2.5 | 0.181661 |
| 2.6 | -0.10942 |
| 2.7 | -0.48431 |
| 2.8 | -0.88575 |
| 2.9 | -1.21564 |
| 3 | -1.40365 |
| 3.1 | -1.48884 |
| 3.2 | -1.44131 |
| 3.3 | -1.16254 |
| 3.4 | -0.77899 |
| 3.5 | -0.4688 |
| 3.6 | 0.030908 |
| 3.7 | 0.352981 |
| 3.8 | 0.797995 |
| 3.9 | 1.12393 |
| 4 | 1.29328 |
| 4.1 | 1.34009 |
| 4.2 | 1.32255 |
| 4.3 | 1.26436 |
| 4.4 | 1.15738 |
| 4.5 | 0.986297 |
| 4.6 | 0.759699 |
| 4.7 | 0.683803 |
| 4.8 | 0.632553 |
| 4.9 | 0.520366 |
| 5 | 0.375769 |
| 5.1 | 0.243973 |
| 5.2 | 0.138067 |
| 5.3 | 0.065831 |
| 5.4 | 0.02898 |
| 5.5 | 0.024243 |
| 5.6 | 0.046924 |
| 5.7 | 0.092147 |
| 5.8 | 0.155019 |
| 5.9 | 0.229962 |
| 6 | 0.31113 |
| 6.1 | 0.391968 |
| 6.2 | 0.465285 |
| 6.3 | 0.523821 |
| 6.4 | 0.563936 |
| 6.5 | 0.587894 |
| 6.6 | 0.602559 |
| 6.7 | 0.589994 |
| 6.8 | 0.56059 |
| 6.9 | 0.519297 |
| 7 | 0.498308 |
| 7.1 | 0.48028 |
| 7.2 | 0.485017 |
| 7.3 | 0.492872 |
| 7.4 | 0.493066 |
| 7.5 | 0.471246 |
| 7.6 | 0.437223 |
| 7.7 | 0.402424 |
| 7.8 | 0.371837 |
| 7.9 | 0.347814 |
| 8 | 0.331227 |
| 8.1 | 0.32189 |
| 8.2 | 0.325757 |
| 8.3 | 0.336703 |
| 8.4 | 0.350507 |
| 8.5 | 0.364169 |
| 8.6 | 0.366615 |
| 8.7 | 0.369932 |
| 8.8 | 0.374709 |
| 8.9 | 0.380306 |
| 9 | 0.385537 |
| 9.1 | 0.389889 |
| 9.2 | 0.393436 |
| 9.3 | 0.396412 |
| 9.4 | 0.399437 |
| 9.5 | 0.403375 |
| 9.6 | 0.4092 |
| 9.7 | 0.416918 |
| 9.8 | 0.425172 |
| 9.9 | 0.423002 |
| 10 | 0.421749 |
| 10.1 | 0.421499 |
| 10.2 | 0.422315 |
| 10.3 | 0.420959 |
| 10.4 | 0.417043 |
| 10.5 | 0.410501 |
| 10.6 | 0.403719 |
| 10.7 | 0.401315 |
| 10.8 | 0.401305 |
| 10.9 | 0.401713 |
| 11 | 0.401741 |
| 11.1 | 0.401143 |
| 11.2 | 0.399957 |
| 11.3 | 0.398351 |
| 11.4 | 0.396548 |
| 11.5 | 0.394071 |
| 11.6 | 0.390856 |
| 11.7 | 0.387314 |
| 11.8 | 0.383993 |
| 11.9 | 0.381502 |
| 12 | 0.380263 |
| 12.1 | 0.380466 |
| 12.2 | 0.382242 |
| 12.3 | 0.389395 |
| 12.4 | 0.396591 |
| 12.5 | 0.400941 |
| 12.6 | 0.405735 |
| 12.7 | 0.410108 |
| 12.8 | 0.413801 |
| 12.9 | 0.416326 |
| 13 | 0.417394 |
| 13.1 | 0.416863 |
| 13.2 | 0.414664 |
| 13.3 | 0.410881 |
| 13.4 | 0.40584 |
| 13.5 | 0.400014 |
| 13.6 | 0.393949 |
| 13.7 | 0.388143 |
| 13.8 | 0.383056 |
| 13.9 | 0.37905 |
| 14 | 0.376329 |
| 14.1 | 0.374892 |
| 14.2 | 0.374567 |
| 14.3 | 0.375128 |
| 14.4 | 0.376361 |
| 14.5 | 0.378085 |
| 14.6 | 0.380142 |
| 14.7 | 0.382381 |
| 14.8 | 0.384649 |
| 14.9 | 0.386804 |
| 15 | 0.388703 |
| 15.1 | 0.390248 |
| 15.2 | 0.391377 |
| 15.3 | 0.392015 |
| 15.4 | 0.392096 |
| 15.5 | 0.391579 |
| 15.6 | 0.390481 |
| 15.7 | 0.388922 |
| 15.8 | 0.387135 |
| 15.9 | 0.385343 |
| 16 | 0.383724 |
| 16.1 | 0.382386 |
| 16.2 | 0.381366 |
| 16.3 | 0.380637 |
| 16.4 | 0.380121 |
| 16.5 | 0.379716 |
| 16.6 | 0.379327 |
| 16.7 | 0.378893 |
| 16.8 | 0.378382 |
| 16.9 | 0.377811 |
| 17 | 0.377238 |
| 17.1 | 0.376679 |
| 17.2 | 0.37618 |
| 17.3 | 0.375772 |
| 17.4 | 0.375492 |
| 17.5 | 0.375349 |
| 17.6 | 0.375311 |
| 17.7 | 0.375319 |
| 17.8 | 0.375307 |
| 17.9 | 0.375221 |
| 18 | 0.37503 |
| 18.1 | 0.37473 |
| 18.2 | 0.374341 |
| 18.3 | 0.373896 |
| 18.4 | 0.373434 |
| 18.5 | 0.372991 |
| 18.6 | 0.372595 |
| 18.7 | 0.372255 |
| 18.8 | 0.371966 |
| 18.9 | 0.371709 |
| 19 | 0.371454 |
| 19.1 | 0.371174 |
| 19.2 | 0.370829 |
| 19.3 | 0.370387 |
| 19.4 | 0.369847 |
| 19.5 | 0.369245 |
| 19.6 | 0.368638 |
| 19.7 | 0.368088 |
| 19.8 | 0.367644 |
| 19.9 | 0.367334 |
| 20 | 0.367161 |
| 20.1 | 0.367108 |
| 20.2 | 0.367135 |
| 20.3 | 0.367196 |
| 20.4 | 0.36724 |
| 20.5 | 0.367229 |
| 20.6 | 0.367144 |
| 20.7 | 0.366992 |
| 20.8 | 0.366794 |
| 20.9 | 0.366552 |
| 21 | 0.366262 |
| 21.1 | 0.365914 |
| 21.2 | 0.365513 |
| 21.3 | 0.365079 |
| 21.4 | 0.364637 |
| 21.5 | 0.364206 |
| 21.6 | 0.363796 |
| 21.7 | 0.363411 |
| 21.8 | 0.36305 |
| 21.9 | 0.362707 |
| 22 | 0.362376 |
| 22.1 | 0.362054 |
| 22.2 | 0.361744 |
| 22.3 | 0.361467 |
| 22.4 | 0.361227 |
| 22.5 | 0.361031 |
| 22.6 | 0.360879 |
| 22.7 | 0.360758 |
| 22.8 | 0.360653 |
| 22.9 | 0.36055 |
| 23 | 0.360446 |
| 23.1 | 0.360327 |
| 23.2 | 0.360187 |
| 23.3 | 0.360018 |
| 23.4 | 0.359818 |
| 23.5 | 0.359593 |
| 23.6 | 0.359354 |
| 23.7 | 0.359116 |
| 23.8 | 0.358894 |
| 23.9 | 0.358703 |
| 24 | 0.358555 |
| 24.1 | 0.358454 |
| 24.2 | 0.358397 |
| 24.3 | 0.358374 |
| 24.4 | 0.358367 |
| 24.5 | 0.358356 |
| 24.6 | 0.358324 |
| 24.7 | 0.35826 |
| 24.8 | 0.358163 |
| 24.9 | 0.358038 |
| 25 | 0.357888 |
| 25.1 | 0.357718 |
| 25.2 | 0.357523 |
| 25.3 | 0.357315 |
| 25.4 | 0.3571 |
| 25.5 | 0.356889 |
| 25.6 | 0.356692 |
| 25.7 | 0.356516 |
| 25.8 | 0.356365 |
| 25.9 | 0.356238 |
| 26 | 0.356129 |
| 26.1 | 0.356032 |
| 26.2 | 0.355937 |
| 26.3 | 0.355835 |
| 26.4 | 0.355722 |
| 26.5 | 0.355593 |
| 26.6 | 0.355448 |
| 26.7 | 0.355298 |
| 26.8 | 0.355147 |
| 26.9 | 0.355001 |
| 27 | 0.354866 |
| 27.1 | 0.354742 |
| 27.2 | 0.354636 |
| 27.3 | 0.354547 |
| 27.4 | 0.354475 |
| 27.5 | 0.354421 |
| 27.6 | 0.354379 |
| 27.7 | 0.354355 |
| 27.8 | 0.354338 |
| 27.9 | 0.354324 |
| 28 | 0.354307 |
| 28.1 | 0.354287 |
| 28.2 | 0.35426 |
| 28.3 | 0.354226 |
| 28.4 | 0.354186 |
| 28.5 | 0.354141 |
| 28.6 | 0.354092 |
| 28.7 | 0.354042 |
| 28.8 | 0.353997 |
| 28.9 | 0.353954 |
| 29 | 0.353939 |
| 29.1 | 0.353938 |
| 29.2 | 0.353933 |
| 29.3 | 0.353923 |
| 29.4 | 0.353907 |
| 29.5 | 0.353884 |
| 29.6 | 0.353856 |
| 29.7 | 0.353821 |
| 29.8 | 0.353778 |
| 29.9 | 0.353725 |
| 30 | 0.353662 |
| 30.1 | 0.353586 |
| 30.2 | 0.353496 |
| 30.3 | 0.35339 |
| 30.4 | 0.353274 |
| 30.5 | 0.353158 |
| 30.6 | 0.353042 |
| 30.7 | 0.352926 |
| 30.8 | 0.352813 |
| 30.9 | 0.352701 |
| 31 | 0.352594 |
| 31.1 | 0.352496 |
| 31.2 | 0.352409 |
| 31.3 | 0.352335 |
| 31.4 | 0.352276 |
| 31.5 | 0.352228 |
| 31.6 | 0.352187 |
| 31.7 | 0.352147 |
| 31.8 | 0.352106 |
| 31.9 | 0.35206 |
| 32 | 0.352009 |
| 32.1 | 0.351954 |
| 32.2 | 0.351899 |
| 32.3 | 0.351845 |
| 32.4 | 0.351797 |
| 32.5 | 0.351755 |
| 32.6 | 0.35172 |
| 32.7 | 0.35169 |
| 32.8 | 0.351664 |
| 32.9 | 0.35164 |
| 33 | 0.351618 |
| 33.1 | 0.351598 |
| 33.2 | 0.35158 |
| 33.3 | 0.351567 |
| 33.4 | 0.351562 |
| 33.5 | 0.351565 |
| 33.6 | 0.35158 |
| 33.7 | 0.351607 |
| 33.8 | 0.351647 |
| 33.9 | 0.351696 |
| 34 | 0.351753 |
| 34.1 | 0.351812 |
| 34.2 | 0.351867 |
| 34.3 | 0.351914 |
| 34.4 | 0.351947 |
| 34.5 | 0.351965 |
| 34.6 | 0.351965 |
| 34.7 | 0.351947 |
| 34.8 | 0.351913 |
| 34.9 | 0.351865 |
| 35 | 0.351805 |
| 35.1 | 0.351735 |
| 35.2 | 0.351658 |
| 35.3 | 0.351577 |
| 35.4 | 0.351492 |
| 35.5 | 0.351406 |
| 35.6 | 0.351318 |
| 35.7 | 0.351231 |
| 35.8 | 0.351145 |
| 35.9 | 0.351063 |
| 36 | 0.350987 |
| 36.1 | 0.350919 |
| 36.2 | 0.350861 |
| 36.3 | 0.350814 |
| 36.4 | 0.350777 |
| 36.5 | 0.350748 |
| 36.6 | 0.350733 |
| 36.7 | 0.350725 |
| 36.8 | 0.350721 |
| 36.9 | 0.35072 |
| 37 | 0.350717 |
| 37.1 | 0.35071 |
| 37.2 | 0.350698 |
| 37.3 | 0.350679 |
| 37.4 | 0.350652 |
| 37.5 | 0.350618 |
| 37.6 | 0.350577 |
| 37.7 | 0.350532 |
| 37.8 | 0.350483 |
| 37.9 | 0.350432 |
| 38 | 0.350381 |
| 38.1 | 0.350332 |
| 38.2 | 0.350286 |
| 38.3 | 0.350246 |
| 38.4 | 0.35021 |
| 38.5 | 0.350182 |
| 38.6 | 0.350163 |
| 38.7 | 0.35015 |
| 38.8 | 0.350144 |
| 38.9 | 0.350143 |
| 39 | 0.350146 |
| 39.1 | 0.350152 |
| 39.2 | 0.350158 |
| 39.3 | 0.350163 |
| 39.4 | 0.350166 |
| 39.5 | 0.350166 |
| 39.6 | 0.350162 |
| 39.7 | 0.350156 |
| 39.8 | 0.350147 |
| 39.9 | 0.350136 |
| 40 | 0.350123 |
| 40.1 | 0.35011 |
| 40.2 | 0.350096 |
| 40.3 | 0.350084 |
| 40.4 | 0.350075 |
| 40.5 | 0.350072 |
| 40.6 | 0.350078 |
| 40.7 | 0.350093 |
| 40.8 | 0.350121 |
| 40.9 | 0.350161 |
| 41 | 0.35021 |
| 41.1 | 0.350265 |
| 41.2 | 0.35032 |
| 41.3 | 0.35037 |
| 41.4 | 0.350408 |
| 41.5 | 0.350433 |
| 41.6 | 0.350442 |
| 41.7 | 0.350434 |
| 41.8 | 0.350411 |
| 41.9 | 0.350374 |
| 42 | 0.350326 |
| 42.1 | 0.350271 |
| 42.2 | 0.35021 |
| 42.3 | 0.350147 |
| 42.4 | 0.350083 |
| 42.5 | 0.350019 |
| 42.6 | 0.349958 |
| 42.7 | 0.3499 |
| 42.8 | 0.349846 |
| 42.9 | 0.349798 |
| 43 | 0.349759 |
| 43.1 | 0.34973 |
| 43.2 | 0.349715 |
| 43.3 | 0.349718 |
| 43.4 | 0.349741 |
| 43.5 | 0.349788 |
| 43.6 | 0.349859 |
| 43.7 | 0.349951 |
| 43.8 | 0.35006 |
| 43.9 | 0.350177 |
| 44 | 0.350293 |
| 44.1 | 0.350398 |
| 44.2 | 0.350482 |
| 44.3 | 0.350539 |
| 44.4 | 0.350565 |
| 44.5 | 0.35056 |
| 44.6 | 0.350526 |
| 44.7 | 0.350468 |
| 44.8 | 0.350393 |
| 44.9 | 0.350308 |
| 45 | 0.350219 |
| 45.1 | 0.350129 |
| 45.2 | 0.35004 |
| 45.3 | 0.349953 |
| 45.4 | 0.349869 |
| 45.5 | 0.349792 |
| 45.6 | 0.349721 |
| 45.7 | 0.349661 |
| 45.8 | 0.349614 |
| 45.9 | 0.349583 |
| 46 | 0.34957 |
| 46.1 | 0.349576 |
| 46.2 | 0.349601 |
| 46.3 | 0.349643 |
| 46.4 | 0.3497 |
| 46.5 | 0.349767 |
| 46.6 | 0.349839 |
| 46.7 | 0.349913 |
| 46.8 | 0.349982 |
| 46.9 | 0.350043 |
| 47 | 0.350094 |
| 47.1 | 0.350134 |
| 47.2 | 0.350162 |
| 47.3 | 0.35018 |
| 47.4 | 0.350188 |
| 47.5 | 0.350188 |
| 47.6 | 0.350182 |
| 47.7 | 0.350173 |
| 47.8 | 0.350161 |
| 47.9 | 0.350148 |
| 48 | 0.350137 |
| 48.1 | 0.350129 |
| 48.2 | 0.350128 |
| 48.3 | 0.350136 |
| 48.4 | 0.350158 |
| 48.5 | 0.350195 |
| 48.6 | 0.350248 |
| 48.7 | 0.350315 |
| 48.8 | 0.350389 |
| 48.9 | 0.350465 |
| 49 | 0.350534 |
| 49.1 | 0.350592 |
| 49.2 | 0.350633 |
| 49.3 | 0.350653 |
| 49.4 | 0.350653 |
| 49.5 | 0.350632 |
| 49.6 | 0.350596 |
| 49.7 | 0.350548 |
| 49.8 | 0.350494 |
| 49.9 | 0.350435 |
| 50 | 0.350375 |
| 50.1 | 0.350314 |
| 50.2 | 0.350253 |
| 50.3 | 0.350191 |
| 50.4 | 0.350129 |
| 50.5 | 0.350068 |
| 50.6 | 0.350007 |
| 50.7 | 0.349948 |
| 50.8 | 0.349891 |
| 50.9 | 0.349837 |
| 51 | 0.349788 |
| 51.1 | 0.349743 |
| 51.2 | 0.349706 |
| 51.3 | 0.349678 |
| 51.4 | 0.349665 |
| 51.5 | 0.349672 |
| 51.6 | 0.349702 |
| 51.7 | 0.349755 |
| 51.8 | 0.349821 |
| 51.9 | 0.349889 |
| 52 | 0.349949 |
| 52.1 | 0.349989 |
| 52.2 | 0.350005 |
| 52.3 | 0.349994 |
| 52.4 | 0.349958 |
| 52.5 | 0.349901 |
| 52.6 | 0.349834 |
| 52.7 | 0.349766 |
| 52.8 | 0.3497 |
| 52.9 | 0.349635 |
| 53 | 0.349569 |
| 53.1 | 0.349502 |
| 53.2 | 0.349434 |
| 53.3 | 0.349364 |
| 53.4 | 0.349294 |
| 53.5 | 0.349225 |
| 53.6 | 0.349156 |
| 53.7 | 0.349089 |
| 53.8 | 0.349025 |
| 53.9 | 0.348967 |
| 54 | 0.348915 |
| 54.1 | 0.348872 |
| 54.2 | 0.348842 |
| 54.3 | 0.348826 |
| 54.4 | 0.348823 |
| 54.5 | 0.348833 |
| 54.6 | 0.348851 |
| 54.7 | 0.348875 |
| 54.8 | 0.348901 |
| 54.9 | 0.348925 |
| 55 | 0.348945 |
| 55.1 | 0.348959 |
| 55.2 | 0.348967 |
| 55.3 | 0.348967 |
| 55.4 | 0.348959 |
| 55.5 | 0.348943 |
| 55.6 | 0.348922 |
| 55.7 | 0.348897 |
| 55.8 | 0.348869 |
| 55.9 | 0.348838 |
| 56 | 0.348804 |
| 56.1 | 0.348767 |
| 56.2 | 0.348727 |
| 56.3 | 0.348684 |
| 56.4 | 0.348639 |
| 56.5 | 0.348595 |
| 56.6 | 0.348553 |
| 56.7 | 0.348517 |
| 56.8 | 0.348487 |
| 56.9 | 0.348466 |
| 57 | 0.348454 |
| 57.1 | 0.348449 |
| 57.2 | 0.348451 |
| 57.3 | 0.348456 |
| 57.4 | 0.348461 |
| 57.5 | 0.348464 |
| 57.6 | 0.348462 |
| 57.7 | 0.348453 |
| 57.8 | 0.348436 |
| 57.9 | 0.348413 |
| 58 | 0.348385 |
| 58.1 | 0.348354 |
| 58.2 | 0.348323 |
| 58.3 | 0.348292 |
| 58.4 | 0.348261 |
| 58.5 | 0.348231 |
| 58.6 | 0.348203 |
| 58.7 | 0.348176 |
| 58.8 | 0.348151 |
| 58.9 | 0.348128 |
| 59 | 0.348107 |
| 59.1 | 0.348088 |
| 59.2 | 0.34807 |
| 59.3 | 0.348052 |
| 59.4 | 0.348033 |
| 59.5 | 0.348013 |
| 59.6 | 0.347992 |
| 59.7 | 0.347968 |
| 59.8 | 0.347942 |
| 59.9 | 0.347915 |
| 60 | 0.347886 |
| 60.1 | 0.347855 |
| 60.2 | 0.347821 |
| 60.3 | 0.347784 |
| 60.4 | 0.347746 |
| 60.5 | 0.347707 |
| 60.6 | 0.347668 |
| 60.7 | 0.347629 |
| 60.8 | 0.347589 |
| 60.9 | 0.34755 |
| 61 | 0.34751 |
| 61.1 | 0.347469 |
| 61.2 | 0.347427 |
| 61.3 | 0.347384 |
| 61.4 | 0.34734 |
| 61.5 | 0.347295 |
| 61.6 | 0.347249 |
| 61.7 | 0.347203 |
| 61.8 | 0.347155 |
| 61.9 | 0.347109 |
| 62 | 0.347062 |
| 62.1 | 0.347017 |
| 62.2 | 0.346973 |
| 62.3 | 0.346931 |
| 62.4 | 0.34689 |
| 62.5 | 0.34685 |
| 62.6 | 0.346811 |
| 62.7 | 0.346771 |
| 62.8 | 0.34673 |
| 62.9 | 0.346688 |
| 63 | 0.346645 |
| 63.1 | 0.346601 |
| 63.2 | 0.346558 |
| 63.3 | 0.346514 |
| 63.4 | 0.346471 |
| 63.5 | 0.346427 |
| 63.6 | 0.346384 |
| 63.7 | 0.34634 |
| 63.8 | 0.346295 |
| 63.9 | 0.34625 |
| 64 | 0.346205 |
| 64.1 | 0.346159 |
| 64.2 | 0.346113 |
| 64.3 | 0.346068 |
| 64.4 | 0.346024 |
| 64.5 | 0.34598 |
| 64.6 | 0.345937 |
| 64.7 | 0.345895 |
| 64.8 | 0.345855 |
| 64.9 | 0.345817 |
| 65 | 0.345779 |
| 65.1 | 0.345742 |
| 65.2 | 0.345705 |
| 65.3 | 0.345654 |
| 65.4 | 0.345585 |
| 65.5 | 0.345516 |
| 65.6 | 0.345459 |
| 65.7 | 0.345415 |
| 65.8 | 0.345379 |
| 65.9 | 0.345349 |
| 66 | 0.345317 |
| 66.1 | 0.345278 |
| 66.2 | 0.345227 |
| 66.3 | 0.345162 |
| 66.4 | 0.345081 |
| 66.5 | 0.344989 |
| 66.6 | 0.344892 |
| 66.7 | 0.344793 |
| 66.8 | 0.344698 |
| 66.9 | 0.344607 |
| 67 | 0.344524 |
| 67.1 | 0.344448 |
| 67.2 | 0.344379 |
| 67.3 | 0.344318 |
| 67.4 | 0.344264 |
| 67.5 | 0.344217 |
| 67.6 | 0.344176 |
| 67.7 | 0.344142 |
| 67.8 | 0.344114 |
| 67.9 | 0.34409 |
| 68 | 0.344071 |
| 68.1 | 0.344056 |
| 68.2 | 0.344043 |
| 68.3 | 0.344032 |
| 68.4 | 0.34402 |
| 68.5 | 0.344007 |
| 68.6 | 0.343991 |
| 68.7 | 0.34397 |
| 68.8 | 0.34394 |
| 68.9 | 0.3439 |
| 69 | 0.343845 |
| 69.1 | 0.343772 |
| 69.2 | 0.343674 |
| 69.3 | 0.343547 |
| 69.4 | 0.343388 |
| 69.5 | 0.343195 |
| 69.6 | 0.34297 |
| 69.7 | 0.342737 |
| 69.8 | 0.34251 |
| 69.9 | 0.342294 |
| 70 | 0.342107 |
| 70.1 | 0.341957 |
| 70.2 | 0.341848 |
| 70.3 | 0.341779 |
| 70.4 | 0.341744 |
| 70.5 | 0.341735 |
| 70.6 | 0.341742 |
| 70.7 | 0.341758 |
| 70.8 | 0.341777 |
| 70.9 | 0.341798 |
| 71 | 0.341822 |
| 71.1 | 0.341849 |
| 71.2 | 0.341879 |
| 71.3 | 0.341913 |
| 71.4 | 0.341949 |
| 71.5 | 0.341984 |
| 71.6 | 0.342016 |
| 71.7 | 0.34204 |
| 71.8 | 0.342052 |
| 71.9 | 0.342047 |
| 72 | 0.342022 |
| 72.1 | 0.341975 |
| 72.2 | 0.341904 |
| 72.3 | 0.341813 |
| 72.4 | 0.341704 |
| 72.5 | 0.341582 |
| 72.6 | 0.341453 |
| 72.7 | 0.341322 |
| 72.8 | 0.341194 |
| 72.9 | 0.341073 |
| 73 | 0.340961 |
| 73.1 | 0.34086 |
| 73.2 | 0.34077 |
| 73.3 | 0.340687 |
| 73.4 | 0.340609 |
| 73.5 | 0.340533 |
| 73.6 | 0.340459 |
| 73.7 | 0.340385 |
| 73.8 | 0.340311 |
| 73.9 | 0.340236 |
| 74 | 0.340161 |
| 74.1 | 0.340086 |
| 74.2 | 0.340012 |
| 74.3 | 0.339939 |
| 74.4 | 0.339868 |
| 74.5 | 0.3398 |
| 74.6 | 0.339735 |
| 74.7 | 0.339675 |
| 74.8 | 0.339619 |
| 74.9 | 0.339568 |
| 75 | 0.339522 |
| 75.1 | 0.339481 |
| 75.2 | 0.339445 |
| 75.3 | 0.339413 |
| 75.4 | 0.339385 |
| 75.5 | 0.339362 |
| 75.6 | 0.33934 |
| 75.7 | 0.339321 |
| 75.8 | 0.339301 |
| 75.9 | 0.339281 |
| 76 | 0.339257 |
| 76.1 | 0.339229 |
| 76.2 | 0.339196 |
| 76.3 | 0.339156 |
| 76.4 | 0.339107 |
| 76.5 | 0.339049 |
| 76.6 | 0.338976 |
| 76.7 | 0.338887 |
| 76.8 | 0.33878 |
| 76.9 | 0.338651 |
| 77 | 0.338501 |
| 77.1 | 0.338331 |
| 77.2 | 0.33815 |
| 77.3 | 0.337967 |
| 77.4 | 0.337795 |
| 77.5 | 0.337643 |
| 77.6 | 0.337518 |
| 77.7 | 0.337425 |
| 77.8 | 0.337362 |
| 77.9 | 0.337329 |
| 78 | 0.337318 |
| 78.1 | 0.337325 |
| 78.2 | 0.337344 |
| 78.3 | 0.337367 |
| 78.4 | 0.337391 |
| 78.5 | 0.337416 |
| 78.6 | 0.337442 |
| 78.7 | 0.337469 |
| 78.8 | 0.337498 |
| 78.9 | 0.337526 |
| 79 | 0.337552 |
| 79.1 | 0.337575 |
| 79.2 | 0.337591 |
| 79.3 | 0.337597 |
| 79.4 | 0.337589 |
| 79.5 | 0.337563 |
| 79.6 | 0.337514 |
| 79.7 | 0.33744 |
| 79.8 | 0.337341 |
| 79.9 | 0.337219 |
| 80 | 0.33708 |
| 80.1 | 0.336929 |
| 80.2 | 0.336773 |
| 80.3 | 0.33662 |
| 80.4 | 0.336477 |
| 80.5 | 0.33635 |
| 80.6 | 0.336241 |
| 80.7 | 0.336152 |
| 80.8 | 0.336081 |
| 80.9 | 0.336026 |
| 81 | 0.335982 |
| 81.1 | 0.335947 |
| 81.2 | 0.335917 |
| 81.3 | 0.33589 |
| 81.4 | 0.335865 |
| 81.5 | 0.335841 |
| 81.6 | 0.335817 |
| 81.7 | 0.335792 |
| 81.8 | 0.335764 |
| 81.9 | 0.335734 |
| 82 | 0.335701 |
| 82.1 | 0.335663 |
| 82.2 | 0.33562 |
| 82.3 | 0.335574 |
| 82.4 | 0.335525 |
| 82.5 | 0.335475 |
| 82.6 | 0.335425 |
| 82.7 | 0.335377 |
| 82.8 | 0.335332 |
| 82.9 | 0.335291 |
| 83 | 0.335256 |
| 83.1 | 0.335224 |
| 83.2 | 0.335196 |
| 83.3 | 0.335172 |
| 83.4 | 0.335149 |
| 83.5 | 0.335127 |
| 83.6 | 0.335105 |
| 83.7 | 0.335081 |
| 83.8 | 0.335054 |
| 83.9 | 0.335024 |
| 84 | 0.334989 |
| 84.1 | 0.334947 |
| 84.2 | 0.334897 |
| 84.3 | 0.334839 |
| 84.4 | 0.33477 |
| 84.5 | 0.334691 |
| 84.6 | 0.334602 |
| 84.7 | 0.334504 |
| 84.8 | 0.334399 |
| 84.9 | 0.334291 |
| 85 | 0.334181 |
| 85.1 | 0.334074 |
| 85.2 | 0.333972 |
| 85.3 | 0.333878 |
| 85.4 | 0.333794 |
| 85.5 | 0.333721 |
| 85.6 | 0.333658 |
| 85.7 | 0.333606 |
| 85.8 | 0.333563 |
| 85.9 | 0.333528 |
| 86 | 0.333498 |
| 86.1 | 0.333473 |
| 86.2 | 0.33345 |
| 86.3 | 0.333428 |
| 86.4 | 0.333407 |
| 86.5 | 0.333386 |
| 86.6 | 0.333363 |
| 86.7 | 0.333339 |
| 86.8 | 0.333311 |
| 86.9 | 0.333279 |
| 87 | 0.333242 |
| 87.1 | 0.3332 |
| 87.2 | 0.333152 |
| 87.3 | 0.333098 |
| 87.4 | 0.333038 |
| 87.5 | 0.332974 |
| 87.6 | 0.332907 |
| 87.7 | 0.332838 |
| 87.8 | 0.332768 |
| 87.9 | 0.332701 |
| 88 | 0.332636 |
| 88.1 | 0.332575 |
| 88.2 | 0.33252 |
| 88.3 | 0.332471 |
| 88.4 | 0.332429 |
| 88.5 | 0.332392 |
| 88.6 | 0.33236 |
| 88.7 | 0.332331 |
| 88.8 | 0.332303 |
| 88.9 | 0.332275 |
| 89 | 0.332244 |
| 89.1 | 0.332212 |
| 89.2 | 0.332175 |
| 89.3 | 0.332134 |
| 89.4 | 0.332088 |
| 89.5 | 0.332035 |
| 89.6 | 0.331974 |
| 89.7 | 0.331904 |
| 89.8 | 0.331824 |
| 89.9 | 0.331732 |
| 90 | 0.33163 |
| 90.1 | 0.331517 |
| 90.2 | 0.331395 |
| 90.3 | 0.331267 |
| 90.4 | 0.331139 |
| 90.5 | 0.331014 |
| 90.6 | 0.330898 |
| 90.7 | 0.330794 |
| 90.8 | 0.330705 |
| 90.9 | 0.330631 |
| 91 | 0.330572 |
| 91.1 | 0.330526 |
| 91.2 | 0.330492 |
| 91.3 | 0.330466 |
| 91.4 | 0.330446 |
| 91.5 | 0.330427 |
| 91.6 | 0.330408 |
| 91.7 | 0.33039 |
| 91.8 | 0.330374 |
| 91.9 | 0.330358 |
| 92 | 0.330343 |
| 92.1 | 0.330326 |
| 92.2 | 0.330306 |
| 92.3 | 0.330282 |
| 92.4 | 0.330252 |
| 92.5 | 0.330215 |
| 92.6 | 0.330171 |
| 92.7 | 0.330119 |
| 92.8 | 0.330058 |
| 92.9 | 0.32999 |
| 93 | 0.329915 |
| 93.1 | 0.329835 |
| 93.2 | 0.329751 |
| 93.3 | 0.329665 |
| 93.4 | 0.329578 |
| 93.5 | 0.329493 |
| 93.6 | 0.32941 |
| 93.7 | 0.329331 |
| 93.8 | 0.329255 |
| 93.9 | 0.329183 |
| 94 | 0.329112 |
| 94.1 | 0.329044 |
| 94.2 | 0.328976 |
| 94.3 | 0.328908 |
| 94.4 | 0.328839 |
| 94.5 | 0.328769 |
| 94.6 | 0.328698 |
| 94.7 | 0.328625 |
| 94.8 | 0.328549 |
| 94.9 | 0.32847 |
| 95 | 0.328389 |
| 95.1 | 0.328306 |
| 95.2 | 0.328222 |
| 95.3 | 0.328137 |
| 95.4 | 0.328054 |
| 95.5 | 0.327975 |
| 95.6 | 0.3279 |
| 95.7 | 0.327831 |
| 95.8 | 0.327769 |
| 95.9 | 0.327714 |
| 96 | 0.327666 |
| 96.1 | 0.327622 |
| 96.2 | 0.327582 |
| 96.3 | 0.327544 |
| 96.4 | 0.327508 |
| 96.5 | 0.327473 |
| 96.6 | 0.327439 |
| 96.7 | 0.327406 |
| 96.8 | 0.327375 |
| 96.9 | 0.327345 |
| 97 | 0.327315 |
| 97.1 | 0.327284 |
| 97.2 | 0.327249 |
| 97.3 | 0.32721 |
| 97.4 | 0.327163 |
| 97.5 | 0.327108 |
| 97.6 | 0.327042 |
| 97.7 | 0.326966 |
| 97.8 | 0.326877 |
| 97.9 | 0.326776 |
| 98 | 0.326663 |
| 98.1 | 0.326538 |
| 98.2 | 0.326404 |
| 98.3 | 0.326265 |
| 98.4 | 0.326122 |
| 98.5 | 0.32598 |
| 98.6 | 0.325842 |
| 98.7 | 0.325711 |
| 98.8 | 0.32559 |
| 98.9 | 0.325482 |
| 99 | 0.325386 |
| 99.1 | 0.325303 |
| 99.2 | 0.325231 |
| 99.3 | 0.325168 |
| 99.4 | 0.325114 |
| 99.5 | 0.325066 |
| 99.6 | 0.325023 |
| 99.7 | 0.324984 |
| 99.8 | 0.324947 |
| 99.9 | 0.324913 |
| 100 | 0.32488 |
| 100.1 | 0.324847 |
| 100.2 | 0.324815 |
| 100.3 | 0.324784 |
| 100.4 | 0.324754 |
| 100.5 | 0.324726 |
| 100.6 | 0.324697 |
| 100.7 | 0.324669 |
| 100.8 | 0.32464 |
| 100.9 | 0.324609 |
| 101 | 0.324576 |
| 101.1 | 0.324539 |
| 101.2 | 0.324497 |
| 101.3 | 0.32445 |
| 101.4 | 0.324398 |
| 101.5 | 0.324341 |
| 101.6 | 0.324279 |
| 101.7 | 0.324215 |
| 101.8 | 0.324151 |
| 101.9 | 0.324088 |
| 102 | 0.324027 |
| 102.1 | 0.323969 |
| 102.2 | 0.323913 |
| 102.3 | 0.323857 |
| 102.4 | 0.323802 |
| 102.5 | 0.323744 |
| 102.6 | 0.323682 |
| 102.7 | 0.323616 |
| 102.8 | 0.323543 |
| 102.9 | 0.323464 |
| 103 | 0.323377 |
| 103.1 | 0.323282 |
| 103.2 | 0.323177 |
| 103.3 | 0.323065 |
| 103.4 | 0.322944 |
| 103.5 | 0.322816 |
| 103.6 | 0.322682 |
| 103.7 | 0.322546 |
| 103.8 | 0.322409 |
| 103.9 | 0.322275 |
| 104 | 0.322147 |
| 104.1 | 0.322028 |
| 104.2 | 0.321921 |
| 104.3 | 0.321826 |
| 104.4 | 0.321746 |
| 104.5 | 0.32168 |
| 104.6 | 0.321627 |
| 104.7 | 0.321586 |
| 104.8 | 0.321557 |
| 104.9 | 0.321537 |
| 105 | 0.321526 |
| 105.1 | 0.321521 |
| 105.2 | 0.321521 |
| 105.3 | 0.321523 |
| 105.4 | 0.321527 |
| 105.5 | 0.321532 |
| 105.6 | 0.321535 |
| 105.7 | 0.321537 |
| 105.8 | 0.321537 |
| 105.9 | 0.321534 |
| 106 | 0.321525 |
| 106.1 | 0.321512 |
| 106.2 | 0.321494 |
| 106.3 | 0.321469 |
| 106.4 | 0.321438 |
| 106.5 | 0.3214 |
| 106.6 | 0.321356 |
| 106.7 | 0.321306 |
| 106.8 | 0.321252 |
| 106.9 | 0.321193 |
| 107 | 0.321132 |
| 107.1 | 0.32107 |
| 107.2 | 0.321009 |
| 107.3 | 0.320949 |
| 107.4 | 0.320891 |
| 107.5 | 0.320836 |
| 107.6 | 0.320784 |
| 107.7 | 0.320735 |
| 107.8 | 0.320689 |
| 107.9 | 0.320644 |
| 108 | 0.3206 |
| 108.1 | 0.320555 |
| 108.2 | 0.32051 |
| 108.3 | 0.320464 |
| 108.4 | 0.320416 |
| 108.5 | 0.320366 |
| 108.6 | 0.320313 |
| 108.7 | 0.320258 |
| 108.8 | 0.320199 |
| 108.9 | 0.320136 |
| 109 | 0.320071 |
| 109.1 | 0.320002 |
| 109.2 | 0.31993 |
| 109.3 | 0.319856 |
| 109.4 | 0.319781 |
| 109.5 | 0.319706 |
| 109.6 | 0.319633 |
| 109.7 | 0.319563 |
| 109.8 | 0.319497 |
| 109.9 | 0.319438 |
| 110 | 0.319385 |
| 110.1 | 0.31934 |
| 110.2 | 0.319303 |
| 110.3 | 0.319274 |
| 110.4 | 0.319251 |
| 110.5 | 0.319235 |
| 110.6 | 0.319224 |
| 110.7 | 0.319217 |
| 110.8 | 0.319214 |
| 110.9 | 0.319212 |
| 111 | 0.319212 |
| 111.1 | 0.319212 |
| 111.2 | 0.319213 |
| 111.3 | 0.319214 |
| 111.4 | 0.319214 |
| 111.5 | 0.319213 |
| 111.6 | 0.319211 |
| 111.7 | 0.319206 |
| 111.8 | 0.319199 |
| 111.9 | 0.319189 |
| 112 | 0.319176 |
| 112.1 | 0.319159 |
| 112.2 | 0.319139 |
| 112.3 | 0.319115 |
| 112.4 | 0.319088 |
| 112.5 | 0.319059 |
| 112.6 | 0.319026 |
| 112.7 | 0.318992 |
| 112.8 | 0.318956 |
| 112.9 | 0.31892 |
| 113 | 0.318885 |
| 113.1 | 0.318851 |
| 113.2 | 0.318819 |
| 113.3 | 0.318788 |
| 113.4 | 0.318761 |
| 113.5 | 0.318736 |
| 113.6 | 0.318713 |
| 113.7 | 0.318693 |
| 113.8 | 0.318676 |
| 113.9 | 0.31866 |
| 114 | 0.318647 |
| 114.1 | 0.318637 |
| 114.2 | 0.318628 |
| 114.3 | 0.318621 |
| 114.4 | 0.318616 |
| 114.5 | 0.318611 |
| 114.6 | 0.318605 |
| 114.7 | 0.318598 |
| 114.8 | 0.318588 |
| 114.9 | 0.318576 |
| 115 | 0.318561 |
| 115.1 | 0.318543 |
| 115.2 | 0.318521 |
| 115.3 | 0.318495 |
| 115.4 | 0.318467 |
| 115.5 | 0.318436 |
| 115.6 | 0.318403 |
| 115.7 | 0.31837 |
| 115.8 | 0.318336 |
| 115.9 | 0.318302 |
| 116 | 0.31827 |
| 116.1 | 0.31824 |
| 116.2 | 0.318213 |
| 116.3 | 0.318189 |
| 116.4 | 0.318168 |
| 116.5 | 0.31815 |
| 116.6 | 0.318134 |
| 116.7 | 0.318121 |
| 116.8 | 0.31811 |
| 116.9 | 0.318101 |
| 117 | 0.318093 |
| 117.1 | 0.318085 |
| 117.2 | 0.318078 |
| 117.3 | 0.31807 |
| 117.4 | 0.318062 |
| 117.5 | 0.318054 |
| 117.6 | 0.318044 |
| 117.7 | 0.318034 |
| 117.8 | 0.318023 |
| 117.9 | 0.318012 |
| 118 | 0.318001 |
| 118.1 | 0.317989 |
| 118.2 | 0.317978 |
| 118.3 | 0.317966 |
| 118.4 | 0.317955 |
| 118.5 | 0.317944 |
| 118.6 | 0.317935 |
| 118.7 | 0.317927 |
| 118.8 | 0.317922 |
| 118.9 | 0.317919 |
| 119 | 0.31792 |
| 119.1 | 0.317925 |
| 119.2 | 0.317935 |
| 119.3 | 0.317949 |
| 119.4 | 0.317968 |
| 119.5 | 0.31799 |
| 119.6 | 0.318016 |
| 119.7 | 0.318045 |
| 119.8 | 0.318074 |
| 119.9 | 0.318103 |
| 120 | 0.31813 |
| 120.1 | 0.318152 |
| 120.2 | 0.318171 |
| 120.3 | 0.318185 |
| 120.4 | 0.318198 |
| 120.5 | 0.31821 |
| 120.6 | 0.318223 |
| 120.7 | 0.318234 |
| 120.8 | 0.318245 |
| 120.9 | 0.318253 |
| 121 | 0.318259 |
| 121.1 | 0.318262 |
| 121.2 | 0.31826 |
| 121.3 | 0.318254 |
| 121.4 | 0.318243 |
| 121.5 | 0.318225 |
| 121.6 | 0.318202 |
| 121.7 | 0.318172 |
| 121.8 | 0.318135 |
| 121.9 | 0.318093 |
| 122 | 0.318046 |
| 122.1 | 0.317994 |
| 122.2 | 0.31794 |
| 122.3 | 0.317884 |
| 122.4 | 0.317829 |
| 122.5 | 0.317776 |
| 122.6 | 0.317726 |
| 122.7 | 0.317681 |
| 122.8 | 0.317643 |
| 122.9 | 0.317612 |
| 123 | 0.317588 |
| 123.1 | 0.31757 |
| 123.2 | 0.317559 |
| 123.3 | 0.317554 |
| 123.4 | 0.317554 |
| 123.5 | 0.317559 |
| 123.6 | 0.317569 |
| 123.7 | 0.317583 |
| 123.8 | 0.317602 |
| 123.9 | 0.317625 |
| 124 | 0.317653 |
| 124.1 | 0.317686 |
| 124.2 | 0.317723 |
| 124.3 | 0.317765 |
| 124.4 | 0.317812 |
| 124.5 | 0.317863 |
| 124.6 | 0.317918 |
| 124.7 | 0.317976 |
| 124.8 | 0.318038 |
konishi
Masashi Imano
Jan 4, 2013, 1:40:00 AM1/4/13
to open...@googlegroups.com
今野です。
> 下記1, 2の設定では解析を実行出来ましたが、3のfvSolutionを用いると必ずinterpolateSchemesのエラーが表示され解析が停止する所迄分かりました事をご報告致します。
interpolateSchemesのエラーメッセージを貼ってください。
--
IMANO Masashi, Ph.D.
> 下記1, 2の設定では解析を実行出来ましたが、3のfvSolutionを用いると必ずinterpolateSchemesのエラーが表示され解析が停止する所迄分かりました事をご報告致します。
interpolateSchemesのエラーメッセージを貼ってください。
2013年1月4日 14:22 Konishi <ryo.kon...@gmail.com>:
konishi--
このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
このディスカッションをウェブ上で閲覧するには、https://groups.google.com/d/msg/openfoam/-/tnkkZDYxHJEJ にアクセスしてください。
このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
このグループから退会するには、openfoam+u...@googlegroups.com にメールを送信してください。
詳細については、http://groups.google.com/group/openfoam?hl=ja からこのグループにアクセスしてください。
IMANO Masashi, Ph.D.
Konishi
Jan 4, 2013, 1:46:30 AM1/4/13
to open...@googlegroups.com
今野先生
お世話に成って居ります。konishiで御座います。
interpolateSchemesのエラーメッセージは以下の様に成って居ります。これはpotentialFoam、simpleFoam共に同じエラーが吐き出されます。また、各行頭の[](分割領域の番号でしょうか)は解析を実行する毎に異なります。
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
relaxation factors:
fields:
{
p 0.05;
}
equations:
{
p 0.05;
U 0.15;
k 0.15;
epsilon 0.15;
omega 0.15;
}
Reading field p
Reading field U
[33]
[33]
[33] --> FOAM FATAL IO ERROR:
[33] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[33]
[33] file: IOstream::interpolationSchemes
[33]
[33] From function dictionary::lookupEntry(const word&, bool, bool) const
[33] in file db/dictionary/dictionary.C at line 396.
[33]
FOAM parallel run exiting
[33]
[29]
[29]
[29] --> FOAM FATAL IO ERROR:
[29] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[29]
[29] file: IOstream::interpolationSchemes
[29]
[29] From function dictionary::lookupEntry(const word&, bool, bool) const
[29] in file db/dictionary/dictionary.C at line 396.
[29]
FOAM parallel run exiting
[29]
[42]
[42]
[42] --> FOAM FATAL IO ERROR:
[42] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[42]
[42] file: IOstream::interpolationSchemes
[42]
[42] From function dictionary::lookupEntry(const word&, bool, bool) const
[42] in file db/dictionary/dictionary.C at line 396.
[42]
FOAM parallel run exiting
[42]
[68]
[68]
[68] --> FOAM FATAL IO ERROR:
[68] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[68]
[68] file: IOstream::interpolationSchemes
[68]
[68] From function dictionary::lookupEntry(const word&, bool, bool) const
[68] in file db/dictionary/dictionary.C at line 396.
[68]
FOAM parallel run exiting
[68]
[43]
[43]
[43] --> FOAM FATAL IO ERROR:
[43] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[43]
[43] file: IOstream::interpolationSchemes
[43]
[43] From function dictionary::lookupEntry(const word&, bool, bool) const
[43] in file db/dictionary/dictionary.C at line 396.
[43]
FOAM parallel run exiting
[43]
[9]
[9]
[9] --> FOAM FATAL IO ERROR:
[9] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[9]
[9] file: IOstream::interpolationSchemes
[9]
[9] From function dictionary::lookupEntry(const word&, bool, bool) const
[9] in file db/dictionary/dictionary.C at line 396.
[9]
FOAM parallel run exiting
[9]
[23]
[23]
[23] --> FOAM FATAL IO ERROR:
[23] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[23]
[23] file: IOstream::interpolationSchemes
[23]
[23] From function dictionary::lookupEntry(const [22wor] d&, bool, bool
[22) const]
[23[22] --] > F OAM i FATAn fiL Ile O ERRdbOR:
/dictio[22] keynary/diwordction inaryterpo.late(CU) at lis unine define396.
d in[23] di
FOActionM parallelary ru "In eOstxiting
[ream::23]
interpolationSchemes"
[22]
[22] file: IOstream::interpolationSchemes
[22]
[22] From function dictionary::lookupEntry(const word&, bool, bool) const
[22] in file db/dictionary/dictionary.C at line 396.
[22]
FOAM parallel run exiting
[22]
[21]
[21]
[21] --> FOAM FATAL IO ERROR:
[21] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[21]
[21] file: IOstream::interpolationSchemes
[21]
[21] From function dictionary::lookupEntry(const word&, bool, bool) const
[21] in file db/dictionary/dictionary.C at line 396.
[21]
FOAM parallel run exiting
[21]
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 43 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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 42 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.
--------------------------------------------------------------------------
[50]
[50]
[50] --> FOAM FATAL IO ERROR:
[50] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[50]
[50] file: IOstream::interpolationSchemes
[50]
[50] From function dictionary::lookupEntry(const word&, bool, bool) const
[50] in file db/dictionary/dictionary.C at line 396.
[50]
FOAM parallel run exiting
[50]
[58]
[58]
[58] --> FOAM FATAL IO ERROR:
[58] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[58]
[58] file: IOstream::interpolationSchemes
[58]
[58] From function dictionary::lookupEntry(const word&, bool, bool) const
[58] in file db/dictionary/dictionary.C at line 396.
[58]
FOAM parallel run exiting
[58]
[5]
[5]
[5] --> FOAM FATAL IO ERROR:
[5] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[5]
[5] file: IOstream::interpolationSchemes
[5]
[5] From function dictionary::lookupEntry(const word&, bool, bool) const
[5] in file db/dictionary/dictionary.C at line 396.
[5]
FOAM parallel run exiting
[5]
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 68 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.
--------------------------------------------------------------------------
[39]
[39]
[39] --> FOAM FATAL IO ERROR:
[39] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[39]
[39] file: IOstream::interpolationSchemes
[39]
[39] From function dictionary::lookupEntry(const word&, bool, bool) const
[39] in file db/dictionary/dictionary.C at line 396.
[39]
FOAM parallel run exiting
[39]
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 39 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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 33 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.
--------------------------------------------------------------------------
[32]
[32]
[32] --> FOAM FATAL IO ERROR:
[32] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[32]
[32] file: IOstream::--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 29 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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 9 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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 21 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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 23 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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 22 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.
--------------------------------------------------------------------------
[13]
[13]
[13] --> FOAM FATAL IO ERROR:
[13] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[13]
[13] file: IOstream::interpolationSchemes
[13]
[13] From function dictionary::lookupEntry(const word&, bool, bool) const
[13] in file db/dictionary/dictionary.C at line 396.
[13]
FOAM parallel run exiting
[13]
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 13 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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 50 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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 58 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.
--------------------------------------------------------------------------
[65]
[65]
[65] --> FOAM FATAL IO ERROR:
[65] keyword interpolate(U) is undefined in dictionary "IOstream::interpolationSchemes"
[65]
[65] file: IOstream::interpolationSchemes
[65]
[65] From function dictionary::lookupEntry(const word&, bool, bool) const
[65] in file db/dictionary/dictionary.C at line 396.
[65]
FOAM parallel run exiting
[65]
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 65 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.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
[mpi::mpi-api::mpi-abort]
MPI_ABORT was invoked on rank 5 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.
--------------------------------------------------------------------------
[ERR.] PLE 0019 plexec One of MPI processes was aborted.(rank=5)(nid=0x02010064)(CODE=1878,793702259720650752,256)
konishi
Masashi Imano
Jan 4, 2013, 2:19:39 AM1/4/13
to open...@googlegroups.com
今野です。
3のfvSolutionを用いると何故以下のエラーが出るかわかりません。
特に、途中でメッセージが乱れているところが不可解です。
これについては情報基盤センターに問い合わせてみてください。
2のfvSolutionで動くのであれば、こちらのほうが収束性が良く実際的なので、
当面これを用いて解析すれば良いと思います。
動くようであれば、fvSchemesでの以下のupwindスキームについては
より高精度なスキームに変更していってください。
divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,omega) Gauss upwind;
なお、細い格子での本計算を実行する前には、研究室マシンでの
計算値と合うことは確認したほうが良いと思います。
以上です。
卒論がんばってください。
2013/1/4 Konishi <ryo.kon...@gmail.com>:
3のfvSolutionを用いると何故以下のエラーが出るかわかりません。
特に、途中でメッセージが乱れているところが不可解です。
これについては情報基盤センターに問い合わせてみてください。
2のfvSolutionで動くのであれば、こちらのほうが収束性が良く実際的なので、
当面これを用いて解析すれば良いと思います。
動くようであれば、fvSchemesでの以下のupwindスキームについては
より高精度なスキームに変更していってください。
divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,omega) Gauss upwind;
計算値と合うことは確認したほうが良いと思います。
以上です。
卒論がんばってください。
2013/1/4 Konishi <ryo.kon...@gmail.com>:
> --
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> このディスカッションをウェブ上で閲覧するには、https://groups.google.com/d/msg/openfoam/-/iy5mUNZUjQsJ
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
Konishi
Jan 4, 2013, 6:28:35 AM1/4/13
to open...@googlegroups.com
今野先生
お世話に成って居ります。konishiで御座います。
情報基盤センターへの問い合わせ、Gauss upwindへの変更、所属研究室による解析値の確認、承知致しました。
長きに渡る厚い御指導、本当に有難う御座いました。目標達成へ向け引き続き邁進致します。
今後とも宜しくお願い申し上げます。
konishi
Reply all
Reply to author
Forward
0 new messages