物理速度公式

Simcenter STAR-CCM+ 中的多孔介质模型使用物理速度公式。此公式考虑了流进入多孔介质时的速度增加。

在下面的方程中,速度指物理速度。关于能量,Simcenter STAR-CCM+ 中的多孔介质模型提供非平衡或平衡建模。非平衡模型不假定液相和固相在相同温度下。它对液相和固相的单独能量传输方程进行求解。另一方面,平衡模型假定液相和固相在相同温度下。它对单个能量传输方程进行求解。

此公式假设孔隙率的阶跃变化只发生在交界面之间。区域内部孔隙率的任何阶跃变化都可能导致非物理结果。

体积分数
体积分数表示由 i 次多孔相所占的体积,与所有多孔相所占体积相对,与孔隙率和网格单元体积相关,如下所示:
α s i = V solid , i ∑ i   V solid , i = V solid , i V ⋅ 1 1 − χ
(1849)

其中:

  • V solid , i 为 i 次多孔固相的体积。
  • V 为网格单元体积。
  • χ 为孔隙率。
连续性方程
∂ ∂ t ∫ V ( ρ χ ) d V + ∮ A ( χ ρ v ) ⋅ d a = ∫ V S u c d V
(1850)

其中:

  • V 为体积
  • ρ 为密度
  • χ 为孔隙率
  • v 为物理速度
  • S u c 为用户自定义的源或汇。
动量方程
∂ ∂ t ∫ V ( χ ρ v ) d V + ∮ A χ ρ v ⊗ v ⋅ d a = − ∮ A χ p I ⋅ d a + ∮ A χ T ⋅ d a + ∫ V χ f b d V + ∫ V χ f p d V + ∫ V S u m d
(1851)

其中:

  • p 为压力。
  • I 为单位矩阵。
  • f b 为体积力矢量。
  • f p 为多孔阻力,其中:
    f p = P v v s + P i | v s | v s
    (1852)

    其中:

    • P v 为粘性阻力张量。
    • v s 为表观速度 v s = χ ⁢   v 。
    • P i 为惯性阻力张量。
  • S u m 为用户自定义的动量源或汇。

使用欧拉多相 (EMP) 模型时,多孔粘性阻力 P v 根据多孔介质的内在渗透率和每个液相的相对渗透率进行计算。请参见具有欧拉多相的多孔介质。

热非平衡能量方程
∂ ∂ t ∫ V ( χ ρ fluid E fluid ) d V + ∮ A χ ρ fluid H fluid v ⋅ d a = − ∮ A χ q fluid ⋅ d a + ∮ A χ T ⋅ v ⋅ d a + ∑ solid phases   ∫ V a s i h s i ( T fluid − T s i ) d V + ∫ V χ f b ⋅ v ⁢     d V + ∫ V S u e d V
(1853)
其中:
  • E fluid 为流体的总能量
  • H fluid 为流体的总焓
  • a s i 为 i 次固相的相互作用面积密度
  • h s i 为 i 次固相的有效导热率。
  • T 为应力张量
  • S u e 为用户自定义的流体能量源或能量汇。
  • T s i 为 i 次固相的固体温度
∂ ∂ t ∫ V ( ( 1 − χ ) α s i ρ s i E s i ) d V = − ∮ A ( 1 − χ ) α s i q s i ⋅ d a + ∫ V a s i h s i ( T s i − T fluid ) d V + ∫ V S u e d V
(1854)

其中:

  • ρ s i 为 i 次固相的密度
  • E s i 为 i 次固相的总能量
  • q s i 为 i 次固相的传导热通量
热平衡能量方程
∂ ∂ t ∫ V ( ρE ) eff ⁢     d V + ∮ A χ ρ fluid H fluid v ⋅ d a = ∮ A k eff ∇ T ⋅ d a + ∮ A χ T ⋅ v ⁢   ⋅   d a + ∫ V χ f b ⋅ v ⁢     d V + ∫ V S u e ⁢     d V
(1855)
( ρ E ) eff = χ ρ fluid E fluid + ( 1 − χ ) ρ solid T fluid
(1856)
k e ff = χ k fluid + ( 1 − χ ) k solid k solid = ∑ i   α s i k s i
(1857)

其中, k eff 为有效导热率,由流体导热率 k fluid 和固体导热率 k solid 计算,其中 k s i 为第 i 个固相的有效导热率。

多孔交界面处的条件

在多孔交界面处,物理速度是不连续的。

下表列出了各种变量的多孔交界面条件:
变量 多孔交界面条件
物理速度 不连续
表观速度 连续
静压力 不连续
总压力 连续