请问在使用turbulentDFSEMInlet入口类型计算过程中,随着计算的进行,后半部分的风场慢慢的变得不是对数廓线了,请问是什么原因,有没有大佬知道。
我的初始边界信息:
计算域:xyz为1011(m)
U:
inlet
{
type turbulentDFSEMInlet;
delta 1;
U
{
type mappedFile;
mapMethod nearest;
}
R uniform (0.8 0 0 0.4 0 0.4);
L uniform 0.004;
d 1;
nCellPerEddy 1;
scale 1;
value $internalField;
}
outlet
{
type advective;
value $internalField;
}
"(topWall|front|back)"
{
type zeroGradient;
}
"(bottomWall|Model)"
{
type noSlip;
}
输入文件:
point:
(
( 0 0 0.01 )
( 0 0 0.05 )
( 0 0 0.1 )
( 0 0 0.15 )
( 0 0 0.2 )
( 0 0 0.25 )
( 0 0 0.3 )
( 0 0 0.35 )
( 0 0 0.4 )
( 0 0 0.45 )
( 0 0 0.5 )
( 0 0 0.55 )
( 0 0 0.6 )
( 0 0 0.65 )
( 0 0 0.7 )
( 0 0 0.75 )
( 0 0 0.8 )
( 0 0 0.85 )
( 0 0 0.9 )
( 0 0 0.95 )
( 0 0 0.99 )
)
U:
(
( 0 0 0 )
( 1.048097372 0 0 )
( 1.499488341 0 0 )
( 1.763535131 0 0 )
( 1.95087931 0 0 )
( 2.096194744 0 0 )
( 2.2149261 0 0 )
( 2.31531203 0 0 )
( 2.402270279 0 0 )
( 2.47897289 0 0 )
( 2.547585713 0 0 )
( 2.609653562 0 0 )
( 2.666317069 0 0 )
( 2.718442442 0 0 )
( 2.766702999 0 0 )
( 2.811632503 0 0 )
( 2.853661247 0 0 )
( 2.893141184 0 0 )
( 2.930363858 0 0 )
( 2.965573488 0 0 )
( 2.992431707 0 0 )
)
计算3600s后,平均风速廓线,U1U2U3U4U5U6U7U8U9的位置为0.1,2,4,5,6,7,8,9,9.9
e595f49f-347a-4180-82e2-b6295cd45a15-图片.png