RANS(Reynolds-Averaged Navier-Stokes) 湍流模型类型

RANS(Reynolds-Averaged Navier-Stokes) 湍流模型有多种不同的类型,除了标准的 kkk-ω 湍流模型,还有其他一些常用的湍流模型。RANS 模型的核心思想是对 Navier-Stokes 方程进行 雷诺平均,通过将流动场的瞬时变量分解为平均量和脉动量,从而简化湍流模拟。不同的 RANS 模型采用不同的方式来描述湍流的封闭方程和湍流应力(雷诺应力),每种模型都有其适用范围和优缺点。

以下是几种常见的 RANS 湍流模型:


1. k-ε 湍流模型

k-ε(k-epsilon)湍流模型是最广泛应用的 RANS 模型之一。该模型基于湍流动能 k 和湍流耗散率 ε 的传输方程。它适用于大多数的工业流动,但在某些复杂流动(如边界层流动、流动分离等)中表现不佳。

方程形式
  • 湍流动能方程:与 k-ω 模型相似,描述了湍流动能 k 的传输过程。
  • 湍流耗散率方程:给出湍流能量的耗散过程,即湍流动能转化为热能的速率。
优点
  • 计算简单,收敛速度较快。
  • 适用于大部分高雷诺数的流动问题。
缺点
  • 对近壁流动和低雷诺数流动的预测不准确。
  • 对流动分离和复杂几何的预测精度较差。

2. k-ω 湍流模型

k-ω 湍流模型是另一种常用的 RANS 模型,与 k-ε 模型类似,但它使用湍流频率 ω 代替了湍流耗散率 ε。该模型特别适用于近壁流动的模拟。

优点
  • 对近壁流动和低雷诺数流动的预测较为准确。
  • 比 k-ε 模型在分离流动的预测上更为有效。
缺点
  • 对自由流(远离壁面的流动)中的湍流预测不准确。
  • 对流动中的湍流动能过度高估。

3. k-ω SST(Shear Stress Transport)模型

k-ω SST(剪切应力传输)模型是对标准 k-ω 模型的改进,它结合了 k-ω 模型在近壁区域的优势和 k-ε 模型在自由流中的优点。该模型通过使用 k-ω 模型近壁,k-ε 模型远离壁面,使得模型对整个流场的适应性更强。

优点
  • 适用于分离流动和复杂几何的湍流模拟。
  • 在预测湍流分离和涡轮叶片流动等工程问题中表现优异。
缺点
  • 相对 k-ε 模型,计算量更大。
  • 对一些高雷诺数流动仍可能不准确。

4. RSM(Reynolds Stress Model,雷诺应力模型)

雷诺应力模型(RSM) 是一种基于雷诺应力方程的 RANS 湍流模型。RSM 直接求解湍流应力(即湍流应力张量),相比于 kkk-ε 和 kkk-ω 模型,它不需要使用假设的涡粘性模型,能够更准确地描述湍流的各向异性。

优点
  • 适用于各向异性流动(例如强剪切流动或旋转流动)和复杂湍流现象。
  • 不需要涡粘性假设,能够更精确地捕捉湍流的真实应力。
缺点
  • 计算复杂度高,求解过程需要更多的计算资源。
  • 需要更多的实验数据来校准模型参数。

5. LES-RANS 混合模型

LES-RANS 混合模型是一种结合大涡模拟(LES)和雷诺平均方法(RANS)的模型。在流体力学中,LES 用于模拟大尺度涡流,通常在湍流模拟中提供更精细的预测,尤其适用于瞬态流动。RANS 用于小尺度湍流的建模,这种混合模型能够在保证计算效率的同时,提高预测精度。

优点
  • 结合了 LES 在大尺度流动中的高精度和 RANS 在小尺度流动中的高效性。
  • 适用于复杂流动的长时间模拟,提供更精确的结果。
缺点
  • 计算需求较高,需要更高的计算资源和复杂的数值方法。

总结:

RANS 湍流模型有多种不同类型,每种模型有其特定的应用场景和优缺点。常见的 RANS 模型包括:

  1. k-ε 湍流模型:计算简单,适用于高雷诺数流动,但在近壁流动和低雷诺数流动中表现较差。
  2. k-ω 湍流模型:适用于近壁流动,能更好地处理低雷诺数流动,但对自由流的不准确性较高。
  3. k-ω SST 模型:结合了 k-ω 和 k-ε 的优点,适用于分离流动和复杂几何,但计算复杂度较高。
  4. 雷诺应力模型(RSM):适用于各向异性流动,能够捕捉更复杂的湍流现象,但计算资源需求较大。
  5. LES-RANS 混合模型:结合 LES 和 RANS 模型的优势,适用于复杂流动问题,但计算需求较高。

根据具体的流动问题选择合适的 RANS 湍流模型非常重要,可以根据流动的复杂性、计算资源的限制以及精度要求来进行选择。

1. DNS(直接数值模拟)中的相关模型

