SEARCH · Engineering Papers

Results for “large displacement”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Search for long-lived heavy neutral leptons with displaced vertices in proton-proton collisions at $\sqrt{\mathrm{s}}$ =13 TeV

A search for heavy neutral leptons (HNLs), the right-handed Dirac or Majorana neutrinos, is performed in final states with three charged leptons (electrons or muons) using proton-proton collision data collected by the CMS experiment at √s = 13 TeV at the CERN LHC. The data correspond to an integrated luminosity of 138 fb -1 . The HNLs could be produced through mixing with standard model neutrinos ν. For small values of the HNL mass (<20 GeV) and the square of the HNL-ν mixing parameter (10 -7 –10 -2 ), the decay length of these particles can be large enough so that the secondary vertex of the HNL decay can be resolved with the CMS silicon tracker. The selected final state consists of one lepton emerging from the primary proton-proton collision vertex, and two leptons forming a displaced, secondary vertex. No significant deviations from the standard model expectations are observed, and constraints are obtained on the HNL mass and coupling strength parameters, excluding previously unexplored regions of parameter space in the mass range 1–20 GeV and squared mixing parameter values as low as 10 -7 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Polarity‐Driven Atomic Displacements at the 2D Mg 2 TiO 4 ‐MgO (001) Oxide Interface for Hosting Potential Interlayer Excitons

Abstract Interlayer excitons in solid‐state systems have emerged as candidates for realizing novel platforms ranging from excitonic transistors and optical qubits to exciton condensates. Interlayer excitons have been discovered in 2D transition metal dichalcogenides, with large exciton binding energies and the ability to form various van der Waals heterostructures. Here, an oxide system consisting of a single unit cell of Mg 2 TiO 4 on MgO (001) is proposed as a platform for hosting interlayer excitons. Using a combination of density functional theory (DFT) calculations, molecular beam epitaxy growth, and in situ crystal truncation rod measurements, it is shown that the Mg 2 TiO 4 ‐MgO interface can be precisely controlled to yield an internal electric field suitable for hosting interlayer excitons. The atoms in the polar Mg 2 TiO 4 layers are observed to be displaced to reduce polarity at the interface with the non‐polar MgO (001) surface. Such polarity‐driven atomic displacements strongly affect electrostatics of the film and the interface, resulting in localization of filled and empty band‐edge states in different layers of the Mg 2 TiO 4 film. The DFT calculations suggest that the electronic structure is favorable for localization of photoexcited electrons in the bottom layer and holes in the top layer, which may bind to form interlayer exciton states.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Identifying Heterogeneous Micromechanical Properties of Biological Tissues via Physics–Informed Neural Networks

The heterogeneous micromechanical properties of biological tissues have profound implications across diverse medical and engineering domains. However, identifying full-field heterogeneous elastic properties of soft materials using traditional engineering approaches is fundamentally challenging due to difficulties in estimating local stress fields. Recently, there has been a growing interest in data-driven models for learning full-field mechanical responses, such as displacement and strain, from experimental or synthetic data. However, research studies on inferring full-field elastic properties of materials, a more challenging problem, are scarce, particularly for large deformation, hyperelastic materials. Here, a physics-informed machine learning approach is proposed to identify the elasticity map in nonlinear, large deformation hyperelastic materials. This study reports the prediction accuracies and computational efficiency of physics-informed neural networks (PINNs) in inferring the heterogeneous elasticity maps across materials with structural complexity that closely resemble real tissue microstructure, such as brain, tricuspid valve, and breast cancer tissues. Further, the improved architecture is applied to three hyperelastic constitutive models: Neo-Hookean, Mooney Rivlin, and Gent. Furthermore, the improved network architecture consistently produces accurate estimations of heterogeneous elasticity maps, even when there is up to 10% noise present in the training data.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Progress on 3D SRF-based architecture for quantum computing

