H
重新画了mesh，模型改成了sst。Re 5w
结果更奇怪了，直接算不下去。下图无量纲时间步
现在mesh的y+平均1.5，max 3，min0.11。应该是没什么问题，质量也很高。
上下cyclicAMI边界，圆形流域，spanwise piD。
没想到一个RANS花了这么多天还没算好。。
有大佬能帮我看看么
k
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0.0021;
boundaryField
{
front
{
type cyclicAMI;
value uniform 0.0021;
}
back
{
type cyclicAMI;
value uniform 0.0021;
}
inlet
{
type fixedValue;
value uniform 0.0021;
}
outlet
{
type inletOutlet;
inletValue internalField;
value $internalField;
}
cylinderwalls
{
type kqRWallFunction;
value uniform 0.0021;
}
}
omega
dimensions [ 0 0 -1 0 0 0 0 ];
internalField uniform 10.9;
boundaryField
{
cylinderwalls
{
type omegaWallFunction;
value $internalField;
}
front
{
type cyclicAMI;
}
back
{
type cyclicAMI;
}
inlet
{
type fixedValue;
value uniform 10.9;
}
outlet
{
type inletOutlet;
inletValue internalField;
value $internalField;
}
}
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;
}
cylinderwalls
{
type nutkWallFunction;
value uniform 0;
}
front
{
type cyclicAMI;
}
back
{
type cyclicAMI;
}
}
p
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
front
{
type cyclicAMI;
value uniform 0;
}
back
{
type cyclicAMI;
value uniform 0;
}
inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value uniform 0;
}
cylinderwalls
{
type zeroGradient;
}
}
U
dimensions [0 1 -1 0 0 0 0];
internalField uniform (7.5 0 0);
boundaryField
{
front
{
type cyclicAMI;
value uniform (7.5 0 0);
}
back
{
type cyclicAMI;
value uniform (7.5 0 0);
}
inlet
{
type fixedValue;
value uniform (7.5 0 0);
}
outlet
{
type zeroGradient;
}
cylinderwalls
{
type fixedValue;
value uniform (0 0 0);
}
}
fvscheme
ddtSchemes
{
default backward;
}
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}
divSchemes
{
default none;
div(phi,U) Gauss limitedLinearV 0.5;
div(phi,k) Gauss limitedLinear 1;
div(phi,omega) Gauss limitedLinear 1;
div(phi,R) Gauss limitedLinear 1;
div(R) Gauss linear;
div(phi,B) Gauss limitedLinear 1;
div(B) Gauss linear;
div(phi,nuTilda) Gauss limitedLinear 1;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
interpolate(U) linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p;
Phi;
}
fvsolution
solvers
{
p
{
solver GAMG;
tolerance 1e-04;
relTol 0;
smoother GaussSeidel;
cacheAgglomeration true;
nCellsInCoarsestLevel 200;
agglomerator faceAreaPair;
mergeLevels 1;
}
pFinal
{
solver GAMG;
tolerance 1e-06;
relTol 0;
smoother GaussSeidel;
cacheAgglomeration true;
nCellsInCoarsestLevel 200;
agglomerator faceAreaPair;
mergeLevels 1;
}
"(U|k|nut|omega)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}
"(U|k|nut|omega)Final"
{
$U;
tolerance 1e-06;
relTol 0;
}
}
PIMPLE
{
nOuterCorrectors 3;
nCorrectors 2;
nNonOrthogonalCorrectors 1;
pRefCell 0;
pRefValue 0;
}
