CFD中文网

@coolhhh 好的好的 真的万分感激

@coolhhh 太感谢大佬了，大佬犹如神助，我的问题目前看来已经得到了解决，被这个问题真的困扰了许久，之前有个帖子大佬也帮了我很多。我还有个问题就是通过这个create patch的方法face的Id得到了更改，cell的id是否也会有变化。

@yzxwh 忘了说了 我用的是pimpleFoam求解器

@李东岳 十分感谢李老师愿意帮助我，我是利用of与外部软件耦合模拟管道运动对多孔介质流场影响的，我这边就模仿movingcone算例单独建了个管道算例，与前面所说的网格略有区别，但问题是一致的。模型如下

给管道赋予x方向0.8m/s的速度，由于我之前与外部软件耦合模拟多孔介质环境，会多一个耦合力源项，在周期性边界条件下能够运行，单纯的流体算例似乎无法运行（如果改成墙边界是可以运行的），可能是不满足守恒定理（这一块我比较浅薄）。文件和模型我已放在下方，再次感谢李老师的帮助。
文件：moving_tunnel.rar
模型：tunnel_model.msh

@李东岳 李老师，我尝试了使用createPatch来生成，结果还是报错，这个方法好像无法解决id为倒序的问题
互为周期性边界两个面的face id顺序如上图所示，我是通过在fluent meshing里面设置好了周期性边界，但是植入到of中周期性边界就直接变成了wall ，随后进行了多次尝试，但都没能解决id倒序的问题，老师是否有更好的想法，望指教。

@李东岳 好的老师

@李东岳 李老师，我是先用fluent mesh画了个网格，然后植入到了of里面用的是fluentMeshToFoam这个功能，原本生成的是wall边界，然后手动改成了cyclic边界。

@李东岳 感谢李老师的答复，这个cyclic映射是按照face的id顺序进行处理的，我在模拟时流体的连通性出现了如下的情况
两个边界面的连通是反的，造成了如图的情况。我看of中导入.msh文件后是自动生成面并进行编号的如下图

那么这个id顺序是否人为可以控制。希望老师能解答！

这是fluent画好的网格植入后得到的左右两个边界条件，我进行了修改，想改成周期性边界如下图 但是在后续的检查中发现，left面的网格顺序从图形上看是从下往上排列的，而right是从上往下排列的，这就导致了如果按id号顺序进行映射的话，二者在物理空间上的映射是倒序的，即left上部的网格面映射了right下部的网格面，left下部的映射了right的上部面，如何能让他们在物理空间上对应解决这个倒序的问题 求大佬指教！！

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

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

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

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

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

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

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

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

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

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

@coolhhh 好的大佬

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

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

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

@coolhhh 谢谢大佬指点，我明早试试看。

@coolhhh 大佬，你指的是每一个时间文件替换掉相应动边界的数据吗，这个我用python实现了，但是我在调试中发现中途改变动边界数据并不会使下一个时步读取这个替换的数据进行计算，除非是终止计算，改controldict的starttime才能重新读取替换后的速度。

@coolhhh 谢谢大佬的答复，我去研究研究timeVaryingMappedFixedValue的源码。