Search NASA⌕ Search

SEARCH · Search NASA

Results for “root dynamics”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Connecting Nitrogen Transformations Mediated by the Rhizosphere Microbiome to Perennial Cropping System Productivity in Marginal Lands

The demand for energy from biofuel production is increasing, prompting concerns about the environmental impact and long-term sustainability of bioenergy cropping systems. These cropping systems will make up much of our future landscapes, and threaten to take the place of food cropping systems. Many life cycle analyses of bioenergy sustainability focus on carbon accrual and budgets, since they want to maximize carbon accrual while producing alternative fuel. Less attention has been given to nitrogen (N) dynamics in these systems. N is the most commonly limiting nutrient for plants, but applying nitrogen fertilizer- as we do for most cropping systems – is harmful to the environment, energetically costly, and produces greenhouse gases. In other words, adding nitrogen by fertilizer bioenergy systems could add to the very problems (climate change) it is trying to address. This is especially true for the areas that are proposed for bioenergy systems: marginal lands. These more degraded lands do not complete with food crops, but do have limited nitrogen. If we are to use these marginal lands for bioenergy, we need to understand the mechanisms regulating nutrient acquisition, and identify ways that bioenergy crops can get nitrogen in sustainable ways. Nutrient acquisition in the soil is performed by microbes in the root zone, or rhizosphere. Microbes can either mineralize nitrogen in the soil (from organic forms) or fix nitrogen from the air, in a process called nitrogen fixation. The goal of our project was thus to understand how the rhizosphere microbiome provides nutrients to bioenergy crops on marginal lands. We focus especially on the process of nitrogen fixation, since it has potential to get “fertilizer for free” that has much less environmental harm. We investigated this goal using sites from the DOE Great Lakes Bioenergy Research Center (GLBRC) in the upper Midwest, and associated lab and ‘omics methods. We group our findings into three major areas. First, we showed that nitrogen fixation, the conversion of N2 gas from the air to ammonium that is usable by plants, is performed in bioenergy soils, and benefits switchgrass crops. While more well-studied in leguminous plants, free-living nitrogen fixation can occur in some systems, and represents a potential opportunity to gain ‘free’ sustainable nitrogen source. We identified the nitrogen fixing bacteria that were most active in providing switchgrass with N, and showed that the drivers of nitrogen fixation occurred at a microscale; it is not well-predicted by bulk variables like soil moisture or plant phenology. Second, we showed that nitrogen fixation is not suppressed by long-term fertilizer. We expected that plentiful nitrogen would reduce the symbiotic relationship between nitrogen fixers and plants, and ‘downregulate’ fixation. We did not find evidence for this, either after long-term fertilizer in the field, or short-term fertilizer in the greenhouse. Finally, we identified the root exudates, carbon compounds that are emitted from the root, that best stimulate nitrogen fixation. We found that carbohydrates were better at stimulating fixation than organic acids. We expected these exudates to be emitted from the plant in periods of high N demand, but we found they are emitted when N is plentiful. This suggests that the stimulation of N fixation by plants is a passive process. Overall, we show that nitrogen fixation has potential to support bioenergy cropping system, and future management could develop ways to maximize it. However, this may not be best achieved via the plant – we found very little evidence of a ‘transactional’ system by which plants are controlling when and where nitrogen fixation is stimulated. It will be better to understand how management practices like planting and fertilizer application affect the microscale soil dynamics, which will determine nitrogen fixation rates.

59 BASIC BIOLOGICAL SCIENCES↗
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 OSTI2020

Disruption of Very-Long-Chain-Fatty Acid Synthesis Has an Impact on the Dynamics of Cellulose Synthase in Arabidopsis thaliana

