Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “maximum”

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

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

At least 253 records · Page 14

High-throughput characterization, correlation, and mapping of leaf photosynthetic and functional traits in the soybean ( Glycine max ) nested association mapping population

We report photosynthesis is a key target to improve crop production in many species including soybean [Glycine max (L.) Merr.]. A challenge is that phenotyping photosynthetic traits by traditional approaches is slow and destructive. There is proof-of-concept for leaf hyperspectral reflectance as a rapid method to model photosynthetic traits. However, the crucial step of demonstrating that hyperspectral approaches can be used to advance understanding of the genetic architecture of photosynthetic traits is untested. To address this challenge, we used full-range (500–2,400 nm) leaf reflectance spectroscopy to build partial least squares regression models to estimate leaf traits, including the rate-limiting processes of photosynthesis, maximum Rubisco carboxylation rate, and maximum electron transport. In total, 11 models were produced from a diverse population of soybean sampled over multiple field seasons to estimate photosynthetic parameters, chlorophyll content, leaf carbon and leaf nitrogen percentage, and specific leaf area (with R 2 from 0.56 to 0.96 and root mean square error approximately <10% of the range of calibration data). We explore the utility of these models by applying them to the soybean nested association mapping population, which showed variability in photosynthetic and leaf traits. Genetic mapping provided insights into the underlying genetic architecture of photosynthetic traits and potential improvement in soybean. Notably, the maximum Rubisco carboxylation rate mapped to a region of chromosome 19 containing genes encoding multiple small subunits of Rubisco. We also mapped the maximum electron transport rate to a region of chromosome 10 containing a fructose 1,6-bisphosphatase gene, encoding an important enzyme in the regeneration of ribulose 1,5-bisphosphate and the sucrose biosynthetic pathway. The estimated rate-limiting steps of photosynthesis were low or negatively correlated with yield suggesting that these traits are not influenced by the same genetic mechanisms and are not limiting yield in the soybean NAM population. Leaf carbon percentage, leaf nitrogen percentage, and specific leaf area showed strong correlations with yield and may be of interest in breeding programs as a proxy for yield. This work is among the first to use hyperspectral reflectance to model and map the genetic architecture of the rate-limiting steps of photosynthesis.

54 ENVIRONMENTAL SCIENCES↗

Direct analysis of the broad-line SN 2019ein: connection with the core-normal SN 2011fe

Type Ia supernovae (SNe Ia) are important cosmological probes and contributors to galactic nucleosynthesis, particularly of the iron group elements. To improve both their reliability as cosmological probes and to understand galactic chemical evolution, it is vital to understand the binary progenitor system and explosion mechanism. The classification of SNe Ia into Branch groups has led to some understanding of the similarities and differences among the varieties of observed SNe Ia. Branch groups are defined by the pseudo-equivalent widths of the two prominent Si ii lines, leading to four distinct groups: core-normal (CN), shallow-silicon, cool, and broad-line (BL). However, partly due to small sample size, little work has been done on the BL group. We perform direct spectral analysis on the pre-maximum spectra of the BL SN 2019ein, comparing and contrasting to the CN SN 2011fe. Both SN 2019ein and SN 2011fe were first observed spectroscopically within 2 d of discovery, allowing us to follow the spectroscopic evolution of both SNe in detail. We find that the optical depths of the primary features of both the CN and BL SNe are very similar, except that there is a Doppler shift between them. We further examine the BL group and show that for nine objects with pre-maximum spectra in the range (–6)–(–2) d with respect to B-maximum all the emission peaks of the Si ii λ6355 line of BL are blueshifted pre-maximum, suggesting a possible classification criterion.

79 ASTRONOMY AND ASTROPHYSICS↗

Machine-learning-based inversion of nuclear responses

A microscopic description of the interaction of atomic nuclei with external electroweak probes is required for elucidating aspects of short-range nuclear dynamics and for the correct interpretation of neutrino oscillation experiments. Nuclear quantum Monte Carlo methods infer the nuclear electroweak response functions from their Laplace transforms. Inverting the Laplace transform is a notoriously ill-posed problem; and Bayesian techniques, such as maximum entropy, are typically used to reconstruct the original response functions in the quasielastic region. In this work, we present a physics-informed artificial neural network architecture suitable for approximating the inverse of the Laplace transform. Utilizing simulated, albeit realistic, electromagnetic response functions, we show that this physics-informed artificial neural network outperforms maximum entropy in both the low-energy transfer and the quasielastic regions, thereby allowing for robust calculations of electron scattering and neutrino scattering on nuclei and inclusive muon capture rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Resonant multiscalar production in the generic complex singlet model in the multi-TeV region

