SEARCH · Search NASA

Results for “Maximum allowable concentrations”

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

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Quantifying the thermal effect and methyl radical production in nanosecond repetitively pulsed glow discharges applied to a methane-air flame

In this work, we investigated non-equilibrium plasma produced by nanosecond repetitively pulsed glow discharges applied across a lean premixed methane-air flame. The flame is stationary, axisymmetric, and laminar. The discharges are applied on the symmetry axis crossing the reactant gases, flame front, and product gases, allowing phase-locked averaged measurements and comparisons with axisymmetric numerical simulations. The thermal effect and methyl radical production are quantified in the discharge in the reactant gas region. One-dimensional, two-beam, hybrid, femtosecond-picosecond, coherent anti-Stokes Raman scattering is used to acquire spatial and temporal profiles of temperature and oxygen-to-nitrogen concentration ratio. Photo-fragmentation laser-induced fluorescence is used to acquire quantitative two-dimensional profiles of methyl radicals in the discharge providing the first quantitative imaging of methyl produced ahead of a flame by plasma-induced methane dissociation. The spatial profiles of temperature and oxygen-to-nitrogen concentration ratio are in steady state, indicating that individual discharges have an insignificant heating effect. Upper and lower bounds of the produced mole fraction of methyl radicals in the plasma are obtained due to uncertainties in the collisional quenching rates of excited state methylidyne radicals in the plasma. The discharges produce a maximum of 600–1100 ppm of methyl radicals upstream of the flame front within 25 ns. This amount is similar to the predicted methyl mole fraction for the flame without plasma and thus represents a significant chemical perturbation to the reactants upstream of the flame front. The produced methyl follows an exponential decay in the first microsecond after the discharge with a decay constant of 8 µs close to the flame, and 0.8 µs further from the flame. The decay then deviates from the exponential curve and the methyl persists for tens of microseconds. The results suggest that for the tested configuration, the thermal effect of individual discharges through fast gas heating is negligible, while active chemical species are produced in large quantities in the reactant gases, upstream of the flame front.

hybrid fs-ps CARS↗

DOE OSTI2022

Mechanically-, Thermally, and Chemically-Robust High-Temperature Ceramic Composites

In an effort to lower the levelized cost of electricity generated by Concentrated Solar Power (CSP) plants, the heat-to-electricity conversion efficiency may be increased by operating such plants with higher turbine inlet temperatures using closed Brayton cycles with high-pressure, supercritical carbon dioxide (sCO 2 ) as the working fluid. For example, if the turbine inlet temperature is raised to ≥750°C, in combination with sCO2 closed Brayton cycles, the relative heat-to-electricity efficiency may be raised by more than 20% (compared to subcritical steam Rankine cycles operating with turbine inlet temperatures ≤550°C). The associated reduction in the cost of dispatchable electricity from such CSP plants (utilizing high-temperature thermal energy storage) would be an important breakthrough on the path towards direct competition with fossil-fuel-based plants. A key barrier to achieve such high-temperature efficient operation has been the limited thermomechanical performance of metal alloys used in compact, primary printed circuit- type heat exchangers (PCHEXs) for heat transfer to high-pressure sCO 2 . The maximum allowed stresses for the use of conventional stainless steels and nickel-based superalloys at high sCO 2 pressures (≥20 MPa) decline rapidly at temperatures ≥550°C. This project has been focused on demonstrating the desired high-temperature properties and attractive manufacturing characteristics of mechanically-, thermally-, and chemically-robust ceramic/metal composites (cermets), in order to allow for the use of such composites in primary PCHEXs for heat transfer to sCO 2 at ≥750°C and ≥20 MPa in CSP plants. Prior work with high-temperature, co-continuous carbide/refractory metal cermets, such as ZrC/W composites, has already demonstrated: i) the attractive combination of high-temperature properties exhibited by such materials (e.g., high values of stiffness, failure strength, and thermal conductivity at 800°C), and ii) cost-effective processes for manufacturing PCHEX plates comprised of such cermets with tailorable channel patterns. However, ZrC/W cermets were not found to be inherently corrosion resistant in sCO 2 at 750°C and 20 MPa. While the application of a Cu coating to ZrC/W surfaces, along with a modest addition (50 ppm) of CO to CO 2 , rendered such composites resistant to corrosion for 1000 h in such sCO 2 -based fluids at 750°C, such additional steps would increase the design complexity of such cermet-based components. An alternative approach (the focus of this project) is to develop mechanically-robust, cost- effective cermets that are inherently resistant to high-temperature oxidation in air and CO 2 . The overall goals of this project have been to demonstrate: i) the oxidation resistance of at least one cermet at 750°C in air and in CO 2 (i.e., with a projected annual corrosion ≤30 μm), ii) the mechanically robust nature of at least one cermet at 750°C (i.e., with an average failure strength ≥200 MPa), and iii) that the selected cermet can be manufactured via low-cost forming techniques.

