SEARCH · Engineering Papers

Results for “model checking”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2023

ENSO-Induced Teleconnection: Process-Oriented Diagnostics to Assess Rossby Wave Sources and Ambient Flow Properties in Climate Models

Climate model fidelity in representing ENSO-induced teleconnection is assessed with process-oriented diagnostics that examine a chain of processes, from equatorial Pacific precipitation to the midlatitude circulation pattern over the Pacific–North American regions. Such processes are rarely addressed during model development. Using an upper-tropospheric divergent level, local vorticity gradient of the ambient zonal flow (∂ 2 U¯/∂y 2 ) and a restoring force for Rossby waves (β * ) are estimated, the equivalent barotropic vorticity equation is solved, and an anomalous Rossby wave source (RWS') quantified. The analysis is applied to AMIP5 and AMIP6 simulations. For a realistic circulation response representation, the hypothesis that models accurately represent the strength and location of RWS', and spatial variations in β * is tested. Compared to AMIP5, in AMIP6 there are clear improvements in representing RWS' and β * . To validate the hypothesis, the analysis identifies two metrics: spatially coherent RWS' in the subtropical North Pacific, and longitudes of negative β * over the western-central North Pacific. By projecting these metrics in two and three-dimensional views, improvements or degradations in model versions are apparent. If a model’s fidelity in representing ∂ 2 U¯/∂y 2 and RWS' are compromised, then radiated Rossby waves are reflected more equatorward, resulting in zonally elongated circulation anomalies over the central North Pacific. Thus, during climate model development, applying this analysis frequently will keep a regular check on the fidelity of the modeled response to anomalous El Niño convection in conjunction with changing model ambient flow dependencies. As a result, this analysis is intended to form a process-oriented diagnostics package, a community contribution to the NOAA Model Diagnostics Task Force.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Bayesian inference of fiber orientation and polymer properties in short fiber-reinforced polymer composites

Herein we present a Bayesian methodology to infer the elastic modulus of the constituent polymer and the fiber orientation state in a short-fiber reinforced polymer composite (SFRP). The properties are inversely determined using only a few experimental tests. Developing composite manufacturing digital twins for SFRP composite processes, including injection molding and extrusion deposition additive manufacturing (EDAM) requires extensive experimental material characterization. In particular, characterizing the composite mechanical properties is time consuming and therefore, micromechanics models are used to fully identify the elasticity tensor. Hence, the objective of this paper is to infer the fiber orientation and the effective polymer modulus and therefore, identify the elasticity tensor of the composite with minimal experimental tests. To that end, we develop a hierarchical Bayesian model coupled with a micromechanics model to infer the fiber orientation and the polymer elastic modulus simultaneously which we then use to estimate the composite elasticity tensor. We motivate and demonstrate the methodology for the EDAM process but the development is such that it is applicable to other SFRP composites processed via other methods. Our results demonstrate that the approach provides a reliable framework for the inference, with as few as three tensile tests, while accounting for epistemic and aleatory uncertainty. Posterior predictive checks show that the model is able to recreate the experimental data well. The ability of the Bayesian approach to calibrate the material properties and its associated uncertainties, make it a promising tool for enabling a probabilistic predictive framework for composites manufacturing digital twins.

36 MATERIALS SCIENCE↗

DOE OSTI2023

IDAES-PSE 2.1.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.1.0 Release Highlights New IDAES Examples Repository Starting with this release, the IDAES examples are developed in the new IDAES/examples repository. Along with many content and usability improvements, the most significant changes are: To install the examples, after installing IDAES, run pip install idaes-examples The idaes get-examples command, previously used for this, has been removed The HTML version is now available at https://idaes-examples.readthedocs.io The previous URL, https://idaes.github.io/examples-pse, will not be updated and may be removed at some point in the future For more details, refer to the resources available at IDAES/examples. Removal of Non-Functional Apps A review of the code in the idaes/apps and idaes/models_extra folders was undertaken, and a number of tools were identified as being outdated or non-functional and no longer supported by their development teams. Due to this, the following tools have been removed: idaes/apps/alamopy_depr (note that the new ALAMOpy interface remains avaialble in idaes/core/surrogates) idaes/apps/helmet idaes/apps/ripe idaes/apps/roundingRegression idaes/models_extra/carbon_capture Pyomo 6.6 This version of IDAES is the first requiring Pyomo 6.6. This version of Pyomo contains multiple internal improvements and refactorings. While for the majority of cases this should have positive or no impact on solvability of IDAES models, we are aware of a small number of models that have been affected as a result of these changes. For more information, refer to the Pyomo 6.6.1 release notes. Other highlights Model Initialization A prototype API for a new approach to initializing IDAES models is now available which makes available some new techniques for initializing models. This is documented in the Initializing Models Reference Guide Modular Properties Framework Support for some transport properties Helmholtz Equation of State properties Better error checking for case where unit models are set to include phase equilibrium but the property package is set to support only a single phase Multi-Stream Contactor model: a new base model for systems involving contacting of two or more streams with mass transfer. This model is intended to be used as the foundation for models such as membrane separators, solvent extraction and other similar processes. This is documented in the Multi-Stream Contactor Reference Guide idaes/models_extra/power_generation report() methods for unit models using Helmholtz equation of state General Code Maintenance Streamlining of dependencies and creation of new optional dependency groupings to support non-core tools General linting of codebase to ensure compliance with most pylint checks Spell checking of all code and doc strings Removal of backward compatibility code for Python 2

IDAES↗

DOE OSTI2025

Correlating processing variables to material properties in recycled polypropylene: A data‐driven approach

Abstract Polypropylene (PP) is one of the most widely used plastics, yet its recycling remains limited, with less than 1% of solid waste PP being reprocessed. Mechanical recycling through extrusion is the most practical method, but inconsistent reprocessing conditions introduce variability in material properties. While temperature, screw speed, and residence time influence the thermomechanical stress applied during reprocessing, there are no standardized guidelines for optimizing these parameters. This study examines how these factors shape the properties of recycled PP, using conditions designed to mimic post‐industrial recycled (PIR) scrap. Residence time was measured using colorimetric tracking and correlated with molecular weight, viscosity, and mechanical properties over multiple extrusion cycles. Data‐driven modeling, including response surface methodology, support vector machines, and artificial neural networks, identified processing temperature as the dominant factor in material degradation, followed by residence time. Mechanical properties remained stable, while viscosity decreased predictably with increasing residence time. By linking reprocessing conditions to property evolution, this study provides a method to optimize processing parameters and reduce variability in recycled PP. These findings help manufacturers improve process control, making recycled PP more predictable for reuse in manufacturing. Highlights Study of PIR‐quality PP without additives or compatibilizers. Residence time analysis shows processing temperature drives PP property changes. Mark‐Houwink enables quick molecular weight checks for quality control. Models predict mechanical and rheological shifts in reprocessing. Optimized processing parameters minimize property degradation in recycling.

Estela‐García, John E. [Polymer Engineering Center↗

DOE OSTI2026

Seawater Acidification and Bubble Plume Dispersion from Accidental Subsea CO 2 Pipeline Rupture: A Multiphase CFD Study

If a CO 2 reservoir or transmission pipeline were to leak, both the surrounding ecology and maritime traffic safety could be put at risk. To better understand and prepare for this risk, multiphase Computational Fluid Dynamics (CFD) models were built in ANSYS Fluent to capture the behavior of a leak once it enters the water. A 3D Eulerian–Eulerian model was used for validation, while a simplified 2D model was applied to simulate conditions at a 50-m depth. The models integrate bubble dynamics, gas holdup, CO 2 dissolution, dissolved species transport, and seawater acidification into a unified CFD framework. Mass transfer was calculated using the Hughmark correlation, and local seawater temperature and salinity were factored in to determine dissociation behavior and the relevant Henry’s Law constant. To confirm the 3D model’s accuracy, results were checked against two experimental datasets: the QICS field study and the Hauser Tank experiments. The team also modeled a hypothetical release scenario at the High Island 10L site and compared the results with earlier published work. The results show that at a depth of 50 m, the surrounding water column can completely absorb a CO 2 release at a rate of 35 kg/s, since the gas dissolves into the seawater as it rises toward the surface. Beyond confirming this mitigation capacity, the simulations shed light on how a leak would actually unfold in the environment, including the shape and movement of the rising bubble plume, how much CO 2 dissolves along the way, and the resulting shifts in seawater pH and pCO 2 . Together, this provides a practical framework for assessing how CO 2 leaks could affect marine environments in the Gulf of Mexico.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2013

Measurements of charmless b-hadron decays at CDF; first evidence for the annihilation $B^0_s \to \pi^+ \pi^-$ decay mode.

In this thesis I report on several results obtained from an accurate analysis of the $Bhh$ data sample: Annihilation decay modes: using 6 fb$^{-1}$ of data collected by the CDF experiment, I obtained the first evidence of the $B_s^0\to \pi^+ \pi^-$ decay mode, with a significance of about 3.7$\sigma$, and I performed the most precise measurement of the Branching Ratio of the $B_d^0\to K^+ K^-$ decay mode at the time, corresponding to an excess of 2$\sigma$ level. These measurements have important consequences as discussed in the previous sections: the comparison between these branching fractions and the theoretical predictions provides unique information for tuning the phenomenological models. Our results prompted intense theoretical work both using pQCD and QCDF approaches to accommodate the discrepancies between predictions and experimental results. CP violation: using the full CDF 10 fb$^{-1}$ of data collected by the CDF experiment, I found evidence of the \textbf{CP violation} of $B_s^0 \to K \pi$ decay mode, with a significance of 3$\sigma$, confirming the result recently obtained at LHCb. I also measured the CP asymmetries of the $B_d^0 \to K \pi$ decay modes, in agreement and of comparable resolution with the b-factories experiments, while the most precise measurement up to date is from LHCb. The precision reached in the Acp($B_s^0 \to K \pi$ ) measurement gives the possibility to perform interesting checks of the Standard Model predictions, explained in the previous sections. I also obtained the world's best measurements of the CP asymmetries in $\Lambda_b^0 \to p\pi$ and $\Lambda_b^0 \to pK$ decay modes. These measurements are unique and represent a step forward in the understanding of the \CP\ violation phenomena in the relatively unexplored $\Lambda_b^0$ sector.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Machine learning models for rat multigeneration reproductive toxicity prediction

Reproductive toxicity is one of the prominent endpoints in the risk assessment of environmental and industrial chemicals. Due to the complexity of the reproductive system, traditional reproductive toxicity testing in animals, especially guideline multigeneration reproductive toxicity studies, take a long time and are expensive. Therefore, machine learning, as a promising alternative approach, should be considered when evaluating the reproductive toxicity of chemicals. We curated rat multigeneration reproductive toxicity testing data of 275 chemicals from ToxRefDB (Toxicity Reference Database) and developed predictive models using seven machine learning algorithms (decision tree, decision forest, random forest, k-nearest neighbors, support vector machine, linear discriminant analysis, and logistic regression). A consensus model was built based on the seven individual models. An external validation set was curated from the COSMOS database and the literature. The performances of individual and consensus models were evaluated using 500 iterations of 5-fold cross-validations and the external validation data set. The balanced accuracy of the models ranged from 58% to 65% in the 5-fold cross-validations and 45%–61% in the external validations. Prediction confidence analysis was conducted to provide additional information for more appropriate applications of the developed models. The impact of our findings is in increasing confidence in machine learning models. We demonstrate the importance of using consensus models for harnessing the benefits of multiple machine learning models (i.e., using redundant systems to check validity of outcomes). While we continue to build upon the models to better characterize weak toxicants, there is current utility in saving resources by being able to screen out strong reproductive toxicants before investing in vivo testing. The modeling approach (machine learning models) is offered for assessing the rat multigeneration reproductive toxicity of chemicals. Our results suggest that machine learning may be a promising alternative approach to evaluate the potential reproductive toxicity of chemicals.

consensus model↗

DOE OSTI2019

A search for cosmological anisotropy using the Lyman alpha forest from SDSS quasar spectra

ABSTRACT The cosmological principle, the combined assumptions of cosmological isotropy and homogeneity, underpins the standard model of big bang cosmology with which we interpret astronomical observations. A new test of isotropy over the redshift range 2 < z < 4 and across large angular scales on the sky is presented. We use the cosmological distribution of neutral hydrogen, as probed by the Ly α forest seen towards distant quasars. The Sloan Digital Sky Survey provides the largest data set of quasar spectra available to date. We use combined information from Data Releases 12 and 14 to select a sample of 142 661 quasars most suitable for this purpose. The scales covered by the data extend beyond post-inflation causality scales, thus probing initial conditions in the early Universe. We identify significant spatially correlated systematic effects that can emulate cosmological anisotropy. Once these systematics have been accounted for, the data are found to be consistent with isotropy, providing an important independent check on the standard model, consistent with results from cosmic microwave background data.

Zavarygin, Evgeny O.↗

DOE OSTI2024

The Atacama Cosmology Telescope: Cosmology from Cross-correlations of unWISE Galaxies and ACT DR6 CMB Lensing

We present tomographic measurements of structure growth using cross-correlations of Atacama Cosmology Telescope (ACT) DR6 and Planck cosmic microwave background (CMB) lensing maps with the unWISE Blue and Green galaxy samples, which span the redshift ranges 0.2 ≲ z ≲ 1.1 and 0.3 ≲ z ≲ 1.8, respectively. We improve on prior unWISE cross-correlations not just by making use of the new, high-precision ACT DR6 lensing maps, but also by including additional spectroscopic data for redshift calibration and by analyzing our measurements with a more flexible theoretical model. We determine the amplitude of matter fluctuations at low redshifts (z ≃ 0.2–1.6), finding S 8 ≡ σ 8 (Ω m /0.3) 0.5 = 0.813 ± 0.021 using the ACT cross-correlation alone and S 8 = 0.810 ± 0.015 with a combination of Planck and ACT cross-correlations; these measurements are fully consistent with the predictions from primary CMB measurements assuming standard structure growth. The addition of baryon acoustic oscillation data breaks the degeneracy between σ 8 and Ω m , allowing us to measure σ 8 = 0.813 ± 0.020 from the cross-correlation of unWISE with ACT and σ 8 = 0.813 ± 0.015 from the combination of cross-correlations with ACT and Planck. These results also agree with the expectations from primary CMB extrapolations in ΛCDM cosmology; the consistency of σ 8 derived from our two redshift samples at z ∼ 0.6 and 1.1 provides a further check of our cosmological model. Our results suggest that structure formation on linear scales is well described by ΛCDM even down to low redshifts z ≲ 1.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Low-density laboratory spectra near the He ii λ 304 line

Aims. Our aim is to interpret the EUV spectra of the solar corona, one hopes for laboratory data of specific chemical elements obtained under coronal conditions. Methods. EUV spectra of He, C, N, O, F, Ne, S, Ar, Fe, and Ni in a 40 Å wide wavelength interval near λ304 were excited in an electron beam ion trap. Results. We observe some two hundred lines about half of which are not yet identified and included in spectral models. Conclusions. Lastly, our data provide a check on the atomic data bases underlying the spectral models that are used to interpret solar corona data. However, a multitude of mostly weak additional lines taken together represent a flux that is comparable to that of various primary lines.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

High-fidelity multiphysics load following and accidental transient modeling of microreactors using NEAMS tools: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The feasibility of modeling microreactors using high-fidelity models with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tools is investigated in this report. Three overarching questions guided this research: can NEAMS tools readily be applied for high-fidelity multiphysics modeling of different types of transients in microreactor designs; how accurate are the results obtained; and are improvements needed in accuracy or user experience of NEAMS tools, especially considering newly developed capabilities? This work builds upon FY-2022 work, and two microreactor concepts considering heat pipe (HP-MR) and gas-cooled (GC-MR) technologies were further analyzed using high-fidelity multiphysics simulations. The NEAMS tools considered and coupled within the MultiApp environment are Griffin for neutronics, BISON for thermo-mechanics, Sockeye for heat pipe modeling (in HP-MR), SAM for 1D Fluid – 3D solid modeling of coolant channels and system modeling of balance of plant components (in GC-MR), and the SWIFT code for hydrogen redistribution in hydride moderator. The Heat Pipe MicroReactor (HP-MR) concept was further analyzed in FY-2023 to demonstrate the stochastic TRISO failure modeling capability in BISON to check operational limits of the TRISO fuel. A new full-core Gas-Cooled MicroReactor (GC-MR) model was developed based on the initial assembly-model used in Y-2022 and used for steady-state and accidental depressurization transient simulations. Accuracy of the simulations performed was assessed through 1) verification analyses completed on the different physics with code-to-code comparison, and 2) validation of the multiphysics simulations based on modeling of the Kilopower Reactor Using Stirling Technology (KRUSTY) experiment. In FY-2023, the mesh and model of KRUSTY was updated to closely match publicly available data, and the neutronic model was verified and validated against experimental control rod worth measurements. The multiphysics model of KRUSTY was developed and used for steady-state analysis and for modeling reactivity insertion transient. The calculated power increase and stabilization agrees well with experimental data following adjustment in fuel thermal expansion coefficient. As an important component of this project, the ANL team gathered experience with a wide range of NEAMS tools: the MOOSE Mesh System, Griffin, BISON, SWIFT, Sockeye, SAM, Workbench, and the MOOSE MultiApp System, and provided assessment of new capabilities. Noteworthy are the user assessment of the “vapor-only” flow model in Sockeye and development of a multiphysics startup transient in HP-MR unit cell for use as tutorial in Sockeye. The full-core GC-MR model was used for assessment of SAM for balance of plant modeling and for demonstrating the SWIFT code capability for hydrogen redistribution modeling in multiphysics transient analyses. In this process, several bugs/issues were identified and reported to developers. Finally, the assembly GC-MR model developed in FY-2022 coupling Griffin, BISON and SAM through flow blockage and rod ejection transients was published to the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB). The Heat Pipe MicroReactor (HP-MR) concept high-fidelity multiphysics coupling of Griffin/BISON/Sockeye in load-following and heat pipe failure transients was also published on the VTB. Those submissions are enabling thorough review of these models as well as wide distribution to industry, regulator, and university users. In this analysis, several new research questions were uncovered, and follow-up analyses are recommended to further improve some models, consider additional transients, and continue development of VTB models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2013