We develop benchmarks for resonant discalar production in the generic complex singlet scalar extension of the Standard Model (SM) with no additional symmetries, which contains two new scalars. These benchmarks maximize discalar resonant production modes at future 𝑝⁢𝑝 colliders: 𝑝⁢𝑝 → ℎ 2 → ℎ 1⁡ ℎ 1 , 𝑝⁢𝑝 → ℎ 2 → ℎ 1 ⁡ℎ 3 , and 𝑝⁢𝑝 → ℎ 2 → ℎ 3 ⁡ℎ 3 , where ℎ 1 is the observed SM-like Higgs boson and ℎ 2,3 are new scalars. The decays ℎ 2 → ℎ 1 ⁡ℎ 3 and ℎ 2 → ℎ 3 ⁡ℎ 3 may be the only way to discover ℎ 3 , leading to a discovery of two new scalars at once. Current LHC and projected future collider (HL-LHC, FCC-ee+HL-LHC, IL⁢C⁢500+HL-LHC) constraints on this model are used to produce benchmarks at the HL-LHC for ℎ 2 masses between 250 GeV and 1 TeV and a future 𝑝⁢𝑝 collider (FCC-hh) for ℎ 2 masses between 250 GeV and 12 TeV. We update the current LHC bounds on the singlet-Higgs boson mixing angle for these benchmarks. As the mass of ℎ 2 approaches the multi-TeV region, certain limiting behaviors of the maximum rates are uncovered due to theoretical constraints on the parameters. These limits, which can be derived analytically, are BR⁡(ℎ 2 → ℎ 1⁡ ℎ 1 ) → 0.25, BR⁡(ℎ 2 → ℎ 3 ⁡ℎ 3 ) → 0.5, and BR⁡(ℎ 2 → ℎ 1⁡ ℎ 3 ) → 0. It can also be shown that the maximum rates of 𝑝⁢𝑝 → ℎ 2 → ℎ 1 ⁡ℎ 1 and 𝑝⁢𝑝 → ℎ 2 → ℎ 3⁡ ℎ 3 approach the same value. Hence, all three ℎ 2 → ℎ 𝑖 ⁡ℎ 𝑗 decays are promising discovery modes for ℎ 2 masses at and below 𝒪⁡(1 TeV), while above 𝒪⁡(1 TeV) the decays ℎ 2 → ℎ 1 ⁡ℎ 1 and ℎ 2 → ℎ 3 ⁡ℎ 3 are more encouraging. We choose benchmark masses for ℎ 3 to produce a large range of decay signatures including multi-𝑏, multivector boson, and multi-SM-like Higgs production. As we will show, the behavior of the maximum rates leads to the surprising conclusion that in the multi-TeV region this model may be discovered in the Higgs quartet production mode via ℎ 2 → ℎ 3 ⁡ℎ 3 → 4⁢ℎ 1 decays before Higgs triple production is observed. The maximum di- and four Higgs production rates are similar in the multi-TeV range.

Extensions of Higgs sector↗

Optimal Power Sharing Speed Compensation in On-road Wireless EV Charging Systems

Dynamic wireless charging of electric vehicles (EV) is an emerging technology with the potential to address range anxiety and reduce the size of batteries or provide charge-sustaining operation. Charging demand for dynamic wireless charging systems (DWCS) varies greatly in response to location-specific traffic behaviors including the number and speed of vehicles. This paper highlights the potential reduction of load variation with speed compensation and simulates opportunities to maximize the number of cars charged concurrently through “power sharing” or altering power output to slower cars while maintaining a maximum delivered energy. An improved sub-minute model for synthetic traffic is proposed to effectively model DWCS load on high speed roadways and a power electronics model is created based on an existing prototype developed by ORNL to investigate the potential for power sharing. Reductions in the average speed of traffic can greatly increase instantaneous DWCS load by as much as 26%, complicating capacity sizing. Parametric studies with a LCC-S power electronics model shows feasible power reduction of more than 20%, while maintaining a maximum achievable efficiency greater than 90%. Speed compensation on a roadway with large speed variation can reduce average power expected by 21% and increase the maximum number of cars charged simultaneously by 30%. The application of power sharing may significantly reduce load variability due to speed, allow for increased car hosting capability, and guarantee a maximum energy delivered.

