Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Adaptive sampling”

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,315 records · Page 73

Searching for fat tails in CRISPR-Cas systems: Data analysis and mathematical modeling

Understanding CRISPR-Cas systems—the adaptive defence mechanism that about half of bacterial species and most of archaea use to neutralise viral attacks—is important for explaining the biodiversity observed in the microbial world as well as for editing animal and plant genomes effectively. The CRISPR-Cas system learns from previous viral infections and integrates small pieces from phage genomes called spacers into the microbial genome. The resulting library of spacers collected in CRISPR arrays is then compared with the DNA of potential invaders. One of the most intriguing and least well understood questions about CRISPR-Cas systems is the distribution of spacers across the microbial population. Here, using empirical data, we show that the global distribution of spacer numbers in CRISPR arrays across multiple biomes worldwide typically exhibits scale-invariant power law behaviour, and the standard deviation is greater than the sample mean. We develop a mathematical model of spacer loss and acquisition dynamics which fits observed data from almost four thousand metagenomes well. In analogy to the classical ‘rich-get-richer’ mechanism of power law emergence, the rate of spacer acquisition is proportional to the CRISPR array size, which allows a small proportion of CRISPRs within the population to possess a significant number of spacers. Our study provides an alternative explanation for the rarity of all-resistant super microbes in nature and why proliferation of phages can be highly successful despite the effectiveness of CRISPR-Cas systems.

59 BASIC BIOLOGICAL SCIENCES↗

Extended Modeling of DOE Sealed Canisters with Updated Chemistry Models

Road-ready and final disposition packaging configurations for the advanced test reactor (ATR) fuel currently specifies storage within helium backfilled DOE sealed standard canisters. The aluminum cladding of the ATR fuel contains an oxyhydroxide layer of boehmite/bayerite that generates hydrogen when subjected to irradiation. Understanding the effect of this hydrogen buildup over time to important for long term storage considerations. Previous modeling efforts have built a coupled CFD-chemical model to simulate the temperature gas phase concentrations within the DOE sealed standard canisters. These models have coupled the temperature conditions to both the gas phase radiolysis chemistry and the surface chemistry associated with the oxyhydroxide layer. The previous iteration of the model utilized constant G-values for the hydrogen generation a 50-year period. This new iteration of the model utilizes new experimental data to update the hydrogen generation rate, as well as increase the simulated time to a 200-year period. Given the half-life assumed for the primary Cs-137 isotope responsible for gamma radiation in the ATR spent fuel, the 200-year period decreases the decay heat and dose rate of the spent fuel by a factor of 100, which combined with the updated chemistry substantially decreases new hydrogen generation. Continued experimental work has identified trends for aluminum surrogate samples with oxyhydroxide layers for tests done at higher dose rates. At low initial doses a fast generation rate of hydrogen occurs which starts to roll over to a lower generation as the total dose applied increase. A small-scale chemical model was built to replicate as mini-canister surrogate system at SRNL as well as for the smaller capsule tests performed at INL. A variety of chemical models to capture this effect were tested, and a back reaction of H radical absorbing onto the surface, or inhibition of the reaction by significant H 2 cover gas were both able to fit both the mini-canister and the small capsule test data. Both kinetic fits were used to generate data from the new 200-year simulation. As additional experiments continue, the kinetic fits may be adjusted to adapt to new data. However, updated data shows that hydrogen atmosphere has little effect on actual generation data, so the model was reverted to use a star-stepped G-value for low-dose and high-dose regions. For the undried fuel case, the three models differ significantly with an end concentration of 1% for back reaction, 4.2% for inhibition, and 13.1% for constant G-value for the nominal scenario. For the dried fuel case, the nominal cases showed end concentrations of 0.23% for back reaction, 1.9% for inhibition, and 4.2% for const G-value for the nominal scenario. For the constant G-value case that is less conservative than the other two, the pressure for undried fuel increases to 1.94 atm over 200 years for the nominal case and 2.11 atm for the high decay heat case. The total hydrogen concentration after 200 years is 4.92% for the low decay heat case and 27.4% for the high decay heat case for undried fuel and is 1.4% and 9.6% for low and high decay heat fuel for the dried fuel case. If small amount of residual air is present, the potential for nitric acid formation of 595, 1568, and 2816 ppm for the lower, nominal, and upper fuel decay heat can occur. At long-term timeframes no significant shift in major species present occurs, so no appreciable amount of oxygen is present in the system. The continued rate of hydrogen generation in the canister occurs at an increase of 0.04% by mole over the final 10 years. Increasing the model range out to 1000 years continues to drop the hydrogen generation rate through decreasing dose rate.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Gravity, Calcium, and Bone: Update, 1989

