Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bulk modulus”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

P-V-T equation of state of boron carbide

We report the P-V-T equation of state measurements of B 4 C to 50 GPa and approximately 2500 K in laser-heated diamond anvil cells. We obtain an ambient temperature, third-order Birch–Murnaghan fit to the P-V data that yields a bulk modulus K 0 of 221(2) GPa and derivative, (dK/dP) 0 of 3.3(1). These were used in fits with both a Mie–Grüneisen–Debye model and a temperature-dependent, Birch–Murnaghan equation of state that includes thermal pressure estimated by thermal expansion (α) and a temperature-dependent bulk modulus (dK 0 /dT). The ambient pressure thermal expansion coefficient (α 0 + α 1 T), Grüneisen γ(V) = γ 0 (V/V 0 ) q and volume-dependent Debye temperature, were used as input parameters for these fits and found to be sufficient to describe the data in the whole P-T range of this study.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of Cr and Al on Elastic Constants of FeCrAl Alloys Investigated by Molecular Dynamics Method

The FeCrAl alloy system is recognized as one of the candidate materials for accident-tolerant fuel (ATF) cladding in the nuclear power industry due to its high oxidation resistance under irradiation and high-temperature environments. The concentrations of Cr and Al have a significant effect on elastic properties of the FeCrAl alloy. In this work, elastic constants C 11 , C 12 , C 44 , bulk modulus and shear modulus of FeCrAl alloy were calculated with molecular dynamics methods. We explored compositions with 1–15 wt.% Cr and 1–5 wt.% Al at temperatures from 0 K to 750 K. The results show that the concentrations of Al and Cr have different effects on the elastic constants. When the concentration of Al was fixed, a decrease in bulk modulus and shear modulus with increasing Cr content was observed, consistent with previous experimental results. The dependence of elastic constants on temperature was also the same as in the experiments. Investigations into elastic properties of defect-containing alloys have shown that vacancies, voids, interstitials and Cr-rich precipitations have different effects on elastic properties of FeCrAl alloys. Investigations of elastic properties of defect-containing alloys have shown that vacancies, void, interstitials and Cr-rich precipitations have different effects on elastic properties of FeCrAl alloys. Therefore, the present results indicate that both the Cr and Al concentrations and radiation defects should be considered to develop and apply the FeCrAl alloy in ATF design.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Deep potential development of transition-metal-rich carbides

In this study, we developed DeepPot-SE type potentials according to the High-Entropy Strategy for Cr 23 C 6 , Fe 23 C 6 , W 23 C 6 , and Mo 23 C 6 systems. Predictive molecular dynamics are then conducted according to various elastic and mechanical properties. Training dataset accuracy is confirmed both numerically and visually. Elastic constants, Poisson Ratio, and bulk modulus are determined and compared to the literature with percent errors ranging between 0.45 and 13.33%. Ground-state lattice constants yield 0.111–1.375% percent error. Linear lattice constant thermal expansion trends are found from 300 to 1500 K as expected. Melting behaviors with clear melting points are observed for each binary. Here, the Birch–Murnaghan Equation of State is suitably fit, providing additional verification of the bulk modulus calculations. Polycrystalline thermal stability is also verified for each binary potential.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Softening of breathing elastic mode and trigonal elastic mode in the disordered pyrochlore magnet NaCaCo 2 F 7

Cobalt pyrochlore fluoride NaCaCo 2 F 7 is a disordered frustrated magnet that is composed of Co 2+ ions with an effective spin-12 magnetic moment and that exhibits spin freezing below $T_f ~$ 2.4 K. Here, we perform ultrasound velocity measurements on a single crystal of the cubic NaCaCo 2 F 7 . The temperature dependence of the bulk modulus (the breathing elastic mode) exhibits Curie-type softening upon cooling below ~20 K down to $T_f$, which is suppressed by the magnetic field. This Curie-type softening should be a precursor to the enhancement of the strength of exchange disorder via the spin-lattice coupling, which causes the spin freezing. In contrast to the magnetic-field-suppressed Curie-type softening in the bulk modulus, the trigonal shear modulus exhibits softening with a characteristic minimum upon cooling, which is enhanced by the magnetic field at temperatures below ~20 K. This magnetic-field-enhanced elastic anomaly in the trigonal shear modulus suggests a coupling of the lattice to the dynamical spin-cluster state. For NaCaCo 2 F 7 , the observed elastic anomalies reveal an occurrence of magnetic-field-induced crossover from an isostructural lattice instability toward the spin freezing to a trigonal lattice instability arising from the emergent dynamical spin-cluster state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

