Search NASA⌕ Search

SEARCH · Search NASA

Results for “Functional expansion”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

A hybrid perturbation-Galerkin technique that combines multiple expansions

A two-step hybrid perturbation-Galerkin method for the solution of a variety of differential equations type problems is found to give better results when multiple perturbation expansions are employed. The method assumes that there is parameter in the problem formulation and that a perturbation method can be used to construct one or more expansions in this perturbation coefficient functions multiplied by computed amplitudes. In step one, regular and/or singular perturbation methods are used to determine the perturbation coefficient functions. The results of step one are in the form of one or more expansions each expressed as a sum of perturbation coefficient functions multiplied by a priori known gauge functions. In step two the classical Bubnov-Galerkin method uses the perturbation coefficient functions computed in step one to determine a set of amplitudes which replace and improve upon the gauge functions. The hybrid method has the potential of overcoming some of the drawbacks of the perturbation and Galerkin methods as applied separately, while combining some of their better features. The proposed method is applied, with two perturbation expansions in each case, to a variety of model ordinary differential equations problems including: a family of linear two-boundary-value problems, a nonlinear two-point boundary-value problem, a quantum mechanical eigenvalue problem and a nonlinear free oscillation problem. The results obtained from the hybrid methods are compared with approximate solutions obtained by other methods, and the applicability of the hybrid method to broader problem areas is discussed.

Geer, James F.↗
NASA NTRS1990

Solar wind proton velocity distributions - Comparison of the bi-Maxwellian based 16-moment expansion with observations

The purpose of this paper is to study the possible types of velocity distributions that can be obtained from the bi-Maxwellian based 16-moment expansion of the distribution function, assuming macroscopic parameter values characteristic of the range of solar wind conditions. While previous studies also took heat flow into account, the theoretical expansions for f and the definitions of the physical moments adopted in these studies were different from those used in this paper. The choice of the 16-moment expansion and corresponding moment definitions was motivated by the fact that this is the correct generalization of the widely-used Maxwellian-based 13-moment expansion to the case where the zeroth-order distribution is a bi-Maxwellian. It is found that most of the features characteristic of solar wind proton distributions can be reproduced with the 16-moment distribution, including the appearance of secondary peaks. It is also shown how each of the physically significant velocity moments affects the shape of the distribution function.

Demars, H. G.↗
NASA NTRS1990

Adaptive pattern recognition by mini-max neural networks as a part of an intelligent processor

In this decade and progressing into 21st Century, NASA will have missions including Space Station and the Earth related Planet Sciences. To support these missions, a high degree of sophistication in machine automation and an increasing amount of data processing throughput rate are necessary. Meeting these challenges requires intelligent machines, designed to support the necessary automations in a remote space and hazardous environment. There are two approaches to designing these intelligent machines. One of these is the knowledge-based expert system approach, namely AI. The other is a non-rule approach based on parallel and distributed computing for adaptive fault-tolerances, namely Neural or Natural Intelligence (NI). The union of AI and NI is the solution to the problem stated above. The NI segment of this unit extracts features automatically by applying Cauchy simulated annealing to a mini-max cost energy function. The feature discovered by NI can then be passed to the AI system for future processing, and vice versa. This passing increases reliability, for AI can follow the NI formulated algorithm exactly, and can provide the context knowledge base as the constraints of neurocomputing. The mini-max cost function that solves the unknown feature can furthermore give us a top-down architectural design of neural networks by means of Taylor series expansion of the cost function. A typical mini-max cost function consists of the sample variance of each class in the numerator, and separation of the center of each class in the denominator. Thus, when the total cost energy is minimized, the conflicting goals of intraclass clustering and interclass segregation are achieved simultaneously.

Szu, Harold H.↗
NASA NTRS1995

An efficient algorithm for choosing scattering directions in Monte Carlo work with arbitrary phase functions

