Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed parallel programming”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

FINAL REPORT - Development of High Pressure Hydrogen Storage Tank for Storage and Gaseous Truck Delivery

The “Development of High Pressure Hydrogen Storage Tanks for Storage and Gaseous Truck Delivery” project [DE-FG36-08GO18062] was initiated on 01 July 2008. Hexagon Lincoln (then Lincoln Composites) received grant funding from the U.S. Department of Energy to support the design and development of an improved bulk hauling and storage solution for hydrogen in terms of cost, safety, weight and volumetric efficiency. The development of this capability required parallel development and qualification of large all-composites pressure vessels, a custom ISO container to transport and store said tanks, and performance of trade studies to identify optimal operating pressure for the system. Qualification of the 250 bar TITAN® module was completed in 2009 with supervision from the American Bureau of Shipping [ABS], and the equipment has been used internationally for bulk transportation of fuel gases since 2010. Phase 1 of the project was successfully completed in 2012 with the issuance of USDOT SP 14951, the special permit authorizing the manufacture, marking, sale and use of TITAN® Mobile Pipeline® equipment in the United States. The introduction of tube trailers with light weight composite tankage has meant that 2 to 3 times as much gaseous fuel can be transported with each trip. This increased hauling efficiency offers dramatically reduced operating costs and has enabled a profitable business model for over-the-road compressed natural gas delivery. The economic drivers of this business opportunity vary from country to country and region to region, but in many places gas distribution companies have realized profitable operations. Additional testing was performed in 2015 to characterize hydrogen-specific operating protocols for use of TITAN® systems in CHG service at 250 bar. This program demonstrated that existing compression and decompression methodologies can efficiently and safely fill and unload lightweight bulk hauling systems. Hexagon Lincoln and U.S. DOE agreed to continue into Phase 2 of the project without pursuing the development of higher pressure capabilities as originally planned. At 250 bar, development of equipment for hydrogen transport is supported by strong activity in the adjacent natural gas transportation sector. Trade studies performed since 2011 indicate optimization of hauling efficiency and system cost for hydrogen transport at about 350 bar (5076 psi). However, due to reduced efficiency of compression of natural gas above 250 bar, 350 bar operation is not an attractive option for natural gas transportation. The CHG market is not developed at this time, and it is difficult to forecast the arrival of significant revenues. On the investment side, the cost to fully qualify a large tank module at 350 bar is estimated at $3MM to $5MM. There is insufficient CHG market definition to support a stand-alone business case for this investment without near term revenue in the adjacent CNG transportation market. Therefore development of a 350 bar TITAN® system was deferred and not pursued under this project. Hexagon Lincoln continues to support the development of tankage and equipment for operation at 350 bar and above; with 700 bar vehicle tanks and 950 bar tanks for ground storage applications. Phase 2 activities were focused on reducing system cost, increasing system capacity, increasing system safety and characterization of polymer material performance specific to hydrogen pressure vessel usage. With the successful launch of TITAN® modules and trailers in natural gas transportation, over 600 units have been produced through the end of 2016, resulting in improved purchasing power for raw materials and manufactured components. This has allowed Hexagon Lincoln to approach the current project goals for system cost. At $590/kg of compressed hydrogen delivered, the system cost of the baseline TITAN® module is below the project’s 2015 target of $730/kg H2 delivered, and very close to the project’s 2020 target of $575/kg H2 delivered. [Based on product pricing in 1Q2017.] Emphasis was placed on configuration of larger capacity systems within the vehicle weights and dimensions allowed on federal and state highways in the United States and other countries. These activities resulted in the design and development of integrated tube trailer systems that have increased delivery capacities by 45%. The hydrogen delivery capacity of our largest system is 845 kg, exceeding the project’s 2015 target of 700 kg H2 delivered. Emerging technologies offering improvement of the safety systems used on the equipment were investigated, with particular focus on improving the reliability and cost of the emergency venting system for fire protection. Finally, investment in our materials laboratory improved detection and characterization of hydrogen-induced damage in polymer materials, supporting the development of operational protocols to avoid damage to pressure vessel liners and valve components.

