Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “correspondence”

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

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

At least 1,081 records · Page 60

Pre-Ionization Controlled Laser Plasma Formation for Ignition Applications (Final Scientific/Technical Report)

This report consists of three grand funded papers. Abstracts for each follow. 1) In this work, we present experimental and modeling studies of air pre-ionization using ultraviolet (UV) laser pulses and its effect on laser breakdown of an overlapped near-infrared (NIR) pulse. Experimental studies are conducted with a 266 nm beam (fourth harmonic of Nd:YAG) for UV pre-ionization and an overlapped 1064 nm NIR beam (fundamental of Nd:YAG), both having pulse duration of 10 ns. Results show that the UV beam produces a pre-ionized volume which assists in breakdown of the NIR beam, leading to reduction in NIR breakdown threshold by factor of >2. Numerical modeling is performed to examine the ionization and breakdown of both beams. The modeled breakdown threshold of the NIR, including assist by pre-ionization, is in reasonable agreement with the experimental results. 2) The present contribution compares the energy absorption, optical emission, temperature, and fluid dynamics of ultraviolet (UV) λ = 266 nm and near infrared (NIR) λ = 1064 nm nanosecond laser induced plasmas in ambient air. For UV pulses at the conditions studied, energy absorption by the plasmas increases relatively gradually with laser pulse energy starting at delivered energy of E~8 mJ. Corresponding measurements of plasma luminosity show that the absorption of UV radiation does not necessarily result in visible plasma emission. For the NIR induced plasmas, the energy absorption profile is far more abrupt and begins at ~55 mJ. In contrast with UV, the absorption of NIR radiation is always accompanied by intense optical emission. The temperatures of both types of plasma have been measured with Rayleigh scattering thermometry (at times after the Thomson signal sufficiently diminishes). The UV plasmas can attain a wider range of temperatures, including lower temperatures, depending on the pulse energy (e.g., T ~400–2000K for E ~7–35 mJ at Δt = 10 μs after the pulse) while the NIR plasmas show only hotter temperatures (e.g., T ~ 12 000 K for E = 75 mJ at Δt = 10 μs after the pulse) as is consistent with the literature. Differences in the fluid dynamics for UV versus NIR pulses are shown with Schlieren imaging. In this work, the contrast in the UV and NIR plasma threshold behavior is attributed to differing roles of avalanche ionization and multiphoton ionization as is also illustrated by a simple numerical model. 3) The present contribution examines the impact of plasma dynamics and plasma-driven fluid dynamics on the flame growth of laser ignited mixtures and shows that a new dual-pulse scheme can be used to control the kernel formation process in ways that extend the lean ignition limit. We perform a comparative study between (conventional) single-pulse laser ignition (λ = 1064 nm) and a novel dual-pulse method based on combining an ultraviolet (UV) pre-ionization pulse (λ = 266 nm) with an overlapped near-infrared (NIR) energy addition pulse (λ = 1064 nm). We employ OH* chemiluminescence to visualize the evolution of the early flame kernel. For single-pulse laser ignition at lean conditions, the flame kernel separates through third lobe detachment, corresponding to high strain rates that extinguish the flame. In this work, we investigate the capabilities of the dual-pulse to control the plasma-driven fluid dynamics by adjusting the axial offset of the two focal points. In particular, we find there exists a beam waist offset whereby the resulting vorticity suppresses formation of the third lobe, consequently reducing flame stretch. With this approach, we demonstrate that the dual-pulse method enables reduced flame speeds (at early times), an extended lean limit, increased combustion efficiency, and decreased laser energy requirements.

Mechanical engineering↗

Low-Temperature Plasma Treatment for Enhanced Recovery of Highly Valued Critical REEs from Coal (Final Report)

