Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity”

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

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

At least 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Analytic formulations for one-dimensional decay of rectangular homoepitaxial islands during coarsening on anisotropic fcc (110) surfaces

Submonolayer homoepitaxial fcc (110) systems display behavior reflecting strong anisotropy at lower temperatures, including one-dimensional decay during Ostwald ripening of rectangular islands maintaining constant width in the (001) direction. To appropriately describe this behavior, we first develop a refined continuum Burton-Cabrera-Frank formalism, which accounts for a lack of equilibration of island shape and importantly also for inhibited incorporation of adatoms at almost-faceted (1¯10) island edges through effective kinetic coefficients. This formalism is shown to describe accurately the adatom diffusion fluxes between islands and thus island evolution for a complex experimental island configuration, as confirmed by matching results from realistic atomistic simulations for this configuration. This approach also elucidates basic dependencies of flux on island geometry and temperature. Second, a further refinement is presented incorporating separate terrace and edge adatom density fields either in a continuum setting or alternatively in a spatially discrete diffusion equation setting. The second approach allows more flexibility and accuracy in accounting for edge-diffusion kinetics including corner rounding, a lack of equilibration of the edge adatom density at (1¯10) island edges, and the effect of rare kinks on (1¯10) island edges. Finally and significantly, it suggests facile two-way corner rounding at the island periphery during island decay, contrasting the previous picture.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structural origin of disorder-induced ion conduction in NaFePO 4 cathode materials

Diffusion in NaFePO 4 can be enhanced through amorphization. Based on computations using DFT and machine learning potentials, we ascribe this phenomenon to the formation of less constrained Na-ion environments upon disordering. Most modern battery technologies depend on solid-state crystalline cathode materials. However, some of these materials are constrained by the low ionic conductivity of their most stable phases. An example of this is maricite (NaFePO 4 ). Interestingly, experiments have shown that maricite can improve its rate capability through disordering (amorphization). However, experimental characterization of amorphous cathode materials remains a major challenge, hindering a clear understanding of the structural origin of the disorder-induced improvement in sodium-ion mobility. To address this, we here employ molecular dynamics simulations by first training a machine learning potential for NaFePO 4 based on the atomic cluster expansion approach and a batch active learning potential parameterization scheme. This potential is then applied to explore the structural and dynamical properties of NaFePO 4 glasses as cathode materials. Specifically, we investigate the effect of glass structure on sodium-ion diffusion, revealing the relative influences of short-range and medium-range order features. We find significant heterogeneity in sodium-ion diffusivity in the glass, with fast-conducting ions residing in less constrained atomic environments with fewer P and Fe neighbors. These more mobile ions are also surrounded by larger ring-type structures. Overall, the results and developed approach present promising avenues for developing high-performance glassy cathodes for next-generation batteries.

Christensen, Rasmus↗
DOE OSTI2019

Dynamics of ionic liquids in the presence of polymer-grafted nanoparticles

We incorporated polymer-grafted nanoparticles into ionic and zwitterionic liquids to explore the solvation and confinement effects on their heterogeneous dynamics using quasi-elastic neutron scattering (QENS). 1-Hexyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (HMIM-TFSI) mixed with deuterated poly(methyl methacrylate) (d-PMMA)-grafted nanoparticles is studied to unravel how dynamic coupling between PMMA and HMIM-TFSI influence the fast and slow diffusion characteristics of the HMIM + cations. The zwitterionic liquid, 1-butyl-3-methyl imidazole-2-ylidene borane (BMIM-BH 3 ) is critically selected and mixed with PMMA-grafted nanoparticles for comparison in this work as its ions do not self-dissociate and it does not couple with PMMA through ion-dipole interactions as HMIM-TFSI does. We find that long-range unrestricted diffusion of HMIM + cations is higher in well-dispersed particles than in aggregated particle systems, whereas the localized diffusion of HMIM + is measured to be higher in close-packed particles. Translational diffusion dynamics of BMIM-BH 3 is not influenced by any particle structures suggesting that zwitterions do not interact with PMMA. In conclusion, this difference between two ionic liquid types enables us to decouple polymer effects from the diffusion of ionic liquids, which is integral to understand the ionic transport mechanism in ionic liquids confined in polymer-grafted nanoparticle electrolytes.

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

Planck 2013 results. XI. All-sky model of thermal dust emission

