Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Comparative framework”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

REDOTHERM (Redox Countercurrent Thermodynamic Limits Model) [SWR-24-88]

REDOTHERM is an open-source, MATLAB-based thermodynamic modeling framework developed to evaluate the performance of redox-active materials for water (H2O) and carbon dioxide (CO2) splitting. It includes models of all major unit operations and supports comparative analysis of different redox-active material candidates. The model is tailored for systems of moving oxide under a parallel/cocurrent flow (PF) and countercurrent flow (CF) configurations. Unvalidated mixed flow reactor (MFR, also known as CSTR) model is also included as an optional addition.

Lidor, Alon [National Renewable Energy Laboratory ↗
DOE OSTI2026

Denudation, solute export, landscape evolution modeling, and geographic information system data for the East River watershed, Colorado, USA (2020-2024)

This data package contains geographic information system (GIS) layers and tabular datasets associated with the study of lithologic controls on denudation, solute export, carbon-scaling relationships, and transient landscape evolution in the East River watershed near Crested Butte, Colorado, USA. The package includes GIS layers used to produce the Figure 2 map, including drainage, hillshade, lithology, sample locations, and basin polygons, together with comma-separated value (CSV) tables and matching CSV data dictionaries. One group of tables reports sample-level and catchment-level information for river-sediment samples analyzed for in situ-produced cosmogenic beryllium-10 (10Be), including sample names, outlet elevations, geographic coordinates, upstream drainage area, rock-type classes, production-rate scaling scheme, analyzed nuclide, catchment-averaged denudation rates, and associated lower and upper analytical uncertainties. Sample and catchment attributes provide the basis for comparing denudation rates across intrusive, shale, sedimentary, and mixed-lithology settings. A second group of tables reports supporting information for landscape-evolution modeling and the mapped geologic framework of the study area. Included files list parameter values and definitions for the two-phase landscape-evolution simulations, summarize full-domain model erosion fluxes and topographic metrics for different simulation configurations, provide a fixed-area carbon-model scaling table, and summarize mapped geologic units within the East River study domain, including geologic code, formation name, lithologic description, mapped area, and lithologic class grouping. Model outputs and geologic summaries support interpretation of transient landscape behavior and its relation to the mapped distribution of shale, intrusive, sedimentary, and surficial units. A third group of tables reports hydrologic and hydrochemical information used to quantify dissolved export from the watershed. Included files provide site-level values for drainage area, mean annual solute export, standard error of annual export, area-normalized solute yield, and equivalent weathering rate for five East River monitoring sites, along with metadata describing the number, sampling cadence, and date range of discharge records and partial and full total dissolved solids observations used in the solute-yield analyses. The package also contains a supplementary daily ion-load time series with daily mean discharge, discharge observation counts, dissolved concentrations, and daily loads for calcium, magnesium, sodium, potassium, chloride, sulfate, nitrate, fluoride, dissolved silica, charge-balance bicarbonate, and total dissolved solids. The package contains GIS files, comma-separated value files (.csv), CSV data dictionaries, a file-level metadata table, a package-tree text file, and a readme text file.

10Be↗
DOE OSTI2020

Modeling Molten Salt Reactor Fission Product Removal with SCALE

Liquid-fuel Molten Salt Reactors (MSRs) allow for continuous online fuel treatment and processing using a variety of subsystems. Among these subsystems, the off-gas system (OGS) continuously removes fission products from the primary fuel salt through a helium sparging process. Accurately modeling fission product behavior is paramount for the development of MSRs; while there are a variety of depletion modeling and simulation tools that allow for removal of material from a defined system, most of them work in a semicontinuous, batch-wise manner. The new tools implemented into TRITON, in the SCALE suite for reactor analysis and design, aim to address this issue by implementing truly continuous material feeds and removals. In addition, material accountancy is enhanced in these tools, which enable tracking of removed material in user-specified mixtures. These continuous depletion tools provide a way to model material transport throughout systems more efficiently and accurately. This report validates the new SCALE MSR modeling capabilities using data from the Molten Salt Reactor Experiment (MSRE). After a framework was established for deriving realistic removal rates based on design specifications, these rates were used by TRITON to continuously remove specific elements while performing burnup calculations. Results were compared to those obtained in MSRE, with a focus on the 135 Xe poison fraction. After these capabilities were validated, the learned concepts and best practices were applied to model the Molten Salt Breeder Reactor (MSBR) and the Molten Salt Demonstration Reactor (MSDR) to study the effects of continuous fission product removal in fuel salt reactivity, as well as the lifetime effects.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

An Overview of Gemma FY20 Verification Activities

Gemma verification activities for FY20 can be divided into three categories: the development of specialized quadrature rules, initial progress towards the development of manufactured solutions for code verification, and automated code-verification testing. In the method-of-moments implementation of the electric-field integral equation, the presence of a Green’s function in the four-dimensional integrals yields singularities in the integrand when two elements are nearby. To address these challenges, we have developed quadrature rules to integrate the functions through which the singularities can be characterized. Code verification is necessary to develop confidence in the implementation of the numerical methods in Gemma. Therefore, we have begun investigating the use of manufactured solutions to more thoroughly verify Gemma. Manufactured solutions provide greater flexibility for testing aspects of the code; however, the aforementioned singularities provide challenges, and existing work is limited in rigor and quantity. Finally, we have implemented automated code-verification testing using the VVTest framework to automate the mesh refinement and execution of a Gemma simulation to generate mesh convergence data. This infrastructure computes the observed order of accuracy from these data and compares it with the theoretical order of accuracy to either develop confidence in the implementation of the numerical methods or detect coding errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

The Kinematically Hot, Extremely Metal-poor C-19 Stellar Stream in DESI DR2

Stellar streams are the result of a host galaxy’s gravitational potential tidally disrupting satellite dwarf galaxies and globular clusters (GCs), causing them to grow leading and trailing tidal tails. The C-19 stellar stream is an extremely metal-poor stellar population, showing chemical abundance patterns characteristic of a GC. However, its large velocity dispersion is difficult to reconcile with a conventional, purely baryonic, disrupting-GC progenitor. Current techniques for stream characterization are primarily applied to Gaia DR3, relying heavily on proper-motion measurements. Using the Dark Energy Spectroscopic Instrument (DESI), which provides radial velocities and metallicities for over 10 million stars and reaches significantly fainter magnitudes than comparable surveys, we employ a mixture model approach to jointly characterize stream populations in proper motions, radial velocities, and metallicities against a Milky Way halo background. By applying this framework to the C-19 stellar stream, we identify a total of 47 spectroscopically confirmed member stars, of which 41 are newly identified and only 6 were previously reported in the literature. In this work, we measure a velocity dispersion of $7.8^{+1.5}_{-1.3}$km s −1 and a mean metallicity of [Fe/H] = $-3.36^{+0.12}_{-0.10}$. We further identify a novel “spur” feature within the stream. We conclude that our measurements are in line with previous works identifying C-19 as a “hot,” metal-poor stream. In forthcoming work, we will apply this approach to many more streams in the DESI footprint, enabling population-level comparisons with predictions from simulations.

Mohammed, Nasser [University of Toronto; Dunlap In↗
DOE OSTI2025

An Evaluation of Representation Learning Methods in Particle Physics Foundation Models

We present a systematic evaluation of representation learning objectives for particle physics within a unified framework. Our study employs a shared transformer-based particle-cloud encoder with standardized preprocessing, matched sampling, and a consistent evaluation protocol on a jet classification dataset. We compare contrastive (supervised and self-supervised), masked particle modeling, and generative reconstruction objectives under a common training regimen. In addition, we introduce targeted supervised architectural modifications that achieve state-of-the-art performance on benchmark evaluations. This controlled comparison isolates the contributions of the learning objective, highlights their respective strengths and limitations, and provides reproducible baselines. We position this work as a reference point for the future development of foundation models in particle physics, enabling more transparent and robust progress across the community.

Chen, Michael [Caltech]↗
DOE OSTI2026

The Kinematically Hot, Extremely Metal-Poor C-19 Stellar Stream in DESI DR2

Stellar streams are the result of a host galaxy's gravitational potential tidally disrupting satellite dwarf galaxies and globular clusters (GCs), causing them to grow leading and trailing tidal tails. The C-19 stellar stream is an extremely metal-poor stellar population, showing chemical abundance patterns characteristic of a globular cluster. However, its large velocity dispersion is difficult to reconcile with a conventional, purely baryonic, disrupting-GC progenitor. Current techniques for stream characterization are primarily applied to Gaia DR3, relying heavily on proper motion measurements. Using the Dark Energy Spectroscopic Instrument (DESI), which provides radial velocities and metallicites for over 10 million stars reaching significantly fainter magnitudes than comparable surveys, we employ a mixture model approach to jointly characterize stream populations in proper motions, radial velocities, and metallicities against a Milky Way halo background. By applying this framework to the C-19 stellar stream, we identify a total of 47 spectroscopically confirmed member stars, of which 41 are newly identified and only 6 were previously reported in the literature. In this work, we measure a velocity dispersion of $7.8^{+1.5}_{-1.3}$ km s$^{-1}$ and a mean metallicity of [Fe/H] = $-3.36^{+0.12}_{-0.10}$. We further identify a novel 'spur' feature within the stream. We conclude that our measurements are in line with previous works identifying C-19 as a 'hot', metal-poor stream. In forthcoming work, we will apply this approach to many more streams in the DESI footprint, enabling population-level comparisons with predictions from simulations.

Mohammed, Nasser [Toronto U.] (ORCID:0009000812240↗
DOE OSTI2026

The Kinematically Hot, Extremely Metal-Poor C-19 Stellar Stream in DESI DR2

Stellar streams are the result of a host galaxy's gravitational potential tidally disrupting satellite dwarf galaxies and globular clusters (GCs), causing them to grow leading and trailing tidal tails. The C-19 stellar stream is an extremely metal-poor stellar population, showing chemical abundance patterns characteristic of a globular cluster. However, its large velocity dispersion is difficult to reconcile with a conventional, purely baryonic, disrupting-GC progenitor. Current techniques for stream characterization are primarily applied to Gaia DR3, relying heavily on proper motion measurements. Using the Dark Energy Spectroscopic Instrument (DESI), which provides radial velocities and metallicites for over 10 million stars reaching significantly fainter magnitudes than comparable surveys, we employ a mixture model approach to jointly characterize stream populations in proper motions, radial velocities, and metallicities against a Milky Way halo background. By applying this framework to the C-19 stellar stream, we identify a total of 47 spectroscopically confirmed member stars, of which 41 are newly identified and only 6 were previously reported in the literature. In this work, we measure a velocity dispersion of $7.8^{+1.5}_{-1.3}$ km s$^{-1}$ and a mean metallicity of [Fe/H] = $-3.36^{+0.12}_{-0.10}$. We further identify a novel 'spur' feature within the stream. We conclude that our measurements are in line with previous works identifying C-19 as a 'hot', metal-poor stream. In forthcoming work, we will apply this approach to many more streams in the DESI footprint, enabling population-level comparisons with predictions from simulations.

Mohammed, Nasser [Toronto U.] (ORCID:0009000812240↗
DOE OSTI2022

Electrocatalytic water oxidation performance in an extended porous organic framework with a covalent alliance of distinct Ru sites

The rational synthesis of durable, earth-abundant efficient electrocatalysts for the oxygen evolution reaction (OER) from water is one of the most important routes for storing renewable energy and minimizing fossil fuel combustion. The prime hurdles for effectively utilizing commercial RuO 2 as (OER) electrocatalysts are its very low stability, catalyst deactivation, and high cost. In this work, we explored a Ru-integrated porous organic polymer (Ru@Bpy-POP) by a facile one-pot Friedel–Crafts alkylation strategy between redox-active (Ru(demob) 3 Cl 2 ) and a carbazole unit, which is composed of unique features including an extended framework unit, isolated active sites, and tunable electrode kinetics. Ru@Bpy-POP can serve as a bridge between a Metal–Organic Framework (MOF) and POP-based catalytic systems with a balanced combination of covalent bonds (structural stability) and open metal sites (single site catalysis). Ru@Bpy-POP, deposited on a three-dimensional nickel foam electrode support, exhibits a promising electrocatalytic OER activity with an ultra-low ruthenium loading compared to a benchmark RuO 2 catalyst, providing an overpotential of about 270 mV to reach 10 mA cm –2 in an alkaline medium. Moreover, a high current density of 248 mA cm –2 was achieved for the Ru@Bpy-POP catalyst at only 1.6 V (vs. RHE), which is much higher than 91 mA cm –2 for commercial RuO 2 . The robust, albeit highly conjugated, POP framework not only triggered facile electro-kinetics but also suppressed aggregation and metallic corrosion during electrolysis. In particular, the benefits of covalent integration of distinct Ru sites into the framework can modulate intermediate adsorption and charge density, which contributes to its exceptional OER activity. All of the critical steps involved in OER are complemented by Density Functional Theory (DFT) calculations, which suggest that electrocatalytic water oxidation proceeds from a closed-shell configuration to open-shell electronic configurations with high-spin states. Finally, these open-shell configurations are more stable than their closed-shell counterparts by 1 eV, improving the overall catalytic activity.

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

Van der Waals density functional study of hydrocarbon adsorption and separation in metal–organic frameworks without open metal sites

Metal–organic frameworks (MOFs) have received significant attention thanks to their promising features in the storage and separation of guest molecules. MOFs without open metal sites are emerging as they are often less susceptible to poisoning compared to those with open metal sites. However, a complete understanding of the binding and gas separation mechanisms in such materials is still missing. In this work, we perform a comparative study of two classes of vanadium-based MOFs without open metal sites: MFM-300-V (III) and MFM-300-V (IV) , as well as MIL-47-V (III) and MIL-47-V (IV) . We employ first-principles van der Waals density functional theory to find the optimal binding conformations and binding energies of a series of small hydrocarbons within the pores of the aforementioned MOFs. Our study provides insight into the host–guest interactions in such MOFs without open metal sites, especially the role played by the bridging hydroxyl group (μ 2 –OH). We conclude that the bridging –OH group acts as a pseudo open metal site in these MOFs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mitigating contamination in LSS surveys: a comparison of methods

ABSTRACT Future large-scale structure surveys will measure the locations and shapes of billions of galaxies. The precision of such catalogues will require meticulous treatment of systematic contamination of the observed fields. We compare several existing methods for removing such systematics from galaxy clustering measurements. We show how all the methods, including the popular pseudo-Cℓ Mode Projection and Template Subtraction methods, can be interpreted under a common regression framework and use this to suggest improved estimators. We show how methods designed to mitigate systematics in the power spectrum can be used to produce clean maps, which are necessary for cosmological analyses beyond the power spectrum, and we extend current methods to treat the next-order multiplicative contamination in observed maps and power spectra, which reduced power spectrum errors from $\Delta \chi ^2_{\rm C_\ell }\simeq 10$ to ≃ 1 in simulated analyses. Two new mitigation methods are proposed, which incorporate desirable features of current state-of-the-art methods while being simpler to implement. Investigating the performance of all the methods on a common set of simulated measurements from Year 5 of the Dark Energy Survey, we test their robustness to various analysis cases. Our proposed methods produce improved maps and power spectra when compared to current methods, while requiring almost no user tuning. We end with recommendations for systematics mitigation in future surveys, and note that the methods presented are generally applicable beyond the galaxy distribution to any field with spatial systematics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Determination of nuclear PDFs using Markov chain Monte Carlo methods

Global QCD analyses of nuclear parton distribution functions (nPDFs) have traditionally relied on the Hessian method for uncertainty estimation. However, the inherent Gaussian approximation and reliance on local curvature often prove insufficient for nPDF fits, which are frequently characterized by limited data constraints and non-Gaussian likelihoods. In this paper, we present the first nPDF determination based on Markov Chain Monte Carlo (MCMC) techniques, implemented within the nCTEQ framework using an adaptive Metropolis-Hastings algorithm. The MCMC approach enables a direct mapping of the posterior distribution and reveals a highly nontrivial parameter-space structure, including multiple modes and pronounced non-Gaussian behavior, particularly for the valence PDFs. We perform the first single-nucleus global analysis of lead PDFs using exclusively lead data and compare it to a multi-nuclei fit employing a standard analytic A dependence. The inclusion of lighter nuclei reduces quark uncertainties and modifies the shape of the lead PDFs, while leaving the gluon distribution largely unaffected. A complementary Hessian analysis exposes systematic limitations of the Gaussian approximation. Our results demonstrate that MCMC methods provide a more reliable framework for uncertainty quantification in nPDF determinations.

Derakhshanian, N. [Institute of Nuclear Physics Po↗
DOE OSTI2025

On the validity and limitations of 1D model for heat and mass transfer performance evaluation in a multilayer binder-free desiccant dehumidifier: isothermal dehumidification with internal cooling

Efficient humidity control is essential for maintaining indoor thermal comfort, yet conventional vapor-compression-based dehumidifiers are energy-intensive. Employing separate sensible and latent cooling through desiccant-coated heat exchangers (DCHEs) combined with evaporative coolers offers energy savings of up to 80 % compared to conventional systems. However, the dehumidification performance of DCHEs remains limited due to the use of polymer binders for coating desiccant materials onto heat exchange surfaces. In our previous study, we developed a multilayer fixed-bed binder-free desiccant dehumidifier (MFBDD) that demonstrated high dehumidification capacity and low pressure drop compared to rotary desiccant wheels. Nevertheless, its potential for further enhancement through internal cooling and the use of step-shaped adsorption isotherms has not been explored. In this study, a physics-based one-dimensional (1D) transient model is developed and validated to capture the coupled heat and mass transfer processes in the MFBDD and extended to simulate internal cooling using a high-capacity composite metal–organic framework, MIL-101/GO-6 (water uptake ≈1.6 g/g within 35–47 % RH). The model enables detailed analysis of local air and bed temperature dynamics and quantifies how internal cooling affects the dehumidification performance under a wide range of operating conditions. Results show that integrating internal cooling and using MIL-101/GO-6 enhance mass adsorbed, moisture removal capacity, and dehumidification effectiveness by 50 %–99 % compared with the M.S. Gel baseline. The study further reveals that achieving near-isothermal operation requires simultaneous enhancement of the convective heat transfer coefficient and heat exchange surface area. In conclusion, this work provides the first detailed physical insight into the interplay between internal cooling and step-shaped isotherms in a binder-free desiccant device and establishes a validated modeling framework for scaling up and system-level performance evaluation of next-generation energy-efficient dehumidification systems.

Heat and mass transfer↗
DOE OSTI2023

Automated Generation of Microkinetics for Heterogeneously Catalyzed Reactions Considering Correlated Uncertainties**

Abstract The study presents an ab‐initio based framework for the automated construction of microkinetic mechanisms considering correlated uncertainties in all energetic parameters and estimation routines. 2000 unique microkinetic models were generated within the uncertainty space of the BEEF‐vdW functional for the oxidation reactions of representative exhaust gas emissions from stoichiometric combustion engines over Pt(111) and compared to experiments through multiscale modeling. The ensemble of simulations stresses the importance of considering uncertainties. Within this set of first‐principles‐based models, it is possible to identify a microkinetic mechanism that agrees with experimental data. This mechanism can be traced back to a single exchange‐correlation functional, and it suggests that Pt(111) could be the active site for the oxidation of light hydrocarbons. The study provides a universal framework for the automated construction of reaction mechanisms with correlated uncertainty quantification, enabling a DFT‐constrained microkinetic model optimization for other heterogeneously catalyzed systems.

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

MIXv2: a long-term mosaic emission inventory for Asia (2010–2017)

The MIXv2 Asian emission inventory is developed under the framework of the Model Inter-Comparison Study for Asia (MICS-Asia) Phase IV and produced from a mosaic of up-to-date regional emission inventories. We estimated the emissions for anthropogenic and biomass burning sources covering 23 countries and regions in East, Southeast and South Asia and aggregated emissions to a uniform spatial and temporal resolution for seven sectors: power, industry, residential, transportation, agriculture, open biomass burning and shipping. Compared to MIXv1, we extended the dataset to 2010–2017, included emissions of open biomass burning and shipping, and provided model-ready emissions of SAPRC99, SAPRC07, and CB05. A series of unit-based point source information was incorporated covering power plants in China and India. A consistent speciation framework for non-methane volatile organic compounds (NMVOCs) was applied to develop emissions by three chemical mechanisms. The total Asian emissions for anthropogenic/open biomass sectors in 2017 are estimated as follows: 41.6/1.1 Tg NO x , 33.2/0.1 Tg SO 2 , 258.2/20.6 Tg CO, 61.8/8.2 Tg NMVOC, 28.3/0.3 Tg NH 3 , 24.0/2.6 Tg PM 10 , 16.7/2.0 Tg PM 2.5 , 2.7/0.1 Tg BC (black carbon), 5.3/0.9 Tg OC (organic carbon), and 18.0/0.4 Pg CO 2 . The contributions of India and Southeast Asia were emerging in Asia during 2010–2017, especially for SO 2 , NH 3 and particulate matter. Gridded emissions at a spatial resolution of 0.1° with monthly variations are now publicly available. This updated long-term emission mosaic inventory is ready to facilitate air quality and climate model simulations, as well as policymaking and associated analyses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Exploring the fusion power plant design space: comparative analysis of positive and negative triangularity tokamaks through optimization

The optimal configuration choice between positive triangularity (PT) and negative triangularity (NT) tokamaks for fusion power plants hinges on navigating different operational constraints rather than achieving specific plasma performance metrics. This study presents a systematic comparison using constrained multi-objective optimization with the integrated FUsion Synthesis Engine (FUSE) framework. Over 200 000 integrated design evaluations were performed exploring the trade-offs between capital cost minimization and operational reliability (maximizing q 95 ) while satisfying engineering constraints including 250 ± 50 MW net electric power, tritium breeding ratio > 1.1, power exhaust limits and an hour flattop time. Both configurations achieve similar cost-performance Pareto fronts through contrasting design philosophies. PT, while demonstrating resilience to pedestal degradation (compensating for up to 40% reduction), are constrained to larger machines (R 0 > 6.5 m) by the narrow operational window between L–H threshold requirements and the research-established power exhaust limit (P sol /R < 15 MW m –1 ). This forces optimization through comparatively reduced magnetic field (∼8 T). NT configurations exploit their freedom from these constraints to access compact, high-field designs (R 0 ~ 5.5 m, B 0 >12 T), creating natural synergy with advancing HTS technology. Sensitivity analyses reveal that PT’s economic viability depends critically on uncertainties in L–H threshold scaling and power handling limits. Notably, a 50% variation in either could eliminate viable designs or enable access to the compact design space. These results suggest configuration selection should be risk-informed: PT offers the lowest-cost path when operational constraints can be confidently predicted, while NT is robust to large variations in constraints and physics uncertainties.

FUSE framework↗
DOE OSTI2025

Quantum annealing-assisted lattice optimization

High Entropy Alloys (HEAs) have drawn great interest due to their exceptional properties compared to conventional materials. The configuration of HEA system is considered a key to their superior properties, but exhausting all possible configurations of atom coordinates and species to find the ground energy state is extremely challenging. In this work, we proposed a quantum annealing-assisted lattice optimization (QALO) algorithm, which is an active learning framework that integrates the Field-aware Factorization Machine (FFM) as the surrogate model for lattice energy prediction, Quantum Annealing (QA) as an optimizer and Machine Learning Potential (MLP) for ground truth energy calculation. By applying our algorithm to the NbMoTaW alloy, we reproduced the Nb depletion and W enrichment observed in bulk HEA. We found our optimized HEAs to have superior mechanical properties compared to the randomly generated alloy configurations. Our algorithm highlights the potential of quantum computing in materials design and discovery, laying a foundation for further exploring and optimizing structure-property relationships.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A constitutive model for elastomers tailored by ionic bonds and entanglements

Over the past decade or two, the concept has emerged of using multiple types of weak interactions simultaneously to enhance the mechanical properties of elastomers. These weak interactions include physical entanglements, hydrogen bonds, metal-coordination bonds, dynamic covalent bonds, and ionic bonds. The combination of entanglements and ionic bonding has been minimally explored and is particularly exciting because of the broad application space for polyelectrolytes. In this work, a constitutive model framework is developed to describe the response of elastomers with both ionic bonds and entanglements. Here we formulate a micromechanical model that couples together chain stretching, ionic bond slipping, and entanglement evolution. The ionic bonds provide toughness by enabling plastic deformation in comparison to covalently crosslinked material and add strength compared to a linear polymer. Evolution of the entanglement density is taken as a key mechanism that can govern stiffness, toughness, and self-recovery in elastomers. The model is used to match bulk polyelectrolytes with different fractions of ionic components under a variety of loading histories. The variations in material parameters are then used to help understand the relative importance of different governing mechanisms in the bulk polymers. We show that the theoretical framework can explain our experimental uniaxial tensile experimental results for polyelectrolytes. This model can help to design better material with high stiffness and toughness. We expect that our model can be extended to explain the mechanical behavior of other polyelectrolytes and other soft materials with a wide range of dynamic bonds.

42 ENGINEERING↗