稍稍看了眼Palabos算例cylinder3d | 动量交换计算drag/lift

这个算例是最近更新的palabos版本才有的算例，它多出了升阻力（英文是drag和lift，如果我没翻译错的话）的计算。

后续是定义了foreceIds这个三维数组的，后续是要求和的，这里subscribe意味登记一下id，简言之就是预定好了这些lattice以后要求和。

lattice.toggleInternalStatistics(true);
    forceIds[0] = lattice.internalStatSubscription().subscribeSum();
    forceIds[1] = lattice.internalStatSubscription().subscribeSum();
    forceIds[2] = lattice.internalStatSubscription().subscribeSum();

如下定义这个圆柱区域时，不是单纯的使用bounceback，而是MomentumExchangeBounceBack，其中有个括号(forceIds)，动量交换的反弹，定义完之后还对整个lattice区域初始化动量交换。

defineDynamics(lattice, lattice.getBoundingBox(),
                   new CylinderShapeDomain3D<T>(cx, cz, radius),
                   new plb::MomentumExchangeBounceBack<T,DESCRIPTOR>(forceIds));
    initializeMomentumExchange (lattice, lattice.getBoundingBox() );

这里就使用getSum进行了求和，forceIds就是之前登记过的区域。

T drag = lattice.getInternalStatistics().getSum(forceIds[0]);
T lift = lattice.getInternalStatistics().getSum(forceIds[2]);

后面还提供了系数。

T drag_coef_factor = 1.0 / (length * util::sqr(vel) * diameter * avgRho);