《微分方程数值解》 课程论文 学生姓名1：许慧卿学号：20144329 学生姓名2：向裕学号：20144327学生姓名3：邱文林学号：20144349学生姓名4：高俊学号：20144305学生姓名5：赵禹恒学号：20144359学生姓名6：刘志刚学号： 20144346 学院：理学院 专业：14级信息与计算科学 指导教师：陈红斌 2017年6 月25日
《偏微分方程数值解》课程论文 《一维热传导方程的差分格式》论文 一、《微分方程数值解》课程论文的格式 1）引言：介绍研究问题的意义和现状 2）格式：给出数值格式 3）截断误差：给出数值格式的截断误差 4）数值例子：按所给数值格式给出数值例子 5）参考文献：论文所涉及的文献和教材 二、《微分方程数值解》课程论文的评分标准 1）文献综述：10分； 2）课题研究方案可行性：10分； 3）数值格式：20分； 4）数值格式的算法、流程图：10分； 5）数值格式的程序：10分； 6）论文撰写的条理性和完整性：10分； 7）论文工作量的大小及课题的难度：10分； 8）课程设计态度：10分； 9）独立性和创新性：10分。 评阅人： - 2 -
一维热传导方程的差分格式 1 引言 考虑如下一维非齐次热传导方程Dirichlet 初边值问题 22(,),u u a f x t t x ??=+?? ,c x d << 0,t T <≤ (1.1) (,0)(),u x x ?= ,c x d ≤≤ (1.2) (,)(),u c t t α= (,)(),u d t t β= 0t T <≤ (1.3) 的有限差分方法, 其中a 为正常数,(,),(),(), ()f x t x t t ?αβ为已知常数, ()(0),c ?α= ()(0).d ?β= 称(1.2)为初值条件, (1.3)为边值条件. 本文将给出(1.1) (1.3)的向前Euler 格式, 向后Euler 格式和Crank Nicolson -格式, 并给出其截断误差和数值例子. 经对比发现, Crank Nicolson -格式误差最小, 向前 Euler 格式次之, 向后Euler 格式误差最大. 2 差分格式的建立 2.1 向前Euler 格式 将区间[,]c d 作M 等分, 将[]0,T 作N 等分, 并记 ()/h d c M =-, /T N τ=, j x c jh =+,0j M ≤≤, k t k τ=,0k N ≤≤. 分别称h 和τ为空间步长和时间步长.用 两组平行直线 j x x =, 0j M ≤≤, k t t =, 0k N ≤≤ 将Ω分割成矩形网格.记{} |0h j x j M Ω=≤≤, {}|0k t k N τΩ=≤≤, h h ττΩ=Ω?Ω. 称() ,j k x t 为结点[1] . 定义h τΩ上的网格函数 {}|0,0k j U j M k N Ω=≤≤≤≤, 其中() ,k j j k U u x t =. 在结点() ,j k x t 处考虑方程(1.1),有
一维热传导方程 一． 问题介绍 考虑一维热传导方程： （1） ,0),(22T t x f x u a t u ≤<+??=?? 其中a 是正常数，)(x f 是给定的连续函数。按照定解条件的不同给法，可将方程（1）的定解问题分为两类： 第一类、初值问题（也称Cauthy 问题）：求具有所需次数偏微商的函数),(t x u ，满足方程（1）（∞<<∞-x ）和初始条件： （2） ),()0,(x x u ?= ∞<<∞-x 第二类、初边值问题（也称混合问题）：求具有所需次数偏微商的函数),(t x u ，满足方程（1）（l x <<0）和初始条件： （3） ),()0,(x x u ?= l x <<0 及边值条件 (4) .0),(),0(==t l u t u T t ≤≤0 假定)(x ?在相应区域光滑，并且在l x ,0=满足相容条件，使上述问题有唯一充分光滑的解。 二． 区域剖分 考虑边值问题（1），（4）的差分逼近。去空间步长N l h /=和时间步长M T /=τ，其中N,M 都是正整数。用两族平行直线： 将矩形域}0;0{T t l x G ≤≤≤≤=分割成矩形网格，网格节点为),(k j t x 。以h G 表示网格内点集合，即位于开矩形G 的网点集合；h G 表示所有位于闭矩形G 的网点集合；h Γ=h G --h G 是网格界点集合。 三． 离散格式 第k+1层值通过第k 层值明显表示出来，无需求解线性代数方程组，这样的格式称为显格式。 第k+1层值不能通过第k 层值明显表示出来，而由线性代数方程组确定，这样的格式称为隐格式。 1. 向前差分格式 （5） ,221 11j k j k j k j k j k j f h u u u a u u ++-=--++τ
有限差分法 finite difference method 用差分代替微分，是有限差分法的基本出发点。是一种微分方程和积分微分方程数值解的方 把连续的定解区域用有限个离散点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。 如何根据问题的特点将定解区域作网格剖分；如何把原微分方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式，为了有实用意义，一个基本要求是它们能够任意逼近微分方程，这就是相容性要求。另外，一个差分格式是否有用，最终要看差分方程的精确解能否任意逼近微分方程的解，这就是收敛性的概念。此外，还有一个重要的概念必须考虑，即差分格式的稳定性。因为差分格式的计算过程是逐层推进的，在计算第n＋1层的近似值时要用到第n层的近似值，直到与初始值有关。前面各层若有舍入误差，必然影响到后面各层的值，如果误差的影响越来越大，以致差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以控制的，就认为格式是稳定的。只有在这种情形，差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。 最常用的方法是数值微分法，比如用差商代替微商等。另一方法叫积分插值法，因为在实际问题中得出的微分方程常常反映物理上的某种守恒原理，一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的差分格式。
变系数线性常微分方程的求解 张慧敏，数学计算机科学学院 摘要：众所周知，所有的常系数一阶、二阶微分方程都是可解的，而变系数 二阶线性微分方程却很难解，至今还没有一个普遍方法。幂级数解法是一个非常有效的方法，本文重点讨论二阶变系数线性常微分方程的解法，从幂级数解法、降阶法、特殊函数法等方面探究了二阶微分方程的解法，简单的介绍了几种高阶微分方程的解法，并讨论了悬链线方程等历史名题。 关键词：变系数线性常微分方程；特殊函数；悬链线方程；幂级数解法 Solving linear ordinary differential equations with variable coefficients Huimin Zhang , School of Mathematics and Computer Science Abstract：As we know, all of ordinary differential equations of first, second order differential equations with constant coefficients are solvable. However, the linear differential equations of second order with variable coefficients are very difficult to solve. So far there is not a universal method. The method of power-series solution is a very efficient method. This article focuses on solving linear ordinary differential equations of second order with variable coefficients, and exploring the solution of in terms of power-series solution, the method of reducing orders, the method of special functions. Also, this paper applies the above methods to solve several linear differential equations of higher order and especially discusses the famous catenary equation. Key words：Linear ordinary differential equations with variable coefficients; Special Functions; catenary equation; Power Series Solution.
文献综述 信息与计算科学 热传导方程差分格式的收敛性和稳定性在实际研究物理问题过程中, 往往能给出问题相应的数学表达式, 但是由于实际物理问题的复杂性, 它的解却一般不容易求出. 由此计算物理应运而生, 计算物理是以计算机为工具, 应用数学的方法解决物理问题的一门应用性学科, 是物理、数学和计算机三者结合的交叉性学科. 它产生于二战期间美国对核武器的研究, 伴随着计算机的发展而发展. 计算物理的目的不仅仅是计算, 而是要通过计算来解释和发现新的物理规律. 这一点它与传统的实验物理和理论物理并无差别, 所不同的只是使用的工具和方法. 计算物理早已与实验物理和理论物理形成三足鼎立之势, 甚至有人提出它将成为现代物理大厦的“栋梁”. 在一个物理问题中一个数值解往往比一个式子更直观, 更有价值. 在实际求解方程时, 除了一些特殊的情况下可以方便地求得其精确解外, 在一般情况下, 当方程或定解条件具有比较复杂的形式, 或求解区域具有比较复杂的形状时, 往往求不到, 或不易求到其精确解. 这就需要我们去寻找方程的近似解, 特别是数值近似解, 简称数值解. 这里主要研究的是热传导方程. 有限差分法是微分方程和积分微分方程数值解的方法. 其基本思想是把连续的定解区域用有限个离散点构成的网格来代替, 这些离散点称作网格的节点;把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似;把原方程和定解条件中的微商用差商来近似, 积分用积分和来近似, 于是原微分方程和定解条件就近似地代之以代数方程组, 即有限差分方程组, 解此方程组就可以得到原问题在离散点上的近似解. 然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解. 热传导的差分法是求解热传导方程的重要方法之一. 对于差分格式的的求解, 我们首先要关注差分格式的收敛性和稳定性. 对于一个微分方程建立的各种差分格式, 为了有实用意义, 一个基本要求是它们能够任意逼近微分方程, 即相容性要求. 一个差分格式是否有用, 就要看差分方程的精确解能否任意逼近微分方程的解, 即收敛性的概念. 此外, 还有一个重要的概念必须考虑, 即差分格式的稳定性. 因为差分格式的计
向后差分格式MATLAB编程： c lear;clc; format short e a=input('请输入系数a的值'); l=input('请输入长度l的值'); M=input('请输入将区间[0,1]等分的个数M '); ot=input('请输入时间增量ot的值'); n=input('请输入运行次数n的值'); ox=1/M; x0=zeros(M+1,1) for ii=1:M x0(ii+1)=ii*ox; end u=sin(pi*x0/l); r=a*ot/(ox)^2; for ii=1:n %数据的输入 B=zeros(M-1,1); A=zeros(M-2,1); C=zeros(M-2,1); S=zeros(M-1,1); for ii=1:M-2 B(ii)=1+2*r;A(ii)=-r;C(ii)=-r; S(ii)=u(ii+1,1); end B(M-1,1)=1+2*r;S(M-1,1)=u(M,1);u(1,2)=0;u(M+1,2)=0; S(1,1)=S(1,1)+r*u(1,2);S(M-1,1)=S(M-1,1)+r*u(M+1,2); %追赶法 S(1)=S(1)/B(1);T=B(1);k=2; while k~=M B(k-1)=C(k-1)/T; T=B(k)-A(k-1)*B(k-1); S(k)=(S(k)-A(k-1)*S(k-1))/T; k=k+1 end k=1; while k~=M-1 S(M-1-k)=S(M-1-k)-B(M-1-k)*S(M-k); k=k+1; end u(2:M,2)=S; u(:,1)=u(:,2); end %计算精确解 for x=0:M
第五次作业（前三题写在作业纸上） 一、用有限差分方法求解一维非定常热传导方程，初始条件和边界条件见说明.pdf 文件，热扩散系数α=const ， 22T T t x α??=?? 1. 用Tylaor 展开法推导出FTCS 格式的差分方程 2. 讨论该方程的相容性和稳定性，并说明稳定性要求对求解差分方程的影响。 3. 说明该方程的类型和定解条件，如何在程序中实现这些定解条件。 4. 编写M 文件求解上述方程，并用适当的文字对程序做出说明。（部分由网络搜索得到，添加，修改后得到。） function rechuandaopde %以下所用数据，除了t 的范围我根据题目要求取到了20000，其余均从pdf 中得来 a=0.00001;%a 的取值 xspan=[0 1];%x 的取值范围 tspan=[0 20000];%t 的取值范围 ngrid=[100 10];%分割的份数，前面的是t 轴的，后面的是x 轴的 f=@(x)0;%初值 g1=@(t)100;%边界条件一 g2=@(t)100;%边界条件二 [T,x,t]=pdesolution(a,f,g1,g2,xspan,tspan,ngrid);%计算所调用的函数 [x,t]=meshgrid(x,t); mesh(x,t,T);%画图，并且把坐标轴名称改为x ，t ，T xlabel('x') ylabel('t') zlabel('T') T%输出温度矩阵 dt=tspan(2)/ngrid(1);%t 步长 h3000=3000/dt;
h9000=9000/dt; h15000=15000/dt;%3000,9000,15000下，温度分别在T矩阵的哪些行T3000=T(h3000,:) T9000=T(h9000,:) T15000=T(h15000,:)%输出三个时间下的温度分布 %不再对三个时间下的温度-长度曲线画图，其图像就是三维图的截面 %稳定性讨论,傅里叶级数法 dx=xspan(2)/ngrid(2);%x步长 sta=4*a*dt/(dx^2)*(sin(pi/2))^2; if sta>0,sta<2 fprintf('\n%s\n','有稳定性') else fprintf('\n%s\n','没有稳定性') error end %真实值计算 [xe,te,Te]=truesolution(a,f,g1,g2,xspan,tspan,ngrid); [xe,te]=meshgrid(xe,te); mesh(xe,te,Te);%画图，并且把坐标轴名称改为xe，te，Te xlabel('xe') ylabel('te') zlabel('Te') Te%输出温度矩阵 %误差计算 jmax=1/dx+1;%网格点数 [rms]=wuchajisuan(T,Te,jmax) rms%输出误差
向前差分格式MATLAB编程： c lear;clc; format short e a=input('请输入系数a的值'); l=input('请输入长度l的值'); M=input('请输入将区间[0,1]等分的个数M '); ot=input('请输入时间增量ot的值'); n=input('请输入运行次数n的值'); ox=1/M; x0=zeros(M+1,1) for ii=1:M x0(ii+1)=ii*ox; end u=sin(pi*x0/l); r=a*ot/(ox)^2; for ii=1:n %数据的输入 B=zeros(M-1,1); A=zeros(M-2,1); C=zeros(M-2,1); S=zeros(M-1,1); for ii=1:M-2 B(ii)=1+2*r;A(ii)=-r;C(ii)=-r; S(ii)=u(ii+1,1); end B(M-1,1)=1+2*r;S(M-1,1)=u(M,1);u(1,2)=0;u(M+1,2)=0; S(1,1)=S(1,1)+r*u(1,2);S(M-1,1)=S(M-1,1)+r*u(M+1,2); %追赶法 S(1)=S(1)/B(1);T=B(1);k=2; while k~=M B(k-1)=C(k-1)/T; T=B(k)-A(k-1)*B(k-1); S(k)=(S(k)-A(k-1)*S(k-1))/T; k=k+1 end k=1; while k~=M-1 S(M-1-k)=S(M-1-k)-B(M-1-k)*S(M-k); k=k+1; end D=(1-2*r)*eye(M-1); temp=r*linspace(1,1,M-2); D=D+diag(temp,1)+diag(temp,-1); S=D*S
一维热传导方程 一． 问题介绍 考虑一维热传导方程： （1） ,0),(22 T t x f x u a t u ≤<+??=?? 其中a 是正常数，)(x f 是给定的连续函数。按照定解条件的不同给法，可将方程（1）的定解问题分为两类： 第一类、初值问题（也称Cauthy 问题）：求具有所需次数偏微商的函数),(t x u ，满足方 程（1）（∞<<∞-x ）和初始条件： （2） ),()0,(x x u ?= ∞<<∞-x 第二类、初边值问题（也称混合问题）：求具有所需次数偏微商的函数),(t x u ，满足方 程（1）（l x <<0）和初始条件： （3） ),()0,(x x u ?= l x <<0 及边值条件 (4) .0),(),0(==t l u t u T t ≤≤0 假定)(x ?在相应区域光滑，并且在l x ,0=满足相容条件，使上述问题有唯一充分光滑 的解。 二． 区域剖分 考虑边值问题（1），（4）的差分逼近。去空间步长N l h /=和时间步长M T /=τ，其中N,M 都是正整数。用两族平行直线： ),,1,0(N j jh x x j === ),,1,0(M k k t t k ===τ 将矩形域}0;0{T t l x G ≤≤≤≤=分割成矩形网格，网格节点为),(k j t x 。以h G 表示网格内点集合，即位于开矩形G 的网点集合；h G 表示所有位于闭矩形G 的网点集合； h Γ=h G --h G 是网格界点集合。 三． 离散格式 第k+1层值通过第k 层值明显表示出来，无需求解线性代数方程组，这样的格式称为
几类二阶变系数常微分方程解法论文
二阶变系数常微分方程几种解法的探讨 胡博（111114109） （湖北工程学院数学与统计学院湖北孝感 432000） 摘要：常系数微分方程是我们目前可以完全解决的一类方程，而求变系数常微分方程的通解是比较困难的，一般的变系数常微分方程目前是还没有通用解法的。本文主要对二阶变系数常微分方程求解进行了探究，利用特解、常数变易法、变量变换等方法求出了某些二阶变系数线性微分方程的通解，并初步归纳了二阶变系数线性方程的求解基本方法及步骤。 关键词：二阶变系数线性微分方程；变换；通解；特解 To explore the solution of some ordinary differential equations of two order variable coefficient Zhang jun(111114128) (School of Mathematics and Statistics Hubei Engineering University Hubei Xiaogan 432000) Abstract:Differential equation with constant coefficients is a class of equations we can completely solve the present general solution, and change coefficient differential equations is difficult, the variable coefficient ordinary differential equation is at present there
一维热传导方程Last revision on 21 December 2020
一维热传导方程 一． 问题介绍 考虑一维热传导方程： （1） ,0),(22T t x f x u a t u ≤<+??=?? 其中a 是正常数，)(x f 是给定的连续函数。按照定解条件的不同给法，可将方程（1）的定解问题分为两类： 第一类、初值问题（也称Cauthy 问题）：求具有所需次数偏微商的函数),(t x u ，满足方程（1）（∞<<∞-x ）和初始条件： （2） ),()0,(x x u ?= ∞<<∞-x 第二类、初边值问题（也称混合问题）：求具有所需次数偏微商的函数),(t x u ，满足方程（1）（l x <<0）和初始条件： （3） ),()0,(x x u ?= l x <<0 及边值条件 (4) .0),(),0(==t l u t u T t ≤≤0 假定)(x ?在相应区域光滑，并且在l x ,0=满足相容条件，使上述问题有唯一充分光滑的解。 二． 区域剖分 考虑边值问题（1），（4）的差分逼近。去空间步长N l h /=和时间步长M T /=τ，其中N,M 都是正整数。用两族平行直线： 将矩形域}0;0{T t l x G ≤≤≤≤=分割成矩形网格，网格节点为),(k j t x 。以h G 表示网格内点集合，即位于开矩形G 的网点集合；h G 表示所有位于闭矩形G 的网点集合；Γ=G --G 是网格界点集合。
三． 离散格式 第k+1层值通过第k 层值明显表示出来，无需求解线性代数方程组，这样的格式称为显格式。 第k+1层值不能通过第k 层值明显表示出来，而由线性代数方程组确定，这样的格式称为隐格式。 1. 向前差分格式 （5） ,22111j k j k j k j k j k j f h u u u a u u ++-=--++τ )(j j x f f =， )(0 j j j x u ??==， 00==k N k u u ， 其中j = 1,2,…，N-1，k = 1,2,…,M-1。以2/h a r τ=表示网比。则方程（5）可以改写为： 易知向前差分格式是显格式。 2. 向后差分格式 （6） ,11111)21(j k j k j k j k j f u ru u u ru τ+=-++-+-+++ )(0 j j j x u ??==， 00==k N k u u ， 其中j = 1,2,…，N-1，k = 1,2,…,M-1，易知向前差分格式是显格式。 3. 六点对称格式（Grank-Nicolson 格式） 将向前差分格式和向后差分格式作算术平均，即得到六点对称格式： （7） 111112)1(2+-+++-++-k j k j k j u r u r u r =j k j k j k j f u r u r u r τ++-+-+112 )1(2 利用0j u 和边值便可逐层求到k j u 。六点对称格式是隐格式，由第k 层计算第k+1层时需解线性代数方程组（因系数矩阵严格对角占优，方程组可唯一求解）。
新疆大学毕业论文(设计) 题目:求解热传导方程的高精度隐式差分格式所属院系：数学与系统科学学院 专业：信息与计算科学
声明 本人郑重声明该毕业论文（设计）是本人在开依沙尔老师指导下独立完成的，本人拥有自主知识产权，没有抄袭、剽窃他人成果，由此造成的知识产权纠纷由本人负责。 声明人（签名）： 年月日 亚库甫江.买买提同学在指导老师的指导下，按照任务书的内容，独立完成了该毕业论文（设计），指导教师已经详细审阅该毕业论文（设计）。 指导教师（签名）： 年月日
新疆大学 毕业论文（设计）任务书 班级：信计07-2 姓名：亚库甫江.买买提论文（设计）题目：求解热传导方程的高精度隐式差分格式 专题：毕业设计 论文（设计）来源：教师自拟 要求完成的内容：学习和掌握一维热传导方程已有的各种差分 格式的基础上，扩散方程对空间变量应用紧 致格式离散，对时间变量应用梯形方法，构 造热传导方程的精度为() 24 τ+数值格式， O h 讨论格式的稳定性，最后数值例子来验证。发题日期：2012 年12月25日完成日期：2012 年5月28 日实习实训单位：数学学院地点：数学学院 论文页数：19页；图纸张数：4 指导教师：开依沙尔老师 教研室主任 院长（系主任）
摘要 本文首先对热传导方程经典差分格式进行复习和讨论，然后热传导方程对空间变量四阶紧致格式进行离散,时间变量保持不变，把一维热传导方程转化为常微分方程组的初值问题, 再利用梯形方法构造热传导方程方程的时间二阶空间四阶精度的一种差分格式，并稳定性进行分析，数值结果与Crank-Nicholson 格式进行比较，数值结果表明, 该方法是有效求解热传导方程的数值计算. 关键词: 热传导方程,高精度紧致格式; 梯形方法；两层隐格式; Crank-Nicolson格式 ABSTRACT This paper first study on some classical finite difference for the heat conduction equation, secondely secondely we apply compact finite difference approximation of fourth order for discretizing spatial derivatives but leave the time variable Continuous. This approach results in a system of ODEs, which can then be used trapezodial formula derived fourth order in space and second order in time unconditionally stable implicit scheme .the stability and local truncation error of the obtained method are analysied. Numerical experiments shows that this method Useful, efficient method for solving diffusion equation Keywords: Heat conduction eqution;Higher- oder compact scheme; Trapezodial formula ;Two- level implict scheme; Crank- Nicolson scheme
有限差分和pde 函数求解一维定态热传导方程 分别用有限差分方法和pde 函数求解一维定态热传导方程，初始条件和边界条件，热扩散系数α=0.00001， 22 T T t x α??=?? （1） 求解过程： 1. 用Tylaor 展开法推导出FTCS 格式的差分方程 首先对T 进行泰勒展开得到如下两式子： 2 3 1231 2 3 ... 232! 3! 2 3 ... 232!3!n n n n n j j j j j n n n n n j j j j j t t T T t x x T T x T T T t t t T T T x x x ++??=+?+ + +??=+?+ + +????????? ? ? ?????? ???? ????????? ? ? ?????? ?? ?? 上述两个方程变换得： ()11223 23...23n n n n n n n j j j j j j j T T T T T t T t T o t t t t t t ++--???? ???????= --=+? ? ? ???????????? （2） 223 123...23n n n n n j j j j j T T T x T x T x x x x --???? ???????= -- ? ? ??????????? ()1232422 342222...3!4!n n n n n n j j j j j j T T T T x T x T x x x x x x +-?? ????????????=--- ? ? ? ??????????????? ()()2112 22 22-n n n j j j T T T T o x x x +--+???=+? ????? （3） 将上述式子（2）（3）代入（1）得：
二维变系数抛物型方程的一个高阶ADI 差分格式 马小霞1，颜晓琳2，陈汝栋2 （1.焦作大学基础部，河南焦作 454003；2.天津工业大学理学院，天津 300387） 摘要：针对二维变系数抛物型方程，构造出了一个高精度、恒稳定的交替方向隐式（ADI ）差分格式，格式的截断误 差阶达O （τ2+h 4）.通过数值实验，验证了理论分析的正确性和差分格式的精确性与有效性. 关键词：抛物型方程；ADI 格式；截断误差；恒稳定中图分类号：O241.82 文献标志码：A 文章编号：１６７１－０２４Ｘ（２０14）０1－００77－０4 A high accuracy ADI difference scheme for solving two-dimension variable coefficients parabolic equation MA Xiao-xia 1，YAN Xiao-lin 2，CHEN Ru-dong 2 （1.Department of Basic Course ，Jiaozuo University ，Jiaozuo 454003，China ；2.School of Science ，Tianjin Polytechnic University ，Tianjin 300387，China ） Abstract ：A high accuracy alternation direction implicit scheme （ADI ）for solving the two-dimensional parabolic equations is presented ，and the scheme is absolutely stable and the truncation error is O （τ2+h 4）.The experiments show the scheme is effective and advantage ，and the theory is right by a numerical example. Key words ：parabolic equation ；ADI difference scheme ；truncation error ；absolutely stable 收稿日期：2013-05-31 基金项目：国家自然科学基金（11071279）；河南省教育厅自然科学基础研究基金（2008B110016） 第一作者：马小霞（1969—），女，硕士，讲师. 通信作者：陈汝栋（1956—），男，教授，硕士生导师.E-mail ：chenrd@https://www.doczj.com/doc/d815748121.html, 天津工业大学学报 ＪＯＵＲＮＡＬＯＦＴＩＡＮＪＩＮＰＯＬＹＴＥＣＨＮＩＣＵＮＩＶＥＲＳＩＴＹ 第33卷第1期２０14年2月 Vol.33No.1February 2014 抛物型方程在处理废料污染、渗透、驱动、海水入侵以及半导体等工程实际问题中有着广泛的应用，因此研究其高精度、高稳定和计算量较小的数值解法具有重要的意义.用有限差分方法研究这类问题的数值方法目前已做了许多工作[1-5].但这些工作大多是对常系数而言的.文献[4]中对二维变系数抛物型方程数值方法仅对系数依赖于一个变量的情况进行了研究，本文的研究是对系数依赖于两个变量的情形进行的.应用Taylor 展开、算子方法[6]以及粘结系数法[7]得到了一个高精度（截断误差阶达O （τ2+h 4））、恒稳定的ADI 格式.格式的建立和稳定性分析都比文献[4]简单得多，文末的数值实验证明了本文理论分析的正确性和所得格式的精确性与有效性. 1差分格式的建立 考虑如下的二维变系数非齐次抛物型方程初边 值问题 鄣u 鄣t =a （x ，y ）鄣2u 鄣x 2+b （x ，y ）鄣2 u 鄣y 2 +f （x ，y ，t ）（x ，y ，t ）∈Ω×（0，T ]（1）u （x ，y ，0）=φ（x ，y ） （x ，y ）∈Ω軍（2） u （x ，y ，t ）=Ψ（x ，y ，t ） （ x ，y ）∈Γ，0，0＜t ≤T （3≤≤≤≤≤≤≤≤≤≤≤≤≤≤≤≤≤≤≤≤≤≤≤≤≤ ）其中：0＜c 1≤a （x ，y ）≤c 2；Γ为Ω的边界. 设τ=Δt =T/N 为时间步长，h =Δx =Δy =1/M 为空间步长，N 、M 均为正整数.u n j ，k 表示在节点（jh ，kh ，n τ）处的网函数值，微分方程问题（1）—（3）的解函数为u （x ，y ，t ），并记u （jh ，kh ，n τ）=u （j ，k ，n ），f n+12 j ，k = 12 （f n +1j ，k +f n j ，k ）.由Taylor 展开式u （j ，k ，n +1）=u （j ，k ，n ）+τ鄣u （j ，k ，n ）+τ22鄣2u （j ，k ，n ）鄣t 2+…=exp （τ鄣鄣t ）u （j ，k ，n ）（4）
用向前差分格式求解二维热传导方程function varargout=liu(varargin) T=1;a=1;h=1/30;dt=1/150; [X,T,Z]=chfenmethed(h,dt,a,T); mesh(X,T,Z(:,:,3)); shading flat; % xlabel('X','FontSize',14); % ylabel('t','FontSize',14); % zlabel('error','FontSize',14); % title('误差图'); function [X,Y,Z]=chfenmethed(h,dt,a,T); %求解下问题 %u_t-a*(u_xx+u_yy)=f(x,y,t) 0 n=length(t); r=a*dt/h^2; [X,Y]=meshgrid(x,y); Z=zeros(m,m,n); U=zeros(m,m,n); for i=1:m for j=1:m U(i,j,1)=d(x(i),y(j)); end end for j=2:n for k=1:m U(1,k,j)=g0(y(k),t(j)); U(m,k,j)=g1(y(k),t(j)); U(k,1,j)=h0(x(k),t(j)); U(k,m,j)=h1(x(k),t(j)); end end for k=2:n for i=2:m-1 for j=2:m-1 微分方程数值解II 主要内容： 第一章有限差分法的理论基础 1. 构造差分格式的主要方法； 2. 差分格式的一般性要求； 3. Lax等价性定理； 4. 差分格式的von Neumann稳定性分析方法； 5. 差分格式的修正方程。 第二章线性抛物型方程的差分方法 1. 扩散方程的显式格式； 2. 扩散方程的隐式格式； 3. 线方法； 4. 多维抛物型方程的ADI方法； 5. 分数步法； 6. Burgers方程的差分法和网格雷诺数。 第三章一维线性双曲型方程的数值方法 1. 线性双曲型系统的特征和Riemann问题； 2. 守恒律的有限体积法； 3. Lax-Friedriches格式、Lax-Wendroff格式、特征线法差分格式； 4. 双曲型方程的迎风格式、CIR格式、Godunov 方法； 5. 二阶Godunov格式、总变差概念及限制器函数； 6. 双曲型方程及变系数双曲型方程的高分辨率(TVD)波传播格式。 第四章一维非线性双曲型守恒律的数值方法 1. 非线性双曲型守恒律的间断解、弱解、熵条件； 2. 标量守恒律的Riemann问题解及Godunov格式； 3. 熵修正、数值粘性、Osher格式及高分辨率波传播格式； 4. 守恒型与Lax-Wendroff定理、离散熵条件、非线性稳定性及收敛性； 5. 典型守恒律方程组的Godunov间断分解方法及Godunov格式； 6. 守恒律方程组的MUSCL格式。 第五章多维双曲型守恒律的高分辨率格式 1. 多维方程组的双曲性； 2．Lax-Wendroff方法、Runge-Kutta推进的半离散方法、维数分裂方法； 3. 标量方程的LW方法、Godunov 格式、方向迎风及角迎风格式； 4. 多维标量方程的高分辨率格式； 5. 多维方程组的高分辨率格式。 第六章双曲型守恒律的其它高分辨率方法 1. ENO与WENO格式； 一维抛物方程的初边值问题 分别用向前差分格式、向后差分格式、六点对称格式，求解下列问题： 22,01,u u a x t x ??=< 在0.05,0.10.2t =和时刻的数值解，并与解析解2 (,)sin()t u x t e x ππ-=进行比较。 1差分格式形式 设空间步长1/h N =, 时间步长0τ>，T M τ=，网比2/r h τ=. （1）向前差分格式 该问题是第二类初边值问题（混合问题），我们要求出所需次数的偏微商的函数 (,)u x t ，满足方程22,01,u u a x t x ??=<。 已知sin x π在相应区域光滑，并且在0,x l =与边值相容，使问题有唯一充分光滑的 解。 取空间步长1/h N =，和时间步长/T M τ=，其中,N M 都是正整数。用两族平行直 线 (0,1,,) j x x jh j N ===L 和 (0,1,,) k t t k k M τ===L 将矩形域 {01,0}G x t =≤≤≥分割成矩形网络，网络格节点为(,)j k x t 。以h G 表示网格内点集合， 即位于矩形G 的网点集合；h G 表示闭矩形G 的网格集合；h h G G -=Γh 是网格界点的集合。 向前差分格式，即 i k j k j k j k j k j f h u u u a u u ++-=--++2 1 112τ （1） KINHIE Fluent菜鸟指南 Fluent 160问 Eric 2009-12-8 如何入门 学习任何一个软件，对于每一个人来说，都存在入门的时期。认真勤学是必须的，什么是最好的学习方法，我也不能妄加定论，在此，我愿意将我三年前入门FLUENT心得介绍一下，希望能给学习FLUENT的新手一点帮助。 由于当时我需要学习FLUENT来做毕业设计，老师给了我一本书，韩占忠的《FLUENT流体工程仿真计算实例与应用》，当然，学这本书之前必须要有两个条件： 第一，具有流体力学的基础，第二，有FLUENT安装软件可以应用。然后就照着书上二维的计算例子，一个例子，一个步骤地去学习，然后学习三维，再针对具体你所遇到的项目进行针对性的计算。不能急于求成，从前处理器GAMBIT，到通过FLUENT进行仿真，再到后处理，如TECPLOT，进行循序渐进的学习，坚持，效果是非常显著的。如果身边有懂得FLUENT的老师，那么遇到问题向老师请教是最有效的方法，碰到不懂的问题也可以上网或者查找相关书籍来得到答案。另外我还有本《计算流体动力学分析》王福军的，两者结合起来学习效果更好。 CFD计算中涉及到的流体及流动的基本概念和术语 理想流体和粘性流体；牛顿流体和非牛顿流体；可压缩流体和不可压缩流体；层流和湍流；定常流动和非定常流动；亚音速与超音速流动；热传导和扩散等。 理想流体（Ideal Fluid）和粘性流体（Viscous Fluid） 流体在静止时虽不能承受切应力，但在运动时，对相邻的两层流体间的相对运动，即相对滑动速度却是有抵抗的，这种抵抗力称为粘性应力。流体所具备的这种抵抗两层流体相对滑动速度，或普遍说来抵抗变形的性质称为粘性。粘性的大小依赖于流体的性质，并显著地随温度变化。实验表明，粘性应力的大小与粘性及相对速度成正比。当流体的粘性较小（实际上最重要的流体如空气、水等的粘性都是很小的），运动的相对速度也不大时，所产生的粘性应力比起其他类型的力如惯性力可忽略不计。此时我们可以近似地把流体看成无粘性的，这样的流体称为理想流体。十分明显，理想流体对于切向变形没有任何抗拒能力。这样对于粘性而言，我们可以将流体分为理想流体和粘性流体两大类。应该强调指出，微分方程数值解II
热传导方程地差分格式
Fluent菜鸟入门(已经排版)
相关主题
文本预览