Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “function as a service”

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 415 records · Page 23

Comprehensive Margin Assessment of the ASME Section III, Division 5, Class A Primary Load Design Rules

This report provides a comprehensive margin assessment of the ASME Section III, Division 5, Subsection HB, Subpart B rules for the design of high temperature reactor components against load controlled stress limits. These rules, often called the Code primary load design rules, provide protection against creep rupture and plastic collapse under steady conditions. The method adopted here is to compare the deterministic ASME design life, according to the primary load design rules, against the expected, statistical service life of the component, considering creep rupture as the relevant failure mode. The report provides the design margin of the ASME rules in terms of the probability of premature failure -- the probability that the component fails in service before reaching the ASME design life. To complete this assessment, the report describes the development of a statistical creep life assessment procedure accounting for variations in the component loading, material creep rate, and the material rupture time. This includes a novel method for correlating creep rupture data using a Gaussian process accounting for heat-to-heat variation in the Larson-Miller time-temperature parameter. The report applies the complete margin assessment process to quantify the design margin inherent in the ASME primary load design rules as a function of temperature for relevant component geometries and materials. The overall conclusion is that the ASME rules are conservative, providing designs with very small probability of premature failure, but some optimization of the Code design rules is possible to provide a more uniform design margin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Deep Learning-Based Weather-Related Power Outage Prediction with Socio-Economic and Power Infrastructure Data

This paper presents a deep learning-based approach for hourly power outage probability prediction within census tracts encompassing a utility company's service territory. Two distinct deep learning models, conditional Multi-Layer Perceptron (MLP) and unconditional MLP, were developed to forecast power outage probabilities, leveraging a rich array of input features gathered from publicly available sources including weather data, weather station locations, power infrastructure maps, socio-economic and demographic statistics, and power outage records. Given a one-hour-ahead weather forecast, the models predict the power outage probability for each census tract, taking into account both the weather prediction and the location's characteristics. The deep learning models employed different loss functions to optimize prediction performance. Our experimental results underscore the significance of socio-economic factors in enhancing the accuracy of power outage predictions at the census tract level.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The Straightening of a River Meander Leads to Extensive Losses in Flow Complexity and Ecosystem Services

To assist river restoration efforts we need to slow the rate of river degradation. This study provides a detailed explanation of the hydraulic complexity loss when a meandering river is straightened in order to motivate the protection of river channel curvature. We used computational fluid dynamics (CFD) modeling to document the difference in flow dynamics in nine simulations with channel curvature (C) degrading from a well-established tight meander bend (C = 0.77) to a straight channel without curvature (C = 0). To control for covariates and slow the rate of loss to hydraulic complexity, each of the nine-channel realizations had equivalent bedform topography. The analyzed hydraulic variables included the flow surface elevation, streamwise and transverse unit discharge, flow velocity at streamwise, transverse, and vertical directions, bed shear stress, stream function, and the vertical hyporheic flux rates at the channel bed. The loss of hydraulic complexity occurred gradually when initially straightening the channel from C = 0.77 to C = 0.33 (i.e., the radius of the channel is three-times the channel width), and additional straightening incurred rapid losses to hydraulic complexity. Other studies have shown hydraulic complexity provides important riverine habitat and is positively correlated with biodiversity. This study demonstrates how hydraulic complexity can be gradually and then rapidly lost when unwinding a river, and hopefully will serve as a cautionary tale.

54 ENVIRONMENTAL SCIENCES↗

Combined Spatially Resolved Characterization of Antireflection and Antisoiling Coatings for PV Module Glass

Characterization of photovoltaic (PV) module materials throughout different stages of service life is crucial to understanding and improving the durability of these materials. Currently the large-scale of PV modules (>1 m 2 ) is imbalanced with the small-scale of most materials characterization tools (= 1 cm 2 ). Furthermore, understanding degradation mechanisms often requires a combination of multiple characterization techniques. Here, we present adaptations of three standard materials characterization techniques to enable mapping characterization over moderate sample areas (=25 cm 2 ). Contact angle, ellipsometry, and UV-Vis spectroscopy are each adapted and demonstrated on two representative samples: a commercial multi-functional coating for PV glass and an oxide combinatorial sample library. Best practices are discussed for adapting characterization techniques for large area mapping and combining mapping information from multiple techniques.

