Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density matrix”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Plasma production and ion heating systems for the Material Plasma Exposure eXperiment

Plasma production and ion heating in the Material Plasma Exposure eXperiment (MPEX), whose design is nearing completion, is accomplished using continuous wave RF power with average power density up to 1.6 MW/m2 at the plasma interface. Plasma is produced using helicon waves coupled through a single helical antenna at 13.56 MHz with power supplied by three 100 kW fixed-frequency RF generators feeding a power combiner network followed by the matching network and launcher. The helicon source can utilize various gasses including hydrogen, deuterium, and helium, with a magnetic field strength in the source region up to 0.2 T, and maximum device |B| of 2.5 T. Power is coupled to ions via ion cyclotron heating at the fundamental resonance using a pair of phased helical antennas operating in the frequency range 4-9 MHz, that launch waves towards the resonance from the antenna region where ω > ωci. ICH power is supplied by a single 500 kW tunable RF transmitter through a 90° power splitter and matching/decoupling network. In the case of both the helicon and ICH systems the antennas are located external to the vacuum, with power transferred through novel water-cooled coaxial vacuum windows consisting of fused quartz outer cylinders and silicon nitride inner cylinders with forced convection water cooling between them. The antenna enclosures are pressurized with dry air to 3 bar absolute for voltage standoff.Several 3-D COMSOL models have been created to simulate the two systems. A model of the helicon region utilizing a cold plasma dielectric tensor with accurate magnetic field and realistic plasma density profiles has been used to calculate the plasma loading/complex antenna input impedance at the launcher feed for various ne and |B| values, for the purpose of estimating power handling. It also incorporates the geometry of all launcher structures relevant to this determination. A still more detailed model of the launcher together with a lossy dielectric plasma surrogate has been used to determine RF electric field values and power losses in the device components.Similar models have been produced to predict the performance and power handling of the ICH launcher. For this device the impedance matrix of the two-element antenna array is calculated using a warm plasma model, necessary to properly determine the wave propagation and absorption.

Goulding, Richard↗
DOE OSTI2021

Isolating the Electrocatalytic Activity of a Confined NiFe Motif within Zirconium Phosphate

Unique classes of active-site motifs are needed for improved electrocatalysis. Herein, the activity of a new catalyst motif is engineered and isolated for the oxygen evolution reaction (OER) created by nickel–iron transition metal electrocatalysts confined within a layered zirconium phosphate matrix. It is found that with optimal intercalation, confined NiFe catalysts have an order of magnitude improved mass activity compared to more conventional surface-adsorbed systems in 0.1 m KOH. Interestingly, the confined environments within the layered structure also stabilize Fe-rich compositions (90%) with exceptional mass activity compared to known Fe-rich OER catalysts. Through controls and by grafting inert molecules to the outer surface, it is evidenced that the intercalated Ni/Fe species stay within the interlayer during catalysis and serve as the active site. After determining a possible structure (wycherproofite), density functional theory is shown to correlate with the observed experimental compositional trends. It is further demonstrated that the improved activity of this motif is correlated to the Fe and water content/composition within the confined space. This work highlights the catalytic enhancement possibilities available through zirconium phosphate and isolates the activity from the intercalated species versus surface/edge ones, thus opening new avenues to develop and understand catalysts within unique nanoscale chemical environments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Magnus method for electronic structure calculations at extreme conditions

We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

general physics↗
DOE OSTI2021

Additive manufacturing of cobalt-based alloy on tool steel by directed energy deposition

