Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained wave”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Improved time-delay lens modelling and H 0 inference with transient sources

ABSTRACT Strongly lensed explosive transients such as supernovae, gamma-ray bursts, fast radio bursts, and gravitational waves are very promising tools to determine the Hubble constant (H0) in the near future in addition to strongly lensed quasars. In this work, we show that the transient nature of the point source provides an advantage over quasars: The lensed host galaxy can be observed before or after the transient’s appearance. Therefore, the lens model can be derived from images free of contamination from bright point sources. We quantify this advantage by comparing the precision of a lens model obtained from the same lenses with and without point sources. Based on Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) observations with the same sets of lensing parameters, we simulate realistic mock data sets of 48 quasar lensing systems (i.e. adding AGN in the galaxy centre) and 48 galaxy–galaxy lensing systems (assuming the transient source is not visible but the time delay and image positions have been or will be measured). We then model the images and compare the inferences of the lens model parameters and H0. We find that the precision of the lens models (in terms of the deflector mass slope) is better by a factor of 4.1 for the sample without lensed point sources, resulting in an increase of H0 precision by a factor of 2.9. The opportunity to observe the lens systems without the transient point sources provides an additional advantage for time-delay cosmography over lensed quasars. It facilitates the determination of higher signal-to-noise stellar kinematics of the main deflector, and thus its mass density profile, which, in turn plays a key role in breaking the mass-sheet degeneracy and constraining H0.

Ding, Xuheng↗
DOE OSTI2024

Hydropower Flexibility and Environmental Tradeoffs Analysis

The importance of hydropower increases as the power grid evolves with the higher variable renewable contribution. As conventional thermal power plants are retired, the importance of hydropower contribution increases to balance the variability of solar and wind generation. However, reservoir water resources are constrained by multiple constraints, and variability of water inflow to the reservoirs creates limitations to dam water releases for power grid needs. Coordinating multiple tools, including water resources, ecological, and technical and economic power grid modeling, informs dam water releases. The case study, the Columbia River Basin multipurpose reservoir project, is operated for hydropower production and many other purposes considering the aquatic habitat of the river basin. Specifically, the river basin fish population is a vital element for the tribal community of the river basin. We integrated a production cost model, a water resource model, and decades of tribal knowledge to analyze the fish-friendly way of operating Columbia hydropower scheduling and grid impacts. We measure power grid impacts for various water resources planning scenarios in terms of total system operating cost, system reliability indicators, changes in wind and solar generation and curtailments, local marginal prices, and revenue for hydropower producers. The study results inform reservoir operating rules decisions from hydropower power producers, system operators, other water users, tribes, environmentalists, and other stakeholders.

Columbia River↗
DOE OSTI2023

Energy Frontier Research Centers: Center for the Computational Design of Functional Layered Materials (CCDM) August 1, 2014 - July 31, 2018; Center for Complex Materials from First Principles (CCM) August 1, 2018 - July 31, 2021 (Final Report)

The mission of the DOE Energy Frontier Research Centers CCDM (2014-2018) and CCM (2018-2021) was to theoretically develop, computationally apply, and experimentally validate electronic structure methods for all materials, with a focus on the complex materials, especially layered and two-dimensional materials, strongly-correlated materials, and liquid water. This was achieved by over 200 published journal articles authored by about 17 senior investigators from physics and chemistry and from theory, computation, and experiment, plus their collaborators. In particular, the Centers confirmed the predictive power of the SCAN (strongly constrained and appropriately normed) density functional, which was constructed to satisfy 17 known exact constraints and several appropriate norms. Without being fitted to real bonded systems, and at a modest computational cost, SCAN correctly predicted covalent, ionic, metallic, hydrogen, and van der Waals bonds in many challenging materials. SCAN gave an improved description of defects in semiconductors, surface properties of metals, seven phases of ice, liquid water, liquid and supercooled silicon, subtle structural distortions in ferroelectrics, formation energies and structural predictions for solids, and critical pressures for structural phase transitions. Perhaps most remarkably, SCAN correctly described some strongly-correlated materials that were previously believed to be beyond the reach of density-functional approximations. SCAN is the only density functional that correctly predicts the band gap closing under chemical doping of the cuprate high-temperature superconducting materials. SCAN also predicts a landscape of competing stripe and magnetic phases in the cuprates. For some materials with some codes, SCAN has convergence problems that are greatly reduced by the CCM-developed r 2 SCAN, without loss of accuracy or rigor. SCAN and r 2 SCAN still make some self-interaction error, which is greatly reduced by the CCDM/ CCM-developed local orbital scaling correction (LOSC). These Centers further proved that the fundamental energy gaps of a solid from an orbital energy difference and from total energy differences are the same for a large class of generalized Kohn-Sham (GKS) functionals, including SCAN and standard hybrid functionals, and that symmetry breaking arises when a dynamic density fluctuation drops to zero frequency. The Centers identified new mechanisms for catalysis in layered materials with ions intercalated between the layers, investigated charge density waves both in model systems and in real layered materials, studied changes of band gap with the number of layers, and explored topological ultrathin films, bent nanoribbons, and defects.

