Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model structural error”

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,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Baryon acoustic oscillations from the complete SDSS-III Lyα-quasar cross-correlation function at z = 2.4

We present a measurement of baryon acoustic oscillations (BAO) in the cross-correlation of quasars with the Lyα-forest flux-transmission at a mean redshift z = 2.40. The measurement uses the complete SDSS-III data sample: 168,889 forests and 234,367 quasars from the SDSS Data Release DR12. In addition to the statistical improvement on our previous study using DR11, we have implemented numerous improvements at the analysis level allowing a more accurate measurement of this cross-correlation. We also developed the first simulations of the cross-correlation allowing us to test different aspects of our data analysis and to search for potential systematic errors in the determination of the BAO peak position. We measure the two ratios D H (z = 2.40)=r d = 9.01 ± 0.36 and D M (z = 2.40)=r d = 35.7 ±1.7, where the errors include marginalization over the non-linear velocity of quasars and the metal - quasar cross-correlation contribution, among other effects. These results are within 1.8σ of the prediction of the flat-ΛCDM model describing the observed CMB anisotropies.We combine this study with the Lyα-forest auto-correlation function (Bautista et al. 2017), yielding D H (z = 2.40)=r d = 8.94 ± 0.22 and D M (z = 2.40)=r d = 36.6 ± 1.2, within 2.3σ of the same flat-ΛCDM model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Ultra-fast computation of electronic spectra for large systems by tight-binding based simplified Tamm-Dancoff approximation (sTDA-xTB)

The computational bottleneck of the extremely fast simplified Tamm-Dancoff approximated (sTDA) time-dependent density functional theory procedure [S. Grimme, J. Chem. Phys. 138, 244104 (2013)] for the computation of electronic spectra for large systems is the determination of the ground state Kohn-Sham orbitals and eigenvalues. This limits such treatments to single structures with a few hundred atoms and hence, e.g., sampling along molecular dynamics trajectories for flexible systems or the calculation of chromophore aggregates is often not possible. The aim of this work is to solve this problem by a specifically designed semi-empirical tight binding (TB) procedure similar to the well established self-consistent-charge density functional TB scheme. The new special purpose method provides orbitals and orbital energies of hybrid density functional character for a subsequent and basically unmodified sTDA procedure. Compared to many previous semi-empirical excited state methods, an advantage of the ansatz is that a general eigenvalue problem in a non-orthogonal, extended atomic orbital basis is solved and therefore correct occupied/virtual orbital energy splittings as well as Rydberg levels are obtained. A key idea for the success of the new model is that the determination of atomic charges (describing an effective electron-electron interaction) and the one-particle spectrum is decoupled and treated by two differently parametrized Hamiltonians/basis sets. The three-diagonalization-step composite procedure can routinely compute broad range electronic spectra (0-8 eV) within minutes of computation time for systems composed of 500-1000 atoms with an accuracy typical of standard time-dependent density functional theory (0.3-0.5 eV average error). An easily extendable parametrization based on coupled-cluster and density functional computed reference data for the elements H–Zn including transition metals is described. The accuracy of the method termed sTDA-xTB is first benchmarked for vertical excitation energies of open- and closed-shell systems in comparison to other semi-empirical methods and applied to exemplary problems in electronic spectroscopy. As side products of the development, a robust and efficient valence electron TB method for the accurate determination of atomic charges as well as a more accurate calculation scheme of dipole rotatory strengths within the Tamm-Dancoff approximation is proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

GPS Altimetry

