twoPhaseEulerFoam中的IATE模型系数
-
各位前辈,经过验证发现,对于IATE模型中表示气泡聚合和破碎的Crc和Cti并未参与计算,因为这两个函数的返回值是tR,而进行赋值计算的是R,经过连续输出发现,赋值计算的R无论是否为0,返回值tR始终为0,因此对于Crc、Cti无论改变为多大的值对于气泡的聚合破碎毫无作用。
在将返回值改为R以后,由于给定算例的连续相速度始终为0,而Crc和Cti计算R表达式中用到了连续相速度,所以还会出现返回值一直为0的情况,若改变算例,连续相速度给一定值后,就会发现此时返回值是tR时,返回值还一直为0,但返回值是R时,发现此时返回值就发生了变化,粒径分布也就发生了变化。因此在想是不是这个程序返回值tR是写错了还是我有哪个地方搞错了,请各位老师指点一下。另外还遇到一个问题就是在linux系统下,将返回值改正为R以后会报错如下:(以上所说的测试是基于Windows版本的绿化版本),希望各位前辈给点建议,这个错误应该怎么处理,不胜感激!--> FOAM FATAL ERROR: different patches for fvPatchField<Type>s From function void Foam::fvPatchField<Type>::check(const Foam::fvPatchField<Type>&) const [with Type = double] in file /opt/openfoam4/src/finiteVolume/lnInclude/fvPatchField.C at line 209. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::operator-=(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:? #3 Foam::diameterModels::IATE::correct() at ??:? #4 Foam::twoPhaseSystem::correct() at ??:? #5 ? at ??:? #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #7 ? at ??:? 已放弃 (核心已转储)
-
-
@李东岳 这个地方我表述错了,但是在算例中即使把Crc、Cti调整到100000000这个量级对计算结果还是毫无影响,即使存在临界数的判断,Crc这个依靠相分数判断的碰撞聚集项肯定会发生变化。还有就是对R赋值,返回却是tR,这个地方对吗,我现在是这个地方不理解,谢谢东岳老师的解答。
-
把这一行
volScalarField R = tR();
改成
volScalarField& R = tR.ref();
再试试
-
怎么样 有进展么
-
东岳老师不好意思,我现在试试,刚看到您的回复贴
-
@李东岳 东岳老师我刚刚试了一下,还是不行,会报同样的错误
-
我测试了一下,可以编译
-
@李东岳 对的,可以编译,但是运行算例还是会报错,和我以前一样,编译通过,运行出错
-
4楼那个错误么,我测试了下运行也没有错误,你是不是动别的东西了
-
那老师方便把您改动好的代码发我一下吗,我对着检查再看看,谢谢老师
-
把这一行
volScalarField R = tR();
改成
volScalarField& R = tR.ref();
就这个 没别的了
-
那最后一行还是return tR吗?
-
-
这种情况,我一会再测试一下Crc和Cti是否参与了计算,谢谢老师百忙之中抽空解答我的疑惑,谢谢
-
有进展么
-
@李东岳 目前是这样的,按照东岳老师您的意见更改完以后,是可行的,Crc和Cti均可以正常参与计算了。但是计算结果和在windows下仅仅返回R的计算结果是不一样的,不知道其中原因。目前正在使用老师您的结果在跑,观察计算结果,十分感谢东岳老师的帮助
第二种是仅仅把return tR; 改成 return R;
-
结果怎么样
-
@李东岳 已经完全正常了,多谢东岳老师一路指点,谢谢。首先能获得结果,然后再获得精准的结果
-
好。
IATE这个目前研究的少,你可以试试。也可以用其他模型对比下。很多模型可以选,比如矩方法、section方法、n传输方法、OPOSPM等
-
好的东岳老师,我再琢磨看看,感谢感谢
-
结果怎样了,有什么进展
-
感谢李老师还一直挂念,目前还在用ITAE计算一些工况水一水毕业论文,计算结果是正常的。您上面说的方法在双欧拉里面可以选嘛,怎么使用还请李老师指点一下
-
完全没问题,目前我对OpenFOAM双欧拉模型中植入DQST算法非常感兴趣。不过自己一直没时间做。如果你读博士的话,或许可以搞一搞。
-
@李东岳 谢谢东岳老师,如果确定读博再深入学习看看能不能尝试下,最后还是很感谢李老师一直的指导,祝愿中文网和OKS课程红红火火