Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Coupling”

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 973 records · Page 54

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗

Initial Investigation of Coupled Electron-Photon Transport with SCEPTRE and RAPTURE (Rev.1.1)

The SCEPTRE1 code models coupled photon and electron transport in material media using discrete ordinates methods. The code can use different solvers in 1D and 3D in parallel, promising users the ability to run problems that are significant both in physics and computational scale. One capability that is perhaps unique at this time in the NNSA complex is SCEPTRE’s coupled photon– electron transport physics in a discrete ordinates approach. For example, incident x–rays can generate a significant electron population, whose transport can be simulated within the same code. In this context, fully–coupled means that photons can produce electrons and electrons can produce photons; uncoupled or no–coupling means photons do not produce electrons. Partial–coupling, an intermediate option in which photons can produce electrons but electrons do not produce photons, is not considered in the present work. The MCNP code can treat this physics with a Monte Carlo approach. MCNP should be viewed as complementary to SCEPTRE in both computational method (Monte Carlo versus discrete ordinates) as well as data libraries. SCEPTRE and MCNP should provide important checks on each other, especially from complementary modeling frameworks of a S n (discrete ordinates) mesh based method versus Monte Carlo particle based method. In order to learn more about SCEPTRE, and how it might inform engineering work at the Laboratory, we used it to carry out a set of initial studies. This was a series of 1D simulations of photon beams hitting Aluminum slabs. Our studies mainly used SCEPTRE directly, and made little use of the closely associated RAPTURE code (more on this below). The goals of these studies were first, to familiarize ourselves with operating the code, and second, to begin to assess the importance of coupled electron–photon transport. To support these studies, we began to develop a set of scripts to initialize inputs, run tools, process outputs, and plot results. An initial version of these scripts has been committed to the Lightning workflow repository. This memo relates the results of those studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Impact of a Regional Drought on Terrestrial Carbon Fluxes and Atmospheric Carbon: Results from a Coupled Carbon Cycle Model

Understanding the underlying processes that control the carbon cycle is key to predicting future global change. Much of the uncertainty in the magnitude and variability of the atmospheric carbon dioxide (CO2) stems from uncertainty in terrestrial carbon fluxes, and the relative impacts of temperature and moisture variations on regional and global scales are poorly understood. Here we investigate the impact of a regional drought on terrestrial carbon fluxes and CO2 mixing ratios over North America using the NASA Goddard Earth Observing System (GEOS) Model. Results show a sequence of changes in carbon fluxes and atmospheric CO2, induced by the drought. The relative contributions of meteorological changes to the neighboring carbon dynamics are also presented. The coupled modeling approach allows a direct quantification of the impact of the regional drought on local and proximate carbon exchange at the land surface via the carbon-water feedback processes.

CO↗

Cost-optimal evaluation of centralized and distributed microgrid topologies considering voltage constraints

Optimal design of hybrid renewable mini-grids requires both economic and power quality considerations. Existing modeling approaches address these considerations via separate or loosely coupled models. Here, we extend REopt—a techno-economic optimization model developed at the National Renewable Energy Laboratory—to consider both within a single model. REopt formulates the design problem as a mixed-integer linear program that solves for a site's optimal technology mix, sizing, and operation to minimize life cycle cost. REopt has traditionally assumed a single node system. In the work presented here, we expand the REopt platform to consider multiple connected nodes with associated voltage constraints. In order to do this, we model power flow using a fixed-point linear approximation method. Additionally, we then use the model to explore design considerations of mini-grids in Sub-Saharan Africa. Specifically, we evaluate under what combinations of transmission line distance and capacity it is technically viable and economically preferable to build multiple isolated mini-grids versus an interconnected, centralized system.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Parametric Studies of Heat Pump Water Heater Using Low GWP Refrigerants

