Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion maps”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A novel design of transitional layer structure between reduced activation ferritic martensitic steels and tungsten for plasma facing materials

Plasma-facing components (PFCs) are among the most critical gaps for fusion energy to establish technical and economic feasibility. Tungsten as a first wall/blanket material in PFCs requires to be integrating with reduced activation ferritic martensitic (RAFM) steels as a structural component. Currently, major drawbacks are the requirement of brazing, the formation of a brittle interface, and a large difference between the coefficients of thermal expansion of tungsten and steel. Here, a novel transitional multilayer structure was designed and investigated to join tungsten and RAFM steels using three interlayers. The composition of each interlayer was selected based on computational thermodynamics and diffusion kinetics to ensure a body-centered cubic (bcc) single-phase structure and prevent the formation of a brittle intermetallic phase region in the temperature range of 600–1150 °C. Although the transitional layer structure was designed for additive manufacturing, spark plasma sintering (SPS) as proof of concept was used to bond the individual layers. Interfaces were investigated using scanning and transmission electron microscopy methods but no layered intermetallic phase was observed. Nanoindentation maps across the interface suggest major hardness differences at the interface between tungsten and the vanadium interlayer, as well as the interface between RAFM steel and the FeCrAl interlayer.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Homoepitaxial growth of 9,10-diphenylanthracene

Organic single crystals are often shown to have significantly improved optoelectronic properties over polycrystalline thin films such as exciton diffusivity and carrier mobility. There is growing interest in incorporating such crystals in organic electronics despite a number of challenges. In this work, the homoepitaxial vapor phase growth modes for 9,10-diphenylanthracene (DPA) homoepitaxy are mapped as a function of growth rate and temperature using in situ ultra-low current reflection high energy electron diffraction (RHEED), X-ray diffraction (XRD), and atomic force microscopy (AFM). Vapor phase deposition was performed on free-standing single crystal DPA substrates, which show clear Kikuchi patterns. At room temperature, a transition from Frank-Van der Merwe layer-by-layer growth to step-flow growth is observed as the deposition rate is decreased, while at lower temperatures a Stranski-Krastanov layer-plus-island growth is observed. This is the first demonstration of true step-flow growth from room temperature vapor phase deposition for organic semiconductors and is reminiscent of traditional semiconductor growth regimes. Accordingly, these results could lead to improved control over growth and doping of organic single crystals, and lead to enhanced single crystal organic optoelectronic applications.

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

GeSn alloys with ∼21% Sn grown by effusion cell molecular beam epitaxy

We report the epitaxial growth of high quality GeSn alloys with Sn compositions up to 21.25 ± 1% by effusion-cell molecular beam epitaxy (MBE). Achieving such a high Sn content with effusion cells is particularly significant, since these sources typically impose stronger radiative heating on the substrate, which is a factor long considered a major obstacle to high-Sn GeSn epitaxy. Contrary to the prevailing assumption that high growth rates are good for suppressing Sn segregation, we demonstrate that carefully controlled low fluxes, combined with stable ultra-low substrate temperatures, enable significant Sn incorporation even under effusion-cell conditions. Structural analyses by high-resolution x-ray diffraction, reciprocal space mapping, atomic force microscopy, and transmission electron microscopy confirm sharp interfaces, high crystallinity, and smooth surfaces at high Sn contents. Furthermore, these findings establish the feasibility of MBE-grown high-Sn GeSn even under diffusion cells, thereby broadening opportunities for Si-compatible photonic and quantum devices.

Alloys↗
DOE OSTI2022

Diffuse back-illumination temperature imaging (DBI-TI), a novel soot thermometry technique