In higher plants, cellulose is synthesized by membrane-spanning large protein complexes named cellulose synthase complexes (CSCs). In this study, the Arabidopsis PASTICCINO2 (PAS2) was identified as an interacting partner of cellulose synthases. PAS2 was previously characterized as the plant 3-hydroxy-acyl-CoA dehydratase, an ER membrane-localized dehydratase that is essential for very-long-chain-fatty acid (VLCFA) elongation. The pas2-1 mutants show defective cell elongation and reduction in cellulose content in both etiolated hypocotyls and light-grown roots. Although disruption of VLCFA synthesis by a genetic alteration had a reduction in VLCFA in both etiolated hypocotyls and light-grown roots, it had a differential effect on cellulose content in the two systems, suggesting the threshold level of VLCFA for efficient cellulose synthesis may be different in the two biological systems. pas2-1 had a reduction in both CSC delivery rate and CSC velocity at the PM in etiolated hypocotyls. Interestingly, Golgi but not post-Golgi endomembrane structures exhibited a severe defect in motility. Experiments using pharmacological perturbation of VLCFA content in etiolated hypocotyls strongly indicate a novel function of PAS2 in the regulation of CSC and Golgi motility. Through a combination of genetic, biochemical and cell biology studies, our study demonstrated that PAS2 as a multifunction protein has an important role in the regulation of cellulose biosynthesis in Arabidopsis hypocotyl.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Silica nanocluster binding rate coefficients from molecular dynamics trajectory calculations

Oxide nanoparticle growth from vapor phase precursors occurs in high temperature aerosol reactors first via the formation of nanoclusters, which are nanometer-scale condensed-phase species composed of 10 1 -10 2 atoms. The binding rate for nanoclusters, defined as the rate at which two nanoclusters collide with and stick to one another, is hence of interest in predicting nanoparticle size distribution functions in an aerosol. We have utilized molecular dynamics (MD) trajectory calculations to determine the homogeneous (equal-sized) and heterogeneous (disparate-sized) binding rate coefficients of SiO2 (silica) nanoclusters composed of 18, 144, and 333 atoms. MD trajectory calculations incorporated all-atom models of nanoclusters using a combined Born–Huggins–Mayer-Lennard-Jones potential model, which accounts for both short range interactions and electrostatic interactions. MD calculations were utilized to determine the binding probability as a function of both initial relative velocity and impact parameter; integration of the binding probability across impact parameter and relative velocity space yields the binding rate coefficient. MD trajectory calculations reveal that the most common type of encounters between nanoclusters leading to binding are grazing collisions, i.e. instances where collision would not occur without induced dipole potential influences. The resulting binding rate coefficients are found to be extremely weakly dependent on system temperature, which is in contrast to the use of rate coefficient models which are the product of a hard-sphere collision rate coefficient and a constant enhancement factor (leading to a rate coefficient proportional to the square root of temperature). Enhancement factors defined from MD trajectory calculations fall in the range of 3–9 as system temperature decreases from 1500 K to 300 K. Furthermore, such large values suggest that potential interactions need to be considered when modeling oxide nanocluster growth in the gas phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Self-consistent hopping theory of activated relaxation and diffusion of dilute penetrants in dense crosslinked polymer networks

We generalize a microscopic statistical mechanical theory of the activated dynamics of dilute spherical penetrants in glass-forming liquids to study the influence of crosslinking in polymer networks on the penetrant relaxation time and diffusivity over a wide range of temperature and crosslink fraction (f n ). Our calculations are relevant to recent experimental studies of a nm-sized molecule diffusing in poly-(n-butyl methacrylate) networks. The theory predicts the penetrant relaxation time increases exponentially with the glass transition temperature, T g (f n ), which grows roughly linearly with the square root of fn due to the coupling of local hopping to longer-range collective elasticity. Moreover, T g is also found to be proportional to a geometric confinement parameter defined as the ratio of the penetrant diameter to the mean network mesh size. The decoupling ratio of the penetrant and Kuhn segment alpha times displays a complex non-monotonic dependence on f n and temperature that is well collapsed based on the variable T g (f n )/T. A model for the penetrant diffusion constant that combines activated relaxation and entropic mesh confinement is proposed, which results in a significantly stronger suppression of mass transport with degree of effective supercooling than predicted for the penetrant alpha time. Furthermore, this behavior corresponds to a new network-based type of “decoupling” of diffusion and relaxation. In contrast to the diffusion of larger nanoparticles in high temperature rubbery networks, our analysis in the supercooled regime suggests that for the penetrants studied the mesh confinement effects are of secondary importance relative to the consequences of crosslink-induced slowing down of activated hopping of glassy physics origin.