08 HYDROGEN↗
DOE OSTI2022

Role of hot electrons in shock ignition constrained by experiment at the National Ignition Facility

Shock ignition is a scheme for direct drive inertial confinement fusion that offers the potential for high gain with the current generation of laser facility; however, the benefits are thought to be dependent on the use of low adiabat implosions without laser–plasma instabilities reducing drive and generating hot electrons. A National Ignition Facility direct drive solid target experiment was used to calibrate a 3D Monte Carlo hot-electron model for 2D radiation-hydrodynamic simulations of a shock ignition implosion. The α = 2.5 adiabat implosion was calculated to suffer a 35% peak areal density decrease when the hot electron population with temperature T ℎ = 55 keV and energy E ℎ = 13 kJ was added to the simulation. Optimizing the pulse shape can recover ~1/3 of the peak areal density lost due to a change in shock timing. Here, despite the harmful impact of laser–plasma instabilities, the simulations indicate shock ignition as a viable method to improve performance and broaden the design space of near ignition high adiabat implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1999

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗
DOE OSTI2023

Inference of the optical depth to reionization τ from Planck CMB maps with convolutional neural networks

The optical depth to reionization, τ, is the least constrained parameter of the cosmological Λ cold dark matter (ΛCDM) model. To date, its most precise value is inferred from large-scale polarized cosmic microwave background (CMB) power spectra from the High Frequency Instrument (HFI) aboard the Planck satellite. These maps are known to contain significant contamination by residual non-Gaussian systematic effects, which are hard to model analytically. Therefore, robust constraints on τ are currently obtained through an empirical cross-spectrum likelihood built from simulations. In this paper, we present a likelihood-free inference of τ from polarized Planck HFI maps which, for the first time, is fully based on neural networks (NNs). NNs have the advantage of not requiring an analytical description of the data and can be trained on state-of-the-art simulations, combining the information from multiple channels. By using Gaussian sky simulations and Planck SRoll2 simulations, including CMB, noise, and residual instrumental systematic effects, we trained, tested, and validated NN models considering different setups. We inferred the value of τ directly from Stokes Q and U maps at ~4° pixel resolution, without computing angular power spectra. On Planck data, we obtained τ NN = 0.0579 ± 0.0082, which is compatible with current EE cross-spectrum results but with a ~30% larger uncertainty, which can be assigned to the inherent nonoptimality of our estimator and to the retraining procedure applied to avoid biases. While this paper does not improve on current cosmological constraints on τ, our analysis represents a first robust application of NN-based inference on real data, and highlights its potential as a promising tool for complementary analysis of near-future CMB experiments, also in view of the ongoing challenge to achieve the first detection of primordial gravitational waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Temperature-induced densification in compressed basaltic glass revealed by in-situ ultrasonic measurements

