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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A Comprehensive Strategy for Grid Forming Control in DC Coupled Photovoltaic and Battery Energy Storage Inverters

This paper presents an integrated DC-DC and DCAC grid-forming control strategy for DC-coupled photovoltaic (PV) plus battery energy storage systems, considering the effect of DC link voltage variations caused by direct PV connections. A power reference algorithm determines power distribution between the PV and battery to the grid while observing device power ratings to prevent the over-rating of components and keep the battery's state of charge within an acceptable range. The simulated utility-scale model in MATLAB/Simulink illustrates its ability against extreme phase angle variation contingencies in the grid while controlled through grid-forming control with a fast dynamic on DC link voltage. The simulation results confirm the effectiveness of the proposed control in integrating PV plus battery configurations with grid forming control and maintaining reliable grid operation under severe grid disturbances.

battery, boost, control, energy storage, grid form↗
DOE OSTI2025

Lens Modeling of STRIDES Strongly Lensed Quasars Using Neural Posterior Estimation

Strongly lensed quasars can be used to constrain cosmological parameters through time-delay cosmography. Models of the lens masses are a necessary component of this analysis. To enable time-delay cosmography from a sample of $\mathcal{O}(10^3)$ lenses, which will soon become available from surveys like the Rubin Observatory’s Legacy Survey of Space and Time and the Euclid Wide Survey, we require fast and standardizable modeling techniques. To address this need, we apply neural posterior estimation (NPE) for modeling galaxy-scale strongly lensed quasars from the Strong Lensing Insights into the Dark Energy Survey (STRIDES) sample. NPE brings two advantages: speed and the ability to implicitly marginalize over nuisance parameters. We extend this method by employing sequential NPE to increase precision of mass model posteriors. We then fold individual lens models into a hierarchical Bayesian inference to recover the population distribution of lens mass parameters, accounting for out-of-distribution shift. After verifying our method using simulated analogs of the STRIDES lens sample, we apply our method to 14 Hubble Space Telescope single-filter observations. We find the population mean of the power-law elliptical mass distribution slope, γ lens , to be $\mathcal{M}_γ$ lens = 2.13 ± 0.06. Our result represents the first population-level constraint for these systems. This population-level inference from fully automated modeling is an important stepping stone toward cosmological inference with large samples of strongly lensed quasars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Impact of Reordering on the LU Factorization Performance of Bordered Block-Diagonal Sparse Matrix

Power engineers rely on computer-based simulation tools to assess grid performance and ensure security. At the core of these tools are solvers for sparse linear equations. When transformed into a bordered block-diagonal (BBD) structure, part of the sparse linear equation solving can be parallelized. This work focuses on using the Schur-complement-based method for LU factorization on BBD matrices, specifically, Jacobian matrices from large-scale systems. Our findings show that the natural ordering method outperforms the default ordering method in computational performance for each block of the BBD matrix. This observation is validated using synthetic 25k-bus and 70k-bus cases, showing a speedup of up to 38% when using natural ordering without permutation. Additionally, the impact of the number of partitions is studied, and the result shows that computational performance improves with more, smaller partitions in the BBD matrices.

BBD matrix↗
DOE OSTI2026

Solar-Tracker Diffuse-Response Algorithm for Balancing Energy Gain and Mechanical Wear in Photovoltaic Systems

Single-axis solar tracking maximizes photovoltaic energy production under clear-sky conditions; however, its effectiveness decreases under cloudy and overcast skies, where diffuse irradiance dominates and the optimal module orientation changes. Conventional tracking algorithms either neglect sky conditions or rely on simplified diffuse-response strategies that may trigger frequent tracker repositioning under variable cloud cover, leading to increased mechanical wear with marginal energy gains. This work proposes an enhanced diffuse-response tracking algorithm that explicitly accounts for both the intensity and temporal persistence of cloudiness. By requiring overcast conditions to persist for a minimum duration before reorienting the tracker to a diffuse-stow position, the proposed approach reduces unnecessary movements while preserving the benefits of diffuse-response operation. The algorithm is evaluated through numerical simulations based on historical meteorological data and validated using field measurements on monofacial and bifacial photovoltaic strings. The results show that the proposed strategy reduces excess tracker movement from 114% to 0.16% while maintaining nearly the same energy yield. Compared to a conventional diffuse-response algorithm, the associated energy reduction is minimal (≈0.17%) relative to the ≈0.37% yield gain observed at the studied location. These findings demonstrate that incorporating cloudiness duration enables a practical compromise between energy performance and tracker durability, particularly for monofacial photovoltaic systems.

14 SOLAR ENERGY↗
DOE OSTI2024

