Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stefan model”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Rapid simulations of hyperspectral near-field images of three-dimensional heterogeneous surfaces – part II

The modeling of the near-field interaction in the scattering-type scanning near-field optical microscope (s-SNOM) is rapidly advancing, although an accurate yet versatile modeling framework that can be easily adapted to various complex situations is still lacking. In this work, we propose a time-efficient numerical scheme in the quasi-electrostatic limit to capture the tip-sample interaction in the near field. This method considers an extended tip geometry, which is a significant advantage compared to the previously reported method based on the point-dipole approximation. Using this formalism, we investigate, among others, nontrivial questions such as uniaxial and biaxial anisotropy in the near-field interaction, the relationship between various experimental parameters (e.g. tip radius, tapping amplitude, etc.), and the tip-dependent spatial resolution. The demonstrated method further sheds light on the understanding of the contrast mechanism in s-SNOM imaging and spectroscopy, while also representing a valuable platform for future quantitative analysis of the experimental observations.

Chen, Xinzhong (ORCID:0000000261033848)↗
DOE OSTI2024

Insights into the laser-assisted photoelectric effect from solid-state surfaces

Photoemission from a solid surface provides a wealth of information about its electronic structure and dynamic evolution. Ultrafast pump-probe experiments offer real-time access to photon-surface interactions and the resulting electron dynamics. Here, we present a femtosecond time-resolved photoelectron spectroscopy study of laser-assisted photoemission (LAPE) from two different metal surfaces, tungsten and platinum. Utilizing synchronized IR laser and x-ray pulses, photoelectron sideband generation up to the sixth order is observed. A significant material-dependent variation of the LAPE response has not been predicted by previous theoretical models, to the best of our knowledge. The observed phenomena are semiquantitatively reproduced by considering distinct dynamic dielectric responses of the two materials. In conclusion, these findings provide a deeper understanding of the LAPE process and insights into the dynamic interplay between optical laser fields and metal surfaces.

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

Low-lying dipole response of Ni 64

In this work, two complementary real-photon scattering experiments were conducted on the proton-magic 64 Ni nucleus to study the dipole response up to its neutron-separation energy of $S_n$ = 9.7 MeV. By combining both measurements, 87 $E$1 and 23 $M$1 transitions were identified above 4.3 MeV. The results of the observed $M$1 transitions were compared to shell-model calculations using two different model spaces. It was found that the inclusion of excitations across the $Z$ = 28 shell gap in the calculations has a large impact. Furthermore, average cross sections for decays to the ground state (elastic transitions) as well as to lower-lying excited states (inelastic decays) were determined. The corresponding $E$1 channel was compared to calculations within the relativistic equation of motion (REOM) framework. Whereas the calculations of highest possible complexity reproduce the fragmentation and overall behavior of the $E$1 average elastic cross section well, the predicted absolute cross sections are approximately twice as high as the experimental upper limits even though the latter also include an estimate of the inelastic-decay channel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A Multi-physics Computational Investigation of Droplet Pathogen Transport Eminating from Synthetic Coughs and Breathing

In response to the global SARS-CoV-2 transmission pandemic, Sandia National Laboratories Rapid Lab-Directed Research and Development (LDRD) COVID-19 initiative has deployed a multi-physics, droplet-laden, turbulent low-Mach simulation tool to model pathogen-containing water droplets that emanate from synthetic human coughing and breathing. The low-Mach turbulent Eulerian/point-particle Lagrangian methodology directly couples mass, momentum, energy, and species to capture droplet evaporation physics that supports the ability to distinguish between droplets that deposit and those that persist in the environment. Additionally, the cough mechanism is modeled as a pulsed-spray with a prescribed log-normal droplet size distribution. Simulations demonstrate direct droplet deposition lengths in excess of three meters while the persistence of a droplet nuclei entrained within a buoyant plume is noted. Including the effect of protective barriers demonstrates effective mitigation of large-droplet transport. For coughs into a protective barrier, jet impingement and large-scale re-circulation can drive droplets vertically and back towards the subject while supporting persistence of droplet nuclei. Simulations in quiescent conditions demonstrate droplet preferential concentrations due to the coupling between vortex ring shedding and the subsequent advection of a series of three-dimensional rings that tilt and rise vertically due to a misalignment between the initial principle vortex trajectory and gravity. These resolved coughing simulations note vortex ring formation, roll-up and breakdown, while entraining droplet nuclei for large distances and time scales.

