Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive models”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Impact Fracture and Fragmentation of Glass via the 3D Combined Finite-Discrete Element Method

A driving technical concern for the automobile industry is their assurance that developed windshield products meet Federal safety standards. Besides conducting innumerable glass breakage experiments, product developers also have the option of utilizing numerical approaches that can provide further insight into glass impact breakage, fracture, and fragmentation. The combined finite-discrete element method (FDEM) is one such tool and was used in this study to investigate 3D impact glass fracture processes. To enable this analysis, a generalized traction-separation model, which defines the constitutive relationship between the traction and separation in FDEM cohesive zone models, was introduced. The mechanical responses of a laminated glass and a glass plate under impact were then analyzed. For laminated glass, an impact fracture process was investigated and results were compared against corresponding experiments. Correspondingly, two glass plate impact fracture patterns, i.e., concentric fractures and radial fractures, were simulated. The results show that for both cases, FDEM simulated fracture processes and fracture patterns are in good agreement with the experimental observations. The work demonstrates that FDEM is an effective tool for modeling of fracture and fragmentation in glass.

42 ENGINEERING↗
DOE OSTI2025

Baseflow Identification via Explainable AI With Kolmogorov‐Arnold Networks

Abstract Hydrological models often involve constitutive laws that may not be optimal in every application. We propose to replace such laws with the Kolmogorov‐Arnold networks (KANs), a class of neural networks designed to identify symbolic expressions. We demonstrate KAN's potential on the problem of baseflow identification, a notoriously challenging task plagued by significant uncertainty. KAN‐derived functional dependencies of the baseflow components on the aridity index outperform their original counterparts; they demonstrate that water availability, rather than potential evapotranspiration, drives baseflow by constraining actual evapotranspiration under arid conditions. On a test set, they increase the Nash‐Sutcliffe efficiency (NSE) by 65%, decrease the root mean squared error by 29%, and increase the Kling‐Gupta efficiency by 34%. This superior performance is achieved while reducing the number of fitting parameters from three to two. Next, we use data from 378 catchments across the continental United States to refine the water‐balance equation at the mean‐annual scale. The KAN‐derived equations based on the refined water balance outperform both the current aridity index model, with up to a 105% increase in NSE, and the KAN‐derived equations based on the original water balance. While the performance of our model and tree‐based machine learning methods is similar, KANs offer the advantage of simplicity and transparency and require no specific software or computational tools. This case study focuses on the aridity index formulation, but the approach is flexible and transferable to other hydrological processes. Plain Language Summary Equations used in hydrologic model are often suboptimal, resulting in reduced prediction accuracy and efficiency. We implemented Kolmogorov‐Arnold networks (KAN), a machine learning algorithm for deriving symbolic formulations, to estimate groundwater recharge and showed that it outperforms an existing state‐of‐the‐art semi‐empirical formulation. In hydrology, Nash‐Sutcliffe efficiency (NSE), root mean squared error (RMSE), and Kling‐Gupta efficiency (KGE) are commonly used to evaluate model performance. Higher NSE and KGE values indicate better performance, while lower RMSE values are preferable. Our results show that NSE increased by 71%, RMSE decreased by 32%, and KGE improved by 25%. In addition, KAN identifies an optimal functional form and can be used to derive new analytical formulas using the prior knowledge. The KAN‐inspired equation outperformed the original formulation and reduced the fitting parameters. Furthermore, we refined the water‐balance equation at the mean‐annual scale and showed that, based on the new water‐balance equation, KAN can derive new formulations that are superior to the original aridity index formulations (up to 105% increase in NSE) and KAN‐derived equations based on the original water balance. These findings highlight the significant potential of KAN to advance the scientific understanding of a wide range of hydrologic processes. Key Points Kolmogorov‐Arnold networks (KANs) enhance interpretability of machine‐learned hydrological models KAN‐derived symbolic formulations outperform state‐of‐the‐art semi‐empirical aridity indices KAN‐identified functional form yields an analytical index with fewer fitting parameters and improved performance

baseflow↗
DOE OSTI2022

Machine learning in nuclear materials research