The principle objective of this project was to develop a novel process flowsheet using low-temperature plasma treatment (LTP) integrated with conventional hydrometallurgical processes to recover rare earth elements (REEs), especially highly valued REEs (i.e. scandium and critical REEs), from coal and coal byproducts. LTP treatment was used to pretreat the solid feed material prior to leaching in an effort to improve the leaching characteristics of the REEs and potentially reduce acid consumption. Initially, a detailed experimental program was performed on a laboratory LTP furnace to identify a set of process parameter values to be used for assessing the enhanced leaching performance achieved on four specific gravity fractions of West Kentucky No. 13 and Fire Clay coarse coal refuse materials. High-temperature oxidation (HTO)) was also evaluated as a means for improving REE leachability on the same fractions of both coal sources. The lixiviants used to assess leachability improvements included: i) untreated water, ii) 0.1 M ammonium sulfate solution and iii) 1.2 M sulfuric acid solution at a temperature of 75°C. The laboratory data was used to develop multiple process flow sheets using HTO only to produce a rare earth oxide mix concentrate or utilizing both LTP and HTO treatment processes to generate a high purity heavy REO mix concentrate and a light REO mix concentrate. A techno-economic analysis was performed on each circuit. Low-temperature plasma oxidation was found to enhance REE leaching characteristics for all specific gravity (SG) fractions of the West Kentucky No. 13 coarse refuse material. The leaching recovery improvement was preferentially achieved for the heavy REEs with higher values obtained when treating the lower density fraction material. For example, the leaching recovery for the heavy REEs in the untreated 1.6 SG Float fraction material was 8% in a 0.1 M ammonium sulfate solution. After LTP treatment for 32 hours, heavy REE recovery increased to 33%. The corresponding recovery value for the light REEs was 10%. These findings indicate a strong affinity of the heavy REEs with the organic material, which is the component being reduced in content during the LTP treatment. Scanning electron micrographs showed the development of a honeycomb structure with increased LTP treatment with higher surface area as confirmed from BET measurements. Higher surface area promotes improved leaching efficiency due to enhance chemical interactions. When using 1.2 M sulfuric acid, the LTP treatment was found to provide a heavy REE recovery increase from 23% to 53% for the 1.6 SG Float fraction which corresponded to an ash content increase from around 20% to 35%. However, LTP treatment provided no improvement in leaching characteristic for the 2.2 SG sink material in the acid solution, which was likely due to the low organic content. The LTP treatment of the Fire Clay SG fractions provided no improvement in REE leaching characteristic, which may be due to the association with REE minerals such as monazite and REE-associated mineral matter such as clays. High temperature oxidation using temperatures in the range of 600-750OC was found to significantly enhance the leaching characteristics of both the light and heavy REEs for all density fractions in both coarse coal refuse materials. Generally, the increase in light REE recovery for the West Kentucky material was from 15% for the untreated material to 90% for the HTO treated material at 600°C. Recovery values dropped significantly when using temperatures below 600°C and above 750°C. The impact of the excessive temperatures was most likely a result of sintering. For the heavy REEs, HTO significantly improved recovery in the lighter SG fractions, which may be due to an affinity with the organic content. Only limited heavy REE recovery improvement was obtained for the denser SG fractions. Another benefit of HTO is the conversion of calcium-based minerals to a very soluble calcium oxide, which is removable using a mild acid solution, and the conversion of Fe minerals to relatively insoluble minerals such as hematite and magnetite. The pregnant leaching solution produced from the leaching of all SG fractions obtained from both coal sources was treated using two stages of solvent extraction procedurally modeling a rougher-cleaner circuit followed by oxalic acid precipitation. Rare earth oxide products having a purity greater than 97% was produced from nearly all of the SG fractions. A novel process flowsheet was developed based on the selective leaching of heavy REEs that was made possible through LTP treatment of West Kentucky No. 13 coarse refuse followed by salt leaching. In the flowsheet, coal refuse material is ground to a top size of 1mm and then treated by LTP. The treated material is mixed into a 0.1 M ammonium sulfate solution and then filtered to produce a heavy REE enriched pregnant leach solution. The heavy REEs are concentrated by two stages of selective precipitation to produce a heavy REE mix oxalate and then roasted to obtain a heavy REO mix concentrate. The solid residues remaining after salt leaching is dried and then subjected to HTO treatment followed by concentration using selective precipitation and roasting to produce a light REO mix concentrate. An alternative circuit was developed using HTO treatment only to treat 100% of the coarse refuse material followed by selective precipitation and roasting to generate a high purity REO mix product. A techno-economic analysis performed for the treatment of West Kentucky No. 13 and Fire Clay coarse refuse based on a hypothetical 500-tph plant revealed that the net present value for all scenarios was significantly negative. Primary reasons include: i) a substantial decrease in the basket value of the final product relative to the feed and ii) high chemical costs relative to the amount of product generated. The basket value reduction was primarily due to little to no scandium recovery due to poor leaching characteristics when using low acid concentrations as necessitated by the low feed grades. Operating costs were primarily influenced by the chemical costs relative to the kilograms of concentrate produced and feed grade. For the Fire Clay coarse refuse material, the total operating cost was estimated to be $354 per kg of REO produced which was higher than the basket value of the REO product ($276/kg). Approximately 37% of the total cost was the sulfuric acid used during the leaching process. A comparative analysis performed to assess the economic viability of using low-temperature plasma furnaces found that market value enhancement for the generation of both heavy REO and light REO concentrates will be needed in addition to significant technology advancements to reduce the cost of plasma furnaces.

01 COAL, LIGNITE, AND PEAT↗

Investigation of Thermolytic Hydrogen Generation Rate in Tank 44 Dissolved Saltcake Samples

