本文详细介绍了三维非稳态导热方程使用Brian ADI格式结合TDMA方法的求解过程，包括沿x、y、z三个方向的迭代计算，并提供了一个具体的算例来说明算法的应用。这种方法适用于处理对时间步长限制严格的三维导热问题，能有效提高计算效率。
基本原理
在上一篇文章，讲述了二维导热问题采用交替方向隐式方法[1]的求解过程。对于二维问题，Peaceman-Rachford ADI格式是无条件稳定的。然而，该格式对于三维问题是条件稳定的，其稳定条件为：
a Δ t [ 1 Δ x 2 + 1 Δ y 2 + 1 Δ z 2 ] ≤ 1.5 a\Delta t[\frac{1}{\Delta x^2}+\frac{1}{\Delta y^2}+\frac{1}{\Delta z^2}]\le 1.5 aΔt[Δx21+Δy21+Δz21]≤1.5
可见该格式对时间步长 Δ t \Delta t Δt的限制很严格。相比之下，Brian ADI格式则是绝对稳定的，故可以选取较大的时间步长，同时每一时间步的计算量较小，可大幅加快计算速度。
下面介绍Brian ADI格式下三维非稳态导热方程的求解过程，求解过程中涉及两个温度中间量： U U U和 V V V。
基于有限体积法的全隐离散方程可表示为：
a P 0 ( T P − T P 0 ) = a E ( T E − T P ) − a W ( T P − T W ) + a F ( T F − T P ) − a B ( T P − T B ) + a N ( T N − T P ) − a S ( T P − T S ) + S Δ x Δ y Δ z a_P^0(T_P-T_P^0)=a_E(T_E-T_P)-a_W(T_P-T_W)+a_F(T_F-T_P)-a_B(T_P-T_B)+a_N(T_N-T_P)-a_S(T_P-T_S)+S\Delta x\Delta y\Delta z aP0(TP−TP0)=aE(TE−TP)−aW(TP−TW)+aF(TF−TP)−aB(TP−TB)+aN(TN−TP)−aS(TP−TS)+SΔxΔyΔz
(1)沿 x x x方向求解：
先取 Δ t / 2 \Delta t/2 Δt/2的时间步长，对x方向做隐式处理，y、z方向取温度场 T 0 T^0 T0的值，求解温度场的中间量 U U U，则全隐离散方程可改成：
a P 0 ( U P − T P 0 ) = a E ( U E − U P ) − a W ( U P − U W ) + a F ( T F 0 − T P 0 ) − a B ( T P 0 − T B 0 ) + a N ( T N 0 − T P 0 ) − a S ( T P 0 − T S 0 ) + S Δ x Δ y Δ z a_P^0(U_P-T_P^0)=a_E(U_E-U_P)-a_W(U_P-U_W)+a_F(T_F^0-T_P^0)-a_B(T_P^0-T_B^0)+a_N(T_N^0-T_P^0)-a_S(T_P^0-T_S^0)+S\Delta x\Delta y\Delta z aP0(UP−TP0)=aE(UE−UP)−aW(UP−UW)+aF(TF0−TP0)−aB(TP0−TB0)+aN(TN0−TP0)−aS(TP0−TS0)+SΔxΔyΔz
整理后得到：
( a P 0 + a E + a W ) U P = a E U E + a W U W + ( a P 0 − a F − a B − a N − a S ) T P 0 + a F T F 0 + a B T B 0 + a N T N 0 + a S T S 0 + S Δ x Δ y Δ z (a_P^0+a_E+a_W)U_P=a_EU_E+a_WU_W+(a_P^0-a_F-a_B-a_N-a_S)T_P^0+a_FT_F^0+a_BT_B^0+a_NT_N^0+a_ST_S^0+S\Delta x\Delta y\Delta z (aP0+aE+
最低0.47元/天 解锁文章
5561
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
