CFD中文网

@fangyuanaza :146: :146: :146:

@田畔的风 是的，感谢大佬，果然可以了。得抓紧学python了，批量后处理急需。
bf4b5b57-9231-415f-b4a7-970bce7e09e5-截图_20231102134831.png

这里有介绍：link text

类似的需求：

How do I add some other package to the container?

The easiest way is to modify your Dockerfile and rebuild. Particularly easy if the software you want to add has an Ubuntu deb package.

里面也有 docker file 的示例。看起来这个比较可行

没有经验，一直在 linux 上跑

@李东岳 感谢老师的回答，抱歉才看到回信，最终解决办法是强制读取边界信息，检测距离边界最近的网格点数据并给出。黄色震荡部分是因为插值了内部网格点和边界值的结果。

@fangyuanaza 解决了这个问题，代码分享如下，以免大家有需要：

const volScalarFields& p = this->db().objectRegistry::lookupObject<volScalarField>("p")

经过Info输出发现，确实是每一步的压力场

@luofq-sysu 第一个问题：是的，我取的是流向速度时均值。
第二个问题，是的，是不同流向位置。我当时做的是一个平板流动，平板上流向的网格113个。150个不是一个必须的值，看你算例的情况跟需求。我的是二维流动，因此没有展向的问题。

请教各位老师~
****目标：***希望给雷诺应力增加高阶项来修改湍流模型，高阶项的形式是固定的（标量对称张量）
***想法：***因为雷诺应力本身就是对称张量，其中的标量是随时变化的，机器学习生成的。如果用include **.H的方式来修改湍流模型，则每次需要编译，比较耗时，（因为增加项会不断变化）。所以想通过在constant文件夹中以.H形式存储标量表达式，从而在turbulenceProperties中读取文件，以系数的形式修改表达式，由于系数是实时读取的，因此只需要编译一次。.H文件中就是一个标量表达式，命名为zetai.H
****问题：****如果仅有zetai.H，内容如下（I1 I2是已经在湍流模型中计算好了的体标量场）：
8b12a85a-ca5c-48ca-8d0a-0e643d6d5861-image.png
并在湍流模型中读取：
fb6a553a-b074-470e-97ef-4ccc9cc304ec-image.png
如何能将zetai中的内容转换为体标量场进行计算？

@李东岳 谢谢李老师的回复~ 教师节快乐~ 是一个标准的channel。但是雷诺数非常高，Re=80 million, 没有8000万网格，是NASA官网上的一个标准算例，独特的地方就是雷诺数非常大

请教一下各位CFDer, 很多算例的设置是无量纲化的，那么有没有热物性参数是怎么无量纲化的介绍呢？或者比如有理想状态下空气热物性参数的无量纲化结果？因为不同的无量纲化过程可能结果是有差异的，比如说压力，有的用rhoU^2来无量纲化，有的用0.5rho*U^2来无量纲化，导致压力相差一倍。

@李东岳 我重新定义了k_, 将其强制变成const，但是还是报错，错误如下：

error: no matching function for call to ‘ddt(const alphaField&, const rhoField&, const volScalarField&)’ + fvc::ddt(alpha, rho, k_) ~~~~~~~~^~~~~~~~~~~~~~~~

可见k_前面已经有const，但是还是不行，可见应该还是要修改alpha与rho的类型？但是我重新赋值定义alpha， rho为volScalarField, 会报初始化错误，不能直接这样改，修改方式如下：
const volScalarField alphat = this->alpha_; 再在fvc中调用fvc::ddt(alphat, rhot, k_)

@李东岳 是的~编译成功了！ 给李老师打call! :high:

@fangyuanaza
:146:

@bestucan 谢谢回复~ 主要是2%的条件放在原来的网格上又能算，原来的网格跟现在的网格只是相当于平移了一小段距离而已。

是我没有解释清楚么？这里的genotype[i]就是输入的字符，因此minus 1.0 2.0 对应三个字符串，从i=0开始，对应minus, i++变为1， 返回第一个递归函数evalGeno对应的i=1,也就是进入else, i++ 为2，numValue为1.0，继续减号，然后到第二个递归evalGeno，此时i为2，所以按程序流程就是1-2，但是程序测试结果是2-1的结果，为什么会这样呢？请各位老师指点一下

已解决。。。先给出ADD的表达式，然后对其进行限制。修改代码为：
屏幕截图 2021-04-07 112950.png

我错了，是$S_{ii}=0$。因为满足不可压缩流体的连续性方程。尴尬。。。:shangxue:

@李东岳 在 跨音速求解器选择 中说：

am求解器本源是打

@李东岳 谢谢李老师！我发现问题了，这个报错其实就是最后的correctADDStress的问题，所以看这个报错提示就是看最后一个信息，错误的原因是我用wclean去移除编译的库的时候其实没有移除，所以后面不管我怎么改，都是调用的原来的库，解决办法就是手动删除openfoam-7的platform中生成的自己编译的库，重新编译就好了:136:

@bestucan 谢谢你~ 我找到问题了，就是ADD范围不对，我用了最值函数限制了就好了，这一点论文里面没有提，很不好找，谢谢你~

@chengan-wang :146: