计算流体力学的发展及应用_魏淑贤

第27卷第2期河北理工学院学报V o l.27N o.2 2005年5月J o u r n a l o f H e b e i I n s t i t u t e o f T e c h n o l o g y M a y.2005

文章编号:100 7-28 29(20 05)02 -0 115 -03

计算流体力学的发展及应用

魏淑贤, 沈跃,黄延军

(石油大学(华东)物理科学与技术学院, 山东东营 2 570 61)

关键词:计算流体力学;发展;应用

摘要:计算流体力学是流体力学的一个分支。它用于求解固定几何形状空间内的流体的动

量、热量和质量方程以及相关的其它方程,并通过计算机模拟获得某种流体在特定条件下的有

关信息,是分析和解决问题的强有力和用途广泛的工具。对计算流体力学的发展和应用进行

了综述并对其发展趋势做了探讨。

中图分类号:O368文献标识码:A

1 计算流体力学的发展

20世纪30年代,由于飞机工业的需要,要求用流体力学理论来了解和指导飞机设计, 当时,由于飞行速

度很低,可以忽略粘性和旋涡,因此流动的模型为L a p l a c e方程,研究工作的重点是椭圆型方程的数值解[ 1]。利用复变函数理论和解的迭加方法来求解析解。随着飞机外形设计越来越复杂,出现了求解奇异边界积分

方程的方法。以后,为了考虑粘性效应,有了边界层方程的数值计算方法,并发展成以位势方程为外流方程, 与内流边界层方程相结合,通过迭代求解粘性干扰流场的计算方法。

同一时期,许多数学家研究了偏微分方程的数学理论, H a d a m a r d,C o u r a n t,F r i e d r i c h s等人研究了偏微分方程的基本特性、数学提法的适定性、物理波的传播特性等问题,发展了双曲型偏微分方程理论。以后, C o u- r a n t,F r i e d r i c h s,L e w y等人发表了经典论文[ 2],证明了连续的椭圆型、抛物型和双曲型方程组解的存在性和唯一性定理,且针对线性方程的初值问题,首先将偏微分方程离散化,然后证明了离散系统收敛到连续系统,最后利用代数方法确定了差分解的存在性;他们还给出了著名的稳定性判别条件:C F L条件。这些工作是差分方法的数学理论基础。20世纪40年代, V o n N e u m a n

n,R i c h t m y e r,H o p f,L a x和其他一些学者建立了非线性双曲型方程守恒定律的数值方法理论, 为含有激波的气体流动数值模拟打下了理论基础。

在2 0世纪50年代,仅采用当时流体力学的方法,研究较复杂的非线性流动现象是不够的, 特别是不能满足高速发展起来的宇航飞行器绕流流场特性研究的需要。针对这种情况, 一些学者开始将基于双曲型方程数学理论基础的时间相关方法用于求解宇航飞行器的气体定常绕流流场问题,这种方法虽然要求花费更多的计算机时,但因数学提法适定,又有较好的理论基础,且能模拟流体运动的非定常过程, 所以在60年代

[3]

这是应用范围较广的一般方法。以后由L a x、K r e i s

s和其他著者给出的非定常偏微分方程差分逼近的稳定性理论,进一步促进了时间相关方法。当时还出现了一些针对具体问题发展起来的特殊算法。

值得一提的是,我国在2 0世纪50年代也开始了计算流体力学方面的研究[3]。我国早期的工作是研究钝头体超声速无粘绕流流场的数值解方法,研究钝头体绕流数值解的反方法和正方法。以后,随着我国宇航事业的发展,超声速、高超声速绕流数值计算方法的研究工作发展很快。对定常欧拉方程数值解的计算方法进行研究,并给出了钝体超声速三维无粘绕流流场的计算结果。

20世纪70年代,在计算流体力学中取得较大成功的是飞行器跨音速绕流数值计算方法的研究。首先

收稿日期:200 4-11-02

作者简介:魏淑贤(19 78-),女, 山东临沂人, 石油大学(华东)物理科学与技术学院硕士研究生。

制 11 6 河 北 理 工 学 院 学 报 第 27卷 是 M u r m a n 和 C o l e 用松弛方法求解位势流小扰动方程 , 数值模拟带激波的跨声速绕流场 , 解决了跨声速绕流

中的混合问题 。 在他们的工作中第一次将迎风格式应用于空气动力学问题的模拟 。 不久以后 J a m e s o n 提出

了旋转格式 , 将穆尔曼 -科勒方法推广于求解三维跨声速绕流的全位势流方程 , 获得成功 [ 3] 。 同一时期 , 我

国开展了采用时间相关方法求解非定常欧拉方程 、可压缩 N -S 方程和简化 N -S 方程的计算方法研究 。

在 差分格式的构造方面 , 提出了求解欧拉方程的特征符号分裂法和三层格式[ 4] 等 。 在可压缩 N -S 方程

