Lesson 05 - Application II - Convection and Diffusion

1、上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分本科生课程本科生课程冶金过程数值模拟冶金过程数值模拟上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 授课内容授课内容上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 目录目录上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 一维稳态对流与扩散

2、一维稳态对流与扩散Sxyzxxyyzz uuu对流产生于流体流动，本次课的任务就是在已知流场（即速度和密度）的情况下，求得的解，即求由于对流和扩散同时存在时某自变量（如温度或浓度）的解。计算流场部分最后表达。对流项和扩散项之间具有不可分割的联系，因此需要把这两项处理成一个单位。从数学角度，对流项不过是一阶导数项，其离散毫无问题。但从物理角度来看，这是最难离散的导数项。这主要与对流作用带有强烈的方向性有关。对流项离散方程的构建是否合适强烈地影响到数值解的准确性、稳定性和经济性。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值

3、对流与换热对流与换热 一维稳态对流与扩散一维稳态对流与扩散考虑无源一维稳态问题，控制方程为：uxxx 其中，u代表x方向的速度。同时连续方程为：0 or const.uxu 该式对可压缩流体和不可压缩流体都成立。换言之，单位时间单位面积上流过的质量流量为常数。当流体为不可压缩流体时，=const.，则有各断面上速度处处相等。区域离散如图。针对如图所示的控制体定解区域对控制方程积分：P(i)W(i-1)E(i+1)(x)w(x)e(x)weuwueewewuuxx 上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对

4、流与换热 一维稳态对流与扩散一维稳态对流与扩散ewewuuxx 对于物理量的一阶导数用中心差商代替，于是有11;22eEPwWP式中，系数1/2为内插因子，表示假设界面处于中心。对于不同的界面位置则要采用其它的内插因子。把上式代入顶式得到： 1122eEPwPWEPWPewewuuxx为了把方程写的更简洁，我们定义两个新的符号：;Fu Dx式中，F称为对流（或流动）强度，其物理意义为单位时间单位面积上流过的流体质量，F越大，流过的质量越多，对流强度越大。F值可正可负，仅取决于流动的方向。D称为扩散传导系数（diffusion conductance），D值具有与F相同的量纲，表示扩散状况，但D

5、值永远为正。当=0时，D值也等于0。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 一维稳态对流与扩散一维稳态对流与扩散FuxPeD这两个变量可以组成贝克来数Pe，它是以网格间距定义的无量纲特征数：当Re=0，意味着没有对流只有扩散；Pe=，表示扩散作用可以忽略，仅有对流；Pe=1，表示对流和扩散的作用相当；当Pe1，表示对流作用远大于扩散作用。引入上述两个变量后，控制方程可以离散表达为： 11;22 where 1122EeeWwwPPEEWWPeewwEWewaDF aDFaaaaDFDFaaFF

6、一个显然的结果是，对流的引入并未改变最终的离散化方程形式，只是使系数计算和导热问题有些差别，所以前面涉及的离散化方程的求解方法依然适用。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 一维稳态对流与扩散一维稳态对流与扩散 11;22 where and ;1122EeeWwwPPEEWWPeewwEWewaDF aDFaaaFu DxaDFDFaaFF1.连续方程u=const表示Fe=Fw，这时aP=aE+aW，说明满足连续方程就确保了系数之和规则成立。但在迭代过程中，中间过程的迭代值可能不保证满足

7、连续方程，所以上式中不能直接省掉Fe和Fw。根据差分格式稳定性的4个基本准则之一：各系数为正，发现必须满足D-(1/2)F0，也就是D|F/2|，或者|Pe|2。2.也就是说，要用中心差分格式，就要保证上面的贝克来数关系，只能选取较小的网格间距，同时还违背了斯卡巴勒准则。3.另外，中心差分还不能处理=0的情况，因为满足连续方程后，aP=0，导致方程无法进行计算。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 目录目录上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7

8、月4日5时27分冶金数值冶金数值 对流与换热对流与换热 对流项的其他离散格式对流项的其他离散格式00,;,LxxL对于对流强度和扩散系数都已知的对流扩散问题，我们可以得到严格解。如求解区域为0 xL，边界条件为其严格解为00exp1exp1LxPeLPe其中，这里的Pe=(uL)/。Pe=100Pe=20Pe=2Pe=1Pe=0Pe=-1Pe=-2Pe=-20LLx00Pe=-100由图中的精确解可以知道，当Pe=0，即纯扩散时，与x的关系式是线性的；对于小Pe而言，变化关系偏离线性不大。因此，在这种情况下采用中心差分格式是可行的。但随着|Pe|的增加，节点间的值越来越为上游节点所影响，当|P

