Radiative transfer equation for graded index medium in cylindrical and spherical coordinate systems

In graded index medium, the ray goes along a curved path determined by Fermat principle, and the curved ray-tracing is very difficult and complex. To avoid the complicated and time-consuming computation of curved ray trajectory, the methods not based on ray-tracing technique need to be developed for the solution of radiative transfer in graded index medium. For this purpose, in this paper the streaming operator along a curved ray trajectory in original radiative transfer equation for graded index medium is transformed and expressed in spatial and angular ordinates and the radiative transfer equation for graded index medium in cylindrical and spherical coordinate systems are derived. The conservative and the non-conservative forms of radiative transfer equation for three-dimensional graded index medium are given, which can be used as base equations to develop the numerical simulation methods, such as finite volume method, discrete ordinates method, and finite element method, for radiative transfer in graded index medium in cylindrical and spherical coordinate systems.     

The use of polarized coordinates for solving the macroscopic equations of collisionless plasma

A method is developed for the iterative solution of an infinite system of moment equations derived from the collisionless Boltzmann (Vlasov) equation. Since in a strong magnetic field the dominant terms of these equations have a complicated form, the iterative solution after a few steps comes up against difficulties connected with the inversion of matrices of high order. It is found that the transformation of moment equations into the natural polarized coordinate system (Buneman O.: Phys. Fluids4 (1961), 669) diagonalizes the dominant matrices in the magnetic fields with straight lines of force exactly and in adiabatic magnetic fields approximately. In addition, this transformation reveals the symmetries of the moment equations which in some cases permit an exact solution to be found. Apart from the first three moment equations (the equation of continuity, the dynamic (Euler) equation, the equation for pressure) the general moment equation of then-th order is also derived in the Cartesian and general curvilinear coordinate system and in the corresponding polarized natural coordinate systems. The paper deals with some special cases of natural polarized coordinate systems belonging to a magnetic field with straight lines of force (plane, cylindrical and spherical geometry) and with curved lines of force (magnetic trap with mirrors).     

笛卡儿坐标下空间转子体系的双波函数描述

本文给出了笛卡儿坐标下空间转子体系的双波函数描述,得到了该坐标下每一个力学量的时间演化方程。因而我们的描述是完备的。经典力学运动方程是我们所得演化方程的经典极限。而通常的量子力学描述是我们描述的统计结果。本文还表明,笛卡儿坐标比球坐标能提供更多的物理内容。 关键词：     

Pomraning-Eddington approximation for radiative transfer in a spherical turbid medium

The Pomraning-Eddington approximation is used to solve the radiative transfer problem for anisotropic scattering in a spherical homogeneous turbid medium with diffuse and specular reflecting boundaries. This approximation replaces the radiative transfer integro-differential equation by a second-order differential equation which has an analytical solution in terms of the modified Bessel function. Here, we calculate the partial heat flux at the boundary of anisotropic scattering on a homogeneous solid sphere. The calculations are carried out for spherical media of radii 0.1, 1.0 and 10 mfp and for scattering albedos between 0.1 and 1.0. In addition, the calculations are given for media with transparent, diffuse reflecting and diffuse and specular reflecting boundaries. Two different weight functions are used to verify the boundary conditions. Our results are compared with those given by the Galerkin technique and show greater accuracy for thick and highly scattering media.     

Pomraning-Eddington approximation for radiative transfer in a spherical turbid medium

Abstract

The Pomraning-Eddington approximation is used to solve the radiative transfer problem for anisotropic scattering in a spherical homogeneous turbid medium with diffuse and specular reflecting boundaries. This approximation replaces the radiative transfer integro-differential equation by a second-order differential equation which has an analytical solution in terms of the modified Bessel function. Here, we calculate the partial heat flux at the boundary of anisotropic scattering on a homogeneous solid sphere. The calculations are carried out for spherical media of radii 0.1, 1.0 and 10 mfp and for scattering albedos between 0.1 and 1.0. In addition, the calculations are given for media with transparent, diffuse reflecting and diffuse and specular reflecting boundaries. Two different weight functions are used to verify the boundary conditions. Our results are compared with those given by the Galerkin technique and show greater accuracy for thick and highly scattering media.     

On vector radiative transfer equation in curvilinear coordinate systems

The differential operator of polarized radiative transfer equation is examined in case of homogeneous medium in Euclidean three-dimensional space with arbitrary curvilinear coordinate system defined in it. This study shows that an apparent rotation of polarization plane along the light ray with respect to the chosen reference plane for Stokes parameters generally takes place, due to purely geometric reasons. Analytic expressions for the differential operator of transfer equation dependent on the components of metric tensor and their derivatives are found, and the derivation of differential operator of polarized radiative transfer equation has been made a standard procedure. Considerable simplifications take place if the coordinate system is orthogonal.     

