Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spherical harmonic method”

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 181 records · Page 10

Seismic Constraints on Interior Solar Convection

We constrain the velocity spectral distribution of global-scale solar convective cells at depth using techniques of local helioseismology. We calibrate the sensitivity of helioseismic waves to large-scale convective cells in the interior by analyzing simulations of waves propagating through a velocity snapshot of global solar convection via methods of time-distance helioseismology. Applying identical analysis techniques to observations of the Sun, we are able to bound from above the magnitudes of solar convective cells as a function of spatial convective scale. We find that convection at a depth of r/R(solar) = 0.95 with spatial extent l < 30, where l is the spherical harmonic degree, comprise weak flow systems, on the order of 15 m/s or less. Convective features deeper than r/R(solar) = 0.95 are more difficult to image due to the rapidly decreasing sensitivity of helioseismic waves.

Hanasoge, Shravan M.↗

The Impact of New Accelerometer Transplant Data (ACH) on GRACE Follow-On Along-Orbit Inter-Satellite Laser Ranging Observations and Monthly Time-Variable Gravity and Mascon Solutions

GRACE-D accelerometer data show significant bias jumps since one month after the launch of the GRACE Follow-On (GRACE-FO) satellites in May 2018, making them inapplicable for correcting GRACE-FO products for non-gravitational accelerations. The GRACE-FO Science Data System (SDS) compensated this issue by transplanting the GRACE-C accelerometer data toward that of GRACE-D. Recently, GRACE-FO SDS implemented an updated transplant method, used in the latest release of GRACE-FO data. Here, we examine the impact of updated accelerometer transplant data (ACH) on GRACE-FO measurements at all levels: (a) Level-1B inter-satellite laser ranging residuals measured along satellite orbit, (b) Level-2 time-variable gravity solutions from all SDS centers (JPL, CSR, and GFZ), and (c) Level-3 mascon solutions. We show that inter-satellite laser ranging residuals are modified at low frequencies below 1 mHz, affecting the along-orbit analysis of large-scale time-variable gravity signals. When mapped into monthly Level-2 spherical harmonic coefficients of geopotential, the low-frequency change in inter-satellite ranging residuals leads to substantial improvement of coefficients associated with resonant orders (i.e., 15, 30, 45, etc.) and C 30 . We also present an improved SLR-derived C 30 which significantly improves the agreement with updated GRACE-FO C 30 at seasonal and interannual timescales. Moreover, we demonstrate the noise reduction in mass change estimates from new GRACE-FO Level-2 data over oceans, Greenland, and Antarctica for all SDS solutions. GRACE-FO mascon solutions show a moderate change in the updated release. Our comprehensive analyses demonstrate high-quality estimates of non-gravitational accelerations by the updated transplant method, resulting in more accurate GRACE-FO time-variable gravity and mass change observations.

Laser ranging interferometer↗

Incompressible fluid ellipsoids in halos. I - The second-harmonic oscillations of the Maclaurin spheroids

The structure and stability of Maclaurin spheroids embedded in rigid uniform-density oblate spheroidal halos are determined by the tensor virial-equation method. These spheroid-halo systems can be thought of as crude fluid analogs of disk galaxies with halos. The halos are assumed to have the same center, the same axis of symmetry, and the same equatorial radius as the Maclaurin spheroids. Only halos with lower eccentricity than the Maclaurin spheroids are considered. The dynamic instability of the toroidal (barlike) modes is suppressed when m, the ratio of the halo mass to Maclaurin spheroid mass, is greater than 3 pi/8 for spherical halos and when m is greater than 1/2 for halos congruent to the Maclaurin spheroids. Intermediate halo flattenings yield intermediate critical m-values. On the other hand, a neutral point of the toroidal modes in the rotating and inertial frames occurs for all m and for all allowed halo flattenings. Growth rates for secular instability beyond the neutral point are calculated, and the eigenfrequencies of all second-harmonic modes are given for select cases. The Ostriker-Peebles (1973) conjecture concerning the stability of disk galaxies against barlike perturbations appears to be incorrect.

