CFDEM里,在进行两相耦合求解时,需要在CFD/constant文件夹下设置coupingProperties文件,本文给出这个文件里一些参数的意义。
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.4 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object couplingProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
//===========================================================================//
// sub-models & settings
modelType "A"; // 耦合模型 有A类何B类
couplingInterval 100;//cfd与dem数据交互的间隔,几个DEM步长交换一次
voidFractionModel divided;//centre,计算提及分数的方式divide用于颗粒何流体网格大小接近的情况
locateModel engine;//engineIB;//定位颗粒在流体网格的模型
meshMotionModel noMeshMotion;//根据DEM网格的移动旋转CFD网格
IOModel basicIO;//basicIO把颗粒的位置速度何直径信息在每个CFD步长写道文件CFD/prco*/time/lagrangian中,sophIO除了以上还写体积分数,作用力。
probeModel off;
dataExchangeModel twoWayMPI;//双向耦合,
averagingModel dense;//dilute;//颗粒很多的情况,流体网格的颗粒速度通过颗粒的系综平均得到,还有dilute用于很稀的情况,一个网格最多只有一个颗粒
clockModel off;//standardClock;//计算程序耗时
smoothingModel off;// localPSizeDiffSmoothing;// constDiffSmoothing; //
forceModels
(
KochHillDrag
gradPForce
viscForce
);
momCoupleModels//动量方程源项耦合方式,这里选择的是隐式
(
implicitCouple
);
//turbulenceModelType "RASProperties";//"LESProperties";//OFversion24x
turbulenceModelType turbulenceProperties;//"LESProperties";//OFversion30x//采用湍流模型的选择,这里与turbulenceProperties保持一致
//===========================================================================//
// sub-model properties
implicitCoupleProps
{
velFieldName "U";
granVelFieldName "Us";
voidfractionFieldName "voidfraction";
}
gradPForceProps
{
pFieldName "p";
voidfractionFieldName "voidfraction";
velocityFieldName "U";
interpolation true;
}
viscForceProps
{
velocityFieldName "U";
interpolation true;
}
KochHillDragProps
{
velFieldName "U";
rhoParticle 1500;
voidfractionFieldName "voidfraction";
interpolation true;
}
engineProps
{
treeSearch true;
}
dividedProps
{
alphaMin 0.1;
scaleUpVol 1.0;//颗粒直径人工变大位置
}
twoWayMPIProps
{
maxNumberOfParticles 10100;
liggghtsPath "../DEM/in.liggghts_resume";
}
// ************************************************************************* //