SEARCH · Search NASA

Results for “Function Driven Optimization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Application of Koopman operator for model-based control of fracture propagation and proppant transport in hydraulic fracturing operation

This work explores the application of the recently developed Koopman operator approach for model identification and feedback control of a hydraulic fracturing process. Controlling fracture propagation and proppant transport with precision is a challenge due in large part to the difficulty of constructing approximate models that accurately capture the characteristic moving boundary and highly-coupled dynamics exhibited by the process. Koopman operator theory is particularly attractive here as it offers a way to explicitly construct linear representations for even highly nonlinear dynamics. The method is data-driven and relies on lifting the states to an infinite-dimensional space of functions called observables where the dynamics are governed by a linear Koopman operator. Here this work considers two problems: (a) fracture geometry control, and (b) proppant concentration control. In both cases, an approximate linear model of the corresponding dynamics is constructed and used to design a model predictive controller (MPC). The manuscript shows that in the case of highly nonlinear dynamics, as observed in the proppant concentration, use of canonical functions in the observable basis fails. In such cases, a priori system knowledge can be leveraged to choose the required basis. The numerical experiments demonstrate that the Koopman linear model shows excellent agreement with the real system and successfully achieves the desired target values maximizing the oil and gas productivity. Additionally, due to its linear structure, the Koopman models allow convex MPC formulations that avoid any issues associated with nonlinear optimization.

42 ENGINEERING↗

DOE OSTI2025

Evaluation of probability density function descriptions for three-component Rayleigh–Taylor mixing

Results from simulations of a three-component Rayleigh–Taylor (RT) mixing problem are presented. These simulations are conducted in heavy–light–heavy and heavy–intermediate–light configurations, and each of these configurations are further considered in high- and low-Reynolds-number regimes. This results in RT-unstable flow with one or both interfaces initially unstable, permitting the influence of problem configuration on the statistical description of three-component RT-driven mixing to be considered. Mass fraction covariances are observed to undergo a sign change through the mixing layer in all four configurations considered. This appears to be unique to the multi-component case and represents another way in which multi-component RT mixing differs from the two-component case. Qualitative and quantitative comparisons of joint and marginal probability density function (PDF) descriptions of species concentration are made. Three-, five-, and six-parameter model PDFs are compared against simulation data to assess how accurately they describe the mixing, and it is found that three-component mixing requires at least a five-parameter model PDF to accurately describe the mixing. Notably, the marginal distributions of three-component mixing do not appear to conform to a beta distribution, representing a departure from the classical two-component RT case. In conclusion, statistical neutrality also appears to influence the optimal choice of model PDF, which is found to be a function of problem configuration.

Large-eddy simulation↗

DOE OSTI2024

The Overlooked Potential of Sulfated Zirconia: Reexamining Solid Superacidity Toward the Controlled Depolymerization of Polyolefins

Closed-loop recycling via an efficient chemical process can help alleviate the global plastic waste crisis. However, conventional depolymerization methods for polyolefins, which compose more than 50% of plastics, demand high temperatures and pressures, employ precious noble metals, and yield complex mixtures of products limited to single-use fuels or oils. Superacidic sulfated zirconia (SZrO) with Hammet Acidity Functions (H 0 ) = -12 (i.e., stronger than 100% H 2 SO 4 ) is an industrially deployed heterogeneous catalyst capable of activating hydrocarbons under mild conditions and is shown to decompose polyolefins at ambient pressure and temperatures near 200 °C. Additionally, confinement of active sites in porous supports is known to radically increase selectivity, coking and sintering resistance, and acid site activity, presenting a possible approach to low-energy polyolefin depolymerization. However, a critical examination of the literature on SZrO led us to a surprising conclusion: despite 40 years of catalytic study, engineering, and industrial use, the surface chemistry of SZrO is poorly understood. Ostensibly spurred by SZrO's impressive catalytic activity, the application-driven study of SZrO has resulted in deleterious ambiguity in requisite synthetic conditions and insufficient characterization of acidity, porosity, and active site structure. Here, this ambiguity has produced significant knowledge gaps surrounding the synthesis, structure, and mechanisms of hydrocarbon activation for optimized SZrO, stunting the potential of this catalyst in olefin cracking and other industrially relevant reactions, such as isomerization, esterification, and alkylation. Toward mitigating these long extant issues, we herein identify and highlight these current shortcomings and knowledge gaps, propose explicit guidelines for characterization of and reporting on characterization of solid acidity, and discuss the potential of pore-confined superacids in the efficient and selective depolymerization of polyolefins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Using Particle Image Velocimetry to Determine Turbulence Model Parameters