Abstract Acoustic velocities of a model basalt glass (64 mol% CaMgSi2O6 + 36 mol% CaAl2Si2O8) were measured along different pressure-temperature (P-T) paths. One set of experiments involved isothermal compression-decompression cycles, performed at temperatures of 300, 641, 823, and 1006 K and pressures up to 12.2 GPa. The other set of experiments involved constant-load heating-cooling cycles at temperatures up to 823 K and pressures up to 7.5 GPa. Both sets of experiments were performed in a multi-anvil apparatus using a synchrotron-based ultrasonic technique. Our results show that the glass compressed isothermally at 300 K (cold-compression) displays anomalously decreasing compressional (VP) and shear (VS) wave velocities with increasing pressure until ~8 GPa. Beyond 8 GPa, both VP and VS start to increase sharply with pressure and irreversible densification of the glass occurred, producing large hysteresis loops of velocities upon decompression. However, for the glass compressed isothermally at increasingly higher temperatures (hot-compression), the velocity minima gradually shift to lower pressures. At temperature close to the glass transition temperature Tg, the velocity minima disappear completely, displaying a monotonic increase of velocities during compression and higher VP and VS during decompression. In addition, constant-load heating-cooling experiments show that velocities generally decrease slightly with increasing temperature, but start to increase once heated above a threshold temperature (~650 K). During cooling the velocities increase almost linearly with decreasing temperature, resulting in higher velocities (~1.5–2.5% higher) when returned to 300 K. This implies that a temperature-induced densification may have occurred in the glass at high pressures. Raman spectra on recovered samples show that the hot-compressed and high-P heated glasses contain distinctly densified and depolymerized structural signatures compared to the initial glass and the cold-compressed glass below the velocity transition pressure PT (~8 GPa). Such densification may be attributed to the breaking of bridging oxygen bonds and compaction in the intermediate-range structure. Our results demonstrate that temperature can facilitate glass densification at high pressures and point out the importance of P-T history in understanding the elastic properties of silicate glasses. Comparison with melt velocity suggests that hot-compressed glasses may better resemble the pressure dependence of velocity of silicate melts than cold-compressed glasses, but still show significantly higher velocities than melts. If the abnormal acoustic behaviors of cold-compressed glasses were used to constrain melt fractions in the mantle low-velocity regions, the melt fractions needed to explain a given velocity reduction would be significantly underestimated at high pressures.

Geochemistry & Geophysics↗
NASA NTRS1994

The Balmer-dominated northeast limb of the Cygnus loop supernova remnant

We present a comprehensive investigation of the Balmar-dominated northeast limb of the Cygnus Loop supernova remnant. Data presented include H alpha (O III), and X-ray images, UV and visible spectrophotometry, and high-resolution spectroscophy. The two relatively bright Balmer-dominated filaments visible on the POSS prints are seen to be part of a very smooth and regular complex of filaments. These filaments mark the current location of the blast wave and are seen to bound the sharply limb-brightened X-ray emission, including the previously reported X-ray, 'halo.' The (O III)/h beta ratio throughout the region is approximately 0.1, except for regions in which the shock is undergoing a transition from nonradiative to incomplete radiative to incomplete radiative conditions. At these locations (O III) emission from the cooling region is quite strong, while collisionally excited Balmer-line emission can be weak because of photoionization of the preshock medium by UV from the nascent cooling region. As a result (O III)/H beta is greater than 100 in some locations. The nonradiative/radiative transition is best studied along the length of the northwestern of the two brightest filaments, where the shock velocity and swept-up column go from approximately 180 km/s and 10(exp 17)/sq cm at one end to approximately 140 km/s and 8 x 10(exp 17)/cm at the other. There are also a number of locations of such incomplete radiative emission where the shock has recently encountered denser regions with characteristic sizes of approximately 10(exp 18) cm. There is a considerable amount of evidence that the shock has decelerated from approximately 400 km/s to less than 200 km/s in the last 1000 yr. We interpret this as the result of the blast wave hitting the wall of a cavity which surround the supernova precursor and succeed in matching a wide range of data with a reflected shock model in which the density ofthe cavity wall is approximately 1.2/cu cm and the density in the interior of the cavity is about 0.09/cu cm. This analysis supports previous suggestions that the Cygnus Loop was a cavity explosion. We opredict that the X-ray shell will continue to brighten for the next 1000 years or so, during which time the shock front will become radiative throughout the northeast region. The similarity of the time scales for post-shock cooling, shock deceleration, and the 1ressure evolution of the region violates assumptions of steady flow, and greatly complicated interpretation of shock specptra in the region. This may be a common and important effect in many supernova remnants. High-resolutoin spectroscopy of the H alpha line profile shows that the broad component is surprisingly narrow (approximately 130 km/s). Consistency between this line width and the shock velocity inferred from spectrophotometry requires either rapid postshock equilibration of T(sub i) and T(sub e), or that the shock front has declerated from 180 km/s to 150 km/s in the last approximately 200 yr. While we cannot unambiguously distinquish between these models, the required deceleration is consistent with extrapolation of the shock decleration rate inferred over the last 1000 years. The width of the narrow component of the h alpha line profile is approximately 3 km/s, which indicates that the preshock medium has been heated to approximately 25,000 K by a shock precursor. The most promising types of precursor are heating either by fast neutrals which overtake the shock or by cosmic rays. In the latter case the observed properties of the precursor constrain the cosmic-ray diffusion coefficient to be between about 10(exp 23) sq cm/s. The characteristics of the precursor seem to be about right to explain the presence of weak (S II) and (N II) emission previously reported by Fesen & Itoh (1985).

