《积云数值模拟》由会员分享，可在线阅读，更多相关《积云数值模拟（39页珍藏版）》请在金锄头文库上搜索。
1、积云数值模拟,一、数值模式中云的微物理过程,动力学方程组中与微物理过程有关的项有： 1）运动方程中的拖带力。 2）热量方程中相变过程的热力作用。 3）水汽方程中相变过程改变了水汽的含量。,微物理过程处理的方法,1）参数化方法：考虑云中粒子的总体特征2）详细描述方法：各种粒子分档,微物理过程参数化,基本概念： 云雨粒子分类：云滴、雨滴、冰晶、霰、雪、雹。,微物理过程参数化,比含量定义,式中，mJ 为某种粒子的质量；ma 空气的质量。类似比湿的定义 优点：在模式计算中，避免了空气密度的计算。,微物理过程参数化,守恒方程,云中粒子的定义,暖云： 云滴（D100）；雨滴（D100） 冷云： 云滴（ D
2、100）；雨滴（D100） 冰晶（ D100）；霰、雪、雹（D100）,粒子间相互转化的物理过程,水汽,云水,雨水,降水,凝结,蒸发,自动转化,碰并,蒸发,冰晶,活化,凝华,升华,雪霰雹,自动转化,碰并,升华,淞附,冻结,融化,淞附,碰并,冻结,融化,粒子谱分布,云滴：c2 雨滴：R0 雹、霰：g0 冰晶：i2,若雨水比含量为qr,单参数法：给出N0R，由qR的含量决定谱型，简单，有误差。 双参数法： N0R不给定，给出总浓度N的预报方程，由N0R和qR共同决定谱型。,粒子下落末速度,单个雨滴,单个霰、雹,群滴下落末速度,质量加权末速度；MDJ为直径D的单个的J类粒子质量；NJ为谱分布,凝结、
3、蒸发过程,1、云滴： 单滴凝结公式谱积分求解的缺陷：水汽含量下降很快，潜热释发增温，导致不饱和。 解决方法：平衡处理法：,平衡处理法,在凝结、蒸发计算前，温度为T* 比湿为q*；凝结水量为后，达到平衡。,有关系式,Teten公式 可得：,平衡处理法,设L/(cp)是小量，对qvs(T*+ L/(cp)做泰勒展开，取一阶近似有：,精确解可以用迭代法求解,平衡处理法,为空气过饱和时通过凝结过程使得空气中水汽达到饱和时所凝结出的水汽量 当0，空气不饱和， 为通过蒸发过程使得空气达到饱和时所需要蒸发的液水量 当云水qc 200 qc10-5,霰、雹结凇增长,干生长：雹块处于较冷，过冷水滴较少的云中，冻
4、结潜热很快散失，形成雹块的不透明层。增长率由碰併水量来决定。（不考虑热量平衡问题）同理可以计算雹对冰晶的碰并率 对与对于霰碰併雨滴，应考到不同雨滴大小霰粒子之间的相对运动 。,霰、雹结凇增长,湿增长：雹块处于较暧且液水丰富的云中，有一部分收集到的液水不会冻结而保持液体状态，在雹块中形成透明层，此时要考虑热量平衡。 假定：在一个格点上，只存在一种增长方式，比较干湿增长率的大小，选取取值小的一种增长方式为当前增长方式 ； 在暖区，还应考虑正在融化的冰雹对云水的碰併过程并且脱落成为雨水的过程,设直径为D的雨滴在直径为Di，浓度为Ni的冰晶空间扫过，单位时间内捕获一个冰晶的概率为：可得：,过冷雨滴碰併
5、冰晶 形成冻滴转为霰,自动转化过程,云滴自然长大到雨滴尺度所形成雨水的过程 同理，冰晶长大到霰和雹,Kessler云雨自动转换公式,经验公式：Acr=K(qc-a)qca, K0; qcb, K0; qib, K=0阈值b根据结果进行调试 转换率 K=10-3exp(T-T0)，与温度有关,核化、冻结与融化,）冰晶的核化 : 自然界中活化冰模浓度与温度的关系（经验、观测）给出 冰晶的初始质量mi0=10-12kg, 相当于D16.3m ）雨滴的冻结 Bigg给出的：水滴自发形成冰相胚胎的机率与水滴的体积成正比 T0时 冰晶瞬时全部融化成云水；霰、雹有有热量平衡过程控制 对单个出，然后对全谱积分
6、。,繁生与破碎,对冰相而言，有三类重要的繁生机制。 a. 冰晶与霰、雨及其它冰晶相碰发生断裂。（浓度的增加） b. 冰粒子在凇附较大云滴时引起破碎冰层脱落。（质量转移） c. 大过冷滴在冻结过程中发生破碎而产生次生冰晶。（质量转移） 雨滴破碎，改变谱型和浓度，无种类间质量转移。,微物理过程的非参数方法,云微物理数值模拟最早开展（1959）的工作，粒子划分为许多档。计算各档粒子浓度的演变。 早期：拉格朗日模式，水滴跟随气块移动 近期：欧拉模式：,水滴大小分档,凝结核0.03m；小云滴12m，雨滴13mm 在相同的尺度间隔内，粒子浓度差别最大可达10个量级，小粒子多，大粒子少 。,标准云雨滴尺度,Berry指数分档,rI=r0exp(I-1)/I0 xI=x0exp3(I-1)/I0 r0, x0第一档的半径和质量 I0由总档数或间隔大小来确定 r0 考虑盐核到雨滴的生长过程，r00.03m,谢谢,已有的积云模式,气泡模式（0维模式） 一维非定常模式 一维半非定常模式 二维柱对称非定常模式 二维面对称非定常模式（流体力学模式） 三维模式完全弹性模式,存在的问题,微物理过程的处理 动力框架结构 数值计算方法 不同问题，用不同模式,