Durisen, R. H.↗

Imaging and Segmenting Grains and Subgrains Using Backscattered Electron Techniques

We present two new methods of processing data from backscattered electron signals in a scanning electron microscope to image grains and subgrains. The first combines data from multiple backscattered electron images acquired at different specimen geometries to (1) better reveal grain boundaries in recrystallized microstructures and (2) distinguish between recrystallized and unrecrystallized regions in partially recrystallized microstructures. The second utilizes spherical harmonic transform indexing of electron backscatter diffraction patterns to produce high angular resolution orientation data that enable the characterization of subgrains. Subgrains are produced during high-temperature plastic deformation and have boundary misorientation angles ranging from a few degrees down to a few hundredths of a degree. Here, we also present an algorithm to automatically segment grains from combined backscattered electron image data or grains and subgrains from high angular resolution electron backscatter diffraction data. Together, these new techniques enable rapid measurements of individual grains and subgrains from large populations.

36 MATERIALS SCIENCE↗

Geomagnetism of earth's core

Instrumentation, analytical methods, and research goals for understanding the behavior and source of geophysical magnetism are reviewed. Magsat, launched in 1979, collected global magnetometer data and identified the main terrestrial magnetic fields. The data has been treated by representing the curl-free field in terms of a scalar potential which is decomposed into a truncated series of spherical harmonics. Solutions to the Laplace equation then extend the field upward or downward from the measurement level through intervening spaces with no source. Further research is necessary on the interaction between harmonics of various spatial scales. Attempts are also being made to analytically model the main field and its secular variation at the core-mantle boundary. Work is also being done on characterizing the core structure, composition, thermodynamics, energetics, and formation, as well as designing a new Magsat or a tethered satellite to be flown on the Shuttle.

Benton, E. R.↗

Narrow Scale Flow and a Weak Field by the Top of Earth's Core: Evidence from Orsted, Magsat and Secular Variation

As Earth's main magnetic field weakens, our magnetic shield against the onslaught of the solar wind thins. And the field strength needed to fend off battering by solar coronal mass ejections is decreasing, just when the delicate complexity of modem, vulnerable, electro-technological systems is increasing at an unprecedented rate. Recently, a working group of distinguished scientist from across the nation has asked NASA's Solid Earth and Natural Hazards program a key question: What are the dynamics of Earth s magnetic field and its interactions with the Earth system? Paleomagnetic studies of crustal rocks magnetized in the geologic past reveal that polarity reversals have occurred many times during Earth s history. Networked super-computer simulations of core field and flow, including effects of gravitational, pressure, rotational Coriolis, magnetic and viscous forces, suggest how this might happen in detail. And space-based measurements of the real, time-varying magnetic field help constrain estimates of the speed and direction of fluid iron flowing near the top of the core and enable tests of some hypotheses about such flow. Now scientists at NASA s Goddard Space Flight Center have developed and applied methods to test the hypotheses of narrow scale flow and of a dynamically weak magnetic field near the top of Earth s core. Using two completely different methods, C. V. Voorhies has shown these hypotheses lead to specific theoretical forms for the "spectrum" of Earth s main magnetic field and the spectrum of its rate of change. Much as solar physicists use a prism to separate sunlight into its spectrum, from long wavelength red to short wavelength blue light, geophysicists use a digital prism, spherical harmonic analysis, to separate the measured geomagnetic field into its spectrum, from long to short wavelength fields. They do this for the rate of change of the field as well.

Voorhies, Coerte V.↗

Ion sheath effects on RF plasma probes - Experimental results in laboratory plasmas