9、e|1时，即强制对流时，节点间的大部分区域的值几乎就是上游节点的值。显然，这个时候如果再用线性分布来近似就会使所得离散方程违反正系数规则，从而导致物理上的不真实解。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 对流项的其他离散格式对流项的其他离散格式针对中心差分的缺陷，考虑到大Pe时节点间的大部分区域的值基本上等于上游节点的值，于是在扩散项仍采用中心差分格式的同时，对流项中的e与w值均取上游节点的值，如0,0,0,0,eePwWeeEwPFuFuP(i)W(i-1)E(i+1)(x)w(x)e(x)

10、weuwue于是就有 max,0 ;max,0 where EeeWwwPPEEWWPEWewaDFaDFaaaaaaFF这种格式的系数可以看到不会有负系数问题，从而保证了任何Pe都能得到物理上真实的解。但一个明显缺陷是：无论小Pe还是大Pe，计算误差都较大。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 对流项的其他离散格式对流项的其他离散格式根据严格解的构造可以定义dJudx 将一维无源稳定严格解代入J的定义式，得到00exp1LJFPe对于控制容积面w和e，分别用W和P以及P和E替代0和L，并用

11、hW替代L，于是式中，J表示由对流流量u和扩散流量-d/dx组成的总流量。于是，控制方程可以改写成dJ/dx=0。exp1WPwwWwJFPeexp1PEeePeJFPe代入J的控制方程离散形式得到 ;expexp1exp1 where ewEWwewPPEEWWPEWewFFaaPePePeaaaaaaFF这个方法可以满足正系数规则，也可以保证一维无源稳态问题的严格解。但这个严格解只对一维无源稳态问题成立，对一般的对流和扩散问题（多维、有源等），将失去严格解的含义。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换

12、热对流与换热 对流项的其他离散格式对流项的其他离散格式-5-4-3-2-1012345012345aE/De=0aE/De=-PeaE/DePe指数格式aE/De=1-Pe/2做指数格式中系数aE/De与Pe的函数关系如右图。发现在Pe分别趋向于-、 、0的时候， aE/De分别趋向于-Pe、0、1-(Pe)/2。混合格式实际上构成了对严格解的包络线，它代表了一种合理近似。混合格式就是由这三条线组成，即：2,22,122,0EeeeEeeeEeePeaDPePeaDPePeaD 代入控制方程并整理可以得到 max,2,0 ;max,2,0 where EeeeWwwwPPEEWWPEWewaF

13、 DFaF DFaaaaaaFF与指数格式一样，混合格式没有要求控制容积面一定要在两节点的正中间，可以在两节点间的任意位置。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 对流项的其他离散格式对流项的其他离散格式设场变量经过对流扩散过程从一维区域的x=0点传输到x=L点，如图。流体密度为=1.0kg/m3，L=1.0m，扩散系数=0.1kg/(ms)。求流速分别为0.1m/s和2.5m/s时，将区域离散成10个节点网格时，的分布。x=0 x=L=1=0u此问题的解析解为：00exp1exp1LuxuL

14、 可采用中心差商离散。x=0123i-1 i i+1n-1x=Lxnxx采用外节点法将区域离散成10等分，每个控制容积长度为0.1m，相应的F=u，D=/x，Fe=Fw=F，De=Dw=D，对所有控制容积成立。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 对流项的其他离散格式对流项的其他离散格式设场变量经过对流扩散过程从一维区域的x=0点传输到x=L点，如图。流体密度为=1.0kg/m3，L=1.0m，扩散系数=0.1kg/(ms)。求流速分别为0.1m/s和2.5m/s时，将区域离散成10个节点网

15、格时，的分布。x=0 x=L=1=0u用中心差商离散，其离散方程为： 11;22 where 1122EeeWwwPPEEWWPeewwEWewaDF aDFaaaaDFDFaaFF针对节点1和10，采用控制容积积分得到（中间过程略）：0000202 for point 122eeePWeEFFDDDDF202 for point 1022WWWLPEEWLLLFFDDDDF129103.050.9502.101.0520.950 for 0.1m/s1.0520.950001.052.950u0.01.00.00.40.