Nuclear materials are often demanded to function for extended time in extreme environments, including high radiation fluxes with associated transmutations, high temperature and temperature gradients, mechanical stresses, and corrosive coolants. They also have a wide range of microstructural and chemical makeups, resulting in multifaceted and often out-of-equilibrium interactions. Machine learning (ML) is increasingly being used to tackle these complex time-dependent interactions and aid researchers in developing models and making predictions, sometimes with better accuracy than traditional modeling that focuses on one or two parameters at a time. Conventional practices of acquiring new experimental data in nuclear materials research are often slow and expensive, limiting the opportunity for data-centric ML, but new methods are changing that paradigm. Here we review high-throughput computational and experimental data approaches, especially robotic experimentation and active learning that is based on Gaussian process and Bayesian optimization. We show ML examples in structural materials (e.g., reactor pressure vessel (RPV) alloys and radiation detecting scintillating materials) and highlight new techniques of high-throughput sample preparation and characterizations, and automated radiation/environmental exposures and real-time online diagnostics. Herein, this review suggests that ML models of material constitutive relations in plasticity, damage, and even electronic and optical responses to radiation are likely to become powerful tools as they develop. Finally, we speculate on how the recent trends of using natural language processing (NLP) to aid the collection and analysis of literature data, interpretable artificial intelligence (AI), and the use of streamlined scripting, database, workflow management, and cloud computing platforms that will soon make the utilization of ML techniques as commonplace as the spreadsheet curve-fitting practices of today.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On learning what to learn: Heterogeneous observations of dynamics and establishing possibly causal relations among them

Abstract Before we attempt to (approximately) learn a function between two sets of observables of a physical process, we must first decide what the inputs and outputs of the desired function are going to be. Here we demonstrate two distinct, data-driven ways of first deciding “the right quantities” to relate through such a function, and then proceeding to learn it. This is accomplished by first processing simultaneous heterogeneous data streams (ensembles of time series) from observations of a physical system: records of multiple observation processes of the system. We determine (i) what subsets of observables are common between the observation processes (and therefore observable from each other, relatable through a function); and (ii) what information is unrelated to these common observables, therefore particular to each observation process, and not contributing to the desired function. Any data-driven technique can subsequently be used to learn the input–output relation—from k-nearest neighbors and Geometric Harmonics to Gaussian Processes and Neural Networks. Two particular “twists” of the approach are discussed. The first has to do with the identifiability of particular quantities of interest from the measurements. We now construct mappings from a single set of observations from one process to entire level sets of measurements of the second process, consistent with this single set. The second attempts to relate our framework to a form of causality: if one of the observation processes measures “now,” while the second observation process measures “in the future,” the function to be learned among what is common across observation processes constitutes a dynamical model for the system evolution.

Sroczynski, David W.↗
DOE OSTI2025

Vegetation biogeography is a main source of uncertainty in modelling the land carbon cycle

The terrestrial biosphere exchanges a large amount of CO 2 with the atmosphere through photosynthesis and respiration, determining the magnitude of land carbon sink and consequently influencing the rate of global warming. The magnitudes of global photosynthesis and respiration, however, vary widely across models (100-200 PgC/year), constituting a key and persistent source of uncertainty in carbon cycle and climate modelling. Here, we argue that the uncertainty in the land carbon cycle modelling is largely attributable to the uncertainty in biogeography – the distribution of plant functional types (PFTs). Using an ensemble of dynamic global vegetation models (DGVMs), we find a strong dependence of total photosynthesis on total area for each PFT. The dependence allows us to reduce the spread of land carbon cycle estimates by ~75% using remote sensing-based PFT maps. We further find that 56 ± 21% of climate-driven changes in global photosynthesis modelled by DGVMs are caused by changes in PFT distribution in the last two decades. Our study identifies vegetation biogeography as a main controlling factor of uncertainty in land carbon cycle modelling and highlights the importance of biogeography-climate interactions in carbon cycle and climate studies.