Lewis, Donovin D.↗

Residual stresses in Cu matrix composite surface deposits after laser melt injection

Abstract Tungsten carbide particles reinforced metal matrix composite (MMC) coatings can significantly improve surface wear resistance owing to their increased surface hardness. However, the presence of macro‐ and micro‐residual stresses in MMC coatings can have detrimental effects, such as reducing service life. In this study, neutron diffraction was used to determine the residual stresses in spherical fused tungsten carbide (sFTC) reinforced Cu matrix composite surface deposits after laser melt injection. We also developed a thermo‐mechanical coupled finite element model to predict residual stresses. Our findings reveal that sFTC/Cu composite deposits produced with a preheating temperature of 400°C have low residual stresses, with a maximum tensile residual stress of 98 MPa in the Cu matrix on the top surface. In contrast, the sFTC/bronze (CuAl10Ni5Fe4) composite deposit exhibits very high residual stresses, with a maximum tensile residual stress in the Cu matrix on the top surface reaching 651 MPa. These results provide a better understanding of the magnitudes and distributions of residual stresses in sFTC‐reinforced Cu matrix composite surface deposits manufactured via laser melt injection.

Zhang, Xingxing↗

Thermomechanical Stress and Creep-Fatigue Analysis of a High-Temperature Prototype Receiver for Heating Particles

This work presents a three-dimensional (3D) thermomechanical model of a prototype-scale enclosed light trapping solar receiver for heating particles. Results of the thermoelastic model are used to estimate receiver lifetime under maximum flux conditions. A computational fluid dynamics (CFD) model is first developed to predict the temperature fields in a multi-panel assembly under steady operating conditions. Solar flux distributions on the receiver are obtained from the software package SolTrace and applied to the 3D thermal model. The subsequent particle heating is captured through a simplified 1D energy balance. Panel reradiation is considered through a surface-to-surface radiation model and natural convection loss to the surrounding air is captured in a representative fluid domain surrounding the receiver. The resulting temperature fields from the CFD analysis are used as inputs for a thermoelastic mechanical model with representative boundary conditions. With the resultant temperature and stress fields, a creep-fatigue damage and lifetime analysis is performed using the linear damage accumulation (LDA) theory. The Manson-Coffin formula and Larson Miller correlation are used to calculate the fatigue and creep, respectively. A maximum damage (corresponding to a 30-year service life) is defined for design assessment. The model was first developed and verified in detail by comparing with published results in the literature (temperature and stress profiles and distributions, and creep/fatigue damage fractions) for tubular solar receivers with supercritical carbon dioxide as the working fluid. It was then implemented to model a planar-cavity receiver with various design parameters. Specifically, three different design geometries are considered, and the results show that a maximum temperature of approximately 1200 K could be reached for each design with the given incident solar flux, with the main difference being the distribution of these temperatures. Preliminary resulting stresses for the small-scale prototype without design optimization vary from 20 MPa to 250 MPa for each design, with the maximum stresses occurring on the front face and concave geometry on the side of the panel. In future work, the developed methodology shown here will be applied to analyze a full-scale (50-150 MWth) receiver.

concentrated solar power↗

Clumped isotope constraints on changes in latest Pleistocene hydroclimate in the northwestern Great Basin: Lake Surprise, California

