SEARCH · Engineering Papers
Results for “diffusion behavior”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers
Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.
Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows
Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.
Ultra-selective CANAL Polymers for Hydrogen-based Membrane Separations after Long-Term Aging
Hydrogen is a critically important molecule in the chemicals and energy industries, driving the need for efficient purification technologies such as membrane-based separations. However, polymer membranes often suffer from physical aging, leading to decreased gas permeability over time. This study leverages the unique aging behavior of contorted microporous polymers, synthesized via catalytic arene-norbornene annulation (CANAL) polymerization, to enhance molecular sieving for H 2 separations. While prior studies demonstrated a 1000% increase in H 2 /CH 4 selectivity with minimal permeability loss over 150 days, we extend this investigation to hyperaged (∼1 year) CANAL polymers and evaluate their performance under various conditions of industrial relevance. We report a remarkable 6100% increase in H 2 /CH 4 selectivity and a 2400% improvement in H 2 /N 2 selectivity, with only a 55% decline in H 2 permeability. These membranes exhibit excellent long-term stability in binary and ternary gas mixtures and at elevated temperatures. Additionally, extended aging enhances H 2 /CO 2 selectivity beyond the 2008 upper bound due to improved size-sieving. High-temperature permeation experiments and grand canonical Monte Carlo simulations reveal that H 2 permeability increases from 570 to 3500 barrer as the temperature rises from 35 °C to 190 °C, while H 2 /CO 2 selectivity declines from 14 to 3 due to the higher diffusion activation energy of CO 2 . Despite this tradeoff, hyperaged CANAL membranes maintain performance beyond the 200 °C upper bound, having selectivity–permeability performance comparable to many state-of-the-art membranes for H 2 /CO 2 separations. As a result, this study highlights the robustness and industrial viability of ultra-selective CANAL polymers for hydrogen purification, contributing to advancements in clean energy and sustainable separation technologies.
Hydrodynamic behavior of wrapped coil water heater tank
In this work, a numerical investigation using computational fluid dynamics (CFD) has been carried out to analyze the thermal behavior of a water tank integrated to a heat pump water heating system (HPWH). The heat pump condenser is a long coil wrapped around the water tank and insulated from outside from which heat is transferred from the heat pump cycle into the water stored inside the tank. The simulated water tank has a capacity of about 170 L with the charging inlet located at the bottom segment of the tank and the discharging outlet located at the top segment of the tank like commercially available water heaters. The tank was simulated during charging/discharging process with isothermal boundary condition at the tank wall to represent the majority of the process for refrigerant condensation inside the heat pump condenser (excluding superheated and sub-cooled section). Three flow rates were selected to represent low, medium and high flow pattern. Both mean and turbulent flow properties were calculated and analyzed. Heat transfer mechanisms inside the tank are quantified by calculating Richardson and Péclet numbers to show the relative contribution of driving forces. The calculations show that advection is dominant over diffusion. Mixing factor is calculated and presented. It is found that the rate of increase of mixing is nonlinear with the increase in flow rate. Mixing increased by a factor of four when the flow rate increased by a factor of two.
First-principles study of hydrogen behavior near W/WC interfaces
Hydrogen (H) solution, migration, and binding to vacancies near the tungsten/tungsten carbide (W/WC) interfaces based on hcp WC or bcc W (W/WCb and Wb/WC) have been studied by first-principles simulations to provide insights into the effect of C on H behavior in W exposed to C deposition in fusion reactors. The H behavior near the W/WC interfaces is related to the tensile or compressive strain due to the lattice match and interfacial effects. Near the W/WCb interface, the H solution energy in hcp WC increases approaching the interface, while it decreases from the W side to the interface. The configuration change of the W slab significantly decreases the H solution energy but has only a modest effect on the H migration energy barrier in comparison with the values in bcc W. The interface enhances the H binding to vacancies in WC layers. Both the interface and configuration adjustment of the W slab reduce the H binding to vacancies on the W side of the interface. Near the Wb/WC interface, H solution energies in both WC and W slabs weakly decrease approaching the interface. The configuration change of the WC slab significantly increases the H solution energy and binding energy to vacancies but decreases the H migration energy barrier in comparison with the values in hcp WC. The H energetics near the W/WC interfaces indicate the likelihood of H localization and diffusion in the W side of the interface, rather than diffusing through the interface into WC layers.
Helium effects on the surface and subsurface evolutions in single-crystalline tungsten
Tungsten (W) has been perceived as one of the most promising plasma facing materials (PFMs) for future fusion reactors. In the past decade, its behavior under irradiation and helium (He) plasma interaction has been extensively studied. However, some key knowledge gaps still exist, such as the influence of crystallographic orientation on the surface and subsurface evolutions. In this work, we focus on the He ion-beam irradiation damage effects in mirror-polished single-crystalline W samples with three different surface planes of {100}, {110} and {111}. Irradiation was performed at room temperature using 40 keV He+ to a fluence of 1 × 10 16 /cm 2 , followed by thermal desorption spectroscopy (TDS) up to ~1920 K. The microstructures of He-irradiated W before and after TDS heat treatment were characterized by scanning and transmission electron microscopy. Subsurface He bubbles were imaged in all irradiated samples, but newly formed <111>-oriented surface grains and surface blisters were only observed in W {100} and {110} starting orientations. These results reveal that radiation damage, He thermal desorption, and surface/subsurface evolution are all strongly dependent on crystallographic orientation. Underlying physical mechanisms are discussed based on ion channeling effects, He-vacancy interactions, and surface diffusion. In conclusion, these findings provide new insights into He effects in W.
Effects of unburnt reaction progress on stretch flame dynamics under elevated temperatures
Flames in practical combustors and engines are inevitably subject to the effects of upstream chemical reaction progress and stretch, induced by elevated thermodynamic conditions and flow non-uniformities, respectively. Recent shock tube experiments and simulation studies on flame propagation have shown that flame speed under engine-relevant conditions can be enhanced with non-negligible upstream chemical reaction progress, especially when low-temperature heat release is involved in the unburnt mixture. On the other hand, depending on the mixture equivalence ratio and diluents, nonequidiffusion (including the non-unity Lewis number effect and the preferential diffusion effect) can couple with the flame stretch to fundamentally affect the flame propagation, which can manifest as either facilitation or suppression. Depending on the transport property of a reacting mixture, there hence can be either inhibition or promotion from unburnt reaction progress and stretch effects on flame propagation. In the current work, through one-dimensional numerical simulations of transient planar and spherical flames of n-heptane/air under elevated thermodynamic conditions, the combined effects of upstream chemical reaction progress and stretch on flame propagation are investigated. Results show that for both lean and rich n-heptane/air mixtures, flame speed can be substantially promoted with reaction progress, while the rich mixtures can exhibit opposite stretch dependence after first-stage ignition. Different definitions of effective Lewis number are adopted to explain the change in Markstein length for spherical flame in lean and rich mixtures reformed by the low temperature chemistry. Further, this work fills an important gap in laminar premixed flame research relevant to practical combustion systems and can provide useful insight into local turbulent flame behaviors and phi-sensitivity in engine combustion.
Revealing the Drag Instability in One-fluid Nonideal Magnetohydrodynamic Simulations of a 1D Isothermal C-shock
C-type shocks are believed to be ubiquitous in turbulent molecular clouds thanks to ambipolar diffusion. We investigate whether the drag instability in 1D isothermal C-shocks, inferred from the local linear theory of Gu & Chen, can appear in nonideal magnetohydrodynamic simulations. Two C-shock models (with narrow and broad steady-state shock widths) are considered to represent the typical environment of star-forming clouds. The ionization-recombination equilibrium is adopted for the one-fluid approach. In the 1D simulation, the inflow gas is continuously perturbed by a sinusoidal density fluctuation with a constant frequency. The perturbations clearly grow after entering the C-shock region until they start being damped at the transition to the post-shock region. We show that the profiles of a predominant Fourier mode extracted locally from the simulated growing perturbation match those of the growing mode derived from the linear analysis. Moreover, the local growth rate and wave frequency derived from the predominant mode generally agree with those from the linear theory. Therefore, we confirm the presence of the drag instability in simulated 1D isothermal C-shocks. We also explore the nonlinear behavior of the instability by imposing larger-amplitude perturbations to the simulation. We find that the drag instability is subject to wave steepening, leading to saturated perturbation growth. Issues concerning local analysis, nonlinear effects, one-fluid approach, and astrophysical applications are discussed.
Accelerating network layouts using graph neural networks
Graph layout algorithms used in network visualization represent the first and the most widely used tool to unveil the inner structure and the behavior of complex networks. Current network visualization software relies on the force-directed layout (FDL) algorithm, whose high computational complexity makes the visualization of large real networks computationally prohibitive and traps large graphs into high energy configurations, resulting in hard-to-interpret “hairball” layouts. Here we use Graph Neural Networks (GNN) to accelerate FDL, showing that deep learning can address both limitations of FDL: it offers a 10 to 100 fold improvement in speed while also yielding layouts which are more informative. We analytically derive the speedup offered by GNN, relating it to the number of outliers in the eigenspectrum of the adjacency matrix, predicting that GNNs are particularly effective for networks with communities and local regularities. Finally, we use GNN to generate a three-dimensional layout of the Internet, and introduce additional measures to assess the layout quality and its interpretability, exploring the algorithm’s ability to separate communities and the link-length distribution. The novel use of deep neural networks can help accelerate other network-based optimization problems as well, with applications from reaction-diffusion systems to epidemics.
Effect of Flexibility in Molecular Simulations of Carbon Dioxide Adsorption and Diffusion in a Cuprous Triazolate Framework
Using fixed atom force fields to model gas adsorption in flexible metal organic frameworks (MOFs) is known to pose difficulties in accurately reproducing and predicting experimental results. This paper studies the difference in accuracy between flexible and fixed atom force fields in reproducing CO 2 adsorption measurements in MAF-2 ([Cu(etz)]∞ (MAF-2, Hetz) 3,5-diethyl-1,2,4-triazole), an NbO-type triazolate scaffold with a bcu cavity system and attached ethyl groups. The flexible force field used to run the hybrid molecular dynamics and grand canonical Monte Carlo calculations were generated using the QuickFF software incorporating van der Waals parameters from the Universal Force Field (UFF) and density derived electrostatic and chemical (DDEC) charges. The fixed atom force field used was composed of UFF van der Waals parameters together with DDEC charges. The calculations were run at 298 K and at pressures of 0.1, 0.3, 0.5, 0.8, and 1 bar. It was observed that for this MOF the rigid force field overpredicted gas adsorption, whereas the flexible force field values closely matched experimental data. In the flexible structure, the freely moving ethyl groups of MAF-2 hindered adsorption, reducing the interaction energy between CO 2 and the N atoms of the triazolate framework as well as reducing the size of the largest cavity diameter. Here, the combination of these factors led to improved prediction of adsorption values with the flexible forcefield as compared to the rigid forcefield, demonstrating the need for modeling MOFs in a way more indicative of their behavior.
Particle charge distributions in the effluent of a flow-through atmospheric pressure low temperature plasma
Atmospheric pressure low-temperature plasmas are often utilized to perform particle synthesis, treatment, and removal. It is well-known that dust particles are highly negatively charged in these plasmas; however, little is known about dust particle charging behavior as particles leave the plasma volume and pass through the spatial afterglow region. In this work, monodisperse particles of various sizes and work functions were introduced into an atmospheric pressure radiofrequency capacitively coupled flow-through plasma. Dust particle electrical mobility distributions downstream of the flow-through plasma were measured utilizing a differential mobility analyzer in conjunction with a condensation particle counter at various gas flow velocities. Charge distributions were determined from the measured electrical mobility distributions. Experiments confirm that particles become less negatively charged, and even net-positively charged after leaving the plasma volume, with a distribution that follows a shifted Boltzmann charge distribution. Additionally, particle charge in the effluent of the flow-through plasma is negligibly dependent on work function but highly size and flow velocity dependent. Larger particles were shown to have a higher magnitude of charge under all studied conditions; however, particle polarity was switchable by varying gas flow velocity. The charging dynamics were simulated utilizing a constant number Monte Carlo model that accounts for electron temperature decay and the transition from ambipolar to free diffusion of electrons and ions in the spatial afterglow. Simulation results also suggest that, at the same flow velocity, larger particles obtain a greater magnitude of charge, negative or positive. Furthermore, the decrease in electron mobility and the difference between ion and electron convective loss rates create an ion-rich region in the plasma effluent that promotes ion–particle collisions and drives particle charge removal and even reversal of polarity. Larger particles more favorably collide with energetic species in these environments, which results in higher charge states.
Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)
In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.
High-Temperature 2D Optical Relaxation Visualizes Enhanced Oxygen Exchange Kinetics at Metal-Mixed Conducting Oxide Interfaces
Solid-state heterointerfaces are of interest for emergent local behavior that is distinct from either bulk parent compound. One technologically relevant example is the case of mixed ionic/electronic conductor (MIEC)–metal interfaces, which play an important role in electrochemistry. Metal–MIEC composite electrodes can demonstrate improved catalytic activity vs single-phase MIECs, improving fuel cell efficiency. Similarly, MIEC surface reaction kinetics are often evaluated using techniques that place metal current collectors in contact with the surface under evaluation, potentially altering the response vs the native surface. Here, techniques enabling direct and local in situ observation of the behavior at and around such heterointerfaces are needed. Here, we develop a spatially resolved optical transmission relaxation (2D-OTR) method providing continuous evaluation of local, high-temperature, controlled atmosphere defect kinetics across a ~1 cm2 sample area simultaneously in a contact-free manner. We apply it to observe the spatial variance of oxygen incorporation and evolution rates at ~525–620 °C, in response to step changes in oxygen partial pressure, on MIEC SrTi 0.65 Fe 0.35 O 3–x films as a function of distance from porous Pt and Au layers. Using this model geometry, we find significant enhancements in kinetics adjacent to the metals that decay over a few millimeter distance. To extract kinetic parameters, we fit the short-term optical data (initial portion of relaxations) with an exponential decay function appropriate for surface-exchange-limited kinetics, yielding apparent surface exchange coefficients (k chem ) with spatial resolution, decreasing with distance from the metal. To understand the kinetic processes governing the complete (long-term) optical relaxations, we performed COMSOL simulations, which demonstrated that a combination of laterally varying k chem and in-plane diffusion controls the observed kinetics over the full time range. Further support for spatially varying kchem comes from demonstrations of changing surface and bulk chemistry vs distance from the metal–MIEC interface, by X-ray photoelectron and optical absorption spectroscopies, respectively. Although microporous Pt and Au are not excellent electrodes in isolation, both metals exert a synergistic effect on the oxygen surface exchange rate in the presence of the mixed conducting film.
Finite temperature properties of uranium mononitride
Uranium mononitride (UN) is a promising nuclear fuel that combines the advantageous properties of readily used UO 2 and uranium alloys, such as high melting temperature and high uranium density, and thermal conductivity, respectively. A better understanding of UN behavior at operating temperatures can be obtained from finite temperature data, such as elastic properties. To get this information, ab initio molecular dynamics (AIMD) simulations were performed at five different temperatures using constant volume (NVT) and constant pressure (NPT) ensembles. Initially, the performance of PBE functional in reproducing experimental crystallographic properties and magnetic ordering is assessed. The finite temperature phonon dispersions are calculated using NVT simulation results, which show a softening of the phonon modes with increasing temperature. The NPT results are used to obtain the thermal expansion of UN and finite temperature electronic properties. The calculated thermal expansion is compared with our measurements using neutron diffraction. Additionally, the temperature dependent elastic properties of UN are evaluated using the strain-stress method in AIMD simulations, indicating that UN becomes softer and more compressible with increasing temperature. Also, the calculated Young’s modulus slope is in very good agreement with the experiment. The finite temperature heat capacity and electronic thermal conductivity are calculated from AIMD simulations, which are in better agreement with the experiment than the heat capacity and thermal conductivity calculated using the structures relaxed at 0K. Finally, the thermal diffusivity from AIMD has opposite temperature dependence compared to experimental results, which we argued comes from the underestimated electronic thermal conductivity.
RELAP5-3D Modeling of High Temperature Test Facility (HTTF) Test PG-26
The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled nuclear reactors with prismatic fuel and reflector blocks. Several tests have been completed, and more are still planned to at the HTTF, including depressurized conduction cooldown (DCC) and pressurized conduction cooldown (PCC) transients. This report analyses test PG-26, a progression of the Double Ended Inlet-Outlet Crossover Duct Break transient that is referred to as a DCC. PG-26 has been performed at the HTTF between May 30 and June 30, 2019. Core initial conditions (i.e., before the DCC started) have been met using low power (<100 kW) and two of ten available electric heaters. The DCC transient was initiated during the 50 th hour of the test. The break valves were opened, and hot helium from the core and cold helium from the reactor cavity simulation tank (RCST) started mixing. The gases flowed in a countercurrent fashion, where the top half of the hot duct contained hot helium that flowed in one direction and cold helium that flowed in the other direction in the bottom half of the duct. After the pressure and density reached equilibrium, the event entered a diffusion mode. The onset of a reverse natural circulation was not observed during the DCC period of the test. Version 4.4.2ie of the RELAP5-3D computer code has been used to model the HTTF PG-26 test, and results have been compared to available high-quality measured data. The model used in this study is the quality-controlled HTTF RELAP5-3D model (HTTF base 2018-04-19 QA), originally developed by P. Bayless. The report includes RELAP5-3D results of the “base calculations” as well as some sensitivities to important uncertain model inputs, such as primary helium mass flow rate, core ceramic thermal properties, as well as heat evacuation and loop friction models. Using the base RELAP5-3D model predicts a countercurrent helium flow in the hot duct observed at the beginning of the DCC, but instead of going into a molecular diffusion mode, the model predicts the onset of natural convection. Increasing friction in the core and hot duct prevents the natural convection from happening in some of the simulations. Although some temperatures are well predicted (and even overpredicted), the general tendency is to underpredict the ceramic and helium temperatures and heat removal rates during the DCC, resulting in many of the assessment findings being in minimal or insufficient agreement with the data. It is worth noting that the described discrepancies between measured data and RELAP5-3D predictions are not RELAP5-3D code limitations. More so, they reflect limitations in boundary condition and thermal property knowledge. While the RELAP5-3D calculations of the test provide some insights into what happens during the transient, and point to missing or potentially uncertain data to which the experimenters can direct their attention, the principal conclusion is that the PG-26 test data are insufficient for a system code assessment.
The Application of Differential Scanning Calorimetry to Investigate Precipitation Behavior in Nickel-Base Superalloys Under Continuous Cooling and Heating Conditions
In this study, a suite of experimental tools and fast-acting, numerical-simulation techniques was used to quantify the precipitation behavior of three nickel-base superalloys: IN-100, LSHR, and 718. Experimental methods comprised differential scanning calorimetry (DSC) to establish the specific heat as a function of temperature and selected direct-resistance heating trials (using a Gleeble ® machine) to obtain samples for microstructural analysis. For the DSC experiments, each alloy was cooled at a prescribed constant rate (between 5 and 20 K/min) after an initial soak/equilibration in the high-temperature, single-phase (supersolvus) temperature regime. On-heating DSC trials beginning at ambient temperature were also performed on alloy 718 in three different starting conditions: super- δ -solvus solution treated and water quenched (denoted as ST), solution treated and aged (STA), and solution treated and overaged (STOA). DSC results, revealing the thermal signatures associated with the kinetics of precipitation of γ ' (IN-100, LSHR) or γ ' and γ " (718), were interpreted using a previously-developed fast-acting routine that treats concurrent nucleation, growth, coarsening, and dissolution. For these simulations, special attention was paid to various thermo-kinetic input parameters including equilibrium solvus-approach curves, bulk free energies of transformation, matrix-precipitate interface energies, and effective diffusivities. For the γ - γ ' superalloys (IN-100 and LSHR), estimates of precipitate volume fraction as a function of temperature from the specific-heat data revealed semi-quantitative agreement with simulation predictions. For the γ - γ '- γ " superalloy (718), simulation predictions of precipitate volume fractions were converted to specific heat as a function of temperature and showed semi-quantitative agreement with the direct measurements.
Studies of Helium Behavior in Uranium Hexafluoride
Several issues have arisen from experimental observations and/or theoretical postulation related to the behavior of helium in uranium hexafluoride (UF 6 ) matrices and how it may impact the ability to use this noble gas to quantify the fill date of the cylinder. A series of experiments was conducted to address several such issues, and the results and implications are reported here. Two issues relate to bias in the helium measurement due to ambient air from extraneous sources such as atmospheric air or from potential holdup in a large, activated alumina chemical trap. These issues are really a question of differential retention of helium versus three other inert gas isotopes (representative of other noble gases found in air) used as a calibration for air in leakage. Discernable bias was not detected in either case. Another issue had to do with helium potentially brought into the cylinder by a liquid UF 6 fill and subsequently released as the material solidifies and/or during storage. The solubility of helium in liquid UF 6 was measured to investigate this theory. The experiment quantified the potential magnitude of the initial helium contribution if present. The final issue was that of retention of helium in solid UF 6 . Investigation of this theory necessitated quantifying both the solubility and the diffusivity of helium in UF 6 . In one experiment, the solubility was measured at cryogenic temperature, allowing quantification of helium in solid UF 6 such as would be found in a cold trap. In a second experiment, transport of helium into and through a macroscopic plug of solid UF 6 was examined. The transport was very slight although apparently detectible. Effectively, only an upper limit to the transport could be determined, but that limit implied that helium holdup in solid UF 6 is a real possibility.