X-Ray Diffraction and Electron Microscopy Studies of the Size Effects on Pressure-Induced Phase Transitions in CdS Nanocrystals

In recent years, investigations of the phase transition behavior of semiconducting nanoparticles under high pressure has attracted increasing attention due to their potential applications in sensors, electronics, and optics. However, current understanding of how the size of nanoparticles influences this pressure-dependent property is somewhat lacking. In particular, phase behaviors of semiconducting CdS nanoparticles under high pressure have not been extensively reported. Therefore, in this work, CdS nanoparticles of different sizes are used as a model system to investigate particle size effects on high-pressure-induced phase transition behaviors. In particular, 7.5, 10.6, and 39.7 nm spherical CdS nanoparticles are synthesized and subjected to controlled high pressures up to 15 GPa in a diamond anvil cell. Analysis of all three nanoparticles using in-situ synchrotron wide-angle X-ray scattering (WAXS) data shows that phase transitions from wurtzite to rocksalt occur at higher pressures than for bulk material. Bulk modulus calculations not only show that the wurtzite CdS nanomaterial is more compressible than rocksalt, but also that the compressibility of CdS nanoparticles depends on their particle size. Furthermore, sintering of spherical nanoparticles into nanorods was observed for the 7.5 nm CdS nanoparticles. Our results provide new insights into the fundamental properties of nanoparticles under high pressure that will inform designs of new nanomaterial structures for emerging applications.

Meng, Lingyao↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2020

First-principles study on structural mechanical and thermodynamic properties of HfMoTaTiZr

This study calculates the phase diagram of HfMoTaTiZr using the TCHEA1 database in Thermo-Calc software. Furthermore, it also investigates the structural and mechanical properties of HfMoTaTiZr by applying first-principles calculations. The elastic and other mechanical properties such as Young’s modulus, pugh ratio, hardness, bulk modulus was derived and discussed. These parameters were found to be in good harmony with previously available experimental reports. Vickers hardness of HfMoTaTiZr predicted from the Tians model agrees with the available experimental result. This study also investigates the thermal coefficient of linear expansion and vibrational heat capacity of the HfMoTaTiZr using the Debye–Grüneisen theory.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Analysis of differential scanning calorimetry data for aged plutonium

Differential scanning calorimetry data for samples of a 52 year old plutonium alloy with 3.3 at. % Ga that were heated beyond the melting point is analyzed using transition state theory to find activation energies for the δ to ε and ε to liquid phase transitions. A Bayesian statistical method involving a Gaussian process model is used to find mean values and confidence intervals for the activation energies. The activation energy for the δ to ε phase transition increases by 3.3 ± 3.8% per decade, relative to the case when all age related plutonium lattice point defects have been removed through annealing. The corresponding increase in activation energy for the ε to liquid transition is shown to be 7.1 ± 1.8% per decade. It is postulated that the change in activation energy with age for both phase transitions is caused, in part, by the accumulation of the same type of lattice point defects associated with the observed increase in elastic bulk modulus over time.

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

Mechanical and Thermal Properties for Uranium and U–6Nb Alloy from First-Principles Theory