Chemistry↗
DOE OSTI2020

Assessment of Low-Altitude Atmospheric Turbulence Models for Aircraft Aeroelasticity

We investigate the dynamic aeroelastic response of large but slow aircraft in low-altitude atmospheric turbulence. To this end, three turbulence models of increasing fidelity, namely, the one-dimensional von Kármán model, the two-dimensional Kaimal model and full three-dimensional wind fields extracted from large-eddy simulations (LES) are used to simulate ambient turbulence near the ground. Load calculations and flight trajectory predictions are conducted for a representative high-aspect-ratio wing aircraft, using a fully coupled nonlinear flight dynamics/aeroelastic model, when it operates in background atmospheric turbulence generated by the aforementioned models. Comparison of load envelopes and spectral content, on vehicles of varying flexibility, shows strong dependency between the selected turbulence model and aircraft aeroelastic response (e.g. 58% difference in the predicted magnitude of the wing root bending moment between LES and von Kármán models). This is mainly due to the presence of large flow structures at low altitudes that have comparable dimensions to the vehicle, and which despite the relatively small wind speeds within the Earth boundary layer, result in overall high load events for slow-moving vehicles. Results show that one-dimensional models that do not capture those effects provide fairly non-conservative load estimates and are unsuitable for very flexible airframe design.

atmospheric turbulence↗
DOE OSTI2025

Effects of shielding gas and spot reduction on humping suppression in high-speed laser welding

High-speed laser welding is a crucial process to enhance the production rate of fuel cell fabrication. However, exceeding a critical welding speed leads to defect formation such as humps and root concavities, which degrade weld quality and mechanical performance. This study employed shielding gas and spot reduction approaches to suppress these defects. The relationships between individual defect types and their influence on weld strength were systematically examined. Moreover, in situ high-speed synchrotron X-ray imaging was employed to comprehend the mechanisms of humping mitigation by shielding gas through analysis of instantaneous keyhole and melt pool dynamics. The observations reveal that shielding gas mitigates humping by shortening the molten pool through enhanced surface heat convection and by altering surface tension to prevent melt flow collisions. The effects of shielding gas and spot reduction were further incorporated into the dimensionless humping index (π h ) that describes the humping formation tendency. These findings provide mechanistic insights and quantitative guidance for defect control strategies in high-speed laser welding.

08 HYDROGEN↗
DOE OSTI2021

Comparative Analysis of ML Techniques for Data-Driven Anomaly Detection, Classification and Localization in Distribution System

High penetration of Distributed Energy Resources (DERs), fundamental load behavior changes, controllable loads, and significant increase in Electrical Vehicles (EVs) lead to complex dynamic behavior of the electric distribution system. Increasing number of components also means more measurements, more data and more data anomalies. Detecting, classifying and localizing these anomalies are important for situational awareness, and at the same time, very challenging given increasing complexity of the system. Highly accurate and high-resolution analytical techniques are needed to support anomaly detection, classification and localization (AD-C-L) for monitoring, root cause analysis and decision making. This paper provides comprehensive review and analysis of the existing spatio-temporal AD-C-L techniques within the distribution system. Challenges for specific problems in AD-C- L have been also discussed in this paper. Existing AD-C- L techniques have been categorized and synthesized for specific merits and limitations of multiple Machine Learning (ML) methodologies using common developed metrics of performance. The comparative analysis is summarized and presented with the open research challenges and path forward for future research needs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Modeling Plant Nutrient Acquisition Strategies Alters Projections of Carbon and Nitrogen Dynamics in Bioenergy Agroecosystems

