Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable geometry”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

XMM-Newton and TESS observations of the highly variable polar V496 UMa

We studied the temporal and spectral behavior of V496 UMa from the optical to the X-ray regimes. Methods. We used archival XMM-Newton and TESS observations obtained in 2017 and 2019, respectively, to perform a spectral and timing analysis of the highly variable magnetic CV. Results. The light curves of the TESS and XMM-Newton satellites reveal a double-humped pattern modulated with the periodicity of 91.058467 ± 0.00001 min. V496 UMa displays a two-pole accretion geometry in the high accretion state. X-ray spectra from these regions are composed of thermal plasma radiation and soft blackbody components with almost identical temperatures and a total accretion rate of $\dot{M}$ = 1.4(8) × 10 -11 M ⊙ yr -1 . The X-ray centers of the humps show longitudinal shifts of -18° and 4°, and shifts around photometric phase zero of -172° and -186°, for the main hump and second hump, respectively. The long-term ZTF light curves reveal high and low accretions states. Low-state ZTF and SDSS photometric data are consistent with an 0.8 M ⊙ white dwarf at 10 000 K and a main-sequence donor star with a spectral type of M5.0 at a Gaia determined distance of 758 pc. Conclusions. V496 UMa is a very bright polar in X-rays when it is in the high state. Due to its unusual geometric structure, mass accretion onto the second accretion pole is interrupted occasionally. This discontinuous behavior does not follow a certain pattern in time and has been observed so far only in the high state. The X-ray light curves display clear evidence of an accretion stream at the photometric phase of φ = 0.81, which does not show up in optical light curves. An accurate period was derived using the combined TESS and XMM-Newton data, which differs by 3.8 σ from published results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Experimental and Numerical Investigation of a Miniature Additively Manufactured Vortex Tube

Ranque–Hilsch vortex tubes are simple devices that can produce a cooling effect using compressed air. A key advantage of vortex tubes is the lack of moving solid parts; however, their efficiencies are relatively low. The present study focuses on the development of a miniature variable-diameter tube using additive manufacturing. A metal-based 3D printing technique was utilized to fabricate this vortex tube monolithically. Computational fluid dynamics simulations employing software STAR-CCM + with a compressible Reynolds-Averaged Navier–Stokes (RANS) approach and the elliptic-blending lag k-epsilon turbulence model have been applied to model thermofluid processes inside the vortex tube, to good agreement with the experiment. In conclusion, a temperature decrease of 13.3 °C and a cooling power of approximately 4 W were experimentally achieved with a pressure ratio of 4 in the air at normal conditions. This result shows promise for the goal of utilizing additive manufacturing to design and build complex-geometry vortex tubes intended for use with cryogenic fluids.

08 HYDROGEN↗
DOE OSTI2026

Solid-State Syntheses, Crystallographic Spatial Disorders, Thermal Behaviors, and Bandgaps of Hybrid Organic-Inorganic Manganese Halides: A2Mn(Cl/Br)4 (A = NH4+, C(NH2)3+, & C3H4N2+)

By systematically optimizing solid-state synthesis via the concerted use of differential scanning calorimetry and in-situ variable-temperature X-ray diffraction, we report the discovery of four new hybrid organic-inorganic manganese halides A2Mn(Cl/Br)4 (A = NH4+, C(NH2)3+, C3H4N2+), with emphasis on how the organic fragment geometry, polarity, and electron-donating properties influence crystallographic spatial disorders, thermal behaviors, and optical band gaps. With a nonpolar tetrahedral cation, as in (NH4)2MnCl4 (P21/c, mP30) or (NH4)2MnBr4 (C2/c, mC180), the direct optical bandgaps (4.36–4.28 eV) are wider than those with an organic nonpolar planar geometry, as in (C(NH2)3)2MnBr4 (P21/c, mP300, 4.07 eV), or an organic polar planar geometry, as in (C3H4N2+)2MnBr4 (I41/a, tI384, 4.05 eV). Furthermore, the organic components affect the spatial arrangement and dimensionality of the resulting inorganic frameworks comprised of Mn(Cl/Br)6 octahedra, with varying distortions, or spatially disordered MnBr4 tetrahedra. These resulting trends, coupled with the systematic investigation into synthesis, and the motivation behind elucidating the nuanced role of the organic cation, altogether expand upon the phase-space of hybrid materials beyond perovskites and demonstrate the potential for a rational design of hybrid materials with tailored optoelectronic functionalities for advanced energy applications.

