Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical visualization”

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 361 records · Page 20

RE-EXAMINING SUNSPOT TILT ANGLE TO INCLUDE ANTI-HALE STATISTICS

Sunspot groups and bipolar magnetic regions (BMRs) serve as an observational diagnostic of the solar cycle. We use Debrecen Photohelographic Data (DPD) from 1974-2014 that determined sunspot tilt angles from daily white light observations, and data provided by Li and Ulrich that determined sunspot magnetic tilt angle using Mount Wilson magnetograms from 1974-2012. The magnetograms allowed for BMR tilt angles that were anti-Hale in configuration, so tilt values ranged from 0 to 360° rather than the more common ±90°. We explore the visual representation of magnetic tilt angles on a traditional butterfly diagram by plotting the mean area-weighted latitude of umbral activity in each bipolar sunspot group, including tilt information. The large scatter of tilt angles over the course of a single cycle and hemisphere prevents Joy's law from being visually identified in the tilt-butterfly diagram without further binning. The average latitude of anti-Hale regions does not differ from the average latitude of all regions in both hemispheres. The distribution of anti-Hale sunspot tilt angles are broadly distributed between 0 and 360° with a weak preference for east-west alignment 180° from their expected Joy's law angle. The anti-Hale sunspots display a log-normal size distribution similar to that of all sunspots, indicating no preferred size for anti-Hale sunspots. We report that 8.4% ± 0.8% of all bipolar sunspot regions are misclassified as Hale in traditional catalogs. This percentage is slightly higher for groups within 5° of the equator due to the misalignment of the magnetic and heliographic equators.

79 ASTRONOMY AND ASTROPHYSICS↗

Value of Nonrigid Registration of Pre-Procedure MR with Post-Procedure CT After Radiofrequency Ablation for Hepatocellular Carcinoma

PurposeTo evaluate the value of pre-radiofrequency ablation (RFA) MR and post-RFA CT registration for the assessment of the therapeutic response of hepatocellular carcinoma (HCC).Materials and MethodsA total of 178 patients with single HCC who received RFA as an initial treatment and had available pre-RFA MR and post-RFA CT images were included in this retrospective study. Two independent readers (one experienced radiologist, one inexperienced radiologist) scored the ablative margin (AM) of treated tumors on a four-point scale (1, residual tumor; 2, incomplete AM; 3, borderline AM; 4, sufficient AM), in two separate sessions: (1) visual comparison between pre-and post-RFA images; (2) with addition of nonrigid registration for pre- and post-RFA images. Local tumor progression (LTP) rates between low-risk (response score, 3–4) and high-risk groups (1–2) were analyzed using the Kaplan–Meier method at each interpretation session.ResultsThe patients’ reassignments after using the registered images were statistically significant for inexperienced reader (p < 0.001). In the inexperienced reader, LTP rates of low- and high-risk groups were significantly different with addition of registered images (session 2) (p < 0.001), but not significantly different in session 1 (p = 0.101). However, in the experienced reader, LTP rates of low- and high-risk groups were significantly different in both interpretation sessions (p < 0.001). Using the registered images, the cumulative incidence of LTP at 2 years was 3.0–6.6%, for the low-risk group, and 18.6–27.8% for the high-risk group.ConclusionRegistration between pre-RFA MR and post-RFA CT images may allow better assessment of the therapeutic response of HCC after RFA, especially for inexperienced radiologists, helping in the risk stratification for LTP.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Experimentally validated high-fidelity simulations of a liquid jet in supersonic crossflow

Here, utilizing recent advancements in computational schemes for compressible, multiphase flows, this work features a parametric study of a pure liquid jet in supersonic crossflow that involves simulating the atomization process for four values of momentum-flux ratio. These simulations are validated against experimental results measured with high-speed X-ray imaging, which confirm the accuracy of the numerical approach. Also, the effect of numerical resolution on some flow behavior is investigated, revealing convergence of the jet shape and surface instability wavelength. Analysis of the resulting sprays includes statistical descriptions of the liquid distribution, liquid structures created through breakup, interfacial instabilities, and dominant flow features. As the flowrate increases, the spray penetrates further, it becomes more disperse, and less liquid impacts the wall, but the droplet size distribution changes little. The wavelength of instabilities on the windward side of the jet diverges from measured trends in subsonic crossflows. In a visualization of the time-averaged flow, counter-rotating vortices are observed along the jet core and in the wake, affecting the process of primary atomization and early droplet trajectories.