42 ENGINEERING↗
DOE OSTI2020

Taking stock of national climate policies to evaluate implementation of the Paris Agreement

Many countries have implemented national climate policies to accomplish pledged Nationally Determined Contributions and to contribute to the temperature objectives of the Paris Agreement on climate change. In 2023, the global stocktake will assess the combined effort of countries. Here, based on a public policy database and a multi-model scenario analysis, we show that implementation of current policies leaves a median emission gap of 22.4 to 28.2 GtCO 2 eq by 2030 with the optimal pathways to implement the well below 2 °C and 1.5 °C Paris goals. If Nationally Determined Contributions would be fully implemented, this gap would be reduced by a third. Interestingly, the countries evaluated were found to not achieve their pledged contributions with implemented policies (implementation gap), or to have an ambition gap with optimal pathways towards well below 2 °C. This shows that all countries would need to accelerate the implementation of policies for renewable technologies, while efficiency improvements are especially important in emerging countries and fossilfuel- dependent countries.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Cloudy-sky contributions to the direct aerosol effect

The radiative forcing of the aerosol–radiation interaction can be decomposed into clear-sky and cloudy-sky portions. Two sets of multi-model simulations within Aerosol Comparisons between Observations and Models (AeroCom), combined with observational methods, and the time evolution of aerosol emissions over the industrial era show that the contribution from cloudy-sky regions is likely weak. A mean of the simulations considered is 0.01±0.1 W m-2. Multivariate data analysis of results from AeroCom Phase II shows that many factors influence the strength of the cloudy-sky contribution to the forcing of the aerosol–radiation interaction. Overall, single-scattering albedo of anthropogenic aerosols and the interaction of aerosols with the short-wave cloud radiative effects are found to be important factors. A more dedicated focus on the contribution from the cloud-free and cloud-covered sky fraction, respectively, to the aerosol–radiation interaction will benefit the quantification of the radiative forcing and its uncertainty range.

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

A pathogenic role for histone H3 copper reductase activity in a yeast model of Friedreich’s ataxia

Disruptions to iron-sulfur (Fe-S) clusters, essential cofactors for a broad range of proteins, cause widespread cellular defects resulting in human disease. A source of damage to Fe-S clusters is cuprous (Cu 1+ ) ions. Since histone H3 enzymatically produces Cu 1+ for copper-dependent functions, we asked whether this activity could become detrimental to Fe-S clusters. Here, we report that histone H3–mediated Cu 1+ toxicity is a major determinant of cellular functional pool of Fe-S clusters. Inadequate Fe-S cluster supply, due to diminished assembly as occurs in Friedreich’s ataxia or defective distribution, causes severe metabolic and growth defects in Saccharomyces cerevisiae. Decreasing Cu 1+ abundance, through attenuation of histone cupric reductase activity or depletion of total cellular copper, restored Fe-S cluster–dependent metabolism and growth. Our findings reveal an interplay between chromatin and mitochondria in Fe-S cluster homeostasis and a potential pathogenic role for histone enzyme activity and Cu 1+ in diseases with Fe-S cluster dysfunction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Broadband X-ray ptychography using multi-wavelength algorithm

