Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Static Program Analysis”

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

A novel transfer learning framework for sorghum biomass prediction using UAV-based remote sensing data and genetic markers

Yield for biofuel crops is measured in terms of biomass, so measurements throughout the growing season are crucial in breeding programs, yet traditionally time- and labor-consuming since they involve destructive sampling. Modern remote sensing platforms, such as unmanned aerial vehicles (UAVs), can carry multiple sensors and collect numerous phenotypic traits with efficient, non-invasive field surveys. However, modeling the complex relationships between the observed phenotypic traits and biomass remains a challenging task, as the ground reference data are very limited for each genotype in the breeding experiment. In this study, a Long Short-Term Memory (LSTM) based Recurrent Neural Network (RNN) model is proposed for sorghum biomass prediction. The architecture is designed to exploit the time series remote sensing and weather data, as well as static genotypic information. As a large number of features have been derived from the remote sensing data, feature importance analysis is conducted to identify and remove redundant features. A strategy to extract representative information from high-dimensional genetic markers is proposed. To enhance generalization and minimize the need for ground reference data, transfer learning strategies are proposed for selecting the most informative training samples from the target domain. Consequently, a pre-trained model can be refined with limited training samples. Field experiments were conducted over a sorghum breeding trial planted in multiple years with more than 600 testcross hybrids. The results show that the proposed LSTM-based RNN model can achieve high accuracies for single year prediction. Further, with the proposed transfer learning strategies, a pre-trained model can be refined with limited training samples from the target domain and predict biomass with an accuracy comparable to that from a trained-from-scratch model for both multiple experiments within a given year and across multiple years.

36 MATERIALS SCIENCE↗

Semiglobal Safety-Filtered Extremum Seeking With Unknown CBFs

We introduce a safe extremum-seeking (Safe ES) algorithm which achieves the minimization of an unknown objective function while ensuring that an unknown, yet measured, control barrier function (CBF) remains above an arbitrarily small negative value for all time. In other words, “practical safety” is maintained during the entire period of convergence to the constrained extremum. Our design is based on quadratic program (QP) CBF style filters for safety, which is applied in an average and estimated sense. Using nonsmooth analysis tools, we guarantee semiglobal practical asymptotic (SPA) stability of the global constrained optimum, practical convergence to the safe set if starting in a condition violating the CBF, and practical safety for all time—semiglobally—if starting in safe set. The safety result of the paper is analogous with modern notions of SPA stability, guaranteeing that, for any small violation of safety, there exist design coefficients which guarantee that such a small violation is not exceeded. The paper outlines a set of sufficient conditions on the barrier and objective functions, and by way of a Lyapunov argument, we demonstrate that nonconvex constrained optimization problems can be solved. We present these results in the setting of a static map and a dynamical system. A simulation example illustrates the results.

97 MATHEMATICS AND COMPUTING↗

FY23 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY23. As a continuation of the modeling effort from FY22, this year’s work focuses on improving the RELAP5-3D model developed previously for two-phase flow simulations. The RELAP5-3D model is updated to more accurately capture the heat loss experienced by the facility. The updated model is compared against experimental data for benchmarking purposes of the RELAP5-3D input model. By correctly accounting for heat loss, the updated RELAP5-3D model can now predict the two-phase baseline case more accurately. The onset and the duration of instability are captured well by the model. Furthermore, analyses are performed to better understand the instability mechanism experienced by the flow where the expansion of the boiling boundary in the chimney is studied in details and the fundamental frequencies of the oscillations are obtained. The updated RELAP5-3D model is further compared against four fault conditions, namely the reduction of riser header inlet flow area, depletion of system inventory, blocked riser channels, and static boiling scenario. For each fault condition, minor modifications and tuning are performed to improve the predictions of the model. The purpose of the analyses is to investigate the capability of RELAP5-3D in predicting complex two-phase flows in possible accident scenarios in actual Reactor Cavity Cooling System (RCCS). Overall, the model is able to capture the behaviors and trends of these fault conditions relatively well. Some discrepancies remain between the experimental data and the predictions, many of which are likely due to the differences in the predicted and experimental vapor generation rate. Future work will focus on the continued development of the current RELAP5-3D input model of the NSTF to both improve the accuracy of the model’s predictive capability and continue supporting the experimental program needs. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Field-Coupled Water Splitting with Metal-Free Donor–Acceptor Covalent Organic-Framework Junctions