To meet stringent emissions regulations on soot emissions, it is critical to further advance the fundamental understanding of in-cylinder soot formation and oxidation processes. Among several optical techniques for soot quantification, diffuse back-illumination extinction imaging (DBI-EI) has recently gained traction mainly due to its ability to compensate for beam steering, which if not addressed, can cause unacceptably high measurement uncertainty. Until now, DBI-EI has only been used to measure the amount of soot along the line of sight, and here, we extend the capabilities of a DBI-EI setup to also measure in-cylinder soot temperature. This proof of concept of diffuse back-illumination temperature imaging (DBI-TI) as a soot thermometry technique is presented by implementing DBI-TI in a single cylinder, heavy-duty, optical diesel engine to provide 2-D line-of-sight integrated soot temperature maps. The potential of DBI-TI to be an accurate thermometry technique for use in optical engines is analyzed. The achievable accuracy is due in part to simultaneous measurement of the soot extinction, which circumvents the uncertainty in dispersion coefficients that depend on the optical properties of soot and the wavelength of light utilized. Analysis shows that DBI-TI provides temperature estimates that are closer to the mass-averaged soot temperature when compared to other thermometry techniques that are more sensitive to soot temperature closer to the detector. Furthermore, uncertainty analysis and Monte Carlo (MC) simulations provide estimates of the temperature measurement errors associated with this technique. The MC simulations reveal that for the light intensities and optical densities encountered in these experiments, the accuracy of the DBI-TI technique is comparable or even better than other established optical thermometry techniques. Thus, DBI-TI promises to be an easily implementable extension to the existing DBI-EI technique, thereby extending its ability to provide comprehensive line-of-sight integrated information on soot.

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

In-Operando Characterization of P-I-N Perovskite Solar Cell Under Reverse Bias

A solar module is often exposed to a variety of intensity and uniformity illumination conditions throughout a day. In the field, objects or dirt can block light from hitting the module creating a partial shading event. When a partial shading event occurs the module power point tracking system will typically push the shaded cells into reverse bias in order to avoid the current loss from the other series connected illuminated cells. When perovskite solar cells are operated in reverse bias, electrochemical reactions can occur, leading to both reversible and irreversible performance losses. Excess holes generated by reverse bias cause and push the oxidation reaction of halide ions from the lattice to neutral halogen interstitial. Small interstitial halogen atoms can diffuse across the device and potentially move from the perovskite layer and into the transport layers, shifting the fermi levels in the devices. In this investigation we utilize cross-sectional kelvin probe atomic force microscopy to map the band energetics of perovskite solar cells during reverse bias. We observe a 150 meV decrease in the work function of the C 60 . This leads to a decrease of 150 mV and 4 mA/cm 2 in V oc and in J sc respectively. Finally, we discuss methods to prevent this irreversible degradation after a perovskite device is held under reverse bias.

fullerene↗
DOE OSTI2024

The Exometabolome of Xylella fastidiosa in Contact with Paraburkholderia phytofirmans Supernatant Reveals Changes in Nicotinamide, Amino Acids, Biotin, and Plant Hormones

Microbial competition within plant tissues affects invading pathogens’ fitness. Metabolomics is a great tool for studying their biochemical interactions by identifying accumulated metabolites. Xylella fastidiosa, a Gram-negative bacterium causing Pierce’s disease (PD) in grapevines, secretes various virulence factors including cell wall-degrading enzymes, adhesion proteins, and quorum-sensing molecules. These factors, along with outer membrane vesicles, contribute to its pathogenicity. Previous studies demonstrated that co-inoculating X. fastidiosa with the Paraburkholderia phytofirmans strain PsJN suppressed PD symptoms. Here, we further investigated the interaction between the phytopathogen and the endophyte by analyzing the exometabolome of wild-type X. fastidiosa and a diffusible signaling factor (DSF) mutant lacking quorum sensing, cultivated with 20% P. phytofirmans spent media. Liquid chromatography–mass spectrometry (LC-MS) and the Method for Metabolite Annotation and Gene Integration (MAGI) were used to detect and map metabolites to genomes, revealing a total of 121 metabolites, of which 25 were further investigated. These metabolites potentially relate to host adaptation, virulence, and pathogenicity. Notably, this study presents the first comprehensive profile of X. fastidiosa in the presence of a P. phytofirmans spent media. The results highlight that P. phytofirmans and the absence of functional quorum sensing affect the ratios of glutamine to glutamate (Gln:Glu) in X. fastidiosa. Additionally, two compounds with plant metabolism and growth properties, 2-aminoisobutyric acid and gibberellic acid, were downregulated when X. fastidiosa interacted with P. phytofirmans. These findings suggest that P. phytofirmans-mediated disease suppression involves modulation of the exometabolome of X. fastidiosa, impacting plant immunity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Insights into the dynamics between viruses and their hosts in a hot spring microbial mat