The primary parameter of a standard k-ϵ model, C μ , was calculated from stereoscopic particle image velocimetry (PIV) data for a supersonic jet exhausting into a transonic crossflow. This required the determination of turbulent kinetic energy, turbulent eddy viscosity, and turbulent energy dissipation rate. Image interrogation was optimized, with different procedures used for mean strain rates and Reynolds stresses, to produce useful turbulent eddy viscosity fields. The eddy viscosity was calculated by a least-squares fit to all components of the three-dimensional strain-rate tensor that were available from the PIV data. This eliminated artifacts and noise observed when using a single strain component. Local dissipation rates were determined via Kolmogorov’s similarity hypotheses and the second-order structure function. The eddy viscosity and dissipation rates were then combined to determine C μ . Considerable spatial variation was observed in C μ , with the highest values found in regions where turbulent kinetic energy was relatively low but where turbulent mixing was important, e.g., along the high-strain jet edges and in the wake region. This suggests that use of a constant C μ in modeling may lead to poor Reynolds stress predictions at mixing interfaces. A data-driven modeling approach that can predict this spatial variation of C μ based on known state variables may lead to improved simulation results without the need for calibration.

42 ENGINEERING↗

DOE OSTI2023

Changing Directions: Influence of Ligand Electronics on the Directionality and Kinetics of Photoinduced Charge Transfer in Cu(I)Diimine Complexes

A key challenge to the effective utilization of solar energy is to promote efficient photoinduced charge transfer, specifically avoiding unproductive, circuitous electron-transfer pathways and optimizing the kinetics of charge separation and recombination. We hypothesize that one way to address this challenge is to develop a fundamental understanding of how to initiate and control directional photoinduced charge transfer, particularly for earth-abundant first-row transition-metal coordination complexes, which typically suffer from relatively short excited-state lifetimes. Here, in this paper, we report a series of functionalized heteroleptic copper(I)bis(phenanthroline) complexes, which have allowed us to investigate the directionality of intramolecular photoinduced metal-to-ligand charge transfer (MLCT) as a function of the substituent Hammett parameter. Ultrafast transient absorption suggests a complicated interplay of MLCT localization and solvent interaction with the Cu(II) center of the MLCT state. This work provides a set of design principles for directional charge transfer in earth-abundant complexes and can be used to efficiently design pathways for connecting the molecular modules to catalysts or electrodes and integration into systems for light-driven catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Layer Time Control for Large Scale Additive Manufacturing Using High Performance Computing

This work proposes to optimize an additive manufacturing AM process to reduce energy and printing cost. The polymer AM process is inherently dependent on the time-temperature history of each layer to maintain geometric tolerances and mechanical integrity. Our preliminary study shows that regression-based layer time control model using thermal images could result in up to 30% build time reduction for simple geometries. This proposed work would use high-performance computing (HPC) to couple the data-driven model with thermal simulation for better predicting layer temperature profiles, improving throughput of large-scale additive manufacturing, and reducing its energy cost. We have developed a method to optimize a layer deposition time (a.k.a. layer time) for large-scale AM via physics-based simulations. A long layer time leads to an over-cooled surface on which a new layer is deposited, and therefore, it may result in a weak bonding or debonding between layers, cracking, or warping. A short layer time leads to a high temperature of the structure due to insufficient cooling, and therefore, the structure may not be stiff enough and may collapse during manufacturing. Therefore, it is important to estimate the optimal layer time in additive manufacturing for a high-quality product. The temperature of a top layer right before deposition is recommended to be slightly higher than the glass temperature of the material. A temperature cooling was approximated to an exponential function of time, and the optimized layer time was obtained based on a target temperature while maintaining a minimal printing time. The material used is carbon fiber-reinforced polycarbonate (CF/PC), and the large-scale deposition system used is LSAM TM from Thermwood Corporation. Three different layer time cases were used for experiments, and a series of thermal images were obtained via an infra-red (IR) camera during the entire AM processes. AM process simulations were performed using a finite element method and the temperature profiles from the simulation were in good agreements with those from experiments. The layer time optimization was performed based on the temperature profiles from the simulations. A layer temperature with the optimal layer time was confirmed as the target temperature through simulation. In addition to the development of a layer time optimization method, we have developed a numerical framework for AM simulation with element activations in sync with toolpath, based on an open source finite element framework, DEAL.II. A major portion of this work was presented at SAMPE 2022 Conference and Exhibition on May 2022, and published in Proceedings of SAMPE 2022.