Advancing metal-free electrocatalysts for hydrogen and oxygen evolution reactions (HER/OER) across acidic and alkaline media requires coordinated control of intermediate binding thermodynamics, interfacial charge delivery, and near-electrode transport dynamics. Here, we design amide-linked benzene–triazine covalent organic frameworks (BTA/TzTA-Hz COFs) and integrate them with carbon nanotubes (CNTs) to form COF–CNT junctions that establish a built-in interfacial electric field. Density functional theory (DFT) and electrostatic potential maps indicate complementary active motifs, with benzene-proximal fragments associated with HER and triazine-proximal motifs associated with OER. CNT integration shifts the contact-potential difference by ≈0.20 V, while operando electrochemical impedance spectroscopy suggests partially separable high-frequency junction-charging and lower-frequency Faradaic/transport responses. A 300 mT static magnetic field lowers the HER and OER overpotentials by tens of millivolts. Under anodic bias, the effective interfacial charging capacitance increases, and Mott–Schottky analysis shows an apparent ∼0.15 V flat-band shift with an essentially unchanged slope. Together, these observations are consistent with field-perturbed interfacial charging and altered bias partitioning. Field-dependent impedance and bubble imaging are consistent with magnetohydrodynamic convection that promotes bubble detachment and near-electrode mass transport for both half-reactions, and they reveal an OER-specific high-frequency perturbation under anodic bias. Under field, the heterostructure reaches an OER onset overpotential of ∼261 mV and requires an overpotential of 366 mV at 10 mA cm –2 in alkaline electrolyte. These results illustrate how reticular-framework chemistry, junction engineering, and both built-in and applied fields can program reactivity through interfacial electrostatics and near-electrode transport in organic-framework electrocatalysts.

