冻结土壤温度场数值模拟的改进.pdfVIP

_______________________________________________________________________________ 冻结土壤温度场数值模拟的改进* 商翔宇 周国庆 别小勇 （221008 江苏徐州 中国矿业大学建筑工程学院） 摘 要：针对土壤冻结过程的数值分析，基于显热容法用有限差分法构造出热传导方程变空间步长的半隐 格式和全隐格式。结合实例分析了FEM 和FDM 方法造成的‘相变遗漏’，提出一种自调节时间步长、半 隐和全隐格式交替使用的方法，使模型求解迭代大为简化，有效避免了‘相变遗漏’引起的温度场分析误 差。该方法能够精确求解冻土温度场，是对传统方法的重要改进，也为精确分析土体冻结及冻胀融沉过程 创造了条件。 关键词：土壤冻结 数值分析 相变遗漏 土壤冻结是一种常见的现象。人工地层冻结技术在土木工程建设中也越来越发挥重要的 作用。要分析天然冻土和人工冻土的冻结工程，必须准确确定土体的温度场。常见的冻胀模 型如Konrad 的分凝冰模型以及Daiyu.Wang 的冻胀模型在计算冻胀量之前都要求准确求出冻 土的温度场。然而，在用有限差分法求温度场过程中，由于时间步长取法不当很容易造成‘相 变遗漏’。本文提出一种自调节时间步长、半隐和全隐格式交替使用的方法，该方法能够有 效防止分析过程中相变遗漏，并且计算速度较之固定时间步长方法有了很大的提高。 1 数学模型 常见的一维土壤冻结过程热传导方程为： ∂T ∂ ∂T ∂θ C L i （1） V (λ ) + ρi t z z ∂t ∂ ∂ ∂ 3 式中： CV ，λ—土壤体积热容量（J/m *℃）和热导率（J/m*s* ℃）； L ，T—融化潜热（J/m3 ）和温度（℃）； 3 3 3 Ρi ，θ—冰密度（kg/m ）和冰含量（m /m ）； i Z ，t—空间坐标（m ）和时间（s ）。 求解上述方程的一种有效方法是显热容法，即在发生相变的一个小的温度范围之内构造 [2] 比热函数，以代替以上方程的右端第二项，只以温度为待求函数 ： ∂T ∂ ∂T * * 　　　　　　　　　　　　　　　C V (λ ) 　　　　　　　　　（2 ） ∂t ∂z ∂z *基金项目：国家973项目（2002CB412704 ）、中国科学院知识创新重大项目（KZCX1-SW-04 ）、中国科学 院冻土工程国家重点实验室开放基金（SKLFSE200304）以及教育部优秀青年教师奖励计划、江苏省“青蓝 工程”新世纪学术带头人专项基金联合资助 作者简介：商翔宇（1977—），男，河南洛阳人，2001年毕业于中国矿业