基于有限差分法的微分方程离散化求解.pdf

《基于有限差分法的微分方程离散化求解.pdf》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、基于有限差分法的微分方程离散化求解◆白红杰白阳椿韩自源(四川警察学院)【摘要】目前偏微分方程数值求解的方法主要有两种，即有限差分是用它来逼近一阶导数。法和有限元方法。本文论述了基于有限差分法的微分方程求解，离2．2二阶偏导数的差分逼近散化过程，并对结果进行了分析。在渗流微分方程中，除含有求知函数的一阶偏导数外，还含有【关键词】有限差分法离散化数值模拟未知函数的二阶偏导数。二阶偏导数扩“／良。相当于对函数一阶偏导数afalf胁)／既。1．前言在x轴上取三点：X—Ax、x、x+AX，如图1所示。令Oul~x：F有限差分法是计算机数值模拟最早采用比较成熟的方法，该磐．．；方法将求解域划分为差分网格，

2、用有限个网格节点代替连续的求x一△—一ix+解域，是一种直接将微分问题变为代数问题的近似数值解法，数学。竺2概念直观，表述简单。有限元方法的基础是变分原理和加权余量图1二阶差商处理方法示意法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在将aF，缸在点X处展开为中心差商每个单元内必改写成由各变量或其导数的节点值与所选用的插值冬一函数组成的线性表达式，借助于变分原理或加权余量法，将微分方—(5)程离散求解。再将缸／觑zF分别在x—Ax／2，x+Ax／2处展成一阶中心差用有限差分法求解偏微分方程必须把连续问题进行离散化，商，得F一l一一为此首先要对求解区域进行离散化。构造离散网格系统的目的

3、在{—一-_—凸_(6)于将表现为非均系统的大尺度用若干可以近似为均匀系统的尺度}～也一i(如网格)表征。构造差分形式就是对参数在一定的离散点中心网一格或块中心网格上离散。其中，离散网格可以是空间离散网格，也x专一：==(7)可以是时间离散网格(即离散时间步长)。平面网格的形式是多种将式(6)、(7)代入(5)得多样的，如矩形网格、柱形网格、多边形网格等。空间离散网格是和⋯垒二一二垄一，、～、边界条件相关联的，一般来说，对于第一类边界条件采用点中心网盘盘一。(8)格较方便，第二、三类边界条件采用块中心网格比较合适。式(8)即是以二阶差商来代替二阶偏导数的公式，这种形成的2．微分方程的离散化．差

4、商称为二阶中心差商。2．1一阶偏导数的差商逼近2．3抛物型方程的离散化方法设有函数u(x，Y，t)，对其自变量X的偏导数可以表示成函数在各类低渗透介质渗流问题中，大多属于不稳定渗流，描述这的差商的极限形式类渗流问题的数学方程不仅含有未知量对空间坐标的偏导数，而t)：lira—u(x+A％y~t)-u(x,y,t)且含有各变量对时间坐标的偏导数。这类方程在数学上属抛物型—，，盘‘△x(1)方程。对这类问题求解时，不仅要求出各变量在空间的分布状态，在(1)式中，当自变量增量充分小时，如果能用比较简单的函数而且还要求出这种分布在时间区域上的变化情况。对数值解法来差商来代替偏导数，即说，就是要求出各变

5、量在不同时刻在各空间点上的数值分布。8uu(x+Ax，Y，t)一“(，)’，f)因为这类方程所描述的是不稳定渗流过程，所以方程中的各～瓦————一(2)个变量都是随时间而变的，而空间差分项本身并没有指明它应该这样就可以把偏微分方程用差分方程代替，从而把难于求解在什么时间点上取值，因此对应于不同的时间取值方法就可以得的偏微分方程化成代数方程组。利用Taylor级数可以说用(2)形到不同形式的差分方程。更主要的是，由于求解这类问题时得到的式的差商来逼近一阶偏听偏信导数，其误差对AX来说是一阶的。是一系列时间点上的参数分布，后一个时间点上的结果依赖于前式(2)是用前差商来代替一阶偏导数即一个时间点的

6、结果，这就产生了求解误差在时间坐标上的传播与Ouu(x,y,t)-u(x-A~,y,t)变化的问题，即所谓的解的稳定性问题。————一(3)取最简单的一维方程进行讨论：同理，后差商也可以用来代替一阶偏导数，且其误差也为0：o≤fo(Ax)。类似地，还可以用中心差商(9)x+等，卜卜等0其初始条件为u(x，O)=)———————一(4)边界条件为，f)：(f)来代替偏导数Ou／Ox。经证明，用中心差商代替一阶偏导数，其，O=uz∞误差为O(Ax2)。这就意味着用中心差商来逼近偏导数，当Ax非式(9)左端是对空间坐标的偏导数，右端是对时间坐标的偏导常小时，是一种比前两种差商逼近更为精确的方法。但是

7、，在具体数。为了对时间的偏导数离散化，取时间坐标上两个时刻的步长为求解一个微分方程时，不仅要考虑差商对导数的逼近精度，更主要At。对方程(9)的右端以任意时刻t“的后差商进行近似，在任意点的是把差分方程作为一个整体考虑。因此，油藏数值模拟中也不总xj得到165万方数据配一2+一任何误差不在以后的进一步计算中被扩大。因为一旦这种误差能一一广一一—(10)被不断扩大而成为越来越大的误差，必然使所得的结