Garcia-Enriquez, Lissette [Univ. of Texas at El Pa↗

A large-scale analysis of bioinformatics code on GitHub

In recent years, the explosion of genomic data and bioinformatic tools has been accompanied by a growing conversation around reproducibility of results and usability of software. However, the actual state of the body of bioinformatics software remains largely unknown. Here, the purpose of this paper is to investigate the state of source code in the bioinformatics community, specifically looking at relationships between code properties, development activity, developer communities, and software impact. To investigate these issues, we curated a list of 1,720 bioinformatics repositories on GitHub through their mention in peer-reviewed bioinformatics articles. Additionally, we included 23 high-profile repositories identified by their popularity in an online bioinformatics forum. We analyzed repository metadata, source code, development activity, and team dynamics using data made available publicly through the GitHub API, as well as article metadata. We found key relationships within our dataset, including: certain scientific topics are associated with more active code development and higher community interest in the repository; most of the code in the main dataset is written in dynamically typed languages, while most of the code in the high-profile set is statically typed; developer team size is associated with community engagement and high-profile repositories have larger teams; the proportion of female contributors decreases for high-profile repositories and with seniority level in author lists; and, multiple measures of project impact are associated with the simple variable of whether the code was modified at all after paper publication. In addition to providing the first large-scale analysis of bioinformatics code to our knowledge, our work will enable future analysis through publicly available data, code, and methods.

97 MATHEMATICS AND COMPUTING↗

Controlling Molecular Structure and Spin with Multiconfigurational Quantum Chemistry (Final Technical Report)

For many first-row transition metal complexes, structure-property relationships can be obtained from high-level molecular geometry optimizations and subsequent electronic structure studies. The project funded under this award utilized newly implemented fully internally contracted (FIC) nuclear gradients for extended multi-state (XMS) complete active space second-order multireference perturbation theory (CASPT2) to explore geometry changes in first-row transition metal coordination complexes. At the start of the project, only one full geometry optimization using FIC-CASPT2 analytical gradients had been reported (J. Chem. Theory Comput., 2016, 12 (8), 3781), and many open-questions regarding the performance and achievable accuracy in applying such computations to larger complexes persisted. Key deliverables in the report include the implementation of the numerical Hessian and subsequent vibrational analysis in the BAGEL program package. This includes both full Hessian vibrational analysis (FHVA) and a partial Hessian vibrational analysis (PHVA). The latter of which allows us to reduce the significant cost of computing the vibrational frequencies in a large molecule by focusing on the modes of highest interest. The gradient code also proved useful in the development of an approach to improve the Hubbard-U correction by removing bias in predicting spin-splitting in Fe(II) complexes via plane-wave density functional theory (DFT). Finally, we showed for three families of complexes (spin-crossover (SCO) complexes, metallocorroles complexes, and those with metal-metal bonds) that CASPT2 can result in good molecular geometries and established best practices in undertaking this work. We are now using this approach in collaborative projects with experimental groups, which would not have been possible at the start of this project. A common theme also arose through this work that static and dynamic correlation must be recovered in a balanced way to yield quantitative results. By systematically varying how both effects were included, we were able to make key chemical insights. This work supported the training of two postdoctoral scholars, one graduate student, and one undergraduate student in applying high-level wavefunction based methods to challenging systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗

Adaptive elasticity policies for staging-based in situ visualization

In situ processing aims to alleviate the growing gap between computation and I/O capabilities by performing data processing close to the data source. In situ processing is widely used to process data generated by multiple data sources, including observation data from edge devices or scientific observational facilities and the simulation data generated by scientific computation on a high-performance computing (HPC) platform. For a scientific workflow that is run on an HPC platform and composed of a simulation program and an in situ data analytics or visualization (abbreviated as ana/vis) task, there is an implicit assumption that the computing resources assigned to the workflow keep static during the workflow execution. However, with the converging trend between the HPC and cloud computing platform, running the in situ ana/vis task in an elastic way is promising to decrease its overhead and improve its resource utilization rate. Resource elasticity represents the ability to change resource configurations such as the number of computing nodes/processes during workflow execution. An elastic job may dynamically adjust resource configurations; it may use a few resources at the beginning and more resources toward the end of the job when interesting data appear. However, it is hard to predict a priori how many computing nodes/processes need to be added/removed during the workflow execution to adapt to changing workflow needs. How to efficiently guide elasticity operations, such as growing or shrinking the number of processes used for in situ analysis during workflow execution, is an open-ended research question. In this article, we present adaptive elasticity policies that adopt workflow runtime information collected during workflow execution to predict how to trigger the addition/removal of processes in order to minimize in situ processing overhead. Taking in situ visualization tasks as an example, we integrate the presented elasticity policies into a staging-based elastic workflow and evaluate its efficiency in multiple elasticity scenarios. Compared with the situation without elasticity or with a static elasticity policy that uses a fixed number of processes for each rescaling operation, the adaptive elasticity policy can save overhead in finding a proper resource configuration and improve resource utilization efficiency. Furthermore, one experiment illustrates that the adaptive elasticity policy saves 41% of core-hours compared with the situation without the resource elasticity.

97 MATHEMATICS AND COMPUTING↗

Creep-fatigue Behavior and Damage Accumulation of a Candidate Structural Material for Concentrating Solar Power Solar Thermal Receiver (Final Technical Report)

Creep-fatigue deformation is an important consideration for a thermal receiver in Concentrating Solar Power systems due to the constant static stress or pressure, diurnal cycling, and elevated service temperatures required for efficient operation. An accurate description of the creep-fatigue behavior, not available for five of the six candidate materials, is important for assessment of preliminary designs. This project has provided a detailed analysis of the creep-fatigue behavior and damage accumulation of Alloy 740H, a candidate structural material for a CSP solar thermal receiver to address a critical knowledge barrier for receiver designs identified in the CSP Gen3 Demonstration Roadmap. This project was divided in to three tasks, split between two general categories: Experimental and design methodology. Tasks 1 and 3, relating to the collection of data on Alloy 740H plate and sheet forms, respectively, were performed at the Idaho National Laboratory and used as input to Task 2, the development of design methodology performed at Argonne National Laboratory. Experimental work covered limited tensile and creep testing to fill gaps in available data that is currently available and to demonstrate how the plate of Alloy 740H tested in this work compared to historical data. A more substantial testing program for fatigue and creep-fatigue was performed, with particular emphasis on creep-fatigue as there is little available data and it is necessary to understand the creep-fatigue interaction and develop related design models. This effort includes the development of rules for the design of solar receiver components against high temperature creep-fatigue and ratcheting failure modes. The ASME Code rules for high temperature nuclear components was used to the basis of the method but adjustments will be made to reflect the generally shorter, diurnal operating cycles of thermal receivers and the relative consequences of failure, comparing nuclear to solar components. This work has provided guidance for the design of components for concentrating power facilities operating at high temperatures and undergoing high, secondary thermal stresses relative to the applied primary pressure stress. The design rules were developed for the design of Generation 3 CSP tubular receivers manufactured from Alloy 740H, but are generally applicable to a wide range of component types undergoing similar loads. Procedural design rules for components to be used in conjunction with the 2019 edition of the ASME Boiler & Pressure Vessel Code and the corresponding design data for Alloy 740H have been developed. A commentary describing the rationale behind the design rules and the data underlying the design material properties, as well as an extensive set of worked sample problems detailing the application of the rules to CSP components were developed as a part of this project, with a separate report covering these aspects in detail and included in the Appendix of this report.

14 SOLAR ENERGY↗

Oak Ridge National Laboratory Annual Sustainability Report 2023

ORNL, managed under contract by UT-Battelle LLC, is DOE’s largest science and energy laboratory and, as such, executes the widest range of mission capabilities. Diverse expertise spans a broad range of scientific and engineering disciplines, enabling research and science achievements to accelerate the delivery of solutions to the marketplace. ORNL supports DOE’s national missions of scientific discovery, clean energy, and security. To execute these activities, ORNL has grown significantly over 80 years of continuous operations, consisting of facilities with commissioning dates ranging from the 1940s to the present—an extraordinary set of distinctive scientific facilities and equipment. The complexities of such a variety of facilities require teamwork among divisions, a wide variety of conservation projects, and creative strategies to achieve the desired energy and water savings. Such a diverse and unique set of major facilities, totaling over 5.5 million square feet, with 6,000 employees, requires an innovative plan to accomplish advancements in operational efficiencies. ORNL is tasked with the management of an extraordinary set of distinctive scientific facilities and equipment for DOE. ORNL is mission-driven, and its mission has grown substantially over the decades. ORNL’s core research capabilities provide broad science and technology support for DOE in the areas of energy, environment, and national security. Currently, ORNL is a world leader in materials, neutron, and nuclear science and engineering, and in high-performance computing and data analytics. ORNL’s vast portfolio of research facilities must be maintained and carefully upgraded to protect the nation’s investment in scientific analysis. The goal of sustainable and resilient operations is to enable more effective execution of ORNL’s science and technology mission. Sustainable operational practices and enhanced resilience strive for excellent results while remaining diligent in energy conservation, environmental stewardship, asset management, and community engagement. The Sustainable ORNL Program (Sustainable ORNL) Continuous improvements in operational and business processes must be integrated into the fabric of the ORNL culture to maximize the return from the investment made in modernizing facilities and equipment. The Sustainable ORNL program promotes the legacy of system-wide best practices, management commitment, and employee engagement that will lead ORNL into a future of efficient, resilient, and sustainable operations. ORNL leadership and Sustainable ORNL champions receive regular status reports on the progress of each project and focus area (i.e., roadmap) and periodic summary reports. More information can be found at the program’s website. The Sustainable ORNL roadmap structure endorses 15 vital roadmaps. The figure below summarizes the current project assignments and demonstrates that each project contributes to the wellbeing of the whole. Continuous employee engagement and regular status reports confirm the ideals of the program. The roadmap structure is not static; as the science mission advances and the needs of the organization evolve, the Sustainable ORNL roadmap structure elements are modified to align with developing priorities. In 2022, Sustainable ORNL made roadmap changes to better align ORNL to support new federal requirements that have been issued.

54 ENVIRONMENTAL SCIENCES↗

Novel experimental probes of QCD in SIDIS and e + e - annihilation (Final Technical Report)

The research addressed with this award seeks to advance our understanding of the structure and dynamics underlying the properties of visible matter. In our current understanding the nucleons (protons and neutrons) are not fundamental but are comprised of quarks and gluons, which are collectively called partons. A static quark picture fails to explain the properties of the nucleons, such as their mass and intrinsic spin, which are thought to emerge dynamically from the quark-gluon interactions via the strong force. In this work, novel observables employing correlations of particles produced in the scattering of high energy electrons off protons at the CLAS12 experiment at Jefferson Lab were analyzed to probe quark-gluon interactions. The ultimate goal of this line of inquiry is to be able to describe the properties of protons and neutrons from first principles, similar to how studying the hydrogen atom has led to the formulation of the theory of Quantum Electrodynamics. Because quarks cannot be observed directly but only as part of more complex composite particles, a smaller, complimentary part of this work was the analysis of particle production in electron-positron annihilation at the Belle II experiment to understand the production of particles from initial quarks. This takes advantage of the fact that in e + e - annihilation the initial quark dynamics is known, unlike in the scattering of nucleons. Several new applications using Machine Learning algorithms for event tagging and reconstruction to support this program were developed as part of this award. In addition, we had a significant role in the development of the physics program for the future Electron-Ion Collider, which is a new collider to be build in the US within the next decade.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Plasma Surface Interactions: Bridging from the Surface to the Micron Frontier through Leadership Class Computing (Final Technical Progress Report)

The performance demands on plasma-facing components (PFCs) in future fusion power plants are beyond the capability of current materials. In defining the plasma-surface interaction (PSI), it is evident that three coupled spatial regions influence PFC materials evolution and performance. These regions consist of (1) the edge and scrape-off layer region of the plasma, (2) the near-surface material response to extreme thermal and particle fluxes under the influence of, and feedback to, the plasma sheath, and (3) the structural material’s response to an intense, 14 MeV peaked neutron spectrum, which produces very high concentrations of transmuted elements through (n,p) and (n,α) reactions and structural material property degradation. The objective of this project is to develop and deploy validated, high-performance simulation tools capable of predicting the performance of tungsten-based plasma-facing components (PFCs) in a burning fusion plasma environment, which includes modeling surface morphology evolution in either erosion or re-deposition regimes, and the recycling of hydrogen species. This requires the development of a leadership-scale computational code to predict PFC behavior that is well integrated to a suite of multiscale modeling techniques to bridge the scales needed to address complex physical and computational issues at the plasma surface interface and the transition region below the surface where neutron damage processes in the bulk dominate material behavior. Successful completion of this project will provide PFC simulation tools to evaluate steady-state performance of tungsten-based PFC and divertor components in burning plasma environments. This will enable the identification of critical experiments to confirm whether practical PFC solutions exist for magnetic fusion energy beyond ITER. More specifically, the research activities within this project focus on two broadly defined research thrusts: 1) Developing a new simulation code, Xolotl-PSI to predict PFC operating lifetime and performance, that is specifically designed to take advantage of leadership class computing facilities; and 2) Integrating and applying discrete particle-based, as well as continuum-based, multiscale modeling techniques to provide scientific discovery of the mechanisms controlling PFC and bulk materials evolution under fusion plasma and 14-MeV neutron exposure; in this thrust, the challenge of scale-bridging between atomistic/microstructural modeling and the continuum-scale of PFC device performance is achieved through multiscale modeling techniques. Throughout the grant period, research at UMass Amherst focused on the atomic-scale analysis of the dynamics and kinetics (transport and reactions) of small mobile helium clusters near surfaces and grain boundaries in plasma-exposed tungsten, as well as the proper parameterization and incorporation of these processes into our SciDAC team’s continuum-level cluster dynamics code (Xolotl). Our work also included characterization of large-scale molecular-dynamics (MD) simulation results of helium implantation into tungsten based on the findings of our analysis and comparisons of Xolotl simulation predictions with the large-scale MD simulation results for Xolotl validation purposes; emphasis was placed on the impact of helium cluster dynamics/kinetics on tungsten surface morphology, near-surface tungsten structure, and helium retention. Furthermore, our work included initiation of a research program on evaluating the impact of structural defects in tungsten due to plasma exposure (such as helium nanobubbles) on tungsten’s thermophysical properties, focusing on lattice thermal conductivity. Our studies contributed to a fundamental understanding of the effects of plasma-surface interactions on the dynamical response of plasma-facing materials in nuclear fusion devices and the development of predictive hierarchical multiscale computational tools for the quantitative description of the material’s dynamical response to plasma exposure.

