这些国内软件包都是直接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贴出来看一下

@李东岳 好滴，谢谢李老师~

发现这个参数匹配好麻烦啊，层数、第一层高度等等之类的，感觉按照计算得到的画一下基本不出来

:136:

CFDEM3.8及OF5.X是用这个reconstructPar -noLagrangian，或者直接在solver中把postpro改成true

@liujunCFD 最后用matlab怎么解决的呢，我目前也遇到这个问题了，没人带，只能自己摸索，所以大佬可以告知一下方法吗？

@李东岳 老师，我有一个问题，压力梯度是指沿着流体流动的方向压力的变化，如果我的流体入口与水平面有45°的夹角，但我的水是沿着x方向进入的，即只有x方向的速度，我的压力梯度应该是把x方向的压力梯度转换到面法向方向上去么？