Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “corrections”

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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Stable isotope equilibria in the dihydrogen-water-methane-ethane-propane system. Part 1: Path-integral calculations with CCSD(T) quality potentials

Isotopic compositions of alkanes are typically assumed to be kinetically controlled, but recently is has been proposed that alkanes can isotopically equilibrate for both C and H isotopes during natural gas generation. Evaluation of this requires knowledge of the isotopic equilibrium between alkanes and other common hydrogen and carbon bearing species. Here, in this study, we calculate isotopic equilibria within and between gaseous dihydrogen (H 2 ), water (H 2 O), methane (CH 4 ), ethane (C 2 H 6 ) and propane (C 3 H 8 ), including isotope fractionation among molecules, clumped isotope effects, as well as among sites of propane (i.e., the site-specific isotope effects) from 0°C to 500°C using a path-integral method paired with high-level descriptions of molecular potentials and the diagonal correction to the Born Oppenheimer approximation. While path-integral calculations with high- level CCSD(T) potentials are available for the isotopic equilibria involving methane, the path-integral calculations for ethane and propane have only been performed based on lower-level descriptions of the molecular potentials. We analyze the relative importance of various approximations that are commonly employed when isotopic equilibria are evaluated. We find that clumped isotope effects can be calculated to the same accuracy using computationally inexpensive combination of the Bigeleisen-Mayer-Urey model with the molecular potential from density functional theory. In contrast, fractionation and site preferences of both deuterium and carbon-13 benefit from the use of the higher level CCSD(T) potentials and accounting for anharmonic effects. Additionally, for fractionation and site preference of deuterium corrections to Born-Oppenheimer approximation can also be important.

03 NATURAL GAS↗
DOE OSTI2021

An investigation into the effects of cyclic strain rate on the high cycle and very high cycle fatigue behaviors of wrought and additively manufactured Inconel 718

Additive Manufacturing (AM) has increasingly been used to fabricate parts in aerospace applications, which may require service lives beyond ten-million cycles due to the imposed high loading frequencies. Understanding the very high cycle fatigue (VHCF) behavior of these additive manufactured (AM) parts is an important step towards their design and qualification processes. This study focuses on the high cycle fatigue (HCF) and VHCF behaviors of both wrought and laser beam-powder bed fusion (LB-PBF) fabricated Inconel 718 in machined/polished surface condition, emphasizing on the influence of test frequency (i.e., cyclic strain rate). Uniaxial, fully-reversed force- and stress-controlled fatigue tests were conducted utilizing a servo-hydraulic and an ultrasonic test system operating at 5 Hz and 20 kHz, respectively, on wrought as well as LB-PBF vertically and diagonally built specimens. Fatigue cracks in the majority of the specimens were found to initiate from intra-granular slip bands near or at the surface, which gives rise to strong anisotropy in fatigue resistance in LB-PBF specimens due to the presence of columnar grains along the build directions. Longer fatigue lives were obtained at 20 kHz, which was ascribed to possibly lower-than-intended stresses applied in the ultrasonic tests. The corrected stress-life fatigue data at 20 kHz were found to converge to the one obtained from conventional testing at 5 Hz, implying no effect of cyclic strain rate on the fatigue behavior of Inconel 718 regardless of the fabrication process. The findings of this work confirm the use of ultrasonic fatigue testing to expedite generation of AM materials data to keep up with the current demand; however, the applied stress may need to be corrected.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Innovative approach to counterfeit and noncompliant refrigerant detection: A cost-effective, portable solution

The increasing prevalence of counterfeit and incompatible refrigerants presents significant risks to Heating, Ventilation, Air Conditioning, and Refrigeration (HVAC&R) systems, including compromised equipment performance, safety hazards, and non-compliance. This article details the development of a novel, cost-effective, and portable detection device designed to accurately verify refrigerants. The device utilizes a controlled gas sampling and analysis system within a sealed chamber, ensuring precise measurements while maintaining safety through a purging mechanism. The system features a high-sensitivity sensor integrated with an onboard control module that analyzes gas composition in real-time, providing feedback within a 2-minute duration. Laboratory validation demonstrated the device’s high accuracy (>95 % based on correct identification of compliant vs. non-compliant blends) in detecting unauthorized refrigerant blends. The projected cost of the product stands at ∼ $150, based on the retail pricing of individual components. Laboratory validation demonstrated the device’s high accuracy (>95 % for composition identification, 100 % rejection of tested counterfeit/incorrect blends) in detecting unauthorized refrigerant blends with a response time <2 min. The device correctly identified authentic R-454A/B/C blends and reliably rejected R-407F and closely related counterfeit mixtures. Key advantages include affordability, ease of use, rapid response time, and compatibility with a wide range of refrigerants. This solution supports compliance with regulatory frameworks, enhances safety in HVAC&R operations, and mitigates the risks associated with counterfeit refrigerants.

