Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inclusion”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Low-temperature carbonization of polyacrylonitrile/graphene carbon fibers: A combined ReaxFF molecular dynamics and experimental study

Graphene inclusion in a polymer matrix is a promising route to significantly enhance the mechanical properties of low-grade carbon fibers (CFs). Using ReaxFF molecular dynamics simulation, the atomistic mechanism leading to this enhancement is investigated. We demonstrate that the graphene edges along with the nitrogen and oxygen functional groups play a catalytic role and act as seeds to expedite alignment of the all-carbon rings, which are starting sites for the growth of graphitic structures. To examine the role of this proposed mechanism that enhances the graphitic structure of PAN/graphene CFs, we discuss the experimental results wherein the PAN/graphene CFs carbonized at 1250 C demonstrate 91% (from 632 to 1207 MPa) increase in strength and 101% (from 88 to 177 GPa) enhancement in Young’s modulus compared to PAN-based CFs carbonized at 1500 C. In conclusion, these enhanced mechanical properties of low-grade carbon fibers achieved via graphene inclusion at decreased carbonization temperature provide a means to realize both energy savings and cost reduction.

42 ENGINEERING↗
DOE OSTI2022

Incorporating PCM-enabled thermal energy storage into 3D printable cementitious composites

This paper delineates the feasibility of incorporating microencapsulated phase change materials (mPCM) into 3D printable cementitious composite materials. A comprehensive experimental program was carried out to evaluate the impacts of mPCM on the printability, microstructures, mechanical and thermal properties of cementitious 3D printing ‘inks’. Results showed that the mPCM affected the printability of the cementitious ink material based on its physical properties (e.g., particle size) and volume loading – at lower volume loadings, mPCM increased the flowability of the cementitious ink material while leading to increased compressive strength and thermal conductivity for the hardened printed material. However, further increase in mPCM dosage led to a decrease in printability and, therefore, decrease in compressive strength and thermal conductivity as compared to the reference mixture. Here the results also showed that the inclusion of mPCM influence the printing parameters. In general, the inclusion of higher volume contents of mPCM necessitates a higher extrusion rate to achieve a desirable extrudability. Lastly, a thermal network model was formulated for 3D printed mPCM charged building components (e.g., wall). The study shows that microencapsulated PCM materials have good potential to be used in 3D printable cementitious mixtures for improving the thermal and energy performance of 3D printed buildings.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Influence of carbon nanotubes on microstructure and corrosion performance of additively manufactured 316L stainless steel

For this study, laser powder bed fusion (LPBF) was performed on feedstock-modified 316L stainless steel powder with 1 and 2 vol% carbon nanotubes (CNT). The corrosion resistance was evaluated following cyclic potentiodynamic polarization tests conducted in 0.6 M NaCl at room temperature, 35 and 50 °C, and immersion tests in 6 % FeCl 3 . The CNT addition has increased the pitting potential and decreased number of pits formed during immersion in FeCl 3 solution, which could be attributed to Mn-Si-O nano inclusions refinement and segregation of chromium around the inclusions. The observed corrosion behavior was correlated with the altered microstructure due to CNT addition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Validating corrosion models: A comparison of governing equations

Experimental validation of Finite Element Method (FEM) models varying electrochemical governing equations, inclusion of chemical reactions, and time on the resultant damage profile for two galvanic couples is explored. Two anode materials (Magnesium AZ31 and Carbon Steel) in contact with a cathode (Stainless Steel 304 L) were modeled in/exposed to NaCl (1 and 0.1 M respectively for the anode materials) for up to one week. The physics approach, inclusion of chemical reactions, and the boundary conditions required to accurately represent the damage profile in FEM models depended on the galvanic couple materials and, ultimately, the corrosion rate. For high rates of corrosion (i.e., magnesium anode), the Nernst-Planck equation with Electroneutrality was sufficient to describe the damage, while, for low rates of corrosion (i.e., carbon steel anode), the Laplace equation was sufficient. In all cases, the most complete governing equation (Nernst-Planck-Poisson Equation) was not necessary to accurately describe the damage. Precipitation reactions in solution also played a critical role in the predicted damage profile, especially for high corrosion rate systems. Finally, for short time periods (< 6 h), the choice of governing equations does not significantly influence damage profile results. Overall, the choice of physics to reduce error in simulations relies on the boundary conditions, geometry, conductivity of the solution, electrochemical potential differences, and time of exposure. The above results are discussed with regard to accuracy and computational savings.

