SEARCH · Engineering Papers

Results for “Conditional rate functions”

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 1,549 records · Page 86

All sources NASA NTRS DOE OSTI

NASA NTRS2024

Aerothermal Analysis of the Rocket Lab Venus Probe Heatshield

This document provides an overview of the aerothermodynamic analyses performed by the Aerothermodynamics Branch at NASA Langley Research Center for the Rocket Lab Venus Probe (RLVP). In addition to defining the baseline heating environment to the heatshield, this document pursues the experimental validation of key physical models at RLVP-relevant conditions. This experimental validation analysis, which captures the model form uncertainty, is used as one of two primary components of the margin assessment, where the other component is the parametric uncertainty. These model form (experimental) and parametric uncertainty components are used to construct a spatial and time varying margin for the heating to the RLVP heatshield. The margin is evaluated as the sum of the parametric and model form uncertainty components. The model form uncertainty is defined as the difference between the RLVP-relevant measurements and their simulations, using the upper limit uncertainty bounds for both the measurements and simulations in the comparisons. The differences in the dominant physics in the stagnation region and flank lead to the separate RLVP-relevant measurements for assessing the model form uncertainty in these two regions. These regions are addressed as follows: Stagnation Region Heating Environment: For the high-temperature stagnation-region, both the radiative heating and impact of blowing on convective heating are significant, while the impacts of turbulence and roughness are negligible. Coupled radiation and ablation LAURA/HARA solutions with ray-tracing provide the radiative heating over the entire vehicle, including the contributions from the Venus atmosphere and ablation species. Non-ablating LAURA simulations provide the convective heating. During the material-response computation typically used for TPS sizing, this non-ablating convective heating is corrected for the impact of ablation using the blowing correction. Coupled ablation LAURA simulations that capture finite-rate sur-face processes show that this blowing correction may be non-conservative over most of the heatshield. This non-conservatism is due to hydrogen recombination in the finite-rate surface model, which tends to increase the coupled ablation convective heating to near the non-ablating values, therefore making any reduction in the non-ablating value through the blowing correction non-conservative. This non-conservatism due to H catalysis is captured in the parametric component of the margin. The best available ground-test measurements that capture the impact of blowing on stagnation region convective heating, at RLVP-relevant conditions, indicate that the current blowing reduction model is non-conservative by up to 20% at RLVP-relevant blowing rates (the coupled blowing simulations were also non-conservative). Because of the relatively low velocity of the ground tests and the non-Venus atmospheric chemistry, these measurements do not capture the chemistry and therefore do not inform the uncertainty due to H catalysis. However, they do capture the fluid mechanics of blowing. The non-conservatism of the blowing correction implied by these measurements is covered by the model form component of the margin, which leads to total margin values over 50%. For the radiative heating, the shock-tube informed bias approach suggests a model form uncertainty of roughly 20%, while the parametric uncertainty analysis suggests values over 100%. The combined stagnation-point radiation margin of over 100% leads to peak margined radiative heating values of over300 W/cm2, which remains small relative to the peak margined convective heating of nearly 2000 W/cm2. Based on this analysis, at the stagnation point, the peak margined heat rate is 2203 W/cm2 and the margined total heat load is 31.5 kJ/cm2 for the current nominal trajectory. Flank Heating Environment: The forebody flank (and near-shoulder) heating environment is dominated by the impact of turbulence, roughness augmentation, and ablation on the convective heating. An extensive collection of ground test measurements with RLVP-relevant turbulence and roughness is studied to show that the maximum difference between the simulated and measured convective heating is 5%. However, with the exception of the Holden measurements from the 1980s, these measurements do not include roughness elements extending into the supersonic region of the boundary layer, which is likely to occur for RLVP (due to the 45 degree sphere-cone geometry). The interaction between the supersonic flow and roughness could cause convective heating augmentation fundamentally different than for locally subsonic flow. Although these Holden measurements are consistent with the other measurements considered, another path was pursued to assure that the rough-ness height extending into supersonic flow does not fundamentally change the roughness augmentation. This additional path was a computational effort to resolve the roughness elements in the CFD grid, so that the interaction be-tween the roughness elements and locally supersonic flow may be simulated in detail. This roughness-resolved CFD simulation is feasible because of the RLVP forebody TPS’s patterned roughness, which may be approximated analytically, and because of the axisymmetric nominal flow field, which allows a narrow surface region to be simulated and therefore make the computational expense feasible. These grid-resolved roughness simulations, which are performed at actual RLVP flight conditions, result in heating augmentation values that are below the design approach for roughness augmentation. This provides evidence that the design approach for RLVP roughness augmentation is sufficient. Based on this analysis, at this flank or near-shoulder location, the peak margined total heat rate is 2088 W/cm2and the margined total heat load is 26.0 kJ/cm2for the current nominal trajectory. Heat flux, shear, pressure and heat transfer coefficient at the RLVP stagnation point and near shoulder location are evaluated for the entire trajectory, and curved fit to a functional form of F=AρB∞UC∞. These simplified relationships for the nominal and margined aerothermal environments are referred to as aerothermal indicators, and presented at the end of this document.

