计算方法 偏微分方程数值解.pptxVIP

1;本章要求;5. 1 偏微分方程简介;5.1 偏微分方程简介;5.1 偏微分方程简介;5.1 微分方程的求解思路;5.2 离散化公式;5.2 离散化公式;5.2 离散化公式推导;几种常见偏微分方程的离散化计算;波动方程求解;波动方程求解;波动方程求解;一维流动热传导方程;将上式进行处理得到： 该式是显式格式。只要保证式中各项系数大于零，一般情况下是稳定的，可以获得稳定的解。 分析上式可以发现，当为了提高数值精度取适当小的Δx 时，最有可能小于零的系数是 uin的系数，若要保证此项系数大于零，此时Δt必须相应地更小,会导致计算量将大大增加，这是显式格式的缺点，为了克服此缺点，下面提出一种隐式格式： 偏微分方程在 点上进行离散化，且对时间的偏微分采用向后欧拉公式得到原偏微分方程的离散化公式：; 从图5-3中可见要由初值及边界条件一排一排推上去是不行的，需解线性方程组,同时添上二边界条件： 正好共有m+2个方程，同时有m+2个变量，就能解出n+1排上各点值。这样，每解一个线性方程组，就可以往上推算一排点的u值，虽然引入了方程组的求解，有可能增加计算量，但由于隐式格式无条件稳定，Δt的取法与Δx 无关，可以少计算许多排节点上的u 值，相应于显式格式来说，最终