求助，使用Comsol计算聚合物沿厚度方向的浓度

hq61a

        按照文献《Photomechanics of light-activated polymers》中的公式(5)-(8)模拟聚合物在光照过程中沿聚合物厚度方向CR(concentration fields of free radicals)、CI(concentration fields of initiators)及I(light intensity)的变化，结果如文献中Fig.5、Fig.6所示。NA, h及v分别为 Avogadro数,Planck常数及入射光频率。
        我使用Comsol模拟这一过程，但结果与文献中的结果有出入，求助于好心人，望好心人能给予指导，谢谢。我的QQ为823641289。