Search NASA⌕ Search

SEARCH · Search NASA

Results for “Using”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Water use and radiation balance of miscanthus and corn on marginal land in the coastal plain region of North Carolina

Abstract Miscanthus is a perennial grass that can yield substantial amounts of biomass in land areas considered marginal. In the Coastal Plain region of North Carolina, marginal lands are typically located in coarse‐textured soils with low nutrient retention and water‐holding capacity, and high erosivity potential. Little is known about miscanthus water use under these conditions. We conducted a study to better understand the efficiency with which miscanthus uses natural resources such as water and radiant energy to produce harvestable dry biomass in comparison to corn, a typical commodity crop grown in the region. We hypothesized that under non‐limiting soil water conditions, miscanthus would have greater available energy and water use rates owing to its greater leaf area, thus leading to greater agronomic yields. Conversely, these effects would be negated under drought conditions. Our measurements showed that miscanthus intercepted more radiant energy than corn, which led to greater albedo (by 0.05), lower net radiation (by 4% or 0.4 MJ m −2 day −1 ), and lower soil heat flux (by 69% or 1.0 MJ m −2 day −1 ) than corn on average. Consequently, miscanthus had greater available energy (by 7% or 0.6 MJ m −2 day −1 ) and water use rates (by 14% or 0.5 mm day −1 ) than corn throughout the growing season on average, which partially confirmed our hypothesis. Greater water use rates and radiation interception by miscanthus did not translate to greater water‐use (1.5 g kg −1 vs. 1.6 g kg −1 ) and radiation‐use (0.9 g MJ −1 vs. 1.1 g MJ −1 ) efficiencies than corn. Compared to literature values, our data indicated that water and radiation availability were not limiting at our study site. Thus, it is likely that marginal land features present at the Coastal Plain region such as low soil fertility and high air temperatures throughout the growing season may constrain agronomic yields even if soil water and radiant energy are non‐limiting.

Carvalho, Henrique D. R.↗
DOE OSTI2026

Optimal Methods for Estimating Cactus Pear Biomass Using Cladode Dimensions of Morphologically Diverse Accessions

Current allometric methods for photosynthetic-stem (cladode) plants, such as cactus pear (Opuntia spp.), require refinement to be used in field settings in which diverse accessions are grown. We analysed cladode dimensional data using 14 accessions representing four species and two hybrids to quantify statistically significant morphological differences among accessions and derived cross-accession models to approximate cladode fresh weight. A Box model using cladode dimensions (e.g., length, width, thickness and diameter) and factorial combinations of these measures (e.g., length*width*thickness*diameter vs. fresh weight) resulted in the highest coefficient of determination (R 2 = 0.95 general fit) across all accessions for estimating fresh weight along with parsimony estimates using the Schwarz–Bayes Criterion (SBC), which assesses the most consistent performance on individual accessions. A Fitting-box modelling approach used the measured cladode area captured using ImageJ (R 2 = 0.93 general fit). Lastly, an Elliptical model used an elliptical approximation for the measured area and performed well over all accessions (R 2 = 0.94 general fit) while avoiding extensive manual measurements. These models meet or exceed the performance of previously published approaches when applied across morphologically diverse accessions, providing efficient tools for nondestructive estimation of cactus pear biomass under the conditions tested.

Opuntia↗
DOE OSTI2024

Recovery of Oxpure 612C-50 Coconut Carbon Particles Using Bump Arrays

