爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 22268|回复: 22

[求助] 理查森数的计算

[复制链接]

新浪微博达人勋

发表于 2013-1-10 16:52:43 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x

有大神计算过梯度理查森数,以及湍流交换系数吗?

我根据梯度理查森数的计算公式,利用上下两气层的温度差、风速差求解,但是由于某些层的风速差值太小,甚至是相等的,导致理查森数太大,虽然从计算公式上看不出有什么问题,但是理论上感觉不合理,且对于较大的理查森数(有的甚至是不存在的值)不适合后续的计算湍流交换系数!对于这个问题该怎么处理呢?

附件贴了理查森数和湍流交换系数的公式
360截图20130110164312245.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-1-21 15:03:39 | 显示全部楼层
本帖最后由 爱的侍者 于 2013-1-21 16:14 编辑

梅子 ~ 你这个Ri定义是哪本书中规定的,我找的定义不是这个呀~!l0是混合长还是特征长度?
探空数据是三维的嘛??
如果取x轴与所在噶偶平均风向一致,那么分母只选取前一项,只需知道另个高度见的位温与风速梯度即可,位温如何计算的?找了一下没找到,根据位温的定义,气体从原有的压强与温度出发,绝热膨胀或压缩到标准压强时(1000mb)的温度,可用来比较不同气压下的气体热状态。你公式中的rd为1吗?另外你的Km和Kh的定义在哪里规定的,我查到的是要根据三个方向风的数据算出摩擦速度,结合特征位温,计算出混合长度,Km的计算需要用到摩擦速度,但是没见过你给的公式,提供出处学习一下~!~谢啦

QQ截图20130121144629.png

密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2013-1-12 15:55:07 | 显示全部楼层
上下两层的风速差值较小,甚至出现上层风速比下层风速小的状况,实际测风资料就会出现这个状况。
你用那两层算的10 70??
有没有程序没,共享学习下嘛~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-1-12 15:55:29 | 显示全部楼层
上下两层的风速差值较小,甚至出现上层风速比下层风速小的状况,实际测风资料就会出现这个状况。
你用那两层算的10 70??
有没有程序没,共享学习下嘛~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-1-12 18:31:33 | 显示全部楼层

10和70 是什么东西,程序倒是有,是我自己编写的,也不知道对不对,我可以从给我整个程序里面截一段给你看看,顺便帮我看一下写得是否正确,谢谢!
i=1
   Ri(i)=g/temp1(1)*((temp1(i)-Tg)/dz+rd)/((u1(i)-0)/dz)**2
   do i=2,ny
      Ri(i)=2*g/(temp1(i)+temp1(i-1))*((temp1(i)-temp1(i-1))/dz+rd)/((u1(i)-u1(i-1))/dz)**2
   enddo
    i=1
   if(Ri(i)<0)then
      km(i)=lo**2*abs((u1(i)-0)/z(i))*(1-16*Ri(i))**0.5   !动力交换系数单位:m^2/s
      kh(i)=(1-16*Ri(i))**0.25*km(i)
   elseif(Ri(i)>0)then
      km(i)=lo**2*abs((u1(i)-0)/z(i))*(1-5*Ri(i))**2
      kh(i)=km(i)
   endif
   do i=2,ny
      if(Ri(i)<0)then
        km(i)=lo**2*abs((u1(i)-u1(i-1))/dz)*(1-16*Ri(i))**0.5
            kh(i)=(1-16*Ri(i))**0.25*km(i)
      elseif(Ri(i)>0)then
        km(i)=lo**2*abs((u1(i)-u1(i-1))/dz)*(1-5*Ri(i))**2
            kh(i)=km(i)
      endif
   enddo
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-1-13 22:05:55 | 显示全部楼层
一颗梅子在南京 发表于 2013-1-12 18:31
10和70 是什么东西,程序倒是有,是我自己编写的,也不知道对不对,我可以从给我整个程序里面截一段给你看 ...

你是用啥数据算的~ 一般这种近地层边界层参数,都是用铁塔数据 。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-1-14 08:53:02 | 显示全部楼层
爱的侍者 发表于 2013-1-13 22:05
你是用啥数据算的~ 一般这种近地层边界层参数,都是用铁塔数据 。

我用的是自己放球的数据!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-1-14 09:09:39 | 显示全部楼层
一颗梅子在南京 发表于 2013-1-14 08:53
我用的是自己放球的数据!

o 哦 知道啦~谢谢回复!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-1-14 09:46:00 | 显示全部楼层
爱的侍者 发表于 2013-1-14 09:09
o 哦 知道啦~谢谢回复!

大家一起讨论问题不用谢的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-1-14 09:47:57 | 显示全部楼层
爱的侍者 发表于 2013-1-14 09:09
o 哦 知道啦~谢谢回复!

大家一起讨论问题,不用谢的!有什么有用资料,大家一起共享啊!
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表