沙尘大气两相流数值模式的理想试验.pdfVIP

湿沉降（干沉降比如植被截留）和化学反应等过程，我们基于两相流（至少一相 为流体的两相物质所组成的流动系统）和质量变化的观点推导出了沙尘大气物理 约束方程组[26] 。理论分析揭示了许多有新意的结论，比如沙尘倾向于加强对流、 加强夹卷、改变沙尘云内部的细微结构等。但是理论分析是基于许多假设基础之 上得到的，实际情况可能有一定的特殊性。本文从数值模拟的角度对此问题进行 了探讨。 1. 沙尘大气两相流数值模式 1.1 沙尘大气两相流方程组 我们推导出的沙尘大气两相流方程组[26]在不计与分子扩散、辐射和化反应学 等有关的项时的形式如下： ∂ρ +∇⋅ρV ρS (1a) ∂t dV 1 =×Ω− − − ∇ (1b) 2V gk SV p dt ρ p ρRT (1c) 2 d c T ⎛ ⎞ ( ) p V v =− ∇⋅ − − + (1d) V S ⎜c T ⎟ Q dt ρ ⎜⎝ v 2 ⎟⎠ ∂qs V q q ( S S) (1e) + ⋅∇ =− s s s ∂t 其中，ρ 、V 、p 、T 和qs 是两相流密度、速度、气压、温度和沙尘组分的质量 R 浓度，cv 、 、Sa 、Ss 和S 是质量定容热容、‘气体常数’以及空气组分、沙尘 组分和混合系统的质量变化率。 一般而言，沙尘暴是大中尺度天气过程，但是也经常发生局地性的沙尘暴， 比如 2010 年 5 月 14 日，青海省格尔木地区出现了局地沙尘暴过程；2010 年 9 月 5 日下午新疆喀什出现局地沙尘暴，能见度不足 100m；2010 年 10 月 9 日， 内蒙古额济纳旗居延海地区出现局地沙尘暴，最小能见度不足百米，瞬时最大风 力在 7—8 级。张智等[27]就讨论了宁夏盐池春季局地沙尘暴特征及其与下垫面的 关系。而且，沙尘天气还包含扬沙天气，比如 2011 年 4 月 19 日20:00-22:00 左 右兰州发生局地扬沙。另外，热力不稳定是沙尘暴的三要素之一。为了把握沙尘 暴与扬沙天气内在的结构特征，采用小尺度的高分辨率的非静力平衡模式是合理 之选。 2 (1a)和(1c)代入(1d)式，注意 V T p 00 Rd / cp 1 、θ ( ) 以及忽略R 随时间的变化， 2c T p p 则得 2 dθ =− θ (2) S dt 引入气压函数