Ptychography is a rapidly developing scanning microscopy which is able to view the internal structures of samples at a high resolution beyond the illumination size. The achieved spatial resolution is theoretically dose-limited. A broadband source can provide much higher flux compared with a monochromatic source; however, it conflicts with the necessary coherence requirements of this coherent diffraction imaging technique. In this paper, a multi-wavelength reconstruction algorithm has been developed to deal with the broad bandwidth in ptychography. Compared with the latest development of mixed-state reconstruction approach, this multi-wavelength approach is more accurate in the physical model, and also considers the spot size variation as a function of energy due to the chromatic focusing optics. Therefore, this method has been proved in both simulation and experiment to significantly improve the reconstruction when the source bandwidth, illumination size and scan step size increase. It is worth mentioning that the accurate and detailed information of the energy spectrum for the incident beam is not required in advance for the proposed method. Further, we combine multi-wavelength and mixed-state approaches to jointly solve temporal and spatial partial coherence in ptychography so that it can handle various disadvantageous experimental effects. The significant relaxation in coherence requirements by our approaches allows the use of high-flux broadband X-ray sources for high-efficient and high-resolution ptychographic imaging.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

An overview of the Western United States Dynamically Downscaled Dataset (WUS-D3)

Abstract. Predicting future climate change over a region of complex terrain, such as the western United States (US), remains challenging due to the low resolution of global climate models (GCMs). Yet the climate extremes of recent years in this region, such as floods, wildfires, and drought, are likely to intensify further as climate warms, underscoring the need for high-quality and high-resolution predictions. Here, we present an ensemble of dynamically downscaled simulations over the western US from 1980–2100 at 9 km grid spacing, driven by 16 latest-generation GCMs. This dataset is titled the Western US Dynamically Downscaled Dataset (WUS-D3). We describe the challenges of producing WUS-D3, including GCM selection and technical issues, and we evaluate the simulations' realism by comparing historical results to temperature and precipitation observations. The future downscaled climate change signals are shaped in physically credible ways by the regional model's more realistic coastlines and topography. (1) The mean warming signals are heavily influenced by more realistic snowpack. (2) Mean precipitation changes are often consistent with wetting on the windward side of mountain complexes, as warmer, moister air masses are uplifted orographically during precipitation events. (3) There are large fractional precipitation increases on the lee side of mountain complexes, leading to potentially significant changes in water resources and ecology in these arid landscapes. (4) Increases in precipitation extremes are generally larger than in the GCMs, driven by locally intensified atmospheric updrafts tied to sharper, more realistic gradients in topography. (5) Changes in temperature extremes are different from what is expected by a shift in mean temperature and are shaped by local atmospheric dynamics and land surface feedbacks. Because of its high resolution, comprehensiveness, and representation of relevant physical processes, this dataset presents a unique opportunity to evaluate societally relevant future changes in western US climate.

Rahimi, Stefan (ORCID:0000000331884462)↗
DOE OSTI2025

HFBTHO-AD: Differentiation of a nuclear energy density functional code

The HFBTHO code implements a nuclear energy density functional solver to model the structure of atomic nuclei. HFBTHO has previously been used to calibrate energy functionals and perform sensitivity analysis by using derivative-free methods. To enable derivative-based optimization and uncertainty quantification approaches, we must compute the derivatives of HFBTHO outputs with respect to the parameters of the energy functional, which are a subset of all input parameters of the code. Here, we use the algorithmic/automatic differentiation (AD) tool Tapenade to differentiate HFBTHO. We compare the derivatives obtained using AD against finite-difference approximation and examine the performance of the derivative computation.

Algorithmic differentiation↗
DOE OSTI2026

Goal-oriented real-time Bayesian inference for linear autonomous dynamical systems with application to digital twins for tsunami early warning

