CFD中文网

我想把图片中的红色框的一部分输出，然后在paraview中做后处理，为了方便，我建立了一个同样的参数（黄色框），想用类似湍流模型里边的方法试一下，但是写不出来b910a92ec6b2db806f3cf7ed6ac073a.png
我还尝试了另一种方法：在createfield.h里面定义了一个参数，这样会写出一个文件，但是写出的文件不会更新。55cdd8962f572210b4507aa805be817.png
其中m是自定义的一个参数453c8ef0-3b5a-481d-90df-3fb52ad4320e-image.png

@bestucan 谢谢，这个问题解决了，需要设置UDS的个数，不论有没有标量方程

箭头是对指针的成员的引用，找到这个类的定义：turbulenceModel.H

epsilon()函数声明在205行，是个虚函数，所以在 turbulenceModel.C 里没有定义，而是在使用其作为模板创建新类的时候再各自定义

在上面那个窗口里搜 epsilon() 函数：
Screenshot from 2020-09-15 16-37-07.png

看你的情况是用的哪种定义，其实里面 epsilon() 函数定义大差不差，
你用的 LES，咋找上层流的了，你找的那个确实是层流的对于 epsilon() 的一个实现，在我发的那个图片里也有

脉动速度都用上了，那应该不是层流了，速度分成脉动速度和平均速度这就是用了湍流模型了

@Stan同学 别邮箱请教了，发布在论坛估计能更快得到解答，祝好

目前 paraview 没有这个功能。

我也遇到这个问题,版本v1912,调用了表面采样,截面穿过了AMI,计算的时候把界面关了,reconstructPar成功了

@一二 都是用matlab做的

最近计算一个方腔流动遇到一些问题，这个流动沿流向（x方向）是周期性边界，yz方向是壁面，流动固定x方向的平均速度，自动求解压力梯度，在采用rans计算时是正常的，可以计算出一个合理的结果

pressure gradient = 0.179097

但是当采用PANS（自己编译的湍流模型）计算的时候，压力梯度会变得非常小，导致计算结果很差，

pressure gradient = 0.0903267

自己编译的模型可以在channel395中正常计算，但是在该算例中出现压力的降低，不知道是什么原因引起的

关于OpenFOAM中模化的量与直接求解的量有一些疑问。在OpenFOAM后处理中会有UPrimeMean和R的出现，这两个量分别代表的是什么的意思？在CFDonline中提到，R给出的是模化的雷诺应力，UPrime2Mean给出的是直接求解的雷诺应力，这种说法正确吗？如果是正确的，那总的雷诺应力就是R+UPrime2Mean？

R gives the modelled Reynolds stress tensor and that Uprime2Mean gives the resolved Reynolds stress tensor.

附上CFDonline链接链接文本

一年多了还是没有解决这个算例，不知道有没有人做过相类似的，重新描述下这个算例，如图所示一个后台阶流动，但是是在图中左下（0,0,0）位置为原点绕Z轴旋转，主要流动参数，
雷诺数$Re_b=\frac{U_b(H-h)}{\nu}=5600$，其中$U_b$是体积平均速度；
旋转数$U_b=\frac{\Omega(H-h)}{U_b}$，$\Omega$是角速度；
台阶前段槽道流部分壁面剪切雷诺数$Re_\tau=\frac{u_\tau(H-h)}{\nu}=360$，这个与槽道流DNS结果是完全一致的，也就是说台阶前是完全发展的槽道流。目前存在以下问题：
1、求解器选择，SRF类还是需要自己修改求解器？
2、合理的进口条件，不旋转的相对好给，但是旋转的case，速度进口边界条件和初始值不知道如何给合理。
3、无量纲化的流向速度为什么会小于1，如下图所示。以我的经验，槽道流里面最大速度比体积平均速度肯定是大的，也就是$\frac{U_{max}}{U_b}$是大于1的，但是为什么在旋转流动里面就会于1呢？
6c444575-46fa-460b-bb57-8d655a396d45-image.png
参考文献
1、DNS原始文献：Turbulent flow over a backward-facing step. Part 1. Effects of anti-cyclonic system rotation
2、图片来自：Development and Application of Rotation and Curvature Correction，是篇博士论文

@fangyuanaza 我也遇到这个问题了，貌似是工具里面判定网格中心点距离的问题，太小的网格就提取不到值

对于稳态问题的收敛，我认为其判定可从两个方面进行，一是根据残差收敛到一定标准如10e-5，二是监测点的速度或者压力是呈现周期性变化的；
目前我计算一个较为复杂的泵的模型，计算采用simpleFoam，湍流模型采用kepsilon，计算过程中，可以实现残差收敛到一定的标准，但是对于检测点，压力速度湍动能等都是呈现上升的趋势。从残差角度，该计算符合我的收敛标准，但是对于检测点却没有出现周期性变化。
我的问题是，这种情况下，该结果可信吗，计算算是手收敛吗？
如果不能算收敛的话，这种情况可能是哪里的原因？该如何改进？

bounded不建议用于瞬态计算哦

是否是真实的流场中湍流长度尺度呢？

Kolmogorov是一种假说：所有雷诺数足够高的湍流中，最小的尺度由一个普适性方程来决定，其参数为粘度和耗散率。这个假说建立在若干的假定基础之上，如各向同性。随后，有了Kolmogorov尺度。所以针对你的问题，至少各大教材都认为Kolmogorov尺度是最小尺度。

这个泰勒湍流尺度能否代表RANS求解中的尺度信息呢？

是的

在大涡模拟中的湍流长度如何让计算？

LES里面的湍流长度通过$\Delta$控制。大于$\Delta$被解析，小于$\Delta$被模化。不确定是否可以通过$\Delta$定义LES的湍流长度。

湍流模型kOmegaSST是低雷诺数模型？在近壁区采用k-omega模型，在湍流区采用k-epsilon模型。
那在计算的时候应该采用什么样的边界条件？
我的case是这样的，网格较密，yplus在1附近，雷诺数为7000，存在两个壁面(wall)
我参考的是researchgate中学者的设定方法，具体方法见链接(https://www.researchgate.net/post/Should_the_RANS_model_kOmegaSST_be_fully_resolved_at_low_Re-number_Is_it_possible_to_use_it_with_a_wall_function)
方法有一定合理性，但是对于我的算例可能存在一定的问题，不知大家在应用SST模型时，有没有什么经验，如何设定边界条件？

是否为容差 matchTolerance 太大导致的衔接不上~

@benqing 谢谢，等我研究研究，不懂的再向您请教。

@西湖冷月 对的，这里说硬盘已经是默认是移动硬盘了，移动硬盘就当个大u盘了，拔插就关机太费劲了…

现在正常了

低雷诺数模型(例如v2f，SST k-w等)在精细网格的情况下可以完全求解壁面的流动，那为什么还要壁面函数(比如v2WallFunction、fWallFunction等)？