@李东岳 在 分享CFDEM+OpenFOAM+LIGGGHTS初次编译安装过程 中说:
ERROR: Invalid dump style (/home/dyfluid/CFDEM/LIGGGHTS/src/output.cpp:589)
代码的最后一行也是同样的问题,可能时安装的时候有点问题,我再找找原因
@李东岳 在 分享CFDEM+OpenFOAM+LIGGGHTS初次编译安装过程 中说:
ERROR: Invalid dump style (/home/dyfluid/CFDEM/LIGGGHTS/src/output.cpp:589)
代码的最后一行也是同样的问题,可能时安装的时候有点问题,我再找找原因
@李东岳 李老师,在用您提供的虚拟机跑算例时碰到如下问题:
dump dmp all custom/vtk 5000 post/liggghts_init_*.vtk id type x y z vx vy vz fx fy fz omegax omegay omegaz radius
ERROR: Invalid dump style (/home/dyfluid/CFDEM/LIGGGHTS/src/output.cpp:589)
$ \nabla\cdot(Ue)=U\cdot\nabla e + e \nabla \cdot U$ ,$\nabla \cdot U \neq 0$
关于这一项$ \nabla\cdot(Ue)$在openfoam中是直接写成div(phi,e)呢?还是div(phi,e)+ediv(U)呢?
目前测试出来差别比较大
@李东岳 从另外一个思路来做,能不能指定内场中一区域的速度分布满足某一个函数关系式呢?
@yhdthu 前辈能够分享一下解决的经验吗?
二维计算域示意图如上图所示,请问大家在openFoam中能否指定计算域中的某一块区域的涡量为某一值吗?
@fubianhanshu 需要调用particleTrack函数哟
有没有大佬可以提点意见呀!
背景:在使用欧拉-拉格朗日算法求解固液两相流动时,使用颗粒追踪的后处理函数,如图所示:
在结果中将会输出颗粒的属性:
现在我希望增加输出项,简单的如颗粒的体积,我修改了kinematicParcelIO.C中的writefields函数:添加了如下代码:
但是计算结果中还是没有颗粒体积的输出?
请问大家有碰到过类似的问题吗?
谢谢各位了!
@bubbleEL 谢谢你
计算区域为旋转区域与非旋转区域,颗粒由非旋转区域入射,但是后期处理颗粒轨迹时候发现,颗粒无法通过动静交界面(cyclicAMI).
在颗粒求解文件中设置交界面为none
请问大家有碰到过这问题吗?有什么好的解决方法吗?
@bubbleEL 我目前也碰到了颗粒无法进入旋转区域的问题,我的交界面设置的也时none,请问楼主是进行了怎样的修改,使得颗粒能够进入旋转区域呢?
请问楼主这个问题解决了吗?
@东岳
我的网格有多个交界面(用的是cyclicAMI),交界面对于拉格朗日的计算会有影响吗?
目前我在采用稳态的欧拉拉格朗日求解器,求解水轮机的磨损问题,颗粒数为50000,网格数为1100w左右。颗粒文件的设置如下:```
code_text
solution
{
active true;
calcFrequency 15;
maxTrackTime 6;
maxCo 0.6;
coupled true;
cellValueSourceCorrection off; //suggested to be off
interpolationSchemes
{
Urel cellPoint; //U for water
rho cell; //rho of water
mu cell; //mu of water
DUcDt cell; //DUcDt of water
}
integrationSchemes
{
U Euler;
}
sourceTerms
{
resetOnStartup false;
schemes
{
U semiImplicit 1;
}
}
}
constantProperties
{
rho0 2650; //the density for the particle
}
subModels
{
particleForces
{
sphereDrag; //drag forces
gravity; // gravity force
pressureGradient //pressure gradient force
{
U Urel;
}
virtualMass //virtual mass force
{
Cvm 0.5;
}
MRF
{
omega (0 0 -104.72);
centreOfRotation (0 0 0);
Ctemp MRFC;
}
}
injectionModels
{
model1
{
type patchInjection;
parcelBasisType fixed;
patchName INLET_SV;
U0 (0 0 0);
nParticle 420;
duration 1;
parcelsPerSecond 50000;
flowRateProfile constant 1;
massFlowRate 0;
//positionsFile "kinematicCloudPositions";
sizeDistribution
{
type fixedValue;
fixedValueDistribution
{
value 0.0002;//m
}
}
}
}
dispersionModel stochasticDispersionRAS;
patchInteractionModel UserLocalInteraction;
//using the tabakoff rebound model
UserLocalInteractionCoeffs
{
patches
(
"(.*)"
{
type rebound;
reboundModel 2;
}
OUTLET_SC
{
type none;
}
INLET_SV
{
type escape;
}
OUTLET_SV
{
type none;
}
INLET_GV
{
type none;
}
OUTLET_GV
{
type none;
}
INLET_RUN
{
type none;
}
OUTLET_RUN
{
type none;
}
INLET_DT
{
type none;
}
INLET_SC
{
type none;
}
OUTLET_DT
{
type escape;
}
);
}
heatTransferModel none;
surfaceFilmModel none;
collisionModel none;
stochasticCollisionModel none;
radiation off;
}
流场迭代15步之后,会进行颗粒的迭代,颗粒计算平均耗时5min。
但每次计算至210步时,颗粒就算不动了(前面的颗粒计算均是正常的)。
我如果把颗粒数降至10000,计算至750步,颗粒就算不动了。
请问一下大家又碰到这类问题吗?有什么比较好的解决办法吗?
@xpqiu 您能够介绍一下上述方法,单看您发的文件不太懂
@xpqiu 手动划分网格吗?