This work introduces an all-sky model of dust emission from the Planck 353, 545, and 857 GHz, and IRAS 100 μm data. Using a modified blackbody fit to the data we present all-sky maps of the dust optical depth, temperature, and spectral index over the 353-3000 GHz range. This model is a good representation of the IRAS and Planck data at 5' between 353 and 3000 GHz (850 and 100 μm). It shows variations of the order of 30% compared with the widely-used model of Finkbeiner, Davis, and Schlegel. The Planck data allow us to estimate the dust temperature uniformly over the whole sky, down to an angular resolution of 5', providing an improved estimate of the dust optical depth compared to previous all-sky dust model, especially in high-contrast molecular regions where the dust temperature varies strongly at small scales in response to dust evolution, extinction, and/or local production of heating photons. An increase of the dust opacity at 353 GHz, τ 353 / N H , from the diffuse to the denser interstellar medium (ISM) is reported. It is associated with a decrease in the observed dust temperature, T obs , that could be due at least in part to the increased dust opacity. We also report an excess of dust emission at H I column densities lower than 10 20 cm -2 that could be the signature of dust in the warm ionized medium. In the diffuse ISM at high Galactic latitude, we report an anticorrelation between τ 353 / N H and T obs while the dust specific luminosity, i.e., the total dust emission integrated over frequency (the radiance) per hydrogen atom, stays about constant, confirming one of the Planck Early Results obtained on selected fields. This effect is compatible with the view that, in the diffuse ISM, T obs responds to spatial variations of the dust opacity, due to variations of dust properties, in addition to (small) variations of the radiation field strength. The indication is that in the diffuse high-latitude ISM τ 353 is not as reliable a tracer of dust column density as we conclude it is in molecular clouds where the correlation of τ 353 with dust extinction estimated using colour excess measurements on stars is strong. To estimate Galactic E(B-V) in extragalactic fields at high latitude we develop a new method based on the thermal dust radiance, instead of the dust optical depth, calibrated to E(B -V) using reddening measurements of quasars deduced from Sloan Digital Sky Survey data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Dynamic behavior of CH 3 NH 3 PbI 3 perovskite twin domains

We report the recent discovery of twin domains in MAPbI 3 perovskites has initiated contentious discussion on the ferroic nature of hybrid perovskites. Ferroelectric polarization is thought to facilitate the dissociation of photoinduced electron-hole pairs, helping to explain the extraordinary photovoltaic performance exhibited by this class of materials. Alternate to ferroelectricity, which has yet to be unambiguously established despite considerable efforts to do so, ferroelasticity was also proposed in these materials. Meanwhile, given the coupling of ionic states and ferroelectricity and the interconnected nature of defect chemistry and ferroelasticity, the electrochemical reactivity can no longer be ignored. In this work, using band excitation piezoresponse force microscopy, we reveal the variation in elasticity between adjacent domains, indicating the ferroelasticity and the difference in the crystallographic states of the twin domain. Moreover, using band excitation contact Kelvin probe force microscopy, we dynamically map the evolution of the twinning structure under electric bias. These results help decipher the effect of the twin domains on ionic mobility and ion diffusion pathways. Combining these results, we reveal the interaction of twin domains and ionic activity in this material. Furthermore, this work provides insights into the twinning structure in MAPbI 3 and its potential effects on the hybrid perovskite optoelectronics.

36 MATERIALS SCIENCE↗
DOE OSTI2017

The arbitrary order mimetic finite difference method for a diffusion equation with a non-symmetric diffusion tensor

Here, we present the arbitrary order mimetic finite difference (MFD) discretization for the diffusion equation with non-symmetric tensorial diffusion coefficient in a mixed formulation on general polygonal meshes. The diffusion tensor is assumed to be positive definite. The asymmetry of the diffusion tensor requires changes to the standard MFD construction. We also present new approach for the construction that guarantees positive definiteness of the non-symmetric mass matrix in the space of discrete velocities. The numerically observed convergence rate for the scalar quantity matches the predicted one in the case of the lowest order mimetic scheme. For higher orders schemes, we observed super-convergence by one order for the scalar variable which is consistent with the previously published result for a symmetric diffusion tensor. The new scheme was also tested on a time-dependent problem modeling the Hall effect in the resistive magnetohydrodynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Interfacial phonon scattering and transmission loss in > 1 µm thick silicon-on-insulator thin films

