Search NASA⌕ Search

SEARCH · Search NASA

Results for “Legendre coefficients”

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.

At least 19 records

Coefficients of associated Legendre functions

The exact coefficients for the explicit forms of the associated Legendre functions Pm for integer values of m,n=0,1,2,...25 are presented in tabular form together with two cross-referenced listings of the zeroes for these functions rounded to five decimal places. The unfactored coefficients and the interger coefficients are presented in adjacent columns for each function. The greatest common factor and divisor have been removed and listed separately for the integer coefficients.

Peasley, Q. D.↗

Neutron elastic and inelastic scattering differential cross sections on carbon

Elastic and inelastic neutron scattering angular distributions were measured on natural carbon samples to confirm existing experimental data and evaluations in the fast neutron region and to guide improvements in resonance parameters, where needed. Sixty-four (n,n') differential cross section measurements were performed at 45 incident neutron energies between 0.5 and 8.0 MeV. Experimental angle-integrated elastic scattering cross sections are consistent with ENDF/B-VIII.0 values with the exception of the region from 3.2 to 4.0 MeV where our results are ~3% higher. In the 3.4 to 3.6 MeV region our differential cross sections are slightly lower at forward angles and somewhat higher at backward angles than the ENDF/B-VIII.0 calculations; however, our results are consistent with measured data from other research groups. The first- through fourth-order elastic scattering Legendre coefficients from fits to these experimental data are consistent with ENDF/B-VIII.0 values across this range of incident neutron energies. Inelastic scattering cross sections were measured at 12 incident neutron energies between 5.6 and 7.8 MeV. Angle-integrated cross sections agree well with the 1978 Perey and ENDF/B-VIII.0 values; however, the angular-distribution shapes deviate somewhat from the ENDF/B-VIII.0 calculations. These data will allow refinement of the resonance parameter description in this difficult energy region.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evaluation and Validation of the n+ 63,65 Cu Cross Sections [Slides]

This presentation discusses the objective of the research which is to update n+ 63,65 Cu cross section evaluations with recently measured data to resolve discrepancies in benchmark performance. It also discusses the models used which are the R -matrix analysis and the analysis of angular distribution coefficients. Additionally, the validation methods that were used are discussed, including the Rez shielding benchmark and the ICSBEP criticality benchmarks. In conclusion, he n+ 63,65 Cu cross sections have been updated via R-matrix analysis up to 100 keV. An increased average capture cross section and the adoption of experimentally based Legendre coefficients lead to improved performance in reactivity benchmarks. In the fast region, the adoption of the JENDL-4.0 cross sections above 4.0 MeV improves the performance in shielding benchmarks. Ultimately, the n 63,65 Cu ENDF files will be submitted to ENDF/B-VIII.1.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Updates to the n+ 63,65 Cu Angular Distributions [Slides]

The performance of the benchmark suite is very sensitive to changes in the angular distribution data. The quasi-differential measurements performed by Blain et al at RPI provide a valuable constraint. Furthermore, ENDF/B-VIII.0 disagrees with their measurements consistently at 300 keV (the transition between RRR and high energy). The next step is to conservatively adjust the Legendre coefficients near 300 keV, validating performance against both critical benchmarks and the RPI quasi-differential measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Principal Component Analysis of Up-The-Ramp Sampled Infrared Array Data

We describe the results of principal component analysis (PCA) of up-the-ramp sampled infrared (IR) array data from the Hubble Space Telescope wide field camera 3 (WFC3 IR), James Webb Space Telescope NIRSpec, and prototype Wide Field Infrared Survey Telescope’s wide field instrument detectors. These systems use, respectively, Teledyne H1R, H2RG, and H4RG-10 near-IR detector arrays with a variety of IR array controllers. The PCA shows that the Legendre polynomials approximate the principal components of these systems (i.e., they roughly diagonalize the covariance matrix). In contrast to the monomial basis that is widely used for polynomial fitting and linearization today, the Legendre polynomials are an orthonormal basis. They provide a quantifiable, compact, and (nearly) linearly uncorrelated representation of the information content of the data. By fitting a few Legendre polynomials, nearly all of the meaningful information in representative WFC3 astronomical datacubes can be condensed from 15 up-the-ramp samples down to 6 compressible Legendre coefficients per pixel. The higher order coefficients contain time domain information that is lost when one projects up-the-ramp sampled datacubes onto two-dimensional images by fitting a straight line, even if the data are linearized before fitting the line. Going forward, we believe that this time domain information is potentially important for disentangling the various nonlinearities that can affect IR array observations, i.e., inherent pixel nonlinearity, persistence, burn in, brighter-fatter effect, (potentially) nonlinear interpixel capacitance, and perhaps others.

Rauscher, Bernard J.↗

Updates and Validation for the n+ 63,65 Cu Cross Sections [Abstract]

The neutron induced total, elastic, and capture cross sections of 63,65 Cu isotopes were selected for evaluation in the resolved and unresolved resonance energy ranges by the National Criticality Safety Program to resolve discrepancies related to benchmark performance. This is especially evident for the series of ZEUS benchmarks in which copper is used as a reflector. Because copper is also used as structural material in both fission and fusion reactors, the need to address benchmark discrepancies linked to nuclear data deficiencies is a task of primary importance. The aim of this work is to describe the steps of evaluation work towards a consistent improvement of the benchmark performance. The R-matrix analysis with the SAMMY code focused on the 63 Cu(n,γ) reaction channel between 100-300 keV coupled to unresolved resonance region parameters up to 650 keV to fit average cross section data from a recent experiment. Due to the high sensitivity of many benchmarks to elastic scattering angular distribution data, especially for the 65 Cu isotope, the impact of these data was tested by generating Legendre coefficients from both resonance parameters and the Hauser-Feshbach model. Guided by the findings of Shaw et al., the performance of the current evaluation for 65 Cu was compared to that of ENDF/B-VII.1 and ENDF/B-VIII.0 by testing the reactivity coefficients corresponding to the validation suite of experimental criticality benchmarks for thermal, intermediate, and fast systems taken from the International Criticality Safety Benchmark Experiments Project Handbook. The benchmark performance is especially sensitive to 63 Cu(n,γ) and 65 Cu elastic scattering for neutron energies in the 100–500 keV region, whereas 100 keV is the upper limit of the resolved resonance region in the ENDF/B-VIII.0 evaluations for 63,65 Cu. The results highlight the need to handle the transition from the resolved resonance region to the high energy region carefully.

07 ISOTOPE AND RADIATION SOURCES↗

A comparison of past and present computational methods for shape analysis of double-shell x-ray radiographs

Implosion symmetry is a key requirement in achieving a robust burning plasma in inertial confinement fusion experiments. In double-shell capsule implosions, we are interested in the shape of the inner shell as it pushes on the fuel. Shape analysis is a popular technique for studying said symmetry during implosion. Combinations of filtering and contour-finding algorithms are studied for their promise in reliably recovering Legendre shape coefficients from synthetic radiographs of double-shell capsules with applied levels of noise. A radial lineout max(slope) method when used on an image pre-filtered with non-local means and a variant of the marching squares algorithm are able to recover p 0 , p 2 , and p 4 maxslope Legendre shape coefficients with mean pixel discrepancy errors of 2.81 and 3.06, respectively, for the noisy synthetic radiographs we consider. Here, this improves upon prior radial lineout methods paired with Gaussian filtering, which we show to be unreliable and whose performance is dependent on input parameters that are difficult to estimate.

47 OTHER INSTRUMENTATION↗

Phase function, backscatter, extinction, and absorption for standard radiation atmosphere and El Chichon aerosol models at visible and near-infrared wavelengths

Tabular values of phase function, Legendre polynominal coefficients, 180 deg backscatter, and extinction cross section are given for eight wavelengths in the atmospheric windows between 0.4 and 2.2 microns. Also included are single scattering albedo, asymmetry factor, and refractive indices. These values are based on Mie theory calculations for the standard rediation atmospheres (continental, maritime, urban, unperturbed stratospheric, volcanic, upper atmospheric, soot, oceanic, dust, and water-soluble) assest measured volcanic aerosols at several time intervals following the El Chichon eruption. Comparisons of extinction to 180 deg backscatter for different aerosol models are presented and related to lidar data.

Whitlock, C. H.↗

Helioseismology from the South Pole: Comparison of 1987 and 1981 Results

Full disk images with 10 arc sec pixels and filtered to a 7 A pass band centered on the Ca II K line were obtained from the geographic South Pole in 1981 and 1987. In 1981, 50hr of essentially uninterrupted data were obtained. In 1987, three such runs were obtained over a period of 325 hours for a duty cycle of 47 percent. The 1987 observations are characterized by a much lower level of solar activity than 1981, a much improved CCD camera, considerably better image stability and a varying amount of instrumental scatter. The 1987 data have a substantially better signal-to-noise ratio than the 1981 data so that oscillations with degrees from 0 to 150 and frequencies from 2 to 7 mHz are well observed. The observations were reduced to spectra in l, m, and v. A comparison of p-mode frequencies measured in 1981 and 1987, and coefficients of Legendre polynomial expansions of frequency shifts caused by solar rotation are presented. The time behavior of systematic frequency shifts which depend upon m but which do not arise from rotation is described.

Jefferies, S. M.↗

Considerations for modeling thin cirrus effects via brightness temperature differences

Brightness temperature difference (BTD) values are calculated for selected Geostationary Operational Environmental Satellite (GOES-6) channels (3.9, 12.7 micrometer) and Advanced Very High Resolution Radiometer channels (3.7, 12.0 micrometer). Daytime and nighttime discrimination of particle size information is possible given the infrared cloud extinction optical depth and the BTD value. BTD values are presented and compared for cirrus clouds composed of equivalent ice spheres (volume, surface area) versus randomly oriented hexagonal ice crystals. The effect of the hexagonal ice crystals is to increase the magnitude of the BTD values calculated relative to equivalent ice sphere (volume, surface area) BTDs. Equivalent spheres (volume or surface area) do not do a very good job of modeling hexagonal ice crystal effects on BTDs; however, the use of composite spheres improves the simulation and offers interesting prospects. Careful consideration of the number of Legendre polynomial coefficients used to fit the scattering phase functions is crucial to realistic modeling of cirrus BTDs. Surface and view-angle effects are incorporated to provide more realistic simulation.

Schmidt, E. O.↗

The Gravitational Potential Due to Uniform Disks and Rings

The gravitational potential of bodies possessing axial symmetry can be expressed as a power series in distance, with the Legendre polynomials as coefficients. Such series, however, converge so slowly in the neighborhood of thin, uniform disks and rings that too many series terms must be summed in order to obtain an accurate field measure. A gravitational potential expression is presently obtained in closed form, in terms of complete elliptic integrals.

Lass, H.↗

Direct solution for thermal stresses in a nose cap under an arbitrary axisymmetric temperature distribution

The design of a nose cap for a hypersonic vehicle is an iterative process requiring a rapid, easy to use and accurate stress analysis. The objective of this paper is to develop such a stress analysis technique from a direct solution of the thermal stress equations for a spherical shell. The nose cap structure is treated as a thin spherical shell with an axisymmetric temperature distribution. The governing differential equations are solved by expressing the stress solution to the thermoelastic equations in terms of a series of derivatives of the Legendre polynomials. The process of finding the coefficients for the series solution in terms of the temperature distribution is generalized by expressing the temperature along the shell and through the thickness as a polynomial in the spherical angle coordinate. Under this generalization the orthogonality property of the Legendre polynomials leads to a sequence of integrals involving powers of the spherical shell coordinate times the derivative of the Legendre polynomials. The coefficients of the temperature polynomial appear outside of these integrals. Thus, the integrals are evaluated only once and their values tabulated for use with any arbitrary polynomial temperature distribution.

Davis, Randall C.↗

Numerical experiments in geomagnetic modeling

Numerical tests were made, using least squares fitting of a spherical harmonic model, to a selection of Magsat data to determine the practical limits of this technique with modern computers. The resulting (M102189) model, whose coefficients were adjusted up to n = 50, was compared with M07AV6, a previous model which used least squares (on vector data) for coefficients up to n = 29, and Gauss-Legendre quadrature (on Z residuals) to adjust the coefficients up to n = 63. For the new least squares adjustment to n = 50 a condition number of 115 was obtained for the solution matrix, with a resulting precision of 11 significant figures. The M102189 model shows a lower and more Gaussian residual distribution than did M07AV6, though the Gaussian envelope fits to the residual distributions, even for the scalar field, gives "standard deviations' never lower than 6 nT, a factor of three higher than the estimated Magsat observational errors. Ionospheric currents are noted to have a significant effect on the coefficients of the internal potential functions.

Cain, Joseph C.↗

Legendre-Tau approximations for functional differential equations

The numerical approximation of solutions to linear functional differential equations are considered using the so called Legendre tau method. The functional differential equation is first reformulated as a partial differential equation with a nonlocal boundary condition involving time differentiation. The approximate solution is then represented as a truncated Legendre series with time varying coefficients which satisfy a certain system of ordinary differential equations. The method is very easy to code and yields very accurate approximations. Convergence is established, various numerical examples are presented, and comparison between the latter and cubic spline approximations is made.

Ito, K.↗

Legendre-tau approximations for functional differential equations

The numerical approximation of solutions to linear retarded functional differential equations are considered using the so-called Legendre-tau method. The functional differential equation is first reformulated as a partial differential equation with a nonlocal boundary condition involving time-differentiation. The approximate solution is then represented as a truncated Legendre series with time-varying coefficients which satisfy a certain system of ordinary differential equations. The method is very easy to code and yields very accurate approximations. Convergence is established, various numerical examples are presented, and comparison between the latter and cubic spline approximation is made.

Ito, K.↗

Solar activity and oscillation frequency splittings

Solar p-mode frequency splittings, parameterized by the coefficients through order N = 12 of a Legendre polynomial expansion of the mode frequencies as a function of m/L, were obtained from an analysis of helioseismology data taken at Big Bear Solar Observatory during the 4 years 1986 and 1988-1990 (approximately solar minimum to maximum). Inversion of the even-index splitting coefficients confirms that there is a significant contribution to the frequency splittings originating near the solar poles. The strength of the polar contribution is anti correlated with the overall level or solar activity in the active latitudes, suggesting a relation to polar faculae. From an analysis of the odd-index splitting coefficients we infer an uppor limit to changes in the solar equatorial near-surface rotatinal velocity of less than 1.9 m/s (3 sigma limit) between solar minimum and maximum.

Woodard, M. F.↗

Extensional flow convecting a reactant undergoing a first order homogeneous reaction and diffusional mass transfer from a sphere at low to intermediate Peclet and Damkohler numbers

Forced convective diffusion-reaction is considered for viscous axisymmetric extensional convecting velocity in the neighborhood of a sphere. For Peclet numbers in the range 0.1 less than or equal to Pe less than or equal to 500 and for Damkohler numbers increasing with increasing Pe but in the overall range 0.02 less than or equal to Da less than or equal to 10, average and local Sherwood numbers have been computed. By introducing the eigenfunction expansion c(r, Theta) = Sum of c(n)(r)P(n)(cos Theta) into the forced convective diffusion equation for the concentration of a chemical species undergoing a first order homogeneous reaction and by using properties of the Legendre functions Pn(cos Theta), the variable coefficient PDE can be reduced to a system of N + 1 second order ODEs for the radial functions c(sub n)(r), n = 0, 1, 2,..., N. The adaptive grid algorithm of Pereyra and Lentini can be used to solve the corresponding 2(N + 1) first order differential equations as a two-point boundary value problem on 1 less than or equal to r less than or equal to r(sub infinity). Convergence of the expansion for a specific value of N can thus be established and provides 'spectral' behavior as well as the full concentration field c(r, Theta).

Shah, N. Y.↗

A note of expansions of functions of velocity in the two-body problem

Fourier expansions based on both the true anomaly and the mean anomaly are obtained for the functions of velocity in the two-body problem; the series of coefficients is written from classical formulae involving associated Legendre polynomials, Gegenbauer polynomials, or Bessel functions. The Fourier expansions are compared with the expansions in powers of eccentricity developed by Broucke (1974) through use of computerized Poisson series manipulation.

Kinoshita, H.↗