Hester, J. Jeff↗
DOE OSTI2023

Uncertainty-quantified phenomenological optical potentials for single-nucleon scattering

Optical-model potentials (OMPs) continue to play a key role in nuclear reaction calculations. However, the uncertainty of phenomenological OMPs in widespread use—inherent to any parametric model trained on data—has not been fully characterized, and its impact on downstream users of OMPs remains unclear. Here we assign well-calibrated uncertainties for two representative global OMPs, those of Koning-Delaroche and Chapel Hill '89, using Markov-chain Monte Carlo for parameter inference. By comparing the canonical versions of these OMPs against the experimental data originally used to constrain them, we show how a lack of outlier rejection and a systematic underestimation of experimental uncertainties contributes to bias of, and overconfidence in, best-fit parameter values. Our updated, uncertainty-quantified versions of these OMPs address these issues and yield complete covariance information for potential parameters. Scattering predictions generated from our ensembles show improved performance both against the original training corpora of experimental data and against a new “test” corpus comprising many of the experimental single-nucleon scattering data collected over the last twenty years. Finally, we apply our uncertainty-quantified OMPs to two case studies of application-relevant cross sections. We conclude that, for many common applications of OMPs, including OMP uncertainty should become standard practice. Furthermore, to facilitate their immediate use, digital versions of our updated OMPs and related tools for forward uncertainty propagation are included as Supplemental Material.

150 ≤ A ≤ 189↗
DOE OSTI2024

Classification of compact objects and model comparison using EOS knowledge

Nuclear theory and experiments, alongside astrophysical observations, constrain the equation of state (EOS) of supranuclear-dense matter. Conversely, knowledge of the EOS allows an improved interpretation of nuclear or astrophysical data. In this article, we use several established constraints on the EOS and the new NICER measurement of PSR J0437-4715 to comment on the nature of the primary companion in GW230529 and the companion of PSR J0514-4002E. We find that, with a probability of ≳84% and ≳68%, respectively, both objects are black holes. These likelihoods increase to above 95% when one uses GW170817’s remnant as an upper limit on the TOV mass. We also demonstrate that the current knowledge of the EOS substantially disfavors high masses and radii for PSR J⁢0030+0451, inferred recently when combining NICER with XMM-Newton background data and using particular hot-spot models. Lastly, we also use our obtained EOS knowledge to comment on measurements of the nuclear symmetry energy, finding that the large value predicted by the PREX-II measurement displays some mild tension with other constraints on the EOS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Uncertainty quantification of optical models in fission fragment deexcitation

Here, we take the first step towards incorporating compound nuclear observables at astrophysically relevant energies into the experimental evidence used to constrain optical models, by propagating the uncertainty in two global optical potentials, one phenomenological and one microscopic, to correlated fission observables using the Monte Carlo Hauser-Feshbach formalism. We compare to a wide range of historic and recent experimental fission measurements, and discuss in detail regions of disagreement. We find that the parametric optical model uncertainty in neutron-fragment correlated observables involving neutron energy is significant. On the other hand, we observe that other experimental features, particularly neutron-fragment correlations near the 132 Sn shell closure and the high energy component of neutron spectra, are unlikely to be explained by the optical potential, and will require further experimental and theoretical effort to explain.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRSDate not supplied

Nancy Grace Roman Space Telescope Observatory Implementation and Challenges