A storage heat pump water heater (HPWH) comprises a vapor compression system and a hot water tank. It can lead to more than 70% energy saving, compared with electric resistance heating. HPWH usually uses a wrapped-tank coil as the water heater. The wrapped coil has a strong interaction with the water tank and impose a new challenge for the engineering design and simulation study. Based on the US DOE/ORNL Heat Pump Design Model (HPDM), i.e., a public-domain, component-based thermal system modelling platform, a segment-to-segment wrapped-tank coil model coupled with a one-dimensional stratified water tank model were developed, which facilitates modeling and design HPWHs having various configurations and components. Additionally, the hardware-based design model was used to evaluate refrigerants having lower global warming potentials to replace R-134a, including R-1234yf, R-1234ze, R-290, R-513A and R-450A, when effects of compressor size, refrigerant charge on 24-hour uniform energy factor (UEF) and first-hour rating (FHR) were investigated.In the drop-in replacements, R-1234ze achieved the best UEF, and R-290 resulted in the lowest UEF. When downsizing the compressor to maximize the UEFs, the UEFs of R-134a, R-290, R1234yf, R-450A and R-513A ascended to the range of 4.4 to 4.53, with R-513A having the highest UEF. The R-1234ze ended up having the lowest UEF of 4.27. Refrigerants having larger volumetric vaporization heat led to higher FHRs if using the fixed-size compressor, that R-290 resulted in the highest FHR.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Roles of irrigation and reservoir operations in modulating terrestrial water and energy budgets in the Indian subcontinental river basins

Water management activities affect the terrestrial water cycle by increasing water supply and hence changing the partitioning of surface water and energy budgets over irrigated agricultural land, as well as modifying streamflow seasonality through reservoir operations (for irrigation, hydropower generation, and other purposes). Indian sub-continental basins are strongly influenced by human activities. However, the dynamics of how anthropogenic activities interact with the terrestrial water cycle have not been sufficiently studied over Indian sub-continental Basins. In this study, we applied a well-evaluated water management model, coupled with a large-scale land surface hydrological model with irrigation and reservoir operation schemes to examine the impacts of irrigation and reservoir operations on terrestrial water budgets over the Indian sub-continental basins. We conducted two simulations: one NATURAL (i.e., no irrigation and no reservoirs) and water management scenario (i.e., with irrigation and with reservoir regulations, referred to as MANAGED) over 18 basins from 1951 to 2012. The water management impacts were evaluated by comparing a number of key water and energy budget terms such as evapotranspiration (ET), total runoff (TR), land surface temperature (LST), reservoir storage, and discharge. Our results show that water management activities have significant (p-value < 0.05) impacts on ET and LST in most of the river basins during the pre-monsoon season (February-May). Comparing to the NATURAL scenario, ET under the MANAGED scenario is significantly increased in most basins during pre- and post-monsoon seasons especially in semi-arid regions, and LST is reduced predominately in the pre-monsoon season when irrigation water demand is high.

Shah, Harsh L.↗

Global-Scale Convergence Obscures Inconsistencies in Soil Carbon Change Predicted by Earth System Models

Soil carbon (C) responses to environmental change represent a major source of uncertainty in the global C cycle. Feedbacks between soil C stocks and climate drivers could impact atmospheric CO 2 levels, further altering the climate. Here, we assessed the reliability of Earth system model (ESM) predictions of soil C change using the Coupled Model Intercomparison Project phases 5 and 6 (CMIP5 and CMIP6). ESMs predicted global soil C gains under the high emission scenario, with soils taking up 43.9 Pg (95% CI: 9.2–78.5 Pg) C on average during the 21st century. The variation in global soil C change declined significantly from CMIP5 (with average of 48.4 Pg [95% CI: 2.0–94.9 Pg] C) to CMIP6 models (with average of 39.3 Pg [95% CI: 23.9–54.7 Pg] C). For some models, a small C increase in all biomes contributed to this convergence. For other models, offsetting responses between cold and warm biomes contributed to convergence. Although soil C predictions appeared to converge in CMIP6, the dominant processes driving soil C change at global or biome scales differed among models and in many cases between earlier and later versions of the same model. Random Forest models, for soil carbon dynamics, accounted for more than 63% variation of the global soil C change predicted by CMIP5 ESMs, but only 36% for CMIP6 models. Although most CMIP6 models apparently agree on increased soil C storage during the 21st century, this consensus obscures substantial model disagreement on the mechanisms underlying soil C response, calling into question the reliability of model predictions.

54 ENVIRONMENTAL SCIENCES↗

A variational approach for dissipative quantum transport in a wide parameter space