This paper describes an efficient Monte Carlo algorithm for choosing a new direction of a photon after a scattering interaction. The algorithm chooses a scattering angle by linear interpolation in a table of the inverse cumulative scattering probability. A Legendre expansion of the phase function makes it easy to apply Clenshaw's algorithm to build the interpolation table. The points in the table are close enough together that linear interpolation is accurate. With a table of 100,000 entries, we can keep the absolute and relative errors in matching the probability distribution below 10(exp -5).

Barkstrom, Bruce R.↗
NASA NTRS2003

Approximating the Helium Wavefunction in Positronium-Helium Scattering

In the Kohn variational treatment of the positronium- hydrogen scattering problem the scattering wave function is approximated by an expansion in some appropriate basis set, but the target and projectile wave functions are known exactly. In the positronium-helium case, however, a difficulty immediately arises in that the wave function of the helium target atom is not known exactly, and there are several ways to deal with the associated eigenvalue in formulating the variational scattering equations to be solved. In this work we will use the Kohn variational principle in the static exchange approximation to d e t e e the zero-energy scattering length for the Ps-He system, using a suite of approximate target functions. The results we obtain will be compared with each other and with corresponding values found by other approximation techniques.

DiRienzi, Joseph↗
NASA NTRS1994

Spherical disharmonics in the Earth sciences and the spatial solution: Ridges, hotspots, slabs, geochemistry and tomography correlations

There is increasing use of statistical correlations between geophysical fields and between geochemical and geophysical fields in attempts to understand how the Earth works. Typically, such correlations have been based on spherical harmonic expansions. The expression of functions on the sphere as spherical harmonic series has many pitfalls, especially if the data are nonuniformly and/or sparsely sampled. Many of the difficulties involved in the use of spherical harmonic expansion techniques can be avoided through the use of spatial domain correlations, but this introduces other complications, such as the choice of a sampling lattice. Additionally, many geophysical and geochemical fields fail to satisfy the assumptions of standard statistical significance tests. This is especially problematic when the data values to be correlated with a geophysical field were collected at sample locations which themselves correlate with that field. This paper examines many correlations which have been claimed in the past between geochemistry and mantle tomography and between hotspot, ridge, and slab locations and tomography using both spherical harmonic coefficient correlations and spatial domain correlations. No conclusively significant correlations are found between isotopic geochemistry and mantle tomography. The Crough and Jurdy (short) hotspot location list shows statistically significant correlation with lowermost mantle tomography for degree 2 of the spherical harmonic expansion, but there are no statistically significant correlations in the spatial case. The Vogt (long) hotspot location list does not correlate with tomography anywhere in the mantle using either technique. Both hotspot lists show a strong correlation between hotspot locations and geoid highs when spatially correlated, but no correlations are revealed by spherical harmonic techniques. Ridge locations do not show any statistically significant correlations with tomography, slab locations, or the geoid; the strongest correlation is with lowermost mantle tomography, which is probably spurious. The most striking correlations are between mantle tomography and post-Pangean subducted slabs. The integrated locations of slabs correlate strongly with fast areas near the transition zone and the core-mantle boundary and with slow regions from 1022-1248 km depth. This seems to be consistent with the 'avalanching' downwellings which have been indicated by models of the mantle which include an endothermic phase transition at the 670-km discontinuity, although this is not a unique interpretation. Taken as a whole, these results suggest that slabs and associated cold downwellings are the dominant feature of mantle convection. Hotspot locations are no better correlated with lower mantle tomography than are ridge locations.

Ray, Terrill W.↗
NASA NTRS2011

Underexpanded Supersonic Plume Surface Interactions: Applications for Spacecraft Landings on Planetary Bodies

