SEARCH · Engineering Papers

Results for “equilibrium modeling”

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 829 records · Page 46

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Light Curves and Event Rates of Axion Instability Supernovae

Abstract It was recently proposed that exotic particles can trigger a new stellar instability that is analogous to the e − e + pair instability if they are produced and reach equilibrium in the stellar plasma. In this study, we construct axion instability supernova (AISN) models caused by the new instability to predict their observational signatures. We focus on heavy axion-like particles (ALPs) with masses of ∼400 keV–2 MeV and coupling with photons of g a γ ∼ 10 −5 GeV −1 . It is found that the 56 Ni mass and the explosion energy are significantly increased by ALPs for a fixed stellar mass. As a result, the peak times of the light curves of AISNe occur earlier than those of standard pair-instability supernovae by 10–20 days when the ALP mass is equal to the electron mass. Also, the event rate of AISNe is 1.7–2.6 times higher than that of pair-instability supernovae, depending on the high mass cutoff of the initial mass function.

Mori, Kanji (ORCID:0000000325951657)↗

DOE OSTI2021

SYK wormhole formation in real time

We study the real time formation of the ground state of two coupled SYK models. This is a highly entangled state which is close to the thermofield double state and can be viewed as a wormhole. We start from a high temperature state, we let it cool by coupling to a cold bath. We numerically solve for the large N dynamics. Our main result is that the system forms a wormhole by going through a region with negative specific heat, taking time that is independent of N. The dynamics is smooth everywhere and it seems to follow the equilibrium thermodynamic configurations of the microcanonical ensemble. Also we comment on the relation between this coupled SYK model and Jackiw-Teitelboim gravity theory with bulk fields.

1/N expansion↗

DOE OSTI2020

The effects of introducing elasticity using different interpolation schemes to the grand potential phase field model

Introducing elastic energy in the phase field method has been shown to influence interfacial energy, depending on the elastic interpolation scheme. This study investigates the impact of the elastic energy when using a grand potential-based phase field method, comparing the result of Khachaturyan’s strain interpolation scheme (KHS) and Voight-Taylor’s elastic energy interpolation scheme (VTS). The KHS model leads to a decrease in the interfacial energy, while the VTS model leads to an increase. The change in interfacial energy is greater with the VTS model than the KHS model, which suggests that the KHS model is more appropriate to limit the artificial impact of the elastic energy on the interfacial energy. When the contribution at the interface is not negligible, it is shown that both the microstructure evolution kinetics and the equilibrium microstructure can be influenced by the choice of the elastic scheme being used. In addition, this paper shows that the grand potential model might not be appropriate when the system requires the introduction of a composition-dependent term in the elastic energy contribution. This limitation is due to the need for an explicit and invertible relation between the total potential and the composition.

36 MATERIALS SCIENCE↗

DOE OSTI2024