Laboratory work dealing with the frequency characteristic of the plasma impedance of spherical and cylindrical electrode systems is reported. The influence of the ion sheath on various features of the impedance characteristic is emphasized. Those features are the series and parallel resonance as well as additional resonances due to the excitation of electroacoustic and cyclotron harmonic waves. The dependence of the series and parallel resonances on dc biasing leads to a method of determining the ion sheath capacity for a cylindrical electrode system. The obtained values agree fairly well with those obtained from a theoretical model for the density and potential distribution in the sheath of a cylindrical sensor aligned with a supersonic plasma flow. The amplitude of resonances due to excitation of longitudinal plasma waves (electroacoustic and cyclotron harmonic) is reduced or even vanishes for sufficiently negative dc bias. Positive bias first leads to an increased amplitude up to a certain dc bias value above which, however, the amplitude decreases again due to electron absorption at the sensor surface.

Kist, R.↗

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↗

A T Matrix Method Based upon Scalar Basis Functions

A surface integral formulation is developed for the T matrix of a homogenous and isotropic particle of arbitrary shape, which employs scalar basis functions represented by the translation matrix elements of the vector spherical wave functions. The formulation begins with the volume integral equation for scattering by the particle, which is transformed so that the vector and dyadic components in the equation are replaced with associated dipole and multipole level scalar harmonic wave functions. The approach leads to a volume integral formulation for the T matrix, which can be extended, by use of Green's identities, to the surface integral formulation. The result is shown to be equivalent to the traditional surface integral formulas based on the VSWF basis.

wave functions↗

A method of defining topographic datums of planetary bodies

Since no known planetary bodies in the solar system other than the earth have large bodies of water, a topographical datum other than a sea-level reference must be used as a zero-elevation reference surface. The present paper discusses the definition of the topographic datums of Mars and the moon in terms of a gravity-level surface. Planetary gravitational field potentials were represented by a spherical harmonic expansion in terms of gravity coefficients measured by planetary orbiters, and the topographical datum was taken as the sum of an arbitrarily selected radius of the mean sphere and the radial deviation from the mean sphere. The datum defined for Mars on the basis of Mariner 9 data can be approximated as a triaxial ellipsoid with semimajor axes of 3394.6 and 3393.3 km and semiminor axis of 3376.3 km, based on a mean radius of 3382.9 km. For the moon, Lunar Orbiter IV tracking and ranging data give a datum approximated by a triaxial figure with semimajor axes 1738.30 and 1738.18 km and semiminor axis 1737.65 km, based on a mean radius of 1738 km. A topographic datum of Venus is also planned based on Pioneer-Venus gravity data.

Wu, S. S. C.↗

On the prediction of velocity fields from redshift space galaxy samples

We present a new method for recovering the underlying velocity field from an observed distribution of galaxies in redshift space. The method is based on a kinematic Zel'dovich relation between the velocity and density fields in redshift space. This relation is expressed in a differential equation slightly modified from the usual Poisson equation and which depends nontrivially on Beta identically equal to Omega(exp 0.6)/b. The linear equation can be readily solved by standard techniques of separation of variables by means of spherical harmonics. One can also include a term describing the 'rocket effect' discussed by Kaiser (1987). From this redshift space information alone, one can generate a prediction of the peculiar velocity field for each harmonic (l, m) as a function of distance. We note that for the quadrupole and higher order moments, the equation is a boundary value problem with solutions dependent on both the interior and exterior mass distribution. However, for a shell at distance r, the dipole, as well as the monopole, of the velocity field in the Local Group frame is fully determined by the interior mass distribution. This implies that the shear of the measured velocity field, when fitted to a dipole distortion, should be aligned and consistent with the gravity field inferred from the well determined local galaxy distribution. As a preliminary application we compute the velocity dipole of distant shells as predicted from the 1.2 Jy IRAS survey compared to the measured velocity dipole on shells, as inferred from a recent POTENT analysis. The coherence between the two fields is good, yielding a best estimate of Beta = 0.6 +/- 0.2.

