Architecture mathématique

Ce chapitre expose l’architecture formelle de la Schubart Master Formula : ce qui est intégré, sur quel domaine, et pourquoi l’expression résultante est dimensionnellement cohérente. Le cadre ayant été motivé physiquement dans les chapitres précédents, l’objectif est ici de rendre pleinement explicite son squelette mathématique, afin que tout travail numérique ou analytique ultérieur repose sur une définition sans ambiguïté.

La représentation physique développée

La Schubart Master Formula canonique est la forme volumique compacte $P(t) = \eta \cdot \int_V \Phi_{\text{eff}}(\mathbf{r},t) \cdot \sigma_{\text{eff}}(E)\,dV$. Pour rendre pleinement explicite son architecture mathématique, ce chapitre travaille avec sa représentation physique développée,

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

Ici $P(t)$ est la puissance électrique utile instantanée [W] ; $\Phi(E,\theta,t)$ est la densité de flux de rayonnement ; $\sigma(E)$ est la section efficace de diffusion ; $\Delta p$ est le transfert de quantité de mouvement dirigée ; $v_{\text{ph}}$ est la vitesse des phonons ; et $\eta_{\text{conv}}$ est le rendement de conversion. La double intégrale porte sur la variable d’énergie $E$ et la variable angulaire $\theta$. C’est la représentation physique développée ; la même expression, accompagnée du glossaire standard des symboles, est reproduite sur la page de l’équation et la référence de notation.

Lire l’intégrande comme un produit d’étapes

L’intégrande se factorise en une chaîne d’étapes physiques, et cette factorisation est l’idée organisatrice de tout le cadre. De gauche à droite : $\Phi$ fournit la ressource incidente ; $\sigma$ fixe l’intensité du couplage de cette ressource au matériau ; $\Delta p$ convertit le couplage en quantité de mouvement dirigée ; $v_{\text{ph}}$ transporte cette quantité de mouvement ; et $\eta_{\text{conv}}$ la rectifie en puissance utile. Parce que ces facteurs apparaissent comme un produit, l’intégrande tout entier s’annule si l’un quelconque des facteurs s’annule — l’expression mathématique de l’argument de « chaîne » développé dans conversion d’énergie hors d’équilibre. La forme produit signifie aussi que les facteurs sont, au premier ordre, ajustables indépendamment : améliorer le transport ne change pas le couplage, et réciproquement, ce qui rend le cadre utile pour isoler les goulets d’étranglement.

La factorisation peut être rendue précise. Seuls deux facteurs portent les variables d’intégration — le flux $\Phi(E,\theta,t)$ dépend de $E$ et de $\theta$, et la section efficace $\sigma(E)$ dépend de $E$ — tandis que $\Delta p$, $v_{\text{ph}}$ et $\eta_{\text{conv}}$ sont, dans la forme de base, traités comme des constantes effectives représentatives du convertisseur. Sous cette hypothèse, les facteurs constants sortent de l’intégrale,

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

de sorte que tout le contenu spectral et angulaire se concentre dans l’intégrale de recouvrement $\iint \Phi\,\sigma\,dE\,d\theta$. Cette séparation est ce qui autorise la méthodologie « isoler le goulet d’étranglement » : on peut étudier le recouvrement environnement–matériau indépendamment du préfacteur de transport et de rectification. Lorsque l’un des trois préfacteurs dépend en réalité de l’énergie ou de l’angle — par exemple $\Delta p(E,\theta)$ issu de la moyenne microscopique — il doit rester à l’intérieur de l’intégrale, et la séparation ci-dessus ne vaut qu’approximativement. Le cadre énonce la forme générale avec tous les facteurs à l’intérieur de l’intégrale précisément pour que cette généralisation ne requière aucun changement de structure.

Différencier la forme factorisée donne la sensibilité de la sortie à chaque étape. Parce que la dépendance est multiplicative, les sensibilités relatives s’ajoutent simplement :

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

Une amélioration de dix pour cent d’un facteur quelconque donne une amélioration de dix pour cent de la sortie, et le facteur le plus mauvais domine le gain atteignable. Cette structure de dérivée logarithmique est l’ossature quantitative de l’argument du goulet d’étranglement.

