Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “isotherm”

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

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

At least 1,009 records · Page 56

Influence of polycarboxylate superplasticizer, citric acid and their combination on the hydration and workability of calcium sulfoaluminate cement

Highlights: • The dispersion efficiency of polycarboxylate superplasticizers drops quickly over time. • Citric acid maintains the dispersion efficiency of PCE over time. • Competitive adsorption between citric acid and PCE occurs. • The combination of citric acid with PCE decreases the initial dispersion efficiency of PCE. PCEs are well known to improve the initial fluidity of CSA. However, their dispersion efficiency drops quickly over time. This issue can be solved by incorporating retarders. In this context, this paper deals with the influence of citric acid, used as a retarder, PCE and their combination on the hydration and workability of CSA. Isothermal calorimetry, XRD and TG analysis were used to describe the hydration process, while workability was characterized with the mini-cone test. Adsorption behavior was investigated using total organic carbon analyzer coupled with ion chromatography. Results show that the introduction of citric acid retained the dispersion efficiency of PCE over time. However, the initial dispersion efficiency of PCE was decreased by citric acid as the latter tend to adsorb first on the surface of cement grains, inhibiting the adsorption of PCE. A dispersion model was proposed to describe the acting mechanism of these admixtures on CSA.

36 MATERIALS SCIENCE↗

Direct observation of C 3 S particle dissolution using fast nano X-ray computed tomography

Tricalcium silicate (C 3 S) occupies 50 % to 70 % of ordinary portland cement (OPC) by mass and it is an important component affecting the hydration of OPC [1], [2], [3], [4], [5], [6], [7]. Generally, the hydration of C 3 S is described by two processes: the dissolution of C 3 S particles and the precipitation of hydration products. While it is understood that the dissolution rates of C 3 S vary with time, more precise measurements are needed to understand this process. Many mechanisms have been proposed to explain the time-evolving dissolution rates of C 3 S [2]. The metastable barrier hypothesis suggests that a thin metastable layer of hydrates forms around the C 3 S particle surface and prohibits the access of grains to the aqueous solution [8], [9], [10], [11], [12]. The slow dissolution step hypothesis suggests that the increased ion concentration from the initial reaction delays the C 3 S dissolution [2], [13], [14], [15], [16], [17]. More recent publications suggest that C 3 S may react differently depending on the existence of crystallographic defects [18], [19], [20]. Etch pits are thought to open on the particle surface during the initial reaction and this contributes to the C 3 S dissolution [21], [22]. As hydrates precipitate and cover these highly reactive surfaces, hydration slows down and the induction period starts [18], [23], [24], [25], [26]. Many experiments have been conducted to test the aforementioned mechanisms. Some hydration studies utilize bulk measurements, such as isothermal calorimetry [27], [28], [29], pore solution analysis [30], quasi-elastic neutron scattering [31], and nuclear magnetic resonance spectroscopy [32], [33]. One limitation of these measurements is that they do not provide direct and detailed information on the individual C 3 S particles. Some other studies utilize imaging techniques, such as scanning electron microscopy (SEM) [34], [35], [36] and transmission electron microscopy (TEM) [37]. However, SEM/TEM cannot track the evolution of individual particles throughout hydration [34], [35], [38], [39], [40] and they do not give insights into the microstructure of materials before hydration [34], [35], [39]. This makes it challenging to draw strong conclusions from only SEM or TEM observations. Synchrotron X-ray tomography techniques have been used more broadly in recent years to study cement hydration. They are not only non-destructive but also able to image a sample in full 3D with resolutions that can reach from micron to nanoscale. Nano computed tomography (nCT) is one technique that has been applied to study cement hydration at the nanoscale [26], [41]. A typical nCT can reach a pixel size from 15 to 65 nm, providing enough detail for observing features <1 μm. However, nCT often takes >0.5 h to finish one scan. This makes the application of this technique on continuous scans for in-situ observations challenging. Fast X-ray computerized tomography (fCT) is another technique that has shown success in studying the time-evolving cement microstructures [20], [42], [43], [44], [45], [46], [47], [48], [49]. Due to the high flux of the X-ray beam from the synchrotron ring, fCT allows a scan to be captured within 1 min at a pixel size of 1 μm. This allows a paste sample to be continuously scanned during the hydration process. However, the micron-sized resolutions limit does not provide detailed insights for particles <5 μm [20], [49]. Fortunately, the combination of nCT and fCT has allowed the development of fast nano X-ray computed tomography (fnCT). fnCT can capture a 3D data set in <2 min at a pixel size of 50 nm. This makes this procedure an exciting method to evaluate hydrating pastes. fnCT collects multiple X-ray radiographs at various rotation angles and generates a 3D model of the scanned sample, which is also referred to as a 3D tomography [50], [51]. In one tomography, the X-ray absorptions of different components (e.g., C 3 S and hydrates) differ as functions of density and chemistry [52], [53]. These X-ray absorption contrasts can be used to extract detailed information about the 3D microstructure [26], [54], [55]. In this paper, fnCT is used to collect time-lapse tomographs of hydrating C 3 S paste from 18 min after mixing to 7 h of hydration. The bulk measurements of anhydrous C 3 S, as well as the microstructural changes of individual C 3 S particles, are directly observed, quantified, and discussed. The dissolution behavior of C 3 S particles at various size scales is systematically analyzed and compared. This work aims to find the relationship between the size of C 3 S particle sizes and their dissolution rates. This provides significant insights into the early-age hydration of C 3 S on length and time scales not previously possible. Because of the magnitude of the data and the substantial amount of observations, this work will solely focus on the change in the anhydrous particles. Changes in the hydration products will be reported in future work.

