平面二维明渠非恒定流的数值模拟
|
摘要
二维明渠水流数学模型的研究与运用己有相当长的历史，许多成果己经应用到水利、环境、港口航道等工程领域中，为工程的规划设计、施工及水域的可持续发展提供了重要的指导作用。随着国民经济的发展，更多的江河水力资源得到了开发利用，对明渠水流数学模型也提出了新的要求，同时也为水流运动基本理论研究的发展提供了一个很好的机遇和条件，以便明渠水流模型能够更好地服务于生产实践。本文利用边界拟合坐标技术生成正交曲线网格，基于此网格系统，建立平面二维明渠非恒定流的数值模型，实现复杂边界明渠水流的数值模拟。主要研究内容如下：
1、采用边界拟合坐标技术建立二维正交曲线网格以克服由于复杂边界而引起的计算困难。提出了一种经济、有效的滑动边界处理方法，提高网格边界上的正交性。
2、利用张量分析的方法，推导出正交曲线坐标下二维浅水方程的具体形式。
3、采用有限差分法对二维浅水方程进行离散，利用交替方向隐格式法(Alternating Direction Implicit Method简称ADI法)实现在计算区域内对二维浅水方程的求解。将得到的数学模型用于长江南通河段进行算例验证，验证结果比较令人满意。针对紊动粘性系数ε这个重要参数，对该算例进行分析验证，验证结果表明模型中考虑ε结果更好。
4、在得到的数学模型中分别加入迭代的思想和迎风格式，以此对模型进行修正。实际算例表明修正是合理的，有助于提高模型的计算精度。
The research and application of two-dimensional open channel flow mathematical model has been a long history, and many researches have been applied to such projects as water conservancy, environment, port and channel, flood defense, which plays an important guiding role in the design and construction for a project and sustainable development of waters. With the development of national economy, more hydraulic power resources have been further exploited and utilized, which makes a higher demand for open channel flow mathematical model. Meanwhile, it also provides a good chance and prerequisite for the development of fundamental research of water current and gives a better service to the production practice. This thesis is going to establish a horizontal two-dimensional mathematical model of unsteady open channel flow under curvilinear coordinate system on the basis of the proposed boundary-fitted coordinate technology and try to obtain the flow fields with complex boundaries. The research details of the thesis are as follows.
1. An orthogonal curvilinear gird has been established by employing the boundary-fitted coordinate technology, which aims to overcome difficulties of numerical simulation resulting from complex boundary. Besides, an economical and effective sliding boundary processing method of improving the orthogonality of the boundary has been proposed.
2. By means of tensor analysis, shallow water wave equations of two dimensional unsteady channel flow in orthogonal curvilinear coordinates have been deducted.
3. In the computational domain, shallow water wave equations of two dimensional unsteady channel flow have been dispersed by finite difference method and solved by the alternating direction implicit technique. The established mathematical model mentioned above is applied to test and verify Nantong section of Yangtze River, the result of which proves to be satisfactory. And if it is testified by the eddy viscosity coefficient, the result turns out to be much better.
4. Iteration method and upwind scheme have been respectively added to the mathematical mode mentioned above to modify the model, and the result of which proves to be reasonable and helpful to improve the computational accuracy of the model.
1、采用边界拟合坐标技术建立二维正交曲线网格以克服由于复杂边界而引起的计算困难。提出了一种经济、有效的滑动边界处理方法，提高网格边界上的正交性。
2、利用张量分析的方法，推导出正交曲线坐标下二维浅水方程的具体形式。
3、采用有限差分法对二维浅水方程进行离散，利用交替方向隐格式法(Alternating Direction Implicit Method简称ADI法)实现在计算区域内对二维浅水方程的求解。将得到的数学模型用于长江南通河段进行算例验证，验证结果比较令人满意。针对紊动粘性系数ε这个重要参数，对该算例进行分析验证，验证结果表明模型中考虑ε结果更好。
4、在得到的数学模型中分别加入迭代的思想和迎风格式，以此对模型进行修正。实际算例表明修正是合理的，有助于提高模型的计算精度。
The research and application of two-dimensional open channel flow mathematical model has been a long history, and many researches have been applied to such projects as water conservancy, environment, port and channel, flood defense, which plays an important guiding role in the design and construction for a project and sustainable development of waters. With the development of national economy, more hydraulic power resources have been further exploited and utilized, which makes a higher demand for open channel flow mathematical model. Meanwhile, it also provides a good chance and prerequisite for the development of fundamental research of water current and gives a better service to the production practice. This thesis is going to establish a horizontal two-dimensional mathematical model of unsteady open channel flow under curvilinear coordinate system on the basis of the proposed boundary-fitted coordinate technology and try to obtain the flow fields with complex boundaries. The research details of the thesis are as follows.
1. An orthogonal curvilinear gird has been established by employing the boundary-fitted coordinate technology, which aims to overcome difficulties of numerical simulation resulting from complex boundary. Besides, an economical and effective sliding boundary processing method of improving the orthogonality of the boundary has been proposed.
2. By means of tensor analysis, shallow water wave equations of two dimensional unsteady channel flow in orthogonal curvilinear coordinates have been deducted.
3. In the computational domain, shallow water wave equations of two dimensional unsteady channel flow have been dispersed by finite difference method and solved by the alternating direction implicit technique. The established mathematical model mentioned above is applied to test and verify Nantong section of Yangtze River, the result of which proves to be satisfactory. And if it is testified by the eddy viscosity coefficient, the result turns out to be much better.
4. Iteration method and upwind scheme have been respectively added to the mathematical mode mentioned above to modify the model, and the result of which proves to be reasonable and helpful to improve the computational accuracy of the model.
引文
[1] 周雪漪．计算水力学[M]．北京：清华大学出版社，1995．
[2] Pieter Wesseling．计算流体力学原理(影印版)[M]．北京：科学出版社，2006．
[3] 郑邦民，赵昕．计算水动力学[M]．武汉：武汉大学出版社，2002．
[4] 刘儒勋，舒其望．计算流体力学的若干新方法[M]．北京：科学出版社，2003．
[5] 槐文信，赵明登，童汉毅．河道及近海水流的数值模拟[M]．北京：科学出版社，2005．
[6] 傅德薰，马延文．计算流体力学[M]．北京：高等教育出版社，2002．
[7] 夏军强，王光谦，吴保生．游荡型河流演变及其数值模拟[M]．北京：中国水利水电出版社，2005．
[8] 冯民权．淹没式温排水试验及三维紊动扩散数值模拟[D]．武汉：武汉水利电力学院，1988．
[9] 曹祖德等．水动力泥沙数值模拟[M]．天津：天津大学出版社，1994．
[10] 姜礼尚等．有限元方法及理论基础[M]．北京：人民教育出版，1979．
[11] 陈景仁．湍流模型及有限分析法[M]．上海：上海交通大学出版社，1988．
[12] C.J.Chen. The finite analytic method in flows and heat transfer[D].The university of Iowa, USA, 1987.
[13] 陈景仁著．流体力学及传热学[M]．北京：国防工业出版社，1987．
[14] 李炜，对流扩散方程与边界层层流流动的混合有限分析[J]，武汉水利电力大学学报，1994．
[15] 刘希云，赵润祥．流体力学中的有限元与边界元方法[M]．上海：上海交通大学出版社，1993．
[16] 谭维炎．计算浅水动力学—有限体积法的应用[M]．北京：清华大学出版社，1998．
[17] 陶文铨．数值传热学[M]．西安：西安交通大学出版社，1988．
[18] 陶文铨．计算传热学的最近进展[M]．西安：西安交通大学出版社．
[19] Patanker，S．V著，张政译．传热与流体流动的数值计算[M]．北京：科学出版社，1992．
[20] Thompson F. Computational Fluid Dynamics[M]. MCGRAW-HILL INTERNATIONAL BOOK COMPANY, 1980.
[21] Peacemen D W, Rachford H H, The Numerical Solution of Parabolic and Elliptic Differential Equations[J], J SIAM, 1995.3:28-41.
[22] Douglas J, Jr, Alternating Direct Method for Three Space Variables[J], NumeriseheMathematic, 1962 (4):41-63.
[23] 刘光宗编．流体力学原理与计算基础[M]．西安：西安交通大学出版社，1988．
[24] Stelling G S. On the Construction of Computational Methods for Shallow Water Flow Problems [D].Delft: Delft University of Technology, the Netherlands, 1983.
[25] Patankar SV, Spalding DB. Calculation Procedure for Heat, Mass and Momentum Transfer in 3 D Flows [J], Int. J.Heat Mass Transfer, 1972 (15): 1806-1987.
[26] Patankar SV. Numerical Heat Transfer and Fluid Flow[M].New York: Me Graw Hill,1980.
[27] Van Doormaal J P.Raithby GD, Enhancement of SIMPLE Method for Predicting Incompressible Fluid Flows [J], Numer Heat Transfer, 1984(7): 147-163.
[28] Raithby GD, Schneider G E. Elliptic Systems:Finite Difference Methods Ⅱ [M].New York:John Wiley & Sons, 1998.
[29] Sheng Y, Shoukri M,Sheng G, Wood P,A modification to the SIMPLE Method for Buoyancy2driven Flows [J],Numer Heat Transfer, Part B, 1998,33 (1):65-78.
[30] Jian Guo Zhou, Velocity2Depvh Coupling in Shallow2WaterFlows[J], J. Hydr. Engrg, ASCE,1995,121 (10):717-724.
[31] HoltM. Numerical Methods in Fluid Dynamics [M].Berlin: Springer-Verlag, 1984.
[32] 王船海，程文辉，河道二维非恒定流场计算方法研究[J]，水利学报，1991(1)：10-17．
[33] 王船海，程文辉，河道二维非恒定流计算[J]，河海大学学报，1987，15(3)：39-53．
[34] 程文辉，王船海，用正交曲线网格及“冻结”法计算河道流速场[J]，水利学报，1988(6)：13-24．
[35] 周建军，林秉南，王连祥，河道平面二维水流数值计算[J]，水利学报，1991(5)：8-18．
[36] 吴修广，沈永明等，非正交曲线坐标下二维水流计算的SIMPLEC算法[J]，水利学报，2003(2)：25-30．
[37] 方春明，全隐式差分法求解河道平面二维恒定水流运动方程[J]，水利学报，1997(4)：42-48．
[38] 刘晓东，华祖林，赵玉萍，基于四叉树网格的GODUNOV型二维水流数值计算模式[J]，河海大学学报，2002，30(6)：7-10．
[39] 李光炽，周晶晏，张贵寿，高桩码头对河道流场影响的数值模拟[J]，河海大学学报，2004，32(2)：217-220．
[40] 邓家泉，二维明渠非恒定水流BGK数值模型[J]，水利学报，2002(4)：1-7．
[41] 胡四一，谭维炎，无结构网格上二维浅水流动的数值模拟[J]，水科学进展，1995，6(1)：1-9．
[42] 王如云，张东生等，曲线坐标网格下二维涌波数值模拟的TVD型格式[J]，水利学报，2002(10)：72-77．
[43] 邵颂东，王光谦，费详俊，平面二维Lagrange-Euler方法及其在水流运动中的应 用[J]，水利学报，1997(6)：34-39．
[44] 程永光，索丽生，二维明渠非恒定流的格子Boltzmann模拟[J]，水科学进展，2003，14(9)：9-14．
[45] 魏文礼，沈永明，二维溃坝洪水波演进的数值模拟[J]，水力学报，2003(9)：43-47．
[46] Thompson J F. Boundary-fitted Coordinate System for Numerical Solution of Partial Differential Equation[J],J ComputerPhys, 1982,47:1-108.
[47] 余利仁，许卫新，张书农，正交贴体坐标系的生成[J]，河海大学学报，1988，16(5)：78-85．
[48] 何国建，汪德爟，刘晓波，定点边界二维正交曲线网格的数值生成方法[J]，河海大学学报，2004，32(2)：140-143．
[49] 揭冠周，介玉新，李广信，模拟自由面渗流的适体坐标变换方法[J]，清华大学学报，2003，43(2)：273-284．
[50] 魏文礼，王德意，刘玉玲，正交曲线网格生成技术的一点改进[J]，武汉水力水电大学学报，2000，33(2)：40-42．
[51] 董耀华，河势贴体河道平面二维正交网格生成方法的研究及应用[J]，长江科学院院报，2001，18(4)：14-17．
[52] 黄炳彬，方红卫，刘斌，复杂边界水流数学模型的斜对角笛卡尔方法[J]，水动力学研究与进展，2003，18(6)：679-685．
[53] 刘玉玲，张宗孝，魏文礼，刘明琴，曲线拟合坐标变换技术在三角形网格生成中的应用[J]，陕西水力发电，2000，16(1)：13-15．
[54] 吴江航，韩庆书．计算流体力学的理论、方法及应用[M]．北京：科学出版社，1988．
[55] W.Kollmann, Computafivnal Fluid Dynamics[J], Henisp Here Publishing CorpoRaition,1980:1-98.
[56] J.F. Thompson et al, TOMCAT-A Code for Numerical Generation of Bounday—Fitted Curvilinear Coordinate system on Fields Containing and Number of Arbitrary Tow-Dimensional Bodies[J],J. of Computational physics, 1977,24(274).
[57] 陈景仁．流体力学及传热学[M]．北京：国防工业出版社，1984．
[58] 陶文铨．数值传热学(第2版)[M]．西安：西安交通大学出版社，2001．
[59] 汪德爟．计算水力学理论与应用[M]．南京：河海大学出版社，1989．
[60] 魏文礼．计算水力学理论及应用[M]．西安：陕西科学技术出版社，2001：135-136．
[61] 黄克智．张量分析[M]．北京：清华大学出版社，2003．
[62] 蔡国琰，刘存芳，马春元．张量运算及其应用[M]．山东：山东大学出版社，2000．
[63] 余天庆．张量分析及演算[M]．湖北：华中理工大学出版社，1996．
[64] 郑邦民，黄克中，魏良琰．流体力学[M]．北京：水利电力出版社，1989．
[65] 李开泰，黄艾香．张量分析及其应用[M]．西安：西安交通大学出版，1983．
[66] 冯潮清，赵愉深，何浩法．矢量与张量分析[M]．北京：国防工业出版社，1986．
[67] 张涤明．计算流体力学[M]．广州：中山大学出版社，1991．
[68] 何杰．潮汐河口汊道治理数值模拟研究[D]．南京：南京水利科学研究所，2006．
[69] 李洪声．河口、海岸泥沙输移数值模拟[D]．大连：大连理工大学，2006．
[70] 何少荃，I．R．Warren，平面二维明渠流动数值计算中涡粘系数项的合理选取[J]，水利学报，1990(4)：1-10．
[2] Pieter Wesseling．计算流体力学原理(影印版)[M]．北京：科学出版社，2006．
[3] 郑邦民，赵昕．计算水动力学[M]．武汉：武汉大学出版社，2002．
[4] 刘儒勋，舒其望．计算流体力学的若干新方法[M]．北京：科学出版社，2003．
[5] 槐文信，赵明登，童汉毅．河道及近海水流的数值模拟[M]．北京：科学出版社，2005．
[6] 傅德薰，马延文．计算流体力学[M]．北京：高等教育出版社，2002．
[7] 夏军强，王光谦，吴保生．游荡型河流演变及其数值模拟[M]．北京：中国水利水电出版社，2005．
[8] 冯民权．淹没式温排水试验及三维紊动扩散数值模拟[D]．武汉：武汉水利电力学院，1988．
[9] 曹祖德等．水动力泥沙数值模拟[M]．天津：天津大学出版社，1994．
[10] 姜礼尚等．有限元方法及理论基础[M]．北京：人民教育出版，1979．
[11] 陈景仁．湍流模型及有限分析法[M]．上海：上海交通大学出版社，1988．
[12] C.J.Chen. The finite analytic method in flows and heat transfer[D].The university of Iowa, USA, 1987.
[13] 陈景仁著．流体力学及传热学[M]．北京：国防工业出版社，1987．
[14] 李炜，对流扩散方程与边界层层流流动的混合有限分析[J]，武汉水利电力大学学报，1994．
[15] 刘希云，赵润祥．流体力学中的有限元与边界元方法[M]．上海：上海交通大学出版社，1993．
[16] 谭维炎．计算浅水动力学—有限体积法的应用[M]．北京：清华大学出版社，1998．
[17] 陶文铨．数值传热学[M]．西安：西安交通大学出版社，1988．
[18] 陶文铨．计算传热学的最近进展[M]．西安：西安交通大学出版社．
[19] Patanker，S．V著，张政译．传热与流体流动的数值计算[M]．北京：科学出版社，1992．
[20] Thompson F. Computational Fluid Dynamics[M]. MCGRAW-HILL INTERNATIONAL BOOK COMPANY, 1980.
[21] Peacemen D W, Rachford H H, The Numerical Solution of Parabolic and Elliptic Differential Equations[J], J SIAM, 1995.3:28-41.
[22] Douglas J, Jr, Alternating Direct Method for Three Space Variables[J], NumeriseheMathematic, 1962 (4):41-63.
[23] 刘光宗编．流体力学原理与计算基础[M]．西安：西安交通大学出版社，1988．
[24] Stelling G S. On the Construction of Computational Methods for Shallow Water Flow Problems [D].Delft: Delft University of Technology, the Netherlands, 1983.
[25] Patankar SV, Spalding DB. Calculation Procedure for Heat, Mass and Momentum Transfer in 3 D Flows [J], Int. J.Heat Mass Transfer, 1972 (15): 1806-1987.
[26] Patankar SV. Numerical Heat Transfer and Fluid Flow[M].New York: Me Graw Hill,1980.
[27] Van Doormaal J P.Raithby GD, Enhancement of SIMPLE Method for Predicting Incompressible Fluid Flows [J], Numer Heat Transfer, 1984(7): 147-163.
[28] Raithby GD, Schneider G E. Elliptic Systems:Finite Difference Methods Ⅱ [M].New York:John Wiley & Sons, 1998.
[29] Sheng Y, Shoukri M,Sheng G, Wood P,A modification to the SIMPLE Method for Buoyancy2driven Flows [J],Numer Heat Transfer, Part B, 1998,33 (1):65-78.
[30] Jian Guo Zhou, Velocity2Depvh Coupling in Shallow2WaterFlows[J], J. Hydr. Engrg, ASCE,1995,121 (10):717-724.
[31] HoltM. Numerical Methods in Fluid Dynamics [M].Berlin: Springer-Verlag, 1984.
[32] 王船海，程文辉，河道二维非恒定流场计算方法研究[J]，水利学报，1991(1)：10-17．
[33] 王船海，程文辉，河道二维非恒定流计算[J]，河海大学学报，1987，15(3)：39-53．
[34] 程文辉，王船海，用正交曲线网格及“冻结”法计算河道流速场[J]，水利学报，1988(6)：13-24．
[35] 周建军，林秉南，王连祥，河道平面二维水流数值计算[J]，水利学报，1991(5)：8-18．
[36] 吴修广，沈永明等，非正交曲线坐标下二维水流计算的SIMPLEC算法[J]，水利学报，2003(2)：25-30．
[37] 方春明，全隐式差分法求解河道平面二维恒定水流运动方程[J]，水利学报，1997(4)：42-48．
[38] 刘晓东，华祖林，赵玉萍，基于四叉树网格的GODUNOV型二维水流数值计算模式[J]，河海大学学报，2002，30(6)：7-10．
[39] 李光炽，周晶晏，张贵寿，高桩码头对河道流场影响的数值模拟[J]，河海大学学报，2004，32(2)：217-220．
[40] 邓家泉，二维明渠非恒定水流BGK数值模型[J]，水利学报，2002(4)：1-7．
[41] 胡四一，谭维炎，无结构网格上二维浅水流动的数值模拟[J]，水科学进展，1995，6(1)：1-9．
[42] 王如云，张东生等，曲线坐标网格下二维涌波数值模拟的TVD型格式[J]，水利学报，2002(10)：72-77．
[43] 邵颂东，王光谦，费详俊，平面二维Lagrange-Euler方法及其在水流运动中的应 用[J]，水利学报，1997(6)：34-39．
[44] 程永光，索丽生，二维明渠非恒定流的格子Boltzmann模拟[J]，水科学进展，2003，14(9)：9-14．
[45] 魏文礼，沈永明，二维溃坝洪水波演进的数值模拟[J]，水力学报，2003(9)：43-47．
[46] Thompson J F. Boundary-fitted Coordinate System for Numerical Solution of Partial Differential Equation[J],J ComputerPhys, 1982,47:1-108.
[47] 余利仁，许卫新，张书农，正交贴体坐标系的生成[J]，河海大学学报，1988，16(5)：78-85．
[48] 何国建，汪德爟，刘晓波，定点边界二维正交曲线网格的数值生成方法[J]，河海大学学报，2004，32(2)：140-143．
[49] 揭冠周，介玉新，李广信，模拟自由面渗流的适体坐标变换方法[J]，清华大学学报，2003，43(2)：273-284．
[50] 魏文礼，王德意，刘玉玲，正交曲线网格生成技术的一点改进[J]，武汉水力水电大学学报，2000，33(2)：40-42．
[51] 董耀华，河势贴体河道平面二维正交网格生成方法的研究及应用[J]，长江科学院院报，2001，18(4)：14-17．
[52] 黄炳彬，方红卫，刘斌，复杂边界水流数学模型的斜对角笛卡尔方法[J]，水动力学研究与进展，2003，18(6)：679-685．
[53] 刘玉玲，张宗孝，魏文礼，刘明琴，曲线拟合坐标变换技术在三角形网格生成中的应用[J]，陕西水力发电，2000，16(1)：13-15．
[54] 吴江航，韩庆书．计算流体力学的理论、方法及应用[M]．北京：科学出版社，1988．
[55] W.Kollmann, Computafivnal Fluid Dynamics[J], Henisp Here Publishing CorpoRaition,1980:1-98.
[56] J.F. Thompson et al, TOMCAT-A Code for Numerical Generation of Bounday—Fitted Curvilinear Coordinate system on Fields Containing and Number of Arbitrary Tow-Dimensional Bodies[J],J. of Computational physics, 1977,24(274).
[57] 陈景仁．流体力学及传热学[M]．北京：国防工业出版社，1984．
[58] 陶文铨．数值传热学(第2版)[M]．西安：西安交通大学出版社，2001．
[59] 汪德爟．计算水力学理论与应用[M]．南京：河海大学出版社，1989．
[60] 魏文礼．计算水力学理论及应用[M]．西安：陕西科学技术出版社，2001：135-136．
[61] 黄克智．张量分析[M]．北京：清华大学出版社，2003．
[62] 蔡国琰，刘存芳，马春元．张量运算及其应用[M]．山东：山东大学出版社，2000．
[63] 余天庆．张量分析及演算[M]．湖北：华中理工大学出版社，1996．
[64] 郑邦民，黄克中，魏良琰．流体力学[M]．北京：水利电力出版社，1989．
[65] 李开泰，黄艾香．张量分析及其应用[M]．西安：西安交通大学出版，1983．
[66] 冯潮清，赵愉深，何浩法．矢量与张量分析[M]．北京：国防工业出版社，1986．
[67] 张涤明．计算流体力学[M]．广州：中山大学出版社，1991．
[68] 何杰．潮汐河口汊道治理数值模拟研究[D]．南京：南京水利科学研究所，2006．
[69] 李洪声．河口、海岸泥沙输移数值模拟[D]．大连：大连理工大学，2006．
[70] 何少荃，I．R．Warren，平面二维明渠流动数值计算中涡粘系数项的合理选取[J]，水利学报，1990(4)：1-10．