The advent of satellite altimetry has greatly improved our ability to observe global ocean circulation. However, the swath of a single, nadir-viewing satellite altimeter is only a few km and the track spacing is several hundred km to resolve the two-dimensional structure of ocean eddies. Our goal is to increase spatial and temporal coverage by monitoring Global Positioning System (GPS) signals reflected from the ocean. A constellation of spacecraft would each carry a GPS receiver capable of recording 8 reflections simultaneously. The reflections are well distributed in azimuth and elevation and can be tracked continuously while the satellite is in view, and another is then acquired, as illustrated below. The diagram depicts a new approach at altimetry measurements where ocean surface reflected GPS signals are simultaneously tracked and processed in a GPS flight receiver in space. The reflected GPS signals from the ocean must be compared precisely with the direct GPS signals in order to infer the characteristics of the ocean from the combined data set. Understanding the features and accuracy of GPS altimetry measurement is crucial to establishing its suitability for oceanography. Preliminary work has enabled us to theoretically model the signal output of the correlator for a variety of system parameters such as wind speed (sea roughness), receiver height, incidence angle, receiver range and Doppler filter bandwidth and antenna gain. Expected signal-to-noise ratio has been estimated from which we have inferred, to a first approximation, the basic receiver gain requirements for a space-based altimeter and the expected range raw error. In 1998, work on a different task led to the extraction of the first reflected GPS signal observed from a spaceborne receiver during the 1995 Space Transportation System-68 (STS-68) Shuttle Radar Laboratory-2 (SRL-2) high resolution synthetic aperture radar mission. Good comparisons with our signal models have been obtained. Having established that only modest signal-to-noise ratios are obtainable unless very high gain antennas are used, we expect that single measurements might not provide an estimate of sea state parameters as accurate as that obtainable with traditional remote sensing instruments, if costs are limited. Therefore spatial and temporal averaging of many measurements is required. Since the receiving satellite tracks do not repeat, measurements in a given area will be used to refine the local solution as a function of time, and will define the spatial resolution. In order to make this instrument viable in space, we need to detect and process many scattered signals. In 1998 we made a first step at understanding what configurations of antenna gain and orientations capture the largest number of viable signals to be used in the subsequent spatial/temporal averaging process. Additional information is contained in the original.

Lichten, Stephen M.↗
NASA NTRS2007

Staggering Inflation To Stabilize Attitude of a Solar Sail

A document presents computational-simulation studies of a concept for stabilizing the attitude of a spacecraft during deployment of such structures as a solar sail or other structures supported by inflatable booms. Specifically, the solar sail considered in this paper is a square sail with inflatable booms and attitude control vanes at the corners. The sail inflates from its stowed configuration into a square sail with four segments and four vanes at the tips. Basically, the concept is one of controlling the rates of inflation of the booms to utilize in mass-distribution properties to effect changes in the system s angular momentum. More specifically, what was studied were the effects of staggering inflation of each boom by holding it at constant length for specified intervals between intervals of increasing length until full length is reached. The studies included sensitivity analyses of effects of variations in mass properties, boom lengths, rates of increase in boom length, initial rates of rotation of the spacecraft, and several asymmetries that could arise during deployment. The studies led to the conclusion that the final attitude of the spacecraft could be modified by varying the parameters of staggered inflation. Computational studies also showed that by feeding back attitude and attitude-rate measurements so that corrective action is taken during the deployment, the final attitude can be maintained very closely to the initial attitude, thus mitigating the attitude changes incurred during deployment and caused by modeling errors. Moreover, it was found that by optimizing the ratio between the holding and length-increasing intervals in deployment of a boom, one could cause deployment to track a desired deployment profile to place the entire spacecraft in a desired attitude at the end of deployment.

Quadrelli, Marco↗
NASA NTRS2010

James Webb Space Telescope Studies of Dark Energy