Lee, Shannon J.↗
DOE OSTI2021

Sculpting crystals one Burgers vector at a time: Toward colloidal lattice robot swarms

Plastic deformation of crystalline materials with isotropic particle attractions proceeds by the creation and migration of dislocations under the influence of external forces. If dislocations are produced and migrated under the action of local forces, then material shape change can occur without the application of surface forces. We investigate how particles with variable diameters can be embedded in colloidal monolayers to produce dislocations on demand. Additionally, we find in simulation that when embedded clusters of variable diameter particles are taken through multiple cycles of swelling and shrinking, large cumulative plastic slip is produced by the creation and biased motion of dislocation pairs in the solid for embedded clusters of particular geometries. In this way, dislocations emitted by these clusters (biased “dislocation emitters”) can be used to reshape colloidal matter. Our results are also applicable to larger-scale swarms of robotic particles that organize into dense ordered two-dimensional (2D) arrangements. We conclude with a discussion of how dislocations fulfill for colloids the role sought by “metamodules” in lattice robotics research and show how successive applications of shear as a unit operation can produce shape change through slicing and swirling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Disentangling ferroelectric domain wall geometries and pathways in dynamic piezoresponse force microscopy via unsupervised machine learning

In this work, domain switching pathways in ferroelectric materials visualized by dynamic piezoresponse force microscopy (PFM) are explored via variational autoencoder, which simplifies the elements of the observed domain structure, crucially allowing for rotational invariance, thereby reducing the variability of local polarization distributions to a small number of latent variables. For small sampling window sizes the latent space is degenerate, and variability is observed only in the direction of a single latent variable that can be identified with the presence of domain wall. For larger window sizes, the latent space is 2D, and the disentangled latent variables can be generally interpreted as the degree of switching and complexity of domain structure. Applied to multiple consecutive PFM images acquired while monitoring domain switching, the polarization switching mechanism can thus be visualized in the latent space, providing insight into domain evolution mechanisms and their correlation with the microstructure.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Geometric and Material Variability Influences Stress States Relevant to Coastal Permafrost Bluff Failure

Scientific knowledge and engineering tools for predicting coastal erosion are largely confined to temperate climate zones that are dominated by non-cohesive sediments. The pattern of erosion exhibited by the ice-bonded permafrost bluffs in Arctic Alaska, however, is not well-explained by these tools. Investigation of the oceanographic, thermal, and mechanical processes that are relevant to permafrost bluff failure along Arctic coastlines is needed. We conducted physics-based numerical simulations of mechanical response that focus on the impact of geometric and material variability on permafrost bluff stress states for a coastal setting in Arctic Alaska that is prone to toppling mode block failure. Our three-dimensional geomechanical boundary-value problems output static realizations of compressive and tensile stresses. We use these results to quantify variability in the loci of potential instability. We observe that niche dimension affects the location and magnitude of the simulated maximum tensile stress more strongly than the bluff height, ice wedge polygon size, ice wedge geometry, bulk density, Young's Modulus, and Poisson's Ratio. Our simulations indicate that variations in niche dimension can produce radically different potential failure areas and that even relatively shallow vertical cracks can concentrate displacement within ice-bonded permafrost bluffs. These findings suggest that stability assessment approaches, for which the geometry of the failure plane is delineated a priori, may not be ideal for coastlines similar to our study area and could hamper predictions of erosion rates and nearshore sediment/biogeochemical loading.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Calibrating hypersonic turbulence flow models with the HIFiRE-1 experiment using data-driven machine-learned models.

In this paper we study the efficacy of combining machine-learning methods with projection-based model reduction techniques for creating data-driven surrogate models of computationally expensive, high-fidelity physics models. Such surrogate models are essential for many-query applications e.g., engineering design optimization and parameter estimation, where it is necessary to invoke the high-fidelity model sequentially, many times. Surrogate models are usually constructed for individual scalar quantities. However there are scenarios where a spatially varying field needs to be modeled as a function of the model’s input parameters. Here we develop a method to do so, using projections to represent spatial variability while a machine-learned model captures the dependence of the model’s response on the inputs. The method is demonstrated on modeling the heat flux and pressure on the surface of the HIFiRE-1 geometry in a Mach 7.16 turbulent flow. The surrogate model is then used to perform Bayesian estimation of freestream conditions and parameters of the SST (Shear Stress Transport) turbulence model embedded in the high-fidelity (Reynolds-Averaged Navier–Stokes) flow simulator, using shock-tunnel data. The paper provides the first-ever Bayesian calibration of a turbulence model for complex hypersonic turbulent flows. We find that the primary issues in estimating the SST model parameters are the limited information content of the heat flux and pressure measurements and the large model-form error encountered in a certain part of the flow.

