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.
Under the MuSiKAL project, we developed a framework for a coastal digital twin (DT) platform capable of integrating diverse data resources, configuring multiscale model simulations, performing SciML‐accelerated predictions, with applications primarily driven by storm surge and heavily rainfall events impacting the Gulf Coast of the U.S.
To maximize the performance of high-speed air-breathing engines, such as dual-mode scramjets, the streamwise location of the shock train leading edge (STLE) is ideally placed as far upstream in the isolator as possible while avoiding engine unstart. Thus, it is of interest to quantify and control the STLE location as the vehicle travels along its flight trajectory. The STLE location is typically quantified using wall static pressure measurements but there are often restrictions on the number and placement of transducers, thus reducing the accuracy and overall capability of STLE detection methods. In this work, the Adaptive Pressure Profile (APP) method is introduced to address such sparsity concerns. This method is data driven and does not heavily rely on prior information about the flow regime or engine model. Instead, the \mname method uses real-time pressure measurements from a small number of transducers to adaptively learn the isolator pressure profile. This adaptively-learned profile is fit to the pressure data at each time instance to estimate STLE location. The \mname method produces accurate estimates even when (1) the STLE location is not bounded by two or more transducers or (2) when the STLE location is between two transducers that are situated several duct heights apart. Data from two direct-connect isolator models are used to evaluate the accuracy of the \mname method and demonstrate its robustness for different back-pressure scenarios and transducer configurations.
To maximize the performance of high-speed air-breathing engines, such as dual-mode scramjets, the streamwise location of the shock train leading edge (STLE) is ideally placed as far upstream in the isolator as possible while avoiding engine unstart. Thus, it is of interest to quantify and control the STLE location as the vehicle travels along its flight trajectory. The STLE location is typically quantified using wall static pressure measurements but there are often restrictions on the number and placement of transducers, thus reducing the accuracy and overall capability of STLE detection methods. In this work, the Adaptive Pressure Profile (APP) method is introduced to address such sparsity concerns. This method is data driven and does not heavily rely on prior information about the flow regime or engine model. Instead, the \mname method uses real-time pressure measurements from a small number of transducers to adaptively learn the isolator pressure profile. This adaptively-learned profile is fit to the pressure data at each time instance to estimate STLE location. The \mname method produces accurate estimates even when (1) the STLE location is not bounded by two or more transducers or (2) when the STLE location is between two transducers that are situated several duct heights apart. Data from two direct-connect isolator models are used to evaluate the accuracy of the \mname method and demonstrate its robustness for different back-pressure scenarios and transducer configurations.
Wetlands are responsible for 20%–31% of global methane (CH 4 ) emissions and account for a large source of uncertainty in the global CH 4 budget. Data-driven upscaling of CH 4 fluxes from eddy covariance measurements can provide new and independent bottom-up estimates of wetland CH 4 emissions. Here, we develop a six-predictor random forest upscaling model (UpCH4), trained on 119 site-years of eddy covariance CH 4 flux data from 43 freshwater wetland sites in the FLUXNET-CH4 Community Product. Network patterns in site-level annual means and mean seasonal cycles of CH 4 fluxes were reproduced accurately in tundra, boreal, and temperate regions (Nash-Sutcliffe Efficiency ~0.52–0.63 and 0.53). UpCH4 estimated annual global wetland CH 4 emissions of 146 ± 43 TgCH 4 y –1 for 2001–2018 which agrees closely with current bottom-up land surface models (102–181 TgCH 4 y –1 ) and overlaps with top-down atmospheric inversion models (155–200 TgCH 4 y –1 ). However, UpCH4 diverged from both types of models in the spatial pattern and seasonal dynamics of tropical wetland emissions. We conclude that upscaling of eddy covariance CH 4 fluxes has the potential to produce realistic extra-tropical wetland CH 4 emissions estimates which will improve with more flux data. To reduce uncertainty in upscaled estimates, researchers could prioritize new wetland flux sites along humid-to-arid tropical climate gradients, from major rainforest basins (Congo, Amazon, and SE Asia), into monsoon (Bangladesh and India) and savannah regions (African Sahel) and be paired with improved knowledge of wetland extent seasonal dynamics in these regions.
Abstract Methods in system identification are used to obtain linear time-invariant state-space models that describe how horizontal averages of temperature and humidity of a large cumulus ensemble evolve with time under small forcing. The cumulus ensemble studied here is simulated with cloud-system-resolving models in radiative–convective equilibrium. The identified models extend steady-state linear response functions used in past studies and provide accurate descriptions of the transfer function, the noise model, and the behavior of cumulus convection when coupled with two-dimensional gravity waves. A novel procedure is developed to convert the state-space models into an interpretable form, which is used to elucidate and quantify memory in cumulus convection. The linear problem studied here serves as a useful reference point for more general efforts to obtain data-driven and interpretable parameterizations of cumulus convection.
A key input in modeling population distribution and flows following conflict events is the inclusion of Internally Displaced Person (IDP) flows between administrative units within a country. These flows are critical for capturing population movement and redistribution driven by current events, particularly conflict. In some cases, IDP destination data are available while origin data is incomplete or unavailable. This creates a gap in understanding where displacement is occurring, limiting the ability to model population redistribution accurately. Without origin data, it is not possible to reallocate population flows or accurately represent where displacement is occurring within the country. This report evaluates whether Armed Conflict Location Event Data (ACLED) political violence event data can be used to estimate IDP origin distributions when direct origin data are unavailable. The approach is validated using historical IDP flow data from Lebanon, where both origins and destinations are observed. Results show that ACLED event distributions strongly correspond to observed IDP-origin patterns, particularly when using cumulative 60-day event windows. The method is most reliable for identifying major origin districts and approximating proportional origin shares. However, it is not intended to reconstruct exact individual displacement flows, but rather to provide a probabilistic spatial allocation of displacement origins.
Two approximations to convective transport have been implemented in an offline chemistry transport model (CTM) to explore the impact on calculated atmospheric CO2 distributions. GlobalCO2 in the year 2000 is simulated using theCTM driven by assimilated meteorological fields from the NASA s Goddard Earth Observation System Data Assimilation System, Version 4 (GEOS-4). The model simulates atmospheric CO2 by adopting the same CO2 emission inventory and dynamical modules as described in Kawa et al. (convective transport scheme denoted as Conv1). Conv1 approximates the convective transport by using the bulk convective mass fluxes to redistribute trace gases. The alternate approximation, Conv2, partitions fluxes into updraft and downdraft, as well as into entrainment and detrainment, and has potential to yield a more realistic simulation of vertical redistribution through deep convection. Replacing Conv1 by Conv2 results in an overestimate of CO2 over biospheric sink regions. The largest discrepancies result in a CO2 difference of about 7.8 ppm in the July NH boreal forest, which is about 30% of the CO2 seasonality for that area. These differences are compared to those produced by emission scenario variations constrained by the framework of Intergovernmental Panel on Climate Change (IPCC) to account for possible land use change and residual terrestrial CO2 sink. It is shown that the overestimated CO2 driven by Conv2 can be offset by introducing these supplemental emissions.
The NASA Soil Moisture Active Passive (SMAP) mission provides observations of L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations at a resolution of ~40 km globally every 2-3 days. These observations are routinely assimilated into the NASA Catchment land surface model to generate the Level-4 Soil Moisture (L4_SM) product, which provides global estimates of surface and root-zone soil moisture, soil temperature, and surface fluxes (among others) at 9-km, 3-hourly resolution with ~2.5-day latency. The Catchment land surface model in the L4_SM algorithm is driven with 0.25°, hourly surface meteorological forcing data from the NASA Goddard Earth Observing System (GEOS) "forward-processing" product. Outside of Africa and the high latitudes, the GEOS precipitation forcing is corrected using the Climate Prediction Center Unified (CPCU) gauge-based, 0.5°, daily precipitation product.Soil moisture estimates from the L4_SM product were previously shown to improve over land model-only estimates that do not benefit from the assimilation of Tb observations, thereby demonstrating the value of assimilating SMAP observations for soil moisture estimation. In this presentation, we further isolate the contribution of the gauge-based precipitation corrections to the skill of the L4_SM soil moisture estimates. Specifically, we compare the skill of the L4_SM soil moisture to that of separate model-only and assimilation estimates obtained without the benefit of the gauge-based precipitation corrections.Preliminary results suggest that the soil moisture skill added by the CPCU-based precipitation corrections primarily depends on the quality of the CPCU precipitation product and is greatest in regions where the CPCU gauge network is dense and reliable. Conversely, in regions where the CPCU product is known to be of poor quality, for example in central Australia, the assimilation of SMAP Tb observations provides the most benefit. The presentation will provide an in-depth evaluation of the soil moisture skill of the model-only and assimilation estimates vs. independent in situ and satellite measurements.
The rapid increase in the volume and variety of terrestrial biosphere observations (i.e., remote sensing data and in situ measurements) offers a unique opportunity to derive ecological insights, refine process‐based models, and improve forecasting for decision support. However, despite their potential, ecological observations have primarily been used to benchmark process‐based models, as many past and current models lack the capability to directly integrate observations and their associated uncertainties for parameterization. In contrast, data assimilation frameworks such as the CARbon DAta MOdel fraMework (CARDAMOM) and its suite of process‐based models, known as the Data Assimilation Linked Ecosystem Carbon Model (DALEC), are specifically designed for model‐data fusion. This review, motivated by a recent CARDAMOM community workshop, examines the development and applications of CARDAMOM, with an emphasis on its role in advancing ecosystem process understanding. CARDAMOM employs a Bayesian approach, using a Markov Chain Monte Carlo algorithm to enable data‐driven calibration of DALEC parameters and initial states (i.e., carbon pool sizes) through observation operators. CARDAMOM's unique ability to retrieve localized model process parameters from diverse datasets—ranging from in situ measurements to global satellite observations—makes it a highly flexible tool for analyzing spatially variable ecosystem responses to environmental change. However, assimilating these data also presents challenges, including data quality issues that propagate into model skill, as well as trade‐offs between model complexity, parameter equifinality, and predictive performance. We discuss potential solutions to these challenges, such as reducing parameter equifinality by incorporating new observations. This review also offers community recommendations for incorporating emerging datasets, integrating machine learning techniques, strengthening collaboration with remote sensing, field, and modeling communities, and expanding CARDAMOM's relevance for localized ecosystem monitoring and decision‐making. CARDAMOM enables a deep, mechanistic understanding of terrestrial ecosystem dynamics that cannot be achieved through empirical analyses of observational datasets or weakly constrained models alone.
Identification and zonal classification of the faults is a decisive factor in the relay’s decision to trip or not. Different types of fault like three-phase, line-to-line-to-ground and single-line-to-ground can occur at various locations in the feeder. These faults are seen as the variation in the instantaneous values of three-phase voltages and currents, i.e., waveforms, that are measured at the relay location. The objective of this work is to develop a machine learning model that can identify a fault and classify it to various protection zones based on measured waveforms. In this work, a data-driven relay based on Convolutional Neural Network (CNN) is proposed for fault identification in distribution feeders with high penetration solar PV. The proposed CNN model takes local current and voltage waveforms as input and classify it into fault, no-fault or a capacitor switching. Further, the CNN also attempts to identify fault zones based on the images of waveforms. The overall testing accuracy of the trained model exceeds 95%.
Computational Fluid Dynamics (CFD) analyses of axisymmetric circular-arc boattail nozzles have been completed in support of NASA's Next Generation Launch Technology Program to investigate the effects of high-speed nozzle geometries on the nozzle internal flow and the surrounding boattail regions. These computations span the very difficult transonic flight regime, with shock-induced separations and strong adverse pressure gradients. External afterbody and internal nozzle pressure distributions computed with the Wind code are compared with experimental data. A range of turbulence models were examined in Wind, including an Explicit Algebraic Stress model (EASM). Computations on two nozzle geometries have been completed at freestream Mach numbers ranging from 0.6 to 0.9, driven by nozzle pressure ratios (NPR) ranging from 2.9 to 5. Results obtained on converging-only geometry indicate reasonable agreement to experimental data, with the EASM and Shear Stress Transport (SST) turbulence models providing the best agreement. Calculations completed on a converging-diverging geometry involving large-scale internal flow separation did not converge to a true steady-state solution when run with variable timestepping (steady-state). Calculations obtained using constant timestepping (time-accurate) indicate less variations in flow properties compared with steady-state solutions. This failure to converge to a steady-state solution was found to be the result of difficulties in using variable time-stepping with large-scale separations present in the flow. Nevertheless, time-averaged boattail surface pressure coefficient and internal nozzle pressures show fairly good agreement with experimental data. The SST turbulence model demonstrates the best over-all agreement with experimental data.
The size and complexity of energy system optimization models have increased significantly in recent years, driven by the availability of high-resolution spatial data. We present complexity reduction and solution methods that enable us to efficiently represent high-resolution spatial data in the network design of large-scale energy systems. We aim to reduce the size and enhance the computational efficiency of network design models without sacrificing solution accuracy. Specifically, we first present how to aggregate highly granular data into larger resolutions without averaging out their specific properties through a composite-curve-based approach and then develop a method to linearly represent these curves. Second, we utilize a general clustering method to determine groups of geographically proximate biomass fields and establish a single transportation arc for all of them, reducing the number of transportation-related variables while maintaining an accurate representation of the system. Finally, we introduce a two-step algorithm that decomposes large-scale network design problems into two smaller, more manageable subproblems. We demonstrate the application of our methods using a case study of switchgrass-to-biofuels network design in the eight states of the U.S. Midwest, using realistic and highly explicit spatial data.
The Deep Underground Neutrino Experiment (DUNE) is a next-generation long baseline neutrino oscillation experiment designed to make precision measurements in a 1.2 2.4 MW neutrino beam, which is directed 1285 km from the Fermi National Accelerator Laboratory (Fermilab) to the Sandford Underground Research Facility (SURF) in South Dakota. Neutrinos are measured at two detector facilities: a near detector located at Fermilab close where the beam is produced and a far detector at SURF. The neutrino beam can be configured to be composed primarily of either muon or anti-muon neutrinos. DUNE measures the disappearance of muon and anti-muon neutrinos and appearance electron and anti-electron neutrinos in the neutrino beam. Measuring these neutrino flavour transitions provides DUNE with sensitivity to the neutrino mass ordering, $\delta_{CP}$, $\theta_{13}$, $\theta_{23}$ and the magnitude of $\Delta m^2_{23}$. The DUNE Precision Reaction Independent Spectrum Measurement (DUNE-PRISM) concept presents a novel way to perform a neutrino oscillation analysis, which has the potential to significantly reduce the impact of large systematic uncertainties in the neutrino interaction model. The PRISM method linearly combines measurements of off-axis neutrino interactions at the DUNE near detector to produce data-driven predictions of the oscillated neutrino event rate spectrum at the far detector. By building an oscillated far detector prediction directly from data, any unknown or poorly modelled neutrino interaction effects will be naturally incorporated into the measurement of the parameters of the neutrino oscillation model. This thesis presents the first complete neutrino oscillation analysis for DUNE using the PRISM method. Details of the methodology are fully explained and the prospects for further improvements to the techniques described are highlighted. The expected impact and relative importance of the neutrino flux, cross section and detector systematic uncertainties are described in detail. Finally, this thesis demonstrates that the PRISM method is capable of performing a measurement of the oscillation parameters that is robust against neutrino interaction modelling errors.
For economic reasons, the nuclear industry is witnessing premature closure of nuclear power plants, despite excellent safety records. Operations and Maintenance (O&M) activities are some of the largest costs in operating legacy light-water plants. By reducing O&M costs, nuclear energy can become more economically competitive with other energy sources. This can be achieved by leveraging machine-learning and artificial intelligence technologies to develop data-driven algorithms to better diagnose potential faults within the system. Improved accuracy of the models can lead to a reduction in unnecessary maintenance, thus reducing costs associated with parts, labor, and unnecessary planned, forced, or extended outages. To address these challenges, the goal of this project is to perform research and development in the area of digital monitoring, i.e., the application of advanced sensor technologies (particularly wireless sensor technologies) and data science based analytic capabilities, to advance online monitoring and predictive maintenance in nuclear plants and improve plant performance (efficiency gain and economic competitiveness). This report summarizes the fiscal year 2020 research progress encompassing (1) different wireless vibration sensor and data indicators used to assess the health of a plant asset; (2) development of diagnostic models for fault detection; and (3) development of prognostic models for estimating the health of the system up to 7 days ahead.
Occupancy schedule is one of the key inputs in Building Energy Modeling (BEM) to reflect the interaction between buildings and occupants. Over the past decades, standardized occupancy schedules, developed mainly by engineering rule-of-thumb, have been widely used in BEM due to its simplicity and lack of real measured occupancy data. However, the BEM community has recognized their association with uncertainty and reliability in simulation results from BEM. This study introduces representative occupancy schedules in the U.S. residential buildings, derived from a large smart thermostat dataset and time-series K-means clustering, and an open-source tool to generate a stochastic residential occupancy schedule. Over 90,000 residential occupancy schedules were estimated from the ecobee Donate Your Data dataset. Then, the representative occupancy schedules were identified through clustering. This study further investigated the impacts of three parameters (day, house type, and state) on residential occupancy schedules. Then, a tool, the Residential Occupancy Schedule Simulator (ROSS), is developed using the representative occupancy schedules derived in this study. Details of this tool are presented in this paper. In conclusion, the derived representative occupancy schedules and the ROSS tool can help improve the energy modeling of residential buildings.
The response of liquid xenon to various types of ionizing radiation has been extensively studied theoretically and experimentally. Recent progress in direct detection dark matter experiments highlights the significance of composite events, where multiple particles interact with xenon simultaneously and generate overlapping ionization signatures. In these events, recombination of electrons and ions associated with different primary particles leads to additional suppression of the ionization signal, introducing a new source of uncertainty in dark matter searches and Migdal effect studies. We developed a model to estimate the recombination enhancement for overlapping low-energy particle interactions. This method, which has minimal dependence on xenon microphysics and is primarily driven by existing experimental data, yields predictions that are consistent with available measurements of composite interactions. Furthermore, we demonstrate that the model predictions are robust against xenon microphysics assumptions.
The Sublimator Driven Coldplate (SDC) is a unique piece of thermal control hardware that has several advantages over a traditional thermal control scheme. The principal advantage is the possible elimination of a pumped fluid loop, potentially increasing reliability and reducing complexity while saving both mass and power. Because the SDC requires a consumable feedwater, it can only be used for short mission durations. Additionally, the SDC is ideal for a vehicle with small transport distances and low heat rejection requirements. An SDC Engineering Development Unit was designed and fabricated. Performance tests were performed in a vacuum chamber to quantify and assess the performance of the SDC. The test data was then used to develop correlated thermal math models. Nonetheless, an Integrated Sublimator Driven Coldplate (ISDC) concept is being developed. The ISDC couples a coolant loop with the previously described SDC hardware. This combination allows the SDC to be used as a traditional coldplate during long mission phases and provides for dissimilar system redundancy
Current generation light water reactors (LWRs), sodium cooled fast reactors (SFRs), small modular reactors (SMRs), and next generation nuclear plants (NGNPs) produce harsh environments in and near the reactor core that can severely tax material performance and limit component operational life. To address this issue, several Department of Energy Office of Nuclear Energy (DOE-NE) research programs are evaluating the long duration irradiation performance of fuel and structural materials used in existing and new reactors. In order to maximize the amount of information obtained from Material Testing Reactor (MTR) irradiations, DOE is also funding development of enhanced instrumentation that will be able to obtain in-situ, real-time data on key material characteristics and properties, with unprecedented accuracy and resolution. Such data are required to validate new multi-scale, multi-physics modeling tools under development as part of a science-based, engineering driven approach to reactor development. It is not feasible to obtain high resolution/microscale data with the current state of instrumentation technology. However, ultrasound-based sensors offer the ability to obtain such data if it is demonstrated that these sensors and their associated transducers are resistant to high neutron flux, high gamma radiation, and high temperature. To address this need, the Advanced Test Reactor National Scientific User Facility (ATR-NSUF) is funding an irradiation, led by PSU, at the Massachusetts Institute of Technology Research Reactor to test the survivability of ultrasound transducers. As part of this effort, PSU and collaborators have designed, fabricated, and provided piezoelectric and magnetostrictive transducers that are optimized to perform in harsh, high flux, environments. Four piezoelectric transducers were fabricated with either aluminum nitride, zinc oxide, or bismuth titanate as the active element that were coupled to either Kovar or aluminum waveguides and two magnetostrictive transducers were fabricated with Remendur or Galfenol as the active elements. Pulse-echo ultrasonic measurements of these transducers are made in-situ. This paper will present an overview of the test design including selection criteria for candidate materials and optimization of test assembly parameters, data obtained from both out-of-pile and in-pile testing at elevated temperatures, and an assessment based on initial data of the expected performance of ultrasonic devices in irradiation conditions