On a general theory of characteristics and the method of invariant imbedding
Invariant imbedding equations for transport problems derived and applied to boundary value and initial value problems
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Invariant imbedding equations for transport problems derived and applied to boundary value and initial value problems
Invariant imbedding method applied to energy dependent neutron radiation penetration problems in thick and heterogeneous shields
Invariant imbedding equations of transport theory, noting imbedding of boundary value problem into initial value problems leads to Cauchy problem
Radiation shield penetration calculation from invariant imbedding method
The invariant imbedding method considered is based on an equation which describes the change in the reflected radiation when an optically thin layer is added to the top of the atmosphere. The equation is used to treat the problem of reflection from a planetary atmosphere as an initial value problem. A fast method is discussed for the solution of the invariant imbedding equation. The speed and accuracy of the new method are illustrated by comparing it with the doubling program published by Hansen and Travis (1974). Computations are performed of the equivalent widths of carbon dioxide absorption lines in solar radiation reflected by Venus for several models of the planetary atmosphere.
Invariant imbedding and perturbation techniques applied to diffuse reflection from spherical shells
The extended boundary condition method (EBCM) and invariant imbedding method (IIM) are two fundamentally different T-matrix methods for the solution of light scattering by nonspherical particles. The standard EBCM is very efficient but encounters a loss of precision when the particle size is large, the maximum size being sensitive to the particle aspect ratio. The IIM can be applied to particles in a relatively large size parameter range but requires extensive computational time due to the number of spherical layers in the particle volume discretization. A numerical combination of the EBCM and the IIM (hereafter, the EBCM+IIM) is proposed to overcome the aforementioned disadvantages of each method. Even though the EBCM can fail to obtain the T-matrix of a considered particle, it is valuable for decreasing the computational domain (i.e., the number of spherical layers) of the IIM by providing the initial T-matrix associated with an iterative procedure in the IIM. The EBCM+IIM is demonstrated to be more efficient than the IIM in obtaining the optical properties of large size parameter particles beyond the convergence limit of the EBCM. The numerical performance of the EBCM+IIM is illustrated through representative calculations in spheroidal and cylindrical particle cases.
Invariant imbedding technique applied to electromagnetic wave propagation through inhomogeneous reentry plasmas
Suboptimal closed loop control of nonlinear systems subject to quadratic performance indices by invariant imbedding concepts and maximum principle
Three terms, ''Waterman's T-matrix method'', ''extended boundary condition method (EBCM)'', and ''null field method'', have been interchangeable in the literature to indicate a method based on surface integral equations to calculate the T-matrix. Unlike the previous method, the invariant imbedding method (IIM) calculates the T-matrix by the use of a volume integral equation. In addition, the standard separation of variables method (SOV) can be applied to compute the T-matrix of a sphere centered at the origin of the coordinate system and having a maximal radius such that the sphere remains inscribed within a nonspherical particle. This study explores the feasibility of a numerical combination of the IIM and the SOV, hereafter referred to as the IIMþSOV method, for computing the single-scattering properties of nonspherical dielectric particles, which are, in general, inhomogeneous. The IIMþSOV method is shown to be capable of solving light-scattering problems for large nonspherical particles where the standard EBCM fails to converge. The IIMþSOV method is flexible and applicable to inhomogeneous particles and aggregated nonspherical particles (overlapped circumscribed spheres) representing a challenge to the standard superposition T-matrix method. The IIMþSOV computational program, developed in this study, is validated against EBCM simulated spheroid and cylinder cases with excellent numerical agreement (up to four decimal places). In addition, solutions for cylinders with large aspect ratios, inhomogeneous particles, and two-particle systems are compared with results from discrete dipole approximation (DDA) computations, and comparisons with the improved geometric-optics method (IGOM) are found to be quite encouraging.
The paper studies thermal microwave emission from an inhomogeneous slab of a random medium, with possible nonuniform absorption, scattering, and temperature profiles, bounded by different dielectrics on both sides. The invariant imbedding method is used to cast the boundary value problem of the radiative transfer equations into an initial value problem at zero slab thickness. As a numerical example, the angular and polarization variations of brightness temperatures for ice over water are considered.
The single-scattering properties (SSPs) of hydrometeors are the fundamental quantities for physics-based precipitation retrievals. Thus, efficient computation of their electromagnetic scattering is of great value. Whereas the semi-analytical T-Matrix methods are likely the most efficient for nonspherical hydrometeors with axial symmetry, they are not suitable for arbitrarily shaped hydrometeors absent of any significant symmetry, for which volume integral methods such as those based on Discrete Dipole Approximation (DDA) are required. Currently the two leading T-matrix methods are the Extended Boundary Condition Method (EBCM) and the Invariant Imbedding T-matrix Method incorporating Lorentz-Mie Separation of Variables (IITM+SOV). EBCM is known to outperform IITM+SOV for hydrometeors with modest aspect ratios. However, in cases when aspect ratios become extreme, such as needle-like particles with large height to diameter values, EBCM fails to converge. Such hydrometeors with extreme aspect ratios are known to be present in solid precipitation and their SSPs are required to model the radiative responses accurately. In these cases, IITM+SOV is shown to converge. An efficient, parallelized C++ implementation for both EBCM and IITM+SOV has been developed to conduct a performance comparison between EBCM, IITM+SOV, and DDSCAT (a popular implementation of DDA). We present the comparison results and discuss details. Our intent is to release the combined ECBM IITM+SOV software to the community under an open source license.
Differential-integral equation for dissipation function and derivation of conservation relationship connecting reflection, transmission and dissipation functions of finite slab
Spectral line formation, assuming frequency redistribution for plane parallel stellar atmosphere containing nonuniform distribution of internal emission sources
A matrix Fredholm integral equation of neuronal networks is transformed into a Cauchy system suited for numerical and analytical studies. A special case is discussed, and a connection with the classical renewal integral equation of stochastic point processes is presented.
The invariant imbedding technique is applied to the problem of resonance line scattering with frequency redistribution. By applying the principles of invariance to the resonance line scattering problem a set of non-linear, integrodifferential equations for the scattering and transmission functions are derived. Thus, an alternate approach to the resonance line scattering problem is presented. Rather than solving for the specific intensity directly the invariant imbedding equations may be solved for the scattering and transmission functions from which the specific intensity may be calculated. The technique is very general and in principle, can accommodate scattering in inhomogeneous atmospheres with arbitrary frequency redistribution and scattering phase functions. Solutions to the invariant imbedding equations in some simple cases are presented. If the scattering is isotropic and completely noncoherent and the atmosphere may be modeled as homogeneous, the scattering and transmission functions may be expressed in terms of generalized H, X and Y functions. This technique is extended to the problems of coupled spectral lines and scattering of radiation from internal sources.
To obtain an exact solution of the radiative-transfer equation in media where both absorption and scattering are significant, the usual approach is to use a numerical method. Three methods are known in the literature: invariant imbedding, eigenvalue-eigenfunction, and matrix doubling. This paper examines the practical application of these methods to the problem of emission from an inhomogeneous (Rayleigh) layer, the effects of layer parameters on the stability. It is found that invariant imbedding is most suitable for computing emission from an inhomogeneous layer with a temperature profile but tends to be unstable as the optical thickness of the layer increases beyond 0.5. On the other hand, the matrix-doubling method is stable for arbitrary optical thickness but is not suitable for handling multilayers. The eigenvalue-eigenfunction method is more stable than the invariant imbedding as optical thickness increases up to 2.0. It also permits temperature profile in the layer, but the computation is much more complicated. It is less stable than the matrix-doubling method when optical thickness is larger than 2.0. In general, the choice of a method is dependent on the nature of the problem.
Book on electromagnetic wave propagation and turbulent media covering invariant imbedding method, turbulence generation, statistical methods of analysis, etc