Influence of thermophoresis on heat and mass transfer under non-Darcy MHD mixed convection along a vertical flat plate embedded in a porous medium in the presence of radiation

The paper presents an investigation of the influence of thermophoresis on MHD mixed convective heat and mass transfer of a viscous, incompressible and electrically conducting fluid along a vertical flat plate with radiation effects. The plate is permeable and embedded in a porous medium. To describe the deviation from the Darcy model the Forchheimer flow model is used. The Rosseland approximation is used to describe the radiative heat flux in the energy equation. The governing partial differential equations are transformed into a system of ordinary differential equations using similarity transformation. The nonlinear ordinary differential equations are linearized by using quasilinearization technique and then solved numerically by using implicit finite difference scheme. The numerical results are analyzed for the effects of various physical parameters such as magnetic parameter Ha, mixed convection parameter Ra _d /Pe _d , Reynolds number Red, radiation parameter R, thermophoretic parameter τ, Prandtl number Pr, and Schmidt number Sc. The heat transfer coefficient is also tabulated for different values of physical parameters.     

Theoretical study of combined radiative conductive and convective heat transfer in semi-transparent porous medium in a spherical enclosure

A new method for the solution of the radiative transfer equation in spherical media based on a modified discrete ordinates method is extended to study radiative, conductive and convective heat transfer in a semi-transparent scattering porous medium. The set of differential equations is solved using the fourth-order Runge-Kutta method. Various results are obtained for the case of combined radiative and conductive heat transfer, as well as for the interaction of those modes with convection. The effects of some radiative properties of the medium on the heat transfer rate are examined.     

Elliptic PDE formulation and boundary conditions of the spherical harmonics method of arbitrary order for general three-dimensional geometries

The inherent complexity of the radiative transfer equation makes the exact treatment of radiative heat transfer impossible even for idealized situations and simple boundary conditions. Therefore, a wide variety of efficient solution methods have been developed for the RTE. Among these solution methods the spherical harmonics method, the moment method, and the discrete ordinates method provide means to obtain higher-order approximate solutions to the equation of radiative transfer. Although the assembly of the governing equations for the spherical harmonics method requires tedious algebra, their final form promises great accuracy for any given order, since it is a spectral method (rather than finite difference/finite volume in the case of discrete ordinates). In this study, a new methodology outlined in a previous paper on the spherical harmonics method (P_N) is further developed. The new methodology employs successive elimination of spherical harmonic tensors, thus reducing the number of first-order partial differential equations needed to be solved simultaneously by previous P_N approximations (=(N+1)²). The result is a relatively small set (=N(N+1)/2) of second-order, elliptic partial differential equations, which can be solved with standard PDE solution packages. General boundary conditions and supplementary conditions using rotation of spherical harmonics in terms of local coordinates are formulated for the general P_N approximation for arbitrary three-dimensional geometries. Accuracy of the P_N approximation can be further improved by applying the “modified differential approximation” approach first developed for the P₁-approximation. Numerical computations are carried out with the P₃ approximation for several new two-dimensional problems with emitting, absorbing, and scattering media. Results are compared to Monte Carlo solutions and discrete ordinates simulations and a discussion of ray effects and false scattering is provided.     

New applications of the homogeneous balance principle

The homogeneous balance principle has been widely applied to the exploration of nonlinear transformation, exact solutions (especially solitary wave solution), dromion and similarity reduction to the nonlinear partial differential equations in mathematical physics. In this paper, we use the homogeneous balance principle to derive B?cklund transformations for nonlinear partial differential equations that have more nonlinear terms and more highest-order partial derivative terms. With the aid of the B?cklund transformations derived here, we could obtain exact solutions to the nonlinear partial differential equations. The Davey－Stewartson equation and the Nizhnik－Novikov－Veselov equation are considered as the examples.     

An Hermitian method for the solution of radiative transfer problems

A new method for the solution of the radiative transfer equation is developed. The resulting system, like the ordinary difference equations is tridiagonal in form and easy to use. Numerical examples are given which show the superiority of the new technique to previous differential and integral equation methods.     

Tau method approximation for radiative transfer problems in a slab medium

An approximate method for solving the radiative transfer equation in a slab medium with an isotropic scattering is proposed. The method is based upon constructing the double Legendre series to approximate the required solution using Legendre tau method. The differential and integral expressions which arise in the radiative transfer equation are converted into a system of linear algebraic equations which can be solved for the unknown coefficients. Numerical examples are included to demonstrate the validity and applicability of the method and a comparison is made with existing results.     

Numerical solution of radiative and conductive heat transfer in concentric spherical and cylindrical media

A new technique is presented to improve the performance of the discrete ordinates method when solving the coupled conduction-radiation problems in spherical and cylindrical media. In this approach the angular derivative term of the discretized one-dimensional radiative transfer equation is derived from an expansion of the radiative intensity on the basis of Chebyshev polynomials. The set of resulting differential equations, obtained by the application of the S_N method, is numerically solved using the boundary value problem with the finite difference algorithm. Results are presented for the different independent parameters. Numerical results obtained using the Chebyshev transform method compare well with the benchmark approximate solutions. Moreover, the new technique can easily be applied to higher-order S_N calculations.     

