无结构网格二维河口海洋数值模式的研制.docVIP

无结构网格二维河口海洋数值模式的研制 【摘要】：数值模拟在科学研究、工程应用与灾害预测分析等方面发挥着重要作用。海洋数值模式是对河流、海洋等水体在计算机上模拟的工具,正得到蓬勃发展。传统的海洋数值模式大多采用矩形结构网格,局部加密及拟合岸线能力有限,往往难以拟合河口、海湾或岛屿等区域复杂的岸线和工程结构。无结构网格能够在任意局部区域对网格加密并拟合岸线,从而克服了结构网格的缺陷。当前国际上知名的无结构网格海洋模式大多来自国外,而国内自主开发的模式尚缺少能达到国际先进水平的。所以需要研制具有完全自主知识产权的无结构网格河口海洋数值模式,以期通过不断改进发展最终达到国际先进水平。本文研制了一套无结构网格二维河口海洋数值模式。该模式采用三角形网格,基于有限体积法求解。在网格中心(重心)求解水位,在网格边的中点同时求解x-方向和y-方向流速的配置。这样的配置既不用水平坐标转换,能使用有限体积法直接对原始方程组进行求解,简洁高效且易于守恒,同时也能够较方便地处理移动潮滩边界。模式采用垂向平均的二维控制方程组,对连续方程、动量方程和盐度方程均采用有限体积法求解,以确保守恒性。求解连续方程和盐度方程时采用三角形网格作为控制体。求解动量方程时则采用围绕边中点的菱形作为控制体,并采用TVD格式处理流速平流项以增强稳定性。设计了两种方法求解正压梯度力,第一种方法为对水位插值,再在菱形控制体中用有限体积法求解,第二种方法为先构造4个以水位计算点为顶点的控制体,用有限体积法求解出水位梯度后,再将所得水位梯度平均。斜压梯度力求法与正压梯度力类似。求解盐度平流项时,根据Darwish和Moukalled在无结构网格下设计物质输运TVD平流项算法的经验,以及Li等对其不足的分析和改进,针对本文模式的无结构网格设计了一种计算上游r值的方法,并以此r值为基础构造了TVD格式。模式在时间上用预估修正法显式求解以增加水动力精度,并且包含了移动潮滩边界。通过在一个理想圆湖试验,验证了第二种正压梯度力求解方法更加精确,计算结果与解析解非常接近。而通过理想往复流数值试验测试TVD盐度平流项格式,测试结果表明该格式成功地消除了数值频散,并有效地降低了数值耗散,但使用Superbee限制函数会造成人为的锋面,故推荐使用Monotonic和Minimode限制函数。对黄浦江进行水动力验证,水位结果与实测结果基本一致。在长江口对模式作了率定验证。模式在不同区域取不同底摩擦曼宁系数值,通过计算2007年8月长江口的流场,对比实测水位、流速和流向率定,并通过计算2004年末2005年初的流场和盐度场,对比实测流速、流向和盐度验证。率定结果整体较好,验证的流速流向较为接近,但盐度偏差略大,这是由于二维模式无法再现三维盐水楔现象,以及盐度变化与初始场等有关。通过各种测试和率定验证,可得出结论：模式的正压梯度力计算和盐度平流项格式均具有较高的精度,可应用于河口与近海区域的数值计算。【关键词】：无结构网格数值模式有限体积法TVD格式数值频散数值耗散 【学位授予单位】：华东师范大学 【学位级别】：博士 【学位授予年份】：2011 【分类号】：P714 【目录】：摘要7-9ABSTRACT9-12第一章绪论12-201.1无结构网格河口海洋数值模式的研究意义12-131.2基于有限体积法的无结构网格海洋模式的国内外研究进展13-171.3本文的工作17-20第二章数值模式的建立20-372.1控制方程组20-222.2网格配置和编号22-242.2.1网格配置22-232.2.2网格编号23-242.3方程离散与求解24-352.3.1求解连续方程24-252.3.2求解动量方程25-272.3.3TVD格式简介27-292.3.4动量方程平流项292.3.5动量方程扩散项29-302.3.6正压梯度力项30-312.3.7斜压梯度力项312.3.8盐度方程求解31-332.3.9预估修正法处理水动力332.3.10移动潮滩边界和边界条件33-352.4模式主程序的计算流程35-37第三章理想试验和简单河道验证37-593.1理想圆湖试验37-493.1.1圆湖模型简介373.1.2本文设计的数值格式的表现37-453.1.3只求法向流速的算法之探讨45-493.2理想盐度平流试验49-553.3黄浦江水动力验证55-59第四章模式在长江口的率定验证59-874.1长江口模式网格和基本设置59-624.2模式率定62-784.3模式验证78-864.4长江口试验小结86-87第五章总结87-89附录89-93附录1.文中常用符号一览表89-92附录2.攻读博士学位期间发表的学术论文92-93参考文献93-97论文的特色和创新点97-98论文不足之处及研究展望98-99致谢99