We present a goal-oriented framework for constructing digital twins with the following properties: (1) they employ discretizations of high-fidelity partial differential equation (PDE) models governed by autonomous dynamical systems, leading to large-scale forward problems; (2) they solve a linear inverse problem to assimilate observational data to infer uncertain model components followed by a forward prediction of the evolving dynamics; and (3) the entire end-to-end, data-to-inference-to-prediction computation is carried out without approximation and in real time through a Bayesian framework that rigorously accounts for uncertainties. Several challenges must be overcome to realize this framework, including the large scale of the forward problem, the high dimensionality of the parameter space, and for a class of problems including those we target, the slow decay of the singular values of the parameter-to-observable map. Here we introduce a methodology to overcome these challenges by exploiting the autonomous structure of the forward model to decompose the solution of the inverse problem into a one-time-only offline phase in which the PDE model is solved a limited number of times (equal to the number of sensors), and an online phase that maps well onto GPUs and computes the parameter inference and prediction of quantities of interest in real time, given observational data. Our ultimate goal is to apply this framework to construct digital twins for subduction zones, including Cascadia, to provide early warning for tsunamis generated by megathrust earthquakes. To this end, we demonstrate how our methodology can be used to employ seafloor pressure observations, along with the coupled acoustic–gravity wave equations, to infer the earthquake-induced spatiotemporal seafloor motion (discretized with $\mathscr{O}$ (10 9 ) parameters) and forward predict the tsunami propagation. We present results of an end-to-end inference, prediction, and uncertainty quantification for a representative test problem with $\mathscr{O}$ (10 8 ) inversion parameters for which goal-oriented Bayesian inference is accomplished exactly and in real time, that is, in a matter of seconds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Non-local transport in radiation-hydrodynamics codes for ICF by efficient coupling to an external Vlasov–Fokker–Planck code

Accurately incorporating non-local transport into radiation-hydrodynamics codes, and indeed any fluid system, has long been elusive. To date, a simplified and accurate theory that can be easily integrated has not been available. This limitation affects modeling in inertial confinement fusion (ICF) and magnetic confinement fusion systems, among others, where non-local transport is well-known to be present. Here, we present a coupling methodology between a full Vlasov–Fokker–Planck (VFP) electron kinetic code and radiation-hydrodynamics (rad-hydro) codes. The VFP code is used to adjust native electron transport in the rad-hydro code, thus enabling improved transport without the need to integrate a full electron VFP solver into the rad-hydro code. This approach necessitates only occasional invocation of the VFP code, reducing computational intensity compared to following the dynamic evolution entirely with the VFP code on fluid time scales. We illustrate that the methodology is more accurate than other simplified methods in thermal decay systems relevant to ICF and can replicate standard theoretical results with high accuracy.

Electronic transport↗
DOE OSTI2025

Probing the Kitaev honeycomb model on a neutral-atom quantum computer

Quantum simulations of many-body systems are among the most promising applications of quantum computers. In particular, models based on strongly correlated fermions are central to our understanding of quantum chemistry and materials problems, and can lead to exotic, topological phases of matter. However, owing to the non-local nature of fermions, such models are challenging to simulate with qubit devices. Here we realize a digital quantum simulation architecture for two-dimensional fermionic systems based on reconfigurable atom arrays. We utilize a fermion-to-qubit mapping based on Kitaev’s model on a honeycomb lattice, in which fermionic statistics are encoded using long-range entangled states. We prepare these states efficiently using measurement and feedforward, realize subsequent fermionic evolution through Floquet engineering with tunable entangling gates interspersed with atom rearrangement, and improve results with built-in error detection. Leveraging this fermion description of the Kitaev spin model, we efficiently prepare topological states across its complex phase diagram and verify the non-Abelian spin-liquid phase by evaluating an odd Chern number. We further explore this two-dimensional fermion system by realizing tunable dynamics and directly probing fermion exchange statistics. Finally, we simulate strong interactions and study the dynamics of the Fermi–Hubbard model on a square lattice. These results pave the way for digital quantum simulations of complex fermionic systems for materials science, chemistry and high-energy physics.

atomic and molecular physics↗
DOE OSTI2025

Programmable simulations of molecules and materials with reconfigurable quantum processors