直接数值模拟(DNS)是一种精确求解 Navier-Stokes 方程 的方法,不依赖于湍流模型,而是通过高分辨率的网格和小时间步长捕捉流体中所有的湍流尺度。尽管 DNS 不使用传统的湍流模型,但在实践中,可能会结合一些优化策略或算法来提高计算效率或适应不同的流动条件。

(1)高性能计算优化
  • 自适应网格生成方法(Adaptive Mesh Refinement, AMR)
    • 自适应网格方法根据流动特性动态调整网格分辨率,减少在不重要区域的计算成本,集中计算流动中的关键区域(如涡旋区域)。AMR 可以显著提高计算效率,尤其是在求解大尺度问题时。
  • 多尺度模拟方法
    • 在 DNS 中使用多尺度方法,通过不同尺度的细致模拟来减少计算量,尤其是在流动中涡旋的尺度范围极广时。多尺度方法将复杂的湍流问题分解为多个可管理的尺度进行计算。
(2)DNS 的应用与优化
  • 混合DNS-LES方法

    • 在 DNS 和 LES 之间进行组合,通常在计算中某些区域采用 DNS 来模拟小尺度涡旋,而在其他区域使用 LES 来捕捉大尺度涡旋。这种方法适用于一些高雷诺数的流动问题,通过减小计算需求同时保留较高的精度。
  • 大规模并行计算(Parallel DNS)

    • 在 DNS 中,通过并行计算来加速模拟过程,使用高性能计算集群(如超级计算机)处理更大规模的流动问题。通过这种方式,DNS 可以应用于更复杂的三维流动和高雷诺数湍流。

2. LES(大涡模拟)中的相关模型

大涡模拟(LES)是一种通过直接模拟大尺度涡旋、小尺度涡旋则通过模型来近似的湍流模拟方法。LES 需要较高的计算资源,尤其是在细节捕捉小尺度涡旋时。为了提高 LES 的准确性和计算效率,发展出了一些衍生模型。

(1)SFS 模型(Subgrid-Scale Model,亚格尺度模型)
  • 基础 SFS 模型
    • 在 LES 中,亚格尺度模型(SFS)用于模拟网格分辨率以下的小尺度涡旋(即亚网格涡旋)。SFS 模型的核心目的是通过合适的模型来逼近小尺度涡旋的贡献。
(2)混合 LES-RANS 模型
  • LES-RANS 混合模型

    • 在 LES 模型的基础上,结合 RANS 模型的优势。在流动的远离壁面区域使用 RANS 模型处理小尺度涡旋,而在壁面附近或流动分离区域使用 LES 来捕捉大尺度涡旋。这种方法结合了 LES 的高精度和 RANS 的计算效率,适用于处理更复杂的流动问题。
  • DES(Detached Eddy Simulation)模型

    • DES 是 LES 和 RANS 的混合模型,在近壁区使用 RANS 模型,而在远离壁面或大尺度涡旋区域使用 LES。这种方法通过在流动区域切换不同的模型来提高计算效率,并获得较好的模拟结果。
(3)动态 LES 方法
  • 动态 LES(D-LES)
    • 动态 LES 方法在进行网格划分时自适应地调整空间分辨率,根据流场的不同特征(如湍流强度、剪切流动等)来选择不同的模型或调节参数,优化计算效率并提高精度。

3. 总结

  • DNS
    • DNS 直接求解 Navier-Stokes 方程,捕捉所有涡旋的动态,是一种最精确的湍流模拟方法,但计算量巨大。通过 自适应网格(AMR)多尺度模拟方法DNS-LES 混合方法,可以在一些高雷诺数流动和大规模问题中提高计算效率。
  • LES
    • LES 通过模拟大尺度涡旋并用亚格尺度模型(SFS)近似小尺度涡旋来实现湍流模拟。常见的 SFS 模型有 Smagorinsky 模型Vreman 模型,并且通过 LES-RANS 混合模型DES 等方法来优化计算效率和准确性。动态 LES 方法则能够根据流动的不同特性自适应调整计算参数。

在实际应用中,选择合适的模型和方法取决于问题的具体要求,包括计算资源、精度需求和流动复杂性。

 

