Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Control variates”

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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Heterogeneous oxidation of amorphous organic aerosol surrogates by O 3 , NO 3 , and OH at typical tropospheric temperatures

Typical tropospheric temperatures render possible phase states of amorphous organic aerosol (OA) particles of solid, semisolid, and liquid. This will affect the multiphase oxidation kinetics involving the organic condensed-phase and gaseous oxidants and radicals. To quantify this effect, we determined the reactive uptake coefficients (γ) of O 3 , NO 3 , and OH by substrate films composed of single and binary OA surrogate species under dry conditions for temperatures from 213 to 313 K. A temperature-controlled coated-wall flow reactor coupled to a chemical ionization mass spectrometer was applied to determine γ with consideration of gas diffusion transport limitation and gas flow entrance effects, which can impact heterogeneous reaction kinetics. The phase state of the organic substrates was probed via the poke-flow technique, allowing the estimation of the substrates' glass transition temperatures. γ values for O 3 and OH uptake to a canola oil substrate, NO 3 uptake to a levoglucosan and a levoglucosan / xylitol substrate, and OH uptake to a glucose and glucose / 1,2,6-hexanetriol substrate have been determined as a function of temperature. We observed the greatest changes in γ with temperature for substrates that experienced the largest changes in viscosity as a result of a solid-to-liquid phase transition. Organic substrates that maintain a semisolid or solid phase state and as such a relatively higher viscosity do not display large variations in heterogeneous reactivity. From 213 to 293 K, γ values of O 3 with canola oil, of NO 3 with a levoglucosan / xylitol mixture, and of OH with a glucose / 1,2,6-hexanetriol mixture and canola oil, increase by about a factor of 34, 3, 2, and 5, respectively, due to a solid-to-liquid phase transition of the substrate. These results demonstrate that the surface and bulk lifetime of the OA surrogate species can significantly increase due to the slowed heterogeneous kinetics when OA species are solid or highly viscous in the middle and upper troposphere. This experimental study will further our understanding of the chemical evolution of OA particles with subsequent important consequences for source apportionment, air quality, and climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Initial PIP-II Beam Current Monitor Fault Case Analyses & Beam Position Monitor Linearity Studies in CST Studio Suite

The use of non-invasive sensors & systems to measure particle beam characteristics is a crucial part of modern accelerator control systems due to their ability to return real time beam data while minimizing negative effects on beam quality. To ensure that one can be reasonably confident these sensors will behave as desired upon be-ing implemented within the beamline, simulations pre-dicting the performance of these sensors under beamline conditions can be used as a valuable tool for checking sensor functionality without a physical test bench. This paper details the design, testing, and results of two sensor models developed using CST Studio Suite soft-ware designed to mimic two sensors to be implemented within the PIP-III beamline: an elliptical, large-aperture beam position monitor (BPM) for which vertical & hori-zontal position signal linearity was analyzed, and an AC current transformer (ACCT) beam current monitor (BCM) used to search for potential fault cases within the BCM and beam pipe flange gaps. Special focus is given to the discovery of linearity variations within the BPM and the use of frequency domain techniques in the BCM fault case analyses.

Rouzky, A. R.↗
DOE OSTI2023

Multi-Timescale Integrated Dynamics and Scheduling for Solar (MIDAS-Solar) (Final Technical Report)

