悬沙输移的三维有限元模型
|
摘要
泥沙问题是河口海岸地区研究的主要内容,对于河口治理、港口规划、航道整治、疏浚等有重要的意义。国内外学者对泥沙问题做了大量的研究,泥沙模型作为研究泥沙问题的主要手段之一,也得到了快速的发展,从一维、平面二维、垂向二维到准三维,不断的进行完善中。尽管在工程实践中也得到了广泛的应用,也只是近似的来模拟泥沙沿水深的变化规律。客观的自然现象中,水流输沙的运动具有明显的三维特性,而且河口海岸地区的潮流、波浪、径流等又增加了泥沙运动的复杂性。这就使得一维、二维泥沙模型的应用受到限制。因此,三维泥沙模型是泥沙研究发展的必然趋势。
本文针对河口海岸地区建立了悬沙输移的三维有限元模型:对三维悬沙对流扩散方程,空间域采用Galerkin有限元法进行离散,时间项采用显示差分法进行离散,把微分方程变为线性方程组,求解采用全选主元Gauss-Jordan消去法求解线性方程组;求解区域采用八结点的六面体等参单位剖分;通过分步积分和Green高斯公式得到满足自由表面和底部的边界条件;对本质边界条件(开边界条件)采用消行修正法进行修正,得到一个线性方程组。本文也介绍了悬沙模型中的各种参数理论的现状,很多是从试验中得到的半经验半理论公式,这就要求我们合理地选取参数来保证我们计算的准确性。
用Fortran语言对该模型编写程序,开发得到了可以计算悬沙模型的通用程序。最后采用有解析解的算例进行了验证,结果表明该模型可以用于海岸河口的悬沙输移模拟。
Sediment transport is an important content in estuary and coastal study, which has more sense in estuary regulation, port and harbor planning, channel improvement and dredging. Scholars around the world do lots of research on sediment transport. With the development of computer technology, sediment transport model is quickly developed. From 1-D to 2-D and 3-D, it is continuous optimization. Although these have been widely applied, in some cases they can’t simulate exactly. In the natural word, sediment movement is being of 3-d, further more, tidal current, wave and runoff is being more and more complicated. So the 1-D and 2-D sediment model are limited applied. Therrfore, the three-dimensional sediment model must be inevitable trend of sediment research.
A three-dimensional mathematical model is presented for the suspended sediment transport using Galerkin finite element method. For the governing equation, i.e. the convection-diffusion equation, the space discrete is performed via the standard Galerkin finite element method and the difference scheme is employed for the temporal discrete The boundary conditions on the free surface and bottom bed are introduced exactly by integrating the diffusion terms using integration by parts and Green's Theorem. The essential boundary conditions can be satisfied by modifying the coefficient matrix, so we get a solving linear equation. Also this paper introduces the present state of kinds of parameter theory in suspended sediment model, many of them is empirical formula, this requires us to reasonable choose parameters to ensure the accuracy of the calculation.
Through the debug, we get the procedure which is writed by Fortran computer language can calculate the suspended sediment transport. Finally, the model is verified by the examples of analytical solution, the results show that the model can be applied in the suspended sediment model of estuary and coastal.
本文针对河口海岸地区建立了悬沙输移的三维有限元模型:对三维悬沙对流扩散方程,空间域采用Galerkin有限元法进行离散,时间项采用显示差分法进行离散,把微分方程变为线性方程组,求解采用全选主元Gauss-Jordan消去法求解线性方程组;求解区域采用八结点的六面体等参单位剖分;通过分步积分和Green高斯公式得到满足自由表面和底部的边界条件;对本质边界条件(开边界条件)采用消行修正法进行修正,得到一个线性方程组。本文也介绍了悬沙模型中的各种参数理论的现状,很多是从试验中得到的半经验半理论公式,这就要求我们合理地选取参数来保证我们计算的准确性。
用Fortran语言对该模型编写程序,开发得到了可以计算悬沙模型的通用程序。最后采用有解析解的算例进行了验证,结果表明该模型可以用于海岸河口的悬沙输移模拟。
Sediment transport is an important content in estuary and coastal study, which has more sense in estuary regulation, port and harbor planning, channel improvement and dredging. Scholars around the world do lots of research on sediment transport. With the development of computer technology, sediment transport model is quickly developed. From 1-D to 2-D and 3-D, it is continuous optimization. Although these have been widely applied, in some cases they can’t simulate exactly. In the natural word, sediment movement is being of 3-d, further more, tidal current, wave and runoff is being more and more complicated. So the 1-D and 2-D sediment model are limited applied. Therrfore, the three-dimensional sediment model must be inevitable trend of sediment research.
A three-dimensional mathematical model is presented for the suspended sediment transport using Galerkin finite element method. For the governing equation, i.e. the convection-diffusion equation, the space discrete is performed via the standard Galerkin finite element method and the difference scheme is employed for the temporal discrete The boundary conditions on the free surface and bottom bed are introduced exactly by integrating the diffusion terms using integration by parts and Green's Theorem. The essential boundary conditions can be satisfied by modifying the coefficient matrix, so we get a solving linear equation. Also this paper introduces the present state of kinds of parameter theory in suspended sediment model, many of them is empirical formula, this requires us to reasonable choose parameters to ensure the accuracy of the calculation.
Through the debug, we get the procedure which is writed by Fortran computer language can calculate the suspended sediment transport. Finally, the model is verified by the examples of analytical solution, the results show that the model can be applied in the suspended sediment model of estuary and coastal.
引文
[1]时钟，周洪强．长江口北槽口外悬沙浓度垂线分布的数学模拟[J].海洋工程，2000，18(3)：57—62
[2]夏君强，王光谦．剖面二维悬移质泥沙输移方程的分步解法[J].长江科学院院报，2000，17(2)：14—17
[3]董文军，李世森，白玉川．三维潮流和潮流输沙问题的一种混合数值模拟及其应用[J].海洋学报，1999，21(2)：108—114
[4]朱志夏．二维水流泥沙数学模型理论及应用[D]．天津大学硕士论文，1995
[5]张修忠，王光谦．浅水流动及输运三维数学模型研究进展[J].水利学报，2002，(7)：1—7
[6]李孟国，蔡东明，张征等．海岸河口二维潮流可视化数学模型[J].海洋通报，2000，19(6)：57—65
[7]白玉川，于天一．分步分层拟三维水流数学模型及其在廉州湾潮流计算中的应用[J]_海洋学报，1998，20(5)：126—135
[8]刘桦，何友声．河口三维流动数学模型研究进展[J].海洋工程，2000，5(2)：87—93
[9]朱建荣，沈焕庭，周健．夏季苏北沿岸流对长江冲淡水扩展影响的[J].华东师范大学学报(自然科学版)，1997，(2)：62—67
[10]窦振兴，杨连武，Ozer，J．渤海三维潮流数值模拟[J]_海洋学报，1993，15(5)：1—15
[11]章本照．流体力学中的有限元方法[M]．北京：机械工业出版社，1986．
[12] Shi, Z., Li, S. S., Peterson, O. S. A vertical moving grid finite-element modeling of tidal flow in the Changjiang Estuary, China [J]. International Journal for Numerical Methods in Fluids, 2003, 43: 115-127
[13]时钟，李世森．垂向二维潮流数值模型及其在长江口北槽的应用[J].海洋通报，2003，22(3)：1—8
[14]吴克田．二维水波问题Galerkin解法的质量集中法[J].计算物理，1990，7(1)：1—6
[15] Zhao, C.G., Li, W.H., Wang J.T. et al. An explicit finite element method for dynamic analysis in three-medium coupling system and its application [J]. ACTA Seismologica Sinica, 2003, 16(3): 272-282
[16]董文军．河口近岸区三维潮流及悬沙扩散的一种分步模型[J]．天津大学学报，1999，32(3)：346—349
[17]董文军，陈虹．迎风有限元法在三维潮流数值模拟中的应用[J]．海洋与湖沼，1997，28(3)：320—327
[18]江春波，邢秀英，张庆海．用分步有限元法求解三维不可压缩流动[J]．清华大学学报(自然科学版)，2004，40(8)：110—113
[19]江春波，徐照明，李秀丽．求解不可压缩流动的分步有限元格式[J]．清华大学学报(自然科学版)，2002，42(2)：278—280
[20]朱志夏．波浪、潮流联合作用下泥沙数学模型的理论研究及其应用[D]．天津大学博士论文，1997
[21]刘玉玲，张宗孝．曲线拟合坐标变换技术在三角形网格生成技术中的应用[J]．陕西水力发电，2000，16(1)：13—15
[22]李褆来，窦希萍，黄晋鹏．长江口边界拟合坐标的三维潮流数学模型[J]．南京水利科学研究院水利水运科学研究，2000，(3)：1—6
[23]张瑞秋．二维水流数学模型的综合方法[D]．天津大学硕士论文，1987．
[24]张廷芳，段小宁．二维潮流数值计算的一种显式模式[J]．大连理工大学学报，1992，32(6)：702—706
[25]郭庆超，何明民．控制体积法在二维潮流计算中的应用[J]．水动力学研究与进展，A辑，1995，10(6)：602—609
[26]陈虹．泥质河口与海岸潮流泥沙数学模型的理论及应用[D]．天津大学博士论文，1997
[27]董文军．感潮泥质河口泥沙数学模型理论研究及应用[D]．天津大学博士论文，1996
[28] O'Connor, B. A., Nicholson, J. A three-dimensional model of suspended participate sediment transport[J]. Coastal Engineering, 1988, 12: 157-174
[29] Caneino, L., Neves, R. Hydrodynamic and sediment suspension modeling in Esturine systems: Part l:Desciption of the numerical models[J]. Journal of Marine Systems, 1999, 22(2): 105-116
[30] Lin, B.L., Falconer, R. Numerical modeling of three-dimensional suspended Sediment for estuarine and coastal waters. Journal of Hydraulic Research[J]. 1996, 34(4): 435-456
[31] Prinos P. Compound open flow with suspended sediments.Advance in Hydro-Science and Enfineering, Part B, USA, The University of Mississippi[J]. 1993, 1: 1206-1214
[32]丁平兴，史峰岩，孔亚珍．波一流共同作用下的三维悬沙扩散方程[J]．科学通报，1999，44(2)：1339—1342
[33]陈虹，李大鸣．三维潮流泥沙运动的一种数值模拟[J]．天津大学学报，1999，32(5)：573—279
[34]李蓓，唐士芳．河口海区开挖航道后三维潮流盐度泥沙数值模拟[J]．水道港口，2000，(4)：36—41
[35]江文胜，孙文心．渤海悬浮颗粒物的三维输运模式I模式[J]．海洋与湖沼，2001，3 1(6)：682—688
[36]江文胜，孙文心．渤海悬浮颗粒物的三维输运模式II模拟结果[J]．海洋与湖沼，2001，32(1)：94—100
[37]董文军，李世森，白玉川．三维潮流和潮流输沙问题的一种混合数值模拟及其应用[J]．海洋学报，1999，21(2)：108—114．
[38]王祥三．河口污染综合预测模型[J]．水利学报，1996，(8)：51—58
[39] Peter K Stansty, Peter M LLoyd. A semi-implicit langranian scheme for 3-D shollow water flow with a two-layer turbulence model[J]. International Journal for Numerical Method in Fluid. 1995, 20(2) : 115-134
[40] Davies AM. The numerical solution of the three dimensional hydrodynamic equations using a B-spline represention of the vertical current profiles, In : Bottom Turbulence[J]. Elservien Publication Corporation, 1997. 1-26
[41] Davies A M. Formulation of a linear 3-D hydrodynamic sea model using a Galerkin Eigen function method [J]. International Journal for Numerical Method in Fluid. 1983(3): 110-126
[42] Wang Sam S. Y. Finite element modeling of 3-D hydrodynamic and sediment transport phenomena (research report) [R]. The University of Mississippi, 1987, 6
[43]刘子龙，王船海．长江口三维水流模拟[J]．河海大学学报，1996(5)：108—110
[44] Thompson J F. Numrical solution of flow problem using body fitted coordinates system. In: Compuation Fluid Dynamics[J]. Hemisphere Publishing Corporation. 1980: 183-256
[45] Chenin Modojorich M. I. Numerical problems in coupling two- and three-dimensional models and turbulence measurements[J]. In : The 2nd S ymposium on refied flow modeling and turbulence Measurement . 1988. 302-305
[46]李孟国．海岸河口泥沙数学模型研究进展[J]．海洋工程．2006，24(1)：139—154
[47]陈国祥，陈界仁，沙捞．巴里．三维泥沙数学模型的研究进展[J]．水利水电科技进展．1998，18(1)：13—19
[48]曹祖德，孔令双，李蓓，李孟国．海岸河口水动力数值模拟的研究方向[J]．中国港湾建设．2002，4：15—18
[49]曹祖德，王运洪．水动力泥沙数值模拟[M]．天津大学出版社，1993
[50]窦希萍，罗肇森．波浪潮流共同作用下的二维泥沙数学模型[J]_水利水运科学研究，1992，(4)：331—336
[51]朱志夏，韩其为，丁平兴．海岸悬沙运移数学模型[J].海洋学报，2002，24(1)：101—107
[52]丁平兴，胡克林，孔亚珍．长江河口波一流共同作用下的全沙数值模拟[J]_海洋学报，2003，25(5)：113—124
[53]陈晓宏，陈永勤，赖国友．珠江口悬浮泥沙迁移数值模拟[J].海洋学报，2003，25(2)：120—126
[54]钱宁，万兆惠．泥沙运动力学[M]．北京：科学出版社，1986
[55]张瑞瑾．河流泥沙动力学[M]．北京：水利电力出版社，1988
[56] Wolanski, E. Transport of sediment in mangrove swamps [J]. Hydroniologia, 1995,295,31-42
[57] Jiang, J.H. An examination of estuarine lutocline dynamics [D]. University of Florida, 1999
[58] Guan, W.B. Transport and deposition of high-concentration suspensions of cohesive sediment in a macrotidal estuary [D]. Hong Kong University of Science and Technology, 2003
[59] Van Leussen, W. Aggregation of particles, settling velocities of mud floes: A review. In: Dronkers, J. and Van Leussen, W(eds): Physical Processes in Estuaries[J]. Springer Verlag, Berlin, Heidelberg, 1988, 347-403
[60] Mehta, A.J. Characterization of cohesive sediment properties and transport process in estuaries. In: Estuarine cohesive sediment dynamics[J]. Lecture Notes in Coastal and Estuarine Studies, No 14, Springer, Berlin, 1986, 290-325
[61] Van Leussen, W. Estuarine macrofiocs and their role in fine-grained sediment transport [D]. Ph.D. thesis, University of Utrecht, February 1994
[62] Malcherek, A., Markofsky, M, Zielke, W. et al. Three dimensional numerical modeling of cohesive sediment transport processes in estuarine environments [R]. Final report to the EC contract MAST2-CT92-0013, LNH(EDF-DER), Chatou, France, 1996
[63] Richardson, J.F., Zaki, W.N. Sedimentation and fluidization, Part I. Transactions of the Institution of Chemical Engineers, 1954, 32: 35-53.
[64]陈沈良，谷国传，张国安．长江口南汇近岸水域悬沙沉降速度估算[J].泥沙研究，2003，(6)：45—51
[65]时钟，朱文蔚，周洪强．长江口北槽口外细颗粒悬沙沉降速度[J].上海交通大学学报，2000，34(1)：18—23
[66]曹祖德，王运洪．水动力泥沙数值模拟[M]．天津：天津大学出版社，1994
[67] Hwang, K.N. Erodibility of fine sediment in wave-dominated environments [D]. M.S. Thesis, University of Florida, Gainesville, Florida, 1989, 158
[68]张东升，蒋勤．江苏北部灌河口悬沙输送数学模型[J].海洋学报，1991，13(1)：125—136
[69]蒋建华，苏纪兰．甬江建闸前后冲淤特性的初步数值模拟[J].海洋学报，1995，17(1)：121—129
[70]朱玉荣，常瑞芳．南黄海辐射沙洲区悬沙潮扩散规律数值研究[J]_青岛海洋大学学报，1998，28(4)：615—621
[71]余明辉，杨国录．平面二维非均匀沙数值模拟方法[J].水利学报，2000，(5)：65—69
[72]辛文杰．潮流、波浪综合作用下河口二维悬沙数学模型[J].海洋工程，1997，15(1)：30—47
[73]张华庆，李华国，岳翠平．海河口潮流泥沙运动数值模拟及清淤积方案研究[J].水动力研究与进展，2002，Ser．A，17(3)：318—326
[74]董文军．永定新河二维输沙有限元数值模拟[J].泥沙研究，1996，(4)：86—94
[75]李孟国，时钟，秦崇仁．伶仃洋三维潮流输沙的数值模拟[J].水利学报，2003，(4)：51—57
[76]匡翠萍．长江口拦门沙冲淤及悬沙沉降规律研究和水流盐度泥沙数学模型[D]．南京水利科学研究院，1993
[77]李孟国．海岸河口“五场”数学模型研究及其应用[R]．天津大学博士后研究工作报告，2004
[78]陈立，明宗富．河流动力学[M]．武汉：武汉大学出版社，2001
[79]窦国仁．再论泥沙起动流速[J]_泥沙研究，1999，(6)：1—9
[80]窦国仁，赵士清，黄亦芬．河道二维全沙数学模型研究[J].水利水运科学研究，1987，(2)：1—12
[81]窦希萍，李裎来，窦国仁．长江口全沙数学模型研究[J].水利水运科学研究，1999，(2)：136—145
[82]沙玉清．泥沙运动学引论[M]．陕西：陕西科学技术出版社，1996
[83]白玉川．潮流和波浪联合输沙的理论研究及其数学模型[D]．天津大学博士论文，1994
[84]曹振铁，胡克林．长江口二维非均匀悬沙数值模拟[J].泥沙研究，2002，(6)：66—73
[85]于清来，窦国仁．高含沙河流泥沙数学模型研究[J].水利水运科学研究，1999，(2)：107—115
[86] Partheniades, E. Erosion and deposition of cohesive soils [J]. Journal of the Hydrology Division, ASCE 91, 1965, (NY1): 105-139
[87] Krone, P. B. Flume studies of the transport in estuary shoaling processes [R]. Final report, Hydraulics Engineering Laboratory, University of Berkeley, California, USA, 1962
[88] Nicholson, J., O'connor, B.A. Cohesive sediment transport model [J]. Journal of Hydraulic Engineering, 1986, 112(7): 621-640
[89]曹祖德，王桂芬．波浪掀沙、潮流输沙的数值模拟[J].海洋学报，1993，15(1)：107—118
[90] Nadaoka, K., Yagi, H., Kamata, H. A simple quasi-3-D model of suspended sediment transport in a nonequilibrium state [J]. Coastal Engineering, 1991, 15:459-474
[91]朱建荣．海洋数值计算方法[M]．北京：海洋出版社，2003
[2]夏君强，王光谦．剖面二维悬移质泥沙输移方程的分步解法[J].长江科学院院报，2000，17(2)：14—17
[3]董文军，李世森，白玉川．三维潮流和潮流输沙问题的一种混合数值模拟及其应用[J].海洋学报，1999，21(2)：108—114
[4]朱志夏．二维水流泥沙数学模型理论及应用[D]．天津大学硕士论文，1995
[5]张修忠，王光谦．浅水流动及输运三维数学模型研究进展[J].水利学报，2002，(7)：1—7
[6]李孟国，蔡东明，张征等．海岸河口二维潮流可视化数学模型[J].海洋通报，2000，19(6)：57—65
[7]白玉川，于天一．分步分层拟三维水流数学模型及其在廉州湾潮流计算中的应用[J]_海洋学报，1998，20(5)：126—135
[8]刘桦，何友声．河口三维流动数学模型研究进展[J].海洋工程，2000，5(2)：87—93
[9]朱建荣，沈焕庭，周健．夏季苏北沿岸流对长江冲淡水扩展影响的[J].华东师范大学学报(自然科学版)，1997，(2)：62—67
[10]窦振兴，杨连武，Ozer，J．渤海三维潮流数值模拟[J]_海洋学报，1993，15(5)：1—15
[11]章本照．流体力学中的有限元方法[M]．北京：机械工业出版社，1986．
[12] Shi, Z., Li, S. S., Peterson, O. S. A vertical moving grid finite-element modeling of tidal flow in the Changjiang Estuary, China [J]. International Journal for Numerical Methods in Fluids, 2003, 43: 115-127
[13]时钟，李世森．垂向二维潮流数值模型及其在长江口北槽的应用[J].海洋通报，2003，22(3)：1—8
[14]吴克田．二维水波问题Galerkin解法的质量集中法[J].计算物理，1990，7(1)：1—6
[15] Zhao, C.G., Li, W.H., Wang J.T. et al. An explicit finite element method for dynamic analysis in three-medium coupling system and its application [J]. ACTA Seismologica Sinica, 2003, 16(3): 272-282
[16]董文军．河口近岸区三维潮流及悬沙扩散的一种分步模型[J]．天津大学学报，1999，32(3)：346—349
[17]董文军，陈虹．迎风有限元法在三维潮流数值模拟中的应用[J]．海洋与湖沼，1997，28(3)：320—327
[18]江春波，邢秀英，张庆海．用分步有限元法求解三维不可压缩流动[J]．清华大学学报(自然科学版)，2004，40(8)：110—113
[19]江春波，徐照明，李秀丽．求解不可压缩流动的分步有限元格式[J]．清华大学学报(自然科学版)，2002，42(2)：278—280
[20]朱志夏．波浪、潮流联合作用下泥沙数学模型的理论研究及其应用[D]．天津大学博士论文，1997
[21]刘玉玲，张宗孝．曲线拟合坐标变换技术在三角形网格生成技术中的应用[J]．陕西水力发电，2000，16(1)：13—15
[22]李褆来，窦希萍，黄晋鹏．长江口边界拟合坐标的三维潮流数学模型[J]．南京水利科学研究院水利水运科学研究，2000，(3)：1—6
[23]张瑞秋．二维水流数学模型的综合方法[D]．天津大学硕士论文，1987．
[24]张廷芳，段小宁．二维潮流数值计算的一种显式模式[J]．大连理工大学学报，1992，32(6)：702—706
[25]郭庆超，何明民．控制体积法在二维潮流计算中的应用[J]．水动力学研究与进展，A辑，1995，10(6)：602—609
[26]陈虹．泥质河口与海岸潮流泥沙数学模型的理论及应用[D]．天津大学博士论文，1997
[27]董文军．感潮泥质河口泥沙数学模型理论研究及应用[D]．天津大学博士论文，1996
[28] O'Connor, B. A., Nicholson, J. A three-dimensional model of suspended participate sediment transport[J]. Coastal Engineering, 1988, 12: 157-174
[29] Caneino, L., Neves, R. Hydrodynamic and sediment suspension modeling in Esturine systems: Part l:Desciption of the numerical models[J]. Journal of Marine Systems, 1999, 22(2): 105-116
[30] Lin, B.L., Falconer, R. Numerical modeling of three-dimensional suspended Sediment for estuarine and coastal waters. Journal of Hydraulic Research[J]. 1996, 34(4): 435-456
[31] Prinos P. Compound open flow with suspended sediments.Advance in Hydro-Science and Enfineering, Part B, USA, The University of Mississippi[J]. 1993, 1: 1206-1214
[32]丁平兴，史峰岩，孔亚珍．波一流共同作用下的三维悬沙扩散方程[J]．科学通报，1999，44(2)：1339—1342
[33]陈虹，李大鸣．三维潮流泥沙运动的一种数值模拟[J]．天津大学学报，1999，32(5)：573—279
[34]李蓓，唐士芳．河口海区开挖航道后三维潮流盐度泥沙数值模拟[J]．水道港口，2000，(4)：36—41
[35]江文胜，孙文心．渤海悬浮颗粒物的三维输运模式I模式[J]．海洋与湖沼，2001，3 1(6)：682—688
[36]江文胜，孙文心．渤海悬浮颗粒物的三维输运模式II模拟结果[J]．海洋与湖沼，2001，32(1)：94—100
[37]董文军，李世森，白玉川．三维潮流和潮流输沙问题的一种混合数值模拟及其应用[J]．海洋学报，1999，21(2)：108—114．
[38]王祥三．河口污染综合预测模型[J]．水利学报，1996，(8)：51—58
[39] Peter K Stansty, Peter M LLoyd. A semi-implicit langranian scheme for 3-D shollow water flow with a two-layer turbulence model[J]. International Journal for Numerical Method in Fluid. 1995, 20(2) : 115-134
[40] Davies AM. The numerical solution of the three dimensional hydrodynamic equations using a B-spline represention of the vertical current profiles, In : Bottom Turbulence[J]. Elservien Publication Corporation, 1997. 1-26
[41] Davies A M. Formulation of a linear 3-D hydrodynamic sea model using a Galerkin Eigen function method [J]. International Journal for Numerical Method in Fluid. 1983(3): 110-126
[42] Wang Sam S. Y. Finite element modeling of 3-D hydrodynamic and sediment transport phenomena (research report) [R]. The University of Mississippi, 1987, 6
[43]刘子龙，王船海．长江口三维水流模拟[J]．河海大学学报，1996(5)：108—110
[44] Thompson J F. Numrical solution of flow problem using body fitted coordinates system. In: Compuation Fluid Dynamics[J]. Hemisphere Publishing Corporation. 1980: 183-256
[45] Chenin Modojorich M. I. Numerical problems in coupling two- and three-dimensional models and turbulence measurements[J]. In : The 2nd S ymposium on refied flow modeling and turbulence Measurement . 1988. 302-305
[46]李孟国．海岸河口泥沙数学模型研究进展[J]．海洋工程．2006，24(1)：139—154
[47]陈国祥，陈界仁，沙捞．巴里．三维泥沙数学模型的研究进展[J]．水利水电科技进展．1998，18(1)：13—19
[48]曹祖德，孔令双，李蓓，李孟国．海岸河口水动力数值模拟的研究方向[J]．中国港湾建设．2002，4：15—18
[49]曹祖德，王运洪．水动力泥沙数值模拟[M]．天津大学出版社，1993
[50]窦希萍，罗肇森．波浪潮流共同作用下的二维泥沙数学模型[J]_水利水运科学研究，1992，(4)：331—336
[51]朱志夏，韩其为，丁平兴．海岸悬沙运移数学模型[J].海洋学报，2002，24(1)：101—107
[52]丁平兴，胡克林，孔亚珍．长江河口波一流共同作用下的全沙数值模拟[J]_海洋学报，2003，25(5)：113—124
[53]陈晓宏，陈永勤，赖国友．珠江口悬浮泥沙迁移数值模拟[J].海洋学报，2003，25(2)：120—126
[54]钱宁，万兆惠．泥沙运动力学[M]．北京：科学出版社，1986
[55]张瑞瑾．河流泥沙动力学[M]．北京：水利电力出版社，1988
[56] Wolanski, E. Transport of sediment in mangrove swamps [J]. Hydroniologia, 1995,295,31-42
[57] Jiang, J.H. An examination of estuarine lutocline dynamics [D]. University of Florida, 1999
[58] Guan, W.B. Transport and deposition of high-concentration suspensions of cohesive sediment in a macrotidal estuary [D]. Hong Kong University of Science and Technology, 2003
[59] Van Leussen, W. Aggregation of particles, settling velocities of mud floes: A review. In: Dronkers, J. and Van Leussen, W(eds): Physical Processes in Estuaries[J]. Springer Verlag, Berlin, Heidelberg, 1988, 347-403
[60] Mehta, A.J. Characterization of cohesive sediment properties and transport process in estuaries. In: Estuarine cohesive sediment dynamics[J]. Lecture Notes in Coastal and Estuarine Studies, No 14, Springer, Berlin, 1986, 290-325
[61] Van Leussen, W. Estuarine macrofiocs and their role in fine-grained sediment transport [D]. Ph.D. thesis, University of Utrecht, February 1994
[62] Malcherek, A., Markofsky, M, Zielke, W. et al. Three dimensional numerical modeling of cohesive sediment transport processes in estuarine environments [R]. Final report to the EC contract MAST2-CT92-0013, LNH(EDF-DER), Chatou, France, 1996
[63] Richardson, J.F., Zaki, W.N. Sedimentation and fluidization, Part I. Transactions of the Institution of Chemical Engineers, 1954, 32: 35-53.
[64]陈沈良，谷国传，张国安．长江口南汇近岸水域悬沙沉降速度估算[J].泥沙研究，2003，(6)：45—51
[65]时钟，朱文蔚，周洪强．长江口北槽口外细颗粒悬沙沉降速度[J].上海交通大学学报，2000，34(1)：18—23
[66]曹祖德，王运洪．水动力泥沙数值模拟[M]．天津：天津大学出版社，1994
[67] Hwang, K.N. Erodibility of fine sediment in wave-dominated environments [D]. M.S. Thesis, University of Florida, Gainesville, Florida, 1989, 158
[68]张东升，蒋勤．江苏北部灌河口悬沙输送数学模型[J].海洋学报，1991，13(1)：125—136
[69]蒋建华，苏纪兰．甬江建闸前后冲淤特性的初步数值模拟[J].海洋学报，1995，17(1)：121—129
[70]朱玉荣，常瑞芳．南黄海辐射沙洲区悬沙潮扩散规律数值研究[J]_青岛海洋大学学报，1998，28(4)：615—621
[71]余明辉，杨国录．平面二维非均匀沙数值模拟方法[J].水利学报，2000，(5)：65—69
[72]辛文杰．潮流、波浪综合作用下河口二维悬沙数学模型[J].海洋工程，1997，15(1)：30—47
[73]张华庆，李华国，岳翠平．海河口潮流泥沙运动数值模拟及清淤积方案研究[J].水动力研究与进展，2002，Ser．A，17(3)：318—326
[74]董文军．永定新河二维输沙有限元数值模拟[J].泥沙研究，1996，(4)：86—94
[75]李孟国，时钟，秦崇仁．伶仃洋三维潮流输沙的数值模拟[J].水利学报，2003，(4)：51—57
[76]匡翠萍．长江口拦门沙冲淤及悬沙沉降规律研究和水流盐度泥沙数学模型[D]．南京水利科学研究院，1993
[77]李孟国．海岸河口“五场”数学模型研究及其应用[R]．天津大学博士后研究工作报告，2004
[78]陈立，明宗富．河流动力学[M]．武汉：武汉大学出版社，2001
[79]窦国仁．再论泥沙起动流速[J]_泥沙研究，1999，(6)：1—9
[80]窦国仁，赵士清，黄亦芬．河道二维全沙数学模型研究[J].水利水运科学研究，1987，(2)：1—12
[81]窦希萍，李裎来，窦国仁．长江口全沙数学模型研究[J].水利水运科学研究，1999，(2)：136—145
[82]沙玉清．泥沙运动学引论[M]．陕西：陕西科学技术出版社，1996
[83]白玉川．潮流和波浪联合输沙的理论研究及其数学模型[D]．天津大学博士论文，1994
[84]曹振铁，胡克林．长江口二维非均匀悬沙数值模拟[J].泥沙研究，2002，(6)：66—73
[85]于清来，窦国仁．高含沙河流泥沙数学模型研究[J].水利水运科学研究，1999，(2)：107—115
[86] Partheniades, E. Erosion and deposition of cohesive soils [J]. Journal of the Hydrology Division, ASCE 91, 1965, (NY1): 105-139
[87] Krone, P. B. Flume studies of the transport in estuary shoaling processes [R]. Final report, Hydraulics Engineering Laboratory, University of Berkeley, California, USA, 1962
[88] Nicholson, J., O'connor, B.A. Cohesive sediment transport model [J]. Journal of Hydraulic Engineering, 1986, 112(7): 621-640
[89]曹祖德，王桂芬．波浪掀沙、潮流输沙的数值模拟[J].海洋学报，1993，15(1)：107—118
[90] Nadaoka, K., Yagi, H., Kamata, H. A simple quasi-3-D model of suspended sediment transport in a nonequilibrium state [J]. Coastal Engineering, 1991, 15:459-474
[91]朱建荣．海洋数值计算方法[M]．北京：海洋出版社，2003