Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large displacement”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 343 records · Page 19

Helium causing disappearance of a/2<111> dislocation loops in binary Fe-Cr ferritic alloys

In this work, single and dual-beam self-ion irradiations were performed at 500°C on ultra-high purity Fe14%Cr alloy to ~0.33 displacements-per-atom (dpa), and 0 or 3030 atomic-parts-per-million (appm) helium/dpa, respectively. Using transmission electron microscopy, we reveal that helium can drastically modify the dislocation loop Burgers vector in Fe-Cr alloys. Helium co-implantation caused complete disappearance of a/2<111> type dislocation loops, and the microstructure consisted of only a<100> loops. Conversely, a/2<111> type loops were predominant without He co-implantation. The total loop density remained largely unaffected. The results strikingly contrast literature asserting that helium stabilizes a/2<111> type loops in bcc Fe alloys, based on low temperature irradiations. Collectively analyzing the results with literature suggest that the small positive interaction between helium and self-interstitial atoms (SIA) in Fe predicted by atomistic simulations maybe insufficient to holistically explain the dislocation loop microstructure development in presence of helium. Helium-SIA positive binding inadvertently implies elevated a/2<111> loop fraction and higher loop densities that the present results contradict. Helium induced high cavity density causing a preferential loss of highly glissile <111> clusters, leaving the matrix saturated with <100> type clusters is proposed as a potential mechanism. Further, the in-situ irradiations combined with Burgers vector analysis strengthened the evidence of Cr-induced dislocation loop mobility reduction that appears to stabilize the a/2<111> type loops and causes higher loop densities in Fe-Cr alloys.

36 MATERIALS SCIENCE↗

Anisotropic Structural Collapse of Mg 3 Sb 2 and Mg 3 Bi 2 at High Pressure

Alloys between Mg 3 Sb 2 and Mg 3 Bi 2 have recently been shown to be exceptional thermoelectric materials due in part to their anomalously low thermal conductivity. In the present study, in situ high-pressure synchrotron X-ray diffraction was used to investigate the structure and bonding in Mg 3 Sb 2 and Mg 3 Bi 2 at pressures up to 50 GPa. Our results confirm prior predictions of isotropic in-plane and out-of-plane compressibility but reveal large disparities between the bond strength of the two distinct Mg sites. Using single-crystal diffraction, we show that the octahedral Mg–Sb bonds are significantly more compressible than the tetrahedral Mg–Sb bonds in Mg 3 Sb 2 , which lends support to prior arguments that the weaker octahedral Mg bonds are responsible for the anomalous thermal properties of Mg 3 Sb 2 and Mg 3 Bi 2 . Further, we report the discovery of a displacive and reversible phase transition in both Mg 3 Sb 2 and Mg 3 Bi 2 above 7.8 and 4.0 GPa, respectively. Furthermore, the transition to the high-pressure structure involves a highly anisotropic volume collapse, in which the out-of-plane axis compresses significantly more than the in-plane axes. Single-crystal diffraction at high pressure was used to solve the monoclinic high-pressure structure ( C 2/ m ), which is a distorted variant of the ambient-pressure structure containing four unique Mg coordination environments.

36 MATERIALS SCIENCE↗

Deciphering the Mechanism of Base-Triggered Conversion of Ammonia to Molecular Nitrogen and Methylamine to Cyanide

Here we report an in-depth investigation into the ammonia oxidation mechanism by the catalyst [Ru III (tpy)(dmabpy)NH 3 ] 3 + ([Ru(NH 3 )] 3 + ). Stoichiometric reactions of [Ru(NH 3 )] 3 + were carried out with exogenous noncoordinating bases to trigger a proposed redox disproportionation reaction, which was followed using variable-temperature NMR spectroscopy. An intermediate species was identified as a dinitrogen-bridged complex using 15 N NMR and Raman spectroscopy on isotopically labeled complexes. This intermediate is proposed to derive from coupling of nitridyl species formed upon sequential redox disproportion reactions. Acetonitrile displaces the dinitrogen bridge to yield free N 2 . DFT calculations support this lower-energy pathway versus that previously reported for ammonia oxidation by the parent [Ru III (tpy)(bpy)NH 3 ] 3 + complex. These experimental and computational results are consistent with the interpretation of redox disproportionation involving sequential hydrogen atom transfer reactions by an amide/aminyl intermediate, [Ru(NH 2 ) – ] + $\Leftrightarrow$ [Ru(NH 2 ) • ] + , formed upon deprotonation of the parent complex. Control experiments employing a large excess of ammonia as a base indicate this new proposed lower-energy pathway contributes to the oxidation of ammonia to dinitrogen in conditions relevant to electrocatalysis. In addition, analogous methylamine complexes, [Ru(NH 2 CH 3 )] 2 + /3 + , were prepared to further test the proposed mechanism. Treating [Ru(NH 2 CH 3 )] 3 + with a base cleanly yields two products [Ru(NH 2 CH 3 )] 2 + and [Ru(CN)] + in an ~3:1 ratio, fully consistent with the proposed cascade of hydrogen atom transfer reactions by an intermediate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Topanga: A kinetic ion plasma code for large-scale ionospheric simulations on magnetohydrodynamic timescales

Topanga is a kinetic ion code developed for simulating large-scale plasma phenomena in the Earth's ionosphere on magnetohydrodynamic timescales. It is a domain-decomposed parallel code that runs on high-performance computing platforms. Features of Topanga include spherical geometry for simplified boundary conditions and computational efficiency; a hybrid plasma model with inertia-less fluid electrons, kinetic ions, and an electric field specified via an Ohm's law; a Maxwell-FDTD (finite difference time domain) plasma model which retains the displacement current in Maxwell's equations and models electron currents in the ionosphere with a tensor conductivity; sponge-layer boundary conditions for absorption of electromagnetic and plasma waves incident on the domain boundaries; and a novel mixed-implicit algorithm for evolving the EM fields inside the Maxwell-FDTD region that is stable over many orders of magnitude in the electron–ion collision frequency. We verify the numerical methods used in Topanga on a pair of test problems. The first test involves modeling a three-dimensional collisionless shock using the hybrid set of equations. The second test involves modeling a spherical TEM mode in vacuum using the Maxwell-FDTD set of equations. Finally, we demonstrate how using the combined set of hybrid and Maxwell-FDTD equations to model the Starfish Prime high-altitude nuclear test recovers a “missing” EM signal on the ground that is not present when using only the hybrid set of equations. The magnitude of this signal in the simulation containing the Maxwell-FDTD region agrees well with the E3a portion of the magnetohydrodynamic electromagnetic pulse from Starfish Prime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generation of random geological models using multi-randomization for machine learning

Generating high-fidelity geological models is essential for advancing machine learning (ML) methods in automated seismic interpretation. For instance, seismic images paired with corresponding fault labels are foundational for ML-based fault detection from seismic migration sections. While several open-access datasets of random geological models exist, open-source tools specifically designed to produce large volumes of such models for ML applications remain scarce. To address this gap, we present RGM (Random Geological Model), an open-source software package for efficiently generating 2D and 3D synthetic geological models tailored for ML workflows. RGM supports the creation of diverse model components, including medium property distributions (P-/S-wave velocities and density), seismic reflectivity images (i.e., synthetic migration sections), relative geological time, and discrete fault attributes such as probability, dip, strike, rake, and displacement. It also accommodates the creation of complex geological features such as salt bodies and unconformities. The model generation algorithm employs a multi-randomization strategy, yielding an effectively infinite-dimensional model space that encompasses a wide range of geological scenarios and associated seismic features. Furthermore, RGM incorporates a method to generate synthetic elastic migration images using analytical elastic reflection coefficients combined with frequency-dependent scaling. This functionality enables the creation of training datasets for ML models that leverage elastic seismic images. RGM is implemented in modern object-oriented Fortran, allowing users to flexibly control statistical parameters governing model variability. We demonstrate the capability, performance, and geological realism of the package through comprehensive 2D and 3D examples.

58 GEOSCIENCES↗

Neutrino masses from low scale partial compositeness

We consider a class of models in which the neutrinos acquire Majorana masses through mixing with singlet neutrinos that emerge as composite states of a strongly coupled hidden sector. In this framework, the light neutrinos are partially composite particles that obtain their masses through the inverse seesaw mechanism. We focus on the scenario in which the strong dynamics is approximately conformal in the ultraviolet, and the compositeness scale lies at or below the weak scale. The small parameters in the Lagrangian necessary to realize the observed neutrino masses can naturally arise as a consequence of the scaling dimensions of operators in the conformal field theory. We show that this class of models has interesting implications for a wide variety of experiments, including colliders and beam dumps, searches for lepton flavor violation and neutrinoless double beta decay, and cosmological observations. At colliders and beam dumps, this scenario can give rise to striking signals involving multiple displaced vertices. The exchange of hidden sector states can lead to observable rates for flavor violating processes such as μ → eγ and μ → e conversion. If the compositeness scale lies at or below a hundred MeV, the rate for neutrinoless double beta decay is suppressed by form factors and may be reduced by an order of magnitude or more. The late decays of relic singlet neutrinos can give rise to spectral distortions in the cosmic microwave background that are large enough to be observed in future experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Eastern Interconnection Wide-Area Oscillation Assessment and Study Report

To meet the government mandates such as Renewable Portfolio Standards (RPS), conventional synchronous generators with rotating mass are being displaced by inverter-based resources (IBRs), leading to high penetration of renewable energy sources (RES). As a result, the inherent properties and characteristics of the transformed generation mix may significantly impact the grid behaviour, and it is of paramount significance that the impact and the consequences of the resource mix change is thoroughly understood to adopt measures to maintain reliable grid operations. Of great concern is the inter-area low-frequency oscillation, which usually propagates through a large region and has a system-wide impact. Such oscillations may lead to unnecessary or inadvertent tripping of generators, that may be simply reacting to the oscillation originating from geographically remote sites. Such tripping of generators can lead to cascading outages, system split and load loss events. Over the years, several system-wide oscillation events have been observed across all three North American interconnections. This project titled Wide-area Oscillation Assessment and Trending Study, sponsored by the Office of Electricity (OE) of the Department of Energy (DOE), aims to conduct the required research to capture and investigate the potential changes in wide-area oscillatory performance of the system as a result of resource mix transition. The Pacific Northwest National Laboratory (PNNL) conducted this research and technical staff of Federal Energy Regulatory Commission (FERC) served as advisors to the project. The primary objective of this study is to assess whether there are significant trends in the power system wide-area oscillatory behaviours, as a result of the generation-mix changes due to the increased penetration of RES in the U.S. Eastern Interconnection (EI). The analysis presented in this report can serve as reference for future grid planning and operation of the EI system. This research conducted in this project evaluates the wide-area oscillatory behaviour of the EI, using a model-based approach for potential future resource mixes, and the measurement-based analysis with 21-month phasor measurement unit (PMU) time series. The oscillation modes in EI system along with the impact of the changes of the generation mix on the frequency and damping ratio (DR) of the oscillations have been studied.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗

Quantitative proton radiography and shadowgraphy for arbitrary intensities

Charged-particle radiography and shadowgraphy data can be directly inverted to obtain a line-integrated transverse Lorentz force or a line-integrated transverse refractive index gradient if intensity modulations due to scattering and absorption are negligible, and angular deflections are small. We develop a new direct-inversion algorithm based on plasma physics and compare it to a new Monge–Ampère code and an existing power diagram code. The measured or source intensity is represented by electrons subject to drag, and the other intensity by fixed ions. The decrease in kinetic plus electrostatic energy determines convergence. The displacement of the electrons from their initial to their equilibrium positions determines the line-integrated force or refractive index gradient. We have implemented two approaches: PIC (particle in cell) and Lagrangian fluid, in 1-D and 2-D. The PIC code works for arbitrary intensities, can work efficiently in parallel, and can make use of existing codes. The Lagrangian code requires less memory and is faster than the PIC code without massively parallel processing, but fails in 2-D for large intensity modulations. The Monge–Ampère code is by far the fastest in 2-D, without massively parallel processing, but fails for intensities with large voids, high contrast ratios and large deflections across the boundaries, and could not obtain the degree of convergence possible with the PIC code. As a result, the power diagram code was by far the slowest and most memory intensive, and failed for large peaks in the measured intensity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Lagrangian perturbation theory in the presence of massive neutrinos