Cyanide is used to leach gold from crushed ore (solid matrix) into a gold cyanide solution. The gold is extracted from the cyanide solution by adsorption onto activated carbon. The gold extraction process occurs when the activated coconut carbon is placed into tanks that contain the gold cyanide solution either in a batch process or into a continuous flow circuit. The coconut carbon is then removed from the solution for gold recovery. The use of the mesofluidic separation system (a deterministic lateral displacement system) can significantly reduce costs and waste generated from the process of removing the coconut carbon from the solution. In the coconut carbon recovery process used at many gold mines in Nevada, the gold is attached to carbon particles that are still in the cyanide liquor. The mesofluidic system can rapidly remove the gold bearing carbon particles from the cyanide liquor quickly and with no operating costs. The benefits include: • Reduce operational costs related to filtration and hydro-cyclones • Reduced acid usage for the elution gold stripping phase as only 25% of the fluids will follow the carbon particles to the final express lane when two mesofluidic systems are used in series This alternative particle removal techniques would be advantageous. A promising technique for removing larger particles from slurries is mesofluidic separation, similar to deterministic lateral displacement arrays or “bump” arrays [1] but operates at much larger flow rates. As described by Pease, et al. [2], “Bump arrays in deterministic lateral displacement devices separate large particles from small particles using arrays of staggered posts. Large particles, defined as those with radii larger than the distance between the edge of a post and the stagnation streamline from the next downstream post, must bump toward one side of the device, whereas particles smaller than this distance slalom from entrance to exit without net lateral displacement.” Unlike filters or sieves, the posts in the arrays that cause separation do not block or occlude particle flow but work because particles go around the posts. Unlike hydrocyclones separation, separation is not driven by particle density, because gravitational forces are unimportant to mesofluidic separation. Burns, et al., and Pease, et al., have shown that particles may be separated from complex suspensions, under turbulent flow conditions, and at industrially important flow rates [3-5]. They have also recently shown that mesofluidic devices may be arranged in series to increase separation performance [6]. In this paper, we evaluate the mesofluidic system for the separation of commercial OxPure GR 612 Charred Coconut shell particles as a proof of principle test for the rapid separation using the 1500 micron cut mesofluidic separator. We first describe the experimental system and conditions. We then present the experimental results. [1] Huang, L.R., E.C. Cox, R.H. Austin, J.C. Sturm. 2004. Continuous particle separation through deterministic lateral displacement. Science, 304, 987–990. https://doi.org/10.1126/science.1094567. [2] Pease L.F., J.A. Bamberger, C.A. Burns, and M.J. Minette. 2021. Large Particle Separation from Non-Newtonian Slurries using Bump Arrays. In Proceedings of the ASME 2021 Fluids Engineering Division Summer Meeting FEDSM2021-65904, V003T08A023. New York, New York: ASME. doi:10.1115/FEDSM2021-65904 [3] Burns C.A., T.G. Veldman, J. Serkowski, R.C. Daniel, X.-Y. Yu, M.J. Minette, L.F. Pease. 2021. Mesofluidic separation versus dead-end filtration. Separation and Purification Technology, 254, 117256. https://doi.org/10.1016/j.seppur.2020.117256. [4] Pease L.F., J.E. Serkowski, T.G. Veldman, J. Williams, X.-Y. Yu, M.J. Minette, J.A. Bamberger, C.A. Burns. 2021. Can Bump Arrays Separate Particles from Turbulent Flows?. In Proceedings of the ASME 2021 Fluids Engineering Division Summer Meeting FEDSM2021-67696, V003T08A024. New York, New York: ASME. doi:10.11

mesofluidic separation, slurry, bump array↗
DOE OSTI2025

Bayesian Inference for the Seismic Moment Tensor Using Regional Waveforms and Teleseismic- P Polarities with a Data-Derived Distribution of Velocity Models and Source Locations

The largest source of uncertainty in any source inversion is the velocity model used in the transfer function that relates observed ground motion to the seismic moment tensor. However, standard inverse procedure often does not quantify uncertainty in the seismic moment tensor due to error in the Green’s functions from uncertain event location and Earth structure. Here, we incorporate this uncertainty into an estimation of the seismic moment tensor using a data-derived distribution of velocity models based on complementary geophysical data sets, including thickness constraints, velocity profiles, gravity data, surface-wave group velocities, and regional body-wave travel times. The data-derived distribution of velocity models is then used as a prior distribution of Green’s functions for use in Bayesian inference of an unknown seismic moment tensor using regional and teleseismic-P waveforms. The use of multiple data sets is important for gaining resolution to different components of the moment tensor. The combined likelihood is estimated using data-specific error models and the posterior of the seismic moment tensor is estimated and interpreted in terms of the most probable source type.

