温盐双扩散系统对流扩散的罚函数迎风有限元模拟

对温盐双扩散系统的对流扩散问题采用了罚函数迎风有限元数值方法模拟。该方法具有稳定性好，二阶精度，可直接求解原始方程等优点，并用一个具有横向受热的方腔加以验证。所得定常结果和前人用速度压力法所做非定常结果的趋势很吻合。在此基础上，应用该方法对具有横向受热，设置障碍物的方腔中的温盐双扩散问题数值计算，得到了温度雷诺数Ｒａ＝１０６及热浮力比Ｎ＝１，３，５，７下的流场、温度场和浓度场的情况，考察了Ｎ与障碍物的影响．