CN111967203A - 一种半解析半数值的大气边界层三维台风风场建模方法 - Google Patents

一种半解析半数值的大气边界层三维台风风场建模方法 Download PDF

Info

Publication number
CN111967203A
CN111967203A CN202010793686.XA CN202010793686A CN111967203A CN 111967203 A CN111967203 A CN 111967203A CN 202010793686 A CN202010793686 A CN 202010793686A CN 111967203 A CN111967203 A CN 111967203A
Authority
CN
China
Prior art keywords
typhoon
boundary layer
wind speed
equation
ground
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010793686.XA
Other languages
English (en)
Other versions
CN111967203B (zh
Inventor
段忠东
杨剑
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Graduate School Harbin Institute of Technology
Original Assignee
Shenzhen Graduate School Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen Graduate School Harbin Institute of Technology filed Critical Shenzhen Graduate School Harbin Institute of Technology
Priority to CN202010793686.XA priority Critical patent/CN111967203B/zh
Publication of CN111967203A publication Critical patent/CN111967203A/zh
Application granted granted Critical
Publication of CN111967203B publication Critical patent/CN111967203B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明涉及一种半解析半数值的大气边界层三维台风风场建模方法，包括以下步骤：S1.根据大气边界层内空气微团的动力学平衡理论，建立相对于地面的空气微团动力学平衡方程；S2.根据混合边界层理论和梯度流理论，得到台风风场近地面的边界条件；S3.求解边界层内相对于地面的台风风速；S4.将垂直风速分量的初始值取为0进行迭代计算，设置收敛条件，取满足收敛条件的风速值作为最终边界层内相对于地面的台风风速。本发明的三维台风风场模型考虑了多种物理过程，包括：水平平流、垂直平流和垂直扩散过程。相比于其他台风模型，该模型计算效率高，能够得到更为真实的三维台风风场结构，且模拟结果与实测情况很接近。

Description

一种半解析半数值的大气边界层三维台风风场建模方法
技术领域
本发明涉及风速解析方法，特别涉及一种半解析半数值的大气边界层三维台风风场建模方法。
背景技术
现有的台风风速计算方法有：
(1)梯度风速折减法。采用风速折减系数把梯度高度风速折算到近地面高度风速，比如10m，通常折减系数在0.6～0.9之间。
梯度风速折减系数法无法描述台风风场的空间结构，只能近似求出近地面风场。近地面风速的大小受折减系数控制，但折减系数随台风的个例、登陆情况、观测地点等不同发生变化，因此很难通过折减系数精确描述台风风场。
(2)边界层平板模型。假设台风边界层高度沿径向不变，建立沿边界层高度平均的动力方程，使用近地风速与平均风速的参数关系，计算近地面台风风场。
边界层平板模型通常假设台风边界层高度不发生变化，而事实上台风边界层高度随着离台风中心距离的增加而增大；而且不同台风的边界层高度也不相同。平板模型的另一个缺点是对控制方程进行了沿高度的平均，忽略了垂直方向上气体的对流，使得平板模型无法正确计算垂直方向上的风速。
(3)沿高度解析的边界层模型。该模型简化了边界层内空气微团的动力学方程，仅考虑了对台风风场影响较大的几个物理过程，该模型可以解析地求出台风风速。
沿高度解析的边界层模型可以解析地计算不同高度处的台风风速，但是该模型简化了很多物理过程，使得其模拟的台风风场结构不准确，特别是对离地面300m以上的风速模拟，其模拟值偏小。目前随着超高层建筑的建设和对台风强降水灾害关注度的提升，该模型的缺陷导致其适用范围受到很大影响。
发明内容
本发明提供一种半解析半数值的大气边界层三维台风风场建模方法，旨在一个能够更精确更高效地模拟三维台风风场的模型。
本发明提供一种半解析半数值的大气边界层三维台风风场建模方法，包括以下步骤：
S1.根据大气边界层内空气微团的动力学平衡理论，建立相对于地面的空气微团动力学平衡方程；
S2.根据混合边界层理论和梯度流理论，得到台风风场近地面的边界条件；
S3.求解边界层内相对于地面的台风风速；
S4.将垂直风速分量的初始值取为0进行迭代计算，设置收敛条件，取满足收敛条件的风速值作为最终边界层内相对于地面的台风风速。
作为本发明的进一步改进，所述步骤S1包括：
根据大气边界层内空气微团的动力学平衡理论，建立相对于地面的空气微团动力学平衡方程：
Figure BDA0002624689700000021
式中，
Figure BDA0002624689700000022
f为科氏参数；k为垂直方向的单位向量；ρ为空气密度；F为摩擦力；
Figure BDA0002624689700000023
Figure BDA0002624689700000024
为相对于地面的台风风速，V为相对于台风中心的风速，Vc为台风移动速度；
Figure BDA0002624689700000025
Figure BDA0002624689700000026
为气压场，pc为台风气压场，pL为大尺度背景气压场，-ΔpL/ρ＝f·k×Vc；
方程(1)可改写为：
Figure BDA0002624689700000027
方程(2)中，V＝vg+v′，其中vg为梯度风速，v′为地表摩擦引起的风速；
对方程(2)进行分解，得到：
Figure BDA0002624689700000031
Figure BDA0002624689700000032
模拟的台风风场是逐小时平均的，可忽略方程(3)和方程(4)中的时间不稳定项
Figure BDA0002624689700000033
和
Figure BDA0002624689700000034
对方程(3)在柱坐标系r,λ,z下展开，vg为vg的切向分量，其径向分量ug＝0假设梯度风场vg为轴对称分布，可得到，
Figure BDA0002624689700000035
其中，
Figure BDA0002624689700000036
式中，p0为台风中心气压；pc为台风气压场；Δp＝1010-p0为台风中心气压差；Rmax为台风最大风速半径；B为气压剖面参数；Rmax和B可通过下式计算：
Rmax＝exp(3.015-6.291×10-5Δp2+0.0337·lat) (7)
Figure BDA0002624689700000037
式中，lat为台风中心的纬度；
对方程(4)在柱坐标系r,λ,z下展开，并忽略高阶项，得到控制方程为，
Figure BDA0002624689700000038
Figure BDA0002624689700000039
式中，Kv为湍流扩散系数，取Kv＝50m2/s；u′,v′和w′是v′在柱坐标下沿径向、切向和高度的分量；在方程(9)和(10)中，u′,v′,w′为未知量，需要通过连续性方程使其闭合，
Figure BDA0002624689700000041
作为本发明的进一步改进，所述步骤S2具体包括：
根据混合边界层理论和梯度流理论，得到近地面边界条件：
Figure BDA0002624689700000042
Figure BDA0002624689700000043
其中，Cd为摩阻力系数，
Figure BDA0002624689700000044
κ＝0.4为卡曼常数，z1＝1m为近地面高度，z0为地表粗糙长度；u和v为
Figure BDA0002624689700000045
在径向和切向的分量，且u＝u′+uc，v＝vg+v′+vc，
Figure BDA0002624689700000046
uc和vc为Vc在径向和切向的分量，其指数形式表示为，
Figure BDA0002624689700000047
Figure BDA0002624689700000048
把方程(12)和(13)展开，并进行简化，得到边界条件为，
Figure BDA0002624689700000049
Figure BDA00026246897000000410
作为本发明的进一步改进，所述步骤S3具体包括S31：
对控制方程(9)和(10)重新表达为：
Figure BDA00026246897000000411
其中，
Figure BDA00026246897000000412
Figure BDA00026246897000000413
式中，i为单位虚数；方程(18)中的解ω可以进行如下形式的傅里叶展开，
Figure BDA0002624689700000051
其中，ak(z)为傅里叶系数，k为傅里叶展开级数；把方程(21)带入方程(18)中，得到：
Figure BDA0002624689700000052
假设方程(22)的解的形式为：
ak(z)＝Akexp(pkz) (23)
式中，Ak和pk为待求未知量；
方程(22)变形为：
Figure BDA0002624689700000053
通过边界条件方程(16)和(17)求解Ak和pk，求得ω，再用公式(19)得出u′和v′。
作为本发明的进一步改进，所述步骤S3还包括步骤S32：
根据方程(19)和(20)，提取出u′和v′的表达式，并带入方程(16)和(17)中得到，
Figure BDA0002624689700000054
Figure BDA0002624689700000055
对于方程(25)和(26)，仅考虑|k|≤1的情况；经过一系列公式运算得到：
当
Figure BDA0002624689700000056
时，
Figure BDA0002624689700000057
Figure BDA0002624689700000058
Figure BDA0002624689700000059
当
Figure BDA0002624689700000061
时，
Figure BDA0002624689700000062
Figure BDA0002624689700000063
其中，
Figure BDA0002624689700000064
Figure BDA0002624689700000065
为pk的共轭复数；
得到u′和v′的表达式为：
当
Figure BDA0002624689700000066
时，
Figure BDA0002624689700000067
v′(r,λ,z)＝Im[A-1exp(p-1·z-iλ)+A0exp(p0·z)+A1exp(p1·z+iλ)] (33)
当
Figure BDA0002624689700000068
时，
Figure BDA0002624689700000069
v′(r,λ,z)＝Im[A′-1exp(p′-1·z-iλ)+A′0exp(p′0·z)+A′1exp(p′1·z+iλ)] (35)
边界层内相对于地面的风速表示为：
u(r,λ,z)＝u′(r,λ,z)+uc (36)
v(r,λ,z)＝v′(r,λ,z)+vg+vc (37)
其中，
Figure BDA00026246897000000610
Figure BDA00026246897000000611
为台风移动方位角；垂直风数w(w′)通过对方程(11)进行数值求解得到：
Figure BDA00026246897000000612
作为本发明的进一步改进，所述步骤S4包括：
S41.取wn的初始值w0＝0，得到δ1＝w0/2Kv＝0，用建立的台风风场模型计算出n＝1时的风速(u1,v1,w1)，其中n为迭代计算的次数；
S42.在n＝2时，取δ2＝w1/2Kv，计算出风速(un,vn,wn)；
S43.设定
Figure BDA0002624689700000071
判断是否max(|V1-V2|)≤0.01或者是否n≥100；若否，则执行步骤S44，若是则执行步骤S45；
S44.取n＝n+1，δn+1＝wn/2Kv后执行重新步骤S42；
S45.最后确定台风风场模型计算得的(un,vn,wn)为边界层风速。
本发明的有益效果是：本发明的三维台风风场模型考虑了多种物理过程，包括：水平平流、垂直平流和垂直扩散过程。相比于其他台风模型，该模型计算效率高，能够得到更为真实的三维台风风场结构，且模拟结果与实测情况很接近。
附图说明
图1是本发明中半解析半数值的边界层台风风场模型求解流程图；
图2是本发明模型模拟结果与H*Wind风场结果对比图；
图3是本发明模型模拟结果与下投式探空仪实测平均风剖面对比图；
图4是本发明模拟台风Hagupit风速风向与实测风速风向对比图；
图5是本发明模拟飓风Irene风速风向与实测风速风向对比图；
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。
如图1所示，本发明一种半解析半数值的大气边界层三维台风风场建模方法，包括以下步骤：
S1.控制方程：根据大气边界层内空气微团的动力学平衡理论，建立相对于地面的空气微团动力学平衡方程：
Figure BDA0002624689700000072
式中，
Figure BDA0002624689700000081
f为科氏参数；k为垂直方向的单位向量；ρ为空气密度；F为摩擦力；
Figure BDA0002624689700000082
Figure BDA0002624689700000083
为相对于地面的台风风速，V为相对于台风中心的风速，Vc为台风移动速度；
Figure BDA0002624689700000084
Figure BDA0002624689700000085
为气压场，pc为台风气压场，pL为大尺度背景气压场，-ΔpL/ρ＝f·k×Vc。
方程(1)可改写为：
Figure BDA0002624689700000086
方程(2)中，V＝vg+v′，其中vg为梯度风速，v′为地表摩擦引起的风速。对方程(2)进行分解，得到：
Figure BDA0002624689700000087
Figure BDA0002624689700000088
本发明模拟的台风风场是逐小时平均的，可忽略方程(3)和方程(4)中的时间不稳定项
Figure BDA0002624689700000089
和
Figure BDA00026246897000000810
对方程(3)在柱坐标系(r,λ,z)下展开，vg为vg的切向分量，其径向分量ug＝0假设梯度风场vg为轴对称分布，可得到：
Figure BDA00026246897000000811
其中，
Figure BDA00026246897000000812
式中，p0为台风中心气压；Δp＝1010-p0为台风中心气压差；Rmax为台风最大风速半径；B为气压剖面参数。Rmax和B可通过下式计算：
Rmax＝exp(3.015-6.291×10-5Δp2+0.0337·lat)， (7)
Figure BDA00026246897000000813
式中，lat为台风中心的纬度。
对方程(4)在柱坐标系(r,λ,z)下展开，并忽略高阶项，可以得到控制方程为，
Figure BDA0002624689700000091
Figure BDA0002624689700000092
式中，Kv为湍流扩散系数，取Kv＝50m2/s；u′,v′和w′是v′在柱坐标下沿径向、切向和高度的分量。在方程(9)和(10)中，u′,v′,w′为未知量，需要通过连续性方程使其闭合，
Figure BDA0002624689700000093
S2.边界条件：根据混合边界层理论和梯度流理论，可以得到近地面边界条件：
Figure BDA0002624689700000094
Figure BDA0002624689700000095
其中，Cd为摩阻力系数，
Figure BDA0002624689700000096
κ＝0.4为卡曼常数，z1＝1m为近地面高度，z0为地表粗糙长度；u和v为
Figure BDA0002624689700000097
在径向和切向的分量，且u＝u′+uc，v＝vg+v′+vc，
Figure BDA0002624689700000098
uc和vc为Vc在径向和切向的分量，其指数形式可表示为，
Figure BDA0002624689700000099
Figure BDA00026246897000000910
把方程(12)和(13)展开，并进行简化，可得到边界条件为，
Figure BDA00026246897000000911
Figure BDA00026246897000000912
S3.求解边界层风速，控制方程(9)和(10)可以重新表达为：
Figure BDA00026246897000000913
其中，
Figure BDA0002624689700000101
Figure BDA0002624689700000102
式中，i为单位虚数。方程(18)中的解ω可以进行如下形式的傅里叶展开，
Figure BDA0002624689700000103
其中，ak(z)为傅里叶系数，k为傅里叶展开级数。把方程(21)带入方程(18)中，得到，
Figure BDA0002624689700000104
假设方程(22)的解的形式为，
ak(z)＝Akexp(pkz)， (23)
式中，Ak和pk为待求未知量。方程(22)可变形为，
Figure BDA0002624689700000105
现在，我们需要通过边界条件(方程(16)和(17))求解Ak和pk，从而求得ω，再用公式(19)得出u′和v′。根据方程(19)和(20)，提取出u′和v′的表达式，并带入方程(16)和(17)中得到，
Figure BDA0002624689700000106
Figure BDA0002624689700000107
对于方程(25)和(26)，我们仅考虑|k|≤1的情况。经过一系列公式运算我们可以得到，
当
Figure BDA0002624689700000108
时，
Figure BDA0002624689700000109
Figure BDA0002624689700000111
Figure BDA0002624689700000112
当
Figure BDA0002624689700000113
时，
Figure BDA0002624689700000114
Figure BDA0002624689700000115
其中，
Figure BDA0002624689700000116
Figure BDA0002624689700000117
为pk的共轭复数。
因此，我们可以得到u′和v′的表达式为：
当
Figure BDA0002624689700000118
时，
Figure BDA0002624689700000119
v′(r,λ,z)＝Im[A-1exp(p-1·z-iλ)+A0exp(p0·z)+A1exp(p1·z+iλ)]； (33)
当
Figure BDA00026246897000001110
时，
Figure BDA00026246897000001111
v′(r,λ,z)＝Im[A′-1exp(p′-1·z-iλ)+A′0exp(p′0·z)+A′1exp(p′1·z+iλ)]。 (35)
最后，边界层内相对于地面的风速可表示为，
u(r,λ,z)＝u′(r,λ,z)+uc， (36)
v(r,λ,z)＝v′(r,λ,z)+vg+vc， (37)
其中，
Figure BDA0002624689700000121
Figure BDA0002624689700000122
为台风移动方位角。垂直风数w(w′)可通过对方程(11)进行数值求解得到，
Figure BDA0002624689700000123
S4.将垂直风速分量的初始值取为0进行迭代计算，设置收敛条件，取满足收敛条件的风速值作为最终边界层内相对于地面的台风风速。
如图1所示，具备为S41.取wn的初始值w0＝0，得到δ1＝w0/2Kv＝0，用建立的台风风场模型计算出n＝1时的风速(u1,v1,w1)，其中n为迭代计算的次数，δn是垂直速度分量w与两倍湍流扩散系数Kv的比值；
S42.在n＝2时，取δ2＝w1/2Kv，计算出风速(un,vn,wn)；
S43.设定
Figure BDA0002624689700000124
判断是否max(|V1-V2|)≤0.01或者是否n≥100；若否，则执行步骤S44，若是则执行步骤S45；
S44.取n＝n+1，δn+1＝wn/2Kv后执行重新步骤S42；
S45.最后确定台风风场模型计算得的(un,vn,wn)为边界层风速。
发明模型模拟结果与H*Wind风场结果对比情况见图2。H*Wind是美国国家海洋和大气管理局(NOAA)飓风研究所，结合飓风模型和实测数据，同化得到的10m高度处1-min平均的台风风场(图2中已转换到1h平均)；该套数据现已转交给RMS公司管理，可从网上免费获取。从图2中的对比可以发现，本发明模型模拟的10m高度处风场与H*Wind数据很接近，尤其是在最大风速和台风风场二维结构上。通过误差分析可以发现，本模型模拟的风速在最大风速半径附近比H*Wind值要偏大，最大偏差可达到20％；在台风外围，本模型模拟的风速比H*Wind值要偏小约5％。造成偏差的原因可能是我们使用的梯度风速(方程(5))与H*Wind数据或实际情况不一致导致，这也指明了下一步风场模型的改进方向。
本发明模型模拟结果与下投式探空仪实测平均风剖面的对比见图3。图3显示了不同平均风速区间的风剖面形状。本模型模拟结果与实测结果很接近，说明本模型能够准确模拟台风平均状态的风剖面。
本发明模型模拟风速风向与实测风速风向对比见图4、图5。台风Hagupit(2008)为登陆我国广东省的台风；飓风Irene(2011)是影响美国墨西哥湾的台风。从图4、图5中可以看出，我们模拟的风速值与实测的风速风向值很接近；说明本模型能合理计算实际台风风场，且该模型的适用性可不受地域限制。
本发明采用半解析半数值的方法计算大气边界层三维台风风场。能够在不牺牲运算速度的条件下实现对台风风场空间结构的合理模拟。本模型考虑了多个物理过程，包括垂直平流和垂直扩展；这使得本模型的模拟结果与实测情况更为接近。本模型在空间上是解析的，即可得到任意高度任意位置处的台风风速。因此本模型不需要使用风速折减系数和台风边界层固定高度假设，克服了之前模型的缺点。
本模型具有较为广泛的应用，包括：
(1)为政府部门提供灾害预警；
(2)为巨灾保险和再保险行业提供高精度的台风风场模型；
(3)为沿海受台风影响区域的输电线路抗风设计及灾害分析提供数据基础；
(4)可结合地形和地表覆盖数据模拟登陆台风的风场结构，为台风影响区域的建筑物及构筑物设计提供抗风标准；
(5)结合台风降水模型，可提供登陆台风的极值降水预测；
(6)结合水动力模型，可提供登陆台风时的极值潮位。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。

Claims (6)

1.一种半解析半数值的大气边界层三维台风风场建模方法，其特征在于，包括以下步骤：
S1.根据大气边界层内空气微团的动力学平衡理论，建立相对于地面的空气微团动力学平衡方程；
S2.根据混合边界层理论和梯度流理论，得到台风风场近地面的边界条件；
S3.求解边界层内相对于地面的台风风速；
S4.将垂直风速分量的初始值取为0进行迭代计算，设置收敛条件，取满足收敛条件的风速值作为最终边界层内相对于地面的台风风速。
2.根据权利要求1所述半解析半数值的大气边界层三维台风风场建模方法，其特征在于，所述步骤S1包括：
根据大气边界层内空气微团的动力学平衡理论，建立相对于地面的空气微团动力学平衡方程：
Figure FDA0002624689690000011
式中，
Figure FDA0002624689690000012
f为科氏参数；k为垂直方向的单位向量；ρ为空气密度；F为摩擦力；
Figure FDA0002624689690000013
Figure FDA0002624689690000014
为相对于地面的台风风速，V为相对于台风中心的风速，Vc为台风移动速度；
Figure FDA0002624689690000015
Figure FDA0002624689690000016
为气压场，pc为台风气压场，pL为大尺度背景气压场，-ΔpL/ρ＝f·k×Vc；
方程(1)可改写为：
Figure FDA0002624689690000017
方程(2)中，V＝vg+v′，其中vg为梯度风速，v′为地表摩擦引起的风速；
对方程(2)进行分解，得到：
Figure FDA0002624689690000018
Figure FDA0002624689690000019
模拟的台风风场是逐小时平均的，可忽略方程(3)和方程(4)中的时间不稳定项
Figure FDA0002624689690000021
和
Figure FDA0002624689690000022
对方程(3)在柱坐标系r,λ,z下展开，vg为vg的切向分量，其径向分量ug＝0假设梯度风场vg为轴对称分布，可得到，
Figure FDA0002624689690000023
其中，
Figure FDA0002624689690000024
式中，p0为台风中心气压；pc为台风气压场；Δp＝1010-p0为台风中心气压差；Rmax为台风最大风速半径；B为气压剖面参数；Rmax和B可通过下式计算：
Rmax＝exp(3.015-6.291×10-5Δp2+0.0337·lat) (7)
Figure FDA0002624689690000025
式中，lat为台风中心的纬度；
对方程(4)在柱坐标系r,λ,z下展开，并忽略高阶项，得到控制方程为，
Figure FDA0002624689690000026
Figure FDA0002624689690000027
式中，Kv为湍流扩散系数，取Kv＝50m2/s；u′,v′和w′是v′在柱坐标下沿径向、切向和高度的分量；在方程(9)和(10)中，u′,v′,w′为未知量，需要通过连续性方程使其闭合，
Figure FDA0002624689690000028
3.根据权利要求2所述半解析半数值的大气边界层三维台风风场建模方法，其特征在于，所述步骤S2具体包括：
根据混合边界层理论和梯度流理论，得到近地面边界条件：
Figure FDA0002624689690000031
Figure FDA0002624689690000032
其中，Cd为摩阻力系数，
Figure FDA0002624689690000033
κ＝0.4为卡曼常数，z1＝1m为近地面高度，z0为地表粗糙长度；u和v为
Figure FDA0002624689690000034
在径向和切向的分量，且u＝u′+uc，v＝vg+v′+vc，
Figure FDA0002624689690000035
uc和vc为Vc在径向和切向的分量，其指数形式表示为，
Figure FDA0002624689690000036
Figure FDA0002624689690000037
把方程(12)和(13)展开，并进行简化，得到边界条件为，
Figure FDA0002624689690000038
Figure FDA0002624689690000039
4.根据权利要求3所述半解析半数值的大气边界层三维台风风场建模方法，其特征在于，所述步骤S3具体包括S31：
对控制方程(9)和(10)重新表达为：
Figure FDA00026246896900000310
其中，
Figure FDA00026246896900000311
Figure FDA00026246896900000312
式中，i为单位虚数；方程(18)中的解ω可以进行如下形式的傅里叶展开，
Figure FDA00026246896900000313
其中，ak(z)为傅里叶系数，k为傅里叶展开级数；把方程(21)带入方程(18)中，得到：
Figure FDA0002624689690000041
假设方程(22)的解的形式为：
ak(z)＝Akexp(pkz) (23)
式中，Ak和pk为待求未知量；
方程(22)变形为：
Figure FDA0002624689690000042
通过边界条件方程(16)和(17)求解Ak和pk，求得ω，再用公式(19)得出u′和v′。
5.根据权利要求4所述半解析半数值的大气边界层三维台风风场建模方法，其特征在于，所述步骤S3还包括步骤S32：
根据方程(19)和(20)，提取出u′和v′的表达式，并带入方程(16)和(17)中得到，
Figure FDA0002624689690000043
Figure FDA0002624689690000044
对于方程(25)和(26)，仅考虑|k|≤1的情况；经过一系列公式运算得到：
当
Figure FDA0002624689690000045
时，
Figure FDA0002624689690000046
Figure FDA0002624689690000047
Figure FDA0002624689690000048
当
Figure FDA0002624689690000049
时，
Figure FDA0002624689690000051
Figure FDA0002624689690000052
其中，
Figure FDA0002624689690000053
Figure FDA0002624689690000054
为pk的共轭复数；
得到u′和v′的表达式为：
当
Figure FDA0002624689690000055
时，
Figure FDA0002624689690000056
v′(r,λ,z)＝Im[A-1exp(p-1·z-iλ)+A0exp(p0·z)+A1exp(p1·z+iλ)] (33)
当
Figure FDA0002624689690000057
时，
Figure FDA0002624689690000058
v′(r,λ,z)＝Im[A′-1exp(p′-1·z-iλ)+A′0exp(p′0·z)+A′1exp(p′1·z+iλ)] (35)
边界层内相对于地面的风速表示为：
u(r,λ,z)＝u′(r,λ,z)+uc (36)
v(r,λ,z)＝v′(r,λ,z)+vg+vc (37)
其中，
Figure FDA0002624689690000059
Figure FDA00026246896900000510
为台风移动方位角；垂直风数w(w′)通过对方程(11)进行数值求解得到：
Figure FDA00026246896900000511
6.根据权利要求5所述半解析半数值的大气边界层三维台风风场建模方法，其特征在于，所述步骤S4包括：
S41.取垂直风速分量wn的初始值w0＝0，得到δ1＝w0/2Kv＝0，用建立的台风风场模型计算出n＝1时的风速(u1,v1,w1)，其中n为迭代计算的次数；
S42.在n＝2时，取δ2＝w1/2Kv，计算出风速(un,vn,wn)；
S43.设定
Figure FDA0002624689690000061
判断是否max(|V1-V2|)≤0.01或者是否n≥100；若否，则执行步骤S44，若是则执行步骤S45；
S44.取n＝n+1，δn+1＝wn/2Kv后执行重新步骤S42；
S45.最后确定台风风场模型计算得的(un,vn,wn)为边界层风速。
CN202010793686.XA 2020-08-10 2020-08-10 一种半解析半数值的大气边界层三维台风风场建模方法 Active CN111967203B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010793686.XA CN111967203B (zh) 2020-08-10 2020-08-10 一种半解析半数值的大气边界层三维台风风场建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010793686.XA CN111967203B (zh) 2020-08-10 2020-08-10 一种半解析半数值的大气边界层三维台风风场建模方法

Publications (2)

Publication Number Publication Date
CN111967203A true CN111967203A (zh) 2020-11-20
CN111967203B CN111967203B (zh) 2024-05-24

Family

ID=73364829

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010793686.XA Active CN111967203B (zh) 2020-08-10 2020-08-10 一种半解析半数值的大气边界层三维台风风场建模方法

Country Status (1)

Country Link
CN (1) CN111967203B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112347715A (zh) * 2020-11-25 2021-02-09 中国人民解放军国防科技大学 一种诊断台风生成正压能量来源的新方法
CN112748255A (zh) * 2020-12-29 2021-05-04 国网河北省电力有限公司保定供电分公司 基于台风模型的输电铁塔台风风速风向计算方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090240352A1 (en) * 2008-03-18 2009-09-24 Powell Mark D Predicting tropical cyclone destructive potential by integrated kinetic energy according to the powell/reinhold scale
KR100974958B1 (ko) * 2009-09-29 2010-08-10 제주대학교 산학협력단 태풍－파랑－해양 접합 시스템과 이를 이용한 태풍 예측 방법
KR20150117972A (ko) * 2014-04-11 2015-10-21 한국해양과학기술원 연안국지 파랑예측 방법 및 시스템
CN107194097A (zh) * 2017-05-27 2017-09-22 中国大唐集团科学技术研究院有限公司 基于风电场气动模拟和风速风向数据的分析方法
CN107273334A (zh) * 2017-06-09 2017-10-20 哈尔滨工业大学深圳研究生院 移动台风边界层三分量风速解析方法
CN107657331A (zh) * 2017-08-16 2018-02-02 哈尔滨工业大学深圳研究生院 基于边界层模型的台风降水解析方法
CN110298115A (zh) * 2019-07-02 2019-10-01 中国气象局上海台风研究所 一种基于简化地形气动参数的风场动力降尺度方法
CN110399625A (zh) * 2019-09-05 2019-11-01 重庆大学 一种改进的台风风场经验模型
CN111310109A (zh) * 2020-03-13 2020-06-19 中铁二院工程集团有限责任公司 一种基于vmd-arma-garch模型的非良态风速建模方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090240352A1 (en) * 2008-03-18 2009-09-24 Powell Mark D Predicting tropical cyclone destructive potential by integrated kinetic energy according to the powell/reinhold scale
KR100974958B1 (ko) * 2009-09-29 2010-08-10 제주대학교 산학협력단 태풍－파랑－해양 접합 시스템과 이를 이용한 태풍 예측 방법
KR20150117972A (ko) * 2014-04-11 2015-10-21 한국해양과학기술원 연안국지 파랑예측 방법 및 시스템
CN107194097A (zh) * 2017-05-27 2017-09-22 中国大唐集团科学技术研究院有限公司 基于风电场气动模拟和风速风向数据的分析方法
CN107273334A (zh) * 2017-06-09 2017-10-20 哈尔滨工业大学深圳研究生院 移动台风边界层三分量风速解析方法
CN107657331A (zh) * 2017-08-16 2018-02-02 哈尔滨工业大学深圳研究生院 基于边界层模型的台风降水解析方法
CN110298115A (zh) * 2019-07-02 2019-10-01 中国气象局上海台风研究所 一种基于简化地形气动参数的风场动力降尺度方法
CN110399625A (zh) * 2019-09-05 2019-11-01 重庆大学 一种改进的台风风场经验模型
CN111310109A (zh) * 2020-03-13 2020-06-19 中铁二院工程集团有限责任公司 一种基于vmd-arma-garch模型的非良态风速建模方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
JEFF KEPERT: "The Dynamics of Boundary Layer Jets within the Tropical Cyclone Core. Part I: Linear Theory", JOURNAL OF THE ATMOSPHERIC SCIENCES, vol. 58, 1 September 2001 (2001-09-01), pages 2469 - 2484 *
JEFFREY D. KEPERT: "Choosing a Boundary Layer Parameterization for Tropical Cyclone Modeling", MONTHLY WEATHER REVIEW, 1 May 2012 (2012-05-01), pages 1427 - 1445 *
REDA SNAIKI 等: "Modeling tropical cyclone boundary layer: Height-resolving pressure and wind fields", JOURNAL OF WIND ENGINEERING & INDUSTRIAL AERODYNAMICS, 31 December 2017 (2017-12-31), pages 18 - 27 *
W. LI 等: "A discussion on influences of turbulent diffusivity and surface drag parameterizations using a linear model of the tropical cyclone boundary layer wind field", ATMOSPHERIC RESEARCH, 30 June 2020 (2020-06-30), pages 1 - 38 *
杨剑: "WRF与CFD嵌套的局地台风风场数值模拟研究", 中国优秀硕士学位论文全文数据库, no. 2, 15 February 2016 (2016-02-15), pages 038 - 397 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112347715A (zh) * 2020-11-25 2021-02-09 中国人民解放军国防科技大学 一种诊断台风生成正压能量来源的新方法
CN112748255A (zh) * 2020-12-29 2021-05-04 国网河北省电力有限公司保定供电分公司 基于台风模型的输电铁塔台风风速风向计算方法

Also Published As

Publication number Publication date
CN111967203B (zh) 2024-05-24

Similar Documents

Publication Publication Date Title
Coffer et al. Simulated supercells in nontornadic and tornadic VORTEX2 environments
JP6896180B2 (ja) 風流検知システム及び風流の速度場を求める方法
Chen et al. A deterministic–stochastic hybrid model of downbursts and its impact on a cantilevered structure
Floors et al. The wind profile in the coastal boundary layer: Wind lidar measurements and numerical modelling
CN114357571B (zh) 建成建筑环境下大气边界层风场特性的反演方法和系统
CN111967203A (zh) 一种半解析半数值的大气边界层三维台风风场建模方法
Wang et al. Data-driven CFD modeling of turbulent flows through complex structures
Bierdel et al. Accuracy of rotational and divergent kinetic energy spectra diagnosed from flight-track winds
US11294063B2 (en) System and method for fast wind flow measurement by LiDAR in a complex terrain
JP7351029B2 (ja) 複雑地形においてＬｉＤＡＲで風の流れの乱流を測定するためのシステムおよび方法
Zaron et al. Bottom topography mapping via nonlinear data assimilation
CN115525978A (zh) 用于飞行器波浪情况水动力分析的速度入口数值造波方法
CN115994496A (zh) 城市公园高分辨率大气co2浓度三维场的数值模拟方法
Hangan et al. Big data from big experiments. The WindEEE Dome
CN112001099A (zh) 大型射电望远镜天线面板脉动风压快速测算方法
Fan et al. A numerical rotating water tank can reproduce the Coriolis effect on the urban heat dome flow
Xiaoyong et al. Dynamical and microphysical retrieval from simulated Doppler radar observations using the 4DVAR assimilation technique
Paiva et al. Numerical simulation of atmospheric boundary layer flow over isolated and vegetated hills using RAMS
Han et al. Comparison of methods for curvature estimation from volume fractions
Wiersema et al. Development of a multiscale modeling framework for urban simulations in the Weather Research and Forecasting model
CN112149302B (zh) 一种基于非均匀下垫面的台风建模方法及相关装置
Casella et al. Dynamic flow analysis using an OpenFOAM based CFD tool: Validation of Turbulence Intensity in a testing site
Gairola Generic numerical tornado model for common interpretation of existing experimental simulators
Kepert A parametric model of tropical cyclone surface winds for sea and land
Oreskovic Numerical investigation of full scale thunderstorm downbursts: A parametric study and comparison to meteorological model

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant