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 199 records · Page 11

Theory of an experiment in an orbiting space laboratory to determine the gravitational constant.

An experiment is discussed for determining the gravitational constant with the aid of an isolated system consisting of an artificial satellite moving around an artificial planet. The experiment is to be conducted in a spherical laboratory traveling in an orbit around the earth. Difficulties due to the gravity-gradient term are considered, and the three-tunnel method proposed by Wilk (1969) is examined. The rotation of the sphere is discussed together with aspects of the reference systems used, the equations of motion of the spacecraft and of the test objects, the field from the earth's gravity gradient at the test object, higher harmonic terms in the gravity gradient force, gravitational effects of the spacecraft itself, and a computer simulation.

Vinti, J. P.↗

Vertical structure of the wind field during the Special Observing Period I of the Global Weather Experiment

The vertical structure of the global atmosphere is analyzed for selected periods of the Special Observing Period I (SOP-I) for the Global Weather Experiment (GWE). The analysis consists of projection of the stream-function and velocity potential at 200 and 850 mb on spherical harmonics and of the wind and height fields on the normal modes of a linearized form of the primitive equations for a basic state at rest. The kinematic vertical structure is discussed in terms of correlation coefficients of the 200 mb and 850 mb winds and analysis of the internal and external normal modes of the primitive equations. The reliability of the results is checked by applying the same analysis methods to data sets obtained from three different institutions: Geophysical Fluid Dynamics Laboratory (GFDL), European Center for Medium Range Weather Forecasting (ECMWF), and Goddard Laboratory for the Atmospheres (GLA). It is found that, on a global basis, vertically reversing circulations are as important as the equivalent barotropic structures. For the verticaly reversing components, the gravity and mixed Rossby-gravity modes have contributions of the same order of magnitude as those of the Rossby modes in tropical latitudes.

Paegle, J. N.↗

Solar oscillations - A method for deriving nonlinear effects

The frequencies of solar oscillations are so closely spaced that nonlinear interactions among modes are probable. The rate of interaction is proportional to an integral involving the eigenfunctions of the interacting modes, which usually are known only numerically. An approximation in which the eigenfunctions are strictly sinusoidal functions of a suitably defined radial variable, and the numerical details of stellar structure are banished to a coefficient in the integrand, are here explored. The physical assumptions are the same as in the asymptotic approximation of p- or g-modes. This method should allow a general investigation as to likely nonlinear interactions and which modes may participate in such interactions. The coupling of two p-modes by possible large-scale internal magnetic fields is introduced as an anisotropic pressure response to a displacement. Pairs of modes differing in frequency by less than the fraction magnetic/thermal pressure are strongly coupled, and energy appears to oscillate slowly between the two associated spherical harmonics. Potentially, upper limits may be derived for internal magnetic fields.

Wentzel, Donat G.↗

A method of evaluating quantitative magnetospheric field models by an angular parameter alpha

