fortran理想流体的平面圆柱绕流程序.pdf

扬州大学 理想流体的平面圆柱绕流 计算水力学 题目： 用 Fortran 语言编写程序解决理想流体的平面圆柱绕流问题，如下图所示。由 于流动的对称性，可以只研究其中的四分之一区域，如图中 abcde 所示。 y e d c =1 2 r a b u =1 x x 2 3.5 3.5 在理想流体的平面运动中，流函数 ψ和势函数 Φ均满足拉氏方程： 2 2 2 2 2 2 0 ， 2 2 0 x y x y 其边界条件如下表所示。 流函数 势函数 在下边界 abc上 0 0 n 在出口边界 cd 上 0 const n 在上边界 de上 2 0 n 在进口边界 ea 上 y y ux 1 n 说明： 是切向流速， 是法向流速。 n n 下面就流函数进行讨论，为便于分析，把边界条件写成： ~ 在 1 上 其中： 1 为具有本质 B、C 的边界 0 在 2 上 2 为具有自然 B、C 的边界 n 水研 2010 级 第 1 页 共 15 页 水利学院 扬州大学 理想流体的平面圆柱绕流 计算水力学 解题步骤： (1)写出Галёркин积分表达式 2 2 dxdy 0