[发明专利]一种反演弱非线性普朗特模式中的涡流热导系数的方法

说明书

本发明公开了一种反演弱非线性普朗特模式中的涡流热导系数的方法，该方法包含以下步骤：1，给定正演模型以及风速和位势温度的观测数据；2，利用观测与模式的差值平方构建来构建代价泛函并得出泛函梯度的计算公式；3，根据经验给出初始猜测值；4，利用计算得出的泛函梯度以及最速下降法进行迭代求解；5，迭代退出的结果即为反演出的涡流热导系数。本方法利用风速和位势温度的观测，采用变分同化的方法来精确反演弱非线性普朗特模式中的涡流热导系数，代价泛函由观测与模式差值的平方构建。

技术领域

本发明属于大气边界层参数反演领域，公开了一种反演弱非线性普朗特模式中的涡流热导系数的方法。

背景技术

大气行星边界层的研究一直是大气科学研究中的热点问题，尤其是近些年的研究表明上下坡风能够影响冰川融化，海平面上升，气候变化等受到广泛关注的问题。但是，非常稳定分层大气边界层的理论和模型的研究还远没有完成。同时，随着数值天气预报模式和气候模式的发展，模式精度越来越高，斜坡地形上空的气流已经能被模拟，但是相应的下边界条件的处理和近地面湍流参数化都还是悬而未决的问题。所以，上下坡风仍然是有活力的研究课题之一。

在上下坡风的研究中，其解析模型的研究是最基础的工作。大多数下坡风的解析模型都是一维的，这些研究热力驱动的简单斜坡流的解析模型中最为经典的是普朗特模式。经典普朗特模式的优缺点已经在Grisogono和Axelsen的文中讨论，主要是在大涡度模拟和有限的观测数据集方面。这些缺点在Mo的文中得到了改进。但是还存在一个很大的缺陷没有改进，大多数下坡风的解析模型都是线性的，并且其涡流扩散系数和导热系数假定为常值。第一，线性模型带来的最大的问题就是无法将近地面强的流诱导的有限振幅的位势温度梯度反馈到环境位势温度梯度的特性体现在热力学方程中。

众所周知，参数反演问题通常是高度非线性的不适定问题，所以参数反演问题难以解决，还有许多理论问题等待解决。然而现实生活中有许多方面需要用到参数反演问题，例如地球科学，生物科学，材料科学，还有大气科学、海洋科学。涡流扩散系数的反演问题就是典型的参数反演问题。涡流扩散系数无论在流体力学还是大气科学中都是一个非常重要的参数，渐变涡流扩散系数K(z)通常采用的是经验参数化方案得出。而实际上涡流扩散系数是一个很难确定的物理参数，它与流体的动力热力结构密切相关。因此给出弱非线性普朗特模式中的涡流扩散系数合理的计算方法无论对于理论的发展还是数值模式的改进都是非常迫切的。

发明内容

发明目的：针对上述现有技术，提出一种反演弱非线性普朗特模式中的涡流热导系数的方法。

技术方案：一种反演弱非线性普朗特模式中的涡流热导系数的方法，包括以下步骤：

步骤一：给定正演模型以及风速和位势温度的观测数据：

(2)正演模型为弱非线性普朗特模式：

其中：g是重力加速度；θ是位势温度与背景场位势温度Θ的偏差；Θ₀是参考温度；α是坡角，α<0代表的是下坡；K(z)和Pr分别是涡流热导系数和湍流普朗特常数；u是风速；z是垂直于斜坡的坐标；z^*代表的是垂直于水平面的方向；ε为常数参数，0≤ε≤0.01；

涡流热导系数K(z)的模型采用其中K_min是K(z)的最小值；K₀是K(z)的最大值；h是最大的K值出现的高度；

步骤二：利用观测与模式的差值平方来构建代价泛函并得出泛函梯度的计算公式：

(3)代价泛函J[K₀,h]为：