42 ENGINEERING↗

DOE OSTI2025

Information-entropy-driven generation of material-agnostic datasets for machine-learning interatomic potentials

In contrast to their empirical counterparts, machine-learning interatomic potentials (MLIAPs) promise to deliver near-quantum accuracy over broad regions of configuration space. However, due to their generic functional forms and extreme flexibility, they can catastrophically fail to capture the properties of novel, out-of-sample configurations, making the quality of the training set a determining factor, especially when investigating materials under extreme conditions. We propose a novel automated dataset generation method based on the maximization of the information entropy of the feature distribution, aiming at an extremely broad coverage of the configuration space in a way that is agnostic to the properties of specific target materials. The ability of the dataset to capture unique material properties is demonstrated on a range of unary materials, including elements with the FCC (Al), BCC (W), HCP (Be, Re and Os), graphite (C), and trigonal (Sb, Te) ground states. MLIAPs trained to this dataset are shown to be accurate over a range of application-relevant metrics, as well as extremely robust over very broad swaths of configurations space, even without dataset fine-tuning or hyper-parameter optimization, making the approach extremely attractive to rapidly and autonomously develop general-purpose MLIAPs suitable for simulations in extreme conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Expediting field-effect transistor chemical sensor design with neuromorphic spiking graph neural networks

Improving the sensitive and selective detection of analytes in a variety of applications requires accelerating the rational design of field-effect transistor (FET) chemical sensors. Achieving high-performance detection relies on identifying optimal probe materials that can effectively interact with target analytes, a process traditionally driven by chemical intuition and time-consuming trial-and-error methods. To address the difficulties in probe screening for FET sensor development, this work presents a methodology that combines neuromorphic machine learning (ML) architectures, specifically a hybrid spiking graph neural network (SGNN), with an enriched dataset of physicochemical properties through semi-automated data extraction using large language models. Achieving a classification accuracy of 0.89 in predicting sensor sensitivity categories, the SGNN model outperformed traditional ML techniques by leveraging its ability to capture both global physicochemical properties and sparse topological features through a hybrid modeling framework. Next-generation sensor design was informed by the actionable insights into the connections between material properties and sensing performance offered by the SGNN framework. Through virtual screening for the detection of per- and polyfluoroalkyl substances (PFAS) as a use case, the effectiveness of the SGNN model was further validated. Density functional theory simulations confirmed graphene as a promising active material for PFAS detection as suggested by the SGNN framework. By bridging gaps in predictive modeling and data availability, this integrated approach provides a strong foundation for accelerating advancements in FET sensor design and innovation.

