Search NASA⌕ Search

SEARCH · Search NASA

Results for “loading simulations”

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 541 records · Page 30

Omega-Limit Sets and Input-To-State Stability in Power Grids with Switching Equilibria

This paper studies a power transmission system with both conventional generators (CGs) and distributed energy assets (DEAs) providing frequency control. We consider an operating condition with demand aggregating two dynamic components: one that switches between different values on a finite set, and one that varies smoothly over time. Such dynamic operating conditions may result from protection scheme activations, external cyber-attacks, or due to the integration of dynamic loads, such as data centers. Mathematically, the dynamics of the resulting system are captured by a system that switches between a finite number of vector fields - or modes -, with each mode having a distinct equilibrium point induced by the demand aggregation. To analyze the stability properties of the resulting switching system, we leverage tools from hybrid dynamic inclusions and the concept of ..omega.. -limit sets from sets. Specifically, we characterize a compact set that is semi-globally practically asymptotically stable under the assumption that the switching frequency and load variation rate are sufficiently slow. For arbitrarily fast variations of the load, we use a level-set argument with multiple Lyapunov functions to establish input-to-state stability of a larger set and with respect to the rate of change of the loads. The theoretical results are illustrated via numerical simulations on the IEEE 39-bus test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗

A Data-Driven Pivot-Point-Based Time-Series Feeder Load Disaggregation Method

The load profile at a feeder-head is usually known to utility engineers while the nodal load profiles are not. However, the nodal load profiles are increasingly important for conducting time-series analysis in distribution systems. Therefore, in this paper, we present a pivot-point based, two-stage feeder load disaggregation algorithm using smart meter data. The two stages are load profile selection (LPS) and load profile allocation (LPA). In the LPS stage, a random load profile selection process is first executed to meet the load diversity requirement. Then, a few pairs of pivot points are selected as the matching targets. After that, a matching algorithm will run repetitively to select one load profile at a time for matching the reference load profile at the pivot points. In the LPA stage, the LPS selected load profiles are allocated to each load node on the feeder considering distribution transformer loading limits, load composition, and square-footage. The proposed method is validated using actual data collected in a North Carolina service area. Finally, simulation results show that the proposed method can generate a unique load shape for each load node while match the shape of their aggregated profile with the actual feeder head load profile.

42 ENGINEERING↗

Aeroelastic Modeling and Full-Scale Loads Measurements for Investigation of Single-Axis PV Tracker Wind-Driven Dynamic Instabilities

While wind tunnel testing and proprietary industry modeling tools have been used for years to design and develop PV tracker systems, wind induced dynamic failures are becoming more prevalent and high visibility. Designers and manufactures have reacted by developing new systems or add-on products to improve system dynamics adding to overall system costs. To better understand the physics and aeroelastic behavior that leads to dynamic instabilities NREL researches have taken up a first of a kind study to both develop open source aeroelastic modeling tools and measure tracker loads on a single axis full-scale tracker system at a high wind site. These modeling tools can simulate the fluid-structure interaction driving torsional instabilities under a wide range of turbulent inflow conditions and stow angles. The results of these simulations reveal the flow features associated with panel rotation and the induced loads in the structure. The experimental load measurements can be used to validate models and provide a quantifiable understanding of tracker dynamics, critical loads paths, identify instability markers, inform resilient design, and suggest the most favorable stow approach.

14 SOLAR ENERGY↗

A physics-informed data-driven low order model for the wind velocity deficit at the wake of isolated buildings

Detailed simulations of the Atmospheric boundary layer (ABL) provides information on wind velocity, direction and turbulence for optimizing the design of wind turbines and wind farms and their operations. Distributed wind turbines have been considered for use in urban or suburban areas for covering energy needs using a sustainable energy solution. Simulations of high Reynolds number turbulent ABL flows past obstacles can be used for the prediction of the velocity deficit and turbulence characteristics in the leeward side of building-like obstacles, allowing for the estimation of power losses and fatigue loads in wind turbines. However, fully three-dimensional simulations of high Reynolds number ABL flows using Large Eddy Simulation (LES) or Reynolds Average Navier-Stokes Simulations (RANS) continue to require large computational/ execution times. As a result, these 3D simulations may not be ideal tools for the rapid assessment of wind energy potential. In the present study, 3D numerical results from wall resolved RANS simulations were used for the development of a low order model that allows the rapid assessment of wind velocity and turbulence in the leeward side of an isolated building.

17 WIND ENERGY↗

EnergyPlus Model Context Protocol Server (EnergyPlus-MCP) v0.1

