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,243 records · Page 69

Time variations of solar UV irradiance as measured by the SOLSTICE (UARS) instrument

An analysis is presented of solar ultraviolet irradiance measurements made by the SOLSTICE spectrometers on the Upper Atmosphere Research Satellite (UARS). Reported observations cover the wavelength interval 119-420 nm, and the analysis discussed here is for the time period 26 Nov 1991 to 31 Dec 1992, during which time solar activity decreased in intensity. At the time of peak activity, the average 27-day variation had a relative amplitude of about 8 percent at Ly-alpha, tailing off to about 0.6 percent at 260 nm. It is shown that over the spectral interval 119-260 nm, the relative 27-day harmonic was about a factor of two larger during the strongly disturbed as compared with the moderately disturbed period.

London, Julius↗

X 17 boson and the H 3 ( p , e + e – ) He 4 and He 3 ( n , e + e – ) He 4 processes: A theoretical analysis

The present work deals with e + – e – pair production in the four-nucleon system. We first analyze the process as a purely electromagnetic one in the context of a state-of-the-art approach to nuclear strong-interaction dynamics and nuclear electromagnetic currents, derived from chiral effective field theory ( χ EFT ) . Next, we examine how the exchange of a hypothetical low-mass boson would impact the cross section for such a process. We consider several possibilities, that this boson is either a scalar, pseudoscalar, vector, or axial particle. The ab initio calculations use exact hyperspherical-harmonics methods to describe the bound state and low-energy spectrum of the A = 4 continuum, and they fully account for initial state interaction effects in the 3 + 1 clusters. While electromagnetic interactions are treated to high orders in the chiral expansion, the interactions of the hypothetical boson with nucleons are modeled in leading-order χ EFT (albeit, in some instances, selected subleading contributions are also accounted for). We also provide an overview of possible future experiments probing pair production in the A = 4 system at a number of candidate facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Distortion Representation of Forecast Errors for Model Skill Assessment and Objective Analysis

We study a novel characterization of errors for numerical weather predictions. In its simplest form we decompose the error into a part attributable to phase errors and a remainder. The phase error is represented in the same fashion as a velocity field and will be required to vary slowly and smoothly with position. A general distortion representation allows for the displacement and a bias correction of forecast anomalies. In brief, the distortion is determined by minimizing the objective function by varying the displacement and bias correction fields. In the present project we use a global or hemispheric domain, and spherical harmonics to represent these fields. In this project we are initially focusing on the assessment application, restricted to a realistic but univariate 2-dimensional situation. Specifically we study the forecast errors of the 500 hPa geopotential height field for forecasts of the short and medium range. The forecasts are those of the Goddard Earth Observing System data assimilation system. Results presented show that the methodology works, that a large part of the total error may be explained by a distortion limited to triangular truncation at wavenumber 10, and that the remaining residual error contains mostly small spatial scales.

Hoffman, Ross N.↗

Infrared parameterization and simple climate models

The geographical distribution of the annual mean and the annual cycle in surface temperature and satellite-observed IR is examined, and the spherical harmonic representation of the data fields is exploited to demonstrate how variances and covariances are dominated by the largest space and time scales. The geographical distribution of the annual cycle in the T and IR fields is explored; the strong imprint of the continents in both is clearly evident. The influence of the cloudiness of seasonal precipitation regimes on the IR annual cycle is also quite striking, especially over the subtropics. Analysis of the data shows that the simple form IR = A+BT (with A = 204 W/sq m and B = 1.93 W/sq m/K) explains 90 percent of the area-weighted variance in the annual mean and annual cycle of the zonally averaged IR field.

Short, D. A.↗

Integrated CFD modeling of gas turbine combustors