14 SOLAR ENERGY↗

Recycling of natural fiber composites: Challenges and opportunities

Natural fibers have been widely used for reinforcing polymers attributed to their sustainable nature, excellent stiffness to weight ratio, biodegradability, and low cost compared with synthetic fibers like carbon or glass fibers. Thermoplastic composites offer an advantage of recyclability after their service life, but challenges and opportunities remain in the recycling of natural fiber reinforced polymer composites (NFRPCs). This article summarized the effects of reprocessing/recycling on the material properties of NFRPCs. The material properties considered include mechanical performance, thermal properties, hygroscopic behavior, viscoelasticity, degradation, and durability. NFRPCs can generally be recycled approximately 4–6 times until their thermomechanical properties change. After recycling 7 times, the tensile strength of NFRPCs can decrease by 17%, and the tensile modulus can decrease by 28%. The mitigation approaches to overcome degradation of material properties of NFRPCs such as adding functional additives and virgin plastics are also discussed. The main challenges in these approaches such as degradation and incompatibility are discussed, and an effort is made to provide a rationale for reprocessing/recyclability assessment. Finally, future applications of NFRPCs such as additive manufacturing and automotive part use are discussed.

36 MATERIALS SCIENCE↗

The Sensitivity of Hyporheic Exchange to Fractal Properties of Riverbeds

Hyporheic exchange in riverbeds is driven by current-bed topography interactions. Because riverbeds exhibit topographic roughness across scales, from individual grains to bedforms and bars, they can exhibit fractal patterns. This study analyzed the influence of fractal properties of riverbed topography on hyporheic exchange. A set of synthetic fractal riverbeds with different scaling statistics was used as inputs to sequentially coupled numerical simulations of turbulent channel flow and hyporheic flow. In the analysis, the maximum power spectrum (dune size) and the fractal dimension (topographic complexity) were considered as independent variables and we then investigated how interfacial fluxes and hyporheic travel times are functionally related to these variables. As the maximum power spectrum increases (i.e., dune height to flow depth ratio), the average interfacial flux increases logarithmically whereas it increases exponentially with an increase in fractal dimension. Hyporheic exchange is more sensitive to additional roughness (larger fractal dimensions) than to bedform size (larger maximum power). Our results imply that fractal properties of riverbeds are crucial to predicting hyporheic exchange. The predictive relationships we propose could be integrated with reduced complexity, large-scale models. Further, they can also be used to design artificial topographies that target hyporheic ecosystem services.

54 ENVIRONMENTAL SCIENCES↗

Modeling Framework for Bulk Electric Grid Impacts from HEMP E1 and E3 Effects (Tasks 3.1 Final Report)

This report presents a framework to evaluate the impact of a high-altitude electromagnetic pulse (HEMP) event on a bulk electric power grid. This report limits itself to modeling the impact of EMP E1 and E3 components. The co-simulation of E1 and E3 is presented in detail, and the focus of the paper is on the framework rather than actual results. This approach is highly conservative as E1 and E3 are not maximized with the same event characteristics and may only slightly overlap. The actual results shown in this report are based on a synthetic grid with synthetic data and a limited exemplary EMP model. The framework presented can be leveraged and used to analyze the impact of other threat scenarios, both manmade and natural disasters. This report d escribes a Monte-Carlo based methodology to probabilistically quantify the transient response of the power grid to a HEMP event. The approach uses multiple fundamental steps to characterize the system response to HEMP events, focused on the E1 and E3 components of the event. 1) Obtain component failure data related to HEMP events testing of components and creating component failure models. Use the component failure model to create component failure conditional probability density function (PDF) that is a function of the HEMP induced terminal voltage. 2) Model HEMP scenarios and calculate the E1 coupled voltage profiles seen by all system components. Model the same HEMP scenarios and calculate the transformer reactive power consumption profiles due to E3. 3) Sample each component failure PDF to determine which grid components will fail, due to the E1 voltage spike, for each scenario. 4) Perform dynamic simulations that incorporate the predicted component failures from E1 and reactive power consumption at each transformer affected by E3. These simulations allow for secondary transients to affect the relays/protection remaining in service which can lead to cascading outages. 5) Identify the locations and amount of load lost for each scenario through grid dynamic simulation. This can be an indication of the immediate grid impacts from a HEMP event. In addition, perform more detailed analysis to determine critical nodes and system trends. 6) To help realize the longer-term impacts, a security constrained alternating current optimal power flow (ACOPF) is run to maximize critical load served. This report describes a modeling framework to assess the systemic grid impacts due to a HEMP event. This stochastic simulation framework generates a large amount of data for each Monte Carlo replication, including HEMP location and characteristics, relay and component failures, E3 GIC profiles, cascading dynamics including voltage and frequency over time, and final system state. This data can then be analyzed to identify trends, e.g., unique system behavior modes or critical components whose failure is more likely to cause serious systemic effects. The proposed analysis process is demonstrated on a representative system. In order to draw realistic conclusions of the impact of a HEMP event on the grid, a significant amount of work remains with respect to modeling the impact on various grid components.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Utility-Scale Shared Energy Storage Deployment: Challenges, Research Gaps, and Opportunities

