CFD中文网

@李东岳 不好意思，李老师，工作原因回复慢了一些，我目前只在尝试稳态simple方法的实现，那就先只考虑边界面 + 的贡献么，我去修改代码再调试一下，非常感谢！

@李东岳 @李东岳 李老师，我可能还是有点绕糊涂了，实际上我是在自己写代码尝试复现有限体积法的simple算法，那么在构建如图压力泊松方程的时候， 我的迷惑点就在于在边界处单元对该项 进行所有面求和时，边界面的Apf该如何处理，是直接不考虑其贡献不加它，还是把其取值为边界处的 + ，还是说是其它的处理方法？

@李东岳 李老师，Ap定义我大概了解了，意思是计算域内部面的Apf是直接插值得到，而实际上Ap本身是对当前vol所有面如图公式求和得到的， ，那么在边界处的Apf也就考虑在边界面处的 和 就行了么。

@李东岳 李老师，那意思是说在如图公式构建压力泊松方程的时候，边界上的Ap并不做出贡献么，公式右侧的外部面求和不用考虑边界面么？我对该过程的理解有限，还希望李老师解答。

@李东岳 李老师您好，我最近也对这个问题感到好奇，看了您的回答，也就是说phiHbyA在边界处直接赋值为边界处的通量phi即可，然后压力泊松方程里还有Ap这一项也要插值到面上，还想请教那么Ap在边界处是如何处理的呢？

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

@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()进行调用，可供参考。