K ∗ (892) ± resonance production in Pb−Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The production of K *⁢(892) ± meson resonance is measured at midrapidity (| y |<0.5) in Pb-Pb collisions at $\sqrt{^sNN}$ = 5.02 TeV using the ALICE detector at the CERN Large Hadron Collider. The resonance is reconstructed via its hadronic decay channel K *⁢(892) ± → K $^0_S$⁢π ± . The transverse momentum distributions are obtained for various centrality intervals in the p T range of 0.4-16 GeV/ c . Measurements of integrated yields, mean transverse momenta, and particle yield ratios are reported and found to be consistent with previous ALICE measurements for K *⁢(892) 0 within uncertainties. The p T -integrated yield ratio 2 K *⁢(892) ± ⁢/(⁢ K + + K - ) in central Pb-Pb collisions shows a significant suppression at a level of 9.3σ⁢ relative to pp collisions. Thermal model calculations result in an overprediction of the particle yield ratio. Although both hadron resonance gas in partial chemical equilibrium (HRG-PCE) and MUSIC + SMASH simulations consider the hadronic phase, only HRG-PCE accurately represents the measurements, whereas MUSIC + SMASH simulations tend to overpredict the particle yield ratio. These observations, along with the kinetic freeze-out temperatures extracted from the yields measured for light-flavored hadrons using the HRG-PCE model, indicate a finite hadronic phase lifetime, which decreases with increasing collision centrality percentile. The p T-differential yield ratios 2 K *⁢(892) ±⁢ /(⁢ K + + K - ) and 2 K *⁢(892) ± ⁢/(⁢π + +π - ) are presented and compared with measurements in pp collisions at $\sqrt{s}$ = 5.02 TeV. Both particle ratios are found to be suppressed by up to a factor of five at p T <2.0 GeV/ c in central Pb-Pb collisions and are qualitatively consistent with expectations for rescattering effects in the hadronic phase. The nuclear modification factor ( R AA ) shows a smooth evolution with centrality and is found to be below unity at p T >8 GeV/ c , consistent with measurements for other light-flavored hadrons. The smallest values are observed in most central collisions, indicating larger energy loss of partons traversing the dense medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Quarkyonic matter equation of state in beta-equilibrium

Quark matter may appear due to a hadronic-quark transition in the core of a hybrid star. Quarkyonic matter is an approach in which both quarks and nucleons appear as quasiparticles in a crossover transition, and provides an explicit realization of early ideas concerning quark matter (e.g., the MIT bag model). This description has recently been employed by McLerran and Reddy to model chargeless (pure neutron) matter with an approach that has the virtue that the speed of sound rises quickly at a neutron-quark transition so as to satisfy observational constraints on the neutron star maximum mass (≳ 2 $M_⊙$) and the radius of a 1.4 $M_⊙$ star ($R_{1.4}$ ≲ 13.5 km). Traditional models involving first-order transitions result in softer pressure-energy density relations that have difficulty satisfying these constraints except with very narrow choices of parameters. Here, we propose a variation of quarkyonic matter involving protons and leptons whose energy can be explicitly minimized to achieve both chemical and beta equilibrium, which cannot be done in the chargeless formulation. Quarkyonic stellar models are able to satisfy observed mass and radius constraints with a wide range of model parameters, avoiding the obligatory fine-tuning of conventional hybrid star models, including requiring the transition density to be very close to the nuclear saturation density. Our formulation fits experimental and theoretical properties of the nuclear symmetry energy and pure neutron matter, and contains as few as three free parameters. This makes it an ideal tool for the study of high-density matter that is an efficient alternative to piecewise polytrope or spectral decomposition methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Adsorption Equilibrium, Kinetics, and Column Breakthrough Data for Aqueous Solutions of Binary-Acid and Ternary-Acid Mixtures of Acetic Acid, Butyric Acid, and Lactic Acid on IRN-78 Ion-Exchange Resin at Initial pH Levels of ∼3–7 and at 25–55 °C

This work is part of an effort to develop thermophysical property data and models supporting adsorptive process development for organic acid separation from a dilute aqueous solution of fermentation broth. It presents systematic experimental measurements for aqueous-phase adsorption equilibrium, kinetics, and column breakthrough for three binary-acid aqueous mixtures (acetic acid + lactic acid, butyric acid + lactic acid, and acetic acid + butyric acid) and one ternary-acid aqueous mixture (acetic acid + butyric acid + lactic acid) on Amberlite IRN-78 ion-exchange resin. The equilibrium measurements covered broad ranges of the initial acid concentration (100–400 mmol/L), initial pH (∼3–7), and temperature (25–55 °C). The equilibrium data for the binary-acid and ternary-acid aqueous mixtures indicate selective adsorption of lactic acid at initial pH levels of ∼3 and 4, while acetic acid and butyric acid are selectively adsorbed at initial pH levels of 5, 6, and 7. The subsequent kinetics and column breakthrough experiments were performed at 200 mmol/L with equimolar ratios, an initial pH of 6, and 25 °C. In conclusion, the measurements provide essential data sets for rigorous thermodynamic modeling and process simulation of adsorptive separation processes for organic acid separation from fermentation broth.

