Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cohesive coupling”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Strain engineering of doped hydrogen passivated silicon quantum dots

Silicon quantum dots are nanomaterials that are attractive candidates for photovoltaic applications. Doping of these materials creates p-n junctions and is important for solar cells. In this work, we present a first-principles study of the coupled influence of doping and strain on the stability, energy gap, Fermi level, electronic density, and density of states of hydrogen-passivated silicon quantum dots. We find that the cohesive energy and the energy gap decrease with increasing quantum dot size and are strongly influenced by strain. Furthermore, the response to strain also depends on the size of the quantum dot and dopant type. We present expressions of cohesive energy and energy gap as power-law of size and polynomial dependence on strain. We also show that the Fermi energy increases with size for pristine and p-type doping but decreases with size for n-type doping. We also discuss the influence of strain and dopant type on the density of states and electron density of the quantum dots.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

A new approach to model geomaterials with heterogeneous properties in thermo-hydro-mechanical coupled problems

The main objective of this article is to present a new approach to model coupled thermo-hydro-mechanical problems considering geomaterials with heterogeneous properties. This approach has been implemented in the software CODE_BRIGHT and it provides the possibility of considering geomaterials with a spatially correlated heterogeneous field of porosity, following a normal distribution. This spatial correlation can be isotropic or anisotropic. An important feature of this approach is that material properties such as intrinsic permeability, thermal conductivity, diffusivity, retention curve, elastic modulus or cohesion are defined as a function of porosity and, thus, they become heterogeneous with spatial correlation and, eventually, anisotropic. A validation exercise and other basic numerical examples have been carried out to illustrate the possibilities of the proposed approach. The results, which have been compared with a homogeneous case, show that considering heterogeneous fields can be relevant in different modelling problems, especially coupled thermo-hydro-mechanical problems.

42 ENGINEERING↗
DOE OSTI2021

Attaining regularization length insensitivity in phase-field models of ductile failure

In this work, a cohesive phase-field model of ductile fracture in a finite-deformation setting is presented. The model is based on a free-energy function in which both elastic and plastic work contributions are coupled to damage. Using a strictly variational framework, the field evolution equations, damage kinetics, and flow rule are jointly derived from a scalar least-action principle. Particular emphasis is placed on the use of a rational function for the stress degradation that maintains a fixed effective strength with decreasing regularization length. The model is employed to examine crack growth in pure mode-I problems through the generation of crack growth resistance (J-R) curves. In contrast to alternative models, the current formulation gives rise to J-R curves that are insensitive to the regularization length. Numerical evidence suggests convergence of local fields with respect to diminishing regularization length as well.

42 ENGINEERING↗
DOE OSTI2023

Pressure Transient Analyses and Poroelastic Modeling of Hydraulic Fracture Dilation for Multiple Injections at the Devine Fracture Pilot Site

Our team has conducted electromagnetic (EM) surveys for the past six years to monitor hydraulic-fracture behavior at the Devine Fracture Pilot Site (DFPS). The sub-horizontal orientation of a shallow hydraulic fracture at the DFPS provides uniform access to the fracture area for interrogation and data collection. Ahmadian et al. (2023) suggested a possible correlation between spatiotemporal changes in the flow rate, bottomhole pressure (BHP), and the observed surface recorded electric field at the DFPS. In this paper, we present the development of poroelastic forward models and pressure transient analyses (PTAs) to support the development of a multiphysics inverse model for these EM surveys. First, we conducted PTAs of the shut-in periods after six injections out of 10 to determine the fracture closure pressure (FCP) or the overburden pressure used in a poroelastic fracture reopening model. Second, we developed a finite-element poroelastic model throughout five injection cycles to include the effect of the cumulative injected volumes due to the previous injections on current fracture dilation in the presence of highly permeable unpropped and propped zones adjacent to the cohesive layer that models fracture reopening. Fracture reopening in this poroelastic model is based on a calibrated traction-separation response using the bottomhole pressure collected in two injection campaigns in 2020 and 2022. We used the outcomes of a previous simulation study of the primary hydraulic-fracturing stimulation to define the dimension of an unpropped fracture zone ahead of the propped fracture area. The PTAs led to FCPs consistent with those obtained using the injection data collected at the DFPS in 2020. Further, these analyses showed that at later injections, the fracture closure occurred at a later time with respect to the shut-in time, inferring the effect of cumulative injected volumes in previous injections. The simulation results show that considering the propped and unpropped fracture zones improves our poroelastic model in predicting the injection-well BHP. The numerical simulation results demonstrate a significant excess pore pressure near the fracture because of the preceding formation loadings by the previous injections. The obtained fracture dilation area and fluid pressure distribution provide a basis to improve the development of a multiphysics inverse model. Furthermore, in an iteratively coupled scheme, this pressure distribution can be introduced into EM models to render a holistic view of the causative mechanisms for the surface signal anomalies.