Abstract During the Last Glacial Maximum (LGM) and subsequent deglaciation, the Great Basin in the southwestern United States was covered by numerous extensive closed-basin lakes, in stark contrast with the predominately arid climate observed today. This transition from lakes in the Late Pleistocene to modern aridity implies large changes in the regional water balance. Whether these changes were driven by increased precipitation rates due to changes in atmospheric dynamics, decreased evaporation rates resulting from temperature depression and summer insolation changes, or some combination of the two remains uncertain. The factors contributing to these large-scale changes in hydroclimate are critical to resolve, given that this region is poised to undergo future anthropogenic-forced climate changes with large uncertainties in model simulations for the 21st century. Furthermore, there are ambiguous constraints on the magnitude and even the sign of changes in key hydroclimate variables between the Last Glacial Maximum and the present day in both proxy reconstructions and climate model analyses of the region. Here we report thermodynamically derived estimates of changes in temperature, precipitation, and evaporation rates, as well as the isotopic composition of lake water, using clumped isotope data from an ancient lake in the northwestern Great Basin, Lake Surprise (California). Compared to modern climate, mean annual air temperature at Lake Surprise was 4.7 °C lower during the Last Glacial Maximum, with decreased evaporation rates and similar precipitation rates to modern. During the mid-deglacial period, the growth of Lake Surprise implied that the lake hydrologic budget briefly departed from steady state. Our reconstructions indicate that this growth took place rapidly, while the subsequent lake regression took place over several thousand years. Using models for precipitation and evaporation constrained from clumped isotope results, we determine that the disappearance of Lake Surprise coincided with a moderate increase in lake temperature, along with increasing evaporation rates outpacing increasing precipitation rates. Concomitant analysis of proxy data and climate model simulations for the Last Glacial Maximum are used to provide a robust means to understand past climate change, and by extension, predict how current hydroclimates may respond to expected future climate forcings. We suggest that an expansion of this analysis to more basins across a larger spatial scale could provide valuable insight into proposed climate forcings, and aid in climate model process depiction. Ultimately, our analysis highlights the importance of temperature-driven evaporation as a mechanism for lake growth and retreat in this region.

Geology↗

Charging and ion ejection dynamics of large helium nanodroplets exposed to intense femtosecond soft xray pulses

Ion ejection from charged helium nanodroplets exposed to intense femtosecond soft X-ray pulses is studied by single-pulse ion time-of-flight (TOF) spectroscopy in coincidence with small-angle X-ray scattering. Scattering images encode the droplet size and absolute photon flux incident on each droplet, while ion TOF spectra are used to determine the maximum ion kinetic energy, Ekin, of He+ j fragments (j = 1-4). Measurements span HeN droplet sizes between N similar to 10(7) and similar to 10(10) (radii R-0 = 78-578 nm), and droplet charges between similar to 9x10(-5) and similar to 3x10(-3) e/atom. Conditions encompass a wide range of ionization and expansion regimes, from departure of all photoelectrons from the droplet, leading to pure Coulomb explosion, to substantial electron trapping by the electrostatic potential of the charged droplet, indicating the onset of hydrodynamic expansion. The unique combination of absolute Xray intensities, droplet sizes, and ion Ekin on an event-by-event basis reveals a detailed picture of the correlations between the ionization conditions and the ejection dynamics of the ionic fragments. The maximum E-kin of He+ is found to be governed by Coulomb repulsion from unscreened cations across all expansion regimes. The impact of ion-atom interactions resulting from the relatively low charge densities is increasingly relevant with less electron trapping. The findings are consistent with the emergence of a charged spherical shell around a quasineutral plasma core as the degree of ionization increases. The results demonstrate a complex relationship between measured ion E-kin and droplet ionization conditions that can only be disentangled through the use of coincident single-pulse TOF and scattering data.

Helium↗

Efficient Probabilistic Prediction and Uncertainty Quantification of Tropical Cyclone–Driven Storm Tides and Inundation

Abstract This study proposes and assesses a methodology to obtain high-quality probabilistic predictions and uncertainty information of near-landfall tropical cyclone–driven (TC-driven) storm tide and inundation with limited time and resources. Forecasts of TC track, intensity, and size are perturbed according to quasi-random Korobov sequences of historical forecast errors with assumed Gaussian and uniform statistical distributions. These perturbations are run in an ensemble of hydrodynamic storm tide model simulations. The resulting set of maximum water surface elevations are dimensionality reduced using Karhunen–Loève expansions and then used as a training set to develop a polynomial chaos (PC) surrogate model from which global sensitivities and probabilistic predictions can be extracted. The maximum water surface elevation is extrapolated over dry points incorporating energy head loss with distance to properly train the surrogate for predicting inundation. We find that the surrogate constructed with third-order PCs using elastic net penalized regression with leave-one-out cross validation provides the most robust fit across training and test sets. Probabilistic predictions of maximum water surface elevation and inundation area by the surrogate model at 48-h lead time for three past U.S. landfalling hurricanes (Irma in 2017, Florence in 2018, and Laura in 2020) are found to be reliable when compared to best track hindcast simulation results, even when trained with as few as 19 samples. The maximum water surface elevation is most sensitive to perpendicular track-offset errors for all three storms. Laura is also highly sensitive to storm size and has the least reliable prediction. Significance Statement The purpose of this study is to develop and evaluate a methodology that can be used to provide high-quality probabilistic predictions of hurricane-induced storm tide and inundation with limited time and resources. This is important for emergency management purposes during or after the landfall of hurricanes. Our results show that sampling forecast errors using quasi-random sequences combined with machine learning techniques that fit polynomial functions to the data are well suited to this task. The polynomial functions also have the benefit of producing exact sensitivity indices of storm tide and inundation to the forecasted hurricane properties such as path, intensity, and size, which can be used for uncertainty estimation. The code implementing the presented methodology is publicly available on GitHub.

