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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS2012

Prediction of BVI Noise for an Active Twist Rotor Using a Loosely Coupled CFD/CSD Method and Comparison to Experimental Data

Numerical predictions of the acoustic characteristics of an Active Twist Rotor (ATR), using two methods to compute the rotor blade aerodynamics and elastic blade motion are compared to experimental data from a wind tunnel test in the NASA Langley Transonic Dynamics Tunnel (TDT) in 2000. The first method, a loosely coupled iterative method, utilizes the Computational Fluid Dynamics (CFD) code OVERFLOW 2 and the Computational Structural Dynamics (CSD) code CAMRAD II. The second method utilizes the CAMRAD II free-wake model only. The harmonic active-twist control to the main rotor blade system is identified with three parameters - harmonic actuation frequency, actuation amplitude, and control phase angle. The resulting aerodynamics and blade motion data from the two methods are then used in the acoustics code PSU-WOPWOP to predict acoustic pressure on a spherical array of equally spaced observers surrounding the rotor. This spherical distribution of pressure is used to compute the sound power level representing baseline and actuated conditions. Sound power levels for three categories of noise are defined as - blade-vortex interaction sound power level (BVIPWL), low frequency sound power level (LFPWL), and overall sound power level, OAPWL. Comparisons with measured data indicate the CFD/CSD analysis successfully captures the trends in sound power levels and the effects of active-twist control at advance ratios of 0.14 and 0.17. The free-wake model predictions show inconsistent sound power levels relative to the trends in the experimental and CFD data. This paper presents the first ever comparison between CFD/CSD acoustic predictions for an active-twist rotor and experimental measurements.

Fogarty, David E.↗

NASA NTRS1975

Downward continuation of gravity information from satellite to satellite tracking or satellite gradiometry in local areas

Integral formulas in the parameter domain are used instead of a representation by spherical harmonics. The neglected regions will cause a truncation error. The application of the discrete form of the integral equations connecting the satellite observations with surface gravity anomalies is discussed in comparison with the least squares prediction method. One critical point of downward continuation is the proper choice of the boundary surface. Practical feasibilities are in conflict with theoretical considerations. The properties of different approaches for this question are analyzed.

Rummel, R.↗

DOE OSTI2022

Full-potential KKR within the removed-sphere method: A practical and accurate solution to the Poisson equation

Here an efficient and accurate generalization of the removed-sphere method (RSM) to solve the Poisson equation for total charge density in a solid with space-filling convex Voronoi polyhedra (VPs) and any symmetry is presented. The generalized RSM avoids the use of multipoles and VP shape functions for cellular integrals, which have associated ill-convergent large, double-internal L sums in spherical-harmonic expansions, so that fast convergence in single-L sums is reached. Our RSM adopts full Ewald formulation to work for all configurations or when symmetry breaking occurs, such as for atomic displacements or elastic constant calculations. The structure-dependent coefficients A L that define RSM can be calculated once for a fixed structure and speed up the whole self-consistent-field procedure. The accuracy and rapid convergence properties are confirmed using two analytic models, including the Coulomb potential and energy. We then implement the full-potential RSM using the Green's function Korringa-Kohn-Rostoker (KKR) method for real applications and compare the results with other first-principle methods and experimental data, showing that they are equally as accurate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1979

The magnetic field of Mercury

Data from Mariner 10 observations of Mercury indicate that there exists an intrinsic magnetic field of the planet, sufficiently strong at present to deflect the solar wind flow around the planet and to form a detached bow shock wave in the super Alfvenic solar wind. Four methods used to analyze the magnetic field data and derive quantitative values for the description of the planetary field include (1) comparison of bow shock and magnetopause relative positions at Mercury to those at Earth; (2) direct spherical harmonic analysis of the data; (3) modeling of the magnetosphere by an image dipole and infinite 2-D current sheet in addition to the planetary field; and (4) scaling of a mathematical model for the terrestrial magnetosphere. The results obtained yield dipole moments ranging from 2.4 to 5.1x10 to the 22d power, with the lower values associated cw cm with certain models using partial quadrupole and octupole terms to improve the least squares fitting of models to observations.

Ness, N. F.↗

NASA NTRS1980

Protostellar formation in rotating interstellar clouds. I - Numerical methods and tests