The Hubble Space Telescope (HST) has contributed significantly to studies of dark energy. It was used to find the first evidence of deceleration at z=1.8 (Riess et al. 2001) through the serendipitous discovery of a type 1a supernova (SN1a) in the Hubble Deep Field. The discovery of deceleration at z greater than 1 was confirmation that the apparent acceleration at low redshift (Riess et al. 1998; Perlmutter et al. 1999) was due to dark energy rather than observational or astrophysical effects such as systematic errors, evolution in the SN1a population or intergalactic dust. The GOODS project and associated follow-up discovered 21 SN1a, expanding on this result (Riess et al. 2007). HST has also been used to constrain cosmological parameters and dark energy through weak lensing measurements in the COSMOS survey (Massey et al 2007; Schrabback et al 2009) and strong gravitational lensing with measured time delays (Suyu et al 2010). Constraints on dark energy are often parameterized as the equation of state, w = P/p. For the cosmological constant model, w = -1 at all times; other models predict a change with time, sometimes parameterized generally as w(a) or approximated as w(sub 0)+(1-a)w(sub a), where a = (1+z)(sup -1) is the scale factor of the universe relative to its current scale. Dark energy can be constrained through several measurements. Standard candles, such as SN1a, provide a direct measurement of the luminosity distance as a function of redshift, which can be converted to H(z), the change in the Hubble constant with redshift. An analysis of weak lensing in a galaxy field can be used to derive the angular-diameter distance from the weak-lensing equation and to measure the power spectrum of dark-matter halos, which constrains the growth of structure in the Universe. Baryonic acoustic oscillations (BAO), imprinted on the distribution of matter at recombination, provide a standard rod for measuring the cosmological geometry. Strong gravitational lensing of a time-variable source gives the angular diameter distance through measured time delays of multiple images. Finally, the growth of structure can also be constrained by measuring the mass of the largest galaxy clusters over cosmic time. HST has contributed to the study of dark energy through SN1a and gravitational lensing, as discussed above. HST has also helped to characterize galaxy clusters and the HST-measured constraints on the current Hubble constant H(sub 0) are relevant to the interpretation of dark energy measurements (Riess et al 2009a). HST has not been used to constrain BAO as the large number of galaxy redshifts required, of order 100 million, is poorly matched to HST's capabilities. As the successor to HST, the James Webb Space Telescope (JWST; Gardner et al 2006) will continue and extend HST's dark energy work in several ways.

Gardner, Jonathan P.↗
DOE OSTI2024

Momentum power spectrum of SDSS galaxies by massE cosmic ruler: 2.1×improvement in measure of growth rate

ABSTRACT Peculiar motion of galaxies probes the structure growth in the universe. In this study, we employ the galaxy stellar mass-binding energy (massE) relation with only two nuisance parameters to build the largest peculiar-velocity (PV) catalogue to date, consisting of 229 890 ellipticals from the main galaxy sample (MGS) of the Sloan Digital Sky Survey (SDSS). We quantify the distribution of the massE-based distances in individual narrow redshift bins (dz = 0.005), and then estimate the PV of each galaxy based on its offset from the Gaussian mean of the distribution. As demonstrated with the Uchuu-SDSS mock data, the derived PV and momentum power spectra are insensitive to accurate calibration of the massE relation itself, enabling measurements out to a redshift of 0.2, well beyond the current limit of z = 0.1 using other galaxy scaling laws. We then measure the momentum power spectrum and demonstrate that it remains almost unchanged if varying significantly the redshift bin size within which the distance is measured, as well as the intercept and slope of the massE relation, respectively. By fitting the spectra using the perturbation theory model with four free parameters, fσ8 is constrained to fσ8 = 0.459$^{+0.068}_{-0.069}$ over Δz = 0.02–0.2, 0.416$^{+0.074}_{-0.076}$ over Δz = 0.02–0.1, and 0.526$^{+0.133}_{-0.148}$ over Δz = 0.1–0.2. The error of fσ8 is 2.1 times smaller than that by the redshift space distortion (RSD) of the same sample. A Fisher matrix forecast illustrates that the constraint on fσ8 from the massE-based PV can potentially exceed that from the stage-IV RSD in late universe (z<0.5).

Shi, Yong (ORCID:0000000286146275)↗
NASA NTRS2002

The Gas Production Rate and Coma Structure of Comet C/1995 01 (Hale-Bopp)