Radiative charge transfer in He{sup +}+ H{sub 2} collisions in the milli- to nano-electron-volt range: A theoretical study within state-to-state and optical potential approaches

The paper presents a theoretical study of the low-energy dynamics of the radiative charge transfer (RCT) reaction He{sup +}({sup 2}S)+H{sub 2}(X{sup 1}{Sigma}{sub g}{sup +}){yields}He({sup 1}S)+H{sub 2}{sup +}(X{sup 2}{Sigma}{sub g}{sup +})+h{nu} extending our previous studies on radiative association of HeH{sub 2}{sup +} [F. Mrugala, V. Spirko, and W. P. Kraemer, J. Chem. Phys. 118, 10547 (2003); F. Mrugala and W. P. Kraemer, ibid. 122, 224321 (2005)]. The calculations account for the vibrational and rotational motions of the H{sub 2}/H{sub 2}{sup +} diatomics and for the atom-diatom complex formation in the reactant and the product channels of the RCT reaction. Continuum states of He{sup +}+ H{sub 2}(v= 0, j= 0) in the collision energy range {approx}10{sup -7}-18.6 meV and all quasi-bound states of the He{sup +}- H{sub 2}(para;v= 0) complex formed in this range are taken into account. Close-coupling calculations are performed to determine rates of radiative transitions from these states to the continuum and quasi-bound states of the He +H{sub 2}{sup +} system in the energy range extending up to {approx}0.16 eV above the opening of the HeH{sup +}+ H arrangement channel. From the detailed state-to-state calculated characteristics global functions of the RCT reaction, such as cross-section {sigma}(E), emission intensity I({nu}, T), and rate constant k(T) are derived, and are presented together with their counterparts for the radiative association (RA) reaction He{sup +}({sup 2}S) +H{sub 2}(X{sup 1}{Sigma}{sub g}{sup +}){yields} HeH{sub 2}{sup +}(X{sup 2}A{sup Prime })+h{nu}. The rate constant k{sup RCT} is approximately 20 times larger than k{sup RA} at the considered temperatures, 0.1 {mu}K-50 K. Formation of rotational Feshbach resonances in the reactant channel plays an important role in both reactions. Transitions mediated by these resonances contribute more than 70% to the respective rates. An extension of the one-dimensional optical potential model is developed to allow inclusion of all three vibrational modes in the atom-diatom system. This three-dimensional optical potential model is used to check to which extent the state-to-state RCT rate constant is influenced by the possibility to access ground state continuum levels well above the opening of the HeH{sup +}+ H arrangement channel. The results indicate that these transitions contribute about 30% to the 'true' rate constant k{sup RCT} whereas their impact on the populations of the vibration-rotational states of the product H{sub 2}{sup +} ion is only minor. Present theoretical rate constant functions k{sup RCT}(T) obtained at different approximation levels are compared to experimental data: 1-1.1 Multiplication-Sign 10{sup -14} s{sup -1} cm{sup 3} at T= 15-35 K and {approx}7.5 Multiplication-Sign 10{sup -15} s{sup -1} cm{sup 3} at 40 K [M. M. Schauer, S. R. Jefferts, S. E. Barlow, and G. H. Dunn, J. Chem. Phys. 91, 4593 (1989)]. The most reliable theoretical values of k{sup RCT}, obtained by combining results from the state-to-state and the optical potential calculations, are between 2.5 and 3.5 times larger than these experimental numbers. Possible sources for discrepancies are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

