积云数值模拟教案资料.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,积云数值模拟,微物理过程参数化,基本概念：,云雨粒子分类：云滴、雨滴、冰晶、霰、雪、雹。,微物理过程参数化,比含量定义,式中，,m,J,为某种粒子的质量；,m,a,空气的质量。类似比湿的定义,优点：在模式计算中，避免了空气密度的计算。,微物理过程参数化,守恒方程,云中粒子的定义,暖云：,云滴（D,100,）；雨滴（D,100,）,冷云：,云滴（D,100,）；雨滴（D,100,）,冰晶（D,100,）；霰、雪、雹（D,100,）,粒子间相互转化的物理过程,水汽,云水,雨水,降水,凝结,蒸发,自动转化,碰并,蒸发,冰晶,活化,凝华,升华,雪霰雹,自动转化,碰并,升华,淞附,冻结,融化,淞附,碰并,冻结,融化,粒子谱分布,云滴：,c,2,雨滴：,R,0,雹、霰：,g,0,冰晶：,i,2,若雨水比含量为q,r,单参数法,：给出,N,0R,，由,q,R,的含量决定谱型，简单，有误差。,双参数法,：,N,0R,不给定，给出总浓度,N,的预报方程，由,N,0R,和q,R,共同决定谱型。,粒子下落末速度,单个雨滴,单个霰、雹,群滴下落末速度,质量加权末速度；M,DJ,为直径D的单个的J类粒子质量；N,J,为谱分布,凝结、蒸发过程,1、云滴：,单滴凝结公式谱积分求解的缺陷：,水汽含量下降很快，潜热释发增温，导致不饱和。,解决方法：,平衡处理法：,平衡处理法,在凝结、蒸发计算前，温度为T*比湿为q*；凝结水量为,后，达到平衡。,有关系式,Teten,公式,可得：,平衡处理法,设,L/(c,p,),是小量,，,对,q,vs,(T*+L/(c,p,),做泰勒展开，取一阶近似有：,精确解可以用迭代法求解,平衡处理法,为空气过饱和时通过凝结过程使得空气中水汽达到饱和时所凝结出的水汽量,当,0,，空气不饱和，,为通过蒸发过程使得空气达到饱和时所需要蒸发的液水量,当云水,q,c,200,q,c,10-5,霰、雹结凇增长,干生长,：雹块处于较冷，过冷水滴较少的云中，冻结潜热很快散失，形成雹块的不透明层。增长率由碰併水量来决定。（不考虑热量平衡问题）,同理可以计算雹对冰晶的碰并率,对与对于霰碰併雨滴，应考到不同雨滴大小霰粒子之间的相对运动,。,霰、雹结凇增长,湿增长,：雹块处于较暧且液水丰富的云中，有一部分收集到的液水不会冻结而保持液体状态，在雹块中形成透明层，此时要考虑热量平衡。,假定,：在一个格点上，只存在一种增长方式，比较干湿增长率的大小，选取取值小的一种增长方式为当前增长方式；,在暖区,，还应考虑正在融化的冰雹对云水的碰併过程并且脱落成为雨水的过程,设直径为D的雨滴在直径为D,i,，浓度为N,i,的冰晶空间扫过，单位时间内捕获一个冰晶的概率为：,可得：,过冷雨滴碰併冰晶形成冻滴转为霰,自动转化过程,云滴自然长大到雨滴尺度所形成雨水的过程,同理，冰晶长大到霰和雹,Kessler,云雨自动转换公式,经验公式,：,A,cr,=K(q,c,-a),q,c,a,K0;q,c,b,K0;,q,i,b,K=0,阈值,b,根据结果进行调试,转换率,K=10,-3,exp(,(T-T,0,),)，与温度有关,核化、冻结与融化,）冰晶的核化:,自然界中活化冰模浓度与温度的关系（经验、观测）给出,冰晶的初始质量m,i0,=10,-12,kg,相当于D,16.3,m,）雨滴的冻结,Bigg给出的：水滴自发形成冰相胚胎的机率与水滴的体积成正比,T0,时,冰晶瞬时全部融化成云水；,霰、雹有有热量平衡过程控制,对单个出，然后对全谱积分。,繁生与破碎,对冰相而言，有三类重要的繁生机制。,a.,冰晶与霰、雨及其它冰晶相碰发生断裂。（浓度的增加）,b.,冰粒子在凇附较大云滴时引起破碎冰层脱落。（质量转移）,c.,大过冷滴在冻结过程中发生破碎而产生次生冰晶。（质量转移）,雨滴破碎,，改变谱型和浓度，无种类间质量转移。,微物理过程的非参数方法,云微物理数值模拟最早开展（1959）的工作，粒子划分为许多档。计算各档粒子浓度的演变。,早期：拉格朗日模式，水滴跟随气块移动,近期：欧拉模式：,水滴大小分档,凝结核0.03,m；小云滴12,m，雨滴,13mm,在相同的尺度间隔内，粒子浓度差别最大可达10个量级，小粒子多，大粒子少,。,标准云雨滴尺度,Berry,指数分档,r,I,=r,0,exp(I-1)/I,0,x,I,=x,0,exp3(I-1)/I,0,r,0,x,0,第一档的半径和质量,I,0,由总档数或间隔大小来确定,r,0,考虑盐核到雨滴的生长过程,，,r,0,0.03,m,谢谢,已有的积云模式,气泡模式（0维模式）,一维非定常模式,一维半非定常模式,二维柱对称非定常模式,二维面对称非定常模式（流体力学模式）,三维模式完全弹性模式,存在的问题,微物理过程的处理,动力框架结构,数值计算方法,不同问题，用不同模式,此课件下载可自行编辑修改，仅供参考！感谢您的支持，我们努力做得更好！谢谢,