matlab解决一类问题,基于MATLAB的一类输运问题的数值分析

第 28 卷 哈尔滨师范大学自然科学学报 Vol． 28，No． 1 2012 第 1 期 NATURAL SCIENCES JOURNAL OF HARBIN NORMAL UNI  VERSITY 基于 MATLAB 的一类输运问题的数值分析* 孙丽男1，张馥菊2 ( 1． 黑河学院;2 哈尔滨师范大学) 【摘 要】利用差分方法对一类输运问题在一维和二维空间中进行数值分析，得到相应数值解，并将数值解与解析解进行比较，且进一步分析了数值方法的有效性． 【关键词】输运问题; MATLAB; 数值分析 收稿日期:2011 －11 －15 * 黑龙江省教育厅科学技术研究项目资助( 11551308) 1 问题提出 在实际生活和科学研究中，经常会遇到当某些分子的物理量在一个截面两边空间分布不均匀时，就会有物理量流从截面的一侧流向另一侧，此类现象，通常称其为输运现象． 例如，热量从温度高的地方向温度低的地方传递; 粒子从浓度高的地方向浓度低的地方传递，前者是热传导现象，后者是扩散现象，二者都是常见的输运现象． 事实上，输运现象在相当广泛的自然现象和日常生活中发生，如中子在核反应堆中的迁移及其所导致的动力学变化; 光子如何从太阳发射和如何穿过地球大气传播到地面的辐射输运; 街道中的车辆输运; 城市中的居民输运; 如此等等． 到目前为止，较为完整的输运理论已经被建立起来，并且已经成为工程技术和物理研究中的重要课题，如土木工程中对非饱和混凝土中氯离子输运问题［1］; 如实验核物理和反应堆物理中光子 －电子耦合的输运问题; 如天体物理、惯性约束聚变研究中的辐射输运问题［2］; 如此等等．该文将研究以下输运问题 ut － a2Δu = － b u， x ∈ Ω，t ＞ 0 u( x，t) = 0， x ∈ Ω，t ≥ 0 u( x，0) = ( x) ， x ∈ 珚{ Ω ( 1) 其中，第一个方程左端第二项为扩散项，右端为对流项，因此称其为对流扩散方程，它是一类基本的运动方程，许多关于物质运输、分子扩散等物理过程及粘体流体流动的数学模型都是由其定解问题刻画的． 这类方程还可以用来描述河流污染、大气污染、核废物污染中污染物质的分布，流体的流动和流体中热传导等众多物理现象． 因此，对其进行研究具有十分重要的理论和 实际应用意义． 对于问题( 1) ，可以采用傅里叶级数法、冲量定理法［3］ 等来求其解析解，但是这些方法都不是万能的，都存在这样或那样的局限性，从而使得要求一些偏微分方程的解析解变得十分困难，有些方程甚至无法求出解析解． 另一方面，即使求出了方程的解析解，由于其复杂的表达形式，解决实际问题时不一定适用． 为此，对一些定解问题进行数值分析十分必要的且非常实用的［4］． 以下，将以 MATLAB 为工具，采取不同的数值方法对问题( 1) 在一维和二维两种情况下进行数值分析并将数值结果与解析结果进行比较分析． 2 一维数值分析 该部分，考虑问题( 1) 在一维空间中的情况 哈尔滨师范大学自然科学学报 2012 年 第 28 卷 ut － a2uxx = － bux，x ∈ ( 0，l) ，t ＞ 0 u( 0，t) = 0，u( l，t) = 0，t ≥ 0 u( x，0) = ( x) ，x ∈［0，1{ ］ ( 2) 以下，首先给出问题( 2) 的几种有效的数值方法，然后利用这些方法进行数值模拟，最后将数值结果与解析结果进行比较． 2． 1 数值方法 对于问题(