CFD中文网

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

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

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

谢谢岳哥，已经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也没解决。。。

在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

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

@coolhhh 发现问题了，因为我的是两相流，需要设定field 为U.water

我在controldict里面加入如下代码想计算vorticity，结果发现只在初始时刻算了一次，后面时刻的值就不变了

vorticity1
{
    type            vorticity;
    libs (fieldFunctionObjects);
    writeControl writeTime;
}

即使我加入执行控制，vorticity依然是不随着计算时间变化，有没有大佬知道这是咋回事？

vorticity1
{
    type            vorticity;
    libs (fieldFunctionObjects);
    executeControl  timeStep;
    executeInterval 1;
    writeControl writeTime;
	
}

@李东岳 加了，就是把杨氏模量调大点就会有堆积了，估计之前的太小了导致没有弹性。。。

谢大佬回复，刚搜了一下貌似跟材料和粒径都有关系

我发现这个youngsModulus的系数貌似跟粒径有关。

1. 当我用粒径为3毫米的时候，我设置系数为1e4粒子并不会发生融合；
2. 但当我用粒径1cm的时候，系数为1e4就会融合了，然后把系数调到1e8就不会融合了。
   但是这个系数越大计算速度越慢。请问这个youngsModulus系数有没有一个设置的标准啊？还是说只能像我现在这样试？

@李东岳 终于找到了问题所在，不知道为什么在这个case里我把young's Modulus设置成了1e-8，改成1e8之后就会有堆积了

谢谢李老师，我试一下然后再来说一下进展

放小了是会好很多，但是我这个scale factor是设置的1呀，这样显示出来的不应该就是实际的粒子大小吗？我设置了粒子直径是1cm，这个计算域的直径也才15cm，一共有2万个粒子，理论上来说应该肯定堆积起来了吧？学生对这个沉积的问题确实是困扰好久了，还望李老师帮忙解答一下，万分感谢啊！

@李东岳 我也想到了这个问题，于是把颗粒直径增大，结果还是一样，貌似这些颗粒都挤在一块了，甚至直接融合在一块了。。。下图是我最新跑的用vector滤镜处理的（中间空出来那一块本来是个injector，在这个测试里我设置成了壁面）

各位大佬好，我用DPMFoam简单算了一个case，颗粒在静止水体中释放，如下图

等颗粒全部沉降后，最后却只沉积在了底面壁面上，如何才能让颗粒像实际一样堆积起来呢？

您好请问这个问题解决了吗？我也遇到了

@李东岳 谢谢李老师答复。我现在就在用-oversubscribe方法跑，但是我不知道这个方法具体是怎么调用的核数，比如是不是我的两个cpu都调用了，还是说只调用了一个cpu里的38个线程（一个cpu总共48线程）

请教各位大佬，我的工作站为双cpu：Intel(R) dual-Xeon(R) Gold 5220R ，每个cpu有24 个cores，总数为48cores，96线程。我在计算时划分了38个processors，但是运行如下命令时出错：mpirun -np 38 interFoam -parallel
错误代码为

There are not enough slots available in the system to satisfy the 38
slots that were requested by the application:

貌似只识别了一个cpu，怎么才能在这两个cpu上进行并行计算呢？

@jasper-0 更新一下目前进度，我尝试在KinematicLookupTableInjection.C里将injectors_设置为MUST_READ_IF_MODIFIED,

template<class CloudType>
Foam::KinematicLookupTableInjection<CloudType>::KinematicLookupTableInjection
(
    const dictionary& dict,
    CloudType& owner,
    const word& modelName
)
:
    InjectionModel<CloudType>(dict, owner, modelName, typeName),
    inputFileName_(this->coeffDict().lookup("inputFile")),
    duration_(this->coeffDict().getScalar("duration")),
    parcelsPerSecond_(this->coeffDict().getScalar("parcelsPerSecond")),
    randomise_(this->coeffDict().getBool("randomise")),
    injectors_
    (
        IOobject
        (
            inputFileName_,
            owner.db().time().constant(),
            owner.db(),
            IOobject::MUST_READ_IF_MODIFIED,
            IOobject::NO_WRITE
        )
    ),

但是在运行solver的时候提示不支持自动重新读取，有没有大佬知道这种情况咋解决呢？

GlobalIOList<kinematicParcelInjectionData> kinematicLookupTableInjection constructed with IOobject::MUST_READ_IF_MODIFIED but GlobalIOList<kinematicParcelInjectionData> does not support automatic rereading.

@oitocfd 谢谢大佬回复，这两天忙着做实验没看到消息。所以我如果只考虑Drag force，那么这个phicForces就不用加了吧？