Abstract Our current knowledge of host–virus interactions in biofilms is limited to computational predictions based on laboratory experiments with a small number of cultured bacteria. However, natural biofilms are diverse and chiefly composed of uncultured bacteria and archaea with no viral infection patterns and lifestyle predictions described to date. Herein, we predict the first DNA sequence-based host–virus interactions in a natural biofilm. Using single-cell genomics and metagenomics applied to a hot spring mat of the Cone Pool in Mono County, California, we provide insights into virus–host range, lifestyle and distribution across different mat layers. Thirty-four out of 130 single cells contained at least one viral contig (26%), which, together with the metagenome-assembled genomes, resulted in detection of 59 viruses linked to 34 host species. Analysis of single-cell amplification kinetics revealed a lack of active viral replication on the single-cell level. These findings were further supported by mapping metagenomic reads from different mat layers to the obtained host–virus pairs, which indicated a low copy number of viral genomes compared to their hosts. Lastly, the metagenomic data revealed high layer specificity of viruses, suggesting limited diffusion to other mat layers. Taken together, these observations indicate that in low mobility environments with high microbial abundance, lysogeny is the predominant viral lifestyle, in line with the previously proposed “Piggyback-the-Winner” theory.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

An Optical Spectrum of the Diffuse Galactic Light from BOSS and IRIS

We present a spectrum of the diffuse Galactic light (DGL) between 3700 and 10,000 Å, obtained by correlating optical sky intensity with far-infrared dust emission. We use nearly 250,000 blank-sky spectra from BOSS/SDSS-III together with IRIS-reprocessed maps from the IRAS satellite. The larger sample size compared to SDSS-II results in a factor-of-2 increase in signal to noise. We combine these data sets with a model for the optical/far-infrared correlation that accounts for self-absorption by dust. The spectral features of the DGL agree remarkably well with the features present in stellar spectra. There is evidence for a difference in the DGL continuum between the regions covered by BOSS in the northern and southern Galactic hemispheres. We interpret the difference at red wavelengths as the result of a difference in stellar populations, with mainly old stars in both regions, but a higher fraction of young stars in the south. There is also a broad excess in the southern DGL spectrum over the prediction of a simple radiative transfer model, without a clear counterpart in the north. We interpret this excess, centered at ~6500 Å, as evidence for luminescence in the form of extended red emission. The observed strength of the 4000 Å break indicates that at most ~7% of the dust-correlated light at 4000 Å can be due to blue luminescence. Our DGL spectrum provides constraints on dust scattering and luminescence, independent of measurements of extinction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Planck 2018 results XII. Galactic astrophysics using polarized dust emission

Observations of the submillimetre emission from Galactic dust, in both total intensity I and polarization, have received tremendous interest thanks to the Planck full-sky maps. In this paper we make use of such full-sky maps of dust polarized emission produced from the third public release of Planck data. As the basis for expanding on astrophysical studies of the polarized thermal emission from Galactic dust, we present full-sky maps of the dust polarization fraction p, polarization angle ψ, and dispersion function of polarization angles S . The joint distribution (one-point statistics) of p and NH confirms that the mean and maximum polarization fractions decrease with increasing NH. The uncertainty on the maximum observed polarization fraction, pmax = 22.0 -1.4 +3.5% at 353 GHz and 80' resolution, is dominated by the uncertainty on the Galactic emission zero level in total intensity, in particular towards diffuse lines of sight at high Galactic latitudes. Furthermore, the inverse behaviour between p and S found earlier is seen to be present at high latitudes. This follows the S ∝ p -1 relationship expected from models of the polarized sky (including numerical simulations of magnetohydrodynamical turbulence) that include effects from only the topology of the turbulent magnetic field, but otherwise have uniform alignment and dust properties. Thus, the statistical properties of p, ψ, and S for the most part reflect the structure of the Galactic magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map S × p, looking for residual trends. While the polarization fraction p decreases by a factor of 3-4 between NH = 10 20 cm -2 and NH = 2 x 10 22 cm -2 , out of the Galactic plane, this product S × p only decreases by about 25%. Because S is independent of the grain alignment efficiency, this demonstrates that the systematic decrease in p with N H is determined mostly by the magnetic-field structure and not by a drop in grain alignment. This systematic trend is observed both in the diffuse interstellar medium (ISM) and in molecular clouds of the Gould Belt. Second, we look for a dependence of polarization properties on the dust temperature, as we would expect from the radiative alignment torque (RAT) theory. We find no systematic trend of S × p with the dust temperature T d , whether in the diffuse ISM or in the molecular clouds of the Gould Belt. In the diffuse ISM, lines of sight with high polarization fraction p and low polarization angle dispersion S tend, on the contrary, to have colder dust than lines of sight with low p and high S . We also compare the Planck thermal dust polarization with starlight polarization data in the visible at high Galactic latitudes. The agreement in polarization angles is remarkable, and is consistent with what we expect from the noise and the observed dispersion of polarization angles in the visible on the scale of the Planck beam. The two polarization emission-to-extinction ratios, R P/p and R S/V , which primarily characterize dust optical properties, have only a weak dependence on the column density, and converge towards the values previously determined for translucent lines of sight. We also determine an upper limit for the polarization fraction in extinction, p V /E(B - V), of 13% at high Galactic latitude, compatible with the polarization fraction p ≈ 20% observed at 353 GHz. Taken together, these results provide strong constraints for models of Galactic dust in diffuse gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Solid-State Calcium-Ion Diffusion in Ca 1.5 Ba 0.5 Si 5 O 3 N 6