Although community energy storage (CES) and behind-the-meter (BTM) energy storage systems have been widely used to offer homeowners and communities a variety of localized benefits, their scalability and grid support functionalities are limited. On the other hand, utility-scale shared energy storage (USES) systems may offer a number of benefits for grid integration, scalability, and economic viability. When compared to BTM and CES alternatives, these large-scale systems provide more storage capacity, more efficient operations, and more economically viable options. The deployment of USES presents opportunities for optimizing grid performance, integrating renewable energy resources, and improving energy security at the community level. However, significant research gaps exist in optimizing the integration and operation of these systems, especially to allocate energy for consumer use, grid services, and enhancing energy resilience. This paper reviews the literature in this regard, focusing on the opportunities, research gaps, and challenges associated with USES deployment. Firstly, the paper provides an overview of USES systems and emphasizes their benefits. Secondly, the key challenges are identified, and research gaps associated with the operation and integration of these systems are highlighted. Lastly, some potential solutions and opportunities that can be adopted to facilitate the rapid deployment and management of USES are presented. Technological, economic, regulatory, and environmental aspects are also discussed in this paper, providing an overview of the current state and future prospects of this technology.

Gautam, Mukesh [BATTELLE (PACIFIC NW LAB)] (ORCID:↗

Effects of hurricane disturbance and increased temperature on carbon cycling and storage of a Puerto Rican forest: a mechanistic investigation of above- and belowground processes (Final Technical Report)

The overall goal of the Tropical Responses for Altered Climate Experiment (TRACE) is to assess the effects of increasing temperature on tropical plant and soil carbon fluxes and storage as the forest recovers from major hurricane disturbance that occurred in September 2017. Ultimately, we aimed to reduce uncertainty and increase confidence with which tropical forests are represented in Earth System models to make more accurate global forecasts of future climate. We focused on both above- and belowground processes and explored temperature controls over critical aspects of carbon and nutrient cycling for tropical plants, soil, and microbes. TRACE is located in a wet tropical forest in the Luquillo Experimental Forest close to the USDA Forest Service Sabana Field Research Station in Luquillo, Puerto Rico. The warming treatment consists of six 4.7 m diameter plots. Three of the plots receive infrared warming and three have the same infrastructure but are not warmed (using ‘dummy’ heaters). Each plot was monitored from 2018-2023 to investigate two major questions: 1) Are there legacy effects of prior warming on forest recovery following hurricane disturbance? 2) Will the trajectory of forest recovery following disturbance be affected by warmer temperatures? Concurrent soil incubation experiments were conducted to enable more controlled mechanistic investigations of temperature response on microbial function. In sum, our goal was to use this novel climate manipulation experiment (the only of its kind in any tropical forest) and once-in-a-lifetime chance to assess how temperature and hurricane disturbance interact to affect coupled biogeochemical cycling in a tropical forest.

54 ENVIRONMENTAL SCIENCES↗

Plutonium-238 Loadout Calorimeter Return To Service

