Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Doubling 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 469 records · Page 26

Gravitational effective field theory islands, low-spin dominance, and the four-graviton amplitude

Abstract We analyze constraints from perturbative unitarity and crossing on the leading contributions of higher-dimension operators to the four-graviton amplitude in four spacetime dimensions, including constraints that follow from distinct helicity configurations. We focus on the leading-order effect due to exchange by massive degrees of freedom which makes the amplitudes of interest infrared finite. In particular, we place a bound on the coefficient of the R 3 operator that corrects the graviton three-point amplitude in terms of the R 4 coefficient. To test the constraints we obtain nontrivial effective field-theory data by computing and taking the large-mass expansion of the one-loop minimally-coupled four-graviton amplitude with massive particles up to spin 2 circulating in the loop. Remarkably, we observe that the leading EFT coefficients obtained from both string and one-loop field-theory amplitudes lie in small islands. The shape and location of the islands can be derived from the dispersive representation for the Wilson coefficients using crossing and assuming that the lowest-spin spectral densities are the largest. Our analysis suggests that the Wilson coefficients of weakly-coupled gravitational physical theories are much more constrained than indicated by bounds arising from dispersive considerations of 2 → 2 scattering. The one-loop four-graviton amplitudes used to obtain the EFT data are computed using modern amplitude methods, including generalized unitarity, supersymmetric decompositions and the double copy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

BAGELS for simultaneous polarization, orbit, and optics control in electron storage rings

We present a new method for minimizing the effects of radiative depolarization in electron storage rings by use of a minimal number of special vertical orbit bumps. The bumps can be used to minimize the effects of radiative depolarization while simultaneously maintaining other common benefits of vertical orbits, e.g., transverse coupling and vertical dispersion control. Because simultaneously optimizing the large number of vertical correctors in a ring is operationally infeasible, we use dimensionality reduction to define a minimal number of the most effective groups of vertical correctors that can be optimized during operation, motivating the name “Best Adjustment Groups for ELectron Spin” (BAGELS). The method is streamlined by using suitable “basis bumps” instead of all individual vertical correctors. We define three types of basis bumps for different purposes: (i) generates no delocalized transverse coupling nor delocalized vertical dispersion, (ii) generates no delocalized vertical dispersion, and (iii) generates no delocalized transverse coupling. BAGELS has been essential in the design of the Electron Storage Ring (ESR) of the Electron-Ion Collider (EIC) and will be beneficial for any polarized electron ring, including FCC-ee. HERA and LEP would have likely benefitted as well. We use BAGELS to significantly increase polarization in the 18 GeV EIC-ESR, beyond what is achievable with conventional methods; in the 1-IP lattice, we nearly double the asymptotic polarization, and in the 2-IP lattice, we more than triple the asymptotic polarization. We also use BAGELS to construct knobs that can be used for global coupling correction, and knobs that generate vertical emittance for beam size matching, all while having minimal impacts on the polarization and orbit/optics. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS↗

Experimental Measurement of Out-of-Time-Ordered Correlators at Finite Temperature

Out-of-time-ordered correlators (OTOCs) are a key observable in a wide range of interconnected fields including many-body physics, quantum information science, and quantum gravity. Measuring OTOCs using near-term quantum simulators will extend our ability to explore fundamental aspects of these fields and the subtle connections between them. Here, we demonstrate an experimental method to measure OTOCs at finite temperatures and use the method to study their temperature dependence. These measurements are performed on a digital quantum computer running a simulation of the transverse field Ising model. Our flexible method, based on the creation of a thermofield double state, can be extended to other models and enables us to probe the OTOC’s temperature-dependent decay rate. Measuring this decay rate opens up the possibility of testing the fundamental temperature-dependent bounds on quantum information scrambling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

ON THE ACCURACY OF MEASUREMENTS MADE UPON FILMS PHOTOGRAPHED BY BAKER-NUNN SATELLITE TRACKING CAMERAS

