Search NASA⌕ Search

SEARCH · Search NASA

Results for “generator parameter estimation”

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,261 records · Page 70

TREAT Neutronics Analysis and Design Support, Part II: Multi-SERTTA-CAL

Experiment vehicle design is necessary in preparation for Transient Reactor Test (TREAT) facility restart and the resumption of transient testing to support Accident Tolerant Fuel (ATF) characterization and other future fuels testing requirements. Currently the most mature vehicle design is the Multi-SERTTA (Static Environments Rodlet Transient Test Apparatuses), which can accommodate up to four concurrent rodlet-sized specimens under separate environmental conditions. Robust test vehicle design requires neutronics analyses to support design development, optimization of the power coupling factor (PCF) to efficiently maximize energy generation in the test fuel rodlets, and experiment safety analyses. In integral aspect of prior TREAT transient testing was the incorporation of calibration experiments to experimentally evaluate and validate test conditions in preparation of the actual fuel testing. The calibration experiment package established the test parameter conditions to support fine-tuning of the computational models to deliver the required energy deposition to the fuel samples. The calibration vehicle was designed to be as near neutronically equivalent to the experiment vehicle as possible to minimize errors between the calibration and final tests. The Multi-SERTTA-CAL vehicle was designed to serve as the calibration vehicle supporting Multi-SERTTA experimentation. Models of the Multi-SERTTA-CAL vehicle containing typical PWR-fuel rodlets were prepared and neutronics calculations were performed using MCNP6.1 with ENDF/B-VII.1 nuclear data libraries; these results were then compared against those performed for Multi-SERTTA to determine the similarity and possible design modification necessary prior to construction of these experiment vehicles. The estimated reactivity insertion worth into the TREAT core is very similar between the two vehicle designs, with the primary physical difference being a hollow Inconel tube running down the length of the calibration vehicle. Calculations of PCF indicate that on average there is a reduction of approximately 6.3 and 12.6%, respectively, for PWR fuel rodlets irradiated under wet and dry conditions. Changes to the primary or secondary vessel structure in the calibration vehicle can be performed to offset this discrepancy and maintain neutronic equivalency. Current possible modifications to the calibration vehicle include reduction of the primary vessel wall thickness, swapping Zircaloy-4 for stainless steel 316 in the secondary containment, or slight modification to the temperature and pressure of the water environment within the primary vessel. Removal of some of the instrumentation within the calibration vehicle can also serve to slightly increase the PCF. Future efforts include further modification and optimization of the Multi-SERTTA and Multi-SERTTA-CAL designs in preparation of actual TREAT transient testing. Experimental results from both test vehicles will be compared against calculational results and methods to provide validation and support additional neutronics analyses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

TREAT neutronics analysis and design support, Part II: Multi-SERTTA-CAL

Experiment vehicle design is necessary in preparation for Transient Reactor Test (TREAT) facility restart and the resumption of transient testing to support Accident Tolerant Fuel (ATF) characterization and other future fuels testing requirements. Currently the most mature vehicle design is the Multi-SERTTA (Static Environments Rodlet Transient Test Apparatuses), which can accommodate up to four concurrent rodlet-sized specimens under separate environmental conditions. Robust test vehicle design requires neutronics analyses to support design development, optimization of the power coupling factor (PCF) to efficiently maximize energy generation in the test fuel rodlets, and experiment safety analyses. An integral aspect of prior TREAT transient testing was the incorporation of calibration experiments to experimentally evaluate and validate test conditions in preparation of the actual fuel testing. The calibration experiment package established the test parameter conditions to support fine-tuning of the computational models to deliver the required energy deposition to the fuel samples. The calibration vehicle was designed to be as near neutronically equivalent to the experiment vehicle as possible to minimize errors between the calibration and final tests. The Multi-SERTTA-CAL vehicle was designed to serve as the calibration vehicle supporting Multi-SERTTA experimentation. Models of the Multi-SERTTA-CAL vehicle containing typical PWR-fuel rodlets were prepared and neutronics calculations were performed using MCNP6.1 with ENDF/B-VII.1 nuclear data libraries; these results were then compared against those performed for Multi-SERTTA to determine the similarity and possible design modification necessary prior to construction of these experiment vehicles. The estimated reactivity insertion worth into the TREAT core is very similar between the two vehicle designs, with the primary physical difference being a hollow Inconel tube running down the length of the calibration vehicle. Calculations of PCF indicate that on average there is a reduction of approximately 6 % for PWR fuel rodlets irradiated under both wet and dry conditions. Changes to the primary or secondary vessel structure in the calibration vehicle can be performed to offset this discrepancy and maintain neutronic equivalency. Current possible modifications to the calibration vehicle include reduction of the primary vessel wall thickness, swapping Zircaloy-4 for stainless steel 316 in the secondary containment, or slight modification to the temperature and pressure of the water environment within the primary vessel. Removal of some of the instrumentation within the calibration vehicle can also serve to slightly increase the PCF. Future efforts include further modification and optimization of the Multi-SERTTA and Multi-SERTTA-CAL designs in preparation of actual TREAT transient testing. Experimental results from both test vehicles will be compared against calculational results and methods to provide validation and support additional neutronics analyses. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Risk Assessment and Management of Uranium Mill Tailings Legacies - 17508