Some of the results of recent short-term flights and ground-based experiments that have contributed new insights into skeletal adaptation, calcium metabolism, and growth processes in 0 g, are highlighted. After 6 months in space, bone demineralization, invariably involving the os calcis, was found not to extend to the lumbar spine in 4 exercising cosmonauts. A flight experiment in the Space Shuttle crew has documented the early events in the calcium endocrine system during spaceflight. On the ground, brief and long-term bed rest studies of healthy volunteers in the head-down tile (HDT) model of weightlessness were completed. The skeleton of the adult male responds more rapidly to unloading than previously recognized. Regional changes in bone density can be quantified in only 30 days, are highly individual, and follow the direction of gravitational forces in the HDT model during inactivity. Bone biopsy results in healthy volunteers after bed rest differ from results in paraplegics from the same sampling site. Flight experiments in growing rats reveal changes in the composition of bone mineral and matrix in the femur postflight that were found to be highly regional and suggestive of an effect of gravity on mineral distribution. These observations may be relevant to the results from an earlier Cosmos flight where artificial gravity in space was found to maintain bone strength, but not to correct the radial growth deficit.

Arnaud, Sara B.↗

Characterizing the Joint Effect of Diverse Test-Statistic Correlation Structures and Effect Size on False Discovery Rates in a Multiple-Comparison Study of Many Outcome Measures

In their 2009 Annals of Statistics paper, Gavrilov, Benjamini, and Sarkar report the results of a simulation assessing the robustness of their adaptive step-down procedure (GBS) for controlling the false discovery rate (FDR) when normally distributed test statistics are serially correlated. In this study we extend the investigation to the case of multiple comparisons involving correlated non-central t-statistics, in particular when several treatments or time periods are being compared to a control in a repeated-measures design with many dependent outcome measures. In addition, we consider several dependence structures other than serial correlation and illustrate how the FDR depends on the interaction between effect size and the type of correlation structure as indexed by Foerstner s distance metric from an identity. The relationship between the correlation matrix R of the original dependent variables and R~, the correlation matrix of associated t-statistics is also studied. In general R~ depends not only on R, but also on sample size and the signed effect sizes for the multiple comparisons.

Feiveson, Alan H.↗

Development of an Unstructured Mesh Code for Flows About Complete Vehicles

This report describes the research work undertaken at the Massachusetts Institute of Technology, under NASA Research Grant NAG4-157. The aim of this research is to identify effective algorithms and methodologies for the efficient and routine solution of flow simulations about complete vehicle configurations. For over ten years we have received support from NASA to develop unstructured mesh methods for Computational Fluid Dynamics. As a result of this effort a methodology based on the use of unstructured adapted meshes of tetrahedra and finite volume flow solvers has been developed. A number of gridding algorithms, flow solvers, and adaptive strategies have been proposed. The most successful algorithms developed from the basis of the unstructured mesh system FELISA. The FELISA system has been extensively for the analysis of transonic and hypersonic flows about complete vehicle configurations. The system is highly automatic and allows for the routine aerodynamic analysis of complex configurations starting from CAD data. The code has been parallelized and utilizes efficient solution algorithms. For hypersonic flows, a version of the code which incorporates real gas effects, has been produced. The FELISA system is also a component of the STARS aeroservoelastic system developed at NASA Dryden. One of the latest developments before the start of this grant was to extend the system to include viscous effects. This required the development of viscous generators, capable of generating the anisotropic grids required to represent boundary layers, and viscous flow solvers. We show some sample hypersonic viscous computations using the developed viscous generators and solvers. Although this initial results were encouraging it became apparent that in order to develop a fully functional capability for viscous flows, several advances in solution accuracy, robustness and efficiency were required. In this grant we set out to investigate some novel methodologies that could lead to the required improvements. In particular we focused on two fronts: (1) finite element methods and (2) iterative algebraic multigrid solution techniques.