58 GEOSCIENCES↗
DOE OSTI2024

Exploring the Nature of f-Element Soft Donor Interactions Using Electronically Tunable Azolate Ionic Liquids

This project was undertaken to advance the understanding of how f-elements interact with moderately soft donors, a heavily investigated yet open question which is of prime importance in spent nuclear fuel processing and fundamental inorganic chemistry. During the course of the project, based on exciting results, a stretch goal was developed to study the hydrolysis products of transuranic actinide metals, a somewhat understudied field even with its significance in nuclear fuel cycle and impacts in environmental chemistry. The stretch goal was to take our serendipitous discovery of an easy route to isolation of crystalline multinuclear ƒ-element hydrolysis products, and apply it to gaining a mechanistic understanding of Pu(III/IV) colloid formation. The simplicity of our techniques should lend themselves to the remote handling required for study of many transuranic elements. We developed several methodologies using azolium azolate chemistry to overcome ƒ-element hydrolysis problems that prohibit the isolation of ƒ-element soft donor complexes and to build a crystallographic library of ƒ-element N-donor complexes as a means to understand the fundamental differences between actinide and lanthanide interactions with moderately soft donor ligands. Our next major endeavor will be to transfer this chemistry from 4ƒ elements to transuranic elements, particularly in the study of hydrolysis of Pu(III/IV). While our work is fundamental in nature, applications of the knowledge we are generating should be felt in such diverse fields as catalysis, separations in general, nuclear waste remediation specifically, and many other applications in f-element magnetic and luminescent properties. The potential ramifications of the consistent and robust formation of hydrolysis controlled hexanuclear lanthanide structures are enormous, with future uses being catalyst formation, higher-nuclearity structure synthesis using our hexanuclear motif as a fundamental building block, Pu waste remediation, separations, and many other potential applications resulting from characteristic magnetic and luminescent properties of lanthanide polynuclear structures. Three synthetic methodologies (direct mixing with variable stoichiometries, use of volatile solvent, metathesis) were developed starting with 7 acidic and 6 basic azoles to obtain ionic liquids suitable for f-element coordination. Proton transfer by acidic/basic azole combination led to suitable low melting salts and two cocrystals. Acid/base reaction of azoles with soft-donor permanent cations of ([X 4444 ][OH] (where [X 4444 ] + = tetrabutylammonium [N 4444 ] + or tetrabutylphosphonium [P 4444 ] + ) with weakly acidic azoles including imidazole, 1,2,3-triazole, 1,2,4-triazole, 5-aminotetrazole, 4,5-dicyanoimidazole, and 2-amino-4,5-dicyanoimidazole) revealed several suitable low-melting salts. Metathesis reactions of Na(azolate) were conducted by first using weakly acidic azoles including 4,5-dicyanoimidazole, 2-amino-4,5-dicyanoimidazole, 5-aminotetrazole, and 1,2,4-triazole to form sodium or lithium salts using group(I) hydroxides in methanolic solutions. The best results were obtained by reacting the basic and acidic azoles directly in 1:1 or 3:1 ratios at elevated temperatures. Twenty-two azole mixtures were identified which are either low melting solids or room temperature liquids. Each of the low melting solids was confirmed to be a new solid phase, each of which is being further characterized. The liquids and solids are anticipated to be ILs, eutectics, or partially ionized systems, all of which will be suitable for the dissolution of f-element salts. Five new synthetic methodologies were developed to finding suitable crystallization conditions for f-element complexation with the goal of finding simple one pot reaction syntheses and crystallization strategies that could be used under the demanding conditions of transuranic chemistries. These synthetic methods yield many new crystalline phases which were studied by single crystal X-ray diffraction.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Modeling graphene sheet growth and dynamical matrix calculations using molecular dynamics