42 ENGINEERING↗
DOE OSTI2021

Bayesian sparse learning with preconditioned stochastic gradient MCMC and its applications

Deep neural networks have been successfully employed in an extensive variety of research areas, including solving partial differential equations. Despite its significant success, there are some challenges in effectively training DNN, such as avoiding overfitting in over-parameterized DNNs and accelerating the optimization in DNNs with pathological curvature. Here, we propose a Bayesian type sparse deep learning algorithm. The algorithm utilizes a set of spike-and-slab priors for the parameters in the deep neural network. The hierarchical Bayesian mixture will be trained using an adaptive empirical method. That is, one will alternatively sample from the posterior using preconditioned stochastic gradient Langevin Dynamics (PSGLD), and optimize the latent variables via stochastic approximation. The sparsity of the network is achieved while optimizing the hyperparameters with adaptive searching and penalizing. A popular SG-MCMC approach is Stochastic gradient Langevin dynamics (SGLD). However, considering the complex geometry in the model parameter space in nonconvex learning, updating parameters using a universal step size in each component as in SGLD may cause slow mixing. To address this issue, we apply a computationally manageable preconditioner in the updating rule, which provides a step-size parameter to adapt to local geometric properties. Moreover, by smoothly optimizing the hyperparameter in the preconditioning matrix, our proposed algorithm ensures a decreasing bias, which is introduced by ignoring the correction term in the preconditioned SGLD. According to the existing theoretical framework, we show that the proposed algorithm can asymptotically converge to the correct distribution with a controllable bias under mild conditions. Numerical tests are performed on both synthetic regression problems and learning solutions of elliptic PDE, which demonstrate the accuracy and efficiency of the present work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Accelerating Additive Manufacturing Process Design for Energy Conversion Materials using In-situ Sensing and Machine Learning

One promising candidate for manufacturing of the bismuth Telluride thermoelectric legs is laser powder bed fusion (LPBF) additive Manufacturing (AM). AM processing parameters highly influence the material properties, however current processing parameter development methods in AM are costly and time consuming. In-situ sensors allow for the capture of physically relevant process information on a layer-by-layer basis and will be used to aide process development. To optimize the AM process for the best thermoelectric performance, process variables, in-situ process sensor data and ex-situ material characterization data are collected. Several different interpretable machine learning (ML) approaches are used, and the performance of each method are assessed. Significant input process variables include laser focus, hatch spacing and laser power. The best performing models are used to determine the manufacturing parameters that maximize the power factor. AM of bismuth telluride material provides the ability to create complex geometries enabling more efficient energy conversion.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Installation and Demonstration of CNC Machining for Mechanical Test Specimens in the IFEL Hot Cell

This report documents the installation and demonstration of computer numerical control (CNC) machining capabilities in the Irradiated Fuels Examination Laboratory hot cell facility at Oak Ridge National Laboratory (ORNL). A modified Tormach PCNC 440 mill was integrated into the hot cell with custom fixturing, fines management, and manipulator-compatible interfaces to enable the fabrication of axial tension test (ATT) and ring tension test (RTT) specimens from irradiated cladding. The first irradiated specimens machined included ATT and RTT geometries harvested from the high-burnup 6XV fuel rod. Dimensional inspections confirmed that machined specimens met the ±0.025 mm tolerance envelope established in prior development; deviations were consistent with expected measurement scatter and inherent specimen variability, such as wall thickness gradients and eccentricity. Comparisons with out-of-cell metrology confirmed that in-cell machining performance aligns with baseline scatter observed under ideal conditions. This work establishes reproducible, end-to-end specimen preparation at ORNL, directly supporting the US Department of Energy’s Accident-Tolerant Fuel program by enabling reliable, traceable mechanical testing of irradiated cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Evaluation and development of flow condensation correlations using the data from low GWP refrigerants in an axial micro-fin aluminum tube