Numerical and experimental investigations of both far-field and near-field supersonic steady jet interactions with a flat surface at various atmospheric pressures are presented in this paper. These studies were done in assessing the landing hazards of both the NASA Mars Science Laboratory and Phoenix Mars spacecrafts. Temporal and spatial ground pressure measurements in conjunction with numerical solutions at altitudes of approx.35 nozzle exit diameters and jet expansion ratios (e) between 0.02 and 100 are used. Data from steady nitrogen jets are compared to both pulsed jets and rocket exhaust plumes at Mach approx.5. Due to engine cycling, overpressures and the plate shock dynamics are different between pulsed and steady supersonic impinging jets. In contrast to highly over-expanded (e <1) and underexpanded exhaust plumes, results show that there is a relative ground pressure load maximum for moderately underexpanded (e approx.2-5) jets which demonstrate a long collimated plume shock structure. For plumes with e much >5 (lunar atmospheric regime), the ground pressure is minimal due to the development of a highly expansive shock structure. We show this is dependent on the stability of the plate shock, the length of the supersonic core and plume decay due to shear layer instability which are all a function of the jet expansion ratio. Asymmetry and large gradients in the spatial ground pressure profile and large transient overpressures are predominantly linked to the dynamics of the plate shock. More importantly, this study shows that thruster plumes exhausting into martian environments possess the largest surface pressure loads and can occur at high spacecraft altitudes in contrast to the jet interactions at terrestrial and lunar atmospheres. Theoretical and analytical results also show that subscale supersonic cold gas jets adequately simulate the flow field and loads due to rocket plume impingement provided important scaling parameters are in agreement. These studies indicate the critical importance of testing and modeling plume-surface interactions for descent and ascent of spacecraft and launch vehicles.

Mehta, M.↗
NASA NTRS1993

Efficient and accurate computation of the incomplete Airy functions

The incomplete Airy integrals serve as canonical functions for the uniform ray optical solutions to several high-frequency scattering and diffraction problems that involve a class of integrals characterized by two stationary points that are arbitrarily close to one another or to an integration endpoint. Integrals with such analytical properties describe transition region phenomena associated with composite shadow boundaries. An efficient and accurate method for computing the incomplete Airy functions would make the solutions to such problems useful for engineering purposes. In this paper a convergent series solution for the incomplete Airy functions is derived. Asymptotic expansions involving several terms are also developed and serve as large argument approximations. The combination of the series solution with the asymptotic formulae provides for an efficient and accurate computation of the incomplete Airy functions. Validation of accuracy is accomplished using direct numerical integration data.

Constantinides, E. D.↗
NASA NTRS1992

Numerical evaluation of the incomplete airy functions and their application to high frequency scattering and diffraction

The incomplete Airy integrals serve as canonical functions for the uniform ray optical solutions to several high frequency scattering and diffraction problems that involve a class of integrals characterized by two stationary points that are arbitrarily close to one another or to an integration endpoint. Integrals of such analytical properties describe transition region phenomena associated with composite shadow boundaries. An efficient and accurate method for computing the incomplete Airy functions would make the solutions to such problems useful for engineering purposes. Here, a convergent series solution form for the incomplete Airy functions is derived. Asymptotic expansions involving several terms were also developed and serve as large argument approximations. The combination of the series solution form with the asymptotic formulae provides for an efficient and accurate computation of the incomplete Airy functions. Validation of accuracy is accomplished using direct numerical integration data.

Constantinides, E. D.↗
NASA NTRS1974

Electromagnetic and scalar diffraction by a right-angled wedge with a uniform surface impedance

The diffraction of an electromagnetic wave by a perfectly-conducting right-angled wedge with one surface covered by a dielectric slab or absorber is considered. The effect of the coated surface is approximated by a uniform surface impedance. The solution of the normally incident electromagnetic problem is facilitated by introducing two scalar fields which satisfy a mixed boundary condition on one surface of the wedge and a Neumann of Dirichlet boundary condition on the other. A functional transformation is employed to simplify the boundary conditions so that eigenfunction expansions can be obtained for the resulting Green's functions. The eigenfunction expansions are transformed into the integral representations which then are evaluated asymptotically by the modified Pauli-Clemmow method of steepest descent. A far zone approximation is made to obtain the scattered field from which the diffraction coefficient is found for scalar plane, cylindrical or sperical wave incident on the edge. With the introduction of a ray-fixed coordinate system, the dyadic diffraction coefficient for plane or cylindrical EM waves normally indicent on the edge is reduced to the sum of two dyads which can be written alternatively as a 2 X 2 diagonal matrix.

