寒区水库冰盖厚度增长数值模拟研究1.pdfVIP

第32卷第1期 渤海大学学报(自然科学版) V01．32No．1 Bohai of Science M盯．20ll 2011年3月 Journal University(NaturalEdition) 寒区水库冰盖厚度增长数值模拟研究 冯景山，白乙拉，李 冰 (渤海大学数理学院，辽宁锦州121013) 摘要：为了研究寒区水库冰盖厚度的增长情况，分析了前人对寒区水库冰盖热力模式 及冰厚增长数值模拟研究的方法，对一维热传导方程和一维单相Steften问题所采用的数值解 法进行了改进，采用隐式的Crank—Nicolson格式进行离散。精度提升到二阶精度。且具有无条 件稳定性。并根据黑龙江省大庆红旗泡水库冰盖的温度、厚度等观测数据，进行了厚度数值模 拟。模拟结果表明：所改进的数值模拟方法及其算法是正确的、有效的。 关键词：冰盖；厚度；数值模拟 中图分类号：0241．82；P338．4 0 引 言 我国北方大部分地区的河流、湖泊、水库冬季常形成冰盖。冰盖的出现改变了水流的水利条件、热力 条件和几何边界条件，继而可能导致冰冻期或开河期的冰塞，甚至产生凌汛等灾害，因此有关冰盏的生 长消融、冰塞的形成等问题的数学模型与数值模拟的研究，一直受到有关学者们的广泛关注。淡水冰与 环境和气候的变化是息息相关的。无论是作为研究淡水冰对气候的影响还是研究寒区水库冰盖生消以 及其本身的热力学特征，前人对寒区水库冰盖热力模式及冰厚增长的数值模拟等研究已有了许多成 果¨。6，。本文主要考虑已知气温影响下冰面温度的条件下。根据寒区水库冰盖表面温度的实测数据．建 立了一维热传导方程作为模型方程，对此自由边界条件的一维抛物型方程的初边值问题(冰盖底郡边界 用非等距步长的方法来处理，数值模拟寒区水库冰盖厚度的增长变化情况。 1 模型的建立 北方地区水库冬季很少有供水任务，因而水库冰盖底部水流缓慢．由机械能转化的热能部分较少。 对冰盖生长影响不大，所以寒冷地区水库冰盖生长主要与气温有关．当气温降到0℃以下时．水由于失 热致使水面开始结冰。随气温不断地降低．为了达到热能平衡，冰盖下的水不断失热由水变成冰，考虑到 冰内温度在垂直方向的梯度变化远远大于水平方向的变化。即热传导主要在垂直方向上进行．因此忽略 水平方向上变化后。寒区水库冰盏热力学模型可由一维热传导方程描述。 收穑日期：2010—12-09． 作者简介：冯景山(1985一)．男．渤海大学硕士研究生．从事应用教学研究工作l 自乙拉(1961一)．男．博士．教授．从事应用教学研究T作． 万方数据 6 渤海大学学报(自然科学版) 第32卷 设冰的比热为c，密度为P，热传导系数为K，取冰表面一点处为坐标原点O，过原点D垂直向下 为坐标轴戈的正向，r(髫，t)表示t时刻算处的温度，冰水交界面为茹=S(t)。由一维热传导方程和一维 单相Stefien问题‘7·8’可得： 冰域方程： 面OT跏窘口_WCp’o菇s㈤p0。 (1) 冰面温度： T(0，t)=T(t)，t20，(r(t)为已知)。 (2) 冰水交界面条件： 罢一Q=A石dS，A等，Q=量，茗=s(t)，‘≥o。 (3) 初始冰厚： s(o)=So。 (4) r(S(t)，1)=0。 (5) 传导系数(w／m·℃)，￡为单位质量冰的熔解潜热(