@aiweimo

k-w方程是低雷诺数方程，本身不需要壁面函数

严格来说，应该是 $k-\omega$ 方程可以一直算到 $y^+<1$，但是，$k-\omega$ 方程也可以结合壁函数来算$y^+>30$ 的网格。相比之下，标准$k-\varepsilon$模型的$\varepsilon$ 方程则只在log-law region之外才严格有效，不可直接用来算很小 $y^+$ 的网格。
OF 里面的那些湍流壁面边界条件，有一点容易引起误解的是他所有的名字都带 wallFunction，而大部分人看到 wall function，想到的都是标准壁函数，也就是 log-law region的 $U^+=\frac{1}{\kappa}\ln(Ey^+)$ 以及层流底层的 $U^+=y^+$ 。 实际情况下原比这个复杂，因为很难保证整个region 都满足 $y^+ > 30$ 或者 $y^+<1$，所以，对于实际问题，湍流模型的壁面处理需要考虑如何处理$1<y^+<30$ 的区域。文献中有很多相关的研究，openfoam 的 omegaWallFunction 其实就是一种，里面构造了一个blending function 把适合 $y^+ > 30$和 $y^+<1$的 omega 公式blend一下，来让omega 在任意 $y^+$ 下都能算出来一个相对合理的值。
所以，OF 里面那些 wallFunction 的边界，可以理解为定义多种不同的 wallTreatment 方式，具体怎么用，需要结合着实际情况来。2 楼提到的就是一种常见的用法。文献中壁面处理的方式太多太多了。。。

至于 $\omega_{wall}$ 指定为零梯度或者固定值，这种做法不太合理。因为 $\omega$ 在壁面附近需要满足某种渐近特性（asymptotic），零梯度或者固定值都不满足这个条件，虽然可以可以算出来结果，但是不太符合物理。

我最近看了下，

替代文字

上面这个方程应该是稳态算法的方程。里面MRF.DDt(U)用于处理$\Omega\times\bfU_I$。

96bfbd01-3d8c-469a-a303-a79d75b4d80c-image.png

但是上面这个方程不是OpenFOAM植入的方程。在瞬态情况下应该是下面这个：
$$
\frac{\p \bfU_I}{\p t}+\nabla\cdot(\bfU_I\bfU_R)+\Omega\times\bfU_I=...
$$
下面的瞬态代码分别对应：

