CFD中文网

需求是 通过probes，随计算时间步实时输出 目标位置的压力值。

现在算一个船体运动响应的动网格算例，用probes监测某几处位置（均位于船体外表面）的压力p。
计划算50s，算到22s的时候计算中止。在22s重新提交，发现由于船体运动，0时刻probes的位置在22s时刻 不在流场域了，所以probes的位置都捕捉不到了。。。
这种情况有什么解决方案吗？
求指教。

@李东岳 解决了，谢谢！

网格划分好了之后，计算报错。貌似是wall边界条件的问题？还是湍流模型的问题呢？
哪位大佬可以点拨一下，不胜感激！

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
[8] 
[8] 
[8] --> FOAM FATAL IO ERROR: 
[8] Entry 'method' not found in dictionary "IOstream.wallDist"
[8] 
[8] 
[8] file: IOstream.wallDist
[8] 
[8]     From function bool Foam::dictionary::readEntry(const Foam::word&, T&, Foam::keyType::option, bool) const [with T = Foam::word]
[8]     in file /home/b408/OpenFOAM/OpenFOAM-v1912/src/OpenFOAM/lnInclude/dictionaryTemplates.C at line 334.
[8] 
FOAM parallel run exiting

@李东岳 好的，我试试这个思路，多谢李老师了！

@李东岳
李老师，我试了下，controlDict文件里是这样写的：

functions
{
    {                                                        // Line 57
        const point A = point[222];
        vector = A[222] - A[222].oldtime();
        Info << vector;
    }
........
}

运行报错：

Create time

[0] 
[0] 
[0] --> FOAM FATAL IO ERROR: 
[0] Found { but expected }
[0] 
[0] file: /home/b408/OpenFOAM/b408-v1912/run_test/DTCHullMoving_t1222_0/system/controlDict at line 57.
[0] 
[0]     From function static bool Foam::entry::New(Foam::dictionary&, Foam::Istream&, Foam::entry::inputMode, int)
[0]     in file db/dictionary/entry/entryIO.C at line 179.
[0] 
FOAM parallel run exiting
[0] 
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD 
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------

@李东岳 这是可以的。我想实现的是 某个位置的位移随每个计算时间步timeStep实时输出。就像controlDict里的function通过probe输出p和U一样。

求助 请教，在包含结构物运动的动网格算例中（例如floatingObject，DTCHullMoving此类），怎样提取结构物表面某一点的位移pointDisplacement？
在controlDict文件中 用Probe尝试过，只能提取U p等，无法提取pointDisplacement。个人猜测原因是：U p是数据记录在网格的场，而pointDisplacement数据记录在节点。
在网上也看了好多关于数据分析和后处理的资料，一直找不到方法。求大神指教！不胜感激！

@五好青年 请问您知道提取位移的方法了吗？遇到了同样的问题，求指导。谢谢！

@cccrrryyy 好的，谢谢！

由于实验室断电，一个正在计算的算例计算中断。瞬态两相流算例，计算时间设置 startTime 为 0 s，endTime为 9 s，输出间隔0.1s，固定时间步长。四核并行。
计算中断时，计算进行到6.56s，输出文件中保留了 0 s 时刻到 6.5 s 时刻的数据。
请问：
我根据输出时间步6.5s的结果数据继续计算（即设置startTime为 6.5 s，endTime为 9 s），这样得到的结果可靠吗？
相比于从0 s 重新计算，两者结果差别大吗？
求大佬指教，谢谢！

@星星星星晴 好的，多谢啦~

@星星星星晴 但是这个报错是运算了一定时间之后出现的，从0时刻开始计算，计算时间设定38s，一般是18s前后出现报错消息，前17s结果良好。

@星星星星晴 谢谢，那应该就是处理器边界条件了...

@cccrrryyy scotch也不行，也会报错......

@cccrrryyy 好的，我试一下，谢谢！

@xiaohsu2006 请问问题解决了吗？我也遇到了类似的问题。谢谢！

@队长别开枪 刚看了一下，设置都没有问题。不过，昨晚试了下，不并行，单核计算的情况下，一切正常，没有出现发散的情况。所以，是并行导致的问题吗？谢谢！

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    note        "mesh decomposition control dictionary";
    object      decomposeParDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

numberOfSubdomains  6;

method          simple;
//method          hierarchical;
// method          simple;
// method          metis;
// method          manual;
// method          multiLevel;
// method          structured;  // does 2D decomposition of structured mesh

simpleCoeffs
{
    n           (6 1 1);
    delta       0.001;
}

请问出现这种情况是什么原因呢？
设置Co <1自动调整时间步长，而Co数还比较小，但是时间步长却已经很小了。如下图所示。

这个时间步之前都很正常（如下图所示）：

请问是什么原因呢？谢谢！