The University of Wisconsin-Madison and NASA-Goddard conducted a comprehensive multi-wavelength observing campaign of coma emissions from comet Hale-Bopp, including OH 3080 A, [O I] 6300 A, H2O(+) 6158 A, H Balmer-alpha 6563 A, NH2 6330 A, [C I] 9850 A CN 3879 A, C2 5141 A, C3 4062 A, C I 1657 A, and the UV and optical continua. In this work, we concentrate on the results of the H2O daughter studies. Our wide-field OH 3080 A measured flux agrees with other, similar observations and the expected value calculated from published water production rates using standard H2O and OH photochemistry. However, the total [O I] 6300 A flux determined spectroscopically over a similar field-of-view was a factor of 3 - 4 higher than expected. Narrow-band [O I] images show this excess came from beyond the H2O scale length, suggesting either a previously unknown source of [O I] or an error in the standard OH + upsilon to O((sup I)D) + H branching ratio. The Hale-Bopp OH and [O I] distributions, both of which were imaged to cometocentric distances greater than 1 x 10(exp 6) km, were more spatially extended than those of comet Halley (after correcting for brightness differences), suggesting a higher bulk outflow velocity. Evidence of the driving mechanism for this outflow is found in the H(alpha) line profile, which was narrower than in comet Halley (though likely because of opacity effects, not as narrow as predicted by Monte-Carlo models). This is consistent with greater collisional coupling between the suprathermal H photodissociation products and Hale-Bopp's dense coma. Presumably because of mass loading of the solar wind by ions and ions by the neutrals, the measured acceleration of H2O(+) down the ion tail was much smaller than in comet Halley. Tailward extensions in the azimuthal distributions of OH 3080 A, [O I], and [C I], as well as a Doppler asymmetry in the [O I] line profile, suggest ion-neutral coupling. While the tailward extension in the OH can be explained by increased neutral acceleration, the [O I] 6300 A and [C I] 9850 A emissions show 13% and less than 200% excesses in this direction (respectively), suggesting a non-negligible contribution from dissociative recombination of CO(+) and/or electron collisional excitation. Thus, models including the effects of photo-and collisional chemistry are necessary for the full interpretation of these data.

Morgenthaler, Jeffrey P.↗
DOE OSTI2020

A Sheaf Theoretical Approach to Uncertainty Quantification of Heterogeneous Geolocation Information

Integration of multiple, heterogeneous sensors is a challenging problem across a range of applications. Prominent among these are multi-target tracking, where one must combine observations from different sensor types in a meaningful and efficient way to track multiple targets. Because different sensors have differing error models, we seek a theoretically justified quantification of the agreement among ensembles of sensors, both overall for a sensor collection, and also at a fine-grained level specifying pairwise and multi-way interactions among sensors. We demonstrate that the theory of mathematical sheaves provides a unified answer to this need, supporting both quantitative and qualitative data. Furthermore, the theory provides algorithms to globalize data across the network of deployed sensors, and to diagnose issues when the data do not globalize cleanly. We demonstrate and illustrate the utility of sheaf-based tracking models based on experimental data of a wild population of black bears in Asheville, North Carolina. A measurement model involving four sensors deployed among the bears and the team of scientists charged with tracking their location is deployed. This provides a sheaf-based integration model which is small enough to fully interpret, but of sufficient complexity to demonstrate the sheaf’s ability to recover a holistic picture of the locations and behaviors of both individual bears and the bear-human tracking system. A statistical approach was developed in parallel for comparison, a dynamic linear model which was estimated using a Kalman filter. This approach also recovered bear and human locations and sensor accuracies. When the observations are normalized into a common coordinate system, the structure of the dynamic linear observation model recapitulates the structure of the sheaf model, demonstrating the canonicity of the sheaf-based approach. However, when the observations are not so normalized, the sheaf model still remains valid.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Flow pattern and void fraction characterization in nitrogen/water flows through a diamond-type triply periodic minimal surface lattice