02 PETROLEUM↗
DOE OSTI2022

ORNL Neutron Cross Section Measurements of 90 Zr NCSP ND-1 Task [Slides]

Luiz Leal's idea to integrate the experimental uncertainties, covariances, and sensitivities in this workflow and centralize the process at one location led to the creation of a cohesive tool for data testing. It has been practiced for more than two decades by the NCS and ND groups at ORNL. The development and use of the ORNL codes TSUNAMI and TSURFER within SCALE made it possible. Today, it has been extended to SAMINT which enables the coupling of differential and integral data evaluation in a continuous-energy framework.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

ORNL Neutron Cross Section Measurements of 90 Zr NCSP ND-1 Task [Slides]

Luiz Leal's idea to integrate the experimental uncertainties, covariances, and sensitivities in the ND workflow and to centralize the process at one location created a cohesive tool for data testing. It has been practiced for more than two decades by the NCS and ND groups at ORNL. The development and use of the ORNL codes TSUNAMI and TSURFER within SCALE made it possible. Today, it has been extended to SAMINT, which enables the coupling of differential and integral data evaluation in a continuous-energy framework.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Electronic structure and magnetism of pristine, defected, and strained $\mathrm{Ti_2N}$ $\mathrm{MX}$ene

Here, from first principles electronic structure calculations, we unravel the evolution of structural, electronic, and magnetic properties of pristine, defected, and strained titanium nitride MXene with different functional groups (-F, -O, -H, and -OH). The formation and cohesive energies reveal their chemical stability. The dynamical stability of Ti 2 N mono-layer is also confirmed by phonon calculations. The MAX phase and defect free functionalized MXenes are metallic except for oxygen terminated (Ti 2 NO 2 ) one which is 100% spin polarized half-metallic ferromagnet. The spin–orbit coupling significantly influences the bare MXene (Ti 2 N) to exhibit Dirac topology and band inversion near the high symmetry directions. The strain effect sways the Fermi level thereby shifting it towards lower energy state under compression and towards higher energy state under tensile strain in Ti 2 NH 2 . The Ti 2 NO 2 exhibits exotic electronic structure not only in pristine but also in strained and defected structures. Its half-metallic nature changes to semi-metallic under 1% compression and it is completely destroyed under 2% compression. In single vacancy defect, its band structure remarkably transforms from half-metallic to semi-conducting with large band gap in 12.5% Ti, weakly semi-conducting in 5.5% Ti, and semi-metallic in 12.5% O. The 25% N defect changes its half-metallic characteristic to metallic. Further, the 12.5% Co substitution preserves its half-metallic character, whereas Mn substitution allows it to convert half-metallic characteristic into weak semi-metallic characteristic preserving ferromagnetism. However, Cr substitution converts half-metallic ferromagnetic state to half-metallic anti-ferromagnetic state. The understanding made here on collective structural stability, and electronic band structure, and magnetic phenomena in novel 2D Ti 2 N derived MXenes open up their possibility in designing them for synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Powering ≈50 µm Motion by a Molecular Event in DNA Crystals

