@youhaoyu 看了本比较通用的书,里面把我第一张图里(1.15)直接写在方程里面的(我感觉这样写比较合适)。那么这个R0和Pg0我还是不理解什么意思,有没有大佬解释解释
youhaoyu
帖子
-
对于RP气泡长大方程疑问 -
对于RP气泡长大方程疑问最近再学习RP方程时候发现方程解释有点不理解。
(1.12)是RP方程,其中(1.13)表示气泡压力等于饱和蒸汽压加上后面一个量,就是后面的Pg0和R0我看大部分的解释都是某一个时间段的平衡的气泡压力和直径,那我理解来看,这两个一个是定值,是一个最开始我们就可以赋予的一个固定值。但是后面(1.14)和(1.15)表示了,随着气泡长大,Pg0和R0在变化????所以这个Pg和R0到底是什么??
第二个问题就是:我是否能够简化一下使用理想气体方程,直接把PBR=Pg0R0? -
Euler-Lagrange的一些解析,sprayFoam@星星星星晴 你好,最近想向DPMFoam添加一些气泡长大方程,比如RP方程,如何去添加啊,有没有什么思路推荐啊,因为本人实在是太不会编程这些东西了。
最近看了篇博士论文,《水质对梢涡空化初生的影响研究》,上面有一些源文件不知道能不能用,意思是在运动包裹源文件中,添加一些颗粒体积更新代码,再植入自己方程什么的。这个思路能行吗? -
拉格朗日求解器合并块问题@Tens 都试了不行,减少并行块数量有点用,但是偶尔也犯病,这玩意我现在看玄学了
-
气泡破裂模型@AppleKiller 感谢,我再看看,以前也接触过这些资料,感觉用欧拉法求解器为基础修改的拉格朗日求解器不大好,没怎么深入去学习。我再看看这些思路,感谢
-
气泡破裂模型@pengdi 十分感谢,我更多的是对于添加代码这个编程逻辑问题一脸茫然,这个如何去学习呢?我目前大致知道可能是去修改一下src里面的一些文件,然后去添加一些方程需要的一些参数,比如方程有的表面张力等等,但是如何去实践在OpenFOAM代码中,我目前不知道怎么去学习,请问您有什么推荐吗
-
气泡破裂模型@AppleKiller 大佬,请问如何去添加RP方程呢,是直接修改源文件的粒子云库吗?请问大佬有什么思路分享么,我太小白了,对于编程这些不大懂
-
气泡破裂模型@pengdi 大佬请问一下,如果在拉格朗日求解器DPMFoam里面加RP方程,应该是怎么一个思路,我太小白了,完全不知道怎么去操作
-
openfoam绘制复杂楔形网格报错@风云5091 什么叫做把轴定义成一个面?我不了解你目前这个方向可能不理解你的这个问题。你可以配点你的模型图我看看
-
openfoam绘制复杂楔形网格报错这种模型用icem挺好用的,直接icem建模就行了,openfoam自带的建模和网格划分我感觉大部分人都不用吧哈哈哈哈,感觉没有界面显示太反人类了(我瞎猜的)
-
拉格朗日中的湍流分散力模型@李东岳 东岳老师您好,我最近发现了一个现象可能造成了我上述提出的问题哈哈哈哈哈,我最开始使用的是DPMFoam默认的的Ergunwenyuforce曳力,我后来换成了sphereforce曳力就能稳定了,而且对于颗粒运动和流场有挺大影响,看来这和我的设置也有一些关系。(最开始图方便就直接使用的默认的模板了,看来还是要根据实际情况来改变)
-
拉格朗日中的湍流分散力模型@李东岳 感谢老师提供思路,我的dispersion模型是默认的none,明天试一下您的思路再去模拟一下看看。这与入射口面积会不会有关?因为我现在这个是简化的模型,真正的模型颗粒入口面积比我现在小得多
-
拉格朗日中的湍流分散力模型@李东岳 东岳老师您好,最近在使用DPMFoam求解器时候发现我的粒子不能稳定,始终是这种摆动状态,与欧拉法算出来的有些区别,以下是我的模型. 造成这种现象的原因是不是和随机游走模型有关?如果我想取消随机游走模型,使得颗粒成稳定的轨迹向上浮动,如何去取消这个模型呢?
-
DPMFoam后处理您好,对于粒子轨迹现在您有什么解决办法吗
-
openFOAM中颗粒压力梯度力于虚拟质量力的计算问题@AppleKiller 您好,请问您做气泡上浮使用什么求解器
-
DPMFoam求解器添加粒子生长程序您好,请问能学习一下您这个的代码么
-
icem非结构网格如何划分block结构化网格非结构化网格按照网格是六面体还是体面体区分好像是不完全正确的。我看一个老师讲解说,结构化和非结构化区分是文件保存格式的区分,结构化是使用索引,非结构化是把每个网格都区分出来了。我看各位大佬上面回答可能是直接讨论的四面体网格,因为我最近也在学openfoam和icem,icem可以通过自动化分成四面体非结构网格,也可以分块化成6面体网格,但是我通常画完6面体网格后,icem有可以转化成非结构网格按钮,将画好的网格转化成非结构体网格(转化后依然是6面体)。
所以如果你要话非结构化网格也可以分块(画了6面体转化一下),如果你说的是画四面体非结构化网格分块的话,这方面知识我就不知道了 -
拉格朗日求解器合并块问题拉格朗日DPMFoam求解器使用过程中,我并行结束后合并不了块,出现以下代码
Create time Reconstructing fields for mesh region0 Time = 0.5 Reconstructing FV fields Reconstructing volScalarField::DimensionedInternalFields kinematicCloud:UCoeff Reconstructing volVectorField::DimensionedInternalFields kinematicCloud:UTrans Reconstructing volScalarFields k.air mu.air alpha.air epsilon.air nut.air p rho.air Reconstructing volVectorFields U.air Reconstructing surfaceScalarFields phi.air Reconstructing point fields No point fields Reconstructing lagrangian fields for cloud kinematicCloud --> FOAM FATAL ERROR: cell, tetFace and tetPt search failure at position (0.855908 3.37988 -0.0348682) for requested cell 78838 If this is a restart or reconstruction/decomposition etc. it is likely that the write precision is not sufficient. Either increase 'writePrecision' or set 'writeFormat' to 'binary' From function void Foam::particle::initCellFacePt() in file /home/zly/OpenFOAM/OpenFOAM-3.0.0/src/lagrangian/basic/lnInclude/particleI.H at line 735. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam::reconstructLagrangianPositions(Foam::polyMesh const&, Foam::word const&, Foam::PtrList<Foam::fvMesh> const&, Foam::PtrList<Foam::IOList<int> > const&, Foam::PtrList<Foam::IOList<int> > const&) at ??:? #3 ? at ??:? #4 __libc_start_main in "/lib64/libc.so.6" #5 ? at ??:? Aborted (core dumped
然后当我想mpirun -np继续运行时候,他就运行不了了,出现以下报错
Create time Create mesh for time = 7 PIMPLE: Operating solver in PISO mode Reading g Reading field U Reading field p Reading/calculating continuous-phase face flux field phic Creating turbulence model Selecting incompressible transport model Newtonian Creating field alphac Constructing kinematicCloud kinematicCloud Constructing particle forces Selecting particle force ErgunWenYuDrag Selecting particle force gravity Constructing cloud functions none Constructing particle injection models Creating injector: model1 Selecting injection model patchInjection Constructing 3-D injection Choosing nParticle to be a fixed value, massTotal variable now does not determine anything. [6] [6] [6] --> FOAM FATAL ERROR: [6] cell, tetFace and tetPt search failure at position (0.856587 3.37816 -0.0300663) for requested cell 78838 If this is a restart or reconstruction/decomposition etc. it is likely that the write precision is not sufficient. Either increase 'writePrecision' or set 'writeFormat' to 'binary' [6] [6] From function void Foam::particle::initCellFacePt() [6] in file /home/zly/OpenFOAM/OpenFOAM-3.0.0/src/lagrangian/basic/lnInclude/particleI.H at line 735. [6] FOAM parallel run aborting [6] [6] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [6] #1 Foam::error::abort() at ??:? [6] #2 ? at ??:? [6] #3 ? at ??:? [6] #4 ? at ??:? [6] #5 ? at ??:? [6] #6 ? at ??:? [6] #7 __libc_start_main in "/lib64/libc.so.6" [6] #8 ? at ??:? -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 6 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. -------------------------------------------------------------------------- ^Cmpirun: abort is already in progress...hit ctrl-c again to forcibly terminate
我认为主要报错原因是这一段话
cell, tetFace and tetPt search failure at position (0.856587 3.37816 -0.0300663) for requested cell 78838 If this is a restart or reconstruction/decomposition etc. it is likely that the write precision is not sufficient. Either increase 'writePrecision' or set 'writeFormat' to 'binary
但是我不理解他意思,是指的这个时间段颗粒在分块的交界线上,所以产生报错现象吗???我感觉of不会出现这种问题把。那出现这种问题有什么补办法么,难道只有修改decomposePar个数,或者把网格画粗一点重新计算么?有没有什么方法能够让他合并成功?
主要为我问题一个是我合并不了块,无法看拉格朗日项,第二个问题是我无法继续进行计算(因为我中途暂停了只是想看一下流场是什么样子的。我能合并流场的块,不能合并拉格朗日项的块很难受)。 -
cell,tecFace and tetPt search failure at position@yuhx 请问您现在解决了么
-
关于MPPICInterFoam@sw001 ssh就是实验室机架式服务器,然后用SSH连接服务器一个软件罢了哈哈哈哈。你这个是电脑自己本身就是Linux系统运行openfoam吗?然后编译好的paraview打开吗?tecplot好像确实无法看openfoam的拉格朗日相。至于这个怎么看相分数,目前你就这样tecplot和paraview混着用吧哈哈哈哈,这些后处理我也不大懂,或者你学学VTK格式,再用paraview打开?