New Solutions of Three Nonlinear Space- and Time-Fractional Partial Differential Equations in Mathematical Physics

Motivated by the widely used ansätz method and starting from the modified Riemann-Liouville derivative together with a fractional complex transformation that can be utilized to transform nonlinear fractional partial differential equations to nonlinear ordinary differential equations, new types of exact traveling wave solutions to three important nonlinear space- and time-fractional partial differential equations are obtained simultaneously in terms of solutions of a Riccati equation. The results are new and first reported in this paper.     

Boundary conditions for differential approximations

Low order spherical harmonic (P-N) approximations are applied to a radiative transfer Marshak wave problem. A modified Milne boundary condition is developed for the P-2 approximation, similar to one suggested earlier for the P-1 approximation. Comparison with exact Monte Carlo results suggests that this modified P-2 method may be an accurate and generally applicable differential approximation to the equation of transfer. The Monte Carlo results presented should be useful for testing other approximate formulations of radiative transfer and validating time dependent numerical solution methods for the equation of transfer.     

Nonlinear,nondispersive wave equations: Lagrangian and Hamiltonian functions in the hodograph transformation

The hodograph transformation is generally used in order to associate a system of linear partial differential equations to a system of nonlinear (quasilinear) differential equations by interchanging dependent and independent variables. Here we consider the case when the nonlinear differential system can be derived from a Lagrangian density and revisit the hodograph transformation within the formalism of the Lagrangian-Hamiltonian continuous dynamical systems.Restricting to the case of nondissipative, nondispersive one-dimensional waves, we show that the hodograph transformation leads to a linear partial differential equation for an unknown function that plays the role of the Lagrangian in the hodograph variables. We then define the corresponding hodograph Hamiltonian and show that it turns out to coincide with the wave amplitude. i.e., with the unknown function of the independent variables to be solved for in the initial nonlinear wave equation.     

Radiative heat transfer between two concentric spheres separated by a two-phase mixture of non-gray gas and particles using the modified discrete-ordinates method

The radiative heat transfer between two concentric spheres separated by a two-phase mixture of non-gray gas and a cloud of particles is investigated by using the combined finite-volume and discrete-ordinates method, named modified discrete-ordinates method (MDOM), which integrates the radiative transfer equation (RTE) over a control volume and a control angle simultaneously like in the finite-volume method (FVM) and treats the angular derivative terms due to spherical geometry as the conventional discrete-ordinates method (DOM). The radiative properties involving non-gray gas and particle behavior are modeled by using the extended weighted sum of gray gases model (WSGGM) with particles. Mathematical formulation and final discretization equations for the RTE are introduced by considering the behavior of a two-phase mixture of non-gray gas and particles in a spherically symmetric concentric enclosure. The present approach is validated by comparing with the results of previous works including gray and non-gray radiative heat transfer. Finally, a detailed investigation of the radiative heat transfer with non-gray gases and/or a two-phase mixture is conducted to examine the dependence of the radiative heat transfer upon temperature ratio between inner and outer spherical enclosure, particle concentration, and particle temperature.     

Photon conservation in scattering by large ice crystals with the SASKTRAN radiative transfer model

The scattering of visible light by ice crystals and dust in radiative transfer models is challenging in part due to the large amount of scattering in the forward direction. We introduce a technique that ensures numerical conservation of photons in any radiative transfer model and that quantifies the integration error associated with highly asymmetric phase functions. When applied to a successive-orders of scatter model, the technique illustrates the high accuracy obtained in numerical integration of molecular and aerosol scattering. As well, a phase function truncation and renormalization technique is applied to scattering by ice crystals with very large size parameters, between 100 and 1000, and the scaled radiative transfer equation is solved with the spherical successive-orders model, SASKTRAN. Since computations shown this work are performed in a fully spherical model atmosphere, the computed radiances are not subject to the discontinuity at the horizon that is inherent in models using a plane–parallel assumption. The methods introduced in this work are of particular interest in modeling limb radiances in the presence of thin cirrus clouds.     

Gröbner bases applied to systems of linear difference equations

In this paper we consider systems of partial (multidimensional) linear difference equations. Specifically, such systems arise in scientific computing under discretization of linear partial differential equations and in computational high-energy physics as recurrence relations for multiloop Feynman integrals. The most universal algorithmic tool for investigation of linear difference systems is based on their transformation into an equivalent basis form. We present an algorithm for this transformation implemented in Maple. The algorithm and its implementation can be applied to automatic generation of difference schemes for linear partial differential equations and to reduction of Feynman integrals. Some illustrative examples are given. The text was submitted by the author in English.