14 SOLAR ENERGY↗

DOE OSTI2020

High-Temperature Heat Transfer Fluid Circulator for Concentrating Solar Power Systems

A high temperature Heat Transfer Fluid circulator compressor with wide turn down ratio aerodynamic performance was designed, fabricated and tested for application with high pressure super Critical CO2 fluid capable of temperatures and pressures to >750°C and 25 MPa respectively. The modular system decoupled the drive motor system from an independent compressor shaft to provide thermal isolation between the high temperature compressor shaft and the lower temperature capable motor elements. While the thermal management approach appears to be viable there are potential losses in the secondary cooling flows that need to be taken into account and fine-tuned during operation in order to achieve desired overall system efficiencies. The system was also shown to be viable and dynamically stable throughout the entire operating speed range, thereby validating the overall system configuration and structural design approach used. Finally, while it was not possible to achieve specific PTC-10 test points, the compressor performance tests provide sufficient data to indicate that the compressor has the ability to achieve the required performance in practice. The two key factors affecting compressor design and cost are the maximum expected temperature and pressure. The maximum allowable temperatures for current state-of-the-art (SOA) motor driven compressor technology is limited to of approximately 150-200°C, but inlet temperatures >750 °C may be encountered. While the optimal efficiency of compression is in the transcritical region, some concepts require inlet pressures to 225 bar.

14 SOLAR ENERGY↗

DOE OSTI2021

Comparison of Deterministic and Statistical Models for Water Quality Compliance Forecasting in the San Joaquin River Basin, California

Model selection for water quality forecasting depends on many factors including analyst expertise and cost, stakeholder involvement and expected performance. Water quality forecasting in arid river basins is especially challenging given the importance of protecting beneficial uses in these environments and the livelihood of agricultural communities. In the agriculture-dominated San Joaquin River Basin of California, real-time salinity management (RTSM) is a state-sanctioned program that helps to maximize allowable salt export while protecting existing basin beneficial uses of water supply. The RTSM strategy supplants the federal total maximum daily load (TMDL) approach that could impose fines associated with exceedances of monthly and annual salt load allocations of up to $1 million per year based on average year hydrology and salt load export limits. The essential components of the current program include the establishment of telemetered sensor networks, a web-based information system for sharing data, a basin-scale salt load assimilative capacity forecasting model and institutional entities tasked with performing weekly forecasts of river salt assimilative capacity and scheduling west-side drainage export of salt loads. Web-based information portals have been developed to share model input data and salt assimilative capacity forecasts together with increasing stakeholder awareness and involvement in water quality resource management activities in the river basin. Two modeling approaches have been developed simultaneously. The first relies on a statistical analysis of the relationship between flow and salt concentration at three compliance monitoring sites and the use of these regression relationships for forecasting. The second salt load forecasting approach is a customized application of the Watershed Analysis Risk Management Framework (WARMF), a watershed water quality simulation model that has been configured to estimate daily river salt assimilative capacity and to provide decision support for real-time salinity management at the watershed level. Analysis of the results from both model-based forecasting approaches over a period of five years shows that the regression-based forecasting model, run daily Monday to Friday each week, provided marginally better performance. However, the regression-based forecasting model assumes the same general relationship between flow and salinity which breaks down during extreme weather events such as droughts when water allocation cutbacks among stakeholders are not evenly distributed across the basin. A recent test case shows the utility of both models in dealing with an exceedance event at one compliance monitoring site recently introduced in 2020.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Defect characterization and charge transport measurements in high-resolution Ni/n-4H-SiC Schottky barrier radiation detectors fabricated on 250 μ m epitaxial layers