42 ENGINEERING↗

FREDA: A Web Application for the Processing, Analysis, and Visualization of Fourier‐Transform Mass Spectrometry Data

The high-resolution measurement capability of Fourier-transform mass spectrometry (FT-MS) has made it a necessity for exploring the molecular composition of complex organic mixtures, like soil, plant, aquatic, and petroleum samples. This demand has driven a need for informatics tools to explore and analyze FT-MS data in a robust and reproducible manner. FREDA is an interactive web application developed to enable spectrometrists to format, process, and explore their FT-MS data without the need for statistical programming expertise. FREDA was built to explore outputs from a molecular identification tool, like CoreMS, and provide a suite of methods to filter data, compute chemical properties of peaks, statistically compare samples and groups of samples, conduct exploratory data analysis, and download the results with a report detailing all steps conducted. To demonstrate the utility of FREDA, an example analysis was conducted using FT-MS data from a soil microbiology study of samples collected in two different soil depths at the Sphagnum bog forest north of Grand Rapids, Minnesota. Differences between the two depths are observed using Kendrick, Gibbs free energy, and van Krevelen plots. G-tests are used to quantify a significant difference between the groups. All analyses and plotting are conducted using only the FREDA application. FREDA is an open-source and readily available web application that allows users to explore and make statistically valid conclusions about their FT-MS data. The application is available online (https://map.emsl.pnnl.gov/app/freda) with a tutorial web series (https://youtu.be/k5HLE2kNSBY?si=yB6sGoyvzxrFf5MP) and freely accessible code on Github (https://github.com/EMSL-Computing/FREDA).

47 OTHER INSTRUMENTATION↗

Classification of Cloud Particle Imagery from Aircraft Platforms Using Convolutional Neural Networks

Abstract A vast amount of ice crystal imagery exists from a variety of field campaign initiatives that can be utilized for cloud microphysical research. Here, nine convolutional neural networks are used to classify particles into nine regimes on over 10 million images from the Cloud Particle Imager probe, including liquid and frozen states and particles with evidence of riming. A transfer learning approach proves that the Visual Geometry Group (VGG-16) network best classifies imagery with respect to multiple performance metrics. Classification accuracies on a validation dataset reach 97% and surpass traditional automated classification. Furthermore, after initial model training and preprocessing, 10 000 images can be classified in approximately 35 s using 20 central processing unit cores and two graphics processing units, which reaches real-time classification capabilities. Statistical analysis of the classified images indicates that a large portion (57%) of the dataset is unusable, meaning the images are too blurry or represent indistinguishable small fragments. In addition, 19% of the dataset is classified as liquid drops. After removal of fragments, blurry images, and cloud drops, 38% of the remaining ice particles are largely intersecting the image border (≥10% cutoff) and therefore are considered unusable because of the inability to properly classify and dimensionalize. After this filtering, an unprecedented database of 1 560 364 images across all campaigns is available for parameter extraction and bulk statistics on specific particle types in a wide variety of storm systems, which can act to improve the current state of microphysical parameterizations.

54 ENVIRONMENTAL SCIENCES↗

Methyl-CpG island-associated genome signature tags

Disclosed is a method for analyzing the organismic complexity of a sample through analysis of the nucleic acid in the sample. In the disclosed method, through a series of steps, including digestion with a type II restriction enzyme, ligation of capture adapters and linkers and digestion with a type IIS restriction enzyme, genome signature tags are produced. The sequences of a statistically significant number of the signature tags are determined and the sequences are used to identify and quantify the organisms in the sample. Various embodiments of the invention described herein include methods for using single point genome signature tags to analyze the related families present in a sample, methods for analyzing sequences associated with hyper- and hypo-methylated CpG islands, methods for visualizing organismic complexity change in a sampling location over time and methods for generating the genome signature tag profile of a sample of fragmented DNA.

59 BASIC BIOLOGICAL SCIENCES↗

Band Excitation in Scanning Probe Microscopy: Recognition and Functional Imaging

Field confinement at the junction between a biased scanning probe microscope s (SPM) tip and solid surface enables local probing of various bias-induced transformations such as polarization switching, ionic motion, or electrochemical reactions to name a few. The nanoscale size of the biased region is smaller or comparable to features like grain boundaries and dislocations, potentially allows for the study of kinetics and thermodynamics at the level of a single defect. In contrast to classical statistically averaged approaches, this allows one to link structure to functionality and deterministically decipher associated mesoscopic and atomistic mechanisms. Furthermore, this type of information can serve as a fingerprint of local material functionality, allowing for local recognition imaging. Here, current progress in multidimensional SPM techniques based on band-excitation time and voltage spectroscopies is illustrated, including discussions on data acquisition, dimensionality reduction, and visualization along with future challenges and opportunities for the field.

Jesse, Stephen↗

Suite of Comprehensive Rapid Analysis Tools for Estimation (SOCRATES) - 19076

Long-term monitoring and remedy implementation require ongoing assessment and analytics, much of which relies on intensive and time-consuming data integration and engineering, data manipulation, and 'big-data' analytics. Many of these analyses, such as groundwater flow direction, contaminant rate attenuation, and non-parametric statistics, are spatial-temporal in nature and involve coalescing disparate data sets from a variety of sources with measurement frequencies from milliseconds to years, as well as total data volumes often exceeding the petabyte scale. Web-based analytical tools are advantageous for ongoing real-time or near-real-time ingestion of data from various sources for general data access and visualization, and also allow for complex, data rich calculations by leveraging expansive networks and virtually unlimited compute power. A web-based framework offers a natural architecture for such endeavors, as access to compute nodes through auto-scaling and dynamic frameworks enables speedy data manipulation and unlimited access to any data set that is made available via an http request. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Z-checker: A framework for assessing lossy compression of scientific data

Because of the vast volume of data being produced by today’s scientific simulations and experiments, lossy data compressor allowing user-controlled loss of accuracy during the compression is a relevant solution for significantly reducing the data size. However, lossy compressor developers and users are missing a tool to explore the features of scientific data sets and understand the data alteration after compression in a systematic and reliable way. To address this gap, we have designed and implemented a generic framework called Z-checker. On the one hand, Z-checker combines a battery of data analysis components for data compression. On the other hand, Z-checker is implemented as an open-source community tool to which users and developers can contribute and add new analysis components based on their additional analysis demands. In this article, we present a survey of existing lossy compressors. Then, we describe the design framework of Z-checker, in which we integrated evaluation metrics proposed in prior work as well as other analysis tools. Specifically, for lossy compressor developers, Z-checker can be used to characterize critical properties (such as entropy, distribution, power spectrum, principal component analysis, and autocorrelation) of any data set to improve compression strategies. For lossy compression users, Z-checker can detect the compression quality (compression ratio and bit rate) and provide various global distortion analysis comparing the original data with the decompressed data (peak signal-to-noise ratio, normalized mean squared error, rate–distortion, rate-compression error, spectral, distribution, and derivatives) and statistical analysis of the compression error (maximum, minimum, and average error; autocorrelation; and distribution of errors). Z-checker can perform the analysis with either coarse granularity (throughout the whole data set) or fine granularity (by user-defined blocks), such that the users and developers can select the best fit, adaptive compressors for different parts of the data set. Z-checker features a visualization interface displaying all analysis results in addition to some basic views of the data sets such as time series. To the best of our knowledge, Z-checker is the first tool designed to assess lossy compression comprehensively for scientific data sets.

97 MATHEMATICS AND COMPUTING↗

Federal Emergency Management Agency Mobile Housing Unit Running Gear Testing (Final Report)

The Federal Emergency Management Agency (FEMA) stages manufactured housing units (MHUs) at key locations to be used during disasters when other housing accommodations cannot be arranged for individuals or families. Due to the unpredictability of natural disasters and the limitations of the manufactured housing industry to produce large numbers of MHUs in very short periods of time, it is necessary for FEMA to store a large number of units for many months or years. With these long storage periods come degradation issues with frames, tires, bearing, and braking systems. These types of issues are also of interest to the U.S. Department of Transportation (DOT) Federal Motor Carrier Safety Administration (FMCSA) because regulated carriers normally transport these units and 49 CFR 393.75 Tires specifies tire loading restrictions for manufactured homes. To compare the running gear (tires, brakes, drums, bearings, axles, springs, and spring mounting hardware) currently specified in the 2014 Rugged Based Performance Requirements (RBPR) and the newly recommended running gear specified in the upcoming revision to the RBPR, the Oak Ridge National Laboratory (ORNL) outlined a Longevity Test Plan, which focused on the assessment of component longevity (failure rates by component type). The Test Plan also studied the forces imparted to the MHU from the running gear from interaction with the roadway for both the running gear in the 2014 RBPR and the running gear specified in the upcoming RBPR revision. Since the assessment of component longevity and failures are also issues of concern for FMCSA, both agencies collaborated to support the research presented in this report. The tests and data collection were conducted at the Navistar Proving Grounds (NPG) in New Carlisle (near South Bend), Indiana, during the May-November 2018 period. Four alternatives for the suspension system of a typical three-bedroom MHU were considered in this test. Alternative 0 was the MHU that was selected from the FEMA staging area in Selma, Alabama. Empirical data had shown that, on average, this alternative experiences one tire failure every 150 miles (or 0.0067 tire failures per mile traveled). Before testing, Alternative 0 was equipped with similar running gear hardware, although in new condition and installed following the manufacturer’s specifications. This became the Baseline Unit, or Alternative 1. Another unit, also selected from the FEMA staging area in Selma, Alabama, was equipped with better axles (i.e., axles that were rated for 7,000lbs as opposed to 6,000lbs for the Baseline Unit) and better tires. This became the Upgrade Unit, or Alternative 2. At the end of the Longevity Test, the Upgrade Unit was mounted with radial tires. This became Alternative 3.Alternative 1 (Baseline Unit) had, on average, one tire failure every 251 miles (or 0.0040 tire failures per mile traveled). So, by simply following the correct manufacturer’s specifications for the assembly of the transportation system, a reduction of 41% in the per mile tire failure was achieved. In the same test, Alternative 2 (Upgrade Unit) showed, on average, one tire failure every 728 miles traveled (or 0.0014 tire failure s per mile traveled). When the Upgrade Unit was mounted with radial tires in the last part of the test, no tire failures were observed for the 2,400 miles tested. That is, in the case of the radial tires (Alternative 3) the tire failure was less than one per 2,400 miles traveled (or less than 0.0004 tire failures per mile traveled). The data collected for the Longevity Test showed that the maximum accelerations were registered in all cases at the center of the MHU in the vertical direction. The overall maximum was registered by the Baseline Unit, Alternative 1. If that value is normalized to 100%, then the Upgrade Unit (Alternative 2) showed, on average, 81% at the same location, and Alternative 3 (radial-tire test) 59%. That is, the upgraded transportation system showed a reduction of 19% and 41% in the vertical accelerations transmitted to MHU at its center when regular and radial tires, respectively, were used. As part of the Test Protocol, NPG conducted a visual inspection of the MHUs and transportation system at the end of each test day. ORNL reviewed all the inspection reports in an attempt to determine if there was any difference in the damage that the Baseline and Upgrade MHUs showed during the test. The Upgrade Unit always showed lower accelerations transmitted to the MHU than the Baseline Unit. In most cases, these differences were statistically significant. The FEMA inspection report, conducted by a FEMA certified inspector at the end of the tests, included a shorter list of damaged items for the Baseline Unit than for the Upgrade Unit, although similar type of damages was described for both units. This may be an indication that, although lower than those of the Baseline Unit, the accelerations transmitted by the Upgrade Unit were above the threshold at which MHU elements start to fail.

42 ENGINEERING↗

Probing and Tuning Strain‐Localized Exciton Emission in 2D Material Bubbles at Room Temperature

In monolayer transition metal dichalcogenides bubbles-nanoscale deformations typically exhibiting a dome-like shape-Excitons are confined by the strain effect, which exhibits extraordinary emission properties, such as single photon generation, enhanced light emission, and spectrally tunable excitonic states. While the strain profiles of these bubbles are extensively studied, this work provides an approach 1) to directly visualize the associated exciton properties in bubbles formed in WSe2 monolayer, revealing an intrinsic emission wavelength shift of ≈40 nm, and 2) actively modify local strain, enabling further exciton emission tuning over a range of 50 nm. These are achieved by emission mapping and nanoindentation using a dielectric near-field probe, which enables the detection of local emission spectra and emission lifetimes within individual bubbles. Statistical analysis of 67 bubbles uncovers an emission wavelength distribution centered around 780 nm. Furthermore, saturation behavior in the power-dependent studies and the associated lifetime change reveal the localized nature of the strain-induced states. These findings provide direct insights into the strain-localized emission dynamics in bubbles and establish a robust framework for non-destructive, reversible, and predictable nanoscale emission control, presenting a potential avenue for developing next-generation tunable quantum optical sources.

