Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized spherical functions”

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 55 records · Page 3

A parametric heat transfer study for cryogenic ball bearings in SSME HPOTP

A numerical modeling is to examine the effects of coolant convective heat transfer coefficient and frictional heating on the local temperature characteristics of a ball element in Space Shuttle Main Engine (SSME) High Pressure Oxidizer Turbopump (HPOTP) bearing. The present modeling uses a control-volume based, finite-difference method to solve the non-dimensionalized heat conduction equation in spherical coordinate system. The dimensionless temperature is found as a function of Biot number, heat flux ratio between the two race contacts, and location in the ball. The current results show that, for a given cooling capability, the ball temperature generally increases almost linearly with the heat input from the race-contacts. This increase is always very high at one of the two contacts. An increase in heat transfer coefficient generally reduces the ball temperature and alleviates the temperature gradient, except for the regions very close to the race contacts. For a 10-fold increase of heat transfer coefficient, temperature decrease is 35 percent for the average over entire ball, and 10 percent at the inner-race contact. The corresponding change of temperature gradient displays opposing trends between the regions immediately adjacent to the contacts and the remaining portion of the ball. The average temperature gradient in the vicinity of both contacts increases approximately 70 to 100 percent. A higher temperature gradient produces excessive thermal stress locally which may be detrimental to the material integrity. This, however, is the only unfavorable issue for an increase of heat transfer coefficient.

Chyu, Mingking K.↗

Wigner functions for nonclassical states of a collection of two-level atoms

The general theory of atomic angular momentum states is used to derive the Wigner distribution function for atomic angular momentum number states, coherent states, and squeezed states. These Wigner functions W(theta,phi) are represented as a pseudo-probability distribution in spherical coordinates theta and phi on the surface of a sphere of radius the square root of j(j +1) where j is the total angular momentum.

Agarwal, G. S.↗

Review of aircraft noise propagation

The current state of knowledge about the propagation of aircraft noise was reviewed. The literature on the subject is surveyed and methods for predicting the most important and best understood propagation effects are presented. Available empirical data are examined and the data's general validity is assessed. The methods used to determine the loss of acoustic energy due to uniform spherical spreading, absorption in a homogeneous atmosphere, and absorption due to ground cover are presented. A procedure for determining ground induced absorption as a function of elevation angle between source and receiver is recommended. Other factors that affect propagation, such as refraction and scattering due to turbulence, which were found to be less important for predicting the propagation of aircraft noise, are also evaluated.

Putnam, T. W.↗

Theory of local thermal instability in spherical systems

The gasdynamical properties of local thermal instability in optically thin astrophysical plasmas as it occurs in spherical accretion and winds is investigated. In a medium characterized by both thermal and hydrostatic equilibrium, if the cooling function is not an explicit function of position and does not display isentropic thermal instability, then isobaric thermal instability by the Field criterion is present if and only if convective instability is present by the Schwarzschild criterion. In this case, thermal overstability cannot occur. Convective instability by the Schwarzschild criterion will also occur in accretion flows locally dominated by external heating or in marginally unbound, radiatively cooling outflows. A very general Lagrangian equation for the development of nonradial thermal instability in flows with spherical symmetry is derived and is solved analytically in certain regimes. The results are applied to cluster X-ray cooling flows.

Balbus, Steven A.↗

Estimating the Subsolar Magnetopause Position from Soft X-Ray Images Using a Low-Pass Image Filter

