快速傅立叶转换和有限元素法求解黏性不可压缩流经过圆柱体之比较.pdfVIP

快速傅立叶转换与有限元素法求解黏性不可压缩 流经过圆柱体之比较 方俊、夏育群、郑仙志 逢甲大学航空工程学系，台中，40724 廖祟钧 逢甲大学机械r-程研究所．台中．40724 【摘要】本研究以两种混合的数值方法求解黏性不可压缩流经圆柱体之瞬 时流场并比较其优劣。针对此类的流场．其统御方程可采用二维圆柱坐标不 可压缩之Navior-Stokes方程式，并简化推导出之压力．速度方程式或涡度．流 线函数方程式，科应用有限差分法或有限元素法等方式求解，近年来已有相 当完整的研究。但是不论有限差分法的迭代技巧及有限元素法的庞人联盘方 程式，均使求解面对严重的效率问题。本研究尝试提出两种混合数值模拟方 式，来求解涡度·流线函数方程式。首先，涡度方程式仍以有限差分法求解。 其次，流线方程式部份则分别以快速傅立叶转换或有限元素法来求解。最后 井将两者之结果比较以了解其优劣。 f关键词】快速傅截叶转换．有限元素法．不可雎缩流，圆柱体，计算流体力 学 1 引言 流体力学方面的研究通常可归纳为理论推导、实验量测观察与数值模拟 等方法。一般来说，藉由实验研究依现场观察及量测所取得归纳的数据最为 可靠、准确。但由于流体现象的复杂，实验常需消耗许多人力、时间及金钱， 甚至受到环境因素的限制，极不利于效率性与经济性的考量。因此在计算机 性能突飞猛进之下，数值模拟的研究开发出非常多元而精确的方法，如有限 差分法(FiniteDifference，FD)、有限体积法(FiniteVolume，FV)及有限元 素法(FiniteElement,FE)等。 图一、非稳态不可压缩流流经圆柱体之实验流场Re=3000，T=5 148 非稳态不可压缩流经过一圆柱体时之流场是数值仿真非常有兴趣的主 题之一，其流场特征如图一所示[121。此类流场之统御方程可采用二维圆柱 坐标不可压缩之Navior-Stokes方程式，并将其推导成压力一速度方程式 function)[6】， (pressure．velocity)[4】或涡度-流线函数方程式(Vorticity-stream Fluid 应用传统之计算流体力学(ComputationalDynamics，CFD)方式，如 有限差分法或有限元素法求解此类方程式，文献均有相当完整的研究【3，4，6， 7．81。不过在求解流线函数方程式时．有限差分法必须运用迭代技巧，而有 限元素法亦须面对庞大联立方程式，其CPU计算时间相当不经济。因此本 研究尝试采用两种混合数值模拟法来求解涡度．流线函数方程式，其中涡度方 Direction 程式仍以有限差分之ADI(AlternatingImplicit)法【l】求解，流线方 程式部份则分别以(一)快速傅立叶转换【6】及(二)有限元素法[2]来求解， 如此可大幅改善其效率问题。另外，本研究并比较两者结果之优劣。 2数值模拟 首先将二维圆柱坐标不可压缩之Navior-Stokes方程式推导成涡度一流线 函数方程式如F： ㈩ 等Ⅵ等+了uo丽a(=去B昙(r等)+上r2丝002] 搿，割+吉等叫 ㈤ l a渺 a∥ 蚱27亩、‰一言 其中f为涡度，矿为流线函数，圆柱半径a，自由流速度为圪，雷诺数 为Re=p比(2口)／∥。无因次化时，长度除以a、速度除以圪、时间除以 口／圪。其边界条件如下： ，=1 y_0，警_o r_∞ f=0，∥=rsin0