Christopher O Johnston↗

DOE OSTI2021

Simulation of PV Variability as a Function of PV Generation and Plant Size: Preprint

The deployment of photovoltaic (PV) systems continues to show significant expansion; however, this growth has brought added attention to issues around the variability of the solar resource. Both spatial and temporal variability exist. Temporal scales can range from the sub-second to multiyear, whereas spatial scales can range from a few meters to tens of kilometers. There are multiple methods described in the literature to quantify PV variability at various spatial and temporal scales. This study focuses on short-term temporal variability and uses similar approaches with the addition of PV plant size a parameter to quantify variability. The method employed here incorporates the normalization of clear- and cloudy-sky conditions and PV plant size to quantify nominal variability metrics. The distribution and fluctuations of these metrics provide relevant information that is useful for system operations. The National Solar Radiation Database (NSRDB) is used to simulate PV variability as a function of PV generation and plant size. Hypothetical but realistic system information at 33 locations is used to model PV generation by feeding NSRDB solar irradiance data to the National Renewable Energy Laboratory’s System Advisor Model (SAM). Over the selected region, it is found that the aggregated ramp rates for the 1-minute data are associated with standard deviations ranging from 0.002–0.055 on a daily basis; however, hourly intervals induce higher aggregated ramp rates than the other timescales. Even though minute-to-minute variations are significant for the 1-minute timescale, the standard deviation aggregated into a daily metric is smaller because of the cancellation of values.

41 EE - Solar Energy Technologies Office (EE-4S)↗

DOE OSTI2015

A priori and a posteriori investigations for developing large eddy simulations of multi-species turbulent mixing under high-pressure conditions

A Direct Numerical Simulation (DNS) database was created representing mixing of species under high-pressure conditions. The configuration considered is that of a temporally evolving mixing layer. The database was examined and analyzed for the purpose of modeling some of the unclosed terms that appear in the Large Eddy Simulation (LES) equations. Several metrics are used to understand the LES modeling requirements. First, a statistical analysis of the DNS-database large-scale flow structures was performed to provide a metric for probing the accuracy of the proposed LES models as the flow fields obtained from accurate LESs should contain structures of morphology statistically similar to those observed in the filtered-and-coarsened DNS (FC-DNS) fields. To characterize the morphology of the large-scales structures, the Minkowski functionals of the iso-surfaces were evaluated for two different fields: the second-invariant of the rate of deformation tensor and the irreversible entropy production rate. To remove the presence of the small flow scales, both of these fields were computed using the FC-DNS solutions. It was found that the large-scale structures of the irreversible entropy production rate exhibit higher morphological complexity than those of the second invariant of the rate of deformation tensor, indicating that the burden of modeling will be on recovering the thermodynamic fields. Second, to evaluate the physical effects which must be modeled at the subfilter scale, an a priori analysis was conducted. This a priori analysis, conducted in the coarse-grid LES regime, revealed that standard closures for the filtered pressure, the filtered heat flux, and the filtered species mass fluxes, in which a filtered function of a variable is equal to the function of the filtered variable, may no longer be valid for the high-pressure flows considered in this study. The terms requiring modeling are the filtered pressure, the filtered heat flux, the filtered pressure work, and the filtered species mass fluxes. Improved models were developed based on a scale-similarity approach and were found to perform considerably better than the classical ones. These improved models were also assessed in an a posteriori study. Different combinations of the standard models and the improved ones were tested. At the relatively small Reynolds numbers achievable in DNS and at the relatively small filter widths used here, the standard models for the filtered pressure, the filtered heat flux, and the filtered species fluxes were found to yield accurate results for the morphology of the large-scale structures present in the flow. Analysis of the temporal evolution of several volume-averaged quantities representative of the mixing layer growth, and of the cross-stream variation of homogeneous-plane averages and second-order correlations, as well as of visualizations, indicated that the models performed equivalently for the conditions of the simulations. The expectation is that at the much larger Reynolds numbers and much larger filter widths used in practical applications, the improved models will have much more accurate performance than the standard one.

