SEARCH · Engineering Papers

Results for “work”

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.

SourceResearch connectionFrom yearThrough year

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

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Benchmark of numerical modeling approaches on the systematic performance evaluation of wave energy converters

Different numerical modeling methods have been developed and applied to evaluate a variety of performance indicators of wave energy converters (WECs), including the power performance, structural loads, levelized cost of energy, etc. Based on the modeling fidelity, the commonly used numerical modeling approaches can be classified as linear modeling, weakly nonlinear modeling and fully nonlinear modeling approaches. Each method differs in accuracy and computational efficiency, making them suitable for different stages of WEC design. However, the selection of modeling approach could significantly impact evaluation outcomes. For instance, simplified linear models may underestimate structural loads or overestimate energy production in some operational conditions, potentially leading to less cost-effective designs. Given the widespread utilization of these models, it is essential to understand the uncertainties brought by them in performance evaluations. This work is dedicated to benchmarking different linear-potential-flow-based numerical models for evaluating the systematic performance of WECs. Three representative numerical modeling approaches are considered in this work, including linear frequency-domain modeling, statistically linearized spectral-domain modeling and Cummins equation-based nonlinear time-domain modeling. A generic point absorber WEC is considered as the research reference in this work, and different sea sites are taken into account. The numerical models are utilized to predict critical performance indicators, including power performance, the annual energy production, the capacity factor, the levelized cost of energy and the PTO fatigue loads. By comparing the results, this work identifies the uncertainties associated with different modeling approaches in evaluating WEC performance.

Fatigue↗

DOE OSTI2025

Degrees of rate control and AutoDiff-driven direct sensitivity analysis in heterogeneous catalysis

Despite the wide application and benefits of the degree of rate control (DRC) analysis, several details remain argued, particularly about the conservation of DRCs at transient (TR) and steady-state (SS) conditions, especially for complex reaction networks. This work argues that previous proofs about the conservation properties of DRCs have been incomplete, and we provide new mathematical proofs at TR and SS conditions. In addition, we use both analytical (automatic differentiation) and numerical (finite difference) approaches to compute DRCs for the case study of ethane hydrogenolysis (EH) over Pt(111). This work confirms that at both TR and SS conditions, the sum of all DRCs, i.e., sum of the degrees of kinetic (DKRC) and thermodynamic rate control (DTRC), is conserved at zero. At SS conditions, the sum of DKRC is conserved at 1 while the sum of DTRC is conserved at −1. In corroboration of previous works, we show that the DTRC for any adsorbate at SS is equal to the product of the species coverage and a constant. In contrast, at TR conditions, the individual sums of both DTRC and DKRC are not conserved and can be any real number, with potential implications for the novel field of dynamic catalysis. Finally, we show that the conventional finite difference (FD) approach, only useful at SS, is prone to inaccuracy and very sensitive to the value of the differential change applied. The optimal differential value also varies significantly with system and rate definition. Consequently, we describe and illustrate in this work the application of the automatic differentiation (AD) approach for the more accurate determination of DRCs at both TR and SS conditions.

Automatic differentiation↗

DOE OSTI2021

Bond of steel-mortar interface interfered by stray current

Highlights: • The focus of this work is to present test results on the bond of steel-mortar interface undergoing stray current; • The bond strength is correlated to the electrochemical response of the steel-mortar interface; • The effects of duration of curing and starting point of stray current are also investigated; • This work also differentiates the effects of stray current on steel-mortar bond versus anodic polarization. The focus of this work is to present test results on the bond of steel-mortar interface undergoing stray current. The bond strength, derived by pull-out tests, is correlated to the electrochemical response of the steel rebar and the properties of the mortar bulk matrix. The effects of curing regimes (in terms of duration of curing) and starting point of stray current are also investigated. It is found that stray current exerts bond degradation of the steel-mortar interface in all investigated cases, irrespective of the presence or absence of a corrodent (Cl{sup −}) in the external medium. For the ease of operation in lab tests, the stray current is generally simulated by anodic polarization, although fundamentally, the stray current effect on the steel surface is composed of both anodic and cathodic polarizations. Hence this work also differentiates the effects of stray current on steel-mortar bond, versus the effects of anodic polarization.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Neural-based time series forecasting of loss of coolant accidents in nuclear power plants