2D materials bubbles↗

External Threat Risk Assessment Algorithm (ExTRAA)

Two risk assessment algorithms and philosophies have been augmented and combined to form a new algorit hm, the External Threat Risk Assessment Algorithm (ExTRAA), that allows for effective and statistically sound analysis of external threat sources in relation to individual attack methods . In addition to the attack method use probability and the attack method employment consequence, t he concept of defining threat sources is added to the risk assessment process. Sample data is tabulated and depicted in radar plots and bar graphs for algorithm demonstration purposes. The largest success of ExTRAA is its ability to visualize the kind of r isk posed in a given situation using the radar plot method.

97 MATHEMATICS AND COMPUTING↗

SU-C-18C-03: Dual-Energy X-Ray Fluoroscopy Imaging System

Purpose: This work studies the clinical utility of dual energy (DE) subtraction fluoroscopy for fiducial-free tumor tracking in lung radiation therapy (RT). Improvement in tumor visualization and quantification of tumor shift within a breathing cycle were analyzed. Methods: Twenty subjects who were undergoing RT for lung cancer were recruited following institutional review board approval. The subjects had a range of tumor sizes, locations in the lungs, and body sizes. An x-ray imaging system was setup with the following components: (a) x-ray tube (Varian G-242, Varian Medical Systems (VMS), CA) (b) flat panel detector (4030CB, VMS, CA) and (c) x-ray generator (EPS 50RF, EMD, Canada). Firmware and software modifications were made to the generator to allow 10 x-ray pulse pairs with alternating low/high kV, 100 ms apart for ∼4s (one breathing cycle). Images were obtained at 4 angles: 0°, 45°, 90° and 135°. Weighted subtraction of a kV-pair image set was used to create a “bone-free” image of the lungs. The 2D tumor-shift in each subtracted image and the 3D shift during a breathing cycle was calculated using all views. Results: The subjects enrolled had the following statistics: average age 62.3±7.1 years, 5 female/15 male, 11 had tumors on the right and 9 on the left and the average tumor size was ∼31.4±10.8 mm. X-ray imaging conditions for the pulse pairs were: 70/120 kVp, 280/221 mA and 65/8 ms. For views where these parameters were insufficient 80/130 kVp, 280/221 mA and 60/12 ms was used. Tumor visibility improved for 0°, 45°, 90° and 135° in 100%, 55%, 75% and 80% of the cases respectively. Tumor shift during a breathing cycle was: 2.4±1.0 mm AP, 2.7±1.4 mm LR and 7.6±4.8 mm IS. Conclusion: DE subtraction fluoroscopy allowed improved visualization and quantification of movement of tumors in the lungs during a breathing cycle. This study was entirely funded by Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗

Realistic simulation of reduced-dose CT with noise modeling and sinogram synthesis using DICOM CT images

Purpose: Reducing the patient dose while maintaining the diagnostic image quality during CT exams is the subject of a growing number of studies, in which simulations of reduced-dose CT with patient data have been used as an effective technique when exploring the potential of various dose reduction techniques. Difficulties in accessing raw sinogram data, however, have restricted the use of this technique to a limited number of institutions. Here, we present a novel reduced-dose CT simulation technique which provides realistic low-dose images without the requirement of raw sinogram data. Methods: Two key characteristics of CT systems, the noise equivalent quanta (NEQ) and the algorithmic modulation transfer function (MTF), were measured for various combinations of object attenuation and tube currents by analyzing the noise power spectrum (NPS) of CT images obtained with a set of phantoms. Those measurements were used to develop a comprehensive CT noise model covering the reduced x-ray photon flux, object attenuation, system noise, and bow-tie filter, which was then employed to generate a simulated noise sinogram for the reduced-dose condition with the use of a synthetic sinogram generated from a reference CT image. The simulated noise sinogram was filtered with the algorithmic MTF and back-projected to create a noise CT image, which was then added to the reference CT image, finally providing a simulated reduced-dose CT image. The simulation performance was evaluated in terms of the degree of NPS similarity, the noise magnitude, the bow-tie filter effect, and the streak noise pattern at photon starvation sites with the set of phantom images. Results: The simulation results showed good agreement with actual low-dose CT images in terms of their visual appearance and in a quantitative evaluation test. The magnitude and shape of the NPS curves of the simulated low-dose images agreed well with those of real low-dose images, showing discrepancies of less than +/−3.2% in terms of the noise power at the peak height and +/−1.2% in terms of the spatial frequency at the peak height. The magnitudes of the noise measured for 12 different combinations the phantom size, tube current, and reconstruction kernel for the simulated and real low-dose images were very similar, with differences of 0.1 to 4.7%. Thep value for a statistical testing of the difference in the noise magnitude ranged from 0.99 to 0.11, showing that there was no difference statistically between the noise magnitudes of the real and simulated low-dose images using our method. The strength and pattern of the streak noise in an anthropomorphic phantom was also consistent with expectations. Conclusions: A novel reduced-dose CT simulation technique was developed which uses only CT images while not requiring raw sinogram data. Our method can provide realistic simulation results under reduced-dose conditions both in terms of the noise magnitude and the textual appearance. This technique has the potential to promote clinical research for patient dose reductions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Towards the distributed burning regime in turbulent premixed flames