Peraire, Jaime↗

Mechanisms for the Crystallization of ZBLAN

A number of research teams have observed that glass forming melts that are solidified in low-g exhibit enhanced glass formation. This project will examine one of these glasses, the heavy metal fluoride glass ZBLAN. A four year ground based research program has been approved to examine the crystallization of ZBLAN glasses with the purpose of testing a theory for the crystallization of ZBLAN glass. The theory could explain the general observations of enhanced glass formation of other glasses melted and solidified in low-g. Fluid flow in 1-g results from buoyancy forces and surface tension driven convection. This fluid flow can introduce shear in undercooled liquids in 1-g. In low-g it is known that fluid flows are greatly reduced so that the shear rate in fluids in low-g are extremely low. It is believed that fluids may have some weak structure in the absence of flow. Even very small shear rates could cause this structure to collapse in response to the shear. A general result would be shear thinning of the fluid. The hypothesis of this research is that: Shear thinning in undercooled liquids increases the rate of nucleation and crystallization of glass forming melts. Shear of the melt can be reduced in low-g enhancing undercooling and glass formation. Samples will be melted and quenched in 1-g under quiescent conditions at a number of controlled cooling rates to determine times and temperatures of crystallization and heated at controlled heating rates to determine kinetic crystallization parameters. Experiments will also be performed on the materials while under controlled vibration conditions and compared with the quiescent experiments in order to evaluate the effect of shear in the liquid on crystallization kinetics. After the experimental parameters are well known, experiments will be repeated under low-g (and 2-g) conditions on the KC-135 aircraft during low-g parabolic maneuvers. The results will determine the effects of shear on crystallization. Our experimental setups will be designed with low-g experiments in mind and will be tested as breadboard low-g experiments. It is very likely that the thermal analysis instrumentation can be adapted to be run in the microgravity glovebox facilities. Critical space experiments may result to test the theory at longer low-g time experiments in space.

Ethridge, Edwin C.↗

Behavioral and Physiological Analyses of Parturition In Pregnant Rats: Insights Derived from Intrauterine Telemetry

During labor and birth, fetuses are exposed to considerable physical stimulation associated with labor contractions and expulsion from the womb These forces are important for the neonates' adaptation to tile extrauterine environment. To further our understanding of the relationship between labor and postpartum outcome, we developed a novel method for measuring intrauterine pressure (IUP) in freely-moving, late pregnant and parturient rats that enables us to make precise, reliable measures of the forces experienced by rat fetuses during parturition. A small (1.25 x 4 cm) telemetric blood pressure sensor was fitted within a fluid-filled balloon, similar in size to a full term rat fetus. On Gestational day (G) 19 of the rats' 22/23 day pregnancy, each dam was anesthetized and a balloon/sensor unit surgically implanted within the uterus following removal of two fetuses. Comparisons were made between sensor-implanted dams (IMPL) and a control conditions: 1) LAP-R, laparotomy with two fetuses removed or 2) LAP-NR, laparotomy with no fetuses removed. IUP signals were sampled at 10s intervals from the IMPL dams during labor and birth. Dams in all three conditions were videorecorded enabling us to analyze the effect of the implant on behavioral expressions of parturition. Contraction frequency, duration, pup-to-pup birth intervals and pup-oriented activities of the dams measured from one hour prior to the first pup birth until the birth of the third pup were unaffected by the sensor implant. Intrauterine telemetry of freely-moving dams offers significant advantages over conventional hardwired IUP measurement techniques. These findings establish and validate intrauterine telemetry as a reliable, non-invasive technique for quantifying pressures associated with parturition.

