双探针法测量月壤导热系数正问题研究.pdfVIP

中国T程热物理学会 传热传质学 学术会议论文 编号：073258 双探针法测量月壤导热系数的正问题研究 孙创岳贤德夏新林 哈尔滨工业大学能源科学与工程学院，哈尔滨15000l 摘要：在考虑月表辐射换热的前提下，数值模拟了由探针和月壤组成的多层介质内二维非稳态传热 过程，通过反演得到的月壤温度响应曲线可求得月壤的导热系数。本文主要探讨加热探针的长度、直 径和加热功率等因素对测量结果的影响，并初步确定了探针间距。受月壤深度方向上温度梯度的作用， 分析了不同初温下响应探针的温度扰动情况。模拟结果表明，由于月壤弱导热性的影响，当探针间距 13mm，加热功率为O．03W抽m时，lOmin的测量时闻内产生8K的温度响应。 关键词：双探针法；月壤；导热系数；数值模拟 1．引言 了解月壤热物性参数和月表辐射特性可为月球热环境的研究提供数据支持，是月球 探测系列科研项目的组成部分。探针法是目前比较常用的物性测量方法之一，具有快速、 简便等特点。原理为通过对被测介质的线热源瞬态加热，获得相应的“温度一时间曲线”， 进而反演得到该介质的热物性参数。已有报道表明，NASA的火星探测计划中已采用探 针法进行火星表面土壤的热物性测量。目前，探针的绝大部分应用还是在地球表面上进 Naidu 行，Dali L．Bristow【23 J等利用多探针对±壤的热物性、含水量等进行分析。谢华清，孟凡凤【4’5l 等则分别采用探针法进行了介质导热系数的测定。 受月壤自身条件约束，采用探针法测量导热系数有其特殊性。由于月昼时间较长， 月表的温度变化明显，应保证探针有良好的适应性；上述文献多采用一维模型进行计算， 而月壤在其深度方向上存在一定温度梯度，为了尽量减小计算误差，本文采用了二维计 算模型。月壤本身为弱导热性t6】对热量的传递产生了抑制作用，导致温度响应具有延迟 性。需保证探针的间距和加热功率才能得到可信的温度响应，但探针的加热功率却同时 受月球车等外部因素制约。因此，合理的设计探针各参数是测量精确的前提保证。 2．计算模型与控制方程 本文数值模拟了由月壤、探针套管及填充材料 组成的多层介质内非稳态导热过程，并考虑了月表 和外界的辐射换热作用。探针插入月壤后即开始加 热，分析不同时刻月壤各处的温度变化情况。 图1为计算区域非均匀网格离散示意图，通过 给出不同位置处不同时刻的温度值，分析探针和月 壤的温度变化情况。测温点的坐标为：半径方向上， Fl。5mm，l。7mm，13mm，18mm，30mm；深度方 向上，x20mm，10mm，30mm，50mm。图中g，。代 表太阳辐射热流，印，为月表自身辐射。 数为常数。建立了圆柱坐标系。卜．二维非稳态多层介数为霉羹多薹皇拿蓄毳鍪架柔宰!墓j 是冀誓墨霁 圊H叫佾划刀小思崮图·网格划分示意图 质的导热模型： 36 K【u，∞ ，糙【u，00 ㈣ 1 万2口 ，’【矿+7石+丽J》加，悟弓詈+引吲啉州岫， 式中，口 ，． 为不同半径处相应材料的热扩散率。 2 m^o ：』五：x∈【叫册d，∈胁】 ‘2’ 丁‘x，，，o’21三 工 ：x∈；fl，，： 。r，．∈ _：j，：； 式 2 为计算初温的选取。瓦 x 为月壤的温度分布，互为探针初始温度值。 I， L J I兰 o：，．∈【o，，i】 f o，x o 出。， 3 、 l一五娑 呸舔～掰r4 o，，，f ：，∈ 巧，oo 式 3 为计算区域上表面的边界条件，探针表面处理为绝热面，月壤与外界存在热交换。 咒为月壤的导热系数，吼为太阳辐射常数，呸为月壤的太阳吸收率，s为月壤红外发 射率，仃为stef抽．Boltzmann常数，r O，r，f 为月表温度值。 z o，， o，x∈ o，‘ ：一^÷ go 4 式 4 表示探针中心处的边界条件，鼋b是探针加热功率，^为探针填充材料的导热系数。 剩余二个边界处理为绝热面，表达式为： II∞ 0 X r 0 5 |I ∞ 0 r 订一缸订一跏 采用物性调和平均法求解不同材料间界面的导热系数： 竺竖：竺墼+箪 61 。 筏。J 筏“ ^?J 对式 1 采用全隐离散方式得到式 7 ，．并对 7 式采用ADI算法求解。 7 口0正≥ 叫：：：，z：u+口0I!。，／+口荔z：+。+匆≥正：：，一。十缆z箩1 3．计算结果及分析 模拟计算中各参数之间是相互联系、相互制约的，在对某一参数进行分析时，通常 假设其它计算条件不变。选取了一个月球昼夜内月表温度最高时的月壤温度作为初值， 一 3．1探针几何尺寸的选择 探针长度的选择与加热功率有关，在加热功率不变的情况下，探针的长度越长，探 针的热容量越大，加热探针表面温升越小，直接影响月壤中温度响应情况。