Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bundle method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Field Test Report Neutron Scintillator Array Dry Storage Cask Scanner FY2024

During two weeks of Field Testing at the Idaho National Laboratory INTEC Cask Farm in July and August 2024, the LLNL Dry Storage Cask Scanner Array was lifted on top of an MC-10 dry storage fuel cask and operated to acquire neutron and gamma-ray data from the 24 fuel bundle positions. Neutron and gamma-ray data acquisition scans across the top of the cask of varying dwell times were performed July 15-18, 2024 and August 19-22, 2024 to evaluate the ability of the scanner data to reveal asymmetries in the fuel positions that reflect asymmetries in the MC-10 cask fuel bundle loading. The MC-10 cask 24 position fuel bundle loading at the INTEC Cask Farm is well documented, including the locations of six empty fuel bundle positions. This loading presents an opportunity to test the ability of the scanner system to detect diversion of spent fuel bundles as well as to validate the MC-10 cask MCNP modeling. The cask scanner array consists of six Stilbene crystal scintillator detectors and a linear actuator frame that moves the six detectors across the MC-10 dry storage cask to obtain data above each of the 24 fuel bundle positions. The detectors are connected to a pulse-shape discrimination data acquisition system capable of generating separate neutron and gamma-ray spectra for each detector and for each scan position. From the prior single detector Field Test in 2021 and iteration with MCNP modeling, the neutron and gamma-ray data were analyzed in multiple energy regions to identify an analysis method that would provide the strongest and most consistent signature of the asymmetric MC-10 cask fuel loading1 . From both the 2021 Field Test and the current Field Test results, the neutron capture gamma-ray count rate around 2.2 MeV provides the strongest signature of the asymmetric MC-10 cask fuel loading and has qualitative agreement with MCNP calculations. Counting all gamma-rays produces a similar signature. Neutrons emerging from the cask top are moderated and captured by the hydrogen in the polyethylene moderator and scintillator detector, producing a 2.2 MeV gamma ray which is seen in the scintillator gamma-ray spectrum. The count rate in the 2.2 MeV gamma-ray region is ~50 c/s, which is ~1000x higher than the ~0.05 n/s rate in the > 4MeV neutron region, and ~50x greater than the ~1 n/s rate in the neutrons > 500 keV region. Analysis of the 2.2 MeV neutron-capture Compton-scattered gamma-rays produces a statistically significant signature of the INTEC Cask Farm MC-10 asymmetric fuel loading. MCNP simulations indicate that the average neutron energy spectrum offers the potential to detect a large asymmetry from several missing bundles as well as individual missing fuel bundles. Testing this feature will require measurements on a cask with single missing elements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Methane monitoring and detection apparatus and methods

A low cost, low power, passive optical methane monitoring system for fixed-position installation at oil and gas production well pads and gathering centers is disclosed. The optical methane monitoring system disclosed can be a scannable field of view Near Infrared (NIR) filter photometer to detect and quantify methane concentration in a two dimensional or a three dimensional grid above and around a facility. A randomized fiber optic bundle is disclosed that can be used to direct the total optical power from a collection lens to two or more isolated optical channels. Band pass filters isolate a desired wavelength range for transmission measurements for the two or more channels. Also disclosed is an absorption algorithm which accounts for variable background spectral intensity as well as correcting for water vapor and overall scattering effects to measure methane concentration for a given field of view.

Lievois, John↗
DOE OSTI2022

Second generation superconducting filaments and cable

A high-temperature superconducting filament and cable, and a method for manufacturing same. The substrate used to grow the superconducting layer is removed, and the exfoliated superconducting layer is coated with a protective layer, and then sliced into narrow strips. The strips are thereafter encapsulated with a conductive metal to provide a high-temperature superconducting filament. The filaments may be bundled together to provide a high-temperature superconducting cable.

Solovyov, Vyacheslav↗
DOE OSTI2026

A methodology for domain overlapping coupling of thermal-hydraulic systems

