Learning continuum-scale models from micro-scale dynamics via Operator Regression.
Abstract not provided.
SEARCH · Engineering Papers
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.
Abstract not provided.
This paper studied the effects of void sizes on the failure behavior of adhesively-bonded materials under global shear via micro-scale computational modeling. The numerical results indicate that large void sizes with increased distances between them at the same void volume fraction can alleviate the reduced joint strength of a weaker adherend-adhesive interface caused by interfacial voids. However, this is not the case for an adherend-adhesive interface with the interfacial tensile strength being sufficiently higher than the adhesive strength, showing the negligible effect on the joint strength due to various void sizes and the amount of interfacial voids. This preliminary investigation provides insightful information for scaling up the void size in the macro-scale computational modeling of joints, and enhances the understanding of the micro-mechanical adhesion in adherend-adhesive interfaces with different levels of bonding.
Published by ECS. Two novel scale-bridging algorithms to model reaction-diffusion transport in porous media are presented. The algorithms are based on direct numerical simulations and couple the information of a micro-scale model, which accounts for the large field of view provided by micro X-ray computed tomography (X-ray CT), and a nano-scale model, which locally resolves transport in the fine structure extracted from nano X-ray CT. The micro-scale model is discretized in the through-plane direction into a 1D grid, where effective properties and internal boundaries are determined based on the results from the nano-scale model. The validated algorithms are used to examine transport of oxygen in precious group metal-free electrodes considering both zero- and first-order kinetics. Unlike conventional methods, the results show that the effective diffusivity is not a passive property but increases in regions where the reaction-rate coefficient is large. The proposed algorithms account for the multiscale coupling of reaction-diffusion transport and material microstructure, thus improving the predictions compared to conventional methods.
Due to a stress redistribution after the excavation of an underground tunnel for radioactive waste disposal, an Ed/DZ (excavation disturbed/damaged zone) will be generated in the near field of the opening, resulting in significant changes in the hydraulic and mechanical properties of the rock mass in the zone. Initially more or less randomly distributed hydrocarbons at grain boundaries in rock salt, which sometimes can only be observed with ultraviolet light, can then be mobilised and migrate at a potentially significant rate towards the excavation. Within the international cooperative project DECOVALEX 2019, the migration mechanism of such fluid inclusions in rock salt is being studied intensively. Specifically, a multi-scale modelling strategy has been developed. A macroscale coupled hydro-mechanical modelling of an underground excavation was performed to determine hydraulic and time-dependent deviatoric stress conditions, by taking into account the rock salt creep behaviour. Under the obtained macro-scale constraints, micro-scale modelling of a pathway dilation along halite grain boundaries was performed using different model strategies: a) coupled hydromechanical modelling with a consideration of hydraulic pressure-induced dilatant deformation, b) nonlinear dynamic model taking account of fluid migration, stress-dependent grain boundary wetting and shear-induced dilatancy of salt, and c) phase-field modelling of flow pathway propagation. The permeability increase resulting from the pathway dilation is estimated to be as high as two orders of magnitude. Based on the permeability determined, a series of pressure buildups measured from a borehole with a high hydrocarbon release rate, a total of 430 build-ups within a monitoring time of 938 days, can be simulated with a macro-scale compressible flow model accounting for different zones around the opening.
This report presents a workflow for advancing fuel performance modeling of SiC composite cladding for light-water reactors by linking microscale, experimental data-informed finite element analysis with rod-scale fuel performance codes such as BISON. The workflow uses X-ray computed tomography (XCT) to capture the actual geometry and processing-induced defects of as-fabricated SiC composite tube specimens, particularly porosity and wall-thickness variations, and converts the segmented XCT volumes into image-based finite element meshes for high fidelity structural analysis.
Various feedstocks from coal, petroleum and bio-based acrylonitrile were investigated for their ability to produce carbon fiber within DOE strength (250 ksi), modulus (25 Msi), elongation (>1%) and cost (≤$5/lb) goals for wider utilization in light-weight vehicles for the automotive industry. The overarching approach was to chemically identify and map the various feedstocks, understand the molecular transformation to produce carbon fiber precursors (mesophase for pitch materials of polyacrylonitrile polymer) and how these impact carbon fiber quality. Predictive models were built to relate the feedstocks and precursors to carbon fiber properties using atomistic and micro-scale modeling within ±15%. Using tow-level carbon fiber epoxy composites, carbon fibers tow-level predictions of elastic components of the carbon fibers were also successfully modeled within ±15%, using finite element analysis to confirm single filament carbon fiber properties. These models were validated with tow-level epoxy composite experiments from carbon fibers produced in this study.
Li-ion battery performance and degradation are strongly correlated with the electrode microstructures and can be modeled at different scales, each with their own limitations. Herein, we compare predictions achieved with a macro- and a micro-scale model, that is, respectively, neglecting or considering the microstructural heterogeneity of the composite electrodes, on virtual numerically generated and real microstructures. While both models are in relative agreement at the low charge rates, differences arise for fast charging scenarios and especially for the real, highly heterogenous, microstructures. The microscale model predicts that electrolyte concentration saturation and depletion, respectively, at the back of the cathode and of the anode are exacerbated, and that lithium plating occurs earlier for real microstructures. The present work also indicates that the mesh representation significantly impacts the microscale model predictions, and consequently that microscale models should add surface area as a parameter to consider explicitly surface roughness. This article is the first of a series, with subsequent entries further investigating in-plane heterogeneities, lithium plating, and the impact of microstructure representativity on model predictions.
This paper describes a new modeling framework for microscopic to reservoir-scale simulations of hydraulic fracturing and production. The approach builds upon a fusion of two existing high-performance simulators for reservoir-scale behavior: the GEOS code for hydromechanical evolution during stimulation and the TOUGH+ code for multi-phase flow during production. The reservoir-scale simulations are informed by experimental and modeling studies at the laboratory scale to incorporate important micro-scale mechanical processes and chemical reactions occurring within the fractures, the shale matrix, and at the fracture-fluid interfaces. These processes include, among others, changes in stimulated fracture permeability as a result of proppant behavior rearrangement or embedment, or mineral scale precipitation within pores and microfractures, at µm to cm scales. In our new modeling framework, such micro-scale testing and modeling provides upscaled hydromechanical parameters for the reservoir scale models. We are currently testing the new modeling framework using field data and core samples from the Hydraulic Fracturing Field Test (HFTS), a recent field-based joint research experiment with intense monitoring of hydraulic fracturing and shale production in the Wolfcamp Formation in the Permian Basin (USA). Below, we present our approach coupling the reservoir simulators GEOS and TOUGH+ informed by upscaled parameters from micro-scale experiments and modeling. We provide a brief overview of the HFTS and the available field data, and then discuss the ongoing application of our new workflow to the HFTS data set.
Abstract Heating, ventilation, and air-conditioning (HVAC) systems are usually an industry’s highest consumer of energy, most of which goes toward space cooling in buildings. Industrial energy-efficiency audits not only benefit manufacturers but also generate significant economic and environmental benefits to localities, states, and the nation. This article analyzes the micro- and macro scale impacts of implementing energy-efficient HVAC systems by integrating the industrial building energy data with the macroeconomic regional economic flow model. Micro-scale data include 10 years of historical energy, cost, and carbon dioxide savings achieved from energy-efficient HVAC implementation offered to manufacturers through industrial energy audits. The data were integrated into the macroeconomic modeling framework to illuminate the cascading regional economic impacts of implementing energy-efficient HVAC recommendations in manufacturing facilities. Results show that if recommendations had been implemented throughout all manufacturers in the region, $656 M energy costs would have been directly saved, 7.8 million metric tons of carbon dioxide emissions would have been avoided, and 4387 jobs could have been created, resulting in a total annual economic impact of $899 M stemming from direct, indirect, and induced impacts. The results offer insight into how industrial energy systems can be designed and provide models for how communities can accomplish a net-zero society.
We model temperature dynamics during Shear Assisted Proccess Extrusion (ShAPE), a solid phase process that plasticizes feedstock with a rotating tool and subsequently extrudes it into a consolidated tube, rod, or wire. Control of temperature is critical during ShAPE processing to avoid liquefaction, ensure smooth extrusion, and develop desired material properties in the extruded products. Accurate modeling of the complicated thermo-mechanical feedbacks between process inputs, material temperature, and heat generation presents a significant barrier to predictive modeling and process design. In particular, connecting micro-structural scale mechanisms of heat generation to macro-scale predictions of temperature can become computationally intractable. In this work we use a neural network (NN) model of heat generation to bridge this gap, by combining it with a simplified model of the temperature dynamics due to conduction and convection to capture the macro scale evolution of temperature. We inform the construction of the NN heat generation model using crystal plasticity simulations at the micro-structural scale to model the effects of process inputs on generation of heat. We achieved close fits of the temperature dynamics model to a diverse experimental data-set. Further, the relationships learned by the NN model between process inputs and heat generation showed qualitative agreement with those predicted by crystal plasticity simulations.
Implementation of High Pressure Die Casting (HPDC) Aluminum (Al) body structures for high volume electrified vehicles (EV) to improve electric efficiency remains a key strategy within many original equipment manufacturer (OEM)s. In addition to high strength for safety requirements, superior Self-Piercing Riveting (SPR) performance is demanded for HPDC Al alloys to be compatible with high volume SPR joining. In this work, it is proposed to extend and validate an existing Contractor finite element multiscale macro-micro modeling approach to quantify the influence of the microstructure of HPDC alloys on the fracture strain/displacement under 3-point bend and clinch testing. The success of this work will allow to replace solution treatment stage with low energy consumption heat treatment (HT) processes, or to design new non heat treatable (NHT) HPDC Al alloys to eliminate HT requirements. Ultimately, this project will facilitate the application of HPDC Al alloys for super-large vehicle structures to significantly reduce vehicle weight, and thus improving energy efficiency. The purpose of this project is to extend and validate an existing finite element code, which is based on the Contractor developed macro-micro multi-scale modeling approach, to numerically simulate the three-point bending and clinch test and study the influences of material microstructural characteristics and phase properties on the rivetability. The macro-micro modeling approach begins with a sample scale model and identify the location, which is mostly prone to failure, the deformation history of the boundaries of that location calculated will be used to drive a microstructure-based sub-models where the material microstructure and microscale properties are considered. Using this approach, the wrap-bending failure for two Al alloys are correctly predicted for the first time. This will start with phase I effort of building a framework of macro-micro three point bending test and clinch test of Al10SiMgMn HPDC alloy in the as-cast and T7 heat treated conditions. Those results will then be validated with experimental test results. The phase II effort will involve the utilization of the knowledge learned in phase I to establish the quantitative correlation between the microstructure characteristics and the riveting performance, which will be further used to guide the optimization of HPDC Al alloy microstructure using heat treatment process to achieve sufficient rivetability to join large thin-wall HPDC alloys.
This paper studied the effects of air void defects on the failure behavior of adhesively-bonded materials under global shear via micro-scale computational modeling. The numerical results indicated that interfacial void defects can largely facilitate interfacial debonding of a weaker adherend/adhesive interface under shear. However, this is not the case for a stronger adherend/adhesive interface, showing the reduction on the nominal shear strength of an adhesive joint is mainly due to internal void defects. The reduced bonding performance due to voids can be improved by leveraging an appropriate surface roughness. This preliminary investigation is a first step to better understand the micro-mechanics of interfacial failure in the adhesion of a roughened/patterned adherend via surface modification(s) and an adhesive, and also shows the importance of minimizing interfacial void defects in particular at a weaker adherend/adhesive interface via different techniques.
A microscale model is presented in this study to simulate electrode kinetics of the oxygen electrode in a solid oxide electrolyzer cell (SOEC). Two mixed ionic/electronic conducting structures are examined for the oxygen producing electrode in this work: single layer porous lanthanum strontium cobalt ferrite (LSCF), and bilayer LSCF/SCT (strontium cobalt tantalum oxide) structures. A yttrium-stabilized zirconia (YSZ) electrolyte separates the hydrogen and oxygen electrodes, as well as a gadolinium doped-ceria (GDC) buffer layer on the oxygen electrode side. Electrochemical reactions occurring at the two-phase boundaries (2PBs) and three-phase boundaries (3PBs) of single-layer LSCF and bilayer LSCF/SCT oxygen electrodes are modeled under various SOEC voltages with lattice oxygen stoichiometry as the key output. The results reveal that there exists a competition in electrode kinetics between 2PBs and 3PBs, but 3PBs are the primary reactive sites for single-layer LSCF oxygen electrode under high voltages. These locations experience the greatest oxygen stoichiometry variations and are therefore the most likely locations for dimensional changes. By applying an active SCT layer over LSCF, the 2PBs become activated to compete with the 3PBs, thus alleviating oxygen stoichiometry variations and reducing the likelihood of dimensional change. This strategy could reduce lattice structural expansion, proving to be valuable for electrode-electrolyte delamination prevention and will be the focus of future work.
Not provided.
Surface evolution due to exposure to harsh environments is of importance in many scientific and technological applications. In plasma-exposed materials, the surface receives charged particles from the plasma, leading to a series of processes that drive the system far from equilibrium and may lead to the severe deterioration of the surface properties. Although surface morphological changes are driven by atomic collisions taking place over picoseconds and nanometers, these processes result in mass loss and redistribution of matter over much larger length and time scales. This necessitates a multi-scale approach capable of capturing the range of processes linking primary atomic collision events with engineering-level surface geometry changes. In this paper, we develop a computational model to simulate the morphological evolution and effective erosion rate of micro-architected tungsten foams during low-energy plasma ion bombardment. Furthermore, the model acts on several length scales, with the energy and angular dependence of the sputtering yield of flat tungsten surfaces determined using the SRIM code based on the binary collision approximation. This information is introduced into a low-fluence, short-term Monte Carlo raytracing model, and further into a high-fluence, long-term particle transport model. In the latter, material particles representing billions of atoms are described in a digitized 3-D representation of the foam structure from X-ray tomography data, and are sputtered off and redeposited using the atomistic information. We show that the redeposition of sputtered atoms leads to partial self-healing in the bottom layers with a sharp reduction in the sputtering coefficient of low density foams, roughly 25% of the solid W value. This is in qualitative agreement with recent experiments on low porosity W structures. At high fluence, the foam structure degrades considerably as there are fewer ligaments available to recapture sputtered atoms and, consequently, the sputtering rate increases again.
Many material and structure failures are initiated at the interface between two distinctive materials, such as the debonding of interfaces of matrix/fiber of composites or bulk matrix/inclusion in metals at micro-scales. Accurate and robust modeling of these interface failures is critical to the design of novel material and structure with improved functionality. This paper presents a robust implementation of a discontinuous Galerkin finite element method in MOOSE based on the Incomplete Interior Penalty Galerkin scheme for modeling material interface behavior. This DG-based interface method naturally accommodates important extrinsic cohesive zone material models and is capable of modeling the failures of material interfaces across which material properties are highly distinctive. We demonstrate a good performance of the proposed method using single fiber and multiple fiber debonding problems.
Abstract not provided.
This paper revealed the detrimental effects of micro-scale air interfacial voids on the debonding at the interface of carbon-fiber-reinforced polyphthalamide (CFRPPA) and thermoset adhesive, representing a weak adherend-adhesive interface, within a dissimilar joint made of an aluminum alloy and a CFRPPA. The reduced lap shear strength of the joint, due to different void area fractions at the CFRPPA-adhesive interface, can be computationally described by using equivalent interfacial properties in the modeling to avoid the explicit modeling of the micro-scale interfacial voids. Such equivalent interfacial properties (e.g., interfacial normal strength, etc.) was found to have a non-linear relationship with respect to interfacial void area fraction as well as lap shear strength. This work has practical applications by utilizing equivalent interfacial properties for the analytical and/or computational design(s) of adhesively bonded joints.