对流扩散方程有限差分方法

对流扩散方程有限差分方法

对流扩散方程有限差分方法

求解对流扩散方程的差分格式有很多种，在本节中将介绍以下3种有限差分格式：中心差分格式、Samarskii格式、Crank-Nicolson型隐式差分格式。

3.1中心差分格式

时间导数用向前差商、空间导数用中心差商来逼近，那么就得到了流扩散方程的显示格式。

处进行Taylor展开: 1)

式的中心差分格式[6]

n 1 n U j U j

n n

U j 1 U j 1 a

2h

n

U j 1

v

n n

2U j U j 1

h2

(3)

若令a h，

n 1 U j

n

U j

Vp，则

h

1 / n

2(U 1

(3)式可改写为

n n

U j 1) (U j 1

2u：n \

U j 1)

(4)

从上式我们看到, 在新的时间层n 1上只包含了一个未知量n

U j

1，它可以由时间层n上的值U；1，U j n，U；1直接计算出来。因此, 中心差分格式是求解对

假定u(x,t)是定解问题的充分光滑的解，将n 1

U j n

U j

U； 1 分别在(X j,t n)

n

U

j

U(X j,t n 1) U(X j,t n) 0( 2)

n

U j 1

U(X j 1,t n) U(X j,t n)

n

U j 1 U(X j 1,t n) U(X j,t n) U n h2 2 U n X j 2 2 X j

U n h2

2

U n

X

j

2 2 X j

代入⑷式，有

T (X j,t n)

n 1

U

j

n

U

j

n n

U j 1 U j 1 a

2h

2U n

h2

n

0()

n

2

a 0(h )

2

U

2

X

n

2

v 0(h )

j

h

h

n

U j 1

0(h3)

0(h3)

n

U j 1

v ---

2

0( h )

显然，当

0, h 0时，T （X j ,t n ） 0，即中心差分格式与定解问题是

相容的。由以上的讨论也可得知，对流扩散方程的中心差分格式的截断误差为

2

O( h )。

对于我们上面构造的差分格式，是否可以直接用于实际计算呢？也就是 说，如果初始

值有误差，在计算过程中误差会不会扩大传播呢？这就是接下来 我们要讨论的是差分方程的稳定性问题。 下面用Fourier 方法来分析中心差分格

式的稳定性。

上式可以改写为

由此得到差分格式（3）的稳定性限制为

n n 2

n

u

u u

a —— v — t j

x

j x

.

2

0( ) (a v) O(h )

令u ；

v n e ikjh ，代入到（

4): 式

n 1 ikjh

v e

v n e ikjh

1

(v n e ik(j 2

1)h

n ik( j 1)h\

v e

)

(v n e

ik(j 1)h

整理得

n 1

一 _

v [1 2

(1 cos kh) i

sin kh]v n

所以该差分格式的增长因子为：

G( ,k) 1 2 (1 coskh) i sin kh

2v n e ikjh v n e ik(j

1)h

)

1 (1 由于 1 coskh 0，

(1 coskh) coskh)[4

所以G （ ,k) 1 (1 4 2(1 coskh)2 2 4 (1 coskh) （即差分格式稳定） 2

coskh) (1 coskh) 0

2

(sin kh)2

2

2

(1 coskh)]

的充分条件为

(2 2

2)—^ 4 2 2 0

1

注意到』

coskh ）

[0,1]，所以上面不等式满足的条件为

(2 2

8

2

) 4 2 2

0， 4 2 2

0。

其模的平方为

G( ,k)

[1 2 (1 coskh)]2

2

(sin kh)2

2v

1

2

， V

N

a

h 2

故有结论：对流扩散方程的中心差分格式是条件稳定的。 根据Lax 等价定理,

我们可以知道，对流扩散方程的中心差分格式是条件收敛的。

3.2 Samarskii 格式

设a>0，先对方程（1）作扰动，得到另一个对流扩散方程 ⑺

2

u v 2 X

的解

用Taylor 级数展开有

再令

用Taylor 级数展开有

2

U n

2

V( 2)j O(h ) X

u(X j ,t n

) a U(X j 」n ) h

(

1 1 八 U(X j 1,t n ) 2U(X j ,t n ) U(X j 1,t n ) 1)v

R h 2

(5)

其中R

对于 n

U j

1

— ha ，当 2v （5）式，构造迎风格式

1 n

U j

h 0时, (5) 式化为（1）式

n

小 n

n

1 U j 1 2U j U j 1

v — 1 R h 2

差分格式（6）称为逼近对流扩散方程的 Samarskii 格式。

首先推导（6）的截断误差。设U （x,t ）是对流扩散方程（1）式的充分光滑

n n U j U j 1 a (6)

T j n

U(X j ,t n 1) U(X j ,t n ) U(X j ,t n ) U(X j 1 ,

t n )

1 U(X j 1,t n ) 2u(X j ,t n ) U V 1 R

(X j 1九) h 2

n U

(X j ,t n1)U (X j ,t n )

j

U(X j i ，t n )

2u(X j ,t n ) U(X j i ，t n )

h 2

2

(U 、n

)

O( h 2)

ah 2

u

）n

2

u、n Rv( 2)j X

2 厂O(h2)

对流扩散方程

徐州工程学院 课程设计报告 课程名称偏微分方程数值解 课题名称对流扩散方程 的迎风格式的推导和求解专业信息与计算科学 班级10信计3 姓名学号 指导教师杨扬 2013年 5 月23 日