Here, this study presents, to the authors knowledge, the first experimental investigation on the void fraction and flow patterns in two-phase flows through a diamond-type Triply Periodic Minimal Surface (TPMS) lattice. An additively manufactured TPMS structure was tested under upward co-current flow of a water/nitrogen mixture. Superficial velocities varied for a total of 42 test conditions (gas: 0.01–2.4 m/s; liquid: 0.01–2.4 m/s; mass flux: 20–2370 kg/m 2 ·s). High-speed video and X-ray imaging enabled time-averaged void fraction measurements and identified six distinct flow regimes which were used to develop a flow pattern map. Comparison of the void fraction data with correlations from literature demonstrated the Rouhani and Axelsson (1970) [49] correlation modified by Steiner (1993) [52] provided the best agreement, which was improved with empirically fit coefficients. This approach predicted the void fraction with errors < ±20% for 64% of the data, with a mean absolute percent deviation (MAPD) of 23%. The measured frictional pressure drop was compared to correlations from literature which captured the observed trends but did not provide good accuracy. The best agreement was found after optimizing the empirical coefficients of the Muller-Steinhagen & Heck (1986) [56] correlation; this approach captured 33% of the data within ±20% with a MAPD of 46% over the full range and captured 71% of the data within ±20% a MAPD of 13.3% at mass fluxes >1100 kg/m 2 -s. These results provide foundational insight into TPMS two-phase flow behavior and inform modeling and design of advanced heat exchange components incorporating TPMS geometries.

42 - ENGINEERING↗
NASA NTRS2016

A Model-Model and Data-Model Comparison for the Early Eocene Hydrological Cycle

A range of proxy observations have recently provided constraints on how Earth's hydrological cycle responded to early Eocene climatic changes. However, comparisons of proxy data to general circulation model (GCM) simulated hydrology are limited and inter-model variability remains poorly characterised. In this work, we undertake an intercomparison of GCM-derived precipitation and P - E distributions within the extended EoMIP ensemble (Eocene Modelling Intercomparison Project; Lunt et al., 2012), which includes previously published early Eocene simulations performed using five GCMs differing in boundary conditions, model structure, and precipitation-relevant parameterisation schemes. We show that an intensified hydrological cycle, manifested in enhanced global precipitation and evaporation rates, is simulated for all Eocene simulations relative to the preindustrial conditions. This is primarily due to elevated atmospheric paleo-CO2, resulting in elevated temperatures, although the effects of differences in paleogeography and ice sheets are also important in some models. For a given CO2 level, globally averaged precipitation rates vary widely between models, largely arising from different simulated surface air temperatures. Models with a similar global sensitivity of precipitation rate to temperature (dP=dT ) display different regional precipitation responses for a given temperature change. Regions that are particularly sensitive to model choice include the South Pacific, tropical Africa, and the Peri-Tethys, which may represent targets for future proxy acquisition. A comparison of early and middle Eocene leaf-fossil-derived precipitation estimates with the GCM output illustrates that GCMs generally underestimate precipitation rates at high latitudes, although a possible seasonal bias of the proxies cannot be excluded. Models which warm these regions, either via elevated CO2 or by varying poorly constrained model parameter values, are most successful in simulating a match with geologic data. Further data from low-latitude regions and better constraints on early Eocene CO2 are now required to discriminate between these model simulations given the large error bars on paleoprecipitation estimates. Given the clear differences between simulated precipitation distributions within the ensemble, our results suggest that paleohydrological data offer an independent means by which to evaluate model skill for warm climates.

Boundary conditions↗
NASA NTRS2017

Program to Optimize Simulated Trajectories II (POST2) Surrogate Models for Mars Ascent Vehicle (MAV) Performance Assessment

The primary purpose of the multiPOST tool is to enable the execution of much larger sets of vehicle cases to allow for broader trade space exploration. However, this exploration is not achieved solely with the increased case throughput. The multiPOST tool is applied to carry out a Design of Experiments (DOE), which is a set of cases that have been structured to capture a maximum amount of information about the design space with minimal computational effort. The results of the DOE are then used to fit a surrogate model, ultimately enabling parametric design space exploration. The approach used for the MAV study includes both DOE and surrogate modeling. First, the primary design considerations for the vehicle were used to develop the variables and ranges for the multiPOST DOE. The final set of DOE variables were carefully selected in order to capture the desired vehicle trades and take into account any special considerations for surrogate modeling. Next, the DOE sets were executed through multiPOST. Following successful completion of the DOE cases, a manual verification trial was performed. The trial involved randomly selecting cases from the DOE set and running them by hand. The results from the human analyst's run and multiPOST were then compared to ensure that the automated runs were being executed properly. Completion of the verification trials was then followed by surrogate model fitting. After fits to the multiPOST data were successfully created, the surrogate models were used as a stand-in for POST2 to carry out the desired MAV trades. Using the surrogate models in lieu of POST2 allowed for visualization of vehicle sensitivities to the input variables as well as rapid evaluation of vehicle performance. Although the models introduce some error into the output of the trade study, they were very effective at identifying areas of interest within the trade space for further refinement by human analysts. The next section will cover all of the ground rules and assumptions associated with DOE setup and multiPOST execution. Section 3.1 gives the final DOE variables and ranges, while section 3.2 addresses the POST2 specific assumptions. The results of the verification trials are given in section 4. Section 5 gives the surrogate model fitting results, including the goodness-of-fit metrics for each fit. Finally, the MAV specific results are discussed in section 6.

