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 145 records · Page 8

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 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.↗

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.↗