Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical balances”

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 721 records · Page 40

Transport phenomena in electrolyte solutions: Nonequilibrium thermodynamics and statistical mechanics

The theory of transport phenomena in multicomponent electrolyte solutions is presented here through the integration of continuum mechanics, electromagnetism, and nonequilibrium thermodynamics. The governing equations of irreversible thermodynamics, including balance laws, Maxwell's equations, internal entropy production, and linear laws relating the thermodynamic forces and fluxes, are derived. Green–Kubo relations for the transport coefficients connecting electrochemical potential gradients and diffusive fluxes are obtained in terms of the flux–flux time correlations. Furthermore, the relationship between the derived transport coefficients and those of the Stefan–Maxwell and infinitely dilute frameworks are presented, and the connection between the transport matrix and experimentally measurable quantities is described. To exemplify the application of the derived Green–Kubo relations in molecular simulations, the matrix of transport coefficients for lithium and chloride ions in dimethyl sulfoxide is computed using classical molecular dynamics and compared with experimental measurements.

42 ENGINEERING↗

Quantifying Groundwater Response and Uncertainty in Beaver‐Influenced Mountainous Floodplains Using Machine Learning‐Based Model Calibration

Abstract Beavers ( Castor canadensis ) alter river corridor hydrology by creating ponds and inundating floodplains, and thereby improving surface water storage. However, the impact of inundation on groundwater, particularly in mountainous alluvial floodplains with permeable gravel/cobble layers overlain by a soil layer, remains uncertain. Numerical modeling across various floodplain structures considers topographic and sediment complexity and multidirectional flow, linking inundation to groundwater response. This study develops a model‐data integration workflow to address uncertainty in groundwater response to beaver‐induced inundations in a mountainous alluvial floodplain in the Upper Colorado River Basin. Uncertain factors include seasonal hydrologic dynamics, hydraulic conductivities, floodplain structures, and meteorological forcings. We employed an ensemble of groundwater models, based on geophysical and hydrologic data, with machine learning‐based calibration using a neural density estimator. This allowed us to quantify the vertical flux from the soil layer to the permeable gravel bed, the down‐valley underflow within the gravel bed, and their ratios. Results show a significant increase in the vertical flux relative to down‐valley underflow, from 2 during dry pond periods to 20 during wet periods, serving as an analogy for conditions without and with beaver ponds. The study highlights the influence of floodplain structure on groundwater storage, water balance, and water quality impacted by beaver ponds. A thick gravel bed layer, with a large down‐valley underflow, minimizes the effect of beaver‐induced inundation on water quality. We emphasize the need for field‐scale measurements of floodplain structure and improved characterization of evapotranspiration changes to reduce uncertainty in groundwater response. Plain Language Summary Beavers change the flow of water in river corridors by creating ponds, expanding wetlands, and flooding floodplains. This increases surface water area, promotes plant growth, and enhances biodiversity. However, the impact of this flooding on groundwater flow is not well understood, especially in mountainous areas with gravel layers where water moves easily beneath soil. In this study, we used numerical modeling to investigate how beaver ponds influence groundwater in a mountainous floodplain of the Upper Colorado River Basin. We adapted a machine learning method to validate our numerical models using multiple field data sets. Our findings show that beaver ponds significantly increase vertical water flow from the soil to the gravel during wet periods, compared to when the ponds are fully drained. The study also highlights the importance of floodplain structure in controlling both water flow in gravel layers along the river direction and vertical flow from the soil to the gravel with the presence of beavers. To reduce uncertainty in groundwater response, we emphasize the need for more field‐scale measurements of floodplain structure, hydraulic properties, and evapotranspiration changes. Key Points Floodplain structures and hydraulic conductivities are important for groundwater response with beaver ponds in mountainous floodplains Large down‐valley underflow in permeability‐stratified floodplains reduces beaver‐induced impacts on groundwater storage and water quality Machine learning‐based model calibration methods are effective for estimating posterior distributions of groundwater model parameters