Plant strategies to acquire nutrients from limited environments help shape ecosystem carbon (C) and nitrogen (N) cycling and response to environmental change. The effects of plant strategies on ecosystem dynamics are largely uncharacterized in bioenergy agroecosystems, where the impacts could determine bioenergy's ability to meet its sustainability goals of storing C and reducing N loss. We used FUN-BioCROP (Fixation and Uptake of Nitrogen-Bioenergy Carbon, Rhizosphere, Organisms and Protection), a plant–microbe interaction model of coupled plant nutrient uptake and soil organic matter decomposition, to simulate the effects of nutrient acquisition strategies on soil microbial activity and ecosystem nutrient cycling in bioenergy feedstocks miscanthus (Miscanthus × giganteus) and sorghum (Sorghum bicolor (L.) Moench). We examined the model's ability to reproduce the relative effects of belowground nutrient uptake on microbial activity using a reanalysis of empirical data showing that miscanthus root exudation provoked a larger soil microbial response than sorghum. From baseline model simulations, we found that the ability of miscanthus to retranslocate N resulted in higher N uptake at a lower C cost than the sorghum/soybean rotation and that soil C and N pools increased under perennial (miscanthus) and decreased under annual (sorghum/soybean) cultivation. The model also predicted that greater root exudation increased soil C accumulation, highlighting the role of roots in forming stable soil C. Overall, the baseline model was unable to reproduce field observations of miscanthus root exudation stimulating microbial activity more than sorghum. To improve the model, we updated the soil microbial parameters in miscanthus to have faster decomposition, a higher C/N ratio, and greater carbon use efficiency. These changes improved the simulated soil microbial response to miscanthus root exudation, supporting the hypothesis that miscanthus soils foster a microbial community that is more responsive to root exudation than that of sorghum.

09 BIOMASS FUELS↗
DOE OSTI2021

Long-term measurements in a mixed-grass prairie reveal a change in soil organic carbon recalcitrance and its environmental sensitivity under warming

Soil respiration, the major pathway for ecosystem carbon (C) loss, has the potential to enter a positive feedback loop with the atmospheric CO2 due to climate warming. For reliable projections of climate-carbon feedbacks, accurate quantification of soil respiration and identification of mechanisms that control its variability are essential. Process-based models simulate soil respiration as functions of belowground C input, organic matter quality, and sensitivity to environmental conditions. However, evaluation and calibration of process-based models against the long-term in situ measurements are rare. Here, we evaluate the performance of the Terrestrial ECOsystem (TECO) model in simulating total and heterotrophic soil respiration measured during a 16-year warming experiment in a mixed-grass prairie; calibrate model parameters against these and other measurements collected during the experiment; and explore whether the mechanisms of C dynamics have changed over the years. Calibrating model parameters against observations of individual years substantially improved model performance in comparison to pre-calibration simulations, explaining 79–86% of variability in observed soil respiration. Interannual variation of the calibrated model parameters indicated increasing recalcitrance of soil C and changing environmental sensitivity of microbes. Overall, we found that (1) soil organic C became more recalcitrant in intact soil compared to root-free soil; (2) warming offset the effects of increasing C recalcitrance in intact soil and changed microbial sensitivity to moisture conditions. Furthermore, these findings indicate that soil respiration may decrease in the future due to C quality, but this decrease may be offset by warming-induced changes in C cycling mechanisms and their responses to moisture conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Three-point energy correlators in hadronic Higgs boson decays

We present the analytic calculation of the leading-order three-point energy correlator (EEEC) in hadronic Higgs decays, including both the gluon-initiated channel H → g g + X and quark-initiated channel H → q q ¯ + X . The phase-space integration is evaluated directly using Mandelstam variables s i j = ( p i + p j ) 2 , and the appearing square roots can be rationalized by either conformal ratios or celestial coordinate variables. Throughout the calculation, we observe the same transcendental function space as in N = 4 super–Yang-Mills theory and e + e − → hadrons . Different infrared limits are also explored using the full analytic result, offering the fixed-order data for EEEC factorization and resummation. Given its nontrivial shape dependence, the EEEC presents an excellent opportunity to explore the dynamics of gluon jets originating from the H → g g decay channel at future lepton colliders. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2020

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗
DOE OSTI2024

