Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “harmonic analysis”

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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

DOE Repository Metadata Profile (DRMP): A Metadata Framework for Advancing Interoperability and AI Readiness Across Scientific Repositories

The Department of Energy (DOE) funds a diverse and distributed ecosystem of repositories that steward scientific data, publications, and software across its research programs, user facilities, and national laboratories. While significant progress has been made in standardizing dataset-level metadata, the metadata describing repositories themselves (their identity, governance, access interfaces, policies, and technical capabilities) remains inconsistent and fragmented across DOE-funded systems. This variability limits discoverability, interoperability, automated validation, and AI-driven analysis, all of which are increasingly essential for modern scientific workflows. To address this gap, the DOE Data Curation Working Group (DCWG) developed the DOE Repository Metadata Profile (DRMP). The DRMP is a practical, community-driven framework that defines how repositories can describe themselves in a consistent, machine-actionable, and scalable manner. The DRMP is not a new metadata schema. Instead, it is a mapping profile and structured element set capturing the essential characteristics of DOE repositories. It harmonizes repository-level metadata across six widely adopted community schemas: RE3Data; DCAT-US v3; Schema.org; Dublin Core; DataCite 4.6; and PREMIS 3.0. This harmonization eliminates reinvention and enables interoperability within DOE and across the broader scientific ecosystem. A core objective of the DRMP is to reduce burden on repositories by allowing them to reuse their existing metadata through a Rosetta-style crosswalk rather than redesigning local implementations. The profile introduces a three-level conformance model that supports incremental adoption: • Level 1 – Minimum Viable Record (MVR): foundational identification elements required for workflows, project registration, and basic repository presence. • Level 2 – Interoperable: structured metadata enabling alignment with national and international discovery systems. • Level 3 – AI-Ready: enhanced provenance, policy transparency, fixity, semantic context, and capabilities that support automated reasoning, model training governance, and machine-assisted curation. To support implementation, the DRMP includes JSON Schema definitions, OpenAPI patterns, and MCP templates that allow repositories to publish machine-readable metadata directly within existing platforms. These resources are modular and lightweight, enabling adoption without major architectural change. Adopting the DRMP enables repositories to: • Enhance discoverability and interoperability by aligning identifiers, classifications, and descriptive elements across widely used schema standards. • Support federated discovery and cross-registration across DOE systems, Data.gov, and international catalogs. • Enable AI agents and workflow orchestration systems to interpret repository-level metadata within the American Science Cloud (AmSC) through Model Context Protocol (MCP)-based context publication. • Demonstrate alignment with DOE’s open science, stewardship, and FAIR data priorities. This guidance represents a community-driven step forward. Through voluntary adoption and continued feedback, the DRMP advances a cohesive, machine-actionable description of DOE repositories that supports FAIR data practices, preparing the infrastructure for AI-enabled research, and strengthening the discoverability and reuse of DOE’s scientific outputs.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
NASA NTRS1983

NASTRAN forced vibration analysis of rotating cyclic structures

Theoretical aspects of a new capability, developed and added to the general purpose finite element program NASTRAN Level 17.7, to conduct forced vibration analysis of turned cyclic structures rotating about their axis of symmetry, are presented. The effects of Coriolis and centripetal accelerations as well as those due to the translational acceleration of the axis of rotation, are included. The equations of motion are first derived for an arbitrary grid point of the cyclic sector finite element model and then extended for the complete model. The equations are solved by four principal steps: (1) transformation of applied loads at frequency-dependent circumferential harmonic components; (2) application of circumferential harmonic-dependent intersegment compatibility constraints; (3) solution of frequency-dependent circumferential harmonic components of displacements; and (4) recovery of frequency-dependent response in various segments of the total structure. Five interrelated examples are presented to illustrate the various features of the development.

Elchuri, V.↗
NASA NTRS1986

The onset of time-dependent convection in spherical shells as a clue to chaotic convection in the earth's mantle