During the last few years, deep learning in neural networks has demonstrated impressive successes in the areas of computer vision, speech and image recognition, text generation, and many others. However, sensitive engineering areas such as nuclear engineering benefited less from these efficient techniques. In this work, deep learning expert systems are utilized to model and predict time series progression of a design-basis nuclear accident, featuring a loss of coolant accident. Two major findings are accomplished in this work. First, the ability to train expert systems with high accuracy, which could help nuclear power plant operators to figure out plant responses during the accident. Second, building fast, efficient, and accurate deep models to simulate nuclear phenomena, which could be valuable to nuclear computational science. In this work, large amount of time series data is obtained from simulation tools by simulating different conditions of the base-case/nominal accident scenario. Four critical outputs/responses are monitored during the accident (e.g. temperature, pressure, break flow rate, water level). Two approaches are adopted in this work. The first approach is to use feedforward deep neural networks (DNN) to fit all time steps and outputs in a single model. The second approach is to use long short-term memory (LSTM) to fit all time steps together for each reactor response separately. Both DNN and LSTM demonstrate very good performance in predicting the test and base-case scenarios, with accuracy as low as 92% and as high as 99%, where these test scenarios are unknown to the expert systems and are not included in the model training. In addition, both approaches demonstrate a significant reduction in computational costs, as the deep expert system is able to accurately predict the accident 100,000 times faster than the original simulation tool. Given sufficient data, the methodology adopted in this study demonstrates that DNN/LSTM expert systems can be used as a decision support system to model advanced time series phenomena within nuclear power plants with high accuracy and negligible computational costs.

42 ENGINEERING↗

DOE OSTI2022

Simulated CO 2 storage efficiency factors for saline formations of various lithologies and depositional environments using new experimental relative permeability data

Saline formations are attractive geologic reservoirs for permanent carbon dioxide (CO 2 ) storage. Here, the U.S. Department of Energy's National Energy Technology Laboratory (DOE-NETL) has worked to develop and refine methods and tools for the calculation of CO 2 storage potential in subsurface reservoirs. DOE-NETL's CO 2 -SCREEN provides an online tool for executing these storage methods. CO 2 storage efficiency terms are input parameters in DOE-NETL's methods and equations embedded in the CO 2 -SCREEN, which assesses pore space available for CO 2 storage. In this work, a modeling workflow was initiated to refine two CO 2 storage efficiency terms - volumetric displacement (E V ) and microscopic displacement (E d ). The models are based on new experimental relative permeability data that are specific to homogenous lithology and depositional environments of key subsurface saline formations targeted for CO 2 storage. In future work, heterogenous features will be added to this initial modeling effort to update efficiency factors as described in DOE-NETL's methods and CO 2 -SCREEN tool. E V accounts for the volume utilized in the reservoir under the areal plume, while E d accounts for saturation values in the plume to assess efficiency of CO 2 storage at the pore scale. The results of this work are significant in that prior values were based on a limited geologically non-specific relative permeability data set that were collected prior to 2009. Specifically, we applied numerical simulations using TOUGH3 models to update CO 2 storage efficiency values for supercritical CO 2 injection into brine-saturated reservoirs for three lithologies (clastics, limestone, dolomite) and six depositional environments (Marginal Marine, Strand Plain, Deltaic Complex Fluvial, Aeolian, Shallow Marine, and Reef) that have a high potential for geologic CO 2 storage. Experimental relative permeability data in cores from these environments were utilized in the models with corresponding rock type/sedimentary environment. Results of this study showed that dolomite followed by limestone generated higher ranges of storage efficiency compared to clastics. The updated values provided a tighter efficiency range for clastics, lower P 10 but higher P 90 range for limestone, and higher P 10 and P 90 for dolomite. In general, tighter reservoirs with relatively low permeability and porosity were associated with higher E V and E d , showing efficient reservoir and pore utilization in these scenarios. High reservoir pressure and temperature associated with increasing depth increased the E V , and high CO 2 injection rates resulted in increases in E V and E d , while the impact of permeability anisotropy was minimal after the 30-year injection period.

