上海海洋大学《海洋要素计算》课件-第11课时 潮流调和常数的计算.pptxVIP

欧拉余流斯托克斯漂流拉格朗日余流其中，时间平均算子欧拉余流和拉格朗日余流的表达式

渤海表层欧拉余流（A），斯托克斯漂流（B）和拉格朗日余流（C）的数值模拟结果资料来源：JMS,2004,44:141-151.ABC

潮流也可以表示为许多分潮流之和（2）（3）其中U0,V0为余流，Ui,ξi为北分流的调和常数，Vi,ηi为东分流的调和常数单纯进行潮流预报，公式（2）足够了，但是调和常数不能直观揭示该地点潮流变化的矢量特征，因此需要引入其它参数。考虑一个分潮的情况其中

潮流椭圆：潮流分量u，v的矢量端划出的轨迹。椭圆长半轴和短半轴就是这个分潮流速可能达到的最大和最小值，分别叫做最大和最小潮流，记作W，w旋转率：k=-w/W，逆时针为正分潮流的椭圆要素：最大分潮流流速W、方向θ、发生的时间τ以及旋转率k决定了分潮流椭圆的基本特征，叫做~

将uv坐标轴旋转一个θ角(最大潮流发生方向)，使其与椭圆的长短轴重合，记为xy坐标轴，则uv在xy上的投影可写为

潮流调和常数和椭圆要素的换算

O（4）

又因为xy轴分别与椭圆的长轴和短轴重合，因此椭圆的参数方程可写为（5）o式中上式y的等号右侧取负号表示旋转率顺时针方向为负

比较（4）和（5）（6a）（5）（6d）（6c）（6b）（4）可以得到

将公式（6a~6d）做如下处理（7a）（7d）（7c）（7b）由潮流椭圆要素计算调和常数的方法求解上面的四元一次方程组即可得到由椭圆要素计算调和常数的公式

计算椭圆要素的Werenski?ld算法（16）（6a）（6b）（6c）（6d）（17）由公式（6a）+（6d）和（6b）-（6c）得到上两式两边分别求平方，相加之后再开方得到