Elasticity, lattice dynamics, and thermal expansion for uranium and U–6Nb alloy (elastic moduli) are calculated from density functional theory that is extended to include orbital polarization (DFT+OP). Introducing 12.5 at.% of niobium, substitutionally, in uranium softens all the cii elastic moduli, resulting in a significantly softer shear modulus (G). Combined with a nearly invariant bulk modulus (B), the quotient B/G increases dramatically for U–6Nb, suggesting a more ductile material. Lattice dynamics from a harmonic model coupled with a DFT+OP electronic structure is applied for α uranium, and the obtained phonon density of states compares well with inelastic neutron-scattering measurements. The Debye temperature associated with the lattice dynamics falls within the range of experimentally observed Debye temperatures and it also validates our quasi-harmonic (QH) phonon model. The QH Debye–Grüneisen phonon method is combined with a DFT+OP electronic structure and used to explore the anisotropic thermal expansion in α uranium. The anomalous negative thermal expansion (contraction) of the b lattice parameter of the α-phase orthorhombic cell is relatively well reproduced from a free-energy model consisting of QH-phonon and DFT+OP electronic structure contributions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Design and Development of Novel Equiatomic Refractory Multi Principal elemental Alloys Based on MoNbTi System for Use in Irradiation Environments

Multi Principal Elemental Alloys (MPEA) have emerged as promising materials for next-generation nuclear reactors due to their exceptional irradiation resistance. Eight equiatomic MPEA based on the MoNbTi system, comprising of elements with low thermal neutron absorption cross-sections were explored using a combined approach employing empirical parameter estimations, and CALPHAD simulations by which the phases and elemental segregation observed in all the alloys in the as-cast state were predicted. Solution heat treatment at 1500°C transformed five alloys into single-phase matrix materials, enhancing homogeneity and reducing hardness. The densities of the alloys ranged between 6.47 to 7.68 g/cm3, hardness between 472 and 656 VHN, Young’s modulus between 142 GPa to 169 GPa, shear modulus between 54 GPa and 62 GPa, bulk modulus between 117 GPa to 194 GPa and Poisson’s ratio between 0.3 to 0.35. The in-situ high temperature Xray diffraction results, differential scanning calorimetry and dilatometry results up to 1000°C suggested the high temperature phase stability of the MPEA. Subsequent ageing heat treatment at 800 and 1000 oC for 96 hours revealed significant secondary phase precipitation in MoNbTiZr, MoNbTiZrV, and MoNbTiCrA. Oxidation studies at 800 oC for 24 hours in air revealed superior oxidation resistance and cubic rate law dependence in Cr containing MPEA especially in MoNbTiCrAl, while severe mass gain resulting in total disintegration and exfoliation in MoNbTiZr and MoNbTiZrV. This comprehensive study underscores the potential of novel MPEA as promising materials for advanced nuclear reactor applications, shedding light on their microstructural control, mechanical properties, thermal stability, and oxidation resistance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

The first experimental evidence for improved nanomechanical properties of calcium silicate hydrate by polycarboxylate ether and graphene oxide

Graphene oxide (GO) efficiently enhances macroscale mechanical properties of cement-based materials. Yet the nanoscale reinforcing mechanism of GO for calcium silicate hydrate (C-S-H), the key binding phase of concrete, remains poorly understood. Moreover, how polycarboxylate ether (PCE)-based superplasticizers affect the nanomechanical properties of C-S-H remains unclear. The intrinsic mechanical properties of the nanocomposites are measured using high-pressure X-ray diffraction. The influences of GO and PCE on the local Ca and Si environments are probed using X-ray absorption spectroscopy. For the first time, we evidence that PCE-induced well-dispersed GO nanoplatelets strengthen the C-S-H basal planes, whereas poorly water-dispersed GO weakly interacts with basal planes and does not strengthen in plane Casingle bondO bonds. Both PCE and GO only slightly strengthen C-S-H's c-axis. Here, the bulk modulus of C-S-H/PCE/GO nanocomposite is 190% than pure C-S-H. Functional groups of PCE and GO prefer interacting with the Casingle bondO sheets of C-S-Hs instead of silicate tetrahedra.

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

Static and shock compression studies of eutectic high-entropy alloy AlCoCrFeNi 2.1 to ultrahigh pressures