Wang, Lijing↗

Annealing Controls Ultrafast Dynamics of Carrier Production in Organic Photovoltaics Incorporating a Nonfullerene Acceptor

Nonfullerene acceptors have emerged as leading candidates in organic photovoltaics pushing past some of the limitations of their fullerene counterparts with absorption extending in the visible and NIR, decreased air-sensitivity, and tunable energy levels. Significant progress has been made in demonstrating the potential of nonfullerene acceptor-based devices to reach high efficiencies; however, a photophysical and mechanistic understanding of charge generation in these nonfullerene acceptors lags behind. In particular, the effects of annealing on carrier production have not been previously examined. Here, we use transient absorption spectroscopy and atomic force microscopy to examine the effect of solvent vapor annealing on a perylene diimide-based nonfullerene acceptor, PDI-DPP-PDI. We find that when mixed with the donor PTB7-Th, the effect of solvent vapor annealing on PDI-DPP-PDI is to reduce geminate recombination, leading to balanced exciton generation and charge separation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of Tethered Polymers on Dynamics of Nanoparticles in Unentangled Polymer Melts

Polymer-tethered nanoparticles (NPs) are commonly added to a polymer matrix to improve the material properties. Critical to the fabrication and processing of such composites is the mobility of the tethered NPs. Here, we study the motion of tethered NPs in unentangled polymer melts using molecular dynamics simulations, which offer a precise control of the grafted chain length N g and the number z of grafted chains per particle. As Ng increases, there is a crossover from particle-dominated to tethered-chain-dominated terminal diffusion of NPs with the same z. The mean squared displacement of loosely tethered NPs in the case of tethered-chain-dominated terminal diffusion exhibits two subdiffusive regimes at intermediate time scales for small z. The first one at shorter time scales arises from the dynamical coupling of the particle and matrix chains, while the one at longer time scales is due to the participation of the particle in the dynamics of the tethered chains. The friction of loosely grafted chains in unentangled melts scales linearly with the total number of monomers in the chains, as the friction of individual monomers is additive in the absence of hydrodynamic coupling. As more chains are grafted to a particle, hydrodynamic interactions between grafted chains emerge. As a result, there is a nondraining layer of hydrodynamically coupled chain segments surrounding the bare particle. Outside the nondraining layer is a free-draining layer of grafted chain segments with no hydrodynamic coupling. The boundary of the two layers is the stick surface where the shear stress due to the relative melt flow is balanced by the friction between the grafted and melt chains in the interpenetration layer. The stick surface is located further away from the bare surface of the particle with higher grafting density.

42 ENGINEERING↗

Arbuscular mycorrhizal fungi equalize differences in plant fitness and facilitate plant species coexistence through niche differentiation

Mycorrhizal fungi are essential to the establishment of the vast majority of plant species but are often conceptualized with contradictory roles in plant community assembly. On the one hand, host-specific mycorrhizal fungi may allow a plant to be competitively dominant by enhancing growth. On the other hand, host-specific mycorrhizal fungi with different functional capabilities may increase nutrient niche partitioning, allowing plant species to coexist. Here, to resolve the balance of these two contradictory forces, we used a controlled greenhouse study to manipulate the presence of two main types of mycorrhizal fungus, ectomycorrhizal fungi and arbuscular mycorrhizal fungi, and used a range of conspecific and heterospecific competitor densities to investigate the role of mycorrhizal fungi in plant competition and coexistence. We find that the presence of arbuscular mycorrhizal fungi equalizes fitness differences between plants and stabilizes competition to create conditions for host species coexistence. Furthermore, our results show how below-ground mutualisms can shift outcomes of plant competition and that a holistic view of plant communities that incorporates their mycorrhizal partners is important in predicting plant community dynamics.

09 BIOMASS FUELS↗

Machine-learned closure of URANS for stably stratified turbulence: connecting physical timescales & data hyperparameters of deep time-series models