Design study of the geometry of the blanking tool to predict the burr formation of Zircaloy-4 sheet

In this work, we investigated factors that influence burr formation for zircaloy-4 sheet used for spacer grids of nuclear fuel roads. Factors we considered are geometric factors of punch. We changed clearance and velocity in order to consider the failure parameters, and we changed shearing angle and corner radius of L-shaped punch in order to consider geometric factors of punch. First, we carried out blanking test with failure parameter of GTN model using L-shaped punch. The tendency of failure parameters and geometric factors that affect burr formation by analyzing sheared edges is investigated. Consequently, geometric factor's influencing on the burr formation is also high as failure parameters. Then, the sheared edges and burr formation with failure parameters and geometric factors is investigated using FE analysis model. As a result of analyzing sheared edges with the variables, we checked geometric factors more affect burr formation than failure parameters. To check the reliability of the FE model, the blanking force and the sheared edges obtained from experiments are compared with the computations considering heat transfer.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

A novel method for the measurement of the von Neumann spike in detonating high explosives

We present detonation wave profiles measured in T2 (97 wt. % TATB) and TX1 (52 wt. % TATB and 45 wt. % HMX) high explosives. The experiments consisted in initiating a detonation wave in a 15 mm diameter cylinder of explosive using an explosive wire detonator and an explosive booster. Free surface velocity wave profiles were measured at the explosive/air interface using a Photon Doppler Velocimetry system. We demonstrate that a comparison of these free surface wave profiles with those measured at explosive/window interfaces in similar conditions allows to bracket the von Neumann spike in a narrow range. For T2, our measurements show that the spike pressure lies between 35.9 and 40.1 GPa, whereas for TX1, it lies between 42.3 and 47.0 GPa. The numerical simulations performed in support to these measurements show that they can be used to calibrate reactive burn models and also to check the accuracy of the detonation products equation of state at low pressure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

