沙尘暴沉降的两相流理论（一）.doc

沙尘暴垂直输运的两相流理论（一） ----气块模式* 巢纪平 刘飞 （1）国家海洋环境预报中心, 北京 100081 （2）中国科学院大气物理研究所东亚区域气候-环境重点实验室， 北京， 100029； 中国科学院研究生院， 北京， 100049 摘要： 沙尘暴已经成为我国西北、华北地区严重的气候灾害，而对沙尘的研究主要集中在起沙问题与沙尘沉降问题，在理论上对前者研究较多，而对后者的研究则比较缺乏。假设空气块中带有沙尘，沙尘颗粒很小，可不计其末速，随空气一起运动，这样由空气和沙尘组成两相流。在两相流情况下，分析了均匀环境和非均匀环境条件沙尘沉降的条件和速率。通过分析得到在均匀背景场下有利于沙尘沉降的条件为大的初始沙尘浓度扰动与小的初始温度扰动。进一步考虑更为真实的环境场，模式给出三种特征沉降区域：1）速降区，沉降速率较快；2）缓降区，沙尘先上升后下降，沉降时间较长；3）波动区，沙尘在空气中不断的上下波动，此时很难判断沙尘何时沉降，取决于地面状况。同时也分析了沙尘不断上扬的背景条件。 关键词 沙尘暴 两相流 沉降 气块模式 引言 沙尘暴的频繁发生己成为我国西北、华北地区严重的气候灾害。关于沙尘暴的源地，发生的天气条件及输送路径等，己有不少研究。在曾庆存等著的“千里黄云－东亚沙尘暴研究” 一书中基本上总结了己有的研究成果并有著者们自身的发展[1]。因此与本文动力学研究无密切关系的文章在此不再引用，均可参阅该书。 * 通讯作者：巢纪平，主要从事热带大气海洋动力学，jipingchao@ 沙尘动力学至少有两个问题需要研究。一个是沙尘暴的生成即起沙问题，显然它涉及下垫面的物理状态，及其上面大气贴地层和边界层的湍流状态，以及向上输送沙尘的条件和过程，相对来讲这个问题的动力学比较复杂[2]。另一个是沙尘的沉降问题，对于该过程的理论分析还比较少。沙尘输送风力的强度，天气情况、起沙过程[3,4]，沙尘输送过程中干湿不同沉降过程[5],云物理混合过程[6]等过程决定了大气中存在不同的沙尘尺度分布。虽然Schütz研究表明沙尘粒径分布最后可能会到达一种平衡态[7]，但是还缺乏充分的试验来证明。最近研究表明不同沙尘粒径的分布在海洋上存在着不同干沉降的关系[8]。但是为了抓住沉降的物理本质，如不考虑由于沙尘而改变的空气中辐射转输的光学性质，只考虑动力学方面，则这个问题比较简单，比如当沙尘颗粒比较小时，它的末速不大，沙尘可以和所在的空气团一起运动，形成了一个两相流的气团，从而改变了气团的密度及相应的阿基米德浮力，使空气团的垂直加速度相应发生改变，也改变了沉降的条件，足够大浓度的沙尘如能沉降到地面即出现沙尘暴。由于是将沙尘大气看作为两相流，本文将不同物理化学性质的沙尘都用初始沙尘浓度来统一区别，比如较大粒径的沙尘在本文中就通过较大的初始沙尘扰动来描述，这有利于把握沙尘沉降的物理本质。 本文从两相流的角度来研究沙尘的沉降和垂直输运，并先从简单的气块法入手，后续工作将进一步发展到气拄挟卷模式，最后给出三维非线性模式的研究。沙尘质量和湍流过程在沙尘沉降过程中有一定的作用，而气块法的局限性在于将气块看成质点，忽略其与环境的交换，所以湍流的过程在本文中不能描述，将在气柱挟卷模式中简单引入，在三维非线性模式中详细讨论。 该简化的气块模型能够给出沙尘沉降的解析解，能够很好的用来讨论沙尘沉降条件与过程，但是正因为是简化的理论模型，为了求得解析解，对真实的物理过程做了极大的简化，所以与复杂的沙尘输送业务模式相比，该理论模型所反映的物理过程比较简单，还有许多不足。但是该理论模式是能够从物理本质上讨论沙尘沉降的一个有力的工具，可以作为发展复杂的沙尘输送业务模式系统的一个判别与检测工具。 两相流的密度变化 设空气和沙尘的密度分别为，混合体的体积为，其中空气占有的体积为，沙尘占有的体积为，总质量和分部质量分别为，并有 (1) 进而有 (2) 式中定义成沙尘的浓度。 密度变化值为 (3) 设空气的热膨胀系数为，对理想气体有，有，为空气温度，上式改写成 (4) 式中，为环境的气候值或平均值。可见混合体的密度不只是决定于空气的温度，也决定于沙尘的浓度。 均匀环境中的气块运动 3．1 基本模型建立 设气块外环境的平均温度和沙尘浓度分别为，气柱和环境之间存在扩散过程，导致气柱内物理量变化，这一过程是线性的，为 (5) (6)