Advances in the growth processes of 4H-SiC epitaxial layers have led to the continued expansion of epilayer thickness, allowing for the detection of more penetrative radioactive particles. In this work, we report the fabrication and characterization of high-resolution Schottky barrier radiation detectors on 250 μm thick n-type 4H-SiC epitaxial layers, the highest reported thickness to date. Several 8 × 8 mm 2 detectors were fabricated from a diced 100 mm diameter 4H-SiC epitaxial wafer grown on a conductive 4H-SiC substrate with a mean micropipe density of 0.11 cm -2 From the Mott–Schottky plots, the effective doping concentration was found to be in the range (0.95–1.85) × 10 14 cm -3 , implying that full depletion could be achieved at ~5.7 kV (0.5 MV/cm at the interface). The current-voltage characteristics demonstrated consistently low leakage current densities of 1–3 nA/cm 2 at a reverse bias of -800 V. This resulted in the pulse-height spectra generated using a 241Am alpha source (5486 keV) manifesting an energy resolution of less than 0.5% full width at half maximum (FWHM) for all the detectors at -200 V. The charge collection efficiencies (CCEs) were measured to be 98–99% with no discernable correlation to the energy resolution. A drift-diffusion model fit to the variation of CCE as a function of bias voltage, revealed a minority carrier diffusion length of ~10 μm. Deep level transient spectroscopy measurements on the best resolution detector revealed that the excellent performance was the result of having ultralow concentrations of the order of 1011 cm-3 lifetime limiting defects—Z 1/2 and EH 6/7 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Projecting Antarctica's contribution to future sea level rise from basal ice shelf melt using linear response functions of 16 ice sheet models (LARMIP-2)

Abstract. The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

High-Temperature Heat Transfer Fluid Circulator for Concentrating Solar Power Systems How Mohawk and Velo3Dare Printing the Future of Energy

14 SOLAR ENERGY↗

DOE OSTI2025

Numerical Investigation of Susceptor-Catalyst Design for Ethylene Generation in Radio Frequency Based Reactors

Ethylene is the most widely produced petrochemical component in the world. Whether reactors are heated directly or indirectly via steam, manufacturers use economies of scale to overcome inherent thermodynamic inefficiencies when burning fossil fuels. While allowing large-scale operations to use alternative sources of energy and raw materials, new methods of supplying energy to chemical reactor systems can reduce the energy waste produced by conventional processes. One viable method for effectively supplying energy to reactor systems is electromagnetic (EM) induction heating. Using the properties of radio frequency (RF) waves, heterogeneous catalysts can be precisely targeted for heating inside reactors. Site-selective heating can greatly lower the energy requirements of the process by supplying heat to reaction sites while reducing needless heat transfer elsewhere. In this study, a microscale model was used to help create guidelines for susceptors and catalysts to improve ethylene production. The oxidative dehydrogenation of ethane is investigated by utilizing several catalysts and potential catalyst/susceptor combinations, with heat provided by an EM susceptor. Having the susceptors and catalyst function separately results in higher gradients in both heat and mass transfer, which drives transport through the catalyst region, while still providing adequate heating for endothermic reactions. However, using a susceptive core with a catalyst covering produces the maximum ethylene concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Displacement Analysis of Geothermal Field Based on PSInSAR And SOM Clustering Algorithms A Case Study of Brady Field, Nevada—USA