Carbon steel↗
DOE OSTI2025

Monitoring pipeline integrity of underground gas storage facilities using membrane-based electrochemical sensors

Effective monitoring of internal corrosion risk is crucial to ensuring the safety and longevity of natural gas pipeline infrastructure. While electrochemical sensors are commonly used to assess corrosion rates and corrosion indicators in aqueous fluids, they are rarely used in gas pipelines as these fluids lack the ionic conductivity needed for electrochemical measurements. The inclusion of ion-conductive membranes into electrochemical sensors can extend their functionality into humidified gas streams, providing critical information about emerging corrosion events that are common during withdrawal season in pipeline systems downstream from underground storage facilities. In parallel, new protective films, like those obtained through cold spray coating, are being developed to protect oil and gas pipelines and recover losses in structural integrity due to corrosion damage. Herein, we demonstrate how membrane-based electrochemical sensors (MBES) can be used to monitor fluid corrosivity by examining their response to changes in water content for a wide range of fluid compositions. It was found that MBES readings were highly sensitive to water content changes with membrane conductivity measurements varying from 10 –6 to 10 –1 S cm -1 , and corrosion rate measurements which varied from 10 –7 to 1 mm y -1 . Electron microscopy confirmed that the self-healing characteristics of metal coating films were still active despite their inclusion into an MBES probe. In conclusion, these findings indicate that membrane-based corrosion monitoring can be expanded to monitor coated-pipeline materials and provide early detection of emerging corrosion upsets relevant to underground gas storage facilities.

Electrochemical sensor↗
DOE OSTI2023

The effect of terrain-influenced winds on fire spread in QUIC-Fire

Here in this manuscript, we describe the implementation of the terrain-following version of QUIC-URB into QUIC-Fire and a demonstration of the impacts of terrain-influenced winds on QUIC-Fire-simulated fire spread. No changes to the underlying QUIC-Fire fire spread algorithm were made other than what was required to correctly account for the inclusion of terrain. This paper summarizes simulations used to understand how the QUIC-URB terrain-influenced winds affect upslope fire behavior without additional changes to the QUIC-Fire fire spread algorithm. Previously published FIRETEC results are compared to simulation results from the modified QUIC-Fire (incorporating terrain-influenced winds) that use the same topographies and fuels. QUIC-Fire results showed overall similar behaviors in terms of how the topographies affected fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all non-flat topographies were over-predicted when compared to FIRETEC. Lateral fire growth showed similar trends with FIRETEC between topographies but did not capture the increase in spread due to a diagonal interface between grassland and forested fuel regions in the test domain. These simulations suggest possible refinements that are necessary to improve QUIC-Fire and thus guide ongoing efforts related to: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Experimental apparatus and methodology to test and quantify thermal performance of micro and macro-encapsulated phase change materials in building envelope applications

Thermal energy storage (TES) is used as a viable technology to shift peak electricity demand caused by the space cooling requirements in buildings. Passive TES is implemented in building envelope via micro and macroencapsulation methods. This study describes the use of a state-of-the-art laboratory to test different PCM inclusions in simplified building walls. A microencapsulated PCM and two macroencapsulated PCMs are tested in a controlled environment to gather data for validation purposes of PCM modelling algorithms in building energy modelling programs. Data indicates that the chamber environment and the heating and cooling system can conduct full-cycle tests of wall panels with PCM inclusions. This study also generates data from parallel tests on 4 wall panels which can be used in building energy modelling programs to validate the PCM modelling algorithms. The cyclic tests also capture the thermal effects of PCMs and complex PCM behaviors like sub-cooling in PCM hydrate-salts.

25 ENERGY STORAGE↗
DOE OSTI2024

Utilization of an Advanced Sensor network to determine fuel heating value and Real-Time net unit heat rate during transient operation

Coal-fired utility boilers are being increasingly used as variable electricity generation to resolve the imbalance in the energy market from the expansion of intermittent renewable energy. The frequent transient operation required to meet residual energy demand has created a challenge for coal-fired units to operate efficiently. This work utilizes an advanced sensor network (ASN) to calculate net unit heat rate (NUHR) of a coal-fired boiler in real time through combustion calculations and statistical correlations to provide the tools for optimizing dynamic operation. Real-time heating values that were necessary to determine fuel input energy to calculate accurate NUHR were found using both fundamental and data-driven methods. Real-time NUHR shows distinct shifts that reflect changes in process conditions that will improve the ability to optimize transient operation. Data-driven heating value correlations had 24% lower root mean square error (RMSE) than the fundamental combustion calculation approach when compared to daily retrospective proximate analysis. Furthermore, the data-driven method RMSE improved by 7% with the inclusion of ASN data. Future work is to validate by comparing unit performance with and without the inclusion of NUHR as a control parameter for the dynamic neural network.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