Saltcake core samples collected from Tank 44 in 2006 were dissolved to provide material for HGR measurements applicable to F-Area dissolved saltcake material. Additionally, characterization was performed on the Tank 44 saltcake material. The following are key results from the Tank 44 saltcake characterization. The Tank 44 Upper Saltcake Composite, corresponding to the 171 to 285 inch tank level, contained by mass approximately 69% sodium nitrate, 11% sodium carbonate, 8% sodium nitrite, smaller amounts of other salts and components, and 9% unquantified (which includes water, water of hydration, oxygen/hydrogen content of oxides and hydroxides, and uncertainty). The Tank 44 Lower Saltcake Composite, corresponding to the 76 to 114 inch tank level, contained by mass approximately 49% sodium carbonate, 18% sodium nitrate, smaller amounts of other salts, at least 8% sludge, and 9% unquantified (see above). The dissolved saltcake contained free hydroxide less than quantifiable (<0.01 M) due to the limited quantity of material that could be removed from the Shielded Cells based on the sample radioactivity. Measurement by pH paper provided an approximate pH of 12. The following are key results from the Tank 44 HGR testing. During boiling at 106.7 °C, HGR for Tank 44 dissolved saltcake without added glycolate was 7.2×10 -8 ft 3 h -1 ga l-1 . During boiling at 106.9 °C, HGR for Tank 44 dissolved saltcake with 1000 mg/L of added glycolate was 8.2×10 -8 ft 3 h -1 gal -1 . ∙ For the test without added glycolate, the first several HGR measurements at 70, 85, and 100 °C gave indication of the release of dissolved hydrogen and should not be used to represent the sustained thermolytic HGR for those temperatures. The measurements at boiling are the best representation of thermolysis in this testing. Carbon dioxide was observed at concentrations up to 6 vol% in the flow-system offgas for the test at boiling. ∙ Methane generation was observed at 100 °C and boiling. Methane concentration in the total gas generated during testing remained well below the lower flammability limit for methane in air. The addition of 1000 mg/L of glycolate did not have a significant impact on the hydrogen generation rates measured during this testing. The low hydroxide concentration in the Tank 44 dissolved saltcake likely influenced the relatively low thermolytic HGR and high carbon dioxide release observations in this testing. Based on the observation that methane was generated or released upon heating SRS radioactive Tank 44 waste samples to 100 °C and above, we recommend gaining a greater understanding of the cause and mechanism of its generation. First, the applicable literature should be reviewed to reveal the thermolytic methane generation mechanisms of possible methane generating species in the SRS CSTF. If warranted, a plan should be developed for simulant tests with methylated siloxanes and other applicable compounds in order to gain a better mechanistic understanding of methane generation in the SRS CSTF.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

H2OI95: A Stand-Alone Fortran Code for Evaluating the IAPWS-95 Equation-of-State Model for Water (Rev. 1)

H2OI95 is a stand-alone Fortran code for evaluating the IAPWS-95 equation-of-state model for the thermodynamic properties of water. It also evaluates the corresponding thermochemical properties of water consistent with the CODATA recommendations. The IAPWS-95 model is based on a model equation for the dimensionless Helmholtz energy for which the primary variables are the inverse reduced temperature (τ = T cr /T ) and reduced density (δ = ρ/ρ cr ). Here T is the absolute temperature (K), ρ is density (kg/m³), and the subscript “ cr ” refers to the critical point of water (647.096 K and 22.064 MPa pressure in this model, for which ρ cr is 322 kg/m³). The code solves four basic types of problems, distinguished by the specified inputs: 1. Temperature (K) and density (ρ kg/m³) or reduced density (δ) 2. Temperature (K) and pressure (MPa). 3. Temperature (K) on the saturation (liquid-vapor equilibrium) curve 4. Pressure (MPa) on the saturation curve. Each type of problem is run using a corresponding input (text) file. All but the first type of problem require iteration. For example, to solve for desired temperature and pressure, the reduced density must be adjusted to give the desired pressure. Iteration is accomplished using the Newton-Raphson method, though the secant method is also used in solving the fourth type of problem. For the last three types of problems, H2OI95 has been used to conduct numerical studies of convergence and the problem of multiple numerical solutions, some of which are not physically valid. Obtaining valid results depends mainly on appropriate choice of starting values for the density. The default values used in H2OI95 appear to consistently lead to generally desired results. With modification (not addressed here), H2OI95 can be used to support SUPCRT92 and similar codes that compute chemical thermodynamic properties of species and reactions over a wide range of temperature and pressure (273.16-1273K and 0-1000 MPa). Version 1.1 of this software calculates an extended set of water properties, including some derivable solely from IAPWS-95, such as the isobaric thermal expansivity and the isothermal compressibility, and still other properties based on supplementary models sanctioned by IAPWS, including the dynamic viscosity, thermal conductivity, surface tension, and static dielectric constant, and various other properties (notably Debye-Hückel parameters) derivable from the dielectric constant and its partial derivatives with respect to temperature and pressure. A coding error in Version 1.0 was identified while testing the implementation of the viscosity model. This error is described in Appendix C. It mainly affects calculation of “higher order” properties such as compressibility and heat capacity near the critical point.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cross section measurements of deuteron electro-disintegration at very high recoil momenta and large 4-momentum transfers (Q^2)

The 2H(e,e'p)n cross sections have been measured at negative 4-momentum transfers of Q^2 = 4.5 +/- 0.5 (GeV/c)^2 and Q^2 = 3.5 +/- 0.5 (GeV/c)^2 reaching neutron recoil (missing) momenta up to p_r ?1.0 GeV/c. The data have been obtained at fixed neutron recoil angles 5 <= theta_nq <= 95 degrees with respect to the 3-momentum transfer q. The new data agree well with the previous data which reached p_r ?550 MeV/c. At theta_nq = 35 and 45 degrees, final state interactions (FSI), meson exchange currents (MEC) and isobar configurations (IC) are suppressed and the plane wave impulse approximation (PWIA) provides the dominant cross section contribution. The new data are compared to recent theoretical calculations, and a significant disagreement for recoil momenta p_r > 700 MeV/c is observed. The experiment was carried out in experimental Hall C at the Thomas Jefferson National Accelerator Facility (TJNAF) and formed part of a group of four experiments that were used to commission the new Super High Momentum Spectrometer (SHMS). The experiment consisted of a 10.6 GeV electron beam incident on a liquid deuterium target which resulted in the break-up of the deuteron into a proton and neutron. The scattered electrons were detected by the SHMS in coincidence with the knocked-out protons detected in the previously existing High Momentum Spectrometer (HMS) and the recoiling neutrons were reconstructed from energy-momentum conservation laws. To ensure that the 2H(e,e'p)n reaction channel was selected, we required the missing energy of the system to be the binding energy of the deuteron (?2.22 MeV). The spectrometers? central angles and momenta were set to measure three central missing momentum settings of the neutron corresponding to p_r = 80, 580 and 750 MeV/c, which required the SHMS central angle and momentum to be fixed and the HMS to be rotated from smaller to larger angles corresponding to the lower and higher missing momentum settings, respectively. The experiment was carried out in a time period of six days with typical electron beam currents of 45-60 uA at about 50% beam efficiency.