42 ENGINEERING↗

DOE OSTI2022

Estimation of the time for steam generator trip due to cyber intrusions

The time required to trip a pressurized water reactor (PWR) by inserting malicious signals into its steam generator (SG) control system has been studied using the Generic PWR (GPWR) Simulator. A semi-analytical model is developed to approximately reproduce the simulator response and understand the dynamics of the control unit. A series of two proportional-integral controllers determines control action according to preset constants, the readings from the feedwater level sensor, and those from feedwater and steam flowrate transmitters. It is observed that the most important factor that determines whether a trip will occur is how much additional water is added to or withheld from the SG over time compared to normal operating conditions. In order to determine the effects of control action on the SG, changes in mass inventory are considered. This approach models the SG water level as a function of mass inventory and has a backward temporal memory. A Python interface is developed for the GPWR framework to automatically simulate different spoofing scenarios and post-process the related data. We observe that the trip times predominantly depend on flow mismatch and/or level errors. Controller parameters, including the integral time and gain constants, either speed up or slow down the rate of progression to a trip setpoint but do not cause a trip by themselves. The reactor can trip on a high-level signal when the reading crosses above 78%, increased from its reference level of 57%, or a low-level reading when it is below 25%. The present results show roughly how long the operators would have to respond to an attack, given a specific set of spoofing signals within the issue space analyzed. Furthermore, we have generated a simple surface by fitting a combination of exponential functions to the data obtained from the GPWR Simulator. In general, trips on a low level have been observed to occur faster than those on a high level.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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 OSTI2016

The mechanism and kinetics of propene ammoxidation over α-bismuth molybdate

Propene ammoxidation over Bi 2 Mo 3 O 12 was investigated to elucidate product (acrylonitrile, acetonitrile, HCN, acrolein, N 2 , etc.) formation pathways. Propene consumption rate is first order in propene and zero order in ammonia (for NH 3 /C 3 H 6 = 0-2) and oxygen (for O 2 /C 3 H 6 ≥ 1.5) partial pressures, with an activation energy (E a = 22 kcal/mol) comparable to that for propene oxidation, suggesting the same rate-limiting step for both reactions. We propose two N-containing species are relevant at ammoxidation conditions: adsorbed NH 3 on surface Bi 3+ ions that reacts with a propene derivative to form products with C-N bonds, and a few metastable M-NH x (M = Mo, Bi; x = 1, 2) groups that are very sensitive to destruction by water, but that are responsible for NH 3 oxidation to N 2 . A proposed reaction mechanism and model that captures the experimental trends in product distribution as a function of partial pressures and temperature are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Impact of hydrate saturation on water permeability in hydrate-bearing sediments

Permeability of the hydrate-bearing sediments vitally affects the hydrate dissociation process as well as the rate and efficiency of gas production. Reported permeability values are observed to be widely are very scattered owing to the dependence on multiple factors such as experimental conditions and test procedures. It is critical that the permeability is measured accurately to enable prediction of long-term gas production using numerical simulation for an economic development of hydrate-bearing reservoirs. In this study, the tetrahydrofuran (THF) hydrates that exhibit pore habits of the pore-filling pattern are formed in sediments, and the water permeability is quantified as a function of hydrate saturation. The results reveal that a decrease in the permeability for an increase in the hydrate saturation. Based on the measured values that are bounded by Kozeny grain-coating and pore-filling models, the fitting parameters of the empirical permeability models are indicated for use in numerical simulation of long-term hydrate dissociation. Wave velocity measurements reveal pore habits of pore-filling pattern of THF hydrates in the sediments with S h < 0.5, and a gradual transition to patchy and load-bearing pore-habits when S h > 0.6. The numerical simulation results obtained using the complementary pore-network model suggest that the THF hydrates in the sediments may form in small clusters with an average patch size of ~4 pores.