Le domaine d’intégration

Le domaine d’intégration mérite de l’attention. L’intégrale en énergie porte sur la plage d’énergies présentes dans le spectre ambiant et auxquelles le convertisseur peut se coupler ; formellement on peut la prendre de zéro à l’infini, étant entendu que l’intégrande est négligeable hors du support de $\Phi(E,\theta,t)\,\sigma(E)$. L’intégrale angulaire porte sur l’angle solide sous-tendu par le flux incident à la surface du convertisseur. Dans de nombreux traitements $\theta$ est un angle polaire et une intégration sur l’angle solide complet introduirait la mesure habituelle $\sin\theta$ ; le cadre garde $\theta$ explicite précisément pour que la mesure angulaire appropriée à une géométrie donnée puisse être insérée sans changer la structure globale. Le point formel important est que les bornes sont fixées par le support physique de l’intégrande, et non choisies arbitrairement.

Dépendance temporelle

La sortie porte un argument temporel explicite parce que le flux en porte un. Le temps entre par $\Phi(E,\theta,t)$ ; les facteurs de couplage et de rectification sont traités comme quasi statiques, c’est-à-dire qu’ils répondent de façon effectivement instantanée aux variations du flux aux échelles de temps d’intérêt. Par conséquent $P(t)$ hérite sa dépendance temporelle de l’environnement, et le cadre décrit naturellement un convertisseur dont la sortie suit un champ ambiant fluctuant. Si l’un des paramètres matériels était lui-même dépendant du temps — par exemple si $\eta_{\text{conv}}$ dérivait avec la température — il pourrait lui aussi recevoir un argument temporel sans altérer la forme de l’intégrale.

L’hypothèse quasi statique peut s’énoncer comme une séparation des échelles de temps. Soit $\tau_{\text{env}}$ le temps caractéristique sur lequel le flux ambiant varie et $\tau_{\text{resp}}$ le temps de réponse interne du convertisseur (fixé par la diffusion et la collecte). La forme instantanée de la formule maîtresse est valide lorsque

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

de sorte qu’à chaque instant le convertisseur a effectivement atteint sa réponse stationnaire au flux courant. Lorsque cela est vérifié, on peut moyenner la sortie sur une fenêtre d’observation $T$ pour obtenir une puissance moyenne,

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

qui est la grandeur que la plupart des expériences rapportent effectivement. Si la séparation venait à échouer — un flux variant plus vite que le convertisseur ne peut répondre — le cadre exigerait une convolution avec un noyau de réponse à la place du produit instantané, une extension qui laisse intacte la structure spectrale–angulaire mais remplace les facteurs quasi statiques par des fonctions mémoire.

Ordre d’intégration et séparabilité

Parce que l’intégrande est un produit d’une partie spectrale–angulaire et d’un ensemble de facteurs quasi constants, le théorème de Fubini garantit que la double intégrale peut être évaluée dans l’un ou l’autre ordre, $\int dE \int d\theta$ ou $\int d\theta \int dE$, dès lors que l’intégrande est absolument intégrable — ce qu’il est, puisque $\Phi$ et $\sigma$ sont non négatifs et que l’intégrande a un support fini. Cette liberté est pratiquement utile : on peut intégrer d’abord sur l’angle pour obtenir un flux spectral effectif $\Phi_{\text{eff}}(E)=\int \Phi(E,\theta,t)\,d\theta$ puis effectuer une seule intégrale de recouvrement spectral,

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

ou, inversement, intégrer d’abord sur l’énergie pour obtenir un diagramme de puissance angulaire. Les deux voies doivent concorder, et vérifier qu’elles le font est un test d’auto-cohérence simple mais efficace pour toute implémentation numérique. Lorsqu’un facteur couple véritablement les deux variables — par exemple une section efficace dépendant de l’angle $\sigma(E,\theta)$ — la séparabilité est perdue et la double intégrale complète doit être conservée, mais l’interversion de l’ordre d’intégration reste valide. C’est une raison de plus pour laquelle le cadre est écrit avec les deux intégrales explicites plutôt que pré-réduites à une forme effective unidimensionnelle.