We develop a Lagrangian Perturbation Theory (LPT) framework to study the clustering of cold dark matter (CDM) in cosmologies with massive neutrinos. We follow the trajectories of CDM particles with Lagrangian displacements fields up to third order in perturbation theory. Once the neutrinos become non-relativistic, their density fluctuations are modeled as being proportional to the CDM density fluctuations, with a scale-dependent proportionality factor. This yields a gravitational back-reaction that introduces additional scales to the linear growth function, which is accounted for in the higher order LPT kernels. Through non-linear mappings from Eulerian to Lagrangian frames, we ensure that our theory has a well behaved large scale behavior free of unwanted UV divergences, which are common when neutrino and CDM densities are not treated on an equal footing, and in resummation schemes that manifestly break Galilean invariance. We use our theory to construct correlation functions for both the underlying matter field, as well as for biased tracers using ConvolutionLPT. Redshift-space distortions effects are modeled using the Gaussian Streaming Model. When comparing our analytical results to simulated data from the Quijote 1 simulation suite, we find good accuracy down to r = 20 Mpc h -1 at redshift z = 0.5, for the real space and redshift space monopole particle correlation functions with no free parameters. The same accuracy is reached for the redshift space quadrupole if we additionally consider an effective field theory parameter that shifts the pairwise velocity dispersion. For modeling the correlation functions of tracers we adopt a simple Lagrangian biasing scheme with only density and curvature operators, which we find sufficient to reach down to r = 20 Mpc h -1 when comparing to simulated halos.

79 ASTRONOMY AND ASTROPHYSICS↗

An optimization-based approach to tailor the mechanical response of soft metamaterials undergoing rate-dependent instabilities

An optimization-based design framework is proposed to tune the response of soft metamaterials involving both geometric instabilities and nonlinear viscoelastic material behavior. Designing the response of soft metamaterials to harness instabilities and undergo large, tailored configuration changes will enable advancements in soft robotics, shock and vibration mitigation, and flexible electronics. In line with the metamaterial concept, the response of these materials is governed to a large extent by the geometric and topological makeup of their small-scale features. However, the link between structure and response is less intuitive for soft metamaterials due to their reliance upon highly nonlinear responses triggered by geometric instabilities. This is further complicated by the effects of viscoelastic relaxation, which recent studies have shown to alter the emergence of instabilities in non-intuitive ways. Here, these effects are accounted for in our framework to achieve various design objectives, including tailored force–displacement response and maximized energy absorption from both geometric and material effects. To fully automate this process, it is essential to have a completely robust equation solver for forward problems involving instabilities and viscoelastic relaxation. We achieve this by casting the search for stable mechanical equilibrium — i.e. the forward problem — as a minimization problem and utilize a trust region algorithm to robustly handle instabilities and follow energetically-favorable equilibrium paths through critical points.

97 MATHEMATICS AND COMPUTING↗

Ionization-driven competitive (recovery) process in pre-damaged KTaO 3 : A brief review

The nuclear (S n ) and electronic (S e ) energy dissipation processes have been considered to be independent and largely uncorrelated, influencing our understanding of ion–solid interaction and damage processes in the last decades. Recently, however, it has become more generally accepted that S n and S e are coupled as they interact both in time and space. To decouple these processes, separating these effects in experiments using sequential dual-beam irradiations have become accepted as the logical path to advance the understanding of complex interactions between S e and pre-existing defects that may be created from displacement events. This experimental approach has been recently applied to studies of KTaO 3 to reveal new insights into this critical research topic. Here, we offer a forward-looking and comprehensive perspective on the fundamental coupling between Se and pre-existing defects in KTaO 3 . The origins behind the competitive two-stage phase transition process leading to damage healing are revealed and discussed. Furthermore, the evidence resulting from synergistic effects is also included for comparison. Additionally, our findings are rationalized using both Se and the ion velocity as key parameters. We highlight how the inelastic thermal spike (i-TS) calculations provide insights into the nature of this coupled process and further confirm that the ion velocity effect governs annealing kinetics. This work emphasizes that through the introduction of a small amount of local disorder in materials, MeV ion irradiation (i.e., not extreme S e ) may also be one additional option in subsequent material modification and functionalization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Experimental Modal Analysis Using Phase Quantities from Phase-Based Motion Processing and Motion Magnification