54 ENVIRONMENTAL SCIENCES↗

A Polarimetric Radar Analysis of Ice Microphysical Processes in Melting Layers of Winter Storms Using S-Band Quasi-Vertical Profiles

Quasi-vertical profiles (QVPs) obtained from a database of U.S. WSR-88D data are used to document polarimetric characteristics of the melting layer (ML) in cold-season storms with high vertical resolution and accuracy. A polarimetric technique to define the top and bottom of the ML is first introduced. Using the QVPs, statistical relationships are developed to gain insight into the evolution of microphysical processes above, within, and below the ML, leading to a statistical polarimetric model of the ML that reveals characteristics that reflectivity data alone are not able to provide, particularly in regions of weak reflectivity factor at horizontal polarization Z H . QVP ML statistics are examined for two regimes in the ML data: Z H ≥ 20 dB Z and Z H < 20 dB Z . Regions of Z H ≥ 20 dB Z indicate locations of MLs collocated with enhanced differential reflectivity Z DR and reduced copolar correlation coefficient ρ hv , while for Z H < 20 dB Z a well-defined ML is difficult to discern using Z H alone. Evidence of large Z DR up to 4 dB, backscatter differential phase δ up to 8°, and low ρ hv down to 0.80 associated with lower Z H (from -10 to 20 dB Z ) in the ML is observed when pristine, nonaggregated ice falls through it. Positive correlation is documented between maximum specific differential phase K DP and maximum Z H in the ML; these are the first QVP observations of K DP in MLs documented at S band. Negative correlation occurs between minimum ρ hv in the ML and ML depth and between minimum ρ hv in the ML and the corresponding enhancement of Z H (Δ Z H = Z H max - Z H rain ).

Meteorology & Atmospheric Sciences↗

Behavior of Tidally Triggered Earthquakes Depends on Fluid Conditions

Seismicity at The Geysers geothermal field and in Oklahoma is heavily influenced by industrial activities related to energy production, though the mechanism in which earthquakes are induced or triggered is different. At The Geysers, much of the seismicity is linked to thermoelastic stresses caused by injecting cold water into hot rocks, whereas in Oklahoma, the seismicity is linked to a reduction of confining stress on faults due to increasing pore pressure resulting from wastewater injections. Here, in this study, we show that these contrasting conditions are also evident in tidally triggered earthquakes. At The Geysers, earthquakes preferentially occur during maximum extensional stress, which does not occur at the same time as maximum shear stress on optimally oriented faults in the regional stress field. In Oklahoma, earthquakes preferentially occur during maximum shear stress on optimally oriented faults rather than maximum extensional stress. The magnitude of tidal extensional stress is naturally much greater than tidal shear stress. However, in a fluid-saturated environment, pore pressure responds to changes in volume, which can counteract or reduce the effect of the applied stress. The difference in behavior at these two sites is indicative of the level of coupling between applied stress and pore pressure, corresponding to unsaturated conditions at The Geysers and high pore pressure in Oklahoma.

58 GEOSCIENCES↗

Inverting infrasound data for the seismoacoustic source time functions and surface spall at the Source Physics Experiments Phase II: Dry Alluvium Geology

This report presents the infrasound data recorded as part of the Source Physics Experiment - Phase 2, Dry Alluvium Geology. This experiment, also known colloquially as DAG, consisted of four underground chemical explosions at the Nevada National Security Site. We focus our analysis on only the fourth explosion (DAG-4) as we determined that this was the only event that produced clear source-generated infrasound energy as recorded by the DAG sensors. We analyze the data using two inversion methods. The first method is designed to estimate the point-source seismoacoustic source time functions, and the second inversion method is designed to estimate the first-order characteristics (e.g. horizontal dimensions and maximum amplitude) of the actual spall surface. For both analysis methods, we are able to fit the data reasonably well, with various assumptions of the source model. The estimated seismoacoustic source appears to be a combination of a buried, isotropic explosion with a maximum amplitude of ~2 x 10 9 Nm and a vertically oriented force, applied to the Earth's surface with a maximum amplitude of 4 x 10 7 N. We use the vertically oriented force to simulate surface spall. The estimated spall surface has an approximate radius of ~40 m with a maximum acceleration magnitude in the range of 0.8 to 1.5 m/s/s. These estimates are approximately similar to the measured surface acceleration at the site.