Stably stratified turbulence (SST), a model that is representative of the turbulence found in the oceans and atmosphere, is strongly affected by fine balances between forces and becomes more anisotropic in time for decaying scenarios. Moreover, there is a limited understanding of the physical phenomena described by some of the terms in the Unsteady Reynolds-Averaged Navier–Stokes (URANS) equations—used to numerically simulate approximate solutions for such turbulent flows. Rather than attempting to model each term in URANS separately, it is attractive to explore the capability of machine learning (ML) to model groups of terms, i.e. to directly model the force balances. We develop deep time-series ML for closure modeling of the URANS equations applied to SST. We consider decaying SST which are homogeneous and stably stratified by a uniform density gradient, enabling dimensionality reduction. We consider two time-series ML models: long short-term memory and neural ordinary differential equation. Both models perform accurately and are numerically stable in a posteriori (online) tests. Furthermore, we explore the data requirements of the time-series ML models by extracting physically relevant timescales of the complex system. We find that the ratio of the timescales of the minimum information required by the ML models to accurately capture the dynamics of the SST corresponds to the Reynolds number of the flow. The current framework provides the backbone to explore the capability of such models to capture the dynamics of high-dimensional complex dynamical system like SST flows.

97 MATHEMATICS AND COMPUTING↗

Genome-Scale Metabolic Modeling Enables In-Depth Understanding of Big Data

Genome-scale metabolic models (GEMs) enable the mathematical simulation of the metabolism of archaea, bacteria, and eukaryotic organisms. GEMs quantitatively define a relationship between genotype and phenotype by contextualizing different types of Big Data (e.g., genomics, metabolomics, and transcriptomics). In this review, we analyze the available Big Data useful for metabolic modeling and compile the available GEM reconstruction tools that integrate Big Data. We also discuss recent applications in industry and research that include predicting phenotypes, elucidating metabolic pathways, producing industry-relevant chemicals, identifying drug targets, and generating knowledge to better understand host-associated diseases. In addition to the up-to-date review of GEMs currently available, we assessed a plethora of tools for developing new GEMs that include macromolecular expression and dynamic resolution. Finally, we provide a perspective in emerging areas, such as annotation, data managing, and machine learning, in which GEMs will play a key role in the further utilization of Big Data.

59 BASIC BIOLOGICAL SCIENCES↗

Non-structural carbohydrates mediate seasonal water stress across Amazon forests

Non-structural carbohydrates (NSC) are major substrates for plant metabolism and have been implicated in mediating drought-induced tree mortality. Despite their significance, NSC dynamics in tropical forests remain little studied. We present leaf and branch NSC data for 82 Amazon canopy tree species in six sites spanning a broad precipitation gradient. During the wet season, total NSC (NSC T ) concentrations in both organs were remarkably similar across communities. However, NSC T and its soluble sugar (SS) and starch components varied much more across sites during the dry season. Notably, the proportion of leaf NSC T in the form of SS (SS:NSC T ) increased greatly in the dry season in almost all species in the driest sites, implying an important role of SS in mediating water stress in these sites. This adjustment of leaf NSC balance was not observed in tree species less-adapted to water deficit, even under exceptionally dry conditions. Thus, leaf carbon metabolism may help to explain floristic sorting across water availability gradients in Amazonia and enable better prediction of forest responses to future climate change.

54 ENVIRONMENTAL SCIENCES↗

Integrating Cambium Marginal Costs into Electric Sector Decisions: Opportunities to Integrate Cambium Marginal Cost Data into Berkeley Lab Analysis and Technical Assistance