Phase-based motion processing and the associated Motion Magnification that it enables has become popular not only for the striking videos that it can produce of traditionally stiff structures visualized with very large deflections, but also for its ability to pull information out of the noise floor of images so that they can be processed with more traditional optical techniques such as digital image correlation or feature tracking. While the majority of papers in the literature have utilized the Phase-based Image Processing approach as a pre-processor for more quantitative analyses, the technique itself can be used directly to extract modal parameters from an image, noting that the extracted phases are proportional to displacements in the image. Therefore, once phases are extracted, they can be fit using traditional experimental modal analysis techniques. This produces a mode “shape” where the degrees of freedom are phases instead of physical motions. These phases can be scaled to produce on-image visualizations of the mode shapes, rather than operational shapes produced by bandpass filtering. Modal filtering techniques can also be used to visualize motions from an environment on an image using the modal phases as a basis for the expansion.

47 OTHER INSTRUMENTATION↗

Metrics and extrapolation of resonant magnetic perturbation thresholds for ELM suppression

This large database study of resonant magnetic perturbation (RMP) edge localized mode (ELM) suppression thresholds in the AUG, DIII-D, EAST, and KSTAR tokamaks details the key strengths and weaknesses of RMP metrics. The RMP ELM suppression database used for this work contains plasma information at the time of transition from ELMing to ELM suppressed states where a clear experimental threshold is identified. The experimental threshold distributions are compared for five metrics: (1) the island overlap width, (2) pedestal top Chirikov overlap, (3) peeling edge displacement, (4) pedestal top resonant drive, and (5) edge dominant mode overlap. The distributions, the regularity of the dependence on RMP coil currents, and the sensitivities of a given metric to equilibrium reconstruction details are compared. The overlap metric proves to be a good compromise between including the appropriate plasma response physics and maintaining a numerical robustness. This quantity does not exhibit clear power-law scalings for projection, but machine learning can assist in predicting thresholds within the existing parameter ranges and providing uncertainty quantification of those predictions. Two new first-principles models, one utilizing a threshold from the non-linear Modified Rutherford equation evaluated at the pedestal top and one utilizing the SLAYER code to calculate the linear tearing threshold from torque balance, offer possible paths to extrapolation beyond the existing database parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗

Twisted symmetric trilayer graphene: Single-particle and many-body Hamiltonians and hidden nonlocal symmetries of trilayer moiré systems with and without displacement field

Here, we derive the Hamiltonian for trilayer moiré systems with the Coulomb interaction projected onto the bands near the charge neutrality point. Motivated by the latest experimental results, we focus on the twisted symmetric trilayer graphene (TSTG) with a mirror symmetry with respect to the middle layer. We provide a full symmetry analysis of the noninteracting Hamiltonian with a perpendicular displacement field coupling the band structure made otherwise of the twisted bilayer graphene (TBG) and the high-velocity Dirac fermions, and we identify a hidden nonlocal symmetry of the problem. In the presence of this displacement field, we construct an approximate single-particle model, akin to the tripod model for TBG, capturing the essence of noninteracting TSTG. We also derive more quantitative perturbation schemes for the low-energy physics of TSTG with displacement field, obtaining the corresponding eigenstates. This allows us to obtain the Coulomb interaction Hamiltonian projected in the active band TSTG wave functions and derive the full many-body Hamiltonian of the system. We also provide an efficient parametrization of the interacting Hamiltonian. Finally, we show that the discrete symmetries at the single-particle level promote the U (2) × U (2) spin-valley symmetry to enlarged symmetry groups of the interacting problem under different limits. The interacting part of the Hamiltonian exhibits a large U (4) × U (4) × U (4) × U (4) symmetry in the chiral limit. Moreover, by identifying a symmetry which we dub spatial many-body charge conjugation, we show that the physics of TSTG is symmetric around charge neutrality.

