Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak transitions”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Temperature-independent giant dielectric response in transitional BaTiO 3 thin films

Ferroelectric materials exhibit the largest dielectric permittivities and piezoelectric responses in nature, making them invaluable in applications from supercapacitors or sensors to actuators or electromechanical transducers. The origin of this behavior is their proximity to phase transitions. However, the largest possible responses are most often not utilized due to the impracticality of using temperature as a control parameter and to operate at phase transitions. This has motivated the design of solid solutions with morphotropic phase boundaries between different polar phases that are tuned by composition and that are weakly dependent on temperature. Thus far, the best piezoelectrics have been achieved in materials with intermediate (bridging or adaptive) phases. But so far, complex chemistry or an intricate microstructure has been required to achieve temperature-independent phase-transition boundaries. Here, we report such a temperature-independent bridging state in thin films of chemically simple BaTiO3. A coexistence among tetragonal, orthorhombic, and their bridging low-symmetry phases are shown to induce continuous vertical polarization rotation, which recreates a smear in-transition state and leads to a giant temperature-independent dielectric response. The current material contains a ferroelectric state that is distinct from those at morphotropic phase boundaries and cannot be considered as ferroelectric crystals. We believe that other materials can be engineered in a similar way to contain a ferroelectric state with gradual change of structure, forming a class of transitional ferroelectrics. Similar mechanisms could be utilized in other materials to design low-power ferroelectrics, piezoelectrics, dielectrics, or shape-memory alloys, as well as efficient electro- and magnetocalorics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dynamical phase transition in a strongly hybridized phonon-triplon chain

In this article, we study a dimerized spin-1/2 chain, such as CuGeO 3 , hosting triplon excitations coupled to optical phonons under weak terahertz laser driving. Both phonons and triplons weakly lose energy into the surrounding baths, forming a non-equilibrium steady state. In the strong phonon-triplon coupling regime, phonons near the two-triplon continuum hybridize strongly with triplons. Using mean-field Lindblad dynamics, we show that strong hybridization induces sharp first-order phase transitions—either single or simultaneous double—in the emission spectrum, mainly due to nonlinearities in the chain and the instantaneous breaking of spatial inversion symmetry, as captured by the phonon displacement. Using mean-field Floquet analysis of harmonic modes in both sectors, we analytically confirm the existence of these phase transitions. Furthermore, we map the complete steady-state phase diagram by varying key control parameters and provide experimentally relevant parameters for observing these transitions in laser-driven CuGeO 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Optical absorption of interlayer excitons in transition-metal dichalcogenide heterostructures

Interlayer excitons, electron-hole pairs bound across two monolayer van der Waals semiconductors, offer promising electrical tunability and localizability. Because such excitons display weak electron-hole overlap, most studies have examined only the lowest-energy excitons through photoluminescence. Here, we directly measured the dielectric response of interlayer excitons, which we accessed using their static electric dipole moment. We thereby determined an intrinsic radiative lifetime of 0.40 nanoseconds for the lowest direct-gap interlayer exciton in a tungsten diselenide/molybdenum diselenide heterostructure. We found that differences in electric field and twist angle induced trends in exciton transition strengths and energies, which could be related to wave function overlap, moiré confinement, and atomic reconstruction. Through comparison with photoluminescence spectra, this study identifies a momentum-indirect emission mechanism. Characterization of the absorption is key for applications relying on light-matter interactions.

42 ENGINEERING↗
DOE OSTI2021

Phase transition mechanism and bandgap engineering of Sb2S3 at gigapascal pressures

Abstract Earth-abundant antimony trisulfide (Sb 2 S 3 ), or simply antimonite, is a promising material for capturing natural energies like solar power and heat flux. The layered structure, held up by weak van-der Waals forces, induces anisotropic behaviors in carrier transportation and thermal expansion. Here, we used stress as mechanical stimuli to destabilize the layered structure and observed the structural phase transition to a three-dimensional (3D) structure. We combined in situ x-ray diffraction (XRD), Raman spectroscopy, ultraviolet-visible spectroscopy, and first-principles calculations to study the evolution of structure and bandgap width up to 20.1 GPa. The optical band gap energy of Sb 2 S 3 followed a two-step hierarchical sequence at approximately 4 and 11 GPa. We also revealed that the first step of change is mainly caused by the redistribution of band states near the conduction band maximum. The second transition is controlled by an isostructural phase transition, with collapsed layers and the formation of a higher coordinated bulky structure. The band gap reduced from 1.73 eV at ambient to 0.68 eV at 15 GPa, making it a promising thermoelectric material under high pressure.

Cui, Zhongxun (ORCID:0000000275044177)↗
DOE OSTI2023