Villareal, J.↗

Hexagonal Lattices of HIV Capsid Proteins Explored by Simulations Based on a Thermodynamically Consistent Model

HIV capsid proteins (CAs) may self-assemble into a variety of shapes under in vivo and in vitro conditions. Here, we employed simulations based on a residue-level coarse-grained (CG) model with full conformational flexibility to investigate hexagonal lattices, which are the underlying structural pattern for CA aggregations. Facilitated by enhanced sampling simulations to rigorously calculate CA dimerization and polymerization affinities, we calibrated our model to reproduce the experimentally measured affinities. Here using the calibrated model, we performed unbiased simulations on several large systems consisting of 1512 CA subunits, allowing reversible binding and unbinding of the CAs in a thermodynamically consistent manner. In one simulation, a preassembled hexagonal CA sheet developed spontaneous curvatures reminiscent of those observed in experiments, and the edges of the sheet exhibited local curvatures larger than those of the interior. In other simulations starting with randomly distributed CAs at different concentrations, existing CA assemblies grew by binding free capsomeres to the edges and by merging with other assemblies. At high CA concentrations, rapid establishment of predominant aggregates was followed by much slower adjustments toward more regular hexagonal lattices, with increasing numbers of intact CA hexamers and pentamers being formed. Our approach of adapting a general CG model to specific systems by using experimental binding data represents a practical and effective strategy for simulating and elucidating intricate protein aggregations.

59 BASIC BIOLOGICAL SCIENCES↗

Musculoskeletal Modeling Component of the NASA Digital Astronaut Project

The NASA Digital Astronaut Project s (DAP) objective is to provide computational tools that support research of the physiological response to low gravity environments and analyses of how changes cause health and safety risks to the astronauts and to the success of the mission. The spaceflight risk associated with muscle atrophy is impaired performance due to reduced muscle mass, strength and endurance. Risks of early onset of osteoporosis and bone fracture are among the spaceflight risks associated with loss of bone mineral density. METHODS: Tools under development include a neuromuscular model, a biomechanical model and a bone remodeling model. The neuromuscular model will include models of neuromuscular drive, muscle atrophy, fiber morphology and metabolic processes as a function of time in space. Human movement will be modeled with the biomechanical model, using muscle and bone model parameters at various states. The bone remodeling model will allow analysis of bone turnover, loss and adaptation. A comprehensive trade study was completed to identify the current state of the art in musculoskeletal modeling. The DAP musculoskeletal models will be developed using a combination of existing commercial software and academic research codes identified in the study, which will be modified for use in human spaceflight research. These individual models are highly dependent upon each other and will be integrated together once they reach sufficient levels of maturity. ANALYSES: The analyses performed with these models will include comparison of different countermeasure exercises for optimizing effectiveness and comparison of task requirements and the state of strength and endurance of a crew member at a particular time in a mission. DISCUSSION: The DAP musculoskeletal model has the potential to complement research conducted on spaceflight induced changes to the musculoskeletal system. It can help with hypothesis formation, identification of causative mechanisms and supplementing small data samples.

Lewandowski, B. E.↗

Hazard Assessment from Storm Tides and Rainfall on a Tidal River Estuary

