En la radiació atmosfèrica , la funció X i Y de Chandrasekhar apareix com les solucions dels problemes que comporten reflexió difusa i transmissió, introduïda per l'astrofísic indi-americà Subrahmanyan Chandrasekhar .[1] [2] [3] [4] [5] La funció X i Y de Chandrasekhar X ( μ ) , Y ( μ ) {\displaystyle X(\mu ),\ Y(\mu )} definida a l'interval 0 ≤ μ ≤ 1 {\displaystyle 0\leq \mu \leq 1} , satisfà la parella d'equacions integrals no lineals:
X ( μ ) = 1 + μ ∫ 0 1 Ψ ( μ ′ ) μ + μ ′ [ X ( μ ) X ( μ ′ ) − Y ( μ ) Y ( μ ′ ) ] d μ ′ , Y ( μ ) = e − τ 1 / μ + μ ∫ 0 1 Ψ ( μ ′ ) μ − μ ′ [ Y ( μ ) X ( μ ′ ) − X ( μ ) Y ( μ ′ ) ] d μ ′ {\displaystyle {\begin{aligned}X(\mu )&=1+\mu \int _{0}^{1}{\frac {\Psi (\mu ')}{\mu +\mu '}}[X(\mu )X(\mu ')-Y(\mu )Y(\mu ')]\,d\mu ',\\[5pt]Y(\mu )&=e^{-\tau _{1}/\mu }+\mu \int _{0}^{1}{\frac {\Psi (\mu ')}{\mu -\mu '}}[Y(\mu )X(\mu ')-X(\mu )Y(\mu ')]\,d\mu '\end{aligned}}} on la funció característica Ψ ( μ ) {\displaystyle \Psi (\mu )} és un polinomi parell μ {\displaystyle \mu } que generalment compleixen la condició
∫ 0 1 Ψ ( μ ) d μ ≤ 1 2 , {\displaystyle \int _{0}^{1}\Psi (\mu )\,d\mu \leq {\frac {1}{2}},} i 0 < τ 1 < ∞ {\displaystyle 0<\tau _{1}<\infty } és el gruix òptic de l'atmosfera. Si la igualtat es compleix en la condició anterior, s'anomena cas conservador , altrament cas no conservador . Aquestes funcions estan relacionades amb la funció H de Chandrasekhar
X ( μ ) → H ( μ ) , Y ( μ ) → 0 com τ 1 → ∞ {\displaystyle X(\mu )\rightarrow H(\mu ),\quad Y(\mu )\rightarrow 0\ {\text{com}}\ \tau _{1}\rightarrow \infty } i també
X ( μ ) → 1 , Y ( μ ) → e − τ 1 / μ com τ 1 → 0. {\displaystyle X(\mu )\rightarrow 1,\quad Y(\mu )\rightarrow e^{-\tau _{1}/\mu }\ {\text{com}}\ \tau _{1}\rightarrow 0.}
Aproximació
Propietats Si X ( μ , τ 1 ) , Y ( μ , τ 1 ) {\displaystyle X(\mu ,\tau _{1}),\ Y(\mu ,\tau _{1})} són les solucions per a un valor particular de τ 1 {\displaystyle \tau _{1}} , llavors les solucions per a altres valors de τ 1 {\displaystyle \tau _{1}} s'obtenen de les següents equacions integro-diferencials ∂ X ( μ , τ 1 ) ∂ τ 1 = Y ( μ , τ 1 ) ∫ 0 1 d μ ′ μ ′ Ψ ( μ ′ ) Y ( μ ′ , τ 1 ) , ∂ Y ( μ , τ 1 ) ∂ τ 1 + Y ( μ , τ 1 ) μ = X ( μ , τ 1 ) ∫ 0 1 d μ ′ μ ′ Ψ ( μ ′ ) Y ( μ ′ , τ 1 ) {\displaystyle {\begin{aligned}{\frac {\partial X(\mu ,\tau _{1})}{\partial \tau _{1}}}&=Y(\mu ,\tau _{1})\int _{0}^{1}{\frac {d\mu '}{\mu '}}\Psi (\mu ')Y(\mu ',\tau _{1}),\\{\frac {\partial Y(\mu ,\tau _{1})}{\partial \tau _{1}}}+{\frac {Y(\mu ,\tau _{1})}{\mu }}&=X(\mu ,\tau _{1})\int _{0}^{1}{\frac {d\mu '}{\mu '}}\Psi (\mu ')Y(\mu ',\tau _{1})\end{aligned}}} ∫ 0 1 X ( μ ) Ψ ( μ ) d μ = 1 − [ 1 − 2 ∫ 0 1 Ψ ( μ ) d μ + { ∫ 0 1 Y ( μ ) Ψ ( μ ) d μ } 2 ] 1 / 2 . {\displaystyle \int _{0}^{1}X(\mu )\Psi (\mu )\,d\mu =1-\left[1-2\int _{0}^{1}\Psi (\mu )\,d\mu +\left\{\int _{0}^{1}Y(\mu )\Psi (\mu )\,d\mu \right\}^{2}\right]^{1/2}.} En casos conservadors , aquesta integral pròpia es redueix a ∫ 0 1 [ X ( μ ) + Y ( μ ) ] Ψ ( μ ) d μ = 1. {\displaystyle \int _{0}^{1}[X(\mu )+Y(\mu )]\Psi (\mu )\,d\mu =1.} Si s'introdueixen les simplificacions x n = ∫ 0 1 X ( μ ) Ψ ( μ ) μ n d μ , y n = ∫ 0 1 Y ( μ ) Ψ ( μ ) μ n d μ , α n = ∫ 0 1 X ( μ ) μ n d μ , β n = ∫ 0 1 Y ( μ ) μ n d μ {\displaystyle x_{n}=\int _{0}^{1}X(\mu )\Psi (\mu )\mu ^{n}\,d\mu ,\ y_{n}=\int _{0}^{1}Y(\mu )\Psi (\mu )\mu ^{n}\,d\mu ,\ \alpha _{n}=\int _{0}^{1}X(\mu )\mu ^{n}\,d\mu ,\ \beta _{n}=\int _{0}^{1}Y(\mu )\mu ^{n}\,d\mu } , llavors tenim una relació que indica ( 1 − x 0 ) x 2 + y 0 y 2 + 1 2 ( x 1 2 − y 1 2 ) = ∫ 0 1 Ψ ( μ ) μ 2 d μ . {\displaystyle (1-x_{0})x_{2}+y_{0}y_{2}+{\frac {1}{2}}(x_{1}^{2}-y_{1}^{2})=\int _{0}^{1}\Psi (\mu )\mu ^{2}\,d\mu .} En el cas conservador , això es redueix a y 0 ( x 2 + y 2 ) + 1 2 ( x 1 2 − y 1 2 ) = ∫ 0 1 Ψ ( μ ) μ 2 d μ {\displaystyle y_{0}(x_{2}+y_{2})+{\frac {1}{2}}(x_{1}^{2}-y_{1}^{2})=\int _{0}^{1}\Psi (\mu )\mu ^{2}\,d\mu } Si la funció característica és Ψ ( μ ) = a + b μ 2 {\displaystyle \Psi (\mu )=a+b\mu ^{2}} , on a , b {\displaystyle a,b} són dues constants, aleshores tenim α 0 = 1 + 1 2 [ a ( α 0 2 − β 0 2 ) + b ( α 1 2 − β 1 2 ) ] {\displaystyle \alpha _{0}=1+{\frac {1}{2}}[a(\alpha _{0}^{2}-\beta _{0}^{2})+b(\alpha _{1}^{2}-\beta _{1}^{2})]} . En els casos conservadors , les solucions no són úniques. Si X ( μ ) , Y ( μ ) {\displaystyle X(\mu ),\ Y(\mu )} són solucions de l'equació original, llavors també són les solucions d'aquestes dues funcions F ( μ ) = X ( μ ) + Q μ [ X ( μ ) + Y ( μ ) ] , G ( μ ) = Y ( μ ) + Q μ [ X ( μ ) + Y ( μ ) ] {\displaystyle F(\mu )=X(\mu )+Q\mu [X(\mu )+Y(\mu )],\ G(\mu )=Y(\mu )+Q\mu [X(\mu )+Y(\mu )]} , on Q {\displaystyle Q} és una constant arbritària.
Referències