This work presents a detailed numerical study of the dynamical behavior of convection in a spherical shell, as applied to mantle convection. From both two-dimensional (120 radial and 360 tangential points) and three-dimensional (60 radial levels and spherical harmonics up to order and degree L = 33, m = 33), it is shown that for a spherical shell (with inner-to-outer radii ratio eta = 0.62) convection becomes time-dependent, with l = 2 dominating, at a Rayleigh number of about 31 times supercritical for a constant-viscosity, base-heated configuration. This secondary instability is characterized by oscillatory time dependence, with higher frequencies involved, at slightly higher Rayleigh numbers. In illustrating the onset of time dependence, the analysis is extended to show that the onset of weak turbulence in spherical-shell convection takes place at about 60 times the critical Rayleigh number via a quasi-periodic mode.

Machetel, Philippe↗
NASA NTRS1956

Lift hysteresis at stall as an unsteady boundary-layer phenomenon

Analysis of rotating stall of compressor blade rows requires specification of a dynamic lift curve for the airfoil section at or near stall, presumably including the effect of lift hysteresis. Consideration of the magnus lift of a rotating cylinder suggests performing an unsteady boundary-layer calculation to find the movement of the separation points of an airfoil fixed in a stream of variable incidence. The consideration of the shedding of vorticity into the wake should yield an estimate of lift increment proportional to time rate of change of angle of attack. This increment is the amplitude of the hysteresis loop. An approximate analysis is carried out according to the foregoing ideas for a 6:1 elliptic airfoil at the angle of attack for maximum lift. The assumptions of small perturbations from maximum lift are made, permitting neglect of distributed vorticity in the wake. The calculated hysteresis loop is counterclockwise. Finally, a discussion of the forms of hysteresis loops is presented; and, for small reduced frequency of oscillation, it is concluded that the concept of a viscous "time lag" is appropriate only for harmonic variations of angle of attack with time at mean conditions other than maximum lift.

Moore, Franklin K↗
NASA NTRS1993

New Mars free-air and Bouguer gravity: Correlation with topography, geology and large impact basins

Free-air and Bouguer gravity anomalies from a 50x50 field (MGM635), derived at the Goddard Space Flight Center, with global topography, geology, and the distribution of large impact basins was compared. The free-air gravity anomalies were derived from re-analysis of Viking Orbiter and Mariner 9 tracking data and have a spatial resolution of 250-300 km. Bouguer anomalies were calculated using a 50x50 expansion of the current Mars topography and the GSFC degree 50 geoid as the equipotential reference surface. Rotational flattening was removed using a moment of inertia of 0.365 and the corrections from Table B2 of Sleep and Phillips. Crustal density and mean density were assumed to be 2.9 and 3.93 gm/cm(sup 3). The spherical harmonic topography used has zero mean elevation, and differs from the USGS maps by about 2 km. Comparisons with global geology use a simplified map with about 1/3 the number of units on the current maps. For correlation with impact basins, the recent compilation by Schultz and Frey was used.

Frey, Herbert↗
NASA NTRS2014

Forced Response Analysis of a Fan with Boundary Layer Inlet Distortion

Boundary layer ingesting propulsion systems have the potential to significantly reduce fuel burn for future generations of commercial aircraft, but these systems must be designed to overcome the challenge of high dynamic stresses in fan blades due to forced response. High dynamic stresses can lead to high cycle fatigue failures. High-fidelity computational analysis of the fan aeromechanics is integral to an ongoing effort to design a boundary layer ingesting inlet and fan for a wind-tunnel test. An unsteady flow solution from a Reynoldsaveraged Navier Stokes analysis of a coupled inlet-fan system is used to calculate blade unsteady loading and assess forced response of the fan to distorted inflow. Conducted prior to the mechanical design of a fan, the initial forced response analyses performed in this study provide an early look at the levels of dynamic stresses that are likely to be encountered. For the boundary layer ingesting inlet, the distortion contains strong engine order excitations that act simultaneously. The combined effect of these harmonics was considered in the calculation of the forced response stresses. Together, static and dynamic stresses can provide the information necessary to evaluate whether the blades are likely to fail due to high cycle fatigue. Based on the analyses done, the overspeed condition is likely to result in the smallest stress margin in terms of the mean and alternating stresses. Additional work is ongoing to expand the analyses to off-design conditions, on-resonance conditions, and to include more detailed modeling of the blade structure.