Attention is given to numerical methods and tests of a series of gravitational hydrodynamics computer codes constructed in order to numerically follow the dynamic collapse of spherically symmetric (1-D), axisymmetric (2-D), and non-axisymmetric (3-D) isothermal interstellar clouds. A spherical harmonic expansion is used to solve the Poisson equation for the gravitational potential. It is shown that the use of explicit donor-cell hydrodynamics on a moving spherical coordinate grid ensures mass and momentum conservation and allows the grid to follow the collapse of the fluid. Finally, the performance of the codes is examined.

Boss, A. P.↗

NASA NTRS1983

On determining the large-scale ocean circulation from satellite altimetry

It is contended that a spherical harmonic expansion of the difference between the altimeter-derived mean sea surface and the geoid estimate should reveal the large-scale circulation of the ocean surface layer when the low-degree terms are examined. Methods based on this principle are proposed and partially demonstrated over the Pacific Ocean with the aid of the mean sea surface derived from the Seasat altimeter and the Goddard Earth Model 9 earth gravity model. The preliminary results reveal a well-defined clockwise gyre in the North Pacific and a much less well defined counterclockwise gyre in the South Pacific. When the dynamic topography thus obtained is compared with Wyrtki's (1975) dynamic topography derived from hydrographic data, the agreement is found to be within the limit of geoid uncertainties and satellite orbital errors.

Tai, C.-K.↗

NASA NTRS1993

Transition operators in electromagnetic-wave diffraction theory. II - Applications to optics

The theory developed by Hahne (1992) for the diffraction of time-harmonic electromagnetic waves from fixed obstacles is briefly summarized and extended. Applications of the theory are considered which comprise, first, a spherical harmonic expansion of the so-called radiation impedance operator in the theory, for a spherical surface, and second, a reconsideration of familiar short-wavelength approximation from the new standpoint, including a derivation of the so-called physical optics method on the basis of quasi-planar approximation to the radiation impedance operator, augmented by the method of stationary phase. The latter includes a rederivation of the geometrical optics approximation for the complete Green's function for the electromagnetic field in the presence of a smooth- and a convex-surfaced perfectly electrically conductive obstacle.

Hahne, G. E.↗

NASA NTRS1989

Radial deformation of the earth by oceanic tidal loading

A high-degree spherical harmonic series is used to compute the radial deformation of the Earth by oceanic tidal loading. By exploiting fast numerical transforms, this approach is found to be much more efficient, but no less accurate, than the traditional Green's function approach. The method is used to derive an atlas of load tide maps for 10 constitutents of the NSWC ocean tide model.

Ray, R. D.↗

DOE OSTI2022

Multiview Incomplete Knowledge Graph Integration with application to cross-institutional EHR data harmonization

Objective: The growing availability of electronic health records (EHR) data opens opportunities for integrative analysis of multi-institutional EHR to produce generalizable knowledge. A key barrier to such integrative analyses is the lack of semantic interoperability across different institutions due to coding differences. We propose a Multiview Incomplete Knowledge Graph Integration (MIKGI) algorithm to integrate information from multiple sources with partially overlapping EHR concept codes to enable translations between healthcare systems. Methods: The MIKGI algorithm combines knowledge graph information from (i) embeddings trained from the co-occurrence patterns of medical codes within each EHR system and (ii) semantic embeddings of the textual strings of all medical codes obtained from the Self-Aligning Pretrained BERT (SAPBERT) algorithm. Due to the heterogeneity in the coding across healthcare systems, each EHR source provides partial coverage of the available codes. MIKGI synthesizes the incomplete knowledge graphs derived from these multi-source embeddings by minimizing a spherical loss function that combines the pairwise directional similarities of embeddings computed from all available sources. MIKGI outputs harmonized semantic embedding vectors for all EHR codes, which improves the quality of the embeddings and enables direct assessment of both similarity and relatedness between any pair of codes from multiple healthcare systems. Results: With EHR co-occurrence data from Veteran Affairs (VA) healthcare and Mass General Brigham (MGB), MIKGI algorithm produces high quality embeddings for a variety of downstream tasks including detecting known similar or related entity pairs and mapping VA local codes to the relevant EHR codes used at MGB. Based on the cosine similarity of the MIKGI trained embeddings, the AUC was 0.918 for detecting similar entity pairs and 0.809 for detecting related pairs. For cross-institutional medical code mapping, the top 1 and top 5 accuracy were 91.0% and 97.5% when mapping medication codes at VA to RxNorm medication codes at MGB; 59.1% and 75.8% when mapping VA local laboratory codes to LOINC hierarchy. When trained with 500 labels, the lab code mapping attained top 1 and 5 accuracy at 77.7% and 87.9%. MIKGI also attained best performance in selecting VA local lab codes for desired laboratory tests and COVID-19 related features for COVID EHR studies. Compared to existing methods, MIKGI attained the most robust performance with accuracy the highest or near the highest across all tasks. Conclusions: The proposed MIKGI algorithm can effectively integrate incomplete summary data from biomedical text and EHR data to generate harmonized embeddings for EHR codes for knowledge graph modeling and cross-institutional translation of EHR codes.