Counterfeit refrigerants↗
DOE OSTI2022

Risk matrix for legacy wells within the Area of Review (AoR) of Carbon Capture & Storage (CCS) projects

The success of CCS depends on the capacity, injectivity, and confinement by the storage medium. Thousands of wells drilled over the past century with the intention to find the trapped oil and gas may penetrate the containment seals. These wells may provide leakage pathways for the CO 2 to escape and contaminate the underground sources of drinking water (USDW) or reach the surface in the worst-case scenario. Identifying the risky wells penetrating the containment seals and predicting their current as well as future well integrity is the most challenging task when limited data is available. Here, this paper proposes a unique methodology for risk assessment of the wells penetrating the containment seals based on the proximity of these wells from the proposed injection location, mechanical integrity, and accessibility of these wells over the lifecycle of the CCS project. This helps in identifying the wells which need immediate attention from the wells that need little to no attention. It also highlights the corrective actions necessary for the success of the CCS project as well as help estimate the approximate cost required to perform the corrective actions. This methodology focuses on all wells (producers, injectors, orphan, abandoned, water, stratigraphic, etc) while the majority of the studies found in literature focused on wells with sustained casing pressure (SCP) reports and cement bond logs (CBL). The proposed risk matrix, if applied to future CCS projects across the globe, will uniformly categorize the wells within the Area of Review (AoR).

58 GEOSCIENCES↗
DOE OSTI2024

Corrigendum to “Cross-sectoral assessment of CO2 capture from U.S. industrial flue gases for fuels and chemicals manufacture” [International Journal of Greenhouse Gas Control 135 (2024) 1-20 / 104137]

The authors regret the inaccuracy in the vertical axis title of Fig. 4 and the distortion in the legend of Fig. 13. Corrections have been made to the vertical axis title and legend of Figs. 4 and 13, respectively. These changes do not affect the analysis, calculations, or results in any way. The authors apologize for any inconvenience caused. The corrected figures are as follows:

Zuberi, M Jibran S↗
DOE OSTI2021

Comprehensive framework for data-driven model form discovery of the closure laws in thermal-hydraulics codes

The two-phase two-fluid model is a basis of many thermal-hydraulics codes used in design, licensing, and safety considerations of nuclear power plants. Thermal-hydraulics codes rely on the closure laws to close the system of conservation equations and describe the interactions between phases. These laws, derived from years of experimental investigations, are semi-empirical correlations that lack generality and have a limited range of applicability. Increase of computational power, availability of new experiments, and development of high-fidelity simulations has increased the number of validation data. The discrepancies between the code predictions and the validation data are a great source of knowledge. Missing physics that are not included in the model but are important for the considered phenomena can be discovered by propagating the information from the experimental results through the model. Furthermore, physics-discovered data-driven model form (P3DM) methodology integrates available integral effect tests and separate effects tests to determine the necessary corrections to the model form of the closure laws. In contrast to existing calibration techniques, the methodology modifies the functional form of the closure laws. Based on the functional form of the correction, the missing physics that were not included in the original model can be discovered. The methodology provides the alternative to the machine learning approach, in which the model is discovered in the form of the intractable black-box relation. In this work, the methodology was applied to the CTF subchannel code to improve the prediction of the two-phase flow phenomena.

42 ENGINEERING↗
DOE OSTI2023

Evaluation of the ring tension test (RTT) for robust determination of material strengths

The ring tension test (RTT) is an experimental method for determining mechanical behavior in a material’s circumferential or hoop direction. It is a crucial test for testing anisotropic materials with tube geometry, such as nuclear fuel cladding or irradiated pipes. Several RTT configurations exist, each with their own advantages and disadvantages. However, this test is significantly more complex than traditional tensile testing and can be especially sensitive to small differences and inconsistencies in the test setup and geometry, ultimately affecting the derived mechanical properties. Previous research has focused on method development, and little work has been done on understanding the subtle differences between an ideal test and experiments, specifically when the tests are performed on highly irradiated materials in hot cells. In this work, a finite element-based investigation of the RTT is conducted to fill this gap. The two most used test configurations are investigated, comparing their ability to determine accurate material strengths through plastic deformation. Several non-ideal conditions and uncontrollable effects which are likely to occur during experimental testing such as machining tolerances, variations of specimen geometry from nominal dimensions, rotation of specimens and fixturing, and other test setup discrepancies are studied. The sensitivity of measured strengths to these conditions is presented. A mechanics-based approach to describing and correcting raw data to determine actual strengths is also included for one of the configurations, resulting in a robust correction method with highly accurate material strength measurements. Finally, based on these analyses, the hemicylindrical mandrel configuration is recommended with a gauge region oriented at a 45° angle.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