Multi-scale coupling has increasingly drawn attention as a promising approach for modeling thermal systems. Traditional system codes provide fast and robust predictions at the plant scale, while high-fidelity computational fluid dynamics (CFD)-based tools resolve localized flow and heat transfer phenomena with greater accuracy. By combining these complementary strengths, co-simulations enable multi-scale analysis that would otherwise be computationally prohibitive for a standalone CFD code. Here, this work introduces a robust and problem-agnostic domain overlapping (DO) coupling between the system thermal-hydraulic (STH) code System Analysis Module (SAM) and the coarse-mesh CFD code Pronghorn. Both applications belong to the Comprehensive Reactor Analysis Bundle (BlueCRAB) code suite, a code suite in active development at the Idaho National Laboratory (INL), tailored for multi-physics analysis of advanced reactors. Unlike previous approaches, BlueCRAB supports an agnostic interface between codes based on different fidelity, while its coupling formulation can address arbitrary flow geometries with multiple inlets and outlets in coupled components. The implemented method leads to consistent pressure drops, enthalpies, and scalar concentrations between coupled SAM and Pronghorn simulations. The methodology is demonstrated through two verification tests, which ensure the numerical consistency and conservation across the codes, and through one validation test against experimental data. The proposed problems explore different physical aspects inherent to thermal systems, with particular attention given to nuclear reactor analysis. These include buoyancy-driven flows, complex flow patterns, and setups with multiple inlets and outlets, representing challenges in advanced reactor applications.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Future Opportunities for LWR Irradiations in US Test Reactors

After several years of relatively low activity in the field of Light Water Reactor (LWR) fuel development, the Department of Energy again began to engage in developing new fuel technologies and irradiation performance data prompted by the Fukushima Daichi nuclear accidents. New competencies for irradiation testing in material test reactors in the United States began to be developed at this time using the Advanced Test Reactor (ATR), High Flux Isotope Reactor (HFIR), Massachusetts Institute of Technology Reactor (MITR), and the Transient Reactor Test Facility (TREAT). Capsules for testing fuel and cladding materials in ATR and HFIR were deployed, a Pressurized Water Reactor (PWR) condition loop for testing fuel rods was established in ATR, cladding corrosion studies were performed using a water loop in MITR, and TREAT pulse testing capabilities were commissioned for fuel rods in water capsules. The more recent and unexpected closure of the Halden Boiling Water Reactor (HBWR) also prompted further investments in Loss of Coolant Accident (LOCA) testing capabilities at TREAT. New configurations of these test devices show further potential in enhanced steam condition control and other investigations are building toward a flowing water loop for testing transient to dryout conditions. The closure of HBWR also prompted a major project currently underway to construct additional water loops in ATR where a novel approach is being pursued to enable Boiling Water Reactor (BWR) conditions. A meaningful collaborative project was awarded to MITR which, amidst an unexpected major overhaul of the reactor, has expanded cladding corrosion test capabilities at MITR. New explorations have led to methods for unique experiments at HFIR including channel box irradiations. New device developments are also bridging toward future potential for instrumented capsule irradiation tests in ATR and HFIR. Finally, a new project referred to as the System Physics Advanced Reactor Critical facility (SPARC) is gaining traction towards a large zero-power reactor able to produce physics validation data for LWR fuel bundle designs with increased enrichment and enhanced absorbers for 24-month operation cycles. This paper provides a brief summary of the status of these irradiation testbed capabilities with an emphasis on current efforts toward future capabilities to obtain new data and maximize the performance potential of LWR fuel technologies.