NREL’s Cambium tool generates forward-looking simulations of marginal wholesale electricity costs associated with NREL’s Standard Scenarios. The scenarios include growing shares of variable renewable energy (VRE, i.e. wind and solar), among other power sector assumptions, between 2018 and 2050. The tool’s primary output—hourly costs at more than 130 balancing areas—could serve as public and transparent data source that supports electric-sector decision-making processes across the U.S. Berkeley Lab conducts a large range of analyses that use historical and forward-looking wholesale electricity prices to inform electric-sector decisions. In this report, Berkeley Lab uses its expertise to evaluate the Cambium cost data. We compare Cambium data with historical wholesale prices for the year 2018 and other modeled prices for the year 2030. We then present eight case studies in which Berkeley Lab researchers use Cambium data to replicate previous analyses based on other price datasets. We describe where primary findings and underlying key price dynamics align or differ, and highlight possible novel insights from the Cambium data. Finally, we qualitatively evaluate the suitability of Cambium costs in ten additional Berkeley Lab studies, though a direct comparison with alternative price data was not feasible at this time. The goal is to inform how electric-sector decision-makers and DOE program offices may be able to use this cohesive dataset, and to highlight what improvements to Cambium may make it even more useful.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Design and Production for Maximum Structural Efficiency With Respect to Fiber Orientation With Increased Understanding of Hybrid Fiber Flow Behavior

Discontinuous fiber-reinforced thermoplastic composites have gained considerable attention in automotive, aerospace, and other industries, due to their high-rate of production combined with their ability to attain complex and intricate shapes. Among other high-rate thermoplastic manufacturing processes, injection-molding is one of the most common manufacturing methods due to fast production and high surface finishing of complex geometries. Fiber orientation in discontinuous fiber composites plays a pivotal role in determining the mechanical, electrical, and thermomechanical properties, underscoring the necessity to comprehend fiber orientation in injection molded parts. Among different fiber types, glass and carbon fibers are most common in the composite industries. The recent trend of hybrid composites comprising both glass fiber (GF) and carbon fiber (CF) is also gaining importance in the automotive industry. Hybrid fiber options allow designers to optimize the balance between glass and carbon fibers by leveraging the high durability and low cost of GF while the strength and lightweight properties of CF. Consequently, comparing the fiber oriented distribution (FOD) of injection molded composites containing GF, CF, and a hybrid of GF/CF is critical to investigating the local mechanical properties of intricate structures for high-end applications. In Phase I of this project, FOD in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for GF, CF, and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varies with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗

Combined selection of the dynamic model and modeling error in nonlinear aeroelastic systems using Bayesian Inference

Here, we report a Bayesian framework for concurrent selection of physics-based models and (modeling) error models. We investigate the use of colored noise to capture the mismatch between the predictions of calibrated models and observational data that cannot be explained by measurement error alone within the context of Bayesian estimation for stochastic ordinary differential equations. Proposed models are characterized by the average data-fit, a measure of how well a model fits the measurements, and the model complexity measured using the Kullback–Leibler divergence. The use of a more complex error models increases the average data-fit but also increases the complexity of the combined model, possibly over-fitting the data. Bayesian model selection is used to find the optimal physical model as well as the optimal error model. The optimal model is defined using the evidence, where the average data-fit is balanced by the complexity of the model. The effect of colored noise process is illustrated using a nonlinear aeroelastic oscillator representing a rigid NACA0012 airfoil undergoing limit cycle oscillations due to complex fluid–structure interactions. Several quasi-steady and unsteady aerodynamic models are proposed with colored noise or white noise for the model error. The use of colored noise improves the predictive capabilities of simpler models.

42 ENGINEERING↗

An Acid-Free, Temperature-Based Cation Contamination Removal Strategy for PEM Water Electrolysis

It is widely understood that the durability and reliability of polymer electrolyte membrane (PEM) water electrolyzers are heavily dependent on feedwater purity, with cation contaminants that originate from incomplete water purification and balance of plant materials significantly harming electrolyzer performance. However, contamination remains a challenge and a common cause of failure at the stack level, indicating the need for strategies to recover the performance of contaminated cells. In this study, we investigate the effects of temperature on the uptake, electrochemical impacts, and removal of contaminant calcium and iron cations. Lower operating temperatures increase the sensitivity of the cell performance to contaminant cations, while also decreasing cation uptake and promoting contaminant removal. Computational charge transfer modelling shows that lower temperature increases the concentration of contaminant at the cathode and facilitates their removal from the cell. By testing single cells under scenarios designed to mimic stack temperature dynamics, we investigate low-temperature operation as an approach to stack-relevant contaminant recovery. Together, these results demonstrate that the low-temperature recovery approach is a promising approach for acid-free contamination recovery for PEM water electrolysis to promote stack reliability and durability.

