计算传热学第2节-第1苏章有限体积法.pptVIP

asdf Sun Jining 2008 @ BUAA * 将整个求解域划分为n个立方体区域，从t1到t2时刻，每立方体能量守恒方程： ((ρcT)Pt2-(ρcT)Pt1)?x?y?z=((λ(?T/?x))e-(λ(?T/?x))w)?y?z?t+SP?x?y?z?t 1 有限体积法 有限体积方法的基本思想 体平均量 每时间步n个未知数 面时平均量 每时间步n-1个未知数 体时平均量 每时间步n个未知数 x z y t2时刻 t1时刻 面时平均量 假设在空间上分段线性分布，在时间上阶梯分布，则温度梯度的面时平均量等于上一时间步两侧中心点的温度差分 QT=((λ(?T/?x))e-(λ(?T/?x))w)?y?z?t ≈(λ(?T/?x))et1-(λ(?T/?x))wt1)?y?z?t ≈(λe(TEt1-TPt1)/δxe-λw(TPt1-TWt1)/δxw) ?y?z?t δxe P E W δxw e w T x P E W e w ?xP T t t2 t1 t3 asdf Sun Jining 2008 @ BUAA * 将整个求解域划分为n个立方体区域，从t1到t2时刻，每立方体能量守恒方程： ((ρcT)Pt2-(ρcT)Pt1)?x?y?z=((λ(?T/?x))e-(λ(?T/?x))w)?y?z?t+SP?x?y?z?t 1 有限体积法 有限体积方法的基本思想 体平均量 每时间步n个未知数 面时平均量 每时间步n-1个未知数 体时平均量 每时间步n个未知数 x z y t2时刻 t1时刻 体时平均量 在空间上一般为已知函数，假设在时间上阶梯分布，则发热功率的体时平均量等于上一时间步的发热功率体平均量 ST=SP?x?y?z?t ≈SPt1?x?y?z?t P E W δxe δxw e w ?xP asdf Sun Jining 2008 @ BUAA * 将整个求解域划分为n个立方体区域，从t1到t2时刻，每立方体能量守恒方程： ((ρcT)Pt2-(ρcT)Pt1)?x?y?z=((λ(?T/?x))e-(λ(?T/?x))w)?y?z?t+SP?x?y?z?t 1 有限体积法 有限体积方法的基本思想 体平均量 每时间步n个未知数 面时平均量 每时间步n-1个未知数 体时平均量 每时间步n个未知数 每个立方体的有限体积方程： ΔUP≈ρPcP(TPt2-TPt1)?x?y?z QT≈(λe(TEt1-TPt1)/δxe-λw(TPt1-TWt1)/δxw)?y?z?t ST≈SPt1?x?y?z?t ρPcP(TPt2-TPt1)?x?y?z≈ (λe(TEt1-TPt1)/δxe-λw(TPt1-TWt1)/δxw)?y?z?t+ SPt1?x?y?z?t ρPcP(TPt2-TPt1)/?t≈ (λe(TEt1-TPt1)/δxe-λw(TPt1-TWt1)/δxw)?x+ SPt1 asdf Sun Jining 2008 @ BUAA * 将整个求解域划分为n个立方体区域，从t1到t2时刻，每立方体能量守恒方程： ((ρcT)Pt2-(ρcT)Pt1)?x?y?z=((λ(?T/?x))e-(λ(?T/?x))w)?y?z?t+SP?x?y?z?t 1 有限体积法 有限体积方法的基本思想 体平均量 每时间步n个未知数 面时平均量 每时间步n-1个未知数 体时平均量 每时间步n个未知数 每个立方体的有限体积方程： ΔUP≈ρPcP(TPt2-TPt1)?x?y?z QT≈(λe(TEt1-TPt1)/δxe-λw(TPt1-TWt1)/δxw)?y?z?t ST≈SPt1?x?y?z?t ρPcP(TPt2-TPt1)?x?y?z≈ (λe(TEt1-TPt1)/δxe-λw(TPt1-TWt1)/δxw)?y?z?t+ SPt1?x?y?z?t ρPcP(TPt2-TPt1)/?t≈ (λe(TEt1-TPt1)/δxe-λw(TPt1-TWt1)/δxw)?x+ SPt1 这种插值处理方法虽然解决了独立方程数目不够的问题，但同时也带来了数值误差！ asdf Sun Jining 2008 @ BUAA * 将整个求解域划分为n个立方体区域，从t1到t2时刻，每立方体能量守恒方程： ((ρcT)Pt2-(ρcT)Pt1)?x?y?z=((λ(?T/?x))e-(λ(?T/?x))w)?y?z?t+SP?x?y?z?t 1 有限体积法 有限体积方法的基本思想 体平均量 每时间步n个未知数 面时平均量 每时间步n-1个未知数 体时平均量 每时间步n个未知数 每个立方体的有限体积方程： ΔUP≈ρPcP(TPt2-TPt1)?x?y