控制方程和边界条件-Read.docVIP

第四章作业 摘要：针对二维电缆导热问题进行了数值计算，采用通用离散形式，加入了初始化条件以提高收敛速度。理论上证明了题目中电流500A是不合理的，若要收敛，电流不得大于496A。计算了不同电流下电缆的温度分布，通过网格考核和结果分析认为数值计算结果是可靠的。 导热问题 如图所示铜芯电缆，电流为500A，内径为10mm，外包材料为聚氯乙烯的厚度为2mm，导热系数为0.15W/(m?K)。电缆左半边为绝热边界条件，右半边为第三类边界条件，空气温度为20℃，绝缘层表面与环境间的符合表面传热系数为10W/(m?K)。铜的电阻率Rt=R01+at-20，R0=1.75×10 要求： （1）使用通用离散形式离散控制方程 （2）通过ADI TDMA算法求解 （3）完成报告（包括控制方程和边界条件、离散过程、边界条件处理、模型验证、网格独立性考核及结果分析）。 1.控制方程和边界条件 1.1 控制方程 一般变截面二维稳态导热问题控制方程的通用形式可表示为 1 该问题宜采用极坐标，形式如下 1 1.2 边界条件 该计算域关于x轴对称，取一半区域计算。由于圆心为奇异点，在中心位置挖去一个半径为mo的小圆面，将该边界设为绝热。 边界条件如下 r=R θ=0 2 方程离散 2.1 通用控制方程离散 鉴于附加源项法的优越性，计算采用内节点法，方程离散采用控制容积积分法，假定控制容积界面上热流密度是均匀的，引入东西尺度系数SX和名义半径R。 扩散项： tt+?ts 源项： t 整理上述结果： a a 2.2 三种坐标系下系数通用表达式 坐标系 直角 圆柱轴对称 极坐标 通用表达式 东西坐标 x x θ X 南北坐标 y r R Y 半径 1 r r R 东西尺度系数 1 1 r SX 3 边界条件处理 边界条件处理采用附加元项法。 3.1 径向边界 r=R （1）r=R，π2≤θ≤π a a （2）r=R，0≤θ≤ a q 整理得： [ a （3）r=mo边界绝热，则T a a 3.2 对称边界 θ=0 （1）θ=π边界绝热，则T a a （2）θ=0边界绝热，则T a a 3.3 内热源 电流通过铜芯产生热量，单位长度的热流量 Q= 电缆横截面上，单位面积内热源表示为 S= 可得 S 4. 题目错误的理论证明与初始化条件的设置 4.1 题目错误的理论证明 分析电缆导热的物理本质，由能量守恒定律，当导热达到稳态时，内热源产生的热量与边界的散热相同。该物理模型中只有电缆右侧边界向外散热，其它边界均为绝热。 单位长度上，内热源产生热量 Q= 内热源产生热量全部从右侧边界流出，热流密度为 q= 第三类边界条件热流密度公式 q 内部边界向外散热时温度逐渐降低，而铜的导热系数很大，近似认为铜线内部温度相同。散热边界温度应低于铜芯温度，设温差为?0，边界温度为T，则铜芯温度为T+?。 由能量守恒 h 整理得 T 方程右端恒大于0，方程左端温度T应大于20℃。若取电流为500A h 方程左端小于0，方程无解，不符合能量守恒定律。不考虑求解的T值是否符合物理学规律，若满足能量守恒定律，必有 h 解得I496A，即电流不得大于496A。 4.2 初始化条件的设置 计算过程中，需要在网格节点上采用ADI TDMA对离散方程进行迭代。提供一个合理的初场，对迭代的收敛速度非常重要。采用4.1中思想，考虑达到稳态时全场温差不会太大，为简化计算，假设计算域内温度相同。取?=0，可以得到全场初始温度 T= 5 网格考核 数值计算中网格划分采用内节点法，先划分界面，下面m×n为x和y方向的界面数。计算仅选择了上半部分，迭代收敛后将计算结果沿y轴镜像后输出。收敛条件为：（1）两次迭代之间全场温度的平均误差小于10-7；（2）若精度不能达到 收敛条件 δ= 括号内数字表示达到收敛所需迭代次数，计算过程中输出迭代次数与温度收敛情况。 题目中没有给出铜的导热系数，查阅相关文献，计算中取为400W/(m·K)。考虑在铜芯与外包层界面处导热系数有统一的计算方法，取此面为界面，通过挖除奇异点的大小mo与Y方向界面数目控制界面的位置。考虑包覆层外侧边界条件不同，不同边界条件采用界面分隔，因此X方向界面数取为奇数。 5.1 9×7网格计算结果（291240） 5.2 19×14网格计算结果（311089） 5.3 49×35网格计算结果（132） 5.4 99×70网格计算结果（401） 5.5 199×140网格计算结果（1090） 可见，网格划分越细，计算结果越精确。当网格划分至49×35时，计算结果与更为细密的网格相差不大，认为得到了网格无关解，且49×35网格迭代次数较少，因此下面验证计算中采用此网格。 6 模型验证 4.1中论证了题目中给出的电流时错误的，并给出了符合物理规律的电流范围。在此范