Cladding hard-surfacing alloys on tool steel is an effective approach to enhance the surface properties of tool steel. In this study, a Co-based alloy was deposited on tool steel by Directed Energy Deposition (DED) following a three-factor three-level design of experiment matrix with varied laser power, scan speed, and powder flow rate. The microstructure of the deposits was characterized using scanning electron microscopy (SEM). The residual stress on the surface of the samples was measured by the X-ray diffraction (XRD) sin2? technique. The parameters that produced promising deposits were used to fabricate samples for tensile test, four-point bending test, Charpy impact test, and hardness measurement. Our result reveals that the processing parameters have a significant role in the residual stress of the coatings. Residual stress reduces with the increase of laser energy density. Cracks were found at samples with energy density below a threshold. Tensile testing of the coating/substrate combined structure reveals fracture at the coatings with an ultimate tensile strength of 633.9 ± 54.7 MPa. The bi-material interface survived the tensile test, indicating a strong interfacial bond. The four-point bending test of coating/substrate laminates shows an ultimate flexure strength of 860.6 ± 36.9 MPa. Cracks initiated from the coatings ignored the interface and penetrated the substrate, suggesting a solid bi-material bond. Charpy impact test shows the absorbed energy of coating/substrate laminates is more than doubled that of the substrate.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Extraction and validation of a new set of CMS pythia8 tunes from underlying-event measurements

New sets of CMS underlying-event parameters ("tunes") are presented for the PYTHIA8 event generator. These tunes use the NNPDF3.1 parton distribution functions (PDFs) at leading (LO), next-to-leading (NLO), or next-to-next-to-leading (NNLO) orders in perturbative quantum chromodynamics, and the strong coupling evolution at LO or NLO. Measurements of charged-particle multiplicity and transverse momentum densities at various hadron collision energies are fit simultaneously to determine the parameters of the tunes. Comparisons of the predictions of the new tunes are provided for observables sensitive to the event shapes at LEP, global underlying event, soft multiparton interactions, and double-parton scattering contributions. In addition, comparisons are made for observables measured in various specific processes, such as multijet, Drell-Yan, and top quark-antiquark pair production including jet substructure observables. The simulation of the underlying event provided by the new tunes is interfaced to a higher-order matrix-element calculation. For the first time, predictions from PYTHIA8 obtained with tunes based on NLO or NNLO PDFs are shown to reliably describe minimum-bias and underlying-event data with a similar level of agreement to predictions from tunes using LO PDF sets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Analysis of diagonal G and subspace W approximations within fully self-consistent GW calculations for bulk semiconducting systems

Fully self-consistent GW (sc-GW) methods are now available to evaluate quasiparticle and spectral properties of various molecular and bulk systems. However, such techniques based on the full matrix of G and W are computationally demanding. Additionally, the routinely used single-shot GW approximation (G 0 W 0 ) has an undesirable dependency on the choice of initial exchange-correlation functional. In the literature, many so-called self-consistent GW methods are based on diagonal approximation of G and low-ranking approximation of W. It is thus worth checking how good such approximations are in comparison with the full matrix method. In this work, we consider AlAs, AlP, GaP, and ZnS as the prototype systems to perform sc-GW calculations by expressing the full G matrix using a plane-wave basis set. We compared our sc-GW results with the diagonal G and subspace W approximated sc-GW results (sc-GW-diagG and sc-GW-subW methods). In the sc-GW-diagG method, interacting G is expanded in the eigenvectors of noninteracting G such that only diagonal elements are retained, whereas the number of eigenmodes is truncated in sc-GW-subW calculations. A systematic analysis of the results obtained from the above techniques is presented. The differences in the quasiparticle band gap between the approximated and the full matrix sc-GW approaches are mostly less than 1.7%, which validates such widely adopted approximations, and also shows how such low-ranking approximation can be used to include higher-order terms such as the vertex correction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Porous and Amorphous Mn x Mo 3 S 13 Chalcogel Electrode for High-Capacity Conversion-Based Lithium-Ion Batteries

