摘要
随着现代计算流体力学(CFD)的发展,对高速度、大容量计算机的需求越来越迫切,计算越来越复杂,规模也呈指数级增长。如此庞大的运算任务单台计算机显然无法胜任。一方面是单处理器性能发展空间有限,另一方面是专用高档并行计算机价格又太昂贵,而且旧有的资源又得不到持续和有效的利用,造成硬件资源空闲和资金浪费。因此,一种能够充分利用现有计算机软硬资源,有效发挥计算资源协同增强作用,同时满足不同专业研究人员高效方便处理的并行计算技术、平台和实现方法就变得十分重要和具有现实意义。
本文在综合国内外关于流体力学并行计算及相关领域研究文献的基础上,搭建了微机并行计算平台,解决了数值模拟对计算资源要求;对并行环境下子区域网格的独立生成和相邻子区域交界面的网格生成进行研究;在原有三维欧拉方程求解器的基础上添加粘性(耗散)项,针对定常问题,考虑粘性的影响,加入无滑移边界条件和两方程湍流模型,得到描述湍流运动的封闭RANS方程;对离心式叶片泵的设计理论进行了深入研究,并利用Visual C++平台编写了基于Pro/E的三维参数化设计和计算软件。最后利用以上得出的计算经验和结论,通过对一个实际工程问题—离心式叶片泵内部流动的数值计算,完成了CFD中湍流高级数值模拟技术的尝试。全文共七章。
第一章为绪论,综述了国内外有关并行计算技术的研究概况,阐述了本文的研究目的、意义与主要研究内容。
第二章对Windows2000操作系统下的并行计算环境进行了深入研究,在此基础上成功搭建了由4台PC组成的局域网并行集群系统,采用消息传递接口MPI作为并行程序解决方案,安装了MPICH并行程序运行环境,用C语言调用MPI消息传递库,并在此基础上编写和实现了流场计算中应用到的并行迭代程序、矩阵二维拓扑和三维拓扑并行程序。
第三章研究了并行环境下子区域网格的独立生成和相邻子区域交界面的网格生成。在网格并行生成过程中,给出了关于边的分类,完善了文献中子区域内网格生成时接受新点及新单元的条件,给出了流场整体网格图以及流场子区域、交界面局部网格图。
第四章对κ-ε两方程模型进行了深入研究,在原有三维欧拉方程求解器的基础上添加粘性(耗散)项,针对定常问题,考虑粘性的影响,加入无滑移边界条件和两方程湍流模型,得到描述湍流运动的封闭RANS方程。
第五章利用Fluent软件在搭建的并行计算平台上对离心式叶片泵的内部流场进行数值模拟,对湍流κ-ε两方程模型进行验证,并得到了一些结论。
第六章通过对离心式叶片泵的优秀水力模型进行了深入的研究,在一元和三元流动理论的基础上给出了叶轮设计中的各种参数计算公式;利用Visual C++平台
With the development of modern computer fluid dynamics, the need for high performance and high capacity becomes more and more urgent., computation becomes more and more complex, scale develops with positive exponent. For such tremendous calculation tasks, a single computer is not able to complete these computations obviously: one is the limited capability of single processor; the other is those special expensive parallel computers are too expensive and formerly resource have to be abandoned, which is waste. Therefore, it is necessary and realistic meaning for us to develop an effective parallel calculating technology, platform and realizing method for different researchers, at the same time we can use present computer software and hardware resources fully. In the basis of the references on a centrifugal pump of two-phase fluids from domestic and oversea, its theory of boundary layer is presented and set up and the modular curve equation of vane is received in this paper. It divides into 7 chapters.The introduction is in chapter one, which summarizes the research results of parallel computation technology in domestic and oversea and illuminates the research aim, significant and main contents.In chapter two, through deep research to windows 2000 operation system, wehave construct a parallel PC cluster-based system with 4 computers; also can extending more computer nodes easily. We use message passing interface (MPI) as a method to program parallel software, then MPICH has been installed and use C langue to call MPI function warehouse. Finally, we have made and realized parallel iterative arithmetic, matrix of two-dimensional and topological arithmetic, three-dimensional topological arithmetic, molecule dynamics simulation and three-dimensional graphics processing parallel programs.In chapter three, we researches the generation of sub domain mesh and the adjacent domain mesh under parallel environment, firstly improve the R lohner's wave front domain-splitting algorithm so that the sub grids and their boundaries are more favorable for grid generation. A new optimization strategy of sub domain's boundary is then presented in order to improve the smoothness of boundaries and the quality of grids. After subdividing the initial meshes by using the above domain-splitting algorithm, we smooth the grid by the parallel method.In chapter four, the NS solver is developed from Euler equation code by adding viscous terms and two-equation turbulence models. In order to get a better result, the improved Jameson second and forth artificial dissipation is used, the implicit residual
引文
1 胡凯,宋京民编著,网络计算新技术,北京:科学出版社,2001.7.pp189-322
2 刘信安,多线程计算在环境与计算化学中的应用,计算机应用化学,2000.17(3).pp 273—279
3 刘信安等,C/S分布式并行处理技术在高速局域网平台上实现非线性环境化学胞映射体系的大规模计算,计算机与应用化学,2002.19(1).pp 163-168
4 陈国良编著,并行计算—结构·算法·编程,北京:高等教育出版社,1999.10
5 金怡镰,高性能计算机的关键技术和发展趋势,中国工程科学,2001.3(6).pp 1-6
6 Gorden Bell等著,马云峰译:高性能计算的现状及发展趋势,高性能计算技术,2003.pp61-64
7 胡峰等,并行计算技术与并行算法综述.电脑与信息技术,1999.5.pp47-59
8 American National Standards Institute. ANSI Technical Committee X3H5. Parallel Processing Model for High-Level Programming Languages, 1993.
9 IEEE. POSIX P1003. 4a: Threads Extension for Portable Operating Systems. Piscataway, NJ: IEEE Press. 1994.
10 Silicon Graphics, IRIS Power C User's Guide, Silicon Graphics Computer Systems, Mountain View, CA, 1989.
11 Tucker, L. W., Mainwaring, a. CMMD: active messages on the CM—5. Parallel Computing, 1994. 20(4). pp481-496
12 Kolawa, A. Parasoft: a comprehensive approach to parallel and distributed computing. In: IEEE Computer Society, ed. Proceedings of the Workshop on Cluster Computing. Los Alamitos, CA: IEEE Press, ]992.
13 Anderson J D,计算流体力学入门[M],北京:清华大学出版社,2002
14 陈国良,并行算法的设计与分析,北京:高等教育出版社,1994.5
15 都志辉编著,高性能计算并行编程技术—MPI并行程序设计,北京:清华大学出版社,2001.8
16 黎康保,黎文楼.用PC机群组构并行超级计算机,计算机工程,2000.26(9).pp 1-3
17 罗海彬等,PC Cluster建立和技术探讨.计算机与应用化学,2002.19(1).pp 116-121
18 孙家褪,张林波等,网络并行计算与分布式编程环境,北京:科学出版社,1997
19 Gropp, William. A high-performance, portable implementation of the MPI message passing interface standard. Parallel Computing, 1996, 22(6). pp789-828
20 Dantas, M. A. R. Efficient scheduling of MPI applications on networks of workstations. Future Generation Computer Systems, 1998, 13(6). pp489-499
21 Owczarz, Wojciech; Zlatev, Zahari. Parallel matrix computations in air pollution modeling, Parallel Computing, 2002, 28(2). pp355—368
22 Wolf, Felix; Mohr, Bernd. Automatic performance analysis of hybrid MPI/OpenMP applications. Journal of Systems Architecture, 2003, 49(10—11). pp421-439
23 Lauria, Mario; Chien, Andrew. MPI-FM: High Performance MPI on Workstation Clusters. Journal of Parallel and Distributed Computing, 1997, 40(1). pp4-18
24 Yu, C. S.; Berlamont, J.; Embrechts, H. Modelling Coastal Sediment Transport on a Parallel Computer. Physics and Chemistry of the Earth, 1998, 23(5-6). pp497-504
25 张兆顺著,《湍流》,国防工业出版社,2002.1
26 朱自强等编著,《应用计算流体力学》,北京航空航天大学出版社,1998.8
27 是勋刚,《湍流》,天津:天津大学出版社,1994
28 傅德熏,《计算空气动力学》,宇航出版社,1994
29 Antti Hellsten, "New Two-Equation Turbulence Model For Aerodynamics Applications", Report A-21
30 Sherrie L. Krist, Robert T. Biedron and Christopher L. Rumsey "CFL3D User's Manual", NASA/TM-1998-208444
31 Abid, R.," Evaluation of Two-Equation Turbulence Models for Predicting Transitional Flows", International Journal of Engineering Science, Vol. 31, pp. 831-840, 1993
32 Wilcox, D., "Turbulence Modeling for CFD", DCW industries, inc., La Canada, California, 1993
33 F. R. Menter, "Two-equation eddy-viscosity turbulence models for engineering applications". AIAA Journal, 32(8): 1598-1605, 1994.—
34 W. Rodi and N. N. Mansour. "Low Reynolds number k—6 modeling with the aid of direct simulation data". Journal of Fluid Mechanics, 250: 509-529, 1993
35 D. C. Wilcox. Reassessment of the scale-determining equation for advanced turbulence models. AIAA Journal, 26(11): 1299-1310, 1988.
36 Girimaji, S., "Fully-Explicit and Self-Consistent Algebraic Reynolds Stress Model," NASA CR-198243, December 1995; also ICASE Report No. 95-82, 1995 and Theoretical and Computational Fluid Dynamics, Vol. 8, pp. 387-402, 1996.
37 Abid, R., Rumsey, C., and Gatski, T., "Prediction of no-equilibrium Turbulent Flows with Explicit Algebraic Stress Models," AMA Journal, Vol. 33, No..11, pp. 2026-2031, November 1995.
38 R. Radespiel and C. Rossow, "A Cell vertex Finite Volume Scheme for the Two-Dimensional Navier-Stokes Equations", IB 129-87140, Institute Entwurfs aerodynamik Braunschweig
39 Abdoi-Hamid, K., Lakshmanan, B., and Carlson, J., "Application of Navier-Stokes Code PAB3D with κ-ε Turbulence Model to Attached and Separated Flows," NASA TP. 3480, 1995
40 何希杰,离心泵理论及设计方法的发展,华北水泵,1993(2):4~7
41 李金生等,国外杂质泵发展概况,全国杂质泵技术交流会论文专辑,1992(7):143~150
42 陈次昌等编著,两相流泵的理论与设计,北京:兵器:工业出版社,1994
43 许洪元,罗先武 主编,磨料固液泵,北京:清华大学出版社,2000.3
44 A.T.特罗斯科兰斯基 & S.拉扎尔基维茨[波兰]编著,耿惠彬译,霍励强校.叶片泵计算与结构,北京:机械工业出版社,1981.11
45 赵振海,离心泵流道中固体颗粒运动规律的研究,硕士学位论文,北京清华大学,1989
46 戴江,离心泵叶轮内固液两相紊流流动规律研究,博士学位论文,北京清华大学,1994
47 姜培正等,离心泵叶轮内液固两相流场的实验研究与理论分析,全国第二届杂质泵及固体物管道水力输送学术讨论会论文集,四川.自贡.1999,7:23~27
48 吴玉林等,杂质泵叶轮中固体颗粒运动规律的实验,清华大学学报(自然科学版),Vol.32.1992(5):52~59
49 吴玉林等,渣浆泵叶轮中固液两相湍流的计算和实验,清华大学学报(自然科学版),Vol.38.1998(1):71~74
50 陶文铨编著,数值传热学[M],西安:西安交通大学出版社,1988
51 李文广著,水力机械流动理论[M],北京:中国石化出版社,2000
52 赵兴艳等,CFD方法在流体机械设计中的应用[[J],流体机械,2000 28(3):22^-25
53 吴玉林等,水轮机转轮内部的三维固液两相紊流计算[[J],工程热物理学报,1997,18(5):580-583
54 曹树良等,三峡机组蜗壳内紊流的三维数值分析[[J],清华大学学报(自然科学版),1997(8):95-98
55 曹树良等,水轮机蜗壳和固定导叶内部三维紊流流场数值研究[J],工程热物理学报,1998(5):586-589
56 曹树良等,用代数应力紊流模型预估水轮机转轮内部三维流场[J],清华大学学报,1998,38(4),113-116
57 周晓泉等,水轮机蜗壳和固定导叶内部流动的数值模拟[[J],清华大学学报(自然科学版),2000(8):93-97
58 陈池,无过载离心泵叶轮内三维不可压湍流场计算,江苏理工大学硕士论文,1999
59 陈洪海,离心泵蜗壳内部三维湍流场数值模拟,江苏大学硕士论文,2002
60 李映,无堵塞泵内部三维固液湍流场数值模拟,江苏理工大学硕士论文,2000.
61 刘厚林,流道式无堵塞泵CAD软件开发及三维不可压湍流计算[D],博士论文,2001.4.
62 朱自强等编著,应用计算流体力学[M],北京:北京航空航天大学出版社,1999
63 陈景仁(美)著,湍流模型及有限分析法[M],上海:上海交通大学出版社,1989
64 M. P. O Brien & R. G. Folsom. The Transportation of Sand in Pipelines. Univ. of Calif. Pub. in Eng. 1937(7): 48
65 L. C. Fairbank. Effect on the characteristics of centrifugal pump. Paper No. 2167, 1942.
66 太田启吉良,理论的考察,机械学会文集,1962,4(66):91~100.
67 W. B. Gregory Pumping Clay slurry through a Four-inch Pipe. Mech. Eng. 1927, (49): 606~616
68 三云,西原,高原.微粒矿石水力输送,機誌,昭和8年12月
69 掘田正雄,性能機論,昭和14年,5(21):1~5
70 佐佐木和郎等,石炭流速时性能变化,礦业會誌,昭和 33年 4(839):287
71 陶文铨著,计算传热学的近代进展[M],北京:科学出版社,2000
72 孙家广等著,计算机图形学(第三版)[M],北京:清华大学出版社,1999
73 王保国等著,叶轮机械跨声速及亚声速流场的计算方法[M],北京:国防出版社,2000
74 金树德等编,现代水泵设计方法[M],北京:兵器工业出版社.1993
75 任静等,基于CFD流态分析技术的水力组转轮改型研究[[J],机械工程学报,2000 36(4),86-90
76 苏铭德,黄素逸,计算流体力学基础[M],北京:清华大学出版社,1997
77 叶正寅,杨永年,非结构网格生成技术中的一种直接提供背景信息的方法,西北工业大学学报,1999,17(1),P1-50
78 叶正寅,杨永年,钟诚文,张仲寅,陈迎春,二维非结构网格生成及Euler方程计算的方法研究,计算物理,第16卷,第6期,1999
79 桑为民,非结构网格多段翼型绕流Eule:及N-S方程的数值模拟,西北工业大学硕士学位论文,1999
80 柳长安,非结构网格生成及二维非定常流场及悬挂物投放数值模拟,西北工业大学硕士学位论文,2000
81 梁强,非定常气动力的N-S方程求解及应用,西北工业大学硕士学位论文,2001
82 王刚,三维非结构粘性网格生成与N-S方程求解,西北工业大学硕士学位论文,2001.
83 朱自强,应用计算流体力学,北京航空航天大学出版社,2001
84 吴子牛,计算流体力学原理,科学出版社,2001
85 张涵信,沈孟育,计算流体力学——差分方法的原理和应用,国防工业出版社,2003
86 M. Dryja and O. B. Widlund, Some domain decomposition algorithms for elliptic problems. Pages 273-291, SIAM, 1992.
87 H. L. de Cougny, M. S. Shephard and C. Ozturan, 'Parallel three-dimensional mesh generaion on distributed memory MIMD Computers: Tech. Rep. Scorec Rep. #7, Rensselaer Polytechnic Institute, 1995.
88 R. Lohner, An adaptive finite element scheme for transient problems in CFD; Comp. Mesh. Appl. Mech. Eng. 61, 323-338, 1987.
89 J. M. Swisshelm and G. M. Johnson, Parallel computation of Euler and N-S flows, Appl Meth. Comp., 321-332, 19, 1986.
90 S. W. Hammond and T. J. Barth, Efficient massively parallel Euler solver for two-dimensional unstructured grids, AIAA J., 947-952, 30, 1992.
91 R. Das, D. J. Mavriplis, S. Gupta and R. Ponusam, The design and implementation of a parallel unstructured Euler solver using software primitives, AIAA-92-0562.