58 GEOSCIENCES↗

DOE OSTI2020

Chespa: Streamlining Expansive Chemical Space Evaluation of Molecular Sets

Thousands of chemical properties can be calculated for small molecules, which can be used to place the molecules within the context of a broader “chemical space.” These definitions vary based on compounds of interest and the goals for the given chemical space definition. Here, we introduce a customizable (i.e., modular) Python module, chespa, built to easily assess different chemical space definitions through cluster-ing of compounds in these spaces and visualize trends of these clusters. To demonstrate this, chespa currently streamlines prediction of vari-ous molecule descriptors (predicted chemical properties, molecular substructures, AI-based chemical space, and chemical class ontology) in order to test 6 different chemical space definitions. Furthermore, we investigated how these varying definitions trend with mass spectrometry (MS)-based observability, i.e., the ability of a molecule to be observed with MS (e.g., as a function of the molecule ionizability), using an example data set from the U.S. EPA's Non-Targeted Analysis Collaborative Trial (ENTACT), where blinded samples had been analyzed previously, providing 1,398 data points. Improved understanding of observability would offer many advantages in small molecule identifica-tion, such as (i) a priori selection of experimental conditions based on suspected sample composition, (ii) the ability to reduce the number of candidate structures during compound identification by removing those less likely to ionize, and, in turn, (iii) a reduced false discovery rate and increased confidence in identifications. Factors controlling observability are not fully understood, making prediction of this property non-trivial and a prime candidate for chemical space analysis. Chespa is available at github.com/pnnl/chespa.

Nunez, Jamie↗

DOE OSTI2015

THE LAUNCHING OF COLD CLOUDS BY GALAXY OUTFLOWS. I. HYDRODYNAMIC INTERACTIONS WITH RADIATIVE COOLING

To better understand the nature of the multiphase material found in outflowing galaxies, we study the evolution of cold clouds embedded in flows of hot and fast material. Using a suite of adaptive mesh refinement simulations that include radiative cooling, we investigate both cloud mass loss and cloud acceleration under the full range of conditions observed in galaxy outflows. The simulations are designed to track the cloud center of mass, enabling us to study the cloud evolution at long disruption times. For supersonic flows, a Mach cone forms around the cloud, which damps the Kelvin–Helmholtz instability but also establishes a streamwise pressure gradient that stretches the cloud apart. If time is expressed in units of the cloud crushing time, both the cloud lifetime and the cloud acceleration rate are independent of cloud radius, and we find simple scalings for these quantities as a function of the Mach number of the external medium. A resolution study suggests that our simulations accurately describe the evolution of cold clouds in the absence of thermal conduction and magnetic fields, physical processes whose roles will be studied in forthcoming papers.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Neuromorphic ionic computing in droplet interface synapses

Ionic devices with memory capabilities can emulate neural functionality, enabling neuromorphic computing and biomedical applications. In this study, we report an ionic spiking synapse based on aqueous droplet interface bilayer assembly. Under stepwise triangular voltages, the device displays coupled memcapacitive-memristive behavior, showing noncrossing pinched hysteretic I-V loops. This hysteretic ion dynamics can be regulated by modifying bilayer components, reconstituting protein channels, or adjusting droplet assembly configuration. Droplet interface synapses (DIS) exhibit fundamental neuromorphic behaviors such as paired-pulse facilitation/depression, spike rate–dependent plasticity, Hebbian learning, and short-term associative learning under classical conditioning. We also used reservoir computing with DIS to implement two learning algorithms: a classification algorithm that recognizes handwritten digits and a reinforcement learning algorithm that learns to play a board game of tic-tac-toe.

Li, Zhongwu [Lawrence Livermore National Laborator↗

DOE OSTI2015

Northeastern Center for Chemical Energy Storage (NECCES)

The chemical reactions that occur in batteries are complex, spanning a wide range of time and length scales from atomic jumps to the entire battery structure. The NECCES team of experimentalists and theorists made use of, and developed new methodologies to determine how model compound electrodes function in real time, as batteries are cycled. The team determined that kinetic control of intercalation reactions (reactions in which the crystalline structure is maintained) can be achieved by control of the materials morphology and explains and allows for the high rates of many intercalation reactions where the fundamental properties might indicate poor behavior in a battery application. The small overvoltage required for kinetic control is technically effective and economically feasible. A wide range of state-of-the-art operando techniques was developed to study materials under realistic battery conditions, which are now available to the scientific community. The team also investigated the key reaction steps in conversion electrodes, where the crystal structure is destroyed on reaction with lithium and rebuilt on lithium removal. These so-called conversion reactions have in principle much higher capacities, but were found to form very reactive discharge products that reduce the overall energy efficiency on cycling. It was found that by mixing either the anion, as in FeOF, or the cation, as in Cu1-yFeyF2, the capacity on cycling could be improved. The fundamental understanding of the reactions occurring in electrode materials gained in this study will allow for the development of much improved battery systems for energy storage. This will benefit the public in longer lived electronics, higher electric vehicle ranges at lower costs, and improved grid storage that also enables renewable energy supplies such as wind and solar.

25 ENERGY STORAGE↗

NASA NTRS1980

Pollutant emissions from flat-flame burners at high pressures

Maximum flame temperatures and pollutant emission measurements for NOx, CO, and UHC (unburned hydrocarbons) are reported for premixed methane air flat flames at constant total mass flow rate over the pressure range from 1.9 to 30 atm and for equivalence ratios from 0.84 to 1.12. For any given pressure, maxima typically occur in both the temperature and NOx emissions curves slightly to the lean side of stoichiometric conditions. The UHC emissions show minima at roughly the same equivalence ratios. The CO emissions, however, increase continually with increasing equivalence ratio. Flame temperature and NOx emissions decrease with increasing pressure, while the opposite is true for the CO and UHC emissions. The NOx data correlate reasonably well as a function of flame temperature only. Four flameholders, differing only slightly, were used. In general, the temperature and emissions data from these four flameholders are similar, but some differences also exist. These differences appear to be related to minor variations in the condition of the flameholder surfaces.

Maahs, H. G.↗

NASA NTRS1998

Small Crack Growth and Fatigue Life Predictions for High-Strength Aluminium Alloys: Experimental and Fracture Mechanics Analysis - Part 1

The small crack effect was investigated in two high-strength aluminium alloys: 7075-T6 bare and LC9cs clad alloy. Both experimental and analytical investigations were conducted to study crack initiation and growth of small cracks. In the experimental program, fatigue tests, small crack and large crack tests A,ere conducted under constant amplitude and Mini-TWIST spectrum loading conditions. A pronounced small crack effect was observed in both materials, especially for the negative stress ratios. For all loading conditions, most of the fatigue life of the SENT specimens was shown to be crack propagation from initial material defects or from the cladding layer. In the analysis program, three-dimensional finite element and A weight function methods were used to determine stress intensity factors and to develop SIF equations for surface and corner cracks at the notch in the SENT specimens. A plastisity-induced crack-closure model was used to correlate small and large crack data, and to make fatigue life predictions, Predicted crack-growth rates and fatigue lives agreed well with experiments. A total fatigue life prediction method for the aluminum alloys was developed and demonstrated using the crack-closure model.

Wu, X. R.↗

NASA NTRS1996

Restoration of plasma volume after 16 days of head-down tilt induced by a single bout of maximal exercise

Seven healthy men performed maximal exercise 24 h before the end of 16 days exposure to 6 degrees head-down tilt (HDT) to test the hypothesis that such an exercise technique could restore plasma volume (PV) at the end of a simulated space mission. Exercise consisted of supine cycling with graded work rates increasing by 16 W/min to volitional fatigue and required an average of 16 min. The experimental protocol was a standard cross-over design in which the order of treatment (exercise or control) was counterbalanced across all seven subjects. PV, fluid intake (ad libitum), urine output, renal function, and hormones associated with fluid homeostasis were measured before HDT, 24 h before the end of HDT just prior to exercise, and at the end of HDT 24 h after exercise. HDT reduced PV by 16% in both control and exercise conditions. Maximal exercise completely restored plasma volume within 24 h to 3.9 +/- 3.2% of pre-HDT levels despite continued HDT. Compared with control, exercise induced a 660-ml larger positive fluid balance because of greater fluid intake and reduced urine volume during the 24 h after exercise. These results suggest that one bout of maximal leg exercise before return from 16 days of spaceflight may be completely effective in stimulating thirst and restoring plasma volume to preflight levels.

Non-NASA Center↗

NASA NTRS2014

Multiphysics Modeling of an Annular Linear Induction Pump With Applications to Space Nuclear Power Systems

An annular linear induction pump (ALIP) that could be used for circulating liquid-metal coolant in a fission surface power reactor system is modeled in the present work using the computational COMSOL Multiphysics package. The pump is modeled using a two-dimensional, axisymmetric geometry and solved under conditions similar to those used during experimental pump testing. Real, nonlinear, temperature-dependent material properties can be incorporated into the model for both the electrically-conducting working fluid in the pump (NaK-78) and structural components of the pump. The intricate three-phase coil configuration of the pump is implemented in the model to produce an axially-traveling magnetic wave that is qualitatively similar to the measured magnetic wave. The model qualitatively captures the expected feature of a peak in efficiency as a function of flow rate.

Kilbane, J.↗

NASA NTRS2016

Thermal Management and Thermal Protection Systems

During my internship in the Thermal Design Branch (ES3), I contributed to two main projects: i) novel passive thermal management system for future human exploration, ii) AVCOAT undercut thermal analysis. i) As NASA prepares to further expand human and robotic presence in space, it is well known that spacecraft architectures will be challenged with unprecedented thermal environments. Future exploration activities will have the need of thermal management systems that can provide higher reliability, mass and power reduction and increased performance. In an effort to start addressing the current technical gaps the NASA Johnson Space Center Passive Thermal Discipline has engaged in technology development activities. One of these activities was done through an in-house Passive Thermal Management System (PTMS) design for a lunar lander. The proposed PTMS, functional in both microgravity and gravity environments, consists of three main components: a heat spreader, a novel hybrid wick Variable Conductance Heat Pipe (VCHP), and a radiator. The aim of this PTMS is to keep electronics on a vehicle within their temperature limits (0 and 50 C for the current design) during all mission phases including multiple lunar day/night cycles. The VCHP was tested to verify its thermal performance. I created a thermal math model using Thermal Desktop (TD) and analyzed it to predict the PTMS performance. After testing, the test data provided a means to correlate the thermal math model. This correlation took into account conduction and convection heat transfer, representing the actual benchtop test. Since this PTMS is proposed for space missions, a vacuum test will be taking place to provide confidence that the system is functional in space environments. Therefore, the model was modified to include a vacuum chamber with a liquid nitrogen shroud while taking into account conduction and radiation heat transfer. Infrared Lamps were modelled and introduced into the model to simulate the sun's rays directly impinging on the system. Heating rate of the lamps were calculated by knowing fraction of emitted energy in a wavelength interval and the filament temperature. This version of the model can be used to predict performance of the system under vacuum with extreme cold or hot conditions. Initial testing of the PTMS showed promise, and the thermal math model predicts even better performance in thermal vacuum testing. ii) Thermal Protection Systems (TPS) are required for vehicles which enter earth's atmosphere to protect from aerodynamic heating caused by the friction between the vehicle and atmospheric gases. Orion's heat shield design has two aspects which needed to be analyzed thermally: i) a small excess of adhesive used to bond the outer AVCOAT layer to the inner composite structure tends to seep from under the AVCOAT and form a small bead in between two bricks of AVCOAT, ii) a silicone rubber with different thermophysical properties than AVCOAT fills the gap between two bricks of AVCOAT. I created a thermal model using TD to determine temperature differences that are caused by these two features. To prevent false results, all TD models must be verified against something known. In this case, the TD model was correlated to CHAR, an ablation modelling software used to analyze TPS. Analyzing a node far from the concerning features, we saw that the TD model data match CHAR data, verifying the TD model. Next, the temperature of the silicone rubber as well as the bead of adhesive were analyzed to determine if they exceeded allowable temperatures. It was determined that these two features do not have a significant effect on the max temperature of the heat shield. This model can be modified to check temperatures at various locations of the heat shield where the composite thickness varies.