Nusser, Adi↗

A minimum assumption approach to MEG sensor array design

Objective. Our objective is to formulate the problem of the magnetoencephalographic (MEG) sensor array design as a well-posed engineering problem of accurately measuring the neuronal magnetic fields. This is in contrast to the traditional approach that formulates the sensor array design problem in terms of neurobiological interpretability the sensor array measurements. Approach. We use the vector spherical harmonics (VSH) formalism to define a figure-of-merit for an MEG sensor array. We start with an observation that, under certain reasonable assumptions, any array of m perfectly noiseless sensors will attain exactly the same performance, regardless of the sensors' locations and orientations (with the exception of a negligible set of singularly bad sensor configurations). We proceed to the conclusion that under the aforementioned assumptions, the only difference between different array configurations is the effect of (sensor) noise on their performance. We then propose a figure-of-merit that quantifies, with a single number, how much the sensor array in question amplifies the sensor noise. Main results. We derive a formula for intuitively meaningful, yet mathematically rigorous figure-of-merit that summarizes how desirable a particular sensor array design is. We demonstrate that this figure-of-merit is well-behaved enough to be used as a cost function for a general-purpose nonlinear optimization methods such as simulated annealing. We also show that sensor array configurations obtained by such optimizations exhibit properties that are typically expected of 'high-quality' MEG sensor arrays, e.g. high channel information capacity. Significance. Our work paves the way toward designing better MEG sensor arrays by isolating the engineering problem of measuring the neuromagnetic fields out of the bigger problem of studying brain function through neuromagnetic measurements.

60 APPLIED LIFE SCIENCES↗

Finite boundary effects on the spherical Rayleigh–Taylor instability between viscous fluids

For the Rayleigh–Taylor unstable arrangement of a viscous fluid sphere embedded in a finite viscous fluid spherical shell with a rigid boundary and a radially directed acceleration, a dispersion relation is developed from a linear stability analysis using the method of normal modes. [Formula: see text] is the radially directed acceleration at the interface. ρ i denotes the density, μ i is the viscosity, and R i is the radius, where i = 1 is the inner sphere and i = 2 is the outer sphere. The dispersion relation is a function of the following dimensionless variables: viscosity ratio [Formula: see text], density ratio [Formula: see text], spherical harmonic mode n, [Formula: see text], [Formula: see text], and the dimensionless growth rate [Formula: see text], where σ is the exponential growth rate. We show that the boundedness provided by the outer spherical shell has a strong influence on the instability behavior, which is reflected not only in the modulation of the growth rate but also in the selection of the most unstable modes that are physically possible. This outer boundary effect is quantified by the relative magnitude of the radius ratio H. We find that when H is close to unity, lower order harmonics are excluded from becoming the most unstable within a vast region of the parameter space. In other words, the effect of H has precedence over the other controlling parameters d, B, and a wide range of s in establishing what the lowest most unstable mode can be. When H ~ 1, low order harmonics can become the most unstable only for s >> 1. However, in the limit when s → ∞, we show that the most unstable mode is n = 1 and derive the dispersion relation in this limit. The exclusion of most unstable low order harmonics caused by a finite outer boundary is not realized when the outer boundary extends beyond a certain threshold length-scale in which case all modes are equally possible depending on the value of B.

74 ATOMIC AND MOLECULAR PHYSICS↗

Dirac oscillator: An alternative basis for nuclear structure calculations