Recent development of theoretical method for dissipative quantum transport has achieved notable progresses in the weak or strong electron-phonon coupling regime. However, a generalized theory for dissipative quantum transport in a wide parameter space had not been established. In this work, a variational polaron theory for dissipative quantum transport in a wide range of electron-phonon coupling is developed. The optimal polaron transformation is determined by the optimization of the Feynman-Bogoliubov upper bound of free energy. The free energy minimization ends up with an optimal mean-field Hamiltonian and a minimal interaction Hamiltonian. Hence, second-order perturbation can be applied to the transformed system, resulting in an accurate and efficient method for the treatment of dissipative quantum transport with different electron-phonon coupling strength. Numerical benchmark calculation on a single site model coupled to one phonon mode is presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The origin of unequal bond lengths in the C 1 B 2 state of SO 2 : Signatures of high-lying potential energy surface crossings in the low-lying vibrational structure

Here the C 1 B 2 state of SO 2 has a double-minimum potential in the antisymmetric stretch coordinate, such that the minimum energy geometry has nonequivalent SO bond lengths. The asymmetry in the potential energy surface is expressed as a staggering in the energy levels of the v' 3 progression. We have recently made the first observation of low-lying levels with odd quanta of v' 3 , which allows us--in the current work--to characterize the origins of the level staggering. Our work demonstrates the usefulness of low-lying vibrational level structure, where the character of the wavefunctions can be relatively easily understood, to extract information about dynamically important potential energy surface crossings that occur at much higher energy. The measured staggering pattern is consistent with a vibronic coupling model for the double-minimum, which involves direct coupling to the bound 2 1 A 1 state and indirect coupling with the repulsive 3 1 A 1 state. The degree of staggering in the v' 3 levels increases with quanta of bending excitation, which is consistent with the approach along the C state potential energy surface to a conical intersection with the 2 1 A 1 surface at a bond angle of ~145°.

74 ATOMIC AND MOLECULAR PHYSICS↗

Numerical Study of Liquid Piston Compression Using Large-Eddy Simulation and Volume-of-Fluid Approach

Efforts to increase the efficiency of residential and commercial air conditioners and heat pumps have demonstrated that the compressor accounts for most of the system’s electrical energy usage. Therefore, the efficiency of this component should be improved to reduce its energy usage. The US Department of Energy’s Oak Ridge National Laboratory developed a near-isothermal liquid piston compressor (LPC) that uses propylene glycol (PG) to compress CO2. This report presents numerical studies of the LPC in which the compression chamber fills with injected PG from the bottom inlet. Numerical simulations were performed using the large-eddy simulation (LES) with the wall-adapting local eddy-viscosity (WALE) subgrid-scale model coupled with the multiphase volume-of-fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfer within the compression chamber. In this effort, effects of boundary conditions applied to the LES-VOF calculations (e.g., no wall, an adiabatic wall, and a wall with a heat flux subscribed) to the overall pressure and temperature of CO2 gas as well as the transient evolution of flow and heat transfer evolution within the compression chamber are investigated and discussed. It was found that the LES calculation with no wall have shown no dynamical flow patterns and the volume-averaged temperature of CO2 increased from 305 to 392.7 K, while LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO2 temperatures. Results of LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO2 gas layers and circulation flow patterns within the chamber that contribute to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor.

Nguyen, Thien D.↗

SU-F-R-10: Selecting the Optimal Solution for Multi-Objective Radiomics Model

Purpose: To develop an evidential reasoning approach for selecting the optimal solution from a Pareto solution set obtained by a multi-objective radiomics model for predicting distant failure in lung SBRT. Methods: In the multi-objective radiomics model, both sensitivity and specificity are considered as the objective functions simultaneously. A Pareto solution set with many feasible solutions will be resulted from the multi-objective optimization. In this work, an optimal solution Selection methodology for Multi-Objective radiomics Learning model using the Evidential Reasoning approach (SMOLER) was proposed to select the optimal solution from the Pareto solution set. The proposed SMOLER method used the evidential reasoning approach to calculate the utility of each solution based on pre-set optimal solution selection rules. The solution with the highest utility was chosen as the optimal solution. In SMOLER, an optimal learning model coupled with clonal selection algorithm was used to optimize model parameters. In this study, PET, CT image features and clinical parameters were utilized for predicting distant failure in lung SBRT. Results: Total 126 solution sets were generated by adjusting predictive model parameters. Each Pareto set contains 100 feasible solutions. The solution selected by SMOLER within each Pareto set was compared to the manually selected optimal solution. Five-cross-validation was used to evaluate the optimal solution selection accuracy of SMOLER. The selection accuracies for five folds were 80.00%, 69.23%, 84.00%, 84.00%, 80.00%, respectively. Conclusion: An optimal solution selection methodology for multi-objective radiomics learning model using the evidential reasoning approach (SMOLER) was proposed. Experimental results show that the optimal solution can be found in approximately 80% cases.

