1　实验方案　Ⅰ　为了比较　3　种自适应方法的实际效果　,在算法　2　的基础上　,设计了如下几组数值实验.　实验　1　:在固定网格　mesh0　上进行计算.　实验　2　:在固定网格　mesh1　上进行计算.　实验　3　:以mesh0　为T0　网格　,采用　Hess　网格自适应.　实验　4　:以mesh0　为T0　网格　,采用　Grad　网格自适应.　实验　5　:以mesh0　为T0　网格　,采用　KGrad　网格自适应.　考虑到方程(1)　的特殊性　,网格变化幅度不宜过大　,因此　,在自适应过程中对网格的规模进行了控制　,实验3、　实验　4　和实验　5　的网格单元数在整个计算过程中都没有超过3　600.　图3记录了实验　3　在t=3108　和t=11185　时刻的光子温度等值线　,其中等值线密集的区域为光子的突　变区.　图3中的等值线呈现出球对称性　,这种现象与模型问题的球对称特性是相适应的.　在实验中　,我们也观　察到实验　4、　实验　5　与实验　3　具有一致的光子温度等值线.　为了比较不同自适应方法下的网格变化情况　,截取了实验　3、　实验　4　和实验　5　两个时刻的网格图　,其结果　见图　4、　图5、　图6.　从这些图形可以看出　:　1)　实验　3　中　,加密产生的新网格落在等值线密集区域的两边.　在图　4(a)　中　,可以观察到等值线密集区域　前的网格被加密　;在图　4(b)　中　,等值线密集在　Ω2　区域　,被加密的却是介质交界面附近的　Ω1　和Ω3　中的单　元.　这说明　,利用　Hess　网格自适应技术　,被加密区域是梯度变化大的区域　,而不完全是梯度值大的区域　,从图　像来看　,加密区域落在了突变区域与正常区域的衔接地带.　3　2　第1期江军等　:求解二维三温辐射热传导方程组的一种网络自适应方法　©　1994-2007　China　Academic　Journal　Electronic　Publishing　House.　All　rights　reserved.　http://www.cnki.net　图3实验　3　的光子温度等值线图　Fig.　3　Photon　temperature　contour　in　Exp.　3　2)　实验　4　中　,梯度值大的区域得到了加密.　图5表现了网格的加密区域与等值线密集区域的一致性.　3)　实验　5　中　,网格加密