Superconducting radio frequency (SRF) cavities are excellent choices for storing and manipulating quantum information as quantum d-level systems (qudits) due to their exceptionally long lifetimes and large accessible Hilbert spaces. A common strategy to manipulate the states is to use a nonlinear element like a transmon. We present preliminary experimental results obtained with cavity displacements and selective number dependent arbitrary phase gates for universal qudit control, and its application towards High-energy physics (HEP) simulations and beyond. We discuss the advantages and challenges associated with building a 3D SRF architecture while maintaining long cavity lifetimes in the presence of lossy components. We show how the system coherence properties can be preserved by carefully engineering to minimize the participation of the long coherence modes in different loss channels, while ensuring sufficient quantum controllability. We further discuss the path towards building multi-qudit systems.

Romanenko, Alexander↗

DOE OSTI2023

Bicrystallography-informed Frenkel–Kontorova model for interlayer dislocations in strained 2D heterostructures

In recent years, van der Waals (vdW) heterostructures and homostructures, which consist of stacks of two-dimensional (2D) materials, have risen to prominence due to their association with exotic quantum phenomena originating from correlated electronic states harbored by them. Atomistic scale relaxation effects play an extremely important role in the electronic scale quantum physics of these systems, providing means of manipulation of these materials and allowing them to be tailored for emergent technologies. We investigate such structural relaxation effects in this work using atomistic and mesoscale models, within the context of twisted bilayer graphene — a well-known heterostructure system that features moiré patterns arising from the lattices of the two graphene layers. For small twist angles, atomic relaxation effects in this system are associated with the natural emergence of interface dislocations or strain solitons, which result from the cyclic nature of the generalized stacking fault energy (GSFE), that measures the interface energy based on the relative movement of the two layers. Here, in this work, we first demonstrate using atomistic simulations that atomic reconstruction in bilayer graphene under a large twist also results from interface dislocations, although the Burgers vectors of such dislocations are considerably smaller than those observed in small-twist systems. To reveal the translational invariance of the heterointerface responsible for the formation of such dislocations, we derive the translational symmetry of the GSFE of a 2D heterostructure using the notions of coincident site lattices (CSLs) and displacement shift complete lattices (DSCLs). The workhorse for this exercise is a recently developed Smith normal form bicrystallography framework. Next, we construct a bicrystallography-informed and frame-invariant Frenkel–Kontorova model, which can predict the formation of strain solitons in arbitrary 2D heterostructures, and apply it to study a heterostrained, large-twist bilayer graphene system. Our mesoscale model is found to produce results consistent with atomistic simulations. We anticipate that the model will be invaluable in predicting structural relaxation and for providing insights into various heterostructure systems, especially in cases where the fundamental unit cell is large and therefore, atomistic simulations are computationally expensive.

2D heterostructures↗

DOE OSTI2021

Floating platform effects on power generation in spar and semisubmersible wind turbines

The design and financing of commercial-scale floating offshore wind projects require a better understanding of how power generation differs between newer floating turbines and well-established fixed-bottom turbines. In floating turbines, platform mobility causes additional rotor motion that can change the time-averaged power generation. In this work, OpenFAST simulations examine the power generated by the National Renewable Energy Laboratory's 5-MW reference turbine mounted on the OC3-UMaine spar and OC4-DeepCWind semisubmersible floating platforms, subjected to extreme irregular waves and below-rated turbulent inflow wind from large-eddy simulations of a neutral atmospheric boundary layer. For these below-rated conditions, average power generation in floating turbines is most affected by two types of turbine displacements: an average rotor pitch angle that reduces power, caused by platform pitch; and rotor motion upwind-downwind that increases power, caused by platform surge and pitch. The relative balance between these two effects determines whether a floating platform causes power gains or losses compared to a fixed-bottom turbine; for example, the spar creates modest (3.1%–4.5%) power gains, whereas the semisubmersible creates insignificant (0.1%–0.2%) power gains for the simulated conditions. Furthermore, platform surge and pitch motions must be analyzed concurrently to fully capture power generation in floating turbines, which is not yet universal practice. Finally, a simple analytical model for predicting average power in floating turbines under below-rated wind speeds is proposed, incorporating effects from both the time-averaged pitch displacement and the dynamic upwind-downwind displacements.