Yero, Carlos↗

Development of Control System Functional Capabilities within the IES Plug-and-Play Simulation Environment

The concept of an Integrated Energy System (IES) is meant to combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework is used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES. A new RAVEN (Risk Analysis Virtual ENvironment) plugin for grid and capacity optimization (HERON) has been developed for optimizing the production variables of the IES given the demand profile. Currently, only the limits that affect the production variables and their corresponding time rates of change are considered (explicit constraints). However, other variables are additionally subject to constraints, but the associated limits are not accounted for (implicit constraints). In particular, for the power dispatch problem, the optimization algorithm takes into account the limits on the electrical power output and the corresponding hourly power variations but does not consider other constraints on process variables whose response effects the service life of the IES. This report describes a scheme that allows accounting for implicit constraints without increasing the size of the optimization problem. The Reference Governor (RG) algorithm is traditionally used for enforcing state and control constraints by modifying the set-point trajectories supplied to the feedback regulators. In our application, the RG is coupled within an iterative loop with the HERON-power dispatcher to generate optimal trajectories that ensure the operational constraints are met. A data-driven procedure to derive a representation of the dynamics of the controlled system was developed. First, the variables that represented the state of the system are selected (PCA-based approach), and then state-space representation matrices are derived from the collected measurements (DMDc algorithm). A preliminary version of the developed workflow based on Linear Time Invariant matrices was assessed by adopting a two-unit test case. More sophisticated versions of this workflow foreseeing the on-line derivation of system matrices will be deployed in FY 2021. Finally, a “plug-and-play” library of controllers and state observers was developed in Dymola. Some aspects of the current configuration of the IES unit components, e.g., the encapsulation of the control schemes into dedicated blocks, are consistent with the “plug-and-play” philosophy. Other features, e.g., the system buses collecting the input and the output variables, are not. For this reason, once listed and described the limits of the current configuration, necessary modifications to the plant model interface are presented. As a test-case, the interfaces of the SES model in the RAVEN-based HYBRID framework were reworked accordingly, and two different control schemes were applied to the same plant model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Oxy-Combustion System Process Optimization

The overall objective of this work is to develop a new chemical absorbent-based, high pressure, CO 2 purification system to remove the residual oxygen that currently contaminates the recovered CO 2 , and to optimize the Pressurized Oxy-Combustion (POxC) process to minimize the Cost of Electricity (COE) generated in this advanced combustion process. TDA developed and validated the performance of the oxygen removal system for CO 2 purification. In collaboration with the Advanced Power & Energy Program (APEP) of University of California, Irvine (UCI), we optimized the POxC process, including thermal management, heat integration, and power cycle optimization using process design and modeling supported with Aspen Plus® process simulations. The techno-economic analysis results indicate that the pressurized oxycombustion coal power plant with Ion-transport membrane (ITM) air separation unit (ASU) (Case 2 – 30.55%) does not show an advantage over a cryogenic ASU (Case 1 – 31.24%) while TDA’s sorption-based ASU (Case 4 - 32.61%) shows a significant advantage over the cryogenic ASU (Case 1 – 331.24%). The specific plant costs show a wide range with a low of $2544/kW for Case 11C (TDA ASU, co-sequestering the SO x , and ultra-supercritical steam cycle) to a high of $2975/kW for Case 12A (cryogenic ASU and sCO 2 cycle). In general, the ITM ASU based cases have lower specific plant costs than the corresponding cryogenic ASU based cases while the TDA ASU based cases show the lowest specific plant costs. The main reason for lowering these costs is the higher overall plant thermal efficiency which decreases the plant cost on a $ per kW basis. Next comparing the cases with different power cycle working fluid conditions in terms of temperature and pressure while all utilizing steam, similar trends as the plant costs may be observed. However, with the supercritical CO 2 (sCO 2 ) cycle, the increase in thermal efficiency of the sCO 2 cycle was not able to offset its increase in plant cost making the plant costs higher than those of the corresponding steam cycle cases. The Cost of Electricity (COE) again shows similar trends as the specific plant costs. The COE for Case 11C at $110.1/MWh is also the lowest, but among all cases that do not co-sequester the SO x , Case 12C (TDA ASU and sCO 2 cycle) has the lowest COE at $\$$117.5/MWh while the highest is for Case 8A (cryogenic ASU and supercritical steam cycle) at $130.4/MWh.

20 FOSSIL-FUELED POWER PLANTS↗

Search for resonance in inclusive and b-tagged dijet mass spectra in proton-proton collision at $\sqrt{s}$ = 13 TeV and constraints on dark matter and other models