The availability of free and high temporal resolution satellite data and advanced SAR techniques allows us to analyze ground displacement cost-effectively. Our aim was to properly define subsidence and uplift areas to delineate a geothermal field and perform time-series analysis to identify temporal trends. A Persistent Scatterer Interferometry (PSI) algorithm was used to estimate vertical displacement in the Brady geothermal field located in Nevada by analyzing 70 Sentinel-1A Synthetic-Aperture Radar (SAR) images, between January 2017 and December 2019. To classify zones affected by displacement, an unsupervised Self-Organizing Map (SOM) algorithm was applied to classify points based on their behavior in time, and those clusters were used to determine subsidence, uplift, and stable regions automatically. Finally, time-series analysis was applied to the clustered data to understand the inflection dates. The maximum subsidence is –19 mm/yr with an average value of –6 mm/yr within the geothermal field. The maximum uplift is 14 mm/yr with an average value of 4 mm/yr within the geothermal field. The uplift occurred on the NE of the field, where the injection wells are located. On the other hand, subsidence is concentrated on the SW of the field where the production wells are located. The coupling of the PSInSAR and the SOM algorithms was shown to be effective in analyzing the direction and pattern of the displacements observed in the field.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Fabrication of Advanced Nanocarbon-Metal Composites for Improved Energy Efficiency

We fabricated nanocarbon metal composites (NCMC) of Al alloys with a process called “electrocharging assisted process” (EAP). This method consists of the application of a high current to a mixture of liquid metal and carbon particles. We investigated aluminum alloys (6061 and 1350) and used activated carbon and graphite powder as the source of carbon. The high current induces the formation of carbon chains and ribbons in the liquid metal. Upon solidification of the metal an epitaxial relation between the carbon nanostructures and the metal lattice is produced. The purpose of the project was to find the parameters during the reaction that would give rise to a high density of nanoribbons extending throughout the metal such that the electrical conductivity and the mechanical strength of the composite increased and that the method could be extended to large scale manufacturing. For this purpose, we designed two reactors for the incorporation of nanocarbon ribbons in aluminum metal. The first reactor was designed to have more control of the region with the high current density. However, there were problems with getting good mixing of the carbon in small volumes. A second reactor was designed with a stirrer that allowed for better mixing of the carbon by the introduction of argon gas through the shaft of the stirrer. NCMC were fabricated with a series of parameters to understand the role of current, time of applied current, type and shape of the cathode electrode and stirring speed of the mixture. We analyzed the composites by Raman scattering to gain information on the crystallite size of the nanocarbon, XRD to obtain the crystal structure of the composite, SEM and TEM to characterize the grain size of the aluminum grains and the quality of the crystal structure. We also measured the electrical conductivity and mechanical properties of selected samples. The nanocomposites showed increase in crystallite size of the nanocarbon with a linear dependence of the crystallite size on the duration of the applied current. There is a minimum current density necessary for the crystallite size to increase compared to the ~ 10 nm size of the activated carbon source used in the fabrication. The electrical conductivity of the nanocomposites increased with crystallite size of the nanocarbon and with the concentration of converted carbon. The maximum increase in electrical conductivity was 5.7% above the baseline for samples with ~ 4 wt % nanocarbon with crystallite size larger than 30 nm. These samples presented a hardness ~ 8% higher than the baseline samples with no carbon.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

Hybrid HEFA-HDCJ Process for the Production of Jet Fuel Blendstocks

