Search NASA⌕ Search

SEARCH · Search NASA

Results for “scalar scattering”

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.

82 records · Page 5

Numerical Simulations of Light Bullets, Using the Full Vector, Time Dependent, Nonlinear Maxwell Equations

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that are currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Kerr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations.

Goorjian, Peter M.↗

Numerical Simulations of Light Bullets, Using The Full Vector, Time Dependent, Nonlinear Maxwell Equations

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that we currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Karr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations.

Goorjian, Peter M.↗

Numerical Simulations of Light Bullets, Using The Full Vector, Time Dependent, Nonlinear Maxwell Equations

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that are currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Kerr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations.

Goorjian, Peter M.↗

Obtaining the Grobner Initialization for the Ground Flash Fraction Retrieval Algorithm

At optical wavelengths and from the vantage point of space, the multiple scattering cloud medium obscures one's view and prevents one from easily determining what flashes strike the ground. However, recent investigations have made some progress examining the (easier, but still difficult) problem of estimating the ground flash fraction in a set of N flashes observed from space In the study by Koshak, a Bayesian inversion method was introduced for retrieving the fraction of ground flashes in a set of flashes observed from a (low earth orbiting or geostationary) satellite lightning imager. The method employed a constrained mixed exponential distribution model to describe the lightning optical measurements. To obtain the optimum model parameters, a scalar function of three variables (one of which is the ground flash fraction) was minimized by a numerical method. This method has formed the basis of a Ground Flash Fraction Retrieval Algorithm (GoFFRA) that is being tested as part of GOES-R GLM risk reduction.

Solakiewicz, R.↗

Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. I: Approach Based on the Far-Field Foldy Equations

In this paper, we revisit, with further enhancements and clarifications, the self-consistent first-principles approach developed previously for deriving the vector radiative transfer theory for a discrete random medium with a sparse concentration of particles. We specifically consider the case of a plane-parallel particulate layer embedded in an otherwise homogeneous unbounded medium. The solution method is based on the far-field Foldy equations, an order-of-scattering expansion for the total field derived under the Twersky approximation, the computation of the coherent field by assuming that the positions of the particles are uncorrelated, and the ladder approximation for the coherency dyadic. The latter yields an integral equation for the diffuse specific coherency dyadic, defined through an angular spectrum repre- sentation for the coherency dyadic, which in turn, gives the vector radiative transfer equation for the diffuse specific intensity column vector. We analyze specifically the computation of the coherent field for inhomogeneous particulate media and multiple species of particles, the continuous extension of the far-field representation to the near field, the Foldy approximation, and the Foldy integral equation for the coherent field. Finally, we discuss the transition from the vector to the scalar radiative transfer equation.

Doicu, Adrian↗

Drift shell splitting by internal geomagnetic multipoles.

Computations on an 80-coefficient model of the earth's field illustrate the 'topography' of the magnetic equatorial surface and the geometry of the drift shells of geomagnetically trapped particles. Individual terms in the spherical harmonic expansion of the geomagnetic scalar potential V(r, theta, phi) are either even or odd in cos theta, where theta = 90 deg denotes the dipole equator. Terms that are even in cos theta tend to 'warp' the equatorial surface, but do not (in first order) distort particle drift shells radially nor split the drift shells of particles having different equatorial pitch angles. Azimuthally asymmetric terms that are odd in cos theta do cause shell splitting in first order. Shell splitting at large L values (neglecting deformation of the earth's field by the solar wind) is found to be dominated by the geomagnetic octupole. At L approximately equal to 1, shell splitting is strongly enhanced by the South American and South African anomalies. When combined with pitch angle diffusion caused by atmospheric scattering, these results may be able to account for anomalous radial diffusion of inner zone electrons.

Roederer, J. G.↗

Investigation of Strain/Vorticity and Large-Scale Flow Structure in Turbulent Nonpremixed Jet Flames

Our study will use the microgravity environment to investigate the underlying flow structure of turbulent nonpremixed round jet flames. In particular, we aim to investigate the large-scale turbulent structure using planar laser Mie scattering (PLMS), and the strain rate and vorticity fields using particle image velocimetry (PIV). This work is motivated by recent studies in our laboratory that have led to several interesting observations of nominally momentum-driven turbulent nonpremixed planar flames. First of all, the organized large-scale turbulent structures that are observed in nonreacting planar jets may be substantially modified or suppressed in nonpremixed planar jet flames. Furthermore, a recent study using PIV and planar laser-induced fluorescence of OH has shown that in transitional and turbulent nonpremixed planar jet flames the presence of the flame seems to greatly influence the underlying vorticity and strain fields, as compared to nonreacting jets. For example, the reaction zones in the jet flames are strongly correlated with regions of high vorticity. A related study has demonstrated that vorticity is not correlated in the same way with either iso-scalar surfaces or scalar dissipation layers in nonreacting planar jets. Furthermore, the relationship between strain and the reaction zone appears to be modified by the presence of high levels of heat release. In particular, the strain rate field in planar jet flames exhibits a preferred direction of principal compressive strain that apparently is related to strong shear across the reaction zone. This preferred direction of strain was not observed in nonreacting jets. One of the major problems encountered when conducting these types of studies is that it is difficult to know to what extent buoyancy influences the results. Therefore, the microgravity environment provides us with an excellent opportunity to explore these issues without the complicating effects of buoyancy. This is particularly the case when studying flames that are transitional between laminar and turbulent states. For example, the strong correlation of vorticity with the reaction zone (discussed above) was observed in both transitional and turbulent planar flames, but the effect was stronger for the transitional case. To date, the reason for the presence of the vorticity-reaction zone correlation is not known, although vorticity production via baroclinic torque is a likely cause. The microgravity environment will allow us to specifically determine whether the vorticity is produced by baroclinic torque resulting from the flame density gradient acting with the hydrostatic pressure gradient.AAdditional details of the planned experiments are provided.