The calorimeter in the transfer area in Building 7920 at the US Department of Energy’s Oak Ridge National Laboratory is used to measure the amount of 238 Pu in the Plutonium Heat Source Oxide. The calorimetry measurement is taken after the oxide is packaged in an ORISC (Oak Ridge Inner Shipping Capsule)/FSO (Fuel Storage Outer) assembly before loading into the Containment Vessel and shipping container. The equipment is used in three specific configurations: zero measurement (nothing in the sample cup), FSO measurement (actual measurement of 238 Pu), and reference calibration measurement (using the reference FSO as a known heat source). System checks were performed in July 2023 in preparation for scheduled material loadout. During these checks, personnel discovered that the reference FSO had experienced an electrical short in the heater cable and was no longer functional. The calorimeter was in the reference calibration measurement configuration at this time. Facility management was notified when the damage was discovered. The data from the Yokogawa recorder were pulled and analyzed for a span of time between September 2022 and July 2023. An unusual spike was found in the sample cup temperature data that started on October 19, 2022, and ended on October 20, 2022. The sample cup reached a temperature of 191.1°C. The maximum allowable temperature for the calorimeter thermopile to reach is 220°C based on the manufacturer statement provided in Appendix B of this report [1]. The calorimeter thermopile and the sample cup thermocouple are physically close to each other; therefore, the data suggest that the calorimeter itself never saw a temperature in excess of 191.1°C. Physical damage present on the cable, tape, and insulation suggests that the temperature reached a higher value in some locations within the sample cup than the data shows. Because of this discrepancy, testing of the calorimeter was required to confirm that the calorimeter is functioning the same as it was before the electrical short. The intent of this report is to assess the functionality of the calorimeter and determine if it is acceptable to use for 238 Pu loadout material measurement. Testing showed that the calorimeter is functionally equivalent to the initial calibration that was performed in the summer of 2021, and thus, the calorimeter was cleared to be used for its intended purposes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Anisotropic thermal diffusivity and conductivity in SiC/SiC tubes studied by infrared imaging and X-ray computed tomography

Traditional flash diffusivity evaluation of thermal diffusivity/conductivity of composite tubes require machining of specimens. For a thin-wall tube, this method can only be used to obtain through-thickness transport property. A novel method to evaluate anisotropic thermal diffusivity in a composite tube has been developed. Braided SiC/SiC composite tubes were subjected to a xenon flash heating pulse. A high-speed, high-sensitivity infrared camera was used to measure surface temperature changes as a function of time and nondestructively detect subsurface defects/damages, such as macroscopic pores. Standard reference material (Pyroceram 9606) and curved SiC/SiC composite tube specimens were used to validate thermal diffusivity obtained from infrared imaging. Unlike the traditional method, there is no need prepare special specimens, and thermal diffusivity values in three orientations are obtained after a single flash. A finite element analysis model based on x-ray computed tomography scans was developed to simulate the heat transfer. In conclusion, this technique is significant in assessing thermal conductivity and inspecting the health of ceramic tubes during and after service.

36 MATERIALS SCIENCE↗

Degradation of impact toughness in cast stainless steels during long-term thermal aging

Cast austenitic stainless steels (CASSs) have been extensively used for the large components of nuclear reactor primary coolant systems. Since the cast steels inevitably contain degradable metastable phases and replacement of the large coolant system components is impractical, the thermal embrittlement of CASS components has been a serious concern in the extended-term operation of nuclear power plants. This study aimed to systematically measure and analyze the effect of long-term thermal aging on the Charpy impact toughness to provide a comprehensive understanding of thermal degradation behavior and a practical aging model to predict the degree of thermal degradation in the cast stainless steels. Additionally, the materials tested in the research include eight CASS alloys (two CF3s, one CF3M, three CF8s, and two CF8Ms) and two reference wrought materials (304L and 316L), in which the nominal δ-ferrite content ranges from ~2% to 33%. These stainless steels have been thermally aged at two light water reactor (LWR) temperatures (290 and 330 °C) and at two accelerated-aging temperatures (360 and 400 °C) for up to 30,000 h; these include both under-aged and over-aged conditions relative to the extended service lifetime (80 years). Charpy impact testing was performed for aged and non-aged specimens, and the impact (absorbed) energy parameters were correlated with a new aging parameter (A). Both the reduction of impact fracture toughness and the shift of ductile-brittle transition temperature were strongly dependent on the δ-ferrite content and degree of thermal aging. A linear relationship was found between the increasing rate of the index transition temperature T 41J and aging parameter A; base on which an empirical model was proposed for prediction of the transition temperature as a simple function of the aging parameter (A) and δ-ferrite content (F). Finally, the critical aging parameter for embrittlement (AC) was evaluated and compared with the existing δ-ferrite content criteria.