The Lunar Environment heliospheric X-ray Imager (LEXI) and Solar wind Magnetosphere Ionosphere Link Explorer (SMILE) missions will image the Earth’s dayside magnetopause and cusps in soft X-rays after their respective launches in the near future, to specify global magnetic reconnection modes for varying solar wind conditions. To support the success of these scientific missions, it is critical to develop techniques that extract the magnetopause locations from the observed soft X-ray images. In this research, we introduce a new geometric equation that calculates the subsolar magnetopause position ( R s ) from a satellite position, the look direction of the instrument, and the angle at which the X-ray emission is maximized. Two assumptions are used in this method: (1) The look direction where soft X-ray emissions are maximized lies tangent to the magnetopause, and (2) the magnetopause surface near the subsolar point is almost spherical and thus R s is nearly equal to the radius of the magnetopause curvature. We create synthetic soft X-ray images by using the Open Geospace General Circulation Model (OpenGGCM) global magnetohydrodynamic model, the galactic background, the instrument point spread function, and Poisson noise. We then apply the fast Fourier transform and Gaussian low-pass filters to the synthetic images to remove noise and obtain accurate look angles for the soft X-ray peaks. From the filtered images, we calculate R 2 and its accuracy for different LEXI locations, look directions, and solar wind densities by using the OpenGGCM subsolar magnetopause location as ground truth. Our method estimates R s with an accuracy of <0.3 R E when the solar wind density exceeds >10 cm -3 . The accuracy improves for greater solar wind densities and during southward interplanetary magnetic fields. The method captures the magnetopause motion during southward interplanetary magnetic field turnings. Consequently, the technique will enable quantitative analysis of the magnetopause motion and help reveal the dayside reconnection modes for dynamic solar wind conditions. This technique will support the LEXI and SMILE missions in achieving their scientific objectives.

Hyangpyo Kim↗

Design of efficient stiffened shells of revolution

A method to produce efficient piecewise uniform stiffened shells of revolution is presented. The approach uses a first order differential equation formulation for the shell prebuckling and buckling analyses and the necessary conditions for an optimum design are derived by a variational approach. A variety of local yielding and buckling constraints and the general buckling constraint are included in the design process. The local constraints are treated by means of an interior penalty function and the general buckling load is treated by means of an exterior penalty function. This allows the general buckling constraint to be included in the design process only when it is violated. The self-adjoint nature of the prebuckling and buckling formulations is used to reduce the computational effort. Results for four conical shells and one spherical shell are given.

Majumder, D. K.↗

Empirical wind model for the middle and lower atmosphere. Part 2: Local time variations

The HWM90 thermospheric wind model was revised in the lower thermosphere and extended into the mesosphere and lower atmosphere to provide a single analytic model for calculating zonal and meridional wind profiles representative of the climatological average for various geophysical conditions. Local time variations in the mesosphere are derived from rocket soundings, incoherent scatter radar, MF radar, and meteor radar. Low-order spherical harmonics and Fourier series are used to describe these variations as a function of latitude and day of year with cubic spline interpolation in altitude. The model represents a smoothed compromise between the original data sources. Although agreement between various data sources is generally good, some systematic differences are noted. Overall root mean square differences between measured and model tidal components are on the order of 5 to 10 m/s.

Hedin, A. E.↗

A Spherical Earth Solution for TOA Lightning Location Retrieval

With the advent of high-speed digital computers, the method of chi square minimization is a highly practical means for analyzing a wide variety of (otherwise intractable) nonlinear inversion problems in applied mathematical physics. Little thought or effort is required to apply the chi square method to obtain quick and reasonable estimates of a solution, and the method offers a means to assess retrieval errors. Because the method is simple and practical it is sometimes hastily applied to problems that can be solved by formal analytic or quasi-analytic means. Presently, Global Atmospherics Inc. (GAI) finds the minimum of a chi square function to analyze time-of-arrival (TOA) and magnetic bearing data derived from the National Lightning Detection Network (NLDN); ellipsoidal Earth geometry is assumed. An analytic solution to this problem has not yet been attained, but the consideration and solving of less general problem statements might eventually lead to a final solution. In the present study, the problem of retrieving lightning ground-strike location on a spherical Earth surface using a network of 4 or more time-of-arrival (TOA) sensors is considered. It is shown that this problem has an analytic solution and therefore does not require the use of nonlinear estimation theory (such as the chi square method mentioned above). The mathematical robustness of the analytic solution is tested using computer-generated lightning sources and simulated TOA measurement errors. A quasi-analytic extension of the spherical Earth solution for an oblate spheroidal Earth geometry is considered in a related study. The incorporation of magnetic bearing information into these analytic solutions would lead to a general and elegant analytic retrieval scheme that would most likely replace the chi square estimation theory currently employed by Global Atmospherics Inc. (GAI).

Koshak, William J,↗

Distortion Representation of Forecast Errors for Model Skill Assessment and Objective Analysis