NASA’s Nancy Grace Roman Space Telescope (Roman), previously referred to as Wide Field Infrared Survey Telescope (WFIRST), was named after Dr. Nancy Grace Roman, an astronomer and NASA pioneer of modern space-based astronomy who is known as the “mother of the Hubble Space Telescope”. Roman is a deep space infrared observatory with a Hubble-sized telescope and wide field of view instrument (greater than 100 times that of Hubble’s) that will conduct a high latitude time-domain survey, a high latitude imaging and spectroscopic survey, and a galactic bulge time-domain survey to characterize dark energy and expand the census of exoplanets in our galaxy while allowing a broad range of astrophysics research. Roman will also demonstrate exoplanet coronagraphy with active wave front control technology and provide general investigator programs for the science community. Roman is finishing the critical design phase and is planning for launch in 2026. It will operate in a quasi-halo orbit about Sun-Earth L2, 1.5 million kilometers from Earth, for a five-year primary mission life. The Observatory features an Optical Telescope Assembly with an existing, repurposed 2.4m primary mirror, a Wide Field Instrument with a focal plane array comprised of 18 HgCdTe near-infrared detectors and a grism, prism and filter elements for imaging and spectroscopy in support of the primary surveys, as well as a Coronagraph instrument technology demonstration with starlight suppression technology for direct imaging and spectroscopy of exoplanets. The telescope is mounted to the Instrument Carrier composite truss structure which also optically meters each instrument, includes a Launch Load and Vibration Isolation System to provide passive isolation of spacecraft jitter sources while also supporting the payload during launch and is attached to the Spacecraft Bus. The Spacecraft also includes a Solar Array Sunshield ,Deployable Aperture Cover, Lower Instrument Sunshade, High Gain Antenna System, and Outer Barrel Assembly. Figure 1 shows an overview of the Roman Observatory. When fully integrated, Roman will be the largest Observatory assembled and tested at NASA’s Goddard Space Flight Center. Figure 1. Roman Observatory Overview Development of scientific satellites is challenging by nature, as the pursuit to broaden scientific knowledge always pushes the boundary of what has come before. The implementation of the Roman mission is a prime example and expected challenges have been augmented by the foundational decision to use the existing telescope components, developed in the early 2000s by another Government agency for a different application. Other unique aspects of the Roman mission, such as its survey nature, the vast amount of data required to meet science objectives, and packaging of the Observatory elements around the existing telescope components, create constrained design spaces that drive competing requirements across Observatory subsystems. Given these challenges, systems engineering has been a critical discipline in balancing implementation decisions for the Roman mission and will continue to play a key role going forward. This paper will discuss details of the Roman Observatory configuration, as well as some of the systems engineering challenges and the decision-making process used to mature the Roman Space Telescope preliminary design to implementation.

Lisa Ml Bartusek↗
DOE OSTI2024

Dynamical ejecta from binary neutron star mergers: Impact of a small residual eccentricity and of the equation of state implementation

Predicting the properties of the matter ejected during and after a neutron star merger is crucial to our ability to use electromagnetic observations of these mergers to constrain the masses of the neutron stars, the equation of state of dense matter, and the role of neutron star mergers in the enrichment of the Universe in heavy elements. Furthermore, our ability to reliably provide such predictions is however limited by a broad range of factors, including the finite resolution of numerical simulations, their treatment of magnetic fields, neutrinos, and neutrino-matter interactions, and the approximate modeling of the equation of state of dense matter. In this manuscript, we study specifically the role that a small residual eccentricity and different implementations of the same equation of state have on the matter ejected during the merger of a 1.3M ⊙ –1.4M ⊙ binary neutron star system. We find that a residual eccentricity e ~ 0.01, as measured ~ 4–6 orbits before merger, causes O(25%–30%) changes in the amount of ejected mass, mainly due to changes in the amount of matter ejected as a result of core bounces during merger. We note that O(1%) residual eccentricities have regularly been used in binary neutron star merger simulations as proxy for circular binaries, potentially creating an additional source of error in predictions for the mass of the dynamical ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Renormalized classical theory of quantum magnets

Here, we derive a renormalized classical spin (RCS) theory for 𝑆 >1/2 quantum magnets by constraining a generalized classical theory that includes all multipolar fluctuations to a reduced CP 1 phase space of dipolar SU(2) coherent states. When the spin Hamiltonian $\hat{ℋ}$(𝑆) is linear in the spin operators $\hat{𝑺}$ 𝑗 for each lattice site 𝑗, the RCS Hamiltonian $\tilde{ℋ}$ cl coincides with the usual classical model ℋ cl = lim 𝑆→∞⁡ $\hat{ℋ}$(𝑆). In the presence of nonlinear terms, however, the RCS theory is more accurate than ℋ cl . For the many materials modeled by spin Hamiltonians with (nonlinear) single-ion anisotropy terms, the use of the RCS theory is essential to accurately model phase diagrams and to extract the correct Hamiltonian parameters from neutron-scattering data.

