CFD中文网

@李东岳 多谢李老师鼓励，还得感谢李老师的讨论答疑和提供的交流平台 。

@ShuanmingZhao 好的，客气了，一起进步。

@ShuanmingZhao 你好，我这个工作已经公开发表了：https://pubs.aip.org/aip/pof/article-abstract/36/10/104111/3315665/Multi-region-coupling-computational-fluid-dynamics?redirectedFrom=fulltext
然后代码也开源了：https://github.com/rararahui/chtPorousMultiRegionFoam
不过我水平有限，仅供参考哈

@李东岳 在 关于进行瞬态计算矩阵迭代过程的数据操作 中说：

是要改很多私有成员变量啥的，来回嵌套。我想想有没有什么好办法。

谢谢李老师关注，确实我本来也以为改一改矩阵求解器就行，结果发现OF里嵌套继承的类太多了，越改越麻，所以才来提问，希望有更好的方法。

@李东岳 在 关于进行瞬态计算矩阵迭代过程的数据操作 中说：

是不是不应该记录时间步信息？

李老师您好，这里是因为我考虑的每个时间步都要输出多个场，所以每个时间步输出数据的filename应当是与当前时间步相关的。所以只在PCG类里做改动无法读取时间步信息的话，无法定义合适的filename，对后处理会带来麻烦。而且没有时间步这种全局变量信息，仅在PCG类里做改动无法定义合适的filename的话，为了不覆盖之前的输出，似乎只能在同一个file里add数据，这种方法感觉是可以记录所有时间步及其迭代步信息的，但是后处理也比较麻烦，需要重新核对时间步，感觉不够优雅，所以我希望有更好的解决方法。个人技术和思考比较浅薄，感谢大佬赐教。

@李东岳 李老师您好，按照您提供的思路我尝试对矩阵求解器的类进行了修改（以PCG为例），我修改了PCG.H和PCG.C，但是在数据输出这块还有点问题想请教一下。我本来想按您说的把场提取到顶层输出，所以想在PCG的solve中多加一个输入用来存储一系列迭代的场，然后再到顶层输出。但我发现这个添加函数输入项的过程牵一发动全身，从lduMatrix、fvMatrix到fvm等一系列的类和源码都需要进行调整，感觉有点杀鸡用牛刀了。但是如果直接在PCG.C中输出好像获取不了时间步的信息，不知道怎么记录输出不同时间步的结果。所以想再次请教一下有没有更好的方法实现这一过程？谢谢。

@李东岳 好的，大概了解了，谢谢李老师，那我这边去尝试一下。

@李东岳 谢谢李老师，这个迭代过程我大概了解了。关于这个输出我也想顺便请教一下，对单个field是用什么命令输出呢，是psi.write(os)这个命令么？

@李东岳 是的李老师，就是这个意思

@李东岳 李老师你好，我这边主要是想做基于深度学习的加速收敛相关的工作，想初步输出各个迭代步的结果用Python和神经网络处理分析一下，然后确定算法后植入到迭代过程。这个后续植入应该直接在求解过程赋值就行，但是目前这个输出每个迭代步的结果不知道怎么实现，感觉需要修改fvMatrix的solve函数？

各位老师好，我最近在OpenFOAM中尝试对瞬态计算矩阵迭代过程中的物理场进行操作，主要希望对各个迭代步的计算结果进行输出、分析和修改，看了相关源码感觉不得要领，希望各位老师不吝赐教应该如何实现。

@李东岳 李老师您好，之前确实没考虑到压力波。不过我这里的泄压过程是缓慢过程，是通过密闭空间周围的多孔介质渗流漏出去的，给的流速很低，大概在10e-4m/s量级，也会有这种问题么？然后考虑到空间内的热对流，我可能更倾向于使用buoyantPimpleFoam求解器？

@李东岳 李老师您好，我尝试了rhoPimpleFoam和buoyantPimpleFoam

各位老师好，我最近在做一个密闭空间充放气过程的数值模拟，压力变化范围大概在常压到MPa级，过程中温度随着压力的升降而升降。但是发现模拟得到的温度变化要比实验数据剧烈得多，实验中温度变化范围较小，但是模拟得到的温度上下限都宽很多，尤其泄压过程温度下降十分剧烈，感觉不合常理。本来猜测可能是状态方程的问题，尝试了OpenFOAM中的理想气体状态方程和PengRobinson状态方程以及热物理模型的一些不同设置，发现结果都是这样，请问老师们有相关情况的研究经验么，希望不吝赐教。

@huangyuhui723 我尝试了可以通过U.boundaryField().oldTime()进行调用，可供参考。

@李东岳 李老师，我尝试了oltTime是可以的。不过还有个新的问题，因为我想做的边界条件是要在之前时刻的基础上自适应调整的，如U(t)=U(t-1)+dU，因此除了上一时刻的内部场，还得获取上一时刻的边界输入U(t-1)，但是边界场好像没有oltTime这个功能，是否需要建立一个全局变量来存储U(t-1)？

@李东岳 在 求问如何在自定义边界条件中获取之前时刻的物理量数值 中说：

如果只需要前两个时刻，可以用p.oldTime()以及p.oldTime().oldTime()

存储scalar就行，不用场

是不是还得看哪个网格点的压力？

谢谢李老师，不用具体哪个网格点，该区域的压力基本上是均匀的，取一个平均值就行，那我去试试看。

@李东岳 李老师，我看了一下，好像没有涉及提取之前时刻压力值的。我只需要之前两个时刻即可，而且该区域压力基本稳定，存储scalar就行，不用场。我有一个思路是在主程序中保存之前时刻的值，然后传递给边界条件，不知道这样是否好实现？

各位老师好，我最近在尝试开发一个入口边界条件，希望根据压力的变化调整进口的质量流率，主要目的是为了使某存在泄露的区域保持压力恒定。我这边的思路是想尝试类似PID控制的算法，基于压力的变化调整入口质量流率，但是不知道有什么方法可以获取之前一个或多个时刻的压力值，并在当前时刻边界条件中使用，希望有经验的老师能指点一下，感谢！

@wangfei9088 感谢老师，不好意思最近没来得及回复。我后面还是尝试了直接写了TEqn，没有去折腾热物理库了。很感谢老师的指点。