Search NASA⌕ Search

SEARCH · Search NASA

Results for “Observing System Simulation”

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 127 records · Page 7

Multiple Coulomb scattering in acrylic of a 221.3 MeV therapeutic proton beam

Measurements of multiple Coulomb scattering (MCS) distributions for 221.3 MeV therapeutic protons are presented using a novel detector system comprised of a thin scintillator, a pellicle mirror, and a digital camera. The MCS distributions were characterized for three acrylic phantoms of varying lengths and for two biological density-equivalent phantoms simulating bone and muscle. Additionally, beam profiles were measured across an energy range of 80.3–221.3 MeV in 20 MeV increments. The observed energy dependence of the photon yields is consistent with the tabulated stopping power values. Finally, the experimental results are benchmarked against Geant4 simulations, demonstrating consistent agreement and validating the capability of the detector system for radiology measurements.

Digital camera↗

The impact of cloud microphysics and ice nucleation on Southern Ocean clouds assessed with single-column modeling and instrument simulators

Abstract. Supercooled liquid clouds are common at higher latitudes (especially over the Southern Ocean) and are critical for constraining climate projections. We take advantage of the Macquarie Island Cloud and Radiation Experiment (MICRE) to perform an analysis of observed and simulated cloud processes over the Southern Ocean in a region and season dominated by supercooled liquid clouds. Using a single-column version of the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS), we compare two different cloud microphysical schemes to ground-based observations of cloud, precipitation, and radiation over a 2.5-month period (1 January–17 March 2017). Both schemes are able to reproduce aspects of the cloud and radiation observations during MICRE to within the uncertainty of the data when the thermodynamic profile is prescribed with relaxation. There are differences in water mass and representation of reflectivity between the schemes. A sensitivity study of the cloud microphysics schemes, one a bulk one-moment scheme and the other a two-moment scheme with prediction of mass and number, indicates that several key processes create differences between the schemes. Surface radiative fluxes and total water path are highly sensitive to the formation and fall speed of precipitation. The prediction of hydrometeor number with the two-moment scheme yields a better comparison with observed reflectivity and radiative fluxes, despite predicting higher liquid water contents than observed. With the two-moment scheme, we are also able to test the sensitivity of the results to the input of liquid cloud condensation nuclei (CCN) and ice nuclei (IN). The cloud properties and resulting radiative effects are found to be sensitive to the CCN and IN concentrations. More CCN and IN increase liquid and ice water paths, respectively. Thus, both the dynamic environment and aerosols, integrated through the cloud microphysics, are important for properly representing Southern Ocean cloud radiative effects.

58 GEOSCIENCES↗

Monthly hydropower generation data for Western Canada to support Western-US interconnect power system studies

Hydroelectric power generation in Western Canada significantly contributes to power grid operations of the North American Western Interconnection through substantial generation, some of which is exported to the United States (U.S.). However, the lack of publicly available hydropower generation datasets poses challenges for future market projections and resource adequacy evaluations. We present a simulation-based monthly power system model-ready hydropower generation dataset for 110 facilities in British Columbia and Alberta from 1981 to 2019. These monthly hydropower generation estimates are developed from integrated hydrologic model simulations of runoff and reservoir-operated streamflow, followed by scaling that considers diversion inflow constraints based on hydropower water license information. To address the lack of comparable hydropower generation records, we conduct step-by-step evaluations for simulated runoff, regulated streamflow, and hydropower generation using available observations or estimates. The presented hydropower dataset aims to enhance the representation of hydropower resources in Western Canada, supporting power grid system studies for the Western Interconnection of the U.S. and Canada.

13 HYDRO ENERGY↗

Synthetic Tuning of Exciton–Phonon Coupling in Janus WS 2(1-x) Se 2x Monolayers Revealed by Resonant Raman Excitation Spectroscopy for Optoelectronic Applications