exploitation activities took place at 61 sites. The great majority of these mining sites are located in central-east Portugal, including the Urgeirica mine that was the most important uranium deposit exploited in the country. The mining activities at this site generated 1.6 million m{sup 3} of radioactive wastes disposed in an open air area originating the largest radioactive tailings pile of the country. The environmental remediation of the Urgeirica site was one of the main priorities of the national plan for the rehabilitation of former uranium mining sites. Site remediation, that included the mine industrial area with several support infrastructures and a processing plant, was executed between 2006 and 2008 (1. phase) and between 2008 and 2014 (2. phase). Nevertheless, the environmental rehabilitation of the Portuguese uranium legacies sites is not yet concluded. One of the last largest mining complexes, the Quinta do Bispo open pit, is scheduled to be rehabilitated only in 2020. This work presents a study to evaluate the environmental impacts from the uranium mill tailings legacies before and after the remediation works, focusing on two different approaches: a deterministic methodology for dose assessment and a probabilistic methodology with the Monte Carlo method to estimate the incremental lifetime risk. Probabilistic distributions of the input parameters were used in the risk calculation approach. The methodology for dose and risk assessments may be applied to other sites that are still under intervention and the results may be used to prioritize the remediation stages. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modern deep neural networks for Direct Normal Irradiance forecasting: A classification approach

The escalating energy demand and the adverse environmental impacts of fossil-fuel use necessitate a shift towards cleaner and renewable alternatives. Concentrated Solar Power (CSP) technology emerges as a promising solution, offering a carbon-free alternative for power generation. The efficiency and profitability of CSP depend on the Direct Normal Irradiance (DNI) component of solar radiation; hence, accurate DNI forecasting can help optimize CSP plants’ operations and performance. The unpredictable nature of weather phenomena, particularly cloud cover, introduces uncertainty into DNI projections. Existing DNI forecasting models use meteorological factors, which are both challenging to estimate numerically over short prediction windows and expensive to model through data at a sufficiently high spatial and temporal resolution. This research addresses the challenge by presenting a novel approach that formulates DNI prediction as a multi-class classification problem, departing from conventional regression-based methods. The primary objective of this classification framework is to identify optimal periods aligning with specific operational thresholds for CSP plants, contributing to enhanced dispatch optimization strategies. We model the DNI classification problem using four advanced deep neural networks – rectified linear unit (ReLU) networks, 1D residual networks (ResNets), bidirectional long short-term memory (BiLSTM) networks, and transformers – achieving accuracies up to 93.5% without requiring meteorological parameters.

14 SOLAR ENERGY↗

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗

Characterization of the porous structures of the green body and sintered biomedical titanium scaffolds with micro-computed tomography

The present research was aimed at gaining an understanding of the porous structure changes from the green body through water leaching and sintering to titanium scaffolds. Micro-computed tomography (micro-CT) was performed to generate 3D models of titanium scaffold preforms containing carbamide space-holding particles and sintered scaffolds containing macro- and micro-pores. The porosity values and structural parameters were determined by means of image analysis. The result showed that the porosity values, macro-pore sizes, connectivity densities and specific surface areas of the titanium scaffolds sintered at 1200 °C for 3 h did not significantly deviate from those of the green structures with various volume fractions of the space holder. Titanium scaffolds with a maximum specific surface area could be produced with an addition of 60–65 vol% carbamide particles to the matrix powder. The connectivity of pores inside the scaffold increased with rising volume fraction of the space holder. The shrinkage of the scaffolds prepared with > 50 vol% carbamide space holder, occurring during sintering, was caused by the reductions of macro-pore sizes and micro-pore sizes as well as the thickness of struts. In conclusion, the final porous structural characteristics of titanium scaffolds could be estimated from those of the green body. - Highlights: •Porous structures of green body and sintered titanium scaffolds was studied. •Porous structures of both samples were quantitatively characterized with micro-CT. •Porous structures of scaffolds could be controlled from the green body. •Shrinkage mechanisms of titanium scaffolds during sintering was established.

