Finite-Slab Reflectance and Transmittance for Henyey-Greenstein Scattering via a First-Passage Transfer Operator

Résumé technique : Réflectance et transmittance de plaques fines pour la diffusion de Henyey–Greenstein via un opérateur de transfert de premier passage

Énoncé du problème
L'article traite du calcul de la réflectance ($R$) et de la transmittance ($T$) pour une plaque à géométrie plane d'épaisseur optique $\tau$ subissant une diffusion de type Henyey–Greenstein (HG). Le système est défini par le paramètre d'asymétrie $g$, l'albédo de diffusion simple $a$ et l'angle d'incidence $\theta_0$. Alors que les approches standard reposent sur l'échantillonnage de Monte Carlo (MC) ou sur des mécanismes de fonctions H ou de l'ordre discret, ce travail cherche une méthode qui soit exacte dans sa formulation (au-delà du modèle de transport lui-même), exempte d'échantillonnage aléatoire, et capable de fournir des résultats pour l'ensemble de l'axe de l'albédo à partir d'une seule évaluation d'opérateur.

Méthodologie : L'opérateur de transfert de premier passage
Le cœur de la méthode est un opérateur de transfert de premier passage agissant sur l'état de la direction de la profondeur $(z, \mu)$, où $z$ est la profondeur optique et $\mu$ est la cosinus de la direction. L'approche réduit la marche aléatoire 3D complète en une chaîne de Markov fermée en $(z, \mu)$ sans rejeter les corrélations azimutales ni invoquer d'approximations de diffusion.

1. Réduction d'état : Les auteurs démontrent que la probabilité d'échappement dépend uniquement de la séquence des pas de profondeur $\{z_k\}$ et des cosinus de direction $\{\mu_k\}$. L'angle azimutal évolue mais n'influence pas les coordonnées transverses pertinentes pour l'échappement en géométrie plane, permettant ainsi de modéliser le processus comme une chaîne de Markov en $(z, \mu)$ uniquement.
2. Construction de l'opérateur :
   • Vol libre : Les photons subissent des vols libres en profondeur à $\mu$ fixe, régis par une distribution exponentielle des longueurs de chemin.
   • Redistribution angulaire : Lors de la collision, la direction est mise à jour via la fonction de phase HG. En effectuant la moyenne sur l'azimut de diffusion uniforme, les auteurs dérivent un noyau de redistribution $\Psi(\mu \to \mu')$ qui dépend uniquement des cosinus de direction entrants et sortants.
   • Frontières absorbantes : L'opérateur est confiné à la plaque $(0, \tau)$ avec des frontières absorbantes. La masse atteignant une frontière est retirée, définissant les lois de réflexion ($P_R(n)$) et de transmission ($P_T(n)$) résolues par ordre pour les photons sortant après exactement $n$ collisions.
3. Pondération de l'albédo : La réflectance et la transmittance totales pour tout albédo $a$ sont obtenues en pondérant les lois résolues par ordre :
   $$R(\tau, g, a, \theta_0) = \sum_{n \ge 1} P_R(n) a^n$$
   $$T(\tau, g, a, \theta_0) = e^{-\tau/\mu_0} + \sum_{n \ge 1} P_T(n) a^n$$
   Cela permet le calcul de $R$ et $T$ pour tout $a \in [0, 1]$ à partir d'une seule évaluation d'opérateur discrétisé.
4. Mise en œuvre numérique : La méthode utilise une grille de profondeur uniforme et une quadrature de Gauss–Legendre pour les angles. L'opérateur est itéré jusqu'à un ordre de coupure $n_{max}$ suffisant pour capturer la masse pour les plaques conservatives ($a=1$).

Résultats clés et validation

• Précision : L'opérateur reproduit les résultats complets de Monte Carlo 3D avec une erreur absolue $\le 1,6 \times 10^{-3}$ sur un large espace de paramètres ($g \in [0, 0,95]$, $\tau \in [0,5, 16]$, $a \in [0,5, 1]$). La conservation de l'énergie ($R+T=1$) est satisfaite à $< 10^{-4}$ pour $a=1$.
• Distributions émergentes : La méthode fournit les distributions angulaires émergentes $R(\mu)$ et $T(\mu)$ sans coût de calcul supplémentaire, correspondant aux données MC avec une erreur $< 1\%$ en réflexion et $\sim 2\%$ en transmission (limitée par le binning MC des pics directionnels pointus).
• Épaisseur et anisotropie : Les résultats confirment que la diffusion vers l'avant ($g > 0$) supprime la réflexion et augmente la transmission. À mesure que l'épaisseur augmente, $R$ tend vers l'unité pour les milieux conservatifs, récupérant la queue diffusive $T \sim 1/\tau$.
• Milieux absorbants : Dans la limite $\tau \to \infty$, l'opérateur récupère l'albédo exact du demi-espace pour les milieux absorbants, en accord avec la solution exacte de la fonction H à $\le 6 \times 10^{-4}$.

Contributions structurelles et importance
L'article souligne que sa principale contribution est structurelle plutôt que simplement computationnelle. Le résultat central est la factorisation exacte de la réflectance de la plaque finie :
$$R_\tau(a) = \sum_{n \ge 1} P_\infty(n) S(n, \tau) a^n$$
où :

• $P_\infty(n)$ est la loi de premier retour en demi-espace (la probabilité de s'échapper du milieu semi-infini au $n$-ième collision).
• $S(n, \tau)$ est le facteur de survie en profondeur, représentant la probabilité qu'un chemin de retour d'ordre $n$ reste plus superficiel que $\tau$.

Cette factorisation établit que le problème de la plaque fine et les statistiques de retour en demi-espace sont « deux faces d'un même opérateur ». Lorsque $\tau \to \infty$, la loi de réflexion de la plaque $P_R(n, \tau)$ converge terme à terme vers la loi de demi-espace $P_\infty(n)$, récupérant le comportement de queue de type Sparre–Andersen/noyau de Cauchy.

Partie II : Notes de travail sur le facteur de survie
L'article inclut des notes préliminaires et empiriques (Partie II) concernant la forme analytique de $S(n, \tau)$.

• Mise à l'échelle (Scaling) : Les données empiriques suggèrent que $S(n, \tau)$ suit un effondrement d'échelle basé sur le chemin moyen libre de transport $\ell^* = 1/(1-g)$ et un décalage de profondeur additif $c_g$.
• Analogie Brownienne : Le facteur de survie semble suivre la fonction de répartition maximale d'une excursion brownienne, modifiée par un décalage dépendant de $g$ et une légère correction de forme.
• Statut : Ces conclusions sont explicitement étiquetées comme « notes de travail » et « empiriques », en attente de dérivation à partir de premiers principes. Elles ne sont pas revendiquées comme des résultats validés du modèle direct de la Partie I.

Limites et portée
Les auteurs notent que la méthode n'est pas destinée à surpasser les codes d'ordre discret optimisés (comme DISORT) pour le calcul d'une radiance unique. Elle est complémentaire, fournissant des lois résolues par ordre et le lien explicite avec les statistiques de retour en demi-espace que les autres solveurs n'exposent pas directement. L'implémentation actuelle fournit des distributions angulaires moyennées en azimut ; la dépendance complète BRDF/BTDF pour une incidence oblique est notée comme une extension directe. L'article conclut que l'unification de l'optique des plaques fines et des statistiques de retour en demi-espace via l'opérateur de premier passage constitue l'avancée conceptuelle majeure.