58 GEOSCIENCES↗

Best Practices for Shielding Analyses of Activated Metals and Spent Resins from Reactor Operation

This report investigated best practices for performing shielding evaluations of Type B waste packages, as defined in 10 CFR Part 71 regulations on packaging and transportation of radioactive material, or packages for which the contents are not defined before loading and may include a broad range of nuclides, geometries, and non-fuel materials. The following non-fuel waste streams were analyzed: Activated metals from decommissioned commercial power reactors, including Type 304 stainless steel, reactor vessel steel, and Inconel, Control blades from boiling water reactors, Neutron-activated corrosion products on surfaces of activated metals, and Spent resins from power plant operations. Measured elemental compositions, including major constituents and impurities, for steel and Inconel samples from commercial power reactors were used in activation calculations to determine radionuclide inventories in activated metals. For a simplified cask model, 60 Co contribution to the total external package dose rate at 30 days after shutdown varied from approximately 60% to 95%, depending on the activated metal, initial cobalt impurity concentration in the metal, and the thickness of the overpack gamma shield. Its maximum contribution to the total external package dose rate of approximately 100% was reached within the time interval of 2 to 5 years after shutdown and was maintained for up to 45 to 60 years after shutdown, depending on material, initial cobalt impurity concentration, and shield thickness. Thereafter, the 60 Co contribution to external package dose rate decreased with increasing decay time. Cobalt-60 is primarily produced by neutron reactions with the cobalt impurity in steel and Inconel. Other important radionuclides in activated metals contributing to package external dose rate are radionuclides with relatively short decay times, including 51 Cr, 59 Fe, 58 Co, and 54 Mn. These radionuclides may be represented as an equivalent 60 Co activity/source because 60 Co gamma ray emissions are bounding in terms of source strength and energy to other important radionuclides identified in the analyzed activated metals. Approaches for modeling the neutron-activated corrosion products that may be attached to activated reactor components were analyzed in this report. It was demonstrated that a surface source is more conservative than a uniform volumetric source for the treatment of neutron-activated corrosion products with respect to external package dose rates. An analysis of the maximum radionuclide loadings reported on spent resins identified the radionuclides 137 Cs, 60 Co, 134 Cs, 65 Zn, and 58 Co as the primary contributors to external package dose rate. For a resin cooled for 3.08 years, the external package dose rate was entirely produced by the reported 137 Cs and 60 Co inventory. The neutron sources from actinides found on spent resins or activated metals produced negligible dose rates and may be ignored in dose rate analyses. Effects of idealized waste material, source geometry, and spatial material/source distributions on external package dose rates were determined based on dose rate results for a simplified cask model under normal conditions of transport. Type 304 stainless steel, zirconium, and aluminum with adjusted mass densities based on a maximum content weight were analyzed for material modeling effects on external package dose rate. These materials produced identical external package dose rates within the statistical uncertainties of the dose rate estimates. Among four different source geometry configurations with homogeneous material of different mass densities, uniform volumetric source distribution, and the same source strength, the geometry configuration with lowest mass density (i.e., minimum self-shielding effect among the four cases) was most conservative. Spatial source distributions that better represent localized peak 60 Co activity values were more conservative than a uniform volumetric source distribution, assuming the same weight and total source strength per package. The increase in external dose rate caused by localized source peaks can be as much as the ratio between source peak activity density to average activity density, depending on the location of the activated metal with peak activity density. Therefore, the shielding analysis may be simplified if localized peak activities and the average activity per package can be measured/determined and documented at the time of cask loading. For simplicity, external package dose rates may be determined based on average source activity and a uniform volumetric source distribution. The dose rate results from that calculation model multiplied by the ratio of peak activity density to average activity density will produce maximum dose rate values for conservative estimates.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Tank 27 Solids Disturbance Testing Using M-Star for the Impact of CSMP Usage on Blend Tank Mixing Effectiveness