Hwang, Y. M.↗
NASA NTRS1977

The decompression of cold neutron star matter

The ejection of cold neutron-star matter is examined, and an attempt is made to determine whether the final composition of this matter may be similar to that normally associated with the hot high-neutron-flux r-process. A semiempirical liquid-drop model is used for the nucleus, and the equilibrium composition of the matter is determined by assuming it to be in its absolute ground state at a given density. Physical mechanisms operating during the expansion are analyzed, and the composition of the ejected matter is found as a function of its density during expansion. The results indicate that it is virtually impossible for deuterium to form, that neutrons can be captured only after beta decay increases the atomic numbers of nuclei, and that no free neutrons can escape. It is concluded that neutron-star ejecta can produce heavy neutron-rich nuclei and may produce somewhat heavier nuclei than a standard r-process.

Lattimer, J. M.↗
NASA NTRS2010

Reflector Surface Error Compensation in Dual-Reflector Antennas

By probing the field on a small subreflector at a minimal number of points, the main reflector surface errors can be obtained and subsequently used to design a phase-correction subreflector that can compensate for main reflector errors. The compensating phase-error profile across the subreflector can be achieved either by a surface deformation or by the use of an array of elements such as patch antennas that can cause a phase shift between the incoming and outgoing fields. The second option is of primary interest here, but the methodology can be applied to either case. The patch array is most easily implemented on a planar surface. Therefore, the example of a flat subreflector and a parabolic main reflector (a Newtonian dual reflector system) is considered in this work. The subreflector is assumed to be a reflector array covered with patch elements. The phase variation on a subreflector can be detected by a small number of receiving patch elements (probes). By probing the phase change at these few selected positions on the subreflector, the phase error over the entire surface can be recovered and used to change the phase of all the patch elements covering the subreflector plane to compensate for main reflector errors. This is accomplished by using a version of sampling theorem on the circular aperture. The sampling is performed on the phase-error function on the circular aperture of the main reflector by a method developed using Zernike polynomials. This method is based upon and extended from a theory previously proposed and applied to reflector aperture integration. This sampling method provides for an exact retrieval of the coefficients of up to certain orders in the expansion of the phase function, from values on a specifically calculated set of points in radial and azimuthal directions in the polar coordinate system, on the circular reflector aperture. The corresponding points on the subreflector are then obtained and, by probing the fields at these points, a set of phase values is determined that is then transferred back to the main reflector aperture for recovering the phase function. Once this function is recovered, the corresponding phase function on the subreflector is calculated and used to compensate for main reflector surface errors. In going back and forth between sub and main reflectors, geometrical (ray) optics is employed, which even though it ignores edge diffraction and other effects, is shown to be accurate for phase recovery.

Jamnejad, Vahraz↗
NASA NTRS1989

Numerically efficient algorithm for model development of high-order systems

A technique for estimating transfer functions in partial fraction expansion form from frequency response data for a high-order system is presented. The problem formulation avoids many of the numerical difficulties associated with high-order polynomials and has the advantage of having the option to fix the camping and frequency of a mode, if known, during the estimation process. The resulting transfer function(s) may be converted to Jordan-Form time domain equations directly. During the implementation of this technique, a frequency and amplitude normalizing window was developed that maximized the efficiency of the optimization algorithm. The combination of estimating the transfer function in factored form, the ability to fix preciously determined parameters and the effectiveness of the normalizing window led to a progressive approach to synthesizing transfer functions from frequency response data for high-order systems.

Parada, L. O.↗
NASA NTRS1989

A test of a modified algorithm for computing spherical harmonic coefficients using an FFT