GARCH-M 模型是 GARCH 模型的一种扩展形式,它可以用来建模多个资产的波动率。其基本结构如下: $$ r_{t,i} = \mu_i + \epsilon_{t,i} \\ \epsilon_{t,i} = \sigma_{t,i} z_{t,i} \\ \sigma^2_{t,i} = \alpha_{i,0} + \sum_{j=1}^q \alpha_{i,j} \epsilon^2_{t-j,i} + \sum_{k=1}^p \beta_{i,k} \sigma^2_{t-k,i} + \sum_{l=1}^m \gamma_{i,l} \sigma^2_{t-l,j} $$ 其中,$r_{t,i}$ 表示第 $i$ 个资产在时间 $t$ 的收益率,$\mu_i$ 表示第 $i$ 个资产的均值,$\epsilon_{t,i}$ 表示第 $i$ 个资产在时间 $t$ 的收益率的误差项,$\sigma_{t,i}$ 表示第 $i$ 个资产在时间 $t$ 的波动率,$z_{t,i}$ 表示第 $i$ 个资产在时间 $t$ 的标准化残差,$\alpha_{i,j}$、$\beta_{i,k}$ 和 $\gamma_{i,l}$ 分别表示第 $i$ 个资产的 ARCH、GARCH 和 MGARCH 参数,$p$、$q$ 和 $m$ 分别表示 GARCH-M 模型的 GARCH 阶数、ARCH 阶数和 MGARCH 阶数。 下面是一个 GARCH-M 模型的例题: 假设有两个资产 $A$ 和 $B$,它们在时间 $t$ 的收益率分别为 $r_{t,A}$ 和 $r_{t,B}$。我们使用 GARCH-M 模型来建模这两个资产的波动率,其中 GARCH 阶数为 $1$,ARCH 阶数为 $1$,MGARCH 阶数为 $1$。已知参数如下: $$ \begin{aligned} & \alpha_{A,0} = 0.01, \quad \alpha_{A,1} = 0.05, \quad \beta_{A,1} = 0.90, \quad \gamma_{A,1} = 0.03 \\ & \alpha_{B,0} = 0.02, \quad \alpha_{B,1} = 0.10, \quad \beta_{B,1} = 0.80, \quad \gamma_{B,1} = 0.05 \\ \end{aligned} $$ 假设在时刻 $t=0$,$r_{0,A} = 0.02$,$r_{0,B} = -0.01$,$\sigma^2_{0,A} = 0.02$,$\sigma^2_{0,B} = 0.03$。求在接下来的 $3$ 个时间点中,$A$ 和 $B$ 的波动率分别为多少。 解: 根据 GARCH-M 模型的公式,我们可以先计算出 $A$ 和 $B$ 在时间 $t=1$ 的波动率: $$ \begin{aligned} \sigma^2_{1,A} & = \alpha_{A,0} + \alpha_{A,1} \epsilon^2_{0,A} + \beta_{A,1} \sigma^2_{0,A} + \gamma_{A,1} \sigma^2_{0,B} \\ & = 0.01 + 0.05 \times 0.02^2 + 0.90 \times 0.02 + 0.03 \times 0.03 \\ & = 0.051 \\ \sigma^2_{1,B} & = \alpha_{B,0} + \alpha_{B,1} \epsilon^2_{0,B} + \beta_{B,1} \sigma^2_{0,B} + \gamma_{B,1} \sigma^2_{0,A} \\ & = 0.02 + 0.10 \times (-0.01)^2 + 0.80 \times 0.03 + 0.05 \times 0.02 \\ & = 0.027 \end{aligned} $$ 接下来,我们可以依次计算出 $A$ 和 $B$ 在时间 $t=2$ 和 $t=3$ 的波动率: $$ \begin{aligned} \sigma^2_{2,A} & = \alpha_{A,0} + \alpha_{A,1} \epsilon^2_{1,A} + \beta_{A,1} \sigma^2_{1,A} + \gamma_{A,1} \sigma^2_{1,B} \\ & = 0.01 + 0.05 \times \left( \frac{r_{1,A} - \mu_A}{\sigma_{1,A}} \right)^2 + 0.90 \times 0.051 + 0.03 \times 0.027 \\ & = 0.058 \\ \sigma^2_{2,B} & = \alpha_{B,0} + \alpha_{B,1} \epsilon^2_{1,B} + \beta_{B,1} \sigma^2_{1,B} + \gamma_{B,1} \sigma^2_{1,A} \\ & = 0.02 + 0.10 \times \left( \frac{r_{1,B} - \mu_B}{\sigma_{1,B}} \right)^2 + 0.80 \times 0.027 + 0.05 \times 0.051 \\ & = 0.040 \\ \sigma^2_{3,A} & = \alpha_{A,0} + \alpha_{A,1} \epsilon^2_{2,A} + \beta_{A,1} \sigma^2_{2,A} + \gamma_{A,1} \sigma^2_{2,B} \\ & = 0.01 + 0.05 \times \left( \frac{r_{2,A} - \mu_A}{\sigma_{2,A}} \right)^2 + 0.90 \times 0.058 + 0.03 \times 0.040 \\ & = 0.064 \\ \sigma^2_{3,B} & = \alpha_{B,0} + \alpha_{B,1} \epsilon^2_{2,B} + \beta_{B,1} \sigma^2_{2,B} + \gamma_{B,1} \sigma^2_{2,A} \\ & = 0.02 + 0.10 \times \left( \frac{r_{2,B} - \mu_B}{\sigma_{2,B}} \right)^2 + 0.80 \times 0.040 + 0.05 \times 0.058 \\ & = 0.046 \end{aligned} $$ 因此,$A$ 在时间 $t=1$、$t=2$ 和 $t=3$ 的波动率分别为 $0.225$、$0.242$ 和 $0.253$,$B$ 在时间 $t=1$、$t=2$ 和 $t=3$ 的波动率分别为 $0.165$、$0.200$ 和 $0.214$。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值