第3章-导热-3非稳态.ppt

* 井筒周围地层的非稳态导热数学模型与求解 初始条件： t0z该段中间位置处原始地层温度： 式中，ts为地表处不受环境影响的温度，℃； m为地温梯度，℃/m； z为该段中间位置处的坐标。 非稳态导热的分析与计算 井筒周围地层的非稳态导热数学模型与求解 边界条件： 在地层和井筒的交界面r=R处： 径向无限远处： 式中，?为地层导热系数； R为交界面处半径； qlz为该段内单位长度井筒通过界面传给地层的热量，W/m； t0z为该段中间位置处的原始地层温度。 非稳态导热的分析与计算 井筒周围地层的非稳态导热数学模型与求解 完整数学模型： 初始条件： 边界条件： 采用Laplace变换可以解得模型的分析解。 非稳态导热的分析与计算 井筒周围地层的非稳态导热数学模型与求解 地层内的温度分布： 其中 式中，J0、J1为零阶、一阶第一类贝塞尔函数； Y0、Y1为零阶、一阶第二类贝塞尔函数； 数学上常用的特殊函数 。 非稳态导热的分析与计算 井筒周围地层的非稳态导热数学模型与求解 地层和井筒交界面处的温度： 式中 关键是如何计算其中积分I。 非稳态导热的分析与计算 解析法很难求出I，而数值积分工作量又很大。 通常利用解析解数据得到简单拟合关系式。 井筒周围地层的非稳态导热数学模型与求解 非稳态导热的分析与计算 无量纲温度： 无量纲时间： 处理地层内非稳态导热的半解析法被国外热采界广泛采用。 Hasan公式 井筒周围地层的非稳态导热数学模型与求解 井筒热损失： 非稳态导热的分析与计算 地层热阻： 单位井筒长度通过地层向无限远处传热的导热热阻。 注意：地层热阻不是常数，与过程的持续时间有关。 过程持续时间越长，地层热阻就越大。 根据过程的持续时间，可得地层热阻的大小; 由于交界面的温度未知，必须结合井筒内的热量传递过程，才可得到井筒内流体和地层间的传热量。 井筒周围地层的非稳态导热／算例 非稳态导热的分析与计算 井筒周围地层的非稳态导热／算例 井筒周围地层的非稳态导热／算例 井筒周围地层的非稳态导热／算例 井筒周围地层的非稳态导热／算例 井筒周围地层的非稳态导热／算例 1.50 井筒周围地层的非稳态导热／算例 井筒周围地层的非稳态导热／算例 井筒周围地层的非稳态导热／算例 井筒周围地层的非稳态导热／算例 1.50 井筒周围地层的非稳态导热／算例 井筒周围地层的非稳态导热／算例 井筒周围地层的非稳态导热／算例 井筒周围地层的非稳态导热／算例 井筒周围地层的非稳态导热／算例 * * * * * * * * * * * * * * 当Bi→0时，可以忽略温度沿空间的变化。 认为物体内任意位置的温度在同一时刻均处于同一温度下； 像该物体原来连续的质量与热容量都汇总到一点上，只有一个温度值那样。 3.4 集总参数分析法 非稳态导热的分析与计算 集总热容系统或称“零维问题” 集总参数法——忽略物体内部导热热阻的简化分析方法 特征：忽略物体温度沿空间的变化（零维） 工程背景： ——金属工件热处理工艺中的加热或冷却过程 ——测温热电偶接点的加热和冷却过程 3.4.1 集总热容系统的温度场 设有任意形状固体 —体积为V，表面积为A —密度ρ、导热系数?和比热容c均为常数 —初始时刻具有均匀的温度t0 —在?＝0时突然将其置于温度为t∞的流体中，流体与物体表面间的表面传热系数为h，且保持不变。假设t0 t∞。 试分析物体的温度随时间的变化规律。 3.4 集总参数分析法 物理模型 实际导热体：?很大，或尺寸l很小，或h极低 导热热阻l /? 1/h 对流热阻，可忽略导热热阻。 Bi较小可看成集总热容系统—— 物理模型。 建立控制方程可采用两种方法： ——基于导热微分方程的方法。忽略了温度沿空间的变化后，可以舍去导热微分方程中的扩散项，将物体与环境间的换热量视为虚拟热源，由导热微分方程可得到集总热容系统的控制方程。 ——基于能量守恒的方法。物理概念清楚、容易理解。 数学模型 3.4.1 集总热容系统的温度场 3.4 集总参数分析法 取整个物体作为控制容积： 无内热源时，物体冷却过程中应满足的关系： 物体表面散热量全部来自于其热力学能的减少。 数学模型 3.4.1 集总热容系统的温度场 3.4 集总参数分析法 初始条件： 非齐次方程 数学描述 3.4.1 集总热容系统的温度场 3.4 集总参数分析法 分离变量、积分： 自然指数函数 集总热容系统的温度分布： 齐次方程 特点：在过程开始阶段，温度变