While Li-ion batteries (LIBs) are a leading energy storage technology, their energy densities are limited by the low capacity of conventional intercalation cathodes, driving interest in high energy-density Li–S batteries that make use of conversion chemistry. Achieving high capacity, reversibility, and cycle stability, and controlling volume changes in conversion batteries during the charge–discharge process, however, remains challenging. Here, in this study, we present a porous, amorphous, sulfide-based Mn x Mo 3 S 13 chalcogel, which concurrently offers high capacity and cycle stability. The solution-processable room temperature synthesized Mn x Mo 3 S 13 (x = 0.25) chalcogel exhibits a local structure that resembles the Mo 3 S 13 cluster with Mn 2+ distributed across the Mo 3 S 13 matrix, as determined by synchrotron X-ray pair distribution function (PDF) and extended X-ray absorption fine structure (EXAFS). Ab initio molecular dynamics (AIMD) simulations reveal that Mn 2+ incorporation shortens the polysulfide chain in the gel matrix compared to the Mo 3 S 13 chalcogel, while forming a coordination environment with disulfide groups, analogous to the experimental findings. A Li/Mn 0.25 Mo 3 S 13 half-cell delivers 897 mAh g –1 capacity during the first discharge and retains 571 mAh g –1 capacity after 100 cycles at a C/3 rate. Distribution of relaxation time (DRT) unveils a stable solid–electrolyte interphase (SEI) formation upon cycling that enables charge–discharge reversibility. Here, the enhanced capacity retention and cycle stability compared to those of the Li/Mo 3 S 13 cell are attributed to the reduced dissolution of active mass into the electrolyte, facilitated by the formation of shorter polysulfide chains within the Mn 0.25 Mo 3 S 13 structure and the strong affinity of Lewis-acidic Mn 2+ for polysulfide anions generated during the charge–discharge process of the Li/Mn 0.25 Mo 3 S 13 cell. Thus, this work illustrates a design principle of material for high-capacity and cycle-stable Li-metal sulfide batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Influence of gadolinium doping on structural, optical, and electronic properties of polymeric graphitic carbon nitride

Polymeric graphitic carbon nitride (gCN) materials have received great attention in the fields of photo and electrocatalysis due to their distinct properties in metal-free systems with high physicochemical stability. Nevertheless, the activity of undoped gCN is limited due to its relatively low specific surface area, low conductivity, and poor dispersibility. Doping Gd atoms in a gCN matrix is an efficient strategy to fine-tune its catalytic activity and its electronic structure. Herein, the influence of various wt% of gadolinium (Gd) doped in melon-type carbon nitride was systematically investigated. Gadolinium-doped graphitic carbon nitride (GdgCN) was synthesized by adding gadolinium nitrate to dicyandiamide during polymerization. The X-ray diffraction (XRD) and transmission electron microscopy (TEM) results revealed that the crystallinity and the morphological properties are influenced by the % of Gd doping. Furthermore, X-ray photoelectron spectroscopy (XPS) studies revealed that the gadolinium ions bonded with nitrogen atoms. Complementary density functional theory (DFT) calculations illustrate possible bonding configurations of Gd ions both in bulk material and on ultrathin melon layers and provide evidence for the corresponding bandgap modifications induced by gadolinium doping.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Quantitative relationships between film morphology, charge carrier dynamics, and photovoltaic performance in bulk-heterojunction binary vs. ternary acceptor blends