Impact of volcanic eruptions on CMIP6 decadal predictions: a multi-model analysis

Abstract. In recent decades, three major volcanic eruptions of different intensity have occurred (Mount Agung in 1963, El Chichón in 1982 and Mount Pinatubo in 1991), with reported climate impacts on seasonal to decadal timescales that could have been potentially predicted with accurate and timely estimates of the associated stratospheric aerosol loads. The Decadal Climate Prediction Project component C (DCPP-C) includes a protocol to investigate the impact of volcanic aerosols on the climate experienced during the years that followed those eruptions through the use of decadal predictions. The interest of conducting this exercise with climate predictions is that, thanks to the initialisation, they start from the observed climate conditions at the time of the eruptions, which helps to disentangle the climatic changes due to the initial conditions and internal variability from the volcanic forcing. The protocol consists of repeating the retrospective predictions that are initialised just before the last three major volcanic eruptions but without the inclusion of their volcanic forcing, which are then compared with the baseline predictions to disentangle the simulated volcanic effects upon climate. We present the results from six Coupled Model Intercomparison Project Phase 6 (CMIP6) decadal prediction systems. These systems show strong agreement in predicting the well-known post-volcanic radiative effects following the three eruptions, which induce a long-lasting cooling in the ocean. Furthermore, the multi-model multi-eruption composite is consistent with previous work reporting an acceleration of the Northern Hemisphere polar vortex and the development of El Niño conditions the first year after the eruption, followed by a strengthening of the Atlantic Meridional Overturning Circulation the subsequent years. Our analysis reveals that all these dynamical responses are both model- and eruption-dependent. A novel aspect of this study is that we also assess whether the volcanic forcing improves the realism of the predictions. Comparing the predicted surface temperature anomalies in the two sets of hindcasts (with and without volcanic forcing) with observations we show that, overall, including the volcanic forcing results in better predictions. The volcanic forcing is found to be particularly relevant for reproducing the observed sea surface temperature (SST) variability in the North Atlantic Ocean following the 1991 eruption of Pinatubo.

Bilbao, Roberto (ORCID:0000000307294980)↗
DOE OSTI2025

Extracted Water Induces Concentration Fluctuations in Model Ternary Liquid–Liquid Extraction System

Organic phase aggregation and phase splitting are two important (and connected) phenomena in liquid–liquid extraction of metal ions. Previously, we demonstrated how structure in binary mixtures of extractant and diluent originates from concentration fluctuations associated with the so-called “third-phase formation” liquid–liquid phase instability. For more complex systems, we sought to understand how extracted aqueous components affect these concentration fluctuations. Here, in this study, we use small-angle X-ray scattering and molecular dynamics simulations to investigate organic phase structural changes upon contact with water in a model ternary system of water/alkane/extractant. Critical exponents were determined from scattering with temperature variation and were consistent with mean field values. Then, by controlling the water content in the organic phase through direct dissolution of fixed quantities at constant temperature, we relate the observed increase in the fluctuation length scale to changes to the spinodal temperature of the third-phase formation phase boundary. We find a roughly linear relationship between water concentration and the spinodal temperature and a similar slope in that relationship for different linear alkane diluents. This suggests that the mechanism by which water extraction impacts the organic phase structure is to enhance concentration fluctuations through the introduction of additional polar extractant-water complexes, which increases the spinodal temperature and thereby decreases the reduced temperature, i.e., the distance to the instability. This explanation for how solutes induce aggregation in liquid–liquid extraction organic phases may extend to more complex systems, explaining the large increases in correlation lengths reported under high acid and metal loading.