08 HYDROGEN↗

DCFC + Hydrogen Station Design Optimization [Slides]

Both direct current fast charging (DCFC) and Hydrogen stations are working to create successful long-term business models; however, analyses for DCFC and hydrogen fueling infrastructure are almost always performed separately. This work provides a detailed exploration of the benefit of integrating DCFC and hydrogen stations to lower the total system cost from load balancing and equipment cost sharing. To achieve this we have adapted the REopt optimization framework to simultaneously optimize the design and operation of integrated DCFC and H2 fueling station. Results indicate that 1) combining hydrogen fueling and DCFC stations can significantly reduce lifetime costs compared to separated stations. 2) Co-location with additional site load reduces DCFC costs, however, integration of DCFC with hydrogen provides an even greater cost reduction. 3) Adding PV to combined stations further reduces the lifetime station cost. 4) Capital investments in station combination today can help reduce the cost of operating DCFC tomorrow and 5) product diversification acts as a hedge against variability and enables a more dynamic response to market changes.

08 HYDROGEN↗

ELM2.1-XGBfire1.0: improving wildfire prediction by integrating a machine learning fire model in a land surface model

Wildfires have shown increasing trends in both frequency and severity across the contiguous United States (CONUS). However, process-based fire models have difficulties in accurately simulating the burned area over the CONUS due to a simplification of the physical process and cannot capture the interplay among fire, ignition, climate, and human activities. The deficiency of burned area simulation deteriorates the description of fire impact on energy balance, water budget, and carbon fluxes in the Earth system models (ESMs). Alternatively, fire models based on machine learning (ML), which capture statistical relationships between the burned area and environmental factors, have shown promising burned area predictions and corresponding fire impact simulation. We develop a hybrid framework (ELM2.1-XGBFire1.0) that integrates an eXtreme Gradient Boosting (XGBoost) wildfire model with the Energy Exascale Earth System Model (E3SM) land model (ELM) version 2.1. A Fortran–C–Python deep learning bridge is adapted to support online communication between ELM and the ML fire model. Specifically, the burned area predicted by the ML-based wildfire model is directly passed to ELM to adjust the carbon pool and vegetation dynamics after disturbance, which are then used as predictors in the ML-based fire model in the next time step. Evaluated against the historical burned area from Global Fire Emissions Database 5 from 2001–2019, the ELM2.1-XGBFire1.0 outperforms process-based fire models in terms of spatial distribution and seasonal variations. The ELM2.1-XGBFire1.0 has proven to be a new tool for studying vegetation–fire interactions and, more importantly, enables seamless exploration of climate–fire feedback, working as an active component of E3SM.

54 ENVIRONMENTAL SCIENCES↗

Machine learning methods for particle stress development in suspension Poiseuille flows

Numerical simulations are used to study the dynamics of a developing suspension Poiseuille flow with monodispersed and bidispersed neutrally buoyant particles in a planar channel, and machine learning is applied to learn the evolving stresses of the developing suspension. The particle stresses and pressure develop on a slower time scale than the volume fraction, indicating that once the particles reach a steady volume fraction profile, they rearrange to minimize the contact pressure on each particle. Here we consider how the stress development leads to particle migration, time scales for stress development, and present a new physics-informed Galerkin neural network that allows for learning the particle stresses when direct measurements are not possible. The particle fluxes are compared with the Suspension Balance Model with good agreement. We show that when stress measurements are possible, the MOR-physics operator learning method can also capture the particle stresses.

97 MATHEMATICS AND COMPUTING↗

Elevating density functional theory to chemical accuracy for water simulations through a density-corrected many-body formalism