Panama Rainforest Changes with Experimental Drought (PaRChED): Initial Effects of Partial Throughfall Exclusion on Soil Dynamics in Lowland Forests Across Variation in Rainfall and Soil Fertility

Changes in rainfall are predicted across tropical regions, with effects on nutrient, water, and carbon cycling. This chapter summarizes results from the first two years of a throughfall exclusion experiment in four lowland Panamanian forests that span a 1,000-mm change in rainfall and variation in soil fertility. Soil respiration (i.e., soil carbon dioxide [CO2 ] flux) declined with throughfall exclusion, with a site*season interaction, and the radiocarbon age of respired carbon was older in exclusion versus control plots. The decline in soil CO2 flux could be related to reduced fine root production and soil microbial biomass. Microbial community composition also changed with throughfall exclusion in infertile soils, and soil nutrients accumulated more in exclusion versus control plots during the dry season. The net effects on soil carbon storage will depend on the relative strengths of these effects over time. Continued research could improve predictions of tropical forest-climate feedbacks with changes in precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Parked aeroelastic field rotor response for a 20% scaled demonstrator of a 13‐MW downwind turbine

Abstract Aeroelastic parked testing of a unique downwind two‐bladed subscale rotor was completed to characterize the response of an extreme‐scale 13‐MW turbine in high‐wind parked conditions. A 20% geometric scaling was used resulting in scaled 20‐m‐long blades, whose structural and stiffness properties were designed using aeroelastic scaling to replicate the nondimensional structural aeroelastic deflections and dynamics that would occur for a lightweight, downwind 13‐MW rotor. The subscale rotor was mounted and field tested on the two‐bladed Controls Advanced Research Turbine (CART2) at the National Renewable Energy Laboratory's Flatiron Campus (NREL FC). The parked testing of these highly flexible blades included both pitch‐to‐run and pitch‐to‐feather configurations with the blades in the horizontal braked orientation. The collected experimental data includes the unsteady flapwise root bending moments and tip deflections as a function of inflow wind conditions. The bending moments are based on strain gauges located in the root section, whereas the tip deflections are captured by a video camera on the hub of the turbine pointed toward the tip of the blade. The experimental results are compared against computational predictions generated by FAST, a wind turbine simulation software, for the subscale and full‐scale models with consistent unsteady wind fields. FAST reasonably predicted the bending moments and deflections of the experimental data in terms of both the mean and standard deviations. These results demonstrate the efficacy of the first such aeroelastically scaled turbine test and demonstrate that a highly flexible lightweight downwind coned rotor can be designed to withstand extreme loads in parked conditions.

17 WIND ENERGY↗
DOE OSTI2025

Modeling the mechanisms of coastal vegetation dynamics and ecosystem responses to changing water levels

Coastal forests are increasingly experiencing mortality due to inundation by fresh- and seawater, leading to their replacement by marshes. These shifts alter vegetation composition, biogeochemical cycling, carbon storage, and hydrology. Using a hydraulically enabled ecosystem demography model (FATES-Hydro), we conducted numerical experiments to investigate the mechanisms behind inundation-driven forest loss and the ecosystem-scale consequences of forest-to-marsh transitions. We compared mortality processes and their effects across broadleaf and conifer trees at two coastal sites – Lake Erie (freshwater) and Chesapeake Bay (saline). Our simulations show that hydraulic failure, driven by root loss under prolonged flooding, is the primary mortality mechanism across both tree types and sites. Forest replacement by marsh reduced ecosystem-scale leaf area index (LAI), gross primary production (GPP), transpiration, and deep soil water uptake in conifer forests, while broadleaf forests experienced smaller changes due to lower initial LAI and greater marsh compensation. Marsh invasion occurred following canopy thinning driven by tree mortality. These findings suggest that, under similar root loss, hydraulic failure dominates coastal tree mortality regardless of species or water type, with denser forests experiencing stronger ecosystem impacts. Our study identifies key mortality mechanisms and offers testable hypotheses for future empirical studies on coastal vegetation change.