Bonnett, Brittany L. [Argonne National Laboratory ↗
DOE OSTI2026

Disentangling cation effects on ion mobility and structure in ionic liquid electrolytes

Ionic liquids (ILs) are low-temperature molten salts, where ion transport is primarily governed by ion–ion interactions. Yet, the impact of organic IL cations on critical electrolyte properties such as ion dissociation and overall transport behavior in lithium-salt-doped ILs remains poorly understood. Moreover, despite their critical role in designing IL-based electrolytes for energy storage applications, ion–ion interactions and ion-specific transport under an applied electrical potential are seldom quantified, largely due to the unique experimental and computational challenges involved. Herein, we compare transport properties obtained using 1 H, 7 Li, and 19 F pulsed-field gradient nuclear magnetic resonance (NMR) and electrophoretic NMR (eNMR) with those measured by electrochemical impedance spectroscopy. Non-equilibrium molecular dynamics (MD) simulations and eNMR confirm the presence of negatively charged [Li(TFSI) n ] (1−n) aggregates that migrate towards the positive electrode, resulting in negative lithium transference numbers. Equilibrium MD simulations reveal a vehicular Li ion transport mechanism facilitated by long-lived aggregates with Li + cations strongly bound to multiple TFSI − anions. Finally, we observe an inverse relationship between the apparent charge of the TFSI − anion in the neat IL, which is dictated by the IL cation, and Li + transport in the salt-doped systems. This highlights the opportunity to tune electrolyte performance by tailoring cation chemistry.

Li-ion batteries↗
DOE OSTI2024

High-Entropy Alloy R&D for Accelerator Beam Window Applications

High-Entropy Alloys are a class of novel material that can offer improved resistance to beam-induced radiation damage and thermal shock. Development of these new alloys to serve as beam windows in multi-megawatt accelerator target applications is ongoing at Fermilab. Currently we are investigating AlCoCrMnTiV alloy systems of 4-6 components for service as beam windows; these compositions are predicted by CALPHAD simulation to have a low density and single-phase BCC crystal structure. The microstructures of these systems are being studied by electron microscopy techniques such as energy dispersive X-ray spectroscopy (EDS) to determine elemental homogeneity and composition, electron backscatter diffraction (EBSD) to quantify grain structure and orientation, and transmission electron microscopy (TEM) to observe defect structures and precipitate formation; nanoindentation is used to probe microstructural mechanical properties. Initial bulk property characterization utilizes differential scanning calorimetry to quantify specific heat capacity (cp), dilatometry to determine the coefficient of thermal expansion (CTE), and time-domain thermoreflectance (TDTR) to measure thermal conductivity (K). A miniature tensile testing apparatus is also being developed to test tensile properties. Post-irradiation examination is currently ongoing for several of these compositions that have been irradiated by low-energy ions to damage levels and at a temperature relevant to beam window applications at future next-generation accelerator facilities. This talk will briefly describe the alloy design and synthesis before going into more depth covering microstructural pre-characterization, and post-irradiation examination results from low-energy ion irradiated specimens. This will be followed by the plans for alloy down selection and future prototypic proton irradiations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

The Role of Deep Convection and Large-scale Circulation in Driving Model Spread in Low Cloud Feedback and Equilibrium Climate Sensitivity

This project aims to advance the understanding of the processes that drive the large uncertainties in climate change projections, use observations to constrain model physics and reduce the inter-model spread in equilibrium climate sensitivity (ECS). There are three major goals: 1) Characterize the representation of the physical pathways that link deep convection, large-scale circulation and low cloud feedback in CMIP6 model simulations and determine the relative contribution of each pathway to the CMIP6 model spread in low cloud feedback and ECS; 2) Use process-oriented diagnostics and multiple observations to evaluate CMIP6 model performance in capturing the observed cloud-circulation relation and deep convection characteristics including convective transition statistics and the bulk properties of mesoscale convective systems (MCSs). Error decomposition in CMIP6 models will be performed. 3) Conduct E3SM short-range hindcasts following the DOE Cloud-Associated Parameterizations Testbed (CAPT) protocol to pinpoint specific model parameters/processes that are crucial to the representation of deep convection, circulation, clouds and the pathways that connect them. We will modify convective parameters in E3SM and analyze the perturbed physics experiments (PPEs) to isolate model parameters that are critical to the uncertainty of ECS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Feasibility of fusion plasma burn control via real-time, sub-divertor neutral gas isotopic and compositional analysis

The ability to provide fusion burn control without requiring physical access through the first wall and fuel breeding blankets, would be vital for any future, magnetically confined fusion power reactor. A multi-sensor, fusion fuel cycle exhaust, neutral gas analysis system on JET, capable of delivering real time data, and accessing only the sub-divertor region, provides an excellent example of such capability. Optimized for and operated during the deuterium–tritium experimental campaigns 2 and 3 (DTE2, DTE3), it is proving valuable for planning to explore fusion reactor burn control in ITER with a comparable diagnostic system called the Diagnostic Residual Gas Analyzer (DRGA). This paper aims to show feasibility of developing model-based controllers for ITER and next generation, reactor-relevant devices, by building both on the empirical experience in JET-DTE2, and on the already emerging experience on developing such models specifically for ITER. The paper begins with a specific use-case from JET-DTE2, pertaining to the observed sensitivity of the fusion neutron yield on the concentration of isotopic helium-3 ( 3 He), with data from one of the high-performance DT shots exhibited with emphasis on the 3 He measurement via the sub-divertor. Then, a first model is developed and then explored with simulations that aim to discover how well the controllers in the model react to either insufficient levels of 3 He or excessive levels of 3 He. The simulations then explore potential impact from a delay in the measurement (or the response) that would be comparable to the ∼1 s, conductance limited response for the ITER DRGA system, currently in its final design. The simulations show that control is feasible, and that its effectiveness is not significantly impacted by such delay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

