@星星星星晴 非常感谢友友的建议！

你这个就是标准的类似6dof的动网格原理。你这个力要推动网格点的移动。你要么看一下SDA、或者6dof那面的处理方法。openfoam这面网格点的运动通过septernion进行。不是传统的力来移动3D坐标。要么你就自己写一个网格边界点受力被移动的在笛卡尔坐标系下的xyz函数。

@吴建民 我前两天看LIGGGHTS已经被替代了吧。我个人经验是 OF5 - OF8 在拉格朗日部分计算速度有很大的提升。

千万级别的颗粒是什么意思，是你要追踪千万级别？还是说一共千万个颗粒，在domain中只有10w个？ 或者你考虑一下 一个parcel代表更多的particle。

如大佬所说，你如果考虑碰撞，什么碰撞模型？deterministic的话根本不可能。

@李东岳 确实没看，主要在网上搜来着🤦‍，感谢老师！我现在研究一下！

@李东岳 嗯嗯，好的，明白了，谢谢老师

动网格一般发散的原因就是网格变形太大了导致网格尺寸太小了

@李东岳 感谢老师的回复。foamListtimes -help 里好像没有cp的选项，但在Allrun里加上两行

{latestTime="$(foamListTimes -latestTime)" cp -r "${latestTime}" newTime}

可以正常运行了。

您好，对于粒子轨迹现在您有什么解决办法吗

@AppleKiller 您好，请问您做气泡上浮使用什么求解器

@杨英狄 老铁这个研究怎么样了。发表没有。

请问还有和这个类似的求解器吗，是可以用于可压缩的，这好像是不可压的

时间步长缩小。这个是不好整。主要网格太大。调试起来都费劲。

可以参考这个写一个 https://github.com/OpenFOAM/OpenFOAM-dev/blob/master/src/meshTools/sets/cellSources/cylinderToCell/cylinderToCell.C

:146: :146: :146:

好的，李老师，谢谢

github上搜一下cfdem就能找到

是的。

@chapaofa 大佬辛苦啦，我跑就是这样，不知道为啥，感觉命令没问题，我换台电脑试试吧。谢谢啦!

@李东岳 谢谢李老师，网格检查过了，没有问题，把cyclicAMI的边界设置为symmetry是可以计算的。刚才按照您的建议试过了，也还是报相同的错误，错误信息如下：

/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2112 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : _6e1fca0e-20220610 OPENFOAM=2112 patch=220610 version=2112 Arch : "LSB;label=32;scalar=64" Exec : dbnsFoam Date : Jun 06 2024 Time : 18:13:40 Host : DESKTOP-H49IU2O PID : 137843 I/O : uncollated Case : /home/duan/JIA_Xiaotong/ICSFoam-master/tutorials/test_5 nProcs : 1 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Steady-state analysis detected pseudoTime: tolerance = 1e-10, relTol = 1e-10 Local timestepping selected Reading thermophysical properties Selecting thermodynamics package { type hePsiThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } AMI: Creating addressing and weights between 18620 source faces and 18620 target faces AMI: Patch source sum(weights) min:0 max:1.84903 average:0.999732 AMI: Patch target sum(weights) min:0.030709 max:1.15763 average:0.994612 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in /lib/x86_64-linux-gnu/libc.so.6 #3 ? at ??:? #4 ? at ??:? #5 Foam::psiThermo::addfvMeshConstructorToTable<Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:? #6 ? at ??:? #7 Foam::psiThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:? #8 ? in ~/OpenFOAM/duan-v2112/platforms/linux64GccDPInt32Opt/bin/dbnsFoam #9 ? in /lib/x86_64-linux-gnu/libc.so.6 #10 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #11 ? in ~/OpenFOAM/duan-v2112/platforms/linux64GccDPInt32Opt/bin/dbnsFoam Floating point exception

老师，有没有可能是因为叶栅模型的旋转半径太大（120m)导致的呢？

@李东岳 非常感谢！

不好意思，下面的问题也许会有点啰嗦。

初始条件 0 文件中的 alpha.phase1 默认是指定主相的边界条件吧，怎么特别指定次相的alpha进出口条件呢？

此外请教一下，当我给出口压力设置定值，出口速度给 codeStream 写的 parabolic 速度分布，入口速度设置为固定值或者 zeroGradient，都会导致计算很快出现 Floating point exception 。有没有什么方法可以特别指定此相的流出速度分布吗？

非常感谢大大的时间和耐心回复。