Arquitetura Matemática

Este capítulo expõe a arquitetura formal da Schubart Master Formula: o que está sendo integrado, sobre qual domínio, e por que a expressão resultante é dimensionalmente consistente. Tendo motivado o arcabouço fisicamente nos capítulos anteriores, o objetivo aqui é tornar seu esqueleto matemático plenamente explícito, de modo que qualquer trabalho numérico ou analítico subsequente repouse sobre uma definição inequívoca.

A representação física expandida

A Schubart Master Formula canônica é a forma de volume compacta $P(t) = \eta \cdot \int_V \Phi_{\text{eff}}(\mathbf{r},t) \cdot \sigma_{\text{eff}}(E)\,dV$. Para tornar sua arquitetura matemática plenamente explícita, este capítulo trabalha com sua representação física expandida,

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

Aqui $P(t)$ é a potência elétrica utilizável instantânea [W]; $\Phi(E,\theta,t)$ é a densidade de fluxo de radiação; $\sigma(E)$ é a seção de choque de espalhamento; $\Delta p$ é a transferência de momento direcionado; $v_{\text{ph}}$ é a velocidade dos fônons; e $\eta_{\text{conv}}$ é a eficiência de conversão. A integral dupla percorre a variável de energia $E$ e a variável angular $\theta$. Esta é a representação física expandida; a mesma expressão, junto com o glossário padrão de símbolos, é reproduzida em a equação e a referência de notação.

Lendo o integrando como um produto de estágios

O integrando se fatoriza em uma cadeia de estágios físicos, e essa fatoração é a ideia organizadora de todo o arcabouço. Lendo da esquerda para a direita: $\Phi$ fornece o recurso incidente; $\sigma$ fixa quão fortemente esse recurso se acopla ao material; $\Delta p$ converte o acoplamento em momento direcionado; $v_{\text{ph}}$ transporta esse momento; e $\eta_{\text{conv}}$ o retifica em potência utilizável. Porque esses fatores aparecem como um produto, todo o integrando se anula se qualquer fator se anular — a expressão matemática do argumento da “cadeia” desenvolvido em conversão de energia em não equilíbrio. A forma de produto também significa que os fatores são, em primeira ordem, ajustáveis de forma independente: melhorar o transporte não altera o acoplamento, e vice-versa, o que é o que torna o arcabouço útil para isolar gargalos.

A fatoração pode ser tornada precisa. Apenas dois fatores carregam as variáveis de integração — o fluxo $\Phi(E,\theta,t)$ depende de $E$ e $\theta$, e a seção de choque $\sigma(E)$ depende de $E$ — enquanto $\Delta p$, $v_{\text{ph}}$ e $\eta_{\text{conv}}$ são, na forma básica, tratados como constantes efetivas representativas do conversor. Sob essa hipótese os fatores constantes saem da integral,

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

de modo que todo o conteúdo espectral e angular fica concentrado na integral de sobreposição $\iint \Phi\,\sigma\,dE\,d\theta$. Essa separação é o que autoriza a metodologia de “isolar o gargalo”: pode-se estudar a sobreposição ambiente–material independentemente do pré-fator de transporte–retificação. Quando algum dos três pré-fatores for de fato dependente da energia ou do ângulo — por exemplo $\Delta p(E,\theta)$ a partir da média microscópica — ele deve permanecer dentro da integral, e a separação acima vale apenas aproximadamente. O arcabouço enuncia a forma geral com todos os fatores dentro da integral precisamente para que essa generalização não exija mudança de estrutura.

Diferenciar a forma fatorada dá a sensibilidade da saída a cada estágio. Como a dependência é multiplicativa, as sensibilidades relativas simplesmente se somam:

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

Uma melhoria de dez por cento em qualquer fator único produz uma melhoria de dez por cento na saída, e o pior fator domina o ganho alcançável. Essa estrutura de derivada logarítmica é a espinha dorsal quantitativa do argumento do gargalo.

O domínio de integração

O domínio de integração merece cuidado. A integral de energia percorre a faixa de energias presentes no espectro ambiental e às quais o conversor pode se acoplar; formalmente pode-se tomá-la de zero a infinito, no entendimento de que o integrando é desprezível fora do suporte de $\Phi(E,\theta,t)\,\sigma(E)$. A integral angular percorre o ângulo sólido subtendido pelo fluxo incidente na superfície do conversor. Em muitos tratamentos $\theta$ é um ângulo polar e uma integração completa de ângulo sólido introduziria a medida usual $\sin\theta$; o arcabouço mantém $\theta$ explícito precisamente para que qualquer medida angular apropriada a uma dada geometria possa ser inserida sem alterar a estrutura geral. O ponto formal importante é que os limites são fixados pelo suporte físico do integrando, não escolhidos arbitrariamente.

Dependência temporal

A saída carrega um argumento temporal explícito porque o fluxo o carrega. O tempo entra através de $\Phi(E,\theta,t)$; os fatores de acoplamento e retificação são tratados como quase estáticos, significando que respondem de modo efetivamente instantâneo a mudanças no fluxo nas escalas de tempo de interesse. Em consequência, $P(t)$ herda sua dependência temporal do ambiente, e o arcabouço descreve naturalmente um conversor cuja saída acompanha um campo ambiental flutuante. Se algum dos parâmetros do material fosse ele próprio dependente do tempo — por exemplo se $\eta_{\text{conv}}$ derivasse com a temperatura — ele também poderia receber um argumento temporal sem alterar a forma da integral.

A hipótese quase estática pode ser enunciada como uma separação de escalas de tempo. Seja $\tau_{\text{env}}$ o tempo característico ao longo do qual o fluxo ambiental varia e $\tau_{\text{resp}}$ o tempo interno de resposta do conversor (fixado pelo espalhamento e pela coleta). A forma instantânea da fórmula mestra é válida quando

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

