如何给定网格内部面上的流动参数？

[我是河滩]

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

[李东岳]

1. 可以对Matrix直接设定cell上的值。类似：

   TEqn.setValues(Cells(), 0.0);
   

2. 可以给方程做一个比较大的源项。

[我是河滩]

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

[D.Benjamin]

@我是河滩 这个图像反应堆，你是做堆物理或热工的？

[我是河滩]

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

[李东岳]

@我是河滩 算啥呢老哥，感觉你啥都搞的样子，很神奇

[我是河滩]

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

[xiaolin]

@我是河滩 老哥请问你这块解决了吗，如何控制网格面通量值为0以达到壁面效果。

[我是河滩]

@xiaolin 没有，给定网格单元值可行。

[xiaolin]

@我是河滩 那这样在与相邻的计算单元交界面还是会有通量进出吧

[我是河滩]

@xiaolin 是的