Rechargeable batteries based on multivalent working ions are promising candidates for next-generation high-energy-density batteries. Development of these technologies, however, is largely limited by the low diffusion rate of multivalent ions in solid-state materials, thereby necessitating a better understanding of the design principles that control multivalent-ion mobility. We report Ca 1.5 Ba 0.5 Si 5 O 3 N 6 as a potential calcium solid-state conductor and investigate its Ca migration mechanism by means of ab initio computations and neutron diffraction. This compound contains partially occupied Ca sites in close proximity to each other, providing a unique mechanism for Ca migration. Nuclear density maps obtained with the maximum entropy method from neutron powder diffraction data provide strong evidence for low-energy percolating one-dimensional pathways for Ca-ion migration. Ab initio molecular dynamics simulations further support a low Ca-ion migration barrier of ~400 meV when Ca vacancies are present and reveal a unique "vacancy-adjacent"concerted ion migration mechanism. This work provides a new understanding of solid-state Ca-ion diffusion and insights into the future design of novel cation configurations that utilize the interactions between mobile ions to enable fast multivalent-ion conduction in solid-state materials.

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

Mapping Current Collection in Cross Section: The Case of Copper-Doped CdTe Solar Cells

For decades, copper has been used to improve the performance of cadmium telluride thin film solar cells. However, it has also been shown to be the main cause of metastability in CdTe. Recently a low activation energy has been reported for the thermal diffusion of Cu in CdTe explaining the ease of motion that it has under moderate temperatures. The community consensus is that copper segregates to the absorber grain boundaries, where it's either beneficial or detrimental to device performance depending on its concentration. Using nanoscale X-ray micrsocopy and a two-dimensional drift-diffusion model we present a preliminary correlation between local copper distribution and electrical performance of a single-junction CdTe/CdS solar cell.

CdTe↗
DOE OSTI2023

Colloidal State Machines as Smart Tracers for Chemical Reactor Analysis

A widely utilized tool in reactor analysis is passive tracers that report the residence time distribution, allowing estimation of the conversion and other properties of the system. Recently, advances in microrobotics have introduced powered and functional entities with sizes comparable to some traditional tracers. This has motivated the concept of Smart Tracers that could record the local chemical concentrations, temperature, or other conditions as they progress through reactors. Herein, the design constraints and advantages of Smart Tracers by simulating their operation in a laminar flow reactor model conducting chemical reactions of various orders are analyzed. It is noted that far fewer particles are necessary to completely map even the most complex concentration gradients compared with their conventional counterparts. Design criteria explored herein include sampling frequency, memory storage capacity, and ensemble number necessary to achieve the required accuracy to inform a reactor model. Cases of severe particle diffusion and sensor noise appear to bind the functional upper limit of such probes and require consideration for future design. The results of the study provide a starting framework for applying the new technology of microrobotics to the broad and impactful set of problems classified as chemical reactor analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Bayesian reduced-order deep learning surrogate model for dynamic systems described by partial differential equations

