[论文解读] On exact solutions and numerics for cold, shallow, and thermocoupled ice sheets

Ed Bueler, Jed Brown|arXiv (Cornell University)|Oct 16, 2006

Cryospheric studies and observations参考文献 10被引用 4

一句话总结

本文为冷、浅冰盖的热耦合冰盖模型提供了精确的解析解，使数值格式的严格测试成为可能。推导了使指定厚度和温度场满足耦合质量与能量平衡方程所需的补偿积累和热源函数，并提供了C语言实现及误差分析，通过温度平流格式中的不稳定性解释了EISMINT模拟中‘辐条’现象的成因。

ABSTRACT

This three section report can be regarded as an extended appendix to (Bueler, Brown, and Lingle 2006). First we give the detailed construction of an exact solution to a standard continuum model of a cold, shallow, and thermocoupled ice sheet. The construction is by calculation of compensatory accumulation and heat source functions which make a chosen pair of functions for thickness and temperature into exact solutions of the coupled system. The solution we construct here is ``TestG'' in (Bueler and others, 2006) and the steady state solution ``Test F'' is a special case. In the second section we give a reference C implementation of these exact solutions. In the last section we give an error analysis of a finite difference scheme for the temperature equation in the thermocoupled model. The error analysis gives three results, first the correct form of the Courant-Friedrichs-Lewy (CFL) condition for stability of the advection scheme, second an equation for error growth which contributes to understanding the famous ``spokes'' of (Payne and others, 2000), and third a convergence theorem under stringent fixed geometry and smoothness assumptions.

研究动机与目标

• 为在冷、浅冰近似下满足热耦合冰盖模型构建精确解析解。
• 通过提供具有指定厚度和温度场的已知解，实现对冰盖模型的严格数值验证。
• 分析有限差分格式在温度方程中的稳定性和收敛性，特别关注CFL条件和误差增长。
• 通过将‘辐条’（EISMINT II模拟中的伪影）与应变加热对温度的局部Lipschitz常数联系起来，解释其成因。
• 提供精确解的参考C语言实现，用于模型验证与测试。

提出的方法

• 通过指定时间与径向相关的厚度 $ H(t,r) = H_s(r) + \phi(r)\gamma(t) $ 构造解，其中 $ \phi(r) $ 为余弦形扰动。
• 假设温度分布 $ T(t,r,z) = T_s(r) \frac{\nu(t,r) + H(t,r)}{\nu(t,r) + z} $，其中 $ \nu $ 由 $ H $ 和 $ T_s $ 推导得出，以满足边界条件。
• 利用具有整数 $ n $ 的Glen流变定律，通过积分恒等式 $ \int e^x x^n dx = p_n(x) e^x $（$ p_n $ 为递归定义的多项式）解析计算速度场 $ \mathbf{U} $。
• 通过将精确场代入控制PDE并求解所需强迫项，推导出必要的源项——补偿积累和热源——$ \Sigma_c $。
• 分析了温度方程的有限差分格式，推导出正确的CFL条件及包含局部应变加热梯度的误差增长不等式。
• C语言实现提供了精确解的参考代码，包括在网格点上对所有场和源项的计算。

实验结果

研究问题

• RQ1 使指定厚度和温度场成为热耦合冰盖方程精确解的补偿积累和热源函数的显式形式是什么？
• RQ2 应变加热对温度的局部Lipschitz常数如何影响有限差分格式中的误差增长？
• RQ3 为何‘辐条’——EISMINT II模拟中的伪影——在较暖冰区出现？这是否可与流动定律的温度依赖性相关联？
• RQ4 此模型中温度方程迎风格式的正确CFL条件为何？
• RQ5 有限差分格式的误差如何随时间增长？在何种条件下可保证收敛？

主要发现

• 通过推导出所需的源项 $ \Sigma_c $ 和 $ M $，构造了精确解‘Test G’，使得给定的 $ H(t,r) $ 和 $ T(t,r,z) $ 满足耦合质量与能量平衡方程。
• 迎风格式的正确CFL条件被推导为 $ \Delta t \leq \min_{ijkl} \left\{ \frac{\Delta z}{|w_{ijkl}|} \right\}^{-1} $，其依赖于垂直速度的大小。
• 格式中的误差增长受一个局部不等式控制，涉及时间步长、局部应变加热和截断误差，其中前者在误差放大中起主导作用。
• 局部Lipschitz常数 $ L_{ijkl} = \Sigma_{ijkl} \frac{Q}{R T_{ijkl}^2} $ 解释了EISMINT II模拟中‘辐条’的出现，因其在从冷冰过渡到暖冰时急剧增加。
• 在严格假设下收敛定理成立：几何形状固定且精确解光滑，确保当 $ \Delta x \to 0 $ 时，有限差分格式收敛。
• C语言实现提供了精确解的完整参考实现，使数值模型可基于已知解析基准进行验证。