Clemens, N. T.↗

Radiative Transfer Modeling of Phytoplankton Fluorescence Quenching Processes

We report the first radiative transfer model that is able to simulate phytoplankton fluorescence with both photo chemical and non-photo chemical quenching included. The fluorescence source term in the inelastic radiative transfer equation is proportional to both the quantum yield and scalar irradiance at excitation wavelengths. The photo chemical and non photo chemical quenching processes change the quantum yield based on the photosynthetic active radiation. A sensitivity study was performed to demonstrate the dependence of the fluorescence signal on chlorophyll a concentration, aerosol optical depths and solar zenith angles. This work enables us to better model the phytoplankton fluorescence, which can be used in the design of new space-based sensors that can provide sufficient sensitivity to detect the phytoplankton fluorescence signal. It could also lead to more accurate remote sensing algorithms for the study of phytoplankton physiology.

Inelastic scattering↗

Electromagnetic Scattering by Discrete Random Media. I: The Dispersion Equation and the Configuration-Averaged Exciting Field

We consider the scattering of a plane electromagnetic wave obliquely incident on a plane-parallel layer of discrete random medium with non-scattering boundaries. We solve the Lax integral equation for the conditional configuration-averaged exciting field coefficients by assuming a special-form solution, that is, by representing the conditional configuration-averaged exciting field coefficients as a linear combination of the coefficients corresponding to an up-going and a down-going wave. This solution representation is supposed to be valid within the whole domain occupied by the particles, even in the close proximity of the boundaries. By balancing the waves with different propagation directions and wavenumbers we derive two homogeneous systems of equations corresponding to the generalized Lorenz–Lorentz law and two inhomogeneous systems of equations corresponding to the generalized Ewald–Oseen extinction theorem. It is shown that (i) the two homogeneous systems of equations of the generalized Lorenz–Lorentz law reduce to a single homogeneous system of equations corresponding to a semi-infinite discrete random medium at normal incidence; (ii) the dispersion equation is direction and polarization independent; and (iii) the two inhomogeneous systems of equations of the generalized Ewald–Oseen extinction theorem can be reduced to two scalar equations by means of the addition theorem for vector spherical harmonics. It is also shown that the same dispersion equation can be obtained without assuming a special-form solution representation in the proximity of the boundaries.

Adrian Doicu↗

Issues and Methods Concerning the Evaluation of Hypersingular and Near-Hypersingular Integrals in BEM Formulations

It is known that higher order modeling of the sources and the geometry in Boundary Element Modeling (BEM) formulations is essential to highly efficient computational electromagnetics. However, in order to achieve the benefits of hIgher order basis and geometry modeling, the singular and near-singular terms arising in BEM formulations must be integrated accurately. In particular, the accurate integration of near-singular terms, which occur when observation points are near but not on source regions of the scattering object, has been considered one of the remaining limitations on the computational efficiency of integral equation methods. The method of singularity subtraction has been used extensively for the evaluation of singular and near-singular terms. Piecewise integration of the source terms in this manner, while manageable for bases of constant and linear orders, becomes unwieldy and prone to error for bases of higher order. Furthermore, we find that the singularity subtraction method is not conducive to object-oriented programming practices, particularly in the context of multiple operators. To extend the capabilities, accuracy, and maintainability of general-purpose codes, the subtraction method is being replaced in favor of the purely numerical quadrature schemes. These schemes employ singularity cancellation methods in which a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. An example of the sin,oularity cancellation approach is the Duffy method, which has two major drawbacks: 1) In the resulting integrand, it produces an angular variation about the singular point that becomes nearly-singular for observation points close to an edge of the parent element, and 2) it appears not to work well when applied to nearly-singular integrals. Recently, the authors have introduced the transformation u(x(prime))= sinh (exp -1) x(prime)/Square root of ((y prime (exp 2))+ z(exp 2) for integrating functions of the form I = Integral of (lambda(r(prime))((e(exp -jkR))/(4 pi R) d D where A (r (prime)) is a vector or scalar basis function and R = Square root of( (x(prime)(exp2) + (y(prime)(exp2) + z(exp 2)) is the distance between source and observation points. This scheme has all of the advantages of the Duffy method while avoiding the disadvantages listed above. In this presentation we will survey similar approaches for handling singular and near-singular terms for kernels with 1/R(exp 2) type behavior, addressing potential pitfalls and offering techniques to efficiently handle special cases.

Fink, P. W.↗