Here, we report on methods and results for a model-based flood hazard assessment we have conducted for the Hudson River from New York City to Troy/Albany at the head of tide. Our recent work showed that neglecting freshwater flows leads to underestimation of peak water levels at up-river sites and neglecting stratification (typical with two-dimensional modeling) leads to underestimation all along the Hudson. As a result, we use a three-dimensional hydrodynamic model and merge streamflows and storm tides from tropical and extratropical cyclones (TCs, ETCs), as well as wet extratropical cyclone (WETC) floods (e.g. freshets, rain-on-snow events). We validate the modeled flood levels and quantify error with comparisons to 76 historical events. A Bayesian statistical method is developed for tropical cyclone streamflows using historical data and consisting in the evaluation of (1) the peak discharge and its pdf as a function of TC characteristics, and (2) the temporal trend of the hydrograph as a function of temporal evolution of the cyclone track, its intensity and the response characteristics of the specific basin. A k-nearest-neighbors method is employed to determine the hydrograph shape. Out of sample validation tests demonstrate the effectiveness of the method. Thus, the combined effects of storm surge and runoff produced by tropical cyclones hitting the New York area can be included in flood hazard assessment. Results for the upper Hudson (Albany) suggest a dominance of WETCs, for the lower Hudson (at New York Harbor) a case where ETCs are dominant for shorter return periods and TCs are more important for longer return periods (over 150 years), and for the middle-Hudson (Poughkeepsie) a mix of all three flood events types is important. However, a possible low-bias for TC flood levels is inferred from a lower importance in the assessment results, versus historical event top-20 lists, and this will be further evaluated as these preliminary methods and results are finalized. Future funded work will quantify the influences of sea level rise and flood adaptation plans (e.g. surge barriers). It would also be valuable to examine how streamflows from tropical cyclones and wet cool-season storms will change, as this factor will dominate at upriver locations.

Hazard assessment↗

A Comparative Study of Physics‐Informed and Data‐Driven Neural Networks for Compound Flood Simulation at River‐Ocean Interfaces: A Case Study of Hurricane Irene

Simulating compound flooding (CF) at the river-ocean interface within large-scale Earth System Models (ESMs) presents significant challenges due to complex interactions between river discharge, storm surge, and tides. This study assesses the comparative advantages of physics-informed and data-driven machine learning (ML) approaches for enhancing local ESM performance. We systematically compare data-driven neural network models (i.e., CNNs, U-Net, Long Short-Term Memory (LSTM), Gated Recurrent Unit), and physics-informed neural network (PINN) models, including vanilla PINN and a finite-difference-based PINN (FD-PINN). Specifically, FD-PINN is introduced to enhance computational efficiency, accelerating vanilla PINNs by ∼6.5 times while improving accuracy. To enhance data-driven model training, a new data-generation approach is developed to sample historical fluvial and coastal flood events, which ensures a robust data set for extreme event prediction. The models are evaluated using a realistic one-dimensional river domain extracted from an ESM's river mesh and the Hurricane Irene event as an independent test case. Results show that FD-PINN achieves accurate predictions with significantly reduced computational costs relative to vanilla PINNs. Among data-driven models, the best overall performance is achieved by a CNN-LSTM hybrid, which balances accuracy and efficiency. While a fully connected CNN (CNN-FC) provides the best accuracy, it incurs high computational cost. Architectures lacking strong temporal modeling tend to underperform on unseen events. These findings highlight the importance of sequence-aware designs for robust generalization. This study reveals the trade-offs between physics-informed and data-driven models and proposes an adaptive hybrid framework for integrating ML into ESMs to enhance local flood simulations.

Earth Systems Modeling↗

Post-Transcriptional Mechanisms Respond Rapidly to Ecologically Relevant Thermal Fluctuations During Temperature-Dependent Sex Determination

An organism’s ability to integrate transient environmental cues experienced during development into molecular and physiological responses forms the basis for adaptive shifts in phenotypic trajectories. During temperature-dependent sex determination (TSD), thermal cues during discrete periods in development coordinate molecular changes that ultimately dictate sexual fate and contribute to patterns of inter- and intra-sexual variation. How these mechanisms interface with dynamic thermal environments in nature remain largely unknown. By deploying thermal loggers in wild nests of the American alligator (Alligator mississippiensis) over two consecutive breeding seasons, we observed that 80% of nests exhibit both male- and female-promoting thermal cues during the thermosensitive period, and of these nests, all exhibited both male- and female-promoting temperatures within the span of a single day. These observations raise a critical question—how are opposing environmental cues integrated into sexually dimorphic transcriptional programs across short temporal scales? To address this question, alligator embryos were exposed to fluctuating temperatures based on nest thermal profiles and sampled over the course of a daily thermal fluctuation. We examined the expression dynamics of upstream genes in the temperature-sensing pathway and find that post-transcriptional alternative splicing and transcript abundance of epigenetic modifier genes JARID2 and KDM6B respond rapidly to thermal fluctuations while transcriptional changes of downstream effector genes, SOX9 and DMRT1, occur on a delayed timescale. Our findings reveal how the basic mechanisms of TSD operate in an ecologically relevant context. We present a hypothetical hierarchical model based on our findings as well as previous studies, in which temperature-sensitive alternative splicing incrementally influences the epigenetic landscape to affect the transcriptional activity of key sex-determining genes.