DPM: A deep learning PDE augmentation method with application to large-eddy simulation

A framework is introduced that leverages known physics to reduce overfitting in machine learning for scientific applications. The partial differential equation (PDE) that expresses the physics is augmented with a neural network that uses available data to learn a description of the corresponding unknown or unrepresented physics. Training within this combined system corrects for missing, unknown, or erroneously represented physics, including discretization errors associated with the PDE's numerical solution. For optimization of the network within the PDE, an adjoint PDE is solved to provide high-dimensional gradients, and a stochastic adjoint method (SAM) further accelerates training. Additionally, the approach is demonstrated for large-eddy simulation (LES) of turbulence. High-fidelity direct numerical simulations (DNS) of decaying isotropic turbulence provide the training data used to learn sub-filter-scale closures for the filtered Navier–Stokes equations. Out-of-sample comparisons show that the deep learning PDE method outperforms widely-used models, even for filter sizes so large that they become qualitatively incorrect. It also significantly outperforms the same neural network when a priori trained based on simple data mismatch, not accounting for the full PDE. Measures of discretization errors, which are well-known to be consequential in LES, point to the importance of the unified training formulation's design, which without modification corrects for them. For comparable accuracy, simulation runtime is significantly reduced. A relaxation of the typical discrete enforcement of the divergence-free constraint in the solver is also successful, instead allowing the DPM to approximately enforce incompressibility physics. Since the training loss function is not restricted to correspond directly to the closure to be learned, training can incorporate diverse data, including experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Bayesian sparse learning with preconditioned stochastic gradient MCMC and its applications

Deep neural networks have been successfully employed in an extensive variety of research areas, including solving partial differential equations. Despite its significant success, there are some challenges in effectively training DNN, such as avoiding overfitting in over-parameterized DNNs and accelerating the optimization in DNNs with pathological curvature. Here, we propose a Bayesian type sparse deep learning algorithm. The algorithm utilizes a set of spike-and-slab priors for the parameters in the deep neural network. The hierarchical Bayesian mixture will be trained using an adaptive empirical method. That is, one will alternatively sample from the posterior using preconditioned stochastic gradient Langevin Dynamics (PSGLD), and optimize the latent variables via stochastic approximation. The sparsity of the network is achieved while optimizing the hyperparameters with adaptive searching and penalizing. A popular SG-MCMC approach is Stochastic gradient Langevin dynamics (SGLD). However, considering the complex geometry in the model parameter space in nonconvex learning, updating parameters using a universal step size in each component as in SGLD may cause slow mixing. To address this issue, we apply a computationally manageable preconditioner in the updating rule, which provides a step-size parameter to adapt to local geometric properties. Moreover, by smoothly optimizing the hyperparameter in the preconditioning matrix, our proposed algorithm ensures a decreasing bias, which is introduced by ignoring the correction term in the preconditioned SGLD. According to the existing theoretical framework, we show that the proposed algorithm can asymptotically converge to the correct distribution with a controllable bias under mild conditions. Numerical tests are performed on both synthetic regression problems and learning solutions of elliptic PDE, which demonstrate the accuracy and efficiency of the present work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A positivity-preserving and conservative intersection-distribution-based remapping algorithm for staggered ALE hydrodynamics on arbitrary meshes

In this study, we introduce new intersection-distribution-based remapping tools for indirect staggered arbitrary Lagrangian-Eulerian (ALE) simulations of multi-material shock hydrodynamics on arbitrary meshes. In addition to conserving momentum and total energy, the three-stage remapper proposed in this work preserves non-negativity of the internal energy. At the first stage, we construct slope-limited piecewise-linear reconstructions of all conserved quantities on zones of the source mesh and perform intersection-based remap to obtain bound-preserving zonal quantities on the target mesh. At the second stage, we define bound-preserving nodal quantities of the staggered ALE discretization as convex combinations of corner quantities. The nodal internal energy is corrected in a way which keeps it non-negative, while providing exact conservation of total energy. At the final stage, we distribute the non-negative nodal internal energy to corners, zones and materials using non-negative weights. Proofs of positivity preservation are provided for each stage. This work is a natural extension of our paper [14] in which a similar intersection-distribution-based remapping procedure was employed. The original version used a nodal kinetic energy fix which did not provably ensure positivity preservation for the zonal internal energy after the final distribution stage. The new algorithm cures this potential drawback by using ‘coordinated’ limiters for piecewise-linear reconstructions, remapping the internal energy to nodes and correcting it before redistribution. The effectiveness of the new nodal fix is illustrated by numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

