|
楼主 |
发表于 2011-6-20 17:31:32
|
显示全部楼层
从开始的地方开始… 1
在我的第一篇博客中,我承诺将描述我在CFD行业的40年历程。如何才能是最好的开始?
红心国王给白兔什么样的忠告呢:“从开始的地方开始吧,一直到末尾,然后停止(译著:《爱丽丝漫游奇境记》中白兔问国王应该如何开始阅读)。” 尽管国王的回答并不是特指博客内容,但是它似乎是个很好的提示。所以,我将从开始的地方开始…
必须承认,我得以进入流体力学界纯属偶然。获得理学士学位后,我原本打算在帝国理工学院润滑学(目前比较通用的名称是摩擦学)获取博士学位。在最后时刻我判断该系没有录取我,但是当时我没有开始找工作,而是绝望地奔走在帝国理工学院寻找其他攻读硕士学位的机会。我记得1968年夏,我在帝国理工学院的机械工程学院围着课程安排的老师以及课题研究的督导们了解情况,他们中,就有Brian Launder,当时他是负责传热系硕士课程的初级讲师。我被他的热情所打动(可见他成功地向我推销了该学科!),于是注册了他的课程,并被承诺以后很有可能获得攻读博士学位的机会。
当时Spalding领导的传热系刚刚在CFD领域起步。Suhas Patankar具有开拓性的博士研究课题(与Spalding教授一同完成)在1965年结束,完成了二维边界层的求解方法, 现名为GENMIX, 这些年一直被广泛应用。大概在1968年,Akshai Runchal 和Mika Wolfstein 主导的博士课题,关于椭圆形循环流动的二维求解方法也已完成(参考文献:1969年, 伦敦,Gosman, Pun, Runchal, Spalding & Wolfstein共同发表的循环流动中的热量与质量传递)。我的CFD入门就是从使用这个软件开始的,为我的硕士项目“热虹吸管轴对称层流的流动和传热”。
之后,SRC出资组建了一个研究组,研究三维抛物型边界层流动的数值求解方法 —— 这是HTS(帝国理工学院传热系)第一次涉足三维流体力学领域 (也应该是所有人首次参与到这个领域)。这个项目开始于1969年。最初的团队成员包括Spalding教授,David Gosman (当时新上任的讲师,我博士课题导师),Larry Caretto(我认为,当时他从教于Berkeley Univeristy,正在休假),Bob Curr(HTS的一位研究人员),Devraj Sharma (PhD在读学生)和我。
与其它研究项目一样,我们从文献调查开始。然而大部分现成的先期工作却不是在HTS完成的,而是由在Los Alamos的Francis Harlow团队完成的,他们研发出了二维瞬态流动求解方法。我记得自己狠学了Harlow和 Welch发表的关于任意表面二维瞬态流动的论文。
与HTS的先期工作不一样,Harlow的求解方法直接求解的是压力以及速度等原始变量。相反,HTS的二维椭圆求解方法,求解的是流函数和旋涡状态(由此可推断压力和速度)等衍生变量。该方法的优点在于,避免了求解原始变量时压力和速度的耦合问题,缺点在于它很难扩展到三维流动。所以,处理我们的三维抛物型流动的研究工作,我们决定像Harlow的团队那样,直接求解压力和速度。
我忘记了研究组里面的具体工作——但是我非常幸运,我研究出了第一个软件,并将之应用到一个简单的三维抛物型流动中 ——一个矩形通道的入口区。马上我就遇到了后来历代CFD使用者都熟悉的问题。最初,自然而然地,我们使用了目前称之为同位网格的技术,即压力和速度在同一位置求解。使用它似乎是显然的,因为它是HTS已经完成的求解方法。有趣的是,尽管Harlow 和 Welch使用的是交错网格,我们在最初却没有想明白为什么他们使用这一技术!
我第一次求解时,马上遇到了经典问题“非耦合压力场”(有时候也被描述为“压力振荡”)。我们考虑到,中间差分方法可以求解动量方程中不可压流的压力,有限体积方程可以满足多个(二维情况下,4个)耦合压力场的求解。我推荐使用交错网格,速度在压力节点中间位置求解(与Harlow 和 Welch的技术相同)。自此,HTS的计算工作都使用该技术,几乎所有 CFD领域的工作都使用这一技术,直到上世纪80年代早期,Rhie 和 Chow提出了“非耦合压力场”问题的求解后,情况才有改变。
在1970左右的时间,Suhas Patenkar结束在IIT Kanpur的短暂工作,回到HTS,加入到我们的研究组。CFD历史将要书写开来!
待续…
更多更新,请登陆 http://blogs.mentor.com/mechanical-cn |
|