Математическая архитектура

Эта глава излагает формальную архитектуру Schubart Master Formula: что интегрируется, по какой области и почему получающееся выражение размерно согласовано. Обосновав концепцию физически в предыдущих главах, цель здесь — сделать её математический скелет полностью явным, чтобы любая последующая численная или аналитическая работа опиралась на однозначное определение.

Расширенное физическое представление

Каноническая Schubart Master Formula — это компактная объёмная форма $P(t) = \eta \cdot \int_V \Phi_{\text{eff}}(\mathbf{r},t) \cdot \sigma_{\text{eff}}(E)\,dV$. Чтобы сделать её математическую архитектуру полностью явной, эта глава работает с её расширенным физическим представлением,

$$P(t) = \iint \Phi(E,\theta,t)\,\sigma(E)\,\Delta p\,v_{\text{ph}}\,\eta_{\text{conv}}\;dE\,d\theta .$$

Здесь $P(t)$ — мгновенная полезная электрическая мощность [Вт]; $\Phi(E,\theta,t)$ — плотность потока излучения; $\sigma(E)$ — сечение рассеяния; $\Delta p$ — направленная передача импульса; $v_{\text{ph}}$ — скорость фононов; а $\eta_{\text{conv}}$ — эффективность преобразования. Двойной интеграл берётся по энергетической переменной $E$ и угловой переменной $\theta$. Это и есть расширенное физическое представление; то же выражение, вместе со стандартным глоссарием символов, воспроизведено в уравнении и справочнике по нотации.

Чтение подынтегрального выражения как произведения стадий

Подынтегральное выражение факторизуется в цепочку физических стадий, и эта факторизация — организующая идея всей концепции. Читая слева направо: $\Phi$ поставляет падающий ресурс; $\sigma$ задаёт, насколько сильно этот ресурс связывается с материалом; $\Delta p$ превращает связь в направленный импульс; $v_{\text{ph}}$ переносит этот импульс; а $\eta_{\text{conv}}$ выпрямляет его в полезную мощность. Поскольку они входят как произведение, всё подынтегральное выражение обращается в нуль, если обращается в нуль любой множитель — математическое выражение аргумента «цепочки», развитого в неравновесном преобразовании энергии. Форма произведения также означает, что множители в первом приближении независимо настраиваемы: улучшение переноса не меняет связь, и наоборот, — именно это делает концепцию полезной для изоляции узких мест.

Факторизацию можно сделать точной. Только два множителя несут переменные интегрирования — поток $\Phi(E,\theta,t)$ зависит от $E$ и $\theta$, а сечение $\sigma(E)$ зависит от $E$ — тогда как $\Delta p$, $v_{\text{ph}}$ и $\eta_{\text{conv}}$ в базовой форме трактуются как эффективные константы, представляющие преобразователь. При этом предположении постоянные множители выносятся из-под интеграла,

$$P(t) = \Delta p\, v_{\text{ph}}\, \eta_{\text{conv}} \iint \Phi(E,\theta,t)\,\sigma(E)\;dE\,d\theta ,$$

так что всё спектральное и угловое содержание сосредоточено в интеграле перекрытия $\iint \Phi\,\sigma\,dE\,d\theta$. Это разделение и есть то, что лицензирует методологию «изолируй узкое место»: можно изучать перекрытие окружения и материала независимо от предфактора переноса и выпрямления. Когда какой-либо из трёх предфакторов на деле зависит от энергии или угла — например $\Delta p(E,\theta)$ из микроскопического усреднения — он должен оставаться под интегралом, и разделение выше справедливо лишь приближённо. Концепция формулирует общую форму со всеми множителями под интегралом именно для того, чтобы это обобщение не требовало изменения структуры.

Дифференцирование факторизованной формы даёт чувствительность выхода к каждой стадии. Поскольку зависимость мультипликативна, относительные чувствительности просто складываются:

$$\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} .$$

Десятипроцентное улучшение любого отдельного множителя даёт десятипроцентное улучшение выхода, и наихудший множитель доминирует над достижимым приростом. Эта логарифмически-производная структура — количественный хребет аргумента об узком месте.

Область интегрирования

Область интегрирования заслуживает внимания. Энергетический интеграл берётся по диапазону энергий, присутствующих в окружающем спектре, с которыми преобразователь может связываться; формально его можно брать от нуля до бесконечности, с пониманием того, что подынтегральное выражение пренебрежимо мало вне носителя $\Phi(E,\theta,t)\,\sigma(E)$. Угловой интеграл берётся по телесному углу, стягиваемому падающим потоком у поверхности преобразователя. Во многих трактовках $\theta$ — полярный угол, и полное интегрирование по телесному углу вводило бы обычную меру $\sin\theta$; концепция держит $\theta$ явным именно для того, чтобы любую угловую меру, подходящую для данной геометрии, можно было подставить без изменения общей структуры. Важный формальный момент в том, что пределы задаются физическим носителем подынтегрального выражения, а не выбираются произвольно.

Зависимость от времени