2-dimensional systems↗

Node-Accessible Zirconium MOFs

High-stability, zirconium-based metal–organic frameworks are attractive as heterogeneous catalysts and as model supports for uniform arrays of subsequently constructed heterogeneous catalysts—for example, MOF-node-grafted metal–oxy and metal–sulfur clusters. For hexa-Zr(IV)-MOFs characterized by nodes that are less than 12-connected, sites not used for linkers are ideally occupied by reactive and displaceable OH/H 2 O pairs. Here, the desired pairs are ideal for grafting the aforementioned catalytic clusters, while aqua-ligand lability renders them effective for exposing highly Lewis-acidic Zr(IV) sites (catalytic sites) to candidate reactants. New single-crystal X-ray studies of an eight-connected Zr-MOF, NU-1000, reveal that conventional activation fully removes modulator ligands, but replaces them with three node-blocking formate ligands (from solvent decomposition) and only one OH/H 2 O pair, not four—a largely overlooked complication that now appears to be general for Zr-MOFs. Here we describe an alternative activation protocol that effectively removes modulators, avoids formate, and installs the full complement of terminal OH/H 2 O pairs. It does so via an unusual isolatable intermediate featuring eight aqua ligands and four non-ligated chlorides—again as supported by single-crystal X-ray data. We find that complete replacement of node-blocking modulators/formate with the originally envisioned OH/OH 2 pairs has striking consequences; here we touch upon just three. First, elimination of unrecognized formate renders aqua ligands much more thermally labile, enabling open Zr(IV) sites to be obtained at lower temperature. Second, in the absence of formate, which otherwise links and locks pairs of node Zr(IV) ions, reversible removal of aqua ligands engenders reversible contraction of MOF meso- and micropores, as evidenced by X-ray diffraction. Third, formate replacement with OH/OH 2 pairs renders NU-1000 ca.10× more active for catalytic hydrolytic degradation of a representative simulant of G-type chemical warfare agents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rethinking PCMI Failures Under Transient Conditions

Under PCMI conditions, the cladding strain is driven by either pellet swelling during normal operation or pellet thermal expansion in transients. The loading in PCMI is displacement controlled because friction between the pellet and the cladding results in a strain compatibility at the pellet/cladding interface. Failure is known to occur under PCMI loading conditions during operation when normal power changes result in release of iodine (and other fission products) causing stress corrosion cracking where the crack-propagation rate is dependent upon the cladding stress but slow relative to strictly mechanical failures. Failures after normal power changes are often referred to as failure by pellet-cladding interaction or PCI to distinguish it from PCMI. Such failures are considered displacement-controlled failures, as noted above, and occur with limited plastic strain. This is different to what is typically thought of as PCMI failure under transient conditions, either Anticipated Operational Occurrences (AOOs) or Reactivity Insertion Accidents (RIAs), which is generally assumed to occur once the cladding reaches a mechanical strain limit. But there are two features in the typical PCMI failure in an RIA that are not consistent with the PCMI assumptions. The first is the ductile shear failure feature in cross section and the second is a large crack opening post failure (e.g., see OI-11 failure in NSRR program). The ductile feature is indicative of a plastic instability that requires very large local strains, and the large crack opening must be the result of a release of energy, both of which are not expected if stress in the cladding was only driven by thermal expansion in the pellet. It is concluded that there must be a mechanism (likely related to rod internal pressure) that eliminates strain compatibility between the pellet and cladding – allowing strain localization to occur in what are known as PCMI transient failures. In some sense, these PCMI failures may actually be early burst failures that occur when thermal expansion causes uniform strains greater than a defined limit and fission gas pressure eliminate the strain compatibility

Cantonwine, Paul↗