Molecular dynamics (MD) has been an incredibly useful tool to model physical processes that were synthesized experimentally but not fully understood. MD, through the use of semi-empirical inter-atomic potentials, has allowed understanding of different physical processes in materials science. Yet as well as providing useful insights into materials science, molecular dynamics has a wider range of usability. In this report, I will be detailing how MD can be used to study graphene formation from a carbon liquid which requires high temperatures and pressures. Beyond this, I will describe the usefulness of MD for understanding the physics for phonon transport quantum sensors. To do this, MD was employed to determine the dynamical matrix by treating atoms as coupled oscillators. An accurate understanding of the dynamical matrix of a system is required to calculate the non-equilibrium Green’s function used to describe the phonon transport within phonon wave-guides. I found that, across multiple pressures and temperatures, randomly placed carbon atoms will show evidence of pent-first formation with semi-empirical models. Density functional theory (DFT), on the other hand, was too computationally expensive to use for full scale MD simulations, but we have the possibility of training a machine learned interatomic potential to approximate DFT for carbon in the environments being studied for pent-first graphene sheet formation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

End-Use Savings Shapes Upgrade Package Documentation: Wall and Roof Insulation and New Windows

Building on the successfully completed effort to calibrate and validate the U.S. Department of Energy's ResStock and ComStock models over the past 3 years, the objective of this work is to produce national data sets that empower analysts working for federal, state, utility, city, and manufacturer stakeholders to answer a broad range of analysis questions. The goal of this work is to develop energy efficiency, electrification, and demand flexibility end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) upgrade measures, or upgrades. "Measures" refers to energy efficiency variables that can be applied to buildings during modeling. An end-use savings shape is the difference in energy consumption between a baseline building and a building with an energy efficiency, electrification, or demand flexibility upgrade applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ComStock is a highly granular, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the commercial building stock across the United States. The baseline model intends to represent the U.S. commercial building stock as it existed in 2018. The methodology and results of the baseline model are discussed in the final technical report of the End-Use Load Profiles project. This documentation focuses on an upgrade package of three end-use savings shapes upgrades - Window Replacement, Exterior Wall Insulation, and Roof Insulation, which we will refer to collectively as the "High Efficiency Envelope" package. More details on the individual upgrades can be found on the ComStock Measures Documentation page. An upgrade package applies two or more EUSS upgrades to a single building model simulation. Since ComStock is a bottom-up physics-based model, an upgrade package will go beyond aggregating or summing the individual upgrade results and produce novel results by simulating interactions between the upgrades. For example, pairing an envelope upgrade with an electrification upgrade would likely result in higher savings results than the sum of these upgrades individually, and the size of the heating, ventilating, and air conditioning (HVAC) equipment may be reduced if the envelope upgrade reduces the loads significantly.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

End-Use Savings Shapes Upgrade Package Documentation: LED Lighting, HP-RTU and ASHP-Boiler

Building on the successfully completed effort to calibrate and validate the U.S. Department of Energy's ResStock and ComStock models over the past 3 years, the objective of this work is to produce national data sets that empower analysts working for federal, state, utility, city, and manufacturer stakeholders to answer a broad range of analysis questions. The goal of this work is to develop energy efficiency, electrification, and demand flexibility end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) upgrade measures, or upgrades. "Measures" refers to energy efficiency variables that can be applied to buildings during modeling. An end-use savings shape is the difference in energy consumption between a baseline building and a building with an energy efficiency, electrification, or demand flexibility upgrade applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ComStock is a highly granular, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the commercial building stock across the United States. The baseline model intends to represent the U.S. commercial building stock as it existed in 2018. The methodology and results of the baseline model are discussed in the final technical report of the End-Use Load Profiles project. This documentation focuses on an upgrade package of three end-use savings shapes upgrades - Light Emitting Diode (LED) Lighting, Heat Pump Rooftop Unit (RTU) (HP-RTU), and Air-Source Heat Pump (ASHP) Boiler, which we will refer to collectively as the "Interior Lighting and Heat Pump" package. More details on the individual upgrades can be found on the ComStock Measures Documentation page. An upgrade package applies two or more EUSS upgrades to a single building model simulation. Since ComStock is a bottom-up physics-based model, an upgrade package will go beyond aggregating or summing the individual upgrade results and produce novel results by simulating interactions between the upgrades. For example, pairing an envelope upgrade with an electrification upgrade would likely result in higher savings results than the sum of these upgrades individually, and the size of the heating, ventilating, and air conditioning (HVAC) equipment may be reduced if the envelope upgrade reduces the loads significantly.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Using Parameter Sweep in WaterTAP to Analyze New Water Treatment Technologies