Two analyses searching for new heavy resonances and interpreting the result with dark matter are presented in this dissertation. The first analysis searches for resonances decaying into dijet final states using the proton-proton collision data at center mass of $\sqrt{s}=13$ TeV with the CMS detector at Large Hadron Collider in2016, corresponding to a total integrated luminosity of 36 fb -1 . No significant evidence for the production of new particles is observed. A simplified model of interactions between quarks and dark matter is used to interpret the result. The second analysis searches for resonances decaying into dijet final states using the proton-proton collision data at center mass of $\sqrt{s}=13$ TeV with the CMS detector at Large Hadron Collider in 2016, 2017, and 2018, corresponding to a total integrated luminosity of 137.5 fb -1 . A b-tagger named DeepJet is used to examine the events with jets containing B hadrons. No evidence of new heavy resonances is observed, as well. The 95% confidence level cross-section upper limit is set on excited b-quark b* and Coloron models. Limits are also set on a simplified model of interactions between b-quarks and dark matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Validation and Demonstration of Control System Functional Capabilities within the IES Plug-and-Play Simulation Environment

The concept of an integrated energy system (IES) is meant to combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework is used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES. The new RAVEN plugin for grid and capacity optimization (HERON) only addresses the limits that affect the production variables and the corresponding rates of variation (explicit constraints). However, other variables are subject to constraints, and the associated limits should be accounted for (implicit constraints). In particular, for the power dispatch problem, the optimization algorithm takes into account the limits on the electrical power output and the corresponding hourly power variations but does not consider other constraints on process variables whose response affects the service life of the IES. This report describes a scheme that allows accounting for implicit constraints without increasing the size of the optimization problem. To obtain a more accurate approximation of the nonlinear dynamic behavior, a parametric version of the dynamic mode decomposition with control (DMDc) algorithm was developed to derive the state-space representation matrices of the IES components at different scheduling parameter. Thanks to this approach, a more accurate approximation of the system response can be obtained, the limits imposed by thermal mechanical implicit constraints can be translated into power dispatch limits, and the feedbacks to HERON power dispatcher can be provided. To assess the developed methodology, a power dispatching test case composed of three power generating and storage units (Balance of Plant, Secondary Energy Source, Thermal Energy Storage) was developed. The power output of each one of the three units was optimized to meet the imposed time-dependent load demand trajectory and to maximize the IES profitability by meeting both the explicit and implicit constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Evaluation of Two-Dimensional to One-Dimensional Site Response for Idaho National Laboratory

We perform two-dimensional (2D) site response analyses accounting for spatial variability of soil properties and subsurface geometry of the Eastern Snake River Plane (ESRP), and quantify their effects on ground surface motion relative to one-dimensional (1D) site response analyses at the Idaho National Laboratory (INL). We first present the development of random field idealizations of the repeated basalt lava flows, heterogeneously inter-layered with sediments, from seismic velocity data collected over four decades in the ESRP. Using realizations of the stochastic fields mapped on 2D deterministic finite element models, we perform 2D viscoelastic and equivalent-linear wave propagation simulations, and quantify the mean and variance of site response aggravation factors, defined as the response spectral ratio of 2D to 1D analyses on the ground surface. Results are shown to be insensitive to the constitutive material behavior considered here, for strains induced by rock outcrop peak ground acceleration (PGA) as high as 0.7g: viscoelastic and equivalent-linear analyses predict peak mean 2D/1D aggravation factor 1.05 at period T=0.075 sec (i.e. the 2D response spectrum is 5% higher than the corresponding 1D at that period, on average), which corresponds to the wavelength of the horizontal correlation length of the random field (50m). For periods longer than the fundamental period of the site (here, T 1 =0.3125 sec), the propagating wavelengths are too long to be affected by the 1D site response and the aggravation factor becomes equal to 1. The standard deviation of the natural logarithms of the 2D/1D aggravation factors is ~0.15 for periods shorter than the fundamental period of the site, and decays thereafter at a steady rate.

58 GEOSCIENCES↗

Investigation of Thermolytic Hydrogen Generation Rate in Tank 44 Dissolved Saltcake Samples

Saltcake core samples collected from Tank 44 in 2006 were dissolved to provide material for HGR measurements applicable to F-Area dissolved saltcake material. Additionally, characterization was performed on the Tank 44 saltcake material. The following are key results from the Tank 44 saltcake characterization. The Tank 44 Upper Saltcake Composite, corresponding to the 171 to 285 inch tank level, contained by mass approximately 69% sodium nitrate, 11% sodium carbonate, 8% sodium nitrite, smaller amounts of other salts and components, and 9% unquantified (which includes water, water of hydration, oxygen/hydrogen content of oxides and hydroxides, and uncertainty). The Tank 44 Lower Saltcake Composite, corresponding to the 76 to 114 inch tank level, contained by mass approximately 49% sodium carbonate, 18% sodium nitrate, smaller amounts of other salts, at least 8% sludge, and 9% unquantified (see above). The dissolved saltcake contained free hydroxide less than quantifiable (<0.01 M) due to the limited quantity of material that could be removed from the Shielded Cells based on the sample radioactivity. Measurement by pH paper provided an approximate pH of 12. The following are key results from the Tank 44 HGR testing. During boiling at 106.7 °C, HGR for Tank 44 dissolved saltcake without added glycolate was 7.2×10 -8 ft 3 h -1 gal -1 . During boiling at 106.9 °C, HGR for Tank 44 dissolved saltcake with 1000 mg/L of added glycolate was 8.2×10 -8 ft 3 h -1 gal -1 . For the test without added glycolate, the first several HGR measurements at 70, 85, and 100 °C gave indication of the release of dissolved hydrogen and should not be used to represent the sustained thermolytic HGR for those temperatures. The measurements at boiling are the best representation of thermolysis in this testing. Carbon dioxide was observed at concentrations up to 6 vol% in the flow-system offgas for the test at boiling. Methane generation was observed at 100 °C and boiling. Methane concentration in the total gas generated during testing remained well below the lower flammability limit for methane in air. The addition of 1000 mg/L of glycolate did not have a significant impact on the hydrogen generation rates measured during this testing. The low hydroxide concentration in the Tank 44 dissolved saltcake likely influenced the relatively low thermolytic HGR and high carbon dioxide release observations in this testing. Based on the observation that methane was generated or released upon heating SRS radioactive Tank 44 waste samples to 100 °C and above, we recommend gaining a greater understanding of the cause and mechanism of its generation. First, the applicable literature should be reviewed to reveal the thermolytic methane generation mechanisms of possible methane generating species in the SRS CSTF. If warranted, a plan should be developed for simulant tests with methylated siloxanes and other applicable compounds in order to gain a better mechanistic understanding of methane generation in the SRS CSTF.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An Overview of Argonne’s Advanced Mobility Technology Laboratory Vehicle Systems Instrumentation and Evaluation Methodology