Woolstenhulme, Nicolas [Idaho National Laboratory ↗
DOE OSTI2026

Enamel nanocrystal misorientation increased with meat-eating and agriculture

Enamel covers teeth, is the hardest tissue in the vertebrate body and has a complex multiscale structure from nanometres to millimetres. The structure comprises thin, long hydroxyapatite (Ca 5 (PO 4 ) 3 OH) nanocrystals, 50–70 nm wide, many micrometres long, parallel and bundled into approximately 5-µm-wide rods. The rods undulate and cross into a microscale ‘decussation pattern’ that toughens enamel by deflecting cracks. However, the crystallographic orientation of enamel nanocrystals is poorly understood. Here we show that the misorientation angle of adjacent nanocrystals varies markedly across 12 primate teeth spanning 9 species, 17.8 million years of evolution and diverse diets. Using a method called Polarization Enabled Large Input of Crystal Angles at the Nanoscale (PELICAN), we compare nanocrystals in the same (pre)molar locations and show that misorientation increases with food hardness in extant and fossil non-human apes and monkeys. We compare misorientation across three major dietary shifts in human evolution: the transition to meat-eating about 2.0–1.5 million years before present, to agriculture (about 12,000 years before present), and the Industrial Revolution (about 250 years before present). We show that over the past 1.6 million years, in the human lineage misorientation increased with time, especially when meat and stone-ground grains were introduced into human diets, but not with the Industrial Revolution. Thus, besides macro-changes, teeth adapted to dietary change at the nanoscale and crystallographically. This observation suggests that misorientation may contribute to enamel’s resilience; thus, bioinspired materials may consider small misorientation angles for added resilience.

biomaterials↗
DOE OSTI2025

Origin of Pulsed Radio Emission from Magnetars

Extended periods of radio pulsations have been observed for six magnetars, displaying characteristics different from those of ordinary pulsars. In this Letter, we argue that radio emission is generated in a closed, twisted magnetic flux bundle originating near the magnetic pole and extending beyond 100 km from the magnetar. The electron–positron flow in the twisted bundle has to carry electric current and, at the same time, experiences a strong drag from the radiation field of the magnetar. This combination forces the plasma into a “radiatively locked” state with a sustained two-stream instability, generating radio emission. We demonstrate this mechanism using novel first-principles simulations that follow the plasma behavior by solving the relativistic Vlasov equation with the discontinuous Galerkin method. First, using one-dimensional simulations, we demonstrate how radiative drag induces the two-stream instability, sustaining turbulent electric fields. When extended to two dimensions, the system produces electromagnetic waves, including superluminal modes capable of escaping the magnetosphere. We measure their frequency and emitted power and incorporate the local simulation results into a global magnetospheric model. The model explains key features of the observed radio emission from magnetars: its appearance after an X-ray outburst, wide pulse profiles, luminosities ∼10 30 erg s −1 , and a broad range of frequencies extending up to ∼100 GHz.

Astronomical simulations↗
DOE OSTI2022

Pressure Drop in Seven-Pin Wire-Wrapped Rod Bundle for the Sodium Cartridge Loop in Versatile Test Reactor

This work studies the hydrodynamics of the seven pin wire-wrapped rod bundle in the sodium cartridge loop for the Versatile Test Reactor (VTR) through scaled water experiments and computational fluid dynamics (CFD) simulations. Here, the scaling analysis is first performed to demonstrate the hydrodynamic similarity between water and sodium flows at the same Reynolds number Re. A separate-effects test facility is designed and constructed based on the scaling analysis. Detailed experimental data on the pressure drop covering a wide range of Re values (1165 to 27 689) are obtained, which are used to evaluate existing correlations for friction factor and to benchmark CFD simulations. The experimentally determined friction factors f Exp agree well with the Upgraded Cheng and Todreas Detailed Correlation and Pacio-Chen-Todreas Detailed Model within ±7% but are significantly underpredicted by Rehme's correlation by 25%. Various CFD near-wall treatment methods are tested using ANSYS Fluent and evaluated by experimental data. It is found that when the recommended wall y + values are met, most of the near-wall treatment methods can give accurate friction factor predictions. The resolved near-wall method (y + ~ 1) with the Shear Stress Transport k - ω turbulence model and the scalable wall functions (y + > 12) with the realizable k - ε turbulence model can predict f Exp within ±5% The standard wall functions (y + > 30) and nonequilibrium wall functions (y + > 30) with the realizable k - ε model can predict f Exp within ±10%.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Automatic Calibration and Health Monitoring of Infrastructure Sensors

Smart transportation infrastructure relies on networks of heterogeneous sensors - cameras, radars, and lidars - continuously monitoring traffic conditions. However, executing the initial spatial calibration of multiple sensors and the subsequent health monitoring presents significant operational challenges. Environmental factors, mechanical vibrations, and gradual drift cause spatial misalignment, degrading fusion performance and tracking accuracy. Traditional calibration approaches require manual intervention with specialized targets or survey equipment, resulting in service interruptions and high maintenance costs. This work presents an automated framework for initial calibration and continuous health monitoring without human intervention or service disruption. Our approach addresses two critical problems: (1) detecting when sensors become miscalibrated during operation, and (2) automatically re-establishing spatial alignment using only operational traffic data. The health monitoring component analyzes measurement innovations - differences between sensor observations and predicted object states - to detect systematic biases indicative of calibration drift. By computing bias magnitude, directional consistency, and rejection rates, the system identifies miscalibrations as small as 0.5 meters. Unlike traditional methods requiring known calibration targets, our diagnostic operates continuously on live traffic observations, enabling early detection before fusion quality degrades. The automatic recalibration algorithm leverages overlapping sensor fields-of-view and temporal correlation of vehicle observations. Using graph-based optimization, the system automatically discovers which sensor pairs observe common regions, estimates pairwise spatial transformations using RANSAC-based robust estimation, and jointly optimizes all sensor poses through bundle adjustment. The framework handles practical deployment challenges, including different sensor sampling rates (1-10 Hz), varying installation positions, unknown orientations, and limited overlap regions (>10%). When approximate sensor positions are available from installation surveys (+/-1m accuracy), the algorithm additionally estimates sensor orientations, refining both position and rotation to sub-meter and sub-degree accuracy. We validate the framework on multi-hour traffic datasets from six heterogeneous sensors with sampling rates ranging from 1 Hz to 10 Hz. Results demonstrate successful calibration even with sparse overlap (<20%) and automatic detection of miscalibrations exceeding 0.8 meters. This work enables a "deploy-and-forget" sensor infrastructure that maintains calibration autonomously, reducing maintenance costs while improving tracking accuracy. The techniques generalize beyond transportation to any multi-sensor monitoring application requiring robust spatial alignment, including smart cities, industrial monitoring, and surveillance systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Universal corner symmetry and the orbit method for gravity

A universal symmetry algebra organizing the gravitational phase space has been recently found. It corresponds to the subset of diffeomorphisms that become physical at corners – codimension-2 surfaces supporting Noether charges. It applies to both finite distance and asymptotic corners. In this paper, we study this algebra and its representations, via the coadjoint orbit method. We show that generic orbits of the universal algebra split into sub-orbits spanned by finite distance and asymptotic corner symmetries, such that the full universal symmetry algebra gives rise to a unified treatment of corners in a manifold. We then identify the geometric structure that captures these algebraic properties on corners, which is the Atiyah Lie algebroid associated to a principal GL(2,R) $\ltimes$ R 2 -bundle. This structure is suggestive of the existence of a novel quantum gravitational theory which would unitarily glue such geometric structures, with spacetime geometries appearing as semi-classical configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Strength distributions of laminated FeNi-based metal amorphous nanocomposite ribbons

Metal Amorphous Nanocomposite (MANC) materials offer low losses at high magnetic switching frequency, enabling high power density motors with increased rotational speed. While MANCs have high strength, they are brittle. The use of motor components such as a rotor consisting of brittle material presents a reliability concern. Here, a promising MANC alloy is subjected to tensile tests and failure is observed with high-speed photography. A method is developed to prepare tensile specimens of laminated MANC and epoxy layers, simulating the stacking of an epoxy-impregnated tape-wound core. Tensile tests are conducted for single layer ribbon and for five- and ten-layer stacks of laminated material with thin layers of thermosetting epoxy. Failure distributions are shown to have increasing Weibull modulus with increasing layer count. The composite MANC material system is modeled using chain-of-bundles models. Using a k-failure model, we show that single ribbon strength distribution data can be used to predict well the failure distribution of laminated stacks. The agreement occurs when the assumed ineffective length, over which load is recovered in a failed layer, is comparable to the observed interlaminar separation length.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Reduced Order Models Generation for HTGRs Pebble Shuffling Procedure Optimization Studies

This report provides an initial study for producing reduced-order models (ROMs) of pebble-bed high temperature gas reactor (HTGR) models for the purposes of design optimization. As an initial study, this work is meant to be exploratory---identifying useful workflows and methods for ROM generation---and not meant to be a catch-all analysis of HTGR ROM generation and usage for optimization. This report summarizes three tasks performed in Fiscal Year 2022: 1) the creation of HTGR model, 2) the sensitivity analysis of model design parameters, and 3) an introduction to ROM generation techniques. The representative HTGR model created in this work is a multiphysics equilibrium-core using the BlueCRAB (comprehensive reactor analysis bundle) reactor analysis application, coupling four physical phenomena: neutronics, streamline depletion, porous flow thermal hydraulics, and pebble heat conduction. Part of the model creation was identifying some design parameters and quantities of interest that are relevant in an optimization analysis and adjustable in the model. The sensitivity analysis utilized a polynomial chaos expansion methodology to compute global sensitivity metrics. This analysis showed that thermal hydraulics parameters and quantities of interest had a relatively small impact on simulation results. Finally, the ROM generation work involved exploring three different ROM methodologies: polynomial regression, a Gaussian process, and artificial neural networks. Using a cross-validation technique to characterize ROM performance, the Gaussian process and single-layer artificial neural networks showed the most promising results. Overall, this study was insightful and the lessons learned will be invaluable for the eventual development of an HTGR design optimization workflow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Multiscale Mechanisms of Twisted Carbon Nanotube Yarns Probed In Situ by Soft X-rays during Tensile Loading