The photographs taken by the Baker-Nunn cameras of the Smithsonian Astrophysical Observatory were formerly measured with Van Biesbroeck goniometers, but are now measured exclusively with Mann two-screw comparators. Instead of measuring point images, we usually have to measure oblong images, either trails or breaks (the latter produced by interruption of the exposure). We can therefore expect that the accuracy of the determination of their positions is more influenced by the personal errors than it is in conventional astrometry. Approximately 800 measurements were made on 34 images of different length, and, from these, the frequency distribution of settings was determined. Then the relationship between this distribution and the length of the image (trail or break) was examined, and the relationship between the "magnitude error" and the length of image was determined. Measurements were made at different angles formed by the trails and the reticle lines, and it was found that the accuracy was not affected by the angle of the setting. Position determinations have been made using reference stars at different distances. On the basis of these results, we can conclude that the accuracy is neither influenced by the distortion of the emulsion nor by the optical distortion within an area of a diameter of S cm (5°.8). We believe that we are justified in using the linear plate constant method, at least within this area. From numerous double and multiple measurements, the following standard errors were determined for a single position.

ASTROMETRY↗

A fast invariant imbedding method for multiple scattering calculations and an application to equivalent widths of CO2 lines on Venus

The invariant imbedding method considered is based on an equation which describes the change in the reflected radiation when an optically thin layer is added to the top of the atmosphere. The equation is used to treat the problem of reflection from a planetary atmosphere as an initial value problem. A fast method is discussed for the solution of the invariant imbedding equation. The speed and accuracy of the new method are illustrated by comparing it with the doubling program published by Hansen and Travis (1974). Computations are performed of the equivalent widths of carbon dioxide absorption lines in solar radiation reflected by Venus for several models of the planetary atmosphere.

Sato, M.↗

Characterization of Interlaminar Crack Growth in Composites with the Double Cantilever Beam Specimen

A project to examine the double cantilever beam specimen as a quantitative test method to assess the resistance of various composite materials to interlaminar crack growth is discussed. A second objective is to investigate the micromechanics of failure for composites with tough matrix resins from certain generic types of polymeric systems: brittle thermosets, toughened thermosets, and thermoplastics. Emphasis is given to a discussion of preliminary results in two areas: the effects of temperature and loading rate for woven composites, and the effects of matrix toughening in woven and unidirectional composites.

Hunston, D. L.↗

The vibrating ribbon problem revisited

A revised formal solution of the vibrating ribbon problem of hydrodynamic stability is presented. The initial formulation of Gaster (1965) is modified by application of the Briggs method and a careful treatment of the complex double Fourier transform inversions. Expressions are obtained in a natural way for the discrete spectrum as well as for the four branches of the continuous spectra. These correspond to discrete and branch-cut singularities in the complex wavenumber plane. The solutions from the continuous spectra decay both upsteam and downstream of the ribbon, with the decay in the upstream direction being much more rapid than that in the downstream direction. Comments and clarification of related prior work are made.

Ashpis, David E.↗

Excitation of Continuous and Discrete Modes in Incompressible Boundary Layers

This report documents the full details of the condensed journal article by Ashpis & Reshotko (JFM, 1990) entitled "The Vibrating Ribbon Problem Revisited." A revised formal solution of the vibrating ribbon problem of hydrodynamic stability is presented. The initial formulation of Gaster (JFM, 1965) is modified by application of the Briggs method and a careful treatment of the complex double Fourier transform inversions. Expressions are obtained in a natural way for the discrete spectrum as well as for the four branches of the continuous spectra. These correspond to discrete and branch-cut singularities in the complex wave-number plane. The solutions from the continuous spectra decay both upstream and downstream of the ribbon, with the decay in the upstream direction being much more rapid than that in the downstream direction. Comments and clarification of related prior work are made.

Ashpis, David E.↗

Nominal 30-M Cropland Extent Map of Continental Africa by Integrating Pixel-Based and Object-Based Algorithms Using Sentinel-2 and Landsat-8 Data on Google Earth Engine