03 NATURAL GAS↗

DOE OSTI2020

Modeling of hydrogen liquefaction using magnetocaloric cycles with permanent magnets

Hydrogen (H 2 ) is promising alternative to replace fossil fuels, but its transport and storage has been challenging. As H 2 fuel cell vehicles are gaining traction, the infrastructure for storing large amounts of liquid H 2 is needed. However, liquid H 2 would suffer from boil-off loss, and traditional vapor compression refrigeration systems would not be able to economically recover the lost H 2 due to the low efficiencies at cryogenic temperature. Magnetocaloric (MC) refrigeration systems could possess much higher coefficient of performance (COP) at cryogenic temperature compared to the vapor compression ones. Previous work on cryogenic MC systems, however, have only focused on large scale applications which use superconducting magnets to provide a large magnetic field but are prohibitively expensive to operate for small scale applications, such as that of a H 2 refilling station. In this work, we model the performance of a MC refrigeration cycle using 1-Tesla permanent magnets for H 2 liquefaction, with the objective of cooling H 2 from 80 K (using liquid nitrogen as the heat sink) to 20 K (boiling point of hydrogen). We evaluate main performance metrics including the total work input to the refrigeration system, COP, total MCM mass in the system, and total volume of the permanent magnets, etc. Our modeling results indicate that such a permanent magnet-based MC cooling system is feasible for small-scale H 2 liquefaction, with projected COP values significantly higher than those of vapor compression systems. In conclusion, this work provides design guidelines for future experimental efforts on permanent magnet MC cooling systems for cryogenic cooling.

08 HYDROGEN↗

DOE OSTI2025

Fuel reid vapor pressure level and ethanol content on stochastic preignition, effects at steady and unsteady engine operation

The present work investigates relations between fuel Reid vapor pressure (RVP) and biofuel (ethanol) content on stochastic preignition (SPI) at both sustained steady-state engine operation and following load transients. This work stems from in-field observations that automotive original equipment manufacturers have observed consistent seasonal increases in United States customer drivability complaints and warranty claims during September and October where SPI is suspected to be responsible. The seasonal timing of these events coincides with the United States seasonal fuel property changeover initiating on September 15 each year, where fuel RVP increases. To explore potential linkage between fuel RVP and SPI the present study employs engine SPI experiments coupled with laboratory spray measurements of fuels with RVPs of 8, 12, and 16 psi in both E10 (10% ethanol) and E25 (25% ethanol) fuels. Engine results are partitioned into fuel RVP and ethanol content effects on SPI in steady-state, sustained high-load engine operation and unsteady-state low- to high-load transitions, where off-engine spray vessel patternation and tip penetration results help to elucidate the observed fuel effects on SPI. A boosted direct-injected, spark-ignition engine was fueled with three market relevant E10 and E25 fuels with RVPs of 8, 12, and 16 to characterize the interplay between winter fuels and abnormal combustion behavior. The steady-state work shows that for high-load, steady-state engine operation, SPI is directly linked to fuel retention, which was found to be dependent on fuel distillation. The unsteady-state engine operation work shows that following low-to high-load transitions, SPI can occur from a memory of fuel property effects at low-load operation. Specifically, the fuel RVP effect on fuel spray collapse at low loads was found to correlate with SPI with a more than 95% confidence interval following low- to high-engine-load transitions. Results suggest that fuel-wall impingement at low-load operation could carry over into high-load transitions and generate SPI events following low- to high-load transitions.

09 BIOMASS FUELS↗

DOE OSTI2020

Determining the impacts of venoarterial extracorporeal membrane oxygenation on cerebral oxygenation using a one-dimensional blood flow simulator

