Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical relativity”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Influence of temperature on accessible pyrolysis pathways of homopolymerized bisphenol A/F epoxies and copolymers

Understanding the thermal pyrolysis of epoxies and their copolymers is important for identifying structural changes resultant from thermal transients, enabling identification of failure modes of high-performance composite materials. This work expands our understanding of the thermal pyrolysis of cured epoxies and the role of temperature and composition on products, pathways, and relative rates. Numerous researchers have explored the pyrolysis of bisphenol A (BPA) epoxy. Significantly fewer have studied bisphenol F (BPF), and copolymers of BPA and BPF have been neglected. In this work, a pyrolysis gas chromatography mass spectrometer (PY-GC-MS) was used to investigate the degradation mechanism of homopolymerized BPA, BPF epoxies and their copolymers. Additionally, for polymer identification, pyrolysis >450 °C resulted in the most extensive fragmentation and is useful for material identification, though lower temperatures show different degradation product profiles that provide mechanistic insight into thermal degradation pathways. Temperature greatly influences the accessible pyrolysis pathways of BPA, revealing dual mechanisms of formation for products p-isopropylphenol and p-isopropenylphenol. For BPF at low pyrolysis temperature the p,p-bisphenol F isomer is produced at significantly lower relative abundance compared to the abundance at higher temperatures tested, but production of the other two isomers changes little with respect to temperature. This suggests the epoxy components consisting of the p,p-bisphenol F isomer have higher thermal stability. Overall, the copolymer epoxies were found to have similar degradation products in stoichiometric distributions. The major exception was the formation of the p,p-bisphenol F isomer, which shows evidence of thermal stabilizing effects from addition of BPA epoxy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Locally purified maximally mixed states at scale: Entanglement pruning and symmetries

Locally Purified Density Operators (LPDOs) are state-of-the-art tensor network ansatze candidates that efficiently represent mixed quantum states at scale. However, given their non-uniqueness, their representational complexity is generally sub-optimal in practical computations. Here, in this work we perform a comprehensive numerical and analytical analysis and resolve this issue in the experimentally relevant limit where noise depolarizes the density operator into a maximally mixed state. To resolve the sub-optimality issue, we analyze two numerical tools, one analytic method, and detail the relations between them. The numerical tools used are fidelity-preserving truncations and isometric gauge transformations leveraging Riemannian optimizations over entropic objective functions. In addition, by invoking the injectivity and symmetry constraints of the maximally mixed LPDO, we also present analytical closed-form expressions for the disentangler and discuss their relation to numerical optimizers. Further, away from the maximally mixed state, our simulations highlight how the truncation threshold smoothly interpolate, as a function of depolarization, between established matrix product results and our new results. Our work shows how, by minimizing the resources required to represent key states of practical interest in experiment, the efficiency of tensor network algorithms can be substantially increased. This paves the path for uncovering tensor network’s fundamental scalability limits and latent potential in representing the wide locus of mixed quantum states that are accessible on near-term quantum devices.