Zwack, M. R.↗
DOE OSTI2021

EMISSA (Exploring Millimeter Indicators of Solar-Stellar Activity): I. The initial millimeter–centimeter main-sequence star sample

Due to their wide wavelength coverage across the millimeter to centimeter (mm–cm) range and their increased sensitivity, modern interferometric arrays facilitate observations of the thermal and non-thermal radiation that is emitted from different layers in the outer atmospheres of stars. We study the spectral energy distribution (S obs (ν)) of main-sequence stars based on archival observations in the mm–cm range with the aim to study their atmospheric stratification as a function of stellar type. The main-sequence stars with significant detection in mm bands were identified in the ALMA Science Archive. These data were then complemented with spectral flux data in the extreme ultraviolet to cm range as compiled from various catalogues and observatory archives. We compared the resultant S obs (ν) of each star with a photospheric emission model (S mod (ν)) calculated with the PHOENIX code. The departures of S obs (ν) from S mod (ν) were quantified in terms of a spectral flux excess parameter (ΔS/S mod ) and studied as a function of stellar type. The initial sample consists of 12 main-sequence stars across a broad range of spectral types from A1 to M3.5 and the Sun-as-a-star as reference. The stars with T eff = 3000–7000 K (F–M type) showed a systematically higher S obs (ν) than S mod (ν) in the mm–cm range. Their ΔS/S mod exhibits a monotonic rise with decreasing frequency. The steepness of this rise is higher for cooler stars in the T eff = 3000–7000 K range, although the single fully convective star (T eff ~ 3000 K) in the sample deviates from this trend. Meanwhile, Sobs(ν) of the A-type stars agrees with S mod (ν) within errors. The systematically high ΔS/S mod in F–M stars indicates hotter upper atmospheric layers, that is, a chromosphere and corona in these stars, like for the Sun. The mm–cm ΔS/S mod spectrum offers a way to estimate the efficiency of the heating mechanisms across various outer atmospheric layers in main-sequence stars, and thereby to understand their structure and activity. We emphasise the need for dedicated surveys of main-sequence stars in the mm–cm range.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Accurate shear estimation with fourth-order moments

ABSTRACT As imaging surveys progress in exploring the large-scale structure of the Universe through the use of weak gravitational lensing, achieving sub-per cent accuracy in estimating shape distortions caused by lensing, or shear, is imperative for precision cosmology. In this paper, we extend the Fourier power function shapelets (FPFS) shear estimator using fourth-order shapelet moments and combine it with the original second-order shear estimator to reduce galaxy shape noise. We calibrate this novel shear estimator analytically to a sub-per cent level-accuracy using the AnaCal framework. This higher order shear estimator is tested with realistic image simulations, and after analytical correction for the detection/selection bias and noise bias, the multiplicative shear bias $|m|$ is below $3\times 10^{-3}$ (99.7 per cent confidence interval) for both isolated and blended galaxies. Once combined with the second-order FPFS shear estimator, the shape noise is reduced by $\sim 35~{{\ \rm per\ cent}}$ for isolated galaxies in simulations with Hyper Suprime-Cam and Vera C. Rubin Observatory Legacy Survey of Space and Time observational conditions. However, for blended galaxies, the effective number density does not significantly improve with the combination of the two estimators. Based on these results, we recommend exploration of how this framework can further reduce the systematic uncertainties in shear due to point spread function leakage and modelling error, and potentially provide improved precision in shear inference in high-resolution space-based images.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Reanalysis of the 𝛽− $\overline{v}$ 𝑒 Angular Correlation Measurement from the aSPECT Experiment with New Constraints on Fierz Interference