This report will provide a general overview of the testing facilities, research equipment, and general testing methodologies Argonne utilizes to conduct vehicle technology evaluations of advanced technology research vehicles. Data is captured from vehicle evaluations that provide powertrain operation and corresponding energy consumption based on a combination of in-depth instrumentation and focused testing. This resulting dataset of hundreds of time-resolved vehicle signals provide a basis for direct analysis and model validation of vehicle specific technologies. Argonne has attempted to standardize the approach to vehicle instrumentation and testing at Argonne, but it should be noted that each vehicle, and the corresponding assessment, remains unique. It is suggested that the reader reference vehicle specific testing reports for an overview of the unique aspects for each test vehicle. Additionally, the datasets for research vehicles are made publicly available through the Advanced Mobility Technology Laboratory’s Downloadable Driving Database (D3) at www.anl.gov/d3.

33 ADVANCED PROPULSION SYSTEMS↗

Structural Health Monitoring of Microreactor Safety Systems Using Convolutional Neural Networks

Microreactors, a class of modular reactors with net power output of less than 20 MWth, have innovative applications in nuclear and nonnuclear industries due to their portability, reliability, resilience, and high capacity factors. In order to operate microreactors on a wider scale, it is essential to bring down maintenance life-cycle costs while ensuring the integrity of operating such systems. Autonomous operations in microreactors using augmented digital-twin (DT) technology can serve as a cost-effective solution by increasing awareness about the system’s health. Structural health monitoring (SHM) is a key component of nuclear DT frameworks. Artificial neural networks can be beneficial to detect degradation in the nuclear safety systems, such as piping equipment systems, by monitoring the sensor data obtained from the plant and its corresponding structures, systems and components. In this report, an SHM methodology is presented which uses convolutional neural networks to determine degraded locations and their corresponding degradation-severity levels at various locations of nuclear piping equipment systems. A simple pipe system, subjected to seismic loads, is selected to design the post-hazard SHM framework. The effectiveness of the proposed SHM methodology is demonstrated by obtaining high accuracy in detecting degraded locations as well as the severity levels.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Graduate Fellowships [Slides]

All fellowship awards and their corresponding internships must be approved by LANL management. Awards and corresponding internship assignments are not guaranteed, even if student has management support.

99 GENERAL AND MISCELLANEOUS↗

Evaluation of Electrical Resistivity Tomography to Monitor the Transport of Past Releases Beneath Tank Farms

Underground storage tanks at the Hanford Site, in southeastern Washington State, hold radioactive waste generated from four decades of plutonium production. The 149 single-shell tanks and the 28 double-shell tanks have all exceeded their initial design life of approximately 25 years. At least 67 tanks are assumed to have leaked in the past, resulting in radioactive releases into the vadose zone. Gamma ray logging within dry monitoring wells is currently the primary method for tracking the migration of leaked tank waste through the vadose zone. While this approach provides an accurate assessment of radioactive contamination, that information is only provided near (within ~1m) the borehole, leaving most of the vadose zone unmonitored, particularly the important region directly beneath the tank. This report describes a numerical study that investigates the feasibility and performance of time-lapse 3D electrical resistivity tomography (ERT) for long-term monitoring of a hypothetical tank waste location and migration through the vadose zone. ERT is a method of remotely imaging the bulk electrical conductivity (EC) of the subsurface, which is significantly impacted by the presence of conductive solid and liquid tank waste. The release of liquid tank wastes increases subsurface fluid conductivity and saturation over time, creating a target to use time-lapse ERT for long-term monitoring. Although the presence of metallic infrastructure can cause ERT interference, recent advancements in ERT data processing enable the deleterious effects of buried metallic infrastructure (e.g. pipes, wellbore casings, tanks) to be removed to better determine the liquid tank waste migration over time. Three hypothetical realistic scenarios were simulated in the ERT evaluation. The first two scenarios assume the same leak amount and rate (i.e., between 1/1/1951 and 12/31/1951 at the rate of 347 m 3 per year) but different leaky tanks. Scenario 1 assumes leaks under tank B-102, which is located on the edge of the B-tank farm and surrounded by a few metallic infrastructure including cased pipes/wells/tanks. Scenario 2 assumes leaks under tank B-108, which is located near the center of the B-tank farm and surrounded by larger amount of metallic infrastructure than B-102. Scenario 3 assumes the same metallic infrastructure as B-102, with a more recent contaminant leak that was simulated to have occurred between 1/1/2018 and 12/31/2023 at a rate of 1.89 m 3 per year. The leak time in Scenarios 1 and 2 corresponds to a historical overfill event in 1951 and Scenario 3 corresponds to a recent found tank leak in 2019. In each scenario, a “true” bulk EC model vs. time reflecting contaminant migration was generated. ERT data was simulated from these “true” bulk EC models and a time-lapse ERT inversion produced “imaged” bulk EC vs. time. Three electrode configurations in two, four and eight boreholes surrounding the leak tank were used in the ERT simulations in each scenario. These borehole configurations were considered logistically feasible and cost-effective for monitoring. The hypothetical ERT boreholes are assumed to have non-metallic casing. By comparing the “imaged” bulk EC with the “true” bulk EC, it was demonstrated that the three configurations of wells used (two, four, and eight wells) were able to successfully monitor the migration of tank leaks through the vadose zone, with bulk EC resolution increasing with the number of down borehole ERT arrays for all scenarios. Therefore, the use of eight boreholes to perform ERT monitoring beneath the tanks provided the best spatiotemporal information.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Big Data For Operation and Maintenance Cost Reduction

