architecture

数学架构

P(t) 主公式的形式结构:被积函数、积分域、量纲一致性,以及向各物理阶段的因子分解。

最后更新:

本章阐述 Schubart Master Formula 的形式架构:被积的是什么、在什么域上积分,以及为何所得表达式在量纲上一致。在前几章从物理上给出框架的动机之后,这里的目标是使其数学骨架完全显式,以便后续任何数值或解析工作都建立在一个无歧义的定义之上。

扩展的物理表示

标准的 Schubart Master Formula 是紧凑的体积形式 P(t)=η⋅∫VΦeff(r,t)⋅σeff(E) dVP(t) = \eta \cdot \int_V \Phi_{\text{eff}}(\mathbf{r},t) \cdot \sigma_{\text{eff}}(E)\,dV。为使其数学架构完全显式,本章采用它的扩展的物理表示,

P(t)=∬Φ(E,θ,t) σ(E) Δp vph ηconv  dE dθ.P(t) = \iint \Phi(E,\theta,t)\,\sigma(E)\,\Delta p\,v_{\text{ph}}\,\eta_{\text{conv}}\;dE\,d\theta .

这里 P(t)P(t) 是瞬时可用电功率 [W];Φ(E,θ,t)\Phi(E,\theta,t) 是辐射通量密度;σ(E)\sigma(E) 是散射截面;Δp\Delta p 是有向动量转移;vphv_{\text{ph}} 是声子速度;ηconv\eta_{\text{conv}} 是转换效率。双重积分遍及能量变量 EE 和角度变量 θ\theta。这就是扩展的物理表示;同一表达式连同标准符号术语表一并复现于方程页和记号参考。

将被积函数读作各阶段之积

被积函数分解为一条物理阶段的链,而这种因子分解是整个框架的组织思想。从左到右读:Φ\Phi 提供入射资源;σ\sigma 设定该资源与材料耦合的强弱;Δp\Delta p 把耦合转化为有向动量;vphv_{\text{ph}} 输运该动量;ηconv\eta_{\text{conv}} 将其整流为可用功率。由于它们以乘积形式出现,只要任一因子消失,整个被积函数就消失——这是在非平衡能量转换中展开的”链条”论证的数学表达。乘积形式还意味着各因子在一阶近似下可独立调节:改善输运不改变耦合,反之亦然,这正是使该框架在孤立瓶颈方面有用的原因。

这种因子分解可以说得更精确。只有两个因子携带积分变量——通量 Φ(E,θ,t)\Phi(E,\theta,t) 依赖于 EE 和 θ\theta,截面 σ(E)\sigma(E) 依赖于 EE——而 Δp\Delta p、vphv_{\text{ph}} 和 ηconv\eta_{\text{conv}} 在基本形式中被视为代表转换器的有效常数。在该假设下,常数因子提出积分之外,

P(t)=Δp vph ηconv∬Φ(E,θ,t) σ(E)  dE dθ,P(t) = \Delta p\, v_{\text{ph}}\, \eta_{\text{conv}} \iint \Phi(E,\theta,t)\,\sigma(E)\;dE\,d\theta ,

于是全部谱与角内容集中于重叠积分 ∬Φ σ dE dθ\iint \Phi\,\sigma\,dE\,d\theta。正是这种分离许可了”孤立瓶颈”的方法学:人们可以独立于输运–整流前因子去研究环境–材料的重叠。当三个前因子中有任何一个实际上与能量或角度相关时——例如来自微观平均的 Δp(E,θ)\Delta p(E,\theta)——它就必须留在积分之内,上述分离只近似成立。框架以全部因子都在积分之内的一般形式陈述,正是为了使这种推广无需改变结构。

对因子分解后的形式求导,给出输出对每一阶段的敏感性。由于依赖是乘性的,相对敏感性简单相加:

dPP=d(Δp)Δp+dvphvph+dηconvηconv+d ⁣(∬Φσ)∬Φσ.\frac{dP}{P} = \frac{d(\Delta p)}{\Delta p} + \frac{dv_{\text{ph}}}{v_{\text{ph}}} + \frac{d\eta_{\text{conv}}}{\eta_{\text{conv}}} + \frac{d\!\left(\iint \Phi\sigma\right)}{\iint \Phi\sigma} .

任一单个因子有百分之十的改善,就使输出有百分之十的改善,而最差的因子主导着可达的增益。这种对数导数结构是瓶颈论证的定量骨干。

积分域

积分域值得仔细对待。能量积分遍及环境谱中存在的、且转换器能够耦合的能量范围;形式上可取从零到无穷,并理解为被积函数在 Φ(E,θ,t) σ(E)\Phi(E,\theta,t)\,\sigma(E) 的支集之外可忽略。角度积分遍及入射通量在转换器表面所张的立体角。在许多处理中 θ\theta 是极角,完整的立体角积分会引入常见的 sin⁡θ\sin\theta 测度;框架保持 θ\theta 显式,正是为了让适合给定几何的任何角度测度都可插入而不改变整体结构。重要的形式要点是:积分限由被积函数的物理支集设定,而非任意选取。

时间依赖

输出带有显式的时间变量,因为通量带有它。时间通过 Φ(E,θ,t)\Phi(E,\theta,t) 进入;耦合与整流因子被视为准静态,意即它们对通量变化的响应在所关心的时间尺度上实际是瞬时的。因此 P(t)P(t) 从环境继承其时间依赖,框架自然地描述一台其输出跟踪起伏环境场的转换器。若材料参数中有任何一个本身随时间变化——例如 ηconv\eta_{\text{conv}} 随温度漂移——它们也可被赋予时间变量而不改变积分的形式。