A major challenge in material design is to couple nanoscale molecular and supramolecular events into desired chemical, physical, and mechanical properties at the macroscopic scale. Here, a novel self-assembled DNA crystal actuator is reported, which has reversible, directional expansion and contraction for over 50 μm in response to versatile stimuli, including temperature, ionic strength, pH, and redox potential. The macroscopic actuation is powered by cooperative dissociation or cohesion of thousands of DNA sticky ends at the designed crystal contacts. The increase in crystal porosity and cavity in the expanded state dramatically enhances the crystal capability to accommodate/encapsulate nanoparticles/proteins, while the contraction enables a “sponge squeezing” motion for releasing nanoparticles. This crystal actuator is envisioned to be useful for a wide range of applications, including powering self-propelled robotics, sensing subtle environmental changes, constructing functional hybrid materials, and working in drug controlled-release systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling Sequences of Earthquakes and Aseismic Slip (SEAS) in Elasto-Plastic Fault Zones With a Hybrid Finite Element Spectral Boundary Integral Scheme

We present a coupled finite element spectral boundary integral framework for modeling sequences of earthquakes and aseismic slip on a 2-D planar rate-and-state fault with off-fault visco-plastic response in the plane strain approximation. The model resolves both slow aseismic deformation and inertia effects during rapid slip. As an application, we perform two sets of simulations with different choices of cohesion to explore the co-evolution of fault slip, bulk plasticity and local stress fields. The first set implements a relatively large value of the cohesion parameter, which results in limiting inelastic strain accumulation to dynamic rupture phases. The second set implements a smaller cohesion, allowing for plastic strain to accumulate in both seismic and aseismic phases. For the first model, our results indicate that the extent and distribution of plastic strain depend on the angle of maximum compressive principal stress. At larger angles, inelastic strain accumulates on the extensional side of a dynamically propagating rupture. At smaller angles, the extent of plasticity is limited to the compressional side of the domain. At smaller cohesion values, off-fault plasticity may occur during aseismic slip, which alters the nucleation characteristics and earthquake sequence pattern. Furthermore, our results at lower cohesion values indicate that plastic strain accumulation may occur in both the extensional and compressional sides of the off-fault bulk even at higher angles of maximum compression. In conclusion, this produces damage patterns that deviate from the traditional off-fault fan-like distribution observed in dynamic rupture simulations and emphasizes the significance of long-term deformation in interpreting observations.

58 GEOSCIENCES↗
DOE OSTI2020

Advancement of Commercial Intumescent Expanding Foams for Deactivation and Decommissioning in the Nuclear Sector - 20198

Florida International University (FIU), in collaboration with The Department of Energy's Office of Environmental Management (DOE-EM), Savannah River National Laboratory (SRNL), and sites across the Savannah River complex, have identified an operational requirement for a fixative technology that is intended to immobilize and/or isolate residual contamination within a 3-dimensional space. Fixation of radiological contamination can reduce worker risk and mitigate potentially hazardous conditions, however nearly every marketed contamination fixative has been found to be flammable; a significant concern in radiological facilities. Coupled with this, industry fixatives are normally used as a thin coating which can present problems when attempting to stabilize irregular geometry or areas that are difficult to access whilst ensuring full coverage. The technical evaluation and advancement of commercial-off-the-shelf (COTS) polyurethane foams has yielded a down-selected candidate that shows potential in meeting the requirements to support deactivation and decommissioning activities. Several performance criteria have been established and tested to progress the technology readiness level towards an active field demonstration (TRL-7). Such criteria include: mechanical failure limits, adhesive and cohesive properties, thermal/fire resilience, determining thermal behavior, ability to immobilize contamination, and a means of non-destructive evaluation of applications. The test scenario examined was targeted towards an application for decommissioning nuclear pipework, in which the down-selected polyurethane foam would act as a barrier to segregate pipework and mitigate the potential for release during cutting, packaging, and storage operations. Testing carried out at SRNL included: mechanical evaluation of tensile, compressive, and adhesion strength by dynamic mechanical analysis (DMA), as well as thermogravimetric analysis (TGA). FIU examined the foam's fixative properties by utilizing phosphorescent europium-dysprosium doped strontium aluminate powder to investigating the extent to which contamination can be immobilized. FIU has also exploited previous successes in the field of intumescent technologies to assess the down-selected foam's tolerance to an extreme fire scenario, while maintaining the ability to effectively mitigate a contamination release. Parallel to this, extensive thermal investigations were carried out to determine the upper boundary of anticipated heat generation during the curing process as heat generation has the potential to compromise rubber parts of contaminated enclosures. These investigations subsequently yielded a promising method for a non-destructive application evaluation by means of infrared thermography. Utilizing the high sensitivity of modern IR cameras, coupled with the heat generated during the curing process of the polyurethane foam, FIU has been exploring the concept of monitoring the external pipe surface for indications of an irregular or abnormal application, thus informing operational decision making. The testing carried out utilized several current 'best fit' ASTM standards, which serve as helpful guidelines for testing, however, a precise definition of the operational parameters and requirements is still necessary. With continued collaboration with SRNL, FIU aims refine said definitions and develop new standards by which this, and other decommissioning technologies, can be accredited by relevant standards based testing. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

