关于openfoam网格转换

Afr1yne

请教下，网格用icem画好转换过来的，还需要blockMesh和snappyHexMesh吗？

李东岳

不需要

Afr1yne

@李东岳李老师我在做一个稳态的算例但是算出来的结果总是很离谱不收敛可以指导下吗

Afr1yne

用的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 其他的都没动

CFD中文网

关于openfoam网格转换