The combined effects of Mg 2+ and Sr 2+ incorporation during CaCO 3 precipitation and crystal growth

Calcite (CaCO 3 ) composition and properties are defined by the chemical environment in which CaCO 3 forms. However, a complete understanding of the relationship between aqueous chemistry during calcite precipitation and resulting chemical and physical CaCO 3 properties remains elusive; therefore, we present an investigation into the coupled effects of divalent cations Sr 2+ and Mg 2+ on CaCO 3 precipitation and subsequent crystal growth. Through chemical analysis of the aqueous phases and microscopy of the resulting calcite phases in compliment with density functional theory calculations, we elucidate the relationship between crystal growth and the resulting composition (elemental and isotopic) of calcite. The results of this experimental and modeling work suggest that Mg 2+ and Sr 2+ have cation-specific impacts that inhibit calcite crystal growth, including: (1) Sr 2+ incorporates more readily into calcite than Mg 2+ (D Sr > D Mg ), and increasing [Sr 2+ ] t or [Mg 2+ ] t increases D Sr ; (2) the inclusion of Mg 2+ into structure leads to a reduction in the calcite unit cell volume, whereas Sr 2+ leads to an expansion; (3) the inclusion of both Mg 2+ and Sr 2+ results in a distribution of unit cell impacts based on the relative positions of the Sr 2+ and Mg 2+ in the lattice. These experiments were conducted at saturation indices of CaCO 3 of ~4.1, favoring rapid precipitation. This rapid precipitation resulted in observed Sr isotope fractionation confirming Sr isotopic fractionation is dependent upon the precipitation rate. We further note that the precipitation and growth of calcite favors the incorporation of the lighter 86 Sr isotope over the heavier 87 Sr isotope, regardless of the initial solution conditions, and the degree of fractionation increases with D Sr . In sum, these results demonstrate the impact of solution environment to influence the incorporation behavior and crystal growth behavior of calcite. These factors are important to understand in order to effectively use geochemical signatures resulting from calcite precipitation or dissolution to gain specific information.

58 GEOSCIENCES↗
DOE OSTI2025

Comparison of the earliest NC and CC planetesimals: Evidence from ungrouped iron meteorites

Isotope anomalies in meteorites reveal a fundamental dichotomy between Non-Carbonaceous- (NC) and Carbonaceous-type (CC) planetary bodies. Until now, this dichotomy is established for the major meteorite groups, representing about 36 distinct parent bodies. Ungrouped meteorites represent an even larger number of additional parent bodies, but whether they conform to the overall NC-CC dichotomy is unknown. Here, the genetics and chronology of 26 ungrouped iron meteorites are considered through nucleosynthetic Mo and radiogenic W isotopic compositions. Secondary cosmic ray-induced modifications of these isotope compositions are corrected using Pt isotope measurements on the same samples. We find that all of the ungrouped irons have Mo isotope anomalies within the range of the major meteorite groups and confirm the NC-CC dichotomy for Mo, where NC and CC meteorites define two distinct, subparallel s-process mixing lines. All ungrouped NC irons fall on the NC-line, which is now precisely defined for 41 distinct parent bodies. The ungrouped CC irons show scatter around the CC-line indicative of small r-process Mo heterogeneities among these samples. These r-process Mo isotope variations correlate with O isotope anomalies, most likely reflecting mixing of CI chondrite-like matrix, chondrule precursors, and Ca-Al-rich inclusions. This implies that CC iron meteorite parent bodies accreted the same nebular components as the later-formed carbonaceous chondrites. The Hf-W model ages of core formation for the ungrouped irons overlap with those of the iron meteorite groups from each reservoir and reveal a narrow age peak at ∼3.3 Ma after Ca-Al-rich inclusions for the CC irons. By contrast, the NC irons display more variable ages, including younger ages indicative of impact-induced melting events, which seem absent among the CC irons. This is attributed to the more fragile and porous nature of the CC bodies, making impact-induced melting on their surfaces difficult. The chemical characteristics of all iron meteorites together reveal slightly more oxidizing conditions during core formation for CC compared to NC irons. More strikingly, strong depletions in moderately volatile elements, typical of many iron meteorite parent bodies, predominantly occur among CC irons, for reasons that remain unclear at present.