To mitigate global warming, the world is transitioning to refrigerants with low global warming potential (GWP). Supporting this shift requires a model that can accurately predict the heat transfer and pressure drop of new refrigerants, crucial for designing efficient heat exchangers. Existing models, however, are largely based on currently deployed refrigerants and primarily developed for unexpanded micro-fin tubes with spiral angles of 6° to 30°. Their applicability to new refrigerants, especially in expanded micro-fin tubes, is uncertain. This study assesses the performance of four well-known condensation models for six emerging refrigerants—R-32, R-454B, R-454C, R-455A, R-1234yf, and R-1234ze(E)—against experimental data. Initially, the Han and Lee (2005) model shows the best prediction accuracy with a mean absolute deviation (MAD) of 22.1 %. To enhance the accuracy of heat transfer models for new refrigerants and geometries with large temperature glides, two approaches are proposed. Here, the first approach applies a simple correction factor, reducing the MAD of the Cavallini et al. (2009) model from 68.2 % to 15.4 %. The second approach uses the variable metric method for minimization, fitting new constants to the data. This optimization results in the Kedzierski and Goncalves (1997) model achieving the highest accuracy, with a MAD of 13.1 %. For pressure drop models, the Cavallini et al. (1997) model is the most accurate with a MAD of 6.4 %, followed by the Haraguchi et al. (1993) model with a MAD of 9.4 %. Due to its simplicity, the Haraguchi et al. (1993) model is a practical option for predicting frictional pressure drop.

Axial micro-fin tubes↗
DOE OSTI2022

MOOSE Reactor Module Meshing Enhancements to Support Reactor Physics Analysis

The U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program develops an integrated suite of advanced reactor physics tools built upon the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Each code generally requires an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The meshing process is often burdensome for the complex geometries present in reactors due to lack of easy-to-use, open-source meshing tools. To address the bottleneck associated with meshing complex geometries found in nuclear reactors, NEAMS initiated the development of the MOOSE Reactor Module starting in FY21. The Reactor Module builds off the existing MOOSE Mesh System to include targeted meshing capabilities such as the ability to generate hexagonal pin cells, assemblies with ducts, rotating control drums, cores, peripheral zones around a core, as well as the automatic labeling (“reporting IDs”) of pin, assembly, and planar zones to simplify post-processing of results. As a Physics Module in MOOSE, the Reactor Module is open-source, available with any MOOSE installation, directly compatible with MOOSE-based tools, and can be invoked from MOOSEbased applications to generate meshes. Functionality from the Reactor Module has been applied to several advanced reactor concepts to demonstrate user workflow improvements and accuracy. The primary objective of the Reactor Module is to improve useability of MOOSEbased tools by streamlining mesh generation and output inspection processes. During FY22, the functionality of the Reactor Module (and accompanying Mesh System) has been expanded based on user needs. First, the Reactor Geometry Mesh Builder capability developed primarily in FY21 has been refactored and merged to the public MOOSE repository. This capability wraps underlying Reactor Module mesh generators into a “Pin – Assembly – Core” workflow appropriate for conventional Cartesian and hexagonal geometries, and notably assigns material IDs during mesh generation stage and generates only the minimal number of blocks needed in order to reduce computational burden. Biasing and boundary layer options have been added to the base mesh generators as required by thermal hydraulics solvers. The reporting ID functionality has been expanded to differentiate ring-wise and azimuthal sectors within a pin for use with depletion algorithms, and VectorPostProcessor and Reporter objects are now available to integrate solution variables across zones based on ID combinations. Functionality to trim hexagonal meshes along the center or periphery has been developed so users may leverage symmetry and reflective boundary conditions to reduce the mesh size. A flexible and powerful tool to fill the space between two sidesets has been introduced to the framework and can be used for transition layers such as stitching two assemblies together with different numbers of pins, or for complex geometries which do not follow conventional Cartesian/hexagonal patterns. Finally, additional verification problems were performed with NEAMS physics tools in complement with existing NEAMS work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Improved Chopper System for a Direct Geometry Spectrometer at the Second Target Station

This report describes a chopper cascade suitable for a cold direct geometry spectrometer at a pulsed neutron source with a low repetition rate such as the Second Target Station (STS) currently being planned at Oak Ridge National Laboratory (ORNL). The chopper cascade offers the following capabilities: selection of a variable number of pulses per frame (flexible pulse suppression), selection of pulses with increasing timing separation within one frame, and an increasing pulse width ∆t within one frame, in order to achieve an even energy resolution over several pulses at different wavelengths across the frame.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Area-Preserving Mapping of 3D Carotid Ultrasound Images Using Density-Equalizing Reference Map