3D, curvilinear, multi-domain CFD analysis is becoming a valuable tool in gas turbine combustor design. Used as a supplement to experimental testing. CFD analysis can provide improved understanding of combustor aerodynamics and used to qualitatively assess new combustor designs. This paper discusses recent advancements in CFD combustor methodology, including the timely integration of the design (i.e. CAD) and analysis (i.e. CFD) processes. Allied Signal's F124 combustor was analyzed at maximum power conditions. The assumption of turbulence levels at the nozzle/swirler inlet was shown to be very important in the prediction of combustor exit temperatures. Predicted exit temperatures were compared to experimental rake data, and good overall agreement was seen. Exit radial temperature profiles were well predicted, while the predicted pattern factor was 25 percent higher than the harmonic-averaged experimental pattern factor.

Fuller, E. J.↗

Variability of accretion disks surrounding black holes: The role of inertial-acoustic mode instabilities

The global nonlinear time-dependent evolution of the inertial-acoustic mode instability in accretion disks surrounding black holes has been investigated. The viscous stress is assumed to be proportional to the gas pressure only, i.e., tau = alphap(sub g). It is found that an oscillatory nonsteady behavior exists in the inner regions of disks (r is less than 10r(sub g) where r(sub g) is the Schwarzschild radius) for sufficiently large alpha(greater than or approximately equal to 0.2) and for mass accretion rates less than about 0.3 times the Eddington value. The variations of the integrated bolometric luminosity from the disk, Delta L/L, are less than 3%. A power spectrum analysis of these variations reveals a power spectrum which can be fitted to a power-law function of the frequency Pis proportional to f(exp -gamma), with index gamma = 1.4-2.3 and a low-frequency feature at about 4 Hz in one case. In addition, a narrow peak centered at a frequency corresponding to the maximum epicyclic frequency of the disk at approximately 100-130 Hz and its first harmonic is also seen. The low-frequency modulations are remarkably similar to those observed in black hole candidate systems. The possible existence of a scattering corona in the inner region of the disk and/or other processes contributing to the power at high frequencies in the inner region of the accretion disk may make the detection of the high-frequency component difficult.

Chen, Xingming↗

Lunar dynamics and selenodesy - Results from analysis of VLBI and laser data

Very long base line interferometry (VLBI) observations of lunar radio transmitters have been combined with data from laser ranging to lunar retroreflectors to estimate simultaneously (1) parameters in models of the lunar orbit and libration and (2) the selenodetic coordinates of the radio transmitters and retroreflectors. For the ratio of the mass of the sun to that of the earth plus moon, a value of approximately 328,900.5 is obtained. The lunar moment-of-inertia ratios beta and gamma are also determined. The values and most of the uncertainties obtained for the third-degree harmonics of the moon's gravity field are comparable to those which have been obtained by others from observations of lunar-orbiting spacecraft. However, the results for the two third-degree gravity coefficients have much less uncertainty than determinations based on laser data alone. For the relative position vectors of the lunar radio transmitters and the retroreflectors, the estimates have uncertainties of about 30 m along the earth-moon direction and about 10 m in each of the two transverse coordinates.

King, R. W.↗

Prediction and measurement of low-frequency harmonic noise of a hovering model helicopter rotor

Far-field acoustic data for a model helicopter rotor have been gathered in a large open-jet, acoustically treated wind tunnel with the rotor operating in hover and out of ground-effect. The four-bladed Boeing 360 model rotor with advanced airfoils, planform, and tip shape was run over a range of conditions typical of today's modern helicopter main rotor. Near in-plane acoustic measurements were compared with two independent implementations of classical linear theory. Measured steady thrust and torque were used together with a free-wake analysis (to predict the thrust and drag distributions along the rotor radius) as input to this first-principles theoretical approach. Good agreement between theory and experiment was shown for both amplitude and phase for measurements made in those positions that minimized distortion of the radiated acoustic signature at low-frequencies.

Aggarawal, H. R.↗

The advanced tokamak path to a compact net electric fusion pilot plant

