各位老师，早上好！小白最近在做openfoam燃烧时发现得到的模拟结果过小，比如我在化学反应机理中添加激发态自由基的反应，模拟得到的激发态自由基的结果量级在10^（-9），这在实验中结果是可信的，因为这种自由基不稳定，很快会跃迁到基态，但是从数值角度怎么解释呢？也有一些基态自由基比如CH,它的量级在10^(-5)。在openfoam中很多方程，它的残差到10^（-6），我们就认为它收敛了。10^(-9)量级的值我们怎么才能够保证其数值模拟结果的可信性呢，而不是计算得到的误差?是对一些不大可能发生的反应乘上了一个很小的系数吗?还是经过什么处理?文献中也只是给出了一个结果，但并没有对结果的可信性作出解释，恳请各位老师解惑。

@Prometheus10 在 OpenFOAM计算圆球绕流过程中，如何输出切向粘性阻力系数和法向粘性阻力系数？或是如何从输出的结果中计算得到？ 中说：

我尝试手动在结果中将Symmetry边界条件改为Wall，然后执行一次wallshearstress；
通过对比文献中气泡的压差阻力系数、法向粘性应力和切向粘性应力结果，输出的结果并不对

是的，因为wallshearstress代码会自动过滤非wall边界，可以手动改为wall。既然计算剪切力的公式都是一样的，这么改一下不会影响计算正确性的。和文献对不上是不是symmetry边界条件的问题呢？（你们这个帖子挺长，我看大意是不是：楔形网格算球noslip正确，但算slip和symmetry和文献对不上，且symmetry效果更好些？）
用整球网格算会如何呢？做一些其它情况的验证呢，比如壁面滑移的方腔流动？

我之前做过低雷诺数下的一阶maxwell速度滑移的圆柱绕流，算出来的受力以及粘性正应力占比和文献比得上。不过你的自由滑移情况我没考虑过，不知道为啥出问题。

@Prometheus10 在 OpenFOAM计算圆球绕流过程中，如何输出切向粘性阻力系数和法向粘性阻力系数？或是如何从输出的结果中计算得到？ 中说：

除此以外，是否只有像您提到的【方案二】基于wallshearstress functionObject进行修改，才能输出“绕流气泡”的粘性剪切力结果。

这个只是方便些，直接输出壁面切向和法向粘性力场，不用在paraview中写公式，但是需要改OpenFOAM代码。得先看明白wallshearstress中咋算的，主要是改calcShearStress函数，不难，慢慢看慢慢改。

能量方程收敛非常慢，导致计算时间增加了数十倍

从你的log来看，能量方程都没有计算了。不是能量方程的问题

这些国内软件包都是直接fork国外github上面的OpenFOAM。因为有时候github访问不上去。就fork到国内。

下列地址可直接使用git clone，例如，打开终端，在终端直接将下面的复制进去：

git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-6.git

即可在本地创建ThirdParty-6文件夹。如果提示你没有git，那么输入下面的即可：

sudo apt install git

正常的话是下图这样：

捕获.JPG

相关git地址： git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2312.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2306.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2206.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2106.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2012.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2006.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v1912.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v1906.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v1712.git

第三方包地址：

git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2306.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2212.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2206.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2106.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2012.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2006.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v1906.git

在下载的时候，需要下载相应的OpenFOAM与相应版本的第三方包

我觉得你这个挺有意思的。初学可以练手用。但是个人觉得算法应该没有问题，出现不一致的情况可能还是哪里出现了小bug。

@BlancheDH 好嘞，我去试试，谢谢大佬:xinxin:

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

@李东岳 我用的v2006，Phi是速度势，用-grad(Phi)计算的速度和输出的U有些差异，尤其是在远场边界。

@YuPeng 我也遇到了相同的问题，请问一下您现在解决了吗

感谢李老师！经过修改，已经可以成功运行了。:146:

@李东岳 就是单核跑也不行，现在解决了，应该是openfoam10里面的pitzDaily哪里要调一下，换成其他算例就能跑了，谢谢东哥

root文件夹是不是没空间了？有一些系统root单独空间划片了。我看你跑在了root下。操作不是很常规

这个问题解决了，主要是因为OpenFOAM在solarCalculator.C中有出纰漏
修改如下并编译就好了，可以模拟在useExternalBeam下给定太阳辐射值（随时间变化）的效果

void Foam::solarCalculator::correctSunDirection() { if (sunDirectionModel_ == mSunDirTracking) { calculateBetaTheta(); calculateSunDirection(); //directSolarRad_ = A_/exp(B_/sin(max(beta_, ROOTVSMALL))); if (sunLoadModel_ == mSunLoadTimeDependent) { directSolarRad_ = directSolarRads_->value(mesh_.time().value()); } else if (sunLoadModel_ == mSunLoadFairWeatherConditions) { directSolarRad_ = A_/exp(B_/sin(max(beta_, ROOTVSMALL))); } } } void Foam::solarCalculator::correctDirectSolarRad() { if (sunLoadModel_ == mSunLoadTimeDependent) { directSolarRad_ = directSolarRads_->value(mesh_.time().value()); } else if (sunLoadModel_ == mSunLoadFairWeatherConditions) { directSolarRad_ = A_/exp(B_/sin(max(beta_, ROOTVSMALL))); } }

@李东岳 好的，李老师，那我自己再修改修改定义步长的代码，感谢李老师，有什么参考的资料吗qaq

半圆柱没有

@李东岳 谢谢李老师的回答，我再找找原因！

@沐风 你上传的作者联系我删除你上传的pdf文件，我将在明天将其进行删除

那个就是发散了。

把log贴出来看一下

@zhangK 客气，祝顺利！