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.
Predictive modeling of watershed to river basin scale biogeochemistry through the use of AI techniques and AI-derived model components and the use of AI to design a prediction system comprised of a hierarchy of multi-fidelity models, including AI driven model/component/parameterization selection.
This work describes the initial development of dynamic system models of the cryogenic moderator system (CMS) of the Spallation Neutron Source (SNS) at ORNL as a part of the ORNL LDRD funded project Building TRANSFORM to Accelerate Digital Twin Applications for Nuclear Systems, LOIS 10563. The goal of the work is to start the dynamic system modeling effort with the end goal of using them for real-time applications as digital twins. The CMS is a cryogenic liquid hydrogen flow loop that provides moderation of the neutrons that are generated by the SNS. For optimal neutron production, the CMS needs to maintain a steady and controlled density of cryogenic hydrogen in the moderator section thus requiring precise temperature and pressure control. Due to the varied time scales and system characteristics, control of the system is complex, and diagnostics are also difficult. Difficulty in accessing the flow loop during operations, limited instrumentation and unknown design details of the equipment combine to make the case for having sophisticated digital twin models of the system. Operationally the CMS also provides a strong use case for digital twins due to the constant need of optimization and for troubleshooting/diagnostics. The large amount of data collected which are freely available for using in building the model and verifying and validating the model also makes it a great candidate for a proof-of-concept for digital twins. The project extends ORNL's capacity of development and implementation of the open-source dynamic system modeling tool TRANSFORM for engineering design and digital twin/real-time applications. Specific system configuration data for the CMS have been gathered and an initial dynamic model was created in the TRANSFORM library using Dymola as the solution platform. Models of increasing complexity are created to demonstrate the need for a multi-layered approach in digital twin modeling depending on the scale and phenomena being focused on. The dynamic modeling is shown to bring the dynamic operational aspects to the design process for systems as well as serve as a digital twin to the hardware and allow for models to be tuned and compared against real time operational data. These aims should help to push forward strategic goals of application of digital twins and increase the impact of ORNL systems modeling capabilities with TRANSFORM/Modelica for various advanced energy systems.
The quasi-biennial oscillation (QBO) is the main mode of variability in the tropical stratosphere, influencing the predictability of other regions in the atmosphere through its teleconnections to the stratospheric polar vortices and coupling to surface tropical and extratropical variability. However, climate and forecasting models consistently underestimate QBO amplitudes in the lower stratosphere, likely contributing to their failure to simulate these teleconnections. One underexplored contributor to model biases is missing representation of ozone-radiative feedbacks, which enhance temperature variability in the lower stratosphere, particularly at periods at and greater than the QBO (>28 months). While previous studies suggest that ozone-radiative feedbacks can impact QBO periods, amplitudes and the associated secondary circulation in the lower stratosphere, the reported impacts differ widely among models and are hard to interpret due to differences in methodology. To this end, here we propose a coordinated experimental protocol – held joint between the Atmospheric Processes and their Role in Climate (APARC) Quasi-Biennial Oscillation Initiative (QBOi) and Chemistry Climate Modeling Initiative (CCMI) activities – which is aimed at assessing the coupling between stratospheric ozone, temperature and the circulation. We use the proposed experiments to define the ozone feedback on the QBO in both present-day and idealized (abrupt quadrupling of carbon dioxide) climates. While primary focus is on the QBO, the proposed protocol also enables analysis of other aspects of ozone-radiative-dynamical coupling in the atmosphere, including impacts on the Brewer-Dobson Circulation and tropospheric eddy-driven jet responses to future climate change. Here we document the scientific rationale and design of the QUOCA Phase 1 experiments, summarize the data request, and give a brief overview of participating models. Preliminary results using the NASA Goddard Institute for Space Studies E2-2 climate model are used to illustrate sensitivities to certain methodological choices.
Despite significant advancements in fenestration technology in the last two decades, the thermal transmittance of fenestration products is still significantly higher than that of walls. This corresponds to 60% of the total energy loss of a modern building envelope through the windows. Hence, further development and improvements of fenestration products are necessary. Increasingly stringent codes and standards for fenestration stimulate industry to work on improved solutions. Thus, it is crucial that assessment techniques are able to account for innovations accurately. The thermal effects of non-continuous hardware in window frames are currently ignored by international rating procedures. A preliminary investigation conducted by our team showed significant performance degradation in two of the three out-opening casement profiles caused by the presence of operating hardware. Frames with the structure made of vinyl and fiberglass consist of many air cavities that are penetrated by operating hardware made of highly conductive materials. In these frames, in order to have an accurate assessment, it may be required to employ three-dimensional modeling due to the convective nature of heat transfer within the cavities. However, in this study, we demonstrate that the three-dimensional (3D) effects of non-continuous hardware can be approximated accurately with simpler two-dimensional (2D) simulations. We then develop a simplified model based on weighted average capable of replacing the time- and computation-intensive 3D simulations with 2D simulations and validate it against market available frames and their corresponding hardware. Validation results show that our approximation technique results in discrepancies lower than 0.05 W/(m2K), or 3% of the total thermal transmittance. Thus, we conclude that simplified 2D simulation models may be used for predicting hardware impact in window frames with reasonable accuracy. As windows and glazing structures are becoming ever better thermally insulated, it is becoming even more important to be able to model the impact of the operating hardware on the total thermal performance in order to design the best windows possible and not let the operating hardware ruin an otherwise well-proven design, which is hence addressed in this study.
This paper reports on a measurement of electron-ion recombination in liquid argon in the ICARUS liquid argon time projection chamber (LArTPC). A clear dependence of recombination on the angle of the ionizing particle track relative to the drift electric field is observed. An ellipsoid modified box (EMB) model of recombination describes the data across all measured angles. These measurements are used for the calorimetric energy scale calibration of the ICARUS TPC, which is also presented. The impact of the EMB model is studied on calorimetric particle identification, as well as muon and proton energy measurements. Accounting for the angular dependence in EMB recombination improves the accuracy and precision of these measurements.
Inhaled radioactive materials can pose a long-term health concern, as the material can be incorporated into the body’s metabolic pathways and remain in organs and tissues for extended durations. During the retention period, the radioactive material may localize in a source organ and irradiate adjacent target organs and tissues. Distribution of these materials changes over time, requiring biokinetic modeling to evaluate their movement through various tissues and organs. The evolving distribution depends on multiple inputs characterizing the inhaled material, such as particle size and size distribution, particle density, aspect ratio, specific radionuclide, the chemical form, and solubility. In addition, biological parameters such as breathing rate, breathing type (nasal or nasal/oral), respiratory system morphometry, tidal volume, functional residual capacity, and anatomical dead space all influence material transport. These aerosol properties and physiological characteristics of the respiratory tract jointly define a range of initial conditions that influence the time-dependent distribution of radioactive material. To evaluate both uncertainty in the initial conditions of inhalation exposure and the final output (committed effective dose) from biokinetic models, a Python-based software tool, Radiological Exposure Dose Calculator (REDCAL), was developed to propagate uncertainty within the human respiratory tract model. Focusing on deposition fraction uncertainty, the primary objective was to characterize the initial activity distribution across respiratory regions as a function of anticipated particle sizes and distributions. The impact of the deposition fraction uncertainty was propagated to committed effective dose coefficients for selected radionuclides in a companion publication. For each particle size, a lognormal distribution, characterized by its geometric mean as defined within ICRP Publication 66, serves as the basis for introducing uncertainty into the physical processes governing deposition in various lung regions. Finally, this study addresses the deposition process and examines how uncertainty in deposition mechanisms affects activity distribution in the airways, ultimately presenting the expected range and standard deviation of deposited activity as a function of particle size.
An approach to increase the value of carbon fiber for wind turbines blades, and other compressive strength driven designs, is to identify pathways to increase its cost-specific compressive strength. A finite element model has been developed to evaluate the predictiveness of current finite element methods and to lay groundwork for future studies that focus on improving the cost-specific compressive strength. Parametric studies are conducted to understand which uncertainties in the model inputs have the greatest impact on compressive strength predictions. Furthermore, a statistical approach is also presented that enables the micromechanical model, which is deterministic, to efficiently account for statistical variability in the fiber misalignment present in composite materials; especially if the results from the hexagonal and square pack models are averaged. The model was found to agree well with experimental results for a Zoltek PX-35 pultrusion. The sensitivity studies suggest that the fiber packing and the interface shear strength have the greatest impact on compressive strength prediction for the fiber reinforced polymer studied here. Based on the performance of the modeling approach presented in this work, it is deemed sufficient for future work which will seek to identify carbon fiber composites with improved cost-specific compressive strength.
The dynamic analysis of fusion power plant (FPP) fuel cycles highlights the challenge of achieving tritium self-sufficiency in future FPPs. While state-of-the-art fuel cycle models offer valuable insights into the necessary design parameters for attaining tritium self-sufficiency, none of these models currently consider the impact of tritium trapping within fuel cycle components. However, detailed analysis of individual components reveals that substantial amounts of tritium can be trapped within the first wall, divertors, and breeding blanket systems, suggesting that tritium trapping may significantly influence the FPP ability to achieve self-sufficiency. The compounded effects of additional tritium traps generated by irradiation effects and component replacements further exacerbate this challenge. The novelty of this work is the integration of an explicit, physics-based model for tritium trapping, evolution of damage-induced traps, and component replacements into a dynamic, system-level model of a fuel cycle. The results show an increase of a factor 10 3 – 10 4 of tritium inventory in the first wall and vacuum vessel of an ARC-class FPP when accounting for the aforementioned phenomena. This, coupled with the replacement of components subject to significant tritium trapping, slows down fuel cycle dynamics, resulting in an extended tritium doubling time (50% increase), higher start-up inventory (30% increase), and higher required tritium breeding ratio (2%–5%) compared to a scenario without tritium trapping.
Fluorinated graphite materials are of interest for an assortment of applications and can be synthesized under a variety of synthetic conditions from many different types of carbon. Due to such variations, structural disorders in the form of defects and polymorphism are often present. Here, we investigate the impact of local structural variations on the C–F bond dissociation energies (BDEs) in carbon-based fluoride materials using density functional theory (DFT) computational methods. Employing fluorographene (FG) cluster models, we determine the impact of different C–F bonding configurations in the core of each platelet on the equilibrium BDEs for each C–F bond. The introduction of structural disorder decreases the first C–F BDE by approximately 1 eV compared to the canonical arrangement of axial C–F bonds ordered as in a network of cyclohexane “chairs”. Variability of calculated BDEs among the different polymorphs decreases upon subsequent F removal. Common structural tendencies of the adiabatic defluorination pathways for each polymorph are identified. In conclusion, our analysis suggests that at F/C ratios near 1.0, disorder in the local structure can play a significant role in the energetics of the initial carbon fluoride defluorination and that the influence of this configurational disorder diminishes with decreasing F/C ratios.
Measurements of carbonyl sulfide (OCS) enable independent estimates of regional stomatal conductance provided that non-stomatal OCS fluxes are well constrained. OCS is taken up through plant leaves, following the same pathway as CO 2 ; in contrast to CO 2 , OCS is irreversibly destroyed in plant leaves and plants do not typically exhibit OCS emissions. Ecosystem uptake of OCS can indicate changes in stomatal opening. Here we present an empirical model to assess the potential impact of soil OCS exchange, the non-Stomatal OCS exchange Empirical Model (SOCSEM, version 0). We created biome-specific response curves characterizing soil OCS exchange and restricted the model design to require only knowledge of soil moisture and surface temperature because remote sensing observations are available for these two features. Comparing the model to field-based chamber observations reveal deviations that can be attributed to missing complexity of the ground surface (having excluded litter and plants without regulated stomata), shortwave radiation, or the soil environment. For agricultural regions with known net emissions, we use remotely-sensed surface temperature data and demonstrate that data resolution can affect anticipated fluxes. We further investigate the influence of regions with unknown soil OCS responses, for example, Arctic tundra. We compare our model to a process-based and respiration-based soil OCS exchange model that has been implemented in a land surface model. Further field study of tropical and arctic ecosystems in conjunction with studies of non-stomatal surfaces in addition to soil (e.g., bryophytes) will increase confidence in applying OCS as a regional tracer for stomatal conductance.
Novel high performance computing systems that feature hybrid architectures require large scale code refactoring to unravel underlying exploitable parallelism. Such redesign can often be accompanied with machine-precision changes as the order of computation cannot always be maintained. For chaotic systems like climate models, these round-off level differences can grow rapidly. Systematic errors may also manifest initially as machine-precision differences. Isolating genuine round off level differences from such errors remains a challenge. Here, we apply two-sample equality of distribution tests to evaluate statistical reproducibility of the ocean model component of US Department of Energy's Energy Exascale Earth System Model (E3SM). A 2-year control simulation ensemble is compared to a modified ensemble as a test case - after a known non-bit-for-bit change in a model component is introduced - to evaluate the null hypothesis that the two ensembles are statistically indistinguishable. To quantify the false negative rates of these tests, we conduct a formal power analysis using a targeted suite of short simulation ensembles. The ensemble suite contains several perturbed ensembles, each with a progressively different climate than the baseline ensemble - obtained by perturbing the magnitude of a single model tuning parameter, the Gent and McWilliams κ, in a controlled manner. The null hypothesis is evaluated for each of perturbed ensembles using these tests. The power analysis informs on the detection limits of the tests for given ensemble size allowing model developers to evaluate the impact of an introduced non-bit-for-bit change to the model.
Both field line curvature (FLC) and electromagnetic ion cyclotron (EMIC) wave scattering are believed to be associated with energetic ion precipitation, but their relative contributions to ionospheric ion precipitation and subsequent effects are still unclear. Here, in this study, by using a kinetic ring current model, we investigate their impact on the ionosphere from two aspects: the global distribution of ion precipitation and resulting ionospheric conductance. Our results show that the intensity and coverage of ion precipitation due to EMIC waves are larger than that due to FLC scattering, while the latter mostly contributes to ion precipitation in outer regions (L > 4–5). We then estimate the conductance with empirical models using the simulated ion precipitation energy flux. When the EMIC wave associated proton precipitation is included, the conductance is significantly enhanced in the dusk-to-midnight sector and has a wide magnetic latitude (MLAT) range from around 52° to 62°, considerably altering the electric potential in the dusk sector, which further influences particle dynamics in the magnetosphere. On the contrary, the proton precipitation caused by FLC scattering only occurs at higher latitudes above MLAT = 60° and the corresponding conductance is slightly enhanced at midnight, with negligible influence on the convective electric potential. Although electron precipitation-associated conductance is predominant globally, the results show that proton precipitation can also play an important role in ionospheric electrodynamics, especially when EMIC wave-associated precipitating proton energy flux exceeds that of electrons in the dusk sector, a region where many subauroral coupling processes occur.
Feedstock attributes of lignocellulosic biomass, such as particle size, compositional makeup, and moisture content, can vary substantially even within pre-processed materials and have a significant effect on conversion in fast pyrolysis-based processes. However, the economic impacts of these attributes are not well understood. To address this, biomass deconstruction phenomena captured with a versatile particle-scale simulation were linked to techno-economic impacts via reduced-order models. Parametric analysis of the particle-scale model, which was validated using literature data, was used in combination with multiple linear regression models to develop correlations between feedstock attributes and yields of pyrolysis oil, gas, and char. Yields were then correlated with the minimum fuel selling price (MFSP) using a techno-economic model, bridging the gap between physics-based biomass conversion simulations and predictions of MFSP for a catalytic fast-pyrolysis process. Empirical correlations derived from the literature regarding the impact of mineral matter (ash) on oil yield were also considered. The model correlations deployed in the integrated framework capture the impacts of variation in feedstock attributes on the MFSP. Variations in ash were shown to have the biggest impact, varying MFSP by -13%/+22% due to catalytic effects and lower relative amounts of convertible lignocellulosic material. It was also found that, if ash can be controlled to low levels, the increased extractives in forest residues can help compensate for some yield losses associated with increased ash. As a result, other inputs considered (particle size, moisture content, and reactor temperature) had relatively negligible effects on process economics within the ranges analyzed considering particle-scale effects alone.
Dust emissions related to anthropogenic activities (i.e., anthropogenic dust) is not represented in most global climate models and its radiative impact remains unassessed. In this study, we develop a new and physically based method to parameterize anthropogenic dust emission from agricultural sources (i.e., agricultural dust, or AD) based on the DOE's Energy Exascale Earth System Model version 1 (E3SMv1). This method relates AD emission to the crop land use fraction in the E3SMv1 land component. Major AD sources simulated by our parameterization include those over Central America, the Sahel, North India, and North China. The annual averaged AD emission with diameter less than 10 μm is 567 Tg yr −1 in present-day (year 2000), which contributes to 13.3% of total dust emission. Model evaluation against satellite and ground-based observations shows that the new parameterization can represent AD emissions and global dust cycle reasonably well. We find that the total dust emission increases by 13.1% (495 Tg yr −1 ) from 1850 to 2000 due to the cropland land use fraction changes without considering the impact from climate changes. This induces a net dust direct effective radiative forcing of −0.041 W m −2 at top of the atmosphere. This dust-induced cooling exceeds 10% of the total anthropogenic aerosol direct effective radiative forcing from 1750 to 2014 estimated by the Intergovernmental Panel on Climate Change Sixth Assessment Report. Our findings indicate an important role of anthropogenic dust in global climate change, which should be included in future climate change assessments.
Abstract This work presents a pseudo-two-dimensional proton exchange membrane fuel cell (PEMFC) model incorporating Nafion ionomer structure–property relationships in the cathode catalyst layer (CL) to capture and explain losses at low Pt loading. Structural data from neutron reflectometry and thin film Nafion conductivity measurements predict variations in the oxygen diffusion coefficient and ionic conductivity with changing CL ionomer thickness and Pt loading. By including these structure–property relationships, predicted polarization curves agree closely with previously published experimental data from cells with Pt loadings between 0.025 and 0.2 mg/cm2. Results demonstrate that structure–property relationships based on physically measurable ionomer and CL properties provide a feasible interpretation of PEMFC CL phenomena for a range of Pt loadings and help explain previously unaccounted-for losses at low Pt. Results also show that simulations must account for surface species coverage variations in order to properly capture the kinetic losses. Finally, results suggest that an increase in ionomer thickness surrounding the C/Pt surfaces may lead to improved cell performance due to improved ionic conductivity.
Land-atmosphere interactions are central to the evolution of the atmospheric boundary layer and the subsequent formation of clouds and precipitation. Existing global climate models represent these connections with bulk approximations on coarse spatial scales, but observations suggest that small-scale variations in surface characteristics and co-located turbulent and momentum fluxes can significantly impact the atmosphere. Recent model development efforts have attempted to capture this phenomenon by coupling existing representations of subgrid-scale (SGS) heterogeneity between land and atmosphere models. Such approaches are in their infancy and it is not yet clear if they can produce a realistic atmospheric response to surface heterogeneity. Here, we implement a parameterization to capture the effects of SGS heterogeneity in the Community Earth System Model (CESM2), and compare single-column simulations against high-resolution Weather Research and Forecasting (WRF) large-eddy simulations (LESs), which we use as a proxy for observations. The CESM2 experiments increase the temperature and humidity variances in the lowest atmospheric levels, but the response is weaker than in WRF-LES. In part, this is attributed to an underestimate of surface heterogeneity in the land model due to a lack of SGS meteorology, a separation between deep and shallow convection schemes in the atmosphere, and a lack of explicitly represented mesoscale secondary circulations. These results highlight the complex processes involved in capturing the effects of SGS heterogeneity and suggest the need for parameterizations that communicate their influence not only at the surface but also vertically.
U.S. critical infrastructure assets are often designed to operate for decades, and yet long-term planning practices have historically ignored climate change. With the current pace of changing operational conditions and severe weather hazards, research is needed to improve our ability to translate complex, uncertain risk assessment data into actionable inputs to improve decision-making for infrastructure planning. Decisions made today need to explicitly account for climate change – the chronic stressors, the evolution of severe weather events, and the wide-ranging uncertainties. If done well, decision making with climate in mind will result in increased resilience and decreased impacts to our lives, economies, and national security. We present a three-tier approach to create the research products needed in this space: bringing together climate projection data, severe weather event modeling, asset-level impacts, and contextspecific decision constraints and requirements. At each step, it is crucial to capture uncertainties and to communicate those uncertainties to decision-makers. While many components of the necessary research are mature (i.e., climate projection data), there has been little effort to develop proven tools for long-term planning in this space. The combination of chronic and acute stressors, spatial and temporal uncertainties, and interdependencies among infrastructure sectors coalesce into a complex decision space. By applying known methods from decision science and data analysis, we can work to demonstrate the value of an interdisciplinary approach to climate-hazard decision making for longterm infrastructure planning.
The body of knowledge on gray box building energy modeling (GBBEM) has been developed over the past few decades and has undergone some important changes recently. Starting with simple methods and simple buildings, the science of GBBEM has grown to encompass complex and more computationally intensive techniques and complex commercial buildings. Numerous works including a recent review have considered model structure and inputs in a nearly systematic way, but no extant work systematically reviews the approaches for GBBEM parameter search initialization and final identification, despite this being arguably the most difficult and impactful part of the modeling process. To this end, we critically review 55 extant works describing advantages, limitations, and domain of applicability of several classes of parameter initialization and optimization techniques specifically for GBBEM. We categorize the classes of methods and analyze the evidence of their applicability for different applications within the field of GBBEM. We find an emerging consensus that initialization of parameter searches for anything other than the simplest building elements is often challenging and sometimes requires a stochastic approach to begin the parameter identification process. After this initial process, faster methods have been used in some cases but often the nature of the problem requires stochastic methods for this portion of the process as well. For less complex systems, more deterministic and efficient methods have been shown to be effective. Finally, we draw conclusions as to the domain of applicability of different classes of initialization and optimization techniques for GBBEEM and offer suggestions for research directions that are likely to prove fruitful.