36 MATERIALS SCIENCE↗

Grid-scale Energy Storage Hazard Analysis & Design Objectives for System Safety

Battery based energy storage systems are becoming a critical part of a modernized, resilient power system. However, batteries have a unique combination of hazards that can make design and engineering of battery systems difficult. This report presents a systematic hazard analysis of a hypothetical, grid scale lithium-ion battery powerplant to produce sociotechnical "design objectives" for system safety. We applied system's theoretic process analysis (STPA) for the hazard analysis which is broken into four steps: purpose definition, modeling the safety control structure, identifying unsafe control actions, and identifying loss scenarios. The purpose of the analysis was defined as to prevent event outcomes that can result in loss of battery assets due to fires and explosions, loss of health or life due to battery fires and explosions, and loss of energy storage services due to non- operational battery assets. The STPA analysis resulted in identification of six loss scenarios, and their constituent unsafe control actions, which were used to define a series of design objectives that can be applied to reduce the likelihood and severity of thermal events in battery systems. These design objectives, in all or any subset, can be utilized by utilities and other industry stakeholders as "design requirements" in their storage request for proposals (RFPs) and for evaluation of proposals. Further, these design objectives can help to protect firefighters and bring a system back to full functionality after a thermal event. We also comment on the hazards of flow battery technologies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Multiphysics Meshfree Degradation Modeling of Energy Storage Materials with Kernel Enrichment

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and ultimately diminishing performance and service life. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), pixel-based meshfree model construction by the reproducing kernel particle method (RKPM) is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. The first kernel enrichment discussed will be the interface modified reproducing kernel (IM-RK) [1, 2], constructed by scaling a smooth kernel function with an interface-distance function to achieve strategic discontinuity types (i.e. weak discontinuities for strain discontinuities and strong discontinuities for cracks) and alleviate Gibbs oscillations near these transition zones. The IM-RK is especially useful for areas in which a known discontinuity-type is expected a priori. The second kernel enrichment to be discussed is a neural network-enhanced reproducing kernel (NN-RK) [3, 4], which is introduced to effectively model non-obvious damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RK is additionally used to inform how crack opening and closure in turn affect the electro-chemo-mechanical responses in the material microstructure. Reference: [1] Wang, Y., Baek, J., Tang, Y. et al. "Support vector machine guided reproducing kernel particle method for image-based modeling of microstructures," Comput Mech 73, 907-942 (2024). https://doi.org/10.1007/s00466-023-02394-9. [2] Susuki, K., Allen, J. & Chen, J. S.. "Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method," Engineering with Computers (2024). https://doi.org/10.1007/s00366-024-02016-9. [3] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, 4422-4454 (2022). https://doi.org/10.1002/nme.7040.

25 ENERGY STORAGE↗

Correction to Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

This errata is to correct the acknowledgment in the original published article. G.B. (partial support), V.-A.G., M.-T.N., and R.R. were supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division. Computer resources were provided by Research Computing at Pacific Northwest National Laboratory and the National Energy Research Scientific Computing Center (NERSC), a U.S. Department of Energy Office of Science User Facility operated under Contract No. DE-AC02-05CH11231. PNNL is operated by Battelle for the U.S. Department of Energy under Contract DE-AC05-76RL01830. G.B. acknowledges the EPSRC for partial support in PhD funding. The authors acknowledge the use of the IRIDIS High Performance Computing Facility (IRIDIS 5) and associated support services at the University of Southampton in the completion of this work. We are grateful to the UK Materials and Molecular Modeling Hub (Thomas HPC) for computational resources, which is partially funded by EPSRC (EP/P020194/1).

continuum solvation, metal/water interfaces, ab in↗

Interlayer fusion bonding of semi-crystalline polymer composites in extrusion deposition additive manufacturing

