5二十五届全国水动力学研讨会暨第十一届全国水动力学学术会议 图4阻力系数曲线(静止圆柱情形) 图5升力系数曲线(静止圆柱情形) 图6考查了不同积分区域下上述圆柱绕流问题升阻力系数计算的值,如图中所示对于不 同半径选取下的计算值,相对误差控制在1.5%以内,验证了实际数值计算中上述升阻力系数 计算公式与控制面Σ的选取无关,也反映了本问题数值计算结果的可靠性。 积分区域选取对CdC_1计算结果影响 静止圆柱绕流t=395.75 积分区域半径 图6不同积分区域升阻力系数计算结果(静止圆柱情形) 3总体流场特征 通过改变2中微分同胚的具体形式,可以获得不同形状的绕流问题的数值模拟,继而对 流场特征进行分析研究。 31水平放置椭圆柱尾迹 本节考察定常均匀来流椭圆柱绕流问题,该椭圆沿流向半轴为12,法向半轴为0.8,微 分同胚、控制方程、边界条件及其差分方法同前文分析,Re=100,计算域网格数为270×200, 时间步长取Mt=0.00125。图7为无量纲时间t405时刻的局部等流函数图,图8为对应的 等涡量云图,可以看到尾流中稳定发展的卡门涡阶,与圆柱相比,其尾迹较窄,涡阶从物体 表面脱落的形态较长。单个涡形态呈扁圆状第二十五届全国水动力学研讨会暨第十一届全国水动力学学术会议 - 4 - 图 4 阻力系数曲线（静止圆柱情形） 图 5 升力系数曲线（静止圆柱情形） 图 6 考查了不同积分区域下上述圆柱绕流问题升阻力系数计算的值，如图中所示对于不 同半径选取下的计算值，相对误差控制在 1.5%以内，验证了实际数值计算中上述升阻力系数 计算公式与控制面  的选取无关，也反映了本问题数值计算结果的可靠性。 图 6 不同积分区域升阻力系数计算结果（静止圆柱情形） 3 总体流场特征 通过改变 2 中微分同胚的具体形式，可以获得不同形状的绕流问题的数值模拟，继而对 流场特征进行分析研究。 3.1 水平放置椭圆柱尾迹 本节考察定常均匀来流椭圆柱绕流问题，该椭圆沿流向半轴为 1.2，法向半轴为 0.8，微 分同胚、控制方程、边界条件及其差分方法同前文分析，Re=100，计算域网格数为 270 200  ， 时间步长取  t 0.00125。图 7 为无量纲时间 t=405 时刻的局部等流函数图，图 8 为对应的 等涡量云图，可以看到尾流中稳定发展的卡门涡阶，与圆柱相比，其尾迹较窄，涡阶从物体 表面脱落的形态较长。单个涡形态呈扁圆状