CFD中文网

如题，根据probes的字母意思解释，是Set of locations to sample（设置位置集），我看了下probes.H和probes.C的源程序，没有看懂，probes.C的部分代码如下：

const label celli = mesh.findCell(location);
  
         elementList_[probei] = celli;
  
         if (celli != -1)
         {
             const labelList& cellFaces = mesh.cells()[celli];
             const vector& cellCentre = mesh.cellCentres()[celli];
             scalar minDistance = great;
             label minFaceID = -1;
             forAll(cellFaces, i)
             {
                 label facei = cellFaces[i];
                 vector dist = mesh.faceCentres()[facei] - cellCentre;//是指面心和体心的插值么？
                 if (mag(dist) < minDistance)
                 {
                     minDistance = mag(dist);
                     minFaceID = facei;
                 }
             }
             faceList_[probei] = minFaceID;
         }

@chengan-wang 里面基于ASHRAE手册中的公式进行的开发。tutorials中的案例能模拟出地面温度随太阳高度角动态变化的过程。你看看先，希望对你有帮助。

@chengan-wang 你好，建议试一下ESI集团的版本OpenFOAM1804及以后的版本由太阳辐射计算模型。

不知道现在是否发布了理论指南的开源文档？

这里有一个，请问你是做什么方向的呢？

建议列表表示下边界条件设置，实在不行可以试试CFX，它有个Openning 出口边界，对不可压缩稳态问题收敛性较好。

你的问题不清楚，这图是从哪来的？initial Gauge Pressure在fluent帮助文档里面说的是给定一个初始压力，而不是实际计算后入口处的压力值，具体你可以查查。另外，我个人理解，入口边界条件只能给定一个，如质量流量、速度、压力。如果能够同时给定质量流量以及压力，那么会不会导致矩阵方程组无解，这块我理解的也不够深。等待大神解答。

ansys meshing越来越傻瓜式了， 没有网格划分基础的话，很多错误就得一遍一遍去试错。

请问这是哪个版本fluent对应的tutorial？

@东岳 2016-Towards-A-Cooler-Singapore.pdf
我看上传文件最大为400k啊。只能上传这个文档了，这是徐明博士早期的文章。如果有兴趣，可以顺着这个文章找他后来的作品。请东岳查收。

谢谢您的解答！目前我做的是流固多场耦合，固体是阀门（上百个零件），流体是阀门及附属管道内区域（开度约10~20°），介质是可压缩空气，目前已知的只有质量流量入口，换算成管道入口平均流速约400~500m/s，若把流体当不可压计算，那么在阀板附近将达到一万多米每秒。如果当做可压缩计算，结果总是发散：我把流体密度设置为idealgas，其他参数恒定，入口用质量流量，入口恒温，压力出口为大气压。目前的发散状态是阀板附近流体温度不断降低到1K以下，超过fluent最低温度限制了。初步怀疑是阀板附近网格问题，当然，这样的设置是否合理，还未经过严格的论证，请问您还有什么建议么？

@youv 好的，难得碰到一起研究OpenFOAM太阳辐射计算模型开发的同仁啊。我这就发你邮箱去，你说的太阳辐射方位及强度开源程序，我也非常感兴趣！我之前用Excel编写了下，正好可以对照下。望共享！

@东岳 据说从ansys17.0以后推出的fluent meshing划分高级网格的功能很强大，我还用着15.0.改天得空试试。另外，前处理模块用的SCDM，据说比之前的Designmodeler模块强大很多。此外，ansys收购了Ensign后处理器。家族越来越庞大了啊。

今天偶然与安世亚太流体分析工程师谈起ansys fluent的tutorials对应原始几何及网格文件的问题，之前一直以为ansys大厂免费发布帮助文档及几何/网格文件，后来才得知需购买正版版权才能拿到。最近项目碰到可压缩气体通过突缩突扩管而超音速，温度方程始终发散的问题，决定从tutorials开始研究管内高超音速可压缩流，奈何竟然没找到几何、网格文件。今天冒昧在此发帖，希望有ansys fluent15.0的tutorials网格文件或者cas、dat文件的同仁能支个招上哪找，或者共享下。先谢谢啦！