Savannah River Mission Completion (SRMC) plans to operate Tank 27as a blend tank for future waste disposal operations. It is desirable to avoid additional worker exposure by using the planned commercial submersible mixing pump (CSMP) installations in Risers B2 and B4 to function as the blend pumps for prolonged blend tank operations. The purpose of this task is to evaluate the impacts to mixing across a range of pump installation heights from 1 to 201 inches above the tank floor. The CSMPs shall maintain adequate blending within a 24 hour maximum timeframe and an adequate nozzle discharge velocity as outlined in U-ESR-G-00030, using the specifications of the planned CSMPs listed in PO SRRA149971 and M-DS-G-00081. Specifications include: maximum pump diameter of D=22.5 in., dual nozzles placed tangentially opposed with nominal inside diameter of 2 in., total flow range (both nozzles combined) of 500 to 1950 GPM and maximum RPM of 1800. Tank 27 is expected to have a fill height maximum of 360". This analysis began by using results from the 2010 testing and applying them to the Tank 27 design to calculate a mixing time as a function of pump flow rate. These results were modified to account for the higher liquid level in Tank 27. Once the liquid level was increased to 360 inches, the results were modified to account for a higher viscosity in Tank 27 than in the testing. Once the viscosity was increased, the effect of changing pump elevation was applied to calculations of the mixing time in Tank 27. Finally, the uncertainty from the testing was included in the analysis to recommend a blend time as a function of pump flow rate and elevation. The conclusion from this study is that the particles in tank 27 will be suspended at both the 600 and 900 gpm flow rates. Only particles above 90 micron will not be suspended in significant amounts when the flow rate is 600 gpm. At 900 gpm even 100 micron particles were thoroughly suspended in the tank. From the series of simulations it has been concluded that it will not be possible to mix tank 27 at the proposed fluid flow rates and not suspend the solid particles at the bottom of the tank.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Geophysical Characterization and Monitoring of the 200 Area Treated Effluent Disposal Facility to Support Permit Renewal; TEDF Geophysical Monitoring

The Treated Effluent Disposal Facility (TEDF), located in the 200 East Area of the Hanford Site, is a site where non-hazardous and non-radioactive liquid wastes are disposed into two state-permitted infiltration basins. In 2016, the Washington State Department of Ecology denied a permit renewal request for TEDF due to the inability to adequately assess the impact of TEDF discharge water on the underlying groundwater quality. TEDF overlies the relatively impermeable Ringold Lower Mud (RLM) unit, whose upper contact lies in the vadose zone approximately 30 m below ground surface and approximately 10 m above the water table boundary. The RLM is assumed to isolate TEDF discharge water from the natural groundwater aquifer, which is monitored using wellbores screened below the RLM. Therefore, samples collected from monitoring wells near TEDF are not considered representative of TEDF discharge water. Rather, TEDF discharge water is assumed to mound on top of the RLM to form a perched aquifer. To support permitting of the TEDF, a new monitoring well is required that can be used to sample water from the presumed perched water aquifer above the RLM. Ideally, the screened section of the well would be located at the peak of the perched water mound(s), which presumably occurs at the point of maximum vertical flux from TEDF to the RLM, or equivalently where the dominant infiltration flow paths reach the RLM. This report describes how time-lapse 3D electrical resistivity tomography (ERT) was used in conjunction with nominal TEDF discharge operations to image the dominant flow paths from each pond to the RLM. Results are summarized in Figure ES.1. Figure ES.1A shows a satellite image of the TEDF overlain by an array of surface ERT electrodes. The solid and dashed black circles denote the zones of maximum vertical flux at the RLM within the south and north ponds, respectively, and presumably the regions where perched water peaks during discharge. Figure ES.1B shows time-lapse difference imaging results approximately 15 days after switching discharge from the south pond to the north pond. Blue iso-surfaces beneath the north pond denote regions of increasing bulk electrical conductivity caused by increasing saturation due to infiltrating water. Red iso-surfaces beneath the south pond denote regions that were previously saturated during the south pond discharge and are now de-saturating, causing a decrease in bulk electrical conductivity. In both cases, the zones of maximum change mark the dominant flow paths to the RLM. Figure ES.1C shows time-lapse imaging results approximately 15 days after switching discharge from the north pond to the south pond. In this case, blue iso-surfaces mark the dominant flow paths to the RLM from the south pond. The red iso-surfaces mark the dominant flow paths that existed during discharge to the north pond. The regions of maximum change in bulk conductivity (due to saturation or desaturation) that mark the primary flow paths are equivalent in both cases and denoted by the dashed solid and black circles. If a perched water zone forms on the RLM, it is likely to mound within or near the dashed circle during discharge to the north pond, and in the solid circle during discharge to the south pond. In other words, if perched water mounding occurs, the circles mark the optimum locations suggested by the ERT imaging for monitoring boreholes to be placed, enabling samples to be collected that are representative of TEDF discharge impacts on groundwater quality.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Seasonal Cost-Benefit Analysis of Automated Distribution Feeder Upgrades with Advanced Mitigation Technologies

