瞬态辐射传输地最小二乘有限元法.pdfVIP

中国。l：程热物理学会 传热传质学 第十一届学术会议 编号：053291 瞬态辐射传输的最小二乘有限元法+ 安巍阮立明谈和平 (哈尔滨工业大学航空航天热物理研究所，黑龙江， 哈尔滨l50001) Tel：045l—mail：Ruanlm@hit．edu．cn 时间步长较大的情况下仍然可以得到光滑无振荡的解。而且，最小二乘法产生的求解系数矩阵是对称 正定的，与GFEM中的系数矩阵相比，仅需要存储一半的非零系数，可以应用许多高效的迭代求解方 法进行求解。为r检验模型，本文研究了一维和二维吸收散射性介质内瞬态辐射传输阿越，其结果与 复杂形体内的瞬态辐射传输。 关键词：微尺度，瞬态辐射传输．最小=乘有限冗法 1引言 众所周知，辐射是一种电磁波．它的传播速度是光速，在通常情况下．都是假设辐 射传输是瞬间达到稳态的。近期，随着超快激光(Ultra-FastLaser)的快速发展及微尺 度系统【1]的研究，半透明介质内的瞬态辐射传输受到了J“泛的关注。短脉冲激光对金属 的加工，医学成像技术，无损伤激光治疗技术，粒子的探测研究，红外遥感技术等许多 方面的进展应用表明：瞬态辐射的传输过程需要理论意义上的深刻理解和精确的数值模 拟。例如，短脉冲激光医学成像技术发展的关键之一就是要有一个有效而准确的前向计 算。如果没有一个好的前向计算，那么反向的分析就不可能提供高精度的结果。除了上 述的新兴技术的应用背景，在传统的辐射传输领域中也需要准确高效的瞬态辐射数值模 型。例如，在天体物理辐射问题的研究中，由于空间尺度的特性决定了必须考虑辐射传 输随时间的变化。 在处理超短脉冲激光引起的辐射传输模拟中，稳态辐射传递方程不可能提供准确的 结果。这是由于描述瞬态辐射波传递的微分方程已经变为了双曲型的偏微分方程，辐射 强度同时是时间，空间．方向，波长四种参数的函数。因此，瞬态辐射传递的模拟研究 成为近十年来光学和计算辐射学的研究热点。近十年中，已经有许多理论研究和数值方 出了瞬态辐射的积分方法，稍后，Tan[4]等也提出了稍有不同的积分方法。他们的方法 是建立在积分型辐射传递方程的基础上的一种拉格朗日观点的方法，具有精度高的特点， 但是这种方法在处理多维问题的各向异性散射时比较复杂。Guo[5]和他们的同事用 和有限体积法[7】已经成功的预测了辐射波随时间的传递。 有限元法作为一种成熟的数值计算方法在力学、传热学等领域有着广泛的应用。它 有着对复杂形体适应性强．在计算单元内近似程度高等优点。从2004年开始，吸收散射 限元，非连续有限元等逐渐被提出用于模拟辐射传输。但这些研究都局限丁二稳态的辐射 资助项Izl：国家自然科学基金资助fN 682 问题，瞬态辐射传递的有限元研究尚未见报道。 正如Tan[4】指出的，准确高效的瞬态辐射传递模拟方法的研究有着两方面的意义： 首先，瞬态辐射模拟对丁|实验有着重要的指导意义。可以分析极短时间间隔内辐射场的 变化情况。其次，从数值模拟的角度来看，由于时间变量的引入，势必加强了辐射模拟 的复杂程度．瞬态模拟也因而成为检验模型的重要手段。因为，在稳态问题中具有相当 准确程度的数值方法可能在瞬态阀题的求解中遇到许多困难，如数值稳定性，伪振荡等。 冈此，当前的研究目的之一也是检验辐射传递的有限元法在瞬态模拟问题的准确性和稳 定性。 的结果经常产生虚假的振荡。这种现象是由于控制方程中一阶导数项即对流项的出现引 起的。对于对流占优的控制方程，有限元的数值求解会产生较大的不稳定性。这种不稳 定性导致了结果存在不真实的伪振荡。由T-GFEM缺乏消除这种不真实的振荡的有效的 迎风阻尼，因而需要其它的数值技术对其进行改良。为了消除非真实的振荡，一些有限 (StreamlinePetrov．Galerkin Upwind 方法虽然准确，但比较复杂。最小二乘有限元法的刚度矩阵是对称正定的。冈而，仅仅 需要存储一半的稀疏矩阵，一些高效的矩阵求解方法可以用来求解这种对称正定系数矩 非对称的，因而需要非常有效的迭代求解方法。而且在多维问题的模拟中需要非常