08 HYDROGEN↗
DOE OSTI2019

Proactive Data Containers for Scientific Storage (Final Report)

Emerging HPC systems are expected to be deployed with an unprecedented level of complexity, due to a deep system memory/storage hierarchy and heterogeneity of the storage hardware. This hierarchy is expected to range from CPU cache through several levels of volatile memory to nonvolatile memory, traditional hard disks, and tape. Simple and efficient methods of data management and movement through this hierarchy is critical for scientific applications using exascale systems. Existing storage system and I/O (SSIO) technologies face severe challenges in dealing with these requirements. POSIX and MPI I/O standards that are the basis for existing I/O libraries and parallel file systems present fundamental challenges in the areas of scalable metadata operations, semantics-based data movement performance tuning, asynchronous operation, and support for scalable consistency of distributed operations. Moving toward new paradigms for SSIO in the extreme-scale era, we have proposed to investigate novel object-based data abstractions and storage mechanisms that take advantage of the deep storage hierarchy and enable proactive automated performance tuning. In order to achieve these overarching goals, we initiated an effort to develop a fundamental new data abstraction, called Proactive Data Containers (PDC). A PDC is a container within a locus of storage (memory, NVRAM, disk, etc.) that stores science data in an object-centric manner. Managing data as objects enables powerful optimization opportunities for data movement and transformations. The R&D focus of this project are: 1) formulation of object-oriented PDCs and their mapping in different levels of the exascale storage hierarchy; 2) efficient strategies for moving data in deep storage hierarchies using PDCs; 3) techniques for transforming and reorganizing data based on application requirements; and 4) novel analysis paradigms for enabling data transformations and user-defined analysis on data in PDCs. Toward achieving these overarching goals, we designed an object-centric application programing interface (API) for HPC, scalable metadata management for object-centric storage systems, and data movement optimizations such as Data Elevator for moving data between two levels of storage devices and TAPIOCA for efficient aggregation of data on compute nodes. We then implemented several components of the PDC system. They include metadata management, data placement services, remote procedure calls, data aggregation, etc. We have put them together into the overall PDC framework.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Proactive Operations and Investment Planning via Stochastic Optimization to Enhance Power Systems’ Extreme Weather Resilience

We present scalable stochastic optimization approaches for improving power systems’ resilience to extreme weather events. We consider both proactive redispatch and transmission line hardening as alternatives for mitigating expected load shed due to extreme weather, resulting in large-scale stochastic linear programs (LPs) and mixed-integer linear programs (MILPs). We solve these stochastic optimization problems with progressive hedging (PH), a parallel, scenario-based decomposition algorithm. Our computational experiments indicate that our proposed method for enhancing power system resilience can provide high-quality solutions efficiently. With up to 128 scenarios on a 2,000-bus network, the operations (redispatch) and investment (hardening) resilience problems can be solved in approximately 6 min and 2 h of wall-clock time, respectively. Additionally, we solve the investment problems with up to 512 scenarios, demonstrating that the approach scales very well with the number of scenarios. Moreover, the method produces high quality solutions that result in statistically significant reductions in expected load shed. Our proposed approach can be augmented to incorporate a variety of other operational and investment resilience strategies, or a combination of such strategies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

HTGR Multiphysics Application Drivers FY26 Updates

