CFD中文网

@a602657563 嗯嗯，好的，谢谢了。我的邮箱是your-we@163.com，我这里有一些太阳辐射方位及强度的开源程序，也可以共享一下。

@a602657563
你好，我最近也在做相关内容，一起探讨一下

ViewFactor模型仅在内部流体介质为透明介质（对辐射而言）、所有表面均为漫灰表面时才可采用的辐射模型。
nFacesInCoarsestLevel的意思是“生成的聚合表面数目的最大值”，也就是说，当网格面聚合后生成的用于角系数求解的“粗”网格数目大于或等于nFacesInCoarsestLevel时，会继续对网格进行聚合；
featureAngle：是指网格聚合过程中，准备参与聚合的两个面的法相向量交角的最大值，若大于这个角度，两个面就不参与聚合了。
所以，如果设置中nFacesInCoarsestLevel和featureAngle如果都太小的话，很可能就没法完成面聚合的计算了，举个例子，一个正方体面网格，如果nFacesInCoarsestLevel=5，featureAngle=10，肯定是没办法完成面聚合的。

前面三个参数（writeViewFactorMatrix 、writeFacesAgglomeration、writePatchViewFactors）为聚合过程中相应结果是否写成文件的开关，对计算没什么影响。

看到你在做太阳辐射模型的开发，我想请教一下，Open FOAM里有没有“给定时间、位置，求解当地太阳辐射强度和入射角度”的模块？先谢过了！

@张某人
我猜你的问题可能是边界附近的内部界面怎么处理，这个问题在陶老师的第二版《数值传热学》第166页，5.6.1.4中有说明。
但是这个界面很少用到，一般是直接用边界上的值了。

if  (potentialFlow.finalNonOrthogonalIter())
{
     phi -= PhiEqn.flux();
}

这段我觉得倒不是因为phi被初始化了，好像上面求解的Phi方程就不是 http://dyfluid.com/potentialFoam.html 里头的方程（6），而应该是
▽(▽P)=▽(U-▽P)
也就是代码中的phi=(U-▽P)，这样子的话也能解释速度边界条件是怎么添加到Phi 方程中的了。但有一点不很明白，就是代码中好像表达的并不完全是phi=(U-▽P)的意思，所以还是有些糊涂。

另，

MRF.makeAbsolute(phi);//此处的makeAbsolute按字面意思是取绝对值，那这个对下文中重构速度场不会产生影响么？

这句话应该不是取绝对值，可能用了多重坐标系，之前有个变换，现在变回来了，但不是很确定是不是这个意思
@东岳 @张某人
上面是我的理解，如果有不对的地方还请直接指出来，一起讨论讨论