用的windAroundBuildings算例，840万网格的二维模型，网格质量检查过了，计算域是14000*4000米。epsilonInlet 0.0003; // Cmu^0.75 * k^1.5 / L ; L =10
dimensions [0 2 -3 0 0 0 0];
internalField uniform $epsilonInlet;
boundaryField
{
INLET
{
type fixedValue;
value uniform $epsilonInlet;
}
OUTLET
{
type zeroGradient;
}
WALL
{
type epsilonWallFunction;
value uniform $epsilonInlet;
}
frontAndBackPlanes
{
type empty;
}
#includeEtc "caseDicts/setConstraintTypes"
}
kInlet 1.5; // approx k = 1.5*(I*U)^2 ; I = 0.1
dimensions [0 2 -2 0 0 0 0];
internalField uniform $kInlet;
boundaryField
{
INLET
{
type fixedValue;
value uniform $kInlet;
}
OUTLET
{
type zeroGradient;
}
WALL
{
type kqRWallFunction;
value uniform $kInlet;
}
frontAndBackPlanes
{
type empty;
}
#includeEtc "caseDicts/setConstraintTypes"
}
nut
internalField uniform 0;
boundaryField
{
INLET
{
type calculated;
value uniform 0;
}
OUTLET
{
type calculated;
value uniform 0;
}
WALL
{
type nutUSpaldingWallFunction;
value uniform 0;
}
frontAndBackPlanes
{
type empty;
}
#includeEtc "caseDicts/setConstraintTypes"
}
p
internalField uniform 0;
boundaryField
{
INLET
{
type zeroGradient;
}
OUTLET
{
type fixedValue;
value uniform 0;
p0 uniform 0;
}
WALL
{
type zeroGradient;
}
frontAndBackPlanes
{
type empty;
}
#includeEtc "caseDicts/setConstraintTypes"
}
Uinlet (10 0 0);
dimensions [0 1 -1 0 0 0 0];
internalField uniform (10 0 0);
boundaryField
{
INLET
{
type fixedValue;
value uniform $Uinlet;
}
OUTLET
{
type zeroGradient;
}
WALL
{
type noSlip;
}
frontAndBackPlanes
{
type empty;
}
#includeEtc "caseDicts/setConstraintTypes"
}
时间步长试过1 0.1 0.01 其他的都没动