60 APPLIED LIFE SCIENCES↗

Causes and Implications of Extreme Atmospheric Moisture Demand during the Record-Breaking 2011 Wildfire Season in the Southwestern United States

In 2011, exceptionally low atmospheric moisture content combined with moderately high temperatures to produce a record-high vapor pressure deficit (VPD) in the southwestern United States (SW). These conditions combined with record-low cold-season precipitation to cause widespread drought and extreme wildfires. Although interannual VPD variability is generally dominated by temperature, high VPD in 2011 was also driven by a lack of atmospheric moisture. The May–July 2011 dewpoint in the SW was 4.5 standard deviations below the long-term mean. Lack of atmospheric moisture was promoted by already very dry soils and amplified by a strong ocean-to-continent sea level pressure gradient and upper-level convergence that drove dry northerly winds and subsidence upwind of and over the SW. Subsidence drove divergence of rapid and dry surface winds over the SW, suppressing southerly moisture imports and removing moisture from already dry soils. Model projections developed for the fifth phase of the Coupled Model Intercomparison Project (CMIP5) suggest that by the 2050s warming trends will cause mean warm-season VPD to be comparable to the record-high VPD observed in 2011. CMIP5 projections also suggest increased interannual variability of VPD, independent of trends in background mean levels, as a result of increased variability of dewpoint, temperature, vapor pressure, and saturation vapor pressure. Increased variability in VPD translates to increased probability of 2011-type VPD anomalies, which would be superimposed on ever-greater background VPD levels. Although temperature will continue to be the primary driver of interannual VPD variability, 2011 served as an important reminder that atmospheric moisture content can also drive impactful VPD anomalies.

54 ENVIRONMENTAL SCIENCES↗

Analysis of Low-Frequency Precipitation Variability in CMIP5 Historical Simulations for Southwestern North America

Drier future conditions are projected for the arid southwest of North America, increasing the chances of the region experiencing severe and prolonged drought. To examine the mechanisms of decadal variability, 47 global climate model historical simulations performed for phase 5 of the Coupled Model Intercomparison Project (CMIP5) were assessed. On average, the CMIP5 models have higher climatological precipitation over the past century in southwestern North America than current instrumental or reanalysis products. The timing of the winter peak in climatological precipitation over California and Nevada is accurately represented. Models with resolutions coarser than 2° show a larger spread in the location and strength of the North American monsoon ridge and subsequent summer precipitation, in comparison with the higher-resolution models. Less than 20% of decadal variability in wintertime precipitation over California is associated with North Pacific sea surface temperature anomalies, a larger proportion than is associated with the tropical forcing but not sufficient for making decadal drought predictions. North American monsoon precipitation is strongly associated with local land temperatures on interannual-to-decadal time scales attributable to evaporative cooling and radiation changes driven by varying cloud cover. Soil moisture in Texas and Oklahoma in April is shown to be positively correlated with monsoon precipitation for the following summer, indicating a potential source of nonoceanic interseasonal persistence in southwestern North American hydroclimate. To make meaningful decadal predictions in the future, it is likely that forecasting will move away from sea surface temperature–driven anomaly patterns, and focus on land surface processes, which can allow persistence of precipitation anomalies via feedbacks.

54 ENVIRONMENTAL SCIENCES↗

Representation of Modes of Variability in Six U.S. Climate Models