The paper introduces an angular parameter, termed alpha, which represents the angular difference between the observed, or model, field and the internal model field. The study discusses why this parameter is chosen and demonstrates its usefulness by applying it to both observations and models. In certain areas alpha is more sensitive than delta-B (the difference between the magnitude of the observed magnetic field and that of the earth's internal field calculated from a spherical harmonic expansion) in expressing magnetospheric field distortions. It is recommended to use both alpha and delta-B in comparing models with observations.

Sugiura, M.↗

Integrated and spectral energy flows of the GLAS GCM

Methods to analyze the generation, transport, and dissipation of energy to study geophysical fluid flows are discussed. Energetics analyses are pursued in several directions: (1) the longitudinal and time dependence on the energy flow to the stratosphere was examined as a function of geographical sector; (2) strong and weak energy flows were correlated by medium range forecasts; (3) the one dimensional spectral results (Fourier services around latitude circles) were extended to spherical harmonics over a global domain; (4) the validity of vertical velocities derived from mass convergence was examined for their effect on the conversion of eddy available potential energy to eddy kinetic energy.

Tennebaum, J.↗

Application of lagged averaged forecasting to medium range forecasting

The lagged average forecast (LAF) method for predicting and compensating for forecast error growth is applied to two 100 day samples of 10 day forecasts at the 500 mb altitude for winter-summer 1980-81. The LAF parameterizes the forecast error growth in order to weight the forecast with regression calculations. A 5 day LAF forecast was generated for a 10 day forecast made with a 1681 d.o.f. model which produced predictions using a truncated spherical harmonic expansion. A 100 day data set was employed, with the first 90 days serving for LAF forecasts for the model forecast. Significant improvements were obtained in the model forecast when the LAF weightings were introduced into the model variables.

Dalcher, A.↗

Mars gravity derived from long-period motion of Mariner 9

The spherical harmonic coefficients in the gravity field of Mars are determined using a two stage filtering process. In the first stage, a weighted least squares orbit determination processor is applied to Mariner 9 Doppler tracking data to estimate Kepler elements and rates. In the second stage, the Kepler element rates are used as input to another weighted least squares processor that estimates gravity coefficients using the Lagrange perturbation equations. The novelty in this method is the fitting of element rates rather than the elements themselves. This procedure decouples the results from the effects of internal spacecraft dynamics.

Ferrari, A. J.↗

Mars gravity field model from Mariner 9, Viking 1 and 2 data

Earth artificial satellite methods are extended by means of two way Doppler data, to the computation of very accurate orbiter trajectories around another planet, and to the determination of its gravity field. It is reported that in the case of Mars, all observations collected by 10 Deep Space Network stations located at three different sites during the Mariner 9 and Viking 1 and 2 missions have been processed and used to compute a full twelfth degree and order spherical harmonic model of the gravitational potential. It is concluded that the aeroid derived from the model shows very large correlations with the Martian topography, raising questions as to the deep structure of the planet which cannot be interpreted on the basis of topographic and isostatic considerations alone.

Balmino, G.↗

Local Lunar Gravity Field Analysis over the South Pole-aitken Basin from SELENE Farside Tracking Data

We present a method with which we determined the local lunar gravity field model over the South Pole-Aitken (SPA) basin on the farside of the Moon by estimating adjustments to a global lunar gravity field model using SELENE tracking data. Our adjustments are expressed in localized functions concentrated over the SPA region in a spherical cap with a radius of 45deg centered at (191.1 deg E, 53.2 deg S), and the resolution is equivalent to a 150th degree and order spherical harmonics expansion. The new solution over SPA was used in several applications of geophysical analysis. It shows an increased correlation with high-resolution lunar topography in the frequency band l = 40-70, and admittance values are slightly different and more leveled when compared to other, global gravity field models using the same data. The adjustments expressed in free-air anomalies and differences in Bouguer anomalies between the local solution and the a priori global solution correlate with topographic surface features. The Moho structure beneath the SPA basin is slightly modified in our solution, most notably at the southern rim of the Apollo basin and around the Zeeman crater

Lunar Gravity↗

Stabilized determination of geopotential coefficients by the mixed hom-BLUP approach

For the determination of geopotential coefficients, data can be used from rather different sources, e.g., satellite tracking, gravimetry, or altimetry. As each data type is particularly sensitive to certain wavelengths of the spherical harmonic coefficients it is of essential importance how they are treated in a combination solution. For example the longer wavelengths are well described by the coefficients of a model derived by satellite tracking, while other observation types such as gravity anomalies, delta g, and geoid heights, N, from altimetry contain only poor information for these long wavelengths. Therefore, the lower coefficients of the satellite model should be treated as being superior in the combination. In the combination a new method is presented which turns out to be highly suitable for this purpose due to its great flexibility combined with robustness.

Middel, B.↗

Precise solution for a finite set of spherical coefficients from equiangular gridded data

An important goal of geodesy is to determine the anomalous potential and its derivatives outside of the earth. Representing the surface anomalies by a series of spherical harmonics is useful since it is then possible to do a term by term solution of Laplace's equation and upward continuation. The problem of finding such a spherical harmonic series for anomaly values given on an equiangular surface grid is addressed. (This is a first step toward the more complicated problem of finding a function such that locally averaged values fit a grid of mean anomalies.) Three approaches to this fitting problem are discussed and compared: the discrete Fourier technique, the discrete integral technique, and a new approach. The peculiar nature of the equiangular grid, with its increasing density of (noisy) data toward the poles, causes each method to exhibit a different type of difficulty. The new method is shown to be practical as well as precise since the numerical conditioning problems which appear can be successfully handled by such well-known techniques as a (simple) Kalman filter.

Zucker, Paul A.↗

Post-glacial relaxation of a viscously stratified compressible mantle

The postglacial relaxation phenomenon of a viscously compressible nonself-gravitating spherical-shell model is investigated. Analytical solutions to an exponentially depth-dependent viscosity are developed with various types of compressible density models such as those with an exponential dependence of the radius and an algebraic-root dependence of the radius. The solutions of a viscously compressible multilayered model with constant thermodynamic properties and viscosity are developed by the propagator matrix method. The results show that inferences of deep mantle viscosity from postglacial rebound would be hampered by mantle compressibility for long-wavelength harmonics because of the smaller excitation of compressible eigenfunctions in the lower mantle. The relaxation times and velocity fields are more sensitive to higher viscosity contrast for the exponentially varying viscosity than for models with discrete jumps in the viscosity structure.

Wu, J.↗

Method for predicting impulsive noise generated by wind turbine rotors

Large wind turbines can generate both broad band and impulsive noises. These noises can be controlled by proper choice of rotor design parameters such as rotor location with respect to the supporting tower, tower geometry and tip speed. A method was developed to calculate the impulsive noise generated when the wind turbine blade experiences air forces that are periodic functions of the rotational frequency. This phenomenon can occur when the blades operate in the wake of the support tower and the nonuniform velocity field near the ground due to wind shear. Results from this method were compared with measured sound spectra taken at locations of one to two rotor diameters from the DOE/NASA Mod-1 wind turbine. The calculated spectra generally agreed with the measured data in both the amplitude of the predominant harmonics and the roll of rate with frequency. Measured sound pressure levels far from the Mod-1 (15 rotor diameters), however, were higher than predicted. Simultaneous measurements in the near and far field indicated that the propagation effects could enhance the sound levels by more than 10 dB above that expected by spherical dispersion. These propagation effects are believed to be due to terrain and atmospheric characteristics of the Mod-1 site.

Viterna, L. A.↗

Stellar convection. III - Convection at large Rayleigh numbers

The results of a numerical, three-dimensional study of convection in a self-gravitating sphere of Boussinesq fluid with a Rayleigh number of 10-billion and a Prandtl number of 1 are presented. The velocity and temperature are computed by using spectral methods in the horizontal and finite-differencing in the radial directions. An eddy viscosity and diffusivity are needed to model the subresolution flow. For Rayleigh numbers much less than 10-billion the flows do not have well-defined inertial ranges, and an eddy viscosity and diffusivity cannot be assigned in this self-consistent manner. By computing the energy spectra as well as the detailed energy budgets as a function of wavenumber, it is shown that for Rs = 10-billion there is an inertial range for the modes corresponding to spherical harmonics with l greater than 6.

Marcus, P. S.↗

Topographic mapping of the moon

Contour maps of the moon have been compiled by photogrammetric methods that use stereoscopic combinations of all available metric photographs from the Apollo 15, 16, and 17 missions. The maps utilize the same format as the existing NASA shaded-relief Lunar Planning Charts (LOC-1, -2, -3, and -4), which have a scale of 1:2,750,000. The map contour interval is 500 m. A control net derived from Apollo photographs by Doyle and others was used for the compilation. Contour lines and elevations are referred to the new topographic datum of the moon, which is defined in terms of spherical harmonics from the lunar gravity field. Compilation of all four LOC charts was completed on analytical plotters from 566 stereo models of Apollo metric photographs that cover approximately 20 percent of the moon. This is the first step toward compiling a global topographic map of the moon at a scale of 1:5,000,000.

Wu, S. S. C.↗

Mars synthetic topographic mapping

Topographic contour maps of Mars are compiled by synthesizing data acquired from various scientific experiments of the Mariner 9 mission, including S-band radio occultation, the UV spectrometer, the IR radiometer, the IR interferometer spectrometer, and television imagery, as well as earth-based radar information. The entire planet is mapped at scales of 1:25,000,000 and 1:5,000,000 using Mercator, Lambert, and polar stereographic map projections. For the computation of map projections, a biaxial spheroid figure is adopted. The semimajor and semiminor axes are 3393.4 and 3375.7 km, respectively, with a polar flattening of 0.0052. For the computation of elevations, a topographic datum is defined by a gravity field described in terms of spherical harmonics of fourth order and fourth degree combined with a 6.1-mbar occultation pressure surface. This areoid can be approximated by a triaxial ellipsoid with semimajor axes of A = 3394.6 km and B = 3393.3 km and a semiminor axis of C = 3376.3 km. The semimajor axis A intersects the Martian surface at longitude 105 deg W. The dynamic flattening of Mars is 0.00525. The contour interval of the maps is 1 km. For some prominent features where overlapping pictures from Mariner 9 are available, local contour maps at relatively larger scales were also compiled by photogrammetric methods on stereo plotters.

Wu, S. S. C.↗

3D reconstruction of an inertial-confinement fusion implosion with neural networks using multiple heterogeneous data sources

3D asymmetries are major degradation mechanisms in inertial-confinement fusion implosions at the National Ignition Facility (NIF). These asymmetries can be diagnosed and reconstructed with the neutron imaging system (NIS) on three lines of sight around the NIF target chamber. Conventional tomographic reconstructions are used to reconstruct the 3D morphology of the implosion using NIS [Volegov et al., J. Appl. Phys. 127, 083301 (2020)], but the problem is ill-posed with only three imaging lines of sight. Asymmetries can also be diagnosed with the real-time neutron activation diagnostics (RTNAD) and the neutron time-of-flight (nToF) suite. Since the NIS, RTNAD, and nToF each sample a different part of the implosion using different physical principles, we propose that it is possible to overcome the limitations of too few imaging lines of sight by performing 3D reconstructions that combine information from all three heterogeneous data sources. This work presents a new machine learning-based reconstruction technique to do just this. By using a simple physics model and group of neural networks to map 3D morphologies to data, this technique can easily account for data of multiple different types. A simple proof-of-principle is presented, demonstrating that this technique can accurately reconstruct a hot-spot shape using synthetic primary neutron images and a hot-spot velocity vector. In particular, the hot-spot’s asymmetry, quantified as spherical harmonic coefficients, is reconstructed to within ±4% of the radius in 90% of test cases. In the future, this technique will be applied to actual NIS, RTNAD, and nToF data to better understand 3D asymmetries at the NIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Heat Flow and Magnetization in the Oceanic Lithosphere

Two aspects of the processing and interpretation of satellite measurements of the geomagnetic field are described. One deals with the extraction of the part of the geomagnetic field that originates from sources in the earth's atmosphere. The other investigates the possibility of using the thermal state of the oceanic lithosphere to further constrain modelling and interpretation of magnetic anomalies. It is shown that some of the magnetic signal in crustal anomaly maps can be an artifact of the mathematical algorithms that have been used to separate the crustal field from the observed data. Strong magnetic anomalies can be distorted but are probably real, but weak magnetic anomalies can arise from leakage of power from short wavelengths, and will also appear in anomaly maps as repetitions of the strong crustal anomaly. The distortion and the ghost anomalies follow the magnetic dip lines in a way that is similar to actual MAGSAT anomaly fields. This phenomenon will also affect the lower degree spherical harmonic terms in the power spectrum of the crustal field. A model of the magnetic properties of the oceanic crust that has been derived from direct measurements of the rock magnetic properties of oceanic rocks is presented. The average intensity of magnetization in the oceanic crust is not strong enough to explain magnetic anomalies observed over oceanic areas. This is the case for both near surface observations (ship and aeromagnetic data) and satellite altitude observations. It is shown that magnetic sources in the part of the upper mantle that is situated above the Curie isotherm, if sufficiently strong, can produce satellite magnetic anomalies that are comparable to MAGSAT data. The method developed for the study of depth to the Curie isotherm and magnetic anomalies can also be used in inverse modelling of satellite magnetic anomalies when the model is to be adjusted with an annihilator.

Hayling, Kjell Lennart↗