A satellite-derived cropland extent map at high spatial resolution (30-m or better) is a must for food and water security analysis. Precise and accurate global cropland extent maps, indicating cropland and non-cropland areas, is a starting point to develop high-level products such as crop watering methods (irrigated or rainfed), cropping intensities (e.g., single, double, or continuous cropping), crop types, cropland fallows, as well as assessment of cropland productivity (productivity per unit of land), and crop water productivity (productivity per unit of water). Uncertainties associated with the cropland extent map have cascading effects on all higher-level cropland products. However, precise and accurate cropland extent maps at high spatial resolution over large areas (e.g., continents or the globe) are challenging to produce due to the small-holder dominant agricultural systems like those found in most of Africa and Asia. Cloud-based Geospatial computing platforms and multi-date, multi-sensor satellite image inventories on Google Earth Engine offer opportunities for mapping croplands with precision and accuracy over large areas that satisfy the requirements of broad range of applications. Such maps are expected to provide highly significant improvements compared to existing products, which tend to be coarser in resolution, and often fail to capture fragmented small-holder farms especially in regions with high dynamic change within and across years. To overcome these limitations, in this research we present an approach for cropland extent mapping at high spatial resolution (30-m or better) using the 10-day, 10 to 20-m, Sentinel-2 data in combination with 16-day, 30-m, Landsat-8 data on Google Earth Engine (GEE). First, nominal 30-m resolution satellite imagery composites were created from 36,924 scenes of Sentinel-2 and Landsat-8 images for the entire African continent in 2015-2016. These composites were generated using a median-mosaic of five bands (blue, green, red, near-infrared, NDVI) during each of the two periods (period 1: January-June 2016 and period 2: July-December 2015) plus a 30-m slope layer derived from the Shuttle Radar Topographic Mission (SRTM) elevation dataset. Second, we selected Cropland/Non-cropland training samples (sample size 9791) from various sources in GEE to create pixel-based classifications. As supervised classification algorithm, Random Forest (RF) was used as the primary classifier because of its efficiency, and when over-fitting issues of RF happened due to the noise of input training data, Support Vector Machine (SVM) was applied to compensate for such defects in specific areas. Third, the Recursive Hierarchical Segmentation (RHSeg) algorithm was employed to generate an object-oriented segmentation layer based on spectral and spatial properties from the same input data. This layer was merged with the pixel-based classification to improve segmentation accuracy. Accuracies of the merged 30-m crop extent product were computed using an error matrix approach in which 1754 independent validation samples were used. In addition, a comparison was performed with other available cropland maps as well as with LULC maps to show spatial similarity. Finally, the cropland area results derived from the map were compared with UN FAO statistics. The independent accuracy assessment showed a weighted overall accuracy of 94, with a producers accuracy of 85.9 (or omission error of 14.1), and users accuracy of 68.5 (commission error of 31.5) for the cropland class. The total net cropland area (TNCA) of Africa was estimated as 313 Mha for the nominal year 2015.

Cropland mapping; cropland areas; 30-m; Landsat-8;↗

Development of a Combined Cohesive and Virtual Crack-Closure Technique Approach to Represent R-Curves

Resistance curve (R-curve) effects due to fiber bridging, crack migration, and other blunting mechanisms are common in composite laminates. These mechanisms can dramatically increase the resistance to crack propagation but can be challenging to simulate. Delamination with R-curves can be analyzed using cohesive zone modeling (CZM) or the virtual crack closure technique (VCCT). Large fracture process zones can be simulated with CZM, but they require highly refined meshes. Coarser meshes can be used with VCCT, but this method is only applicable to small fracture process zones with R-curves defined as functions of position. Therefore, a technique with the computational efficiency of VCCT and the natural ability of cohesive elements to represent large fracture process zones is desirable. An approach is proposed that starts as a CZM in which the cohesive traction separation law (TSL) is separated into high-strength (HS) and low-strength (LS) components. The HS part, which is responsible for the mesh requirements of the CZM analysis, is replaced by VCCT. The combined method is evaluated by analyzing the response of double cantilevered beam specimens. Two specimen layup configurations, [0/90/90/0] 3s and [0/90/90/0] 9s , are evaluated. The R-curve response from each configuration is determined and applied to the constitutive properties of models built with CZM and with VCCT. The results and computational efficiency of the CZM, VCCT, and combined approaches are compared. The results are indicative that a combined cohesive/VCCT approach can enable progressive failure analyses to retain the computational efficiency of VCCT with the ability of the cohesive elements to capture R-curve effects.

Virtual Crack-Closure Technique↗

Linear-scaling quadruple excitations in local pair natural orbital coupled-cluster theory

Here, we present a fast, asymptotically linear-scaling implementation of the perturbative quadruples energy correction in coupled-cluster theory using local natural orbitals. Our work follows the domain-based local pair natural orbital (DLPNO) approach previously applied to lower levels of excitations in coupled-cluster theory. Our DLPNO-CCSDT(Q) algorithm uses converged doubles and triples amplitudes from a preceding DLPNO-CCSDT computation to compute the quadruples amplitude and energy in the quadruples natural orbital (QNO) basis. We demonstrate the compactness of the QNO space, showing that more than 95% of the (Q) correction can be recovered using relatively loose natural orbital cutoffs, compared to the tighter cutoffs used in pair and triples natural orbitals at lower levels of coupled-cluster theory. We also highlight the accuracy of our algorithm in the computation of relative energies, which yields deviations of sub-kJ mol −1 in relative energy compared to the canonical CCSDT(Q). Timings are conducted on a series of growing linear alkanes (up to 10 carbons and 608 basis functions) and water clusters (up to 49 water molecules and 2842 basis functions) to establish the asymptotic linear-scaling of our DLPNO-(Q) algorithm.