EnergyPlus-MCP is the first open-source Model Context Protocol server specifically designed for EnergyPlus building energy simulation. This innovative software enables AI assistants and other applications to interact programmatically with EnergyPlus through a standardized, secure interface, eliminating traditional technical barriers in building energy modeling. The software provides specialized tools across five functional domains: server management, model configuration and loading, comprehensive building component inspection, systematic model modification, and simulation execution with results visualization. Key features include automated HVAC system discovery and topology mapping, advanced schedule analysis, intelligent model validation, and interactive visualization capabilities. EnergyPlus-MCP's layered architecture ensures robust separation between protocol communication and domain expertise, enabling scalable deployment across organizations, educational institutions, and research teams. Unlike direct LLM approaches that suffer from inconsistent results and security gaps, EnergyPlus-MCP provides validated, reliable interactions while maintaining scientific rigor. This democratizes sophisticated building energy analysis, making EnergyPlus accessible to broader audiences through conversational interfaces and streamlined workflows.

Li, Han [Lawrence Berkeley National Laboratory (LB↗

Exploring the parameter space of MagLIF implosions using similarity scaling. III. Rise-time scaling

Magnetized liner inertial fusion (MagLIF) is a z-pinch magneto-inertial-fusion concept studied at the Z Pulsed Power Facility of Sandia National Laboratories. Two important metrics characterizing current delivery to a z-pinch load are the peak current and the current-rise time, which is roughly the time interval to reach the peak current. It is known that, when driving a z-pinch load with a longer current-rise time, the performance of the z-pinch decreases. However, a theory to understand and quantify this effect is still lacking. Here, in this paper, we utilize a framework based on similarity scaling to analytically investigate the variations in the performance of MagLIF loads when varying the current-rise time, or equivalently, the implosion timescale. To maintain similarity between the implosions, we provide scaling prescriptions of experimental input parameters defining a MagLIF load and derive the expected scaling laws for stagnation conditions and for various performance metrics. We compare predictions of the theory to 2D numerical simulations using the radiation, magneto-hydrodynamic code hydra. For several metrics, we find acceptable agreement between the theory and simulations. Our results show that the voltage φload near the MagLIF load follows a weak scaling law φload ∝$t_{φ}^{-0.12}$ with respect to the characteristic timescale t φ of the voltage source, instead of the ideal φload ∝$t_{φ}^{-1}$ scaling. This occurs because the imploding height of the MagLIF load must increase to preserve end losses. As a consequence of the longer imploding liners, the required total laser preheat energy and delivered electric energy increase. Overall, this study helps understand the trade-offs of the MagLIF design space when considering future pulsed-power generators with shorter and longer current-rise times.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reverse scaling of a bonded-sphere DEM model: Formulation and application to lignocellulosic biomass microstructures

We explore scaling laws for adapting a bonded-sphere discrete element method (BS-DEM) model developed for woody structural mechanics at the millimeter scale to model the mechanics of realistic lignocellulosic biomass microstructures. Two scaling approaches, i.e., the reverse coarse graining (RCG) and equivalent bulk behavior (EBB), are proposed based on the classical mechanics principles and assessed in the single-particle compression and rectangular cuboid block tension tests. The EBB approach is recommended for BS-DEM models with general 3D geometries and is applied to simulate the microindentation test on a realistic microscale pinewood specimen. Simulations are performed to elucidate the impact of specimen thickness and loading position on the specimen’s force–deformation behavior. Furthermore, the range of Young’s modulus obtained from the calibrated BS-DEM simulations can match the experimental measurements. This is the first-of-its-kind study that has explored scale-bridging modeling approaches in biomass micromechanics and has proposed solutions based on BS-DEM models in the microscale.

59 BASIC BIOLOGICAL SCIENCES↗

Grid Impact Analysis of Heavy-Duty Electric Vehicle Charging Stations

This paper presents a grid impact analysis of heavy-duty electric vehicle (EV) charging stations. In this work, heavy-duty EVs have battery capacities high enough to provide a range of 250-500 miles on a single charge, such as long-haul trucks. Heavy-duty EVs will require extremely fast charging rates to reduce charging time and will induce very high charging loads (at the multiple-megawatt scale) if several vehicles charge at the same time. Therefore, analysis is needed to understand the impact of charging station loads on the electric power grid and set the baseline for developing mitigation plans and necessary system upgrades. We develop a systematic procedure to analyze the potential impact of the placement of charging stations on the grid. Charging load is modeled using a DC fast-charging station model. A voltage load sensitivity matrix approach is leveraged to investigate the challenges of placing charging stations on the feeder. Given the charging load profiles and suggested charging station locations, time-series simulations are performed on various connection points on the feeder to understand the impact. The analysis is performed on both the IEEE 34-bus system and a realistic feeder from California. Initial mitigation solutions are developed based on insights from this analysis.

47 OTHER INSTRUMENTATION↗

Sensitivity of Dislocation-GB interactions to simulation setups in atomistic models

Dislocation-grain boundary (GB) interactions play an integral role in strengthening of crystalline materials, and are dependent on external loading conditions and atomic arrangements at the grain boundary. While molecular dynamics (MD) simulations can provide critical insights into relationships between these parameters and the exact interaction, the outcomes are sensitive to the computational setup. Here, in this work, we explore the effect of computational setup (system size, GB orientation and boundary conditions) on stress-controlled dislocation-grain boundary (DGB) reactions using MD simulations for a well-studied copper symmetric $\langle$1 1 2$\rangle$ tilt boundary. The study demonstrates that the DGB reactions (pinning, absorption, transmission, etc.) are sensitive to the system size and the orientation of the GB to the applied stress. Additionally, the current work shows that there is a critical system size for the dislocation-grain boundary reaction stresses to converge, one which was found to be much larger than the setups in comparable MD studies. This was attributed to the specific stress states and the influence of surfaces, fixed or free, as these computational setups are modified. Finally, a stress-controlled setup with minimal model artifacts is examined, and the effect of coupled Schmid and non-Schmid stress components on the dislocation-GB reactions are discussed.

36 MATERIALS SCIENCE↗

Analysis of Grid-Forming Inverter Controls for Grid-Connected and Islanded Microgrid Integration

Autonomous grid-forming (GFM) inverter testbeds with scalable platforms have attracted interest recently. In this study, a self-synchronized universal droop controller (SUDC) was adopted, tested, and scaled in a small network and a test feeder using a real-time simulation tool to operate microgrids without synchronous generators. We presented a novel GFM inverter control adoption to better understand the dynamic behavior of the inverters and their scalability, which can impact the distribution system (DS). This paper provides a steady-state and transient analysis of the GFM power inverter controller via simulation to better understand voltage and frequency stabilization and ensure that the critical electric loads are not affected during a prolonged power outage. The controllers of the GFM inverter are simulated in HYPERSIM to examine voltage and frequency fluctuations. This analysis includes assessing the black start capability for photovoltaic microgrids, both grid-connected and islanded, during transient fault conditions. The high photovoltaic PV penetration levels open exciting opportunities and challenges for the DS. The GFM inverter control demonstrated appropriate response times for synchronization, connection, and disconnection to the grid. The DS has become more resilient and independent of fossil fuels by increasing the penetration of inverter-based distributed energy resources (DERs).

Ward, Laura (ORCID:0000000345501581)↗

Compressible pairwise interaction extended point-particle model for force prediction of shock-particle bed interaction

We propose a pairwise influence framework for the complex unsteady compressible particle-laden flow problem by accounting for the scattered hydrodynamic waves emitting from neighboring particles in a Euler-Lagrange simulation. It has been observed from particle-resolved (PR) simulations of randomly dispersed particle beds under a loading shock that the compressible pseudoturbulence dominates the flow system even after the primary shock has passed, which causes fluctuations observed in the forces experienced by the particle. Moreover, the fact that each particle exists in the vicinity of a random arrangement of other particles modifies the time history of the drag force experienced by each particle during and after the passage of the shock. First, the scattering flow field due to an incoming shock interacting with a single sphere is constructed using an analysis of the flow in the acoustic limit. Then we examine the validity of the compressible Maxey-Riley-Gatignol force model by comparing the force prediction against a PR simulation of two interacting particles for various particle arrangements and incoming shock strength. Subsequently, the neighboring influences are stored as a library of maps that can be used readily in the calculation of the perturbation force. Lastly, the pairwise interaction assumption is evaluated by comparing the force predicted with the model with PR simulations of a randomly packed particle bed of 10% volume fraction for both water and air as the fluid medium for an incoming shock Mach number 1.22. With a considerably lower cost for the implementation of the model compared to PR simulations, it is verified that the model is reasonably accurate in pinpointing particles whose peak force is significantly larger or smaller than the mean drag but also to capture the prolonged fluctuations after the initial shock.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Coupling Microstructural Evolution Simulations to Material Property Degradation Predictions for Plasma-Facing Materials

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500 C to 1500 C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 D/m^2-s. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗

Connect microstructure evolution to property degradation with validated simulation

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500$^o$C to 1500$^o$C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 $\frac{D}{m^2s}$. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗

COHORT: Coordination of Heterogeneous Thermostatically Controlled Loads for Demand Flexibility

Demand flexibility is increasingly important for power grids. Careful coordination of thermostatically controlled loads (TCLs) can modulate energy demand, decrease operating costs, and increase grid resiliency. We propose a novel distributed control framework for the Coordination Of HeterOgeneous Residential Thermostatically controlled loads (COHORT). COHORT is a practical, scalable, and versatile solution that coordinates a population of TCLs to jointly optimize a grid-level objective, while satisfying each TCL’s end-use requirements and operational constraints. To achieve that, we decompose the grid-scale problem into subproblems and coordi- nate their solutions to find the global optimum using the alternating direction method of multipliers (ADMM). The TCLs’ local problems are distributed to and computed in parallel at each TCL, making COHORT highly scalable and privacy-preserving. While each TCL poses combinatorial and non-convex constraints, we characterize these constraints as a convex set through relaxation, thereby making COHORT computationally viable over long planning horizons. After coordination, each TCL is responsible for its own control and tracks the agreed-upon power trajectory with its preferred strategy. In this work, we translate continuous power back to discrete on/off actuation, using pulse width modulation. COHORT is generalizable to a wide range of grid objectives, which we demonstrate through three distinct use cases: generation following, minimizing ramping, and peak load curtailment. In a notable experiment, we validated our approach through a hardware-in-the-loop simulation, including a real-world air conditioner (AC) controlled via a smart thermostat, and simulated instances of ACs modeled after real-world data traces. During the 15-day experimental period, COHORT reduced daily peak loads by an average of 12.5% and maintained comfortable temperatures.

demand response↗

Hybrid power plants: An effective way of decreasing loss-of-load expectation

Diversifying variable renewable resources by combining wind, solar photovoltaic, and battery assets in a hybrid power plant can increase renewable energy usage efficiency and improve system flexibility, particularly in distributed energy systems. However, the resilience impact of these systems, particularly outage mitigation, can be difficult to quantify due to uncertainty in resource, energy demand, and outage occurrence. Here, this study outlines a framework to quantify the incremental benefit of hybrid power plant assets for reducing loss-of-load expectation during random outage events. Hybrid power plant performance during outages (considering varying duration and severity) is simulated using a Monte Carlo methodology to reflect uncertainty associated with renewable resource, load demand, and outage timing. Results demonstrate the additional incremental value from increasingly hybrid designs, in which relative capacities of wind, solar photovoltaic, and storage assets contribute to lower loss-of-load expectation than the constituent technologies would alone. The value of added wind or solar capacity increases as the plant composition approaches an equal split. The value of added battery capacity depends on the outage duration and severity, but the first 50 MWh of added storage capacity is the most valuable for reducing the loss-of-load expectation for all plant designs.

14 SOLAR ENERGY↗

A Scalable and Distributed Algorithm for Managing Residential Demand Response Programs using Alternating Direction Method of Multipliers (ADMM)

For effective engagement of residential demand-side resources and to ensure efficient operation of distribution networks, we must overcome the challenges of controlling and coordinating residential components and devices at scale. To overcome this challenge, we present a distributed and scalable algorithm with a three-level hierarchical information exchange architecture for managing the residential demand response programs. First, a centralized optimization model is formulated to maximize community social welfare. Then, this centralized model is solved in a distributed manner with alternating direction method of multipliers (ADMM) by decomposing the original problem to utility-level and house-level problems. The information exchange between the different layers is limited to the primary residual (i.e., supply-demand mismatch), Lagrangian multipliers, and the total load of each house to protect each customer’s privacy. Simulation studies are performed on the IEEE 33 bus test system with 605 residential customers. The results demonstrate that the proposed approach can save customers’ electricity bills and reduce the peak load at the utility level without much affecting customers’ comfort and privacy. Finally, a quantitative comparison of the distributed and centralized algorithms shows the scalability advantage of the proposed ADMM-based approach, and it gives benchmarking results with achievable value for future research works.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Image-Based Digital Twin for Assessing the Coupled Electro-Chemo-Mechanical Behavior of Li-Ion Batteries

In energy storage materials, strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking during charge/discharge cycling, resulting in reduced performance and service life. In this work, a digital twin is created to investigate the performance of a Li-ion battery cathode and simulate degradation accumulation. Pixel-based model construction is used to represent the complex material geometries from microstructural images supplied by the National Renewable Energy Laboratory (NREL). Because of the expected large deformation and crack opening, the reproducing kernel particle method (RKPM), a meshfree method with discretization at the image pixels, is used to approximate the field variables: electrostatic potential, concentration, and displacement. An interface modified reproducing kernel (IM-RK) is 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 discrete nature of the images' pixel points is employed throughout the model and approximation construction. The mechanical model is verified using an image-based microstructure under tensile loading. A transient electro-chemo-mechanical coupled simulation is performed to evaluate the potential micro-cracking induced degradation of the battery cathode material subjected to charge/discharge cycling.

battery degradation↗