Magnetism and quantum melting in moiré-material Wigner crystals

Recent experiments have established that semiconductor-based moir´e materials can host incompressible states at a series of fractional moir´e-miniband fillings. These states have been identified as generalized Wigner crystals in which electrons localize on a subset of the available triangular-lattice moir´e superlattice sites. In this article, we use momentum-space exact diagonalization to investigate the many-body ground state evolution at rational fillings from the weak-hopping classical lattice gas limit, in which only spin degrees-of-freedom are active at low energies, to the strong-hopping metallic regime where the Wigner crystals melt. We specifically address the nature of the magnetic ground states of the generalized Wigner crystals at fillings ν = 1/3 and ν = 2/3.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Vortex guidance and transport in channeled pinning arrays

Here, we numerically examine vortices interacting with pinning arrays where a portion of the pinning sites have been removed in order to create coexisting regions of strong and weak pinning. The region without pinning sites acts as an easy-flow channel. For driving in different directions with respect to the channel, we observe distinct types of vortex flow. When the drive is parallel to the channel, the flow first occurs in the pin free region followed by a secondary depinning transition in the pinned region. At high vortex densities there is also an intermediate plastic flow phase due to the coupling between the weak and strong pinning regions. For driving applied perpendicular to the channel, we observe a jammed phase in which vortices accumulate on the boundary of the pinned region due to the vortex-vortex repulsion, while at higher drives the vortices begin to flow through the pinning array. For driving at an angle to the channel, depending on the filling we observe a drive-induced reentrant pinning effect as well as negative differential mobility which occurs when vortices move from the unpinned to the pinned portion of the sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Effects of anisotropy and disorder on the superconducting properties of niobium