A combined experimental and computational analysis of failure mechanisms in open-hole cross-ply laminates under flexural loading

In this work, integrated experimental tests and computational modeling are proposed to investigate the failure mechanisms of open-hole cross-ply carbon fiber reinforced polymer (CFRP) laminated composites. In particular, we propose two effective methods, which include width-tapered double cantilever beam (WTDCB) and fixed-ratio mixed-mode end load split (FRMMELS) tests, to obtain the experimental data more reliably. We then calibrate the traction-separation laws of cohesive zone model (CZM) used among laminas of the composites by leveraging these two methods. The experimental results of fracture energy, i.e. G Ic and G Tc , obtained from WTDCB and FRMMELS tests are generally insensitive to the crack length thus requiring no effort to accurately measure the crack tip. Moreover, FRMMELS sample contains a fixed mixed-mode ratio of G IIc /G Tc depending on the width taper ratio. Examining comparisons between experimental results of FRMMELS tests and failure surface of B–K failure criterion predicted from a curve fitting, good agreement between the predictions and experimental data has been found, indicating that FRMMELS tests are an effective method to determine mixed-mode fracture criterion. In addition, a coupled experimental-computational modeling of WTDCB, edge notched flexure, and FRMMELS tests are adopted to calibrate and validate the interfacial strengths. Finally, failure mechanisms of open-hole cross-ply CFRP laminates under flexural loading have been studied systematically using experimental and multi-scale computational analyses based on the developed CZM model. The initiation and propagation of delamination, the failure of laminated layers as well as load-displacement curves predicted from computational analyses are in good agreement with what we have observed experimentally.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Large-scale protein level comparison of Deltaproteobacteria reveals cohesive metabolic groups

Abstract Deltaproteobacteria, now proposed to be the phyla Desulfobacterota, Myxococcota, and SAR324, are ubiquitous in marine environments and play essential roles in global carbon, sulfur, and nutrient cycling. Despite their importance, our understanding of these bacteria is biased towards cultured organisms. Here we address this gap by compiling a genomic catalog of 1 792 genomes, including 402 newly reconstructed and characterized metagenome-assembled genomes (MAGs) from coastal and deep-sea sediments. Phylogenomic analyses reveal that many of these novel MAGs are uncultured representatives of Myxococcota and Desulfobacterota that are understudied. To better characterize Deltaproteobacteria diversity, metabolism, and ecology, we clustered ~1 500 genomes based on the presence/absence patterns of their protein families. Protein content analysis coupled with large-scale metabolic reconstructions separates eight genomic clusters of Deltaproteobacteria with unique metabolic profiles. While these eight clusters largely correspond to phylogeny, there are exceptions where more distantly related organisms appear to have similar ecological roles and closely related organisms have distinct protein content. Our analyses have identified previously unrecognized roles in the cycling of methylamines and denitrification among uncultured Deltaproteobacteria. This new view of Deltaproteobacteria diversity expands our understanding of these dominant bacteria and highlights metabolic abilities across diverse taxa.

Langwig, Marguerite V. (ORCID:0000000202472816)↗
DOE OSTI2026

V2XConnect: Harmonizing the Landscape of Bidirectional Charging Codes, Standards, and Communication Protocols