The hydrotreatment of bio-oil derived from the pyrolysis and biocrude from hydrothermal liquefaction of lignocellulosic materials to produce hydrocarbons faces significant technological challenges, mainly due to the high reactivity and poor thermal stability of bio-oil, resulting in the formation of large quantities of coke. This problem has been addressed by existing PNNL patents with a two-step hydrotreatment technology in which the bio-oil is first stabilized with a noble hydrogenation metal (often Pt or Ru). Then, in the second step, the bio-oil is deoxygenated with a Ni-Mo or Co-Mo sulfide catalyst. The main problem with this approach is that the Pt/Ru catalysts deactivate easily in the presence of S or other impurities, which are commonly present in pyrolysis oils. In this project, we explored technological solutions to mitigate coke formation, avoiding the use of Pt/Ru catalysts. Our strategy is based on three actions: (1) Bio-oil stabilization in the presence of alcohols. In this project, we studied the stabilization with butanol. (2) the use of a cosolvent to solubilize the bio-oil. Because coke formation reactions are second-order reactions a reduction in the concentration of reactive bio-oil molecules. In this case, we used yellow greases as a co-solvent. (3) Separation of bio-oil reactive fractions. In this project, we studied the removal of water-soluble fractions. Our batch co-hydrotreatment studies confirmed that the addition of butanol and methanol and the blend with lipids effectively contributed to mitigating coke formation (reducing coke yield to about 1 wt.% %). The removal of sugars did not have a noticeable effect on the overall coke yield, suggesting that coke precursors are present in all bio-oil fractions. Our analytical work suggests that they may be concentrated in the water-insoluble/CH 2 Cl 2 insoluble fractions of pyrolysis oils. Although the technological strategies tested resulted in significant coke reductions, the levels achieved were not sufficiently low to ensure a reliable operation in continuous, fixed-bed trickle-bed reactors. Long runs of more than 100 hours (maximum: 255 h) of co-processing time on stream were achieved in a continuous 40 mL reactor. When the same test was conducted in a larger 400 mL reactor, pressure drop increases associated with coke formation were observed. This increase in coke formation could be due to larger temperature gradients in the bed. Hydrodeoxygenation tests in moving bed reactors and using more active hydrogenation catalysts (for example Ni) could lead to more reliable operations. Unfortunately, our team did not have access to such experimental setups. The technoeconomic analysis suggests that although alcohol use is an effective means to reduce coke formation, the use of alcohol increases production cost. Thus, its use needs to be minimized. A delicate balance needs to be found between the use of technological solutions that allow the reliable operation of the system (stabilization with Ni catalysts, use of small quantities of solvents, processing in moving bed reactors) with a tolerable level of coke formation for the hydrodeoxygenation reactor used and that result in minimum production costs.

09 BIOMASS FUELS↗

DOE OSTI2021

Predicting Thermal Excursions During in situ Oxidative Regeneration of Packed Bed Catalytic Fast Pyrolysis Catalyst

Ex situ catalytic fast pyrolysis (CFP) uses a secondary reactor to upgrade biomass pyrolysis vapors to stabilized CFP oils with reduced oxygen content. In one configuration, the secondary reactor is operated as a packed-bed swing reactor system which allows coke-deactivated beds to be decarbonized in situ while other beds remain online for vapor upgrading. In situ decarbonization must be done carefully to avoid irreversible deactivation and/or physical degradation of catalyst pellets. Given that packed bed reactors are well known to have poor heat transfer characteristics, this is a critical issue impacting scaleability and commercial viability of the technology. To predict thermal excursions during regeneration, finite element computational models have been built to assist in scaling up oxidative decarbonization of a Pt/TiO2 CFP catalyst (0.5 mm spheres) from a bench scale packed bed with 100 g of catalyst to a pilot scale packed bed with 2 kg of catalyst and internal cooling tubes. Based on transient measurements of outlet temperature and effluent CO2 concentration, and using an assumed coke profile and activation energy, this paper demonstrates that specific combinations of effective thermal conductivity and wall heat transfer coefficient can fit bench scale oxidative regeneration data equally well. For the upscaled 2 kg bed, four bench-scale “best fit” parameter pairs give different predictions for location and magnitude of thermal excursions, with the maximum computed bed temperature gradients ranging from 30 °C cm-1 to as high as 3000 °C cm-1. The larger the fraction of heat removal by conduction through the cooling tubes, the greater the differences between the parameter pairs. The modelling results presented in this paper cast doubt on the industrial viability of the proposed combination of catalyst, bed and regeneration process, and point to the need for alternate reactor designs. However, there is considerable uncertainly in some of the key model parameters. The reliability of model predictions can be increased by adding more temperature measurements at key bed locations, testing additional variations in process conditions, performing careful bed dissections to determine the true coke profile, and perhaps most importantly, directly measuring the effective thermal conductivity of the catalyst pellets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

A validated thermal and biological model for predicting algal productivity in large scale outdoor cultivation systems