magnetic anisotropy↗
NASA NTRS2020

Mineralogy and Geochemistry of Sedimentary Rocks and Eolian Sediments in Gale Crater, Mars: A Review After Six Earth Years of Exploration with Curiosity

The Mars Science Laboratory Curiosity rover arrived at Mars in August 2012 with a primary goal of characterizing the habitability of ancient and modern environments. Curiosity was sent to Gale crater to study a sequence of ∼3.5 Ga old sedimentary rocks that, based on orbital visible and near- to short-wave infrared reflectance spectra, contain secondary minerals that suggest deposition and/or alteration in liquid water. The sedimentary sequence in the lower slopes of Mount Sharp in Gale crater preserves a dramatic shift on early Mars from a relatively warm and wet climate to a cold and dry climate, based on a transition from smectite-bearing strata to sulfate-bearing strata. The rover is equipped with instruments to examine the sedimentology and identify compositional changes in the stratigraphy. The Chemistry and Mineralogy (CheMin) instrument is one of two internal laboratories on Curiosity and includes a transmission X-ray diffractometer (XRD) and X-ray fluorescence (XRF) spectrometer. CheMin measures loose sediment samples scooped from the surface and drilled rock powders, and the XRD provides quantitative mineralogy to a detection limit of ∼1 wt.% for crystalline phases. Curiosity has traversed>20 km since landing and has primarily been exploring an ancient lake environment fed by streams and groundwater. Of the 19 drilled rock samples analyzed by CheMin as of sol 2300 (January 2019), 15 are from fluvio-lacustrine deposits that comprise the Bradbury and Murray formations. Most of these samples were drilled from units that did not have a clear mineralogical signature from orbit. Results from CheMin demonstrate an astounding diversity in the mineralogy of these rocks that signifies geochemical variations in source rocks, transportation mechanisms, and depositional and diagenetic fluids. Most detrital igneous minerals are basaltic, but the discovery in a few samples of abundant silicate minerals that usually crystallize from evolved magmas on Earth remains enigmatic. Trioctahedral smectite and magnetite at the base of the section may have formed from low-salinity pore waters with a circumneutral pH in lake sediments. A transition to dioctahedral smectite, hematite, and Ca-sulfate going up section suggests a change to more saline and oxidative aqueous conditions in the lake waters themselves and/or in diagenetic fluids. Perhaps one of the biggest mysteries revealed by CheMin is the high abundance of X-ray amorphous materials (15–73 wt.%) in all samples drilled or scooped to date. CheMin has analyzed three modern eolian sands, which have helped constrain sediment transport and mineral segregation across the active Bagnold Dune Field. Ancient eolian sandstones drilled from the Stimson formation differ from modern eolian sands in that they contain abundant magnetite but no olivine, suggesting that diagenetic processes led to the alteration of olivine to release Fe(II) and precipitate magnetite. Fracture-associated halos in the Stimson and the Murray formations are evidence for complex aqueous processes long after the streams and lakes vanished from Gale crater. The sedimentology and composition of the rocks analyzed by Curiosity demonstrate that habitable environments persisted intermittently on the surface or in the subsurface of Gale crater for perhaps more than a billion years.

E.B. Rampe↗
DOE OSTI2026

Asymptotic-preserving semi-implicit finite volume scheme for extended magnetohydrodynamics

A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such that it allows the algorithm to retain the use of ideal MHD Riemann solvers and the constrained transport method to preserve divergence-free magnetic fields. Incorporation of electron inertia and displacement current introduces additional numerical stiffness which motivates a semi-implicit FV scheme that re-formulates the XMHD model as a relaxation system. The equations are then advanced in time using an explicit 2nd-order Runge–Kutta scheme with operator splitting applied to the implicit source term updates at each sub-stage. For additional numerical stability, a density-dependent slope limiter is implemented to increase flux diffusivity at low density regions where non-ideal effects become significant. The algorithm is subsequently implemented in a scalable adaptive mesh refinement (AMR) framework. As the new algorithm retains many aspects of the ideal MHD formulations, it asymptotes naturally to the ideal MHD limit. Moreover, it shows promising results at the resistive and Hall MHD limits. This is verified against reference test problems for ideal, resistive and Hall MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

A Compact Dual-Band Bowtie Antenna for RF and ISM bands Operation