Carotid atherosclerosis is a focal disease at the bifurcations of the carotid artery. To quantitatively monitor the local changes in the vessel-wall-plus-plaque thickness (VWT) and compare the VWT distributions for different patients or for the same patients at different ultrasound scanning sessions, a mapping technique is required to adjust for the geometric variability of different carotid artery models. In this work, we propose a novel method called density-equalizing reference map (DERM) for mapping 3D carotid surfaces to a standardized 2D carotid template, with an emphasis on preserving the local geometry of the carotid surface by minimizing the local area distortion. The initial map was generated by a previously described arc-length scaling (ALS) mapping method, which projects a 3D carotid surface onto a 2D non-convex L-shaped domain. A smooth and area-preserving flattened map was subsequently constructed by deforming the ALS map using the proposed algorithm that combines the density-equalizing map and the reference map techniques. This combination allows, for the first time, one-to-one mapping from a 3D surface to a standardized non-convex planar domain in an area-preserving manner. Evaluations using 20 carotid surface models show that the proposed method reduced the area distortion of the flattening maps by over 80% as compared to the ALS mapping method.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Polymorphism and Phase Control in Dion–Jacobson 2D 3-(Aminomethyl)piperidinium-Based Metal Iodide Perovskites

Two-dimensional halide perovskites exhibit rich structural diversity and tunable optoelectronic properties, making them promising materials for energy, sensing, and photonic applications. In this work, we report the structural mapping of four distinct polymorphs, γ (P2 1 /c), β (P4/mbm), α (P4/mmm), and δ (Cmcm) in the two-dimensional iodoplumbate, iodostannate, and iodogermanate perovskite so-called Dion-Jacobson series (3AMP)MI4 (M = Sn, Pb, Ge), where 3AMP is 3-(aminomethyl)piperidinium. The phases exhibit systematic evolution in octahedral distortion, lattice symmetry, and metal–halide geometry, enabling structural control over optoelectronic properties. Notably, the α-phase of (3AMP)SnI 4 represents a rare, ambient-stable, high-symmetry structure for Sn-based perovskites, without a phase transition down to 100 K. Variable-temperature single-crystal diffraction, powder x-ray diffraction (PXRD), and calorimetry reveal metal- and temperature-dependent polymorph interconversions, including the emergence of long range supercell reflections in Pb-rich compositions at low temperature. Optical spectroscopy and photoelectron yield spectroscopy confirm band gap tunability and band alignment trends, highlighting symmetry-dependent shifts and anomalous band gap bowing in mixed-metal systems, verified by electronic structure calculations. Calculations additionally indicate that the higher symmetry phases have reduced electron and hole effective masses compared to the lower symmetry phases.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Deep Learning enabled spectral energy conversion for in situ exposure measurements

A detector-specific deep learning (DL) approach is presented for spectra-to-exposure conversion using large-format sodium iodide (NaI(Tl)) detectors deployed for in situ environmental radiation measurements in emergency response scenarios. Accurate determination of exposure from NaI spectra is challenging due to poor energy resolution, partial energy absorption, and the strong sensitivity of traditionally deployed analytical conversion methods to calibrated source geometry and pre-deployment assumptions. Here, to address these limitations, a multi-layer perceptron model was trained on a hybrid in situ /Monte Carlo dataset constructed to span a broad range of photon energies, spatial extents, and realistic deployment variability, representative of general in situ emergency response conditions. The DL model was evaluated against commonly fielded analytical approaches under matched simulation conditions, including a single-factor method, a G-function method, and a modeled pressurized ion chamber (PIC) baseline. This study was intentionally computational in scope to enable controlled, like-for-like comparisons between conversion techniques while minimizing confounding real-world variability. Comparison to the modeled PIC provides contextual benchmarking and is not intended as a field inter-comparison with deployed instruments. Across the evaluated 20 keV to 3 MeV energy range, the DL approach consistently exhibited higher accuracy and reduced variance relative to the analytical methods against a deterministically calculated exposure. This may indicate improved robustness to spectral complexity without reliance on source-, geometric-, or spectral region-specific optimization. While results do not represent real-world validation, the presented work demonstrates that deep learning may effectively learn the nonlinear detector response-to-exposure relationship for asymmetric NaI(Tl) detectors and offers a promising pathway for improving in situ exposure estimation using spectroscopic systems already integrated into initial real-time emergency response operations.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