Proper assessment of the sustainability of algal products requires an understanding of algal growth rate with regional and temporal resolution. This study developed a bulk growth model that utilizes geospatial data, reactor geometry inputs, and a maximum of six strain-specific parameters to predict algal productivity with temporal and regional resolution. The model requires local weather data to calculate the time-resolved culture temperature which is combined with a biological model based on maximum theoretical carbon fixation to predict the time-resolved algal concentration. Validation of the model against experimental data illustrates an average ac- curacy of 0.27% ± 5.32% when comparing algal concentration at harvest over 45 days of cultivation data for Galdieria sulphuraria 5587.1, and Galdieria sulphuraria Soos grown in photobioreactors (PBRs) in Mesa, AZ (33.4152° N, 111.8315° W). For open raceway ponds (ORPs), an accuracy of 0.9% ± 2.35% was achieved when comparing algal concentration at harvest over 50 days of cultivation data for Chlorella vulgaris, Desmodesmus intermedius, Nannochloropsis oceanica grown in Atlanta, GA (33.7490° N, 84.3880° W) and Mesa, AZ (33.4152° N, 111.8315° W). The validation trials spanned three seasons (winter, spring, and summer) to ensure model ac- curacy throughout the year, and concentration between harvests was modeled with hourly resolution. The validated model is used to assess the productivity at several case-study locations in the United States with Chlorella vulgaris, grown in open raceway ponds and Galdieria sulphuraria grown in vertical flat panel photo- bioreactors. In conclusion, the model has been published as an open-source tool, in an effort to increase the fidelity of future studies that rely on outdoor algae cultivation for feedstock production and allow for time resolved results for future techno-economic and life cycle assessments.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Reduced Temperature Cesium Removal from AP-101 Using Crystalline Silicotitanate

The Tank Side Cesium Removal (TSCR) system, currently operational under Washington River Protection Solutions LLC (WRPS), sends initial low-activity Hanford waste tank supernate feed to the Hanford Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Facility. In addition to entrained solids removal from the supernate, the primary goal of TSCR is to remove cesium-137 ( 137 Cs) by ion exchange, allowing contact handling of the liquid effluent product at the WTP. Crystalline silicotitanate (CST) ion exchange media, manufactured by Honeywell UOP, LLC (product IONSIV™ R9140-B), was selected as the ion exchange media at TSCR. Laboratory-scale ion exchange processing using TSCR prototypic unit operations continues to contribute toward WRPS establishing accurate process flowsheets for the individual feed campaigns planned for TSCR. This report describes the small-scale ion exchange testing with 14.0 L of diluted and filtered supernate from tank 241-AP-101 (AP-101DF) at 16 °C (62 °F) to demonstrate processing at temperature conditions that are more prototypic of what the TSCR system could experience during colder seasons of the year. Since CST Cs capacity increases with decreasing contact temperature, testing at the lower operating temperature will help to predict the maximum 137 Cs loading onto the CST in the TSCR system. One of the waste acceptance criteria (WAC) for the WTP Low-Activity Waste Facility is that the waste must contain less than 3.18×10 -5 Ci 137 Cs per mole of Na. For the AP-101DF tank waste to meet this criterion, only 0.144% of the influent 137 Cs concentration may be delivered to the WTP; this requires a Cs decontamination factor of 694. Testing with AP-101DF matched TSCR prototypic operations where a lead-lag configuration was used until the lag column reached the WAC limit, then a polish column was brought online for continued processing in a lead-lag-polish column configuration. Feed was processed at 1.9 bed volumes (BVs) per hour; the flowrate, in terms of contact time with the CST bed, matched the expected flowrate at TSCR. The Cs-decontaminated product was retained for vitrification testing (to be reported separately). The lead column reached 62% Cs breakthrough after processing ~1400 BVs of feed; the 50% Cs breakthrough occurred at 1250 BVs. Testing compared to previous AP-107 testing at 16 °C showed ~80 BV increases in volume processed to reach the WAC limit for both lead and lag columns. A similar slope in breakthrough curves for both tests indicates similar kinetic behavior, with variations in feed matrices (Na and Cs concentrations) likely responsible for the deviations in reaching the WAC limit. The Cs effluent from the lag column reached the WAC limit after processing 875 BVs. Anticipating this breakthrough point, the polish column was preemptively installed at 770 BVs. Cs breakthrough from the lag column began at 300 BVs, reaching 5.32×10 0 µCi/mL, or 5.6 % Cs breakthrough, after processing all 1400 BVs of feed. Table S.1 and Figure S.1 summarize the observed column performance and relevant Cs loading characteristics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Simulations of short-term exposure to NO 2 and PM 2.5 to inform capture efficiency standards

