SEARCH · Engineering Papers
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.
In this work, we propose MD Loader, a hybrid in-memory data loader for distributed deep neural networks. MDLoader introduces a model-driven performance estimator to automatically switch between one-sided and collective communication at runtime.
Inverse problems constrained by partial differential equations (PDEs) play a critical role in model development and calibration. In many applications, there are multiple uncertain parameters in a model that must be estimated. However, high dimensionality of the parameters and computational complexity of the PDE solves make such problems challenging. A common approach is to reduce the dimension by fixing some parameters (which we will call auxiliary parameters) to a best estimate and use techniques from PDE-constrained optimization to estimate the other parameters. In this article, hyper-differential sensitivity analysis (HDSA) is used to assess the sensitivity of the solution of the PDE-constrained optimization problem to changes in the auxiliary parameters. Foundational assumptions for HDSA require satisfaction of the optimality conditions which are not always practically feasible as a result of ill-posedness in the inverse problem. Here we introduce novel theoretical and computational approaches to justify and enable HDSA for ill-posed inverse problems by projecting the sensitivities on likelihood informed subspaces and defining a posteriori updates. Our proposed framework is demonstrated on a nonlinear multiphysics inverse problem motivated by estimation of spatially heterogeneous material properties in the presence of spatially distributed parametric modeling uncertainties.
This project demonstrates, for the first time, the application of simulation-based inference (SBI) techniques to tune neutrino–nucleus interaction models. Using a mock dataset based on the MicroBooNE tuning of the GENIE event generator, our approach employs a Neural Posterior Estimator (NPE) with Masked Autoregressive Flows (MAF) to infer the posterior distributions of key GENIE parameters directly from simulated histograms. The workflow provides a scalable and amortized framework for performing likelihood-free inference in high-dimensional parameter spaces, offering a pathway to more efficient and uncertainty-aware model tuning for next-generation neutrino experiments such as DUNE and SBND.
Inorganic membrane-based systems are a promising technology for precombustion CO2 capture with simultaneous H2 production. State-of-the-art packages for high temperature and pressure service consist of multiple tube membrane bundles prepared in a "candle filter" configuration, in which the membrane tubes are open at one end and sealed at the other. This configuration is used for practical reasons, specifically the need to minimize problems due to thermal expansion mismatch between the ceramic tube bundle and the steel housing. However, the primary technical problem with the candle filter format for commercial-scale installations is the inability to purge the tube side (typically the permeate side), a feature that is crucial for high H2 recovery. In this study, the focus is to design and fabricate the first dual-end open full ceramic multiple tube membrane bundle that enables tube side (permeate) purge for gas separation applications. An additional key feature of this design is the simplified module layout, as the membrane bundles can be installed end-to-end with tube side (permeate) flow directly from one bundle to the next. This layout simplifies the membrane to housing seals and yields significant improvement in membrane packing density. Detailed focus areas in our studies include: (i) Materials development and preparation of the tube-to-tube sheet potting for the dual-ended bundle; (ii) the sealing and optimal module configuration design to minimize membrane stress upon module mounting; (iii) the demonstration, via the fabrication of CMS and Pd-alloy membranes supported on full-size, dual-ended ceramic support bundles, of the first example of a purgeable ceramic membrane and module; and (iv) development of a CFD model of the membrane module for calculation of feed flow distribution, and for use in scale-up, and capital cost estimating. The CFD model was validated using experimental data with the multi-tubular membrane system, employing He/N2 as a model gas mixture (surrogate for H2/CO2), and has been shown to be quite accurate. Employing the model, we are able to study the effects of operating pressure and temperature, feed and sweep gas flow rates, and the choice of membrane tube configuration on system performance.
The purpose of this study is to reduce hematite (Fe2O3) using hydrogen (H2), carbon monoxide (CO), and methane (CH4) gases using microwave irradiation. Conventional temperature programmed reduction of iron phase transformation under H2, CO, and CH4 agree with the literature. Furthermore, the activation energy for the iron phase transformation in a conventional reactor was in the order of H2>CO>CH4, suggesting that CH4 was the best reductant to initiate the reduction at a much lower temperature than CO, and CO was significantly better than H2 atmosphere. Due to the uneven temperature distribution, the apparent activation energy of microwave reduction under H2 was approximately 1/4th when compared to the conventional reduction in H2. X-ray diffraction (XRD) analysis showed a mixed oxide phase in microwave reduction as opposed to the clear phase transformation in conventional studies. To combat the temperature measurement issues in microwave studies, fiber optic sensors and Forward Looking InfraRed (FLIR) sensors were used to acquire the actual temperature distribution in the catalyst bed. Therefore, a reasonably accurate activation energy estimation was achieved in the catalyst bed. Combining these with the characterization and quantification of the different iron phases in microwave reduction will provide a blueprint for future microwave-assisted gasification studies.
Inorganic membrane-based systems are a promising technology for precombustion CO2 capture with simultaneous H2 production. State-of-the-art packages for high temperature and pressure service consist of multiple tube membrane bundles prepared in a "candle filter" configuration, in which the membrane tubes are open at one end and sealed at the other. This configuration is used for practical reasons, specifically the need to minimize problems due to thermal expansion mismatch between the ceramic tube bundle and the steel housing. However, the primary technical problem with the candle filter format for commercial-scale installations is the inability to purge the tube side (typically the permeate side), a feature that is crucial for high H2 recovery. In this study, the focus is to design and fabricate the first dual-end open full ceramic multiple tube membrane bundle that enables tube side (permeate) purge for gas separation applications. An additional key feature of this design is the simplified module layout, as the membrane bundles can be installed end-to-end with tube side (permeate) flow directly from one bundle to the next. This layout simplifies the membrane to housing seals and yields significant improvement in membrane packing density. Detailed focus areas in our studies include: (i) Materials development and preparation of the tube-to-tube sheet potting for the dual-ended bundle; (ii) the sealing and optimal module configuration design to minimize membrane stress upon module mounting; (iii) the demonstration, via the fabrication of CMS and Pd-alloy membranes supported on full-size, dual-ended ceramic support bundles, of the first example of a purgeable ceramic membrane and module; and (iv) development of a CFD model of the membrane module for calculation of feed flow distribution, and for use in scale-up, and capital cost estimating. The CFD model was validated using experimental data with the multi-tubular membrane system, employing He/N2 as a model gas mixture (surrogate for H2/CO2), and has been shown to be quite accurate. Employing the model, we are able to study the effects of operating pressure and temperature, feed and sweep gas flow rates, and the choice of membrane tube configuration on system performance.
The aerosol chemical speciation monitor (ACSM) was developed to adapt the technology of the aerosol mass spectrometer (AMS) to routine, long-term, standalone monitoring. The calculation of particulate mass concentration from ACSM data requires the measurement of the response of the instrument to aerosol of specific size and composition as well as assumptions about the instrument response based on laboratory measurements and field experience acquired over more than two decades of operation of AMS and a decade of operation of the ACSM. Three parameters in the concentration calculations that are particularly important are the NO3 response factor (RFNO3), relative ionization efficiency (RIE), and the collection efficiency (CE). The values of RFNO3 and RIEs are determined from calibration, however the jump calibration method previously used in calibration can result in errors in the RIE for sulfate. This has been corrected by implementing a continuous calibration method. The default collection efficiency is 0.5. This has been shown to result in mass loadings that do not agree with mass determinations from other instruments because of effects of composition on the vaporization of the particles. The previous work of investigators addressing this issue is discussed. After preliminary work on ACSM and scanning mobility particle sizer (SMPS) data from the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility Southern Great Plains (SGP) observatory collected in late 2016 and 2017 produced a parameterization of composition dependent collection efficiency very different from the results of previous studies, SMPS data were examined and we determined that there was significant mass that the instrument did not capture because the particles with diameters larger than 465 nm are not counted by this instrument. Data for the ultra-high-sensitivity aerosol spectrometer (UHSAS), SMPS, and ACSM are available for nearly all of 2019. The data from UHSAS and SMPS collected in 2019 were compared. We found that the UHSAS data has particle counts and total volumes significantly less than measured by the condensation particle counter (CPC) and SMPS. The SMPS data were extended by fitting the average volume distribution with a log normal curve and using this relationship to estimate a mass value for the SMPS over extended diameter range. The extended SMPS mass values result in a CDCE parameterization that is in better agreement with the results of other investigators, but is still different from other formulations. The working group recommends that the ACSM data be processed with a collective efficiency (CE)=1, that this be documented clearly in the metadata, and the use of the default CE of 0.5 or a formulation of composition-dependent collection efficiencies (CDCE) chosen by the user should be implemented based on the ammonium nitrate mass fraction. This is clearly necessary for the wintertime SGP ACSM data because of the high nitrate concentrations.
The Deep Operator Network (DeepONet) framework is a different class of neural network architecture that one trains to learn nonlinear operators, i.e., mappings between infinite-dimensional spaces. Traditionally, DeepONets are trained using a centralized strategy that requires transferring the training data to a centralized location. Such a strategy, however, limits our ability to secure data privacy or use high-performance distributed/parallel computing platforms. To alleviate such limitations, in this paper, we study the federated training of DeepONets for the first time. That is, we develop a framework, which we refer to as Fed-DeepONet, that allows multiple clients to train DeepONets collaboratively under the coordination of a centralized server. To achieve Fed-DeepONets, we propose an efficient stochastic gradient-based algorithm that enables the distributed optimization of the DeepONet parameters by averaging first-order estimates of the DeepONet loss gradient. Then, to accelerate the training convergence of Fed-DeepONets, we propose a moment-enhanced (i.e., adaptive) stochastic gradient-based strategy. Finally, we verify the performance of Fed-DeepONet by learning, for different configurations of the number of clients and fractions of available clients, (i) the solution operator of a gravity pendulum and (ii) the dynamic response of a parametric library of pendulums.
In this work, generalized polynomial chaos (gPC) expansion for land surface model parameter estimation is evaluated. We perform inverse modeling and compute the posterior distribution of the critical hydrological parameters that are subject to great uncertainty in the Community Land Model (CLM) for a given value of the output LH. The unknown parameters include those that have been identified as the most influential factors on the simulations of surface and subsurface runoff, latent and sensible heat fluxes, and soil moisture in CLM4.0. We set up the inversion problem in the Bayesian framework in two steps: (i) building a surrogate model expressing the input–output mapping, and (ii) performing inverse modeling and computing the posterior distributions of the input parameters using observation data for a given value of the output LH. The development of the surrogate model is carried out with a Bayesian procedure based on the variable selection methods that use gPC expansions. Our approach accounts for bases selection uncertainty and quantifies the importance of the gPC terms, and, hence, all of the input parameters, via the associated posterior probabilities.
In this paper, we present a flexibility estimation mechanism for buildings’ thermostatically controlled loads (TCLs) to enable the distribution level consumption of the majority of solar photovoltaic (PV) generation by local building TCLs. The local consumption of PV generation provides several advantages to the grid operation as well as the consumers, such as reducing the stress on the distribution network, minimizing voltage fluctuations and two-way power flows in the distribution network, and reducing the required battery storage capacity for PV integration. This would result in increasing the solar PV generation penetration levels. The aims of this study are twofold. First, spectral (frequency) analyses of solar PV power generation together with the power consumption of multiple building TCLs (such as heating, ventilation, and air conditioning (HVAC) systems, water heaters, and refrigerators) are performed. These analyses define the bandwidth over which these TCLs can operate and also describe the PV generation frequency bandwidth. Such spectral analyses, in frequency domain, can help identify the flexible components of PV generation that can be consumed by the various TCLs through optimal building load utilization. Second, a quadratic optimization problem based on model predictive control is formulated to allow consuming most of the low and medium frequency content of the PV power locally by building TCLs, while maintaining occupants’ comfort and TCLs’ physical constraints. The solution to the proposed optimization problem is achieved using optimal control strategies. Numerical results show that most of the low and medium frequency content of the PV generation can be consumed locally by building TCLs. The remaining high-frequency content of the PV generation can then be stored/offset using energy storage systems.
Abstract Observations of core-collapse supernovae (CCSNe) reveal a wealth of information about the dynamics of the supernova ejecta and its composition but very little direct information about the progenitor. Constraining properties of the progenitor and the explosion requires coupling the observations with a theoretical model of the explosion. Here we begin with the CCSN simulations of Couch et al., which use a nonparametric treatment of the neutrino transport while also accounting for turbulence and convection. In this work we use the SuperNova Explosion Code to evolve the CCSN hydrodynamics to later times and compute bolometric light curves. Focusing on Type IIP SNe (SNe IIP), we then (1) directly compare the theoretical STIR explosions to observations and (2) assess how properties of the progenitor’s core can be estimated from optical photometry in the plateau phase alone. First, the distribution of plateau luminosities ( L 50 ) and ejecta velocities achieved by our simulations is similar to the observed distributions. Second, we fit our models to the light curves and velocity evolution of some well-observed SNe. Third, we recover well-known correlations, as well as the difficulty of connecting any one SN property to zero-age main-sequence mass. Finally, we show that there is a usable, linear correlation between iron core mass and L 50 such that optical photometry alone of SNe IIP can give us insights into the cores of massive stars. Illustrating this by application to a few SNe, we find iron core masses of 1.3–1.5 M ⊙ with typical errors of 0.05 M ⊙ . Data are publicly available online on Zenodo: doi: 10.5281/zenodo.6631964 .
Abstract Modern surveys of gravitational microlensing events have progressed to detecting thousands per year, and surveys are capable of probing Galactic structure, stellar evolution, lens populations, black hole physics, and the nature of dark matter. One of the key avenues for doing this is the microlensing Einstein radius crossing time (t E) distribution. However, systematics in individual light curves as well as oversimplistic modeling can lead to biased results. To address this, we developed a model to simultaneously handle the microlensing parallax due to Earth's motion, systematic instrumental effects, and unlensed stellar variability with a Gaussian process model. We used light curves for nearly 10,000 OGLE-III and -IV Milky Way bulge microlensing events and fit each with our model. We also developed a forward model approach to infer the t E distribution by forward modeling from the data rather than using point estimates from individual events. We find that modeling the variability in the baseline removes a source of significant bias in individual events, and the previous analyses overestimated the number of t E > 100 day events due to their oversimplistic model ignoring parallax effects. We use our fits to identify the hundreds filling a regime in the microlensing parameter space that are 50% pure of black holes. Finally, we have released the largest-ever catalog of Markov Chain Monte Carlo parameter estimates for microlensing events.
An improved fine-scale soil moisture (SM) dataset at 1 km grid spacing, covering much of the eastern continental US, was generated by assimilating 9 km Soil Moisture Active Passive (SMAP) SM data into the v4.0.1 Noah-MP land surface model. With 12 ensemble members, the assimilation was carried out using the ensemble Kalman filter algorithm within NASA's Land Information System. The SM analysis for 2016 was fully validated against in situ observations from four different networks and compared with four other existing datasets. Results indicate that this SM analysis surpasses other datasets in top-layer SM distribution, including a machine-learning-based product, despite all SM estimates being less heterogeneous than observed. The analysis of anomalous errors suggests that large similarity in intrinsic errors is likely due to overlapping data sources among the selected SM datasets. More detailed evaluations were performed over two geographic areas. The observations collected by the Atmospheric Radiation Measurement facility in Oklahoma suggest that soil temperature and surface heat fluxes are concurrently simulated with good accuracy. Investigation into the 2016 southeastern US drought response further indicates drier conditions and higher evapotranspiration estimates compared to GLEAMv4.1. Notably, large errors are associated with grids having clay soil textures, underscoring the need for refined model treatments for specific soil types to further improve SM estimates. The dataset is publicly available on Zenodo at https://doi.org/10.5281/zenodo.14370563 (Tai et al., 2024).
Context. While supernova remnants (SNRs) have been considered the most relevant Galactic cosmic ray (CR) accelerators for decades, core-collapse supernovae (CCSNe) could accelerate particles during the earliest stages of their evolution and hence contribute to the CR energy budget in the Galaxy. Some SNRs have indeed been associated with TeV γ -rays, yet proton acceleration efficiency during the early stages of an SN expansion remains mostly unconstrained. Aims. The multi-wavelength observation of SN 2023ixf, a Type II supernova (SN) in the nearby galaxy M 101 (at a distance of 6.85 Mpc), opens the possibility to constrain CR acceleration within a few days after the collapse of the red super-giant stellar progenitor. With this work, we intend to provide a phenomenological, quasi-model-independent constraint on the CR acceleration efficiency during this event at photon energies above 100 MeV. Methods. We performed a maximum-likelihood analysis of γ -ray data from the Fermi Large Area Telescope up to one month after the SN explosion. We searched for high-energy, non-thermal emission from its expanding shock, and estimated the underlying hadronic CR energy reservoir assuming a power-law proton distribution consistent with standard diffusive shock acceleration. Results. We do not find significant γ -ray emission from SN 2023ixf. Nonetheless, our non-detection provides the first limit on the energy transferred to the population of hadronic CRs during the very early expansion of a CCSN. Conclusions. Under reasonable assumptions, our limits would imply a maximum efficiency on the CR acceleration of as low as 1%, which is inconsistent with the common estimate of 10% in generic SNe. However, this result is highly dependent on the assumed geometry of the circumstellar medium, and could be relaxed back to 10% by challenging spherical symmetry. Consequently, a more sophisticated, inhomogeneous characterisation of the shock and the progenitor’s environment is required before establishing whether or not Type II SNe are indeed efficient CR accelerators at early times.
The Sodium Ion Expansion Power Block for Distributed CSP was a three-plus-one-year effort under the Concentrating Solar Power: Advanced Projects Offering Low LCOE Opportunities (CSP: APOLLO) funding program within the U.S. Department of Energy Solar Energy Technologies Office. The primary objective of this project is to develop a dual-stage modular sodium thermal electrochemical converter (Na-TEC) heat engine power block, which can be potentially integrated with either a small-scale dish solar or large-scale heliostats and parabolic trough CSP. Na-TEC is a heat engine that generates electricity through the isothermal expansion of sodium ions. The Na-TEC is a closed system that can theoretically achieve conversion efficiencies above 45% when operating between thermal reservoirs at 1150 K and 550 K. However, thermal designs have confined previous single-stage devices to thermal efficiencies below 20%. To mitigate some of these limitations, we consider dividing the isothermal expansion into two stages; one at the evaporator temperature (1150 K) and another at an intermediate temperature (650 K –1050 K). This dual-stage Na-TEC takes advantage of regeneration and reheating, and could be amenable to better thermal management. In light of this, we first designed and developed a thermo-electrochemical model, and thermodynamically demonstrated how the dual-stage device can improve the efficiency by up to 8% points over the best performing single-stage device. We also established an application regime map for the single- and dual-stage Na-TEC in terms of the power density and the total thermal parasitic loss. Moreover, a thermal design of an axisymmetric dual-stage Na-TEC is developed to guide the scale-up and fabrication of sub-components of prototype module. A reduced-order finite-element model is used in conjunction with a Na-TEC thermodynamic model that was developed to determine the total parasitic heat loss of this dual-stage design. A number of simplifications are applied in the reduced-order model to decrease the computational time while maintaining acceptable accuracy. According to this analysis, a maximum efficiency of 29% and a maximum power output of 125 W can be achieved. Ultimately, we were able to demonstrate thermal efficiency improvements of the Na-TEC heat engine from 19% up to 40.3%, in a dual-stage (non-optimized) prototype module that we designed, fabricated, and tested with high temperature stage at 923 K. Furthermore, a cost-performance analysis for this improved dual-stage design was carried out for distributed-CSP systems. A high-level techno-economic analysis (TEA) explores four scenarios where a Na-TEC is used as the heat engine for a distributed-CSP system. Overnight capital cost and levelized cost of electricity (LCOE) are estimated for a system lifetime of 30 years, revealing that overnight capital costs in a range from $3.57 to $17.71 per We are feasible, which equate to LCOEs from 6.9 to 17.2 cents/kWh e -1 . This analysis makes a significant contribution by concurrently quantifying the efficiency and unit costs for a range of multistage configurations, and demonstrating that a Na-TEC may be a promising alternative to Stirling engines for distributed-CSP systems at residential scale of 1–5 kW e .
High penetrations of distributed photovoltaics (PV) could cause adverse grid impacts, such as voltage violations. The recent development in inverter technologies provides the opportunity to develop control systems to realize effective PV governance and thus to improve dynamic PV hosting capacity for distribution grids. In this paper, a novel distributed energy resource management system (DERMS) solution is proposed by adopting the real-time optimal power flow (OPF) approach for coordinated control of the distributed PV inverters. The proposed approach eliminates the dependence on load knowledge via measurement feedback correction, and it can be implemented in real time. One challenge is that the technique is sensitive to the data availability and integrity of voltage measurements. Therefore, a decentralized DERMS approach is developed by leveraging the concepts of state estimation. The framework and the effectiveness of the solution approach are numerically demonstrated on a real distribution feeder in Southern California.
The protection, control, and monitoring of the power grid is not possible without accurate measurement devices. As the percentage of renewable energy sources penetrating the existing grid infrastructure increases, so do uncertainties surrounding their effects on the everyday operation of the power system. Many of these devices are sources of high-frequency transients. These transients may be useful for identifying certain events or behaviors otherwise not seen in traditional analysis techniques. Therefore, the ability of sensors to accurately capture these phenomena is paramount. In this work, two commercial-grade power system distribution sensors are investigated in terms of their ability to replicate high-frequency phenomena by studying their responses to three events: a current inrush, a microgrid “close-in”, and a fault on the terminals of a wind turbine. Kernel density estimation is used to derive the non-parametric probability density functions of these error distributions and their adequateness is quantified utilizing the commonly used root mean square error (RMSE) metric. It is demonstrated that both sensors exhibit characteristics in the high harmonic range that go against the assumption that measurement error is normally distributed.
We present that a sodium thermal electrochemical converter (Na-TEC) generates electricity directly from heat through isothermal expansion of sodium ions across a beta"-alumina solid-electrolyte. This heat engine has been considered for use with conventional concentrating solar power (CSP) systems before. However, unlike previous single-stage devices, the improved design uses two stages with an interstage reheat, allowing more economical and efficient conversion up to 29% at a hot side temperature of 850 °C. Herein, a cost-performance analysis for this improved design assesses opportunities for distributed-CSP in the context of micro-combined heat and power systems. A high-level techno-economic analysis (TEA) is presented that explores four scenarios where a Na-TEC is used as the heat engine for a distributed-CSP system. Overnight capital cost and levelized cost of electricity (LCOE) are estimated for a system lifetime of 30 years, revealing that overnight capital costs in a range from $3.57 to $17.71 per W e are feasible, which equate to LCOEs from 6.9 to 17.2 cents kWh e -1 . This analysis makes a significant contribution by concurrently quantifying the efficiency and unit costs for a range of multistage configurations, and demonstrating that a Na-TEC may be a promising alternative to Stirling engines for distributed-CSP systems at residential scale of 1–5 kW e .