新手小白求教，有没有能够通过读取向量文件数据获取动边界速度的功能

[yzxwh]

@coolhhh 大佬，我把main函数主体放进了python中，利用python可以控制pimplefoam来进行求解，并且可以将每个时步文件的动边界速度pointmotionU进行强制的实时修改，这个时候如果我想让pimplefoam能够时刻读取新的pointmotionU是否就是将pointmotionU的读取和写入功能放在runtime里面？像下面的截图里面一样？

[coolhhh]

@yzxwh

1. 我没用过pointVectorField定义场，这是定义边界场吗，测试是否可以运行？
   我看这个帖子的#9楼，是通过获取边界的patchIndex调用的边界场。
2. 类似pisoFoam这些求解器，定义场只要一开始定义一次即可，把定义场写到循环前面。每个时刻就会自动读取时间文件夹中的数据
3. pointmotionU是自己直接生成文件的吗，那就不需要AUTO_WRITE，可以设置为NO_WRITE。是否用的并行计算，自己强制改pointmotionU，在不同processor文件夹中修改数据是否挺麻烦的？

[yzxwh]

@coolhhh 回答下大佬的问题：1.pointVectorField是定义的边界场；
2.像piso这些求解器，定义场一开始定义一次后，人为调整最新时刻的边界条件，后续的计算不会读取新的边界条件，一直是初始定义的边界条件；
3..在dynamicmeshdict文件里面设置motionSolver velocityLaplacian这个运动求解器就需要初始化pointmotionU，后续计算中会直接生成pointmotionU；应该不是并行计算，这块我比较生疏，没有专门去设置过。

[coolhhh]

@yzxwh 这块我没接触过，说下自己的理解：

1. 如果不人为调整，后面自动生成的pointmotionU，不同时刻的结果是不一样，还是一直等同于初始0文件中设定的值？
2. 大致看了下程序，velocityLaplacianFvMotionSolver是继承velocityMotionSolver，pointMotionU_在velocityMotionSolver.C中已经定义，在主函数中就不需要重复定义了。即使主函数定义了，读取的数值也没传递到定义的类中，也没起到作用。

3. velocityLaplacianFvMotionSolver.C里边有个updateCoeffs()，正常是通过这里更新数值。如果自动生成的pointmotionU一直都是一样的，考虑在updateCoeffs()前强制增加读取文件试试

pointmotionU_.read();

4. 后面计算网格量大的话，单进程估计计算太慢。不过可以先把单进程调试通，再考虑如何改并行程序。

[yzxwh]

@coolhhh 万分感谢大佬的指点，回答一下第一点：后面生成的pointmotionU的确是不一样的，我可能表达有误，由于圆管的运动是通过别的软件生成决定的，所以圆管的运动数据是实时传输到of里面的。of里面的动边界相当于我另一个软件的圆管，但他的运动完全由那边所决定，我的想法是在运行过程中修改最新时步的pointmotionU文件最下面几行的边界条件代码，举个例子 如第一张图所示这是当前最新时刻的边界速度，我通过python强行修改了这个边界速度如下图
我是想让下一步能够以修改的边界条件进行运算，以此类推，往后每一步都会修改边界速度，然后下一步以新的边界速度运算，不过我的尝试没有成功，pointmotionU边界的读取似乎只与初始时刻相关。
关于大佬说的剩下三点我再好好理解理解，再次感谢大佬的悉心指点。

[coolhhh]

@yzxwh 确实可能MUST_READ只是在初始启动时候会读取文件而已，大概率需要手动增加pointmotionU_.read();到程序中，实现每个时间步调用时候先读取数据

[yzxwh]

@coolhhh 好的大佬

[yzxwh]

@coolhhh 大佬，我这边在velocityLaplacianFvMotionSolver.C添加了读取文件的代码pointmotionU_.read();文件的确是每个时步都被读取了，但是读取的边界条件一直是修改之前的。我还采用赋值的方法强行给pointmotionU_赋值，代码如下

强制读取当前时刻修改后的pointmotionu的数据存到temppointmotionu中然后再赋给pointmotionu，两个结果都在运算中打印出来了，结果发现，temppointmotionu是成功读取修改后数据的，但是赋值后的pointmotionu依旧是修改前的数据。是否是该类型无法被直接赋值，望指点。

[coolhhh]

@yzxwh 最后输出语句如果直接写pointMotionU_ ，输出结果会等于tempPointMotionU吗？不能直接等于应该会报错的，可以再新建一个tempPointMotionU2，直接幅值tempPointMotionU2 = tempPointMotionU；看看输出结果

Info << "Reading pointMotionU_ from tempPointMotionU from file: " << pointMotionU_ << endl;

[yzxwh]

@coolhhh 大佬，我已经尝试过了 直接幅值tempPointMotionU2 = tempPointMotionU；是可以的，但是pointMotionU_ 还是没办法被赋值到，同时也不存在报错的问题。这是由于pointMotionU_ 是存在于velocityMotionsolver中的protected类型所导致的吗？
下面是velocityMotionsolver的C文件构造函数中初始化列表进行初始化的。

新手小白不太明白，望指点。

[coolhhh]

@yzxwh
意思是赋值完，直接 Info << pointMotionU_ 也结果有误？
你是把赋值pointMotionU_ 写在velocityMotionsolver.C文件中吗？还是写在这个函数外边？关注变量的作用域范围，可以把velocityMotionsolver理解为一个子函数，在子函数定义的pointMotionU_为局部变量，如果在这个函数外边重新定义变量赋值的话，并没有修改子函数中的值

[yzxwh]

@coolhhh 大佬，赋值之前与赋值之后的PointMotionU_输出是一样的。我是在velocityMotionSolver的类的public类型中添加了一个成员函数readPointmotionU（）；

并在其中进行了tempPointmotionU的读取与对PointmotionU_的赋值操作，如下：

然后我在velocitylaplacianfvmotionsolver的c文件中调用了这个函数：

我这样算是在子函数中定义的变量吗？

[coolhhh]

@yzxwh

1. 看起来没问题，我也按照你的问题写了个简单例子测试，输出结果是正确的。

#include "fvCFD.H"

class Parent {
protected:
    int protectedVar = 5;

public:
    void modifyVariable() {
        int tmp = 10;
        int tmp2 = tmp;
        Parent::protectedVar = tmp;
        Info << "tmp2 = " << tmp2 << nl << endl;
        Info << "Parent::protectedVar = " << Parent::protectedVar << nl << endl;
    }
};

class Child : public Parent {
public:
    void accessParentFunction() {
        Parent::modifyVariable();
        Info << "Child::protectedVar = " << Parent::protectedVar << nl << endl;
    }
};


int main(int argc, char *argv[])
{

    #include "setRootCase.H"
    Child child;
    child.accessParentFunction();

    Info<< "End\n" << endl;

    return 0;
}

---

2. 看到velocityMotionsolver.H函数中有个返回引用函数：pointVectorField& pointMotionu()，试试能不能创建一个引用，然后修改引用的值，看看结果会不会变化。

pointVectorField& pointMotionURef = pointMotionU();
pointMotionURef = tempPointMotionU;
Info << "pointMotionURef = " << pointMotionURef << endl;
Info << "pointMotionU = " << pointMotionU_ << endl;

[yzxwh]

@coolhhh 大佬，我这边已经把问题解决了，如图所示，添加了一段delete代码删掉了旧对象，再重新分配，就可以进行赋值了。

大佬的悉心指点对我颇有助益，万分感谢大佬在这段时间的指导与帮助！！！

[coolhhh]

@yzxwh 厉害厉害，我也学习了

[李东岳]