Adsorption↗

DOE OSTI2022

Nonequilibrium effects and transverse spherocity in ultrarelativistic proton-nucleus collisions

We investigate the effects of nonequilibrium dynamics in small colliding systems by comparing a nonequilibrium transport approach, Parton-Hadron-String Dynamics (PHSD), with a (2+1)-dimensional viscous hydrodynamic model, VISHNew. Focusing on p + Pb collisions at energies available at the Large Hadron Collider, we extract the initial conditions for the hydrodynamic model from PHSD, in order to reduce the impact of the early out-of-equilibrium dynamics and focus on the traces of nonequilibrium in the ensuing medium evolution. We find that in the transport approach quantities like energy density and bulk viscous pressure are highly inhomogeneous on the transverse plane during the whole evolution, whereas the hydrodynamic simulations dissolve more efficiently the initial spatial irregularities, still keeping a high degree of inhomogeneity due to the smaller size and lifetime of the medium produced in p + Pb collisions with respect to heavy-ion reactions. As a first step that will help to identify the impact of these nonequilibrium effects on final observables in proton-nucleus collisions, we perform an analysis of the transverse spherocity, an event-shape observable able to distinguish between jetty and isotropic configurations of transverse momenta. We found that the spherocity distribution in PHSD is slightly shifted towards the isotropic limit with respect to the hydrodynamic result. Even though this dissimilarity is partially due to the difference in the final charged particle production, it mainly comes from the different description within the two frameworks of the medium produced in small colliding systems. This finding supports the idea that multidifferential measurements, such as those based on event categorization according to multiplicity and spherocity, are useful to study final-state observables in ultrarelativistic proton-nucleus collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Prediction of transport in the JET DTE2 discharges with TGLF and NEO models using the TGYRO transport code

Abstract The JET Deuterium-Tritium-Experiment Campaign 2 (DTE2) has demonstrated the highest-ever fusion energy production. To forecast the transport dynamics within these discharges, the TGLF and NEO models within the TGYRO transport code were employed. A critical development in this study is the new quasilinear transport model, TGLF-SAT2, specifically designed to resolve discrepancies identified in JET deuterium discharges. This model accurately describes the saturated three-dimensional (3D) fluctuation spectrum, aligning closely with a database of nonlinear CGYRO turbulence simulations, thereby enhancing the predictive accuracy of TGYRO simulations. In validating against the JET DTE2 discharges across two primary operating scenarios, TGYRO effectively predicted the temperature profiles within a broad radial window ( ρ ∼ 0.2–0.85), though with minor ion temperature discrepancies near the core. However, a consistent underprediction of electron density profiles by 20% across the simulation domain was noted, indicating areas for future refinement. To achieve a self-consistent steady-state solution based on the JET DTE2 discharges, an integrated modeling workflow TGYRO-STEP within the OMFIT framework was introduced. This workflow iterates among the core transport, the pedestal pressure and the MHD equilibrium, ultimately yielding a converged solution that significantly reduces dependence on experimental boundary conditions for temperature and density profiles. The integrated simulation results show negligible differences in electron density and temperature profiles compared to standalone TGYRO modeling, while the ion temperature profile is lower due to the updated boundary condition in TGYRO-STEP. The application of the TGYRO-STEP workflow to JET DTE2 discharges serves as a crucial test to validate its robustness and highlights its limitations, providing valuable insights for its potential future application in ITER and Fusion Power Plant deuterium and tritium prediction modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Applying deep learning methods to develop new models of molecular charge transfer, nonadiabatic dynamics, and nonlinear spectroscopy in the condensed phase