The increasing deployment of distributed solar photovoltaics (DPV) to meet clean energy goals can trigger adverse grid operation issues, such as voltage excursions and the violation of thermal loading constraints of the power delivery elements (e.g., lines and transformers) on the evolving electricity infrastructure. Such integration issues would require distribution upgrades with associated costs to mitigate them and to maintain reliable and resilient grid operating conditions. Traditional distribution network upgrade approaches use a specific single snapshot analysis that is overly conservative. This study considers a multi-time point analysis to capture both moderate (probable bounds) and extreme grid operating conditions using time points such as minimum load with minimum photovoltaics (PV), maximum load with maximum PV, maximum load with minimum PV, and minimum load with maximum PV. Further, this study investigates seasonal variation impacts and associated distribution upgrade costs for a spring season case (March, representing a low load and high PV scenario) and a summer case (July, representing a high load and high PV scenario). Such seasonal analysis will allow system operators to characterize upgrade requirements and associated costs across various periods. Because the spatial distribution of DPV can impact upgrade and associated costs, this study investigates three common DPV deployment scenarios - randomly deployed, close to the substation, and far from the substation - at different penetration levels. Apart from spatial distribution impacts, this project evaluates the techno-economic impacts of the nodal photovoltaic penetration factor (NPPF) for generating the various DPV deployment scenarios at increasing penetration levels. This project investigates the impact of varying nodal PV-to-load ratios using conservative and extreme NPPF values of 3 and 10, respectively. This study investigates the deployment of traditional infrastructure upgrade strategies, such as installing new voltage regulating equipment, transformers and lines replacements, and the activation of advanced inverter functionality (e.g., autonomous volt/VAR) in expanding PV hosting capacity. Existing DPV systems are assumed to operate with the legacy unity power factor, and we considered the possibility of retrofitting such systems with the activation of volt/VAR control as integration standards and regulations continue to evolve to allow such functions. The cost-benefit analysis metrics used in study include distribution upgrade costs, average cost per watt of the upgrade cost, average marginal cost per watt of the upgrade cost, and power losses.

14 SOLAR ENERGY↗

Parametric and Sensitivity Analysis of a Steam Generator Model Using Python and Machine-Learning Tools

For this study, we used Python and machine-learning tools to perform a comprehensive parametric and sensitivity analysis on a steam generator (SG) model. (The Python model was based on a previously completed MATLAB framework for the Holtec SMR-160 SG.) We investigated the influence of various input parameters (e.g., heat transfer coefficient [HTC], Nusselt number, and heat exchanger effectiveness) on the system’s output. With machine-learning tools such as the Risk Analysis Virtual Environment (RAVEN), which was developed at Idaho National Laboratory, we were then able to perform an automated analysis of the SG inputs’ effect on the HTC. The analysis results give valuable insights into the performance and optimization of SG systems. We found the inlet mass flow rate (MFR) to have the greatest impact on the HTC, followed closely by the inlet temperature, and then pressure. Shifting of the input parameters causes the location of the maximum HTC along the SG length to change incrementally. The cold leg (CL) MFR was also found to impact the HTC magnitude as well as the location of the maximum HTC. At between 0.4–0.9 of the total SG length, the input parameters experience maximum impact on the HTC, leading us to suggest that sensors be efficiently placed on the SG so as to closely and effectively monitor thermal-hydraulic properties during reactor operation. We also found that the sensitivity data calculated manually agrees with the RAVEN – based data, confirming the same range of maximum sensitivity. However, the RAVEN-based analysis showed that cold leg pressure and hot leg temperature have a greater impact on the heat transfer coefficient than the mass flow rate, implying that a manual sensitivity study taking only two samples is not accurate.

20 FOSSIL-FUELED POWER PLANTS↗