Ferreira, Rodrigo Pires [Univ. of Chicago, IL (Uni↗

DOE OSTI2024

An integrated high-throughput robotic platform and active learning approach for accelerated discovery of optimal electrolyte formulations

Solubility of redox-active molecules is an important determining factor of the energy density in redox flow batteries. However, the advancement of electrolyte materials discovery has been constrained by the absence of extensive experimental solubility datasets, which are crucial for leveraging data-driven methodologies. In this study, we design and investigate a highly automated workflow that synergizes a high-throughput experimentation platform with a state-of-the-art active learning algorithm to significantly enhance the solubility of redox-active molecules in organic solvents. Our platform identifies multiple solvents that achieve a remarkable solubility threshold exceeding 6.20 M for the archetype redox-active molecule, 2,1,3-benzothiadiazole, from a comprehensive library of more than 2000 potential solvents. Significantly, our integrated strategy necessitates solubility assessments for fewer than 10% of these candidates, underscoring the efficiency of our approach. Our results also show that binary solvent mixtures, particularly those incorporating 1,4-dioxane, are instrumental in boosting the solubility of 2,1,3-benzothiadiazole. Beyond designing an efficient workflow for developing high-performance redox flow batteries, our machine learning-guided high-throughput robotic platform presents a robust and general approach for expedited discovery of functional materials.

25 ENERGY STORAGE↗

DOE OSTI2026

A Review on Direct Air Capture of Carbon Dioxide: Sorbent Materials, Process Engineering, Industrial Scale-Up, and Future Perspectives

The relentless accumulation of anthropogenic greenhouse gases has driven atmospheric carbon dioxide concentrations to approximately 426 ppm, necessitating the aggressive deployment of negative-emission technologies to achieve net zero by 2050. Direct air capture (DAC) offers a scalable, location-independent approach to atmospheric carbon removal; however, it is fundamentally constrained by the significant thermodynamic barriers associated with capturing CO 2 from ultradilute ambient conditions, requiring minimum thermodynamic energy inputs substantially higher than those for postcombustion point sources. This comprehensive review critically examines the technological landscape of DAC, focusing on the interdependent triad of sorbent material design, contactor engineering, and regeneration thermodynamics. We evaluate the fundamental boundaries of adsorption, emphasizing that an optimal adsorption enthalpy and isosteric heat of adsorption must balance the high CO 2 uptake capacity with the energetic penalties of sorbent regeneration. A systematic, comparative analysis of state-of-the-art sorbents is presented, encompassing mesoporous silicas, zeolites, carbon-based materials (CBMs), metal−organic frameworks (MOFs), porous organic polymers (POPs), and polymeric membranes. Special attention is devoted to surface functionalization strategies, particularly amine grafting and impregnation, which transition capture mechanisms from physisorption to chemisorption to enhance selectivity under ambient moisture and low partial pressures. Furthermore, we assess the operational merits of various reactor configurations, including gas−solid, gas−liquid, and membrane contactors, alongside regeneration cycles such as temperature, vacuum, pressure, and moisture swing adsorption. Finally, the review bridges fundamental materials science with industrial application by chronicling the scale-up milestones of pioneering entities and providing a strategic roadmap for advancing DAC technology readiness levels toward global deployment.

Adsorption↗

DOE OSTI2025

Process intensification approach to enhancing heat and mass transfer during drying: Ultrasonic (US) assisted drying of paper and board

Drying of paper and board is conventionally achieved through alternating conduction (steam-heated cylinders) and pocket convection (heated air over the paper web surface). These conventional drying systems rely heavily on steam from fossil fuels, resulting in inefficiencies, high energy usage, and thermal losses due to surface-driven mechanisms. Here, to address these challenges, an experimental system, with in-situ drying characteristics measurements, was developed to investigate process intensification using ultrasonic-based dewatering—a volumetric, pressure-driven acoustic energy system—integrated with conventional drying. The objectives of this study are to assess the impact of ultrasonics (US) on dewatering; compare performances to conventional drying systems; identify improvements in drying rate and energy use as a function of moisture content; and gain potential insights on heat and mass transfer mechanisms during US-assisted drying. US performance was evaluated across frequencies, power levels, pulp types, and basis weights. Results show that improvements to ultrasonic applications in conjunction with convection were 30-43% in drying rate and 20-35% in drying time over continuous and intermittent applications. When combined with conduction and convection, ultrasonics yielded up to 20% improvement in both rate and time and up to 20% reduction in energy consumption. Observations support a hypothesis of extension of the constant rate period due to improved capillary flow at higher moisture content and enhancing vapor diffusion and boundary layer disruption at lower moisture contents during falling rate period. These findings will inform future modeling, simulation, design and optimization of advanced drying systems.

42 ENGINEERING↗

DOE OSTI2023

Experimental Verification of An Electro-Hydraulic Actuation System Driven by An Integrated Electro-Hydraulic Unit

This paper proposes an electro-hydraulic actuator (EHA) system, and a novel-designed electro-hydraulic unit (EHU) consisting of a fixed-displacement hydraulic gear machine and a variable-speed electric machine. The novel EHU design integrates an electric machine and a hydraulic machine in a single housing, targeting optimal compactness, power density, and component reduction. The EHA system features an open circuit design, where hydraulic hoses connect the EHU with a tank, a valve manifold, and the hydraulic cylinder. In this way, the proposed EHA technology can be used to implement distributed hydraulic actuation in a vehicle without requiring changes to the hydraulic actuators or at the overall layout of the hydraulic components with respect to the original vehicle design. A dedicated test rig is developed to verify the performance of the proposed EHA system. Here, the efficiency of the EHA is measured in a steady state and under a realistic duty cycle of a commercial compact loader. The efficiency of the EHA system based on the measurements on the test rig can reach 54% with up to 20 kN load and 6 kW power level. Furthermore, the EHU is implemented on the reference machine to drive the boom function. Preliminary results on the performance of the EHU show an efficiency of up to 68% under different loading conditions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

High-throughput phase-field simulations and machine learning of resistive switching in resistive random-access memory

Metal oxide-based Resistive Random-Access Memory (RRAM) exhibits multiple resistance states, arising from the activation/deactivation of a conductive filament (CF) inside a switching layer. Understanding CF formation kinetics is critical to achieving optimal functionality of RRAM. Here a phase-field model is developed, based on materials properties determined by ab initio calculations, to investigate the role of electrical bias, heat transport and defect-induced Vegard strain in the resistive switching behavior, using MO 2- x systems such as HfO 2- x as a prototypical model system. It successfully captures the CF formation and resultant bipolar resistive switching characteristics. High-throughput simulations are performed for RRAMs with different material parameters to establish a dataset, based on which a compressed-sensing machine learning is conducted to derive interpretable analytical models for device performance (current on/off ratio and switching time) metrics in terms of key material parameters (electrical and thermal conductivities, Vegard strain coefficients). These analytical models reveal that optimal performance (i.e., high current on/off ratio and low switching time) can be achieved in materials with a low Lorenz number, a fundamental material constant. This work provides a fundamental understanding to the resistive switching in RRAM and demonstrates a computational data-driven methodology of materials selection for improved RRAM performance, which can also be applied to other electro-thermo-mechanical systems.

36 MATERIALS SCIENCE↗

DOE OSTI2021

A Hybrid Energy System Workflow for Energy Portfolio Optimization

This manuscript develops a workflow, driven by data analytics algorithms, to support the optimization of the economic performance of an Integrated Energy System. The goal is to determine the optimum mix of capacities from a set of different energy producers (e.g., nuclear, gas, wind and solar). A stochastic-based optimizer is employed, based on Gaussian Process Modeling, which requires numerous samples for its training. Each sample represents a time series describing the demand, load, or other operational and economic profiles for various types of energy producers. These samples are synthetically generated using a reduced order modeling algorithm that reads a limited set of historical data, such as demand and load data from past years. Numerous data analysis methods are employed to construct the reduced order models, including, for example, the Auto Regressive Moving Average, Fourier series decomposition, and the peak detection algorithm. All these algorithms are designed to detrend the data and extract features that can be employed to generate synthetic time histories that preserve the statistical properties of the original limited historical data. The optimization cost function is based on an economic model that assesses the effective cost of energy based on two figures of merit: the specific cash flow stream for each energy producer and the total Net Present Value. An initial guess for the optimal capacities is obtained using the screening curve method. The results of the Gaussian Process model-based optimization are assessed using an exhaustive Monte Carlo search, with the results indicating reasonable optimization results. The workflow has been implemented inside the Idaho National Laboratory’s Risk Analysis and Virtual Environment (RAVEN) framework. The main contribution of this study addresses several challenges in the current optimization methods of the energy portfolios in IES: First, the feasibility of generating the synthetic time series of the periodic peak data; Second, the computational burden of the conventional stochastic optimization of the energy portfolio, associated with the need for repeated executions of system models; Third, the inadequacies of previous studies in terms of the comparisons of the impact of the economic parameters. The proposed workflow can provide a scientifically defendable strategy to support decision-making in the electricity market and to help energy distributors develop a better understanding of the performance of integrated energy systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2022

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗

DOE OSTI2025

Dynamic Analysis of a Six-Cable Parallel Robot for Automated Panelized Building Retrofits

Cable-Driven Parallel Robots (CDPRs) are highly suitable for automated panelized building retrofits, thanks to their compact footprint and high payload-to-weight ratio. A common CDPR configuration featuring eight cables, where the anchors form a rectangular prism in front of the building facade, offers a large wrench feasible workspace and good control versatility. However, installing upper anchor points requires additional support structures, such as towers or beams, increasing setup complexity and posing logistical challenges in construction settings. To mitigate these challenges, we propose a six-cable CDPR model specifically designed for automated panelized building retrofits. Although the feasible workspace is limited, our analysis shows that the proposed CDPR adequately covers the critical areas required for panel installation. To validate that the six-cable system can effectively transport the end effector to the desired installation pose, we calculated optimal trajectories based on a constrained dynamic model. The simulation results of the six-cable CDPR demonstrate promising potential for automated panelized building retrofits, effectively balancing simplicity, cost-effectiveness, and functionality.

Liu, Yifang [ORNL] (ORCID:0000000190817417)↗

DOE OSTI2025

JAX-CanVeg: A Differentiable Land Surface Model

Land surface models consider the exchange of water, energy, and carbon along the soil-canopy-atmosphere continuum, which is challenging to model due to their complex interdependency and associated challenges in representing and parameterizing them. Differentiable modeling provides a new opportunity to capture these complex interactions by seamlessly hybridizing process-based models with deep neural networks (DNNs), benefiting both worlds, that is, the physical interpretation of process-based models and the learning power of DNNs. Here, we developed a differentiable land model, JAX-CanVeg. The new model builds on the legacy CanVeg by incorporating advanced functionalities through JAX in the graphic processing unit support, automatic differentiation, and integration with DNNs. We demonstrated JAX-CanVeg's hybrid modeling capability by applying the model at four flux tower sites with varying aridity. To this end, we developed a hybrid version of the Ball-Berry equation that emulates the water stress impact on stomatal closure to explore the capability of the hybrid model in (a) improving the simulations of latent heat fluxes (LE) and net ecosystem exchange (NEE), (b) improving the optimization trade-off when learning observations of both LE and NEE, and (c) benefiting a multi-layer canopy model setup. Our results show that the proposed hybrid model improved the simulations of LE and NEE at all sites, with an improved optimization trade-off over the process-based model. Additionally, the multi-layer canopy set benefited hybrid modeling at some sites. Anchored in differentiable modeling, our study provides a new avenue for modeling land-atmosphere interactions by leveraging the benefits of both data-driven learning and process-based modeling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Comprehensive Land Use and Environmental Stewardship

Comprehensive Land Use and Environmental Stewardship (CLUES) Report serves as a summary document of the land use and environmental stewardship activities occurring on the Idaho National Laboratory (INL) Site (Desert) and the Research and Education Campus (REC) within Idaho Falls. Land and facility use planning and decisions at the INL Site are guided by a comprehensive planning process in accordance with the United States (U.S.) Department of Energy (DOE) Order 430.1C, “Real Property Asset Management,” which states "Establish a data-driven, risk-informed, performance-based approach to the life-cycle management of real property assets that aligns the real property portfolio with DOE mission needs; acquire, manage, positively account for, and dispose of real property assets in a safe, secure, cost-effective, and sustainable manner; and ensure the real property portfolio is appropriately sized, aligned, and in the proper condition to support efficient mission execution." Land use planning, like Campus Master Planning and 5-Year Facility Planning, provides a means for better, more sustainable use of the INL Site in a coordinated effort to ensure current and future mission needs are met, including acquisition, recapitalization, maintenance, disposition, real property utilization and long-term stewardship. This document and all functions of INL are guided by the DOE Vision for INL and the INL Mission and builds on the baseline established in the FY 2015 CLUES Report. However, it delivers a revised structure with focus on new resource management zones, which provide organization of key management considerations and access restrictions to optimize land use and environmental stewardship. A new set of thematic Guiding Principles presented in this report provide the intent and sustainable management direction to protect the INL Site natural environment. They demonstrate that the mission and vision of DOE and INL can be realized with inclusion of first of a kind technology, private sector development, and globally recognized testing and demonstration. This CLUES Report encourages comprehensive management decision-making with additional resource discussions for sustainability, the built environment, and the ecological landscape at the INL Site. Enhanced resource considerations and trends emphasize the importance of air, land use, environmental, subsurface, and cultural resources. The INL Site supports exceptional and interdependent resources. A diversity of bat species are accommodated by culturally significant caves distributed around the Site. Big game species like elk, deer, and moose seek vegetation communities which are protected by long-term stewardship and monitoring programs.

99 GENERAL AND MISCELLANEOUS↗