Vibrational stress↗
NASA NTRS1987

A new hydrodynamic analysis of double layers

A genuine two-fluid model of plasmas with collisions permits the calculation of dynamic (not necessarily static) electric fields and double layers inside of plasmas including oscillations and damping. For the first time a macroscopic model for coupling of electromagnetic and Langmuir waves was achieved with realistic damping. Starting points were laser-produced plasmas showing very high dynamic electric fields in nonlinear force-produced cavitous and inverted double layers in agreement with experiments. Applications for any inhomogeneous plasma as in laboratory or in astrophysical plasmas can then be followed up by a transparent hydrodynamic description. Results are the rotation of plasmas in magnetic fields and a new second harmonic resonance, explanation of the measured inverted double layers, explanation of the observed density-independent, second harmonics emission from laser-produced plasmas, and a laser acceleration scheme by the very high fields of the double layers.

Hora, Heinrich↗
NASA NTRS2011

Suzaku Observations of the HMXB 1A 1118-61

We present broad band analysis of the Be/X-ray transient lA 1118-61 by Suzaku at the peak of its 3rd observed outburst in January 2009 and 2 weeks later when the source flux had decayed by an order of magnitude. The continuum was modeled with a cutoffpl model as well as a compTT model, with both cases requiring an additional black body component at lower energies, We confirm the detection of a cyclotron line at approx. 55 keV and discuss the possibility of a first harmonic at approx. 110 keV. Pulse profile comparisons show a change in the profile structure at lower energies, an indication for possible changes in the accretion geometry, Phase resolved spectroscopy in the outburst data show a change in the continuum throughout the pulse period. The decrease in the CRSF centroid energy also indicates that the viewing angle on the accretion column is changing throughout the pulse period. Subject headings: X-rays: stars - X-rays: binaries - stars: pulsars: individual (lA 1118-61) - stars: magnetic fields

Suchy, Slawomir↗
DOE OSTI2021

OVERVIEW OF RECENT PU-238 PRODUCTION ACTIVITIES AT IDAHO NATIONAL LABORATORY

The Plutonium Fuel Services (PFS) program at Idaho National Laboratory (INL) is active in the qualification of irradiation targets containing Np-237 for irradiation in the Advanced Test Reactor (ATR) to produce Pu-238 for future NASA missions. INL qualified and loaded 7 targets in ATR’s South Flux Trap (SFT) for cycle 169A, which occurred in Spring 2021. This program was reinitiated after two baseline production targets in three positions validated significant production of Pu-238 [ref.1]. The validation model was followed by the PFS-1 experimental test in the ATRC (Critical) facility [ref.2]. This paper outlines the progress and status of the PFS program. The qualification effort, safety analysis, hardware status, and future activities for qualification of an updated target design for use in the ATR will be discussed.

36 MATERIALS SCIENCE↗
NASA NTRS2019

Stability and Spin-Orbit Resonance Analysis of Low Altitude Martian Orbits

Orbit stability has been thoughtfully studied in various celestial bodies. In particular the focus has been placed on the Moon, due to the unstable natural perturbations of its gravity field. The increasing interest in Mars orbiters brings the question of the likelihood of natural decay in low altitude regimes. This paper studies the change in shape of low altitude Mars orbits by carrying out large sets of numerical high fidelity simulations. Results showed that various configurations of the orbital elements gave perturbations that result-ed in unstable orbits. The paper also studies the potential causes of the observed unstable regions. First by taking a close look at zonal and tesseral harmonics to find the implications of Mars mass concentrations of the used gravity fields, and second by computing theoretical spin-orbit resonances to study their implications in the stability at low altitudes.

gravity perturbations↗
DOE OSTI2025

Solution of the linear wave-particle kinetic equation for global modes of arbitrary frequency in a tokamak