的求 解中 , 计算方法有了很大进展 , 先后提出了开关函数法 、调解因子方法 、紧致迎风格式 、推进迭代法 、无波

动无 自由参数的耗散格式 、界值为限格式和耗散比拟方法等[ 3] 。 这些研究工作进一步改进了计算方法精

度 , 提 高了求解效率 , 且对流场激波的数值模拟有较高的分辨能力 。 而且这些研究成果使得我们在计算流体力学 的差分方法研究工作中初步形成了自己的特点 。

进入 20 世纪 8 0年代以后 , 计算机硬件技术有了突飞猛进的发展 , 千万次机 、亿次机逐渐进入人们的实 践活动范围 。随着计算方法的不断改进和数值分析理论的发展高精度数值模拟已不再是天方夜谭 。 同时随 着人类生产实践活动的不断发展 , 科学技术的日新月异 , 一大批高新技术产业对计算流体力学提出了新的要 求 , 同时也为计算流体力学的发展提供了新的机遇 。实践与理论的不断互动 , 形成计算流体力学的新热点 、 新动力 , 从而推动计算流体力学不断向前发展 。

首先 , 在计算模型方面 , 又提出了一些新的模型 , 如新的大涡模拟模型 、考虑壁面曲率等效应的新的湍流

模式 、新的多相流模式 、新的飞行器气动分析与热结构的一体化模型等 [ 5] 。这就使得计算流体力学的计算

模型由最初的 E u l e r 和 N -S 方程 , 扩展到包括湍流 、两相流 、化学非平衡 、太阳风等问题研究模型在内的多

个模型 [ 6] 。其中以考虑更多流动机制 , 如各向异性的非线性 (应力 /应变关系 )湍流研究为重点 。研究结果 再次证明 , 万能的湍流模型还不存在 , 重要的是如何在模型精度和计算量上较好地取得折中 ;也有学者从更 高层次研究湍流模型问题 , 由湍流流动中速度不可微 , 怀疑 N -S 方程的有效性 , 进而提出以积分方程为基

础的数学模型[ 6] 。

其次 , 在计算方法方面 , 又提出了一些新的计算方法 , 如新的遗传算法 、无网格算法 、新型高精度紧致格

式 、气动计算的新变分原理 、结构 /非结构混合网格新技术 、新型动网格技术等等[ 5] 。 目前计算方法的研究

集中在高精度格式方法 , 即追求三阶精度以上 , 其中又以解决真正实际问题 。 除此之外 , 计算方法研究还涉 及带限制器的高阶插值 、谱方法 、拉格朗日方法 , 时 -空守恒元方法等等 。 将其它方法引进传统的计算流体

力学也是现阶段的重要成果之一[ 6] , 其中特别值得一提的是将基因算法与传统计算流体力学结合在一起 ,

在域分裂和最优化设计等许多方面显示出了良好的应用前景 。在算法分析上 , 除传统的精度 、稳定性 、收敛 性等方面的分析 , 还有更深层次的数值动力学分析 , 即将数值方法看成是动力系统来进行分析 , 揭示了许多 奇异的数值现象 。

再次 , 在研究成果方面 , 英国 M .A .L e s d h z i a e r 关于湍流模型 、美国 H .C .Y e e 关于计算不确定性 、

日本学 者的玻耳兹曼方程解流动问题 、德国的 E .v o n L a v a n t e 关于使用并行计算机进行发动机气缸流场涡和激波的 非定常流动模拟等等 [ 6]

, 都有较新的学术思想 , 较高的学术水平 。

目前 , 计算流体力学研究的热点是 :研究计算方法 , 包括并行算法和各种新型算法 ;研究涡运动和湍流 , 包括可压和不可压湍流的直接数值模拟 、大涡模拟和湍流机理 ;研究网格生成技术及计算机优化设计 ;研究 计算流体力学用于解决实际流动问题 , 包括计算生物力学 、计算声学 、微型机械流动 、多相流及涡轮机械流动

的数值模拟等[ 7] 。 2 计算流体力学的应用

计算流体力学的应用已经从最初的航空航天领域不断地扩展到船舶

[ 8] 、海洋 [ 9] 、化学

[ 10] 、铸造 [ 11] 、 冷 [ 12] 、工业设计 、城市规划设计[ 13] 、建筑消防设计 [ 14] 、汽车 [ 15] 等多个领域 。 近几年来计算流体力学在全机

流场计算 、旋翼计算 、航空发动机内流计算 、导弹投放 、飞机外挂物 、水下流体力学 、汽车等 [ 6] 方面获得广泛

应用 。 这表明计算流体力学在解决工程实际问题方面具有重要的应用价值 。下面仅以在汽车领域的应用为 例 , 介绍计算流体力学应用于工程实际中的速度和深度 。 20 世纪 80 年代初期才开始有计算流体力学应用 于汽车领域的论文发表 , 经过短短的二十余年 , 其应用已涉及到汽车车身设计 、汽车内部空间的空调与通风 、