谢谢回复！你的解释很清楚了！根据我的调查，OpenFOAM中没有你所说的“给定时间、位置，求解当地太阳辐射强度和入射角度”的模块，之前见新加坡国立大学-未来城市研究中心的徐明（华东理工大学本科）博士开发过相关模块，计算过城市热环境。但是后来据说他去搞金融去了。他发表的一些文献我有下载，如果你也打算开发太阳辐射计算模型，我们可以交流下，资料可以分享给你。

希望有更多的同仁一起讨论辐射模型的问题。

谢谢东岳！你推荐的这篇文章我迟迟没有下载下来，未及时回复，见谅！目前，对该文件中的字典有了初步认识，辐射方程可以采用三种方法：P1模型、fvDOM模型和ViewFactor模型，其中ViewFactor模型即角系数模型，即建立面与面之间几何和位置对应关系的数学表达，为了避免每个网格都要建立关系而导致的计算量过大的问题，需要将一个面上紧邻的网格形成表面束，新的表面束需要指定聚集的网格数目（nFacesInCoarsestLevel），而新形成的表面束也会有新的中心，从而与其他各个表面束中心形成一个特征角度（featureAngle），至于这里的特征角度是上限还是下限或者其他什么意思就不得而知了。

最近在从事太阳辐射模型的开发，涉及到一个chtMultiReigonSimpleFoam的算例，我用的是OF3.0.1版本，案例里面constant/bottomAir/viewFactorDict中有一段这样的文本：

writeViewFactorMatrix     true;
writeFacesAgglomeration   true;
writePatchViewFactors     false;

bottomAir_to_heater
{
    nFacesInCoarsestLevel     40;
    featureAngle              24;
}

我有尝试过设置新的模型的nFacesInCoarsestLevel 和featureAngle值，如果设置不当，在生成角系数时无法通过，目前是随便试试出来的，现在想知道这两个参数的具体含义以及前面三个参数（writeViewFactorMatrix 、writeFacesAgglomeration、writePatchViewFactors）的用法。在CFD-online上也没有过多的解释。请问哪位同仁能解答下困惑么？感激不尽！

谢谢东岳的回答，您发的是热辐射的一般理论，可能在OpenFOAM中也是依据“吸收率+反射率+透射率=1 ”的原理进行设置的，但是问题的关键是无法确定constant-radiationProperties中的字段是如何与以上三个参数对应的。此外，我还注意到，在tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation/constant/heater/radiationProperties中有这样的字段：

radiation       on;
radiationModel  opaqueSolid;
absorptionEmissionModel constantAbsorptionEmission;
constantAbsorptionEmissionCoeffs
{
    absorptivity    absorptivity [ 0 -1 0 0 0 0 0 ] 0.0;  //opaque
    emissivity      emissivity [ 0 -1 0 0 0 0 0 ] 0.1;
    E               E [ 1 -1 -3 0 0 0 0 ] 0;
}
scatterModel    none;
sootModel none;

这里，解释说absorptivity指定的参数为0时表示opaque（不透明），但是若absorptivity翻译为“吸收率”，那固体不透明介质吸收率不可能为0。因此，这里很矛盾。此外，这里的absorptivity的单位为m-1，教科书中吸收率应该是一个无量纲数。因此比较费解。

此外，值得一提的是，该案例system/bottomAir中，与辐射相关的边界条件如Qr、G、IDefault等有默认壁面边界的emissionMode lookup； emissivity uniform 1.0；这里与constant中的定义的emissivity又有什么区别呢?

cfd-online上的同仁也有这样的困惑。我贴下链接，以供参考：
https://www.cfd-online.com/Forums/openfoam-programming-development/135502-understanding-marshak-boundary-condition-radiation.html

案例位置：tutorials/heatTransfer/hotRoomRadiation/constant/radiationProperties中的

constantAbsorptionEmissionCoeffs
{
    absorptivity    absorptivity    [0 -1 0 0 0 0 0] 0.5;
    emissivity      emissivity      [0 -1 0 0 0 0 0] 0.5;
    E               E               [1 -1 -3 0 0 0 0] 0;
}

这里面emissivity是发射率，但是作为灰体辐射，发射率等于吸收率。那absorptivity又是指代什么呢？我在cfd-online上看到有的同仁说absorptivity其实指代的是透射率，若透射率为0，吸收率/发射率为0.7，那么反射率可以自动确定为0.3。不知道这样理解对不对。（PS：这里absorptivity翻译为中文是“吸收率”的意思，这点着实让人疑惑）