On the basis of revisions of some of the systematic errors, we reanalyzed the electron-antineutrino angular correlation (𝑎 coefficient) in free neutron decay inferred from the recoil energy spectrum of the protons which are detected in 4⁢𝜋 by the aSPECT spectrometer. With 𝑎=−0.104 02⁢(82), the new value differs only marginally from the one published in 2020. Here, the experiment also has sensitivity to 𝑏, the Fierz interference term. From a correlated (𝑏,𝑎) fit to the proton recoil spectrum, we derive a limit of 𝑏 =−0.0098⁢(193) which translates into a somewhat improved 90% confidence interval region of −0.041 ≤ 𝑏 ≤ 0.022 on this hypothetical term. Tighter constraints on 𝑏 can be set from a combined analysis of the PERKEO III (𝛽 asymmetry) and aSPECT measurement which suggests a finite value of 𝑏 with 𝑏 (𝑐) =−0.0181 ± 0.0065 deviating by 2.82⁢𝜎 from the standard model.

Beta decay↗
DOE OSTI2022

Structural stability of thin overhanging walls during material extrusion additive manufacturing of thermoset-based ink

Recent developments have enabled material extrusion additive manufacturing of thermoset-based composite inks on the large scale. In addition, printing out-of-plane components is of broad interest to the polymer material extrusion community. Here, we address some of the challenges associated with both large-scale and out-of-plane thermoset material extrusion additive manufacturing by studying the height at which thin overhanging walls collapse. Walls at a range of overhang angles were printed until they collapsed. An optical camera captured the profile of each wall throughout the print, allowing the collapse height to be identified and the geometric fidelity to the programmed angle to be evaluated. Using previously measured rheological properties, predictive models were generated to approximate the collapse height and profile of the deflected walls. First, an analytical model was created to predict the height at which the walls would yield. The analytical model assumes the walls exhibit a perfectly linear profile; however, experiments proved this assumption to be false. Therefore, a finite element simulation was developed to account for the elastic deflection that occurs during printing. The finite element simulation predicts both the yield height and the deflected profile after the deposition of each layer. For the properties of the thermoset ink used here, the yield height predicted by the analytical model and finite element simulation are virtually identical. These predictions match experimental data reasonably well, but minor errors are observed. Accounting for the fully plastic moment appears to explain the small mismatch between experimental data and predictions. Additionally, the finite element simulation provides an excellent prediction of the deflected profile before the wall begins to collapse. Finally, by demonstrating that the collapse height and deflected profile of thin overhanging walls can be predicted, this work illustrates how the soft viscoelastic properties of thermoset-based composite inks limit the scale of a key feature required to print some nonplanar components. It also provides a basis to tailor in-process curing systems to suppress deflection and collapse of thin overhanging walls.

36 MATERIALS SCIENCE↗
DOE OSTI2017

On the adiabatic representation of Meyer-Miller electronic-nuclear dynamics