We report results for the superconducting transition temperature and anisotropic energy gap for pure niobium based on Eliashberg’s equations and electron and phonon band structures computed from density functional theory. The electronic band structure is used to construct the Fermi surface and calculate the Fermi velocity at each point on the Fermi surface. The phonon bands are in excellent agreement with inelastic neutron scattering data. The corresponding phonon density of states and electron–phonon coupling define the electron–phonon spectral function, α 2 F(p, p'; ω), and the corresponding electron–phonon pairing interaction, which is the basis for computing the superconducting properties. The electron–phonon spectral function is in good agreement with existing tunneling spectroscopy data except for the spectral weight of the longitudinal phonon peak at ℏω LO = 23 meV. We obtain an electron–phonon coupling constant of λ = 1.057, renormalized Coulomb interaction μ* = 0.218, and transition temperature T c = 9.33 K. The corresponding strong-coupling gap at T = 0 is modestly enhanced, Δ 0 = 1.55 meV, compared to the weak-coupling BCS value $Δ^{\text {wc}}_{0}$ = 1.78 k B T c = 1.43 meV. The superconducting gap function exhibits substantial anisotropy on the Fermi surface. We analyze the distribution of gap anisotropy and compute the suppression of the superconducting transition temperature using a self-consistent T-matrix theory for quasiparticle-impurity scattering to describe niobium doped with non-magnetic impurities. We compare these results with experimental results on niobium SRF cavities doped with nitrogen impurities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Resonant Soft X-Ray Scattering from Stripe-Ordered La 2 - x Ba x Cu O 4 Detected by a Transition-Edge Sensor Array Detector

Resonant soft x-ray scattering (RSXS) is a leading probe of valence-band order in materials best known for detecting charge-density-wave order in copper-oxide superconductors. One of the biggest limitations on the RSXS technique is the presence of a severe fluorescence background, which, like the RSXS cross section itself, is enhanced under resonant conditions. This background prevents the study of weak signals such as diffuse scattering from glassy or fluctuating order that is spread widely over momentum space. Recent advances in superconducting transition-edge-sensor (TES) detectors have led to major improvements in energy resolution and detection efficiency in the soft x-ray range. Here, we perform a RSXS study of stripe-ordered La 2 - x Ba x Cu O 4 at the Cu L 3 / 2 edge (932.2 eV) using a TES detector with 1.5-eV resolution, to evaluate its utility for mitigating the fluorescence-background problem. We find that, for suitable degree of detuning from the resonance, the TES rejects the fluorescence background, leading to a five to ten times improvement in the statistical quality of the data compared to an equivalent, energy-integrated measurement. We conclude that a TES introduces a promising method to reducing background in RSXS studies and may lead to discoveries in materials exhibiting valence-band order that is fluctuating or glassy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Photon-mediated entanglement scheme between a ZnO semiconductor defect and a trapped Yb ion

We propose an optical scheme to generate an entangled state between a trapped ion and a solid state donor qubit through which-path erasure of identical photons emitted from the two systems. The proposed scheme leverages the similar transition frequencies between In donor bound excitons in ZnO and the 2 P 1/2 to 2 S 1/2 transition in Yb+. The lifetime of the relevant ionic state is longer than that of the ZnO system by a factor of 6, leading to a mismatch in the temporal profiles of emitted photons. A detuned cavity-assisted Raman scheme weakly excites the donor with a shaped laser pulse to generate photons with a 0.99 temporal overlap to the Yb+ emission and partially shift the emission of the defect toward the Yb+ transition. Here, the remaining photon shift is accomplished via the dc Stark effect. We show that an entanglement rate of 2.1 × 10 4 s –1 and an entanglement fidelity of 94% can be attained using a weak excitation scheme with reasonable parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Turbulence regimes in the nocturnal roughness sublayer: Interaction with deep convection and tree mortality in the Amazon

We investigated the influence of seasonality and proximity to the forest canopy on nocturnal turbulence regimes in the roughness sublayer of a Central Amazon forest. Since convective systems of different scales are common in this region, we also analyzed the effect of extreme wind gusts (propagated from convective downdrafts) on the organization of the turbulence regimes, and their potential to cause the mortality of canopy trees. Our data include high-frequency winds, temperature and ozone concentration at different heights during the dry and wet seasons of 2014. In addition, we used critical wind-speed data derived from a tree-winching experiment and a modeling study conducted in the same study site. Two different turbulence regimes were identified at three heights above the canopy: a weakly stable (WS) and a very stable regime (VS). The threshold wind speeds that mark the transition between turbulence regimes were larger during the dry season and increased as a function of the height above the canopy. The turbulent fluxes of sensible heat and momentum during the WS accounted for 88% of the entire nighttime flux. Downdrafts occurred only in the WS and favored a fully coupled state of wind flow along the canopy profile. The destructive potential of winds was four times higher than on nights without downdrafts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Cooper minima of atoms in a plasma environment

Recent efforts have extended the calculations of photon ionization and associated Cooper minima to atoms in a plasma environment simulated by screening the atomic potential. Furthermore, in this paper previously unrecognized features of bound electron photon excitations in screened atomic potentials are examined. First, the Cooper minima location as a function of photon energy is shown to be weakly dependent on the screening, which is explained by stationary perturbation theory. Second, the screening transfers transition probability from high to low photon energies relative to the isolated atom result.

Cooper minima↗
DOE OSTI2023

The gallium anomaly

In order to test the end-to-end operations of gallium solar neutrino experiments, intense electron-capture sources were fabricated to measure the responses of the radiochemical SAGE and GALLEX/GNO detectors to known fluxes of low-energy neutrinos. Such tests were viewed at the time as a cross-check, given the many tests of 71 Ge recovery and counting that had been routinely performed, with excellent results. However, the four 51 Cr and 37 Ar source experiments yielded rates below expectations, a result commonly known as the Ga anomaly. As the intensity of the electron-capture sources can be measured to high precision, the neutrino lines they produce are fixed by known atomic and nuclear rates, and the neutrino absorption cross section on 71 Ga is tightly constrained by the lifetime of 71 Ge, no simple explanation for the anomaly has been found. To check these calibration experiments, a dedicated experiment BEST was performed, utilizing a neutrino source of unprecedented intensity and a detector optimized to increase statistics while providing some information on counting rate as a function of distance from the source. The results BEST obtained are consistent with the earlier solar neutrino calibration experiments, and when combined with those measurements, yield a Ga anomaly with a significance of approximately 4𝜎, under conservative assumptions. But BEST found no evidence of distance dependence and thus no explicit indication of new physics. In this review we describe the extensive campaigns carried out by SAGE, GALLEX/GNO, and BEST to demonstrate the reliability and precision of their experimental procedures, including 71 Ge recovery, counting, and analysis. We also describe efforts to define uncertainties in the neutrino capture cross section, which now include estimates of effects at the ≲ 0.5% level such as radiative corrections and weak magnetism. With the results from BEST, an anomaly remains even if one retains only the transition to the 71 Ge ground state, whose strength is fixed by the known lifetime of 71 Ge. We then consider the new-physics solution most commonly suggested to resolve the Ga anomaly, oscillations into a sterile fourth neutrino, 𝜈 𝑒 → 𝜈 𝑠 . We find such a solution generates substantial tension with several null experiments, owing to the large mixing angle required. While this does not exclude such solutions – the sterile sector might include multiple neutrinos as well as new interactions – it shows the need for more experimental constraints, if we are to make progress in resolving the Ga and other low-energy neutrino anomalies. We conclude by consider the role future low-energy electron-capture sources could play in this effort.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Terahertz Nano-Imaging of Electronic Strip Heterogeneity in a Dirac Semimetal

Emerging topological semimetals offer promise of realizing topological electronics enabled by terahertz (THz) current persistent against impurity scattering. Yet most fundamental issues remain on how to image nanoscale conductivity inhomogeneity. Here in this paper we show noninvasive and contactless conductivity mapping at THz-nm limit of electronic heterogeneity and nanostrip junctions in a Dirac material ZrTe 5 . A clear Dirac Fermion density transition, manifested as the exclusive THz conductivity contrast, is quantitatively analyzed and profiled on both sides of the junction. This also allows the determination of variable junction width of ~25–220 nm, depending on the THz conductivity contrast of adjacent strips. The unique THz-nm contrast is absent in mid-infrared nano-imaging measurements since topological semimetals with small Fermi pockets exhibit a better matching of their plasma frequency and scattering rate to the THz spectral region. The first-principles calculations provide two compelling implications: the conductivity nanocontrast can be induced by a small anisotropic strain, even less than 0.5%, due to an extreme strain sensitivity in ZrTe 5 ; A nanoscale topological phase transition is realized across some junctions induced by the strain, between strong topological insulators (TIs) and weak TIs/Dirac semimetals (DSMs).

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

Mechanochemical synthesis and structural analysis of trivalent lanthanide and uranium diphenylphosphinodiboranates

Phosphinodiboranates (H 3 BPR 2 BH 3 – ) are a class of borohydrides that have merited a reputation as weakly coordinating anions, which is attributed in part to the dearth of coordination complexes known with transition metals, lanthanides, and actinides. We recently reported how K(H 3 BP t Bu 2 BH 3 ) exhibits sluggish salt elimination reactivity with f-metal halides in organic solvents such as Et 2 O and THF. Here we report how this reactivity appears to be further attenuated in solution when the t Bu groups attached to phosphorus are exchanged for R = Ph or H, and we describe how mechanochemistry was used to overcome limited solution reactivity with K(H 3 BPPh 2 BH 3 ). Grinding three equivalents of K(H 3 BPPh 2 BH 3 ) with UI 3 (THF) 4 or LnI 3 (Ln = Ce, Pr, Nd) allowed homoleptic complexes with the empirical formulas U(H 3 BPPh 2 BH 3 ) 3 (1), Ce(H 3 BPPh 2 BH 3 ) 3 (2), Pr(H 3 BPPh 2 BH 3 ) 3 (3), and Nd(H 3 BPPh 2 BH 3 ) 3 (4) to be prepared and subsequently crystallized in good yields (50–80%). Single-crystal XRD studies revealed that all four complexes exist as dimers or coordination polymers in the solid-state, whereas 1 H and 11 B NMR spectra showed that they exist as a mixture of monomers and dimers in solution. Treating 4 with THF breaks up the dimer to yield the monomeric complex Nd(H 3 BPPh 2 BH 3 ) 3 (THF) 3 (4-THF). XRD studies revealed that 4-THF has one chelating and two dangling H 3 BPPh 2 BH 3 – ligands bound to the metal to accommodate binding of THF. In contrast to the results with K(H 3 BPPh 2 BH 3 ), attempting the same mechanochemical reactions with Na(H 3 BPH 2 BH 3 ) containing the simplest phosphinodiboranate were unsuccessful; only the partial metathesis product U(H 3 BPH 2 BH 3 )I 2 (THF) 3 (5) was isolated in poor yields. Here, despite these limitations, our results offer new examples showing how mechanochemistry can be used to rapidly synthesize molecular coordination complexes that are otherwise difficult to prepare using more traditional solution methods.

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

Kohn-Sham density functional theory of Abelian anyons

In this work, we develop a density functional treatment of noninteracting Abelian anyons, which is capable, in principle, of dealing with a system of a large number of anyons in an external potential. Comparison with exact results for a few particles shows that the model captures the behavior qualitatively and semiquantitatively, especially in the vicinity of the fermionic statistics. We then study anyons with statistics parameter 1+1/n, which are thought to condense into a superconducting state. An indication of the superconducting behavior is the mean-field result that, for uniform density systems, the ground-state energy increases under the application of an external magnetic field independent of its direction. Our density functional theory based analysis does not find that to be the case for finite systems of anyons, which can accommodate a weak external magnetic field through density transfer between the bulk and the boundary rather than through transitions across effective Landau levels, but the “Meissner repulsion” of the external magnetic field is recovered in the thermodynamic limit as the effect of the boundary becomes negligible. We also consider the quantum Hall effect of anyons, and show that its topological properties, such as the charge and statistics of the excitations and the quantized Hall conductance, arise in a self-consistent fashion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

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↗
DOE OSTI2026

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↗
DOE OSTI2024

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↗