Chondrite components↗
DOE OSTI2021

A mechanistic damage model for solid oxide fuel cell ceramic materials - Part I: Constitutive modeling

A multiscale mechanistic damage model for solid oxide fuel cell (SOFC) ceramics that combines micromechanics of stiffness reduction due to material porosity change and microcracking with a continuum damage mechanics (CDM) description for the evolution of microcracks up to fracture is developed in this work. The model also accounts for volumetric swelling of the anode due to redox cycling. Porosity change and swelling during redox combined with normal operating thermomechanical loads could increase the stress, strain and damage distributions in SOFC stacks leading them to failure. The proposed model involves three governing parameters: (i) the porosity in terms of the pore or void volume fraction, (ii) material swelling magnitude prescribed at a given loading step, and (iii) the damage variable that describes microcracking caused by thermomechanical loads and swelling. Pores and microcracks in the ceramics are modeled as randomly distributed ellipsoidal inclusions with negligible stiffness by an Eshelby-Mori-Tanaka formulation combined with an inclusion orientation distribution method. Microcracking damage evolution is described by a CDM formulation. Volumetric swelling is treated in a similar way to thermal expansion in the constitutive relations. After validation and numerical checks, the damage model is used to analyze a simplified SOFC positive-electrode/electrolyte/negative-electrode (PEN) structure subjected a redox cycle in addition to thermomechanical loads experienced during normal operation of a SOFC stack.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Consequences of product inhibition in the quantification of kinetic parameters

While the potential for product inhibition in catalytic reactions is well known, the impact of neglected inhibition on measured kinetic parameters is often overlooked. The presence of product inhibition, most often caused by the competitive adsorption of products with reactants on catalytic active sites, is difficult to determine a priori for an arbitrary catalytic system. The significance of product inhibition relies on the concentration of the products, their adsorption thermodynamics on catalytically relevant sites, and process parameters such as temperature and pressure. When inhibition is significant, however, apparent activation energies and reaction orders vary from the differential-reactor apparent activation energy by a factor of 1/(1 - δ), where δ is the total inhibition order (e.g., the factor (1 - δ) = 1.6 for a system with product inhibition of -0.6 order). This is illustrated here with the kinetics of NO oxidation over Cu ion clusters (Cu x O y ) in Cu-SSZ-13, for which the product NO 2 inhibits the forward reaction. Furthermore, in the presence of inhibition, when only reactants are fed to a flow reactor or placed in a batch reactor, there is often no practical conversion that is low enough to guarantee differential behavior. Inclusion of products in the feed solves this problem, allowing accurate determination of kinetic parameters such as apparent reaction orders and activation energies. We also demonstrate that evaluation of the necessity of co-feeding products to assure measurement of differential-reactor data in a given catalytic system is straightforward from a plot of the log of the rate (or conversion) versus the log of the space time in a flow reactor or elapsed time in a batch reactor. We encourage inclusion of this test in all kinetic analyses that are reasonably approximated by power law rate expressions.

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

Growth interface study of CdTeSe crystals grown by the THM technique

The Traveling Heater Method (THM) is a widely accepted technique to grow high-quality detector grade CdTe and CdZnTe ingots, especially for X-ray and gamma-ray detector applications. Unlike melt-growth techniques, the growth interface for ingots produced by the THM technique consists of two different compounds, viz. the solidified region (just below the interface) consists of the compound to be grown (e.g., CdZnTe for CdZnTe growth), while the molten zone just above the growth interface consists of Te-rich CdZnTe. Thus, optimization of the growth parameters is critical to obtain a clean growth interface with the required shape and presence of minimal Te-rich inclusions. Here, in this study, CdTeSe ingots were grown employing the THM technique to investigate the growth interface. Both macroscopic and microscopic behavior of the interface were studied. The study revealed that a microscopically smooth growth interface can be achieved by optimizing the growth parameters, which is essential for obtaining high-quality, inclusion-free ingots.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

Ruthenium solubility and its impact on the crystallization behavior and electrical conductivity of MoO 3 -containing borosilicate-based model high-level nuclear waste glasses