17 WIND ENERGY↗

DOE OSTI2024

Influences of lidar scanning parameters on wind turbine wake retrievals in complex terrain

Abstract. Scanning lidars enable the collection of spatially distributed measurements of turbine wakes and the estimation of wake properties such as magnitude, extent, and trajectory. Lidar-based characterizations, however, may be subject to distortions due to the observational system. Distortions can arise from the resolution of the measurement points across the wake, the projection of the winds onto the beam, averaging along the beam probe volume, and intervening evolution of the flow over the scan duration. Using a large-eddy simulation and simulated measurements with a virtual lidar model, we assess how scanning lidar systems may influence the properties of the retrieved wake using a case study from the Perdigão campaign. We consider three lidars performing range-height indicator sweeps in complex terrain, based on the deployments of lidars from the Danish Technical University (DTU) and German Aerospace Center (DLR) at the Perdigão site. The unwaked flow, measured by the DTU lidar, is well-captured by the lidar, even without combining data into a multi-lidar retrieval. The two DLR lidars measure a waked transect from different downwind vantage points. In the region of the wake, the observation system reacts to the smaller spatial and temporal variations of the winds, allowing more significant observation distortions to arise. While the measurements largely capture the wake structure and trajectory over its 4–5 D extent, limited spatial resolution of measurement points and volume averaging lead to a quicker loss of the two lobes in the near wake, smearing of the vertical bounds of the wake (< 30 m), wake center displacements up to 10 m, and dampening of the maximum velocity deficit by up to a third. The virtual lidar tool, coupled with simulations, provides a means for assessing measurement capabilities in advance of measurement campaigns.

17 WIND ENERGY↗

DOE OSTI2024

How Does Vertical Wind Shear Influence Updraft Characteristics and Hydrometeor Distributions in Supercell Thunderstorms?

Abstract Vertical wind shear is known to affect supercell thunderstorms by displacing updraft hydrometeor mass downshear, thereby facilitating the storms’ longevity. Shear also impacts the size of supercell updrafts, with stronger shear leading to wider, less dilute, and stronger updrafts with likely greater hydrometeor production. To more clearly define the role of shear across different vertical layers on hydrometeor concentrations and displacements relative to supercell updrafts, a suite of idealized numerical model simulations of supercells was conducted. Shear magnitudes were systematically varied across the 0–1, 1–6, and 6–12 km AGL layers, while the thermodynamic environment was held fixed. Simulations show that as shear magnitude increases, especially from 1 to 6 km, updrafts become wider and less dilute with an increase in hydrometeor loading, along with an increase in the low-level precipitation area/rate and total precipitation accumulation. Even with greater updraft hydrometeor loading amid stronger shear, updrafts are more intense in stronger shear simulations due to larger thermal buoyancy owing to wider, less dilute updraft cores. Furthermore, downshear hydrometeor displacements are larger in environments with stronger 1–6-km shear. In contrast, there is relatively less sensitivity of hydrometeor concentrations and displacements to variations in either 0–1- or 6–12-km shear. Results are consistent across free tropospheric relative humidity sensitivity simulations, which show an increase in updraft size and hydrometeor mass with increasing free tropospheric relative humidity owing to a reduction in entrainment-driven dilution for wider updrafts in moister environments. Significance Statement Rotating thunderstorms, known as supercells, are able to persist for multiple hours. One common explanation is that large changes in wind speed and/or direction with height, or shear, transport rain/hail away from supercell updrafts, supporting their maintenance. The strong shear within supercell environments, however, may also lead to greater rail/hail amounts, thereby leading to weaker storms due to this extra mass of water/ice within updrafts. Furthermore, the impact of shear across different height layers on supercell rain/hail characteristics has not been thoroughly investigated. In this study, computer simulations of supercells were conducted to determine that shear occurring between 1 and 6 km above ground level has a large impact on rain/hail distribution in supercells and that stronger shear in this layer leads to wider/stronger supercells with greater rain/hail accumulations at the surface. Additionally, some of the extra mass of water/ice is transported farther away from updrafts due to the stronger environmental storm-relative winds.

