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; }

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

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

这些辐射系数都是和 wavenumber 相关, 在RT E中需要这些带量纲参数.
Thanks