We propose a reduced-order deep-learning surrogate model for dynamic systems described by time-dependent partial differential equations. This method employs space–time Karhunen–Loève expansions (KLEs) of the state variables and space-dependent KLEs of space-varying parameters to identify the reduced (latent) dimensions. Subsequently, a deep neural network (DNN) is used to map the parameter latent space to the state variable latent space. An approximate Bayesian method is developed for uncertainty quantification (UQ) in the proposed KL-DNN surrogate model. The KL-DNN method is tested for the linear advection–diffusion and nonlinear diffusion equations, and the Bayesian approach for UQ is compared with the deep ensembling (DE) approach, commonly used for quantifying uncertainty in DNN models. It was found that the approximate Bayesian method provides a more informative distribution of the PDE solutions in terms of the coverage of the reference PDE solutions (the percentage of nodes where the reference solution is within the confidence interval predicted by the UQ methods) and log predictive probability. The DE method is found to underestimate uncertainty and introduce bias. For the nonlinear diffusion equation, we compare the KL-DNN method with the Fourier Neural Operator (FNO) method and find that KL-DNN is 10% more accurate and needs less training time than the FNO method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

3D hybrid fluid-particle jet simulations and the importance of synchrotron radiative losses

Context. Relativistic jets in active galactic nuclei are known for their exceptional energy output, and imaging the synthetic synchrotron emission of numerical jet simulations is essential for a comparison with observed jet polarization emission. Aims. Through the use of 3D hybrid fluid-particle jet simulations (with the PLUTO code), we overcome some of the commonly made assumptions in relativistic magnetohydrodynamic (RMHD) simulations by using non-thermal particle attributes to account for the resulting synchrotron radiation. Polarized radiative transfer and ray-tracing (via the RADMC-3D code) highlight the differences in total intensity maps when (i) the jet is simulated purely with the RMHD approach, (ii) a jet tracer is considered in the RMHD approach, and (iii) a hybrid fluid-particle approach is used. The resulting emission maps were compared to the example of the radio galaxy Centaurus A. Methods. We applied the Lagrangian particle module implemented in the latest version of the PLUTO code. This new module contains a state-of-the-art algorithm for modeling diffusive shock acceleration and for accounting for radiative losses in RMHD jet simulations. The module implements the physical postulates missing in RMHD jet simulations by accounting for a cooled ambient medium and strengthening the central jet emission. Results. We find a distinction between the innermost structure of the jet and the back-flowing material by mimicking the radio emission of the Seyfert II radio galaxy Centaurus A when considering an edge-brightened jet with an underlying purely toroidal magnetic field. We demonstrate the necessity of synchrotron cooling as well as the improvements gained when directly accounting for non-thermal synchrotron radiation via non-thermal particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Electroluminescence Characterization of Recombination in Back Junction Silicon Heterojunction Test Structures: Role of the Inversion Layer

Electroluminescence (EL) characterization technique is used to quantify recombination properties in back junction silicon (Si) heterojunction (HJ) solar cell test structures with interdigitated back contact (IBC) emitters and metal contacts. These test structures are three terminal devices that allow the selective biasing and injection of minority carriers into spatially isolated emitter strips having different widths. This allows the study of minority carriers' current transport, flow, and recombination in these regions, which can, then, be applied to full IBC cells. EL is used to verify that emission is proportional to the minority carriers' concentration in the corresponding regions of the device, thus providing information relevant to minority carriers' collection when compared with laser beam induced current technique in considerably less time and with a higher spatial resolution. Current transport through the inversion layer that occurs due to the band bending at the n crystalline Si and p amorphous Si interface is discussed relying on the distance that minority carriers travel in different configurations and on their diffusion length. Laser defect spots (LDS) were, then, introduced in the emitter regions to analyze the spatial profile of minority carriers as they recombine at the LDS due to their high surface recombination velocity. Ideality factors in different configurations are spatially mapped from EL before and after introducing LDS and are used to show the effect of the inversion layer on the different recombination mechanisms. Current transport through the inversion layer is studied at different biasing conditions and correlated with current transport after LDS.

42 ENGINEERING↗
DOE OSTI2022

Impact of Concurrent Solubilization and Fines Migration on Fracture Aperture Growth in Shales during Acidized Brine Injection