Piecing together carbon nanotubes (CNTs) into assemblies has so far failed to achieve the same elite strength performance metrics as individual CNTs, highlighting a critical deficiency in understanding the effects that the processing of individual nanostructures have on the performance of their derived macroscale assemblies, thereby hindering the development of a process-structure-performance map for these materials. Here, in this work, we propose a method to decouple the distribution orientation of nanoscale tortuosity and the microscale twist of CNT dry-spun yarns under applied loads via in situ soft X-ray probing at high energy (1200 eV) and low energy (280 eV), respectively. With this decoupling enabled by in situ soft X-ray scattering, we acquired a deeper understanding of the deformation mechanisms of these yarns. We found that for untreated yarns, the twist angle of collective CNT bundles at the macroscale is more sensitive to applied stress than the nanoscale alignment distribution. We also found that increasing nominal twist densities of yarns as well as increased strengthening via plasma treatments and polymer infiltration act to decrease the yarns’ sensitivity to realignment at the nanoscale and prevent failure by the slip mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Design of multi-scale protein complexes by hierarchical building block fusion

A systematic and robust approach to generating complex protein nanomaterials would have broad utility. We develop a hierarchical approach to designing multi-component protein assemblies from two classes of modular building blocks: designed helical repeat proteins (DHRs) and helical bundle oligomers (HBs). We first rigidly fuse DHRs to HBs to generate a large library of oligomeric building blocks. We then generate assemblies with cyclic, dihedral, and point group symmetries from these building blocks using architecture guided rigid helical fusion with new software named WORMS. X-ray crystallography and cryo-electron microscopy characterization show that the hierarchical design approach can accurately generate a wide range of assemblies, including a 43 nm diameter icosahedral nanocage. The computational methods and building block sets described here provide a very general route to de novo designed protein nanomaterials.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Overview of the Instrumentation for the Dark Energy Spectroscopic Instrument

