CFD中文网

用多重网格算法求解线性方程组Ax=b，算法求解过程中按理只需要给定系数粗细网格的系数矩阵就可以了，但在看算法理论时，都涉及到粗细网格的位置关系的判定来确定插值矩阵、限定矩阵。这个就需要在算法上导入网格具体的物理位置了吧。这样好像不对（按理是只要给定系数矩阵A，b,就能求解Ax=b）对于这点一直想不明白:mihu:

@一颗鸭蛋

尝试过docker 没弄明白 当时M1还没适配docker 后来也懒得弄了
我是用parallel 装的ubuntu arm 自己编译的，跑的起来完全没问题。
跑分里面的M1 就是我的mini

@李东岳谢谢东岳老师

1.JPG 捕获.JPG
就是这几项在方程好像没有它们 它们的作用是什么？

谢谢东岳哥

刚刚尝试了下将double型变量赋值给一个新的体标量场，然后比较两个体标量场的大小就可以通过编译了:chitang:

b =mesn.C().component(vector::Y)/a;//b被声明为一个体标量场记录距离壁面的无量纲距离，a为特征长度 Prt = 1.25;//Prt也是被声明为一个体标量场，这里先把Prt当做了一个中间变量，Prt为我研究中的湍流普朗特数，其对流体温度梯度影响很大，应当是一个函数而非常数。 if( b <= Prt ) { Prt = 0.45; } else { Prt=4.75; }

今晚先算用编译好的求解器算一遍对不对，结果好的话，明后天再放上来。

我的上述计算中没有用到估算的方法，那论文里的是不是因为边界层网格采用的估算的行为，所以才会导致friction Reynolds numbe因人而异，出现微小的差别呢！

就是从定义来求解，k还好说，epsilon和omega牵涉到一个特征长度的问题，这个取起来就不容易了，只能多调试。

@dzw05 耦合计算和单域计算的流体域是使用上面提到的模型求解，求解器则是不同的，一个是壁面导电耦合求解器，一个壁面是绝缘求解器。我的耦合也只是壁面和流体的电势耦合。

封闭方腔自然对流

这个不应该出现问题啊，网格很规则，有实验结果可以对比么？还是自己算着玩的。有实验数据对比我可以给你调调，网格少的话

公式就是那几个，你也知道的

@马乔 您好，我参照你的方法将transportProperties代码写进去，但是很奇怪的是，报错显示，request for dictionary transportPropertie ,而可获得的字典里是transportProperties,两者总是存在单复数的差别（就像你贴的代码里一样）。本来想着直接用令锑度边界条件绕开这个错误，但是最近发现还是得用alphat壁面函数。真是很奇怪呀。

最近在家呆的闲，又没有实验室的资源，，想看看磁流体的流固耦合的代码和书最好中文的，，，求推荐啊～

感谢鸭蛋分享

问题解决了，是因为我将Omega设置为零了，所以在计算的时，湍流模型在算湍流黏性系数是nut=k/Omega 出现除零错误。

@Samuel-Tu 嗯嗯

我当时是由于湍流模型虽然编译通过但其实还是有问题的，具体就是H文件申明的虚函数divDevRhoReff，在C文件中没有具体的实现。你可以根据H文档最后几个虚函数（virtual函数）再好好对照一下看看C文档中是否都有对应的具体实现。其次我是2.3.1版本。

求解器编译成功了，，，
将Options中动态链接库中的：-lincompressibleRASModel 去掉就编译通过了。

感谢楼主分享

@星星星星晴
谢谢你，我已经把这个错误解决了，是因为我在files文档的第一行kOmega231/kOmega231.C应该写成kOmega231.C

@东岳 在 在OF2.3.0版本如何编译低版本下的湍流模型？ 中说：

模型计算公式拔出来，自

懂了，，，谢谢东岳老师!