Acidic hydraulic fracturing fluid chemically and physically alters shale rock fabric during injection and shut-in, creating a “reaction-altered zone” along the fracture faces. To better characterize the variable thickness and composition of this reaction-altered zone under advective flow, we take a coupled experimental and modeling approach. A fluidic cell, with six fiducial markers, is first fabricated to keep the rock sample in place during the core floods and to allow image alignment of acquired images. Then, we conduct a series of reactive core floods in a clay-rich siliceous Wolfcamp shale sample with 10 wt % carbonate, using a synthetic fracturing fluid under no confining stress and at room temperature. High-resolution computed tomography (CT) scans are periodically conducted to observe the spatial alteration of the fracture network. We then perform scanning electron microscopy (SEM-EDS) on the two orthogonal surfaces (fracture surface and freshly cut profile face) to generate high-resolution elemental maps that show the change in mineralogy, both with distance along a given flow path along a fracture surface and with depth from the fracture surface into the shale matrix. These results are contrasted against a two-dimensional (2D) advection-diffusion reaction model developed previously for batch reactions between shale and synthetic fracturing fluids. The model simulates the geochemical interaction occurring at the fracture/matrix interface and penetrating into the shale matrix during the reactive core flood. Both model and experimental results show that the acidic brine is neutralized during the core flood, corresponding to an increase in fracture aperture as a function of fluid volume injected with the greatest change near the inlet. SEM-EDS scans reveal significant dissolution of carbonates on the fracture surface without pyrite oxidation. The reactive transport model indicates that carbonate depletion into the shale interior should be observable, yet SEM-EDS shows no discernible loss of carbonate in the orthogonal profile face. The combination of these observations suggests an additional fracture evolution mechanism in the reactive system, i.e., fines migration. We show that fines migration enhances the access of fracturing fluid to the matrix resulting in a more pronounced fracture widening. Finally, we conclude that coupled mineral dissolution and fines migration govern fracture aperture growth during acidized brine injection. In this work, we effectively show the underlying risk of relying solely on models that do not include an important (transport) process that can alter the system significantly and propose a combined chemomechanical mechanism for fracture evolution appropriate for this shale mineralogy.

58 GEOSCIENCES↗
DOE OSTI2024

Fundamentals of Scalable Nanoparticle Assembly in Engineering Polymres

Systematically improving the properties of polymer nanocomposites is contingent on effectively controlling nanoparticle (NP) dispersion; further, it is necessary to attain this at scale if these methods are to impact practical applications. While most past work has focused on amorphous polymers, ≈ 70% of all polymers sold annually are semicrystalline. Inspired by this apparent disconnect, our work in the current DOE funded grant and this renewal proposal employs a novel handle, polymer crystallization, to organize NPs into the various amorphous zones of the lamellar semicrystalline morphology (interlamellar, interfibrillar and interspherulitic domains) of stereoregular polymers. We hypothesize that the isothermal crystallization rate of a polymer host, G, and NP diffusivity, are the two relevant control parameters in this context. In the case where the NPs are miscible in the melt, they remain spatially well-dispersed when the polymer is crystallized at a rate that is fast relative to NP diffusion. In contrast, monolayer-thick NP sheets with intersheet spacing of ~10-100nm form in slowly crystallized samples. This sheet stacking is remarkably similar to that seen in natural materials, e.g., Nacre. There is also ordering of the NPs at the larger interfibrillar scale. Importantly, the relative fraction of engulfed, interlamellar and interfibrillar NPs is determined by G. Further, the modulus of these materials can increase by a factor of 2-3 due to the organization state of the NPs, while leaving fracture toughness unaffected. We have established the broad applicability of this approach in two polymers, polyethylene (PE) and polyethyleneoxide (PEO). To achieve this goal we developed a method for grafting long PE chains from NPs with controlled graft density and molecular weight, giving us the ability to tailor their dispersion in the PE melt prior to crystallization. In this renewal proposal, we will develop the fundamental, quantitative understanding that will allow us to independently tailor polymer crystallization rate, NP diffusivity and NP agglomeration state in semicrystalline polymers so as to facilely control NP ordering and hence properties. A program that tightly integrates experiment and theory is proposed, and a key outcome of this research will be a three-dimensional design map for different regimes of NP organization (and hence properties) that can be achieved by polymer crystallization. A second independent outcome is the role of directional solidification and how it competes with the rates of polymer crystallization and NP mobility to affect the organization of the NP.

36 MATERIALS SCIENCE↗