The phase space and stellar populations of cluster galaxies at z ∼ 1: simultaneous constraints on the location and timescale of satellite quenching

We investigate the velocity versus position phase space of z ∼ 1 cluster galaxies using a set of 424 spectroscopic redshifts in nine clusters drawn from the GCLASS survey. Dividing the galaxy population into three categories, that is, quiescent, star-forming, and poststarburst, we find that these populations have distinct distributions in phase space. Most striking are the poststarburst galaxies, which are commonly found at small clustercentric radii with high clustercentric velocities, and appear to trace a coherent 'ring' in phase space. Using several zoom simulations of clusters, we show that the coherent distribution of the poststarbursts can be reasonably well reproduced using a simple quenching scenario. Specifically, the phase space is best reproduced if these galaxies are quenched with a rapid timescale (0.1 <τ {sub Q} < 0.5 Gyr) after they make their first passage of R ∼ 0.5 R {sub 200}, a process that takes a total time of ∼1 Gyr after first infall. The poststarburst phase space is not well reproduced using long quenching timescales (τ {sub Q} > 0.5 Gyr) or by quenching galaxies at larger radii (R ∼ R {sub 200}). We compare this quenching timescale to the timescale implied by the stellar populations of the poststarburst galaxies and find that the poststarburst spectra are well-fit by a rapid quenching (τ {sub Q} = 0.4{sub −0.4}{sup +0.3} Gyr) of a typical star-forming galaxy. The similarity between the quenching timescales derived from these independent indicators is a strong consistency check of the quenching model. Given that the model implies satellite quenching is rapid and occurs well within R {sub 200}, this would suggest that ram-pressure stripping of either the hot or cold gas component of galaxies are the most plausible candidates for the physical mechanism. The high cold gas consumption rates at z ∼ 1 make it difficult to determine whether hot or cold gas stripping is dominant; however, measurements of the redshift evolution of the satellite quenching timescale and location may be capable of distinguishing between the two.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Constraints on dark matter annihilations from diffuse gamma-ray emission in the Galaxy