Photon- and field-induced charge transfer has central importance in the generation and storage of electricity, the novel properties of materials, photo-induced catalysis, and electro-optic activity (e.g., photovoltaic cells, fuel cells, and organic chromophores for use in optical fibers and light-emission diodes). These non-equilibrium electronic and chemical transformations are probed by ultrafast, nonlinear spectroscopies. Accurate simulations play a crucial role in our ability to understand, optimize, and control these transformations. This project applies modern deep learning and machine learning (ML) methods to dramatically improve models of electronic dynamics, electronic-nuclear dynamics, and spectroscopic measurements for improved simulations of chemistry in complex environments, far from equilibrium phenomena, and processes in extreme environments, such as materials exposed to strong or resonant fields. This project develops accurate neural net models that go beyond predictive capability to also provide new insight into the fundamental physics underlying electron and nuclear dynamics. To achieve its objectives, this project explores and develops customized versions of high-capacity deep learning algorithms/models. These techniques are developed with an emphasis on fundamental chemical insight, not just predictive accuracy, to assist the development of the next generation of quantum simulation methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Using dynamic mode decomposition to predict the dynamics of a two-time non-equilibrium Green’s function

Computing the numerical solution of the Kadanoff–Baym equations, a set of nonlinear integral differential equations satisfied by the two-time Green's functions derived from many-body perturbation theory for a quantum many-body system away from equilibrium, is a challenging task. Recently, we have successfully applied dynamic mode decomposition (DMD) to construct a data driven reduced order model that can be used to extrapolate the time-diagonal of a two-time Green's function from numerical solutions of the KBE within a small time window. In this paper, we extend the previous work and use DMD to predict off-diagonal elements of the two-time Green's function. We partition the two-time Green's function into a number of one-time functions along the diagonal and subdiagonals of the two-time window as well as in horizontal and vertical directions. We use DMD to construct separate reduced order models to predict the dynamics of these one-time functions in a two-step procedure. We extrapolate along diagonal and several subdiagonals within a subdiagonal band of a two-time window in the first step. In the second step, we use DMD to extrapolate the Green's function outside of the sub-diagonal band. In conclusion, we demonstrate the efficiency and accuracy of this approach by applying it to a two-band Hubbard model problem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

ReaxFF reactive force field model enables accurate prediction of physiochemical and mechanical properties of crystalline and amorphous shape‐memory polyurethane

Abstract We investigate reactive force field (ReaxFF) prediction of physiochemical and mechanical properties, and stimuli response behavior of a shape‐memory polyurethane (SMPU) rigid segment. We used SMPUs as a platform to link the rigid/hard segment domain crystallinity with properties that can also be experimentally measured. Specifically, we developed models for crystalline and amorphous forms of a common rigid segment of SMPUs, for example, 4,4′‐diphenylmethane diisocyanate (MDI) with n‐butanediol (BDO), denoted as MDI‐BDO. ReaxFF molecular dynamics simulations with equilibrium dynamics are performed by controlling mass, temperature, and pressure or volume to predict the structural and physiochemical properties of crystalline and amorphous rigid segments of SMPU systems, followed by simulating their XRD and FTIR patterns, respectively. Non‐equilibrium ReaxFF simulations are implemented by using uniaxial box deformation technique to study response‐behavior of crystalline SMPU under tensile loading and during stress‐relaxation. The model and simulation scheme presented in this study demonstrates a robust theoretical pathway toward application‐dictated, on‐demand design of high‐performance SMPUs in applications such as soft robotics.

Afroz, Mohammad M.↗

DOE OSTI2020

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. In this study, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

EF-Hand Battle Royale: Hetero-ion Complexation in Lanmodulin

