@bestucan 感谢您的回复，抱歉之前不知道为什么没看到。好几年过去了，今天刚看到这个之前问的问题，我那个卡住的原因是：（拿好笔记本记好了，这个问题值100万）我刚开始不会用openfoam，然后硬盘满了，所以算不下去了。

@coolhhh 在 有人有用过snappyHexMesh生成特别薄的边界层吗 中说：

@tyhero123 不太清楚无法加上边界层具体类型。我之前是画方柱绕流的网格，方柱的边界层要薄，计算域底面厚点，方柱拐角处、方柱与计算域底面衔接处的边界层无法生成问题。通过下面方法实现生成：

snappyHexMeshDict中用explicitFeatureSnap explicitFeatureSnap true; implicitFeatureSnap false; meshQualityDict中的maxBoundarySkewness设置大点，几个控制网格质量参数设置小点 maxBoundarySkewness 10; maxInternalSkewness 4; minDeterminant 1e-30; minFaceWeight 1e-30; minVolRatio 1e-30;

步骤1和2只是解决方柱拐角处边界层衔接，计算域底面边界层厚度和方柱的边界层厚度设置一样薄，方柱与计算域底面衔接处的边界层才能生成

生成很薄的边界层网格质量其实很差

好的谢谢！我先试试看能生成怎么样的网格。

直接输入paraFoam能打开吗？我以前没装好的时候，提示要输入paraFoam -builtin才能打开软件

用moveMesh程序，可以旋转

tmp<fv::convectionScheme<scalar> > mvConvection ( fv::convectionScheme<scalar>::New ( mesh, fields, phi, mesh.divScheme("div(phi,Yi_h)") ) ); { reaction->correct(); dQ = reaction->dQ(); label inertIndex = -1; volScalarField Yt(0.0*Y[0]); forAll(Y, i) { if (Y[i].name() != inertSpecie) { volScalarField& Yi = Y[i]; fvScalarMatrix YiEqn ( fvm::ddt(rho, Yi) + mvConvection->fvmDiv(phi, Yi)

这个看起来并不会导致区别。有可能cfd-online的用户没做特别仔细的判定吧。

你这个进出口边界有点问题 最后能使用固定值速度进口 固定值压力出口测试一下 瘦脸之后再调试别的边界条件
另外我看你的网格 最好网格质量能提高一下 进口那个地方看起来很奇怪

谢谢老师解决了，原因是使用了不同的版本of来回切换导致的错误，使用同一版本执行命令就没问题了

代码方面你照猫画虎就行呀，它这个是给scalar场添加源项，你给速度场U添加就行了，该变的地方变一变。我用的可能和这个不一样。

codedSource { type coded; selectionMode all; fields (h); //对哪一个场添加源项，你用 U name sourceTime; //名字随便起一个 codeAddSup #{ const Time& time = mesh().time(); //获取时间 const scalarField& V = mesh_.V(); //获取网格的体积 scalarField& heSource = eqn.source(); //获取h的方程的源项的引用，之后在这个scalarField上面加上具体的量就可以了 // Start time const scalar startTime = 2.0; // Retrieve the x component of the cell centres const scalarField& cellx = mesh_.C().component(0);//获取网格的x坐标 // Only apply when we have reached the start time if (time.value() > startTime) { // Apply the source forAll(cellx, i) { // cell volume specific source heSource[i] += 1e5*sin(200*cellx[i])*V[i]; //forAll所有网格的循环，加上这个源项 }; } #}; }

@dxl 在 reconstructPar合并问题 中说：

并行计算结束之后先用的

你并行计算了什么？snappyhexmesh 还是求解器？

目前想用一个类似sprayFoam的求解器，算一个旋流燃烧的算例。在喷入燃料前需要先通入空气，对冷态流场进行计算。但在冷态流场计算时，出现了部分高温、高压的小块出现，同时时间步变的很慢，无法继续正常计算，不知道各位有没有什么建议和指导:140:
18f84d67-316a-46dd-af4d-c6fcaeab642c-image.png
右边绿色区块为inlet，流入空气，仅有流向速度为20m/s（后方圆柱长度约为20mm）。紫色平面包括叶片为wall。
a72a1134-5beb-4e4c-aa34-a9f1909ab6c0-image.p ng
截取了某时间点的温度分布，出现明显高温小块（同时高压），且后续会持续上涨，无法正常计算。

9366c4d5-ce39-46bc-85e5-8dcadbb12850-image.png
哪位老师可以解释一下这个 slhs 1 srhs 2 glhs 0 grhs 4是什么意思，我左右两边质量已经守恒，开启的是质量比不是摩尔比，为什么还是质量不守恒呢

在foam-extend4.1的ibm里折腾了两三个月，发现用foam-extend4.1的浸没边界法计算升阻力系数产生震荡是其本身的问题。具体原因可以参考这篇论文Sources of spurious force oscillations from an immersed boundary
method for moving-body problems
这篇论文提到了三个解决方法：1、在力源项后面再增加一个质量源项以抑制动边界的力震荡 2、增大时间步长 3、将动边界所处位置的网格加密
目前我只是用了后面两种方法可以有效的减少其震荡，但震荡还是会存在，如果想获得一条平滑曲线，可能需要与第一种方法结合。第一种方法涉及到一些算法问题，鉴于目前的水平有限还不会在openfoam中实现，希望有大佬能够实现并给大伙解惑。

谢谢 学到啦 还有一个问题 采用高分辨率网格可以弥补精度不够高的不足嘛

@da 您好我不是做这方面的，谈一下自己的理解，不一定正确。涡方法感觉和谱合成法一样，都有自己的局限性，可以参考罗银师兄的博士论文。设置没有波动，很可能是随着流动距离的增加，涡流衰减掉了，尝试加密一下网格或者用其余方法试一下。

用paraview打开sparyFoam求解的aachenBomb案例文件，结果有这样的报错：
ERROR: In C:\bbd\8fbebceb\build\superbuild\paraview\src\VTK\IO\Geometry\vtkOpenFOAMReader.cxx, line 7891
vtkOpenFOAMReaderPrivate (0000022DE32F20E0): Error reading line 20 of E:\OpenFOAM\OpenFOAM-8\sprayFoam\aachenBomb\5e-05/lagrangian/sprayCloud/positions: Expected punctuation token ')', found f
我用的是Windows环境下的paraview5.11.0。路过的大佬，帮忙看看。

@李东岳 在 momentumLookupTableInjection这个粒子注入模型一直无法正确注入到相应的位置 中说：

那个帖子很老了，我看是用在openfoam2.0.0的，相当老了
我还没测试过你这么玩，感觉你可以把算例发到论坛试试。不过我马上要做国际航班了。够呛能给你debug了近期。

谢谢！我已经解决了，自己编了一个injection模型。

您好 您后来展向平均怎么做的呢

@李东岳 在 最终残差曲线 中说：

看不太出来纵坐标到了多少。不过CFD残差经常这个样子。按道理来说如果最开始就是从openfoam起来的，都不太会关注残差。老铁是不是有商软的使用经验。

是的，李老师（捂脸）

@李东岳 找到了 https://vimeo.com/148670249

@xuanze 3374769460@qq.com