Abstract Density functional theory (DFT) has been extensively used to model the properties of water. Albeit maintaining a good balance between accuracy and efficiency, no density functional has so far achieved the degree of accuracy necessary to correctly predict the properties of water across the entire phase diagram. Here, we present density-corrected SCAN (DC-SCAN) calculations for water which, minimizing density-driven errors, elevate the accuracy of the SCAN functional to that of “gold standard” coupled-cluster theory. Building upon the accuracy of DC-SCAN within a many-body formalism, we introduce a data-driven many-body potential energy function, MB-SCAN(DC), that quantitatively reproduces coupled cluster reference values for interaction, binding, and individual many-body energies of water clusters. Importantly, molecular dynamics simulations carried out with MB-SCAN(DC) also reproduce the properties of liquid water, which thus demonstrates that MB-SCAN(DC) is effectively the first DFT-based model that correctly describes water from the gas to the liquid phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pulsed radio frequency power modulation for profile control in cryogenic reactive ion etching of silicon gratings

Cryogenic silicon etching with SF 6 /O 2 plasmas enables vertical sidewalls and smooth surfaces for high–aspect-ratio structures, yet precise control of etched profiles (sidewall angle, lateral etching, and linewidth retention) remains challenging under continuous-wave excitation due to coupling among radical fluxes, ion energy distributions, and in situ passivation dynamics. In this work, we investigated pulsed radio frequency bias delivery as a process knob for profile control during cryogenic silicon etching. Systematic experiments on 500 nm-wide, 1 μm-pitch gratings etched to 0.8–1.1 μm depth at 0.5–2 kHz pulsing frequency and 25%–75% duty cycle revealed strong dependencies of etch rate, linewidth, and sidewall angle on modulation conditions. We proposed a conceptual kinetic picture combined with feature-scale electrostatic analysis to understand the alternation of ion-driven etching and passivation recovery between on/off segments and local electric field shaping and modulation of ion trajectories within the trench. Moderate duty cycles (50%–75%) near 1 kHz achieved a balance between anisotropy and etch rate. These results demonstrate pulsed-bias modulation can be utilized to tune directionality and sidewall preservation in cryogenic SF 6 /O 2 etching, which are relevant to photonic, metrology, and quantum device fabrication.

36 MATERIALS SCIENCE↗

Simulated wildfire burned area over the CONUS during 2001-2020

Wildfires have shown increasing trends in both frequency and severity across the Contiguous United States (CONUS). However, process-based fire models have difficulties in accurately simulating the burned area over the CONUS due to a simplification of the physical process and cannot capture the interplay among fire, ignition, climate, and human activities. The deficiency of burned area simulation deteriorates the description of fire impact on energy balance, water budget, and carbon fluxes in the Earth System Models (ESMs). Alternatively, machine learning (ML) based fire models, which capture statistical relationships between the burned area and environmental factors, have shown promising burned area predictions and corresponding fire impact simulation. We develop a hybrid framework (ML4Fire-XGB) that integrates a pretrained eXtreme Gradient Boosting (XGBoost) wildfire model with the Energy Exascale Earth System Model (E3SM) land model (ELM). A Fortran-C-Python deep learning bridge is adapted to support online communication between ELM and the ML fire model. Specifically, the burned area predicted by the ML-based wildfire model is directly passed to ELM to adjust the carbon pool and vegetation dynamics after disturbance, which are then used as predictors in the ML-based fire model in the next time step. Evaluated against the historical burned area from Global Fire Emissions Database 5 from 2001-2020, the ML4Fire-XGB model outperforms process-based fire models in terms of spatial distribution and seasonal variations. Sensitivity analysis confirms that the ML4Fire-XGB well captures the responses of the burned area to rising temperatures. The ML4Fire-XGB model has proved to be a new tool for studying vegetation-fire interactions, and more importantly, enables seamless exploration of climate-fire feedback, working as an active component in E3SM.

Liu, Ye↗