@Cr 十分感谢，问题已经解决了:xinxin:

@Cr 原来是在这里啊，感谢感谢！

好的，谢谢您

@Cr 你好 我是把消波边界调整为2D消波边界了 然后可以计算了。

但是提取出来的坝前波面结果很奇怪，我的坝体振幅在三个方向均只有0.4m（初始水深设置的是50m），而水深却会突然变成0。提取出的波面时程曲线如下，暂时没有找到是什么原因导致的这种现象，还请各位老师指教。微信图片_20201009170300.png

@hotairballon1997 和steady Rans结果不同，因为每个时间步的结果不同。

看来还是有只用动网格而不用自适应的。这种状况应该用在变形的铁板、高层建筑随风摆之类的在一定范围内变形的。这样网格即使没有自适应也没关系。一般用来做应力分析之类的。

查看了下网格发现网格变形为负了

你的这个初始网格就不好。在圆形边界上没有用O形网格。畸变本就很严重，然后一变形，还不是震荡变形，是一直一个方向走，网格又没有自适应，然后就负网格了。

我同时想在INLET（火箭喷管口）设置燃气射流环境，目前改了 0/p 0/U 0/T 下的边界条件，并在constant/thermophysicalProperties 设置了燃气参数，但不知道如何指定在INLET上不断喷射燃气。

最简单的把边界设定成固定压力梯度。要是还要计算燃烧过程就不懂了:mihu:

@小考拉

不忙不忙，就是不太懂:mihu:
但是探索不懂的东西的方法都是相通的。你要解决的问题我没解决过，所以要解决花费的时间比较多。思路就是看代码。你也在看，但我们看的版本对不上，开源软件有个特点，太自由了，所以版本间的差距特别大。

在你的代码中找相关定义的方法：

grep -R 'volTensorField' ./

先用 cd 命令 切换到源代码所在的目录 ~/OpenFOAM/OpenFOAM-version/
然后运行上面的那条命令。-R 是递归查找，包括子文件夹，'volTensorField' 是查找的内容，区别大小写， ./代表当前路径。

另外，在 ~/OpenFOAM/OpenFOAM-version/里查，查的太多。可以分别去这个文件夹下的src或者applications下面查。如果还是太多，就把输出重定向到一个文件再看。不过没颜色不太好看。

grep -R 'volTensorField' ./ > grep.log

Thanks. I created two meshes with the same body size but different type of meshes,namely structured and unstructured, the results are different. Note that both of two cases are grid independent. Do you have any idea about that ? since I saw a thread you opened in cfd online.

您好，请问您最后解决了嘛，怎么样输出每个基元反应的速率呀

@King_RVM 没错的问题就在乱用了tmp指针，我已经解决了。关于tmp指针，我在网上搜索到过问题，祝好！

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

请教大家：
提取场内一个面的物理量 Function Obj 的方法
谢谢

另外：
本人原帖 “湍流模型问题请教” 得各位相助， @闻久STU 君 "回访"，李博 @东岳 赐"西语句"，为不在原贴用回谢灌水，在此一并谢过

谢谢，学习到啦！

这个求解器是不可压缩的求解器，他单独求解温度方程，没有用用到热物理模型库，所以就不支持这个边界条件，fluidThermo是热物理模型库的方法。
你可以去看看buoyantPimpleFoam或buoyantSimpleFoam下面的算例，这些算例的constant文件夹下都有thermophysicalProperties这个字典，这些求解器也都用到了热物理模型库求解焓方程，并且支持这个边界条件。

或许可以将代码改为用不可压缩求解器求解的边界上的温度和导热率来获得，边界上的传热量。例如：

if(db().foundObject<volScalarField>("T")) { const volScalarField& T = db().lookupObject<volScalarField>("T"); const fvPatchScalarField& Tp = T.boundaryField()[patchi]; qDot = kappa*Tp.snGrad(); }

其中kappa需要根据你的物性和边界上的温度场来计算获得

@hotairballon1997 算吧 不过不是map，我个人觉得不要轻易用map这个人东西。。因为网格密度不同，map出来的场可能是有问题的

亲测实用:140:

已经解决。我的makeTurbulenceModel.C文件中namespace Foam的
typedef RASModel<EddyDiffusivity<turbulenceModel>> RASfluidThermoCompressibleTurbulenceModel;
typedef LESModel<EddyDiffusivity<turbulenceModel>>
LESfluidThermoCompressibleTurb
应该改成
typedef RASModel<EddyDiffusivity<fluidThermoCompressibleTurbulenceModel>> RASfluidThermoCompressibleTurbulenceModel;
typedef LESModel<EddyDiffusivity<fluidThermoCompressibleTurbulenceModel>> LESfluidThermoCompressibleTurb

@Samuel-Tu 好的，多谢提醒~cfdonline有这两个边界条件差异的讨论。
https://www.cfd-online.com/Forums/openfoam-programming-development/180412-slip-symmetry-boundary-conditions.html

@zhangxc0223 我不清楚你怎么做的，不过我刚发现可以在paraview中可以显示粒子ID，随着时间步推移，好像每个粒子的ID号是不变的，这样就可以通过后处理输出某个固定ID的各个时间步的信息了，可以交流下