Meteorology & Atmospheric Sciences↗

DOE OSTI2022

Towards a Unified Setup to Simulate Mid‐Latitude and Tropical Mesoscale Convective Systems at Kilometer‐Scales

Abstract Mesoscale convective systems (MCSs) are the main source of precipitation in the tropics and parts of the mid‐latitudes and are responsible for high‐impact weather worldwide. Studies showed that deficiencies in simulating mid‐latitude MCSs in state‐of‐the‐art climate models can be alleviated by kilometer‐scale models. However, whether these models can also improve tropical MCSs and whether we can find model settings that perform well in both regions is understudied. We take advantage of high‐quality MCS observations collected over the Atmospheric Radiation Measurement (ARM) facilities in the US Southern Great Plains (SGP) and the Amazon basin near Manaus (MAO) to evaluate a perturbed physics ensemble of simulated MCSs with 4 km horizontal grid spacing. A new model evaluation method is developed that enables to distinguish biases stemming from spatiotemporal displacements of MCSs from biases in their reflectivity and cloud shield. Amazon MCSs are similarly well simulated across these evaluation metrics than SGP MCSs despite the challenges anticipated from weaker large‐scale forcing in the tropics. Generally, SGP MCSs are more sensitive to the choice of model microphysics, while Amazon cases are more sensitive to the planetary boundary layer (PBL) scheme. Although our tested model physics combinations had strengths and weaknesses, combinations that performed well for SGP simulations result in worse results in the Amazon basin and vice versa. However, we identified model settings that perform well at both locations, which include the Thompson and Morrison microphysics coupled with the Yonsei University (YSU) PBL scheme and the Thompson scheme coupled with the Mellor‐Yamada‐Janjic PBL scheme.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Electric field distribution in a non-self-sustained RF discharge with ionization generated by Ns discharge pulses

Abstract Electric field in a capacitively coupled, non-self-sustained RF discharge nitrogen plasma with external ionization generated by high-voltage ns pulses has been measured by ps Electric field induced second harmonic generation (EFISH). The measurements are made both in the bulk of the plasma and in the sheaths, using the absolute calibration by the Laplacian field between two plane electrodes. The results are compared with the kinetic modeling calculations. The RF electric field in the sheaths is significantly higher compared to that in the plasma, due to the displacement of the electrons by the drift oscillations and the resultant plasma self-shielding. However, the kinetic modeling predictions indicate that the electron impact ionization in the sheaths is largely ineffective, due to the low electron density. The reduction of the electric field in the plasma by the self-shielding in the sheaths is moderate, such that the energy is coupled to the plasma by the below-breakdown RF field. The peak RF field in the plasma is in the range of 15–25 Td, indicating the efficient vibrational excitation of N 2 by electron impact. The modeling predictions suggest that the targeted vibrational excitation of molecular species in a non-self-sustained RF discharge with external ionization is scalable to high pressures, electron densities, and discharge powers. The present approach can be extended to the vibrational excitation of other molecular species where the vibrational relaxation is relatively slow, such as CO, CO 2 , and H 2 .

Physics↗

DOE OSTI2022

Understanding the tool influence function during sub-aperture belt-on-wheel glass polishing

The tool influence function (TIF) during sub-aperture belt-on-wheel polishing has been evaluated as a function of various process conditions (belt use/wear, dwell time, displacement, belt velocity, and wheel modulus and diameter) on fused silica glass workpieces using C e O 2 polishing media. TIF spots are circular or elliptical in shape with a largely flat bottom character. Furthermore, the volumetric removal rate varies significantly with belt use (or wear), stabilizing after ~ 15 m i n of use. A modified Preston model, where the pressure dependence is adjusted using a different scaling of the wheel modulus ( E w 0.5 ), largely predicts the volumetric removal rate over the range of process conditions evaluated. The relatively high volumetric removal rate of 30 - 60 m m 3 / h using a fixed C e O 2 -in-resin-host belt offers a rapid, and hence, more economical, initial polish of aspheric and freeform optics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Optimizing Grain Boundary Structures with LAMMPS Using Evolutionary Algorithms