Scattering of phonons at boundaries of a crystal (grains, surfaces, or solid/solid interfaces) is characterized by the phonon wavelength, the angle of incidence, and the interface roughness, as historically evaluated using a specularity parameter p formulated by Ziman [Electrons and Phonons (Clarendon Press, Oxford, 1960)]. This parameter was initially defined to determine the probability of a phonon specularly reflecting or diffusely scattering from the rough surface of a material. The validity of Ziman's theory as extended to solid/solid interfaces has not been previously validated. Here, to better understand the interfacial scattering of phonons and to test the validity of Ziman's theory, we precisely measured the in-plane thermal conductivity of a series of Si films in silicon-on-insulator (SOI) wafers by time-domain thermoreflectance (TDTR) for a Si film thickness range of 1–10 μm and a temperature range of 100–300 K. The Si/SiO 2 interface roughness was determined to be 0.11±0.04nm using transmission electron microscopy (TEM). Furthermore, we compared our in-plane thermal conductivity measurements to theoretical calculations that combine first-principles phonon transport with Ziman's theory. Calculations using Ziman's specularity parameter significantly overestimate values from the TDTR measurements. We attribute this discrepancy to phonon transmission through the solid/solid interface into the substrate, which is not accounted for by Ziman's theory for surfaces. The phonons that are specularly transmitted into an amorphous layer will be sufficiently randomized by the time they come back to the crystalline Si layer, the effect of which is practically equivalent to a diffuse reflection at the interface. Finally, we derive a simple expression for the specularity parameter at solid/amorphous interfaces and achieve good agreement between calculations and measurement values.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Structural and dielectric properties of CuO-doped SrTiO{sub 3} ceramics

Copper-doped SrTiO{sub 3} (ST) ceramic powders were processed by solid-state route diffusion (SRD) bulk preparation technique. We reported the effect of Cu{sup +2} ions on the dielectric response of ST and it established the substantial increase in dielectric constant (ε{sub r}) than undoped ST from 303K-673K and low loss (tanδ) for good dielectric applications. In respect of the electrical properties ac- activation energies were computed of range 0.089-0.319eV using lnσ Vs 1/T plots. The microstructure was examined with grain sizes 6-9µm of uniform distribution by field emission scanning electron microscope (FESEM). Some additional phases SrCu{sub 3}Ti{sub 4}O{sub 12} and TiO{sub 2} rutiles were detected by X-ray diffraction technique.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Phase-space entropy cascade and irreversibility of stochastic heating in nearly collisionless plasma turbulence

We consider a nearly collisionless plasma consisting of a species of “test particles” in one spatial and one velocity dimension, stirred by an externally imposed stochastic electric field—a kinetic analog of the Kraichnan model of passive advection. The mean effect on the particle distribution function is turbulent diffusion in velocity space—known as stochastic heating. Accompanying this heating is the generation of fine-scale structure in the distribution function, which we characterize with the collisionless (Casimir) invariant C 2 ∝ ∫ ∫ d x d v 〈 f 2 〉 —a quantity that here plays the role of (negative) entropy of the distribution function. We find that C 2 is transferred from large scales to small scales in both position and velocity space via a phase-space cascade enabled by both particle streaming and nonlinear interactions between particles and the stochastic electric field. We compute the steady-state fluxes and spectrum of C 2 in Fourier space, with k and s denoting spatial and velocity wave numbers, respectively. In our model, the nonlinearity in the evolution equation for the spectrum turns into a fractional Laplacian operator in k space, leading to anomalous diffusion. Whereas even the linear phase mixing alone would lead to a constant flux of C 2 to high s (towards the collisional dissipation range) at every k , the nonlinearity accelerates this cascade by intertwining velocity and position space so that the flux of C 2 is to both high k and high s simultaneously. Integrating over velocity (spatial) wave numbers, the k -space ( s -space) flux of C 2 is constant down to a dissipation length (velocity) scale that tends to zero as the collision frequency does, even though the rate of collisional dissipation remains finite. The resulting spectrum in the inertial range is a self-similar function in the ( k , s ) plane, with power-law asymptotics at large k and s . Our model is fully analytically solvable, but the asymptotic scalings of the spectrum can also be found via a simple phenomenological theory whose key assumption is that the cascade is governed by a “critical balance” in phase space between the linear and nonlinear timescales. We argue that stochastic heating is made irreversible by this entropy cascade and that, while collisional dissipation accessed via phase mixing occurs only at small spatial scales rather than at every scale as it would in a linear system, the cascade makes phase mixing even more effective overall in the nonlinear regime than in the linear one. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The spread of technological innovations: effects of psychology, culture and policy interventions