36 MATERIALS SCIENCE↗

The LSST DESC DC2 Simulated Sky Survey

Here, we describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep area of approximately 300 deg 2 , as well as a deep drilling field of approximately 1 deg 2 . We simulate 5 yr of the planned 10 yr survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the data set to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic test bed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time domain cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Implementation and Verification of Cable Bending Stiffness in MoorDyn: Preprint

Cable bending stiffness capabilities have been added to MoorDyn to enable modeling of dynamic power cables. The relatively large motions experienced by floating wind turbines and wave energy converters pose a challenge for power cables, whose internal components provide significant bending resistance and are sensitive to deformation. The behavior and associated design considerations of power cables in these highly dynamic applications make coupled analysis relevant for design. MoorDyn is a lumped-mass mooring dynamics model that is used for floating wind energy simulation (often coupled with OpenFAST) and floating wave energy converter simulation (often coupled with WEC-Sim). Its existing line elasticity formulation considers axial stiffness only. To properly capture the dynamics of power cables, a cable bending stiffness model has been added that approximates cable curvature based on the difference in tangent vectors of adjacent line elements. The resulting bending moment is realized by applying forces on adjacent nodes, leaving the underlying lumped-mass formulation unchanged. The bending stiffness implementation is verified in static conditions against analytical solutions and then in a dynamic power cable scenario in comparison with OrcaFlex. The dynamic scenario uses prescribed motions and includes wave loadings on the cable. Results indicate correct implementation of bending stiffness and show close agreement with OrcaFlex.