Addressing pertinent and perplexing questions regarding why nonfullerene acceptors (NFAs) promote higher power conversion efficiencies (PCEs) than traditional fullerenes and how photoactive bulk heterojunction (BHJ) film morphology, charge photogeneration, and recombination dynamics dictate solar cell performance have stimulated many studies of polymer solar cells (PSCs), yet quantitative relationships remain limited. Better understanding in these areas offers the potential to advance materials design and device engineering, afford higher PCEs, and ultimate commercialization. Here we probe quantitative relationships between BHJ film morphology, charge carrier dynamics, and photovoltaic performance in model binary and ternary blend systems having a wide bandgap donor polymer, a fullerene, and a promising NFA. We show that optimal PC 71 BM incorporation in a PBDB-TF:ITIC-Th binary system matrix retains the original π-face-on orientation, ITIC-Th crystallinity and BHJ film crystallite dimensions, and reduces film upper surface ITIC-Th segregation. Such morphology changes together simultaneously increase hole (μ h ) and electron (μ e ) mobilities, facilitate light-activated ITIC-Th to PC 71 BM domain electron delocalization, reduce free charge carrier (FC) bimolecular recombination (BR) within PBDB-TF:ITIC-Th mixed regions, and increase FC extraction pathways viaPBDB-TF:PC 71 BM pairs. The interplay of these effects yields significantly enhanced inverted cell short-circuit current density (J SC ), fill factor (FF), and PCE. Unexpectedly, we also find that excessive PC 71 BM in the PBDB-TF:ITIC-Th binary system alters the PBDB-TF orientation to π-edge-on, increases large scale PC 71 BM-rich aggregations and BHJ upper surface PC 71 BM composition. Further, these morphology changes increase parasitic decay processes such as intersystem crossing from photoexcited PC 71 BM, compromising the J SC , FF, and PCE metrics. ITIC-Th X-ray diffraction reveals a unique sidechain-dominated molecular network with previously unknown sidechain-end group stacking, rationalizing the STEM and GIWAXS results, photophysics, and the high μ e . DFT computation reveals charge transfer networks within ITIC-Th crystallites, supporting excited-state electron delocalization from ITIC-Th to PC 71 BM. This structure–property understanding leads to a newly reported NFA blend with PCE near 17%.

14 SOLAR ENERGY↗
DOE OSTI2020

Capturing non-Gaussianity of the large-scale structure with weighted skew-spectra

The forthcoming generation of wide-field galaxy surveys will probe larger volumes and galaxy densities, thus allowing for a much larger signal-to-noise ratio for higher-order clustering statistics, in particular the galaxy bispectrum. Extracting this information, however, is more challenging than using the power spectrum due to more complex theoretical modeling, as well as significant computational cost of evaluating the bispectrum signal and the error budget. To overcome these challenges, several proxy statistics have been proposed in the literature, which partially or fully capture the information in the bispectrum, while being computationally less expensive than the bispectrum. One such statistics are weighted skew-spectra, which are cross-spectra of the density field and appropriately weighted quadratic fields. Using Fisher forecasts, we show that the information in these skew-spectra is equivalent to that in the bispectrum for parameters that appear as amplitudes in the bispectrum model, such as galaxy bias parameters or the amplitude of primordial non-Gaussianity. We consider three shapes of the primordial bispectrum: local, equilateral and that due to massive particles with spin two during inflation. In conclusion, to obtain constraints that match those from a measurement of the full bispectrum, we find that it is crucial to account for the full covariance matrix of the skew-spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Towards a self-consistent analysis of the anisotropic galaxy two- and three-point correlation functions on large scales: application to mock galaxy catalogues

We establish a practical method for the joint analysis of anisotropic galaxy two- and three-point correlation functions (2PCF and 3PCF, respectively) on the basis of the decomposition formalism of the 3PCF using tripolar spherical harmonics. We perform such an analysis with MultiDark-Patchy mock catalogues to demonstrate and understand the benefit of the anisotropic 3PCF. We focus on scales above $80\, h^{-1}\, {\rm Mpc}$, and use information from the shape and the baryon acoustic oscillation (BAO) signals of the 2PCF and 3PCF. We also apply density field reconstruction to increase the signal-to-noise ratio of BAO in the 2PCF measurement, but not in the 3PCF measurement. In particular, we study in detail the constraints on the angular diameter distance and the Hubble parameter. We build a model of the bispectrum or 3PCF that includes the non-linear damping of the BAO signal in redshift space. We carefully account for various uncertainties in our analysis including theoretical models of the 3PCF, window function corrections, biases in estimated parameters from the fiducial values, the number of mock realizations to estimate the covariance matrix, and bin size. The joint analysis of the 2PCF and 3PCF monopole and quadrupole components shows a $30{{\ \rm per\ cent}}$ and $20{{\ \rm per\ cent}}$ improvement in Hubble parameter constraints before and after reconstruction of the 2PCF measurements, respectively, compared to the 2PCF analysis alone. Furthermore, this study clearly shows that the anisotropic 3PCF increases cosmological information from galaxy surveys and encourages further development of the modelling of the 3PCF on smaller scales than we consider.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Electrical and structural characterization of nano-carbon–aluminum composites fabricated by electro-charging-assisted process