Technological innovations drive the evolution of human societies. The success of innovations depends not only on their actual benefits but also on how potential adopters perceive them and how their beliefs are affected by their social and cultural environment. To deepen our understanding of socio-psychological processes affecting the new technology spread, we model the joint dynamics of three interlinked processes: individual learning and mastering the new technology, changes in individual attitudes towards it, and changes in individual adoption decisions. We assume that the new technology can potentially lead to a higher benefit but achieving it requires learning. We posit that individual decision-making process as well as their attitudes are affected by cognitive dissonance and conformity with peers and an external authority. Individuals vary in different psychological characteristics and in their attitudes. We investigate both transient dynamics and long-term equilibria observed in our model. We show that early adopters are usually individuals who are characterized by low cognitive dissonance and low conformity with peers but are sensitive to the effort of an external authority promoting the innovation. We examine the effectiveness of five different intervention strategies aiming to promote the diffusion of a new technology: training individuals, providing subsidies for early adopters, increasing the visibility of peer actions, simplifying the exchange of opinions between people, and increasing the effort of an external authority. We also discuss the effects of culture on the spread of innovations. Finally, we demonstrate that neglecting the cognitive forces and the dynamic nature of individual attitudes can lead to wrong conclusions about adoption of innovations. Our results can be useful in developing more efficient policies aiming to promote the spread of new technologies in different societies, cultures and countries.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Using Active Learning to Rapidly Develop Machine Learned Diffusion Coefficients of CO 2 Conversion Reagents in Metal–Organic Frameworks

Here, we used a combined molecular dynamics/active learning (AL) approach to create machine learning models that can predict the diffusion coefficient of epichlorohydrin and chloropropene carbonate, the reactant and product of a common CO 2 cycloaddition reaction, in metal–organic frameworks (MOFs). Nanoporous MOFs are effective catalysts for the cycloaddition of CO 2 to epoxides. The diffusion rates within nanoporous catalysts can control the rate of reaction as the reactants and products must diffuse to the active sites within the MOF and then out of the nanoporous material for reusability. However, the diffusion process is routinely ignored when searching for new materials in catalytic applications. Here we verified improvement during the AL process by consistently tracking metrics on the same groups of MOFs to ensure consistency. Metal identity was found to have little impact on diffusion rates, while structural features like pore limiting diameter act as a threshold where a minimum value is needed for high diffusion rates. We identified the MOFs with the highest epichlorohydrin and chloropropene carbonate diffusion coefficients which can be used for further studies of reaction energetics.

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

Comprehensive studies of response characteristics of organic photodetectors based on rubrene and C{sub 60}

We studied the transient response characteristics of organic photodetectors composing of high mobility materials of rubrene and C{sub 60}, respectively, as donor and acceptor. It was found that the response speed was limited by the delay time of both exciton diffusion and transit, and an anomalous phenomenon that the bandwidth decreases as the reverse bias increases was found for the first time. The investigation of frequency dependence at different device structures and light excitations demonstrated that the light absorption of rubrene prevents the photodetector from obtaining a fast response. With the help of magnetic field effect study, it was clearly elucidated that the slow diffusion time of the long lifetime triplet excitons generated from singlet fission in rubrene limited the bandwidth of the device. Moreover, the simulation of the response of photodetector under transient and steady state by exciton transport-diffusion equation showed that the exciton dissociation efficiency in rubrene increases more quickly than that in C{sub 60}, which should account for the negative dependence of bandwidth on the reverse bias in rubrene-based device.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Optical polarization and spectral properties of the hydrogen-poor superluminous supernovae SN 2021bnw and SN 2021fpl