Recent advances in γ-ray cosmic ray, infrared and radio astronomy have allowed us to develop a significantly better understanding of the galactic medium properties in the last few years. In this work using the DRAGON code, that numerically solves the CR propagation equation and calculating γ-ray emissivities in a 2-dimensional grid enclosing the Galaxy, we study in a self consistent manner models for the galactic diffuse γ-ray emission. Our models are cross-checked to both the available CR and γ-ray data. We address the extend to which dark matter annihilations in the Galaxy can contribute to the diffuse γ-ray flux towards different directions on the sky. Moreover we discuss the impact that astrophysical uncertainties of non DM nature, have on the derived γ-ray limits. Such uncertainties are related to the diffusion properties on the Galaxy, the interstellar gas and the interstellar radiation field energy densities. Light ∼ 10 GeV dark matter annihilating dominantly to hadrons is more strongly constrained by γ-ray observations towards the inner parts of the Galaxy and influenced the most by assumptions of the gas distribution; while TeV scale DM annihilating dominantly to leptons has its tightest constraints from observations towards the galactic center avoiding the galactic disk plane, with the main astrophysical uncertainty being the radiation field energy density. In addition, we present a method of deriving constraints on the dark matter distribution profile from the diffuse γ-ray spectra. These results critically depend on the assumed mass of the dark matter particles and the type of its end annihilation products.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

