基于有限差分法的冻结过程数值分析.docxVIP

基于有限差分法的冻结过程数值分析 在中国东北部和西北部的寒冷地区，土木工程方兴不断，人工土壤冷冻技术在土木工程施工中发挥着越来越重要的作用。要分析自然冻土和人工冻土的冰冻工程，必须正确规定土壤温度变化规律。例如，konrad分凝冰模型和daiyuwa的冻胀模型在计算冻胀前必须正确计算冻面的温度场。然而，在使用有限差分法求解温度场的过程中，由于时间间隔不正确，很容易导致“异常日期”。在这项工作中，我们提出了一种基于自适应时间间隔、半隐藏和全隐藏格式的交替使用方法。这种方法可以有效。 防止传统计算过程中出现的相变遗漏, 并且计算速度较之固定时间步长方法有了很大的提高. 1 显热容法求解 常见的一维土壤冻结过程热传导方程为 , (1) 式中:CV,λ为土壤体积热容量 (J/ (m3·K) ) 和热导率 (J/ (m·s·K) ) ;L,T为融化潜热 (J/m3) 和温度 (℃) ;ρi,θi为冰密度 (kg/m3) 和冰含量 (m3/m3) ;z,t为空间坐标 (m) 和时间 (s) . 求解上述方程的一种有效方法是显热容法, 即在发生相变的一个小的温度范围之内构造比热函数, 以代替以上方程的右端第二项, 只以温度为待求函数 , (2) 其中: 式中:C*V为等效热容;λ*为等效导热系数;Cf,λf分别为已冻土的体积热容和导热系数;Cu,λu分别为未冻土的体积热容和导热系数;Tf-ΔT≤T≤Tf+ΔT为假设小的相变温度范围. 2 相变影响下的半隐格式 方程 (2) 是与时间和空间都有关系的发展方程, 采用有限差分法能对其有效求解.在常用的有限差分格式中, 显格式精度最差, 隐格式次之, 而Crank-Nicolson格式精度较高;但是从解的稳定性角度分析, 则隐格式的稳定性最好, 因此本文采用隐格式进行求解.根据方程系数的取法不同, 隐格式又有全隐格式和半隐格式之分.全隐格式公式如下 式中j,n分别为空间和时间坐标标记. 所谓半隐格式是将式 (3) 中系数C*v和λ*的时间步n+1换成n, 通过半隐格式只需求解线性方程组即可.事实上, 半隐格式和全隐格式对本方程的求解精度影响不大可忽略不计.作者也分别编制了两种不同格式的求解程序进行了验证. 由式 (2) , (3) 可见, 相变的影响在模拟中是通过在假设相变范围的等效热容来实现的, 因此要全面考虑相变的影响, 必须保证在每一计算时间步内所有发生相变的空间单元的等效热容能够包含相变潜热项.而土体空间单元的等效热容是否应该包含相变潜热项, 则是通过该单元的温度值来判断的.由式 (2) 知, 只有其温度落在Tf-ΔT≤T≤Tf+ΔT范围之内, 该单元的相变影响才能被考虑.这样, 就不免会出现‘相变遗漏’的情况.比如空间某单元在某一时间步之内温度由T1Tf+ΔT变为T2Tf-ΔT, 显然该单元在这一时间步内发生了相变, 但是根据T1或T2都无法判断出该单元发生相变这一事实, 因此单元等效热容也就没能包含本应包含大相变潜热项.由于半隐格式是根据时间步初的温度值、全隐格式是根据时间步末的温度值来判断等效热容是否包含相变潜热项的, 因此对于半隐和全隐格式, 一旦相应的温度判据即时间步初的温度值和时间步末的温度值没能落在Tf-ΔT≤T≤Tf+ΔT范围之内, 该单元的相变影响才能被考虑, 否则就会出现下述的‘相变遗漏’现象. 例如在任意一个计算时间步内, 假设计算空间上某一点在时间步初的温度值为T1, 在时间步末的温度值为T2.如果时间步长控制不好, 不同的差分格式就会有不同的‘相变遗漏’.采用半隐格式, 若T1Tf+ΔT而Tf-ΔT≤T2≤Tf+ΔT(时段初未冻而时段末处于过渡期) 则相变作用就会遗漏;采用全隐格式, 若Tf-ΔTT1Tf+ΔT而T2Tf-ΔT相变作用也会遗漏.但是如果T1Tf+ΔT且T2Tf-ΔT则不论采用半隐或是全隐格式相变作用都会被遗漏掉. 可见如果仅仅采用一种差分格式, 那么无论是半隐格式还是全隐格式, 对于本文所述问题都会不可避免地存在‘相变遗漏’这一事实, 因此是不正确和不可取的.由于冻结过程伴随着大量相变潜热的释放, 所以原来采用固定格式算法得出温度将会偏低, 而冻结锋面推进速度则会偏快. 3 速度速的数值分析步骤 为了避免数值分析中相变作用的遗漏, 必须有效地控制时间步长.本文根据半隐格式和全隐格式的不同对方程的求解精度影响不大这一特点, 提出了自调节时间步长、半隐和全隐格式交替使用的方法 (Modified FDM以下简称MFDM) , 其中半隐格式直接求解线性方程组即可, 全隐格式则要采用迭代法求解.由于本文提出的方法采用的是稳定的计算格式, 而且采用了自调节时间步长的技术, 所以该方法具有稳定和计算速度快的特点. 采用该方法进行数值分析的