Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel programming patterns”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Solid-propellant rocket motor internal ballistics performance variation analysis, phase 5

The results of research aimed at improving the predictability of internal ballistics performance of solid-propellant rocket motors (SRM's) including thrust imbalance between two SRM's firing in parallel are presented. Static test data from the first six Space Shuttle SRM's is analyzed using a computer program previously developed for this purpose. The program permits intentional minor design biases affecting the imbalance between any two SMR's to be removed. Results for the last four of the six SRM's, with only the propellant bulk temperature as a non-random variable, are generally within limits predicted by theory. Extended studies of internal ballistic performance of single SRM's are presented based on an earlier developed mathematical model which includes an assessment of grain deformation. The erosive burning rate law used in the model is upgraded and made more general. Excellent results are obtained in predictions of the performances of five different SRM's of quite different sizes and configurations. These SRM's all employ PBAN type propellants with ammonium perchlorate oxidizer and 16 to 20% aluminum except one which uses carboxyl terminated butadiene binder. The only non-calculated parameters in the burning rate equations that are changed for the different SRM's are the zero crossflow velocity burning rate coefficients and exponents. The results, in general, confirm the importance of grain deformation. The improved internal ballistic model makes practical development of an effective computer program for application of an optimization technique to SRM design which is also demonstrated. The program uses a pattern search technique to minimize the difference between a desired thrust-time trace and one calculated based on the internal ballistic model.

Sforzini, R. H.↗
NASA NTRS2003

Afferent innervation patterns of the saccule in pigeons

The innervation patterns of vestibular saccular afferents were quantitatively investigated in pigeons using biotinylated dextran amine as a neural tracer and three-dimensional computer reconstruction. Type I hair cells were found throughout a large portion of the macula, with the highest density observed in the striola. Type II hair cells were located throughout the macula, with the highest density in the extrastriola. Three classes of afferent innervation patterns were observed, including calyx, dimorph, and bouton units, with 137 afferents being anatomically reconstructed and used for quantitative comparisons. Calyx afferents were located primarily in the striola, innervated a number of type I hair cells, and had small innervation areas. Most calyx afferent terminal fields were oriented parallel to the anterior-posterior axis and the morphological polarization reversal line. Dimorph afferents were located throughout the macula, contained fewer type I hair cells in a calyceal terminal than calyx afferents and had medium sized innervation areas. Bouton afferents were restricted to the extrastriola, with multi-branching fibers and large innervation areas. Most of the dimorph and bouton afferents had innervation fields that were oriented dorso-ventrally but were parallel to the neighboring reversal line. The organizational morphology of the saccule was found to be distinctly different from that of the avian utricle or lagena otolith organs and appears to represent a receptor organ undergoing evolutionary adaptation toward sensing linear motion in terrestrial and aerial species.

NASA Discipline Neuroscience↗
NASA NTRS1994

Effect of hexylene glycol-altered microtubule distributions on cytokinesis and polar lobe formation in fertilized eggs of Ilyanassa obsoleta

Some effects of gravity on early morphogenesis are correlated with microtubule locations within cells. During first cleavage in Ilyanassa obsoleta embryos, a transitory polar lobe constriction forms and then relaxes, allowing the polar lobe to merge with one daughter cell. If the polar lobe is equally divided or removed, morphogenesis is severely disrupted. To examine microtuble locations during early Ilyanassa development, eggs were fixed and stained for polymerized alpha-tubulin during first cleavage. The mitotic apparatus assembles at the animal pole. The cleavage furrow forms between the asters, constricting to a stabilized intercellular bridge encircling midbody-bound microtubules, whereas the polar lobe constriction forms below and parallel to the spindle, constricting to a transitory intercellular bridge encircling no detectable microtubules. At metaphase an alpha-tubulin epitope is distributed throughout the spindle, whereas a beta-tubulin epitope is present predominantly in the asters. Incubation in hexylene glycol, a drug that increases microtubule polymerization, during mitosis causes the polar lobe constriction to tighten around polymerized alpha-tubulin and remain stably constricted. If hexylene glycol is removed, alpha-tubulin staining disappears from the polar lobe constriction, which relaxes, whereas microtubules remain in the cleavage furrow, which remains constricted. These observations suggest that asymmetric distribution of microtubules affects early Ilyanassa cleavage patterns, and that continued presence of microtubules extending through an intercellular bridge is important for stabilization of the bridge constriction prior to completion of cytokinesis. These data provide the basis for further analysis of the role of microtubules in possible microgravity disruptions of Ilyanassa development.

NASA Discipline Cell Biology↗
DOE OSTI2021

Restoration Strategy for Active Distribution Systems Considering Endogenous Uncertainty in Cold Load Pickup

Cold load pickup (CLPU) phenomenon is identified as the persistent power inrush upon a sudden load pickup after an outage. Under the active distribution system (ADS) paradigm, where distributed energy resources (DERs) are extensively installed, the decreased outage duration can induce a strong interdependence between CLPU pattern and load pickup decisions. In this paper, we propose a novel modelling technique to tractably capture the decision-dependent uncertainty (DDU) inherent in the CLPU process. Subsequently, a two-stage stochastic decision-dependent service restoration (SDDSR) model is constructed, where first stage searches for the optimal switching sequences to decide step-wise network topology, and the second stage optimizes the detailed generation schedule of DERs as well as the energization of switchable loads. Further, to tackle the computational burdens introduced by mixed-integer recourse, the progressive hedging algorithm (PHA) is utilized to decompose the original model into scenario-wise subproblems that can be solved in parallel. The numerical test on modified IEEE 123-node test feeders has verified the efficiency of our proposed SDDSR model and provided fresh insights into the monetary and secure values of DDU quantification.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Production Application Performance Data Streaming for System Monitoring

In this article, we present an approach to streaming collection of application performance data. Practical application performance tuning and troubleshooting in production high-performance computing (HPC) environments requires an understanding of how applications interact with the platform, including (but not limited to) parallel programming libraries such as Message Passing Interface (MPI). Several profiling and tracing tools exist that collect heavy runtime data traces either in memory (released only at application exit) or on a file system (imposing an I/O load that may interfere with the performance being measured). Although these approaches are beneficial in development stages and post-run analysis, a systemwide and low-overhead method is required to monitor deployed applications continuously. This method must be able to collect information at both the application and system levels to yield a complete performance picture. In our approach, an application profiler collects application event counters. A sampler uses an efficient inter-process communication method to periodically extract the application counters and stream them into an infrastructure for performance data collection. We implement a tool-set based on our approach and integrate it with the Lightweight Distributed Metric Service (LDMS) system, a monitoring system used on large-scale computational platforms. LDMS provides the infrastructure to create and gather streams of performance data in a low overhead manner. We demonstrate our approach using applications implemented with MPI, as it is one of the most common standards for the development of large-scale scientific applications. We utilize our tool-set to study the impact of our approach on an open source HPC application, Nalu. Our tool-set enables us to efficiently identify patterns in the behavior of the application without source-level knowledge. We leverage LDMS to collect system-level performance data and explore the correlation between the system and application events. Finally, we demonstrate how our tool-set can help detect anomalies with a low latency. We run tests on two different architectures: a system enabled with Intel Xeon Phi and another system equipped with Intel Xeon processor. Our overhead study shows our method imposes at most 0.5% CPU usage overhead on the application in realistic deployment scenarios.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

An Integrated View of Hydraulic Induced Fracture Geometry in Hydraulic Fracture Test Site 2

Until relatively recently, most of our understanding about the geometry of hydraulically induced fractures in unconventional reservoirs has been inferred by patterns observed via microseismic monitoring. Important information, such as fracture-azimuth, fracture-length, height-growth and other geometric data has been determined from microseismic “event clouds”, which in turn has been used to interpret a stimulated reservoir volume (SRV). It is important to note that relatively limited microseismic results have been acquired in conjunction with other diagnostics and validated in a consistent manner. The Hydraulic Fracture Test Site 2 (HFTS2) in the Permian Delaware Basin provides a unique opportunity to compare the frac geometry interpretations derived from multiple frac diagnostic tools and at multiple scales. The integrated view that emerges from all the frac diagnostic tools is that in HFTS2 the geometries of hydraulically induced fractures are not random and with a highly complex branching architecture. In fact, the fracs are mostly vertical, parallel plana domains, and both vertically and horizontally asymmetrical. The hydraulically induced fractures have a very consistent azimuth with a strike that matches SHmax. At the distance corresponding to the spacing between nearby wells in HFTS2, ~ 660 ft, the pattern of strain interceptions observed from Low Frequency Distributed Acoustic Sensing (LF-DAS) shows far-field dimensions that are consistent with the stage length dimensions at the stimulated well. In general, larger stages with higher number of clusters create larger number of interceptions and wider stimulated intervals than that of shorter stages with the same cluster to cluster spacing and fewer clusters. The hydraulically induced fractures also show a clear tendency to preferentially grow upwards as shown by the microseismic data, by the relative lower number of LF-DAS interceptions in nearby deeper wells and finally by the observed LF-DAS strain patterns in the vertical observation well. In the near field, the geometries of the fracture domains are also consistent. As shown by the distributed strain measurements obtained during production via the monitoring of Rayleigh Frequency Shift (RFS), each cluster has a separate, non-overlapping, frac-zone-domain that are generally centered around the locations of each perforation cluster. Overall, the hydraulically induced fractures in HFTS2 are interpreted as occurring in swarms associated with individual clusters, each with its own unique geometry but with similar predictable geometric and propagation tendencies. The comprehensive frac diagnostic program in HFTS2, the well-pad layout that permits the investigation of the Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/1-21URTC/D011S018R003/2477511/urtec-2021-5396-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5396 2 impact that nearby producing wells has on frac geometry, together with an innovative Design of Experiment (DofE) that includes sequential-fracing, consistent cluster-to-cluster distances and few single-perforation cluster stages, has allowed us to gain a unique perspective about the geometry of the induced hydraulic fractures on this pad. Finally, the HFTS2 dataset has also provided us with a better understanding of the applicability and limitations of the different diagnostic tools.

58 GEOSCIENCES↗
NASA NTRS1989

Evolution and extinction in the marine realm: some constraints imposed by phytoplankton

The organic and mineralized remains of planktonic algae provide a rich record of microplankton evolution extending over nearly half of the preserved geological record. In general, Phanerozoic patterns of phytoplankton radiation and extinction parallel those documented for skeletonized marine invertebrates, both augmenting and constraining thought about evolution in the oceans. Rapidly increasing knowledge of Proterozoic plankton is making possible the recognition of additional episodes of diversification and extinction that antedate the Ediacaran radiation of macroscopic animals. In contrast to earlier phytoplankton history, the late Mesozoic and Cainozoic record is documented in sufficient detail to constrain theories of mass extinction in more than a general way. Broad patterns of diversity change in planktonic algae show similarities across the Cretaceous-Tertiary and Eocene-Oligocene boundaries, but detailed comparisons of origination and extinction rates in calcareous nannoplankton, as well as other algae and skeletonized protozoans, suggest that the two episodes were quite distinct. Common causation appears unlikely, casting doubt on monolithic theories of mass extinction, whether periodic or not. Studies of mass extinction highlight a broader class of insights that paleontologists can contribute to evolutionary biology: the evaluation of evolutionary change in the context of evolving Earth-surface environments.

NASA Discipline Number 52-30↗
NASA NTRS1992

The back photoreaction of the M intermediate in the photocycle of bacteriorhodopsin: mechanism and evidence for two M species

The back photoreaction of the M intermediate in the photocycle of bacteriorhodopsin is investigated both for the native pigment and its D96N mutant. The experimental setup is based on creating the M intermediate by a first pulse, followed by a (blue) laser pulse which drives the back photoreaction of M. Experiments are carried out varying the delay between the two pulses, as well as the temperature over the -25 degrees C-20 degrees C range. It is found that the kinetic patterns of the M back photoreaction change with time after the generation of this intermediate. The data provide independent evidence for the suggestion of a photocycle mechanism based on two distinct M intermediates. They are thus in keeping with the consecutive model of Varo and Lanyi (Biochemistry 30, 5016-5022; 1991), although they cannot exclude other models such as those based on branched or parallel cycles. More generally, we offer a "photochemical" approach to discriminating between intermediate stages in the photocycle which does not depend on spectroscopic and/or kinetic data. While markedly affecting the rate of the M --> N transition in the photocycle, the rate of the thermal step in back photoreaction of M, at both room and low temperatures, is not significantly affected by the D96N mutation. It is proposed that while Asp 96 is the Schiff-base protonating moiety in the M --> N transition, another residue (most probably Asp 85) reprotonates the Schiff base following light absorption by M.

NASA Discipline Exobiology↗
DOE OSTI2014

Damsel: A Data Model Storage Library for Exascale Science

Computational science applications have been described as having one of seven motifs (the “seven dwarfs”), each having a particular pattern of computation and communication. From a storage and I/O perspective, these applications can also be grouped into a number of data model motifs describing the way data is organized and accessed during simulation, analysis, and visualization. Major storage data models developed in the 1990s, such as Network Common Data Format (netCDF) and Hierarchical Data Format (HDF) projects, created support for more complex data models. Development of both netCDF and HDF5 was influenced by multi-dimensional dataset storage requirements, but their access models and formats were designed with sequential storage in mind (e.g., a POSIX I/O model). Although these and other high-level I/O libraries have had a beneficial impact on large parallel applications, they do not always attain a high percentage of peak I/O performance due to fundamental design limitations, and they do not address the full range of current and future computational science data models. The goal of this project is to enable exascale computational science applications to interact conveniently and efficiently with storage through abstractions that match their data models. The project consists of three major activities: (1) identifying major data model motifs in computational science applications and developing representative benchmarks; (2) developing a data model storage library, called Damsel, that supports these motifs, provides efficient storage data layouts, incorporates optimizations to enable exascale operation, and is tolerant to failures; and (3) productizing Damsel and working with computational scientists to encourage adoption of this library by the scientific community. The product of this project, Damsel library, is openly available for download from http://cucis.ece.northwestern.edu/projects/DAMSEL. Several case studies and application programming interface reference are also available to assist new users to learn to use the library.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Association Mapping of Cell Wall Synthesis Regulatory Genes and Cell Wall Quality in Switchgrass

Inefficient conversion of biomass to biofuels is one of the main barriers for biofuel production from such materials. Approximately half of polysaccharides in biomass remain unused by typical biochemical conversion methods. Conversion efficiency is influenced by the composition and structure of cell walls of biomass. Grasses such as wheat, maize, and rice, as well as dedicated perennial bioenergy crops, like switchgrass, make up ~55% of biomass that can be produced in the United States. Grass cell walls have a different composition and patterning compared with dicotyledonous plants, including the well-studied model plant, Arabidopsis. This project identified genetic determinants of cell wall composition in grasses using both naturally occurring genetic variation of switchgrass and gene network reconstruction and functional assays in rice. In addition, the project linked functional data in rice and other species to switchgrass improvement efforts through curation of the most abundant class of regulators in the switchgrass genome. Characterizing natural diversity of switchgrass for variation in cell wall composition and properties, also known as quality, provides an unbiased avenue for identifying biologically viable diversity in switchgrass cell walls. To characterizing natural diversity, this project generated cell wall composition and enzymatic deconstruction data for ~450 genotypes of the Switchgrass Southern Association Collection (SSAC), a diverse collection composed of 36 switchgrass accessions from the southern U.S. distribution of switchgrass. Comparing these data with other measures of cell wall quality for the same samples demonstrated the complementary nature of the diverse characterization platforms now being used for biomass characterization. Association of the composition data with ~3.2K single nucleotide variant markers identified six significant single nucleotide variant markers co-associated with digestibility and another compositional trait. These markers might be used to select switchgrass genotypes with improved composition in breeding programs for biofuel and forage production. Because the SSAC continues to be characterized by collaborators in the bioenergy community, the data generated will be used to identify additional markers in higher resolution genotyping data to approach identifying the genes and alleles that cause natural variation in switchgrass cell wall quality. For example, these markers can be surveyed in the 2100-member Oklahoma Southern and Northern Lowland switchgrass collections that this project also characterized. An orthogonal approach to biodiversity studies, using comparative functional genomics permits systematic querying of how much regulatory information is likely to be transferable from dicots to grasses and use of accumulated functional genomics resources for better-characterized grass species, such as rice, itself a biomass source in global agriculture and in certain regions. The project generated and tested a number of specific hypotheses regarding cell wall transcription factors and enzymes of grasses. To aid identification of cell wall regulators, the project assembled a novel, highdepth and -quality gene association network using a general linearized model scoring system to combine rice gene network data. Using known or putative orthologs of Arabidopsis cell wall biosynthesis genes and regulators, the project pulled from this network a cell wall sub-network that includes 96 transcription factors. Reverse genetics of a co-ortholog of the Arabidopsis MYB61 transcription factor in rice revealed that this regulatory node has evolved the ability to regulate grass-specific cell wall synthesis enzymes. A transcription factor with such activity has not been previously characterized to our knowledge, representing a major conclusion of this work. Changes in gene expression in a protoplast-based assay demonstrated positive or negative roles in cell wall regulation for eleven other transcription factors from the rice gene network. Eight of fifteen (53%) of these have not previously been examined for this function. Some of these may represent novel grass-diverged cell wall regulators, while others are likely to have this function across angiosperms. A parallel effort of this project to expand knowledge of enzymes that have evolved to function in grass cell wall synthesis, revealed that a grass-diverged enzyme in rice, OsAT 5 , ferulates monolignols that are naturally incorporated into grass cell walls. This finding opens potential natural selection avenues for improving biomass composition for downstream processing by weak base pretreatment. Thus, this project has significantly expanded knowledge of cell wall synthesis and regulation in rice, information that can be used in reverse genetics and synthetic biology approaches to re-engineer cell walls for improved production of biofuel and high-value products. To lay the foundation for translating these results directly for switchgrass improvement, the project employed a comparative phylogenetic analysis of the major group of cell wall transcription factors that have been found to function in cell wall regulation, the R 2 R 3 MYBs. This analysis concluded that known cell wall regulators are largely conserved across switchgrass, rice, maize, poplar, and Arabidopsis. This interpretation is also largely consistent with the gene network analysis described above, though both approaches provide evidence that some co-orthologs of Arabidopsis regulators have diminished or increased in importance based on gene expression patterns. Also, several clades containing dicot cell wall regulators have expanded, consistent with the evolution of new cell wall regulators. This latter result is supported by functional analysis of the R 2 R 3 MYB protein SWAM 1 in a collaboration between this project and the DOE-funded group of Dr. S. Hazen at the University of Massachusettes. The curation of the switchgrass genome through this project provides specific targets for future engineering of switchgrass cell wall regulation and may also facilitate identification of regulators that underlie the molecular markers that are genetically linked to differences in cell wall quality. With the goal of spurring further research and technological developments in lignocellulosic biofuel production, this work has been communicated to the bioenergy and cell wall communities though various presentations and publications. To date, three manuscripts have been published, two others are near to publication, three others are in an advanced state, and two to four more are likely to be written based on analyses still in progress. In addition, project participants have presented thirteen posters and talks at regional, national, and international meetings about aspects of this project. In sum, the work supported by this funding has made and communicated significant progress in identifying the genes that grasses use for cell wall synthesis and regulation, information that will be used by project participants and others to improve the efficiency of conversion of lignocellulosic biomass to biofuels.

09 BIOMASS FUELS↗
DOE OSTI2014

Parallel-scanning tomosynthesis using a slot scanning technique: Fixed-focus reconstruction and the resulting image quality

Purpose: Parallel-scanning tomosynthesis (PS-TS) is a novel technique that fuses the slot scanning technique and the conventional tomosynthesis (TS) technique. This approach allows one to obtain long-view tomosynthesis images in addition to normally sized tomosynthesis images, even when using a system that has no linear tomographic scanning function. The reconstruction technique and an evaluation of the resulting image quality for PS-TS are described in this paper. Methods: The PS-TS image-reconstruction technique consists of several steps (1) the projection images are divided into strips, (2) the strips are stitched together to construct images corresponding to the reconstruction plane, (3) the stitched images are filtered, and (4) the filtered stitched images are back-projected. In the case of PS-TS using the fixed-focus reconstruction method (PS-TS-F), one set of stitched images is used for the reconstruction planes at all heights, thus avoiding the necessity of repeating steps (1)–(3). A physical evaluation of the image quality of PS-TS-F compared with that of the conventional linear TS was performed using a R/F table (Sonialvision safire, Shimadzu Corp., Kyoto, Japan). The tomographic plane with the best theoretical spatial resolution (the in-focus plane, IFP) was set at a height of 100 mm from the table top by adjusting the reconstruction program. First, the spatial frequency response was evaluated at heights of −100, −50, 0, 50, 100, and 150 mm from the IFP using the edge of a 0.3-mm-thick copper plate. Second, the spatial resolution at each height was visually evaluated using an x-ray test pattern (Model No. 38, PTW Freiburg, Germany). Third, the slice sensitivity at each height was evaluated via the wire method using a 0.1-mm-diameter tungsten wire. Phantom studies using a knee phantom and a whole-body phantom were also performed. Results: The spatial frequency response of PS-TS-F yielded the best results at the IFP and degraded slightly as the distance from the IFP increased. A visual evaluation of the spatial resolution using the x-ray test pattern indicated that the resolution was 1.8 lp/mm at the IFP and 1.2 lp/mm at heights of −100 and 100 mm from the IFP. The authors demonstrated that a spatial resolution of 1.2–1.8 lp/mm could be obtained within heights of 200 mm of the IFP. The slice sensitivity varied between 11.1 and 13.8 mm for heights between −50 and 100 mm, and there was no critical change in the slice sensitivity within a height range of 150 mm around the IFP. The phantom results demonstrated that tomosynthesis and long-view images could be reconstructed. Conclusions: PS-TS-F provides tomosynthesis images while using low-cost systems that have no tomographic scanning function, such as tableside-controlled universal R/F systems or universal radiographic systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2018

The OpenSE Cookbook: A Practical, Recipe Based Collection of Patterns, Procedures, and Best Practices for Executable Systems Engineering for the Thirty Meter Telescope

The OpenSE Cookbook is an open-sourced collection of patterns, procedures, and best practices targeted for systems engineers who seek guidance on applying model-based and executable systems engineering (MBSE) using SysML. Its content has emerged from the system level modeling effort on the European Framework Program 6 (FP6) and the Thirty Meter Telescope (TMT). The TMT MBSE approach applied the Executable Systems Engineering Method (ESEM) and the open-source Engineering Environment (OpenMBEE) to specify, analyze, and verify requirements of TMT’s Alignment and Phasing System (APS) and the Narrow Field Infrared Adaptive Optics System (NFIRAOS). In these applications, implicit dependencies are made explicit in a formal model through the use of ESEM, OpenMBEE, and SysML modeling constructs. The value proposition for applying this MBSE approach was to establish precise requirements and fine-grained traceability to system designs, and to verify key requirements beginning early in development. The integration of ESEM and the OpenMBEE tooling infrastructure (providing linked-data and web-operability) is a significant added value for the MBSE approach. The APS is responsible for the overall pre-adaptive optics wavefront quality, using starlight to measure wavefront errors and align the TMT optics. In the formally integrated and executable SysML model, simulations are performed to analyze the impact of changed requirements and verify specified constraints for various operational scenarios. The APS team used several modeling patterns to capture information such as the requirements, the operational scenarios, involved subsystems and their interaction points, the estimated or required time durations, and the mass and power consumption. Adaptive optics systems are designed to sense real-time atmospheric turbulence and correct the telescope’s optical beam to remove its effect. The system model for the adaptive optics operational modes was developed to capture sequence behaviors and operational scenarios to run Monte-Carlo simulations for verifying acquisition time, observing efficiency, and operational behavior requirements. The model is particularly useful for investigating the effect of parallelization, identifying interface issues, and re-ordering sequence acquisition tasks. A former version of the Cookbook (which is now updated to MBSE challenges, goals, and lessons learned) included modeling guidelines and conventions for all system aspects, hierarchy levels, and views, which were developed during for the Active Phasing Experiment (APE), an opto-mechatronical system technology demonstrator for the Extremely Large Telescope (ELT). The Cookbook utilizes the above mentioned system models as real-world case-studies to demonstrate and document the applications of the recipes, providing also instructional examples and addressing the available tooling support. The Cookbook is accompanied by a number of SysML models and aodel libraries which facilitate model authoring and maintenance. The Cookbook covers the different aspects of Systems Engineering such as management of Requirements, Design (behavior and structure), Interfaces, Interdisciplinary Integration, Analysis, Trade Studies, and Technical Resources. This paper presents the background, motivation, architecture, and highlights some key content of the Cookbook. For example, interface management, error budget management, requirements verification, Monte Carlo driven analysis, and timing analysis of operational scenarios. The paper discusses how the capabilities of OpenMBEE contributed significantly to the adoption of executable systems engineering.

Brower, Eric↗
DOE OSTI2021

Phasor-Measurement-Unit-Based Data Analytics Using Digital Twin and PhasorAnalytics Software

A major objective of this project was to apply GE’s commercial machine learning and data analytics toolsets to large-scale, real-world, anonymized Phasor Measurement Unit (PMU) datasets in order to extract signatures, correlated and/or causal factors, and precursor patterns associated with significant power system phenomena. The project had a particular emphasis on extraction of insights relevant to asset health monitoring, real-time load modeling and cybersecurity monitoring. Additionally, the team was directed to undertake a comprehensive data quality analysis for the provided datasets and encouraged to estimate the ‘machine-learning readiness’ of the datasets by documenting any major obstacles to the application of commercial machine learning algorithms. To accomplish the aforementioned objectives, the project team’s work centered around the identification of key event signatures and application of the identified event signatures for event detection and event classification. The industry-validated, semi-supervised machine learning strategy employed for event signature identification involved several major tasks, including data-preprocessing, generation of an overabundance of features, normal data identification, normality modeling, and event signature identification through a methodical, quantitative ranking of features in order of relevance to each studied event type. Throughout the project, data quality issues and mitigation techniques were investigated. In this report, insights are provided regarding the readiness of the provided synchrophasor datasets for application of machine learning and data analytics. The methodologies employed for this technical strategy are summarized in this report. With regards to data preprocessing and feature generation, the provided Training and Test Datasets were ingested into GE’s big data environment. Subsequently, the team applied bad data cleansing and data imputation scripts, event detection scripts, and application programming interfaces (APIs) to the datasets for convenient data access. The project team completed development and validation of dozens of physics-based, statistics-based and transformation-based feature functions used for the extraction of over 60 synchrophasor features. Using a new parallel feature generation technology developed on this project, over 60 features have been rapidly generated for the full two years’ worth of Training and Test Dataset data associated with both the Eastern and Western interconnects. Even accommodating for temporal down-sampling inherent to the feature extraction procedure, this parallel feature generation activity resulted in a massive feature set with a storage requirement approximately equal to that of the raw training dataset itself. With regards to normal data identification and normality modeling, a normality model was built using the feature data extracted from the Training Dataset and iteratively refined subsequent to incremental adjustments and expansions of the Training Dataset feature data. With respect to event characterization and signature identification, an event signature identification pipeline was developed and used in conjunction with the normality model to identify over 15 event signatures for key event categories within the Training Dataset. The identified event signatures were used to characterize hundreds of key events in terms of relative severity, duration, and location of the event. An investigation was undertaken to identify correlated and causal factors involved in transformer events. A separate investigation into temporal trends in ring-down analysis results was undertaken to determine possible associations between system dynamics and various other factors such as loading, season or year. To validate the identified event signatures, additional work was undertaken to develop signature-based anomaly detection and classification tools suitable for convenient application to the synchrophasor datasets. The anomaly detection and classification tools, suitable for online application, were then applied to the entirety of the Eastern Interconnect Training and Test Datasets. Performance of the event detection and classification tools was evaluated upon receipt of the Test Dataset event logs (i.e., the labels for events contained in the Test Dataset), and promising results were obtained despite several challenges (documented herein) associated with application of supervised or semi-supervised machine learning methods to large-scale, anonymized datasets. Finally, the detection and classification tools were used to detect, classify, and characterize thousands of new events not included in the original event logs provided by the DOE within both the Training and Test Datasets.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

IAEA Coordinated Research Project on HTGR Physics, Thermal Hydraulics, and Depletion Uncertainty Analysis: PHISICS/RELAP5-3D Results for the Phase III Coupled Core Exercises

This report details the Parallel and Highly Innovative Simulation for INL Code System (PHISICS)/Reactor Excursions and Leak Analysis Program (RELAP5)-3D results obtained for the neutronics stand-alone core exercises defined for Phase II of the International Atomic Energy Agency (IAEA) Coordinated Research Project (CRP) on high-temperature gas cooled reactor (HTGR) uncertainty analysis in modeling (UAM). The Phase II models and results are linked to the earlier Standardized Computer Analyses for Licensing Evaluation (SCALE)/Sampler/New ESC-based Weighting Transport (NEWT) data generated for the lattice physics (lattice) stage Phase I of the CRP. The focus of this report is the propagation of cross-section uncertainties from the lattice to the core calculation (core) phase of the prismatic 350 megawatt (MW) General Atomics (GA) Modular HTGR (MHTGR). The Phase II neutronics exercises consist of two sub-cases (II-2a and II-2b) specified for the fresh and mixed core loading of the MHTGR-350 prismatic high-temperature reactor design. Exercises II-2a and II-2b are performed without temperature feedback using isothermal temperature profiles. For the nominal lattice physics (i.e., best-estimate) model, a series of fresh and depleted single block and super cell SCALE/KENO and SCALE/NEWT models were developed at Idaho National Laboratory (INL) to assess the impact of spectral effects on the few-group cross sections. A combination of these lattice cells were used in a super cell calculation to construct 8-group cross-section libraries for use in the INL-developed PHISICS/RELAP5-3D code suite. The block-level homogenized libraries were assigned to the 220 fuel zones utilized in the MHTGR-350 PHISICS core model, according to the fresh and mixed fresh/depleted core loading pattern. In order to evaluate the effect of uncertainty propagation, several sets consisting of 1,000 perturbed NEWT lattice and PHISICS/RELAP5-3D core calculations were performed using the SCALE 6.2 ENDF/B-VII.1 cross-section co-variance data and the stochastic sampling module Sampler in SCALE 6.2.0. The Figures of Merit (FOM) investigated include the mean and standard deviation uncertainty values obtained for the eigenvalues, control rod worth values, block-averaged power profiles, power axial offsets, and isothermal temperature feedback coefficients. A sensitivity study on the main nuclide reaction contributors to the FOM uncertainties has not yet been performed since the SCALE/PHISICS/RAVEN coupling with ENDF/B-VII co-variance data is still being developed. It was found that uncertainties in the cross sections lead to uncertainties in the local power values of up to 1.5%. The variances between the fresh and mixed/depleted core models are not as large as the axial variations between the top and bottom within each core model. The core eigenvalue uncertainty is approximately 0.5%, which is significant in comparison to typical eigenvalue difference targets (less than ~0.1% for LWRs). The power peaking and axial offset uncertainties are smaller (0.10% – 0.15%) and probably below the level of significance in typical core design metrics. The use of an 8-group energy structure slightly under-estimates some of the uncertainties, but in absolute terms, the uncertainty results are reasonably close to the 26-group results. The 26- and 8-group mean and nominal results differ by a larger degree, and if it can be afforded in terms of calculation resources, the 26-group structure should be used for HTGR calculations. The 8-group structure can, however, still be useful for reasonably accurate uncertainty assessments of transient cases where computational time is important. It was furthermore found that the use of a super cell during the cross-section generation process for better representation of the spectral environment at the core periphery did not lead to significant changes in the core eigenvalues, axial offset, and power peaking uncertainties. The impact of the super cell cross sections on the nominal and mean eigenvalues is, however, important enough (0.6% – 1.1% lower) to justify the use of these super cells for HTGR core simulation to represent the softer spectral environment in the core peripheral regions for nominal best-estimate calculations. The use of super cells for HTGR uncertainty studies is, therefore, not strictly necessary, but the improvements in nominal results could justify their use.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Improving ICARUS Track Reconstruction Algorithms

The ICARUS experiment is part of the Short-Baseline Neutrino (SBN) program at Fermilab. The main goal of the experiment is to investigate the possibility of sterile neutrinos in the O(1 eV) mass region and provide clarification of the anomaly detected from the Liquid Scintillator Neutrino Detector (LSND) and MiniBooNE experiments. The ICARUS-T600 detector is a Liquid Argon Time Projection Chamber (LAr-TPC), that can provide excellent 3D imaging and calorimetric reconstruction of any ionizing particles. This detection technique allows a detailed study of neutrino interactions, spanning a wide energy spectrum (from a few keV to several hundreds of GeV). The detector consists of two identical adjacent modules, filled with a total of 760 tons of ultra-pure liquid argon. Each module houses two LAr-TPCs separated by a common cathode with a maximum drift distance of 1.5 m, equivalent to about 1 ms drift time for the nominal $500$ V/m electric drift field. The anode is made of three parallel wire planes positioned 3 mm apart, where the stainless-steel wires are oriented on each plane at a different angle with respect to the horizontal direction ($+60^\degree$,$-60^\degree$,$0^\degree$). The first two planes (Induction 1 and Induction 2) provide a non-destructive charge measurement, whereas the ionization charge is fully collected by the last collection plane. In total, 53248 wires with a 3 mm pitch and length up to 9 m are installed in the detector. In the first stage of the reconstruction, segments of waveforms corresponding to physical signals (hits) are searched for in the deconvolved wire waveform with a threshold-based hit-finding algorithm. Each hit is then fitted with a Gaussian, whose area is proportional to the number of drift electrons generating the signal. In the second stage of the reconstruction, hits are passed as input to Pandora, a framework software composed of different pattern recognition algorithms, that performs a 3D reconstruction of the full image recorded in the collected event, including the identification of interaction vertices and tracks and showers inside the TPC. These are organized into a hierarchical structure (called slice) of particles generated starting from a primary interaction vertex. In some cases, related to the inefficiencies in the hit detection or excessive deflection of the particle trajectory, Pandora breaks the particle's track into two or more smaller pieces and considers each piece as an independent track. We studied this phenomenon focusing on primary muons from ν_μ CC interactions contained in a single module with a track at least 20 cm long, to exclude delta rays. The study determined that about $7-8\%$ of the muon tracks are broken. Approximately $80\%$ of the times, Pandora assigns all segments of the track to the same slice (intra-slice track split), while in the remaining $20\%$ of the cases, one of the segments is associated with another slice (extra-slice track split). To mitigate this phenomenon, we designed an algorithm that detects and stitches the tracks broken by Pandora for the intra-slice split. In Monte Carlo simulations, the algorithm showed an efficiency exceeding $80\%$ and a purity exceeding $93\%$.

Ricci, Alessandro Maria [Pisa U.; INFN, Pisa] (ORC↗
DOE OSTI2017

CO 2 Capture by Cold Membrane Operation with Actual Power Plant Flue Gas

The main objective of the project was to develop a post-combustion CO 2 capture process based on the hybrid cold temperature membrane operation. The CO 2 in the flue gas from coal fired power plant is pre-concentrated to >60% CO 2 in the first stage membrane operation followed by further liquefaction of permeate stream to achieve >99% CO 2 purity. The aim of the project was based on DOE program goal of 90% CO 2 capture with >95% CO 2 purity from Pulverized Coal (PC) fired power plants with $40/tonne of carbon capture cost by 2025. The project moves the technology from TRL 4 to TRL 5. The project involved optimization of Air Liquide commercial 12” PI-1 bundle to improve the bundle productivity by >30% compared to the previous baseline (DE-FE0004278) using computational fluid dynamics (CFD) modeling and bundle testing with synthetic flue gas at 0.1 MWe bench scale skid located at Delaware Research and Technology Center (DRTC). In parallel, the next generation polyimide based novel PI-2 membrane was developed with 10 times CO 2 permeance compared to the commercial PI-1 membrane. The novel PI-2 membrane was scaled from mini-permeator to 1” permeator and 1” bundle for testing. Bundle development was conducted with a Development Spin Unit (DSU) installed at MEDAL. Air Liquide’s cold membrane technology was demonstrated with real coal fired flue gas at the National Carbon Capture Center (NCCC) with a 0.3 MWe field-test unit (FTU). The FTU was designed to incorporate testing of two PI-1 commercial membrane bundles (12” or 6” diameter) in parallel or series. A slip stream was sent to the next generation PI-2 membrane for testing with real flue gas. The system exceeded performance targets with stable PI-1 membrane operation for over 500 hours of single bundle, steady state testing. The 12” PI-1 bundle exceeded the productivity target by achieving ~600 Nm3/hr, where the target was set at ~455 Nm3/hr at 90% capture rate. The cost of 90% CO 2 capture from a 550 MWe net coal power plant was estimated between 40 and $45/tonne. A 6” PI-1 bundle exhibited superior bundle performance compared to the 12” PI-1 bundle. However, the carbon capture cost was not lower with the 6” PI-1 bundle due to the higher bundle installed cost. A 1” PI-1 bundle was tested to compare bundles with different length / diameter ratios. This bundle exhibited the lowest performance due to the different fiber winding pattern and increased bundle non-ideality. Several long-term and parametric tests were conducted with 3,200 hours of total run-time at NCCC. Finally, the new PI-2 membrane fiber was tested at a small scale (1” modules) in real flue gas and exhibited up to 10 times the CO 2 permeance and slightly lower CO 2 /N 2 selectivity as the commercial PI-1 fiber. This corresponded to a projected 4 - 5 times increase in the productivity per bundle and a potential cost reduction of $3/tonne for CO2 capture, as compared with PI-1. An analytical campaign was conducted to trace different impurities such as NOx, mercury, Arsenic, Selenium in gas and liquid samples through the carbon capture system. An Environmental, Health and Safety (EH&S) analysis was completed to estimate emissions from a 550 MWe net power plant with carbon capture using cold membrane. A preliminary design and cost analysis was completed for 550 tpd (~25 MWe) plant to assess the capital investment and carbon capture cost for PI-1 and PI-2 membrane solutions from coal fired flue gas. A comparison was made with an amine based solution with significant cost advantage for the membrane at this scale. Additional preliminary design and cost analysis was completed between coal, natural gas and SMR flue gas for carbon capture at 550 tpd (~25 MWe) plant.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Distributed out-of-memory NMF on CPU/GPU architectures

We propose an efficient distributed out-of-memory implementation of the non-negative matrix factorization (NMF) algorithm for heterogeneous high-performance-computing systems. The proposed implementation is based on prior work on NMFk, which can perform automatic model selection and extract latent variables and patterns from data. In this work, we extend NMFk by adding support for dense and sparse matrix operation on multi-node, multi-GPU systems. The resulting algorithm is optimized for out-of-memory problems where the memory required to factorize a given matrix is greater than the available GPU memory. Memory complexity is reduced by batching/tiling strategies, and sparse and dense matrix operations are significantly accelerated with GPU cores (or tensor cores when available). Input/output latency associated with batch copies between host and device is hidden using CUDA streams to overlap data transfers and compute asynchronously, and latency associated with collective communications (both intra-node and inter-node) is reduced using optimized NVIDIA Collective Communication Library (NCCL) based communicators. Benchmark results show significant improvement, from 32X to 76x speedup, with the new implementation using GPUs over the CPU-based NMFk. Good weak scaling was demonstrated on up to 4096 multi-GPU cluster nodes with approximately 25,000 GPUs when decomposing a dense 340 Terabyte-size matrix and an 11 Exabyte-size sparse matrix of density 10 -6 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Evaluation of Graph Analytics Frameworks Using the GAP Benchmark Suite

The analysis of connected data is an increasingly important application in high-performance computing. Such analyses can reveal fraudulent patterns in financial transactions, optimize telecommunications networks, predict information flow in social networks, etc. However, the landscape of graph analytics is highly diverse. Graph algorithms stress processor architectures differently, and no one graph can represent all topologies. Consequently, no single approach or framework is expected to be optimal for all graph analytics problems. To help make sense of this diverse landscape, we evaluated four approaches to graph analytics: GraphBLAS, Galois, BGL17, GraphIt; and compare them against hand-tuned implementations that take advantage of hardware features on our test platform. Graph- BLAS formulates graph analytics as sparse linear algebra. Galois provides syntactic constructs for data parallelism over irregular data structures. BGL17 is a generic C++ template library for implementing graph algorithms. GraphIt provides a domain- specific language to describe and optimize graph algorithms. We use the GAP Benchmark Suite to establish baseline performance and guide the side-by-side evaluation of each framework. GAP consists of 30 tests: six graph analytics algorithms (breadth- first search, single-source shortest path, PageRank, betweenness centrality, connected components, and triangle counting) run on five graphs, each with different topological characteristics (e.g., high diameter, skewed degree distribution, high average degree). High-performance reference implementations are included for each benchmark algorithm. Because a graph can be loaded into memory a number of ways (e.g., flat file on disk, compressed sparse format, data frames, retrieved from SQL or NoSQL databases), our evaluation focused on computational performance rather than I/O. Our results show the relative strengths of each framework.

Graph algorithms, Benchmarking, shared-memory prog↗