Solar photovoltaic (PV) installations have experienced unprecedented growth in the United States. PV will become not only an energy producer but also a necessary provider of ancillary services at multiple timescales. Conventional methods to simulate power systems operations - such as long-term production simulation (which typically considers schedules from hours to minutes by using an optimization framework) and short-term transient studies (which simulate dynamics from seconds to sub-seconds using state variables and differential equations) - are not sufficient for studying the multiple-timescale variation of solar generation and its impact on system reliability. Long-term system economics and short-term system dynamics are highly coupled, particularly when the penetration level of renewable generation is extremely high, because the uncertainty and variability of solar generation will impact both power system steady-state and dynamic performance. This project helps meet and exceed the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy Solar Energy Technologies Office goal of systems integration by directly addressing this stability and reliability challenge for power grid planning and operation. We have developed a temporally comprehensive, closed-loop simulation model, named Multi-timescale Integrated Dynamics and Scheduling (MIDAS), that seamlessly simulates power system operations from economic scheduling (day-ahead to hours) to dynamic response analysis (seconds to sub-seconds). For schedules with very high levels of inverter-based resources (IBRs), up to and including 100%, the stability of grid controls has been evaluated through electromagnetic transient (EMT) simulations and power-hardware-in-the-loop (PHIL) simulations of key transient events at key schedule points. Specifically, MIDAS provides: 1) a closed-loop simulation framework for simulating timescales from economic scheduling to dynamic stability analysis; 2) machine learning-based stability assessment; 3) EMT modeling and analysis for large-scale power systems; 4) MIDAS PHIL test bed. We worked with Hawaii Electric Companies to apply the MIDAS study framework to a Maui grid study. The entire island's transmission system was modeled in detail - from a yearly scheduling model, to a second-level frequency dynamic model, down to a sub-second-scale EMT model to address critical stability issues. The project demonstrated how MIDAS can help system planners and operators assess system reliability and stability while the power grid is marching toward a high-renewable, high-IBR future. In this Maui grid study, we found that 100% instantaneous IBR operation is achievable in EMT simulation and PHIL testing, and grid planners and operators might need new analysis/simulation tools to assess grid reliability and stability in the scheduling stage. MIDAS will bring Maui and other systems closer to 100% clean and stable energy futures. (In this study, we examined transient stability. Other topics necessary for 100% IBR operation, such as protection and resource adequacy, were not examined.)

100% Renewables↗
DOE OSTI2021

Interlayer Triplet Energy Transfer in Dion-Jacobson 2D Perovskites Containing Naphthalene Diammonium Cations

The ability to tune the heterojunction energy level alignment in layered 2D perovskites is well documented, with both early works and more recent studies reporting interfacial triplet energy transfer in certain pairings of lead halide with small conjugated organic spacer cations. Such 2D perovskites are promising systems for sensitized molecular phosphorescence, which has varied applications in optoelectronics, solar energy conversion, and more. However, even for a given organic cation core, variations in the cation structure (e.g. location of the ammonium substituent, length of the alkyl linker, and valency of the cation) can dramatically alter the organic layer morphology and electronic coupling between the inorganic and organic layers. While interlayer triplet energy transfer has been studied in Ruddlesden-Popper 2D perovskites containing monovalent naphthalene cations, the photophysical properties of their Dion-Jacobson analogue, formed using divalent naphthalene cations, have not been reported. In this study, we examine interlayer energy transfer in a series of mixed-halide 2D perovskites formed by divalent naphthalene cations. We find that the sensitized phosphorescence in these compounds is dominated by naphthalene triplet excimer emission, but when the lead halide exciton level is tuned near resonance with the triplet level of naphthalene, emission from the naphthalene triplet monomer competes with triplet excimer formation. Interlayer energy transfer in these compounds is further supported by ultrafast transient absorption spectroscopy. Ultimately, the ability to gain control over interlayer interactions in layered 2D perovskites through cation design will help uncover new functions and applications for these materials.

2D↗
DOE OSTI2026

Using Best Basis Inventory Data to Direct Strategies for Real Time Monitoring of Hanford High Level Waste – 26226

The potential to accelerate the processing of low- and high-level tank waste by applying real-time monitoring (RTM) of chemical and physical properties has prompted research into the suitability of multiple analytical methods for that purpose. The broad variety of waste stream properties and the large number of analytes of interest (as evidenced by Waste Acceptance Criteria (WAC) and Process Control Limit (PCL) lists) lead to an overwhelming set of possible analytical scenarios. This report describes the use of Best Basis Inventory (BBI) data to find the most relevant analytical targets for the specific case of monitoring the blending of High Level Waste from multiple tanks prior to introduction into a vitrification facility. Campaigns for blending this waste to minimize the risk of exceeding WACs and PCLs have been proposed. However, the predicted compositions of the blended materials do not incorporate any uncertainties that may be associated with the representativeness of the waste layer samples or the laboratory analyses that generated the BBI data. Also, they do not include any uncertainty associated with the precision of collecting highly specific fractions of the layers during a blending campaign or any inhomogeneities that may exist in those layers. Monte Carlo methods are used to apply uncertainties to the compositions of the individual layers specified in the campaign recipes. The resulting variations in the compositions of the blended materials allow estimation of the risks of exceeding WACs and PCLs for each campaign. A critical subset of WACs/PCLs – NOx, NaK, AlFeZr, and S – are especially at risk of being exceeded in multiple campaigns. These analytes should be the focus of instrument development. We also have extracted the expected solid/supernate distribution for these analytes, which establishes important performance criteria for individual analytical methods. The BBI data also permits an understanding of the different chemical forms in which the analytes appear. Thus, the need to establish instrumental sensitivity to these forms can be gainfully addressed. Although concentrating on one specific application – the blending of tank waste - this approach should be generalizable for the analysis of other possible RTM applications for waste processing.

