openfoam中cyclic周期性边界的问题
-
@Tong 好的!谢谢,您说的文章我也看过,不过涉及到具体细节他没有将,OF2006中已经实现了petsc2foam了。我正在看源码,在这个网址上:https://develop.openfoam.com/modules/external-solver/-/blob/develop/src/petsc4Foam/utils/petscLinearSolverContext.H
我最终的想法是实现一个自己的预处理器或者矩阵求解器然后继承到OF中。 -
@Micro
foam-extend 里面好像也有不少新的求解器,或许可以参考一下
https://github.com/Shadow-fax/foam-extend-4.1/tree/93fffffdb7453736d95f5535f9c8eed1054e48f6/src最近发现对于全耦合求解的矩阵(foam-extend里的coupledMatrix)的确需要针对物理问题的新预条件器来加速收敛。
-
@Micro
个人观点,线性方程组求解来说,感觉目前主要用的是克雷洛夫子空间方法+多重网格方法这两个大类,已经能在理论上比较好的解决求解方程组了(复杂度大概在O(N^2)~O(NlogN))。当然对于不同线性代数库对这些方法的实现效率有高有底(针对体系结构和并行的优化可能不同),但是感觉如果优化到位大概也就是到那个程度了,在同一套系统上不会有数量级的差距。
我之前做实验发现,对于有些多物理量耦合在一起的分块矩阵,通用的预条件器好像对降低条件数作用比较有限(主要也是由于物理上的特殊性造成的),针对某一类物理问题,似乎需要一些特殊的预条件器来有效的降低条件数,听说预条件这方面的论文也挺多的,不过比较偏数学,我不太懂。 -
@Tong 老师,您好。我想请教您一下,interface是针对周期边界、多求解域(如旋转机械)、并行边界这类复杂边界的吗,对于常规边界条件有影响吗,我理解的是边界条件的离散系数已经添加到lduMatrix中的对角线系数和源项中了。在测试中发现initMatrixInterfaces和updateMatrixInterfaces对于X和B都没有影响啊
-
@Micro 请问老师,您这边搞清楚interface并行边界如何处理的了吗
-
@dzw05 请问老师,您这边清楚interface并行边界如何处理吗?我想从openfoam中拿到完整矩阵数据,但是并行的话只能拿到各个线程对应的子矩阵,如何知道并行边界的关系将其组装成一个完整的矩阵
32/34