[论文解读] Radiative Transfer For Variable 3D Atmospheres

F. Golse, F. Hecht+3 more|arXiv (Cornell University)|Aug 11, 2022

Thermal Radiation and Cooling Technologies被引用 1

一句话总结

本文提出了一种快速的、O(N log N)的数值方法，用于在可变大气中进行三维辐射传输计算，采用H-矩阵压缩和有限元空间离散化。该方法可在一台M1笔记本电脑上，仅用不到5分钟时间，对包含50,000个网格点和680个频率带的复杂地形中，实现对温室气体吸收（CO2、CH4）引起的温度差异的高精度模拟，其核心为积分方程公式与纳维-斯托克斯方程及能量平衡方程的耦合。

ABSTRACT

To study the temperature in a gas subjected to electromagnetic radiations, one may use the Radiative Transfer equations coupled with the Navier-Stokes equations. The problem has 7 dimensions; however with minimal simplifications it is equivalent to a small number of integro-differential equations in 3 dimensions. We present the method and a numerical implementation using an H-matrix compression scheme. The result is a very fast: 50K physical points, all directions of radiation and 680 frequencies require less than 5 minutes on an Apple M1 Laptop. The method is capable of handling variable absorptioN and scattering functionS of spatial positions and frequencies. The implementation is done using htool, a matrix compression library interfaced with the PDE solver freefem++. Applications to the temperature in the French Chamonix valley is presented at different hours of the day with and without snow / clouds and with a variable absorption taken from the Gemini measurements. The result is precise enough to assert temperature differences due to increased absorption in the vibrational frequency subrange of greenhouse gasses.

研究动机与目标

• 开发一种数值高效的求解方法，用于在空间和光谱上具有变化吸收与散射特性的三维大气中进行三维辐射传输计算。
• 克服传统方法（如SHDOM和蒙特卡洛方法）在高频率和方向分辨率下的高计算成本问题。
• 实现对由于吸收系数微小变化（如CO2和CH4在振动频带中的变化）引起的温度变化的精确模拟。
• 在真实地形（法国上阿尔卑斯山查马尼克斯地区）上验证该数值工具，使用高分辨率网格和真实的气象数据。
• 证明该方法在标准硬件上具备可扩展性和实用性，适用于与气候相关的模拟。

提出的方法

• 通过将辐射传输方程在立体角上取平均，推导出积分方程公式，将原本7维的辐射传输问题简化为3维系统。
• 在积分方程上采用定点迭代格式，数学上已证明其单调且收敛，确保了数值的稳定性和准确性。
• 采用H-矩阵压缩技术，将密集矩阵运算的计算成本从O(N²)降低至O(N log N)，从而实现高效的矩阵-向量乘法。
• 在FreeFem++中实现有限元方法，使用非结构化三角形网格对复杂三维地形进行建模，并基于真实地形数据生成高精度网格。
• 以静止、雷诺平均形式求解纳维-斯托克斯方程，湍流通过有效粘度模型进行模拟，并与包含辐射加热项的能量方程相耦合。
• 通过压缩的H-矩阵计算辐射源项，对混合边界-内部矩阵条目进行特殊处理，从而高效求解积分方程。

实验结果

研究问题

• RQ1 能否开发一种快速、可扩展的数值方法，用于在复杂三维大气中实现具有可变吸收与散射特性的三维辐射传输计算？
• RQ2 H-矩阵压缩在多大程度上能降低辐射传输模拟的计算成本，同时不损失精度？
• RQ3 该方法能否分辨由CO2和CH4在振动频带中的吸收引起的微小温度差异？
• RQ4 在具有可变吸收特性的实际三维大气条件下，该方法在积分方程上的定点迭代格式是否稳定且收敛？
• RQ5 该方法能否在使用非结构化网格的复杂地形（如法国阿尔卑斯山）上，应用于真实地形，包括显著的垂直梯度？

主要发现

• 该方法在一台苹果M1笔记本电脑上，仅用不到5分钟时间，即可完成包含50,000个物理点、所有辐射方向和680个频率带的完整三维辐射传输模拟。
• H-矩阵压缩将计算成本从O(N²)降低至O(N log N)，使大规模问题的高效处理成为可能。
• 数值解的精度足以分辨由CO2和CH4在振动频带中吸收增强所引起的温度差异，该结果基于Gemini卫星测量数据建模。
• 该方法成功模拟了查马尼克斯山谷的真实温度分布，清晰展示了有无积雪、云层及吸收系数变化时的显著热响应差异。
• 积分方程上的迭代格式被证明具有单调性和收敛性，且可通过上下解获得数值误差界。
• 该代码已公开发布于GitHub，可扩展至更大计算域和超级计算机，可在数分钟内处理数百万个网格节点。