第八讲应用2ApplicationII.pptVIP

上海大学——冶金工程专业本科生课程 吴永全 * 冶金过程数值模拟Numerical Modelling of Metallurgical Processing 主讲：吴永全 上海大学，材料科学与工程学院，材料工程系，冶金工程教研室 本科生课程——冶金过程数值模拟 * 冶金数值—— 授课内容 授课内容 绪论 数值求解方法 冶金过程数值模拟 数学描述 绪论 数值求解方法 冶金过程数值模拟 数学描述 导热问题数值模拟 对流与扩散数值模拟 导热问题数值模拟 对流与扩散数值模拟 流场计算简介 * 冶金数值—— 对流与换热—— 目录 一维稳态对流与扩散 1 对流项的其他离散格式 2 多维对流和扩散 3 虚假扩散 4 一维稳态对流与扩散 1 * 冶金数值—— 对流与换热—— 一维稳态对流与扩散 积累项 对流项 扩散项 源项 对流产生于流体流动，本次课的任务就是在已知流场（即速度和密度）的情况下，求得φ的解，即求由于对流和扩散同时存在时某自变量φ（如温度或浓度）的解。计算流场部分最后表达。对流项和扩散项之间具有不可分割的联系，因此需要把这两项处理成一个单位。 从数学角度，对流项不过是一阶导数项，其离散毫无问题。但从物理角度来看，这是最难离散的导数项。这主要与对流作用带有强烈的方向性有关。对流项离散方程的构建是否合适强烈地影响到数值解的准确性、稳定性和经济性。 * 冶金数值—— 对流与换热—— 一维稳态对流与扩散 考虑无源一维稳态问题，控制方程为： 其中，u代表x方向的速度。同时连续方程为： 该式对可压缩流体和不可压缩流体都成立。换言之，单位时间单位面积上流过的质量流量为常数。当流体为不可压缩流体时，ρ=const.，则有各断面上速度处处相等。 区域离散如图。 针对如图所示的控制体定解区域对控制方程积分： P(i) W(i-1) E(i+1) (Δx)w (Δx)e (Δx) w e uw ue * 冶金数值—— 对流与换热—— 一维稳态对流与扩散 对于物理量φ的一阶导数用中心差商代替，于是有 式中，系数1/2为内插因子，表示假设界面处于中心。对于不同的界面位置则要采用其它的内插因子。把上式代入顶式得到： 为了把方程写的更简洁，我们定义两个新的符号： 式中，F称为对流（或流动）强度，其物理意义为单位时间单位面积上流过的流体质量，F越大，流过的质量越多，对流强度越大。F值可正可负，仅取决于流动的方向。D称为扩散传导系数（diffusion conductance），D值具有与F相同的量纲，表示扩散状况，但D值永远为正。当Γ=0时，D值也等于0。 * 冶金数值—— 对流与换热—— 一维稳态对流与扩散 这两个变量可以组成贝克来数Pe，它是以网格间距定义的无量纲特征数： 当Re=0，意味着没有对流只有扩散；Pe=∞，表示扩散作用可以忽略，仅有对流；Pe=1，表示对流和扩散的作用相当；当Pe1，表示对流作用远大于扩散作用。 引入上述两个变量后，控制方程可以离散表达为： 一个显然的结果是，对流的引入并未改变最终的离散化方程形式，只是使系数计算和导热问题有些差别，所以前面涉及的离散化方程的求解方法依然适用。 * 冶金数值—— 对流与换热—— 一维稳态对流与扩散 问题的引出 连续方程ρu=const表示Fe=Fw，这时aP=aE+aW，说明满足连续方程就确保了系数之和规则成立。但在迭代过程中，中间过程的迭代值可能不保证满足连续方程，所以上式中不能直接省掉Fe和Fw。根据差分格式稳定性的4个基本准则之一：各系数为正，发现必须满足D-(1/2)F≥0，也就是D≥|F/2|，或者|Pe|≤2。 也就是说，要用中心差分格式，就要保证上面的贝克来数关系，只能选取较小的网格间距，同时还违背了斯卡巴勒准则。 另外，中心差分还不能处理Γ=0的情况，因为满足连续方程后，aP=0，导致方程无法进行计算。 寻求其它方法或者离散格式 * 冶金数值—— 对流与换热—— 目录 一维稳态对流与扩散 1 对流项的其他离散格式 2 多维对流和扩散 3 虚假扩散 4 一维稳态对流与扩散 1 对流项的其他离散格式 2 * 冶金数值—— 对流与换热—— 对流项的其他离散格式 对流扩散问题的严格解（解析解） 对于对流强度和扩散系数都已知的对流扩散问题，我们可以得到严格解。如求解区域为0≤x≤L，边界条件为 其严格解为 其中，这里的Pe=(ρuL)/Γ。 由图中的精确解可以知道，当Pe=0，即纯扩散时，φ与x的关系式是线性的；对于小Pe而言，变化关系偏离线性不大。因此，在这种情况下采用中心差分格式是可行的。但随着|Pe|的增加，节点间的φ值越来越为上游节点所影响，当|Pe|1时，即强制对流时，节点间的大部分区域的φ值几乎就是上游节点的φ值。显然，这个时候如果再用线性分布来近似就会使所得离散方程违反正