tmp<fvVectorMatrix> tUEqn ( fvm::ddt(U) //第1项 + fvm::div(phi, U) //第2项 + MRF.DDt(U) //第3项 + turbulence->divDevReff(U) == fvOptions(U) );

:chitang:

@Samuel-Tu 在 local index to global index 中说：

但是写入solver的时候调用不了这个函数。说没有定义此函数，如何调用这个函数呢。

大佬您好，请问您解决这个问题了吗？我也想问一下怎么调用这个函数，感恩！

大佬您好，为啥 那请问局部索引应该怎么显示呢？谢谢！

多谢！已更新

@李东岳 感谢李老师的回复，谢谢

@李东岳 好的，谢谢东岳老师！

@李东岳 岳 东岳老师好，上面公式S其实表达的意思不一样的，我这里用错了，现在存在的困难是怎么求解这个压力方程值p'，感觉好难写出来只关于p'的方程，再进行求解，希望能得到东岳老师和各位前辈的指点，谢谢！
111.png

你这个最好用wedge网格来计算就不会出现这个问题了

@王慧博 在 矩阵求解器 中说：

LUsolve

您好。请问LUsolve是OF自带的函数吗？ 还有一个问题，为啥LUsolve可以直接调用而不是通过对象访问成员函数的方式调用呢？难道LUsolve是一个全局函数吗？谢谢大神！！！

主要输OF中的残差是线性方程组求解的残差，不是时间步迭代的残差，您说对吗？李老师。另外，是不是只有稳态问题我们才会关注线性方程组求解的残差残差，瞬态问题的residual 有意义吗？

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

@李东岳 WALE不需要vanDriest吧？

@huntersong 哇，好详细，我去试试， 十分感谢！

@wangfei9088 哇，太谢谢了。我也考虑过是单位的问题，但还未在代码中找到数据验证猜想，非常感谢大佬！

在dargmodel.C这个衍生类中没有mesh这个成员对象，不知道该如何引入

@xpqiu 感谢您的答疑
原来如此，是两个box相当于框定了一个水域区范围间的意思。
zoneID场的值设置明白了，还想问一下zoneID这个场是不是因为重叠网格所以要设置的，具体作用是什么呢？

我大部分情况下模拟结果都是这种
上面呢种是因为激波前移叠加了？
不太懂啊
Snipaste_2021-04-02_21-04-11.png

@bestucan 好的，谢谢老师的耐心解答！

@sr0922msc 我没有尝试过运行4.0之前的版本，即使是4.0版本，我在运行tutorials案例的时候，也遇到了一些问题，有些可以运行，有些不可以

因为这不是C:mihu:
这是C++。
你在这里修改，但是其他过程在其他地方修改了。而且这个赋值八成也是不成功的。给类成员赋值一般都由成员函数完成。

改起来还是挺麻烦的。:chitang:

掩耳盗铃法：mkdir那些不存在的文件

偷梁换柱法：把of4里这里的代码复制过来

根治：修改代码中调用of的接口，不要调用这些已经不存在或者改名的模块

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

@fubianhanshu 好的，谢谢您，我去试试

在找资料的过程中发现了关于sampleDict中可用type的说明，包括uniform, face, midPoint, midPointAndFace, curve, cloud（之前在苏军伟前辈的博客中看到过相关介绍，但当时并无curve类型）。

可以看到是有针对曲线的类型(curve)，但是没有给出具体的做法，后面会尝试用这个提取一下曲面上的压力看看。
sampling type1.png

@sungda 是的.而且必须用初始化列表

@Zhoujq 好的，谢谢你的回复，我去调整一下。

@warnerchang
好的谢谢你
@李东岳

gnuplot 配合grep

@mark 今天李老师在课上说，如果进出口速度都设为零梯度的话，进出口的压力好像是都需要给定的

@浪迹天大 好的，谢谢你。

@CFDngu 我也在纠结怎么得到波面数据，请问您知道如何得到了吗？

好的，谢谢老师

@队长别开枪 这个是一篇2016年的硕士论文《基于S-CLSVOF方法的固体火箭发动机装药燃面计算》

1.下图中ReImSum按我的理解，就是将一个复数的实部与虚部相加，为什么要这样做呢？为什么不直接取一个复数的实部？
还有下面fft逆变换中的K及forceGen是什么意思呢？K好像是和Kmesh相关，但是forceGen没找着，希望有大神可以解答一下
2.在piso的修正循环中，rAU是什么意思呢？按我的理解，是不是就是第二张图中的$\frac{1}{ A_{p}}$?

Image 77.png

Image 78.png

@蒙恩 我也遇到了相同情况，请问您解决了吗

@cfd_上弦月 已经解决了，点击左侧的选项''Zone Style''，确定内部流场对应的zone（以上面的图片为例，我这里是Step 1737 Incr 0）。选择''Surfaces''标签，然后在''Surfaces to Plot''处右键选择All surfaces，内部的流体域就能显示出来了。

@李东岳 谢谢东岳老师，我再琢磨琢磨

我觉得你要不就真的试试：
（1）确实做一个循环流化床，从出口引入一个管道到入口，然后直接算；
（2）能不能统计每个时步出口截面出颗粒的位置，速度，然后在入口设置patch颗粒入口

感谢您的回复

@李东岳 好的！我先看看这个有没有

@李东岳 确实如李老师所说，将固定壁面边界条件Omega的值设置为1e-6后，能算了！！！

谢谢东岳哥

@李东岳 噢噢，好吧，可能是我浏览器的问题，不过符号函数应该是少了等号
Image 54.png

这么刺激的吗，大佬:135: 这个硬盘才买了不到3个月:xinlei:

@xpqiu
但是有些资料里面是有这样的设置的

@bestucan 我决定直接放图上去了，不用GCI了，主要是因为GCI到底应该限制在多少以下没有定论，同时GCI也需要求解结果满足单调性，偏偏我这个满足不了……

你看下你是不是要欧拉和拉格朗日耦合，如果要耦合的话，如果需要考虑颗粒体积的话，2维就不成了。如果颗粒体积可以忽略，如果耦合效应很小，可以不考虑z方向。

感谢老师的指点，我自己去试一试！！:xinxin: :xinxin:

@bestucan 谢谢老师，我已经解决拉，是有一些物理特性文件没拷贝。