Incorporating nano-carbon phases into metal-matrix composites is a promising strategy for simultaneously enhancing electrical conductivity and mechanical properties of metals. Here, we describe the manufacture of novel nano-carbonaluminum composites by an electro-charge-assisted process (EAP) that show 5.6% ± 1.7% increase in electrical conductivity compared to the base metal alloy. The source of nano-carbon that we used in this study is activated carbon with particle size less than 100 nm. The enhancement is attributed to nano-graphitic structures that extend through the lattice of the metal. Through electron transfer from the metal to the nanostructures the electron density at the interface of nano-crystalline graphite and the metal lattice increases thereby enhancing the bulk electrical conductivity. We identify the important fabrication parameters of the EAP for a reaction system employing a tapered graphite cathode. Here, a high current density of 100 A/cm 2 causes ionization and crystallization of the carbon in the liquid metal. The increase in electrical conductivity of the composite are directly related to the incorporation of the nanocrystalline carbon in the metal lattice. The superior performance of these nano-carbon aluminum composites makes them promising candidates for power transmission lines and other applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Spontaneous Raman–LIF–CO–OH measurements of species concentration in turbulent spray flames

This paper presents new measurements of species concentrations, temperature and mixture fraction in selected regions of a turbulent ethanol spray flame. The line-Raman–LIF–COsingle bondOH setup developed at the Sandia's Combustion Research Facility is utilised to probe regions of a spray flame where laser breakdown of liquid droplets is avoided and the remaining interferences can be corrected. The spray flame is stabilised on the piloted Sydney needle spray burner, where axial translation of the liquid injecting needle in the air-blast stream can transition the spray from dilute to dense. The solution to obtaining successful measurements is found to be multifaceted and includes: the appropriate selection of flame conditions; high sensitivity of the Raman detection system permitting reduced laser energies; development of a pre-processing algorithm to reject strong droplet interferences; and application of the hybrid matrix inversion method combined with wavelet denoising to account for interference corrections and noise at the very low signal levels obtained. Unique and necessary for the successful measurements reported in this paper, a pre-processing algorithm is outlined that removes data points corrupted with strong interferences from droplets. These interferences arise from a range of sources, but the most intense are due to the laser interaction with surrounding mist or liquid fragments, such that measurements near the jet centreline are corrupted and hence discarded. Reliable measurements of mixture fraction, temperature obtained from the sum of the species number densities, and species mole fractions are reported for regions in the flames sufficiently far from the centreline. The paper demonstrates the feasibility of the judicious use of Raman scattering in turbulent spray flames, the results of which will be extremely useful for validating numerical simulations.

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

In situ transmission electron microscopy observations of CaCO 3 crystallization onto polysaccharide-coated nanoparticles