Lascola, Robert [Savannah River National Laborator↗
DOE OSTI2023

Disruption event characterization and forecasting in tokamaks

Disruption prediction and avoidance is a critical need for next-step tokamaks, such as ITER. Disruption Event Characterization and Forecasting (DECAF) research fully automates analysis of tokamak data to determine chains of events that lead to disruptions and to forecast their evolution allowing sufficient time for mitigation or complete avoidance of the disruption. Disruption event chains related to local rotating or global magnetohydrodynamic (MHD) modes and vertical instability are examined with warnings issued for many off-normal physics events, including density limits, plasma dynamics, confinement transitions, and profile variations. Along with Greenwald density limit evaluation, a local radiative island power balance theory is evaluated and compared to the observation of island growth. Automated decomposition and analysis of rotating tearing modes produce physical event chains leading to disruptions. A total MHD state warning model comprised of 15 separate criteria produces a disruption forecast about 180 ms before a standard locked mode detector warning. Single DECAF event analyses have begun on KSTAR, MAST, and NSTX/-U databases with thousands of shot seconds of device operation using from 0.5 to 1 × 10 6 tested sample times per device. An initial multi-device database comparison illustrates a highly important result that plasma disruptivity does not need to increase as βN increases. Global MHD instabilities, such as resistive wall modes (RWMs), can give the briefest time period of warning before disruption compared to other physics events. In an NSTX database with unstable RWMs, the mode onset, loss of boundary and current control, and disruption event warnings are found in all cases and vertical displacement events are found in 91% of cases. An initial time-dependent reduced physics model of kinetic RWM stabilization created to forecast the disruption chain predicts instability 84% of the time for experimentally unstable cases with a relatively low false positive rate. Instances of the disruption event chain analysis illustrate dynamics including H–L back transitions for rotating MHD and global RWM triggering events. Furthermore, disruption warnings are issued with sufficient time before the disruption (on transport timescales) to potentially allow active profile control for disruption avoidance, active mode control, or mitigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

From tides to seasons: How cyclic tidal drivers and plant physiology interact to affect carbon cycling at the terrestrial-estuarine boundary (Final technical report)

Coastal ecosystems are among the most biologically and biogeochemically active and diverse systems on Earth. Because they act as important linkages between terrestrial ecosystems and the open ocean, their incorporation in Earth system models (ESMs) is critical to predict coastal and global responses to environmental changes. However, they vary greatly in the magnitude of tides and the volume and timing of freshwater input from land, making it challenging to model the major biogeochemical reactions that control productivity and greenhouse gas emissions across coastal terrestrial aquatic interfaces (TAIs). Our overall objective was to improve mechanistic process understanding and modeling of tidal wetland hydro-biogeochemistry in coastal TAIs. We established a new flux tower site (Ameriflux US-PLo) in the oligohaline part of the Parker River to continuously monitor ecosystem-scale carbon fluxes under temporally varying salinity conditions. The site is co-located with long-term monitoring plots of the Plum Island Ecosystems LTER project. We installed wells and redox sensors in the marsh interior and creek bank, established biomass monitoring plots and deployed novel optode sensors in both locations. We used this data to parameterize plant-mediated transport in PFLOTRAN and tested the impact of soil heterogeneity on porewater constituents and gas fluxes. We collected observations of root oxygen release with a novel planar optode system in the field. Flux data collected during the measurement period encompasses a large variation in salinity ranging from drought to record precipitation years. We developed a method to extract functional relationships from the flux data using artificial neural networks, identifying salinity thresholds for CH 4 fluxes. Finally, we are using the coupled ELM-PFLOTRAN model to test the impact of antecedent hydrological conditions on the salinity-CH 4 flux relationship. This grant contributed to the professional development of one postdoc, three research assistants and one graduate student. The sensor data has been shared with external collaborators.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Effects of Cinnamomum camphora coppice planting on soil fertility, microbial community structure and enzyme activity in subtropical China

Cinnamomum camphora ( C. camphora ) is a broad-leaved evergreen tree cultivated in subtropical China. Currently, the use of C. camphora clonal cuttings for coppice management has become popular. However, the effects of C. camphora coppice planting on soil abiotic and biotic variances remained unclear. In this study, we collected soil from three points in the seven-year C. camphora coppice planting land: under the tree canopy (P15), between trees (P50), and abandoned land (Control) to investigate the effects of C. camphora coppice planting on soil fertility, microbial community structure and enzyme activity. The results revealed that C. camphora coppice planting significantly increased soil fertility in the point under the tree canopy (P15) and point between trees (P50), and P15 had more significant effects than P50. Meanwhile, in P15 and P50, soil bacterial, fungal alpha-diversity were improved and microbial community structures were also changed. And the changes of soil organic carbon and total nitrogen promote the transformation of soil bacterial, fungal community structures, respectively. In addition, C. camphora coppice planting significantly ( p < 0.05) increased soil urease (UE), polyphenol oxidase, and peroxidase activities, while significantly decreased soil ACP activity. This study demonstrated that the C. camphora coppice planting could improve soil fertility in subtropical China, which promoted the transformation of soil microbial community from oligotrophs ( K -strategist) to copiotrophs ( r -strategist). Thus, this work can provide a theoretical basis for soil nutrient variation and productive management of C. camphora coppice plantation in subtropical China.

Sun, Luyuan↗
DOE OSTI2026

Spatially resolved phase transformation mapping in 410 stainless steel during additive manufacturing visualized via real-time infrared data

Here, in this study, we demonstrate that spatially resolved cooling curves derived from real-time infrared (IR) thermography during additive manufacturing (AM) can capture spatial variations in phase transformation temperatures through cooling curve analysis (CCA). Using this approach, we show that during laser hot-wire deposition of 410 stainless steel (410SS), the martensite start temperature (M s ) evolves dynamically throughout the build. The M s temperature is spatially nonuniform, ranging from 185 °C to 348 °C, with the lowest values toward the build center and higher values toward the upper region of the deposit. In the lower portion of the build, no M s inflection is detected via CCA, consistent with transformation occurring earlier during thermal cycling followed by tempering during subsequent thermal cycles. These trends in M s are corroborated by characterizing the microstructure by electron backscatter diffraction (EBSD). Traditionally, M s is assumed to be constant, and a single interpass temperature is applied during both deposition and residual stress modeling. Our results demonstrate that IR-derived cooling curves provide a route to spatially and temporally resolved transformation temperature tracking for dynamic interpass control and improved residual-stress modeling.

Additive manufacturing↗
DOE OSTI2024

Modeling the effect of wind speed and direction shear on utility‐scale wind turbine power production

Abstract Wind speed and direction variations across the rotor affect power production. As utility‐scale turbines extend higher into the atmospheric boundary layer (ABL) with larger rotor diameters and hub heights, they increasingly encounter more complex wind speed and direction variations. We assess three models for power production that account for wind speed and direction shear. Two are based on actuator disc representations, and the third is a blade element representation. We also evaluate the predictions from a standard power curve model that has no knowledge of wind shear. The predictions from each model, driven by wind profile measurements from a profiling LiDAR, are compared to concurrent power measurements from an adjacent utility‐scale wind turbine. In the field measurements of the utility‐scale turbine, discrete combinations of speed and direction shear induce changes in power production of −19% to +34% relative to the turbine power curve for a given hub height wind speed. Positive speed shear generally corresponds to over‐performance and increasing magnitudes of direction shear to greater under‐performance, relative to the power curve. Overall, the blade element model produces both higher correlation and lower error relative to the other models, but its quantitative accuracy depends on induction and controller sub‐models. To further assess the influence of complex, non‐monotonic wind profiles, we also drive the models with best‐fit power law wind speed profiles and linear wind direction profiles. These idealized inputs produce qualitative and quantitative differences in power predictions from each model, demonstrating that time‐varying, non‐monotonic wind shear affects wind power production.

Energy & Fuels↗
DOE OSTI2021

EXPRES. II. Searching for Planets around Active Stars: A Case Study of HD 101501

By controlling instrumental errors to below 10 cm s{sup −1}, the EXtreme PREcision Spectrograph (EXPRES) allows for a more insightful study of photospheric velocities that can mask weak Keplerian signals. Gaussian processes (GP) have become a standard tool for modeling correlated noise in radial velocity data sets. While GPs are constrained and motivated by physical properties of the star, in some cases they are still flexible enough to absorb unresolved Keplerian signals. We apply GP regression to EXPRES radial velocity measurements of the 3.5 Gyr old chromospherically active Sun-like star, HD 101501. We obtain tight constraints on the stellar rotation period and the evolution of spot distributions using 28 seasons of ground-based photometry, as well as recent Transiting Exoplanet Survey Satellite data. Light-curve inversion was carried out on both photometry data sets to reveal the spot distribution and spot evolution timescales on the star. We find that the >5 m s{sup −1} rms radial velocity variations in HD 101501 are well modeled with a GP stellar activity model without planets, yielding a residual rms scatter of 45 cm s{sup −1}. We carry out simulations, injecting and recovering signals with the GP framework, to demonstrate that high-cadence observations are required to use GPs most efficiently to detect low-mass planets around active stars like HD 101501. Sparse sampling prevents GPs from learning the correlated noise structure and can allow it to absorb prospective Keplerian signals. We quantify the moderate to high-cadence monitoring that provides the necessary information to disentangle photospheric features using GPs and to detect planets around active stars.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Issue Summary of INL Phase IV Transient Results for IAEA CRP on HTGR UAM Benchmark

This report details the Parallel and Highly Innovative Simulation for Idaho National Laboratory (INL) Code System (PHISICS)/Reactor Excursions and Leak Analysis Program (RELAP5)-3D results obtained for the transient core exercises defined for Phase IV of the International Atomic Energy Agency (IAEA) Coordinated Research Project (CRP) on high-temperature gas cooled reactor (HTGR) uncertainty analysis in modeling (UAM). The Phase III models and results are linked to the earlier Standardized Computer Analyses for Licensing Evaluation (SCALE)/Sampler/New ESC-based Weighting Transport (NEWT) data generated for the lattice physics (lattice) stage Phase I of the CRP. The focus of this report is the Uncertainty/Sensitivity Assessment (U/SA) of the prismatic modular high-temperature gas cooled reactor (MHTGR)-350 design, and specifically for Exercises IV-1 and IV-2 of the benchmark: the Control Rod Withdrawal (CRW) and Pressurised Loss of Cooling (PLOFC) events. The statistical U/SA methodology is implemented and demonstrated using the RAVEN code, based on perturbed cross-section libraries obtained from the SCALE/Sampler sequence. Uncertainties in nuclear data (cross-sections and the average number of neutrons produced per fission, 235U[¯v ]) lead to standard deviations (uncertainties of one s) of approximately 0.5% in the core eigenvalues of the MHTGR-350 and core models. For the coupled neutronics/thermal fluid model, local power density uncertainties up to 3.6% were observed in the colder regions of the core, while the local maximum fuel temperature uncertainties reached 1.5% for the models that included thermal fluid uncertainties. The addition of thermal fluid uncertainties dominated the impacts of nuclear data uncertainties in all cases. The main contributors to uncertainties in the power density and fuel temperatures during the transients were uncertainties in the reactor operating conditions (total power, inlet mass flow rate and inlet gas temperature). Variations in the bypass flows did not have significant impact on any of the output variables. For the nuclear data uncertainties it was found that the 235U(¯v ) / 235U(¯v ) covariance produced the largest sensitivities in terms of its impact on the eigenvalue and peak reactor power. It was also observed that the impact of any nuclear data uncertainties on the maximum fuel temperature was much less significant that the impact on eigenvalue and power. Another important finding was that although the use of eight or more energy groups is recommended for best-estimate HTGR simulation, two-group models produced acceptable uncertainty and sensitivity results for most FOMs. Since the statistical U/SA methodology is computationally expensive, and most transient solver requirements will scale directly with the number of energy groups, two energy groups could be used by HTGR developers during the early stages of design when larger uncertainty margins can be tolerated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Signature Analysis Utilizing a Dynamic Molten Salt Reactor Model for MC&A

Moving from traditional fixed fuel nuclear reactor systems to a mobile, dynamic fuel system that has dissolved special nuclear material in a molten salt is a paradigm shift in several respects. One major consideration is how to develop effective nuclear material controls and accounting for these novel reactor systems. The Molten Salt Reactor Experiment is one example of a critical molten salt system and provided a significant reference library based on the documented effort. But it was low thermal power that was not intended to reflect a commercial scale electricity production design. Therefore, to facilitate and assist vendors with domestic licensing considerations, research is underway to identify methods that could be used for domestic safeguards approaches for these novel reactor systems. This research presents the results of a signature analysis of data generated from three simulated scenarios using the molten salt demonstration reactor model defined in the Transient Simulation Framework of Reconfigurable Modules. Each scenario provides 1 hour isotope inventories over a 180 day period. The scenarios investigated provide test cases to examine if direct gamma-ray spectroscopy of the fuel can be used to identify changes comparing a base case (no reactivity control and fixed fission contribution) to an insertion of reactivity (10 pcm no change in fission composition) and a change in fission composition. The analysis demonstrates that monitoring the total count rate in a highly collimated high-resolution photon energy spectrum is sensitive to perturbations imposed into the reactor model. The total photon count rate changes ≈2% for the fission composition change and ≈4.5% for the reactivity insertion compared to the base case. However, both scenarios show an increase in the total photon count rate. The total photon count rate can be used to identify changes due to power (number of fissions) but not due to a change in the material undergoing fission. To distinguish between the two cases of increased power, the photon spectrum would require an intensive analysis technique. A photon peak count rate ratio analysis could be used to identify changes in the fissile material fission generation in the core through identification of a static peak that shows little variation to the source of fission and a highly varying peak. The photon peak strength will ultimately be determined by the isotope’s fission yield. A preliminary analysis investigating the coupling of the isotope’s fission yield and its concentration in the fuel salt derived from the modeling has been performed. A ratio analysis of the photon peak count rates of 140 La to 99 Mo demonstrated that the reactivity insertion creates a distinct difference in the ratio compared to the fission composition change scenario.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Tuneable phase behaviour and glass transition via polymerization-induced phase separation in crosslinked step-growth polymers

We report once limited to chain-growth polymerizations, fine control over polymerization-induced phase separation (PIPS) has recently been demonstrated in rubber-toughened thermoset materials formed through step-growth polymerizations. The domain length scales of these thermoset materials can be elegantly tuned by utilizing a binary mixture of curing agents (CAs) that individually yield disparate morphologies. Importantly, varying the composition of the binary mixture affects characteristics of the materials such as glass transition temperature and tensile behavior. Here, we establish a full phase diagram of PIPS in a rubber-toughened epoxy system tuned by a binary CA mixture to provide a robust framework of phase behaviour. X-Ray scattering in situ and post-PIPS is employed to elucidate the PIPS mechanism whereby an initial polymerization-induced compositional fluctuation causes nanoscale phase separation of rubber and epoxy components prior to local chain crosslinking and potential macrophase separation. We further demonstrate the universality of this approach by alternatively employing binary epoxy or binary rubber mixtures to achieve broad variations in morphology and glass transitions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Correlation and co-localization of QTL for stomatal density, canopy temperature, and productivity with and without drought stress in Setaria

Abstract Mechanistic modeling indicates that stomatal conductance could be reduced to improve water use efficiency (WUE) in C4 crops. Genetic variation in stomatal density and canopy temperature was evaluated in the model C4 genus, Setaria. Recombinant inbred lines (RILs) derived from a Setaria italica×Setaria viridis cross were grown with ample or limiting water supply under field conditions in Illinois. An optical profilometer was used to rapidly assess stomatal patterning, and canopy temperature was measured using infrared imaging. Stomatal density and canopy temperature were positively correlated but both were negatively correlated with total above-ground biomass. These trait relationships suggest a likely interaction between stomatal density and the other drivers of water use such as stomatal size and aperture. Multiple quantitative trait loci (QTL) were identified for stomatal density and canopy temperature, including co-located QTL on chromosomes 5 and 9. The direction of the additive effect of these QTL on chromosome 5 and 9 was in accordance with the positive phenotypic relationship between these two traits. This, along with prior experiments, suggests a common genetic architecture between stomatal patterning and WUE in controlled environments with canopy transpiration and productivity in the field, while highlighting the potential of Setaria as a model to understand the physiology and genetics of WUE in C4 species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Latent Neural ODE for Integrating Multi-Timescale Measurements in Smart Distribution Grids

Under a smart grid paradigm, there has been an increase in sensor installations to enhance situational awareness. The measurements from these sensors can be leveraged for real-time monitoring, control, and protection. However, these measurements are typically irregularly sampled. These measure-ments may also be intermittent due to communication bandwidth limitations. To tackle this problem, this paper proposes a novel latent neural ordinary differential equations (LODE) approach to aggregate the unevenly sampled multivariate time-series measurements. The proposed approach is flexible in performing both imputations and predictions while being computationally efficient. Simulation results on IEEE 37 bus test systems illustrate the efficiency of the proposed approach.

multi time-scale measurements↗
DOE OSTI2023

Seasonal and tidal variations in hydrologic inputs drive salt marsh porewater nitrate dynamics

Salt marshes remove terrestrially derived nutrients en route to coasts. While these systems play a critical role in improving water quality, we still have a limited understanding of the spatiotemporal variability of biogeochemically reactive solutes and processes within salt marshes. We implemented a high-frequency sampling system to monitor sub-hourly nitrate (NO$^{–}_{3}$) concentrations in salt marsh porewater at Elkhorn Slough in central California, USA. We instrumented three marsh positions along an elevation gradient subjected to different amounts of tidal inundation, which we predicted would lead to varied biogeochemical characteristics and hydrological interactions. At each marsh position, we continuously monitored porewater NO$^{–}_{3}$ concentrations at depths of 10, 30, and 50 cm and porewater levels measured at 70 cm depth over seven deployments of ~10 days each that spanned seasonal wet/dry periods common to Mediterranean climates. We quantified tidal event hysteresis between NO$^{–}_{3}$ and water level to understand how NO$^{–}_{3}$ concentrations and sources fluctuate across tidal cycles. In dry periods, the NO$^{–}_{3}$-porewater level relationship indicated that the NO$^{–}_{3}$ source was likely estuarine surface water that flooded the transect during high tides and the salt marsh was a NO$^{–}_{3}$ sink. In wet periods, the NO$^{–}_{3}$-porewater level relationship suggested the salt marsh was a source of NO$^{–}_{3}$. Finally, these findings suggest that tidal and seasonal hydrologic fluxes together control NO$^{–}_{3}$ porewater dynamics and export and influence ecological processes in coastal environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Elastic strain engineering of lattice thermal conductivity of silicon: An ab-initio study

Silicon (Si) is the most essential material in the semiconductor industry. It is important to manage the thermal properties of crystalline Si. Elastic strain engineering (ESE) has proven to be an effective tool in controlling the electrical conductivity of Si in strained-silicon technology; its effects on the thermal conductivity of silicon, therefore, warrants careful investigation. The ESE effect is much more pronounced for nanostructured materials due to the ultralarge elastic strains (on the order of 10%) achievable at the nanoscale. In this work, the lattice thermal conductivity (κ L ) of Si under hydrostatic, biaxial, and uniaxial strain states is studied with ab-initio simulations, and the values of strain-dependent κ L compare well with experimental results and existing molecular dynamics simulations. To understand the mechanisms of strain-modulated κ L , the phonon bands, scattering rate, and Grüneisen parameters of phonon modes are computed. It is shown that strain can significantly change the anharmonicity of the crystal system, thus changing phonon scattering rates and κ L . Our results demonstrate that ESE can reduce silicon κ L by up to approximately 90%. Furthermore, uniaxial and biaxial strains can induce highly anisotropic thermal conductivity in Si, with relative variations up to 58.5% and 14.5%, respectively.

Anisotropy in Thermal Conductivity↗