Zhou, Doudou↗

NASA NTRS2001

A Spectral Algorithm for Solving the Relativistic Vlasov-Maxwell Equations

A spectral method algorithm is developed for the numerical solution of the full six-dimensional Vlasov-Maxwell system of equations. Here, the focus is on the electron distribution function, with positive ions providing a constant background. The algorithm consists of a Jacobi polynomial-spherical harmonic formulation in velocity space and a trigonometric formulation in position space. A transform procedure is used to evaluate nonlinear terms. The algorithm is suitable for performing moderate resolution simulations on currently available supercomputers for both scientific and engineering applications.

Shebalin, John V.↗

NASA NTRS1980

Equivalent source modeling of the main field using MAGSAT data

The software was considerably enhanced to accommodate a more comprehensive examination of data available for field modeling using the equivalent sources method by (1) implementing a dynamic core allocation capability into the software system for the automatic dimensioning of the normal matrix; (2) implementing a time dependent model for the dipoles; (3) incorporating the capability to input specialized data formats in a fashion similar to models in spherical harmonics; and (4) implementing the optional ability to simultaneously estimate observatory anomaly biases where annual means data is utilized. The time dependence capability was demonstrated by estimating a component model of 21 deg resolution using the 14 day MAGSAT data set of Goddard's MGST (12/80). The equivalent source model reproduced both the constant and the secular variation found in MGST (12/80).

Source record↗

NASA NTRS1973

Mars gravity derived from the long-period motion of Mariner 9

The spherical harmonic coefficients in the Mars gravity field are determined using a two-stage filtering process. In the first stage, a weighted least-squares orbit determination processor is applied to Mariner 9 radio data to estimate the mean Kepler elements of each orbit. A finite difference method is then used to generate orbital element rates. In the second stage, the gravity coefficients are determined using another weighted least-squares processor which fits the long period Lagrange perturbation equations to the estimated rates. The resulting gravity potential and associated equipotential surface is presented and error sources discussed. Gravity field results are applied to various physical properties of Mars, such as moments of inertia.-

Ferrari, A. J.↗

NASA NTRS2013

Covariance Analysis Tool (G-CAT) for Computing Ascent, Descent, and Landing Errors

G-CAT is a covariance analysis tool that enables fast and accurate computation of error ellipses for descent, landing, ascent, and rendezvous scenarios, and quantifies knowledge error contributions needed for error budgeting purposes. Because GCAT supports hardware/system trade studies in spacecraft and mission design, it is useful in both early and late mission/ proposal phases where Monte Carlo simulation capability is not mature, Monte Carlo simulation takes too long to run, and/or there is a need to perform multiple parametric system design trades that would require an unwieldy number of Monte Carlo runs. G-CAT is formulated as a variable-order square-root linearized Kalman filter (LKF), typically using over 120 filter states. An important property of G-CAT is that it is based on a 6-DOF (degrees of freedom) formulation that completely captures the combined effects of both attitude and translation errors on the propagated trajectories. This ensures its accuracy for guidance, navigation, and control (GN&C) analysis. G-CAT provides the desired fast turnaround analysis needed for error budgeting in support of mission concept formulations, design trade studies, and proposal development efforts. The main usefulness of a covariance analysis tool such as G-CAT is its ability to calculate the performance envelope directly from a single run. This is in sharp contrast to running thousands of simulations to obtain similar information using Monte Carlo methods. It does this by propagating the "statistics" of the overall design, rather than simulating individual trajectories. G-CAT supports applications to lunar, planetary, and small body missions. It characterizes onboard knowledge propagation errors associated with inertial measurement unit (IMU) errors (gyro and accelerometer), gravity errors/dispersions (spherical harmonics, masscons), and radar errors (multiple altimeter beams, multiple Doppler velocimeter beams). G-CAT is a standalone MATLAB- based tool intended to run on any engineer's desktop computer.

