设置大气边界层

要设置大气边界层 (ABL),可将地面指定为粗糙壁面边界。对于标准 K-Epsilon 湍流模型,可使用大气边界层模型来设置入口边界分布。如果要使用其他湍流模型,则必须使用标准分布(例如,场函数)指 ABL(例如风速)的入口值。

网格注意事项

根据灵敏度研究 [349],建议在为 ABL 模拟设置网格时使用以下值:

  • 确保 y > 0.2 r y 10 y 0 ,其中 y 是指从壁面到网格单元形心的最短距离, r 是指粗糙度高度, y 0 是指地面。
  • 确保在 Richards Hoxley 模型 [349] 适用的边界层的对数区域中, y + 30

要使用标准 K-Epsilon 模型对 ABL 进行建模:

  1. 对于物理连续体,停用自动选择推荐模型,并按顺序选择下列模型:
    组合框 模型
    时间 定常
    材料 气体
    流体 耦合流体分离流
    状态方程 恒密度
    粘滞态 湍流
    湍流 雷诺平均纳维-斯托克斯
    雷诺平均湍流 K-Epsilon 湍流
    K-Epsilon 湍流模型 标准 K-Epsilon
    K-Epsilon 高 y+ 壁面处理 高 y+ 壁面处理
    可选模型 大气边界层
  2. 选择连续体 > [物理连续体] > 模型 > 大气边界层节点,然后设置下列属性:
    属性 设置
    参考速度 给定高度的风速 ,Eqn. (1375) 中。 u h 默认值为 10.0 m/s
    参考高度 采用风速的高度,,Eqn. (1375) 中。 h 默认值为 6.0 m
    粗糙度长度 用于对地面上的流动障碍建模的大小,,Eqn. (1372) 中。 y 0 默认值为 0.01 m
    冯·卡门常数 ABL 速度分布 Eqn. (1374) 中使用的无量纲系数 。 κ 默认值为 0.42
    有关详细信息,请参见大气边界层模型参考
  3. 选择 [物理连续体] > 模型 > 标准 K-Epsilon 节点,然后设置下列属性:
    属性 设置
    Cmu 根据入口速度分布设置值。典型值介于 0.030.11 之间。
    C1e 使用以下方程设置值:
    1. EQUATION_DISPLAY
    c ϵ 1 = c ϵ 2 κ 2 c μ σ ϵ
    (258)

    有关详细信息,请参见 K-Epsilon 模型参考

  4. 要设置湍流粘度比,选择 [物理连续体] > 初始条件 > 湍流粘度比节点,将方法设为常数,将设为 1E10
  5. 使用本示例中的值将地面指定为粗糙壁面边界:
    y 0 = 0.01 m E = 8.75 C = 0.25 并且 y y 0 时,计算得出的粗糙度高度 r 为 0.35 m。
    1. 选择区域 > [区域] > 边界 > [地面边界] 节点
    2. 类型设为壁面,然后设置以下属性:
      节点 属性 设置
      物理条件
      壁面指定 方法 粗糙
      物理值
      标准壁面函数 E 8.75
      壁面粗糙度参数 C 0.25
      RplusRough (R+) 90 以上,因此 Eqn. (1377) 有效。
      粗糙度限制器 停用,因为 Eqn. (1379) 中的假设 有效。 y y 0
      粗糙度高度 价值 0.35 m
  6. 定义流入边界:
    1. 选择边界 > [速度入口边界]节点,然后将类型设为速度入口
    2. 编辑[速度入口边界]节点,然后设置下列属性:
      节点 属性 设置
      物理条件
      湍流指定 方法 K+Epsilon
      物理值
      速度幅值 方法 大气边界层
      湍流耗散率
      湍流动能
      大气边界层模型提供流体速度 u (在 Eqn. (1372) 中)、湍流耗散率 ε (在 Eqn. (1374) 中)和湍动能 k (在 Eqn. (1373) 中)。
如果要对 ABL 模拟使用标准 K-Omega 湍流模型,则选择 [物理连续体] > 模型 > 标准 K-Omega节点,然后设置以下属性:
属性 设置
β* 根据入口速度分布设置值,其中 β * c μ
β 使用以下方程设置值:
2. EQUATION_DISPLAY
β = β * 1.25
(259)
Alpha 使用以下方程设置值:
3. EQUATION_DISPLAY
α = 1 1.25 σ ω κ 2 β *
(260)
Sigma_w 0.5(默认)

有关详细信息,请参见 K-Omega 模型参考