Grain boundary structure optimization is an important part of materials modeling. Current methods for grain boundary structure optimization involve inefficient, time-consuming processes that do not fully explore the interface parameter space. Evolutionary algorithms have recently been demonstrated to be effective at determining both stable and metastable grain boundary interface structures. In this work, we demonstrate the use of GBOpt, a grain boundary structure optimization software designed to use the Large-scale Atomic/Molecular Massively Parallel Simulation (LAMMPS) software to efficiently determine grain boundary structures. We demonstrate that a only a few manipulations, namely atom insertion, atom removal, and relative grain displacement, are sufficient to explore much of the grain boundary structure parameter space. The efficacy of this approach is demonstrated on an FCC Ni system, and a BCC Fe system. The computational cost is compared against the gamma-surface sampling approach to demonstrate performance improvement.

Evolutionary algorithms↗

DOE OSTI2021

Novel 2D velocity estimation method for large transient events in plasmas

Dynamics of fast transient events are challenging to be analyzed with high time resolution. Such events can occur in fusion plasmas such as the filaments during edge-localized modes (ELMs). Here, we present a robust method—the spatial displacement estimation—for estimating the displacements of structures with fast dynamics from high spatial and time resolution imaging diagnostics [e.g., gas-puff imaging (GPI)] with sampling time temporal resolution. First, a background suppression method is shown, which suppresses the slowly time-evolving and spatially non-uniform background in the signal. In the second step, a two-dimensional polynomial trend subtraction method is presented to tackle the remaining polynomial order trend in the signal. After performing these pre-processing steps, the spatial displacement of the propagating structure is estimated from the two-dimensional spatial cross-correlation coefficient function calculated between consecutive frames. The method is tested for its robustness and accuracy by simulated Gaussian events and spatially displaced random noise. An example application of the method is presented on propagating ELM filaments measured by the GPI system on the National Spherical Torus Experiment spherical tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Bio-inspired bistable piezoelectric energy harvester for powering animal telemetry tags: Conceptual design and preliminary experimental validation

Here, this paper presents the conceptual design, preliminary experimental validation, and performance evaluation of a novel bio-inspired bi-stable piezoelectric energy harvester for self-powered animal telemetry tags. The overall conceptual design, which includes a bio-inspired attachment and a bi-stable piezoelectric energy harvester, is introduced firstly with a specific application example of marine fish tracking. The self-powered telemetry tag can be externally deployed on fish (dorsal fin) to monitor fish habitats, population, and underwater environment. Inspired by the Venus flytrap's rapid shape transition, a bi-stable piezoelectric energy harvester is developed to scavenge energy from fish maneuvering and the surrounding fluid flow for a sustainable power supply. The bistability of the harvester is characterized by the measured force-displacement curve and double potential wells. A bluff body is integrated to the free end of the bistable piezoelectric energy harvester to enhance the structure-fluid interaction for the large-amplitude snap-through vibrations and higher voltage output. Controlled laboratory experiments are conducted in a water tank on the bio-inspired bi-stable piezoelectric energy harvester using a servo motor system to simulate fish swing motion at various conditions to evaluate the power generation performance. The preliminary underwater experimental results demonstrated that the proposed bio-inspired bi-stable piezoelectric energy harvester could effectively convert fish swing motions into electricity. The device collected 17.25 mJ of energy over 130 s under a peak-to-peak swing angle of 30° at 1.5 Hz in the capacitor charging experiments.

Animal tracking↗

DOE OSTI2022

Epitaxial intergrowths and local oxide relaxations in natural bixbyite Fe 2– x Mn x O 3