准静态假设可表述为时间尺度的分离。设 τenv\tau_{\text{env}} 为环境通量变化的特征时间,τresp\tau_{\text{resp}} 为转换器的内部响应时间(由散射和收集设定)。主公式的瞬时形式在

τresp≪τenv,\tau_{\text{resp}} \ll \tau_{\text{env}} ,

时有效,使得在每一瞬间转换器都实际达到了对当前通量的稳态响应。当此成立时,可在观测窗口 TT 上对输出求平均以得到平均功率,

⟨P⟩=1T∫0TP(t)  dt,\langle P \rangle = \frac{1}{T}\int_0^T P(t)\;dt ,

这是多数实验实际报告的量。若该分离失效——通量变化快于转换器所能响应——则框架需要以与响应核的卷积取代瞬时乘积,这一扩展不触及谱–角结构,但以记忆函数取代准静态因子。

积分次序与可分离性

由于被积函数是一个谱–角部分与一组准常数因子之积,富比尼定理保证双重积分可按任一次序求值,∫dE∫dθ\int dE \int d\theta 或 ∫dθ∫dE\int d\theta \int dE,只要被积函数绝对可积——它确实如此,因为 Φ\Phi 和 σ\sigma 非负且被积函数具有有限支集。这种自由在实践中有用:可以先对角度积分得到有效谱通量 Φeff(E)=∫Φ(E,θ,t) dθ\Phi_{\text{eff}}(E)=\int \Phi(E,\theta,t)\,d\theta,再做单个谱重叠积分,

P(t)=Δp vph ηconv∫0∞Φeff(E,t) σ(E)  dE,P(t) = \Delta p\, v_{\text{ph}}\, \eta_{\text{conv}} \int_0^\infty \Phi_{\text{eff}}(E,t)\,\sigma(E)\;dE ,

或者反过来先对能量积分以得到一个角功率图样。两条路线必须一致,核验它们是否一致,是对任何数值实现的一项简单却有效的自洽检验。当某个因子确实耦合了两个变量时——例如一个与角度相关的截面 σ(E,θ)\sigma(E,\theta)——可分离性丧失,必须保留完整的双重积分,但积分次序的互换仍然有效。这是框架以两个积分都显式书写、而非预先坍缩为一个有效一维形式的又一原因。

量纲一致性

此类框架唯有量纲健全才可信,故值得核验。通量密度 Φ\Phi 的单位为能量每面积、每时间、每立体角、每能量——即当能量密度按每焦耳表示时为 W m−2 sr−1 J−1\mathrm{W\,m^{-2}\,sr^{-1}\,J^{-1}}。截面 σ\sigma 贡献 m2\mathrm{m^{2}},动量转移 Δp\Delta p 贡献 kg m s−1\mathrm{kg\,m\,s^{-1}},声子速度 vphv_{\text{ph}} 贡献 m s−1\mathrm{m\,s^{-1}},而 ηconv\eta_{\text{conv}} 无量纲。积分测度 dE dθdE\,d\theta 贡献能量乘立体角,J sr\mathrm{J\,sr}。在框架的任何具体实现中,这些单位之积,连同将几何因子与功率相联系的任何比例常数,必须化简为瓦特。框架要求这一化简成立;具体应用必须提供使这笔账恰好闭合的常数。这种纪律——坚持在信任任何数字之前先让单位配平——使形式体系不至于漂入不可证伪的领域。

正性与有界性

两条结构性质使我们确信公式行为良好。第一,正性:通量 Φ\Phi 按定义非负,截面 σ\sigma 非负,vphv_{\text{ph}} 为正,ηconv\eta_{\text{conv}} 落于 [0,1][0,1]。唯一可能带符号的因子是有向动量转移 Δp\Delta p,当净推力朝向收集轴时为正,逆着它时则为负。在恰当取向的转换器中 Δp≥0\Delta p \ge 0,故被积函数非负且 P(t)≥0P(t) \ge 0:框架不会报告负的输送功率,正如物理常识所要求。第二,有界性:因为 ηconv≤1\eta_{\text{conv}} \le 1,输送功率绝不能超过跨越边界的全部可整流功率,

P(t)  ≤  ∬Φ(E,θ,t) σ(E) Δp vph  dE dθ,P(t) \;\le\; \iint \Phi(E,\theta,t)\,\sigma(E)\,\Delta p\,v_{\text{ph}}\;dE\,d\theta ,

而后者本身又受入射通量所界。没有任何参数配置能让输出超过环境所供给的资源。正性与有界性合在一起,正是分别对应于”无自发反向电流”和”无超出可用资源的输出”的数学保证——这是一个可信的能量转换模型绝不能违背的两件事。

从结构到求值

架构既定,余下的工作便是诠释与求值。每个因子的含义与单位在变量定义中编目;勾画平衡的控制面在系统边界中给定;逐步的完整构造见完整推导;而由该结构导出的支配关系汇集于支配关系。

延伸阅读

  • R. Courant and D. Hilbert, Methods of Mathematical Physics, Vol. 1 (Interscience, 1953). 此处所用积分算子与可分离性论证的标准参考。
  • G. B. Arfken, H. J. Weber, and F. E. Harris, Mathematical Methods for Physicists, 7th ed. (Academic Press, 2013). 多维积分、测度与量纲推理。
  • H. B. Callen, Thermodynamics and an Introduction to Thermostatistics, 2nd ed. (Wiley, 1985). 论量纲一致性与能量平衡表达式的单位。
  • J. M. Ziman, Electrons and Phonons (Oxford University Press, 1960). 输运前因子 vphv_{\text{ph}} 的物理依据。
↑