计算结尾如下:
边界条件:
流体:
epsilon 边界为
INLET
{
type inletOutlet;
inletValue uniform 0.02547;
value uniform 0;
}
OUTLET
{
type inletOutlet;
inletValue uniform 0.02547;
value uniform 0;
}
FRONT_AND_BACK
{
type zeroGradient;
}
UP_AND_DOWN
{
type zeroGradient;
}
INTERFACE_INTERFACE
{
type epsilonWallFunction;
value uniform 0.02547;
}
k 的边界如下:
internalField uniform 0.03375;
boundaryField
{
INLET
{
type inletOutlet;
inletValue uniform 0.03375;
value uniform 0;
}
OUTLET
{
type inletOutlet;
inletValue uniform 0.03375;
value uniform 0;
}
FRONT_AND_BACK
{
type zeroGradient;
value uniform 0;
}
UP_AND_DOWN
{
type zeroGradient;
value uniform 0;
}
INTERFACE_INTERFACE
{
type kqRWallFunction;
value uniform 0.03375;
}
nut 的边界如下:
INLET
{
type calculated;
//value uniform 0;
}
OUTLET
{
type calculated;
// inletValue uniform 0;
//value uniform 0;
}
FRONT_AND_BACK
{
type zeroGradient;
value uniform 0;
}
UP_AND_DOWN
{
type zeroGradient;
value uniform 0;
}
INTERFACE_INTERFACE
{
type nutkWallFunction;
value uniform 0;
}
nutilda 边界为:
INLET
{
type zeroGradient;
}
OUTLET
{
type zeroGradient;
}
FRONT_AND_BACK
{
type zeroGradient;
}
UP_AND_DOWN
{
type zeroGradient;
}
INTERFACE_INTERFACE
{
type zeroGradient;
}
p 的边界为:
INLET
{
type zeroGradient;
value uniform 0;
}
OUTLET
{
type fixedValue;
value uniform 0;
}
FRONT_AND_BACK
{
type zeroGradient;
}
UP_AND_DOWN
{
type zeroGradient;
}
INTERFACE_INTERFACE
{
type zeroGradient;
value uniform 0;
}
pointU 边界为:
INLET
{
type fixedValue;
value uniform (0 0 0);
}
OUTLET
{
type zeroGradient;
value uniform (0 0 0);
}
FRONT_AND_BACK
{
type fixedValue;
value uniform (0 0 0);
}
UP_AND_DOWN
{
type fixedValue;
value uniform (0 0 0);
}
INTERFACE_INTERFACE
{
type fixedValue;
value uniform (0 0 0);
}
INLET
{
type kqRWallFunction;
value uniform ( 0 0 0 0 0 0 );
}
OUTLET
{
type kqRWallFunction;
value uniform ( 0 0 0 0 0 0 );
}
FRONT_AND_BACK
{
type kqRWallFunction;
value uniform ( 0 0 0 0 0 0 );
}
UP_AND_DOWN
{
type kqRWallFunction;
value uniform ( 0 0 0 0 0 0 );
}
INTERFACE_INTERFACE
{
type kqRWallFunction;
value uniform ( 0 0 0 0 0 0 );
}
U 的边界为:
internalField uniform (3 0 0);
boundaryField
{
INLET
{
type fixedValue;
value uniform (3 0 0);
}
OUTLET
{
type zeroGradient;
value uniform (0 0 0);
}
FRONT_AND_BACK
{
type fixedValue;
value uniform (0 0 0);
}
UP_AND_DOWN
{
type fixedValue;
value uniform (0 0 0);
}
INTERFACE_INTERFACE
{
type movingWallVelocity;
value uniform (0 0 0);
}
固体的边界如下:
d 文件为
OUTERWALL
{
type tractionDisplacement;
traction uniform ( 0 0 0 );
pressure uniform 0;
value uniform (0 0 0);
}
INNERWALL
{
type tractionDisplacement;
traction uniform ( 0 0 0 );
pressure uniform 0;
value uniform (0 0 0);
}
UP_AND_DOWN
{
type fixedValue;
value uniform (0 0 0);
}
dd 文件为:
OUTERWALL
{
type tractionDisplacementIncrement;
traction uniform ( 0 0 0 );
pressure uniform 0;
value uniform (0 0 0);
}
INNERERWALL
{
type tractionDisplacementIncrement;
traction uniform ( 0 0 0 );
pressure uniform 0;
value uniform (0 0 0);
}
UP_AND_DOWN
{
type fixedValue;
value uniform (0 0 0);
}
pointD 为:
OUTERWALL
{
type calculated;
value uniform (0 0 0);
}
INNERWALL
{
type calculated;
value uniform (0 0 0);
}
UP_AND_DOWN
{
type fixedValue;
value uniform (0 0 0);
}
pointdd 为:
dimensions [0 1 0 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
OUTERWALL
{
type calculated;
value uniform (0 0 0);
}
INNERWALL
{
type calculated;
value uniform (0 0 0);
}
UP_AND_DOWN
{
type fixedValue;
value uniform (0 0 0);
}
controlDict 为:
application fsiFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 0.415;
deltaT 2.5e-4;
writeControl runTime;
writeInterval 5e-4;
purgeWrite 0;
writeFormat ascii;
writePrecision 8;
writeCompression uncompressed;
timeFormat general;
timePrecision 6;
runTimeModifiable yes;
adjustTimeStep yes;
maxCo 0.01;
maxdeltaT 5e-4;
functions
(
pointHistory
{
type pointHistory;
functionObjectLibs
(
"libpointHistory.so"
);
region solid;
refHistoryPoint (0 0.005 0.025);
}
);
边界条件就是这样,可能存在的问题有 边界错误,k epsilon值错误,松弛因子 等还有计算过程中提醒壁面函数边界条件没有G这个参数等等。
另外目前还需要获取升力阻力及系数的提取数据函数。
同样的条件实心管的流固耦合,也就是圆柱绕流却收敛了