数值计算的误差分析(一维显式格式)

1、帕克莱数

首先介绍一下何为帕克莱数(Pe数):佩克莱数(Peclet number,简称Pe数)是一个在连续传输现象研究中经常用到的一个无量纲数。其物理意义为对流速率与扩散速率之比,其中扩散速率是指在一定浓度梯度驱使下的扩散速率。

Pe=v/u=vL/α 其中v为特征速度,L为特征长度,α为特征扩散系数.u=α/L

由上式可以看出帕克莱数越大,对流运输占比越大。

对流-扩散方程的离散过程中,扩散项二阶导数的中心差分离散格式有优良的物理特性和计算精度。

2、正系数法则

正系数法则的提出是为了离散格式的稳定,避免出现随着迭代的进程,节点计算值的波动,这个与离散格式的稳定性有关。

非正系数下的计算错误

3、初值问题的稳定性

    非稳态问题中,计算是由着初始条件一步步的向前推进,初始条件给的误差以及各个时层计算引入的误差,如果随着计算时间的推移初始条件引入的误差不断的被放大则物理问题的解会被破坏,得不到正确解,此时初始问题不具有稳定性。初始问题的稳定性是一个离散格式的固有属性。这一思想也可以应用在离散格式对待误差上,格式的稳定性也是看误差的随时间的变化。

    误差随时间的传递问题不介绍了,主要是由n个方程(显式格式)组合成的矩阵A,既是求解的方程组的矩阵也是误差传递的矩阵(边值为0条件下)。

    误差矢量在计算过程中的传递,可以将误差矢量分解成若干个不同频率的分量的迭加组成,离散点的误差既是各个分量的迭加结果,将其中一个分量进行分析[1]。

初值稳定性-1
初值稳定性-2

4、离散格式的迁移性

中心差分格式与迎风格式的差异体现在对对流项的差分上,中心差分格式以本节点前后的节点数值进行差分,迎风(背风)格式以来流方向的节点值与本节点的数值进行差分。

守恒方程中对流项的中心差分格式不具有迁移性,即它使得扰动向两边扩散,且进行微分方程离散时不能保证正系数法则,迎风差分格式时具有迁移性,扰动随着对流的方向传递,且在满足一定条件可以保证正系数法则。Pe数比较大,可以看出对流作用比较剧烈时,且网格划分受到限制时,慎用中心差分格式。

离散方程还需要具有守恒性:1、控制方程守恒;2、物理量及其导数是连续的。

一阶迎风格式稳定性好,看起来结果合理,但是误差大,基本被限制。

差分格式稳定性分析的目的是为了找到允许的最大时间步长。

库朗数Courant number实际上是指时间步长和空间步长的相对关系。

c=\frac{u\Delta t}{\Delta x}

Courant数决定时间步长与网格尺寸,Pe数决定离散格式?

[1] 陶文铨编著.数值传热学.西安:西安交通大学出版社[M],2001.59-60.

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • 何谓数值分析? 众所周知现实生活中科学技术上面的问题大多数都能够通过建立对应的数学模型把实际问题转化成一个数学问题...
    罗泽坤阅读 8,905评论 0 14
  • 非守恒控制方程的推导 将瞬态项和对流项中的物理量从微分符号中移除 有限差分法 向前(i+1,i)、向后(i,i-1...
    Jonah_Ren阅读 615评论 0 0
  • 计算机控制系统 学习思路 信号变换——系统模型的描述(数学描述、性能分析)——控制器的设计(基于传递函数的间接、直...
    Invoke_r阅读 2,988评论 0 0
  • 要说明白CFD计算中的人工粘性,便没有办法绕开其他几个概念,他们分别是: 截断误差 舍入误差 修正方程 先说舍入误...
    唉那个什么阅读 3,858评论 1 1
  • 我好想你 昨夜星辰昨夜风像你 天边树若荠蒙蔽思你 眼前人是心上人是你 东风破我嘲笑我亦你 我的心好想你
    野原大阅读 189评论 0 0