We describe a powerful and generalized parameter sweep tool in this report that was originally developed to analyze the performance of existing and novel water treatment models being developed in WaterTAP. Since WaterTAP is built upon IDAES and Pyomo, the parameter sweep tool can be used to systematically explore and debug the behavior of most Pyomo and IDAES numerical models. In order to enable meaningful analyses, the parameter sweep tool has been designed with the following features: 1) Model flexibility: The parameter sweep tool does not enforce any restrictions on the types of models that can be used with it. As long as a Pyomo model can be solved and the parameter is active and mutable, the tool only needs functions that describe how to run the model, the sweep parameters, and the output quantities of interest. 2) Flexible sampling: The parameter sweep tool has inbuilt functions to generate samples from a random distribution or a multidimensional Euclidean space. Furthermore, the users have to ability to supply samples generated from a tool of their choice. 3) Multiple sweep types: A user can choose from one of 3 types of parameter sweeps depending on their needs. 4) Detailed outputs: Outputs generated by the parameter sweep tool can be stored in detailed H5 file or user-friendly CSV files for post processing. 5) Parallel computing: The parameter sweep supports shared and distributed memory parallel computing to enable the use of high performance computers (HPC) for large-scale analyses. 6) Modular: The parameter sweep tool is self-contained and can easily be integrated within an outer-loop analysis or as desired by the user. 7) Ease of use: The tool is well documented and a simple sweep can be easily executed by following the online documentation in a few lines of code. We demonstrate the use of the parameter sweep tool on a simple water treatment system from the WaterTAP repository and show its parallel scaling performance on an Apple laptop and NREL's Eagle HPC. The parameter sweep tool is actively being used with models currently being developed within WaterTAP and we expect its use to grow beyond it to other IDAES and Pyomo models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Alpha-Quartz Plastic Strength Investigation Via Diffraction Experiments on Novaculite Using A D-Dia and Elastic Plastic Self- Consistent Interpretation [Thesis}

X-Ray Tomography. Porosity is commonly measured using mercury injection (MI) or water immersion porosimeter (WIP). Both MI and WIP utilize pressure to fill open voids with mercury or water respectively. A measure of the volume change of the sample is then used to estimate the percentage of open voids. However, for the purpose of rheological study it is important to get an accurate measure of open and closed voids. X-ray tomography has been selected as it offers a three-dimensional view into the sample that can quantify all pores limited only by the voxel size which is in the micron range for this study. Radiographs for tomography were collected at the Material Science and Technology Division at Los Alamos National Lab using the Carl Zeiss Xradio 520 instrument and Scout-and-Scan version 16.1 operating software. Two samples were imaged, the starting material and the deformed sample SiO2_65. 3001 radiographs were taken of the starting material with a 6 second exposure time using a 4x objective lens. The x-ray beam was set to 60 kilovoltage peak (kVp) and 5 watts. 1901 radiographs were taken of SiO2_65 with a 25 second exposure time using a 10x objective lens. The x-ray beam was set to 80 kVp and 7 watts. Radiograph files were analyzed by Brian Patterson using Avizo. Void and inclusion volumes were output by voxel sized (1.03 μm) slices, used to calculate a total percentage volume for the starting material.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Using SPHINX Accelerator to Measure Heat Transfer in Metal Foils [Slides]