Hasnain, Aqib↗

DOE OSTI2021

Quantifying Properties for a Mechanistic, Predictive Understanding of Aqueous Impact on Aging of Medium and Low Voltage AC and DC Cabling in Nuclear Power Plants

To address the gap in knowledge in understanding degradation in relevant service conditions, this project aimed to develop a mechanistic, predictive model of medium and low voltage cable failure based on the primary environmental degradation parameters of aqueous immersion time, temperature, and the oxidation extent. To do so, we took a two-fold approach toward evaluating degradation as related to the aqueous condition. First, we evaluated the chemical, mechanical, and electrical properties of polymers that comprise the cabling insulation under varied aqueous conditions at different temperatures. Secondly, we evaluated the performance of insulation materials under similar accelerated aging conditions in addition to the dielectric breakdown of the cabling under submersion conditions. Thermal oxidation and immersion of LDPE thin films were done using a Parr vessel and the extent of oxidation was monitored using the carbonyl index from ATR-FTIR spectrum, a ratio quantifying the carbonyl functionality produced on the LDPE film surface. Increasing the temperature, oxidation time, and oxygen pressure increased the water-vapor permeability and the carbonyl content. Thermogravimetric analysis of the films confirmed that there is indeed an initial weight loss of the smaller molecular weight molecules and volatiles. At higher temperatures, there is a secondary mass loss. Up to 70°C, the mass loss curves were similar until 80°C, when the aged films decreased at a much greater rate with respect to increasing temperature. This was seen in the 10, 50, 90% weight loss temperatures. The mechanical properties of sheet PP and HDPE as well as injection molded LDPE, HDPE, HDPP dog bones were plotted with respect to the time spent in the Parr vessel during thermal oxidation. For the sheet PP and HDPE, the UTS and the modulus of elasticity decreased while the elongation increased. The injection-molded LDPE, HDPE, and HDPP dog bones experienced a similar trend where the changes in the properties were within experimental uncertainty. This analysis of LDPE in dry and immersive oxygen-rich environments is an important baseline for future experiments looking at other degradation mechanisms. A predictive aging model for low-density polyethylene insulative cable housings was developed. This model revealed a combination of physical and chemical processes which lead to an increase in permeability and a decrease in strength of the insulator. In creating this predictive model, a gap of knowledge in the literature was revealed in two areas: understanding how the crystallinity of a polymer changes with time and visualizing the predicted pores formed through the insulator which leads to failure. Developing and using an adapted ATR-FTIR method, crystallinity was monitored and compared to bulk crystallinity found via DSC. Inhomogeneity in crystallinity changes were observed but the FTIR method but was found to be less reliable. Pore visualization was found utilizing electrical impedance spectroscopy saturated aged polyethylene films with synthesized citrate-capped gold nanoparticles, and chloroauric acid precursor. results indicate a decreasing impedance of the polyethylene films caused by the transport of ions through the film. SEM imaging was then performed for elemental analysis of the films and counter electrodes for the presence of ions. Chloride and gold ions were detected; however, no nanoparticles were found. This gives an estimated pore size of at least 0.3 nm through the aged film. Cyclic submergence of polyethylene and polypropylene in aqueous solutions of copper sulfate and Harrison’s solution were examined. It was discovered that aging in these mixed conditions and at 90°C did cause a small, yet significant increase to tensile strength for both unaged polyethylene and polypropylene. The PE and PP in cycled and submerged conditions did not have tensile strengths significantly different from dry-aged after 16 weeks. For both solutions, it was determined that capacitance increases with both water tree depth and cable temperature and is relatively independent of changing water tree AR. As for resistance, there is no apparent change between depth percentages of 10 and 80 for both solutions. This is because the water tree has not yet entered the conductor and thus there is no shorting current flowing through the tree yet. Between 80 and 100% water tree depth, it is observed that for both solutions there is a rapid decrease in resistance because the tree is now impacting the conductor shield and allowing shorting current to flow through it. The relationship between resistance and temperature was found to be inversed. As temperature increases, resistance decreases for both water tree solutions with distilled water having the most apparent difference across the range of temperatures simulated. With regards to water tree AR, resistance was found to be relatively independent at depth percentages of 10, 50, and 70, but there was some relationship at 100% depth for both solutions. At this depth, as water tree AR increased, resistance was found to also increase. Through the examination of voltage and electric field distribution plots, it was confirmed for EPR cables that increasing water tree depth leads to increasing distortion. The increase in shorting current at 100% depth for aqueous copper sulfate compared to distilled water was also visualized. Lastly, the relationship between shorting current and temperature for narrow water trees was visualized and validated using the distribution plots. This study has led to a better understanding of the effect of cable temperature, water tree solution, and geometry of the tree region on cable degradation. Furthermore, this study has documented simulation results of water tree degradation in EPR MV cables, which had been previously lacking in information. Future work in this area will add frequency into the mix and examine what effect it has on the rate of cable degradation as a result of water treeing. To compliment the work being carried by UMD on different polymer chemistries, harvested medium voltage (MV) cables were selected for accelerated aging study at ORNL. These medium voltage cables are representative of the vintages and materials that are currently in use at existing nuclear power plants (NPPs). These cables were immersed in water at 90°C and energized for a period of two years at elevated voltage. Partial discharge was measured during this period of time to track potential degradation. Unfortunately, due to the absence of failure, the extent of integration between the UMD techniques and the harvested MV insulation was limited. However, based on the findings and from the UMD research in the previous sections on polyethylene and polypropylene, follow-on characterization for harvested MV cable insulation should focus on sample preparation to take advantage of the UMD techniques to better understand mechanistic degradation in MV cable insulations. This would include: 1.) Utilization of solutions and impedance measurement insulation to track permeability in harvested and accelerated aging insulation samples, 2.) Adaptation of gold nanoparticle synthesis for electrical impedance spectroscopy with supporting SEM characterization to study pore formation in harvested and accelerated aged insulation samples, & 3.) ATR-FTIR crystallinity with multiple temperature dependent harvested MV cable insulation under different air and submerged accelerated aging. In addition, low voltage dielectric spectroscopy and high voltage dissipation factor, tan δ, could be implemented as a condition monitoring tool for harvested MV cable insulation in submerged environments.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2017

