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

Technische Samenvatting: Reflectie en Transmissie van een Eindige Slab via een First-Passage Transfer Operator voor Henyey–Greenstein Verstrooiing

Probleemstelling
Het artikel behandelt de berekening van reflectie ($R$) en transmissie ($T$) voor een plane-parallel slab met een optische dikte $\tau$ die onderhevig is aan Henyey–Greenstein (HG) verstrooiing. Het systeem wordt gedefinieerd door de asymmetrieparameter $g$, de single-scattering albedo $a$ en de invalshoek $\theta_0$. Hoewel standaardbenaderingen vertrouwen op Monte Carlo (MC) bemonstering of discrete-ordinate/H-functie mechanica, streeft dit werk naar een methode die exact is in formulering (voorbij het transportmodel zelf), vrij is van willekeurige bemonstering en in staat is om resultaten te leveren voor de gehele albedo-as vanuit een enkele operator-evaluatie.

Methodologie: De First-Passage Transfer Operator
De kern van de methode is een first-passage transfer operator die werkt op de diepte-richting toestand $(z, \mu)$, waarbij $z$ de optische diepte is en $\mu$ de richtingscosinus. De benadering reduceert de volledige 3D random walk tot een gesloten Markov-keten in $(z, \mu)$ zonder azimuthale correlaties te verwerpen of diffusie-approximaties aan te roepen.

1. Toestandsreductie: De auteurs tonen aan dat de ontsleppingskans alleen afhangt van de sequentie van dieptestappen $\{z_k\}$ en richtingscosinussen $\{\mu_k\}$. De azimuthale hoek evolueert, maar beïnvloedt de transversale coördinaten die relevant zijn voor plane-parallel ontsnapping niet, waardoor het proces als een Markov-keten in $(z, \mu)$ alleen gemodelleerd kan worden.
2. Constructie van de Operator:
   • Vrije Vlucht: Fotonen ondergaan vrije vluchten in de diepte bij een vaste $\mu$, gestuurd door een exponentiële distributie van padlengtes.
   • Angulaire Herverdeling: Bij een botsing wordt de richting bijgewerkt via de HG-fasefunctie. Door te middelen over de uniforme verstrooiingsazimuth, leiden de auteurs een herverdelingskern $\Psi(\mu \to \mu')$ af die enkel afhangt van de inkomende en uitgaande richtingscosinussen.
   • Absorberende Grenzen: De operator is begrensd binnen de slab $(0, \tau)$ met absorberende grenzen. Massa die een grens bereikt, wordt verwijderd, wat de orde-opgeloste reflectie ($P_R(n)$) en transmissiewetten ($P_T(n)$) definieert voor fotonen die ontsnappen na exact $n$ botsingen.
3. Albedo Weging: De totale reflectie en transmissie voor elke albedo $a$ worden verkregen door de orde-opgeloste wetten te wegen:
   $$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$$
   Dit maakt de berekening van $R$ en $T$ voor alle $a \in [0, 1]$ mogelijk vanuit een enkele gediscretiseerde operator-evaluatie.
4. Numerieke Implementatie: De methode maakt gebruik van een uniforme diepte-grid en Gauss–Legendre kwadratuur voor hoeken. De operator wordt geïtegreerd tot een afkaporde $n_{max}$ die voldoende is om de massa te vangen voor conservatieve slabs ($a=1$).

Belangrijkste Resultaten en Validatie

• Nauwkeurigheid: De operator reproduceert volledige 3D Monte Carlo-resultaten tot binnen een absolute fout van $\le 1,6 \times 10^{-3}$ over een breed parameterrumte ($g \in [0, 0,95]$, $\tau \in [0,5, 16]$, $a \in [0,5, 1]$). Energiebehoud ($R+T=1$) wordt voldaan tot $< 10^{-4}$ voor $a=1$.
• Emergente Distributies: De methode levert de emergente angulaire distributies $R(\mu)$ en $T(\mu)$ zonder extra computationele kosten, waarbij de MC-data wordt benaderd met $<1\%$ fout in reflectie en $\sim 2\%$ in transmissie (beperkt door de MC-binning van scherpe forward peaks).
• Dikte en Anisotropie: Resultaten bevestigen dat forward scattering ($g > 0$) reflectie onderdrukt en transmissie versterkt. Naarmate de dikte toeneemt, nadert $R$ de eenheid voor conservatieve media, waarbij de diffusieve $T \sim 1/\tau$ staart wordt hersteld.
• Absorberende Media: In het limiet $\tau \to \infty$ herstelt de operator de exacte semi-infinite albedo voor absorberende media, overeenkomend met de exacte H-functie oplossing tot $\le 6 \times 10^{-4}$.

Structurele Bijdragen en Betekenis
Het artikel benadrukt dat de primaire bijdrage structureel is en niet louter computationeel. Het centrale resultaat is de exacte factorisatie van eindige slab reflectie:
$$R_\tau(a) = \sum_{n \ge 1} P_\infty(n) S(n, \tau) a^n$$
waarbij:

• $P_\infty(n)$ de half-space first-return law is (de waarschijnlijkheid om de semi-infinite media te verlaten bij de $n$-de botsing).
• $S(n, \tau)$ de diepte-overlevingsfactor is, die de waarschijnlijkheid vertegenwoordigt dat een retourpad van orde $n$ minder diep blijft dan $\tau$.

Deze factorisatie vestigt dat het eindige slab probleem en de half-space return statistieken "twee zijden van één operator" zijn. Wanneer $\tau \to \infty$, convergeert de slab reflectiewet $P_R(n, \tau)$ term-voor-term naar de half-space wet $P_\infty(n)$, waarbij de Sparre–Andersen/Cauchy-kernel staart wordt hersteld.

Deel II: Werknotities over de Overlevingsfactor
Het artikel bevat ook voorlopige, empirische notities (Deel II) met betrekking tot de analytische vorm van $S(n, \tau)$.

• Scaling: Empirische data suggereren dat $S(n, \tau)$ een scaling collapse volgt gebaseerd op de transport gemiddelde vrije weglengte $\ell^* = 1/(1-g)$ en een additieve diepte-offset $c_g$.
• Brownse Analogie: De overlevingsfactor lijkt de cumulatieve distributiefunctie van een Brownse excursie maximum te volgen, gemodificeerd door een $g$-afhankelijke offset en een milde vormcorrectie.
• Status: Deze bevindingen worden expliciet gelabeld als "werknotities" en "empirisch", wachtend op afleiding vanuit eerste principes. Ze worden niet geclaimd als gevalideerde resultaten van het model in Deel I.

Beperkingen en Reikwijdte
De auteurs merken op dat de methode niet bedoeld is om geoptimaliseerde discrete-ordinate codes (zoals DISORT) te overtreffen voor het berekenen van een enkele radiantie. In plaats daarvan is het complementair, waarbij het orde-opgeloste wetten en de expliciete verbinding met half-space statistieken biedt die andere solvers niet direct blootleggen. De huidige implementatie biedt azimuthaal gemiddelde angulaire distributies; volledige BRDF/BTDF afhankelijkheid voor schuine inval wordt genoteerd als een eenvoudige uitbreiding. Het artikel concludeert dat de unificatie van slab optica en half-space return statistieken via de first-passage operator de belangrijkste conceptuele vooruitgang is.