有限体积法边界条件处理问题

spacewarship

    编写有关用有限体积法三维CFD求解程序,在显示处理边界条件时遇到诸多问题,在此提出来急切希望有哪位大侠能指点迷津:
   1.首先有限体积法是求解的是网格中心点的物理量值,那么边界处的各个物理量是不是一定要给出或者程序也要计算出来,但是比如对称边界处物理量值无法显示给出,那该怎么计算
   2.为了处理边界的方便,有的程序常常在边界外引入虚网格,一般虚网格设置几层,适用于什么边界
   3.在处理无穿透壁面和黎曼远场边界时,需要采用内场外插法,那么边界处的法向和另外两个切向如何给定
   先谢了!

spacewarship

  1.但是要求解内部结点,就得确定边界上的通量,这个时候不也要预先知道边界上的物理量值吗,
  2.即使先求出来内结点的值,再怎么计算边界上的值,比如对称边界
   我对有限体积法也是初学,还请多多指教了,非常谢谢了!

anyone

有限体积法就是面上通量和体积内变化的守恒关系. 边界上的通量不知道得话, 就不封闭了.边界条件常见有3种:
1. 给定值,
2. 给定通量
3. 值和通量混在一起.
给定通量的最容易, 直接把通量放到源项就行了. 对称边界就是给定通量,而且通量为0.

anyone

[这个贴子最后由anyone在 2004/12/15 04:06am 第 1 次编辑]

下面引用由dpshaka在 2004/12/14 02:23pm 发表的内容：
第一类边界不用处理
第二类 和第三类，用附加源项法处理。就统统可在内部接点求借了。

如果第一类边界不处理的话, 一定要把边界点放到离散后的代数方程里. 现在的有限体积法一般不这样做了. 代数方程里只有每个单元的中心点.为了引入边界条件, 第一类边界也要处理. 在有限体积法第一类边界比第二类还麻烦一点, 一般要同时加源项和系数.