拉格朗日一个bug

李东岳

上述代码为OpenFOAM-5.x下的结果，经测试OpenFOAM-4.x结果正确。拉格朗日这面玩的这么少么？不仅是Cloud<basicKinematicParcel>，Cloud<basicKinematicCollidingParcel>也有这个问题。

4.x:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       Cloud<basicKinematicCollidingParcel>;
    location    "5/lagrangian/kinematicCloud";
    object      positions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

850
(
(0.1652375049 0.5 0.5) 0
(0.3082808448 0.5 0.5) 0
(0.09319717789 0.5 0.5) 0
(0.1760049582 0.5 0.5) 0
(0.1925717587 0.5 0.5) 0
(0.3543820922 0.5 0.5) 0

5.x:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       Cloud<basicKinematicCollidingParcel>;
    location    "5/lagrangian/kinematicCloud";
    object      positions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

850
(
(0.3261957614 0.3369021193 0.3369021193 0) 0 9 2//什么东西这是？。。
(0.4001988533 0.2999005734 0.2999005734 0) 0 9 2
(0.1530216921 0.4234891539 0.4234891539 0) 0 9 2
(0.3338597466 0.3330701267 0.3330701267 0) 0 9 2
(0.3488235028 0.3255882486 0.3255882486 0) 0 9 2
(0.4730656805 0.2634671597 0.2634671597 0) 0 9 2
(0.2571585867 0.3714207066 0.3714207066 0) 0 9 2
(0.4849837396 0.2575081302 0.2575081302 0) 0 9 2
(0.1467023563 0.4266488218 0.4266488218 0) 0 9 2
(0.2706319485 0.3646840257 0.3646840257 0) 0 9 2
(0.3640142443 0.3179928778 0.3179928778 0) 0 9 2
(0.2279260728 0.3860369636 0.3860369636 0) 0 9 2

yuan

这个跟颗粒追踪算法的更新有关吧，在 Openfoam 5.0 的发布说明里有提到（OpenFOAM Barycentric Tracking），更新之后颗粒输出的格式与之前版本不兼容。

在 $FOAM_SRC/lagrangian/basic/particle/particleIO.C 中可以看到保存颗粒位置信息时的区别（writePosition 函数）：

之前版本中为

os  << position_ << token::SPACE << celli_;

在 5.0 版中为

os  << coordinates_
     << token::SPACE << celli_
     << token::SPACE << tetFacei_
     << token::SPACE << tetPti_;

变量的定义在 $FOAM_SRC/lagrangian/basic/particle/particle.H 中可以找到。

李东岳

有道理！确实输出的是

os  << coordinates_
     << token::SPACE << celli_
     << token::SPACE << tetFacei_
     << token::SPACE << tetPti_;

5.x下面怎么后处理呢？paraview提示

ERROR: In /home/buildslave/dashboards/buildbot/paraview-pvbinsdash-linux-shared-release_superbuild/build/superbuild/paraview/src/VTK/IO/Geometry/vtkOpenFOAMReader.cxx, line 8285
vtkOpenFOAMReaderPrivate (0x5e89d60): Error reading line 20 of /home/dyfluid/OpenFOAM_tutorials/5.x/lagrangian/DPMFoam/Goldschmidt/0.01/lagrangian/kinematicCloud/positions: Expected punctuation token ')', found 0

我看5.0把position位置矢量写到obj文件里面去了

李东岳

um, 看来得学习学习barycentric坐标这个玩意

yuan

@李东岳我这里后处理可以用 paraFoam，也可以先运行 foamToVTK 后用 paraview。没见过这个错误，可能是因为我编译 OpenFOAM 时顺便编译了 paraview ？

李东岳

嗯，有可能，我没编译paraFoam，用的是paraview。

马乔

5.X用了粒子所在剖分四面体内的权重作为局部坐标，这样省去了很多计算量，但是对后处理不友好了，因为之前一直都是输出的全局坐标。

浪迹天大

请问有好的方式处理这个问题吗？

Yufeng Wei

@浪迹天大今天遇到了这个问题，解决方法就是楼上说的，用paraFoam 打开

yfclark

@Yufeng-Wei 要么改用基金会版本，基金会版本专门做了优化，自动把质心坐标转化为三维坐标，要么试一试这个https://github.com/blueCFD/lagrangianExtraFunctionObjects

yuhx

我也遇到了这个问题，查了一圈~发现FoamToVTK比较好使，可以输出保存颗粒数据，但是一定注意scv文件不要用excel打开，否则会自动删减数据！（这是什么神奇的操作不清楚）

myjimp

@yuhx 你好，转化为VTK显示不出粒子,是什么情况？

JqBu

用blueCFD转化了粒子坐标后（https://github.com/blueCFD/lagrangianExtraFunctionObjects/tree/OF5x）突然遇到了另一个报错看起来像是转化过程中浮点数出错了后面找了代码发现是调用了regIOobject::write(cloud_.size())这个函数，然而我只在of里找到了 virtual bool write(const bool write = true) const这个函数看起来不太对的样子

请问有人遇到过吗请问如何解决？

疏影横斜水清浅

在OpenFOAM5.X之后的版本中都改了这个position文件的记录方式，那如果想得到粒子的位置信息从而得到轨迹，是不是需要使用p.position函数得到坐标矢量，存到新的矢量变量中，这样就可以得到粒子坐标轨迹了

CFD中文网

拉格朗日一个bug