扩散问题的有限体积法.doc

流体力学数值方法 第六讲 扩散问题的有限体积法 扩散问题的有限体积法 ◆通用形式流动与传热问题守恒形式的输运方程 ? ρφ → ( ) + ρ div( U ? t φ ) = div( Γ grad ) φ + S φ 瞬变项 对流项 扩散项 源项 ▼在应用有限体积法（控制容积法）进行数值求 解时，通常首先将通用公式在一个容积上进行积 分，将微分方程转化为积分方程，然后采用不同 的近似方式在控制容积的边界上对积分项进行处 理，从而得到不同的差分格式。 扩散问题的有限体积法 ◆有限体积法求解过程 ? ρφ → ( ) ∫ ∫ ∫ ∫ ? + = Γ + dV div( U )dV div( grad )dV S ρ φ φ φ t CV CV CV CV dV 高斯定理 r r ∫ = ∫ ? div(a)dV n CV A r adA → → r ∫ ∫ div(ρ φ) = ? (ρφ ) U dV n U dA CV A r ∫ Γ = ∫ ? Γ div ( φ ) ( φ ) grad dV n grad CV A dA ? → r r ∫( ) + ∫ ? = ∫ ? Γ + ∫ ρφ (ρφ ) ( φ) dV n U dA n grad dA S φ ? t CV A A CV dV 扩散问题的有限体积法 ◆有限体积法输运方程的物理意义 ? → r r ∫( ) + ∫ ? = ∫ ? Γ + ∫ ρφ (ρφ ) ( φ) dV n U dA n grad dA S φ ? t CV A A CV dV 总的变 外法线方 内法线方 源项引 化率 向的对流 通量 向的扩散 通量 起的的 增加率 物理意义 物理意义 因对流而 扩散而引 引起的净 起的净增 减少量 加量 扩散问题的有限体积法 ◆稳态输运方程 r → r ∫ ? = ∫ ? Γ + ∫ n (ρφ ) ( φ) φ U dA n grad dA S A A CV dV ◆非稳态输运方程 ? r → r ∫ ∫( ) ∫∫ ∫∫ ∫ ∫ ρφ + ? (ρφ ) = ? (Γ φ) + dV dt n U dAdt n grad dAdt S φ ? t Δ Δ Δ Δ t CV t A t A t CV dV dt 扩散问题的有限体积法 ◆稳态纯扩散 (Γ )+ = 0 div grad φ S φ ∫ (Γ + ∫ = div φ φ grad )dV S dV CV CV 0 ∫ A r n φ φ ? ∫ (Γ + = grad )dA S dV CV 0 一维稳态扩散问题的有限体积法 ◆一维稳态纯扩散方程 d ?Γ φ S d ? + ? ? dx ? dx ? = 0 ◆节点划分（P点） 有限体积法的第一步是把求解域划为离散的控制容积。 一维稳态扩散问题的有限体积法 ▼控制容积的取法 方法A：一种是把控制容积的界面放在相邻2个节点中间（先划分节点） 方法B：一种是把控制容积的中心节点放在控制容积的几何中心（先划 分控制容积） 一维稳态扩散问题的有限体积法 ◆方程的离散 φ φ φ ?Γ ∫ ? ? ? ?Γ ?Γ d d d d ∫ S V dV SdV A A + = ? + Δ ? ? ? ? ? ? dx ? ? ? ? ? ? dx dx dx CV CV e w = 0 ?Γ ? A ? ? d φ φ ?φ = Γ ? A E P e e δ dx ?e x PE ?Γ ? A ? ? d φ φ ?φ = Γ A ? P W w w δ dx ?w x WP Γ w = Γ W + 2 Γ P Γ e = Γ E + 2 Γ P 中心差分格式 S Δ = + φ V S u S P P 一维稳态扩散问题的有限体积法 ◆方程的离散 φ ?φ φ ?φ Γ A ? Γ + + φ e S S e δ u P P E P P W A w w δ x x PE WP = 0 ? ? ? ? Γ φ A ? + Γ ? ? Γ ? ? Γ ? A A A + ? ?φ = ? ?φ +? e S e S w w e e w w ? ? ? ? δ δ δ δ x ? x x ? P P E W u x ? ? ? PE WP PE WP a φ = φ + φ + P a a S P E E W W u a = + ? P a a S E W P Γ A a = w w W δx WP Γ A a = e e E δx PE 一维稳态扩散问题的有限体积法 ◆方程离散的步骤 首先将微分方程在控制容积上进行积分，利用高斯定理把体 积分转化为控制容积边界界面上的面积分，然后通过对界面 上的参数的近似而得到最终的离散方程。 对界面上的有关参数