Boussalis, Dhemetrios↗

NASA NTRS1989

Downward continuation of the free-air gravity anomalies to the ellipsoid using the gradient solution and terrain correction: An attempt of global numerical computations

The formulas for the determination of the coefficients of the spherical harmonic expansion of the disturbing potential of the earth are defined for data given on a sphere. In order to determine the spherical harmonic coefficients, the gravity anomalies have to be analytically downward continued from the earth's surface to a sphere-at least to the ellipsoid. The goal is to continue the gravity anomalies from the earth's surface downward to the ellipsoid using recent elevation models. The basic method for the downward continuation is the gradient solution (the g sub 1 term). The terrain correction was also computed because of the role it can play as a correction term when calculating harmonic coefficients from surface gravity data. The fast Fourier transformation was applied to the computations.

Wang, Y. M.↗

NASA NTRS1981

Some new methods in geomagnetic field modeling applied to the 1960 - 1980 epoch

The utilization of satellite and surface data together permitted the incorporation of a solution for the anomaly field at each observatory. The residuals of the observatory measurements to such models is commensurate with the actual measurment accuracy. Incorporation of the anomaly estimation enabled the inclusion of stable time derivatives of the spherical harmonic coefficients up to the third derivative. A spherical harmonic model is derived with degree and order 13 in its constant and first time derivative terms, six in its second derivative terms and four in its third derivative terms.

Langel, R. A.↗

NASA NTRS1994

Venus spherical harmonic gravity model to degree and order 60

The Magellan and Pioneer Venus Orbiter radiometric tracking data sets have been combined to produce a 60th degree and order spherical harmonic gravity field. The Magellan data include the high-precision X-band gravity tracking from September 1992 to May 1993 and post-aerobraking data up to January 5, 1994. Gravity models are presented from the application of Kaula's power rule for Venus and an alternative a priori method using surface accelerations. Results are given as vertical gravity acceleration at the reference surface, geoid, vertical Bouguer, and vertical isostatic maps with errors for the vertical gravity and geoid maps included. Correlation of the gravity with topography for the different models is also discussed.

Konopliv, Alex S.↗

NASA NTRS1963

The World Magnetic Survey

The mathematical and graphical description of the earth's main field has been, and is, a "data limited" problem. The World Magnetic Survey (WMS) is an endeavor to minimize this limitation by rapidly and comprehensively blanketing the earth with magnetic field measurements. Satellite surveys, which will play a key role in the W MS, are the principal topic of this paper. Existing magnetic field descriptions, the expected results from new surveys, and the methods of obtaining these results with the POGO satellite are emphasized. It is anticipated on the basis of extrapolation from Vanguard 3 results and other considerations that a factor of 10 improvement will be obtained. This means that the average errors of 1 to 3 percent now present in field charts and spherical harmonic descriptions should be reduced to 0.1 to 0.3 percent as a result of the survey.

MATHEMATICS↗

NASA NTRS1973

Geodesic synchrotron radiation

This paper presents the results and methods of computing the high-frequency radiation emitted by freely falling particles moving in circular geodesic orbits in a spherically symmetric gravitational field. The high-frequency radiation, to which the methods of this paper apply, is the principal part of radiated energy only in the case of a particle moving in a highly relativistic, and therefore unstable, circular geodesic. The geodesic synchrotron radiation emitted in this case shows excitation of high-frequency harmonics and a narrow angular distribution. A Green's-function solution of the scalar wave equation is obtained using WKB methods. For application to relativistic circular orbits, a parabolic WKB approximation is required and yields solutions in terms of parabolic cylinder functions.

Breuer, R. A.↗