SEARCH · Engineering Papers
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.
A first order geometrical optics treatment of holograms combined with the generation of interference fringes by two point sources is used to describe reference fringe formation in non-diffuse dual-hologram interferometry.
The International Space Station (ISS) is currently equipped with a complex, heavy, and power consuming system that recovers approximately 50% of O2 from metabolic CO2. Future long duration human missions to the Moon and Mars will necessitate a sustainable and highly efficient metabolic oxygen recovery system capable of yielding a minimum of 75% O2 recovery. A Macrofluidic Electrochemical Reactor (MFECR) technology development effort is currently underway at NASA Marshall Space Flight Center (MSFC) to significantly increase current metabolic O2 recovery efficiency, expand mission sustainability, and reduce complexity of the system. The novel design combines CO2 conversion to O2 along with C2H4 as byproduct and water electrolysis (currently conducted in two separate units) into a single compact unit that runs at standard conditions and is theoretically capable of generating O2 with a theoretical maximum metabolic CO2 conversion of 73% while consuming less than metabolic water. This paper presents a comprehensive multi-physic 3D model developed at MSFC on CO2 conversion to O2 and C2H4 at standard conditions via MFECR. The 3D spatial domain of the model is a replica of the actual MFECR’s 3D drawing generated for the MFECR fabrication and operated to recover O2 from CO2 yielding C2H4 as byproduct. Electrochemical (EC) physics that includes EC multicomponent reaction mechanisms, mass transport, and electrical current density distributions is coupled in the model with all the other physics phenomena involved in the MFECR’s process, such as two-phase flow, free and porous fluid regimes, multicomponent mass transfer, heat transfer, and DC electrical current generation along with Joule heating effect. The EC reaction sections of the MFECR consists of two porous gas diffusion electrodes (GDE) and an electrolyte serpentine channel sandwiched in the middle. The CO2 feeds the cathode serpentine channel and part of the O2 product is fed back to the anode serpentine chamber. An alkaline solution feeds the electrolyte serpentine chamber wetting the GDEs of both, the anode and cathode allowing the OH- ionic transport between them. The EC reactions in the cathode’s GDE yield C2H4 from CO2 and H2 from water while the EC reaction in the anode’s GDE yields O2. The authors will present in this paper the validation of the model using experimental data and the utilization of the validated model in building a reliable simulator that will not only assist the authors on the MFECR design but also the optimization of its operation in the ISS and future spatial human missions.
We demonstrate that it is technologically possible to obtain graphite-like films which can be used for the creation of a thermoelectric generator. The proposed technology ensures uniformity and rather small thickness of these films and allows them to be formed on a diamond-like film substrate with deposited contacts at acceptable interface. Measurements show that the electron-phonon drag effect in this system ensures thermo-emf values about 100 times as large as those provided by the diffusion process. Arrangement of the graphite-like material on a diamond-like film substrate also favors increase in the thermo-emf, which is a manifestation of the electron–ballistic phonon drag. It is established that conditions necessary for the creation of a thermoelectric generator can be achieved based on the proposed carbon nanostructures.
A kinematic model is presented for the solar cycle which includes not only the transport of magnetic flux by supergranular diffusion and a poleward bulk flow at the sun's surface, but also the effects of turbulent diffusion and an equatorward 'return flow' beneath the surface. As in the earlier models of Babcock and Leighton, the rotational shearing of a subsurface poloidal field generates toroidal flux that erupts at the surface in the form of bipolar magnetic regions. However, such eruptions do not result in any net loss of toroidal flux from the sun (as assumed by Babcock and Leighton); instead, the large-scale toroidal field is destroyed both by 'unwinding' as the local poloidal field reverses its polarity, and by diffusion as the toroidal flux is transported equatorward by the subsurface flow and merged with its opposite hemisphere counterpart. The inclusion of meridional circulation allows stable oscillations of the magnetic field, accompanied by the equatorward progression of flux eruptions, to be achieved even in the absence of a radial gradient in the angular velocity. An illustrative case in which a subsurface flow speed of order 1 m/s and subsurface diffusion rate of order 10 sq km/s yield 22-yr oscillations in qualitative agreement with observations.
The annihilation of dark matter particles is expected to yield a broad radiation spectrum via the production of Standard Model particles in astrophysical environments. In particular, electrons and positrons from dark matter annihilation produce synchrotron radiation in the presence of magnetic fields. Galaxy clusters are the most massive collapsed structures in the universe, and are known to host ∼μG-scale magnetic fields. They are therefore ideal targets to search for, or to constrain the synchrotron signal from dark matter annihilation. In this work, we use the expected sensitivities of several planned surveys from the next generation of radio telescopes to predict the constraints on dark matter annihilation models which will be achieved in the case of non-detections of diffuse radio emission from galaxy clusters. Specifically, we consider the Tier 1 survey planned for the Low Frequency Array (LOFAR) at 120 MHz, the Evolutionary Map of the Universe (EMU) survey planned for the Australian Square Kilometre Array Pathfinder (ASKAP) at 1.4 GHz, and planned surveys for Aperture Tile in Focus (APERTIF) at 1.4 GHz. We find that, for massive clusters and dark matter masses ≲100 GeV, the predicted limits on the annihilation cross section would rule out vanilla thermal relic models for even the shallow LOFAR Tier 1, ASKAP, and APERTIF surveys.
This study aims to investigate the effect of the potential window on heat generation in carbon-based electrical double layer capacitors (EDLCs) with ionic-liquid (IL)-based electrolytes using in operando calorimetry. The EDLCs consisted of two identical activated-carbon electrodes with either neat 1-butyl-1-methylpyrrolidinium bis(trifluoromethane-sulfonyl)imide ([Pyr 14 ][TFSI]) electrolyte or 1.0 m [Pyr 14 ][TFSI] in propylene carbonate (PC) as electrolyte. The instantaneous heat generation rate at each electrode was measured under galvanostatic cycling for different potential windows ranging from 1 to 4 V. First, the heat generation rates at the positive and negative electrodes differed significantly in neat IL owing to the differences in the ion sizes and diffusion coefficients. However, these differences were minimized when the IL was diluted in PC. Second, for EDLC in neat [Pyr 14 ][TFSI] at high potential window (4 V), a pronounced endothermic peak was observed at the beginning of the charging step at the positive electrode owing to TFSI - intercalation in the activated carbon. On the other hand, for EDLC in 1.0 m [Pyr 14 ][TFSI] in PC at potential window above 3 V, an endothermic peak was observed only at the negative electrode owing to the decomposition of PC. Third, for both neat and diluted [Pyr 14 ][TFSI] electrolytes, the irreversible heat generation rate increased with increasing potential window and exceeded Joule heating. This was attributed to the effect of potential-dependent charge redistribution resistance. Additionally, a further increase in the irreversible heat generation rate was observed for the largest potential windows owing to the degradation of the PC solvent. Finally, for both types of electrolyte, the reversible heat generation rate increased with increasing potential window because of the increase in the amount of ion adsorbed/desorbed at the electrode/electrolyte interface.
A large scale (2 in.-wide test vanes at 50-in. inlet diameter), low speed (100-200 ft/sec) steady flow radial cascade wind tunnel for diffuser studies was designed, built and tested. The apparatus was shown to provide flow angles (from radial) in the range of 58-72 deg with suitable spanwise profiles. The novel flow angle control mechanism was shown to work, but measured flow angles were somewhat smaller than expected. The mechanism for controlling case-wall boundary layer profile did not behave as predicted. Attempts to predict the generated flow both analytically and with computational codes are compared with initial measurements. The low speed apparatus will be used to obtain detailed data for diffuser design and analysis code verification, and to provide experience toward the design of a high speed device.
The paper demonstrates that the broad band of whistler-mode waves observed within the high density torus surrounding Io is consistent with electron cyclotron generation. Cyclotron resonant instability of Jovian energetic electrons is enhanced due to the lower resonant electron energy within the equatorial high density plasma torus surrounding the orbit of Io. The higher energy resonant electron scattering and the corresponding energetic electron lifetimes indicate that an efficient local acceleration process is required to replenish the precipitating relativistic electrons. Calculated energy deposition into the Jovian atmosphere should provide a dominant source of middle atmospheric ionization and excite a continuous band of diffuse auroral emission. It is suggested that the diffuse Jovian aurora should be influenced by the variable volcanic activity on Io which is thought to be an important source of plasma, since the cyclotron scattering process is strongly influenced by the ambient equatorial thermal plasma density.
Polymer networks with dynamic covalent bonds have been investigated for their self-healing ability, recyclability, and potential as more sustainable materials. Recent results have indicated that in some cases, bond exchange can enhance the transport of penetrants in dense networks, pointing to their potential for separations of membranes. Here, imine dynamic bonds in ethylene oxide (EO) networks with precise linker lengths were synthesized to investigate the transport of N,N′-bis(2,5-di-tert-butylphenyl)-3,4,9,10-perylenedicarboximide (BTBP), a large, anisotropic dye molecule. Networks with mesh sizes smaller than, comparable to, and greater than the size of the penetrant axes were investigated to probe the effects of bond exchange and network confinement on transport. Mesh sizes, which ranged from 0.5 to 1.62 nm, were determined from shear rheology, glass transitions by calorimetry, and probe diffusion coefficients by fluorescence recovery after photobleaching. Permanent networks with identical EO chain lengths were prepared as control samples, and up to a 3 orders of magnitude increase in diffusion coefficient is observed in the dynamic systems for short linkers containing 13 backbone atoms. The longest linkers with 71 backbone atoms show no difference between the permanent and dynamic networks. Linkers shorter than 11 backbone atoms, corresponding to a mesh size smaller than the penetrant small axis, diffusion is no longer observable on the experimental time scale, indicating a sharp cutoff attributed to the precise linkers and narrow mesh size distribution. The dynamic imine exchange time scales were compared to the diffusive hopping times of penetrants and indicate that exchange can occur during a diffusive displacement. Furthermore, these findings provide insights into the factors affecting penetrant transport in dense polymers and inspire the development of next-generation selective polymer membranes.
Conventional rate theory often uses the mean field concept to describe the effect of inhomogeneous microstructures on the evolution of radiation induced defect and solute/fission product segregation. However, the spatial and temporal evolution of defects and solutes determines the formation and spatial distribution of radiation-induced second phase such as precipitates and gas bubbles/voids, especially in materials with complicated microstructures and subject to high dose radiation. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed to investigate the effect of inhomogeneous thermodynamic and kinetics properties of defects on diffusion and accumulations of solute A in AB binary alloys. Four independent concentrations: atom A, interstitial A, interstitial B, and vacancy on [A, B] sublattice are used as field variables to describe temporal and spatial distribution and evolution of defects and solute A. The independent concentrations of interstitial A and interstitial B allow to describe their different generation rates, thermodynamic and kinetic properties, and release the assumptions of interstitial generation and sink strength used in the conventional rate theory. Microstructure and concentration dependent chemical potentials of defects are used to calculate the driving forces of defect diffusions. With the model, the effects of defect chemical potentials and mobilities on the RIS in polycrystalline AB model alloys have been simulated. The results demonstrated the model capability in predicting defect evolution in materials with inhomogeneous thermodynamic and kinetic properties of defects. The model can be extended to materials with complicated microstructures such as a wide range of grain size distribution, coating structure and multiphases as well as radiation-induced precipitation subject to severe radiation damage.
Recent advancements in the characterization of irradiated U-Pu-Zr fuels have revealed complexities that challenge existing understanding of constituent redistribution. Traditionally, models have proposed three concentric regions within the fuel, each characterized by distinctive phases and porosity. However, through detailed analysis of high burnup U-Pu-Zr, we discovered the presence of four distinct constituent redistribution regions. Particularly novel is the observation of significant Pu redistribution, a previously unreported phenomenon that necessitates a reevaluation of current models. This work aims to delve deeper into these findings, seeking to correlate mesoscale measurements of thermal diffusivity and respective thermal conductivity with the phases present in each redistribution region. To achieve this objective, we employed mesoscale thermoreflectance methods using the unique, Idaho National Laboratory (INL) developed, Thermal Conductivity Microscope (TCM) at INL’s Irradiated Materials Characterization Laboratory. The TCM employs two tightly focused lasers: one for heating to generate periodic thermal waves in the substrate, and another spatially separated probe laser to detect changes in the optical reflectivity of the gold-coated substrate resulting from thermal wave diffusion. We conducted several thermal diffusivity measurements within each region of constituent redistribution of a U-19Pu-10Zr fuel pin cross section irradiated to 11 at. % burnup. The TCM measurement positions strategically aligned with transmission electron microscopy (TEM) lift-out locations previously collected from the fuel sample. Complementary microstructural analysis techniques such as optical and scanning electron microscopy (OM/SEM), electron probe microanalysis for chemical compositions, and TEM-based selective area electron diffraction (SAED) analysis for crystallographic insights into each phase were also utilized. This comprehensive approach allowed us to correlate local thermal diffusivity data with microstructural characteristics, enabling the computation of local thermal conductivity at each position. The significance of this contribution lies in its pioneering use of the TCM for ternary fuel mesoscale examination, shedding light on the previously overlooked effects of Pu redistribution on local thermal conductivity. By informing current models capturing constituent redistribution and heat transfer, our findings pave the way for more accurate predictions of metallic fuel performance. Moreover, this work sets the stage for future comparisons with similar TCM examinations on U-19Pu-10Zr fuels at ultra-low burnup, facilitating a comprehensive understanding of thermal property changes across different burnup levels. Ultimately, our study not only enriches our understanding of the thermophysical properties of individual redistribution regions within U-Pu-Zr fuel but also offers valuable insights for the design and operational parameters of proposed next-generation fast reactors.
The lithium–sulfur (Li–S) battery is a promising next-generation energy storage technology because of its high theoretical energy and low cost. Extensive research efforts have been made on new materials and advanced characterization techniques for mechanistic studies. However, it is uncertain how discoveries made on the material level apply to realistic batteries due to limited analysis and characterization of real high-energy cells, such as pouch cells. Evaluation of pouch cells (>1 A h) (instead of coin cells) that are scalable to practical cells provides a critical understanding of current limitations which enables the proposal of strategies and solutions for further performance improvement. Herein, we design and fabricate pouch cells over 300 W h kg -1 , compare the cell parameters required for high-energy pouch cells, and investigate the reaction processes and their correlation to cell cycling behavior and failure mechanisms. Spatially resolved characterization techniques and fluid-flow simulation reveal the impacts of the liquid electrolyte diffusion within the pouch cells. We found that catastrophic failure of high-energy Li–S pouch cells results from uneven sulfur/polysulfide reactions and electrolyte depletion for the first tens of cycles, rather than sulfur dissolution as commonly reported in the literature. The uneven reaction stems from limited electrolyte diffusion through the porous channels into the central part of thick cathodes during cycling, which is amplified both across the sulfur electrodes and within the same electrode plane. A combination of strategies is suggested to increase sulfur utilization, improve nanoarchitectures for electrolyte diffusion and reduce consumption of the electrolytes and additives.
To decommission deactivated gaseous diffusion enrichment facilities, miles of contaminated pipe must be measured. The current method requires thousands of manual measurements, repeated manual data transcription, and months of manual analysis. The Pipe Crawling Activity Measurement System (PCAMS), developed by Carnegie Mellon University and in commissioning for use at the DOE Portsmouth Gaseous Diffusion Enrichment Facility, uses a robot to measure Uranium-235 from inside pipes and automatically log the data. Radiation measurements, as well as imagery, geometric modeling, and precise measurement positioning data are digitally transferred to the PCAMS server. On the server, data can be automatically processed in minutes and summarized for analyst review. Measurement reports are auto-generated with the push of a button. A database specially-configured to hold heterogeneous data such as spectra, images, and robot trajectories serves as archive. This paper outlines the features and design of the PCAMS Post-Processing Software, currently in commissioning for use at the Portsmouth Gaseous Diffusion Enrichment Facility. The analysis process, the analyst interface to the system, and the content of auto-generated reports are each described. Example pipe-interior geometric surface models, illustration of how key report features apply in operational runs, and user feedback are discussed. (authors)
Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.
In this work, we investigated the effect of proton irradiation on the performance of long wavelength infrared InAs/GaSb photodiodes (λ{sub c} = 10.2 μm), based on the complementary barrier infrared detector design. We found that irradiation with 68 MeV protons causes a significant increase of the dark current from j{sub d} = 5 × 10{sup −5} A/cm{sup 2} to j{sub d} = 6 × 10{sup −3} A/cm{sup 2}, at V{sub b} = 0.1 V, T = 80 K and fluence 19.2 × 10{sup 11 }H{sup +}/cm{sup 2}. Analysis of the dark current as a function of temperature and bias showed that the dominant contributor to the dark current in these devices changes from diffusion current to tunneling current after proton irradiation. This change in the dark current mechanism can be attributed to the onset of surface leakage current, generated by trap-assisted tunneling processes in proton displacement damage areas located near the device sidewalls.
A method of depositing by chemical vapor deposition a modified platinum aluminide diffusion coating onto a superalloy substrate comprising the steps of applying a layer of a platinum group metal to the superalloy substrate; passing an externally generated aluminum halide gas through an internal gas generator which is integral with a retort, the internal gas generator generating a modified halide gas; and co-depositing aluminum and modifier onto the superalloy substrate. In one form, the modified halide gas is hafnium chloride and the modifier is hafnium with the modified platinum aluminum bond coat comprising a single phase additive layer of platinum aluminide with at least about 0.5 percent hafnium by weight percent and about 1 to about 15 weight percent of hafnium in the boundary between a diffusion layer and the additive layer. The bond coat produced by this method is also claimed.
Diffusion flames are extensively encountered in many domestic and industrial processes. Even after many decades of research, a complete understanding of the diffusion flame structure is not available. The structure and properties of the flames are governed by the mixing (laminar or turbulent), chemical kinetics, radiation and soot processes. Another important phenomenon that affects flame structure in normal gravity is buoyancy. The presence of buoyancy has long hindered the rational understanding of many combustion processes. In gas jet diffusion flames, buoyancy affects the structure of the shear layer, the development of fluid instabilities, and formation of the coherent structures in the near nozzle region of the gas jets. The buoyancy driven instabilities generate vorticial structures outside the flame resulting in flame flicker. The vortices also strongly interact with the small-scale structures in the jet shear layer. This affects the transitional and turbulence characteristics of the flame. For a fundamental understanding of diffusion flames it is essential to isolate the effects of buoyancy. This is the primary goal of the experiments conducted in microgravity. Previous investigations, have shown dramatic differences between the jet flames in microgravity and normal gravity. It has been observed that flames in microgravity are taller and more sooty than in normal gravity. The fuels used in these experiments were primarily hydrocarbons. In the absence of buoyancy the soot resides near the flame region, which adversely affects the entrainment of reactants. It is very important to eliminate the interference of soot on flame characteristics in microgravity. The present work, therefore, focuses on the changes in the flame structure due to buoyancy without the added complexities of heterogeneous reactions. Clean burning hydrogen is used as the fuel to avoid soot formation and minimize radiative losses. Because of the low luminosity of hydrogen flames, we use rainbow schlieren deflectometry for visualization. The visualized images are digitized for quantification.The work reported here is divided into three sections; rainbow schlieren deflectometry (RSD), microgravity experiments and sub-atmospheric pressure experiments. The first section demonstrates the application of RSD for quantitative measurements in non-reacting and reacting flow systems. A computational effort to complement the experimental work is also included. In the second section, the experiments conducted at the 2.2s NASA Lewis Drop tower facility are described. The experiments were conducted to study the behavior of laminar, transitional and turbulent hydrogen flames in microgravity. The ability of RSD technique to provide quantitative data is highlighted. The final section deals with the sub-atmospheric pressure tests, which demonstrate that buoyancy in hydrogen diffusion flames can be scaled with pressure at normal gravity.
Replacement reactions commonly alter the multiscale pore structures of rocks during fluid-rock interactions. Analysis of these processes in various model fluid-rock systems during controlled laboratory experiments provides insights into the origins of microstructures found in natural materials. This study focused on understanding the effects of initial starting material permeability and resultant differences in transport pathways on porosity and mineralogical changes during limestone dolomitization.A series of replacement experiments (32–317 days in duration) have been conducted in which 1.59 cm (5/8 in.) diameter cores of two different limestones were reacted with saturated MgCl2 solutions at 200 °C. The Texas Cream (Austin Chalk) is a high-porosity, high-permeability limestone, whereas both the porosity and permeability of the Carthage Marble (Burlington Limestone) are relatively low. Altered limestones were imaged using scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM-EDX), Time-of-Flight Secondary Ion Mass Spectrometry (ToF–SIMS) and electron microprobe analysis (EMPA). A representative grain boundary of the low-porosity limestone was targeted for a focused ion beam (FIB) lift-out and characterized using transmission electron microscopy (TEM). These results were coupled with analyses of radial changes in the porosity distribution of the core derived from X-ray and neutron small- and ultra-small angle scattering ((U)SANS/(U)SAXS).The high-porosity/permeability limestone showed a four times faster bulk replacement rate than the lower-porosity/permeability material, and a different mechanism of porosity development. For the low-porosity limestone, a two-stage replacement occurred, with the reacted region of the core consisting of an inner rim in which the limestone was replaced by two calcite-dolomite solid solutions, and an outer rim in which the dolomite was replaced by magnesite. Elongated pores formed along grain boundaries at the initial limestone/dolomite reaction interface, and additional nanometer-scale porosity was formed at the secondary magnesite replacement rim. Grain boundaries were identified as preferential pathways for transport leading to dolomitization and a grain boundary diffusion rate was calculated based on microstructural characterization. In contrast, replacement in the high-porosity limestone was accompanied by porosity generation through replacement of individual grains by dolomite throughout the sample and, in longer runs, magnesite in outer parts of the core. These observations emphasize that both the mechanisms of the replacement reaction and the microstructure and chemistry of the replaced product are contingent on the initial structure of the starting material.