扩散方程中源项线性化问题的研究.doc
豆丁网精品论文扩散方程中源项线性化问题的研究张敏，胡文斌，李磊，许彬，卢学山南京理工大学动力工程学院，南京 (210094)摘要：用微分方法对扩散方程源项进行线性化处理。以此保证数值计算的收敛性和精确性。 在二个热传导算例的基础上，证明和展示本方法的可行性和实用性。关键词：方程源项，线性化，扩散方程扩散方程是自然界中描述物理过程三大类型方程之一，它最常用的形式是对导热扩散过 程的描述。对于扩散方程中源项的数值离散，在数值模拟计算中是十分重要的。在这方面美 国学者 Patankar1-2和新加坡学者 Chai3-4做了许多工作。本文对扩散基本方程源项的线性化进行系统的描述，并通过两个算例，给出在非结构化网格中，对源项的具体操作过程。1方程源项线性化的基本概念和关系式非稳态的扩散方程或导热方程，对一个标量物理变量 可写成， c p= + S(1.1)txi xi 其中， S 是单位体积中的净源项，是对应于变量 的扩散系数。对于方程（1.1）中的源项，可以分成两个部分，S = Sc + S p p(1.2)其中，Sc 和 S p 为常数项。在有限容积方法中，对于一个控制容积 P，式（1.2）为对标量 p的线性方程。当源项 S 为一个非线性函数时， Sc 和 S p 两量也将变成标量 p 的函数值，此时在数值计算中，我们不得不通过叠代来更新它们的值。在进行 Sc 和 S p 线性化中，有许多方法可以选择，这其中有一种最佳的方法形式如下式所示，ppS = S * + ( dS )* (d * )(1.3)在这里，假设 S 为 的一阶可微函数。星号*表示当前时刻的估计值，如初始值或猜测值。比较方程（1.2）和（1.3），可得到，S = S * ( dS )* *(1.4a)cd pS = ( dS )*pd(1.4b)对于 S p 在此有严格的要求，它必须小于零；否则 S p 设为零，即整个源项为 Sc 。下面通过 例子来解释这方面的处理技巧。1、当 S = 7 3 ，显然有 SC = 7, SP = 32、当 S = 6 + 5 ，如果我们选取 SC = 6, SP = 5 ，将得到 SP 为正值，这是我们不希望出现CPP的。因此另一种选择为， S= 6 + 5 * , S= 0 ，这其中， SC 必须不断叠代更新。3、当 S = 2 3 3 ，我们采用公式(1.4)有，S = S * ( dS )* * = 2 3 * 3 + 9 * 3 = 2 + 6 * 3(1.5a)cdSpPPP= ( dS )* = 9 * 2(1.5b)pd P这样我们确保了 SP 为负值。当然我们还有另外的选择。P4 、当 S = 2 + 3 3 ，如果我 们还 采用公式 (1.4) 将 会导致 S 为正值 ，所以 有CPPS = 2 + 3 *3 , S= 0 。对于方程源项为自变量函数时，它的线性化对收敛求解离散方程十分重要。所谓的线性 化就是采用公式(1.4)保证 SP 不为正值。当然，在保证 SP 不为正值时，源项的处理是有多种 可能的。下面结合实际，给出两个算例。2算例分析算例 1. 二维稳态扩散问题：假设有温度扩散方程（极坐标），1 kr T + sin(2 ) = 0(2.1a)r r r 因为, 方程源项为与自变量 T 无关的函数，我们可取，Sc = sin(2 ),S p = 0(2.1b)图 1 为温度场计算结果。同时，图 2-4 给出了按下列源项计算的结果。如果读者愿意从审美 的视角出发，我们不难找到科学的艺术性5-6。sin(3 ),cos(2 ),cos(3 ) .(2.1c)图 1 不同源项温度场分布算例 2. 二维非稳态扩散问题：实心圆柱体，初始温度为零度，在 0 < t 5s 的时间段内， 物体内产生热量，单位体积发热率为 q(t ) = t ；在 r = R 的边界处，温度保持为零度。非稳态导热方程和精确解7-8为， c T = 1 kr T + t(2.2a)p tr r r T (r, t) = 2 e 2 tJ0 ( m r) t q( )e 2 t d （ 0 < t 5s ）源项的处理有，R m =1mmJ1 ( m R) 0(2.2b)Sc = t p,S p = 0(2.2c)图 2 圆域的网格划分采用显式格式，时间步长 = 0.0005s ，计算结果如图 3 所示。时间步长再增大就有 可能发散。采用 C-N 格式和全隐格式计算了 t = 1s 和 t = 5s 时的温度分布，时间步长 = 0.0005s ，计算结果如图 4 和图 5 所示。表 1 给出了三种格式计算的平均误差。由图中的温度分布可以看出，数值解与精确解基本重合，在趋近圆心的计算结果吻合得不好，对于该处计算误差比较大。表 1 三种格式计算的平均误差时间 (s)显式格式C-N 格式全隐格式13.43%2.89%2.34%54.67%3.67%3.12%110.750.750.50.50.250.25YY00-0.25-0.25-0.5-0.5-0.75-0.75-1-1-0.500.51X-1-1-0.500.51X(a) 1 秒时的温度分布(b) 5 秒时的温度分布图 3 采用显式格式计算得到的温度分布图110.750.750.50.130.50.250.25YY00-0.25-0.25-0.5-0.5-0.75-0.75-1-1-0.500.51X-1-1 -0.5 0 0.5 1X(a) 1 秒时的温度分布(b) 5 秒时的温度分布图 4 采用 C-N 格式计算得到的温度分布图1 10.750.750.50.50.250.25YY0 0-0.25-0.25-0.5-0.5-0.75-0.75-1-1 -0.5 0 0.51X-1-1 -0.5 0 0.5 1X(a) 1 秒时的温度分布(b) 5 秒时的温度分布图 5 采用全隐格式计算得到的温度分3 结果及讨论对于抛物型和椭圆型扩散方程中的源项进行线性化的处理，在数值分析中十分重要，因 为它直接影响到计算的收敛性和精确性。本文通过系统地介绍源项线性化的处理过程，和给出具体算例，证明和展示了此方法的 可行性和实用性。参考文献1Patankar, S. V. Computation of Conduction and Duct Flow Heat Transfer. Innovative Research, Inc, 1991. 2Patankar，SV Numerical Heat Transfer and Fluid FlowNew York，Hemisphere Publishing，1981.3Chai, J. C., Moder, J. P., and Karki, K. C., 2001, "A Procedure for View Factor Calculation Using theFinite-Volume Method," Numerical Heat Transfer, Part B, Vol. 40, No. 1, pp. 23-35.4Chai, J. C., Zhang, M., Moder, J. P., and Patankar, S. V., 2001, "Conduction Heat Transfer Calculations Using Structured and Unstructured Grids," ICHMT Symposium CHT'01 - Advances in Numerical Heat Transfer, Vol. 1, pp. 519-526. (Palm Cove, Cairns, Queens land, Australia).5Zhang M Modeling of Radiative Heat Transfer and Diffusion Processes Using UnstructuredGridPhDDissertation，2000，Tennessee Technological University，USA6Kakac, S. and Yener, Y., Heat Conduction, Third Edition, Taylor & Francis, Publisher, 1993.7商立英非稳态导热问题在非结构化网格中的计算研究南京理工大学硕士学位论文，20068MN奥齐西克 著，俞昌铭 译热传导M北京：高等教育出版社，1984A Study of the Linear Source Terms for Diffusion EquationsZhang Min，Hu Wenbin，Li Lei，Xu Bin，Lu XueshanSchool of Power Engineering，Nanjing University of Science & Technology，Nanjing (210094)AbstractUnsteady/steady state diffusion equations solved using a cell-based finite volume method in unstructuredgrid. The linear source terms were adopted. The two heat conduction examples including the different source terms were given. All above not only demonstrate the correctness of the method and program, but also provide valuable reference for calculation of many source terms.Keywords：source terms，linear method，diffusion equations