ABSTRACT New optical photometric, spectroscopic, and imaging polarimetry data are combined with publicly available data to study some of the physical properties of the two hydrogen-poor superluminous supernovae (SLSNe) SN 2021bnw and SN 2021fpl. For each SLSN, the best-fitting parameters obtained from the magnetar model with Modular Open-Source Fitter for Transients do not depart from the range of parameter obtained on other SLSNe discussed in the literature. A spectral analysis with SYN++ shows that SN 2021bnw is a W type, fast evolver, while SN 2021fpl is a 15bn type, slow evolver. The analysis of the polarimetry data obtained on SN 2021fpl at four epochs (+1.8, +20.6, +34.1, and +43.0 d, rest frame) shows >3σ polarization detections in the range of 0.8–1 per cent. A comparison of the spectroscopy data suggests that SN 2021fpl underwent a spectral transition a bit earlier than SN 2015bn, during which, similarly, it could have underwent a polarization transition. The analysis of the polarimetry data obtained on SN 2021bnw does not show any departure from symmetry of the photosphere at an empirical diffusion time-scale of ≈2 (+81.1 d rest frame). This result is consistent with those on the sample of W-type SLSN observed at empirical diffusion time-scale ≤ 1 with that technique, even though it is not clear the effect of limited spectral windows varying from one object to the other. Measurements at higher empirical diffusion time-scale may be needed to see any departure from symmetry as it is discussed in the literature for SN 2017egm.

Astronomy & Astrophysics↗
DOE OSTI2015

Experimental study on effects of geologic heterogeneity in enhancing dissolution trapping of supercritical CO 2

Abstract Dissolution trapping is one of the primary mechanisms that enhance the storage security of supercritical carbon dioxide (scCO 2 ) in saline geologic formations. When scCO 2 dissolves in formation brine produces an aqueous solution that is denser than formation brine, which leads to convective mixing driven by gravitational instabilities. Convective mixing can enhance the dissolution of CO 2 and thus it can contribute to stable trapping of dissolved CO 2 . However, in the presence of geologic heterogeneities, diffusive mixing may also contribute to dissolution trapping. The effects of heterogeneity on mixing and its contribution to stable trapping are not well understood. The goal of this experimental study is to investigate the effects of geologic heterogeneity on mixing and stable trapping of dissolved CO 2 . Homogeneous and heterogeneous media experiments were conducted in a two‐dimensional test tank with various packing configurations using surrogates for scCO 2 (water) and brine (propylene glycol) under ambient pressure and temperature conditions. The results show that the density‐driven flow in heterogeneous formations may not always cause significant convective mixing especially in layered systems containing low‐permeability zones. In homogeneous formations, density‐driven fingering enhances both storage in the deeper parts of the formation and contact between the host rock and dissolved CO 2 for the potential mineralization. On the other hand, for layered systems, dissolved CO 2 becomes immobilized in low‐permeability zones with low‐diffusion rates, which reduces the risk of leakage through any fault or fracture. Both cases contribute to the permanence of the dissolved plume in the formation.

Agartan, Elif↗
DOE OSTI2024

Validating a microphysical prognostic stratospheric aerosol implementation in E3SMv2 using observations after the Mount Pinatubo eruption

Abstract. This paper describes the addition of a stratospheric prognostic aerosol (SPA) capability – developed with the goal of accurately simulating sulfate aerosol formation and evolution in the stratosphere – in the Department of Energy (DOE) Energy Exascale Earth System Model, version 2 (E3SMv2). The implementation includes changes to the four-mode Modal Aerosol Module microphysics in the stratosphere to allow for larger particle growth and more accurate stratospheric aerosol lifetime following the Pinatubo eruption. E3SMv2-SPA reasonably reproduces stratospheric aerosol lifetime, burden, aerosol optical depth, and top-of-atmosphere flux when compared to remote sensing observations. E3SMv2-SPA also has close agreement with the interactive chemistry–climate model CESM2-WACCM (Community Earth System Model version 2–Whole Atmosphere Community Climate Model) – which has a more complete chemical treatment – and the observationally constrained, prescribed volcanic aerosol treatment in E3SMv2. Global stratospheric aerosol size distributions identify the nucleation and growth of sulfate aerosol from volcanically injected SO2 from both major and minor volcanic eruptions from 1991 to 1993. The modeled aerosol effective radius is consistently lower than satellite and in situ measurements (max differences of ∼ 30 %). Comparisons with in situ size distribution samples indicate that this simulated underestimation in both E3SMv2-SPA and CESM2-WACCM is due to overly small accumulation and coarse-mode aerosols 6–18 months post-eruption, with E3SMv2-SPA simulating ∼ 50 % of the coarse-mode geometric mean diameters of observations 11 months post-eruption. Effective radii from the models and observations are used to calculate offline scattering and absorption efficiencies to explore the implications of smaller simulated aerosol size for the Pinatubo climate impacts. Scattering efficiencies at wavelengths of peak solar irradiance (∼ 0.5 µm) are 10 %–80 % higher for daily samples in models relative to observations through 1993, suggesting higher diffuse radiation at the surface and a larger cooling effect in the models due to the smaller simulated aerosol; absorption efficiencies at the peak wavelengths of outgoing terrestrial radiation (∼ 10 µm) are 15 %–40 % lower for daily samples in models relative to observations, suggesting an underestimation in stratospheric heating in the models due to the smaller simulated aerosol. These potential biases are based on aerosol size alone and do not take into account differences in the aerosol number. The overall agreement of E3SMv2-SPA with observations and its similar performance to the well-validated CESM2-WACCM makes E3SMv2-SPA a viable alternative to simulating climate impacts from stratospheric sulfate aerosols.