Auxiliary functions↗

Ground and excited states of even-numbered Hubbard ring at half-filling: comparison of the extended Gutzwiller approach with exact diagonalization

It remains a great challenge in condensed matter physics to develop a method to treat strongly correlated many-body systems with balanced accuracy and efficiency. We introduce an extended Gutzwiller (EG) method incorporating a manifold technique, which builds an effective manifold of the many-body Hilbert space, to describe the ground- and excited-state properties of strongly correlated electrons. Here we systematically apply an EG projector onto the ground and excited states of a non-interacting system. Diagonalization of the true Hamiltonian within the manifold formed by the resulting EG wavefunctions gives the approximate ground and excited states of the correlated system. To validate this technique, we implement it on even-numbered fermionic Hubbard rings at half-filling with periodic boundary conditions, and compare the results with the exact diagonalization (ED) method. The EG method is capable of generating high-quality ground and low-lying excited state wavefunctions, as evidenced by the high overlaps of wavefunctions between the EG and ED methods. Favorable comparisons are also achieved for other quantities including the total energy, the double occupancy, the total spin and the staggered magnetization. With the capability of accessing the excited states, the EG method can capture the essential features of the one-electron removal spectral function that contains contributions from states deep in the excited spectrum. Finally, we provide an outlook on the application of this method on large extended systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Sensitivity of Simulations of Double-detonation Type Ia Supernovae to Integration Methodology

Abstract We study the coupling of hydrodynamics and reactions in simulations of the double-detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator/Strang splitting versus the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation time step to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Estimation of regional heat flux from scanning radiometer measurements

The paper presents a statistical method for extracting regional information on the radiative heat flux at the top of the atmosphere from a satellite-borne scanning radiometer. The radiative heat flux distribution is denoted by a double Fourier series with coefficients to be determined from measurements. The method determines, in the sense of minimum variance, the best linear estimator of the total flux over a region. The heat flux distribution is represented by a truncated series which is unbiased.

Bess, T. D.↗

Energy calibration of germanium detectors for the Majorana Demonstrator

The Majorana Demonstrator was a search for neutrinoless double-beta decay (0νββ) in the 76 Ge isotope. It was staged at the 4850-foot level of the Sanford Underground Research Facility (SURF) in Lead, SD. The experiment consisted of 58 germanium detectors housed in a low background shield and was calibrated once per week by deploying a 228 Th line source for 1 to 2 hours. The energy scale calibration determination for the detector array was automated using custom analysis tools. We describe the offline procedure for calibration of the Demonstrator germanium detectors, including the simultaneous fitting of multiple spectral peaks, estimation of energy scale uncertainties, and the automation of the calibration procedure.

47 OTHER INSTRUMENTATION↗

High-temperature fatigue in metals

The presentation focuses primarily on the progress we at NASA Lewis Research Center have made. The understanding of the phenomenological processes of high temperature fatigue of metals for the purpose of calculating lives of turbine engine hot section components is discussed. Improved understanding resulted in the development of accurate and physically correct life prediction methods such as Strain-Range Partitioning for calculating creep fatigue interactions and the Double Linear Damage Rule for predicting potentially severe interactions between high and low cycle fatigue. Examples of other life prediction methods are also discussed.

Gary R. Halford↗

High-temperature fatigue in metals - A brief review of life prediction methods developed at the Lewis Research Center of NASA

The presentation focuses primarily on the progress we at NASA Lewis Research Center have made. The understanding of the phenomenological processes of high temperature fatigue of metals for the purpose of calculating lives of turbine engine hot section components is discussed. Improved understanding resulted in the development of accurate and physically correct life prediction methods such as Strain-Range partitioning for calculating creep fatigue interactions and the Double Linear Damage Rule for predicting potentially severe interactions between high and low cycle fatigue. Examples of other life prediction methods are also discussed. Previously announced in STAR as A83-12159

Halford, G. R.↗