COMPASS-FME Synoptic Sites Level 2 Sensor Data v2-1

This is the version 2-1 Level 2 (L2) data release for COMPASS-FME environmental sensors located at our synoptic field sites. COMPASS-FME is studying sites in two distinct regions, the Chesapeake Bay and the Western Lake Erie Basin. We established the network at seven "synoptic" (observational) sites along the Chesapeake Bay and Lake Erie coastlines, collectively generating over three million observations per month, to track and comprehend environmental changes where land and water intersect. Additionally, the two regions provide an interesting contrast of saltwater and freshwater coasts that allow us to differentiate the impacts of inundation and coastal water chemistries in two nationally important coastal systems. Level 2 (L2) data consist of sensor observations from the COMPASS-FME synoptic sites, TEMPEST, and DELUGE. Compared to the L1 data, these are more consistent (always 15-minute timestamps for the entire year); better QA/QC’d (out of bounds, out of service, and extreme outlier values are removed); and more complete, with a gap-filled time series available alongside the main observations, and additional derived (calculated) variables. L2 data are intended to be rapidly and easily usable in analyses and simulations. However, algorithmic outlier identification always carries the risk of removing valid data, and Level 1 data may be more suitable for analyses that focus on variability or extreme events. This dataset includes: - An overall dataset README file that describes the current version, gives citation and contact information, etc. - Site- and year-specific folders, each holding variable-specific Parquet (a high performance, space efficient format; see https://parquet.apache.org) data files for each site and plot in that year. - Metadata files within each site-year folder provide full information on data units, expected ranges, contact information, detailed flood times, as well as a general description of the site. - Environmental sensor types that appear in the data files include weather (ClimaVUE50, CS, RM Young, and LI instruments in the graphs below); soil conditions (TEROS12); soil redox state (Redox); groundwater variables (AquaTROLL200 and AquaTROLL600); open water sondes (Exo); tree sap velocity (Sapflow); and system voltage and state (Datalogger). Data are reported every 15 minutes. Data files are in Apache Parquet, a high performance, space efficient format for tabular data. These files can be read using R's `arrow` package (https://arrow.apache.org/docs/r/), with similar tools available in other languages. Please see v2-1 L2 Sensor Package QStart.pdf for detailed information on data package structure, temporal coverage, and versioning.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗
DOE OSTI2025

Liquid crystal torons in Poiseuille-like flows

Three-dimensional (3D) simulations of the structure of liquid crystal (LC) torons, topologically protected distortions of the LC director field, under material flows are rare but essential in microfluidic applications. Here, we show that torons adopt a steady-state configuration at low flow velocity before disintegrating at higher velocities, in line with experimental results. Furthermore, we show that under partial slip conditions at the boundaries, the flow induces a reversible elongation of the torons, also consistent with the experimental observations. These results are in contrast with previous simulation results for 2D skyrmions under similar flow conditions, highlighting the need for a 3D description of this LC soliton in relation to its coupling to the material flow. These findings pave the way for future studies of other topological solitons, like hopfions and heliknotons, in flowing soft matter systems.

Science & Technology - Other Topics↗
DOE OSTI2024

Low Temperature CO2 Capture from an NGCC Flue Gas Using a Magnetically Stabilized, Inductively Heated Fluidized Bed Reactor

The lab sorbent testing finalized the sorbent composition and operating conditions for the integrated unit. The bench unit design included COMSOL simulation to finalize the coil design and select initial operating conditions of the power supply system. Fluidized bed testing with the induction heating verified that a heating efficiency of up to 90% can be achieved. The integrated bench unit operations provided valuable data confirming the feasibility of using induction heating for the desorption of CO2. Major challenge observed was the limited working capacity of the sorbent, requiring additional development in future projects. In addition, a large-scale reactor demonstration is required to mitigate the risk of uniform temperature distribution and heating efficiency over a large cross-sectional area when using induction heating.

Tong, Andrew↗
DOE OSTI2024

MatCal Users Guide: Release 1.3.0

Any continuum mechanics model will require three components: (1) a discretized geometry of the boundary value problem being studied, (2) the partial differential equations to be solved, and (3) the initial conditions and boundary conditions for the problem. To describe material behavior in these computational models, material models contribute to (2) the underlying equations and, occasionally, to (3) the initial conditions for the simulation. These material models can exhibit a mathematical form that is empirically based, based on first principles, or developed from both empirical observations and known physics. In general, these models are meant to represent a class of materials with well understood behavior. As a result, material models have parameters that must be tuned or calibrated so that the model response matches characterization data available for the specific material it is intended to represent when used to simulate a specific system. For simple models, such as isotropic, linear elastic materials in solid mechanics, this calibration process can be a simple analytical calculation directly extracting the parameters from experimental measurements. For complex models that have many inputs and require many characterization datasets to adequately identify the material behavior, the model calibration process can require an inverse problem approach where an optimization is performed to tune the model parameters to the available data.

36 MATERIALS SCIENCE↗
DOE OSTI2024

HighNESS conceptual design report: Volume II. The NNBAR experiment.

A key aim of the HighNESS project for the European Spallation Source is to enable cutting-edge particle physics experiments. This volume presents a conceptual design report for the NNBAR experiment. NNBAR would exploit a new cold lower moderator to make the first search in over thirty years for free neutrons converting to anti-neutrons. The observation of such a baryon-number-violating signature would be of fundamental significance and tackle open questions in modern physics, including the origin of the matter-antimatter asymmetry. This report shows the design of the beamline, supermirror focusing system, magnetic and radiation shielding, and anti-neutron detector necessary for the experiment. A range of simulation programs are employed to quantify the performance of the experiment and show how background can be suppressed. For a search with full background suppression, a sensitivity improvement of three orders of magnitude is expected, as compared with the previous search. Civil engineering studies for the NNBAR beamline are also shown, as is a costing model for the experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗
DOE OSTI2025

Interactions between land use, fires, and dust as drivers of global climate change (Final Technical Report)

This award to UC Irvine and Cornell University led to the publication of 18 papers, three of which are highly cited. According to ISI Web of Science, the body of work from this grant has been cited over 476 times by May 2025. Several additional publications remain in review. This grant contributed to several significant science advances, including an improved understanding of the influence of fire aerosols on global photosynthesis (Xu et al., 2021), dust-driven teleconnections between Africa and South America (Li et al., 2021), the role of fire aerosols in triggering phytoplankton blooms in the Arctic (Ardyna et al., 2022), and wildfire responses to heat waves (Li et al., in review). Several other studies are being written up for publication, and will acknowledge funding support. These efforts include an analysis in which the magnitude of tropical fire emissions from land use change is constrained by comparing E3SMv3 simulations with atmospheric observations. This work has revealed that carbon cycle model estimates of land use change may considerably overestimate the global tropical deforestation flux. The project also contributed to several technical advances, including a new process-based and scale-aware desert dust emission scheme for global climate models (Leung et al., 2023), an improved representation of dust and its influence on radiation in E3SM (Feng et al., 2022), an improved representation of the dust parameterization in the Community Earth System Model (Li et al., 2022), and an improved coupling of wildfires with atmospheric chemistry in E3SM v3.

58 GEOSCIENCES↗
DOE OSTI2025

Size-resolved process understanding of stratospheric sulfate aerosol following the Pinatubo eruption

Stratospheric sulfate aerosol produced by volcanic eruptions plays important roles in atmospheric chemistry and the global radiative balance of the atmosphere. The simulation of stratospheric sulfate concentrations and optical properties is highly dependent on the chemistry scheme and microphysical treatment. In this work, we implemented a sophisticated gas-phase chemistry scheme (full chemistry, FC) and a 5-mode version of the Modal Aerosol Module with Prognostic Stratospheric Aerosol (MAM5-PSA) for the interactive treatment of stratospheric sulfate aerosol in the Department of Energy's Energy Exascale Earth System Model version 2 (E3SMv2) model to better simulate the chemistry-aerosol feedback following the Pinatubo eruption, and to compare it against a simulation using simplified chemistry (SC) and the default 4-mode version of the Modal Aerosol Module (MAM4). MAM5-PSA experiments were found to better capture the stratospheric sulfate burden from the eruption of the volcano to the end of 1992 as compared to the High-resolution Infrared Sounder (HIRS) observations, and the formation of sulfate in MAM5-PSA with FC (with an additional OH replenishment reaction) was significantly faster than in MAM4 with FC. Analyses of microphysical processes indicate that more sulfate aerosol mass was generated in total in FC experiments than in SC experiments. MAM5-PSA performs better than MAM4 in simulation of aerosol optical depth (AOD); AOD anomalies from the MAM5-PSA experiment have better agreement with observations. The simulated largest changes in global mean net radiative flux at the top of the atmosphere following the eruption were about −3 W m −2 in MAM5-PSA experiments and roughly −1.5 W m −2 in MAM4 experiments.

AEROSOL↗