内科大冶金传输原理课程程序设计上机指导及说明书.docVIP

冶金传输原理课程程序设计上机指导及说明书

冶金传输原理数值计算部分共计8学时，包含有：

1.速度边界层计算，2学时

2.二维稳态导热计算，2学时

3.一维非稳态导热计算（隐式差分方程），2学时

4.二维非稳态导热计算（显式差分方程），2学时

边界层微分方程数值解法

1、目的和意义

通过该程序设计，让学生掌握常微分方程初值问题数值解的求解方法及程序设计。

2、方程的推导

根据普朗特微分方程和布拉休斯解，可以得到如下方程：

运用常微分方程的数值解法进行求解。为了方便，在推导时用函数y代替，自变量x代替，将上述方程改写为：

令，则上列方程可以变为：

利用四阶标准Longe-Ku