The California building code requires all new or renovated residential dwelling units to have kitchen exhaust ventilation to manage air pollutants and moisture generated during cooking. Current performance requirements are specified as a minimum airflow and maximum sound level. This report presents an analysis to support consideration of adding a capture efficiency requirement to the code. The analysis uses a physics-based simulation model to calculate air pollutant concentrations in homes that result from cooking and the inflow of outdoor air, mediated by loss and removal process including deposition, dwelling unit ventilation, and the use of a range hood with varying levels of capture efficiency. Calculated pollutant concentrations are compared to relevant, health-based guidelines. The analysis considers the highest 1-hour concentration of nitrogen dioxide (NO 2 ), which is a product of natural gas combustion and the highest 24-hour concentrations of fine particulate matter (PM 2.5 ) which is emitted in substantial quantities from frying, broiling and grilling among other cooking activities. For NO 2 , the analysis considers cooking of a dinner for 3-4 persons consisting of pasta, meat sauce, a par-boiled vegetable, and baked garlic bread. For PM 2.5 , the analysis considers a day in which breakfast emitted particles at the 80th percentile and lunch and dinner entailed particle emissions at the 50th percentile of cooking emission events reported in the literature. Model simulations were performed to determine the level of range hood capture efficiency that will allow these cooking scenarios to occur in the vast majority (>99%) of new homes being built in California while maintaining pollutant concentrations below the health-based guidelines, if the range hood is used throughout cooking. All homes were assumed to have dwelling unit ventilation at the rate required in the building code. Simulation model input parameters were specified using a Monte Carlo approach to represent a range of housing characteristics, outdoor conditions, and indoor pollutant dynamics. Simulation results suggest that requiring a minimum capture efficiency of at least 70% is needed to avoid unacceptably high NO 2 (1-h average concentration of 100 ppb or higher) and at least 60% to avoid unacceptably high PM 2.5 (24-h average of 25 μg/m³ or higher). These results were driven by multi-family homes, which have a smaller volume of air to dilute any pollutants not captured and removed at the cooking area, resulting in higher concentrations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Addressing Low-Cost Methane Sensor Calibration Shortcomings with Machine Learning

Quantifying methane emissions is essential for meeting near-term climate goals and is typically carried out using methane concentrations measured downwind of the source. One major source of methane that is important to observe and promptly remediate is fugitive emissions from oil and gas production sites but installing methane sensors at the thousands of sites within a production basin is expensive. In recent years, relatively inexpensive metal oxide sensors have been used to measure methane concentrations at production sites. Current methods used to calibrate metal oxide sensors have been shown to have significant shortcomings, resulting in limited confidence in methane concentrations generated by these sensors. To address this, we investigate using machine learning (ML) to generate a model that converts metal oxide sensor output to methane mixing ratios. To generate test data, two metal oxide sensors, TGS2600 and TGS2611, were collocated with a trace methane analyzer downwind of controlled methane releases. Over the duration of the measurements, the trace gas analyzer’s average methane mixing ratio was 2.40 ppm with a maximum of 147.6 ppm. The average calculated methane mixing ratios for the TGS2600 and TGS2611 using the ML algorithm were 2.42 ppm and 2.40 ppm, with maximum values of 117.5 ppm and 106.3 ppm, respectively. A comparison of histograms generated using the analyzer and metal oxide sensors mixing ratios shows overlap coefficients of 0.95 and 0.94 for the TGS2600 and TGS2611, respectively. Overall, our results showed there was a good agreement between the ML-derived metal oxide sensors’ mixing ratios and those generated using the more accurate trace gas analyzer. This suggests that the response of lower-cost sensors calibrated using ML could be used to generate mixing ratios with precision and accuracy comparable to higher priced trace methane analyzers. This would improve confidence in low-cost sensors’ response, reduce the cost of sensor deployment, and allow for timely and accurate tracking of methane emissions.