59 BASIC BIOLOGICAL SCIENCES↗

Indirect decentralized learning control

The new field of learning control develops controllers that learn to improve their performance at executing a given task, based on experience performing this specific task. In a previous work, the authors presented a theory of indirect learning control based on use of indirect adaptive control concepts employing simultaneous identification and control. This paper develops improved indirect learning control algorithms, and studies the use of such controllers in decentralized systems. The original motivation of the learning control field was learning in robots doing repetitive tasks such as on an assembly line. This paper starts with decentralized discrete time systems, and progresses to the robot application, modeling the robot as a time varying linear system in the neighborhood of the nominal trajectory, and using the usual robot controllers that are decentralized, treating each link as if it is independent of any coupling with other links. The basic result of the paper is to show that stability of the indirect learning controllers for all subsystems when the coupling between subsystems is turned off, assures convergence to zero tracking error of the decentralized indirect learning control of the coupled system, provided that the sample time in the digital learning controller is sufficiently short.

Longman, Richard W.↗

Construction of Polarimetric Radar-Based Reference Rain Maps for the Iowa Flood Studies Campaign

The Global Precipitation Measurement (GPM) Mission Iowa Flood Studies (IFloodS) campaign was conducted in central and northeastern Iowa during the months of April-June, 2013. Specific science objectives for IFloodS included quantification of uncertainties in satellite and ground-based estimates of precipitation, 4-D characterization of precipitation physical processes and associated parameters (e.g., size distributions, water contents, types, structure etc.), assessment of the impact of precipitation estimation uncertainty and physical processes on hydrologic predictive skill, and refinement of field observations and data analysis approaches as they pertain to future GPM integrated hydrologic validation and related field studies. In addition to field campaign archival of raw and processed satellite data (including precipitation products), key ground-based platforms such as the NASA NPOL S-band and D3R Ka/Ku-band dual-polarimetric radars, University of Iowa X-band dual-polarimetric radars, a large network of paired rain gauge platforms, and a large network of 2D Video and Parsivel disdrometers were deployed. In something of a canonical approach, the radar (NPOL in particular), gauge and disdrometer observational assets were deployed to create a consistent high-quality distributed (time and space sampling) radar-based ground "reference" rainfall dataset, with known uncertainties, that could be used for assessing the satellite-based precipitation products at a range of space/time scales. Subsequently, the impact of uncertainties in the satellite products could be evaluated relative to the ground-benchmark in coupled weather, land-surface and distributed hydrologic modeling frameworks as related to flood prediction. Relative to establishing the ground-based "benchmark", numerous avenues were pursued in the making and verification of IFloodS "reference" dual-polarimetric radar-based rain maps, and this study documents the process and results as they pertain specifically to efforts using the NPOL radar dataset. The initial portions of the "process" involved dual-polarimetric quality control procedures which employed standard phase and correlation-based approaches to removal of clutter and non-meteorological echo. Calculation of a scale-adaptive KDP was accomplished using the method of Wang and Chandrasekar (2009; J. Atmos. Oceanic Tech.). A dual-polarimetric blockage algorithm based on Lang et al. (2009; J. Atmos. Oceanic Tech.) was then implemented to correct radar reflectivity and differential reflectivity at low elevation angles. Next, hydrometeor identification algorithms were run to identify liquid and ice hydrometeors. After the quality control and data preparation steps were completed several different dual-polarimetric rain estimation algorithms were employed to estimate rainfall rates using rainfall scans collected approximately every two to three minutes throughout the campaign. These algorithms included a polarimetrically-tuned Z-R algorithm that adjusts for drop oscillations (via Bringi et al., 2004, J. Atmos. Oceanic Tech.), and several different hybrid polarimetric variable approaches, including one that made use of parameters tuned to IFloodS 2D Video Disdrometer measurements. Finally, a hybrid scan algorithm was designed to merge the rain rate estimates from multiple low level elevation angle scans (where blockages could not be appropriately corrected) in order to create individual low-level rain maps. Individual rain maps at each time step were subsequently accumulated over multiple time scales for comparison to gauge network data. The comparison results and overall error character depended strongly on rain event type, polarimetric estimator applied, and range from the radar. We will present the outcome of these comparisons and their impact on constructing composited "reference" rainfall maps at select time and space scales.

