想用双欧拉方法模拟水气两相流,请问twoPhaseEulerFoam中可以加入PBM模型来模拟气泡的聚合和破碎吗?该如何操作呢
Latest posts made by Jasper 0
-
twoPhaseEulerFoam中怎么加入PBM模型
-
RE: 拉格朗日粒子浓度与网格大小的关系
@李东岳 李老师您好,还想请教一下,因为我的粒子直径相对于网格来说较大,所以用DPMFoam计算时发散了,取了几个时间步后的压力分布图发现有一个点压强突然变的很小,如下图:
我的粒子浓度相对来说应该是比较大的
但当我把kinematicCloudProperties里的alphamax从0.9变成0.7之后就可以正常计算了,请问这是什么原因呢?通过调整alphamax来进行计算可以吗?谢谢您 -
RE: 拉格朗日粒子浓度与网格大小的关系
@李东岳 但是我的模拟中粒子是通过kinematicCloudPositions来定义的,所以计算域内粒子数量是固定的。而且我的网格大小是6mm,粒子直径是3mm,对于这种情况是不是网格尺寸对粒子浓度影响较大啊?
-
拉格朗日粒子浓度与网格大小的关系
各位大佬好,
本人最近在想是不是拉格朗日粒子的浓度是与网格大小相关的?因为在openfoam中浓度定义为:theta[celli] += p.nParticle()*p.volume(); theta.primitiveFieldRef() /= mesh_.V();
当网格变小的时候,浓度不就自然而然变大了吗?如果是这样的话,该如何跟实验进行验证呢?
-
RE: MPPICInterFoam 中如何通过alphac得到alpha1和alpha2?
@zhe 在 MPPICInterFoam 中如何通过alphac得到alpha1和alpha2? 中说:
那么alpha1才应该是整个cell里的alpha1啊
没看懂这是什么意思。举个例子,比如alphac=1-theta=0.8, alpha1=0.5,那么phase1所占整个cell的体积分数就是0.8*0.5=0.4,phase2也是0.4,0.4+0.4+0.2=1,在MPPICInterFoam中粒子是考虑了体积分数的
-
kinematicCloud.SU(Uc)的返回值
大家好,我想请教一下调用momentum source函数的时候,返回值是哪一个?在KinematicCloudI.H中有这个函数的定义:
template<class CloudType> inline Foam::tmp<Foam::fvVectorMatrix> Foam::KinematicCloud<CloudType>::SU(volVectorField& U, bool incompressible) const { if (debug) { Pout<< "UTrans min/max = " << min(UTrans()).value() << ", " << max(UTrans()).value() << nl << "UCoeff min/max = " << min(UCoeff()).value() << ", " << max(UCoeff()).value() << endl; } dimensionSet dim(dimForce); if (incompressible) { dim.reset(dimForce/dimDensity); } if (solution_.coupled()) { if (solution_.semiImplicit("U")) { volScalarField::Internal Vdt(mesh_.V()*this->db().time().deltaT()); if (incompressible) { Vdt.dimensions() *= dimDensity; } return UTrans()/Vdt - fvm::Sp(UCoeff()/Vdt, U) + UCoeff()/Vdt*U; } else { tmp<fvVectorMatrix> tfvm(new fvVectorMatrix(U, dim)); fvVectorMatrix& fvm = tfvm.ref(); fvm.source() = -UTrans()/(this->db().time().deltaT()); return tfvm; } } return tmp<fvVectorMatrix>::New(U, dim); }
在调用这个函数的时候,如果我用semiImplicit格式的话,返回的是
return UTrans()/Vdt - fvm::Sp(UCoeff()/Vdt, U) + UCoeff()/Vdt*U;
还是
return tmp<fvVectorMatrix>::New(U, dim);
呢?
我认为返回的应该是后者吧 tmp<fvVectorMatrix>::New(U, dim);
我想知道它是什么意思呢?