This report summarizes FY26 progress under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program's high-temperature gas-cooled reactor (HTGR) application driver work, covering a wide range of activities such as code validation and multi-physics code assessment. 1) A detailed SAM model of the High-Temperature Engineering Test Reactor (HTTR) was developed using a unique-block grouping approach, with an extended parallel thermal network method to capture block-to-block conduction and radiation heat transfer, and applied to steady-state simulations of the HTTR 30~MW and 9~MW cases. 2) In another activity, SAM's newly implemented multi-component gas flow model was validated against the Natural convection Shutdown heat removal Test Facility (NSTF) argon ingress experiment, correctly capturing the density-driven suppression and thermal recovery of natural circulation observed when argon is introduced into the air-cooled Reactor Cavity Cooling System (RCCS) loop. 3) For the OECD/NEA High Temperature Test Facility (HTTF) benchmark, we co-led the international benchmark activities as well as the OECD/NEA final benchmark report to be released at the end of this year. 4) Finally, the coupled Griffin-SAM modeling capability for pebble-bed HTGRs was advanced by verifying the Griffin neutronics solution against Serpent Monte Carlo for a realistic non-uniform temperature distribution, resolving several deficiencies in the SAM-to-Griffin temperature transfer scheme, and enabling distinct fuel kernel, moderator, and coolant temperatures for cross section feedback. These new features were demonstrated in a PBR load-following transient.

Lee, Alvin↗
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 OSTI2015

SU-E-T-778: Use of the 2D MatriXX Detector for Measuring Scanned Ion Beam Parameters

Purpose: The quality assurance (QA) procedure has to check the most relevant beam parameters to ensure the delivery of the correct dose to patients. Film dosimetry, which is commonly used for scanned ion beam QA, does not provide immediate results. The purpose of this work is to answer whether, for scanned ion beam therapy, film dosimetry can be replaced with the 2D MatriXX detector as a real-time tool. Methods: MatriXX, equipped with 32×32 parallel plate ion-chambers, is a commercial device intended for pre-treatment verification of conventional radiation therapy.The MatriXX, placed at the isocenter, and GAFCHROMIC films, positioned on the MatriXX entrance, were exposed to 131.44 MeV proton and 221.45 MeV/u Carbon-ion beams.The OmniPro-I’mRT software, applied for the data taking of MatriXX, gives the possibility of acquiring consecutive snapshots. Using the NI LabVIEW, the data from snapshots were logged as text files for further analysis. Radiochromic films were scanned with EPSON scanner and analyzed using software programs developed in-house for comparative purposes. Results: The field dose uniformity, flatness, beam position and beam width were investigated. The field flatness for the region covering 6×6 cm{sup 2} square field was found to be better than 2%. The relative standard deviations, expected to be constant over 2×2, 4×4 and 6×6 pixels from MatriXX measurement gives a uniformity of 1.5% in good agreement with the film results.The beam center position is determined with a resolution better than 200 µm for Carbon and less than 100 µm for proton beam.The FWHM determination for a beam wider than 10 mm is satisfactory, whilst for smaller beams the determination is uncertain. Conclusion: Precise beam position and fast 2D dose distribution can be determined in real-time using MatriXX detector. The results show that MatriXX is quick and accurate enough to be used in charged-particle therapy QA.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

PLATO (Parallel Load Assignment Tool): A Parallel Workload Partitioner for Particle-based Methods

A key first step in many large-scale high-performance scientific computing applications is the distribution of the computational domain across all participating processors of a parallel computing platform such as a shared memory machine, cluster or a tightly-coupled supercomputer. The primary goal of such a partitioning step is to ensure that the computational work performed by each process remains load-balanced throughout the execution of the parallel application. In this project, a general-purpose parallel partitioning tool, code-named PLATO, was designed and implemented. PLATO assumes that the workload is a set of points in an n-dimensional metric space. These points are abstractions of computational units that may represent point-like computational units such as celestial objects in astrophysics, atoms in materials sciences and molecules in biology/chemistry, or may represent more abstract computational units such as grid points in finite difference methods, grid elements in finite element methods or n-dimensional points in the feature space of a host of other applications. PLATO supports several key and useful capabilities, namely: (a) multiple load partitioning strategies, (b) support for arbitrary user-defined processor topology, (c) pre-fetching data from neighboring processors based on user-defined parameters, and (d) support for arbitrary number and types of attributes associated with each point in the computational domain. PLATO is implemented in the C programming language using Message Passing Interface (MPI) for inter-processor communications and its parallel performance tested on up to 48 processors of a shared memory machine. This report includes a detailed performance study of PLATO using parallel speedup and iso-efficiency as the two main metrics of parallel performance. In addition, a user manual and DOXYGEN generated code documentation have been included for completeness.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2013