The purpose of this research is to develop a first-of-a-kind framework for integrating Big Data capability into the daily activities of our current fleet of nuclear power plants. Big Data is traditionally defined as data sets with high volume, velocity, and heterogeneity, and the existing Big Data analytics capabilities are now widely popular in fields such as finance, weather, e-commerce, healthcare and sports. In the nuclear industry, while the volume and velocity of data may present computational challenges for existing analytics capabilities, data heterogeneity are seen to present the major challenge. This research project mainly focuses on incorporating the wide range of data heterogeneities in nuclear power plants into an integrated Big Data Analytics capability. The primary end-product of this project is a Big Data framework that is capable of dealing with the large volume and heterogeneity of the data found in nuclear power plants to extract timely and valuable information on equipment performance. The framework can generate system insights that are actionable relations between measurable impacts and the corresponding maintenance action plans and enable optimization of plant operation and maintenance based on the extracted information. The developed framework is capable of handling heterogeneous data including both image data and time-series sensor data. Specifically, this developed framework includes the following components. The first component is an overarching maintenance ontology which includes system insights required by maintenance optimization. The maintenance ontology interacts with other components in the developed framework. The second component handles Piping & Instrumentation Diagram (P&ID) data. It can be used to extract system components and their relations automatically from the P&IDs. This extracted information is stored in the first component, i.e., maintenance ontology, and is also used as input to the third component, i.e., a tool for generating the fault tree for the corresponding system. The generated fault tree in turn is stored in the ontology for assessing risk that is used as a criterion in maintenance policy optimization. The fourth component is a tool for inferring the parameters in the Markov degradation model for a nuclear system. It uses basic information from the ontology. The fifth component is a tool for assessing the degradation level using sensor measurement data, for example, pressure, flowrate. This tool can be used for determining corrective maintenance actions. The results obtained from components four and five are returned to the ontology. The sixth component of the framework is a tool for optimizing the maintenance policy for a nuclear system of interest. It takes certain basic information from the ontology, e.g., costs of maintenance actions and system failures, as input, and returns the optimal maintenance policy to the ontology. This tool can be used for determining predictive maintenance actions. A set of experiments have also been conducted to verify the algorithms developed in this project for nuclear system degradation monitoring. The experiments are based on four solenoid valves, similar to the ones used in nuclear power plants. The analyses based on the experimental data using two algorithms, i.e., the Randomized Window Decomposition (RWD) algorithm and the particle filtering algorithm, and the results are introduced in the report. The Big Data framework developed in this project can be used as a support tool in daily activities of plant operation and maintenance and will reduce current costs while maintaining or improving safety levels. Overall, the project will not only benefit existing reactors, however it will open new frontiers to realize the long overdue value of Big Data Analytics in the nuclear sphere.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Wood Decomposition: Understanding Processes Regulating Carbon Transfer to Soil Carbon Pools Using FACE Wood at Multiple Scales

