发现错误了，是文件名字的问题，文件名本来是大写字母，我写成小写了。警告对结果没影响，或者使用fin函数进行输入，就没有警告了。

@东岳 总压减小是因为速度变小，还是静压变小，或二者都有？无粘不可压的时候就是动能与压力势能的转化，是守恒的。

@东岳 伯努利方程 速度增大，压力变小的结论，就是基于总压不变，分离后就变了。

不可压计算，不求解能量方程，温度不变，产生分离后，分离区的总压要小于外面流场的总压，减小的总压去哪里了？

@东岳 运动物体的计算，cfd的方法都很慢啊。

@D-Benjamin 不是，是个二维圆柱。

@东岳 setValues实现的是给定网格单元中心的值，源项的话我再想想，谢谢。

@yfclark 目前想用笛卡尔网格来进行计算，贴体网格自然好处理。但以后打算计算边界的运动，不涉及边界的变形，就只是平动和转动。Foam extend 4.0中的浸入边界法比较合适，但它采取的是给定网格单元中心值的办法，而不是直接给定内部面上的值。给定网格中心值，相邻网格面上的值还是正常进行插值计算。给定面上的值，就是要直接给定速度，得到phi（对流项）及通量，给定速度梯度，得到面上应力（laplace项扩散项），大致是这样的思路。

@东岳 确是可能不太符合物理，或者说是对于内部面围成的封闭区域，比如一个长方形或者一个圆，其内部的网格单元不参与求解，此时内部面类似于外部面的边界，上面的值就需要给定。但OpenFoam仍然将其视为内部面处理，应该如何给定面上的值啊？

采用笛卡尔网格进行计算，而不是贴体网格，圆外的网格单元与圆的边界为内部面，阶梯状的。

外部面上的值可以通过边界条件给定，内部面上的值一般均是利用相邻网格单元的值进行插值，比如内部面上的层流粘性系数和湍流粘性系数。通量phi的计算需要面上的速度，也是利用相邻网格单元的速度进行插值再计算。如果现在需要人为给定内部面上的流动参数，并进行流动方程的求解，有什么方法可以实现吗？

离散格式，求解方法，湍流边界条件等都可能造成差异，你说的不够详细