Experiments on the SPHINX accelerator at Sandia National Laboratories (SNL) looked at electron beam heating of metal foils used in Bremsstrahlung flash x-ray production. Titanium and tantalum foils were exposed to electron beam energies of up to 2MeV to determine surface heating. As part of the experiment, the x-ray dose profiles were characterized using a scintillating panel (EJ-230) placed both perpendicular and parallel to the photon beam. This was followed by direct electron beam imaging using a fast-gated ICCD camera and a fused silica Cherenkov witness plate. Once characterized, the foil was directly measured using a mid-wave (3-5 micron) infrared camera (MWIR) to determine the peak surface temperatures and spatially-resolved thermal profiles. Additionally, an IR pyrometer was used in a single wavelength mode (2.3 µm) to determine localized surface temperatures as well. Results are compared with calculations of surface temperatures obtained from e-beam energy depositions combined with the specific heats of the foils. This is the first attempt to directly measure the surface temperature profile on a short-pulse (10 nanoseconds) flash x-ray source using MWIR imaging. These data are directly relevant for Monte Carlo/Integrated Tiger Series electron/photon transport code results being used in the design of future Bremsstrahlung sources such as CREST.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Multiplexing and Demultiplexing Signals for Radiography Application Using the Discrete Fourier Transform

Our goal is to develop an X-ray phase-contrast imaging system that can provide excellent soft tissue contrast of phase, attenuation, and small-angle scatter. We propose to replace the common system of G0, G1, and G2 gradings with a biprism array to replace the G1 grading and introduce a novel X-ray tube designed to replace the motion of the phase stepping grading G2. The proposed X-ray tube uses temporal multiplexing to provide simultaneous virtual “electronic phase stepping.” In this work the discrete Fourier transform is used to separate from the composite measurement individual X-ray phase contrast measurements sampled at different frequencies. The method performs a discrete Fourier transform of a composite refence sequence to obtain using the frequency amplitudes calibration factors needed to extract the X-ray phase contrast measurement amplitudes from the composite image. The composite reference sequence is the sum of the individual sequences, at different frequencies, with amplitudes of one. The method takes the discrete Fourier transform of this composite reference sequence; whereby, the amplitude of each frequency component is compared with the total sum of its stand-alone sequence amplitude. A calibration factor is determined so that the amplitude of this composite reference frequency times the calibration factor must equal the total sum of the sequence amplitude—the zero-frequency amplitude of the discrete Fourier transform of its stand-alone sequence. To demultiplex the composite measured signal these calibration factors are multiplied by the amplitudes of the frequency components of the discrete Fourier transform of the composite X-phase-contrast measurement to obtain the amplitude of each frequency encoded measurement. Using these calibration factors, we demonstrate with the discrete Fourier transform in Mathematica the extraction of individual images from a composite image that one would expect obtaining from our proposed new X-ray phase contrast imaging system. We then demonstrate as an example how using images from X-ray phase contrast data one can calculate phase, attenuation and the dark field images using grading phase step data supplied to use from Microworks, GmbH in Karlsruhe, Germany.

42 ENGINEERING↗
DOE OSTI2024

Precision Agriculture using Networks of Degradable Analytical Sensors (PANDAS) (Final Technical Report)

