WRF耦合风力机参数化问题

dueler

搬砖也要学气象 发表于 2020-9-15 11:58
我理解的是，打开这一项是为了对应MYNN方案，我们用QKE表现2 X TKE，这样可以更明显的展现风机产生的TKE。 ...

不打开这一项，WRF输出的结果也是QKE呀。文献中说的打开这一项与不开这一项的结果对比，TKE的数值变小了。我目前的理解是打开这一项可能使风力机后本来很剧烈的掺混被抹平了不少，因此输出TKE就减小了。

搬砖也要学气象

我理解的不是这样，我认为你用的是MYNN而不是MYJ，这一项本身就应该打开，来配合QKE输出，我不知道你读的是哪篇文献，我不懂了。我按照我的理解打开了，效果还好。另外，你也可以调一下你的积分步长再试试。你的积分步长我个人感觉不太科学。还有你的eta层在1km以内是否还可以加密些呢。以上仅供参考。愿意的话，酌情试试。如果你按照你的思路彻底解决了，可以再交流。

dueler

搬砖也要学气象 发表于 2020-9-15 20:41
我理解的不是这样，我认为你用的是MYNN而不是MYJ，这一项本身就应该打开，来配合QKE输出，我不知道你读的是 ...

我觉得大气稳定度的影响挺大的，换了一个大气稳定的区域尾流就比较明显了。如果按照那篇文献中说的，不打开尾流就会更明显，当然那篇文章也是在稳定大气条件下模拟的。

搬砖也要学气象

dueler 发表于 2020-9-15 21:46
我觉得大气稳定度的影响挺大的，换了一个大气稳定的区域尾流就比较明显了。如果按照那篇文献中说的，不打 ...

老铁，咱现在用的这个不考虑尾流，只是中尺度风机模型，尾流有专门的fitch尾流模块。那个是计算流体力学了，你看的文章和气象学区别较大。这个我很确定。另外，我的看法的依据是WRF手册，对bl_mynn_tkeadvect是这么说的，(default) off; does not advect tke in MYNN scheme (default)。我觉得可以理解为QKE的开关。你不妨试试，也许就可以了。看你自己

dueler

搬砖也要学气象 发表于 2020-9-15 22:21
老铁，咱现在用的这个不考虑尾流，只是中尺度风机模型，尾流有专门的fitch尾流模块。那个是计算流体力学 ...

中尺度风机模型不就是指的Fitch参数化模型么？bl_mynn_tkeadvect这一项是在WRF3.8中新加的，之前WRF的风力机模拟只需要打开windfarm_opt就可以啦

搬砖也要学气象

dueler 发表于 2020-9-16 17:20
中尺度风机模型不就是指的Fitch参数化模型么？bl_mynn_tkeadvect这一项是在WRF3.8中新加的，之前WRF的风 ...

中尺度风机模型当然是FITCH参数化模型，但不是你说的尾流，你觉得中尺度模型还能考虑尾流吗？那个是计算流体力学了。fitch有中尺度参数化模型，也有尾流模型。这个得清楚。之前是之前，你用的新版，你就该打开，而且你的步长设置，也不太科学，最好再看看手册和更多文献。我就建议到这里，兄弟你自己掌握。

dueler

搬砖也要学气象 发表于 2020-9-17 10:35
中尺度风机模型当然是FITCH参数化模型，但不是你说的尾流，你觉得中尺度模型还能考虑尾流吗？那个是计算 ...

好的，多谢

xinshou6

他可能看的使用 WRF-Fitch 对湖区风电场风力发电机尾流效应特征的数值模拟，Fitch 等[19–20]将一个风场参数化模块——Fitch 尾流模型(Fitch Wake Model)添加到WRF 中。将该尾流模型与MYNN 方案结合, 可实现在中尺度天气模型中在线耦合模拟风电场风场特征与尾流效应。实质上, Fitch 模型是将风机涡轮的拖曳系数引入WRF 中的风速和湍流动能的计算中, 将风场解析为大气动量的汇和湍流动能产生的源[21]在 WRF 中, 借助趋势量(tendency)的传递,Fitch 模型对风速、动量和湍流动能等要素进行数据流入和流出的双向计算。建议 Fitch 模型的适用分辨率为低于5 倍扇叶直径。以上摘自它。字面看是加完了的------我不是大佬 我也不会 有空看看怎么算----每次编辑都回退 New in Version 3.3, and in this version it only works with 2.5 level MYNN PBL option (bl_pbl_physics=5), and updated in V3.6。这个是只有5才能用吗 不是说更新了都能用吗？前面说边界层方案要改，QKE开关要打开是什么意思？ windfarm_opt不用全开吗？位置确定在范围里面？用的是默认的参数吗 我觉得一般默认就应该有效果了吧。。要不把参数功率什么的调大看看？du(i,k,j) = du(i,k,j)-.5*u(i,k,j)*thrcof*speed*tarea/dz(i,k,j)*ec 要算的是这个不？Fitch 模型基于风机涡轮推力系数计算风机从大气中提取的动能(kinetic energy, KE)及风速的函数湍流动能 TKE (turbulence kinetic energy)。因此, Fitch 适用的边界层参数化是能够产生湍流动能项输出的 MYNN 方案。由于在高于 5 倍扇叶直径(70 m 左右)的空间分辨率下, 必须再考虑复杂的扇叶旋转, 因此, 为了简化,Porte-Agel 等[22]建议 Fitch 模型的适用分辨率为低于5 倍扇叶直径。------10.11 Registry.EM_COMMON:package   fitchscheme    windfarm_opt==1 ,./phys/module_pbl_driver.F:  windfarm_select: SELECT CASE(windfarm_opt)，CASE (fitchscheme)，  IF (PRESENT(id) .AND.     PRESENT(z_at_w) ) THEN，     CALL wrf_debug(100,'in phys/module_wind_fitch.F'),CALL dragforce,module_wind_fitch.F:Represents kinetic energy extracted by wind turbines and turbulence (TKE) they produce at model levels within the rotor area.   CALL dragforce(id & &,z_at_w,u,v & &,dx,dz,dt,qke & &,du,dv & &,windfarm_opt,power & &,ids,ide,jds,jde,kds,kde & &,ims,ime,jms,jme,kms,kme & &,its,ite,jts,jte,kts,kte)，du(i,k,j) = du(i,k,j)-.5*u(i,k,j)*thrcof*speed*tarea/dz(i,k,j)*ec 是不只要windfarm_opt=1就能算？

dueler

xinshou6 发表于 2020-9-29 20:35
他可能看的使用 WRF-Fitch 对湖区风电场风力发电机尾流效应特征的数值模拟，Fitch 等[19–20]将一个风场参 ...

我上边的问题大佬有什么见解么