de modo que a cada instante o conversor atingiu efetivamente sua resposta estacionária ao fluxo corrente. Quando isso vale, pode-se promediar a saída sobre uma janela de observação $T$ para obter uma potência média,

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

que é a grandeza que a maioria dos experimentos de fato relata. Caso a separação falhe — um fluxo variando mais rápido do que o conversor consegue responder — o arcabouço exigiria uma convolução com um núcleo de resposta no lugar do produto instantâneo, uma extensão que deixa a estrutura espectral–angular intacta mas substitui os fatores quase estáticos por funções de memória.

Ordem de integração e separabilidade

Como o integrando é um produto de uma parte espectral–angular e um conjunto de fatores quase constantes, o teorema de Fubini garante que a integral dupla pode ser avaliada em qualquer ordem, $\int dE \int d\theta$ ou $\int d\theta \int dE$, sempre que o integrando for absolutamente integrável — o que ele é, já que $\Phi$ e $\sigma$ são não negativos e o integrando tem suporte finito. Essa liberdade é praticamente útil: pode-se integrar sobre o ângulo primeiro para obter um fluxo espectral efetivo $\Phi_{\text{eff}}(E)=\int \Phi(E,\theta,t)\,d\theta$ e então executar uma única integral de sobreposição espectral,

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

ou, inversamente, integrar sobre a energia primeiro para obter um padrão angular de potência. As duas rotas devem concordar, e verificar que concordam é um teste de autoconsistência simples, porém eficaz, para qualquer implementação numérica. Quando um fator genuinamente acopla as duas variáveis — por exemplo uma seção de choque dependente do ângulo $\sigma(E,\theta)$ — a separabilidade se perde e a integral dupla completa deve ser mantida, mas a troca da ordem de integração permanece válida. Esta é outra razão pela qual o arcabouço é escrito com ambas as integrais explícitas, em vez de pré-colapsadas em uma forma efetiva unidimensional.

Consistência dimensional

Um arcabouço desse tipo só é crível se for dimensionalmente sólido, de modo que vale a pena verificar. A densidade de fluxo $\Phi$ tem unidades de energia por área, por tempo, por ângulo sólido, por energia — isto é, $\mathrm{W\,m^{-2}\,sr^{-1}\,J^{-1}}$ quando a densidade de energia é expressa por joule. A seção de choque $\sigma$ contribui com $\mathrm{m^{2}}$, a transferência de momento $\Delta p$ contribui com $\mathrm{kg\,m\,s^{-1}}$, a velocidade dos fônons $v_{\text{ph}}$ contribui com $\mathrm{m\,s^{-1}}$, e $\eta_{\text{conv}}$ é adimensional. A medida de integração $dE\,d\theta$ contribui com energia vezes ângulo sólido, $\mathrm{J\,sr}$. Em qualquer realização concreta do arcabouço o produto dessas unidades, junto com qualquer constante de proporcionalidade que relacione os fatores geométricos à potência, deve reduzir-se a watts. O arcabouço exige que essa redução valha; uma aplicação específica deve fornecer as constantes que fecham a contabilidade exatamente. Essa disciplina — insistir em que as unidades se equilibrem antes de qualquer número ser confiado — é o que impede que o formalismo derive para território infalseável.

Positividade e limitação

Duas propriedades estruturais nos asseguram que a fórmula é bem comportada. Primeiro, a positividade: o fluxo $\Phi$ é não negativo por definição, a seção de choque $\sigma$ é não negativa, $v_{\text{ph}}$ é positiva, e $\eta_{\text{conv}}$ está em $[0,1]$. O único fator que pode carregar um sinal é a transferência de momento direcionado $\Delta p$, que é positiva quando o empurrão líquido é em direção ao eixo de coleta e seria negativa para um empurrão contra ele. Em um conversor devidamente orientado $\Delta p \ge 0$, de modo que o integrando é não negativo e $P(t) \ge 0$: o arcabouço não pode relatar potência entregue negativa, como o bom senso físico exige. Segundo, a limitação: porque $\eta_{\text{conv}} \le 1$, a potência entregue nunca pode exceder a potência total retificável que cruza a fronteira,

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

a qual é ela própria limitada pelo fluxo incidente. Não há configuração dos parâmetros que permita à saída exceder o recurso fornecido pelo ambiente. A positividade e a limitação juntas são as garantias matemáticas que correspondem, respectivamente, a “nenhuma corrente reversa espontânea” e “nenhuma saída além do recurso disponível” — as duas coisas que um modelo crível de conversão de energia nunca deve violar.

Da estrutura à avaliação

Com a arquitetura fixada, o trabalho restante é interpretação e avaliação. O significado e as unidades de cada fator são catalogados em definições das variáveis; a superfície de controle através da qual se traça o balanço é especificada em fronteiras do sistema; a construção completa passo a passo é dada em derivação completa; e as relações reguladoras que decorrem da estrutura são reunidas em relações reguladoras.

Leitura adicional

• R. Courant e D. Hilbert, Methods of Mathematical Physics, Vol. 1 (Interscience, 1953). Referência padrão para os operadores integrais e os argumentos de separabilidade usados aqui.
• G. B. Arfken, H. J. Weber e F. E. Harris, Mathematical Methods for Physicists, 7ª ed. (Academic Press, 2013). Integração multidimensional, medidas e raciocínio dimensional.
• H. B. Callen, Thermodynamics and an Introduction to Thermostatistics, 2ª ed. (Wiley, 1985). Sobre consistência dimensional e as unidades de expressões de balanço de energia.
• J. M. Ziman, Electrons and Phonons (Oxford University Press, 1960). Fundamentação física para o pré-fator de transporte $v_{\text{ph}}$.