(19)中华 人民共和国 国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202111594201.5 (22)申请日 2021.12.24 (65)同一申请的已公布的文献号 申请公布号 CN 113987841 A (43)申请公布日 2022.01.28 (73)专利权人 苏州浪潮智能科技有限公司 地址 215168 江苏省苏州市吴中经济开发 区郭巷街道官浦路1号9幢 (72)发明人 冉鑫　 (74)专利代理 机构 北京市万慧达律师事务所 11111 代理人 陈晓磊 (51)Int.Cl. G06F 30/20(2020.01) G06F 119/08(2020.01)(56)对比文件 CN 111986733 A,2020.1 1.24 CN 107942375 A,2018.04.20 CN 110275733 A,2019.09.24 US 2007230135 A1,2007.10.04 审查员 朱琳玲 (54)发明名称 一种求解界面处声子热输运的方法和存储 介质 (57)摘要 本发明涉及一种求解界面处声子热输运的 方法和存储介质， 方法包括： 初始化设定参数及 系统参数； 基于所述系统参数， 获取声子偏差强 度； 基于线性波尔兹曼方程和声子界面条件， 迭 代更新所述声子偏差强度， 并基于任一所述声子 偏差强度获得对应的温度、 伪温度与热流； 获取 当前所述声子偏差强度对应的温度与前一次所 述声子偏 差强度对应的温度的差值， 并将所述差 值定义为相对误差； 将所述相对误差与预设的收 敛误差进行对比， 若所述相对误差不大于所述收 敛误差， 则输出所述当前所述声子偏差强度对应 的温度、 伪温度与热流。 本发明的方法可以直接、 准确、 高效率地对半导体内界面处声子的热输运 进行求解， 且最大程度地保证求解结果的准确 性。 权利要求书4页 说明书17页 附图4页 CN 113987841 B 2022.04.19 CN 113987841 B 1.一种求 解界面处声子热输运的方法， 其特 征在于： 所述方法包括以下步骤： 步骤S1、 初始化设定参数及系统参数； 步骤S2、 基于所述系统参数， 获取声子偏差强度； 步骤S3、 基于线性波尔兹曼方程和声子界面条件， 迭代更新所述声子偏差强度， 基于任 一所述声子偏差强度与温度的基本定义的比值， 获得任一所述声子偏差强度对应的温度， 基于任一所述声子偏差强度与伪温度的基本定义的比值， 获得任一所述声子偏差强度对应 的伪温度， 基于任一所述声子偏差强度进行积分求和， 获得任一所述声子偏差强度对应的 热流； 步骤S4、 获取当前所述声子偏差强度对应的温度与前一 次所述声子偏差强度对应的温 度的差值， 并将所述差值定义 为相对误差； 步骤S5、 将所述相对误差与预设的收敛误差进行对比， 若所述相对误差不大于所述收 敛误差， 则输出 所述当前 所述声子偏差强度对应的温度、 伪温度与热流。 2.根据权利要求1所述的求解界面处声子热输运的方法， 其特征在于： 所述方法还包 括： 步骤S6、 获取所述迭代更新所述声子偏差强度的计算时长， 并定义所述计算时长为第 一计算时长； 步骤S7、 若所述第一计算时长小于预设的总计算时长， 则继续执行步骤S3、 S4、 S5、 S6和 S7； 反之， 则输出当前 所述声子偏差强度对应的温度、 伪温度与热流。 3.根据权利要求1或2所述的求 解界面处声子热输运的方法， 其特 征在于： 所述设定参数包括： 空间步长、 系统尺度、 时间步长、 总计算时长和收敛误差； 所述系统参数包括： 温度和伪温度； 所述声子偏差强度包括： 基于声子强度获取得到的第 一偏差强度和基于声子伪平衡强 度获取得到的第二偏差强度。 4.根据权利要求3所述的求解界面处声子热输运的方法， 其特征在于： 所述步骤S2包括 基于声子强度， 获取 得到所述第一偏差强度； 所述获取 得到所述第一偏差强度具体包括以下步骤： 步骤S21、 设定所述温度的初始值 为T0， 基于所述温度的初始值， 获取 得到声子强度为： 式中， 为声子强度， 为声子群速度， 为约化普朗克常数， 为声子角频率， p为 声子分支， 是声子分布， D为声子态密度， 表示声子分支为 p且声子角频率为 的声子态密度； 步骤S22、 基于所述声子强度， 获取 得到声子参 考平衡强度为： 权　利　要　求　书 1/4 页 2 CN 113987841 B 2式中， 为声子参考平衡强度， 为考虑常数型参考温度 的玻色‑爱因斯坦 分布； 步骤S23、 基于所述声子强度和所述声子参 考平衡强度， 获得第一偏差强度： 式中， 为第一偏差强度， 为声子强度， 为声子参 考平衡强度。 5.权利要求4所述的求解界面处声子热输运的方法， 其特征在于： 所述步骤S2还包括基 于声子伪平衡强度， 获取 得到所述第二偏差强度； 所述获取 得到所述第二偏差强度具体包括以下步骤： 步骤S24、 设定所述伪温度的初始值与所述温度的初始值相等均为T0， 基于所述伪温度 的初始值， 获取 得到声子伪平衡强度为： 式中， 为声子伪平衡强度， 是考虑初始伪温度 的玻色‑爱因斯坦分 布； 步骤S25、 基于所述声子强度和所述声子伪平衡强度， 获得第二偏差强度： 式中， 为第二偏差强度。 6.根据权利要求5所述的求 解界面处声子热输运的方法， 其特 征在于： 所述线性玻尔兹曼 方程为： 式中， 为第一偏差强度， 为第二偏差强度， 为声子弛豫时间， 表示 声子角频率， p表示声子分支， 表示声子群速度， 表示声子分支为 p且声子角频 率为 的声子群速度； 所述界面条件为： 式中， A、 B分别表示两种材料， 为界面法向向量与群速度的夹角， 界面法向正方向由权　利　要　求　书 2/4 页 3 CN 113987841 B 3