Abstract Physics-based simulations project a compact net electric fusion pilot plant with a nuclear testing mission is possible at modest scale based on the advanced tokamak concept, and identify key parameters for its optimization. These utilize a new integrated 1.5D core-edge approach for whole device modeling to predict performance by self-consistently applying transport, pedestal and current drive models to converge fully non-inductive stationary solutions, predicting profiles and energy confinement for a given density. This physics-based approach leads to new insights and understanding of reactor optimization. In particular, the levering role of high plasma density is identified, which raises fusion performance and self-driven ‘bootstrap currents’, to reduce current drive demands and enable high pressure with net electricity at a compact scale. Solutions at 6–7 T, ∼4 m radius and 200 MW net electricity are identified with margins and trade-offs possible between parameters. Current drive comes from neutral beam and ultra-high harmonic (helicon) fast wave, though other advanced approaches are not ruled out. The resulting low recirculating power in a double null configuration leads to a divertor heat flux challenge that is comparable to ITER, though reactor solutions may require more dissipation. Strong H-mode access (x2 margin over L–H transition scalings) and ITER-like heat fluxes are maintained with ∼20%–60% core radiation, though effects on confinement need further analysis. Neutron wall loadings appear tolerable. The approach would benefit from high temperature superconductors, as higher fields would increase performance margins while potential for demountability may facilitate nuclear testing. However, solutions are possible with conventional superconductors. An advanced load sharing and reactive bucking approach in the device centerpost region provides improved mechanical stress handling. The prospect of an affordable test device which could close the loop on net-electric production and conduct essential nuclear materials and breeding research is compelling, motivating research to validate the techniques and models employed here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of acoustic radiation in a jet flow environment

The analysis of sound fields from arbitrary source distributions in terms of Legendre and spherical Hankel functions is well known. The purpose of this paper is to extend this classical method of analysis to environments such as jet flows where flow and flow gradients are inherently present. The wave-equation governing the radiation of sound in such an environment is derived. The steady state flow and flow gradients in the axial and transverse directions appear as coefficients in the terms of the wave-equation. A semi-numerical method is used to solve the wave-equation in terms of modified spherical harmonics yielding the phase velocities and the directivities of an infinite set of modes. The directivity of each mode is obtained in terms of modified Legendre functions by numerical integration. Some results of these directivity and phase-velocity calculations are presented for a limited number of frequency and flow parameters. Both convective and shear refraction are shown to be important.

Mungur, P.↗

Control Design Strategies to Enhance Long-Term Aircraft Structural Integrity

Over the operational lifetime of both military and civil aircraft, structural components are exposed to hundreds of thousands of low-stress repetitive load cycles and less frequent but higher-stress transient loads originating from maneuvering flight and atmospheric gusts. Micro-material imperfections in the structure, such as cracks and debonded laminates, expand and grow in this environment, reducing the structural integrity and shortening the life of the airframe. Extreme costs associated with refurbishment of critical load-bearing structural components in a large fleet, or altogether reinventoring the fleet with newer models, indicate alternative solutions for life extension of the airframe structure are highly desirable. Increased levels of operational safety and reliability are also important factors influencing the desirability of such solutions. One area having significant potential for impacting crack growth/fatigue damage reduction and structural life extension is flight control. To modify the airframe response dynamics arising from command inputs and gust disturbances, feedback loops are routinely applied to vehicles. A dexterous flight control system architecture senses key vehicle motions and generates critical forces/moments at multiple points distributed throughout the airframe to elicit the desired motion characteristics. In principle, these same control loops can be utilized to influence the level of exposure to harmful loads during flight on structural components. Project objectives are to investigate and/or assess the leverage control has on reducing fatigue damage and enhancing long-term structural integrity, without degrading attitude control and trajectory guidance performance levels. In particular, efforts have focused on the effects inner loop control parameters and architectures have on fatigue damage rate. To complete this research, an actively controlled flexible aircraft model and a new state space modeling procedure for crack growth have been utilized. Analysis of the analytical state space model for crack growth revealed the critical mathematical factors, and hence the physical mechanism they represent, that influenced high rates of airframe crack growth. The crack model was then exercised with simple load inputs to uncover and expose key crack growth behavior. To characterize crack growth behavior, both "short-term" laboratory specimen test type inputs and "long-term" operational flight type inputs were considered. Harmonic loading with a single overload revealed typical exponential crack growth behavior until the overload application, after which time the crack growth was retarded for a period of time depending on the overload strength. An optimum overload strength was identified which leads to maximum retardation of crack growth. Harmonic loading with a repeated overload of varying strength and frequency again revealed an optimum overload trait for maximizing growth retardation. The optimum overload strength ratio lies near the range of 2 to 3 with dependency on frequency. Experimental data was found to correlate well with the analytical predictions.