Выход несёт явный аргумент времени, потому что его несёт поток. Время входит через $\Phi(E,\theta,t)$; множители связи и выпрямления трактуются как квазистатические, то есть они реагируют практически мгновенно на изменения потока на интересующих временных масштабах. Следовательно, $P(t)$ наследует свою зависимость от времени от окружения, и концепция естественно описывает преобразователь, чей выход отслеживает флуктуирующее окружающее поле. Если бы какие-либо из материальных параметров сами зависели от времени — например если бы $\eta_{\text{conv}}$ дрейфовал с температурой — им тоже можно было бы придать аргумент времени, не меняя формы интеграла.

Квазистатическое предположение можно сформулировать как разделение временных масштабов. Пусть $\tau_{\text{env}}$ — характерное время, за которое меняется окружающий поток, а $\tau_{\text{resp}}$ — внутреннее время отклика преобразователя (задаваемое рассеянием и сбором). Мгновенная форма основной формулы справедлива, когда

$$\tau_{\text{resp}} \ll \tau_{\text{env}} ,$$

так что в каждый момент преобразователь фактически достиг своего установившегося отклика на текущий поток. Когда это выполнено, можно усреднить выход по окну наблюдения $T$ и получить среднюю мощность,

$$\langle P \rangle = \frac{1}{T}\int_0^T P(t)\;dt ,$$

которая и есть величина, фактически сообщаемая большинством экспериментов. Если разделение нарушается — поток меняется быстрее, чем преобразователь может откликнуться — концепция потребовала бы свёртки с ядром отклика вместо мгновенного произведения, расширение, которое оставляет спектрально-угловую структуру нетронутой, но заменяет квазистатические множители функциями памяти.

Порядок интегрирования и разделимость

Поскольку подынтегральное выражение есть произведение спектрально-угловой части и набора квазипостоянных множителей, теорема Фубини гарантирует, что двойной интеграл можно вычислять в любом порядке, $\int dE \int d\theta$ или $\int d\theta \int dE$, всякий раз, когда подынтегральное выражение абсолютно интегрируемо — а оно таково, поскольку $\Phi$ и $\sigma$ неотрицательны, а подынтегральное выражение имеет конечный носитель. Эта свобода практически полезна: можно сначала проинтегрировать по углу, получив эффективный спектральный поток $\Phi_{\text{eff}}(E)=\int \Phi(E,\theta,t)\,d\theta$, а затем выполнить единственный интеграл спектрального перекрытия,

$$P(t) = \Delta p\, v_{\text{ph}}\, \eta_{\text{conv}} \int_0^\infty \Phi_{\text{eff}}(E,t)\,\sigma(E)\;dE ,$$

или, наоборот, сначала проинтегрировать по энергии, получив угловую диаграмму мощности. Два пути должны согласовываться, и проверка того, что это так, — простой, но эффективный тест на самосогласованность любой численной реализации. Когда множитель действительно связывает две переменные — например угловозависимое сечение $\sigma(E,\theta)$ — разделимость теряется и полный двойной интеграл должен быть сохранён, но перестановка порядка интегрирования остаётся справедливой. Это ещё одна причина, по которой концепция записана с обоими явными интегралами, а не сведена заранее к эффективной одномерной форме.

Размерная согласованность

Концепция такого рода заслуживает доверия, только если она размерно обоснована, так что стоит это проверить. Плотность потока $\Phi$ имеет единицы энергии на площадь, на время, на телесный угол, на энергию — то есть $\mathrm{W\,m^{-2}\,sr^{-1}\,J^{-1}}$, когда плотность энергии выражена на джоуль. Сечение $\sigma$ вносит $\mathrm{m^{2}}$, передача импульса $\Delta p$ вносит $\mathrm{kg\,m\,s^{-1}}$, скорость фононов $v_{\text{ph}}$ вносит $\mathrm{m\,s^{-1}}$, а $\eta_{\text{conv}}$ безразмерна. Мера интегрирования $dE\,d\theta$ вносит энергию, умноженную на телесный угол, $\mathrm{J\,sr}$. В любой конкретной реализации концепции произведение этих единиц, вместе с любой константой пропорциональности, связывающей геометрические множители с мощностью, должно сводиться к ваттам. Концепция требует, чтобы это сведение выполнялось; конкретное применение должно поставить константы, замыкающие бухгалтерию точно. Эта дисциплина — настаивать на том, чтобы единицы сходились прежде, чем какому-либо числу можно доверять, — и есть то, что удерживает формализм от дрейфа в нефальсифицируемую область.

Положительность и ограниченность

Два структурных свойства убеждают нас в том, что формула ведёт себя хорошо. Во-первых, положительность: поток $\Phi$ неотрицателен по определению, сечение $\sigma$ неотрицательно, $v_{\text{ph}}$ положительна, а $\eta_{\text{conv}}$ лежит в $[0,1]$. Единственный множитель, который может нести знак, — это направленная передача импульса $\Delta p$, которая положительна, когда чистый толчок направлен к оси сбора, и была бы отрицательной при толчке против неё. В правильно ориентированном преобразователе $\Delta p \ge 0$, так что подынтегральное выражение неотрицательно и $P(t) \ge 0$: концепция не может сообщить отрицательную выдаваемую мощность, как и требует физический смысл. Во-вторых, ограниченность: поскольку $\eta_{\text{conv}} \le 1$, выдаваемая мощность никогда не может превысить полную выпрямляемую мощность, пересекающую границу,

$$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). Физическое обоснование предфактора переноса $v_{\text{ph}}$.