如何重新读取已有的时间步里的场

[VamosZhi]

在一个自定义solver里，需要在特定的时间步时，重新读取当前时间步文件夹里的场（可以是速度场或者压力场），我使用了：

volScalarField p
(
    IOobject
    (
        "p",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ_IF_MODIFIED,
        IOobject::NO_WRITE
    ),
    mesh
);

但实际上MUST_READ_IF_MODIFIED是不起作用的，直接在solver里写“p.read()” 也是不管用的，值并没有发生变化，还是一开始初始时间步从time文件中读取的值。在CFDonline上我看到了解释：
https://www.cfd-online.com/Forums/openfoam-programming-development/122189-runtime-modification-volscalarfield.html
和IOobject有关系。但是不知道大家有没有更简单的办法。谢谢

[李东岳]

还是一开始初始时间步从time文件中读取的值。

你是不是controlDict里面读取的是0文件夹的

[VamosZhi]

@李东岳 李老师，controlDict读取的也是latestTime （比如是t1）, 然后他读取了t1时间步的场，比如一个force source term ‘f’。现在我在solver里想定义，每运行ndeltaT的时间步时，从新的（t1+ndeltaT）time文件中读取新的 ‘f’场， ‘f’场的值事先准备好了，放在（t1+ndeltaT，t1+2ndeltaT，t1+3n*deltaT，.......）里面。问题是，读了t1时间步的场后，后面的怎么都读不进去，仍然保持t1时候的值不变，不管是设置了“MUST_READ_IF_MODIFIED”，或者强制‘f.read()’。

volScalarField f
(
    IOobject
    (
        "f",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ_IF_MODIFIED,
        IOobject::NO_WRITE
    ),
    mesh
);

[李东岳]

你的意思是，求解器还没有求解的时候，自己定义了一系列的时间序列文件夹0 1 2 3 4，求解器运行到4的时候，去读4文件夹下的文件？

[VamosZhi]

@李东岳 嗯嗯对的，是这个意思。李老师您看看有什么解决办法

[coolhhh]

参考这个帖子新手小白求教，有没有能够通过读取向量文件数据获取动边界速度的功能，看看读取之前先把p的值删掉，再执行 p.read()效果咋样

[VamosZhi]

@coolhhh 谢谢！这确实是一个可行的办法，我看了那个帖子，可以通过删除后重新赋值来解决。我试了一下，直接敲代码“delete p”,“delete U”啥的是不行的，不知道其对应的对象是怎么表示的（像那条帖子一样）。我上网搜了一下，还没搜到。不知道您有什么想法

[coolhhh]

@VamosZhi 试下下面这种方式，看看能不能更新读取数值？

        // 清空 p 的值
        p = dimensionedScalar("zero", p.dimensions(), 0.0);

        // 重新读取 p 
        p.read();

[VamosZhi]

@coolhhh 嗯嗯，刚刚试了，不行。这样子的话，p的值就为0了，再read()一下还是0，读不进来。
我在想，能不能借鉴从constant/transportProperties中读取nu的方法：

dimensionedScalar nu
(
    "nu",
    dimViscosity,
    transportProperties.lookup("nu")
);

这样子是可以实时更新的。但是不知道如何在runTime.timeName()中lookup场的值，这个代码怎么写还没搞明白，一直报错。
这样写是报错的:

p = dimensionedScalar("p", p.dimensions(), runTime.timeName().lookup("p"));
log文件：
interPhaseChangeFoam_justPp.C:120:72: error: ‘class Foam::word’ has no member named ‘lookup’
       p = dimensionedScalar("p", p.dimensions(), runTime.timeName().lookup("p")

[coolhhh]

@VamosZhi 我是在icoFoam尝试的，没有复现出你的问题，不知道为啥你的代码p.read()不生效

[VamosZhi]

@coolhhh 嗯嗯，确实奇怪，我这边是控制每迭代特定时间步读入一次：

各种方法都尝试过了。输出log文件：

比如我自定义的ddtp还是0,我试了p,也是一样的。

[coolhhh]

@VamosZhi 还真不知道怎么解决，有没有一个最小的测试算例复现这个问题的，试试看怎么处理

[VamosZhi]

@coolhhh 好的，谢谢！我再试试看