Radar↗

The GGCMI Phase 2 Experiment: Global Gridded Crop Model Simulations Under Uniform Changes in CO2, Temperature, Water, and Nitrogen Levels (Protocol Version 1.0)

Concerns about food security under climate change motivate efforts to better understand future changes in crop yields. Process-based crop models, which represent plant physiological and soil processes, are necessary tools for this purpose since they allow representing future climate and management conditions not sampled in the historical record and new locations to which cultivation may shift. However, process-based crop models differ in many critical details, and their responses to different interacting factors remain only poorly understood. The Global Gridded Crop Model Intercomparison (GGCMI) Phase 2 experiment, an activity of the Agricultural Model Intercomparison and Improvement Project (AgMIP), is designed to provide a systematic parameter sweep focused on climate change factors and their interaction with overall soil fertility, to allow both evaluating model behavior and emulating model responses in impact assessment tools. In this paper we describe the GGCMI Phase 2 experimental protocol and its simulation data archive. A total of 12 crop models simulate five crops with systematic uniform perturbations of historical climate, varying CO2, temperature, water supply, and applied nitrogen (“CTWN”) for rainfed and irrigated agriculture, and a second set of simulations represents a type of adaptation by allowing the adjustment of growing season length. We present some crop yield results to illustrate general characteristics of the simulations and potential uses of the GGCMI Phase 2 archive. For example, in cases without adaptation, modeled yields show robust decreases to warmer temperatures in almost all regions, with a nonlinear dependence that means yields in warmer baseline locations have greater temperature sensitivity. Inter-model uncertainty is qualitatively similar across all the four input dimensions but is largest in high-latitude regions where crops may be grown in the future.

GGCMI↗

Development of Direct Injection/Ionization Mass Spectrometry Methods for Whole Molecule Characterization

The objective of this work is to adapt ambient ionization mass spectrometry (AMS) techniques for the rapid analysis of intact uranium complexes, stable strontium, cerium, and explosive compounds. The methods used were “soft ionization” techniques, which facilitate the detection of whole molecule complexes. The soft ionization mass spectrometry (MS) techniques that were investigated include paper spray ionization (PSI), matrix-assisted ionization (MAI), electrospray ionization (ESI) and direct analysis in real time (DART). For the first time, PSI-MS methods were successfully developed for whole molecule uranium-containing analytes (uranyl acetate, uranyl nitrate, and uranyl-tributylphosphate complexes). This was also the first demonstration of uranium complex detection and characterization and one of the few examples of inorganic analysis using MAI techniques. Proof of concept experiments also putatively identified matrix-derived ions and ion complexes that have not previously been described in the literature. Additionally, PSI-MS on cotton swipe samples doped with a multi-element standard containing μg levels of U, Bi, Pb, Cd, Fe, and Zn were directly analyzed without purification, representing a major improvement over existing methods. Both PSI and MAI methods demonstrated limits of detection (LODs) in the 10 - 100’s ng for various uranyl species within a range of 10’s ppm - 100’s ppb, dependent on analytical method and analyte species. AMS methods were also developed for other inorganics, including Ce and Sr, and organic explosive residues to address specific challenges in environmental monitoring and forensics. Further refinement and qualification of the AMS techniques developed within this effort would lead to significant cost reduction and timeliness by facilitating the triage and queueing of samples for subsequent more sensitive and time-consuming analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Bridging the length scales on mechanical property evaluation (Final Report)

