如题...

好吧！我先自己回复下自己。在稀疏矩阵存储中，通用有两种格式，如下矩阵

A = [ a00 a01 0 0 0 a11 a12 0 a20 0 a22 0 a30 0 0 a33];

引用一篇博士论文 CRS scheme stores the non-zero elements the matrix in left-to-right and top-to-bottom order in vector v (row-wise storge). 存储格式如下

v = [a00, a01, a11,a12,a20,a22,a30,a33] cI = [0,1,1,2,0,2,0,3] rI = [0,1,1,2,0,2,0,3]

但是在openfoam中，不是这么存储的
openfoam是一种 LDU 格式 对角线 lowerAddr() upperAddr()

A = [ d0 u0 u1 0 0 0 l0 d1 u2 0 0 0 0 l1 l2 d2 u4 u5 . . . ] l = [l0 l1 l2 ...] d = [d0 d1 d2 d3 d4 ...] u = [u0 u1 u2 u3 ...] L = [0 0 1 1 2 2 3 ...] U = [1 2 2 3 3 4 4 ..]

具体可参见lduAddressing.H 头文件
因为计算流体力学中，大部分矩阵都是对称（因为一个internal face对应两个体）除了特别的边界条件外，基本上都是对称矩阵。而矩阵对角元素就更好求了，大部分情况直接把该行相加就得到主对角元素。这种sparse matrix设计是符合cfd计算原理的，但是我的问题还是没解决。看到相关资料分享一下。
稀疏矩阵也分好多种，openfoam中的renumberMesh就是尽量保证在是多对角阵。
谢谢！

@东岳 好的，了解了！！谢谢您

在多区域（液体和固体换热）计算中，某一个物理量（温度T）在两个区域间的interface是一个连续的量，即没有热阻，这个interface边界条件该怎么设置？
在OpenFOAM自带的多区域算例中，这种耦合面用的边界条件是turbulentTemperatureCoupledBaffleMixed，这个其实是根据边界面的值是根据边界上的值来计算传热系数的，里面写的很麻烦，有点看不懂。有没有哪位老师指点一下？多谢多谢~

如题，OpenFOAM中有简单的命令可以快速实现吗？试过refineWallLayer，只能添加一层网格。

@东岳 第一次发帖求助，就能得到东岳老师的回复，真的非常感谢！
（1）是根据laplacianFOAM稍微改了下的求解器，但是如果是探针probe的话是可以运算出结果的。那cellSource不运行function的话，我需要怎么改下求解器吗
（2）还想请教下，算例边界条件设置的问题，左图是压力，右图是温度。
左图压力上下边界条件均为zero flux，右侧问题上下右边界条件为zero conductive flux。对于zero flux和zero conductive flux，边界条件是不是应设置为symmetryPlane，还是zeroGradient。我计算出的结果差别很大，不知道问题出在哪里，所以想确定下是不是我边界条件设置的不对。
case.JPG

@悬铃神木 谢谢。

@气固两相流 fluidisedBed.zip
是一个TFM算例的设置，烦请大神们帮忙看看算例有什么问题，不胜感激

@random_ran 1243.jpg cyclic边界得到的压力梯度恒为零，跟我所需要的连续的周期性边界条件不太一样

@东岳 帖子好像沉掉了，再次求助李老师。主要就是顶部边界条件不知道怎么设置，顶部是定压，也就是P是恒定的，但是流体可以自由出入。而方程中求解的是prgh，所以不知道prgh如何设置边界条件了，可能还需要与出口的速度边界条件保持一致性

@史浩 恩恩，是的，vof文件都得拷进去，感谢感谢

对 :duang:

可以直接忽略icoFoam，直接用pisoFoam，icoFoam是个demo类的求解器，

@myler 用这个代码是将最大值输出至cellSource1的文件夹吗？我也把这个代码添加至了controlDict里面，但是并没有看到最大值的输出结果啊

@zhangdong 谢谢

可以没问题，直接用吧 :duang:

@bestucan

@quotation
People sometimes ask me if it is a sin in the Church of Emacs to use
vi. Using a free version of vi is not a sin; it is a penance. So happy
hacking.
@author Richard Stallman
@end quotation

在wall边界里找到了constantAlphaContactAngle类型，这个类型能同时满足无滑移边界吗？

大家好，我在openfoam中的动网格设置遇到了一些问题，d.png 首先在motionSolverLibs 不知道填写什么，我发现填写不同的Lib里面的可用的求解器有所不同，但是也有相同的，请问相同的部分他们达到的求解效果是一样的么？
我在openfoamwiki中查到solver solidBody可以指定物体的运动。但是求解效果不好，然后solver rigidBodyMotion和sixDofRigbodyMotion我感觉是差不多的。有没有关于这些动态网格求解器的资料，想看看他们具体是解决什么样的问题的？谢谢大家。

@Jacobian 我尝试了您的方法，其实不需要设置rho.water、rho.air。我设置了水的密度随高度按照双曲函数变化，初始时刻以及边界密度是没有问题的。但是运算开始后出现问题，密度分布立即按照thermalphysicalProperties设置变化，也就是说按照目前的方法无法控制密度与其他物理场的联系。实际上我的密度场是可以随着流场变化的，因此我觉得应该添加一个浓度输运方程，然后重新定义热物性，使密度随之变化。