Janus monolayers, such as WSSe, have broken out-of-plane symmetry and an intrinsic dipole moment, impacting exciton transport, lifetime, and phonon interactions while imbuing piezoelectric, photocatalytic, and Rashba spin-splitting properties to transition metal dichalcogenides (TMDs). The new properties of this atomically thin material can be used for optoelectronic device applications. As TMDs are converted into Janus monolayers, e.g., top selenization of WS2 to WSSe, the bandgap and structure smoothly evolve, impacting not only the formation of excitons but also their complex interactions with different phonon modes. Resonant Raman excitation profiles (REPs) are uniquely well-suited to reveal both excitonic transitions and exciton–phonon coupling. Here, the resonant REPs of $A^{'}_{1}$ WS 2 and A 1 WSSe modes are measured to understand the strength of their coupling with the A, B, and C excitonic bands of a WS2 monolayer throughout its stepwise transformation into Janus WSSe by pulsed laser deposition (PLD) of energetic selenium species. In situ Raman spectroscopy during deposition is used to controllably prepare stable intermediate Janus structures, WS 2(1-x) Se 2x (0 ≤ x ≤ 0.5), for ex situ measurement of their resonant REPs. As x increases, REPs reveal not only pronounced excitonic bands that gradually shift toward lower photon energies but also strong, mode-selective exciton–phonon coupling. First-principles resonant Raman simulations independently predict this spectral behavior and are shown capable of matching the spectrally broadened, experimentally observed REP profiles in this model system, indicating their strong predictive capability for future experiments. The combination of controlled synthesis, REP characterization, and predictive theory employed here demonstrates a powerful pathway to understand and ultimately tune exciton–phonon interactions for future quantum optical devices.

Janus monolayers↗

LDES-Sizing-and-siting-model (Power grid test cases for long-duration energy storage (LDES) siting) [SWR-25-75]

This repository contains code and data for performing a siting analysis of long-duration energy storage (LDES) in the 5-bus and RTS systems using the Sienna suite developed by the National Renewable Energy Laboratory for production cost modeling (PCM). This analysis involves moving the LDES component to different buses in the system, running a simulation, and considering the production cost of the simulation. Different system configurations are also analyzed with these scripts (such as moving load or renewable dispatch generators to different buses) to observe the impacts the system configuration has on optimal siting. This repository contains three sub directories discussed below. Both the 5-bus and RTS systems have two different initial configurations for the renewable energy components in them, one that is predominantly PV-driven and one that is predominantly wind-driven. Sienna suite can be found here: https://github.com/NREL-Sienna

Cole, David [University of Wisconsin]↗

Development of a Transferable Density-Functional Tight-Binding Model for Organic Molecules at the Water/Platinum Interface

A computationally efficient and transferable approach for modeling reactions at metal/water interfaces could significantly accelerate our understanding and ultimately the development of new catalytic transformations, particularly in the context of the emerging field of biomass conversion. Here, we present a parametrization of Pt–X (X = H, O, C) density-functional tight-binding (DFTB) for addressing this need. We first constructed Pt–H, Pt–O, and Pt–C repulsive potential splines. These pairwise parameters were then augmented to include many-body interactions using the Chebyshev Interaction Model for Efficient Simulation (ChIMES). We compare the geometrical and energetic performances of both DFTB and DFTB/ChIMES methods with DFT reference data across a variety of organic molecules at the platinum surface from nanoparticles to single-crystal surfaces. DFTB shows limited transferability between extended crystal surfaces and small nanoparticles. This transferability is significantly improved through the introduction of three-body interactions with Pt in DFTB/ChIMES, which provides consistent results across various systems, with reductions in the RMSD from around 30 kcal/mol in DFTB to around 10 kcal/mol. We demonstrate the stability and reliability of the obtained parameters by performing metadynamic simulations for the adsorption of phenol on Pt(111). We observe that DFTB itself is undersolvating the surface, leading to only one or two chemisorbed water molecules in a c(4 × 6) unit cell. In contrast, DFTB/ChIMES leads to a coverage of about 0.5 ML and successfully captures the chemisorbed mode of phenol at both the solid/liquid and the solid/gas interfaces. Furthermore, in agreement with experimental measurements, the adsorption at the solid/liquid interface is significantly weaker than that at the solid/gas interface. As a result, we highlight that even with DFTB, where we can accumulate dynamics for more than 1 ns for a given system, the simulations are not fully converged.

Adsorption↗

Phase-field modeling of thermally-grown oxide and damage evolution in environmental barrier coatings

Silicon carbide-based ceramic matrix composites protected by environmental barrier coatings (EBCs) present a promising materials solution for next-generation gas turbines. Developing more robust and efficient EBCs is therefore of significant technological importance. During the service in high-temperature oxidative environments, there is a thermally grown oxide (TGO) layer, spontaneously formed in the EBC system. TGO is recognized as a critical factor for the degradation and failure of EBCs, yet the detailed mechanisms of TGO growth and its effect on EBC failure remain unclear. In this study we develop a comprehensive chemo-mechano-phase-field model to simulate growth of the TGO in EBCs, factoring in creep and deformation, and especially the cracking behaviors. The volume expansion due to TGO growth and the resulting large inelastic deformation are addressed by using our recently developed, so-called incremental realization of inelastic deformation (IRID) algorithm, in combination with an adapted Hu-Chen spectral solver for elasticity. Simulations of TGO growth are performed considering different growth modes of TGOs determined mainly by the ratio of oxidant permeability in the topcoat to that in the TGO itself. Large-scale three-dimensional (3D) simulations are performed to model the formation of interconnecting vertical/channel cracks (often called ‘mud cracks’). The simulated crack morphology are in excellent agreement with the experimental observations from the literature. The simulations also provide insights into the cracking of EBCs and its dependence on the structure and constituent properties of the coating system. Furthermore, these results demonstrate the developed damage model can be a useful tool for design of more durable EBCs.

Cracking↗

An analytical model of “Electron-Only” magnetic reconnection rates

Abstract “Electron-only” reconnection, which is both uncoupled from the surrounding ions and much faster than standard reconnection, is arguably ubiquitous in turbulence. One critical step to understanding the rate in this novel regime is to model the outflow speed that limits the transport of the magnetic flux, which is super ion Alfvénic but significantly lower than the electron Alfvén speed based on the asymptotic reconnecting field. Here we develop a simple model to determine this limiting speed by taking into account the multiscale nature of reconnection, the Hall-mediated electron outflow speed, and the pressure buildup within the small system. The predicted scalings of rates and various key quantities compare well with fully kinetic simulations and can be useful for interpreting the observations of NASA’s Magnetospheric-Multiscale (MMS) mission and other ongoing missions.

Physics↗

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗

Lattice Effective Field Theory Simulations of Nuclei

Lattice effective field theory applies the principles of effective field theory in a lattice framework where space and time are discretized. Nucleons are placed on the lattice sites, and the interactions are tuned to replicate the observed features of the nuclear force. Monte Carlo simulations are then employed to predict the properties of nuclear few- and many-body systems. Here, we review the basic methods and several theoretical and algorithmic advances that have been used to further our understanding of atomic nuclei.

nuclear lattice effective field theory↗

The interplay of excitonic delocalization and vibrational localization in optical lineshapes: A variational polaron approach

The dynamics of molecular excitonic systems are complicated by a competition between electronic coupling (which drives delocalization) and vibrational-electronic (vibronic) interactions (which tend to encourage electronic localization). A particular challenge of molecular systems is that they typically possess a large number of independent vibrations, with frequencies often spanning the entire spectrum of relevant electronic energy gaps. Recent spectroscopic observations and numerical simulations on a water-soluble chlorophyll-binding protein (WSCP) reveal a transition between two regimes of vibronic behavior, a Redfield-like regime in which low-frequency vibrations respond to a delocalized excitonic state, and a Förster-like regime where high-frequency vibrations act as incoherent excitations on individual pigments. Although numerical simulations can reproduce these effects, there is a need for a simple, systematic theory that accurately describes the smooth transition between these two regimes in experimental spectra. Here we address this challenge by generalizing the variational polaron transform approach of [Bloemsma et al., Chem. Phys. 481, 250 (2016)] to include arbitrary bath densities for systems with or without symmetry. We benchmark this theory against both numerical matrix-diagonalization methods and experimental 77 K fluorescence spectra for two WSCP variants, obtaining quite satisfactory agreement in both cases. Here, we apply this theory to offer an explanation for the large loss in apparent electronic coupling in the WSCP Q57K mutant and to examine the likely impact of the interplay between excitonic delocalization and vibrational localization on vibrational sideband shapes and apparent coupling strengths in high-resolution optical spectra for chlorophyll-protein complexes such as WSCP.

14 SOLAR ENERGY↗

Impacts of sea ice leads on sea salt aerosols and atmospheric chemistry in the Arctic

Abstract. The processes contributing to Arctic cold-season (November–April) sea salt aerosols (SSAs) remain uncertain. Observations from coastal Alaska suggest that emissions from open leads in sea ice, which are not included in climate models, may play a dominant role. Their Arctic-wide significance has not yet been quantified. Here, we create an emission parameterization of SSAs from leads by combining satellite data of lead area (the Advanced Microwave Scanning Radiometer–Earth Observation System (AMSR-E) product) and a chemical transport model (GEOS-Chem) to quantify pan-Arctic SSA emissions from leads during the cold season from 2002 to 2008 and to predict their impacts on atmospheric chemistry, evaluating the results of our simulated SSAs against in situ observations. The AMSR-E product detects large leads with certainty (> 3 km in size), and, hence, our study is limited to quantifying emissions from large leads. Lead emissions vary seasonally and interannually. Simulated total monthly SSA emissions increase by 1.1 %–1.8 % (≥60° N latitude) and 5.6 %–7.5 % (≥75° N) for the 2002–2008 cold seasons. SSA concentrations primarily increase at the location of leads, where standard model concentrations are low. GEOS-Chem overestimates SSA concentrations at Arctic sites compared to ground observations, even when lead emissions are not included, suggesting underestimation of SSA sinks and/or uncertainties in SSA emissions from blowing snow and the open ocean. Multi-year monthly mean surface bromine atom (Br) concentrations increase by 2.8 %–8.8 % due to SSAs from leads for the 2002–2008 cold seasons. Changes in ozone concentrations are negligible. While leads contribute < 10 % to Arctic-wide SSA emissions in the years 2002–2008, these emissions occur in regions of low background aerosol concentrations. Leads may increase in frequency under future climate change, which could increase SSA emissions from leads.

Emme, Erin J. (ORCID:0000000268406352)↗

CoLoRe-2LPT: Lyman-$α$ mock catalogues for the validation of DESI cosmological analyses

The Lyman-$α$ (Ly$α$) forest has become a crucial probe for studying the large-scale structure of the universe at high redshift ($z > 2$), providing powerful constraints on Baryon Acoustic Oscillations (BAO) and the full-shape (FS) clustering of matter. As a key ingredient for upcoming BAO and FS analyses, we present a new generation of fast cosmological Ly$α$ mocks based on second-order Lagrangian perturbation theory (2LPT). These new mocks significantly improve upon previous log-normal approaches, both at accurately capturing small scale clustering and at recovering the non-linear broadening of the BAO peak. They are able to reproduce Ly$α$ statistics within $10\%$ of the latest DESI measurement; including the Ly$α$ bias and the redshift-space distortion $β$ parameter, mean transmitted flux, and 1D power spectrum. The corresponding quasar (QSO) clustering is also improved with respect to previous approaches, calibrated against high-resolution Abacus simulations, recovering the observational QSO linear bias to less than $5\%$ and improving redshift-space distortions via 2LPT velocities and the addition of Fingers-of-God effects. Furthermore, these mocks incorporate high column density systems and metal lines, allowing us to explore the effects and systematics induced by these astrophysical contaminants. This new set of mocks has been key for enhancing the modeling and validation of the DESI DR2 Ly$α$ full shape cosmological analysis. This work provides a physically motivated and computationally efficient tool for simulating current and next-generation Ly$α$ surveys and validating FS and BAO analysis.

Bernal, M.F. Ruiz-Herrera [Madrid, CIEMAT] (ORCID:↗

Fractional AC Josephson effect in a topological insulator proximitized by a self-formed superconductor

A lateral Josephson junction in which the surface of a 3D topological insulator serves as the weak link should support topologically protected excitations related to Majorana fermions. The resulting 4π-periodic current-phase relationship could be detected under high-frequency excitation by the suppression of odd Shapiro steps. Here, in this study, we demonstrate such devices through the self-formation of a Pd-Te superconducting layer from a telluride topological insulator, and observe suppressed first and third Shapiro steps. Other devices, including those where the Pd-Te layer is bolstered by an additional Al layer, show no suppression of Shapiro steps, a difference supported by simulations. Though we rule out the known trivial causes of suppressed Shapiro steps in our devices, we nevertheless argue that corroborating measurements and disorder-aware theoretical descriptions of these systems are needed before confidently claiming the observation of Majorana states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Influences of lidar scanning parameters on wind turbine wake retrievals in complex terrain

Abstract. Scanning lidars enable the collection of spatially distributed measurements of turbine wakes and the estimation of wake properties such as magnitude, extent, and trajectory. Lidar-based characterizations, however, may be subject to distortions due to the observational system. Distortions can arise from the resolution of the measurement points across the wake, the projection of the winds onto the beam, averaging along the beam probe volume, and intervening evolution of the flow over the scan duration. Using a large-eddy simulation and simulated measurements with a virtual lidar model, we assess how scanning lidar systems may influence the properties of the retrieved wake using a case study from the Perdigão campaign. We consider three lidars performing range-height indicator sweeps in complex terrain, based on the deployments of lidars from the Danish Technical University (DTU) and German Aerospace Center (DLR) at the Perdigão site. The unwaked flow, measured by the DTU lidar, is well-captured by the lidar, even without combining data into a multi-lidar retrieval. The two DLR lidars measure a waked transect from different downwind vantage points. In the region of the wake, the observation system reacts to the smaller spatial and temporal variations of the winds, allowing more significant observation distortions to arise. While the measurements largely capture the wake structure and trajectory over its 4–5 D extent, limited spatial resolution of measurement points and volume averaging lead to a quicker loss of the two lobes in the near wake, smearing of the vertical bounds of the wake (< 30 m), wake center displacements up to 10 m, and dampening of the maximum velocity deficit by up to a third. The virtual lidar tool, coupled with simulations, provides a means for assessing measurement capabilities in advance of measurement campaigns.

17 WIND ENERGY↗

Adsorption Hysteresis Under Control: Tuning Host–Guest Interactions via a Genetic Algorithm

Mesoporous adsorbent materials offer a large volumetric capacity; however, cyclic adsorption/desorption processes in these systems often suffer from hysteresis and may require a significant pressure swing to access this capacity. To mitigate hysteresis, a proposed strategy is to include nucleation sites on the walls of the mesoporous material to facilitate droplet and bubble formation, lowering the free energy barriers to the respective phase transitions. It is unclear, however, what combination of adsorbate− adsorbent interactions and spatial patterning would be beneficial for a given application, considering that improvements to some sorption properties may come at the expense of other attributes. To understand these interconnected observables, we examine two model systems, planar-slit and cylindrical pores with tunable interaction sites, using GPU-accelerated transition matrix Monte Carlo simulations. The simulations provide a free energy map of the pressure−adsorption space in a matter of minutes, which we use to track adsorption isotherm characteristics as a function of adsorbent properties. We then leverage the rapid acquisition of simulation data to construct a genetic algorithm to iteratively modify interaction sites of the slit-pore wall to minimize the hysteresis of this system without sacrificing uptake. We find that the adsorption branch of the isotherm is easily modulated via the average host−guest interaction strength, but desorption is only adjustable if there is a suitable bubble nucleation site. Within the context of a slit-pore system, we identify relative interaction strengths and patch sizes required to gain control over both branches of the hysteresis loop.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron

Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. Here, we further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.

Ab-initio molecular dynamics↗

Demography, dynamics and data: building confidence for simulating changes in the world's forests

Vegetation demographic models (VDMs) are advanced tools for simulating forest responses to climate and land-use changes, and are essential for projecting carbon cycling and large-scale forest management strategies. Despite their increasing incorporation into Earth System Models, VDMs differ in their demographic assumptions, with no prior quantitative comparison of their performance. We benchmarked nine VDMs against observational data from boreal, temperate and tropical sites, assessing their accuracy in predicting tree growth, carbon turnover, biomass stocks and size distributions. Models were simulated under consistent climate conditions with postdisturbance recovery monitored for at least 420 yr. Postdisturbance carbon recovery trajectories showed significant variability while remaining within observational ranges. Initial regrowth rates varied substantially (0.03-0.60, 0.18-0.70 and 0.35-1.10 kgCm-2 yr-1 for boreal, temperate and tropical sites, respectively), influenced by each model's initial forest state. Models captured mature forest carbon content but showed compensating effects between overestimated growth and underestimated mortality rates. This first multi-model benchmarking identifies growth and mortality rates as critical calibration targets and highlights the need to refine postdisturbance establishment conditions for model development. We outline specific benchmarking variables needed to improve predictions of forest responses to environmental change.

demographic vegetation model benchmarking↗