16、0.91.0 Analytic u=0.1m/sx, m(x)=2.7183-exp(x)/1.7183上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 对流项的其他离散格式对流项的其他离散格式设场变量经过对流扩散过程从一维区域的x=0点传输到x=L点，如图。流体密度为=1.0kg/m3，L=1.0m，扩散系数=0.1kg/(ms)。求流速分别为0.1m/s和2.5m/s时，将区域离散成10个节点网格时，的分布。x=0 x=L=1=0u用中心差商离散，其离散方程为： 11;22 wh

17、ere 1122EeeWwwPPEEWWPeewwEWewaDF aDFaaaaDFDFaaFF针对节点1和10，采用控制容积积分得到（中间过程略）：0000202 for point 122eeePWeEFFDDDDF202 for point 1022WWWLPEEWLLLFFDDDDF1219207.250.7506.503.2540.7503.2540.750003.254.740123452.750.753.51.7510.7501.7510.7501.7510.7501.750.2500.01.00.00.40

18、.0.91.02.6 Analytic x=0.05m x=0.2mx, m=1+1-exp(25x)/(7.21010)上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 对流项的其他离散格式对流项的其他离散格式设场变量经过对流扩散过程从一维区域的x=0点传输到x=L点，如图。流体密度为=1.0kg/m3，L=1.0m，扩散系数=0.1kg/(ms)。求流速分别为0.1m/s和2.5m/s时，将区域离散成10个节点网格时，的分布。x=0 x=L=1

19、=0u用中心差商离散，其离散方程为： 11;22 where 1122EeeWwwPPEEWWPeewwEWewaDF aDFaaaaDFDFaaFF针对节点1和10，采用控制容积积分得到（中间过程略）：0000202 for point 122eeePWeEFFDDDDF202 for point 1022WWWLPEEWLLLFFDDDDF1234540.53.533.50.5033.50.5033.50.50340123453.50002.52.5002.52.5002.53.500.01.00.0

20、1.02.7 Analytic Difference Upwind Mixedx, m上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 对流项的其他离散格式对流项的其他离散格式1.中心差分格式的解，在小|Pe|时和严格解很符合，但当|Pe|=2的附近，误差明显增大。当|Pe|2后，结果均已超出边界值，违反了物理上的真实性。2.上风格式对任何Pe都能得到物理上真实解，但在整个Pe范围内都有显著误差。这是由于在小Pe时，对流项用了上游节

21、点值计算；而在大Pe时，扩散项仍用中心差分格式离散所致。3.混合格式所得结果在|Pe|2时和中心差分格式的相同，|Pe|2时比上风格式有很大改进，因为此时将扩散项取成了零，因此和严格解符合很好。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 目录目录上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 多维对流和扩散多维对流和扩散对于二维对流和扩散问题，通用方程的二维形式为：uuSxyxxyy 写成总流量形式 wh

22、ere xyxyJuJJxSxyJuy PNSEWnswe(x)w+(x)e-x(x)n-(x)s+yJwJeJsJn将上式在如图所示的控制容积上积分，得到enneennewsswwsswxyyxxyyxyxxyyxxyyxJJd dydxdxdyddydxdSdxdydxy 00enwsxyPPPPxyd dydxx y ,neswyxxx ex wyxJdxdydJJyx ,enwsxyyy ny sxyJdydxdJJxy neswyxCPPyxSdxdydSSx y 上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值

23、对流与换热对流与换热 多维对流和扩散多维对流和扩散对于二维对流和扩散问题，通用方程的二维形式为：uuSxyxxyy 写成总流量形式 where xyxyJuJJxSxyJuy PNSEWnswe(x)w+(x)e-x(x)n-(x)s+yJwJeJsJn将上式在如图所示的控制容积上积分，得到00PPPPewnsCPPx yJJJJSSx y 式中,源项已经线性化，S=SC+SPP；并假定控制容积内的和均可用P点的P和P表示；J是控制容积面上的总流量，如Je=Jdy。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热

24、对流与换热 多维对流和扩散多维对流和扩散针对连续方程：0uuxyPNSEWnswe(x)w+(x)e-x(x)n-(x)s+yJwJeJsJn也在控制容积上进行积分，得到：00PPewnsx yFFFF 式中，F是通过控制容积面上的流量，且假定e、w、n、s的流量分别代表其各自的F（对流强度），即：;ewsnewsnFuy Fuy Fux Fux将离散的连续方程与控制方程合并，得到： 00PPPeePwwPnnPssPCPPx yJFJFJFJFSSx y 采用混合格式代入，得到：00000max 0,2,;max 0,2,max 0,2,;max 0,2,+ where ;EeeeWwwwN