Traditionally, bowtie antennas have been known to exhibit wide impedance characteristics, omnidirectional radiation patterns, and linear polarization. There is a broad range of applications from medical imaging, archaeological survey, and Ground Penetrating Radar (GPR) to trackers and sensor networks where wideband bowtie antenna designs are required for their operations. Broadening the bandwidth of bowtie antenna requires widening the flare angle of the bowtie arms, which consequently results in a large surface area that may not be suitable for space-constrained applications. Moreover, drawback attributes of the wideband bowtie designs feature inconsistent radiation pattern across the bandwidth and low signal-to-noise (SNR) ratios. As it is known, the SNR would be improved in dual- or multi-band antennas due to their reduced bandwidth. To this end, dual-band/multi-band antennas are preferred over wideband antennas in applications where more than a single frequency of interest is present. Previously, a two-port double-dipole elements was reported, whose arms were orthogonally interleaved to facilitate operation in both the standard WLAN frequency bands (J. M. Steyn and et. al, Progress in Electromagnetic Research, Vol. – 10 pp. 151-161). Even though the antenna is not very compact it does exhibit good cross-polarization, moderate gain in both frequency bands. Another dual-band bowtie antenna which excites two bands using a single transmission line was reported (Wen Chao Zheng and et. al, IEEE Trans. Antennas propag., 2014). The design was compact and did not need multi-port feeding network. In this paper, a dual-band compact bowtie antenna operating at 900 MHz (RF band) and 2.45 GHz (ISM band) using a single excitation port is introduced. It is printed on a 1.54mm thick dielectric substrate (εr = 3.38). The antenna consists of two sets of bowtie arms, a microstrip transmission line to feed the bowtie arms, and a ground plane acting as a reflector to partially reduce the back radiation. One of the bowtie arms of each frequency is printed on the top layer and the other arm, which is mirror imaged, is printed on the bottom layer of the substrate. The length of the bowtie controls the resonance frequency of the antenna and the flare angle controls the bandwidth of the antenna. The microstrip transmission line, connected to a 50 Ω SMA probe, feeds the bowtie antenna. The compact antenna can be used for both RF and ISM band applications. The bowtie arms at the lower frequency band are miniaturized by elongating their electrical lengths. The influence of miniaturizing the bowtie arms and the supporting partial ground plane is observed in the reduced peak gain and degraded front-to-back ratio. These are partly neutralized using four quarter-wave choke-slots in the ground plane with two on each side of the feeding transmission line. The proposed antenna is numerically investigated and finalized by the finite-element based full-wave EM solver, ANSYS HFSS. The miniaturization has reduced the ground plane size by ~45% and the arms size by ~41%. In addition to the size reduction benefits, the antenna shows reasonable peak gain and front-to-back ratio in both the bands. The corresponding results will be presented and discussed at the conference.

Saininad Naik↗
DOE OSTI2025

Direct nonparametric multimessenger constraints on the equation of state of cold dense nuclear matter

We utilize the now substantial amount of astrophysical observations of neutron stars (NSs), along with perturbative quantum chromodynamics (pQCD) calculations at high density, to directly constrain the NS equation of state (EOS). To this end, we construct nonparametric EOS priors by using Gaussian processes trained on 75 EOSs, which include models with either hadrons, hyperons, or quarks at high densities. We create a prior using the full EOS sample (model agnostic), and one prior for each EOS family to test model discrimination. We introduce a novel inference approach, which allows the simultaneous sampling of intrinsic and extrinsic parameters of binary NS mergers, as well as a nonparametric equation of state. We showcase this method in a Bayesian updating scheme by first performing a complete analysis of the binary NS merger event GW170817 with minimal assumptions, and sequentially adding information from x-ray and radio NS observations, along with pQCD calculations. Besides providing standard constraints, such as the pressure at twice nuclear saturation density 𝑝⁡(2⁢𝜌 sat ) = 4.3$^{+0.6}_{−0.6}$ × 10 34 dyne/cm 2 , at 95% confidence level, for the model agnostic prior, our methodology shows how the choice of EOS families used in conditioning changes the inferred astrophysical properties of the EOS, namely tidal deformability and maximum supported NS mass. We find hyperonic priors predicting higher tidal deformabilities for a 1.4⁢𝑀 ⊙ NS, and hadronic priors being preferred by the considered astrophysical data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