An explicit, energy-conserving particle-in-cell scheme

We present an explicit temporal discretization of particle-in-cell schemes for the non-relativistic Vlasov equation that results in exact energy conservation when combined with an appropriate spatial discretization. The scheme is inspired by a simple, second-order explicit scheme that conserves energy exactly in the Eulerian context. We show that direct translation to particle-in-cell does not result in strict conservation, but derive a simple correction based on an analytically solvable optimization problem that recovers conservation. While this optimization problem is not guaranteed to have a real solution for every particle, we provide a correction that makes imaginary values extremely rare and still admits $\mathcal{O}$(10 –12 ) fractional errors in energy for practical simulation parameters. We present the scheme in both electrostatic – where we use the Ampère formulation – and electromagnetic contexts. With an electromagnetic field solve, the field update is most naturally linearly implicit, but the more computationally intensive particle update remains fully explicit. Here, we also show how the scheme can be extended to use the fully explicit leapfrog and pseudospectral analytic time-domain (PSATD) field solvers. The scheme is tested on standard kinetic plasma problems, confirming its conservation properties.

Energy conservation↗
DOE OSTI2025

Multi-material ALE remap with interface sharpening using high-order matrix-free finite element methods

The arbitrary Lagrangian-Eulerian (ALE) technique involves remapping field quantities from a Lagrangian mesh to an optimized mesh in a conservative, accurate and bounds-preserving manner. For methods based on arbitrary order finite elements, as described in a reference, material volume fractions are advected in pseudo-time using flux-corrected transport (FCT) without any form of interface reconstruction. In practice, this can lead to excessive propagation of small volume fractions throughout the domain. In addition, this method requires assembly of a global advection matrix to compute the bounds-preserving low-order FCT solution. In this work, we introduce a new approach for ALE remap using a high-order matrix-free technique which incorporates a flux modification to sharpen material interfaces in a conservative manner. Our approach begins with computing a bounds-preserving low-order solution to the ALE remap equations at the element level. We then compute a sharp interface solution (not guaranteed to be bounds-preserving) which comes from solving an augmented version of the ALE remap equations with a conservative flux modification which acts to sharpen material volume fractions based on their gradients and transport directions. Using the sharp interface solution, we make global corrections to the bounds-preserving solution while maintaining preservation of bounds. By blending with the sharpened solution at the global level we are able to globally conserve mass without hindering the remap pseudo-time step. This new interface-aware ALE remap method is based entirely on partial assembly techniques where globally assembled matrix operators are no longer needed, resulting in a globally matrix-free FCT method for multi-material, multi-field ALE remap with high performance on GPU architectures. We present results of our new remap method on 1D, 2D and 3D benchmarks and describe the algorithmic tailoring for GPU architectures that was developed.