Ding, Junyan [Occidental College, Los Angeles, CA ↗
DOE OSTI2024

Development of Steady-State and Dynamic Mass and Energy Constrained Neural Networks for Distributed Chemical Systems Using Noisy Transient Data

The paper presents the development of algorithms for mass and energy constrained neural network models that can exactly conserve the overall mass and energy of distributed chemical process systems, even though the noisy transient data used for optimal model training violate the same. In contrast to approximately satisfying mass and energy balance constraints of a system by soft penalization of objective function, algorithms have been developed for solving equality-constrained nonlinear optimization problems, thus providing the guarantee of exactly satisfying the system mass and energy conservation laws. For developing dynamic mass-energy constrained network models for distributed systems, hybrid series and parallel dynamic-static neural networks have been leveraged. The developed algorithms for solving both the training and forward problems are validated using both steady-state and dynamic data in the presence of various noise characteristics. The developed data-driven algorithms are flexible to exactly satisfy mass and energy balance constraints for dynamic chemical processes if the system holdup information is available. The proposed network structures and algorithms are applied to the development of data-driven lumped and distributed models of an adiabatic superheater/reheater system, a nonisothermal continuous stirred tank reactor, as well as an electrically heated plug-flow reactor system where one form of energy gets transformed to another. It has been observed that the mass-energy constrained neural networks yield a root mean squared error of <1% with respect to the system truth for the case studies evaluated in this work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Kinetics of Magnetic Skyrmion Crystal Formation from the Conical Phase

The particle-like magnetic skyrmion or skyrmion lattice (SkX) formation has promoted strong application and fundamental science interests. Despite extensive research, the kinetic of the SkX development is much less understood because of the ultrafast spin rotation and high sensitivity to external perturbations. In this work, using in situ Lorentz transmission electron microscopy, we successfully measured the dynamics of SkX formation from the conical phase with precise control of both the temperature and the magnetic field. We discovered that the Avrami equation can accurately describe the transition process with an initial Avrami constant around 1, suggesting that the rate-limiting step for the quasiparticle lattice formation is one-dimensional heterogeneous nucleation of individual skyrmions. A modified Arrhenius rate law is established, with an energy barrier that has a square-root dependence on temperature and a quadratic dependence on the magnetic field. This study paves the way toward precise and predictable manipulation of topological spin structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Topographical Controls on Hillslope-Scale Hydrology Drive Shrub Distributions on the Seward Peninsula, Alaska

Observations indicate shrubs are expanding across the Arctic tundra, mainly on hillslopes and primarily in response to climate warming. However, the impact topography exerts on hydrology, nutrient dynamics, and plant growth can make untangling the mechanisms behind shrub expansion difficult. We examined the role topography plays in determining shrub expansion by applying a coupled transect version of a mechanistic ecosystem model (ecosys) in a tundra hillslope site in the Seward Peninsula, Alaska. Modeled biomass of the dominant plant functional types agreed well with field measurements (R 2 = 0.89) and accurately represented shrub expansion over the past 30 years inferred from satellite observations. In the well-drained crest position, canopy water potential and plant nitrogen (N) uptake was modeled to be low from plant and microbial water stress. Intermediate soil water content in the mid-slope position enhanced mineralization and plant N uptake, increasing shrub biomass. The deciduous shrub growth in the mid-slope position was further enhanced by symbiotic N 2 fixation primed by increased root carbon allocation. The gentle slope in the poorly drained lower-slope position resulted in saturated soil conditions that reduced soil O 2 concentrations, leading to lower root O 2 uptake and lower nutrient uptake and plant biomass. Additionally, a simulation that removed topographical interconnectivity between grid cells resulted in (1) a 28% underestimate of mean shrub biomass and (2) over or underestimated shrub productivity at the various hillslope positions. Our results indicate that land models need to account for hillslope-scale coupled surface and subsurface hydrology to accurately predict current plant distributions and future trajectories in Arctic ecosystems.

54 ENVIRONMENTAL SCIENCES↗