拉格朗日粒子和动网格不兼容吗？

Jasper 0

在interIsoFoam基础上加入了拉格朗日粒子，编译成功后在网格不变的情况下运行没有问题，但是加入动网格后，计算在进行一段时间后突然终止，没有报错，就是停在了Solving 3-D cloud kinematicCloud这里：

Courant Number mean: 4.05392e-05 max: 0.899987
Interface Courant Number mean: 3.47277e-06 max: 0.47951
deltaT = 4.6765e-06
Time = 0.102143

PIMPLE: iteration 1
Selected 4 cells for refinement out of 208887.
Refined from 208887 to 208915 cells.
Selected 0 split points out of a possible 6252.
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 6.35461e-06, No Iterations 25
GAMGPCG:  Solving for pcorr, Initial residual = 0.0776576, Final residual = 7.22028e-06, No Iterations 15
time step continuity errors : sum local = 4.77276e-15, global = 1.316e-17, cumulative = 8.05896e-11
Evolving kinematicCloud

Solving 3-D cloud kinematicCloud

重新跑了好多次，每次都是停在这里，有大神知道怎么解决吗？

李东岳

你这个最好从代码里面去debug一下，看看那一行hang了

Jasper 0

谢谢岳哥，已经debug过了，发现是KinematicCloud.C里运行到 cloud.motion(cloud, td);这一行就停了，定义在这里：

template<class CloudType>
template<class TrackCloudType>
void Foam::KinematicCloud<CloudType>::motion
(
    TrackCloudType& cloud,
    typename parcelType::trackingData& td
)
{
    td.part() = parcelType::trackingData::tpLinearTrack;
    CloudType::move(cloud, td, solution_.trackTime());

    updateCellOccupancy();
}

有点debug不动了，现在还是毫无头绪。问了chatgpt也没解决。。。

李东岳

CloudType::move(cloud, td, solution_.trackTime());

看这里面哪一行有问题

余正东

我出现过很多次这种情况，你使用的AMI吧，使用NCC动网格就能解决了

Jasper 0

@余正东谢谢大佬，我用的AMI，不过我用的是openfoam.com版本，貌似NCC没有植入。。。

Jasper 0

这玩意好像可以停了后重新接着跑，再停再跑，最终也能到需要计算的时间，但是这样真麻烦啊，啥时候是个头啊唉

Jasper 0

debug根本没法进行，上一次停的时间步这一次直接跑过去了，然后过会又停了，一直等着也太浪费时间了。。。

余正东

@Jasper-0 建议更换org的openfoam11+版本即可，改动没有想象的那么大

CFD中文网

拉格朗日粒子和动网格不兼容吗？