The Meyer-Miller (MM) classical vibronic (electronic + nuclear) Hamiltonian for electronically non-adiabatic dynamics—as used, for example, with the recently developed symmetrical quasiclassical (SQC) windowing model—can be written in either a diabatic or an adiabatic representation of the electronic degrees of freedom, the two being a canonical transformation of each other, thus giving the same dynamics. Although most recent applications of this SQC/MM approach have been carried out in the diabatic representation—because most of the benchmark model problems that have exact quantum results available for comparison are typically defined in a diabatic representation—it will typically be much more convenient to work in the adiabatic representation, e.g., when using Born-Oppenheimer potential energy surfaces (PESs) and derivative couplings that come from electronic structure calculations. The canonical equations of motion (EOMs) (i.e., Hamilton’s equations) that come from the adiabatic MM Hamiltonian, however, in addition to the common first-derivative couplings, also involve second-derivative non-adiabatic coupling terms (as does the quantum Schrödinger equation), and the latter are considerably more difficult to calculate. This paper thus revisits the adiabatic version of the MM Hamiltonian and describes a modification of the classical adiabatic EOMs that are entirely equivalent to Hamilton’s equations but that do not involve the second-derivative couplings. The second-derivative coupling terms have not been neglected; they simply do not appear in these modified adiabatic EOMs. This means that SQC/MM calculations can be carried out in the adiabatic representation, without approximation, needing only the PESs and the first-derivative coupling elements. Here, the results of example SQC/MM calculations are presented, which illustrate this point, and also the fact that simply neglecting the second-derivative couplings in Hamilton’s equations (and presumably also in the Schrödinger equation) can cause very significant errors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Modest set of simulations of local-type primordial non-Gaussianity in the DESI era

A measurement of a primordial non-Gaussianity (PNG) signal through late- or early-Universe probes has the potential to transform our understanding of the physics of the primordial Universe. While large-scale structure observables in principle contain vital information, interpreting these measurements is challenging due to poorly understood astrophysical effects. Luckily, N-body simulations, such as the abacuspng set presented in this study, consisting of 9 boxes, each of size L box = 2 Gpc/h and particle mass of 1.01 × 10 10 M ⊙ /h, provide a viable path forward. As validation, we find good agreement between the simulations and our expectations from one-loop perturbation theory (PT) and the "separate universe" method for the matter bispectrum, matter power spectrum and the halo bias parameter associated with PNG, b φ . As a science application, here we investigate the link between halo assembly bias and b φ for halo properties known to play a vital role in accurately predicting galaxy clustering: concentration, shear (environment), and accretion rate. We find a strong response for all three parameters, suggesting that the connection between b φ and the assembly history of halos needs to be taken into account by future PNG analyses. We further perform the first study of the b φ parameter from fits to early DESI data of the luminous red galaxy (LRG) and quasistellar object (QSO) samples and comment on the effect on f NL constraints for the allowed galaxy-halo models (note that σ[f NL ]∝ $\frac{σ[b_φ]}{b_φ}$). We find that the error on f NL is 15, 8, 7 for the LRGs at z = 0.5 and z = 0.8 and QSOs at z = 1.4, respectively, suggesting that a thorough understanding of galaxy assembly bias is warranted so as to perform robust high-precision analysis of local-type PNG with future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Anisotropic yield models for lattice unit cell structures exploiting orthotropic symmetry

We report numerical homogenization enables efficient computational analysis and design of multiscale structures made of micro-architected materials including lattice unit cells. However, predicting yield is nontrivial because it requires accurate and efficient predictions of the maximum stress inside the homogenized unit cells. To address this challenge, we develop a macroscale anisotropic yield function for micro-architected materials. The yield function depends on the three-dimensional macroscale stress state and the parameters describing a family of micro-architectures, such as the radii of the struts in a lattice unit cell. To ensure accuracy, we determine the maximum stress using three-dimensional continuum finite-element analysis. To ensure efficiency, we construct surrogate models from the aforementioned high-fidelity results for yield prediction. To reduce simulation costs and surrogate modeling complexity, we leverage orthotropic symmetry commonly found in lattice unit cells. In this paper, we provide a thorough presentation of group representations and the systematic procedure to exploit orthotropic domain symmetry in homogenization and surrogate modeling. We illustrate the use of linear homogenization results to predict yield in specific unit cells without further simulations. We furthermore show that surrogate modeling presents a viable option for anisotropic yield prediction in a continuously parametrized family of micro-architectures. More specifically, despite the dimensionality and degree of nonlinearity in the maximum micro von Mises stress within the unit cell, the surrogate models can predict it with less than 5% error at least 90% of the time. Moreover, the largest under-prediction error, which is more critical than the over-prediction error, is typically less than 10%.

42 ENGINEERING↗

Refine your search to explore more results.