03 NATURAL GAS↗

DOE OSTI2024

Strain-concentration for fast, compact photonic modulation and non-volatile memory

A critical figure of merit (FoM) for electro-optic (EO) modulators is the transmission change per voltage, d T / d V . Conventional approaches in wave-guided modulators maximize d T / d V via a high EO coefficient or longer light-material interaction lengths but are ultimately limited by material losses and nonlinearities. Optical and RF resonances improve d T / d V at the cost of spectral non-uniformity, especially for high- Q optical cavity resonances. Here, we introduce an EO modulator based on piezo-strain-concentration of a photonic crystal cavity to address both trade-offs: (i) it eliminates the trade-off between d T / d V and waveguide loss—i.e., enhancement of the resonance tuning efficiency d v c / d V for the fixed EO coefficient, waveguide length, and cavity Q —and (ii) at high DC strains it exhibits a non-volatile (NV) cavity tuning Δ v c ,NV for passive memory and programming of multiple devices into resonance despite fabrication variations. The device is fabricated on a scalable silicon nitride-on-aluminum nitride platform. We measure d v c / d V =177±1MHz/V, corresponding to Δ v c =40±0.32GHz for a voltage spanning ±120V with an energy consumption of δ U /Δ v c =0.17nW/GHz. The modulation bandwidth is flat up to ω BW,3dB /2 π =3.2±0.07MHz for broadband DC-AC and 142±17MHz for resonant operation near a 2.8 GHz mechanical resonance. Optical extinction up to 25 dB is obtained via Fano-type interference. Strain-induced beam-buckling modes are programmable under a “read-write” protocol with a continuous, repeatable tuning range of 5±0.25GHz, allowing for storage and retrieval, which we quantify with mutual information of 2.4 bits and a maximum non-volatile excursion of 8 GHz. Using a full piezo-optical finite-element-model (FEM) we identify key design principles for optimizing strain-based modulators and chart a path towards achieving performance comparable to lithium niobate-based modulators and the study of high strain physics on-chip.

Wen, Y. Henry (ORCID:0009000685423628)↗

DOE OSTI2023

Reduced Order Model to Predict Dispersion of Flammable Refrigerant into a Space

As the HVAC&R industry mobilizes to deploy more low-GWP refrigerants, relevant standards are being continually reviewed and updated. Those include the general safety standards ISO 5149 and ASHRAE 15, and the equipment standards IEC and UL. The standards systematically set the allowable maximum amount of refrigerant that should be used in different equipment types and different applications. To do so, they rely on predictions of how a leaked refrigerant mass will disperse into a space. Dispersion characteristics, such as total flammable volume and its residence time, determine the risk associated with the presence of the flammable refrigerant. The standards have included provisions for the use of flammable refrigerants for approximately two decades. They relied on limited analytical analyses and test cases in their development. Dispersion of a refrigerant into a space is complex. Computational fluid dynamics (CFD) are the most accurate in predicting a given problem. However, CFD is computationally expensive and requires specialized expertise and resources and is not suitable for use by standards development working group as prediction tool. This paper presents the development of a reduced order model (ROM) that predicts the key dispersion characteristics relevant to the dispersion of a leaked refrigerant into a space for any combination of input variables. The inputs are the refrigerant release height, the total released refrigerant mass and its release flow rate, the refrigerant molecular weight, the ventilation flow rate, the floor area and height of the space, recirculation air flow rate, and the tightness of the space. The outputs are histograms of volume fraction of the room in prescribed concentration bins and the total mass of the refrigerant in each bin normalized by the total refrigerant charge at 13 prescribed simulation time stamps between 1 and 900 seconds. The ROM is constructed from a set of CFD simulations with carefully chosen combinations of input parameters. The selection if done using a multidimensional sparse grid which is a generalization of the classical tensor approach but offers additional flexibility and thus can be more carefully tuned towards a specific model. The tuning is done to improve the accuracy, measured in the difference between the output values of the ROM and the CFD model, while minimizing the computational cost, measured in number of CFD simulations which is orders of magnitude more expensive than the processing the training data.

Edwards, Dean↗