Gangapuram, Amit Jamadagni [Oak Ridge National Lab↗
DOE OSTI2022

Reflection and transmission of electromagnetic pulses at a planar dielectric interface -- theory and quantum lattice simulations

There is considerable interest in the application of quantum information science to advance computations in plasma physics. A particular point of curiosity is whether it is possible to take advantage of quantum computers to speed up numerical simulations relative to conventional computers. Many of the topics in fusion plasma physics are classical in nature. In order to implement them on quantum computers it will require couching a classical problem in the language of quantum mechanics. Electromagnetic waves are routinely used in fusion experiments to heat a plasma or to generate currents in the plasma. The propagation of electromagnetic waves is described by Maxwell equations with an appropriate description of the plasma as a dielectric medium. Before advancing to the tensor dielectric of a magnetized plasma, this paper considers electromagnetic wave propagation in a one-dimensional inhomogeneous scalar dielectric. The classic theory of scattering of plane electromagnetic waves at a planar interface, separating two different dielectric media, leads to Fresnel equations for reflection and transmission coefficients. In contrast to plane waves, this paper is on the reflection and transmission of a spatially confined electromagnetic pulse. Following an analytical formulation for the scattering of a Gaussian pulse, it is deduced that the maximum transmission coefficient for a pulse is $\sqrt{n_2/n_1}$ times that for a plane wave; the incident and transmitted pulses propagating in dielectric media with refractive indices $n_1$ and $n_2$, respectively. The analytical theory is complemented by numerical simulations using a quantum lattice algorithm for Maxwell equations. The algorithm, based on the Riemann-Silberstein-Weber representation of the electromagnetic fields and expressed in term of qubits, is an interleaved sequence of entangling operators at each lattice site and unitary streaming operators which transmit information from one site to an adjacent lattice site. Besides substantiating results from the theory for Gaussian pulses, numerical simulations show their validity for non-Gaussian pulses. Apart from their time-asymptotic forms, the simulations display an interplay between the incident, reflected, and transmitted pulses in the vicinity of the transition region between two dielectric media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Neutron stars and the dense matter equation of state

The past two decades have witnessed tremendous progress in understanding the properties of neutron stars, their maximum mass and radii, and the properties of the dense matter in their cores, made possible by electromagnetic observations of neutron stars and the detection of gravitational waves from their mergers. These observations have provided novel constraints on neutron-star structure that are intimately related to the properties of dense neutron-rich matter described by the nuclear equation of state. Nevertheless, constraining the equation of state over the wide range of densities probed by astrophysical observations is still challenging, as the physics involved is broad and the system spans many orders of magnitude in density. Here, theoretical approaches to calculate and model the neutron-star equation of state in various regimes of densities are reviewed, and the related consequent properties of neutron stars are discussed. How the equation of state at low densities can be calculated from nuclear interactions that are constrained and benchmarked by nuclear experiments is described. Neutron-star observations, with a particular emphasis on information provided by gravitational-wave signals and electromagnetic observations, are reviewed. Finally, future challenges and opportunities in the field are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Robustness of neutron star merger simulations to changes in neutrino transport and neutrino-matter interactions

Binary neutron star mergers play an important role in nuclear astrophysics: Their gravitational wave and electromagnetic signals carry information about the equation of state of cold matter above nuclear saturation density, and they may be one of the main sources of r-process elements in the Universe. Neutrino-matter interactions during and after merger impact the properties of these electromagnetic signals and the relative abundances of the produced r-process elements. Existing merger simulations are, however, limited in their ability to realistically model neutrino transport and neutrino-matter interactions. Here, we perform a comparison of the impact of the use of state-of-the art two-moment or Monte Carlo transport schemes on the outcome of merger simulations, for a single binary neutron star system with a short-lived neutron star remnant [(5–10) ms]. We also investigate the use of different reaction rates in the simulations. While the best transport schemes generally agree well on the qualitative impact of neutrinos on the system, differences in the behavior of the high-density regions can significantly impact the collapse time and the properties of the hot tidal arms in this metastable merger remnant. The chosen interaction rates and transport algorithm as well as recent improvements by Radice et al. to the two-moment algorithms can all contribute to changes at the (10–30)% level in the global properties of the merger remnant and outflows. Furthermore, the limitations of previous moment schemes fixed by Radice et al. also appear sufficient to explain the large difference that we observed in the production of heavy-lepton neutrinos in a previous comparison of Monte Carlo and moment schemes in the context of a low-mass binary neutron star system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Characterization, Structure, and Reactivity of Hydroxyl Groups on Metal‐Oxide Cluster Nodes of Metal–Organic Frameworks: Structural Diversity and Keys to Reactivity and Catalysis

Among the most stable metal–organic frameworks (MOFs) are those incorporating nodes that are metal oxide clusters with frames such as Zr 6 O 8 . This review is a summary of the structure, bonding, and reactivity of MOF node hydroxyl groups, emphasizing those bonded to nodes containing aluminum and zirconium ions. Hydroxyl groups are often present on these nodes, sometimes balancing the charges of the metal ions. They arise during MOF syntheses in aqueous media or in post-synthesis treatments. They are identified with infrared and 1 H nuclear magnetic resonance spectroscopies and characterized by their reactivities with polar compounds such as alcohols. Terminal OH, paired µ 2 -OH, and aqua groups on nodes are catalytic sites in numerous reactions. Relatively unreactive hydroxyl groups (such as isolated µ 2 -OH groups) may replace reactive groups and inhibit catalysis; some node hydroxyl groups (e.g., µ 3 -OH) are mere spectators in catalysis. There are similarities between MOF node hydroxyl groups and those on the surfaces of bulk metal oxides, zeolites, and enzymes, but the comparisons are mostly inexact, and much remains to be understood about MOF node hydroxyl group chemistry. It is posited that understanding and controlling this chemistry will lead to tailored MOFs and improved adsorbents and catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Advances in Machine and Deep Learning for Modeling and Real-time Detection of Multi-Messenger Sources

We live in momentous times. The science community is empowered with an arsenal of cosmic messengers to study the universe in unprecedented detail. Gravitational waves, electromagnetic waves, neutrinos, and cosmic rays cover a wide range of wavelengths and timescales. Combining and processing these datasets that vary in volume, speed, and dimensionality requires new modes of instrument coordination, funding, and international collaboration with a specialized human and technological infrastructure. In tandem with the advent of large-scale scientific facilities, the last decade has experienced an unprecedented transformation in computing and signal-processing algorithms. The combination of graphics processing units, deep learning, and the availability of open source, high-quality datasets has powered the rise of artificial intelligence. This digital revolution now powers a multibillion dollar industry, with far-reaching implications in technology and society. In this chapter, we describe pioneering efforts to adapt artificial intelligence algorithms to address computational grand challenges in multi-messenger astrophysics. We review the rapid evolution of these disruptive algorithms, from the first class of algorithms introduced in early 2017 to the sophisticated algorithms that now incorporate domain expertise in their architectural design and optimization schemes. We discuss the importance of scientific visualization and extreme-scale computing in reducing time-to-insight and obtaining new knowledge from the interplay between models and data.

Artificial Intelligence↗
DOE OSTI2021

blocks_3d: software for general 3d conformal blocks

We introduce the software blocks_3d for computing four-point conformal blocks of operators with arbitrary Lorentz representations in 3d CFTs. It uses Zamolodchikov-like recursion relations to numerically compute derivatives of blocks around a crossing-symmetric configuration. It is implemented as a heavily optimized, multi-threaded, C++ application. We give performance benchmarks for correlators containing scalars, fermions, and stress tensors. As an example application, we recompute bootstrap bounds on four-point functions of fermions and study whether a previously observed sharp jump can be explained using the “fake primary” effect. We conclude that the fake primary effect cannot fully explain the jump and the possible existence of a “dead-end” CFT near the jump merits further study.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Overview of high-density QCD studies with the CMS experiment at the LHC

We review key measurements performed by CMS in the context of its heavy ion physics program, using event samples collected in 2010-2018 with several collision systems and energies. These studies provide detailed macroscopic and microscopic probes of the quark-gluon plasma (QGP) created at the LHC energies, a medium characterized by the highest temperature and smallest baryon-chemical potential ever reached in the laboratory. Numerous observables related to high-density quantum chromodynamics (QCD) were studied, leading to some of the most impactful and qualitatively novel results in the 40-year history of the field. Using a dedicated high-multiplicity trigger in the first pp run, CMS discovered that small collision systems can exhibit signs of collectivity, a generic phenomenon with significant implications and presently understood to affect essentially all soft physics processes. This observation opened new paths to understand how fluidity and plasma properties emerge in QCD matter as a function of system size. Measurements of jet quenching have reached a completely new level of detail by directly assessing, for the first time, the medium modification of parton showers, as opposed to simply observing leading hadrons or di-hadrons. The first fully reconstructed beauty hadron and heavy-flavor jet nuclear modifications were also measured. The large size of the event samples, the precision of the measurements, and the extension of the probed kinematical phase space, allowed many other hard probes of the QGP medium to be explored in detail, leading to multiple groundbreaking findings. In particular, the seminal measurements of bottomonium suppression patterns answer fundamental questions that have been actively pursued, both theoretically and experimentally, by the community since the mid-1980s. We conclude by outlining the opportunities offered by the continuation of this physics program at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

State-of-the-art short-term electricity market operation with solar generation: A review

With the accelerating penetration of solar energy in energy systems, market operations concerning the solar-associated intermittency are widely discussed. How to correctly model the solar generation in the market and solve the uncertainty-based operation problems call for solutions. Unlike other renewable resources, solar power can be more easily applied to the demand side through behind-the-meter installations in a distributed manner with significant future extensibility. The corresponding operational rules in the markets are largely being revised. Over the past two decades, with the development of deregulated power markets, numerous studies related to uncertainty-based market operations tend to transform the market from a deterministic and centralized framework to a stochastic and decentralized one. As the relevant literature is most recent, this paper reviews their contributions to electricity market operations with solar integration and underlines the research directions of this topic and ongoing industry practices. A thorough review of the electricity markets worldwide with solar energy is provided. A variety of proposed mathematical solutions to the problem are also addressed.

14 SOLAR ENERGY↗
DOE OSTI2024

Haasttse‐baad Tessera Ring Complex: A Valhalla‐Type Impact Structure on Venus?

Abstract Venus preserves ∼1,000 impact craters, yet to date no impact basins larger than 300 km in diameter—common in the oldest terrains on Mercury, Mars and the Moon—are recognized on Venus. The tessera terrain is Venus' oldest recognized terrain. We describe a ∼1,500 km‐diameter concentric ring‐graben complex preserved on Haasttse‐baad Tessera, Venus that we identify as the Haasttse‐baad Tessera Ring Complex (HTRC). Based on geologic relations and numerical modeling, we propose that the HTRC may represent a Valhalla‐type multiring impact basin formed late during the evolution of its host ribbon‐tessera terrain (rtt). Formation of Valhalla‐type impact basins could involve a unique three‐layer target rheology with a thin elastic layer above a low viscosity layer above a deep strong layer. This multi‐layer rheological sandwich is consistent with crustal rheology previously proposed for the formation of Venus' rtt. If the HTRC is a Valhalla‐type impact basin, it would be Venus' oldest, and currently largest, impact structure, providing a rare window into Venus' ancient past and with implications for early crustal processes on Venus.

58 GEOSCIENCES↗
DOE OSTI2025

Response of Atmospheric River Width and Intensity to Aquaplanet Warming: A Detection Algorithm‐ and Background Moisture‐Independent Approach

The width of an atmospheric river (AR) is an important parameter when evaluating its impact. Although previous research suggests ARs will widen with global warming, a precise response has been muddled by the large sensitivity of width to a diverse set of AR detection techniques (ARDTs). Here, we propose a methodology that removes the influence of the ARDT by modeling AR‐integrated vapor transport (IVT) profiles as idealized exponential curves with free parameters given by background IVT, intensity above background IVT, and profile width. Notably, our definition for AR profile width does not include any explicit numerical thresholds, relative or absolute, for IVT. We apply our approach to a series of idealized aquaplanet experiments, first with a baseline sea surface temperature (SST) distribution, and then with +2K, +4K, and +6K uniform warming, so as to determine the contributions of each free parameter to AR width. We also apply our approach to high‐resolution atmosphere‐only models forced with SSTs modified to preserve historical variability but following projected warming over 2016–2050. Our results show that contributions to impacts‐relevant AR widening comes primarily from enhancements in background IVT and AR intensity, as opposed to from dynamic width changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Machine learning assisted derivation of minimal low-energy models for metallic magnets

Abstract We consider the problem of extracting a low-energy spin Hamiltonian from a triangular Kondo Lattice Model (KLM). The non-analytic dependence of the effective spin-spin interactions on the Kondo exchange excludes the use of perturbation theory beyond the second order. We then introduce a Machine Learning (ML) assisted protocol to extract effective two- and four-spin interactions. The resulting spin model reproduces the phase diagram of the original KLM as a function of magnetic field and single-ion anisotropy and reveals the effective four-spin interactions that stabilize the field-induced skyrmion crystal phase. Moreover, this model enables the computation of static and dynamical properties with a much lower numerical cost relative to the original KLM. A comparison of the dynamical spin structure factor in the fully polarized phase computed with both models reveals a good agreement for the magnon dispersion even though this information was not included in the training data set.

Sharma, Vikram (ORCID:0000000189105519)↗
DOE OSTI2025

Black hole-neutron star binaries near neutron star disruption limit in the mass regime of event GW230529

In May 2023, the LIGO-Virgo-KAGRA (LVK) Collaboration reported the likely black hole-neutron star (BHNS) merger GW230529_181500. The signal was observed with high significance in only one detector, limiting constraints on the black hole spin and motivating our study of disruption in this mass regime. That event is expected to be the merger of a 2.5–4.5 M ⊙ primary with a secondary compact object of mass between 1.2–2.0 M ⊙ . This makes it the first BHNS merger with a significant potential for the production of electromagnetic (EM) counterparts, and provides further evidence for compact objects existing within the suspected lower mass gap. To produce post-merger EM transients, the component of the black hole spin aligned with the orbital angular momentum must be sufficiently high, allowing the neutron star to be tidally disrupted. The disrupting BHNS binary may then eject a few percent of a solar mass of matter, leading to an observable kilonova driven by radioactive decays in ejecta, and/or a compact-binary gamma-ray burst (cbGRB) resulting from the formation of an accretion disk and relativistic jet. Determining which mergers lead to disruption of the neutron star is necessary to predict the prevalence of EM signals from BHNS mergers, yet most BHNS simulations so far have been performed far from the minimum spin required for tidal disruption. Here, we use the Spectral Einstein Code to explore the behavior of BHNS mergers in a mass range consistent with GW230529_181500 close to that critical spin, and compare our results against the mass remnant model currently used by the LVK Collaboration to predict the probability of tidal disruption. Furthermore, our numerical results reveal the emergence of non-zero accretion disks even below the predicted NS disruption limit, of low mass but capable of powering cbGRBs. Our results also demonstrate that the remnant mass model underpredicts the disk mass for the DD2 equation of state, while they are within expected modeling errors for SFHo. The disruption limit itself, however, is not found to significantly differ from the predictions of the analytical model, unless remnant masses M rem ≲ 0.001 M ⊙ prove interesting observationally. In all of our simulations, any kilonova signal would be dim and most likely dominated by post-merger disk outflows.

GW230529_181500↗
DOE OSTI2026

Investigating the impact of higher-order phase transitions in binary neutron-star mergers

In this paper we investigate quark deconfinement in neutrons stars and their mergers, focusing on the effects of higher orders for the phase transition between hadronic and quark matter. The different descriptions we use to describe matter microscopically contain varying particle degrees of freedom, including nucleons, hyperons, Delta baryons, and light and strange quarks. We use tabulated equations of state from the CompOSE database in which the quark deconfinement phase transition is described as being first order, and then smooth it out by introducing a percolation, replacing the single first-order phase transition with two transitions of second or third order. We then perform binary neutron-star merger simulations using these new equations of state, focusing on groups of binaries with the same single-star mass, radius, and tidal deformability, but different equations of state. Here, we go on to discuss differences in their evolution, and the ramifications for interpreting future gravitational wave observations and the potential to learn about dense matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Isospin-symmetry implications for nuclear two-body distributions and short-range correlations

Here, we study the implications of isospin symmetry in nuclear systems on two-body distributions and derive relations between different densities or momentum distributions of a given nucleus or of different nuclei in the same isospin multiplet. A connection between neutrinoless double beta decay transition of mirror nuclei and two-body densities is also obtained. The relations are numerically verified by means of quantum Monte Carlo calculations. We discuss the relevance of these results to the analysis of nuclear short-range correlations and also show that isospin symmetry allows us to extract information about the spectator nucleons when a short-range correlated pair is formed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