36 MATERIALS SCIENCE↗

Production efficiencies of sound waves in the intracluster medium driven by AGN jets

Feedback from active galactic nuclei (AGNs) is believed to be the most promising solution to the cooling flow problem in cool-core clusters. Dissipation of sound waves is considered as one of the possible heating mechanisms; however, its relative contribution to heating remains unclear. To estimate the energy budget for heating, we perform 3D hydrodynamic simulations of AGN jet injections in a Perseus-like cluster and quantify the amount of energy stored in the forms of weak shocks and waves. For this work, e find that, for a single jet injection with typical parameters in cool-core clusters, ${\sim}9{{\ \rm per\ cent}}$ of the total jet energy is stored in compressional waves (including both shocks and waves). However, due to the destructive effects among randomly phased waves as well as the dissipation of shock energies, in our simulations including self-regulated AGN feedback, no more than $3{{\ \rm per\ cent}}$ of the total injected energy goes into compressional waves. We further separate the energy contribution from shocks and waves and find that, for a single outburst, the shocks can only contribute to ${\sim}20{-}30{{\ \rm per\ cent}}$ of the total compressional energy in the inner radii and quickly dissipate away. In the self-regulated case where shocks are repeatedly generated, shocks completely dominate over sound waves in the inner region and can still provide ${\sim}40{-}50{{\ \rm per\ cent}}$ of the total compressional energy at outer radii. Our results suggest that the production of sound waves is not as efficient as what was previously found, and thus sound wave dissipation may be a subdominant source of heating in cool-core clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

Adiabatic Gravitational Waveform Model for Compact Objects Undergoing Quasi-Circular Inspirals Into Rotating Massive Black Holes

We present bhpwave: a new python-based, open-source tool for generating the gravitational waveforms of stellar-mass compact objects undergoing quasicircular inspirals into rotating massive black holes. These binaries, known as extreme-mass-ratio inspirals (EMRIs), are exciting mHz gravitational wave sources for future space-based detectors such as the Laser Interferometer Space Antenna (LISA). Relativistic models of EMRI gravitational wave signals are necessary to unlock the full scientific potential of mHz detectors, yet few open-source EMRI waveform models exist. Thus we built bhpwave, which uses the adiabatic approximation from black hole perturbation theory to rapidly construct gravitational waveforms based on the leading-order inspiral dynamics of the binary. In this work, we present the theoretical and numerical foundations underpinning bhpwave. We also demonstrate how bhpwave can be used to assess the impact of EMRI modeling errors on LISA gravitational wave data analysis. In particular, we find that for retrograde orbits and slowly spinning black holes we can mismodel the gravitational wave phasing by as much as ∼10 radians without significantly biasing EMRI parameter estimation.

Zachary Nasipak↗

Volatile Organic Compounds and Meteorological Conditions in the Missouri Ozark AmeriFlux (MOFLUX) Site, 2023