Hypersonic Inflatable Aerodynamic Decelerator (HIAD) Technology Development Overview

The successful flight of the Inflatable Reentry Vehicle Experiment (IRVE)-3 has further demonstrated the potential value of Hypersonic Inflatable Aerodynamic Decelerator (HIAD) technology. This technology development effort is funded by NASA's Space Technology Mission Directorate (STMD) Game Changing Development Program (GCDP). This paper provides an overview of a multi-year HIAD technology development effort, detailing the projects completed to date and the additional testing planned for the future. The effort was divided into three areas: Flexible Systems Development (FSD), Mission Advanced Entry Concepts (AEC), and Flight Validation. FSD consists of a Flexible Thermal Protection Systems (FTPS) element, which is investigating high temperature materials, coatings, and additives for use in the bladder, insulator, and heat shield layers; and an Inflatable Structures (IS) element which includes manufacture and testing (laboratory and wind tunnel) of inflatable structures and their associated structural elements. AEC consists of the Mission Applications element developing concepts (including payload interfaces) for missions at multiple destinations for the purpose of demonstrating the benefits and need for the HIAD technology as well as the Next Generation Subsystems element. Ground test development has been pursued in parallel with the Flight Validation IRVE-3 flight test. A larger scale (6m diameter) HIAD inflatable structure was constructed and aerodynamically tested in the National Full-scale Aerodynamics Complex (NFAC) 40ft by 80ft test section along with a duplicate of the IRVE-3 3m article. Both the 6m and 3m articles were tested with instrumented aerodynamic covers which incorporated an array of pressure taps to capture surface pressure distribution to validate Computational Fluid Dynamics (CFD) model predictions of surface pressure distribution. The 3m article also had a duplicate IRVE-3 Thermal Protection System (TPS) to test in addition to testing with the Aerocover configuration. Both the Aerocovers and the TPS were populated with high contrast targets so that photogrammetric solutions of the loaded surface could be created. These solutions both refined the aerodynamic shape for CFD modeling and provided a deformed shape to validate structural Finite Element Analysis (FEA) models. Extensive aerothermal testing has been performed on the TPS candidates. This testing has been conducted in several facilities across the country. The majority of the testing has been conducted in the Boeing Large Core Arc Tunnel (LCAT). HIAD is continuing to mature testing methodology in this facility and is developing new test sample fixtures and control methodologies to improve understanding and quality of the environments to which the samples are subjected. Additional testing has been and continues to be performed in the NASA LaRC 8ft High Temperature Tunnel, where samples up to 2ft by 2ft are being tested over representative underlying structures incorporating construction features such as sewn seams and through-thickness quilting. With the successful completion to the IRVE-3 flight demonstration, mission planning efforts are ramping up on the development of the HIAD Earth Atmospheric Reenty Test (HEART) which will demonstrate a relevant scale vehicle in relevant environments via a large-scale aeroshell (approximately 8.5m) entering at orbital velocity (approximately 7km/sec) with an entry mass on the order of 4MT. Also, the Build to Print (BTP) hardware built as a risk mitigation for the IRVE-3 project to have a "spare" ready to go in the event of a launch vehicle delivery failure is now available for an additional sub-orbital flight experiment. Mission planning is underway to define a mission that can utilize this existing hardware and help the HIAD project further mature this technology.

Hughes, Stephen J.↗
DOE OSTI2017

MPI-FAUN: An MPI-Based Framework for Alternating-Updating Nonnegative Matrix Factorization