Polysaccharides and proteoglycans are widely associated with the organic matrix at sites of CaCO 3 biomineralization, and previous studies indicate that these macromolecules may confer greater roles in mineral nucleation than previously recognized. This investigation uses in situ liquid-phase transmission electron microscopy (LP-TEM) to observe CaCO 3 nucleation onto aminated silica (SiO 2 –NH 3 + ) nanoparticles treated with a layer of chitosan (near-neutral derivative of chitin) or heparin (a carboxylated and highly sulfated glycosaminoglycan). In the absence of polysaccharides, few CaCO 3 particles formed and exhibited mobility. However, the SiO 2 –NH 3 + nanoparticles were enveloped in a region of higher mass density relative to the bulk solution, suggesting the development of a local solute-rich environment that surrounds the charged NH 3 + groups. The heparin- or chitosan-coated silica particles also exhibited regions of higher mass density around the nanoparticles. In the presence of these polysaccharide coatings, we observed the nucleation of abundant CaCO 3 particles whereby the polyanionic heparin promoted more nucleation than the weakly cationic chitosan. Many crystallites appeared to form at the polysaccharide–TEM cell membrane–solution interface, further indicating interfacial and macromolecule-specific control on crystallization. The combined results demonstrate that chitosan and heparin have an appreciable effect on the timing, size, and location of CaCO 3 nucleation compared to the polysaccharide-free nanoparticles.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Chirality reversal at finite magnetic impurity strength and local signatures of a topological phase transition

Here, we study the honeycomb lattice with a single magnetic impurity modeled by adding imaginary next-nearest-neighbor hopping 𝑖⁢ℎ on a single hexagon. This Haldane defect gives a topological mass term to the gapless Dirac cones and generates chirality. For a small density of defects, Neehus et al. [Phys. Rev. Lett. 135, 126604 (2025)] found that the system's chirality reverses at a critical ℎ 𝑐 ≈ 0.95 associated with an unexpected tricritical point of Dirac fermions at zero defect density. We investigate this zero-density limit by analyzing a single defect and computing two experimentally relevant measures of chirality: (1) orbital magnetization via local Chern marker, a bulk probe of all occupied states; and (2) electronic currents of low-energy states. Both probes show a chirality reversal at a critical ℎ 𝑐 ≈ 0.9–1.0. Motivated by this consistency, we propose a defect-scale toy model whose low-energy states reverse their chirality at ℎ$^{'}_{c}$ ≈ 0.87. Remarkably, the same pair of zero-energy bound states also generates the critical point ℎ 𝑐 in the full impurity projected T-matrix. Our results show how the chirality reversal produced by an impurity can be observed either in local probes or in the global topology, and suggest a possible role of the microscopic defect structure at the critical point.

Chern insulators↗
DOE OSTI2022

Helium interaction with solutes and impurities in neutron-irradiated nanostructured ferritic alloys: A first principles study

Density functional theory calculations are performed to explore the binding between He and alloying solutes, impurities, and transmutation products expected in neutron irradiated nanostructured ferritic alloys (NFAs), here 14YWT is taken as an example. Elements that exhibit significant binding (attraction) with an interstitial He are Y (binding energy = 0.46 eV), Mg (0.32), O (0.33), Ti (0.16), and C (0.15). Those that provide significant binding to a substitutional He are O (1.44), Y (1.24), N (0.73), H (0.56), Mg (0.52), Ti (0.34), Si (0.34), C (0.33), Al (0.32), Ni (0.26), Ta (0.23), and Mn (0.16). The presence of these elements in Fe matrix could reduce the transport of He towards oxide particles, dislocations, and internal boundaries, and could promote He bubble nucleation in the matrix. For convenience, we compile existing binding energy data of He with He n and He n V (He-vacancy) clusters. Dissociation pathway analysis reveals that, in general, the most likely dissociation of a He n V cluster is by a sequential emission of individual He atoms. Furthermore, larger bubbles are more prone to dissociation than smaller ones. In addition, we estimate the binding energy (segregation energy) of He in bulk Y 2 Ti 2 O 7 (YTO) single crystal, YTO/Fe interface, and YTO particle embedded in Fe, with respect to interstitial He in Fe, from existing formation energies of He in these structures. We also compile available data of He binding with Fe self-interstitial atom (SIA), SIA clusters, and edge and screw dislocations. Note that given the absence of DFT data, the binding with SIA clusters and dislocations are gathered from simulations with empirical potentials. Finally, the data presented in this paper is important to inform multiscale simulations of He bubble accumulation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling of Magnetic Properties of Magnetorheological Elastomers Using JA Hysteresis Model

