对流扩散方程.pdfVIP

徐 州 工 程 学 院 课 程 设 计 报 告 课程名称 偏微分方程数值解 课题名称 对流扩散方程 的迎风格式的推导和求解 专 业 信息与计算科学 班 级 10 信计 3 姓名学号 指导教师 杨 扬 2013 年 5 月 23 日 一、实验目的： 进一步巩固理论学习的结果，学习双曲型对流扩散方程的迎风格式的构造 方法，以及稳定的条件。从而进一步了解差分求解偏微分方程的一些基本概念， 掌握数值求解偏微分方程的基本过程。在此基础上考虑如何使用 Matlab 的软件 进行上机实现，并针对具体的题目给出相应的数值计算结果。 二、实验题目 ： u t a 1u x b 1u xx f , 0 x 1, 0 t 1 u ( x , 0 ) exp( x / 2 ); u ( 0 , t ) exp( t ); u (1, t ) exp( 1 / 2 t ); 其中 a1=1,b1=2, f ( x , t ) exp( x / 2 t ) 。 用迎风格式求解双曲型对流扩散方程 , 观差分解对真解的敛散性 ( u exp(x / 2 t ) 三、实验原理： 1、用迎风格式求解双曲型对流扩散方程，迎风格式为： n 1 n n n n n n u u u u u 2 u u j j j j 1 j 1 j j 1 n a 1 b1 2 f j ( a 1 0) h h n 1 n n n n n n u u u u u 2 u u j j j 1 j j 1 j j 1 n a 1 b1 2 f j ( a 1 0) h h 2 若令 r a 1 * / h , b1 * / h , 则迎风格式可整理为： n 1 n n n n u j (1 r 2 ) u j ( r ) u j 1 u j 1 f j ( a1 0 ) 1 n 1 n n n n u j (1 r 2 ) u j ( r ) u j 1 u j 1 f j ( a1 0 ) 2 2、稳定条件： h 2 /( 2 * b 1 a 1 * h ), （ a 1 0 ）