We study a novel characterization of errors for numerical weather predictions. In its simplest form we decompose the error into a part attributable to phase errors and a remainder. The phase error is represented in the same fashion as a velocity field and will be required to vary slowly and smoothly with position. A general distortion representation allows for the displacement and a bias correction of forecast anomalies. In brief, the distortion is determined by minimizing the objective function by varying the displacement and bias correction fields. In the present project we use a global or hemispheric domain, and spherical harmonics to represent these fields. In this project we are initially focusing on the assessment application, restricted to a realistic but univariate 2-dimensional situation. Specifically we study the forecast errors of the 500 hPa geopotential height field for forecasts of the short and medium range. The forecasts are those of the Goddard Earth Observing System data assimilation system. Results presented show that the methodology works, that a large part of the total error may be explained by a distortion limited to triangular truncation at wavenumber 10, and that the remaining residual error contains mostly small spatial scales.

Hoffman, Ross N.↗

General relativity and satellite orbits - The motion of a test particle in the Schwarzschild metric

The motion of a satellite of negligible mass about a massive spherically symmetric object in a space of Schwarzschild geometry is calculated by using ordinary Newtonian methods for an inverse-square gravitational field acted upon by an inverse-cube disturbing function. The disturbing function is expressed in terms of the Keplerian elements of the orbit and substituted in the Lagrange planetary equations. The equations can be integrated to calculate the displacement in position of the satellite due to the relativistic potential. For example, the Beacon Explorer C satellite is found to be displaced by about 17.4 cm after one revolution. For Mercury the advance of the perihelion is about 85 km after one orbit, and the maximum periodic displacement is about 13 km.

Rubincam, D. P.↗

Predictions of dose from electrons in space

The objective of the project is to develop a general-purpose, user-friendly computerized database and code package, for the PC as well as larger computers, which can be used for the routine prediction of the absorbed dose from incident electrons and their secondary bremsstrahlung (and from incident protons) as functions of the thickness of aluminum shielding in space. The assumption of homogeneous aluminum shields and of isotropic incident fluxes (at least in a time-averaged sense) allows for the rather reliable conversion of doses in slabs to those in other simple bodies, such as spherical and cylindrical solids and shells. On such a basis, depth-dose data for monoenergetic incident radiation can be generated once-and-for-all from accurate transport calculations, and this database can then be used repeatedly in rapid dose predictions for arbitrary radiation spectra and for a variety of spacecraft sizes and shapes, without recourse to the very time-consuming Monte Carlo calculations. This project entails a thorough updating, extension, and refinement of our earlier SHIELDOSE package, with the goal of a more reliable, fool-proof, and general system.

Seltzer, Stephen M.↗

Retrieving the Balanced Winds on the Globe as a Generalized Inverse Problem

A generalized inverse technique is applied to retrieve two types of balanced winds that characterize the large-scale dynamics of the atmosphere: rotational winds based upon the linear balance equation, and divergent winds based upon the vorticity budget equation. Both balance equations are singular at or near the equator. The balance equations are transformed in spherical harmonic function space to an under-determined system, for which the scale-weighed least-squares solution consists of a sum of principal and singular components. The principal components represent the response to the source function for the regular eigenmodes, while the singular components are determined by the projection of an independent measurement on the singular eigenmodes. The method was tested with the NCEP/NCAR reanalysis data in which a quasi-balance condition exists. A realistic balanced wind field is retrievable when the singular components are computed based upon the reanalyzed wind data.

Lu, Huei-Iin↗

Analysis of axisymmetrical vibration of a partially liquid-filled elastic sphere by the method of Green's function

The longitudinal dynamic behavior of launch vehicles is largely determined by their huge liquid propellant masses which are spring supported by the elastic tank walls. Throughout much of the powered flight time, these masses constitute a high percentage of the entire vehicle mass and, therefore, may dominate the fundamental modes of the vehicle. In this report, a spherical container is considered. The analysis is based on a Galerkin approach, in the course of which a second-order differential equation must be solved. The solution has been obtained by the method of Green's function. This method is favorable because it displays the manner in which the analysis can be extended to partially liquid-filled general shells of revolution. The computer programs currently available for partially liquid-filled propellant tanks are based on the finite element methods and result in analytical models having as many as several hundred degrees of freedom. The method applied in this report results in a model having less than 10 degrees of freedom as can be shown by numerical evaluation. Therefore, it will be possible to analyze propellant tanks using much less computer time with comparable accuracy.