50 EE - Wind and Water Power Program - Water (EE-4↗

Nevada National Security Site 2018 Data Report: Groundwater Monitoring Program, Area 5 Radioactive Waste Management Site

This report is a compilation of the groundwater and waste cell leachate sampling results from the Area 5 Radioactive Waste Management Site (RWMS) at the Nevada National Security Site, Nye County, Nevada. Groundwater samples from the aquifer immediately below the Area 5 RWMS have been collected and analyzed and static water levels have been measured in this aquifer since 1993. These groundwater data are evaluated for evidence of effects on the aquifer related to the Area 5 RWMS. Leachate from one lined mixed waste cell at the Area 5 RWMS (Cell 18) has been sampled since the cell opened in 2011. These leachate data are analyzed for hazardous contaminants to determine appropriate leachate handling and disposal. A second lined mixed waste cell at the Area 5 RWMS (Cell 25) began receiving waste in August 2018, but no leachate samples were collected or analyzed during 2018. This report includes data from 2014 through 2018 but focuses on the 2018 data. During 2018, groundwater samples were collected and static water levels were measured at three wells surrounding the Area 5 RWMS. Groundwater samples were collected at wells UE5PW-1, UE5PW-2, and UE5PW-3 on March 6 and August 14, 2018. Static water levels were measured at each of these wells on March 5, June 4, August 13, and October 25, 2018. Groundwater monitoring at the Area 5 RWMS complies with Title 40 Code of Federal Regulations (CFR) Part 264.97, General ground-water monitoring requirements , and 40 CFR 264.98, Detection monitoring program. Groundwater samples were analyzed for indicators of contamination (pH, specific conductance [SC], total organic carbon, total organic halides, and tritium) and, beginning in 2017, toxicity characteristic metals (arsenic, barium, cadmium, chromium, lead, mercury, selenium, and silver). General water chemistry parameters (calcium, iron, magnesium, manganese, potassium, sodium, sulfate, chloride, fluoride, bicarbonate, and silicate) were also measured. Results from samples collected in 2018 are within the limits for each parameter established by the Nevada Division of Environmental Protection in Resource Conservation and Recovery Act (RCRA) Permit NEV HW0101, Revision 6. Data from the shallow aquifer indicate there has been no measurable impact to the uppermost aquifer from the Area 5 RWMS, and there were no significant changes in measured groundwater parameters compared to previous years. Leachate from lined mixed waste Cell 18 drains into a sump and is pumped into an above ground tank. Samples are collected from the tank when the leachate volume approaches the 3,000-gallon tank capacity. During 2018, leachate samples were collected on March 1, 2018 and August 1, 2018. All leachate analysis results are below the regulatory levels for toxicity characteristic contaminants (40 CFR 261.24) and below the investigation levels for polychlorinated biphenyls (PCBs), tritium, pH, and SC specified in RCRA Permit NEV HW0101, Revision 6. No quantifiable PCB levels were detected in any leachate sample. Based on the leachate analysis results, RCRA Permit NEV HW0101 allowed all leachate collected from Cell 18 to be pumped from the collection tank and used for dust control at Cell 18.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Nanosecond Gated CMOS Camera (NSGCC) ICD (Rev. 2.1)

The Ultra-Fast X-ray Imager (UXI) program is an ongoing effort at Sandia National Laboratories to create high speed, multi-frame, time-gated Read Out Integrated Circuits (ROICs), and a corresponding suite of photodetectors to image a wide variety of High Energy Density (HED) physics experiments on both Sandia’s Z-Machine and LLNL’s National Ignition Facility (NIF). Several cameras have been designed over the length of the program; one of the most recent is the Icarus, which is an improvement on past imagers (Furi and Hippogriff). A second sensor that can be connected is the Daedalus sensor. The Icarus is a 1024 × 512-pixel array with either 25 μm or 8 µm spatial resolution containing four frames of storage per pixel and has improved timing generation and distribution components while achieved 2 ns time gating. The Daedalus sensor is also a 1024 x 512-pixel array with 25 µm special resolution containing three frames of storage per pixel and has an increased set of features for a wider variety of applications from interlacing of rows in each frame to configurability of all shutters. See Section 14 for details regarding the Icarus implementation of the firmware and Section 15 for details regarding the Daedalus implementation. Due to the unique test environments UXI sensors are targeted for, full custom hardware was required to physically mount an Icarus or Daedalus sensor, manage its various functions, and read out pixel data for transfer to a host computer. Beyond experimental functionality, the hardware also needed to accommodate sensor characterization requirements. Lawrence Livermore National Laboratory’s ‘Version 4.0 Board’ was the result of these efforts. It mounts all the components required to fully utilize the Icarus and Daedalus sensors including analog to digital converters to convert pixel data and various system voltages to digital form for readout and analysis, DAC channels for remote configuration of critical bias voltages, static random-access memories to buffer pixel data, RS422 and Gigabit Ethernet communications for remote access, and an FPGA to tie these components together. This document describes the FPGA electrical interfaces in detail to allow the reader a greater understanding of the device, and to facilitate implementation of custom software to control and manage it. The Version 4.0 Board is a continuation of the Nano-second Gated CMOS hardware design that retains much of the functionality of the Version 1.0 Board while adding features including a DAC instead of digital potentiometers, as well as sensors for pressure and radiation. The Version 4.0 board is intended for applications requiring tight form-factor enclosures. It is composed of two stacking boards; one holds the FPGA and regulators to power the various components of the board while the other contains the mating connector to the sensor, image-readoff ADCs, the DAC, and other components.

42 ENGINEERING↗

Comparison of two software systems for quantification of myocardial blood flow in patients with hypertrophic cardiomyopathy

Backgorund: Quantification of myocardial blood flow (MBF) by positron emission tomography (PET) is important for investigation of angina in hypertrophic cardiomyopathy (HCM). Several software programs exist for MBF quantification, but they have been mostly evaluated in patients (with normal cardiac geometry), referred for evaluation of coronary artery disease (CAD). Software performance has not been evaluated in HCM patients who frequently have hyperdynamic LV function, LV outflow tract (LVOT) obstruction, small LV cavity size, and variation in the degree/location of LV hypertrophy. Aim: We compared results of MBF obtained using PMod, which permits manual segmentation, to those obtained by FDA-approved QPET software which has an automated segmentation algorithm. Methods: {sup 13}N-ammonia PET perfusion data were acquired in list mode at rest and during pharmacologic vasodilation, in 76 HCM patients and 10 non-HCM patients referred for evaluation of CAD (CAD group.) Data were resampled to create static, ECG-gated and 36-frame-dynamic images. Myocardial flow reserve (MFR) and MBF (in ml/min/g) were calculated using QPET and PMod softwares. Results: All HCM patients had asymmetric septal hypertrophy, and 50% had evidence of LVOT obstruction, whereas non-HCM patients (CAD group) had normal wall thickness and ejection fraction. PMod yielded significantly higher values for global and regional stress-MBF and MFR than for QPET in HCM. Reasonably fair correlation was observed for global rest-MBF, stress-MBF, and MFR using these two softwares (rest-MBF: r = 0.78; stress-MBF: r = 0.66.; MFR: r = 0.7) in HCM patients. Agreement between global MBF and MFR values improved when HCM patients with high spillover fractions (> 0.65) were excluded from the analysis (rest-MBF: r = 0.84; stress-MBF: r = 0.72; MFR: r = 0.8.) Regionally, the highest agreement between PMod and QPET was observed in the LAD territory (rest-MBF: r = 0.82, Stress-MBF: r = 0.68) where spillover fraction was the lowest. Unlike HCM patients, the non-HCM patients (CAD group) demonstrated excellent agreement in MBF/MFR values, obtained by the two softwares, when patients with high spillover fractions were excluded (rest-MBF: r = 0.95; stress-MBF: r = 0.92; MFR: r = 0.95). Conclusions: Anatomic characteristics specific to HCM hearts contribute to lower correlations between MBF/MFR values obtained by PMod and QPET, compared with non-HCM patients. These differences indicate that PMod and QPET cannot be used interchangeably for MBF/MFR analyses in HCM patients.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Implementation and Verification of Cable Bending Stiffness in MoorDyn

The relatively large motions experienced by floating wind turbines and wave energy converters pose a challenge for power cables, whose internal components provide significant bending resistance and are sensitive to deformation. The behavior and associated design considerations of power cables in these highly dynamic applications make coupled analysis relevant for design. Bending stiffness capabilities have recently been added to the lumped-mass mooring dynamics model MoorDyn to enable simulation of dynamic power cables. MoorDyn is a common modeling choice for floating wind energy simulation (often coupled with OpenFAST) and floating wave energy converter simulation (often coupled with WEC-Sim) but the model’s previous line elasticity formulation only considered axial stiffness. To properly capture the dynamics of power cables, a bending stiffness model has been added that approximates cable curvature based on the difference in tangent vectors of adjacent elements. The resulting bending moment is realized by applying forces on adjacent nodes, enabling cable modeling while leaving the underlying lumped-mass formulation unchanged. In this paper, the new bending stiffness implementation is verified in static conditions against analytical solutions and then in a dynamic power cable scenario in comparison with the commercial simulator OrcaFlex. The dynamic scenario uses prescribed motions and includes wave loadings on the cable. Results indicate correct implementation of bending stiffness and show close agreement with OrcaFlex.

50 EE - Wind and Water Power Program - Water (EE-4↗