Newman, Brett A.↗

Higher Order Multipole Analysis of Beam Crabbing Mode at 197 MHz for Electron-Ion Collider

In the local crabbing scheme for Electron Ion Collider interaction region, beam rotation is initiated by one set of crab cavities at a high beta region with near 90 deg phase advance to the interaction point (IP), and ceased by an identical set of cavities on the symmetrically opposite side. The fundamental mode of all crab cavities is designed to apply a momentum kick to the particle bunches to complete this half cycle of rotation in the horizontal plane, meanwhile the electromagnetic (EM) field can be expanded into a series of harmonics. Despite the lowest dipole mode, higher harmonics cause instabilities during beam operation and loss of particles. We used SimTrack, which is the 6-d symplectic element-by-element tracking code developed at Brookhaven National Laboratory in 2015, to study the threshold of each higher-order harmonic multipolar mode with the criteria of maintaining a dynamic aperture of 6 σ or greater for the Hadron Storage Ring (HSR). In this paper, we showed the threshold for each higher harmonic multipole and their comparison with the multipole strength for the current crab cavity design.

43 PARTICLE ACCELERATORS↗

Jarzynski Equality for Conditional Stochastic Work

It has been established that the inclusive work for classical, Hamiltonian dynamics is equivalent to the two-time energy measurement paradigm in isolated quantum systems. However, a plethora of other notions of quantum work has emerged, and thus the natural question arises whether any other quantum notion can provide motivation for purely classical considerations. In the present analysis, we propose the conditional stochastic work for classical, Hamiltonian dynamics, which is inspired by the one-time measurement approach. This novel notion is built upon the change of expectation value of the energy conditioned on the initial energy surface. As main results, we obtain a generalized Jarzynski equality and a sharper maximum work theorem, which account for how non-adiabatic the process is. Our findings are illustrated with the parametric harmonic oscillator.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Saturn gravity results obtained from Pioneer 11 tracking data and earth-based Saturn satellite data

Improved gravity coefficients for Saturn, its satellites and rings are calculated on the basis of a combination of Pioneer 11 spacecraft Doppler tracking data and earth-based determinations of Saturn natural satellite apse and node rates. Solutions are first obtained separately from the coherent Doppler tracking data obtained for the interval from August 20 to September 4, surrounding the time of closest approach, with the effects of solar plasma on radio signal propagation taken into account, and from secular rates for Mimas, Enceladus, Tethys, Dione, Rhea and Titan determined from astrometric data by Kozai (1957, 1976) and Garcia (1972). Combination of the data by the use of the Pioneer solution and corresponding unadjusted covariance matrix as a priori information for a secular rate analysis results in values for the total ring mass of essentially zero at a standard error level of 1.7 x 10 to the -6th Saturn masses, a ratio of solar mass to that of the Saturn system of 3498.09 + or - 0.22, masses of Rhea, Titan and Iapetus of 4.0 + or - 0.9, 238.8 + or - 3, and 3.4 + or - 1.3 x 10 to the -6th Saturn masses, respectively, and second and fourth zonal harmonics of 16,479 + or - 18 and -937 + or - 38, respectively. The harmonic coefficients are noted to be important as boundary conditions in the modeling of the Saturn interior.

Null, G. W.↗