一、实验目的： 进一步巩固理论学习的结果，学习双曲型对流扩散方程的迎风格式的构造 方法，以及稳定的条件。从而进一步了解差分求解偏微分方程的一些基本概念，掌握数值求解偏微分方程的基本过程。在此基础上考虑如何使用Matlab 的软件进行上机实现，并针对具体的题目给出相应的数值计算结果。 二、实验题目： ?? ? ??-=-==<<<<+=+);2/1exp(),1();exp(),0();2/exp()0,(10,10,11t t u t t u x x u t x f u b u a u xx x t 其中a1=1,b1=2, ) 2/exp(),(t x t x f --=。 用迎风格式求解双曲型对流扩散方程,观差分解对真解的敛散性()2/exp(t x u -= 三、实验原理： 1、用迎风格式求解双曲型对流扩散方程，迎风格式为： ) 01(21 1 )01(2112 1 1112 1 11 1<++-=-+->++-=-+--+++-+-+a f h u u u b h u u a u u a f h u u u b h u u a u u n j n j n j n j n j n j n j n j n j n j n j n j n j n j n j n j τ τ 若令,/*1,/*12h b h a r τμτ== 则迎风格式可整理为： > <<++-+-+=><>++++--=-+++-+2)01()()21(1)01()()21(111111a f u u r u r u a f u u r u r u n j n j n j n j n j n j n j n j n j n j τμμμτμμμ2、稳定条件： ） （） （01),*11*2/(01),*11*2/(2 2<-≤>+≤a h a b h a h a b h ττ（*） 四、数值实验的过程、相关程序及结果： 本次的实验题目所给出的边界条件是第一边界条件，直接利用所给的边界条件，我们可以给出界点处以及第0层的函数值，根据a1的正负性，使用相应的<1>或者<2>式，求出其他层的函数值。误差转化成图的形式，并输出最大值。 针对三种不同的输入对应输出结果 ：

扩散方程稳态扩散与非稳态扩散

扩散方程稳态扩散与非稳态扩散 1．稳态扩散下的菲克第一定律（一定时间内，浓度不随时间变化dc/dt=0） 单位时间内通过垂直于扩散方向的单位截面积的扩散物质流量（扩散通量）与该面积处的浓度梯度成正比 即J=－D（dc/dx） 其中D：扩散系数，cm2/s，J：扩散通量，g/cm2〃s ，式中负号表明扩散通量的方向与浓度梯度方向相反。 可见，只要存在浓度梯度，就会引起原子的扩散。 x轴上两单位面积1和2，间距dx，面上原子浓度为C1、C2 则平面1到平面2上原子数n1=C1dx ，平面2到平面1上原子数n2=C2dx 若原子平均跳动频率f, dt时间内跳离平面1的原子数为 n1f〃dt 跳离平面2的原子数为n2fdt，但沿一个方向只有1/2的几率，则单位时间内两者的差值即扩散原子净流量。 令，则上式 2．扩散系数的测定：

其中一种方法可通过碳在γ-Fe中的扩散来测定纯Fe的空心园筒，心部通渗碳气氛，外部为脱碳气氛，在一定温度 下经过一定时间后，碳原子从内壁渗入，外壁渗出达到平衡，则为稳态扩散单位时单位面积中碳流量： A：圆筒总面积，r及L：园筒半径及长度，q：通过圆筒的碳量 则： 即： 则： q可通过炉内脱碳气体的增碳求得，再通过剥层法测出不同r处的碳含量，作出C-lnr曲线可求得D。 第一定律可用来处理扩散中浓度不因时间变化的问 3．菲克第二定律：解决溶质浓度随时间变化的情况，即dc/dt≠0

两个相距dx垂直x轴的平面组成的微体积，J1、J2为进入、流出两平面间的扩散通量，扩散中浓度变化为，则单元体积中溶质积累速率为 （Fick第一定律） （Fick第一定律） （即第二个面的扩散通量为第一个面注入的溶质与在这一段距离内溶质浓度变化引起的扩散通量之和） 若D不随浓度变化，则 故： 4．Fick第二定律的解：很复杂，只给出两个较简单但常见问题的解 a. 无限大物体中的扩散

第三章 一维扩散方程

第三章 一维扩散方程 本章讨论一维扩散方程。首先，从随机过程中的一维扩散方程的讨论可直接得到扩散方程的解。然后对非齐次和各类边值问题相应的扩散方程作了讨论。讨论的方程类型 （1）直线上的齐次和非齐次扩散方程： 2,,0 (,0)() t xx u c u x t u x x ??=-∞<<∞>? =?；（利用随机过程的理论得到结论，再直接验证） (,),,0 (,0)() t xx u ku f x t x t u x x ?-=-∞<<∞>?? =?；（算子方法，与常微分方程类比） （2）半直线上的扩散方程0,0,0(,0)(),(0,)0t xx u ku x t u x x u t ?-=<<∞>?? =??=? ；（其它非齐次边界等） 对扩散方程理论方面的探讨：最大（最小）值原理。由此证明方程解的唯一性和稳定性。 §3.1全直线上的扩散方程 首先讨论随机过程中的扩散过程。设想粒子在一维直线上作连续随机游动（Brown 运动），满足性质：在t ?时间内位移转移概率为均值为0，方差为2 t σ?的正态分布。在时刻t 处于x 的概率密度记为(,)Pr(())u x t dx X t x dx ==。则 2 ()2(,)(,)x y t u x t t u y t dy σ-∞ -?-∞+?=?， 或 2 2 (,)(,)y u x t t u x y t dy ∞ -+?= +? 2222 1 [(,)(,)(,)()]2 y x xx u x t u x t y u x t ty o t dy σ∞ - = ++?+?? 21 (,)(,)()2 xx u x t u x t t o t σ=+?+? 因此， 2 2 t xx u u σ= 。 可见：一维Brown 运动的状态概率密度满足扩散方程。 从随机过程的角度，可直接写出状态概率密度： 22()2(,)(,0)y x t u x t e u y dy σ-∞ - = ?。 所以，有如下定理。 定理 扩散方程2,,0 (,0)() t xx u c u x t u x x ??=-∞<<∞>?=?的解为

扩散方程的差分解法

扩散方程的差分解法 在研究热传导过程、扩散过程、边界层现象时，我们常常遇到抛物型方程，这类方程中最典型、最简单的就是热传导方程。热传导方程中的自变量中包括时间t ，它是描述一种随时间变化的物理过程，即所谓不定常现象。这类问题的基本定解问题应是初值问题，即在初始时刻（t=0）时给定定解条件，求解t>0时的解。 本文主要运用有限差分法对一维扩散方程进行求解，并对差分解的适定性、相容性、收敛性及稳定性进行分析，同时与解析解进行对比。 1.扩散方程 一维扩散方程为： 22u u t x α??=?? （1） 式中，u 为因知量，α为扩散系数，x 为坐标，t 为时间。 其定解条件如下： 初始条件： (,0)() 0x u x f x L =≤≤ （2） 边界条件： 12(0,)() , (,)()u t f t u L t f t == （3） 一般假定函数()f x ，1()f t ，2()f t 满足连接条件，即1(0)(0) f f =，2()(0) f L f =。 2.有限差分法 有限差分法是数值计算解微分方程古老的方法之一，也是系统化地、数值地求解数学物理方法的方程。其控制方程中的导数用离散点上函数值的差商代替。 差分格式可以分为显格式和隐格式。所谓显格式是指在任一结点上因变量在新是时间层上的值可以通过之前的时间层上相邻结点变量的值显式解出来。由于这些层的变量值是已知的，当时间向前推进时，空间点上的新的变量值就只需逐点计算就行了，因此显格式计算起来比较省事。隐格式则是指任一结点上变量在新的时间层的值，不能通过之前的时间层上相邻结点的值显式解出来，它不仅与之前的时间层上的已知值有关，而且也与新时间层的相邻结点的变量值有关。因而一个差分方程常常包括几个相邻结点上的未知数，未知数的个数取决于格式的构成形式。为了解出这些未知数需要联立新的方程，而每引进一个新的方程往往又同时引进了新的未知数。因此，隐格式总是伴随着求解巨大的代数方程组。隐格式的主要缺点是计算工作量大，因而不如显格式计算得快，但这只是就时间步长一样的情况而言的。隐格式的主要优点是时间步长可以比显格式能够采用的最大步长大很多。显格式的时间步长受到稳定性条件的限制，而隐格式则几乎不受限制。 3.方程的离散 3.1 显格式 采用时间前差及第n 时间层的空间中心差，得一维扩散方程的显格式解： 111 2 2()n n n n n j j j j j u u u u u t x α ++---+=?? （4） 即 111(2) n n n n n j j j j j u u r u u u ++-=+-+ （5）

对流扩散方程引言

对流扩散方程的定解问题是指物质输运与分子扩散的物理过程和黏性流体流动的数学模型，它可以用来描述河流污染、大气污染、核污染中污染物质的分布，流体的流动和流体中热传导等众多物理现象。关于对流扩散方程的求解很也备受关注，因此寻找一种稳定实用的数值方法有着重要的理论与实际意义。 求解对流扩散方程的数值方法有多种，尤其是对流占优扩散方程，这些方法有迎风有限元法，有限体积法，特征有限体积法，特征有限差分法和特征有限元法，广义差分法，流线扩散法，以及这些方法与传统方法相结合的方法如迎风广义差分法，迎风有限体积法有限体积——有限元法等这些方法数值求解效果较好，及有效的避免了数值震荡，有减少了数值扩散，但是一般计算量偏大 近年，许多研究者进行了更加深入的研究，文献提出了对流扩散方程的特征混合元法，再次基础上，陈掌引入了特征间断混合元方法，还有一些学者将特征线和有限体积法相结合，提出了特征有限体积元方法（非线性和半线性），于此同时迎风有限元也得到较大的发展，胡建伟等人研究了对流扩散问题的Galerkin 部分迎风有限元方法和非线性对流扩散问题的迎风有限元，之后又有人对求解发展型对流扩散问题的迎风有限元法进行了理论分析 有限差分法和有限元是求解偏微分方程的常用数值方法，一般情况下考虑对流占优的扩散方程，当对流项其主导作用时，其解函数具有大梯度的过渡层和边界层，导致数值计算困难，采用一般的有限元或有限体积方法虽然具有形式上的高精度，不能解决数值震荡的问题，虽然我们不能简单的将对流占优扩散方程看做对流方程，但由于次方程中含有一阶不对称的导数，对流扩散方程仍会表现出“对流效应”，从而采用迎风格式逼近，尽量反应次迎风特点，此格式简单，克服了锋线前沿的数值震荡，计算结果稳定，之前的迎风格式只能达到一阶精度，我们采用高精度的广义迎风格式，此格式是守恒的，精度高，稳定性好，具有单调性，并且是特征线法的近似，有效的避免了锋线前沿的数值震荡。 有限体积是求解偏微分方程的新的离散技术，日益受到重视。有限体积与有限差分、有限元法最大的区别及优点在于有限体积将求解区域内的计算转化到控制体积边界上进行计算，而后二者均是直接（或间接）在域内计算，故有限体积有着明显的物理涵义，在很大程度上减少计算工作量又能满足计算精度要求，加快收敛速度。由于此方法讲散度的积分化为子域边界积分后子啊离散，数值解满足离散守恒，而且可以采用非结构网格，所以在计算物理特别是计算流体力学领域上有限体积有广阔的前景。 间断Galerkin（DG）方法是在1973年，Reed和Hill在求解种子迁移问题时，针对一阶双曲问题的物理特点提出的。之后C.Johnson，G.R.Richter等人对双曲问题的DG方法做了进一步的研究，并且得到了该机的误差分析结果，由于这种方法具有沿流线从“上游”到“下游”逐层逐单元计算的显示求解的特点，并且可以进行并行计算，所以被广泛应用于各类方程的求解。最近Douglas等人在{25}中处理二阶椭圆问题时，得到DG方法的有限元空间不需要满足任何连续性条件，因此空间构造简单，具有较好的局部性和并行性。DG发展的一个重要方面是对对流占优扩散方程的应用。G.R.Richter等在1992年提出利用DG方法求解定长对流扩散问题 近年DG方法有了新的发展，其中YeXiu提出间断体积元方法备受人们关注，2004年，她将有限体积法与DG相结合，提出了椭圆问题的间断有限体积法，此方法解除了逼近函数在跨越边界上连续的限制，之后更多的研究者应用到Stokes问题，抛物问题，双曲问题，并得到了较好的结果，该方法不但继承了有限体积的高精度计算简单及保持物理间局部守恒等优点，而且有限元空间无需满足任何连续性要求，空间构造简单，有较好的局部和并行性。 当对流扩散方程中的对流项占主导地位时，方程具有双曲方程的特点，这是由于对流扩散方程中的非对称的对流项所引起的迎风效应使对流扩散方程的数值求解更困难，用传统的中心差分法和标准的有限元求解会差生数值的震荡，从而使数值模拟失真，为了克服这一困难，早在20世纪50年代，就有人提出了迎风思想，由于使用迎风技巧可以有效的消除数值解不稳定性，因此吸引了众多学者的关注，从1977年，Tabata等人就针对对流扩散方程提出了三角形网格上的迎风格式{42,38}，并进行了深入的研究，梁栋基于广义差分法，提出并分析了一类建立在三角网格上的广义迎风差分格式，袁益让2001年就多层渗流方程组合系统提出并分析了迎风分数步长差分方法，以上均是讨论的线性对流扩散问题，胡建伟等通过引入质量集中算子，构造并分析了一类基于三角网格的质量集中型的部分有限元方法处理线性和非线性对流扩散问

一维扩散方程的差分法matlab

一维扩散方程的差分法 m a t l a b Company number【1089WT-1898YT-1W8CB-9UUT-92108】

一维扩散方程的有限差分法 ——计算物理实验作业七 陈万 题目： 编程求解一维扩散方程的解 取1.0,1.0,1.0,10,0.1,0,1,1,0,1,1max 0222111======-=====τh D t a c b a c b a 。输出t=1,2,...,10时刻的x 和u(x),并与解析解u=exp(x+0.1t)作比较。 主程序： % 一维扩散方程的有限差分法 clear,clc; %定义初始常量 a1 = 1; b1 = 1; c1 = 0; a2 = 1;b2 = -1; c2 = 0; a0 = 1.0; t_max = 10; D = 0.1; h = 0.1; tao = 0.1; %调用扩散方程子函数求解 u = diffuse_equation(a0,t_max,h,tao,D,a1,b1,c1,a2,b2,c2); 子程序1： function output = diffuse_equation(a0,t_max,h,tao,D,a1,b1,c1,a2,b2,c2) % 一维扩散方程的有限差分法，采用隐式六点差分格式(Crank-Nicolson) % a0: x 的最大值 % t:_max: t 的最大值

% h: 空间步长 % tao: 时间步长 % D：扩散系数 % a1,b1,c1是（x=0）边界条件的系数；a2,b2,c2是（x=a0）边界条件的系数 x = 0:h:a0; n = length(x); t = 0:tao:t_max; k = length(t); P = tao * D/h^2; P1 = 1/P + 1; P2 = 1/P - 1; u = zeros(k,n); %初始条件 u(1,:) = exp(x); %求A矩阵的对角元素d d = zeros(1,n); d(1,1) = b1*P1+h*a1; d(2:(n-1),1) = 2*P1; d(n,1) = b2*P1+h*a2; %求A矩阵的对角元素下面一行元素e e = -ones(1,n-1);

fick定律扩散方程

扩散方程 扩散方程稳态扩散与非稳态扩散 1．稳态扩散下的菲克第一定律（一定时间内，浓度不随时间变化dc/dt=0） 单位时间内通过垂直于扩散方向的单位截面积的扩散物质流量（扩散通量）与该面积处的浓度梯度成正比 即J=－D（dc/dx） 其中D：扩散系数，cm2/s，J：扩散通量，g/cm2·s ，式中负号表明扩散通量的方向与浓度梯度方向相反。 可见，只要存在浓度梯度，就会引起原子的扩散。 x轴上两单位面积1和2，间距dx，面上原子浓度为C1、C2 则平面1到平面2上原子数n1=C1dx ，平面2到平面1上原子数n2=C2dx 若原子平均跳动频率f, dt时间内跳离平面1的原子数为n1f·dt 跳离平面2的原子数为n2fdt，但沿一个方向只有1/2的几率，则单位时间内两者的差值即扩散原子净流量。 令，则上式 2．扩散系数的测定： 其中一种方法可通过碳在γ-Fe中的扩散来测定纯Fe的空心园筒，心部通渗碳气氛，外部为脱碳气氛，在一定温度

下经过一定时间后，碳原子从内壁渗入，外壁渗出达到平衡，则为稳态扩散单位时单位面积中碳流量： A：圆筒总面积，r及L：园筒半径及长度，q：通过圆筒的碳量 则： 即： 则： q可通过炉内脱碳气体的增碳求得，再通过剥层法测出不同r处的碳含量，作出C-lnr曲线可求得D。 第一定律可用来处理扩散中浓度不因时间变化的问 3．菲克第二定律：解决溶质浓度随时间变化的情况，即dc/dt≠0 两个相距dx垂直x轴的平面组成的微体积，J1、J2为进入、流出两平面间的扩散通量，扩散 中浓度变化为，则单元体积中溶质积累速率为 （Fick第一定律） （Fick第一定律） ，，， （即第二个面的扩散通量为第一个面注入的溶质与在这一段距离内溶质浓度变化引起的扩散通

对流扩散方程有限差分方法.

对流扩散方程有限差分方法 求解对流扩散方程的差分格式有很多种，在本节中将介绍以下3种有限差分格式：中心差分格式、Samarskii 格式、Crank-Nicolson 型隐式差分格式。 3.1 中心差分格式 时间导数用向前差商、空间导数用中心差商来逼近，那么就得到了（1）式的中心差分格式]6[ 2 1 11 1122h u u u v h u u a u u n j n j n j n j n j n j n j -+-+++-=-+-τ （3） 若令 h a τ λ=，2h v τ μ=，则（3）式可改写为 )2()(2 111111 n j n j n j n j n j n j n j u u u u u u u -+-+++-+--=μλ （4） 从上式我们看到，在新的时间层1+n 上只包含了一个未知量1 +n j u ，它可以由时间层n 上的值n j u 1-，n j u ，n j u 1+直接计算出来。因此，中心差分格式是求解对 流扩散方程的显示格式。 假定),(t x u 是定解问题的充分光滑的解，将1 +n j u ，n j u 1+，n j u 1-分别在),(n j t x 处 进行Taylor 展开： )(),(),(211ττO t u t x u t x u u n j n j n j n j +??? ?????+==++ )(2),(),(3 22211 h O x u h x u h t x u t x u u n j n j n j n j n j +????????+????????+==++ )(2),(),(3 22211 h O x u h x u h t x u t x u u n j n j n j n j n j +????????+????????-==-- 代入(4)式，有 2 111 1122),(h u u u v h u u a u u t x T n j n j n j n j n j n j n j n j -+-+++---+-= τ )()()(2222 h O v x u v h O a x u a O t u n j n j n j ?-????????-?+????????++????????=τ )()()(222h O v a O x u v x u a t u n j n j n j ?-++????????-??? ?????+????????=τ

基于Peclet数判别法的一维对流扩散方程分类研究

基于Peclet 数判别法的一维对流扩散方程分类研究 摘要：采用Peclet 数的绝对值大小来判别一维对流扩散方程为对流占优型或是扩散占优型方程，运用三种隐式差分格式—中心隐式格式、对流C-N 型格式和扩散C-N 格式，对不同Peclet 数的算例进行离散和求解。然后，将计算区域中所有节点的解析解与数值解表示成矩阵形式，并求解出它们的矩阵2范数之后作比较，两者越接近则代表差分格式精度越高。通过比较得出了当方程Peclet 数的绝对值小于等于0.5时，方程为扩散占优型方程。在离散方法选取方面，针对扩散项的离散可以采用更高精度的差分格式，如扩散C-N 格式；当Peclet 数的绝对值大于等于20时，方程为对流占优型方程。此时，针对对流项可以采用更高精度的差分格式，如对流C-N 格式；当Peclet 数的绝对值介于0.5与20之间时，无法用Peclet 数判断方程类型,不过可以选择折衷的离散格式减小误差，如中心隐式格式。 关键字：一维对流扩散方程 Peclet 数判别法 有限差分方法 数值模拟 MR(2010)主题分类号：39A14；65M06 中图分类号：O242.2 文献标识码： A 1.引言 一维对流扩散方程是描述流体流动和传热问题的一类线性化模型方程。土壤、大气等多孔介质中水、盐分、温度以及污染物质的对流扩散问题都会遇到此类方程。在一维对流扩散方程的求解过程中，反映流体对流和扩散两种物理作用的分别是对流项和扩散项。所以，根据方程中对流项还是扩散项占主导作用，通常可将方程分为对流占优型和扩散占优型两类方程。然而，要想得到精确度较高的数值结果，这两种类型方程的离散方法不能采用相同的离散格式。因此，需要有一种判别方法来判断方程的类型，关于对流占优型和扩散占优型方程的判别方法一直是近年来研究的热点问题。这对研究不同类型的方程使用合适的差分格式进行离散具有实际的意义。由于Peclet 数的绝对值表示了对流作用相对扩散作用的大小，即绝 大，扩散所起的作用就可以忽略。反之，当Peclet 数为零时，方程就为纯扩散方程。本文选用一维定解非稳态对流扩散方程为例，通过考察Peclet 数的绝对值大小来对方程进行分类，方程一般形式如下： 2(,),,0 122(1)(,0)()(,)(),(,)()12(,) u u u a f x t x x x t t x x u x g x u x t t u x t t u u x t υ?φ???? ?? ?? ????+=+≤≤≥???==== 其中a 和υ分别代表对流项系数和扩散项系数。假定求解区间长度为s ， Peclet 数的绝对值

剖面二维非恒定悬移质泥沙扩散方程的数值方法.

剖面二维非恒定悬移质泥沙扩散方程的数值方法< 1 引言 数学模拟方法正在成为研究河流泥沙问题的重要手段。目前，一维数学模型发展较成熟，已广泛应用于模拟长河段的长期变形，但它只能给出河段平均冲淤深度的沿程变化，如需了解短河段的河床变形细节，则要采用二维以至三维数学模型。不论是一维数学模型还是平面二维维数学模型，都不能反映含沙量沿垂线的分布状况，并忽略了含沙量沿垂线分布对垂线平均含沙量变化过程的影响。要解决这类问题，必须建立剖面二维数学模型。这种模型主要通过解剖面二维泥沙扩散方程来研究悬移质泥沙沿水深的分布及含沙量的变化过程，对水电站进口和其它引水工程的引水口高程的确定都能提供较好的数值模拟。 泥沙扩散方程实际上是一个变系数的二阶线性偏微分方程，这样的方程在各种复杂边界条件下求解是极为困难的。求扩散方程的解析解在数学上存在着难以克服的困难，往往只能通过对方程的简化，才能得到一些简单边界条件下的解析解，在这方面， A.A.Kalinske 、野满隆治、 W.E.Dobbins 、俞维强、张启舜、韦直林 [ ２ ] 等都做了有益的尝试；求扩散方程的数值解曾经因为缺乏高效率的计算工具而难以实现，直到 60 年代后，随着计算机的广泛应用，在各种复杂边界条件下求扩散方程的数值解不但成为可能，而且得到迅速的发展，在这方面，曹志先、崔侠 [ ４ ] 等做了大量工作，取得了很多成果。 数值方法相对于解析方法在求解偏微分方程上有着明显的优势，即简单灵活、计算方便快捷，但要寻找一种精度高、稳定性好、计算方便的差分格式也并非易事。本文拟在前人研究的基础上着重讨论剖面二维泥沙扩散方程的数值解问题，希望能提供一种精度高、稳定性好、计算方便的数值解。 2 基本方程 剖面二维泥沙扩散方程的形式为

一维对流扩散方程的稳定性条件推导

一维稳态对流扩散方程稳定性条件的推导 姓名： 班级：硕5015 学号： 2015/12/15

证明： 一维稳态对流扩散方程： 22u x x φφρ??=Γ?? 采用控制容积积分法，对上图P 控制的容积作积分，取分段线性型线，对均分网格可得下列离散方程： ()()()()()()()()11112222e w e w P E W e w e w w w e e u u u u x x x x φρρφρφρδδδδ??????ΓΓΓΓ+-+=-++????????????????记：()()()()1122e w P e w w e a u u x x ρρδδΓΓ=+-+ ()()12 e E e e a u x ρδΓ=- ()()12w W w w a u x ρδΓ= + 定义通过界面的流量u ρ记为F ,界面上单位面积扩散阻力的倒数x δΓ记为D ,则原式简化为： P P E E W W a a a φφφ=+ 12 E e e a D F =- 12 W w w a D F =+ ()P E W e w a a a F F =++- 令 u x F Pe D ρδ==Γ 则 1111222 E W P Pe Pe φφφ????-++ ? ?????=

当Pe 大于2以后，数值解出现了异常；P φ小于其左右邻点之值，在无源项情 况下是不可能的。因为当2Pe >时系数12 E e e a D F =-小于零，即右边点的通过对流及扩散作用对中间点所产生的影响是负的，这会导致物理上产生不真实的解，所以2u x Pe ρδ=≤Γ 证毕。

二维扩散方程

一、用有限差分法求解二维扩散方程的初边值问题 该问题的精确解为1 ()2 (,,)x y t u x y t e +-= 二、用下列差分格式编程计算，并比较计算速度、精度、稳定性。 1． 古典显式格式： 1,,1,,1,,1,,1 2 2 22n n n n n n n n j l j l j l j l j l j l j l j l u u u u u u u u h h τ ++-+---+-+= + 2． 交替方向隐式格式（P-R 格式）： 1111 2222 ,,1,,1,,1,,1 22 1111 1111 2222,,1,,1,,1,,122 222 222 n n n n n n n n j l j l j l j l j l j l j l j l n n n n n n n n j l j l j l j l j l j l j l j l u u u u u u u u h h u u u u u u u u h h ττ+++++-+-+++++++++-+-?--+-+?=+?????--+-+?=+??? 3． 局部一维格式： 1111 2222 ,,1,,1,1,,1,221111 11112222 ,,,1,,1,1,,122 2211()222 2211()222 n n n n n n n n j l j l j l j l j l j l j l j l n n n n n n n n j l j l j l j l j l j l j l j l u u u u u u u u h h u u u u u u u u h h ττ+++++-+-++++ +++++-+-?--+-+?=+?????--+-+?=+??? 4． 预测-校正格式： 第一步：1111 4444 ,,1,,1,2 11111 24222 ,,,1,,1222 22 n n n n n j l j l j l j l j l n n n n n j l j l j l j l j l u u u u u h u u u u u h ττ+++++-++++++-?--+?=?????--+?=??? ， 122 ()2221 ()2 11 (1)2211 (1)223,(0,1,01) 2(,,0)(0,1)(0,,),(1,,)(01,01) (,0,),(,1,)(01,01)x y t x y y t y t x t x t u u u e x y t t x y u x y e x y u y t e u y t e y t u x t e u x t e x t +-+-+--+-?????-+=-<<<≤??????=<

一维扩散方程的差分法matlab

一维扩散方程的有限差分法 ——计算物理实验作业七 陈万 题目： 编程求解一维扩散方程的解 取1.0,1.0,1.0,10,0.1,0,1,1,0,1,1max 0222111======-=====τh D t a c b a c b a 。输出t=1,2,...,10时刻的x 和u(x),并与解析解u=exp(x+0.1t)作比较。 主程序： % 一维扩散方程的有限差分法 clear,clc; %定义初始常量 a1 = 1; b1 = 1; c1 = 0; a2 = 1;b2 = -1; c2 = 0; a0 = 1.0; t_max = 10; D = 0.1; h = 0.1; tao = 0.1; %调用扩散方程子函数求解 u = diffuse_equation(a0,t_max,h,tao,D,a1,b1,c1,a2,b2,c2); 子程序1：

function output = diffuse_equation(a0,t_max,h,tao,D,a1,b1,c1,a2,b2,c2) % 一维扩散方程的有限差分法，采用隐式六点差分格式 (Crank-Nicolson) % a0: x的最大值 % t:_max: t的最大值 % h: 空间步长 % tao: 时间步长 % D：扩散系数 % a1,b1,c1是（x=0）边界条件的系数；a2,b2,c2是（x=a0）边界条件的系数 x = 0:h:a0; n = length(x); t = 0:tao:t_max; k = length(t); P = tao * D/h^2; P1 = 1/P + 1; P2 = 1/P - 1;

一维对流扩散方程的数值解法

一维对流扩散方程的数值解法 对流-扩散方程是守恒定律控制方程的一种模型方程，它既是能量方程的表示形式，同时也可以认为是把压力梯度项隐含到了源项中去的动量方程的代表。因此，以对流-扩散方程为例，来研究数值求解偏微分方程的相容性、收敛性和稳定性具有代表性的意义。 1 数学模型 本作业从最简单的模型方程，即一维、稳态、无源项的对流扩散方程出发，方程如下： 22, 02f f f U D x t x x ???+=≤≤??? (1) 初始条件 (),0sin(2)f x t A kx π== (2) 解析解 ()()()224,sin 2Dk t f x t e A k x Ut ππ-=- (3) 式中，1,0.05,0.5,1U D A k ==== 函数(3)描述的是一个衰减波的图像，如图1所示 t=0 t=0.5 t=1 图1 函数()()()224,sin 2Dk t f x t e k x Ut ππ-=- 的图像（U=1，D=0.05，k=1） 2 数值解法 2.1 数值误差分析 在网格点(),i n 上差分方程的数值解n i f 偏离该点上相应的偏微分方程的精确解 (),f i n 的值，称为网格节点上的数值误差。 当取定网格节点数21N =时，观察差分方程的解与微分方程的解在不同时间步长下的趋近程度，其中时间步长分别取值0.05,0.025,0.0125,0.0005t ?=。

（a ）21,0.05N t =?= （b ）21,0.025N t =?= （c ）21,0.0125N t =?= （d ）201,0.0005N t =?= 图2 数值误差随步长的变化情况 从图2的(a)~(d)可以定性的看出，数值误差与步长的大小有关。在满足稳定性条件的前提下，数值误差随着时间步长的减小而减小，同时，图（d ）表示增大网格的分辨率也有助于减小网格误差。 为了对数值误差有一个定量的认识，接下来取定时间步长为0.0005t ?=，分别算出 11,21,41,61,81,101,121,161N =时，指标E =1所示。 表1 不同网格节点数下指标E 的值

物理分析课程设计---一维扩散方程求解

课程设计报告 课程名称：核反应堆物理分析题目：一维扩散方程求解院系：核科学与工程学院班级： 学号： 姓名： 指导教师： 成绩： 教师签名： 日期：2011 年6月日

目录 摘要 (1) 课程设计的目的与要求 (1) 设计正文 (1) 课程设计总结或结论 (3) 参考文献 (4)

摘要和关键词 摘要 这个设计用微分方程的差分数值求解方法，运用MATLAB编程计算出一维扩散方程中子通量密度的离散解。 关键词：一维扩散方程 一．课程设计的目的与要求 学习使用微分方程的数值解法（差分方法）来近似求解一维扩散方程，掌握差分方法的核心思想，熟练使用matlab数据处理，origin绘图软件。通过给定的微分方程及边界条件，计算平板型，圆柱形，球形反应堆中子通量密度分布。 二．设计正文 通过查找有关资料，根据二阶线性微分方程 ○1 转换为差分方程的一般公式 其中 ○2 h为给定步长, 我们把原方程化简为○3

对比方程○1和○3得出 ○4 把○4代入○2 等式右端向量 差分方程其实就是一个线性方程组，此线性方程组的系数矩阵为： 则有 这是一个三对角阵，故可用追赶法解式○3。 下面通过matlab程序来计算变换后的差分方程的解。 所编程序如下： clear; N=input('请输入参数：'); alpha=input('请输入alpha值：'); if alpha==0 rmax=input('请输入平板的厚度：'); f0=input('请输入平板中心的中子通量密度：'); elseif alpha==1 rmax=input('请输入堆芯半径：'); f0=input('请输入圆柱中心的中子通量密度：'); elseif alpha==2 rmax=input('请输入堆芯半径：'); f0=input('请输入球形中心的中子通量密度：'); end

对流扩散方程背景

对流扩散方程背景 提出一种隐格式用于解决二维时间依赖的Burgers型方程。迎风单边差分格式被用于对流项离散；对扩散项用二阶中心差分格式离散。我们建立了全隐的数值有限差分格式，分析了无条件稳定性和严格推导了收敛性，在空间是二阶收敛的和时间一阶收敛的。给出数值结果验证理论正确性。 关键词：隐格式，单边差分逼近，Burgers方程，稳定性，收敛阶 对流扩散方程背景 对流扩散方程描述黏性流体的动力学行为，这在许多工程应用中发挥了重要作用。对流占优型扩散方程一般具有对流比扩散的系数大得多的特点，通常数值模拟具有一定难度，因为一方面，扩散系数比传输速度小，并且在另一方面，由于数值扰动容易出现边界层现象。许多格式已用于这些问题的模拟，并有大量成功的数值方法[1-3]。通过离散方法来解决对流扩散问题时，一般运用标准Galerkin有限元方法求解，但此方法会导致非物理特征扰动。为了解决这类缺陷，几种稳定的有限元方法已经在[4]中被提出了。 我们感兴趣的是建立非耗散方法来克服数值扰动，并有鲁棒性和二阶精度，尤其是对Burgers问题。Burgers问题通常被认为非线性流体的流动和扰动的经典模型。在二维非线性的情况下，可以描述对流和扩散的现象，Burgers方程代表一种最基本的非线性模型方程。从一个数值格式出发研究是相当有趣的，因为Burgers已出现在众多的流体方程中[5-7]。并已经由霍普夫-科尔计算出多种组合的初始条件和边界条件下的结果[10,11]。此外，对于非线性Burgers方程的解析解也可以通过Homotop Perturbation法[12]得到。 众所周知，单独的选择一种基本差分格式如中心差分或者迎风格式，来计算纯对流式的方程，扩散项通常只是中心近似。而解决问题的关键在于对流方面构造稳定的离散结构来克服数值扰动。虽然单边差分近似格式已经被提出了30年之久[13]，人们却很少关注他们在计算流动问题。一阶或者二阶单边迎风有限差分

对流方程及算法介绍

1 引言 2 对流方程及算法介绍 2.1对流方程的概述 对流：是指由于流体的宏观运动，从而使流体各部分之间发生相对位移，冷热流体相互掺混所引起的热量传递过程。 对流仅发生在流体中，对流的同时必伴随有导热现象。 人们研究对流扩散方程，主要的研究对象是流体在流动过程中，流体所携带的某种物质的物理量的变化规律，例如传热过程中温度的变化规律或者溶解于流体中溶质的物质浓度等物理量的变化规律。这些变化通常包括对流、扩散以及由于某种物理或者化学的因素而引起的物理量的自身衰减或增长。 最简单的一维对流扩散方程形如(2-1)式: （2-1） 其中C 是常数，它属于双曲型方程，可以被用来描述流体的运动等物理现象。 2.2水对流现象的简易演示 2.2.1 基本步骤 用两只相同的小烧杯，各装上冷水，再如图1所示插入长短两根吸管，虹吸管由普通化学实验用玻璃管在酒精灯上加热弯成，一根查到被子底部，一根只插入水的表面，再在右杯中滴入几滴墨水并搅拌均匀，现在开始用酒精灯加热左边的烧杯，一段时间后就可以明显的看到染了颜色的水从右杯源源不断的流入左杯，左杯的水源源不断的流入右杯，最后两杯水都变成了墨水的颜色，与此同时用手摸右边的杯子，右边的水也热了起来，这就是冷热水发生了对流的缘故。 2.2.2 实验注意事项 短虹吸管只插入水的表面，不能过深。 玻璃管宜选壁较厚一些的，这样绝热性好一些，效果也好一些。 0=??+??x u C t u

2.2.3 实验原理分析 对左边的水杯用酒精灯加热，水受热密度变小开始上升，右边水杯的 冷水从下边的吸管流向左边的水杯进行补充，左边水杯的热水从上边的吸管流向右边的水杯，这样一会儿两杯水都变成墨水的颜色了[1]。 在冷水里面掺热水也是一样的道理，在不搅拌的情况下，最后水温基本都是一个温度，这就是水的对流，除了水的对流还有刮风是空气的对流，气压高的一方向气压低的一方补充空气，这就形成了对流，就会产生风；还有冬天在家里开空调，形成空气对流，最后整个房间的温度都升了起来。 2.3对流方程及其现有算法 1.针对常系数对流扩散方程，我们利用指数变换， 构造四阶紧致差分格式。 2.针对一维变系数对流扩散方程，将其转化为扩散方程，并构造四阶紧致差分格式。 3.对于常系数二维对流扩散方程，构造出四阶紧致差分方程，以及特殊的变系数 对流扩散方程的四阶紧致差分格式。 4.针对一维常系数对流扩散方程 和一维变系数对流扩散方程，分别构造了几种基于线性和双线性插值 的特征差分格式。 5.针对二维对流扩散方程 ，构造了几种基于线性和双线性插值的特征差分格式。 2.3影响物理量?的三个过程 用),,,(t z y x ??=来表示流体中单位体积的流体所携带的某种物理量，它可以是流体的质量或温度。流体的温度可以用?来表示，流体的密度ρ也可以用? ),(22t x f x u x u a t u +??=??+??ε2 2),(x u x u t x a t u ??=??+??εf y u x u a y u q x u p t u =??+??-??+??+??)(2222) ,()()()(2222y x f y u y q x u x p y u x u =??+??+??+??22x u x u v t u ??=??+??ε) ,()()()(22t x f x x a x u x b x u x c u =??-??+??) ,())(,(),(),(),(222221y x f y u x u y x a y u y x b x u y x b t u y x c =??+??-??+??+??

一维抛物线方程混合边界问题matlab求解(一维扩散方程)

在差分方法中分为前差、后差、中心差以及显式、隐式、中心式。这些概念分别对应着对空间和时间的离散。 针对以上抛物线方程的求解方法，分别采用向前、向后、对称六点和三行式进行计算。 一维抛物线的一般方程为 此题为混合边界问题 1、前差的一般格式为 将τ 成为网比，记做r，则差分格式可以写成 ?2

将上标为k+1的放在一边，k的放在一边，这样就可以写成矩阵形式。 根据已知条件便可以求解。（下面是前差的matlab代码） r=10; x=0:0.1:1; t=0:0.1:1; U=[]; U(:,1)=0; U(:,11)=0; U(1,:)=sin(pi.*x); for i=2:11 %%行 for j=2:10 %% 列 U(i,j)=r*U(i-1,j-1)+(1-2*r)*U(i-1,j)+r*U(i-1,j+1); end end figure; surf(x,t,U); 例题中并没有给出t的具体值，这里取了1，如同正方形laplace方程一样，纯粹为了好求。上面两个for循环代替了矩阵的作用。如果想改成矩阵，很简单，自行解决。

2、后差的一般格式为 这里请注意前差后差的不同点在于右式中上标的变化。 附上matlab代码： h=0.1; k=0.1; N=1/h; M=1/0.1; r=k/h^2; for i=1:N-1 B(i)=sin(i/N*pi); end U(:,1)=B; A=diag(ones(1,N-1)*(1+2*r))-diag(ones(1,N-2)*r,1)-diag(ones(1,N-2)*r, -1); for i=2:M+1 U(:,i)=A\U(:,i-1); end Z=[zeros(11,1),U',zeros(11,1)]; [X,Y]=meshgrid(linspace(1,0,11),linspace(1,0,11)); surf(X,Y,Z);

对流扩散方程.

A

对流扩散方程的求解 对流扩散问题的有效数值解法一直是计算数学中重要的研究内容，求解对流扩散方程的数值方法主要是有限差分法(FDM)、有限元法(FEM)、有限体积法(FVM)、有限解析法(FAM)、边界元法(BEM)、谱方法(SM) 等多种方法。但是对于对流占优问题，用通常的差分法或有限元法进行求解将出现数值震荡。 为了克服数值震荡，80年代，J.Douglas,Jr.和T.F.Russell 等提出特征修正技术求解对流扩散占优的对流扩散问题，与其它方法相结合，提出了特征有限元方法、特征有限差分方法、特征混合元方法；T.J.Hughes和A.Brooks提出过一种沿流线方向附加人工黏性的间断有限元法，称为流线扩散方法(SDM)。有限差分法、有限元法、有限体积法是工程应用中的主要方法。 对流扩散方程的特点 对流扩散方程右端第一项为扩散项，左端第二项则是对流项。由于其方程本身的特点，给建立准确有效的数值求解方法带来一定的困难。对流和扩散给流体中由流体携带的某种物理量的变化过程，可以通过一个无量纲的特征参数(Peclet数)来描述，Peclet数Pe的定义为：Pe=|ν|L/D。这里v是来流速度，L是特征长度，D是物质的扩散系数。如果Pe数较小，即对流效应相对较弱，这类问题中，扩散占主导地位，方程是椭圆型或抛物线型；

如果Pe数较大，即溶质分子的扩散相对于流体速度而言是缓慢的，这类问题中，对流占优，方程具有双曲型方程的特点。 对于对流占优问题的求解，采用常规的Galerkin有限元方法，为了避免求解结果产生数值振荡，获得稳定解，则应使每个单元的局部Peclet数，Peh=|ν|h/D≤2，这里h为单元的最大尺寸，|v|为单元中的最大速度分量值。因此，用本文方法求解对流占优对流扩散问题，要得到稳定解，则要通过加密有限元网格来实现。