The present study focuses on investigating the solubility of RuO 2 in a borosilicate-based model high-level waste glass and understanding its impact on the crystallization behavior and electrical conductivity of the resulting vitrified waste forms. The solubility limit of RuO 2 in the investigated glass composition has been determined to be 460 ppm by weight. The higher concentration of RuO 2 results in sub-micron sized Ru-rich inclusions in the glassy matrix, which eventually agglomerate to form needle-like and polyhedral RuO 2 crystals. It is observed that RuO 2 selectively promotes the crystallization of the rare-earth apatite phase over the powellite phase. The as-synthesized RuO 2 -containing glasses exhibit semiconducting behavior with a similar level of electrical conductivity below the glass transition. Here, the percolation of non-uniformly distributed RuO 2 inclusions may result in a formation of short-range low-resistive conduction pathways in the host glass matrix leading to an apparent metallic-like behavior of selected thin samples with the highest ruthenium content.

36 MATERIALS SCIENCE↗
DOE OSTI2020

On residual stress analysis and microstructural evolution for stainless steel type 304 spent nuclear fuel canisters weld joint: Numerical and experimental studies

Chloride-induced stress corrosion cracking (CISCC) in the weldments of austenitic stainless steel canisters is one of the primary safety concerns during the dry storage of used nuclear fuel at Independent Spent Fuel Storage Installations in coastal areas. In order to evaluate the CISCC behavior in the canister, a 3D sequentially coupled thermo-mechanical finite element model was built to simulate the residual stresses induced by two intersecting longitudinal and circumferential multi-pass welds in austenitic stainless steel mockup canister. Weld-induced residual stresses from simulation were in good agreement with the experimental measurement results by deep-hole drilling and contour methods. Through-wall high tensile axial and hoop stresses were observed in both the fusion zone and heat affected zone (HAZ) in the longitudinal and circumferential weld, respectively. On the other hand, the microstructure within the weldments of 304 stainless steel mockup was characterized using optical microscopy, scanning electron microscopy, electron back-scattered diffraction and energy dispersive X-ray spectroscopy methods. A few amount of MnS inclusions in both the base metal and HAZ, and a lot of (Mn,Si) x O y inclusions in the fusion zone were found. Meanwhile, there is some α'-martensite in the base metal, while the martensite was reversed in the HAZ. In the fusion zone, there is a lot of δ-ferrite. The martensite reversion phenomenon was also verified by the combination of Feritscope test and thermal simulation, which could predict a martensite reversion temperature as 400C°. The welding residual stresses and microstructure analysis results could offer the instructive information for the prediction of CISCC behavior of the spent nuclear fuel canisters.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Global and local texture development during initial plastic deformation of cold-pilgered Zircaloy-4 tubing

In this work we present a thorough characterization of the texture development during the initial step of a cold pilgering process of a Zircaloy-4 TREX tube, on which the axial cross section is reduced by 80% over a 40cm length with an overall Q-ratio of 2.34, followed by a recrystallization thermal treatment. Texture across a tapered specimen representing this working zone was characterized by neutron diffraction and high energy synchrotron x-ray diffraction at 28 positions (4 locations along its length and 7 locations across its thickness with 0.3 mm spatial resolution), and on the deformed end after the thermal treatment. Results show that the starting material (TREX tube) presents a texture gradient through its wall thickness, with the c poles aligned to the hoop direction at the outer surface (with radial and hoop Kearns factors: f R =0.36, f H = 0.50), and gradually rotating around the axial direction to become a quasi-radial orientation (~35 o out of the radial direction) at the inner surface ( f R =0.43, f H =0.35). The overall texture becomes more radial due to the cold pilgering, but the texture gradient remains during the whole deformation process, being always more radial at the inner surface, ( f R =0.53, f H =0.37 at 80%) than at the outer surface ( f R =0.41, f H =0.50at 80%). The most remarkable changes in texture produced by the rolling process is a 30 o rotation of the grains around their c-axis, producing a marked increase of the <11 2 ¯ 0> // axial fiber texture at expenses of the original <10 1 ¯ 0> // axial fiber texture observed in the TREX tube. The observed changes in texture during deformation were correctly explained by visco-plastic self-consistent (VPSC) modelling. This was achieved after the inclusion of the basal <a> slip system, which plays a central role to explain the 30 o rotation around the c-axis during deformation. The texture evolution was correctly described after the inclusion of a drag effect to grain rotation due to the interaction with neighbouring grains. The through thickness texture differences are greatly reduced by the recrystallization thermal treatment, with similar Kearns factors at the inner surface ( f R =0.50, f H =0.28) and the outer surface ( f R =0.48, f H =0.39). The predominance of the <10 1 ¯ 0> // axial fiber texture is recovered after recrystallization.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Factors influencing mode choice of adults with travel-limiting disability