Three-dimensional numerical simulations of canonical statistically steady, statistically planar turbulent flames have been used in an attempt to produce distributed burning in lean methane and hydrogen flames. Dilatation across the flame means that extremely large Karlovitz numbers are required; even at the extreme levels of turbulence studied (up to a Karlovitz number of 8767) distributed burning was only achieved in the hydrogen case. In this case, turbulence was found to broaden the reaction zone visually by around an order of magnitude, and thermodiffusive effects (typically present for lean hydrogen flames) were not observed. In the preheat zone, the species compositions differ considerably from those of one-dimensional flames based a number of different transport models (mixture averaged, unity Lewis number and a turbulent eddy viscosity model). The behaviour is a characteristic of turbulence dominating non-unity Lewis number species transport, and the distinct limit is again attributed to dilatation and its effect on the turbulence. Peak local reaction rates are found to be lower in the distributed case than in the lower Karlovitz cases but higher than in the laminar flame, which is attributed to effects that arise from the modified fuel-temperature distribution that results from turbulent mixing dominating low Lewis number thermodiffusive effects. Finally, approaches to achieve distributed burning at realisable conditions are discussed; factors that increase the likelihood of realising distributed burning are higher pressure, lower equivalence ratio, higher Lewis number and lower reactant temperature.

42 ENGINEERING↗