The Dilts (1985) algorithm for computing the spherical harmonic expansion coefficients for a function on a sphere, on the basis of a two-dimensional FFT, is presently modified, tested, and found to eliminate problems of overflow and large storage requirements associated with the encounter of harmonic degree values greater than 16. Results from timing tests show the Dilts program to be impractical, however, for the computation of spherical harmonic expansion coefficients for large harmonic degree values.

Elowitz, Mark↗
NASA NTRS1994

An Efficient Spectral Method for Ordinary Differential Equations with Rational Function Coefficients

We present some relations that allow the efficient approximate inversion of linear differential operators with rational function coefficients. We employ expansions in terms of a large class of orthogonal polynomial families, including all the classical orthogonal polynomials. These families obey a simple three-term recurrence relation for differentiation, which implies that on an appropriately restricted domain the differentiation operator has a unique banded inverse. The inverse is an integration operator for the family, and it is simply the tridiagonal coefficient matrix for the recurrence. Since in these families convolution operators (i.e. matrix representations of multiplication by a function) are banded for polynomials, we are able to obtain a banded representation for linear differential operators with rational coefficients. This leads to a method of solution of initial or boundary value problems that, besides having an operation count that scales linearly with the order of truncation N, is computationally well conditioned. Among the applications considered is the use of rational maps for the resolution of sharp interior layers.

Coutsias, Evangelos A.↗
NASA NTRS1994

Thermal expansion and elastic anisotropy in single crystal Al2O3 and SiC reinforcements

In single crystal form, SiC and Al2O3 are attractive reinforcing components for high temperature composites. In this study, the axial coefficients of thermal expansion and single crystal elastic constants of SiC and Al2O3 were used to determine their coefficients of thermal expansion and Young's moduli as a function of crystallographic orientation and temperature. SiC and Al2O3 exhibit a strong variation of Young's modulus with orientation; however, their moduli and anisotropies are weak functions of temperature below 1000 C. The coefficients of thermal expansion exhibit significant temperature dependence, and that of the non-cubic Al2O3 is also a function of crystallographic orientation.

Salem, Jonathan A.↗
NASA NTRS1988

Investigating a hybrid perturbation-Galerkin technique using computer algebra

A two-step hybrid perturbation-Galerkin method is presented for the solution of a variety of differential equations type problems which involve a scalar parameter. The resulting (approximate) solution has the form of a sum where each term consists of the product of two functions. The first function is a function of the independent field variable(s) x, and the second is a function of the parameter lambda. In step one the functions of x are determined by forming a perturbation expansion in lambda. In step two the functions of lambda are determined through the use of the classical Bubnov-Gelerkin method. The resulting hybrid method has the potential of overcoming some of the drawbacks of the perturbation and Bubnov-Galerkin methods applied separately, while combining some of the good features of each. In particular, the results can be useful well beyond the radius of convergence associated with the perturbation expansion. The hybrid method is applied with the aid of computer algebra to a simple two-point boundary value problem where the radius of convergence is finite and to a quantum eigenvalue problem where the radius of convergence is zero. For both problems the hybrid method apparently converges for an infinite range of the parameter lambda. The results obtained from the hybrid method are compared with approximate solutions obtained by other methods, and the applicability of the hybrid method to broader problem areas is discussed.

Andersen, Carl M.↗
NASA NTRS1977

A recursively formulated first-order semianalytic artificial satellite theory based on the generalized method of averaging. Volume 1: The generalized method of averaging applied to the artificial satellite problem

A recursively formulated, first-order, semianalytic artificial satellite theory, based on the generalized method of averaging is presented in two volumes. Volume I comprehensively discusses the theory of the generalized method of averaging applied to the artificial satellite problem. Volume II presents the explicit development in the nonsingular equinoctial elements of the first-order average equations of motion. The recursive algorithms used to evaluate the first-order averaged equations of motion are also presented in Volume II. This semianalytic theory is, in principle, valid for a term of arbitrary degree in the expansion of the third-body disturbing function (nonresonant cases only) and for a term of arbitrary degree and order in the expansion of the nonspherical gravitational potential function.

Mcclain, W. D.↗