Non-negative matrix factorization (NMF) is the problem of determining two non-negative low rank factors W and H, for the given input matrix A, such that A≈WH. NMF is a useful tool for many applications in different domains such as topic modeling in text mining, background separation in video analysis, and community detection in social networks. Despite its popularity in the data mining community, there is a lack of efficient parallel algorithms to solve the problem for big data sets. The main contribution of this work is a new, high-performance parallel computational framework for a broad class of NMF algorithms that iteratively solves alternating non-negative least squares (NLS) subproblems for W and H. It maintains the data and factor matrices in memory (distributed across processors), uses MPI for interprocessor communication, and, in the dense case, provably minimizes communication costs (under mild assumptions). The framework is flexible and able to leverage a variety of NMF and NLS algorithms, including Multiplicative Update, Hierarchical Alternating Least Squares, and Block Principal Pivoting. Our implementation allows us to benchmark and compare different algorithms on massive dense and sparse data matrices of size that spans from few hundreds of millions to billions. We demonstrate the scalability of our algorithm and compare it with baseline implementations, showing significant performance improvements. The code and the datasets used for conducting the experiments are available online.

2D↗
NASA NTRS1997

Effects of Fuel Preheat on Soot Formation in Microgravity Laminar Diffusion Flames

Nonbuoyant flames offer themselves as an attractive and promising platform to gain a better understanding of soot mechanisms. The effects of buoyancy can be eliminated temporarily in drop towers which sustain brief intervals of reduced gravity-typically lower than 10(exp -3)g- extending up to several seconds at a time. Microgravity facilities have been employed to show that nonbuoyant flames are longer, wider and sootier than their normal-gravity counterparts. Sunderland et al. recently verified the existence of smoke point in laminar nonbuoyant flames. As reported, microgravity flames operating above their smoke point displayed a blunt tip and much broader soot-containing regimes in comparison to their buoyant counterparts. Mortazavi et al. established that residence times in microgravity laminar jet diffusion flames with Re=0(100) tend to be proportional to burner diameter and inversely proportional to burner exit velocity. This offers the capability to alter residence times in nonbuoyant laminar jet diffusion flames when varying the burner exit diameters and velocities. Megaridis et al. presented a quantitative definition of the soot-field structure within laminar microgravity jet diffusion flames which operated well above their smoke point. The experimental methodology involved a full-field laser-light extinction technique and jet diffusion flames of nitrogen-diluted (50% vol.) acetylene fuel burning in quiescent air at atmospheric pressure. The work was conducted at the 2.2s drop tower of the NASA Lewis Research Center (NASA-LeRC). Parallel work on 1-g flames was also presented in (6) to facilitate comparisons on the effect of gravity on the soot fields. As reported, the soot spatial distributions in 0-g flames did not change in a detectable manner after 1s within a typical 2.2s experiment. During that period, the soot field was shown to sustain a pronounced annular structure throughout the luminous nonbuoyant-flame zone. The maximum soot volume fraction measured at 0-g was nearly a factor of two higher than that at 1-g, thus confirming the enhanced sooting tendency of nonbuoyant flames. Greenberg and Ku presented a similar study and reported trends that matched those of for the 50% (vol.) nitrogen-diluted acetylene fuel. Furthermore, they examined pure acetylene flames and reported similar trends with respect to the influence of gravity on maximum soot volume fractions and flame cross-section-averaged soot loadings. Both studies clearly demonstrated the improved spatial resolution of microgravity flames compared to their normal-gravity counterparts. The current study evaluates the influence of moderate fuel preheat on soot formation within 0-g laminar gas jet diffusion flames. While fuel temperature variations have little influence on residence times in 1-g, they have a much more significant effect in 0-g. The primary objective of this program is to quantify this effect and its consequences on sooting by comparing soot volume fraction distributions under preheated and unpreheated-fuel conditions. Furthermore, the current work aims at expanding the limited soot database available for nonbuoyant flames. Soot fields in such flames can be used to perform additional tests of recently developed soot sub-models which have the potential to become powerful predicting tools in combustion design.

Konsur, Bogdan↗
DOE OSTI2024

Using Parameter Sweep in WaterTAP to Analyze New Water Treatment Technologies

