[求助]simple算法中的难点——边界条件的处理

zlf81

    目前在学simple算法，可是老碰到问题，特别是边界条件的处理上，弄不懂。好头痛，希望大家帮帮忙！
    书上说——当为第一类边界条件的问题，由于边界温度为已知值，使求解区域仅限于内部节点。
    也就是说，对于第一个内节点的ap*Tp＝ae*Te+aw*Tw+b来说，Tw已知，这样的话，aw*Tw要不要作为源项的Sc部分？？
    书上说——当为第二、三类边界条件，可以采用附加源项法，但是是处理多维非稳态导热问题，如果多了udT/dx+vdT/dy项，即有对流项，附加源项法，是不是也成立啊？？
    希望大家对于这个问题能多说几句。
    谢谢！！！

zlf81

    有没有人知道啊？
    请大家帮帮忙，怎么处理这几种情况？？
    谢谢！

dongshanyan

你先针对你的问题，确定好使用哪一种边界条件。我个人认为直接使用第一二类边界条件就没必要用附加源项法。另，你对湍流的壁面函数的边界条件搞清楚吗，如果清楚的话可以大家一起讨论

zlf81

我计算固液对流传热的能量方程，使用第三类边界条件q=h(T-Tw)。
我看书上说，第二、三类边界条件要用附加源项法，因为对于第一个内节点的离散方程要用到边界节点的温度值，而该值并不知道，因此把边界条件代入消去边界节点的未知温度Tw，得到Sc，ad和Sp，ad，再令aw＝0（编程时令边界节点的热导率为0即可），从而使方程仅限于内节点的求解，然后根据第一个内节点的值计算边界节点。求得整个温度场。
但是我在编程时，令导热率为0后，aw就变为0，这不是第二类绝热边界条件的处理方法吗，计算的结果也比较接近绝热情况。
应该还是在某个地方的理解上弄错了，但又弄不清楚。
望大家指点一下！！谢谢！！

zlf81

    谢谢！ 明白了！
    但是在学习的时候，发现书上处理的第三类边界条件为q=h(T-Tw)。
    而我要处理的第三类边界条件中还含有辐射项，即
    q=h(T-Tw)+StefanBoltzmann*emissivity*(T^4-Ta^4)

    这样就不知道怎么办了，有T的4次方？
    不知道你知道不？望指教！谢谢

zlf81

我怎么没想到，还是学的不透彻啊！  万分感谢！
还想和你讨论一下：速度边界条件的问题。麻烦了！
附加源项法对于动量方程的速度边界条件应该也是通用的把？
比如：有速度边界条件为du/dy=-(dT/dx)(dr/dt)，这是温度变化引起表面张力引起速度变化。
由于右边项中不含变量u，所以可以把它看成du/dy=q，（类似于热边界条件的第二类边界条件：边界上受热流密度q=kdT/dy的作用），故属于第二类边界条件。可按第二类边界条件得出附加源项Sc,ad。
不知道我的理解对不对，望再次指点！谢谢！^_^