The development of small-scale mechanical testing in combination with microstructural investigation is of great interest to the nuclear materials community for both materials development and monitoring applications. Dramatically reducing the sample sizes to reduce radioactivity and obtaining mechanical properties of irradiated samples is truly intriguing. Moreover, such studies promise a range of benefits including cost reduction, fundamental insight in structure-property relationships, increased statistics on less sample material, and reinvestigation of prior irradiated and tested reactor samples while simultaneously enabling the generation of mechanical test data on ion beam irradiated materials with limited penetration depths. Small scale materials testing on sub-sized samples has been studied for several decades, though it has only been after the development of micro-testing based on Focused Ion Beam (FIB) sample manufacturing in that orders of magnitude smaller samples could really be investigated in a quantitative manner. In recent years, small scale mechanical testing techniques at a number of length scales has been developed for both unirradiated and irradiated (ion and neutron) materials. Technological advances made in this field have enabled ex-situ and in-situ transmission electron microscopy (TEM) and scanning electron microscopy (SEM) examination, thus leading to more accurate measurements as well as additional mechanistic information. A recent review of the benefits of these techniques show that these techniques are at a stage to tackle multi-scale ranges of materials investigations and can be utilized to obtain fundamental science-based understanding of nuclear materials. Considering the tremendous advances made, one can see how small-scale mechanical testing techniques combined with modeling can enable true small scale to bulk scale mechanical property correlations. However, for the engineering community to adapt this approach fully, one needs to demonstrate that a) that these techniques can produce results with high fidelity and reproducibility, b) generate engineering stress-strain data that one can utilize to understand bulk behavior, and c) generate new insight into relevant phenomena fostering the true understanding of radiation damage and microstructure in materials for nuclear applications. It is the objective of this proposal to bridge the length scale between macro- and micro- scale mechanical testing of unirradiated and irradiated materials. This involves the development and demonstration of procedures for multi-scale mechanical testing that enable high fidelity reproducibility of data and the generation bulk property data from small-scale mechanical tests. Through this, the proposal aims to enhance the confidence in the obtained data at the smaller length scales and enhance the insight provided from these techniques for bulk scale applications on both unirradiated and irradiated materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Adaptive Thresholding and Parameter Estimation for PPM

A method of adaptive setting of a threshold level for the detection of pulses in a pulse-position modulation (PPM) free-space optical communication system has been developed. In simplified terms, it is desirable to set a threshold value high enough to greatly reduce the probability (PFA as defined below) of erroneously detecting noise as signal pulses but not so high as to greatly reduce the probability (PD as defined below) of detecting any signal pulses that may be present along with noise. In the present method, the threshold level is varied with time, in response to changing conditions in the optical-communication channel, in an effort to maintain a balance between the aforesaid competing requirements. An integral part of this adaptation scheme is a scheme for estimating key parameters of the optical-communication channel in particular, parameters that describe the fading and total attenuation in the channel, and parameters that characterize spreading of pulses by atmospheric and other effects. The method can be implemented by software processing of digitized optoelectronic-detector output, and has been tested by computational simulation. In the first stage of processing by this method, the digitized values of the detector output during noise-only time slots of received PPM symbols are averaged to obtain a background level. This background level is subtracted from the detector output in the hope of reducing or eliminating the noise component in the remaining signal. (This background level should not be confused with the detection threshold, which is computed in the last stage of processing.) Next, the remaining signal - in effect, a vector of pulse samples - is normalized by dividing it by its L1 norm (in general, the L1 norm of a vector is defined as the sum of absolute magnitudes of its orthogonal components).

Arabshahi, Payman↗