The high-entropy alloy with composition AlCoCrFeNi 2.1 , additively manufactured with the laser powder-bed fusion technique, has a far-from-equilibrium BCC/FCC eutectic nanolamellar structure. We studied the high-pressure response of this alloy under both static compression and high-strain rate shock compression. The response to static compression using a diamond anvil cell was studied at pressures up to 302 GPa with synchrotron x-ray diffraction at the advanced photon source. The high-pressure FCC-only phase of the EHEA previously observed by Pope et al. [AIP Adv. 13, 035124 (2023)] is found to be stable up to the highest pressure achieved in this study with a volume compression of V/V 0 = 0.587 at ambient temperature. The shock experiments were performed by using GEKKO XII lasers at the Institute of Laser Engineering, Osaka University. The principal Hugoniot equation-of-state of the EHEA was measured up to a pressure of 515 GPa and a compression of V/V 0 = 0.613. Additionally, the thermal equation of state of the EHEA was measured up to 6.2 GPa and 1623 K using a large-volume Paris–Edinburgh cell to obtain the temperature dependence of bulk modulus and thermal expansion coefficients. The melting temperature for EHEA AlCoCrFeNi 2.1 at a pressure of 5.6 GPa was measured to be 1648 ± 25 K. These results can be used to refine stochastic (or special) quasi-random structure (SQS) models for high-pressure high-temperature behavior of high-entropy alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Shear strength measurements and hydrostatic compression of rhenium diboride under high pressures