Precision agriculture, where sensing of soil, environment and crop conditions are used to precisely synchronize inputs (such as water and fertilizer) to crop needs enhances input use efficiency. This can improve yields and farm profitability while mitigating environmental losses, improving soil carbon content and substantially decreasing energy use for food, feed and fuel crops. Unfortunately, farmers are not yet able to harness the full potential of these management technologies as there is a lack of available management information, and there is therefore a need for sensors that are able to economically measure spatio-temporal variability in soil and crop properties of extremely heterogeneous farm fields precisely at high resolution and at low cost. Real-time, in-situ monitoring of agricultural soil conditions is today carried out using devices that limit the total number of nodes that can be used economically to typically one per acre or less. Higher spatio-temporal resolution sensing would enable more precise agricultural input optimization, with significant benefits to the farmer and the environment. In order to address this issue, this project focused on developing additively manufactured, biodegradable, soil sensors with predicted costs of < $\$$1 per unit to monitor crop inputs (such as water and fertilizer) that predictably, harmlessly degrade away into the soil when no longer needed. These sensor nodes should be easy to place, accurately and continuously monitor soil and crop conditions for an entire season, be read remotely using existing farm equipment, require no ongoing maintenance, not impede farm operations and produce no persistent waste. This approach could enable a >100× increase in information density over current solutions for precision farming of row and other crops, and lead to significant reductions in input energy use and provide increased yield for biofuel crops. Over the course of this project the team at the University of Colorado Boulder, University of California Berkeley, and Colorado State University/Kansas State University investigated a wide range of printable biodegradable electronic materials and sensor designs for determining soil moisture and soil nitrate concentration. These efforts expanded the available materials set for printed soil degradable electronic materials, particularly for conductors, enabling high conductivity and stability. Printed soil moisture and nitrate sensors with suitable sensitivity and selectivity were developed and characterized. Low power and passive wireless electronic systems were integrated with the soil sensors, and testing was carried out with completed sensors to understand their functionality under agricultural conditions. Additionally, other sensor types enabled by the biodegradable materials set created during this project, such as soil microbial activity sensors, were also developed and demonstrated. Project outputs include 10 peer reviewed publications, 4 patent applications, 21 technical presentations, 3 PhD thesis, 10 media reports, 8 additional grants worth over $\$$6M, and the formation of 3 start-up companies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Using Machine Learning to Improve Thermostability of MHETase

Protein engineering is a field which utilizes proteins as tools, which has many useful applications in medicine, industry, biofuels and more.1 One such protein is MHETase, which is a protein that plays an important function in the degradation of polyethylene terephthalate (PET) plastics, which are commonly used in water and soda bottles.2 However, these proteins are adapted to work in specific conditions, and may not satisfy the desired properties that a new application would desire, or could be improved. For instance, a more thermostable MHETase would be more effective in the plastic degradation conditions.3 To make these desired changes, the primary structure of the protein is mutated, but there are many possible mutations and positions to mutate to make with the 20 canonical amino acids. Therefore, to narrow down the possibilities and to make the process of finding a thermostable MHETase variant, we used sequence design tools that are grounded in machine learning to find mutations that would improve thermostability of MHETase.4 In particular, we used the tools Protein MPNN and Fireprot to design a more thermostable MHETase enzyme. We then compiled these mutations into a library and grew these proteins using bacteria colonies, and measured their effectiveness using a fluorescent protein marker. Thermostable proteins and their marker would fold correctly and fluorescence would be seen, but if neither folded correctly then there would be no marker detected. We grew these proteins in bacteria and then intend to use these methods to evaluate their thermostability.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Using Machine Learning to Improve Thermostability of MHETase

Protein engineering is a field which utilizes proteins as tools, which has many useful applications in medicine, industry, biofuels and more. One such protein is MHETase, which is a protein that plays an important function in the degradation of polyethylene terephthalate (PET) plastics, which are commonly used in water and soda bottles.2 However, these proteins are adapted to work in specific conditions, and may not satisfy the desired properties that a new application would desire, or could be improved. For instance, a more thermostable MHETase would be more effective in the plastic degradation conditions.3 To make these desired changes, the primary structure of the protein is mutated, but there are many possible mutations and positions to mutate to make with the 20 canonical amino acids. Therefore, to narrow down the possibilities and to make the process of finding a thermostable MHETase variant, we used sequence design tools that are grounded in machine learning to find mutations that would improve thermostability of MHETase.4 In particular, we used the tools Protein MPNN and FireProt to design a more thermostable MHETase enzyme. We then compiled these mutations into a library and grew these proteins using bacteria colonies, and measured their effectiveness using a fluorescent protein marker. Thermostable proteins and their marker would fold correctly and fluorescence would be seen, but if neither folded correctly then there would be no marker detected. We grew these proteins in bacteria and then intend to use these methods to evaluate their thermostability.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Evaluation of DNA Extraction Efficiency in Diverse Algae Strains Using Commercial Kits and Lysis Approaches