Venoarterial extracorporeal membrane oxygenation (VA-ECMO) is a mechanical system that provides rapid and short-term support for patients with cardiac failure. In many patients, pulmonary function is also impaired, resulting in poorly-oxygenated cardiac outflow competing against well-oxygenated VA-ECMO outflow, a condition known as North-South syndrome. North-South syndrome is a primary concern because of its potential to cause cerebral hypoxia, which has a critical influence on neurological complications often seen in this patient population. In order to reduce ischemic neurological complications, it is important to understand how clinical decisions regarding VA-ECMO parameters influence blood oxygenation. In our work, we studied the impacts of flow rate and cannulation site on oxygenation using a one-dimensional (1D) model to simulate blood flow. Our model was initially tested by comparing blood flow results to those observed from experimental work in VA-ECMO patients. The 1D model was combined with a two-phase flow model to simulate oxygenation. Additionally, the influence of various other clinician-tunable parameters on oxygenation in the common carotid arteries (CCAs) were tested, including, blood viscosity, cannula position within the insertion artery, heart rate, and systemic vascular resistance (SVR), as well as geometrical changes such as arterial radius and length. Our results indicated that blood oxygenation to the brain strongly depended on the cannula insertion site and the VA-ECMO flow rate with a weaker but potentially significant dependence on arterial radius. During femoral cannulation, VA-ECMO flow rates greater than ~4.9 L/min were needed to perfuse the CCAs. However, axillary and central cannulation began to perfuse the CCAs at significantly lower flow (~1 L/min). These results may help explain the incidence of cerebral hypoxia in this patient population and the common need to change cannulation strategies during treatment to address this clinical problem. While this work describes patient-averaged results, determining these relationships between VA-ECMO parameters and cerebral hypoxia is an important step towards future work to develop patient-specific models that clinicians can use to improve outcomes.

1D blood flow↗

DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Investigation of γ - ( U , Z r ) structural properties and its interfacial properties with liquid sodium using ab initio molecular dynamics

In this study, the elastic properties, structural parameters, sound velocity, and Debye temperature of γ–(U,Zr) were computed using ab initio molecular dynamics (AIMD) at temperatures between 1000 K and 1400 K and for Zr content between 0 at. % and 100 at. %. UZr is used as a metallic fuel for Sodium Fast Reactors (SFRs). The study of the mechanical and thermal behavior of these alloys leads to a better data-informed fuel design. The bulk modulus, shear modulus, Young's modulus, and Poisson's ratio were calculated from the elastic constants and their dependence on Zr content and temperature was investigated, comparing the results with previous computational work and the available experimental data in the literature. Interfacial properties between UZr (up to 32 at. % which typically exists in nuclear fuel) and liquid sodium are also of interest due to the presence of a sodium bond between the fuel and the cladding in metallic nuclear fuel. The interfacial energy between γ–(U,Zr) and liquid sodium, the surface tension of liquid sodium, and the work of adhesion were computed at different temperatures and Zr concentrations. It was demonstrated that γ–(U,Zr) is completely wetted by liquid sodium at all the investigated temperatures and Zr concentrations. Finally, this work provides the basis for the determination of interfacial resistances in SFRs and their implementation into heat transfer fuel performance simulations, which will be the subject of future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

New Synthesis Strategies to Improve Co-Free LiNi0.5Mn0.5O2 Cathodes: Early Transition Metal d0 Dopants and Manganese Pyrophosphate Coating

In this work, we report solution-based doping and coating strategies to improve the electrochemical performance of the Co-free layered oxide cathode LiNi0.5Mn0.5O2 (LNMO). Small amounts of d0 dopants (e.g., Mo6+and Ti4+, 0.5-1 at. %) increase the cathode’s specific capacity, cycling stability, and rate capability. More specifically, a Mo-doped cathode with the nominal composition LiNi0.495Mn0.495Mo0.01O2 achieves a high reversible capacity of 180 mAh/g at 20 mA/g with good retention at higher rates (e.g., 120 mAh/g at 100 mA/g). Effects of 1 at.% Mo dopant on the cathode structure were studied using a suite of characterization tools including X-ray diffraction (XRD), Raman spectroscopy, and transmission electron microscopy (TEM). These measurements demonstrate that Mo6+ reduces Li+/Ni2+cation mixing and mitigates undesirable phase transformations near the cathode surface during cycling. This work also reports the use of an inorganic Mn2P2O7 coating which enhances cycling stability, presumably through formation of a stable cathode electrolyte interphase (CEI) layer. Overall, the synthesis approaches reported herein are quite general and can potentially be expanded to other high voltage LIB cathodes.