42 ENGINEERING↗

A perspective on challenges and opportunities in developing high entropy-ultra high temperature ceramics

The emergence of a new class of high entropy ultra-high temperature ceramics (HE-UHTCs) is expected to possess exceptionally superior mechanical, oxidation, and erosion resistance as compared to conventional UHTCs. The improvement in the properties of HE-UHTCs is attributed to the high configurational entropy of the multi-component system stabilized into a single-phase, but with little discussion on entropy contributions. This perspective article presents a viewpoint on the critical scientific question regarding the actual role entropy plays in stabilizing and enhancing the thermo-mechanical properties of a multi-component system. The high throughput first-principle calculations, which can provide the entropy formation ability (EFA) of thermodynamically stable HE-UHTCs, can prove to be valuable toolset to this puzzle. HE-UHTCs show enhancement in the mechanical properties much higher than the predicted rule of the mixture, which cannot traditionally be explained by a mere solid solution formation effect. Additionally, isothermal oxidation studies of HE-UHTCs limited to furnace testing state remarkable oxidation resistance which is attributed to the passivity of complex, layered structure, and composition of the oxide scale. However, whether all the constituents in the HE-UHTC system react simultaneously or preferential oxidation takes place needs to be addressed. To address all these issues, an integrated computational and experimental approach has been postulated towards designing a thermal protection system (TPS) for re-entry applications. Succeeding the approach will provide new opportunities in the compositional space of HE-UHTCs with tailored properties.

36 MATERIALS SCIENCE↗

Rich methane oxidation on Pt/Pd/Al 2 O 3 : Steady state performance, multiplicity features, and spatial patterns

Performance features are reported for the catalytic sub-stoichiometric (rich) oxidation of methane on a Pt + Pd/ Al 2 O 3 washcoated monolith over a range of feed temperatures and O 2 /CH 4 for both dry and wet feeds. Near isothermal steady state multiplicity and co-existing states with different methane conversions are presented. The methane conversion vs. O 2 feed concentration consists of a higher conversion regime at lower concentrations and a lower conversion, inhibited regime at higher concentration. Here, the conversion maximum is accompanied by an overlap of the two regimes in the form of a clockwise hysteresis loop. Spatial concentration profiles for the dry feed reveal an upstream oxidation zone and downstream stream reforming zone. Measurements extending beyond the catalyst imply the co-existence of active and inhibited states. Addition of water promotes oxidation and steam reforming while the multiplicity and nonuniform states persist. Underlying mechanistic aspects responsible for the inhibition, rate multiplicity and co-existing states are discussed.