A gamma-ray pulsar timing array constrains the nanohertz gravitational wave background

We report after large galaxies merge, their central supermassive black holes are expected to form binary systems. Their orbital motion should generate a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background use pulsar timing arrays, which perform long-term monitoring of millisecond pulsars at radio wavelengths. We used 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95% credible limit on the GWB characteristic strain of 1.0 × 10 -14 at a frequency of 1 year –1 . The sensitivity is expected to scale with t obs , the observing time span, as t$^{-13/6}_{obs}$. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Contrasting galaxy formation from quantum wave dark matter, ψDM, with ΛCDM, using Planck and Hubble data

The newly established luminosity functions (LFs) of high-z galaxies at 4 ≲ z ≲ 10 can provide a stringent check on dark matter models that aim to explain the core properties of dwarf galaxies. The cores of dwarf spheroidal galaxies are understood to be too large to be accounted for by free streaming of warm dark matter without overly suppressing the formation of such galaxies. Here we demonstrate with cosmological simulations that wave dark matter, ψDM, appropriate for light bosons such as axions, does not suffer from this problem, given a boson mass of m{sub ψ}⩾1.2×10{sup −22} eV (2σ). In this case, the halo mass function is suppressed below ∼10{sup 10} M{sub ⊙} at a level that is consistent with the high-z LFs, while simultaneously generating the kiloparsec-scale cores in dwarf galaxies arising from the solitonic ground state in ψDM. We demonstrate that the reionization history in this scenario is consistent with the Thomson optical depth recently reported by Planck, assuming a reasonable ionizing photon production rate. We predict that the LF should turn over slowly around an intrinsic ultraviolet luminosity of M{sub UV}≳−16 at z≳4. We also show that for galaxies magnified >10× in the Hubble Frontier Fields, ψDM predicts an order of magnitude fewer detections than cold dark matter at z≳10 down to M{sub UV}∼−15, allowing us to distinguish between these very different interpretations for the observed coldness of dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