Factors affecting the bacterial community composition and heterotrophic production of Columbia River estuarine turbidity maxima

Estuarine turbidity maxima (ETM) function as hotspots of microbial activity and diversity in estuaries, yet, little is known about the temporal and spatial variability in ETM bacterial community composition. To determine which environmental factors affect ETM bacterial populations in the Columbia River estuary, we analyzed ETM bacterial community composition (Sanger sequencing and amplicon pyrosequencing of 16S rRNA gene) and bulk heterotrophic production (3H-leucine incorporation rates). We collected water 20 times to cover five ETM events and obtained 42 samples characterized by different salinities, turbidities, seasons, coastal regimes (upwelling vs. downwelling), locations, and particle size. Spring and summer populations were distinct. All May samples had similar bacterial community composition despite having different salinities (1–24 PSU), but summer non-ETM bacteria separated into marine, freshwater, and brackish assemblages. Summer ETM bacterial communities varied depending on coastal upwelling or downwelling conditions and on the sampling site location with respect to tidal intrusion during the previous neap tide. In contrast to ETM, whole (>0.2 μm) and free-living (0.2–3 μm) assemblages of non-ETM waters were similar to each other, indicating that particle-attached (>3 μm) non-ETM bacteria do not develop a distinct community. Brackish water type (ETM or non-ETM) is thus a major factor affecting particle-attached bacterial communities. Heterotrophic production was higher in particle-attached than free-living fractions in all brackish waters collected throughout the water column during the rise to decline of turbidity through an ETM event (i.e., ETM-impacted waters). However, free-living communities showed higher productivity prior to or after an ETM event (i.e., non-ETM-impacted waters). This study has thus found that Columbia River ETM bacterial communities vary based on seasons, salinity, sampling location, and particle size, with the existence of three particle types characterized by different bacterial communities in ETM, ETM-impacted, and non-ETM-impacted brackish waters. Taxonomic analysis suggests that ETM key biological function is to remineralize organic matter.

Microbiology↗