The isotropic harmonic oscillator supplemented by a strong spin-orbit interaction has been the cornerstone of nuclear structure since its inception more than seven decades ago. In this paper we introduce—or rather re-introduce—the “Dirac oscillator,” a fully relativistic basis that has all the desired attributes of the ordinary harmonic oscillator while naturally incorporating a strong spin-orbit coupling. To assess the power and flexibility of the Dirac oscillator basis in the solution of nuclear structure problems within the framework of covariant density-functional theory. Here, self-consistent calculations of binding energies and ground-state densities for a selected set of doubly magic nuclei are performed using the Dirac oscillator basis and are then compared against results obtained with the often-used Runge-Kutta method. Results obtained using the Dirac oscillator basis reproduce with high accuracy those derived using the Runge-Kutta method and suggest a clear path for a generalization to systems with axial symmetry. Although the harmonic oscillator with spin-orbit corrections has been the staple of the nuclear shell model since the beginning, the Dirac oscillator is practically unknown among the nuclear physics community. In this paper we illustrate the power and flexibility of the Dirac oscillator and suggest extensions to the study of systems without spherical symmetry, as required in constrained calculations of nuclear excitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A stochastic-dynamic model for global atmospheric mass field statistics

A model that yields the spatial correlation structure of atmospheric mass field forecast errors was developed. The model is governed by the potential vorticity equation forced by random noise. Expansion in spherical harmonics and correlation function was computed analytically using the expansion coefficients. The finite difference equivalent was solved using a fast Poisson solver and the correlation function was computed using stratified sampling of the individual realization of F(omega) and hence of phi(omega). A higher order equation for gamma was derived and solved directly in finite differences by two successive applications of the fast Poisson solver. The methods were compared for accuracy and efficiency and the third method was chosen as clearly superior. The results agree well with the latitude dependence of observed atmospheric correlation data. The value of the parameter c sub o which gives the best fit to the data is close to the value expected from dynamical considerations.

Ghil, M.↗

Mapping experiment with space station

Mapping of the Earth from space stations can be approached in two areas. One is to collect gravity data for defining topographic datum using Earth's gravity field in terms of spherical harmonics. The other is to search and explore techniques of mapping topography using either optical or radar images with or without reference to ground central points. Without ground control points, an integrated camera system can be designed. With ground control points, the position of the space station (camera station) can be precisely determined at any instant. Therefore, terrestrial topography can be precisely mapped either by conventional photogrammetric methods or by current digital technology of image correlation. For the mapping experiment, it is proposed to establish four ground points either in North America or Africa (including the Sahara desert). If this experiment should be successfully accomplished, it may also be applied to the defense charting systems.

Wu, S. S. C.↗

Reflection of solar radiation by a cylindrical cloud

Potential applications of an analytic method for computing the solar radiation reflected by a cylindrical cloud are discussed, including studies of radiative transfer within finite clouds and evaluations of these effects on other clouds and on remote sensing problems involving finite clouds. The pattern of reflected sunlight from a cylindrical cloud as seen at a large distance has been considered and described by the bidirectional function method for finite cloud analysis, as previously studied theoretically for plane-parallel atmospheres by McKee and Cox (1974); Schmetz and Raschke (1981); and Stuhlmann et al. (1985). However, the lack of three-dimensional radiative transfer solutions for anisotropic scattering media have hampered theoretical investigations of bidirectional functions for finite clouds. The present approach permits expression of the directional variation of the radiation field as a spherical harmonic series to any desired degree and order.

Smith, G. L.↗

Spherical harmonic analysis of steady photospheric flows. II

The use of the spherical harmonic functions to analyze the nearly steady flows in the solar photosphere is extended to situations in which B(0) the latitude at disk center, is nonzero and spurious velocities are present. The procedures for extracting the rotation profile and meridional circulation are altered to account for the seasonal tilt of the sun's rotation axis toward and away from the observer. A more robust and accurate method for separating the limb shift and meridional circulation signals is described. The analysis procedures include the ability to mask out areas containing spurious velocities (velocity-like signals that do not represent true flow velocities in the photosphere). The procedures are shown to work well in extracting the various flow components from realistic artificial data with a broad, continuous spectrum for the supergranulation. The presence of this supergranulation signal introduces errors of a few m/s in the measurements of the rotation profile, meridional circulation, and limb shift from a single Doppler image.

Hathaway, David H.↗