Radiative transfer 2925371 224381823 2008-07-08T16:35:51Z ResidueOfDesign 5557416 Plural of "Medium" The term '''radiative transfer''' refers to the physical phenomena of energy transfer in the form of electromagnetic radiation. The propagation of radiation through a medium is affected by [[Absorption (electromagnetic radiation)|absorption]], [[Emission (electromagnetic radiation)|emission]] and [[scattering]] processes. The '''equation of radiative transfer''' describes these interactions mathematically. Equations of radiative transfer have application in wide variety of subjects including optics, astrophysics, atmospheric science, and remote sensing. Analytic solutions to the radiative transfer equation (RTE) exist for simple cases but for more realistic media with complex multiple scattering effects numerical methods are required <ref name="chandrasekhar"> {{cite book | author=S. Chandrasekhar | title=Radiative Transfer | publisher=Dover Publications Inc. | year=1960 | isbn=0-486-60590-6 | pages=393}} </ref> <ref name="lenoble"> {{cite book | author=Jacqueline Lenoble | title=Radiative Transfer in Scattering and Absorbing Atmospheres: Standard Computational Procedures | publisher=A. Deepak Publishing | year=1985 | isbn= 0-12-451451-0 | pages=583}} </ref>. == Definitions == The fundamental quantity which describes a field of radiation is the '''spectral intensity'''. If we think of a very small area element in the radiation field, there will be radiation energy flowing through that area element. The flow can be completely characterized by the amount of energy flowing per unit time per unit [[solid angle]], the direction of the flow, and the wavelength interval being considered ([[polarization]] will be ignored for the moment). In terms of the spectral intensity, <math>I_\nu</math>, the energy flowing across an area element of area <math>da\,</math> located at <math>\mathbf{r}</math> in time <math>dt\,</math> in the solid angle <math>d\Omega</math> about the direction <math>\hat{\mathbf{n}}</math> in the frequency interval <math>\nu\,</math> to <math>\nu+d\nu\,</math> is :<math>dE_\nu = I_\nu(\mathbf{r},\hat{\mathbf{n}},t) \cos\theta ~ d\nu da d\Omega dt</math> where <math>\theta</math> is the angle that the unit direction vector <math>\hat{\mathbf{n}}</math> makes with a normal to the area element. The units of the spectral intensity are seen to be energy/time/area/solid angle/frequency. In MKS units this would be W&middot;m<sup>-2</sup>·sr<sup>-1</sup>·Hz<sup>-1</sup> (watts per square-metre-steradian-hertz). == The equation of radiative transfer== The equation of radiative transfer simply says that as a beam of radiation travels, it loses energy to the atmosphere by absorption and gains energy by atmospheric emission, and redistributes energy by scattering. The differential form of the equation for radiative transfer is: :<math>\frac{dI_\nu}{ds}=j_\nu - \alpha_\nu I_\nu + \iint \sigma_{\nu}(\Omega,{\nu'}) I_{\nu'} d\nu' d\Omega </math> where <math>j_\nu</math> is the spectral emission coefficient which accounts for emission and scattering, <math>\alpha_\nu</math> is the spectral absorption coefficient, and <math>\sigma_\nu</math> is the scattering coefficient. The emission coefficient has units of energy/time/volume/solid angle/frequency. The amount of energy contributed by a volume element <math>dV</math> is :<math>dE_\nu = j_\nu(\mathbf{r},\hat{\mathbf{n}},t) ~ d\nu dV d\Omega dt</math> The emission coefficient may be a function of the intensity if scattering occurs. The absorption coefficient <math>\alpha</math> is the infinitesimal fraction by which the intensity is reduced as the beam travels an infinitesimal distance. It has units of 1/length and the fraction by which a beam's intensity is reduced when traveling a distance <math>ds</math> is just <math>\alpha(\mathbf{r},\hat{\mathbf{n}},t)\,ds</math>. == Solutions to the equation of radiative transfer == Solutions to the equation of radiative transfer form an enormous body of work. The differences however, are essentially due to the various forms for the emission and absorption coefficients. If scattering is ignored, then a general solution in terms of the emission and absorption coefficients may be written: :<math>I_\nu(s)=I_\nu(s_0)e^{-\tau(s_0,s)}+\int_{s_0}^s j_\nu(s') e^{-\tau(s',s)}\,ds'</math> where <math>\tau(s_1,s_2)</math> is the [[optical depth]] of the atmosphere between <math>s_1</math> and <math>s_2</math>: :<math>\tau(s_1,s_2) \ \stackrel{\mathrm{def}}{=}\ \int_{s_1}^{s_2} \alpha_\nu(s)\,ds</math> === Local thermodynamic equilibrium === A particularly useful simplification of the equation of radiative transfer occurs under the conditions of [[local thermodynamic equilibrium]] (LTE). In this situation, the atmosphere consists of massive particles which are in equilibrium with each other, and therefore have a definable temperature. The radiation field is not, however in equilibrium and is being entirely driven by the presence of the massive particles. For an atmosphere in LTE, the emission coefficient and absorption coefficient are functions of temperature and density only, and are related by: :<math>\frac{j_\nu}{\alpha_\nu}=B_\nu(T)</math> where <math>B_\nu(T)</math> is the [[black body]] intensity at temperature ''T''. The solution to the equation of radiative transfer is then: :<math>I_\nu(s)=I_\nu(s_0)e^{-\tau(s_0,s)}+\int_{s_0}^s B_\nu(T(s'))\alpha_\nu(s') e^{-\tau(s',s)}\,ds'</math> Knowing the temperature profile and the density profile of the atmospheric components will be enough to calculate a solution to the equation of radiative transfer. === The Eddington approximation === The '''Eddington approximation''' can be used to obtain the intensity in a plane-parallel atmosphere with isotropic frequency-independent scattering. It assumes that the intensity is a linear function of <math>\mu=\cos\theta</math>. i.e. :<math>I_\nu(\mu,z)=a(z)+\mu b(z)</math> where <math>z</math> is the normal direction to the slab atmosphere. Note that expressing angular integrals in terms of <math>\mu</math> simplifies things because <math>d\mu=-\sin\theta d\theta</math> appears in the [[Jacobian]] of integrals in [[Spherical coordinate system|spherical coordinates]]. Extracting the first few moments of the intensity with respect to <math>\mu</math> yields :<math>J_\nu=\frac{1}{2}\int^1_{-1}I_\nu d\mu = a</math> :<math>H_\nu=\frac{1}{2}\int^1_{-1}\mu I_\nu d\mu = \frac{b}{3}</math> :<math>K_\nu=\frac{1}{2}\int^1_{-1}\mu^2 I_\nu d\mu = \frac{a}{3}</math> Thus the Eddington approximation is equivalent to setting <math>K_\nu=1/3J_\nu</math>. Higher order versions of the Eddington approximation also exist, and consist of more complicated linear relations of the intensity moments. This extra equation can be used as a closure relation for the truncated system of moments. Note that the first two moments have simple physical meanings. <math>J_\nu</math> is the isotropic intensity at a point, and <math>H_\nu</math> is the flux through that point in the <math>z</math> direction. The radiative transfer through an isotropically scattering atmosphere at local thermal equilibrium is given by :<math>\mu \frac{dI_\nu}{dz}=- \alpha_\nu (I_\nu-B_\nu) + \sigma_{\nu}(J_\nu -I_\nu)</math> Integrating over all angles yields :<math>\frac{dH_\nu}{dz}=\alpha_\nu (B_\nu-J_\nu)</math> Premultiplying by <math>\mu</math>, and then integrating over all angles gives :<math>\frac{dK_\nu}{dz}=-(\alpha_\nu+\sigma_\nu)H_\nu</math> Substituting in the closure relation, and differentiating with respect to <math>z</math> allows the two above equations to be combined to form the radiative diffusion equation :<math>\frac{d^2J_\nu}{dz^2}=3\alpha_\nu(\alpha_\nu+\sigma_\nu)(J_\nu-B_\nu)</math> This equation shows how the effective optical depth in scattering-dominated systems may be significantly different from that given by the scattering opacity if the absorptive opacity is small. == See also == * [[Absorption]] * [[Atomic line spectra]] * [[Beer-Lambert law]] * [[Emission]] * [[List of atmospheric radiative transfer codes]] * [[Scattering]] * [[Radiative transfer equation and diffusion theory for photon transport in biological tissue]] == Further reading == * {{cite book | author=S. Chandrasekhar | title=Radiative Transfer | publisher=Dover Publications Inc. | year=1960 | isbn=0-486-60590-6 | pages=393}} * {{cite book | author=Jacqueline Lenoble | title=Radiative Transfer in Scattering and Absorbing Atmospheres: Standard Computational Procedures | publisher=A. Deepak Publishing | year=1985 | isbn= 0-12-451451-0 | pages=583}} * {{cite book | author=Dimitri Mihalas, Barbara Weibel-Mihalas | title=Foundations of Radiation Hydrodynamics | publisher=Dover Publications, Inc. | year=1984 | isbn=0-486-40925-2}} * {{cite book | author=George B. Rybicki, Alan P. Lightman | title=Radiative Processes in Astrophysics | publisher=Wiley-Interscience | year=1985 | isbn=0-471-82759-2}} * {{cite book | author=G. E. Thomas and K. Stamnes | title=Radiative Transfer in the Atmosphere and Ocean | publisher=Cambridge University Press. | year=1999 | isbn=0-521-40124-0}} == References == <references/> [[Category:Radiometry]] [[Category:Electromagnetic radiation]] [[Category:atmospheric radiation]] [[de:Strahlungstransport]] [[fr:Transfert de rayonnement]] [[it:Trasferimento radiativo]] [[pl:Transport promieniowania]] [[vi:Phương trình truyền xạ]]