This work focuses on the evolution of interlayer fracture toughness properties of fiber-reinforced, semi-crystalline polymers in the extrusion deposition additive manufacturing (EDAM) process. Further, this work bridges the gap between the additive process conditions (time-temperature history) and the effective layer-to-layer fracture properties developed within a printed component. This is the first step to predict delamination that can occur during printing, during cooling to room temperature after printing, and during service performance of an additively manufactured geometry. A phenomenological model is developed for fusion bonding of semi-crystalline polymer matrix composites by coupling the interdiffusion of polymer chains with the evolution of polymer crystallinity. While the interdiffusion is captured by reptation theory of polymer dynamics, the evolution of crystallinity is modeled by phenomenological crystallization kinetics and crystal melting dynamics. Further, a methodology is developed to determine the critical strain energy release rate, G IC of the interlayer interface and experiments are conducted utilizing the double cantilever beam fracture test geometry. In conclusion, predictions of GIC as a function of thermal history are compared with experiments.

36 MATERIALS SCIENCE↗

Influence of shrub willow buffers strategically integrated in an Illinois corn-soybean field on soil health and microbial community composition

We report soil serves many important ecological functions and is an integral part of our existence as a society. However, concerns for soil health are growing globally, in part due to the negative impacts of agricultural management on soil resources. The production of perennial bioenergy crops on marginal land in row-crop production systems is one solution that could improve land-use efficiency and address the sustainability of cropland management. Because the relationship between crop management and the environment is complex, more research is needed to evaluate the potential benefits perennial bioenergy crop production has on soil health, as well as other ecosystem services. In this study, shrub willow buffers were strategically integrated into a corn-soybean cropping system with the main objective of reducing nitrate-N leaching from grain crop production while producing biomass for bioenergy. Two buffer systems (defined by landscape positions) were included for comparison, one on marginal land with exposure to nitrate-N leaching from upslope grain (southern plots) and one on fertile soils with less nitrate-N leaching potential (northern plots). Evaluation of soil (chemistry, bulk density, microbial community) and shrub willow vegetation properties ( fine roots, leaf litter decomposition, and nutrient uptake dynamics), showed that landscape position plays an important role in (1) the dynamics of soil chemical properties, (2) shrub willow's influence and productivity, and (3) the provision of additional ecosystem services such as reductions in nitrous oxide emissions and nitrate-N leaching. In addition, the combination of crop type and landscape position (N-grain, N-willow, S-grain, and S-willow) influenced the species composition of the soil microbial community, resulting in unique and identifiable communities. These results highlight the potential application of shrub willow buffers for ecosystem service provision and support of ecosystem processes; however, understanding the relationship between the microbial community, crop type, and landscape is important for understanding the sustainability of the design.

54 ENVIRONMENTAL SCIENCES↗

Kernel Enriched Meshfree Multiphysics Degradation Modeling of Energy Storage Materials

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and ultimately diminishing performance and service life. With microstructural images supplied by the National Laboratory of the Rockies (NLR), pixel-based meshfree model construction by the reproducing kernel particle method (RKPM) is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. The first kernel enrichment discussed will be the interface modified reproducing kernel (IM-RK) [1, 2], constructed by scaling a smooth kernel function with an interface-distance function to achieve strategic discontinuity types (i.e. weak discontinuities for strain discontinuities and strong discontinuities for cracks) and alleviate Gibbs oscillations near these transition zones. The IM-RK is especially useful for areas in which a known discontinuity-type is expected a priori. The second kernel enrichment to be discussed is a neural network-enhanced reproducing kernel (NN-RK) [3, 4], which is introduced to effectively model non-obvious damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RK is additionally used to inform how crack opening and closure in turn affect the electro-chemo-mechanical responses in the material microstructure. References: [1] Wang, Y., Baek, J., Tang, Y. et al. "Support vector machine guided reproducing kernel particle method for image-based modeling of microstructures," Comput Mech 73, 907-942 (2024). https://doi.org/10.1007/s00466-023-02394-9. [2] Susuki, K., Allen, J. & Chen, J. S.. "Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method," Engineering with Computers (2024). https://doi.org/10.1007/s00366-024-02016-9. [3] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, 4422-4454 (2022). https://doi.org/10.1002/nme.7040.

97 MATHEMATICS AND COMPUTING↗