Vargas, Arturo [Lawrence Livermore National Labora↗
DOE OSTI2025

Climate-extreme modeling framework for sustainable flood management in the Arabian Peninsula

Evaluating extreme precipitation events (EPEs) is essential for building climate-resilient water management strategies, but it remains a major challenge in ungauged basins. Using the 26,070 km 2 Wadi al-Rummah basin in central Saudi Arabia as a case study, we developed an alternative, reliable, cost-effective satellite-based framework that combines empirically derived EPE thresholds, imagery-calibrated 2D hydrodynamic modeling, GRACE water-storage diagnostics, and bias-corrected CMIP6 projections to assess flood hazards and recharge potential under current and future climate scenarios in ungauged basins. The integrated approach and the resulting findings followed four key steps: (1) Identified a 22.5 mm EPE threshold, the 80th percentile of 3-day GPM/IMERG rainfall (2000–2024), aligned with flood-triggering events (Nov 2018: 23–28 mm; Apr, 2023: 42 mm); (2) Developed and calibrated a RiverFlow2D model using Sentinel-2 and PlanetScope imagery for the November 2018 flood, accurately reproducing flood depth and extent (RMSE ≤0.31 m; fuzzy-Dice ≥0.91), and estimating runoff (41 %), infiltration (25 %), and evaporation (34 %); (3) Independently validated the model with the April 2023 event (RMSE ≤0.35 m; fuzzy-Dice ≥0.86); (4) Conducted climate projections (2025–2100) from five bias-corrected NEX-GDDP CMIP6 models that revealed a 34 % increase in EPE intensity under SSP2-4.5 and 48 % under SSP5-8.5 scenarios, relative to 20th-century baselines. Our findings indicate that while intensifying extremes in the 21st century increase flood risk, the results highlight the potential for episodic recharge if effective retention strategies are employed, and offer a transferable model for climate-informed planning in data-scarce arid regions.

CMIP6↗
DOE OSTI2023

Examining the effects of soil entrainment during nuclear cloud rise on fallout predictions using a multiscale atmospheric modeling framework

Current operational models for nuclear cloud rise over land were developed and validated using observations from shallow-buried or surface detonations, where lofted soil quickly mixed with fission products from the detonation. These models poorly predict fallout from elevated detonations near the fallout-free height of burst (FFHOB), where interactions with the ground are limited and the mixing of fission products and lofted soil is incomplete. Fallout-free is a misnomer at this HOB, as fallout was observed in these cases, but was below the levels of concern, especially off-grounds of the nuclear test site. To correctly characterize and model fallout from detonations near the FFHOB, models must be developed which can capture the stratified nature of the particle and activity-size distributions within the cloud. Previously, it was shown that the Weather Research and Forecasting (WRF) model can accurately simulate nuclear cloud rise for airbursts with little to no ground interactions (Arthur et al., 2021). That work is expanded here by (1) using a radiation-hydrodynamics code to improve the fireball initialization in WRF, (2) further developing an aerosol package from WRF-Chem to simulate lofted soil, and (3) combining the WRF cloud rise simulations with the operational models used at the National Atmospheric Release Advisory Center (NARAC) for fallout modeling. Using this combination of codes, the Upshot-Knothole Grable detonation, which was just below the FFHOB, is simulated from seconds after detonation through cloud rise and fallout, and results are compared to historical test data. Here, the results show improved prediction of dose rate and highlight the need to correctly characterize the entrainment of material into the cloud and the subsequent mixing of fission products with entrained material.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

A improved equation of state for Xe gas bubbles in gamma U-Mo fuels

A monolithic fuel design based on a U-Mo alloy has been selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). An issue with U-Mo monolithic fuel is the large amount of swelling that takes place during operation. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic properties into mesoscale and continuum level fuel performance modeling codes. However, the thermodynamic properties of the fission gas bubbles (such as the relationship among bubble size, equilibrium Xe concentration, and bubble pressure) are not well known. This work studies Xe bubbles in ?U-Mo from a diameter of 3 nm up to 8.5 nm and from 400 K up to 700 K. The energetic relationship of Xe bubbles with regard to voids and Xe substitutional atoms is described. The transition is also determined for when a bubble becomes over-pressurized. Finally, an equation of state is fit to the pressure as a function of molar volume and temperature.

Beeler, Benjamin↗
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 OSTI2022

Anisotropic temperature-dependent elastic constants and thermal conductivities of TRISO particle coatings

Tristructural isotropic (TRISO) particles show great promise as a candidate fuel for use in several next-generation high-temperature nuclear reactor designs due to their structural integrity and fuel performance at high temperatures and burnups. Computational codes exist that can simulate TRISO fuel performance characteristics and failure probabilities under extreme conditions which require knowledge of the TRISO coatings’ thermophysical properties. The thermophysical descriptions of the TRISO particle's layers (i.e., buffer, pyrolytic carbon, and silicon carbide) currently used in fuel performance codes, however, assume that many of these properties are constant with respect to temperature or texture. In this paper, we obtain the full elastic stiffness tensors of the carbon and silicon carbide layers, which have transversely isotropic symmetry. Additionally, the calculated elastic properties of each of the coating layers are in remarkable agreement with the current models at room temperature and correct orientations. Additionally, the textured 3C-SiC layer was found to exhibit novel auxetic behavior above 1500 °C. The anisotropic temperature-dependent thermal conductivities of the carbon layers were calculated using acoustical Grüneisen-Debye theory which are in excellent agreement with current models at room temperature and correct orientations. These texture- and temperature-dependent relationships can be incorporated into the thermophysical description of TRISO particles in order to more accurately model fuel performance and failure probabilities under extreme conditions in forthcoming high-fidelity computational simulations.

36 MATERIALS SCIENCE↗