Co-free cathodes, d0 cation, pyrophosphate coating↗

DOE OSTI2022

Monitoring and flaw detection during wire-based directed energy deposition using in-situ acoustic sensing and wavelet graph signal analysis

The goal of this work is to detect flaw formation in the wire-based directed energy deposition (W-DED) process using in-situ sensor data. The W-DED studied in this work is analogous to metal inert gas electric arc welding. The adoption of W-DED in industry is limited because the process is susceptible to stochastic and environmental disturbances that cause instabilities in the electric arc, eventually leading to flaw formation, such as porosity and suboptimal geometric integrity. Moreover, due to the large size of W-DED parts, it is difficult to detect flaws post-process using non-destructive techniques, such as X-ray computed tomography. Accordingly, the objective of this work is to detect flaw formation in W-DED parts using data acquired from an acoustic (sound) sensor installed near the electric arc. To realize this objective, we develop and apply a novel wavelet integrated graph theory approach. The approach extracts a single feature called graph Laplacian Fiedler number from the noise-contaminated acoustic sensor data, which is subsequently tracked in a statistical control chart. Using this approach, the onset of various types of flaws are detected with a false alarm rate less-than 2%. This work demonstrates the potential of using advanced data analytics for in-situ monitoring of W-DED.

42 ENGINEERING↗

DOE OSTI2026

Improving creep resistance of Al-0.27Zr-0.08Sn (wt%) via cold swaging while maintaining high electrical conductivity

This study investigates the effect of mechanical cold work (0, 50, and 90% cross-sectional area reduction via swaging) on the creep properties of cast aluminum alloys with high electrical conductivity including high-purity Al (HP-Al), Al-0.15Zr (Al-Zr), and Al-0.27Zr-0.08Sn (Al-Zr-Sn, wt%). The Al-Zr alloy is strengthened by Zr in solid solution while the Al-Zr-Sn alloy is additionally strengthened by Al3Zr nanoprecipitates (8 nm diameter). After 90% cold swaging, the alloys exhibit grains elongated along the swaging direction with their widths ranging between 8 and 152 μm perpendicular to the swaging direction. Creep testing at 200 °C shows that the minimum creep rate of all 90% swaged alloys shows high sensitivity to stress, which can be modeled via a threshold stress σth. Increasing swaging magnitude from 0 to 50 to 90% in Al-Zr-Sn improves the creep resistance without reducing electrical conductivity. The formation of subgrains, increased dislocation density, and columnar grain structure in swaged alloys all enhance creep resistance. HP-Al swaged to 90% exhibits much lower creep resistance (σth = 18 MPa) due to subgrain coarsening and the absence of precipitates compared to Al-Zr (σth = 40 MPa) and Al-Zr-Sn (σth = 40 MPa), which maintain stable subgrain structures and benefit from solid solution strengthening and Al3Zr nanoprecipitates, respectively. Although Al3Zr precipitates are known to stabilize the deformed microstructure by pinning grain-boundaries, this work demonstrates that Zr in solid solution alone can effectively stabilize the deformed microstructure resulting in enhanced creep resistance. This effect of Zr solute on stabilizing grain boundaries for creep performance was not previously well-defined in the literature. Despite similar creep performances of Al-Zr and Al-Zr-Sn, the latter shows 70% higher room-temperature microhardness and slightly improved electrical conductivity (56%IACS (International Annealed Copper Standard) vs. 57%IACS, respectively). This work provides new insights into creep mechanisms of cold-worked Al–Zr alloys that will guide the design of heat-resistant Al conductors for high-demand applications.