This data set contains measurements of atmospheric components and meteorological conditions of central Missouri, United States during the summer of 2023. Data serve to examine the impact of meteorological conditions relevant to future climate on the emission and transformation of volatile organic compounds (VOCs). During the field campaign, researchers were also able to incorporate opportunistic analyses of the long-range transport of smoke plumes generated from extreme forest fire activities in Canada. VOC measurements were conducted at the Missouri Ozark AmeriFlux (MOFLUX) site (latitude 38.7441, longitude −92.2000) using a proton transfer reaction time of flight mass spectrometer (PTR-ToF-MS 6000 X2). The sampling campaign was conducted during the summer of 2023 (2023-06-25 to 2023-08-12) with measurements being taken at high temporal resolution (1 hour). Particular VOCs analyzed include: Methanol, Acetonitrile, Acetone, Isoprene, Methylvinyl ketone (MVK) and methacrolein (MACr), Benzene, Toluene, Catechol, and Monoterpene. Meteorological parameters included in this dataset were collected from a nearby Columbia Regional Airport (~10 km). Global solar radiation data were measured at a weather site in Ashland, MO, 5.22 km from the MOFLUX tower. The data were accessed using the MesoWest online website (https://mesowest.utah.edu/) provided by the Department of Atmospheric Sciences, University of Utah. Discussion of the data providers, database, and dissemination were highlighted in prior studies (Horel et al., 2002a; Horel et al., 2002b) . Smoke mixing ratios (in mg m−3) were estimated from the High-Resolution Rapid Refresh (HRRR) 3 km weather model for Missouri at 6-hour intervals (Dowell et al., 2022). This dataset contains three data files in comma separate (*.csv) format. Additional metadata are provided: three data dictionaries and a file-level metadata file in comma separate (*.csv) format and a user guide in PDF (*.pdf) format.

Acetone↗

Constraining Carbon and Nutrient Flows in Soil With Ecological Stoichiometry

Soil organic matter (SOM) is central to soil carbon (C) storage and terrestrial nutrient cycling. New data have upended the traditional model of stabilization, which held that stable SOM was mostly made of undecomposed plant molecules. We now know that microbial by-products and dead cells comprise unexpectedly large amounts of stable SOM because they can become attached to mineral surfaces or physically protected within soil aggregates. SOM models have been built to incorporate the microbial to mineral stabilization of organic matter, but now face a new challenge of accurately capturing microbial productivity and metabolism. Explicitly representing stoichiometry, the relative nutrient requirements for growth and maintenance of organisms, could provide a way forward. Stoichiometry limits SOM formation and turnover in nature, but important nutrients like nitrogen (N), phosphorus (P), and sulfur (S) are often missing from the new generation of SOM models. In this synthesis, we seek to facilitate the addition of these nutrients to SOM models by (1) reviewing the stoichiometric bias—the tendency to favor one element over another—of four key processes in the new framework of SOM cycling and (2) applying this knowledge to build a stoichiometrically explicit budget of C, N, P, and S flow through the major SOM pools. By quantifying the role of stoichiometry in SOM cycling, we discover that constraining the C:N:P:S ratio of microorganisms and SOM to specific values reduces uncertainty in C and nutrient flow as effectively as using microbial C use efficiency (CUE) parameters. We find that the value of additional constraints on stoichiometry vs. CUE varies across ecosystems, depending on how precise the available data are for that ecosystem and which biogeochemical pathways are present. Moreover, because CUE summarizes many different processes, stoichiometric measurements of key soil pools are likely to be more robust when extrapolated from soil incubations to plot or biome scale estimates. Our results suggest that measuring SOM stoichiometry should be a priority for future empirical work and that the inclusion of new nutrients in SOM models may be an effective way to improve precision.

54 ENVIRONMENTAL SCIENCES↗

Within-Host Phenotypic Evolution and the Population-Level Control of Chronic Viral Infections by Treatment and Prophylaxis

Chronic viral infections can persist for decades spanning thousands of viral generations, leading to a highly diverse population of viruses with its own complex evolutionary history. We propose an expandable mathematical framework for understanding how the emergence of genetic and phenotypic diversity affects the population-level control of those infections by both non-curative treatment and chemo-prophylactic measures. Our frameworks allows both neutral and phenotypic evolution, and we consider the specific evolution of contagiousness, resistance to therapy, and efficacy of prophylaxis. We compute both the controlled and uncontrolled, population-level basic reproduction number accounting for the within-host evolutionary process where new phenotypes emerge and are lost in infected persons, which we also extend to include both treatment and prophylactic control efforts. We used these results to discuss the conditions under which the relative efficacy of prophylactic versus therapeutic methods of control are superior. Finally, we give expressions for the endemic equilibrium of these models for certain constrained versions of the within-host evolutionary model providing a potential method for estimating within-host evolutionary parameters from population-level genetic sequence data.

sensitivity analysis↗

End-To-End Uncertainty Quantification with Analytical Derivatives for Design Under Uncertainty

Uncertainty quantification (UQ) is a rapidly growing and evolving discipline, especially within the aerospace community. Performing analysis with UQ can provide decision makers with a wealth of information about a candidate design. However, the value of UQ is fully realized when the information gained during UQ analysis is leveraged in a feedback loop of a design optimization process, often referred to as design under uncertainty. Although design under uncertainty can be a powerful risk mitigation technique, there are a number of roadblocks that prevent its implementation. Two primary factors are computational costs and added complexity of the analysis. High fidelity simulations on the order tens of uncertain variables quickly become computationally infeasible. Also, implementing UQ into an existing multidisciplinary design and optimization (MDO) process often requires extensive knowledge of the UQ methods and careful treatment of the problem formulation. The objective of this work is to address these two primary roadblocks and enable practitioners to efficiently perform design under uncertainty with limited knowledge of the UQ discipline. Methods outlined in this paper demonstrate MDO incorporating UQ into the design process, leveraging an analytic derivative tool chain through the entire optimization. The proposed approach leverages machine learning techniques to generate a differentiable confidence interval output from polynomial chaos models. This technique, coupled with the incorporation of analytical derivatives through the Polynomial Chaos Expansion (PCE) process, eliminates the need to estimate derivatives which are usually obtained from finite difference, complex step, or similar methods. Developing a differentiable confidence interval allows mixed uncertainty problems (both epistemic and aleatory) to be modeled. Without such modeling, these problems cannot accurately predict objective functions containing statistical quantities such as mean and variance. The addition of analytic derivatives to a polynomial chaos-based UQ method decreases the computational costs of performing design under uncertainty by orders of magnitude in comparison with methods such as complex step. The method and codes developed are modular in nature and are a drop-in solution for design under uncertainty within existing MDO problems. A low-fidelity analytical multidisciplinary optimization under uncertainty for a wing design in OpenMDAO is detailed in this paper. This demonstration case will include both objective functions and constraints which are influenced by uncertain parameters.

Ben D Phillips↗

End-To-End Uncertainty Quantification with Analytical Derivatives for Design Under Uncertainty

Uncertainty quantification (UQ) is a rapidly growing and evolving discipline, especially within the aerospace community. Performing analysis with UQ can provide decision makers with a wealth of information about a candidate design. However, the value of UQ is fully realized when the information gained during UQ analysis is leveraged in a feedback loop of a design optimization process, often referred to as design under uncertainty. Although design under uncertainty can be a powerful risk mitigation technique, there are a number of roadblocks that prevent its implementation. Two primary factors are computational costs and added complexity of the analysis. High fidelity simulations on the order tens of uncertain variables quickly become computationally infeasible. Also, implementing UQ into an existing multidisciplinary design and optimization (MDO) process often requires extensive knowledge of the UQ methods and careful treatment of the problem formulation. The objective of this work is to address these two primary roadblocks and enable practitioners to efficiently perform design under uncertainty with limited knowledge of the UQ discipline. Methods outlined in this paper demonstrate MDO incorporating UQ into the design process, leveraging an analytic derivative tool chain through the entire optimization. The proposed approach leverages machine learning techniques to generate a differentiable confidence interval output from polynomial chaos models. This technique, coupled with the incorporation of analytical derivatives through the Polynomial Chaos Expansion (PCE) process, eliminates the need to estimate derivatives which are usually obtained from finite difference, complex step, or similar methods. Developing a differentiable confidence interval allows mixed uncertainty problems (both epistemic and aleatory) to be modeled. Without such modeling, these problems cannot accurately predict objective functions containing statistical quantities such as mean and variance. The addition of analytic derivatives to a polynomial chaos-based UQ method decreases the computational costs of performing design under uncertainty by orders of magnitude in comparison with methods such as complex step. The method and codes developed are modular in nature and are a drop-in solution for design under uncertainty within existing MDO problems. A low-fidelity analytical multidisciplinary optimization under uncertainty for a wing design in OpenMDAO is detailed in this paper. This demonstration case will include both objective functions and constraints which are influenced by uncertain parameters.

Ben Phillips↗

KGML-ag: a modeling framework of knowledge-guided machine learning to simulate agroecosystems: a case study of estimating N<sub>2</sub>O emission using data from mesocosm experiments

Abstract. Agricultural nitrous oxide (N2O) emission accounts for a non-trivial fraction of global greenhouse gas (GHG) budget. To date, estimating N2O fluxes from cropland remains a challenging task because the related microbial processes (e.g., nitrification and denitrification) are controlled by complex interactions among climate, soil, plant and human activities. Existing approaches such as process-based (PB) models have well-known limitations due to insufficient representations of the processes or uncertainties of model parameters, and due to leverage recent advances in machine learning (ML) a new method is needed to unlock the “black box” to overcome its limitations such as low interpretability, out-of-sample failure and massive data demand. In this study, we developed a first-of-its-kind knowledge-guided machine learning model for agroecosystems (KGML-ag) by incorporating biogeophysical and chemical domain knowledge from an advanced PB model, ecosys, and tested it by comparing simulating daily N2O fluxes with real observed data from mesocosm experiments. The gated recurrent unit (GRU) was used as the basis to build the model structure. To optimize the model performance, we have investigated a range of ideas, including (1) using initial values of intermediate variables (IMVs) instead of time series as model input to reduce data demand; (2) building hierarchical structures to explicitly estimate IMVs for further N2O prediction; (3) using multi-task learning to balance the simultaneous training on multiple variables; and (4) pre-training with millions of synthetic data generated from ecosys and fine-tuning with mesocosm observations. Six other pure ML models were developed using the same mesocosm data to serve as the benchmark for the KGML-ag model. Results show that KGML-ag did an excellent job in reproducing the mesocosm N2O fluxes (overall r2=0.81, and RMSE=3.6 mgNm-2d-1 from cross validation). Importantly, KGML-ag always outperforms the PB model and ML models in predicting N2O fluxes, especially for complex temporal dynamics and emission peaks. Besides, KGML-ag goes beyond the pure ML models by providing more interpretable predictions as well as pinpointing desired new knowledge and data to further empower the current KGML-ag. We believe the KGML-ag development in this study will stimulate a new body of research on interpretable ML for biogeochemistry and other related geoscience processes.

54 ENVIRONMENTAL SCIENCES↗

Modeling Non-homologous End Joining

Non-homologous end joining (NHEJ) is the dominant DNA double strand break (DSB) repair pathway and involves several NHEJ proteins such as Ku, DNA-PKcs, XRCC4, Ligase IV and so on. Once DSBs are generated, Ku is first recruited to the DNA end, followed by other NHEJ proteins for DNA end processing and ligation. Because of the direct ligation of break ends without the need for a homologous template, NHEJ turns out to be an error-prone but efficient repair pathway. Some mechanisms have been proposed of how the efficiency of NHEJ repair is affected. The type of DNA damage is an important factor of NHEJ repair. For instance, the length of DNA fragment may determine the recruitment efficiency of NHEJ protein such as Ku [1], or the complexity of the DNA breaks [2] is accounted for the choice of NHEJ proteins and subpathway of NHEJ repair. On the other hand, the chromatin structure also plays a role of the accessibility of NHEJ protein to the DNA damage site. In this talk, some mathematical models of NHEJ, that consist of series of biochemical reactions complying with the laws of chemical reaction (e.g. mass action, etc.), will be introduced. By mathematical and numerical analysis and parameter estimation, the models are able to capture the qualitative biological features and show good agreement with experimental data. As conclusions, from the viewpoint of modeling, how the NHEJ proteins are recruited will be first discussed for connection between the classical sequential model [4] and recently proposed two-phase model [5]. Then how the NHEJ repair pathway is affected, by the length of DNA fragment [6], the complexity of DNA damage [7] and the chromatin structure [8], will be addressed

Li, Yongfeng↗

The DESI-Lensing Mock Challenge: large-scale cosmological analysis of 3x2-pt statistics

The current generation of large galaxy surveys will test the cosmological model by combining multiple types of observational probes. Realising the statistical promise of these new datasets requires rigorous attention to all aspects of analysis including cosmological measurements, modelling, covariance and parameter likelihood. In this paper we present the results of an end-to-end simulation study designed to test the analysis pipeline for the combination of the Dark Energy Spectroscopic Instrument (DESI) Year 1 galaxy redshift dataset and separate weak gravitational lensing information from the Kilo-Degree Survey, Dark Energy Survey and Hyper-Suprime-Cam Survey. Our analysis employs the 3x2-pt correlation functions including cosmic shear and galaxy-galaxy lensing, together with the projected correlation function of the spectroscopic DESI lenses. We build realistic simulations of these datasets including galaxy halo occupation distributions, photometric redshift errors, weights, multiplicative shear calibration biases and magnification. We calculate the analytical covariance of these correlation functions including the Gaussian, noise and super-sample contributions, and show that our covariance determination agrees with estimates based on the ensemble of simulations. We use a Bayesian inference platform to demonstrate that we can recover the fiducial cosmological parameters of the simulation within the statistical error margin of the experiment, investigating the sensitivity to scale cuts. This study is the first in a sequence of papers in which we present and validate the large-scale 3x2-pt cosmological analysis of DESI-Y1.

79 ASTRONOMY AND ASTROPHYSICS↗

Three-Dimensional Aeroelastic and Aerothermoelastic Behavior in Hypersonic Flow

The aeroelastic and aerothermoelastic behavior of three-dimensional configurations in hypersonic flow regime are studied. The aeroelastic behavior of a low aspect ratio wing, representative of a fin or control surface on a generic hypersonic vehicle, is examined using third order piston theory, Euler and Navier-Stokes aerodynamics. The sensitivity of the aeroelastic behavior generated using Euler and Navier-Stokes aerodynamics to parameters governing temporal accuracy is also examined. Also, a refined aerothermoelastic model, which incorporates the heat transfer between the fluid and structure using CFD generated aerodynamic heating, is used to examine the aerothermoelastic behavior of the low aspect ratio wing in the hypersonic regime. Finally, the hypersonic aeroelastic behavior of a generic hypersonic vehicle with a lifting-body type fuselage and canted fins is studied using piston theory and Euler aerodynamics for the range of 2.5 less than or equal to M less than or equal to 28, at altitudes ranging from 10,000 feet to 80,000 feet. This analysis includes a study on optimal mesh selection for use with Euler aerodynamics. In addition to the aeroelastic and aerothermoelastic results presented, three time domain flutter identification techniques are compared, namely the moving block approach, the least squares curve fitting method, and a system identification technique using an Auto-Regressive model of the aeroelastic system. In general, the three methods agree well. The system identification technique, however, provided quick damping and frequency estimations with minimal response record length, and therefore o ers significant reductions in computational cost. In the present case, the computational cost was reduced by 75%. The aeroelastic and aerothermoelastic results presented illustrate the applicability of the CFL3D code for the hypersonic flight regime.

McNamara, Jack J.↗

Towards a Circular Economy for PET Bottles in the US

The United States generates the most plastic waste of any country. Along with that GHG emissions from the global plastic economy are expected to increase to 15% of the global carbon budget by 2050. It is imperative that plastic recycling is made a reality to reduce both plastic pollution in the environment and GHG emissions. A portfolio of end-of-life strategies must be implemented to minimize environmental impacts and retain valuable plastic material, but it is challenging to compare options that generate products with different utility and lifetime. Plastic use reduction, reuse and recycling are thus increasingly important, but making informed policy and research decisions within this space can be challenging given the diverse range of available solutions. The novel analysis framework, Plastic Parallel Pathways Platform (4P) has been equipped with consequential life cycle assessment, techno-economic analysis, and a plastic circularity indicator to estimate the greenhouse gas (GHG) emissions, circularity, and cost of polyethylene terephthalate (PET) down-cycling to lower-quality resin, closed-loop recycling to food-grade PET bottles, up-cycling to fiber-reinforced plastic (FRP), and conversion to non-plastic products (electricity, oil) on a United States economy-wide basis. Integrating system dynamics into this robust plastics model that already incorporates techno-economics, circularity, and environmental impacts will enable identification of key bottlenecks between manufacturers, waste sorters, and reclaimers that currently prevent rapid decarbonization of the plastics economy. System dynamics (SD) explore the evolution of activities and technologies based on changed macro parameters such as plastic demand and supply, market shifts, and cross-sectoral interactions. This project particularly aims to explore the interplay between, waste collection, plastic waste sorting, recycling, and manufacturing, as well as the effect of plastic bale quality and plastic reuse initiatives on the surrounding process stages. This functionality will facilitate combinatory analysis in which a portfolio of end-of-life pathways are assessed simultaneously, with the exact makeup of that portfolio affected by parameters such as technology scales, resource constraints, and waste mitigation efforts. Integrating SD with the 4P framework enables analyzing the effect of increased revenue and reinvestment into improving process efficiencies, sorting and collection quantities. Through that, market effects of increased recycled resin availability can be studied for the plastics systems model for the US. The results will help identify technical or economic bottlenecks that currently limit efforts to decarbonize the U.S. plastics economy.

circular economy↗