Efficient DNA extraction is essential for accurately monitoring microalgae communities in large-scale cultivation systems such as raceway ponds and wastewater ponds. Traditional phenol chloroform extracts are a staple in microbiology but are obsolete for routine sampling due to its high toxicity reagents and time intensive setups. Commercial DNA extraction kits are more favorable for the microbes found in these ponds, but lack specific kits made for these communities. Little is known about which kits perform the best, leading researchers to use a variety of different kits with inconsistent results. This project compared one precipitation based commercial kit (Lucigen Masterpure) and five wash based kits (Monarch, Zymo Quick-DNA, and three Qiagen DNeasy kits) using four brackish algae strains to determine which methods yield the greatest quantity and quality of genomic DNA. Extractions were evaluated using the manufacturers protocol, and additional pretreatment options were administered before a single kit to compare its potential in being added routinely before extractions. Pretreatment options included both cryogenic freeze-thawing and heat incubation using enzymes. DNA was quantified using Qubit fluorometry and NanoDrop purity ratios. Overall, the Qiagen PowerWater kit provided the highest DNA yield and purity, but at a significantly higher cost then the precipitation-based kit (MasterPure). It was also noted that while the precipitation-based kit was significantly cheaper, provided similar results, it took significantly more time to complete a single run. Cryogenic pretreatment (6x cycles) increased average DNA yields by up to 80%, whereas enzymatic pretreatment most improved purity ratios without substantially improving quantity. The results suggest that it may be more cost and time efficient to use Qiagen kits with the addition of lysis pretreatments to procure better results. Future works includes developing a better system to efficiently collect multi variable data, and to upscale to artificial polycultures using similar methodologies alongside sequencing to confirm kit results.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Simulation of Physics-Based 0-10Hz Strong Motion Using High Performance Computing Supporting Refinements to Regional Ground Motion Models for the Central Eastern US

In collaboration with the U.S. Nuclear Regulatory Commission (NRC) the LLNL has developed a computationally efficient simulation platform designed to perform physics-based ground motion simulations for crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using high-performance computing. The main objective of the earthquake simulations was to use synthetic ground motion to provide constrains to refinements of existing ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-10Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model (1994). In simulations performed during the first stage of this project the Saikia’s velocity model demonstrated better performance in modelling high frequency regional wave propagation for the CEUS region recorded during the Mw5.0 November 7, 2016, Cushing Oklahoma (Taylor et al., 2017), and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes. The proposed regional 3D model includes random perturbations to the 1D background model using the stochastic scheme of Pitarka and Mellors (2021). In addition, validation analysis of the rupture generator and regional wave propagation models, using comparisons with different GMMs for Mw6.5 and Mw7.0 scenario earthquakes in the CEUS region resulted in a very good match between the simulated and empirical ground motion models. For the purposes of seismic hazard assessment at the existing and planned nuclear power plants, NRC is interested in studies aimed at improving the current ground motion models (GMM) for both Stable Continental Regions (SCR) in the Central and Eastern US and Active Crustal Regions (ACR) in the Western US. Due to lack of recorded data, these improvements require synthetic data for short fault distances and large magnitude earthquakes for which the existing recorded data is not enough to uniquely constrain the GMMs. The need for simulations and strong motion data is especially critical for the CEUS region where we do not have recorded data from potentially large damaging earthquakes with moment magnitudes 6.0 and higher. In this the project, we focused on 10Hz simulations of Mw7.0 scenario earthquakes with strike slip and thrust faulting mechanisms. We used more than 50 Mw7.0 earthquake rupture scenarios to investigate the ground motion uncertainty due to unknown earthquake rupture parameters, in particular, the slip distribution, rupture velocity, and faulting mechanism, and their implication on ground motion amplification due to forward rupture directivity effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

GeoThermalCloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.

15 GEOTHERMAL ENERGY↗