Cohérence dimensionnelle

Un cadre de ce genre n’est crédible que s’il est dimensionnellement sain, il vaut donc la peine de le vérifier. La densité de flux $\Phi$ a les unités d’énergie par surface, par temps, par angle solide, par énergie — c’est-à-dire $\mathrm{W\,m^{-2}\,sr^{-1}\,J^{-1}}$ lorsque la densité d’énergie est exprimée par joule. La section efficace $\sigma$ contribue $\mathrm{m^{2}}$, le transfert de quantité de mouvement $\Delta p$ contribue $\mathrm{kg\,m\,s^{-1}}$, la vitesse des phonons $v_{\text{ph}}$ contribue $\mathrm{m\,s^{-1}}$, et $\eta_{\text{conv}}$ est sans dimension. La mesure d’intégration $dE\,d\theta$ contribue énergie fois angle solide, $\mathrm{J\,sr}$. Dans toute réalisation concrète du cadre, le produit de ces unités, avec la constante de proportionnalité qui relie les facteurs géométriques à la puissance, doit se réduire à des watts. Le cadre exige que cette réduction soit vérifiée ; une application spécifique doit fournir les constantes qui font que la comptabilité se ferme exactement. Cette discipline — exiger que les unités s’équilibrent avant qu’aucun nombre ne soit cru — est ce qui empêche le formalisme de dériver vers un terrain non falsifiable.

Positivité et bornitude

Deux propriétés structurelles nous rassurent sur le bon comportement de la formule. D’abord, la positivité : le flux $\Phi$ est non négatif par définition, la section efficace $\sigma$ est non négative, $v_{\text{ph}}$ est positive, et $\eta_{\text{conv}}$ appartient à $[0,1]$. Le seul facteur qui peut porter un signe est le transfert de quantité de mouvement dirigée $\Delta p$, qui est positif lorsque la poussée nette est dirigée vers l’axe de collecte et serait négatif pour une poussée à contresens. Dans un convertisseur correctement orienté $\Delta p \ge 0$, donc l’intégrande est non négatif et $P(t) \ge 0$ : le cadre ne peut pas rapporter une puissance délivrée négative, comme le sens physique l’exige. Ensuite, la bornitude : parce que $\eta_{\text{conv}} \le 1$, la puissance délivrée ne peut jamais excéder la puissance totale rectifiable traversant la frontière,

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

qui est elle-même bornée par le flux incident. Aucune configuration des paramètres ne permet à la sortie d’excéder la ressource fournie par l’environnement. La positivité et la bornitude sont ensemble les garanties mathématiques qui correspondent, respectivement, à « aucun courant inverse spontané » et « aucune sortie au-delà de la ressource disponible » — les deux choses qu’un modèle crédible de conversion d’énergie ne doit jamais violer.

De la structure à l’évaluation

L’architecture étant fixée, le travail restant est l’interprétation et l’évaluation. La signification et les unités de chaque facteur sont cataloguées dans définitions des variables ; la surface de contrôle à travers laquelle le bilan est établi est spécifiée dans frontières du système ; la construction pas à pas complète est donnée dans dérivation complète ; et les relations directrices qui découlent de la structure sont rassemblées dans relations directrices.

Pour aller plus loin

• R. Courant et D. Hilbert, Methods of Mathematical Physics, Vol. 1 (Interscience, 1953). Référence standard pour les opérateurs intégraux et les arguments de séparabilité utilisés ici.
• G. B. Arfken, H. J. Weber et F. E. Harris, Mathematical Methods for Physicists, 7e éd. (Academic Press, 2013). Intégration multidimensionnelle, mesures et raisonnement dimensionnel.
• H. B. Callen, Thermodynamics and an Introduction to Thermostatistics, 2e éd. (Wiley, 1985). Sur la cohérence dimensionnelle et les unités des expressions de bilan énergétique.
• J. M. Ziman, Electrons and Phonons (Oxford University Press, 1960). Fondement physique du préfacteur de transport $v_{\text{ph}}$.