Coello ramirez, Ismael [Northwestern University,]↗

DOE OSTI2021

A General Strategy for Batch Development of High-Performance and Cost-Effective Sodium Layered Cathodes

High-performance and low-cost transition metal (TM) layered oxides using earth abundant elements are promising cathodes for Na-ion batteries. However, it is challenging to obtain desired materials because the large Na size, different Na occupations and various layer stacking sequences multiply the complication in determining the structure of a given composition and exacerbate uncertainty to the structure-property correlation. In this work, we use the attainment of desired NaxMnyNiZTM1-y-zO2-based cathode materials as an example to demonstrate a general roadmap for batch development of sodium layered cathodes towards practical applications. A synthesis phase diagram of NaxMnyNi1-yO2 was created for pre-screening and rational selection of the platform material of P2/O3-structured Na0.85Mn0.6Ni0.4O2. Cationic potential was leveraged in elemental substitution to further promote the material structural stability and electrochemical performance. Several cost-effective O3 and P2/O3 hybrid cathode materials have been obtained, all of which demonstrate excellent performance. In particular, the Na0.85Mn0.5Ni0.4Ti0.1O2 delivers a high specific capacity of ~130 mAh/g between 2-4 V and 91% retention after 500 cycles. The work discovers multiple materials as high-performance and cost-effective Na-ion battery cathodes and offers critical guidance to the rational design of future layered cathode materials.

Xiao, Biwei↗

DOE OSTI2026

Nuclear Structure and Decay Data for A=169 Isobars

Experimental data pertaining to all nuclei with mass number A=169 (Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Hf, Ta, W, Re, Os, Ir, Pt) have been evaluated. Level schemes from both radioactive decay and reaction studies are presented, along with associated tables of experimental data and adopted properties for levels and γ rays. The present evaluation for A=169 supersedes the 2008 evaluation, 2008Ba31, by C.M. Baglin. A few highlights of this evaluation: More extensive work on ε decay from 169W is needed and new experimental work will be required to resolve a discrepancy between the J π values deduced for a 180-keV level in 169Ta based on extensive band structure from (HI,xnγ) work (J π =1/2−) and TDPAD measurements (J=5/2). Low lying states of 169Os were studied via fine structure of 173Pt α decay in 2014ThZZ. The Eαs feeding the g.s. of 169Os in 2008Ba31 are separated well into two consistent groups to feed the g.s. and the newly proposed state at 34.84 keV. Based on the studies of 2014ThZZ and 2021Zh52, the g.s. spin-parity assignment of 169Os has been proposed to be (7/2−) from (5/2−). The 169Ir g.s. half-life and alpha emission branching reported in 2012Th13 from 173Au α decay measurements are preferred over the values in 2005Sc22. The reported half-life value in 2005Sc22 for 169Ir g.s. is discrepant and the research work was carried out in the same lab of 2012Th13.

Basunia, M Shamsuzzoha↗

DOE OSTI2022

Simulation of natural circulation cartridge loop experiments and application to molten salt reactors

This work uses the TRAC/RELAP Advanced Computational Engine (TRACE) thermal hydraulics (TH) code to model natural circulation cartridge loop experiments previously conducted at Oak Ridge National Laboratory (ORNL) using water and compares the simulated and experimental results. TRACE is also used to characterize natural circulation in the cartridge loop vehicle using FLiNaK as the working fluid. The experimental vehicle is a buoyancy-aided, annular cartridge loop, referred to as a thermosyphon, and is designed to aid in qualifying liquid–fueled and/or liquid–cooled irradiation experiments for the Versatile Test Reactor (VTR), which is currently being designed in the United States. Out-of-pile water experiments have been conducted using the cartridge and the Thermosyphon Test Loop facility at ORNL, and future experiments are anticipated that would use other molten salt surrogates as the working fluid, followed by eventual insertion of a similar cartridge into VTR. Additionally, this work aims to determine how well TRACE can replicate the natural convection conditions that were observed experimentally; this serves as an initial step for validating the modeling tool for design and safety calculations to support future irradiation experiments in VTR. Initial predictions of potential experiments were made using FLiNaK as the natural circulation fluid to demonstrate the relevance of the cartridge design to molten salt reactors (MSRs). Results from this study indicate that TRACE can accurately capture natural convection phenomena in the thermosyphon and that several design changes to the current cartridge vehicle are necessary to achieve hydraulic conditions similar to those expected in MSRs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Sensitivity analysis applied to SiC failure probability in TRISO modeled with BISON