We describe a powerful and generalized parameter sweep tool in this report that was originally developed to analyze the performance of existing and novel water treatment models being developed in WaterTAP. Since WaterTAP is built upon IDAES and Pyomo, the parameter sweep tool can be used to systematically explore and debug the behavior of most Pyomo and IDAES numerical models. In order to enable meaningful analyses, the parameter sweep tool has been designed with the following features: 1) Model flexibility: The parameter sweep tool does not enforce any restrictions on the types of models that can be used with it. As long as a Pyomo model can be solved and the parameter is active and mutable, the tool only needs functions that describe how to run the model, the sweep parameters, and the output quantities of interest. 2) Flexible sampling: The parameter sweep tool has inbuilt functions to generate samples from a random distribution or a multidimensional Euclidean space. Furthermore, the users have to ability to supply samples generated from a tool of their choice. 3) Multiple sweep types: A user can choose from one of 3 types of parameter sweeps depending on their needs. 4) Detailed outputs: Outputs generated by the parameter sweep tool can be stored in detailed H5 file or user-friendly CSV files for post processing. 5) Parallel computing: The parameter sweep supports shared and distributed memory parallel computing to enable the use of high performance computers (HPC) for large-scale analyses. 6) Modular: The parameter sweep tool is self-contained and can easily be integrated within an outer-loop analysis or as desired by the user. 7) Ease of use: The tool is well documented and a simple sweep can be easily executed by following the online documentation in a few lines of code. We demonstrate the use of the parameter sweep tool on a simple water treatment system from the WaterTAP repository and show its parallel scaling performance on an Apple laptop and NREL's Eagle HPC. The parameter sweep tool is actively being used with models currently being developed within WaterTAP and we expect its use to grow beyond it to other IDAES and Pyomo models.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2000

Global Interoperability of High Definition Video Streams Via ACTS and Intelsat

In 1993, a proposal at the Japan.-U.S. Cooperation in Space Program Workshop lead to a subsequent series of satellite communications experiments and demonstrations, under the title of Trans-Pacific High Data Rate Satellite Communications Experiments. The first of which is a joint collaboration between government and industry teams in the United States and Japan that successfully demonstrated distributed high definition video (HDV) post-production on a global scale using a combination of high data rate satellites and terrestrial fiber optic asynchronous transfer mode (ATM) networks. The HDV experiment is the first GIBN experiment to establish a dual-hop broadband satellite link for the transmission of digital HDV over ATM. This paper describes the team's effort in using the NASA Advanced Communications Technology Satellite (ACTS) at rates up to OC-3 (155 Mbps) between Los Angeles and Honolulu, and using Intelsat at rates up to DS-3 (45 Mbps) between Kapolei and Tokyo, with which HDV source material was transmitted between Sony Pictures High Definition Center (SPHDC) in Los Angeles and Sony Visual Communication Center (VCC) in Shinagawa, Tokyo. The global-scale connection also used terrestrial networks in Japan, the States of Hawaii and California. The 1.2 Gbps digital HDV stream was compressed down to 22.5 Mbps using a proprietary Mitsubishi MPEG-2 codec that was ATM AAL-5 compatible. The codec: employed four-way parallel processing. Improved versions of the codec are now commercially available. The successful post-production activity performed in Tokyo with a HDV clip transmitted from Los Angeles was predicated on the seamless interoperation of all the equipment between the sites, and was an exciting example in deploying a global-scale information infrastructure involving a combination of broadband satellites and terrestrial fiber optic networks. Correlation of atmospheric effects with cell loss, codec drop-out, and picture quality were made. Current efforts in the Trans-Pacific series plan to examine the use of Internet Protocol (IP)-related technologies over such an infrastructure. The use of IP allows the general public to be an integral part of the exciting activities, helps to examine issues in constructing the solar-system internet, and affords an opportunity to tap the research results from the (reliable) multicast and distributed systems communities. The current Trans- Pacific projects, including remote astronomy and digital library (visible human) are briefly described.

Hsu, Eddie↗
NASA NTRS1997

Determination of Stress Intensity Factor Distributions for "Interface" Cracks in Incompressible, Dissimilar Materials

