任意曲线坐标系下二维潮流数学模型.pdfVIP

任意曲线坐标系下二维潮流数学模型 澍军 (上海航道勘察设计研究院，上海市200120) 摘要 奉文应用肯限差分法，柱任意曲线坐标下对长江口地区进行了_二绺潮流数学模型计算．计算表明： 该方法嘲格划分任意性火．能较好地拟台边界．边界条件处理简单。计算教率高，计算值与实测值吻台较 好．该方法可在实际[程中加以运用。 关键词数学模型 涮流 曲线坐标 差分法 I前言 海岸、河口地区水域宽广，水深相对较浅 水流运动适宜采用二维模型。只在有特殊 要求对局部水域场合下．才需进行三维模拟。 因此二维数值模拟仍是研究海岸、河口地区 潮汐水流运动的主耍方法之一。 目前海岸、河口地区数值模拟计算经常采用的方法有有限元法、边界单元法与有限差 分法等。有限元法与边界单元法一般适用于方程简单、边界复杂情况，矩形喇格有限差分 标变换技术，将不规则的物理区域变换为正交的规则均匀的计算区域，在计算区域上求解 椭圆方程，得出物理平面上的正交网格位置，将相应方程变换到计算区域上．在计算区域 上进行求解。该方法为解决直角坐标系系下有限差分法边界吻台差的雉点开辟了新的研究 途径。但这一方法对边界较不规则且多岛屿的海岸、河口地区要做到既能较好地拟台边界． 又要保证物理域内网格处处正交，则正交变换技术比较繁杂。 本文采用任意曲线坐标，物理域内网格无须正交，流速分量沿曲线法向，建立了二维 潮流数学模型，对长江口地区潮流运动进行了模拟计算。 2 曲线坐标系下二维潮流数学模型 2．1 曲线坐标系下二维潮流运动基本方程及定解条件 直角坐标系F--：维潮流运动基本方程为： 堂+旦幽+曼幽：o Ot 8x 却 (2) 氨 av 。 。良 塑Ot+u警+y詈一∥+g篆+型垃C2H=。 33S一 詈+u豢+y罾+∥+g芳+警一oo， 其中f为潮位。U、矿为j、Y向速度分量，Ⅳ为总水深，c为谢才系数·科氏系数 f=2口sin驴(脚为地球地转角速度，妒为计算区域平均地理纬度)· 方程(1)～(3)定解条件为： 初始条件： fG，Y，≠】，。=厶G，Y) (4) U(x，Y，fl。=0 (5) 矿G，Y，，】。=0 (6) 边界条件： 在水边界上：“x，Y，f)=厶b，Y，t) (7) 在固边界上： V。=0 (8) y．为圃边界法向流速分量。 每 “=∽【，一却矿)／石 (9) v=ky—yfU)／痧 (10) 在叩、f坐标法向分量a 将(9)、(10)式代入方程(1)一(6)并利用如下关系式： 篆=嚣六+嚣仉=专(嚣蜘一茜，0 ㈣， 面。面￡+瓦巩2了l面即一琵斯j爹=嚣￡+鸶巩=；(荔即一嚣x0 nz，““