大气边界层的数值模拟方法研究：修正的壁面函数.docxVIP

大气边界层的数值模拟方法研究：修正的壁面函数

一、引言

大气边界层作为大气圈与下垫面相互作用的重要区域，其流动特性的准确模拟对于天气预报、空气污染扩散研究以及风能利用等领域都具有至关重要的意义。在大气边界层的数值模拟中，壁面函数法因其计算效率高、对网格要求相对较低等优点而被广泛应用。然而，传统的壁面函数在处理复杂下垫面条件和非平衡流动时，往往存在一定的局限性，导致模拟结果与实际情况存在偏差。因此，对壁面函数进行修正，以提高大气边界层数值模拟的精度，成为当前研究的热点之一。

二、大气边界层的基本特性

（一）大气边界层的结构

大气边界层通常可以分为近地面层和埃克曼层。近地面层是大气边界层的最底层，其厚度一般为几十米到几百米，该层内的流动主要受地面摩擦力的影响，风速随高度的变化遵循对数律。埃克曼层位于近地面层之上，其厚度通常为1-2公里，该层内的流动不仅受地面摩擦力的影响，还受到科里奥利力的作用，风速和风向随高度会发生明显的变化。

（二）大气边界层的流动特性

大气边界层的流动具有强烈的湍流特性，湍流运动使得动量、热量和水汽等物理量在垂直方向上进行强烈的交换。此外，大气边界层的流动还受到下垫面粗糙度、地形以及大气稳定度等因素的影响，使得流动特性变得更加复杂。

三、传统壁面函数的基本原理及局限性

（一）传统壁面函数的基本原理

传统壁面函数基于半经验公式，将壁面附近的湍流流动与主流区的流动联系起来。其基本思想是在壁面附近的粘性底层和对数律层之间建立一种关系，通过求解壁面切应力和湍流动能等参数，来模拟壁面附近的流动特性。

（二）传统壁面函数的局限性

对复杂下垫面的适应性差：传统壁面函数通常假设下垫面是均匀的、光滑的，而实际大气边界层中的下垫面往往是复杂多样的，如城市建筑群、森林等，这些复杂下垫面会导致壁面附近的流动特性发生明显变化，传统壁面函数难以准确描述这种变化。

处理非平衡流动的能力不足：在大气边界层中，由于受到天气系统、地形等因素的影响，流动往往处于非平衡状态，如强风条件下的流动、地形引起的绕流等。传统壁面函数基于平衡流动假设，在处理非平衡流动时，容易出现模拟误差。

对大气稳定度的敏感性：大气稳定度是影响大气边界层流动特性的重要因素之一。传统壁面函数在处理不同稳定度条件下的流动时，往往不能准确反映湍流交换系数的变化，从而导致模拟结果与实际情况存在偏差。

四、修正的壁面函数方法

（一）基于下垫面粗糙度的修正

为了提高壁面函数对复杂下垫面的适应性，研究人员提出了基于下垫面粗糙度的修正方法。该方法通过引入下垫面粗糙度参数，来调整壁面附近的湍流交换系数和风速廓线。例如，在城市大气边界层模拟中，可以将城市建筑群的粗糙度等效为一个粗糙度长度，然后通过修正的壁面函数来模拟城市下垫面附近的流动特性。

（二）考虑非平衡流动的修正

针对传统壁面函数处理非平衡流动能力不足的问题，研究人员提出了考虑非平衡流动的修正方法。该方法通过引入非平衡参数，来调整壁面附近的湍流动能和耗散率等参数。例如，在强风条件下的流动模拟中，可以通过修正的壁面函数来考虑流动的非平衡性，从而提高模拟结果的准确性。

（三）基于大气稳定度的修正

为了准确反映大气稳定度对壁面附近流动特性的影响，研究人员提出了基于大气稳定度的修正方法。该方法通过引入大气稳定度参数，来调整壁面附近的湍流交换系数和风速廓线。例如，在稳定大气条件下，湍流交换系数较小，风速廓线呈现出较为平缓的变化；而在不稳定大气条件下，湍流交换系数较大，风速廓线呈现出较为陡峭的变化。通过修正的壁面函数，可以更好地模拟不同稳定度条件下的大气边界层流动特性。

五、修正的壁面函数在大气边界层数值模拟中的应用

（一）城市大气边界层模拟

将修正的壁面函数应用于城市大气边界层模拟中，可以更好地模拟城市建筑群对大气流动的影响。例如，在模拟城市空气污染扩散时，修正的壁面函数可以更准确地描述城市下垫面附近的风速和湍流特性，从而提高空气污染扩散模拟的精度。

（二）复杂地形大气边界层模拟

在复杂地形大气边界层模拟中，修正的壁面函数可以更好地处理地形对大气流动的影响。例如，在模拟山区风场时，修正的壁面函数可以考虑地形引起的绕流和湍流增强等现象，从而提高风场模拟的准确性。

（三）不同大气稳定度条件下的模拟

将修正的壁面函数应用于不同大气稳定度条件下的大气边界层模拟中，可以更好地反映大气稳定度对流动特性的影响。例如，在模拟夜间稳定大气边界层时，修正的壁面函数可以更准确地描述湍流交换系数的减小和风速廓线的变化，从而提高模拟结果的可靠性。

六、结论与展望

（一）结论

传统壁面函数在处理复杂下垫面条件和非平衡流动时存在局限性，对大气稳定度的敏感性也较高。

修正的壁面函数通过考虑下垫面粗糙度、非平衡流动和大气稳定度等因素，提高了大气边界层数值模拟的精度。

修正