The present study was undertaken in order to develop test methods and procedures for measuring the variation of the stress intensity factor through the thickness in bimaterial specimens containing cracks within and parallel to the bond line using the frozen stress photoelastic method. Since stress freezing materials are incompressible above critical temperature, and since thick plates are to be employed which tend to produce a state of plane strain near the crack tip, the interface near tip fracture equations reduce to the classic form for homogeneous materials. Moreover, zero thickness interfaces do not exist when materials are bonded together. It was decided early on that it would be important to insure a uniform straight and accurate crack tip region through the thickness of the body to reduce scatter in the SIF distribution through the thickness. It was also observed that rubberlike materials which were desired to be modeled exhibited significant tip blunting prior to crack extension and that some blunting of the tip would provide a more realistic model. It should be noted that, in normal stress freezing photoelastic work, it is considered good practice to avoid utilizing data near bond lines in photoelastic models due to the bond line stresses which inevitably develop when two parts are bonded together. Thus, the present study involves certain exploratory aspects in deviating from standard practice in stress freezing work. With the above ideas in mind, several different test methods were investigated and are described in the following sections and appendices. The geometry selected for the program was a thick, edge cracked specimen containing a bond line.

Smith, C. W.↗
DOE OSTI2019

Cataloging the visible universe through Bayesian inference in Julia at petascale

A key task in astronomy is to locate astronomical objects in images and to characterize them according to physical parameters such as brightness, color, and morphology. This task, known as cataloging, is challenging for several reasons: many astronomical objects are much dimmer than the sky background, labeled data is generally unavailable, overlapping astronomical objects must be resolved collectively, and the datasets are enormous – terabytes now, petabytes soon. In this work, we infer an astronomical catalog from 55 TB of imaging data using Celeste, a Bayesian variational inference code written entirely in the high-productivity programming language Julia. Using over 1.3 million threads on 650,000 Intel Xeon Phi cores of the Cori Phase II supercomputer, Celeste achieves a peak rate of 1.54 DP PFLOP/s. Celeste is able to jointly optimize parameters for 188 M stars and galaxies, loading and processing 178 TB across 8192 nodes in 14.6 min. To achieve this, Celeste exploits parallelism at multiple levels (cluster, node, and thread) and accelerates I/O through Cori’s burst buffer. Julia’s native performance enables Celeste to employ high-level constructs without resorting to hand-written or generated low-level code (C/C++/Fortran) while still achieving petascale performance.

Regier, Jeffrey↗
NASA NTRS2007

Modeling of Radiative Transfer in Protostellar Disks

This program implements a spectral line, radiative transfer tool for interpreting Spitzer Space Telescope observations by matching them with models of protostellar disks for improved understanding of planet and star formation. The Spitzer Space Telescope detects gas phase molecules in the infrared spectra of protostellar disks, with spectral lines carrying information on the chemical composition of the material from which planets form. Input to the software includes chemical models developed at JPL. The products are synthetic images and spectra for comparison with Spitzer measurements. Radiative transfer in a protostellar disk is primarily affected by absorption and emission processes in the dust and in molecular gases such as H2, CO, and HCO. The magnitude of the optical absorption and emission is determined by the population of the electronic, vibrational, and rotational energy levels. The population of the molecular level is in turn determined by the intensity of the radiation field. Therefore, the intensity of the radiation field and the population of the molecular levels are inter-dependent quantities. To meet the computational challenges of solving for the coupled radiation field and electronic level populations in disks having wide ranges of optical depths and spatial scales, the tool runs in parallel on the JPL Dell Cluster supercomputer with C++ and Fortran compiler with a Message Passing Interface. Because this software has been developed on a distributed computing platform, the modeling of systems previously beyond the reach of available computational resources is possible.

VonAllmen, Paul↗
DOE OSTI2027

Nodal capacity expansion planning with flexible large-scale load siting