Zhao, Ruiying [National Univ. of Singapore (Singap↗
DOE OSTI2025

Tensile Modeling PVC Gels for Electrohydraulic Actuators

Polyvinyl chloride (PVC)-dibutyl adipate (DBA) gels are a fascinating dielectric elastomer actuator showing promise in soft robotics. When actuated with high voltages, the gel deforms towards the anode. A recent application of PVC gels in electrohydraulic actuators motivates elastic and hyperelastic constitutive relationships for tensile loading modes. PVC gels with plasticizer-to-polymer weight ratios of 2:1, 4:1, 6:1, and 8:1 w/w were evaluated. PVC gels exhibit a linear elastic region up to 25% strain. The elastic modulus decreased with increasing plasticizer content from 288.8 kPa, 56.1 kPa, 24.7 kPa, to 11 kPa. Poisson’s ratio also decreased with increasing plasticizer content from 0.42, 0.43, 0.39, to 0.35. We suggest that the decrease in polymer concentration facilitates a weakly interconnected polymer network susceptible to chain slippage that hinders the network response, thus lowering Poisson’s ratio. Our work suggests that PVC gels can be treated as isotropic and incompressible for large strains and hyperelastic modeling; however, highly plasticized gels tend to act less incompressible at small strains. The power scaling law between the elastic modulus and plasticizer weight ratio showed high agreement, making the elastic modulus deterministic for any plasticizer content. The Neo–Hookean, Mooney–Rivlin, Yeoh, Gent, Ogden, and extended tube hyperelastic constitutive models are investigated. The Yeoh model shows the highest feasibility when evaluated up to 3.5 stretch, showing a maximum normalized root-mean-square-error of 6.85%. Together, these findings establish a constitutive basis for PVC-DBA gels, incorporating small strain elasticity, large strain non-linear behavior, and network analysis while providing suggestive insight into the network structure required for accurately modeling the EPIC.

Polymer Science↗
DOE OSTI2024

Simulating the effect of grain structure and porosity on creep for powder bed fusion 316H

This report details results from modeling studies focused on identifying the potential factors influencing creep performance of 316H stainless steel manufactured using Laser Powder Bed Fusion (LPBF), with a focus on experimentally-observed and predicted differences between the response of the AM material and conventionally manufactured wrought 316H. The studies presented here systematically look at the role of grain structure, porosity and texture in the anisotropic creep behavior of the material under conditions expected in high temperature advanced nuclear reactors. A physics-based Crystal Plasticity Finite Element modeling approach has allowed us to look at these factors in isolation and study their role in the material's long-term deformation behavior. An update to the precipitation model used in the constitutive framework for 316H is also presented here, with the aim of improving the accuracy of our results from this modeling effort. This work is a step in the direction of gaining mechanistic understanding of the long-term deformation behavior of LPBF 316H, helping us model its long term performance and reducing the qualification time for the material.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Magnesium abundances in cool metal-polluted white dwarfs

ABSTRACT The accretion of rocky material is responsible for the presence of heavy elements in the atmospheres of a large fraction of white dwarf stars. Those objects represent a unique opportunity to infer the bulk composition of exoplanetesimals. This chemical characterization requires the use of detailed atmosphere models to determine the elemental abundances at the photospheres of white dwarfs. In this work, we use a state-of-the-art model atmosphere code to reanalyse the first large survey of metal-polluted white dwarfs for which abundances are found for multiple elements. We show that the improved constitutive physics of our models lead to systematically higher Mg abundances than previous analyses. We find an average $\log \, {\rm Mg/Ca}$ number abundance ratio of 1.5. This value is significantly above the reference abundance for chondrites, which is expected as current diffusion models predict that for the cool helium-atmosphere white dwarfs of our sample, Mg should remain in the atmosphere longer than Ca. This helps resolve a recently identified Mg depletion problem, where the planetesimals accreted by white dwarfs were reported to be Mg-deficient compared to the expected composition of their planetary systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Modeling the viscoplastic behavior of a semicrystalline polymer

In this study, a complex constitutive relation is identified using inverse modeling with the nominal mechanical response as sole experimental input. The methodology is illustrated for a semicrystalline thermoplastic in the presence of strain localization at finite deformations. The experimental database includes cylindrical tensile bars, compression pins and round notched bars loaded at strain rates spanning up to five decades and temperatures below and above T g . The data is organized into a calibration set and a validation set. The response of tensile specimens is determined using finite element analyses and a two-phase constitutive relation for semicrystalline polymers that accounts for temperature- and rate-sensitive plastic flow, pressure-sensitivity, small-strain softening and large-strain orientational hardening of the amorphous phase, along with the evolution of crystallinity. The large number of constitutive parameters is identified using an optimization tool coupled with the finite element solver and the calibration set from experiments. The methodology is shown to be successful in predicting the response of round notched bars and replicating the effects of temperature and strain rate on the severity of necking in tensile bars. The proposed model identification strategy is both simple and effective in comparison with other elaborate methods that attempt to access intrinsic behavior directly from high-fidelity experimental measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2022

MULTISTEP AND CONTINUOUS PHYSICS-INFORMED NEURAL NETWORK METHODS FOR LEARNING GOVERNING EQUATIONS AND CONSTITUTIVE RELATIONS

In this work, we investigate the applicability and relative merit of discrete and continuous versions of physics-informed neural network (PINN) methods for learning unknown governing equations or constitutive relations in a nonlinear dynamical system. In the case of unknown dynamics, entire right-hand-side (RHS) equations of the ordinary differential equations are unknown. In the case of unknown constitutive relations, however, the RHS equations are known up to the specification of constitutive relations (that may depend on the state of the system). We use a deep neural network to model unknown governing equations or constitutive relations. The discrete PINN approach combines classical multistep discretization methods for dynamical systems with neural-network-based machine learning methods. On the other hand, the continuous versions utilize deep neural networks to minimize the residual function for the continuous governing equations. We use the case of a fedbatch bioreactor system to study the effectiveness of these approaches and discuss conditions for their applicability. Our results indicate that the accuracy of the trained neural network models is much higher for the cases where we only have to learn a constitutive relation instead of all dynamics. This finding corroborates the well-known fact from scientific computing that building as much structural information as is available into an algorithm can enhance its efficiency and/or accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

An experimental, theoretical, and kinetic modeling study of post-flame oxidation of ammonia

The post -flame oxidation rate of ammonia was investigated in a novel atmospheric pressure flow reactor at temperatures of 1280 ± 16 K and as a function of residence time and mixture composition (1-10% O 2 , dry and moist). The experimental results, as well as selected data from literature, were analyzed using an updated detailed chemical kinetic model. The medium temperature, very lean conditions enhance the importance of reactions of the nitroxyl (HNO) intermediate. High-level theory was used to calculate the rate constant for HNO + NH 2 , indicating that this step is significantly faster than values used in literature. Furthermore, a trajectory based approach was used to determine collision efficiencies for selected bath gases for HNO + M. The experimental results show that the NH 3 oxidation rate increases with temperature and O 2 concentration, while the presence of water vapor slightly inhibits reaction. Formation of NO and N 2 O was strongly promoted at higher levels of O 2 . Modeling results agreed well with the measurements, except at the lowest level of O 2 . The predicted oxidation rate of NH 3 was shown to result from a delicate balance between chain branching and terminating steps involving NH 2 , H 2 NO, and HNO. Recent theoretical work on reactions of these species by Klippenstein and coworkers and Stagni et al. was instrumental in improving modeling predictions. After initiation, NO reached a pseudo -steady-state level, where the pathways to NO were largely balanced by the NH 2 + NO reaction. Nitric oxide was partly oxidized to NO 2 , with the NH 2 + NO 2 reaction responsible for most of the N 2 O formation. Novelty and significance statement: This study provides the first detailed kinetic analysis of the lean postflame oxidation of ammonia, based on time -resolved flow reactor data in a novel reactor. In addition to the post -flame oxidation rate of ammonia, data for formation of NO and N 2 O were compared with modeling predictions. The medium temperature, very lean conditions enhance the importance of reactions of the HNO and H 2 NO intermediates. Inclusion in the model of results from recent high-level theoretical work, including present calculations for HNO + NH 2 and HNO + M, was crucial for capturing the observed behavior. It is argued that the post -initiation steady-state NH 3 oxidation rates constitute important data for model validation, along with ignition delays and laminar flame speeds.

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

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗
DOE OSTI2020

Influence of crystal structure on constitutive anisotropy of silica sand at particle-scale

In spite of recent advanced geotechnical research regarding the influence of particle-scale morphology and fabric on the constitutive behavior of silica sands, there is still no published literature that investigated the influence of crystal structure of sand particles on their constitutive behavior. The main constituent of silica sand particles is quartz mineral, which has a trigonal crystal system and hexagonal lattice. This paper investigates the influence of quartz crystal structure on the constitutive response of silica sand particles using 3D X-ray diffraction (3DXRD), synchrotron micro-computed tomography (SMT), and 3D finite element (FE) analysis. Here, unconfined uniaxial compression experiments on single particles of synthetic silica cubes exposed the notion of constitutive anisotropy caused by the crystal structure of quartz. Since the silica cubes have simple regular geometry (1 mm) 3 , they were loaded in different directions with respect to their crystal local planes. The crystal structure was also identified for a natural silica sand specimen composed of 2705 particles and loaded in confined uniaxial compression. Physical limitations of the experimental measurements made a strong case for the use of 3D FE analysis to study how the change in crystal local orientation of individual sand particles can fundamentally produce anisotropy in the constitutive response of natural sand specimens. The 3D FE analysis was validated to accurately model the crystal-based constitutive anisotropy in silica particles using the results of the silica cube experiments. The analysis was then executed on 3D FE meshes that closely matched the physical shape and fabric of the 2705 particles in the natural sand specimen. Both the experiments and 3D FE analysis showed that the crystal structure of quartz essentially causes a directional anisotropy in the constitutive behavior of silica particles.

58 GEOSCIENCES↗
DOE OSTI2023

Bentonite buffer under high temperature: laboratory experiments and coupled process modeling

Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.

Zheng, Liange↗
DOE OSTI2023

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗
DOE OSTI2023

A Physics-Based Data-Driven Approach for Modeling of Environmental Degradation in Elastomers

Abstract Elastomers are now commonly used in a number of industries, including aerospace, structure, transportation, shipbuilding, and automotive, due to their excellent workability, formability, and flexibility. During their activity, elastomers are subjected to harsh environmental conditions, which decreases their resilience. False predictions made early in their lives can have major financial and environmental implications. Elastomers’ performance and properties, such as strength, durability, and density, are influenced by chemical changes in these materials, known as degradation, which occurs over time. This process can alter the morphology of a polymer matrix as well as cause chain scission and cross-linking, resulting in different behaviors than that of the unaged material. To demonstrate the effect of thermaloxidative aging on the mechanical behavior of elastomers, several experimental and theoretical models have been proposed. In view of the large volume of experimental data available on micro-structural evolution in the course of aging, we propose a physics-based data-driven approach to overcome the shortcomings of both phenomenological and micro-mechanical models. This work presents a novel thermodynamically consistent, multiagent machine-learned model for predicting the constitutive behavior of cross-linked elastomers during environmental aging, such as thermo-oxidative and hydrolytic aging for various states of deformation. Single mechanism degradation changes the polymer matrix over time where it is causing chain scission, reduction of cross-links, and morphology change. To capture the idealized Mullins effect and permanent set due to the effect of single aging mechanisms on nonlinear mechanical responses of elastomers, we propose a data-driven model for simulating inelastic elements in a polymer matrix. By using a sequential order reduction, we were able to reduce the 3D stress-strain tensor mapping problem to a small number of super-constrained 1D mapping problems. To systematically classify such mapping problems into a few categories, an assembly of multiple replicated conditional neural network learning agents (L-agents) is used based on our recent work. Each category is represented by a different type of agent. The effect of deformation history, aging time, and aging temperature is captured by this model. The model is validated using a broad collection of data, ranging from our experimental results to data from the literature. In addition, thermodynamic consistency and frame independence are investigated. The most significant achievements of this model are its precision, simplicity, and prediction of inelasticity under various states of deformation. The model’s accuracy and simplicity make it a good option for commercial and industrial applications. Conveniently, due to the model modular nature, it can be expanded in the future to include viscoelasticity and non-isotropic formation for better precision.

Ghaderi, Aref↗
DOE OSTI2020

Physics-Informed Deep Neural Networks for Learning Parameters and Constitutive Relationships in Subsurface Flow Problems

In this paper, we present a physics informed deep neural network (DNN) method for estimating parameters and unknown physics (constitutive relationships) in partial differential equation (PDE) models. We use PDEs in addition to measurements to train DNNs to approximate unknown parameters and constitutive relationships as well as states. The proposed approach increases the accuracy of DNN approximations of partially known functions when a limited number of measurements is available and allows for training DNNs when no direct measurements of the functions of interest are available. We employ physics informed DNNs to estimate the unknown space-dependent diffusion coefficient in a linear diffusion equation and an unknown constitutive relationship in a non-linear diffusion equation. For the parameter estimation problem, we assume that partial measurements of the coefficient and states are available and demonstrate that under these conditions, the proposed method is more accurate than state-of-the-art methods. For the non-linear diffusion PDE model with a fully unknown constitutive relationship (i.e., no measurements of constitutive relationship are available), the physics informed DNN method can accurately estimate the non-linear constitutive relationship based on state measurements only. Finally, we demonstrate that the proposed method remains accurate in the presence of measurement noise.

42 ENGINEERING↗
DOE OSTI2024

Improved biomass feedstock materials handling and feeding engineering data sets, design methods, and modeling/simulation tools

Forest Concepts led a project to develop a set of tools which enable feedstock handing equipment designers to better understand and model the flowability of bulk particulate biomass materials. The work products included improved flowability mathematical models, data sets used to populate the models, and two new laboratory devices – a biomass-scale true cubical triaxial tester and a biomass-scale gas pycnometer. This work was funded in part by the US Department of Energy under contact DE-EE0008254.

09 BIOMASS FUELS↗