Magnetorheological elastomers (MREs) are composite materials that consist of magnetically permeable particles in a nonmagnetic polymeric matrix. Under the influence of an external magnetic field, a reversible deformation change occurs in the mechanical properties of these materials. Due to their coupled magnetomechanical response, these materials have been found suitable for a range of applications including tunable vibration absorbers, sensors, and actuators. Notably, improvement of such devices are prerequisites to efficient energy conversion systems, hence the need to understand further the MRE technology. The Jiles-Atherton (JA) theory takes into consideration the magneto-coupling experienced by effective domains in a magnetic material. Algorithm based on the theory yields five model parameters; saturation magnetization (M s ), domain density (a), domain coupling (α), loss coefficient (k), and reversibility (c). Using JA theory, model parameters were calculated and linked to the physical attributes of Fe powder and isotropic MRE. The results show that the calculated Ms for the MRE is reasonably related to that of the Fe powder by a factor of the particle's volume fraction used in the MRE. Further, the calculated k, a, and α provided support for the reduced pinning factor, domain density, and increased domain coupling in the MRE due to the changes in the domain structure between the two materials. From the calculated JA parameters, finite-element modeling (FEM) of the MRE hysteresis loop was performed. The analysis showed that the modeled magnetic properties including coercivity, remanence, and coordinates of the hysteresis loop tip vary with geometric position.

42 ENGINEERING↗
DOE OSTI2021

Reactivity of CO 2 with Utica, Marcellus, Barnett, and Eagle Ford Shales and Impact on Permeability

We report that in order to reduce greenhouse gas emissions while recovering hydrocarbons from unconventional shale formations, processes that make use of carbon dioxide to enhance oil recovery while storing carbon dioxide (CO 2 ) should be considered. Here, we examine samples from three shale basins across the United States (Utica and Marcellus Shales in the Appalachian Basin, Barnett Shale in the Bend Arch-Ft. Worth Basin, and Eagle Ford in the Western Gulf Basin) to address the following questions: (1) do changes from reaction with CO 2 and fluids at the micrometer and nanometer scale alter flow pathways and, in turn, impact hydrocarbon production, CO 2 storage, and seal integrity and (2) can CO 2 or fluid reactivity be predicted based on physical or chemical properties of shale formations? Experiments were conducted at 40 °C and 10.3 MPa to characterize the interaction between CO 2 and shale using X-ray diffraction (XRD), carbon and sulfur analysis, in situ Fourier transform infrared spectroscopy (FT-IR), feature relocation scanning electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), mercury (Hg) intrusion porosimetry, and Brunauer-Emmett-Teller (BET) surface area and pore size analysis coupled with density functional theory (DFT) methods. Changes in mechanical, physical, and flow properties of shale cores due to CO 2 exposure were addressed using a New England Research Autolab 1500 and Xenon X-ray computed tomography (CT) scanning. Results showed that CO 2 did not promote significant reactivity with the shale if water was not present; only shales with swelling clays or residual interstitial pore water reacted with dry CO 2 to promote reactivity in shale. When water was added as a reactant, CO 2 formed carbonic acid and reacted with the shale to dissolve carbonate pockets, etched and pitted the shale matrix surfaces, and increased the microporosity and decreased nanoporosity. Porosity and permeability increased appreciably in core shale samples after exposure to CO 2 saturated fluid due to dissolution of carbonate. Shale mechanical properties were not altered. Trends were not observed that could tie CO 2 or fluid reactivity to physical or chemical properties of the shale formations at the basin scale from the samples we examined. However, if the shale contained significant amounts of carbonate and water was available to react with the CO 2 , pore sizes were altered in the matrix and permeability and porosity increased.

04 OIL SHALES AND TAR SANDS↗