计算差分方程的收敛点_PDE有限差分方法(1)——全显格式与全隐格式

4a17c26890b48a4e2c409c259185b25c.png

今天开始更新偏微分方程的有限差分方法!!!!!预计两周三更(考试期间会更得少一点)

参考书:
(1) J.W. Thomas - Numerical Partial Differential Equations_ Finite Difference Methods (1995, Springer)
(2) 张强《偏微分方程的有限差分方法》科学出版社,2019年1月版。
(3) K. W. Morton, D. F. Mayers - Numerical solution of partial differential equations (2005, Cambridge University Press)

本次内容:

cc201ceb4d23e2336a49277c083b9b8e.png

0 (HD)问题

是给定的终止时刻, 考虑一维热传导方程的边值问题(HD)

其中扩散系数

是给定常数,
是已知源项. 相应的定解条件是初值

和边值

其中

是已知函数.

根据偏微分方程理论,模型问题(HD)是适定的, 即真解

唯一存在, 且连续依赖于定解条件. 在适当的条件下, 真解充分光滑的, 其任意阶导数都是连续有界的.

模型问题的真解可以用公式准确表达, 比如分离变量法、Duhamel原理. 但是推导过程繁琐, 公式通常是级数或者卷积, 不适合大量位置的函数值计算; 扩散系数是常数, 无法推广.

所以,有必要考虑求模型问题的数值解。在这里,主要介绍有限差分方法来求数值解。

1 初边值问题基本实现过程

基本实现过程:离散网格

微分方程的离散
定解条件的离散.

一、离散网格

通常采用最简单的矩形结构. (接下来的内容大多数都采用它)

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
有限分初学者必备】如何根据问题的特点将定解区域作网格剖分;如何把原微分方程离散化为方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确,还需从理论上分析方程组的性态,包括解的唯一性、存在性和格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种格式,为了有实用意义,一个基本要求是它们能够任意逼近微分方程,这就是相容性要求。另外,一个格式是否有用,最终要看方程的精确解能否任意逼近微分方程的解,这就是收敛性的概念。此外,还有一个重要的概念必须考虑,即格式的稳定性。因为格式计算过程是逐层推进的,在计算第n+1层的近似值时要用到第n层的近似值 ,直到与初始值有关。前面各层若有舍入误,必然影响到后面各层的值,如果误的影响越来越大,以致格式的精确解的面貌完全被掩盖,这种格式是不稳定的,相反如果误的传播是可以控制的,就认为格式是稳定的。只有在这种情形,格式在实际计算中的近似解才可能任意逼近方程的精确解。关于格式的构造一般有以下3种方法。最常用的方法是数值微分法,比如用商代替微商等。另一方法叫积分插值法,因为在实际问题中得出的微分方程常常反映物理上的某种守恒原理,一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的格式
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值