25、nnnSsssPPEEWWNNSSPPCPPPEWNSPPaDFFaDFFaDFFaDFFaaaaabax ybSx yaaaaaaaSx y 式中“0”表示时间的上一时间层的已知值。各对流项F的定义已知，相应的扩散传导系数D则由下式定义：;ewnsewnsewnsyyxxDDDDxxyy 上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 多维对流和扩散多维对流和扩散00000max 0,2,;max 0,2,max 0,2,;max 0,2,+ where ;EeeeWwwwNnnnSsssPPEEW

26、WNNSSPPCPPPEWNSPPaDFFaDFFaDFFaDFFaaaaabax ybSx yaaaaaaaSx y 上述方程是采用混合格式进行推导得到的，相应的系数表达是混合格式。如果写成更一般的格式，可以得到：00000max,0 ;max,0max,0 ;max,0;EeeWwwewNnnSssnsPPCPPPEWNSPPaD APeFaD APeFaD APeFaD APeFax ybSx yaaaaaaaSx y 式中的A(|Pe|)称为贝克来函数，它代表着对流扩散差分格式的形式。下表给出不同格式下的计算形式。格式A(|Pe|)表达式格式A(|Pe|)表达式中心差分1-0.5|Pe

27、|乘方格式max0, 1-0.1|Pe|5上风格式1指数格式|Pe|/exp(|Pe|)-1混合格式max0, 1-0.5|Pe|上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 多维对流和扩散多维对流和扩散二维问题的推导同样适用于三维情况，这里直接给出三维对流扩散问题的离散化方程：00000max,0max,0max,0max,0max,0max,0EeeeWwwwNnnnSsssTtttBbbbPPCPPPEWNSTBPPaD APeFaD APeFaD APeFaD APeFaD APeFaD A

28、PeFax y zbSx y zaaaaaaaaaSx y z +PPEEWWNNSSTTBBaaaaaaab式中 ,;,;,;,eeewwweewwnnnsssnnsstttbbbttbbFuy z Dy zxFuy z Dy zxFuz x Dz xyFuz x Dz xyFux y Dx yzFux y Dx yz 流率（对流强度）F及扩散传导系数D分别定义为贝克来数仍然定义为：Pe=F/D。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 目录目录上海大学上海大学冶金工程冶金工程专业本科生课程专

29、业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 虚假扩散虚假扩散在处理有流动存在的问题时，经常会遇到所谓的（或人工扩散）问题，比如因为上风格式中的系数aW和aE等系数比中心差分格式相应的系数大|F|/2，这相当于上风格式在真实的广义扩散系数中增加了一个大小为|F|/2的虚假的扩散系数。虚假扩散一般产生于下列三种情况：1.离散化时，非稳态项或对流项采用一阶截断误差的格式，为此可以采用二阶上风格式（又称为QUICK格式（quadratic upstream interpolation for convertive kinematics）等高阶离散化格式

30、加以避免；2.流动方向与网格线呈倾斜交叉状的多维对流扩散问题；3.建立差分格式时没有考虑非常数的源项的影响。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 虚假扩散虚假扩散QUICK格式是一种对流项的二阶上风格式，是英国Leonard于1979年提出的用于计算控制容积界面值的二阶插值格式。它利用控制容积界面两侧的紧挨着的邻近点和处于上风侧的一个远邻近节点，共三个节点的值来进行插值计算。如图所示。认为，界面w处的对流项中的w值除邻近节点P、W影响外，还受其上游的远邻近节点WW影响，而对于反向流动的情况，

31、则受节点E的影响。于是有PWWWEEEweii+1i+2i-1i-2uwue2631 for 028888PWPWWWwWWPWWwu2631 for 028888WPWPEwPPWEwu而扩散项中的可采用上述三点构造的拟合曲线在界面处的斜率计算，也可采用中心差商格式计算。上海大学上海大学冶金工程冶金工程专业本科生课程专业本科生课程 吴永全吴永全2022年7月4日5时27分冶金数值冶金数值 对流与换热对流与换热 虚假扩散虚假扩散以一维稳态对流问题为例。针对一维对流扩散问题控制容积积分方程PWWWEEEweii+1i+2i-1i-2uwue631631888888ePEWewWPWWwewuuFAFA ewewuuxx 其中，当uw0,ue0时，式中对流项采用QUICK格式计算，有扩散项采用二阶中心差分格式，整理得到最后的通式6131,8888+ where WwweEeeWWwPPEEWWWWWWPWEWWewaDFF aDF aFaaaaaaaaFF 上海大学上海大学