Simulations of quantum chemistry and quantum materials are believed to be among the most important applications of quantum information processors. However, realizing practical quantum advantage for such problems is challenging because of the prohibitive computational cost of programming typical problems into quantum hardware. Here we introduce a simulation framework for strongly correlated quantum systems represented by model spin Hamiltonians that uses reconfigurable qubit architectures to simulate real-time dynamics in a programmable way. Our approach also introduces an algorithm for extracting chemically relevant spectral properties via classical co-processing of quantum measurement results. We develop a digital–analogue simulation toolbox for efficient Hamiltonian time evolution using digital Floquet engineering and hardware-optimized multi-qubit operations to accurately realize complex spin–spin interactions. As an example, we propose an implementation based on Rydberg atom arrays. In addition, we show how detailed spectral information can be extracted from the dynamics through snapshot measurements and single-ancilla control, enabling the evaluation of excitation energies and finite-temperature susceptibilities from a single dataset. To illustrate the approach, we show how to use the method to compute key properties of a polynuclear transition-metal catalyst and two-dimensional magnetic materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Establishing structure-sensitivity of ceria reducibility: real-time observations of surface-hydrogen interactions

The first layer of atoms on an oxide catalyst provides the first sites for adsorption of reactants and the last sites before products or oxygen are desorbed. We employ a unique combination of morphological, structural, and chemical analyses of a model ceria catalyst with different surface terminations under an H 2 environment to unequivocally establish the effect of the last layer of atoms on surface reduction. (111) and (100) terminated epitaxial islands of ceria are simultaneously studied in situ allowing for a direct investigation of the structure–reducibility relationship under identical conditions. Kinetic rate constants of Ce 4+ to Ce 3+ transformation and equilibrium concentrations are extracted for both surface terminations. Unlike the kinetic rate constants, which are practically the same for both types of islands, more pronounced oxygen release, and overall higher reducibility were observed for (100) islands compared to (111) ones. The findings are in agreement with coordination-limited oxygen vacancy formation energies calculated by density functional theory. The results point out the important aspect of surface terminations in redox processes, with particular impact on the catalytic reactions of a variety of catalysts.

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

TDCOSMO - XVI. Measurement of the Hubble constant from the lensed quasar WGD 2038–4008

Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant (H0). The TDCOSMO Collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038−4008 using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain DΔ t = 1.68−0.38+0.40 Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat ΛCDM cosmology, we constrain the Hubble constant to be H0 = 65−14+23 km s−1 Mpc−1. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory’s Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.Key words: gravitational lensing: strong / cosmological parameters / distance scale⋆ Corresponding author; kcwong19@gmail.com.⋆⋆ NHFP Einstein fellow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Simulating the Multi-epoch Direct Detection Technique to Isolate the Thermal Emission of the Non-transiting Hot Jupiter HD187123b

We report the 6.5σ detection of water from the hot Jupiter HD187123b with a Keplerian orbital velocity K {sub p} of 53 ± 13 km s{sup −1}. This high-confidence detection is made using a multi-epoch, high-resolution, cross-correlation technique, and corresponds to a planetary mass of 1.4{sub −0.3}{sup +0.5} M {sub J} and an orbital inclination of 21° ± 5°. The technique works by treating the planet/star system as a spectroscopic binary and obtaining high signal-to-noise, high-resolution observations at multiple points across the planet’s orbit to constrain the system’s binary dynamical motion. All together, seven epochs of Keck/NIRSPEC L-band observations were obtained, with five before the instrument upgrade and two after. Using high-resolution SCARLET planetary and PHOENIX stellar spectral models, we were able to drastically increase the confidence of the detection by running simulations that could reproduce, and thus remove, the nonrandom structured noise in the final likelihood space well. The ability to predict multi-epoch results will be extremely useful for furthering the technique. Here, we use these simulations to compare three different approaches to combining the cross correlations of high-resolution spectra and find that the Zucker log(L) approach is least affected by unwanted planet/star correlation for our HD187123 data set. Furthermore, we find that the same total signal-to-noise ratio (S/N) spread across an orbit in many, lower S/N epochs rather than fewer, higher S/N epochs could provide a more efficient detection. This work provides a necessary validation of multi-epoch simulations, which can be used to guide future observations and will be key to studying the atmospheres of farther separated, non-transiting exoplanets.

79 ASTRONOMY AND ASTROPHYSICS↗