The utility grid is constantly evolving with new generation sources and diverse loads from backup generators, backup batteries, rooftop photovoltaics (PV), and flexible charging loads. Each of these technologies presents consumers with an opportunity to work with their local utility to promote grid stability, reliability, and affordability. Through a wide range of programs, utilities are incentivizing the use of these resources to benefit the grid. However, especially for devices with inverters that supply power back to the grid, it is important that these assets comply with local grid codes to ensure human safety, power quality, and grid reliability. One of the most promising of these new technologies is the bidirectional electric vehicle (EV) coupled with bidirectional electric vehicle supply equipment (EVSE). This report introduces the technologies that enable bidirectional charging and how these systems operate. It then outlines the necessary interconnection codes and industry standards that govern the operations of grid-tied inverters, including those in bidirectional charging. A key challenge is then outlined from the complex landscape of communication protocols and proprietary systems designed to coordinate these grid assets. Finally, a single gateway solution is proposed to harmonize communications between all assets across an entire site for a cohesive response to grid codes and utility programs that is flexible enough to support a range of market options.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗
DOE OSTI2026

Ionic-content-driven restructuring of spirobisindane ionene networks: implications for mechanics, self-healing, and gas transport

Polymers of intrinsic microporosity (PIMs) offer exceptional gas permeability but remain brittle and susceptible to physical aging, limiting their durability in separation applications. Here, we introduce a reconfigurable microporous polymer network that uniquely integrates permanent PIM microporosity with autonomous, intrinsic self-healing driven by imidazolium-based ionic motifs. Spirobisindane units generate the intrinsic free-volume architecture, while an imidazolium-containing polyamide ionene supplies dynamic ionic and hydrogen-bonding interactions that reorganize under mild activation. Incorporation of imidazolium-based ionic liquids further tunes cohesion, mobility, and densification, enabling the network to relax, re-associate, and retain microporosity without structural collapse. Through a comprehensive multiscale approach combining spectroscopy, scattering, thermal and mechanical characterization with all-atom molecular dynamics and density functional theory calculations, we elucidate how ionic content, as a single control parameter that reshapes free-volume distributions, modulates local coordination environments, and governs relaxation and healing kinetics. At intermediate ionic loadings, the networks achieve rapid, repeatable self-healing while maintaining CO$_2$ selectivity, demonstrating an optimal balance between segmental mobility and structural integrity. By establishing how hierarchical ionic interactions couple structure, dynamics, and transport in microporous ionene networks, this work provides generalizable design rules for adaptive soft-matter systems that require simultaneous mechanical resilience, reconfigurability, and selective gas transport.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗
DOE OSTI2023

Meshfree Multiphysics Damage Modeling of Li-ion Battery Materials

Through repeated charging and discharging cycles, the electrodes of a Li-ion battery experience significant swelling and contraction due to the movement of lithium, also known as intercalation. This research focuses specifically on the chemo-mechanical cracking in the cathode. Because cathode particles are comprised of many randomly-oriented grains, which have highly anisotropic material properties, the expansion and contraction is very non-uniform. As a result, stress concentrations tend to form between grains, which necessitates the modeling of crack propagation largely along grain boundaries and material interfaces. Chemo-mechanical damage models are generally simulated by one of two methods: the cohesive zone model (CZM) and the continuous damage model (CDM). The CZM is more accurate at capturing the sharp discontinuities of a crack but is very computationally expensive and intractable for large-scale models as a result. Conversely, the CDM is easily computed but not well-suited to easily allow for discontinuous field variables, which are inherent across a crack. This study aims to improve the CDM's ability to capture discontinuous cracks. Current versions of the CDM use the finite element method (FEM), which is one of the most widely used approaches for spatial discretization. This research investigates the use of the Reproducing Kernel Particle Method (RKPM), a meshfree method, for spatial discretization and aims to achieve a coupled chemo-mechanical crack propagation model that enhances accuracy while maintaining high computational efficiency. This model features a fully coupled, iterative electrochemistry solution, which informs the meshfree damage model of impending crack formation. To gain a better understanding of how the meshfree model compares to the current FEM model in capturing cathode crack propagation, further investigation is needed.

damage modeling↗