Dead wood comprises approximately 20% of forest biomass. While it is widely recognized for its role in ecosystem functions relative to habitat and nutrient cycling, there are considerable uncertainties associated with its role in the forest carbon (C) cycle. The principal uncertainty is whether dead wood contributes to the soil C pool. With little supporting evidence, the presumption is that C in wood is respired back to the atmosphere through decomposition processes. The inability to attribute forest soil C to its source components has been a barrier to addressing this important forest carbon cycle question. The U.S. Dept. of Energy’s Free Air CO 2 Enrichment (FACE) experiment produced wood that has a distinct δ 13 C signature that can be used to trace its fate through organic matter decomposition processes. Using two hardwood species (aspen, birch) from the Rhinelander FACE site and pine from the Duke FACE site, the FACE Wood Decomposition Experiment (FWDE) was established in 2011 to provide the foundation for addressing wood decomposition and its associated processes at the continental scale. Plots were installed on nine experimental forests representing a range in bioclimatic forest conditions within the continental U.S. Each identical plot has logs of the three species incubating on the soil surface and pine logs suspended in the air, representing standing dead wood. Abiotic conditions are monitored continuously, and soil and forest conditions characterized. Logs on the sites have been measured every two years to assess decay. This project was initiated in 2016 with specific objectives to: (a) determine whether dead wood C contributes to the soil C pool, (b) assess the fungal community composition mediating the decay, (c) assess the interactions of termites with microbial decay, and (d) develop a mechanistic model to simulate wood decomposition. The principal field sampling was conducted in 2017; correspondingly, the results reported here represent six years of decomposition. Wood decomposition varied from 13 to 75 % loss across the bioclimatic gradient. Sites in the southeastern U.S. coastal plain and piedmont had fastest decomposition rates for each of the species. In contrast sites in the Rocky Mountains had the slowest decay rates, which were associated with cold temperatures. Standing deadwood decomposes at a much slower rate than logs on the soil surface; a response due to differences in wood moisture content. Wood carbon was measured in the mineral soil across each of the sites; accordingly, this is the first study to affirm that a portion of the wood C does become part of the soil C pool. Averaged across the FWDE sites, approximately 1.4% of the soil C in the upper 10 cm of the mineral soil directly beneath the log was derived from the deadwood; within 50 cm of the log, wood C accounted for an average of 2.1%. These results demonstrate that the contribution of deadwood to the soil C pool is a soil volume surrounding the log; and it’s likely a larger volume than what we sampled. We also hypothesized that there would be more wood C in the mineral soil when termites were present. However, that hypothesis wasn’t supported by measurements on the Santee Experimental Forest termite exclusion (SEF-TX) plots under FACE logs that incubated for 6 years. Across all the sites, the majority of the wood C was transferred to the atmosphere. Approximately 7,000 fungal species were sequenced. The fungal communities were different among each of the sites suggesting local colonization; interestingly, within sites the three species also exhibited differences in fungal community composition. White rot was the dominant fungal form in pine and hardwood, which yields less lignin in the wood residue. Brown rot occurred in approximately 20% of the pine and aspen logs, and it did not occur in birch logs. The occurrence of brown rot was only weakly correlated to site factors. We conducted two additional experiments to assess the effects of termites on wood decomposition and the associated interactions with the fungal community. Pine and birch logs were incubated on the Santee Experimental Forest termite exclusion plots in upland and riparian settings. There was no difference in decay between areas where termites were able colonize and plots where they had been excluded. Highly variable foraging by termites was a primary contributor the result. Correspondingly, fungal communities did not vary on logs incubated among the two treatments, but they did reflect site and species differences. The second experiment was conducted on each of the nine FWDE sites where pine and birch wood blocks were incubated for two years with and without a wrapping designed to exclude termites and other arthropods. The effects of termites on wood decomposition could be confirmed on two of the four sites known to have subterranean termites. On the two southeastern sites, termites contributed approximately 8% to mass loss in pine blocks over two years. A major factor clouding the role of termites in wood decomposition is the inability to assess the time, duration, and intensity of the foraging with the logs. Another uncertainty is related to whether a piece of wood will be foraged. A process-based wood decomposition model was developed representing biological processes mediated by fungi, bacteria, termites, beetles, as well as physical and chemical processes. The Coarse Wood Decomposition Model (CWDM) was parameterized based on a thorough review of the literature and it was validated using data from the FWDE. In the initial setup, the model considers softwoods and hardwoods groups, as opposed to individual species. Sensitivity analyses incorporating climate data from 89 locations in North America demonstrated that CWDM was stable from the tropics to the boreal zone. Model calibration and validation showed it to perform very well. The CWDM is poised for additional testing with subsequent FWDE data as well as other data sets, and it is available for incorporation into forest biogeochemical models or DOEs Earth System Modeling framework.

54 ENVIRONMENTAL SCIENCES↗

Automation of FARM from Alpha Phase to Beta Phase

Integrated energy systems (IES) combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework and the RAVEN plugin for grid and capacity optimization (HERON) are used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES, by respecting the limits on the production variables and the corresponding rates of variation (explicit constraints). Besides, there are other process variables whose evolution needs to be bounded to avoid damaging the components (e.g., condensers, heat exchangers, steam generators, etc.) or degrading the process efficiency (e.g., electrolysis in the hydrogen production process). To avoid violating these latter limits (implicit constraints), a proof-of-concept HERON validator based on Feasible Actuator Range Modifier (FARM-Alpha) was developed by Argonne National Laboratory in January 2021. This FARM-Alpha validator calculates the evolution of process variables on whom the implicit constraints are placed, and then provides feedback to HERON dispatcher to adjust the power setpoints of three IES components, i.e., Balance of Plant, Secondary Energy Source, and Thermal Energy Storage, so as to meet both the explicit and implicit constraints. FARM-Alpha was designed to assess the performance of FARM as a HERON validator only, i.e., the list of components and implicit operational constraints were hard-coded within the source code. The lack of flexibility of the corresponding software structure does not allow the deployment in production environment. This report describes the development and the implementation of an enhanced version of the FARM-based validator (FARM-Beta), which ensures more flexibility for the end user in modeling multiple IES configurations and scenarios. Several test cases of the power dispatch problem were then selected to demonstrate the capabilities offered by FARM-beta. The test cases illustrate the efficiency of the closed-loop optimization scheme and the capability to calculate set-point trajectories satisfying both explicit and implicit constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