The lanmodulin (LanM) protein has emerged as an effective means for rare earth element (REE) extraction and separation from complex feedstocks without the use of organic solvents. Whereas the binding of LanM to individual REEs has been well characterized, little is known about the thermodynamics of mixed metal binding complexes (i.e., heterogeneous ion complexes), which limits the ability to accurately predict separation performance for a given metal ion mixture. In this paper, we employ the law of mass action to establish a theory of perfect cooperativity for LanM-REE complexation at the two highest-affinity binding sites. The theory is then used to derive an equation that explains the nonintuitive REE binding behavior of LanM, where separation factors for binary pairs of ions vary widely based on the ratio of ions in the aqueous phase, a phenomenon that is distinct from single-ion-binding chemical chelators. We then experimentally validate this theory and perform the first quantitative characterization of LanM complexation with heterogeneous ion pairs using resin-immobilized LanM. Importantly, the resulting homogeneous and heterogeneous constants enable accurate prediction of the equilibrium state of LanM in the presence of mixtures of up to 10 REEs, confirming that the perfect cooperativity model is an accurate mechanistic description of REE complexation by LanM. We further employ the model to simulate separation performance over a range of homogeneous and heterogeneous binding constants, revealing important insights into how mixed binding differentially impacts REE separations based on the relative positioning of the ion pairs within the lanthanide series. In addition to informing REE separation process optimization, these results provide mathematical and experimental insight into competition dynamics in other ubiquitous and medically relevant, cooperative binding proteins, such as calmodulin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Ringing Artifact Reduction Method for Ultrasound Reconstruction Using Multi-Agent Consensus Equilibrium

Non-destructive characterization of multi-layered structures that can be accessed from only a single side is important for applications such as well-bore integrity inspection. Existing methods related to Synthetic Aperture Focusing Technique (SAFT) rapidly produce acceptable results but with significant artifacts. Recently, ultrasound model-based iterative reconstruction (UMBIR) approaches have shown significant improvements over SAFT. However, even these methods produce ringing artifacts due to the high fractional-bandwidth of the excitation signal.In this paper, we propose a ringing artifact reduction method for ultrasound image reconstruction that uses a multi-agent consensus equilibrium (RARE-MACE) framework. Our approach integrates a physics-based forward model that accounts for the propagation of a collimated ultrasonic beam in multi-layered media, a spatially varying image prior, and a denoiser designed to suppress the ringing artifacts that are characteristic of reconstructions from high-fractional bandwidth ultrasound sensor data. We test our method on simulated and experimental measurements and show substantial improvements in image quality compared to SAFT and UMBIR.

Alanazi, Abdulrahman M↗

DOE OSTI2023

Revealing the dynamics of ultrarelativistic non-equilibrium many-electron systems with phase space tomography

The description of physical processes with many-particle systems is a key approach to the modeling of numerous physical systems. For example in storage rings, where ultrarelativistic particles are agglomerated in dense bunches, the modeling and measurement of their phase-space distribution is of paramount importance: at any time the phase-space distribution not only determines the complete space-time evolution but also provides fundamental performance characteristics for storage ring operation. Here, we demonstrate a non-destructive tomographic imaging technique for the 2D longitudinal phase-space distribution of ultrarelativistic electron bunches. For this purpose, we utilize a unique setup, which streams turn-by-turn near-field measurements of bunch profiles at MHz repetition rates. To demonstrate the feasibility of our method, we induce a non-equilibrium state and show that the phase-space distribution microstructuring as well as the phase-space distribution dynamics can be observed in great detail. Our approach offers a pathway to control ultrashort bunches and supports, as one example, the development of compact accelerators with low energy footprints.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Preliminary Result of Neoclassical Impurity Transport Affected by ICRH with a 4-D Fokker-Planck Code

To understand the role of ion cyclotron resonance heating (ICRH) in mitigating tungsten accumulation, we evaluate the impact of ICRH-induced poloidal asymmetries on neoclassical impurity transport using a newly developed 4-D Fokker–Planck code (FP4D). Unlike conventional approaches that separate RF heating and transport modelling, FP4D solves both fast-ion distribution formation and impurity transport self-consistently without bounce-averaging. Benchmarking against the NEO code confirms the accuracy of the neoclassical physics in FP4D. Using FP4D, we analyze anisotropic and non-Maxwellian minority distributions under ICRH and their effects on the equilibrium potential and tungsten transport. Results show that anisotropy induced poloidal potential variations can reduce inward tungsten flux, depending on the modelling approach. Our findings highlight the importance of resolving the full velocity and poloidal structure of minority species in predicting impurity transport under strong RF heating conditions.