We propose explicitly incorporating large-scale load siting into a stochastic nodal power system capacity expansion planning model that concurrently co-optimizes generation, transmission, and storage expansion. The potential operational flexibility of some of these large loads is also taken into account by considering them as consisting of a set of tranches with different reliability requirements, which are modeled as a constraint on expected served energy across operational scenarios. We implement our model as a two-stage stochastic mixed-integer optimization problem with cross-scenario expectation constraints. To overcome the challenge of scalability, we build upon existing work to implement this model on a high performance computing platform and exploit scenario parallelization using an augmented Progressive Hedging Algorithm. The algorithm is implemented using the bounding features of mpisppy, which have shown to provide satisfactory provable optimality gaps despite the absence of theoretical guarantees of convergence. We test our approach and assess the value of this proactive planning framework on total system cost and reliability metrics using realistic testcases geographically assigned to San Diego and South Carolina, with datacenter and direct air capture facilities as large loads.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2015

The impact of aqueous medium on zeolite framework integrity

Understanding the zeolite framework stability in aqueous phase is crucial to develop stable catalysts. Al K–edge, extended X–ray absorption fine structure and 27Al MAS NMR spectroscopies in combination with DFT calculations have been used to monitor both qualitative and quantitative structural changes of two well–characterized samples with BEA structure. The effects of various properties on stability were explored, including Al concentration, Al distribution, particle size and structural defects. As the samples were degraded by treatment in hot liquid water, the local structure about the Al T–site remained mostly intact, including the Al–O–Si angles and bond distances, while the nano–scale crystalline structure as measured by XRD and TEM was disrupted. The combined data suggest a three–step mechanism in which, initially, the HBEA framework crystallinity decreases via hydrolysis of T–O bonds along polymorph stacking faults and inter–grain boundaries in a mode similar to crack propagation in glass. With prolonged exposure, amorphization occurs via hydrolysis of surface Si–OH groups propagating inward through the zeolite lattice. In parallel, cracks propagate within the crystalline micro–domains along paths through specific T–O–T groups. Authors thank B. W. Arey (PNNL) for HIM measurements, T. Huthwelker for support during Al XAFS measurements at the Swiss Light Source (PSI, Switzerland) and M. Y. Hu (PNNL) for support during NMR experiments. This work was supported by the U. S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences. NMR experiments were performed at the Environmental Molecular Science Laboratory, a national scientific user facility sponsored by the DOE’s Office of Biological and Environmental Research, and Physical Science Laboratory both located at Pacific Northwest National Laboratory (PNNL). PNNL is a multi–program national laboratory operated for DOE by Battelle under Contract No. DE–AC05–76RL01830.

Environmental Molecular Sciences Laboratory↗
NASA NTRS1997

Mid-Infrared Observational and Theoretical Studies of Star Formation and Early Solar Systems

The first 2 years of this program were used to make mid-IR observations of regions of star formation in the Orion nebula with the UCSD mid-IR camera at the UCSD/University of Minnesota telescope at Mt. Lemmon. These observations attempted to make the first systematic study of an extended region, known to have newly forming stars, and expected to have complex mid-IR emission. We discovered, to our surprise, that most of the thermal emission originated from extended sources rather than from point sources. This interesting observation made the analysis of the data much more complex, since the chop/nod procedures used at these wavelengths produce a differential measurement of the emission in one region compared to that in the adjacent region. Disentangling complex extended emission in such a situation is very difficult. In parallel with this work we were also observing comets in the thermal infrared, the other component of the original proposal. Some spectacular data on the comet Swift-Tuttle was acquired and published. A changing jet structure observed over a 2 week period is described. The rotation period of the comet can be measured at 66 hours. The size of the nucleus can also be estimated (at 30 km) from the observed excess flux from the nucleus. These data have lead to the development of models describing the action of dust particles of differing sizes and composition leaving the nucleus. The spatial distribution of the predicted IR emission has been compared to the observed jet structures, leading to estimates of both particles sizes, relative amounts of silicate vs organic grains, and the amounts of dust emitted in the jets vs isotopic emission.

Jones, Barbara↗