[论文解读] Efficient simulation of multidimensional phonon transport using energy-based variance-reduced Monte Carlo formulations

Jean-Philippe M. Péraud, Nicolas G. Hadjiconstantinou|DSpace@MIT (Massachusetts Institute of Technology)|Sep 18, 2011

Thermal properties of materials参考文献 28被引用 18

一句话总结

本文提出了一种基于能量的、方差缩减的蒙特卡洛方法，用于在微米和纳米尺度器件中高效且精确地模拟多维声子输运。通过采用将平衡态与偏离部分分离的控制变量公式，并模拟能量包而非单个声子，该方法实现了精确的能量守恒，显著降低了统计不确定性，从而在计算成本恒定的情况下实现了对任意微弱温度信号的模拟——该方法在多孔硅及热导率谱分析中的应用已得到验证。

ABSTRACT

We present a new Monte Carlo method for obtaining solutions of the Boltzmann equation for describing phonon transport in micro and nanoscale devices. The proposed method can resolve arbitrarily small signals (e.g. temperature differences) at small constant cost and thus represents a considerable improvement compared to traditional Monte Carlo methods whose cost increases quadratically with decreasing signal. This is achieved via a control-variate variance reduction formulation in which the stochastic particle description only solves for the deviation from a nearby equilibrium, while the latter is described analytically. We also show that simulating an energy-based Boltzmann equation results in an algorithm that lends itself naturally to exact energy conservation thereby considerably improving the simulation fidelity. Simulations using the proposed method are used to investigate the effect of porosity on the effective thermal conductivity of silicon. We also present simulations of a recently developed thermal conductivity spectroscopy process. The latter simulations demonstrate how the computational gains introduced by the proposed method enable the simulation of otherwise intractable multiscale phenomena.

研究动机与目标

• 解决传统蒙特卡洛方法在模拟声子输运中微弱温度信号时计算成本过高的问题。
• 通过基于能量的玻尔兹曼公式，实现精确的能量守恒，从而提高模拟保真度。
• 通过采用控制变量方差缩减技术，降低低信号区域的统计不确定性。
• 使以往传统方法难以处理的复杂多尺度声子输运现象得以模拟。
• 通过多孔硅和热导率谱分析的应用验证该方法。

提出的方法

• 该方法采用控制变量公式，其中平衡态可解析求得，仅对偏离平衡的部分进行随机模拟。
• 采用基于能量包而非单个声子的玻尔兹曼方程，实现精确的能量守恒。
• 偏离部分蒙特卡洛方法仅处理非平衡分量，从而将统计噪声降低至与信号大小的平方成正比，优于传统方法。
• 该公式利用了微弱信号对应于小偏离平衡态的特性，从而实现对这些偏离的低方差估计。
• 该算法实现了基于频率和极化依赖的群速度与弛豫时间，其参数源自德拜和爱因斯坦模型。
• 该方法引入瞬态边界条件，并使用海维赛德函数追踪声子在时空中的发射历史。

实验结果

研究问题

• RQ1 如何使声子输运的蒙特卡洛模拟在任意微弱温度信号下保持高效？
• RQ2 将玻尔兹曼方程表述为能量包形式，能在能量守恒和统计精度方面带来多大提升？
• RQ3 通过控制变量实现的方差缩减，在低信号声子输运问题中能将计算成本降低到何种程度？
• RQ4 孔隙率在纳米尺度硅中如何影响其有效热导率？
• RQ5 所提出的方法能否模拟复杂、瞬态且多尺度的热现象，如热导率谱分析？

主要发现

• 所提出的方法在任意微弱温度信号下保持恒定的计算成本，与传统蒙特卡洛方法的二次方增长形成鲜明对比。
• 基于能量的公式实现了精确的能量守恒，与以往方法中人为的能量守恒相比，显著提升了模拟保真度。
• 多孔硅的模拟结果表明，有序排列的孔洞比随机分布的孔洞更显著地降低有效热导率，原因在于增强了声子的弹道散射。
• 由于计算效率高，该方法成功模拟了热导率谱分析中的瞬态响应，这一问题此前在多尺度模拟中难以处理。
• 偏离部分的公式设计使得统计不确定性随信号减小而降低，从而实现了对低振幅温度梯度的高精度结果。
• 瞬态一维弹道输运的解析解已得到验证，其与推导表达式 ΔT(x,t) = ½(1−x/Vgt)H(1−x/Vgt)ΔTl + ½(1−(L−x)/Vgt)H(1−(L−x)/Vgt)ΔTr 在小温度变化假设下完全一致。