空调房间内的三维非均温气流组织的数值研究
|
摘要
本论文应用计算流体力学的理论和方法模拟三维较复杂条件下空调房间室内空气流场的湍流流动。研究在强制对流的作用下，室内气流的流动结构、速度、温度分布以及换热规律。在对数学模型及物理模型进行理论分析的基础上，对含湍流双方程模型及非均温自然对流的耦合方程组采用了联立迭代解方法，并使用了对流-扩散方程系数的计算法、SIMPLEC算法、交错网格法、ADI线迭代法、块修正技术等，边界条件采用了的k-ε模型结合壁面函数的方法。并将数值模拟的结果采用计算机图形学的直观方式来表示多种气流组织方案的气流流型，分析讨论其气流分布规律、特点，及其适用场合。计算的结果证实了利用k-ε湍流模型对非均温的房间的气流组织进行数值计算是可行的。这一数值模拟方法也能帮助设计人员预测其气流设计方案的效果。这对于保证良好的房间空调系统气流组织设计方案、提高室内空气品质、减少建筑物的能耗，对建筑物热环境的节能设计及室内空气质量的评价都有重要的意义。
This paper is an application of Computational Fluid Dynamics in the turbulent flow of three-dimensional air-conditioned room. The airflow structure, velocity, temperature and the law of heat exchange on the condition of the forced convection are studied. On the basis of the theoretical analysis of physical and mathematical model, the iterative method is adopted to the conjugated governing equations, and the numeration of the coefficient to the convection-diffusion equation, SIMPLEC algorism, staggered grid method, ADI line iteration method, and block correction technique are used, and so on. K-ε turbulence model linked with the wall function are applied to the boundary condition. The results of numerical simulation for many different kinds of airflow are expressed in streamline chart of airflow structure through the computer graphics, and the law and character of the airflow distribution and the application are discussed. The results validate that the K-ε turbulence model in the numerical computation to the airflow structure of thermal non-uniform room is feasible. At the same time, it is elicited that the numerical simulation also helps man to predict the effect of the airflow designing. Its application is of great importance to the increase of the level of the designing project of airflow organization in the air-conditioned room, the increase of the air quality in the room, the decrease of the energy expensed in the building, the energy saving design to the heat environment of the building and the evaluation of the air quality in it.
This paper is an application of Computational Fluid Dynamics in the turbulent flow of three-dimensional air-conditioned room. The airflow structure, velocity, temperature and the law of heat exchange on the condition of the forced convection are studied. On the basis of the theoretical analysis of physical and mathematical model, the iterative method is adopted to the conjugated governing equations, and the numeration of the coefficient to the convection-diffusion equation, SIMPLEC algorism, staggered grid method, ADI line iteration method, and block correction technique are used, and so on. K-ε turbulence model linked with the wall function are applied to the boundary condition. The results of numerical simulation for many different kinds of airflow are expressed in streamline chart of airflow structure through the computer graphics, and the law and character of the airflow distribution and the application are discussed. The results validate that the K-ε turbulence model in the numerical computation to the airflow structure of thermal non-uniform room is feasible. At the same time, it is elicited that the numerical simulation also helps man to predict the effect of the airflow designing. Its application is of great importance to the increase of the level of the designing project of airflow organization in the air-conditioned room, the increase of the air quality in the room, the decrease of the energy expensed in the building, the energy saving design to the heat environment of the building and the evaluation of the air quality in it.
引文
1．戴有为．对流传热分析．第一版．北京：兵器工业出版社．1989
2．陶文铨．数值传热学．第一版．西安：西安交通大学出版社．1988
3．陶文铨．计算流体力学与传热学．北京：中国建筑工业出版社．1992
4．帕坦卡[美]．传热与流体流动的数值计算．第一版．北京：科学出版社．1984
5．赵学端．粘性流体力学．第一版．北京：机械工业出版社．1983
6．潘文全．流体力学基础(上册)．第一版．北京：机械工业出版社．1982
7．潘文全．流体力学基础(下册)．第一版．北京：机械工业出版社．1982
8．尤国钊．计算流体力学．第一版．南京：南京理工大学．1994
9．马铁犹．计算流体力学．第一版．北京：航空学院出版社．1986．6
10．第六届全国计算流体力学会议论文集 科学出版社．1992
11．陈景仁．湍流模型及有限分析法．第一版．上海：上海交通大学出版社．1992
12．郭宽良．计算传热学．第一版．合肥：中国科学技术大学出版社．1988
13．章熙民．任泽霈．梅飞鸣．王中铮．传热学．第二版．北京：中国建筑工业出版社．1992
14．赵荣义．范存养．薛殿华．钱以明．空气调节．第三版 北京：中国建筑工业出版社．1994
15．周亨达．工程流体力学．第二版．北京．北京钢铁学院 1988
16．孙一坚．工业通风．第三版．北京．中国建筑工业出版社．1994
17．陆金甫．偏微分方程的数值解法．第一版．北京：清华大学出版社 1987
18．同济大学数学教研室．高等数学(上册)．第三版．上海：高等教育出版社．1995
19．同济大学数学教研室．高等数学(下册)．第三版．上海：高等教育出版社．1995
20．金练，欧阳嘀，张洁等．暖卫 通风 空调技术手册．第一版．北京：中国建筑工业出版社．2000．7
21．陆耀庆．暖通空调设计指南．第一版．北京：中国建筑工业出版社．1987
22．孙一坚．工业通风．第三版．北京：中国建筑工业出版社出版，1994
23．赵荣义，范存养，薛殿华，钱以明．空气调节．第三版．中国建筑工业出版社，1994
24．吴元炜．暖通空调新技术．第一版．北京：中国建筑工业出版社，1999
25．黄翔，连之伟，哈文．空调工程应用．第一版．北京：科学出版社，1999
26．王启杰．对流传热传质分析西安．第一版．西安：西安交通大学出版社．1991
27．陆耀庆．HAvC暖通空调设计指南．第一版．北京：中国建筑工业出版社，1996
28．刘旭，冯玉琪．实用空调技术精华．第一版．北京：人民邮电出版社，2001
29．潘云钢．高层民用建筑空调设计，第一版．北京：中国建筑工业出版社，1990
30．何耀东，何青．中央空调，第一版．北京：冶金工业出版社，1998
31．龚光彩．汤广发．吕文瑚．高大空间横向隔断气流的计算机仿真．湖南大学报．1998．25(4)：71-76
32．黄维平．高大门厅的空调送回风口设计探讨．制冷．1997．3：51-53
33．龚光彩．汤广发．高大空间纵向隔断气流的数值模拟．湖南大学学报(社科版)．1999．26(1)：66-71
34．何冰．空调房间气流组织的三维湍流计算．西安交通大学硕士论文
35．周绍荣，程运林，汤莉．高大空间气流组织的研究．湖南大学学报(自然科学版)．1998．25(2)：70-74
36．范存养．办公室下送风空调方式的应用．暖通空调．1997．27(4)：30-39
37．李湘洲，赵志江．日本的暖通空调及建筑节能现状．节能技术．1995．5：15-18
38．陈爱玲．空调房间气流组织的初步探讨．流体机械．1996．24(9)：51-55
39．王汉青，汤广发．三维紊流室内气流数值解．空气动力学报．1988．4(6)
40．魏学孟，樊洪明，张维功．失流洁净室数值模拟与实验研究．暖通空调．1997．27．(2)：11-14
41．李俊．空调的节能设计．住宅科技．1996．7：15-17
42．谭洪卫．计算流体动力学在建筑环境工程上的应用．暖通空调．1999，29(4)：31-36
43．龚光彩．CFD技术在暖通空调制冷工程中的应用．暖通空调．1999，29(6)：25-27
44．黄晨，李美玲，左涛等．大空间建筑室内垂直温度分布的计算方法．制冷学报．1999．3：25-30
45．陈在康．暖通CAD技术的开发和应用．暖通空调．1993．6
46．郭春信．室外采风口和排风口与系统的节能．1993．6
47．陈宾．日本大空间空调现状．暖通空调．1994．3
48．高养田．空调新风系统设计若干问题分析．暖通空调．1993．2：70-74
49．熊螯魁．不可压流场数值模拟方法综述．武汉交通科技大学学报．1999．10：523-528
50．张维，汤广发．风对建筑绕流流动的数值方法．空气动力学学报．1990．12：437-443
51．杨秉俭，蔡临宁．用代数应力模型模化型腔充填过程中的紊流流动．力学学报．1998．7：475-481
52．汤广发．通风的有效性与室内空气品质．暖通空调．2000．5：20-23
53．黄以主．二维室内通风气流的数值解．南京建筑工程学院学报．1992．4：20-26
54．侯淑芬等．平战结合的高层建筑地下室通风与排烟设计．暖通空调，2001(1)：45～47
55．方祥位，李建中，魏文礼．紊流数值模拟的现状及发展．陕西水利发电．2000．3：16-19
56．是勋刚．湍流直接数值模拟的进展与前景．水动力学研究与进展．1992，A辑7(1)
57．马民仁，连之伟．置换通风几个问题的研究．暖通空调．2000．4：18-22
58．马民仁．论下部送风空调节能及其适用条件．暖通空调．1983．13(3)
59．王中铮，谷操，郭新川等．一种改进的SIMPLE算法．天津大学学报．1995．5：424-428
60．周鹏，李强民．置换通风与冷却顶板．暖通空调．1998．5：1-5
61．李世民．热车间自然通风的简化计算．暖通空调．1993．6：42-43
62．弓南，胡岚．地下汽车库通风排烟不同方式的比较．暖通空调，2001(2)：47～49
63．林其昌，吴时晶．新汽车库防火规范与高层建筑地下汽车库通风排烟设计．暖通空调．2001(6)：66～68
64．王宜义．空调房间气流组织的数值模拟．暖通空调．1993，5：34-37
65．周绍荣，楚惠东，陈汉平．混合对流作用下通风小室内流动换热的数值模拟．建筑热能 2000．1：14-17
66．张建忠．室内二维气流组织数值模拟．南京建筑工程学院学报．1993，3：10-16
67．张维功．失流洁净室室内流场的模型实验．哈尔滨建筑工程学院学报．1993．26(2)：77-81
68．丁卿．空调的节能与环保．江西建筑．1999．6：31-32
69．李先庭，江亿．用计算流体力学方法求解通风房间的空气年龄．清华大学学报(自然科学版)．1998．38(5)：28-31
70．陈斌，汤广发．三维紊流室内空气流动动态数值计算．湖南大学学报．1988．2
71．全国暖通空调制冷2000年学术文集．2000．10
72．全国暖通空渊制冷2000年学术年会．2000．10
73．江苏省制冷学会第四届年会论文集．南京：中国制冷学会．2000．8
74．2000年全国暖通空调年会专题交流与研讨．暖通空调，2001(1)：4～6
75．贺伟进．复杂几何边界流动数值计算．南京理工大学硕士论文．1999
76．杨含离．三维入水发射装置优化设计．南京理工大学硕士论文．1993
77．张向明．应用气动力原理研究拉格朗日问题．南京理工大学硕士论文．1991
78．乔学勇．复杂物体三维瞬态温度场边界元算法．南京理工大学硕士论文．1997
79．梅启元．侧送气流组织的数值模拟．南京理工大学硕士论文．2001
80. Brereton, G. J and Mankbadi, R. (1993).A rapid distortion theory turbulence model for developed unsteady wall-bounded flows. NASA-TM 106249
81. Christo, F. C., Masri A, R., and Nebot, E, M., (1996). An artificial neural network implementation of chemistry with pdf simulation of CO_2 flames. Combust Flames, 106,406
82. Sharath S. Girimaji. Lower-Dimensional manifold (algebraic) representation of Reynolds stress closure equations. Theoretical and Computational Fuild Dynamics. 259-281
83. Patankar S V, Spalding D B, A calculation procedure for heat, mass, and momentum transfer in three-dimensional flows[J].Int J Heat Mass Transfer, 1972,15:1787-1806
84. M F Brunk. Cooling ceilings-an opportunity to reduce energy costs by way of radiant cooling. ASHRAE Transactions, 1993. Vol 99
85. U Busweiler. Air conditioning with a combination of radiant cooling, displacement ventilation and desiccant cooling. ASHRAE Transactions, 1993. Vol 99
86. Lee, S.C (Univ of Missouri-Rolla, MO, USA); Sauer, HJ.J.r. Computer Simulation of turbulent mixing between outdoor air and return air in HVAC system. Eng. Div Publ FED, 1996 (2):237
87. H. M. Mathisen, Case studies of displacement ventilation in public halls, ASHRAE Transactions, 1989, vol. 97, Part 2
88. Ackerman. Mark Y (Univ of Alberta, Edmonton Can). Hoyme. Clifford L. Controlling the flow of combustion air and fresh air in housing. , 1996, 102(2)
89. Chow, Tin-Tai (City Unit of Hong Kong, Hong Kong), Ip, Fred, Dunn, Alan, Tse, wai leung. Numerical modeling of thermal behavior of fluid conduct flow with transport delay. ASHRAE Trans, 1996,102(2): 45-51.
90. T. Cziesla, G. Biswas, H. Chattopqdhyay and NiKMitra. Large-eddy simulation of flow and heat transfer in an impinging slot jet. International Journal of Heat and Fluid Flow, 22(2001): 500-8-508
91. Kirkpatrick, Allan T. (Colorado State Univ, first Collins. CO. OSA). Kenyon, Allen E, flow characteristics of three-dimensional wall jets. ASHRAE Trans V104 nPt 1B 1998 Proceedings of the ASHARE Winter Meeting, Part 2(of 2). San Francisco. CA. USA. ASHARE, Atlanta. GA, USA. P 1755-1762
92. HengSun and Suyi Huang. Simulation of Wind Flow Around a Building with a K-ε Model. Theoretical and Computational Fluid Dynamics, 2001(14): 283-292
93. Lauder, B. E., Recece, G. J and Rodi, W., J. Fluid Mech., 68, Parts(1975), 537-566
94. Fu, S.,Huang, P. G., Lauder, B. E, and Leschziner, M., ASME, J. of Fluid Eng.,110, (1988),216-274
95. Patankar, S. V., Sparrow, E. M., and Ivanovie, H., Int J. Heat Mass Transfer, 24, (1978),269-274
96. Humphrey, J, A. C.,Whitelaw, J.H and Yee, G., J, Fluid Mech.,103(1981), 443-463
2．陶文铨．数值传热学．第一版．西安：西安交通大学出版社．1988
3．陶文铨．计算流体力学与传热学．北京：中国建筑工业出版社．1992
4．帕坦卡[美]．传热与流体流动的数值计算．第一版．北京：科学出版社．1984
5．赵学端．粘性流体力学．第一版．北京：机械工业出版社．1983
6．潘文全．流体力学基础(上册)．第一版．北京：机械工业出版社．1982
7．潘文全．流体力学基础(下册)．第一版．北京：机械工业出版社．1982
8．尤国钊．计算流体力学．第一版．南京：南京理工大学．1994
9．马铁犹．计算流体力学．第一版．北京：航空学院出版社．1986．6
10．第六届全国计算流体力学会议论文集 科学出版社．1992
11．陈景仁．湍流模型及有限分析法．第一版．上海：上海交通大学出版社．1992
12．郭宽良．计算传热学．第一版．合肥：中国科学技术大学出版社．1988
13．章熙民．任泽霈．梅飞鸣．王中铮．传热学．第二版．北京：中国建筑工业出版社．1992
14．赵荣义．范存养．薛殿华．钱以明．空气调节．第三版 北京：中国建筑工业出版社．1994
15．周亨达．工程流体力学．第二版．北京．北京钢铁学院 1988
16．孙一坚．工业通风．第三版．北京．中国建筑工业出版社．1994
17．陆金甫．偏微分方程的数值解法．第一版．北京：清华大学出版社 1987
18．同济大学数学教研室．高等数学(上册)．第三版．上海：高等教育出版社．1995
19．同济大学数学教研室．高等数学(下册)．第三版．上海：高等教育出版社．1995
20．金练，欧阳嘀，张洁等．暖卫 通风 空调技术手册．第一版．北京：中国建筑工业出版社．2000．7
21．陆耀庆．暖通空调设计指南．第一版．北京：中国建筑工业出版社．1987
22．孙一坚．工业通风．第三版．北京：中国建筑工业出版社出版，1994
23．赵荣义，范存养，薛殿华，钱以明．空气调节．第三版．中国建筑工业出版社，1994
24．吴元炜．暖通空调新技术．第一版．北京：中国建筑工业出版社，1999
25．黄翔，连之伟，哈文．空调工程应用．第一版．北京：科学出版社，1999
26．王启杰．对流传热传质分析西安．第一版．西安：西安交通大学出版社．1991
27．陆耀庆．HAvC暖通空调设计指南．第一版．北京：中国建筑工业出版社，1996
28．刘旭，冯玉琪．实用空调技术精华．第一版．北京：人民邮电出版社，2001
29．潘云钢．高层民用建筑空调设计，第一版．北京：中国建筑工业出版社，1990
30．何耀东，何青．中央空调，第一版．北京：冶金工业出版社，1998
31．龚光彩．汤广发．吕文瑚．高大空间横向隔断气流的计算机仿真．湖南大学报．1998．25(4)：71-76
32．黄维平．高大门厅的空调送回风口设计探讨．制冷．1997．3：51-53
33．龚光彩．汤广发．高大空间纵向隔断气流的数值模拟．湖南大学学报(社科版)．1999．26(1)：66-71
34．何冰．空调房间气流组织的三维湍流计算．西安交通大学硕士论文
35．周绍荣，程运林，汤莉．高大空间气流组织的研究．湖南大学学报(自然科学版)．1998．25(2)：70-74
36．范存养．办公室下送风空调方式的应用．暖通空调．1997．27(4)：30-39
37．李湘洲，赵志江．日本的暖通空调及建筑节能现状．节能技术．1995．5：15-18
38．陈爱玲．空调房间气流组织的初步探讨．流体机械．1996．24(9)：51-55
39．王汉青，汤广发．三维紊流室内气流数值解．空气动力学报．1988．4(6)
40．魏学孟，樊洪明，张维功．失流洁净室数值模拟与实验研究．暖通空调．1997．27．(2)：11-14
41．李俊．空调的节能设计．住宅科技．1996．7：15-17
42．谭洪卫．计算流体动力学在建筑环境工程上的应用．暖通空调．1999，29(4)：31-36
43．龚光彩．CFD技术在暖通空调制冷工程中的应用．暖通空调．1999，29(6)：25-27
44．黄晨，李美玲，左涛等．大空间建筑室内垂直温度分布的计算方法．制冷学报．1999．3：25-30
45．陈在康．暖通CAD技术的开发和应用．暖通空调．1993．6
46．郭春信．室外采风口和排风口与系统的节能．1993．6
47．陈宾．日本大空间空调现状．暖通空调．1994．3
48．高养田．空调新风系统设计若干问题分析．暖通空调．1993．2：70-74
49．熊螯魁．不可压流场数值模拟方法综述．武汉交通科技大学学报．1999．10：523-528
50．张维，汤广发．风对建筑绕流流动的数值方法．空气动力学学报．1990．12：437-443
51．杨秉俭，蔡临宁．用代数应力模型模化型腔充填过程中的紊流流动．力学学报．1998．7：475-481
52．汤广发．通风的有效性与室内空气品质．暖通空调．2000．5：20-23
53．黄以主．二维室内通风气流的数值解．南京建筑工程学院学报．1992．4：20-26
54．侯淑芬等．平战结合的高层建筑地下室通风与排烟设计．暖通空调，2001(1)：45～47
55．方祥位，李建中，魏文礼．紊流数值模拟的现状及发展．陕西水利发电．2000．3：16-19
56．是勋刚．湍流直接数值模拟的进展与前景．水动力学研究与进展．1992，A辑7(1)
57．马民仁，连之伟．置换通风几个问题的研究．暖通空调．2000．4：18-22
58．马民仁．论下部送风空调节能及其适用条件．暖通空调．1983．13(3)
59．王中铮，谷操，郭新川等．一种改进的SIMPLE算法．天津大学学报．1995．5：424-428
60．周鹏，李强民．置换通风与冷却顶板．暖通空调．1998．5：1-5
61．李世民．热车间自然通风的简化计算．暖通空调．1993．6：42-43
62．弓南，胡岚．地下汽车库通风排烟不同方式的比较．暖通空调，2001(2)：47～49
63．林其昌，吴时晶．新汽车库防火规范与高层建筑地下汽车库通风排烟设计．暖通空调．2001(6)：66～68
64．王宜义．空调房间气流组织的数值模拟．暖通空调．1993，5：34-37
65．周绍荣，楚惠东，陈汉平．混合对流作用下通风小室内流动换热的数值模拟．建筑热能 2000．1：14-17
66．张建忠．室内二维气流组织数值模拟．南京建筑工程学院学报．1993，3：10-16
67．张维功．失流洁净室室内流场的模型实验．哈尔滨建筑工程学院学报．1993．26(2)：77-81
68．丁卿．空调的节能与环保．江西建筑．1999．6：31-32
69．李先庭，江亿．用计算流体力学方法求解通风房间的空气年龄．清华大学学报(自然科学版)．1998．38(5)：28-31
70．陈斌，汤广发．三维紊流室内空气流动动态数值计算．湖南大学学报．1988．2
71．全国暖通空调制冷2000年学术文集．2000．10
72．全国暖通空渊制冷2000年学术年会．2000．10
73．江苏省制冷学会第四届年会论文集．南京：中国制冷学会．2000．8
74．2000年全国暖通空调年会专题交流与研讨．暖通空调，2001(1)：4～6
75．贺伟进．复杂几何边界流动数值计算．南京理工大学硕士论文．1999
76．杨含离．三维入水发射装置优化设计．南京理工大学硕士论文．1993
77．张向明．应用气动力原理研究拉格朗日问题．南京理工大学硕士论文．1991
78．乔学勇．复杂物体三维瞬态温度场边界元算法．南京理工大学硕士论文．1997
79．梅启元．侧送气流组织的数值模拟．南京理工大学硕士论文．2001
80. Brereton, G. J and Mankbadi, R. (1993).A rapid distortion theory turbulence model for developed unsteady wall-bounded flows. NASA-TM 106249
81. Christo, F. C., Masri A, R., and Nebot, E, M., (1996). An artificial neural network implementation of chemistry with pdf simulation of CO_2 flames. Combust Flames, 106,406
82. Sharath S. Girimaji. Lower-Dimensional manifold (algebraic) representation of Reynolds stress closure equations. Theoretical and Computational Fuild Dynamics. 259-281
83. Patankar S V, Spalding D B, A calculation procedure for heat, mass, and momentum transfer in three-dimensional flows[J].Int J Heat Mass Transfer, 1972,15:1787-1806
84. M F Brunk. Cooling ceilings-an opportunity to reduce energy costs by way of radiant cooling. ASHRAE Transactions, 1993. Vol 99
85. U Busweiler. Air conditioning with a combination of radiant cooling, displacement ventilation and desiccant cooling. ASHRAE Transactions, 1993. Vol 99
86. Lee, S.C (Univ of Missouri-Rolla, MO, USA); Sauer, HJ.J.r. Computer Simulation of turbulent mixing between outdoor air and return air in HVAC system. Eng. Div Publ FED, 1996 (2):237
87. H. M. Mathisen, Case studies of displacement ventilation in public halls, ASHRAE Transactions, 1989, vol. 97, Part 2
88. Ackerman. Mark Y (Univ of Alberta, Edmonton Can). Hoyme. Clifford L. Controlling the flow of combustion air and fresh air in housing. , 1996, 102(2)
89. Chow, Tin-Tai (City Unit of Hong Kong, Hong Kong), Ip, Fred, Dunn, Alan, Tse, wai leung. Numerical modeling of thermal behavior of fluid conduct flow with transport delay. ASHRAE Trans, 1996,102(2): 45-51.
90. T. Cziesla, G. Biswas, H. Chattopqdhyay and NiKMitra. Large-eddy simulation of flow and heat transfer in an impinging slot jet. International Journal of Heat and Fluid Flow, 22(2001): 500-8-508
91. Kirkpatrick, Allan T. (Colorado State Univ, first Collins. CO. OSA). Kenyon, Allen E, flow characteristics of three-dimensional wall jets. ASHRAE Trans V104 nPt 1B 1998 Proceedings of the ASHARE Winter Meeting, Part 2(of 2). San Francisco. CA. USA. ASHARE, Atlanta. GA, USA. P 1755-1762
92. HengSun and Suyi Huang. Simulation of Wind Flow Around a Building with a K-ε Model. Theoretical and Computational Fluid Dynamics, 2001(14): 283-292
93. Lauder, B. E., Recece, G. J and Rodi, W., J. Fluid Mech., 68, Parts(1975), 537-566
94. Fu, S.,Huang, P. G., Lauder, B. E, and Leschziner, M., ASME, J. of Fluid Eng.,110, (1988),216-274
95. Patankar, S. V., Sparrow, E. M., and Ivanovie, H., Int J. Heat Mass Transfer, 24, (1978),269-274
96. Humphrey, J, A. C.,Whitelaw, J.H and Yee, G., J, Fluid Mech.,103(1981), 443-463