流动传热燃烧七的数值计算.pptVIP

： ： ： ： j Δ x 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 n Δ t n j 1 2 3 4 5 6 7 8 9 10 0 0 0.5 1 ε 1.0 2 1/2ε 0 1/2ε 1.5 3 1/4ε 0 1/4ε 0 1/4ε 2.0 4 1/8ε 0 3/8ε 0 3/8ε 0 1/8ε 2.5 5 0 1/4ε 0 3/8ε 0 1/4ε 0 1/16ε 0 : : j Δ x 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 n Δ t n j 1 2 3 4 5 6 7 8 9 10 0 0 0.5 1 ε 1.0 2 ε - ε ε 1.5 3 ε -2 ε 3 ε -2 ε ε 2.0 4 ε -3 ε 6 ε -7 ε 6 ε -3 ε ε 2.5 5 -4 ε 10ε -16 ε 19 ε -16 ε 10 ε -4 ε ε 2.Von Neumann方法（傅里叶级数法、分离变量法） 优点：寻求一般稳定性判据最常用的方法，使用方便 可靠。 局限性：只能在常系数的线性初值问题中求稳定的充分必要条件，在实际变系数的非线性问题以及各种不同边界条件的问题中，应用受到限制。 ： ： ： §1.1.6 收敛性与稳定性的关系，Lax等价定理 收敛性和稳定性是两个不同的概念，分别属于离散化过程和离散代数方程的求解过程。 收敛性分析比稳定性分析更困难，而收敛性与稳定性有密切联系。 Lax等价定理：对于一个适定的线性偏微分方程的初值问题和它的一个具有相容性的差分方程来说，稳定性是收敛性的必要和充分条件。 适定的初值问题：解存在、唯一、并且连续依赖于初值。 §1.1.7 非线性问题的有限差分法 一般求解方法： 1°给出所有结点上T的试探值 2°用试探值计算离散方程中的系数 3°求解名义上的线性方程组，得到新的T值 4°用新的T值作为较好的试探值→2° 5°若T值不再产生大的变化→收敛，否则发散 §1.1.8 关于对流项的迎风差分（upwind difference scheme） 1.迎风差分的基本思想 迎风—迎着来流（从上游）获取信息，来构造对流项的离散格式。 ? 2.对流项离散格式的迁移性 如果对流项的某种离散格式能使扰动仅沿着流动方向传递，则称此离散格式具有迁移特性。 3.讨论 1°迁移性是对流项离散格式的一个重要的物理特性，凡是由不具有迁移性的对流项离散格式组成的离散方程，在计算中可能出现解的振荡，因而只是有条件稳定。 2°虽然对流项的中心差分为二阶精度而迎风差分为一阶精度，但就它们对物理特性的模拟而言，迎风差分要比中心差分更合理。所以在求解实际物理问题时，只注意差分格式的精度是不够的。 3°在数值解不出现振荡的范围内，对流项采用中心差分比采用一阶迎风格式得到的解精度更高；但是，当对流作用十分强烈、网格数又受到限制时，采用中心差分会出现解的振荡，而采用一阶迎风格式却始终可以得到物理上看起来合理的解（尽管可能有较大误差）。 4°数值计算实践已表明，一阶迎风格式可能使计算结果产生比较严重的误差（数值扩散或假扩散），但从一阶迎风格式发展而来的思想已被用于构造具有二阶或三阶精度的迎风格式。 由于一阶迎风格式优良的数值稳定性，它往往用于初步计算，起过渡作用。 第一章：离散化方法 离散化就是将原微分方程在整个求解域内所考虑的物理量求解问题，转化为对求解域内有限个位置（网格结点）上物理量的求解问题。用各网格结点上定义的离散状态量来代替在求解域内连续分布的状态量。 ?这样，将原来的微分方程连同初、边界条件一起，离散成为一组代数方程。 对于给定的微分方程，对应有不同的离散化方程，它们之间的差别是由于所取的分布以及推导离散化方程的方法不同造成的。 §1.1有限差分法 差分法的基本思想是用差商代替微商。 用差分法离散微分方程的具体方法有很多种，对同一微分方程可以建立不同的差分方程。 §1.1.1差分方程的建立 以一维问题为例，设x3-x2=x2-x1=Δx 由Taylor级数展开可得： 则 可以表示为： 236 则 可以表示为： §1.1.2网格划分 离散过程包括两个方面：①微分方程的离散化； ②求解域的离散化。 设微分方程有两个自变量x，t 。 微分方程的解应是x和t的函数u(x,t)，而差分法（数值方法）求的是若干个离散点（x=j Δx，t=nΔt ）上u(x,t)的近似值 因此，首先要在求解域上指定一系列离散点。 一般取时间步长Δt =const.，空间步长Δx =const. 对于二维和三维问题，不要求Δx = Δy= Δ