[转载]一个问题

zhaoliqie

在导热问题中，我们采用了附加源项法对内点法形成的网格进行边界条件的处理，在哪
里我们将边界点处的gam＝0（二，三类边界），这是合理的，
因为这时没有流动，gam＝0与aE＝0是等价的，但是在有流动的情况下好像不是这样：
对于一维稳态流动和扩散问题，有下式
Je－Jw＝（Sc＋SpYp）dx （Y为通用变量phi）
其中， J为对流和扩散的合并项。如果假设在E边界处的Je＝J已知，则有
Je-aW（Yw－Yp）＝（Sc＋SpYp）dx
（aW－Spdx）Yp＝aW*Yw+(Sc-Je/dx)*dx (Je为流出控制体，取为负）
这样，为了使得aE＝0，令gam＝0，将使得De＝0，而这里完全可能有Fe<0，而aE=De*A(
|pe|)+max(0,-Fe)不为0（如：Y 为混合物中的一个组份分数，该组份的J由主体对流的一部分和扩散组成，这时完全可能由于其扩散作用大于对流，而使得J>0，Fe<0），即 这样将会使得边界点仍对内点有影响。
在陶老师分析出流边界的速度时，将相邻的速度都提升到一个正水平，也是基于这个考
虑，使得Fe>0,这样gam＝0,就可以保证aE为0），
请问在此如何采用附加源项法对上述的边界进行处理呢，当然要保持程序的通用性，即只是在gamsor中对gam进行处理？