We compare the performance of several modes of variability across six U.S. climate modeling groups, with a focus on identifying robust improvements in recent models [including those participating in phase 6 of the Coupled Model Intercomparison Project (CMIP)] compared to previous versions. In particular, we examine the representation of the Madden–Julian oscillation (MJO), El Niño–Southern Oscillation (ENSO), the Pacific decadal oscillation (PDO), the quasi-biennial oscillation (QBO) in the tropical stratosphere, and the dominant modes of extratropical variability, including the southern annular mode (SAM), the northern annular mode (NAM) [and the closely related North Atlantic Oscillation (NAO)], and the Pacific–North American pattern (PNA). Where feasible, we explore the processes driving these improvements through the use of “intermediary” experiments that utilize model versions between CMIP3/5 and CMIP6 as well as targeted sensitivity experiments in which individual modeling parameters are altered. We find clear and systematic improvements in the MJO and QBO and in the teleconnection patterns associated with the PDO and ENSO. Some gains arise from better process representation, while others (e.g., the QBO) from higher resolution that allows for a greater range of interactions. Our results demonstrate that the incremental development processes in multiple climate model groups lead to more realistic simulations over time.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the Impacts of Land-Cover Change on the Hydrologic Response to Hurricane Ida in the Lower Mississippi River Basin

Abstract The Lower Mississippi River basin (LMRB) has experienced significant changes in land cover and is one of the most vulnerable regions to hurricanes in the United States. Here, we study the impacts of land-cover change on the hydrologic response to Hurricane Ida in LMRB. By using an integrated surface–subsurface hydrologic model, Energy Exascale Earth System Model (E3SM) Land Model coupled with the three-dimensional ParFlow subsurface flow model (ELM-ParFlow), we simulate the effects of land-cover change on the flood volume and peak timing induced by rainfall from Hurricane Ida. The results show that land-cover changes from 1850 to 2015, which resulted in a smoother surface and less vegetation, exacerbated both flood peak time and volume induced by Hurricane Ida. The effects of land-cover changes can be decomposed into two mechanisms: a smoother surface routes more water faster to a watershed outlet and less vegetation allows more water to contribute to surface runoff. By comparing scenarios in which the two mechanisms were isolated, we found that changes in soil moisture due to vegetation cover change have more dominant effects on floods in the southern part and changes in Manning’s coefficient have the largest effect on floods in the northern part of the LMRB. The study provides important insights into the complex relationship between land-use, land-cover, and hydrologic processes in coastal regions.

54 ENVIRONMENTAL SCIENCES↗

Depletion, Chemical Reaction and Transport in High Burnup Nuclear Fuel

We have developed a formulation and computational model for oxygen transport in Light Water Reactor (LWR) uranium dioxide fuel. The overall model couples the burnup simulation isotopic composition with a thermochemistry model of the fuel phase, and oxygen transport using the driving forces from the thermodynamic calculations. The diffused oxygen is accounted for in the thermochemistry model in order to establish consistent material and thermodynamic conditions in the fuel undergoing burnup. The model has been implemented in the nuclear fuel performance code Bison. The formulation and the model have been successfully demonstrated on fuel burnup data from the open literature. The developed capability enables consideration of complex chemical composition of irradiated fuels beyond available burnup models in Bison.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Mercury Model of the Molly-G Fast Burst Reactor (FBR)

In order to design neutron irradiation experiments at a given nuclear reactor, one requires a representative neutron spectrum for the configuration of the reactor that will exist during the experiment. This usually implies modeling the operation of the reactor using a neutron transport code. Since the modeling of pulsed neutron reactors is quite complex, requiring coupled modeling of neutronics and thermo-structural response, the typical methodology employs analysis of the neutron spectrum which is obtained from a static keff calculation. While this time-independent Eigenfunction is not truly representative of that during a reactor pulse, we postulate that it is a reasonable representation for our purposes. The Nuclear Survivability (NS) Program at LLNL is considering fielding such experiments at the Molly-G (‘Molybdenum Godiva-II’) Fast Burst Reactor (FBR) at the White Sands Missile Range (WSMR). The goal of these experiments is to (a) test / proof / calibrate diagnostics, and (b) test the efficacy of neutron shield configurations, before fielding at other nuclear reactors, such as the Annular Core Research Reactor (ACRR) at SNL/NM. After attempting to secure either drawings or a copy of an existing MCNP model of the Molly-G FBR from various institutions, it was decided to develop a new model for use with the Mercury Monte Carlo transport code. The model described herein will be used to design these experiments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