03 NATURAL GAS↗

Balanced Fast-SpaciMS capillary configurations provide practically noninvasive channel-average measurements in catalytic monoliths

Spatially resolved capillary-inlet mass spectrometry (SpaciMS) provides a detailed picture of the spatiotemporal evolution of reaction network in catalytic monoliths. In the present work, we combine the SpaciMS experiments with a newly developed non-isothermal 3D CFD model for heterogeneously catalyzed reactive flows, including diffusion and permeation through the coated catalyst and channel wall. We explore how the capillary size and sampling rate can be balanced to minimize the impact on probed-channel conversions and provide species concentration measurements representative of the free channel average. In all studied configurations, the balanced sampling-rate fraction is noticeably higher than the corresponding capillary occlusion fraction. For a typical 350-micron capillary that occupies 11% of the channel cross- section, the balanced sampling rate represents 43% of the channel flow. Furthermore, the present work shows that together, a balanced SpaciMS configuration and numerical simulation provide accurate channel-averaged information of a monolithic catalyst under realistic operation conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Intrinsic kinetics of water-inhibited ultra-lean methane oxidation over PtPd-Mg/ θ -Al 2 O 3 catalyst

Here, this study develops intrinsic methane oxidation kinetics for ultra-lean methane conditions in the presence of water over a highly active and stable PtPd–Mg/θ-Al 2 O 3 catalyst. Comprehensive laboratory experiments were conducted over a wide range of methane concentrations (150–1200 ppm CH 4 ), water contents (1–5% H 2 O), and industrially relevant space velocities (80,000 ≤ GHSV ≤ 160,000 h -1 ). These systematic experiments informed a two-dimensional, axisymmetric, multiscale reactor model that was used to develop and validate methane oxidation kinetics under practically relevant conditions, including non-isothermal operation and high conversion regimes. Combined experimental and modeling results revealed significant intraparticle diffusion resistance and transport-induced reaction exotherm at elevated temperatures, which limited catalyst utilization despite high intrinsic activity. These transport effects were explicitly incorporated into the reactor model, enabling accurate estimation of intrinsic kinetic parameters without reliance on conventional effectiveness-factor corrections. The resulting kinetic model successfully captured both kinetically controlled and mass-transfer-limited regimes and reliably predicted CH 4 conversion across broad ranges of temperature, methane concentration, and water content

Heat and mass transfer limitations↗

Mobility of heavy rare earth elements in magmatic-hydrothermal systems: Experimental determination of DyPO 4 solubility in supercritical fluids