Abstract The Dark Energy Spectroscopic Instrument (DESI) embarked on an ambitious 5 yr survey in 2021 May to explore the nature of dark energy with spectroscopic measurements of 40 million galaxies and quasars. DESI will determine precise redshifts and employ the baryon acoustic oscillation method to measure distances from the nearby universe to beyond redshift z > 3.5, and employ redshift space distortions to measure the growth of structure and probe potential modifications to general relativity. We describe the significant instrumentation we developed to conduct the DESI survey. This includes: a wide-field, 3.°2 diameter prime-focus corrector; a focal plane system with 5020 fiber positioners on the 0.812 m diameter, aspheric focal surface; 10 continuous, high-efficiency fiber cable bundles that connect the focal plane to the spectrographs; and 10 identical spectrographs. Each spectrograph employs a pair of dichroics to split the light into three channels that together record the light from 360–980 nm with a spectral resolution that ranges from 2000–5000. We describe the science requirements, their connection to the technical requirements, the management of the project, and interfaces between subsystems. DESI was installed at the 4 m Mayall Telescope at Kitt Peak National Observatory and has achieved all of its performance goals. Some performance highlights include an rms positioner accuracy of better than 0.″1 and a median signal-to-noise ratio of 7 of the [O ii ] doublet at 8 × 10 −17 erg s −1 cm −2 in 1000 s for galaxies at z = 1.4–1.6. We conclude with additional highlights from the on-sky validation and commissioning, key successes, and lessons learned.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