Edge counts for the auxiliary pair graph within the graphical unitary group approach

Closed-form expressions are presented for the numbers of edges in the auxiliary pair graphs (APGs) associated with non spin-orbit and spin-orbit Shavitt graphs for full configuration interaction expansions. A Shavitt graph is a visual representation of a configuration state function expansion space constructed via the graphical unitary group approach (GUGA). An APG is an organisational aid and a programmatic tool generated from a Shavitt graph. The number of edges in an APG determines bounds on the computational scaling as a function of the total numbers of electrons, orbitals, and spin multiplicities. The edge counts extend a suite of Shavitt graph statistics based on these functional parameters. The derivation and the presentation of the formulas for the edge counts has been assisted by the bra-ket interchange symmetry and the particle-hole interchange symmetry in the GUGA formalism. Furthermore, these symmetry operators produce one-to-one correspondences between various sets of edges, and this yields identities among some edge count formulas. There are 208 possible edge types. Of these, some do not contribute to two-electron operators, some are related by bra-ket interchange symmetry, and some are related by particle-hole interchange symmetry. For the remaining unique edge types, explicit expressions are derived for the numbers of edges.

74 ATOMIC AND MOLECULAR PHYSICS↗

City Scale Modeling with OpenStudio

Assessing the impact of energy efficiency technologies at a district or city scale is of great interest to local governments, real estate developers, utility companies, and policymakers. This paper describes a flexible framework that can be used to create and run district and city scale building energy simulations. The framework is built around the new OpenStudio City Database (CityDB). Building footprints, building height, building type, and other data can be imported from public records or other sources. Missing data can be inferred or assigned from a statistical sampling of other datasets. Once all required data is available, OpenStudio Measures are used to create starting point energy models and to model energy efficiency measures for each building. Together this framework allows a user to pose several scenarios such as 'what if 30% of the commercial retail buildings added rooftop solar' or 'what if all elementary schools converted to ground source heat pumps' and then visualize the impacts at a district or city scale. This paper focuses on modeling existing building stock using public records. However, the framework is capable of supporting the evaluation of new construction, district systems, and the use of proprietary data sources.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