Dysprosium (Dy) is part of the heavy rare earth elements (REE) and a critical component in the production of high-performance super magnets used in wind turbines, electric vehicles, and other green technologies. The HREE phosphate, xenotime, is commonly associated with hydrothermal systems, where the mobility, fractionation, and deposition of REE are controlled by the stability of aqueous complexes and solubility of REE minerals which depend on temperature, pressure, and fluid composition (salinity, pH, and ligand availability). Here, we conducted a series of synthetic DyPO 4 solubility experiments in NaCl-HCl-bearing aqueous solutions (pH of 2 at 25 °C) from 500 to 700 °C and 1.5 kbar using externally heated pressure vessels. DyPO 4 displays a strong prograde solubility with temperature across all studied fluid salinities (0–0.5 m NaCl). The DyPO₄ solubility isotherms at 500 °C, 600 °C, and 700 °C show a complex relationship with salinity where DyPO₄ solubility is highest (39.2–1552 ppm at 500–700 °C) in 0 m NaCl solutions, decreases strongly (21–331 ppm at 500–700 °C) in the 0.06 m NaCl solutions and remains relatively constant in the 0.25 m (42.8–248 ppm) and 0.5 m NaCl (40.2–124 ppm) solutions. Thermodynamic modeling suggests that pH increases as a function of temperature and salinity, and controls solubility due to the stability of chloride complexes (i.e., DyCl +2 and DyCl 2 + ) at 500–600 °C, transitioning to hydroxyl complexes (i.e., Dy(OH) 3 0 ) dominating at 700 °C and relatively alkaline pH. Comparison with NdPO 4 solubility data indicates that the light REE in monazite is significantly more soluble than the heavy REE in xenotime in saline fluids at high temperature (500–700 °C), implying efficient light/heavy REE fractionation in supercritical fluids. The lower solubility of DyPO 4 relative to NdPO 4 suggests that at high temperatures, acidic (pH 2–4), low-salinity (∼5 wt% NaCl) fluids will favor the precipitation of xenotime over monazite as it will cool from 700 °C to below 500 °C. These observations align with paragenetic mineral assemblages in NYF pegmatites (e.g., Baveno, Italy), where xenotime commonly immobilizes heavy REE during the high-temperature hydrothermal stages.

58 GEOSCIENCES↗

Reduction and bacterial adsorption of dissolved mercuric ion by indigenous bacteria at the Oak Ridge Reservation site

Mercury exists in various forms in the environment and the indigenous bacteria mediated processes have the potential to be used for mercury remediation. In this study, two mixed cultures of indigenous bacteria at the Oak Ridge Reservation site (i.e., ORR soil culture and ORR sediment culture) were selected to study the microbial mediated mercuric reduction under an aerobic condition as well as mercury adsorption onto bacterial surfaces. PCR analysis was performed to provide insights into the microbial community. The mercuric volatilizing experiment demonstrated the mercuric reducing capacity for both ORR cultures, in which the Pseudomonas genus was the dominating Hg 0 producer. The investigation of the impact of the sole carbon source revealed the energy-dependent characteristics of the mercuric reduction in this study. Namely, the mercuric reduction was nearly not impacted by the type of carbon source but positively related to the energy that a unit amount of substrate could provide. The study also indicated that the mercury adsorption competed with the reduction. According to the fitting of the Langmuir isotherm, the ORR soil culture was found to have a higher mercury adsorption capacity (i.e., 67.5 mg Hg/g dry biomass) than the ORR sediment culture (i.e., 53.1 mg Hg/g dry biomass). The negative correlation between the reduced mercury mass and adsorbed mercury mass was identified for both ORR cultures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A fully light-driven approach to separate carbon dioxide from emission streams

Carbon capture from industrial point sources is an essential component of the global effort to mitigate climate risks. However, traditional approaches require significant energy input—often provided, counterproductively, by fossil fuel combustion. Using sunlight directly as the energy source would significantly improve the energy efficiency of carbon capture processes. Herein, we report the first fully visible-light-driven CO 2 separation system, in which carbon capture is achieved via the photoenolization/addition reaction of inexpensive 2-methylbenzophenone with CO 2 , and CO 2 release is realized through an intramolecular photodecarboxylation reaction. Finally, this system operates isothermally, works with natural sunlight, and facilitates CO 2 removal from natural gas flue emissions, providing a blueprint for other non-thermal chemical separations.