High yield production of ultrathin fibroid semiconducting nanowire of Ta 2 Pd 3 Se 8

Immediately after the demonstration of the high-quality electronic properties in various two dimensional (2D) van der Waals (vdW) crystals fabricated with mechanical exfoliation, many methods have been reported to explore and control large scale fabrications. Comparing with recent advancements in fabricating 2D atomic layered crystals, large scale production of one dimensional (1D) nanowires with thickness approaching molecular or atomic level still remains stagnant. Here, we demonstrate the high yield production of a 1D vdW material, semiconducting Ta 2 Pd 3 Se 8 nanowires, by means of liquid-phase exfoliation. The thinnest nanowire we have readily achieved is around 1 nm, corresponding to a bundle of one or two molecular ribbons. Transmission electron microscopy (TEM) and transport measurements reveal the as-fabricated Ta 2 Pd 3 Se 8 nanowires exhibit unexpected high crystallinity and chemical stability. Our low-frequency Raman spectroscopy reveals clear evidence of the existing of weak inter-ribbon bindings. The fabricated nanowire transistors exhibit high switching performance and promising applications for photodetectors.

1D semiconductor↗
DOE OSTI2021

Rubisco production in maize mesophyll cells through ectopic expression of subunits and chaperones

Abstract C4 plants, such as maize, strictly compartmentalize Rubisco to bundle sheath chloroplasts. The molecular basis for the restriction of Rubisco from the more abundant mesophyll chloroplasts is not fully understood. Mesophyll chloroplasts transcribe the Rubisco large subunit gene and, when normally quiescent transcription of the nuclear Rubisco small subunit gene family is overcome by ectopic expression, mesophyll chloroplasts still do not accumulate measurable Rubisco. Here we show that a combination of five ubiquitin promoter-driven nuclear transgenes expressed in maize leads to mesophyll accumulation of assembled Rubisco. These encode the Rubisco large and small subunits, Rubisco assembly factors 1 and 2, and the assembly factor Bundle sheath defective 2. In these plants, Rubisco large subunit accumulates in mesophyll cells, and appears to be assembled into a holoenzyme capable of binding the substrate analog CABP (carboxyarabinitol bisphosphate). Isotope discrimination assays suggest, however, that mesophyll Rubisco is not participating in carbon assimilation in these plants, most probably due to a lack of the substrate ribulose 1,5-bisphosphate and/or Rubisco activase. Overall, this work defines a minimal set of Rubisco assembly factors in planta and may help lead to methods of regulating the C4 pathway.

Plant Sciences↗
DOE OSTI2020

A Parallel Computing Infrastructure for Building Energy Simulation

In order to study grid-interactive efficient buildings, Pacific Northwest National Laboratories (PNNL) needs an infrastructure for urban-scale building energy modeling. Such an infrastructure should be fast, scalable, and easy-to-use. Given a set of data from the Energy Information Administration’s Commercial Building Energy Consumption Survey (CBECS) and tool to translate survey data into simulation inputs, this project aimed to conduct the simulation of the entire dataset in parallel. Before running the simulations, the necessary software was bundled into a container for use on the PNNL supercomputing network. Then, the parallel simulation workflow was designed using GNU Make, a file creation software, and submitted to a supercomputing partition which could run hundreds of simulations simultaneously. The EnergyPlus simulations output hourly electric meter data for each CBECS sample, which represents the electricity consumption of similar commercial buildings across the United States. Analyzing and visualizing the meter data is important to the future of the work, and this project wrote code to make common analysis methods simple, fast, and accessible. Moving forwards, the model will need to be expanded to include data from other sources and its accuracy will need to be improved and eventually validated.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