Brown, Hunter York (ORCID:0000000218147874)↗
DOE OSTI2024

The ab initio non-crystalline structure database: empowering machine learning to decode diffusivity

Non-crystalline materials exhibit unique properties that make them suitable for various applications in science and technology, ranging from optical and electronic devices and solid-state batteries to protective coatings. However, data-driven exploration and design of non-crystalline materials is hampered by the absence of a comprehensive database covering a broad chemical space. In this work, we present the largest computed non-crystalline structure database to date, generated from systematic and accurate ab initio molecular dynamics (AIMD) calculations. We also show how the database can be used in simple machine-learning models to connect properties to composition and structure, here specifically targeting ionic conductivity. These models predict the Li-ion diffusivity with speed and accuracy, offering a cost-effective alternative to expensive density functional theory (DFT) calculations. Furthermore, the process of computational quenching non-crystalline structures provides a unique sampling of out-of-equilibrium structures, energies, and force landscape, and we anticipate that the corresponding trajectories will inform future work in universal machine learning potentials, impacting design beyond that of non-crystalline materials. In addition, combining diffusion trajectories from our dataset with models that predict liquidus viscosity and melting temperature could be utilized to develop models for predicting glass-forming ability.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fast estimation of reaction rates in spherical and non-spherical porous catalysts

Here we present a methodology for modeling multi-step reaction rates in porous catalyst particles for use in CFD-DEM and two fluid models. Single-step effectiveness factors based on a Thiele modulus, while useful, cannot accurately capture the cascading reaction systems common in high temperature vapor-phase chemical reactors like fluidized catalytic cracking units and catalytic biomass fast pyrolysis systems. Instead, multi-step effectiveness vectors derived from steady-state solutions to the governing reaction-diffusion equations are needed. Solutions for various catalyst shapes are presented, including spheres, cylinders, and prisms. Computational challenges inherent in repeated evaluation of reaction rates with diffusion limitations are discussed, and an efficient implementation based on pre-computed lookup tables is proposed and demonstrated on a simulation of a fluidized bed reactor. Open-source code is provided for the compilation of reaction rate tables for use in ODE, DEM, and two-fluid models.

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

Strengthening of Calcite Assemblages Through Chemical Complexation Reactions

Abstract Utilization of subsurface reservoirs, including fluid extraction/injection, can induce stress changes and modify in‐situ chemical equilibrium causing subcritical fracturing and deformation. Here, we show that chemical complexation reactions at the newly created fracture surfaces within the crack tip limit fracturing during consolidation of granular calcite assemblages. Previously, we showed that analogous chemical complexation reactions can increase in‐situ fracture toughness in pre‐fractured single crystal calcite. To test the chemical complexation effect under subsurface conditions, we consolidated samples under increasing hydrostatic pressure to induce widespread intergranular fracturing in the presence of fluid containing common anions. Measured consolidation, acoustic emissions, and microfracturing demonstrate that deformation correlates with complexation affinity for calcium in calcite and aqueous anions. We observed that deformation is lowest in the presence of strongly complexing anions (e.g., sulphate), provided sufficient anion diffusion rates. This understanding is important for predicting reservoir effective strength in the presence of complex aqueous fluids.

Choens, R. C.↗