The 1983 June 15 occultation by Neptune. II - The oblateness of Neptune

The oblateness and radius of Neptune were determined from an analysis of photoelectric observations of the June 15, 1983 occultation by Neptune at six stations, combined with the results of the Apr. 7, 1968 Neptune occultation of BD - 17 deg 4388. The oblateness is 0.0191 + or - 0.0017 and the equatorial radius is 25,268 + or - 12 km at the level probed by the occultation. The results are consistent with recent determinations of Neptune's rotation period and J2 (the second-order gravitational harmonic coefficient) and suggest that Neptune is far less centrally condensed than Uranus. Temperature profiles for Neptune's upper atmosphere were also derived from these data.

French, R. G.↗

Triggering of longitudinal combustion instabilities in solid rocket motors: Nonlinear combustion response

Pulsed oscillations in solid rocket motors are investigated with emphasis on nonlinear combustion response. The study employs a wave equation governing the unsteady motions in a two-phase flow, and a solution technique based on spatial- and time-averaging. A wide class of combustion response functions is studied to second-order in fluctuation amplitude to determine if, when, and how triggered instabilities arise. Conditions for triggering are derived from analysis of limit cycles, and regions of triggering are found in parametric space. Based on the behavior of model dynamical systems, introduction of linear cross-coupling and quadratic self-coupling among the acoustic modes appears to be the manner in which the nonlinear combustion response produces triggering to a stable limit cycle. Regions of initial conditions corresponding to stable pulses were found, suggesting that stability depends on initial phase angle and harmonic content, as well as the composite amplitude, of the pulse.

Wicker, J. M.↗

Performance Improvements for the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. Griffin includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. A code assessment performed in FY-20 identified two significant issues with the transport solvers in Griffin: first, the primary heterogeneous SN (discrete ordinates) transport solver based on continuous finite element methods required significant mesh refinement and higher memory usage compared to solvers based on the method of characteristic for equivalent accuracy. Second, the homogeneous PN (spherical harmonics expansion) transport solver did not adequately support polynomial refinement, which is a feature usually required for problems with spatial homogenization and pronounced streaming, typical in fast or gas-cooled reactor systems. To address the first issue, the development effort focused on the more promising discontinuous finite element method (DFEM)-based SN transport solver in Griffin. The addition of an asynchronous parallel transport sweeper and coarse mesh finite difference (CMFD) acceleration have rendered a superior heterogeneous SN transport capability for multiphysics problems that requires far less computing resources in terms of both CPU time and memory usage. This is demonstrated with typical thermal- and fast-spectrum reactor benchmark problems, including 2D Transient Reactor Test, 3D Advanced Burner Test Reactor (ABTR), and 2D and 3D Empire microreactor. For the second issue, the development effort focused on a new transport solver based on the hybrid finite element PN method (HFEM-PN), equivalent to the variational nodal method, as well as a new diffusion solver based on HFEM-Diffusion. This solver is intended for homogenized domains with multiphysics coupling (i.e., supports mesh displacement, seamless temperature feedback, etc.). Initial calculations with the HFEM-Diffusion implementation show very good parallel efficiency for the residual evaluations with the 2D ABTR benchmark. A future development effort will be centered on further improvements to the CMFD, HFEM-PN, and DFEM diffusion solvers to ensure Griffin meets performance and software quality assurance requirements for advanced reactor design and analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Geodetic and oceanographic results from the analysis of 1 year of Geosat data

The radial error model of Engelis (1987) is extended using the altimeter data obtained by the Geosat Exact Repeat Mission for simultaneous measurements of the stationary sea surface topography and geopotential coefficients. It is shown that the proper utilization of the Geosat altimeter data significantly improved the quality of field modeling at wave numbers by providing, simultaneously, the solution for improved orbits and the oceanographic information on the global scale. A new set of Geosat geophysical data records with the improved orbits was created, which also contains the values of the geoid undulation computed from a recent spherical harmonic expansion to degree 360.

Denker, Heiner↗