Jeong, Yunho [Hanyang Univ., Seoul (Korea, Republi↗

DOE OSTI2024

Uncovering grain and subgrain microstructure at the scale of additive manufacturing melt tracks with a scalable cellular automaton solidification model

Metal additive manufacturing, characterized by rapid solidification, yields refined grains with a distinctive cellular subgrain microstructure that plays a pivotal role in determining material properties. Due to the significant computational expense demanded to simulate the required physics with submicron spatial resolution, their numerical simulations have been limited to proof-of-concept studies to either 2D or small subregions of a melt pool. In this study, an open-source, scalable, solidification code, muMatScale, based on the cellular automaton method, has been developed to predict the grain and the underlying subgrain microstructure over an entire melt pool. The model incorporates flexible parallelization schemes, utilizing MPI and OpenMP GPU Offloading, in addition to appropriate multi-physics specific to non-equilibrium rapid solidification in AM. The impact of nucleation parameters on grain microstructures was investigated with a focus on grain size variations and morphology transitions. With selected nucleation parameters, the simulation predicted the grain size, subgrain morphology, crystallographic orientation, and microsegregation aligned with experimental measurements. The model demonstrates that epitaxial grain growth is a dominant factor at the melt pool boundary, influencing grain size variation under different grain sizes in the build plate while maintaining consistent primary dendrite arm spacing under identical thermal conditions. Here, the highly efficient numerical model enables large-scale simulations with a spatial resolution of 100 nm or less, unveiling unprecedented insights into thermal and solutal diffusion driven grain growth, and the subgrains with microsegregation within grains in 3D across scales. muMatScale will enable the linking of submicron length-scale microstructure to part-level material behavior by investigating fundamental solidification problems at the intercellular scale in many-track and many-layer builds.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A reactive electrochemomechanical theory for growth and remodeling of polyelectrolyte hydrogels and application to dynamic polymerization of $\text{DNA}$ hydrogels

Here this study develops a framework for growth and remodeling of active polyelectrolyte hydrogels that accounts for effects of compositional changes on the mechanical response. By developing a reactive electrochemomechanical theory, thermodynamical constraints upon reactive and remodeling processes are elucidated within a general framework that allows any number of chemical reactions to evolve the response of the gel and transfer mass and charge between constituents. Fully coupled, nonlinear constitutive relations are adopted for molar fluxes, allowing exploration of effects including cross-diffusion, electrophoresis, and electro-osmosis. A robust finite element implementation is developed in the open source FEBio software (febio.org) by exploiting an equivalence between electrochemomechanics and mixture theory. The implementation is verified against analytical solutions for free swelling, and a proper reduction to a prior chemomechanical theory is demonstrated for neutral gels swollen only by a solvent with no solutes. The theory and implementation are then applied to model the tunable large swelling achieved through dynamic polymerization of DNA crosslinkers seen in our recently developed experimental hydrogel system (Cangialosi et al., 2017). A novel constitutive model for reaction-driven evolution of the locking stretch λ L in a non-Gaussian mechanical free energy was developed, where the increasing concentration of DNA crosslinkers makes further swelling energetically favorable. With a single free parameter, excellent agreement was found between measured and predicted equilibrium swelling ratios. This study demonstrated the ability to extend the electrochemomechanical framework to include chemical reactions and composition-aware constitutive models, and showed that development of reactive models allows simulation of complex dynamic polymerization phenomena not treated before. The theoretical frame here can be further expanded in scope to incorporate additional non-ideal and nonlinear phenomena.

42 ENGINEERING↗