The scattering pattern of a crystal obeys the symmetry of the crystal structure through the corresponding Laue group. This is usually also true for the diffuse scattering, containing information about disorder, but here a case is reported where the diffuse scattering is of lower symmetry than the parent crystal structure. The mineral bixbyite has been studied by X-ray and neutron scattering techniques since 1928 with some of the most recent studies characterizing the low-temperature transition to a magnetically disordered spin-glass state. However, bixbyite also exhibits structural disorder, and here single-crystal X-ray and neutron scattering is used to characterize the different modes of disorder present. One-dimensional rods of diffuse scattering are observed in the cubic mineral bixbyite, which break the expected symmetry of the scattering pattern. It is shown that this scattering arises from epitaxial intergrowths of the related mineral, braunite. The presence of this disorder mode is found to be directly observable as well-defined residuals in the average structure refined against the Bragg diffraction. An additional three-dimensional diffuse scattering component is observed in neutron scattering data, which is shown to originate from the substitutional disorder on the Fe/Mn sites. This occupational disorder gives rise to local relaxations of the oxide sublattice, and the pattern of oxide displacements can be rationalized based on crystal-field theory. The combined use of neutron and X-ray single-crystal scattering techniques highlights their great complementarity. In particular, the large sample requirements for neutron scattering experiments prove to be an obstacle in solving the intergrowth disorder due to several growth orientations, whereas for X-ray scattering the one-dimensional nature of the intergrowth disorder renders solving this a more tractable task. On the other hand, the oxide relaxations cannot be resolved using X-rays due to the low Mn/Fe contrast. By combining the two approaches both types of disorder have been characterized.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Optimization of COMPASS-U Support Structure

The global support structure for Compass Upgrade (COMPASS-U) employs many bolted and pinned components to improve the manufacturability of the whole structures. This is a departure from other designs used for highly loaded tokamaks, like Alcator C-Mod which used large forgings. Modeling of the numerous connections and interfaces was a challenge. This article will describe how to optimize the C-frame support structure with the goal to eliminate differential displacements at those interfaces and, importantly, minimize toroidal mutual displacement at the sliding joint of the toroidal field (TF) coils. The global support structure consists of 16 C-frames, and it can be simplified as a 1/8th cyclic symmetry model due to the eight vertical rods mounting that connect bottom CS and PF1-3 assemblies to the C-frames. Structure variants with different connection webs between two neighboring C-frames have been analyzed. The preloading of bolts can be split into several preload substeps to simulate the real assembly process of the COMPASS-U support structure. Additionally, the TF crown and the wedge plate can be further improved to become a solid part to better hold the TF coils.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Bio-inspired Bistable Piezoelectric Energy Harvester for Powering Animal Telemetry Tags: Concept Design and Preliminary Experimental Validation

This paper presents the concept design, preliminary experimental validation, and performance evaluation of a novel bio-inspired bi-stable piezoelectric energy harvester for self-powered animal telemetry tags. The overall concept design, which includes a bio-inspired attachment and a bi-stable piezoelectric energy harvester, is introduced firstly with a specific application example of marine fish tracking. The self-powered telemetry tag can be externally deployed on fish (dorsal fin) to monitor fish habitats, population, and underwater environment. Inspired by the Venus flytrap’s rapid shape transition, a bistable piezoelectric energy harvester is developed to scavenge energy from fish maneuvering and the surrounding fluid flow for a sustainable power supply. The bistability of the harvester is characterized by the measured force-displacement curve and double potential wells. A bluff body is integrated to the free end of the bistable piezoelectric energy harvester to enhance the structure-fluid interaction for the large-amplitude snap-through vibrations and higher voltage output. Controlled laboratory experiments are conducted in a water tank on the bio-inspired bi-stable piezoelectric energy harvester using a servo motor system to simulate fish swing motion at various conditions to evaluate the power generation performance. The preliminary underwater experimental results demonstrated that the proposed bio-inspired bi-stable piezoelectric could effectively converter fish swing motions into electricity. 17.25 mJ energy was obtained within 130 s under the peak-to-peak swing angle of 30o at 1.5 Hz in the capacitor charging experiments.

Qian, Feng↗

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↗