Introduction: Despite the plethora of research devoted to analyzing the impact of disability on travel behavior, not enough studies have investigated the varying impact of social and environmental factors on the mode choice of people with disabilities that restrict their ability to use transportation modes efficiently. This research gap can be addressed by investigating the factors influencing the mode choice behavior of people with travel-limiting disabilities, which can inform the development of accessible and sustainable transportation systems. Additionally, such studies can provide insights into the social and economic barriers faced by this population group, which can help policymakers to promote social inclusion and equity. Method: This study utilized a Random Parameters Logit model to identify the individual, trip, and environmental factors that influence mode selection among people with travel-limiting disabilities. Here, using the 2017 National Household Travel Survey data for New York State, which included information on respondents with travel-limiting disabilities, the analysis focused on a sample of 8,016 people. In addition, climate data from the National Oceanic and Atmospheric Administration were integrated as additional explanatory variables in the modeling process. Results: The results revealed that people with disabilities may be inclined to travel longer distances walking in the absence of suitable accommodation facilities for other transportation modes. Furthermore, people were less inclined to walk during summer and winter, indicating a need to consider weather conditions as a significant determinant of mode choice. Moreover, low-income people with disabilities were more likely to rely on public transport or walking. Conclusion: Based on this study’s findings, transportation agencies could design infrastructure and plan for future expansions that is more inclusive and accessible, thus catering to the mobility needs of people with travel-limiting disabilities.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Reaction path model of the formation of abiotic immiscible hydrocarbon fluids in subducted carbonated serpentinites, Lanzo Massif (Western Italian Alps)

Fluids generated from subducted slabs participate in the cycling of deep carbon in the crust and upper mantle. In these fluids, aqueous carbon species vary in oxidation state between +IV and -IV depending on whether the fluids are oxidizing or reducing, respectively. Most studies of subduction-zone fluids have focused on oxidized carbon species. However, recent studies of natural samples have demonstrated the occurrence of deep, reducing fluids, generated in both subducted oceanic upper mantle and crustal rocks. CH 4 -H 2 -rich fluid inclusions in subducted carbonated serpentinites have demonstrated the existence of abiotic, immiscible, hydrocarbon fluids at upper mantle conditions. To investigate the formation of such immiscible hydrocarbon fluids during the evolution of subducted carbonated serpentinites, we used equilibrium constants from the Deep Earth Water model to carry out predictive chemical mass transfer modeling to simulate the alteration reactions. A novel feature of the models was the inclusion of an immiscible hydrocarbon fluid containing six components (CH 4,f , C 2 H 6,f , C 3 H 8,f , isoC 4 H 10,f , CO 2,f , H 2,f ). This feature enabled prediction of the formation of a separate immiscible fluid in equilibrium with aqueous species and minerals. We developed a predictive reaction path model of invasive H 2,f reacting with carbonated serpentinites and interstitial aqueous fluids for comparison with the natural samples from the Lanzo Massif, western Italian Alps. Over a range of temperatures and pressures, immiscible hydrocarbon fluids formed in association with altered mineral assemblages. Reaction progress caused the transformation of carbonated serpentinites and the formation of clinopyroxene, brucite, graphite, and hydrocarbon fluids, along with changes of pH, logfO 2 , and aqueous species. CH 4,f was the most abundant hydrocarbon species in all the models. The overall results at 2.0 GPa and 400 to 450 °C were consistent with the natural samples from the Lanzo Massif. Interestingly, large amounts of H 2 O formed due to oxidation of H 2 . More hydrocarbons and H 2 O formed in models with lower fluid/rock mass ratios or with more reactant H 2 . Models at different pressure and temperature conditions showed similar results with some variation in the relative stabilities of aragonite, graphite and olivine solid solution, and associated differences in mineral sequences, hydrocarbon fluids, values of aqueous species, and the final log fO 2 and pH. As a result, our models strongly support the laboratory and field evidence that reduction of carbonated serpentinites by infiltrating H 2 fluids can cause the formation of immiscible, abiotic hydrocarbon fluids in subduction zones.

58 GEOSCIENCES↗