The linear response of a plasma to perturbations of arbitrary frequency and wavelength is derived for any axisymmetric magnetized toroidal plasma. An explicit transformation to action-angle coordinates is achieved using orthogonal magnetic coordinates and the Littlejohn Lagrangian, establishing the validity of this result to arbitrary order in normalized Larmor radius. The global resonance condition for compressional modes is clarified in more detail than in previous works, confirming that the poloidal orbit-average of the cyclotron frequency gives the desired result at lowest order in Larmor radius. The global plasma response to the perturbation at each resonance is captured by a poloidal and gyroaverage of the perturbing potential. A “global gyroaveraging” of the potential is a natural by-product of this analysis which takes into account the changing of the magnetic field over an orbit. The resonance condition depends on two arbitrary integers which completely separately capture the effects poloidal non-uniformity and finite Larmor radius in generating sidebands. We learn that poloidal sidebands generated for compressional modes are dominated by the change in gyrofrequency over the orbit, which is very different to shear modes where the gyrofrequency only contributes via a finite Larmor radius effect. This increases the number of bounce harmonics required to compute the linear drive, giving a more complicated resonance map. An example calculation is given comparing resonance of shear and compressional modes in a published DIII-D case.

Compressional↗
DOE OSTI2021

Modeling Framework for Evaluating Grid Disturbances

Recent research in the literature proposed parameter-based functional forms to quantify the impact of disturbances on electric power grids. This modeling was based on a concept of generalized grid disturbances that harmonized reliability and resilience. Three main stages—avoid, react, and recover—were identified, and to support practicality, were intended to be independently planned. An objective of the grid disturbance model was to fairly analyze the potential impacts of novel mitigations like transactive energy systems (TES) on grid reliability and resilience. Quantitative evaluation of a common performance metric (e.g., area under the curve of percentage customers online before, during, and after an event) over time was shown as an effective tool for such analysis. Such quantification may also enable system planners to assess the impact of appropriate mitigative actions for improving grid performance against potential disturbances. This is the focus of this report, i.e., presenting an evaluation modeling framework that can be utilized to evaluate grid performance against disturbances, under a variety of existing and possible new mitigative actions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1981

A time Fourier analysis of zonal averaged ozone heating rates

A time-Fourier analysis is presented for the yearly variation of the zonal averaged ozone heating rates in the middle atmosphere based on a model study. The ozone heating rates are determined by utilizing two-dimensional ozone distributions, the altitude and latitude, and by including the effect of the curved earth's atmosphere. In addition, assumptions are introduced to the yearly variations of the ozone distributions due to the lack of sufficient existing ozone data. Among other results, it is shown that the first harmonic component indicates that the heating rates are completely out of phase between the northern and southern hemispheres. The second Fourier component shows a symmetric pattern with respect to the equator, as well as five distinct local extreme values of the ozone heating rate. The third harmonic component shows a pattern close to that of the first component except in the regions above 70 deg between 45-95 km in both hemispheres.

Wang, P.-H.↗
DOE OSTI2025

Development and investigation of efficient resonance ionization mass spectrometry schemes of gadolinium

Resonance ionization mass spectrometry of gadolinium can be used for nuclear forensics and to further the understanding of stellar nucleosynthesis but has been used only a handful of times due to the high laser power required and interference from non-resonant ionization of molecules of other elements. Herein we present the development of two novel resonance ionization spectroscopy schemes for gadolinium that provide improvements in isotopic fractionation and ionization efficiency, respectively, opening new applications for gadolinium analysis. The schemes are demonstrated and compared in a mixed sample of gadolinium and neodymium.

and nuclear chemistry↗
NASA NTRS1991

The ionizing radiation environment of LDEF prerecovery predictions