Here, a sensitivity analysis on the failure probability of the Silicon Carbide (SiC) layer in tristructural isotropic (TRISO) nuclear fuel during transient conditions predicted by the BISON fuel performance code is performed. The principal goal of the analysis is to understand the most important parameters dictating SiC failure behavior in BISON during Reactivity Initiated Accidents (RIAs). SiC brittle fracture probability is modeled using Weibull statistics. A total of seven inputs related to SiC failure has been selected for the analysis, including the Weibull statistics parameters, elastic moduli for SiC and Pyrolitic Carbon (PyC) and SiC stress-free temperature. A 1D TRISO BISON model has been established for various reactivity insertions performed at the Nuclear Safety Research Reactor (NSRR). The principal advantage associated with the 1D TRISO model developed in this work is its computational efficiency. The Sobol variance decomposition method is used, and the sensitivity indices are presented for eight different values of the energy deposition. The results show that the two most important parameters impacting the predicted SiC failure probability are the Weibull modulus and the characteristic stress, and a co-variance amongst these parameters is obtained for low reactivity insertions. An additional new finding of this work is that the relative importance of Weibull parameters depends on the energy deposition, and thus reactivity, regime. For low energy depositions, two parameters are of influence on SiC failure probability, while for high energy depositions only one parameter impacts failure probability results. Moreover, optimization of the Weibull modulus and characteristic stress is performed by minimizing the RMSE between BISON failure probability predictions and experimental failure fractions for each energy deposition. This work also demonstrates the validity of the NSRR tests BISON simulations and of the respective sensitivity analysis results as conservative, yet indicative for slower transients characterized by lower deposited energies. Such verification is achieved through the partial extension of the analysis to a group Control Rod Withdrawal reproduced from a previous study. Another novel result of this analysis is that no single set of Weibull parameters can reproduce all reactivity insertion experimental failure results, which is related to the intrinsic nature of SiC failure properties, and a new range for the parameters is proposed to produce a failure probability envelope that encompasses the experimental fractions. Additionaly, this work proposes a new approach for future failure analysis with BISON, consisting in the use of Weibull parameters ranges, rather than fixed sets, along with failure envelopes generation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Insights on keto-hydroperoxide formation from O 2 addition to the beta-tetrahydrofuran radical

Cyclic ethers provide an interesting case study of low-temperature oxidation chemistry, especially in relevance to biofuels. A recent experimental study revealed new questions regarding the low-temperature oxidation mechanism of tetrahydrofuran (THF) concerning the formation of keto-hydroperoxides. In particular, keto-hydroperoxides originating from the THF-β radical were not captured accurately by current literature models, motivating this work to calculate the energetics of the first and second O 2 -addition pathways for THF radicals. Electronic structure calculations at the CCSD(T)/cc-pV∞Z//M06-2X/cc-pVTZ level of theory were used to generate potential energy surfaces for the α-C and β-C THF radicals and subsequent pathways to the formation of the keto-hydroperoxide isomers. These are the first theoretical calculations of the second O 2 -addition radical pathways for the THF-β radical. Results from the theorical work provided further insight into the low-temperature oxidation of THF. This included identifying the pathways most likely to form the keto-hydroperoxide isomers observed in prior experimental work; and detecting that the shortcomings in prior models are likely due to uncertainties in R + THF abstraction reaction rates. Here, these conclusions will motivate future work for accurate THF kinetic model development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