Glaser, R. F.↗

Radiative transfer in a sphere illuminated by a parallel beam - An integral equation approach

The problem of multiple scattering of nonpolarized light in a planetary body of arbitrary shape illuminated by a parallel beam is formulated using the integral equation approach. There exists a simple functional whose stationarity condition is equivalent to solving the equation of radiative transfer and whose value at the stationary point is proportional to the differential cross section. The analysis reveals a direct relation between the microscopic symmetry of the phase function for each scattering event and the macroscopic symmetry of the differential cross section for the entire planetary body, and the interconnection of these symmetry relations and the variational principle. The case of a homogeneous sphere containing isotropic scatterers is investigated in detail. It is shown that the solution can be expanded in a multipole series such that the general spherical problem is reduced to solving a set of decoupled integral equations in one dimension. Computations have been performed for a range of parameters of interest, and illustrative examples of applications to planetary problems as provided.

Shia, R.-L.↗

Radiative transfer in spherical dust shells using a generalized two-stream Eddington approximation

Application of a generalized two-stream Eddington approximation to the problem of radiative transfer in extended, spherically symmetric dust shells is presented. It is assumed that the radiation field can be characterized by the mean intensity, the flux, and a positionally dependent direction cosine specifying the division into two solid-angle ranges. The direction cosine is not specified a priori and is a function of the geometry, opacity, and emissivity in the dust shell. A multiple-grain-size multiple-temperature-distribution dust shell is postulated in which isotropic and anisotropic scattering as well as absorption and thermal reemission are allowed. A program has been developed that solves for the multiple temperature distributions by applying the constraint of radiative equilibrium to each grain size, and then calculates emergent fluxes. Results of one such calculation are presented for a model dust shell having a maximum optical depth (approximately 41) in the visible, clearly showing large optical extinction and a moderate infrared excess.

Haisch, B. M.↗

Transition regime droplet growth and evaporation - An integrodifferential variational approach

A variational principle based on the integrodifferential form of the Boltzmann equation, which allows very general forms of the collision integral and of the boundary operator, has been used to compute the mass flux to a spherical particle. For simplicity the problem has been restricted to the one-speed, constant cross-section approximation, the black sphere problem of neutron transport theory. Using the Hilbert expansion of the solution far from the sphere as a trial function leads to a rational expression in the inverse Knudsen number for the extrapolation distance, which is simply related to the number flux to the sphere. A five-term expansion gives results that are no more than 6 percent in error when compared to an accurate numerical solution. The technique can be applied to other forms of the collision model and other boundary conditions.

Cipolla, J. W., Jr.↗

Tectonics of planetary loading - A general model and results

The tectonics of planetary loading is investigated using an analytical model for determining the stresses in an arbitrarily thick spherical shell due to an idealized axisymmetric load. The model includes the flat plate and thin shell membrane approximations as end members, and makes it possible to determine the nature of the transition between them. Using this model, the stress states and the resulting tectonic patterns due to an idealized exponential load are determined as functions of five dimensionless parameters: the ratio of the lithospheric thickness to the planetary radius; the decay width of the load; the 'support parameter', which is the ratio of the buoyancy to the flexural support; the angular distance from the load center; and the normalized radial distance from the planet center.

Janes, D. M.↗

Rigorous Justification of the Localized Approximation to the Beam-Shape Coefficients in Generalized Lorenz-Mie Theory: On-Axis Beams - Part 1

Generalized Lorenz-Mie theory describes electromagnetic scattering of an arbitrary light beam by a spherical particle. The computationally most expensive feature of the theory is the evaluation of the beam-shape coefficients, which give the decomposition of the incident light beam into partial waves. The so-called localized approximation to these coefficients for a focused Gaussian beam is an analytical function whose use greatly simplifies Gaussian-beam scattering calculations. A mathematical justification and physical interpretation of the localized approximation is presented for on-axis beams.

Lock, James A.↗