The Long Duration Exposure Facility (LDEF) was exposed to several sources of ionizing radiation while in orbit. The principal ones were trapped belt protons and electrons, galactic cosmic rays, and albedo particles (protons and neutrons) from the atmosphere. Large solar flares in 1989 may have caused a small contribution. Prior to the recovery of the spacecraft, a number of calculations and estimates were made to predict the radiation exposure of the spacecraft and experiments. These were made to assess whether measurable radiation effects might exist, and to plan the analysis of the large number of radiation measurements available on the LDEF. Calculations and estimates of total dose, particle fluences, linear energy transfer spectra, and induced radioactivity were made. The principal sources of radiation is described, and the preflight predictions are summarized.

Watts, John W., Jr.↗
DOE OSTI2023

Topological Equivalence of VSI and CSI Commutation Cells and Its Application to Switching Resonance Analysis and Damper Design

This paper proposes the concept of the topological equivalence between the commutation cells of a voltage source inverter (VSI) and a current source inverter (CSI) using duality and equivalent circuit transformation techniques. A switching resonance of CSI is analyzed and an effective three-phase RC damper is derived from the dc-link RC damper of VSIs based on the topological equivalence of the VSI and CSI commutation cells. Furthermore, a new three-phase RC damper design framework utilizing the topological equivalence is proposed. Here, experimental results verify that the proposed three-phase RC damper design for CSI significantly reduces the switching resonance and improves the CSI’s output waveform total harmonics distortions (THDs).

42 ENGINEERING↗
DOE OSTI2023

Optimal Frequency-domain Analysis for Spacecraft Time Series: Introducing the Missing-data Multitaper Power Spectrum Estimator

While the Lomb–Scargle periodogram is foundational to astronomy, it has a significant shortcoming: the variance in the estimated power spectrum does not decrease as more data are acquired. Statisticians have a 60 yr history of developing variance-suppressing power spectrum estimators, but most are not used in astronomy because they are formulated for time series with uniform observing cadence and without seasonal or daily gaps. Here we demonstrate how to apply the missing-data multitaper power spectrum estimator to spacecraft data with uniform time intervals between observations but missing data during thruster fires or momentum dumps. The F-test for harmonic components may be applied to multitaper power spectrum estimates to identify statistically significant oscillations that would not rise above a white noise–based false alarm probability. Multitapering improves the dynamic range of the power spectrum estimate and suppresses spectral window artifacts. We show that the multitaper–F-test combination applied to Kepler observations of KIC 6102338 detects differential rotation without requiring iterative sinusoid fitting and subtraction. Significant signals reside at harmonics of both fundamental rotation frequencies and suggest an antisolar rotation profile. Next we use the missing-data multitaper power spectrum estimator to identify the oscillation modes responsible for the complex "scallop-shell" shape of the K2 light curve of EPIC 203354381. We argue that multitaper power spectrum estimators should be used for all time series with regular observing cadence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Window convolution of the galaxy clustering bispectrum

In galaxy survey analysis, the observed clustering statistics do not directly match theoretical predictions but rather have been processed by a window function that arises from the survey geometry including the sky footprint, redshift-dependent background number density and systematic weights. While window convolution of the power spectrum is well studied, for the bispectrum with a larger number of degrees of freedom, it poses a significant numerical and computational challenge. In this work, we consider the effect of the survey window in the tripolar spherical harmonic decomposition of the bispectrum and lay down a formal procedure for their convolution via a series expansion of configuration-space three-point correlation functions, which was first proposed by Sugiyama et al. (2019). We then provide a linear algebra formulation of the full window convolution, where an unwindowed bispectrum model vector can be directly premultiplied by a window matrix specific to each survey geometry. To validate the pipeline, we focus on the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) luminous red galaxy (LRG) sample in the South Galactic Cap (SGC) in the redshift bin 0.4 ≤ z ≤ 0.6. We first perform convergence checks on the measurement of the window function from discrete random catalogues, and then investigate the convergence of the window convolution series expansion truncated at a finite of number of terms as well as the performance of the window matrix. This work highlights the differences in window convolution between the power spectrum and bispectrum, and provides a streamlined pipeline for the latter for current surveys such as DESI and the Euclid mission.

79 ASTRONOMY AND ASTROPHYSICS↗