120 重庆理工大学学报 外部闭曲线划分为60个单元，边界总剖分为76 r(1,0) -10≤x≤10，y=5 (u,)= 单元，由Matlab剖分在矩形区域内确定1521个 (0,0)x=±10，ory=-5 点，除去圆内部余1492个有效点，由Fortran程序 计算矩形内单位圆以外的点的流速并画出流 计算出流速后输出到文本文件，耗费时间以秒为 线图。 单位，再由Matlab程序读取画出流线图，耗费时间 在此算例中，边界剖分与计算流程与上述两 也是以秒为单位，与其他方法计算结果相对比精 例完全一致，计算、画图耗用时间以及结果精确度 度较高，画出的流线图与实际情况符合很 也完全一致。不同点在于，Matlab画图中放大后 好（图2）。 将清晰的出现2个小的漩涡。漩涡如图4所示。 42924.10 图2算例1流线图 算例2单位圆流速为(0,0)，矩形区域 [-5,5,-10,10]边界上的速度为 r(0,0) y=±5 25 (u,)= 20 (1-y2/25,0)x=±10 15 计算矩形内单位圆以外的点的流速并画出流 10 线图。 在此算例中，边界剖分与计算流程与算例1 完全一致，计算、画图耗用时间以及结果精确度也 完全一致（图3）。 15 图4算例3的漩涡图 0 2468 通过以上数值算例可以看出，该种方法精度 图3算例2流线图 较高，计算结果符合客观事实。 算例3单位圆流速为(0,0)，矩形区域 [-5,5,-10,10]边界上的速度为 (下转第126页) ?1994-2017 China Academic Journal Electronic Publishing House.All rights reserved.http://www.cnki.net外部闭曲线划分为 60 个单元，边界总剖分为 76 单元，由 Matlab 剖分在矩形区域内确定 1 521个 点，除去圆内部余 1 492个有效点，由 Fortran 程序 计算出流速后输出到文本文件，耗费时间以秒为 单位，再由 Matlab 程序读取画出流线图，耗费时间 也是以秒为单位，与其他方法计算结果相对比精 度较 高，画出的流线图与实际情况符合很 好( 图 2) 。 图 2 算例 1 流线图 算 例 2 单 位 圆 流 速 为 ( 0，0 ) ，矩 形 区 域 ［－ 5，5，－ 10，10］边界上的速度为 ( ) u，v = ( ) 0，0 y = ± 5 1 － y 2 ( ) /25，0 x = ± { 10 计算矩 形 内 单 位 圆 以 外 的 点 的 流 速 并 画 出 流 线图。 在此算例中，边界剖分与计算流程与算例 1 完全一致，计算、画图耗用时间以及结果精确度也 完全一致( 图 3) 。 图 3 算例 2 流线图 算 例 3 单 位 圆 流 速 为 ( 0，0 ) ，矩 形 区 域 ［－ 5，5，－ 10，10］边界上的速度为 ( ) u，v = ( ) 1，0 － 10 ≤ x ≤ 10，y = 5 ( ) 0，0 x = ± 10，or y = － { 5 计算矩 形 内 单 位 圆 以 外 的 点 的 流 速 并 画 出 流 线图。 在此算例中，边界剖分与计算流程与上述两 例完全一致，计算、画图耗用时间以及结果精确度 也完全一致。不同点在于，Matlab 画图中放大后 将清晰的出现 2 个小的漩涡。漩涡如图 4 所示。 图 4 算例 3 的漩涡图 通过以上数值算例可以看出，该种方法精度 较高，计算结果符合客观事实。 ( 下转第 126 页) 120 重庆理工大学学报 