Ahmad, Bayu I. Z. [Cornell Univ., Ithaca, NY (Unit↗

Inelastic equation of state for solids

In this study, a complete inelastic equation of state (IEOS) for solids is developed based on a superposition of thermodynamic energy potentials. The IEOS allows for a tensorial stress state by including an isochoric hyperelastic Helmholtz potential in addition to the zero-kelvin isotherm and lattice vibration energy contributions. Inelasticity is introduced through the nonlinear equations of finite strain plasticity which utilize the temperature dependent Johnson–Cook yield model. Material failure is incorporated into the model by a coupling of the damage history variable to the energy potentials. The numerical evaluation of the IEOS requires a nonlinear solution of stress, temperature and history variables associated with elastic trial states for stress and temperature. The model is implemented into the ALEGRA shock and multi-physics code and the applications presented include single element deformation paths, the Taylor anvil problem and an energetically driven thermo-mechanical problem.

42 ENGINEERING↗

An attention-based neural ordinary differential equation framework for modeling inelastic processes

To preserve strictly conservative behavior as well as model the variety of dissipative behavior displayed by solid materials, we propose a significant enhancement to the internal state variable-neural ordinary differential equation (ISV-NODE) framework. In this data-driven, physics-constrained modeling framework internal states are inferred rather than prescribed. The ISV-NODE consists of: (a) a stress model dependent on observable deformation and inferred internal state, and (b) a model of the evolution of the internal states. The enhancements to ISV-NODE proposed in this work are multifold: (a) a partially input convex neural network stress potential provides polyconvexity in terms of observed strain while leaving the inferred state unconstrained, and (b) an internal state flow model uses common latent features to inform novel attention-based gating and drives the flow of internal state only in dissipative regimes. We demonstrated that this architecture can accurately model dissipative and conservative behavior across an isotropic, isothermal elastic-viscoelastic-elastoplastic spectrum with three exemplars, while maintaining fundamental principles by design.

97 MATHEMATICS AND COMPUTING↗

Non-equilibrium low-temperature plasma-assisted combustion of iso-octane: Perturbing pyrolysis and oxidation kinetics

Here, in this study, a plasma-coupled flow reactor facility is used to examine the effects of non-equilibrium low-temperature plasmas on perturbing the pyrolysis and oxidation kinetics of iso-octane. Experiments were performed in highly dilute reactive mixtures of nitrogen, at near isothermal conditions for temperatures ranging from 523 K to 1203 K. Experiments cumulatively demonstrated enhanced chemical reactivity with the plasma for temperatures below 900 K, and a lowering of the hot-ignition temperature. Detailed kinetic insight was derived from a 0D plasma-coupled kinetic model, utilizing a constructed mechanism that combined both plasma-specific chemistry and the neutral combustion chemistry. For pyrolysis conditions, the model displayed relatively good agreement with fuel consumption and the formation of most intermediates compared to the experimental data, demonstrating the model is able to accurately predict primary radical formation from the plasma directly interacting with the fuel. Enhanced reactivity was attributed to collisional quenching of excited-states of N 2 with fuel, which led to efficient fuel fragmentation and enhancement of the H-radical flux. For oxidation conditions, the model displayed satisfactory agreement with the experiments. Model predictions were able to accurately predict fuel consumption and most intermediate speciation data for T > 800 K, but most discrepancies were towards T < 800 K in particular with oxygenated intermediates. In the presence of oxygen, plasma effects were predominantly spent on efficient enhancement of O- and H-radical fluxes, leading to further fuel fragmentation and initiation of both the OH- and HO 2 -radical pools. Subsequent reactivity of iso-octane was then dictated by the response of the temperature-dependent neutral chemistry. At low-temperatures (T = 643 K), enhanced fuel radicals and O 2 -additon chemistry lead to the formation of oxygenated species, while at intermediate temperatures (T = 843 K) net decrease in OH-radical reactivity led to an increase in hydrocarbon speciation. Near the self-ignition threshold (T = 1163 K), radicals generated by high-temperature branching reactions dominate the oxidation process and effectively ignition. This study ultimately demonstrated that the enhancement of radicals afforded by the plasma causes a deviation in known understanding of iso-octane kinetics in some regards and warrants future studies to reconcile these discrepancies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Dependence of predicted bulk properties of hexagonal hydroxyapatite on exchange-correlation functional

Hydroxyapatite (HA) is the main component of human bones and teeth. HA has also been widely applied in various technological fields due to its unique properties. Reliable computational simulations for this critical material often require input parameters obtained from the first-principles density functional theory (DFT) calculations with appropriate exchange–correlation functionals. Previous DFT calculations are insufficient to verify the reliabilities of various functionals, particularly as they fail to assess predictions for multiple properties of the material. In this paper, we first select 18 different functionals to calculate geometric, elastic, electronic, and thermodynamic properties of hexagonal HA bulk crystal. We find that the results from optB86b-vdW and optB88-vdW functionals with dispersion corrections have the overall best agreement with available experimental data. Then, we choose optB88-vdW functional, as well as PBE functional without dispersion corrections as a comparison, to perform extensive first-principles DFT phonon calculations under the quasiharmonic approximation. Various thermodynamic properties (including phonon contributions to internal energy, entropy, and Helmholtz free energy) and thermal parameters (including thermal expansions of volume, thermal expansion coefficients, heat capacities, isothermal bulk moduli, etc.) versus temperature are consequently obtained. By comparing these quantities, we report that the results from optB88-vdW functional have significantly better agreement with available experimental data than those from PBE functional although the latter has been widely used in previous DFT calculations for HA-based material systems.

36 MATERIALS SCIENCE↗

A simple and practical process modeling methodology for pressure swing adsorption

Although many dynamic models exist for the design and simulation of pressure swing adsorption (PSA) processes, these models involve the solution of a complex system of coupled partial differential equations. Process engineers need a simple, practical, and yet robust short-cut model that helps decide whether to implement a PSA system in a process flowsheet. This work presents a “virtual” moving bed modeling methodology that considers only mass and energy balances and adsorption isotherms to describe the cyclic steady state behavior of PSA systems. Similar to tray efficiencies in distillation calculations, adsorption efficiencies are further introduced to account for system “non-ideality.” Finally, a lab-scale air separation system is used to illustrate the application of this modeling methodology.

42 ENGINEERING↗

Sulfuric Acid Improves the Reactivity of Zeolites via Dealumination

The demand for supplementary cementitious materials (SCMs), like fly ash, slag, and silica fume, is projected to rise over the next few decades. However, global decarbonization efforts, such as decommissioning of coal-fired power plants, will limit the supply of SCMs. Therefore, alternative SCM sources are a critical need for the cement and concrete industry. However, the high crystallinity of zeolites limits their pozzolanic reactivity. The goal of this work was to improve the reactivity of two common zeolites, namely clinoptilolite and SSZ-13, by treating them with sulfuric acid. The treatment enhanced zeolite reactivity via dealumination and decreased particle sizes, which led to increased reactivity with calcium hydroxide, as indicated by a 42.0% increase in reaction heat measured by isothermal calorimetry. Here, we elucidate mechanistic understanding that smaller particle sizes induced by acid treatment accelerated peak heat evolution, while free aluminates from dealumination increased heat via direct reaction with aqueous Ca 2+ .

36 MATERIALS SCIENCE↗

Assessing the hygrothermal performance of bio-based materials in building wall systems

Building envelope systems are crucial in regulating thermal and moisture exchange between interior and exterior environments, accounting for approximately 28 % of building energy consumption in the United States with walls being the primary contributors. Improper selection of building envelope materials can lead to moisture-related issues, reduced resilience, and compromised durability. Hygrothermal performance assessment is a key factor in efficient building design. As such, improving the energy and hygrothermal performance of opaque wall materials, through careful assessment of material choices, is essential to enhancing building resilience, lowering energy costs, and improving occupant comfort. As the building industry seeks new strategies to reduce material energy intensity, bio-based materials emerge as a promising solution. However, their long-term hygrothermal performance in building envelope systems remains underexplored. To fill this gap, this study evaluates the hygrothermal behavior of 13 bio-based materials in residential wall systems across three U.S. climate zones. Laboratory experiments were performed to measure material properties such as density, thermal conductivity, moisture transmission, and sorption isotherms. These data were integrated into the WUFI® simulation tool to assess wall hygrothermal performance in Houston, Baltimore, and Chicago. A three-phase modeling approach was used: (1) baseline residential walls with oriented strand board (OSB) and gypsum board; (2) replacing OSB with bio-based materials; and (3) replacing drywall with bio-based materials. Results showed that the evaluated bio-based materials maintained acceptable moisture thresholds of ≤ 16 % across all climates, confirming their viability as an alternative for current sheathing materials. Furthermore, this study provides a foundation for future research and innovation in material science on the use of certain bio-based materials in high-performance, low energy use residential construction. Ultimately, providing critical data, offering a database of bio-based material properties, and supplying a simulation-based approach will help designers make informed decisions for future efficient building practices.

Bio-based materials↗

Modelling uranium dioxide corrosion under repository conditions: A pore-scale study of the chemical and thermal processes

This paper investigates the corrosion behaviors of spent uranium dioxide (UO 2 ) fuel when placed in geological repositories. We performed pore-scale reactive transport simulations of the UO 2 corrosion process in a defective fuel rod with different orientations of breach or fracture on its clad. It is found that the corrosion rate has strong dependency on the breach orientation. The highest corrosion rate of UO 2 is calculated when the angle between the groundwater flow direction and the clad fracture reaches 180°. UO 2 corrosion simulations with different flow rates are compared. The results show the higher flow rate accelerates the corrosion of UO 2 fuel. The effect of pH on the corrosion process is also determined. It is indicated that higher pH can limit the corrosion of UO 2 by reducing the reaction rate. The dependence of reaction rate and time on reactive surface area is explored. Spent fuel with lower surface area demonstrates longer lifetime under corrosion conditions. A numerical model coupling thermal conduction and chemical reactions is developed to assess the impact of temperature on the process of UO 2 corrosion. The results show that higher temperature leads to larger corrosion rates for UO 2 . Here, the predicted reaction rates are higher in comparison with the isothermal results. Hence, our results can help improve the fundamental understanding of UO 2 corrosion in geological repositories for long-term storage of spent nuclear fuels and provide guidance for the safe operations and selection of appropriate repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling of High-Temperature Corrosion of Zirconium Alloys Using the eXtended Finite Element Method (X-FEM)

Oxidation modeling in modern nuclear fuel performance codes is currently limited by the lack of coupling with mechanics, thus preventing proper description of how high-temperature oxidation impacts mechanical properties. This is mostly due to the fact that the finite difference formalism adopted in corrosion models is incompatible with the direct coupling with mechanics in the finite element modeling employed in modern nuclear fuel performance codes. In this study, a physically based zirconium alloy corrosion model called the Coupled-Current Charge Compensation (C4) model, which was initially developed for operating temperature conditions, has been updated to include high-temperature corrosion in order to provide additional critical information (e.g., oxygen concentration profile) under loss-of-coolant accident (LOCA) conditions—information lacking in existing empirical models. The C4 model was implemented in the MOOSE finite-element framework developed at Idaho National Laboratory, enabling it to be used in the BISON nuclear fuel performance code based on the MOOSE framework. To precisely track the different interfaces at a relatively low computational cost, the eXtended Finite Element Method (X-FEM) was applied in MOOSE. The model’s results were compared to those of existing empirical models as well as metallographic analysis of high-temperature oxidized Zircaloy-4 coupons. Oxygen diffusivities in the a and ß phases resulting from this comparison closely agree with those found in the literature. The C4 model implemented with X-FEM in MOOSE now has the capability to accurately predict oxide, oxygen-stabilized a, and prior ß phase layer growth kinetics under isothermal exposure at high temperature (1000–1500°C). Furthermore, in contrast with the empirical models, the C4 model accounts for the finite thickness of the fuel cladding. It can predict the oxygen concentration profile evolution through the whole cladding, enabling evaluation of the remaining ductile thickness—a crucial variable for modeling the mechanical behavior of the fuel cladding under LOCA. Furthermore, this implementation allows direct coupling with mechanics, at a low computing cost, using finite-element-based nuclear fuel performance codes such as BISON.

36 MATERIALS SCIENCE↗