Shear strength measurements have been carried out on rhenium diboride, ReB 2 , to a pressure of 74 GPa using a Radial X-ray Diffraction (R-XRD) technique in a diamond anvil cell using platinum as an internal x-ray pressure standard. The R-XRD result has provided a unique insight into the deformation of hexagonal crystal lattice under non-hydrostatic compression and variation of shear strength with increasing pressure. From R-XRD data, we have estimated hydrostatic component of compression to determine an equation of state of rhenium diboride yielding a bulk modulus of K 0 = 366 ± 25 GPa with a pressure derivative $K^{'}_{0}$ = 4.3 ± 0.5 in good agreement with hydrostatic density functional theory calculations. Here, the average lower bound of shear strength (τ) from various diffraction planes was then calculated using the measured interplanar d-spacing (d m ) and hydrostatic component of d-spacing (d p ) to be shown to approach 6.7 ± 0.4 GPa at 70 GPa. Our results show that the anisotropic compression effects observed in ReB 2 under hydrostatic compression are correlated to electronic structure changes under compression as predicted by theoretical calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermoelastic properties of tungsten at simultaneous high pressure and temperature

In this work, the compressional ( P ) and shear wave velocities ( S ) and unit cell volumes (densities) of polycrystalline tungsten (W) have been measured simultaneously up to 10.5 GPa and 1073 K using ultrasonic interferometry in conjunction with x-ray diffraction and x-radiography techniques. Thermoelastic properties of W were derived using different methods. We obtained the isothermal bulk modulus K T 0 = 310.3(1.5) GPa, its pressure derivative K ' T 0 = 4.4(3), its temperature derivative at constant pressure ( ∂ K T / ∂ T ) P = – 0.0138 ( 1 ) GPa K – 1 and at constant volume ( ∂ K T / ∂ T ) V = – 0.0050 GPa K – 1 , the thermal expansion α ( 0, T ) = 1.02(27) × 10 –5 + 7.39(3.2) × 10 –9 T (K –1 ), as well as the pressure derivative of thermal expansion ( ∂ α / ∂ P ) T = – 1.44 ( 1 ) × 10 – 7 K – 1 GP a – 1 based on the high-temperature Birch–Murnaghan equation of state (EOS), the Vinet EOS, and thermal pressure approach. Finite strain analysis allowed us to derive the elastic properties and their pressure/temperature derivatives independent of the choice of pressure scale. A least-squares fitting yielded K S 0 = 314.5(2.5) GPa, K S 0 ' = 4.45(9), (∂ K S /∂ T ) P = – 0.0076(6) GPa K –1 , G 0 = 162.4(9) GPa, G 0 ' = 1.8(1), (∂ G /∂ T ) P = – 0.0175(9) GPa K –1 , and α 298 K = 1.23 × 10 – 5 K – 1 . Fitting current data to the Mie–Grüneisen–Debye EOS with derived θ 0 = 383.4 K yielded γ 0 = 1.81 ( 6 ) and q = 0.3 . The thermoelastic parameters obtained from various approaches are consistent with one another and comparable with previous results within uncertainties. Our current study provides a complete and self-consistent dataset for the thermoelastic properties of tungsten at high P – T conditions, which is important to improve the theoretical modeling of these materials under dynamic conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High-Pressure Phase Transition of Metastable Wurtzite-Like CuInSe 2 Nanocrystals

Ternary I–III–VI 2 semiconductors, such as CuInSe 2 , exhibit diverse polymorphs with unique structural characteristics and optoelectronic properties. This study investigates the pressure-induced phase transitions of metastable wurtzite-like CuInSe 2 nanocrystals. Using a combination of synchrotron X-ray diffraction, pair distribution function analysis, and density functional theory calculations, we reveal a transition from cation-ordered wurtzite-like (Pmc2 1 ) to cation-disordered NaCl-like (Fm3̅m) structures at 7.7 GPa. The cation-disordered NaCl-like phase persists upon decompression. Bulk modulus calculations highlight size-dependent deviations from bulk material behavior. These findings deepen our understanding of phase stability in colloidal I–III–VI 2 semiconductor nanocrystals, with implications for tailoring functional materials under extreme conditions.

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

Structure, Stability, and Electronic Properties of Boron Suboxide: A Density Functional Theory Study

Boron suboxide (B 6 O) is a boron-rich compound derived from the α-rhombohedral boron lattice with extreme hardness and unusual semiconducting properties. In this work, density functional theory (DFT) was used to show that unit cell volume, mechanical strength, band gaps, and thermodynamic stabilities of B 6 O were influenced by the interstitial elements and point defects at the icosahedral sites. While the hexagonal unit cell volume (HUCV) varies with interstitial occupancy, it is the icosahedral defect that weakens the intrinsic bulk modulus of B 6 O. Using the hybrid HSE functional, we confirmed that the perfect B 6 O bulk is a p-type semiconductor with a direct band gap of 2.8 eV. Furthermore, by screening α-boron compounds systematically, we found that a simple octet rule may offer a consistent explanation for the variations in the computed electronic structures. The formation free energies calculated over a wide range of temperatures (0–2500 K) and pressures (0–80 GPa) predict that formations of interstitial defects become favorable only at higher temperatures (ca. 1800 K) in bulk B 6 O lattices. The nudged elastic band (NEB) method was employed to identify the minimum energy pathways for the diffusions of dislocated B and O atoms. The diffusion of icosahedral B atoms has an energy barrier of 0.16 eV. More complex B diffusion paths involving the reorganization of icosahedral boron atoms incur higher barriers (>1 eV). In contrast, the diffusion of interstitial O atoms is facile with a barrier of 0.4 eV. Lastly, successive O insertions into the α-B lattice were performed using DFT to generate a basic understanding of the oxidation process. Furthermore, these calculations provide fundamental atomistic insights into the growth of B 6 O crystals and control of their point defects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Origin of Heating-Induced Softening and Enthalpic Reinforcement in Elastomeric Nanocomposites

Molecular simulations demonstrate that the enthalpic softening of elastomeric nanocomposites upon heating can arise naturally from a Poisson’s ratio mismatch between elastomer and nanoparticle networks, providing a more parsimonious explanation for this phenomenon than the widely accepted interpretation based on glassy interparticle bridging. Despite a century of use, the mechanism of nanoparticle-driven mechanical reinforcement of elastomers is unresolved. Here, a major hypothesis attributes it to glassy interparticle bridges, supported by an observed inversion of the variation of the modulus E(T) on heating – from entropic stiffening in elastomers to enthalpic softening in nanocomposites. Here, molecular simulations reveal that elastomer enthalpic softening can instead emerge from a competition over the preferred volumes between elastomer and nanoparticulate networks. A theory for this competition accounting for softening of the bulk modulus on heating predicts the simulated E(T) inversion, suggesting that reinforcement is driven by a volume-competition mechanism unique to cocontinuous systems of soft and rigid networks.

Biopolymers↗