Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “I&C Modernization”

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 595 records · Page 33

Artificial Intelligence for Natural Gas Utilities: A Primer

Modern natural gas utilities face numerous challenges and competing priorities from various stakeholders. Policymakers, customers, and advocacy groups want to see gas utilities improve performance on safety, reliability, resilience, affordability, and environmental stewardship. State utility regulators — public utility commissions — are responsible for overseeing utility performance, ensuring that ratepayer funds are being spent in the public interest, and aligning utility goals with public goals. The use of new technologies is critical to enabling cost-effective performance on these attributes. Artificial intelligence (AI) is a widely used term among utilities and regulators, but the term means different things to different stakeholders, and it is often used to describe data analytics approaches that fall short of the formal definition of AI, which is: “…the ability of a machine to receive inputs and produce a behavior or reaction similar to that of an intelligent human being.” AI (and related tools, techniques, and technologies) can help utilities solve current and emerging challenges. By combining customer and system data with analytical tools and technologies, AI can augment human decision-makers by assisting in identifying problems and events before they occur, enabling resources to be more efficiently directed across utility infrastructure. The intended primary audience for this primer is state regulators, although utilities and other stakeholders might also find it useful and relevant to improve their awareness of AI. The objectives of this primer are to: (a) offer a set of broadly applicable definitions for AI and related terms, allowing regulators, utilities, and other stakeholders to speak the same language; (b) discuss how AI is currently being implemented in the gas utility sector; and (c) understand the challenges affecting AI solutions and how tools might be implemented in the future. This primer fits within NARUC’s goals of providing impartial information to improve the ability of public utility commissions to regulate in the public interest. As such, this primer does not seek to recommend AI over any other investment, nor does it endorse any particular vendor, product, or approach. It does seek to prepare state regulators to oversee AI investments by sharing information about the current landscape of commercially available tools. To these ends, the primer is organized as follows: Section I discusses the current environment in which natural gas utilities operate and how AI, when thoughtfully designed and implemented, can enable utilities to achieve performance goals; Section II offers definitions of AI and related terms within the data analytics discipline; Section III provides three current opportunities for which AI can offer solutions: replacing aging gas distribution infrastructure, preventing excavator damage to gas distribution infrastructure, and improving energy efficiency programs. This section discusses each problem statement in detail. Second, Section III includes a discussion of how costs and benefits of investments to solve each problem are measured. And third, this section offers real-world examples of utility implementation of AI solutions; Section IV discusses challenges with implementing AI, both from the perspective of utilities and regulators; Section V suggests areas in which AI could feasibly be implemented in the near future; Finally, Section VI offers concluding thoughts and areas for further research.

03 NATURAL GAS↗

Impacts of sulfonic acids on fungal manganese oxide production

Microbial manganese (Mn) oxidation plays a critical role in Mn(III/IV) oxide formation in modern day environmental systems. These oxidation processes and resulting biominerals are sensitive to aqueous conditions, such as pH and dissolved constituent concentrations. With Mn and sulfur (S) biogeochemical cycling closely linked in many environmental systems, and dissolved organic sulfur comprising a substantial pool of total sulfur in several environments, the impact of dissolved organic sulfur compounds on Mn redox processes is important to consider. Sulfonic acids, environmentally ubiquitous organosulfur compounds, play substantial roles in S cycling in many natural and contaminated systems. Research to assess the effects of these abundant sulfonic acids on Mn biogeochemical cycling, microbial Mn oxidation processes, and Mn biominerals is needed for understanding and predicting the impact of coupled S and Mn biogeochemical cycles, particularly in environments with dynamic redox gradients or in anthropogenically contaminated systems. Further, with research on microbial and abiotic Mn oxidation processes often using aminosulfonic acids to control pH, understanding the impact of such sulfonic acids on microbial Mn oxidation processes is critical. Several recent studies found that commonly used zwitterionic N substituted aminosulfonic acids, known as Good’s buffers, such as HEPES and MES, can alter abiotic birnessite sheet structures. Here we investigate the impact of two sulfonic acids with broad applications to natural and contaminated sites as well as laboratory settings (HEPES and MES) on fungal Mn oxidation relative to a carbonate buffer and a buffer-free control by three Ascomycete fungi known to oxidize Mn(II): Stagonospora sp. SRC1lsM3a, Paraphaeosphaeria sporulosa AP3s5–JAC2a, and Plectosphaerella cucumerina DS2psM2a2. Structural analyses of the products show that sulfonic acids promote Mn oxidation by P. cucumerina, producing hexagonally symmetric phyllomanganates analogous to hexagonal birnessite or c–disordered H + birnessite [(Ca,Na,K)(Mn 4+ ,Mn 3+ ,$\square$)O 2 nH 2 O], with solid–associated Mn(II) bound to vacancy sites and biomass, while in their absence almost all Mn remains as either aqueous Mn(II) or solid–associated Mn(II) bound to biomass. In contrast, sulfonic acids exert the opposite effect on Mn oxidation by P. sporulosa, with their presence suppressing Mn(II) oxidation to Mn(IV), likely leading to the formation of mycogenic bixbyite (Mn 3+ 2 O 3 ) while the buffer–free control forms a poorly crystalline phyllo- or tectomanganate. Meanwhile, all treatments exert a minimal effect on Mn(II) uptake from solution and Mn oxidation with Stagonospora sp., with all experimental systems and controls forming poorly crystalline, hexagonally symmetric phyllomanganates. The fact that the sulfonic acids here studied exert similar effects on Mn oxidation, but substantially different effects for each fungus, suggests they affect Mn oxidation via mechanistically similar pathways that are likely dependent on interactions with fungal exudates (which vary from species to species) or specific fungal Mn oxidation processes. Interestingly, for all fungi, MES increases Mn(III) in the resulting biominerals, while the carbonate buffer consistently decreases Mn(III). Further, these results clearly demonstrate that sulfonic acids not only alter Mn oxide structures, as has been previously noted in abiotic studies, but can interfere with Mn oxidation reactions themselves, highlighting the incredible sensitivity of both Mn oxide structures and the Mn oxidation process to the aqueous environment.

58 GEOSCIENCES↗

A hybrid optoelectronic Mott insulator

The coupling of electronic degrees of freedom in materials to create “hybridized functionalities” is a holy grail of modern condensed matter physics that may produce versatile mechanisms of control. Correlated electron systems often exhibit coupled degrees of freedom with a high degree of tunability which sometimes lead to hybridized functionalities based on external stimuli. However, the mechanisms of tunability and the sensitivity to external stimuli are determined by intrinsic material properties which are not always controllable. A Mott metal-insulator transition (MIT) is technologically attractive due to the large changes in resistance, tunable by doping, strain, electric fields, and orbital occupancy but not, in and of itself, controllable with light. Here, an alternate approach is presented to produce optical functionalities using a properly engineered photoconductor/strongly correlated hybrid heterostructure. This approach combines a photoconductor, which does not exhibit an MIT, with a strongly correlated oxide, which is not photoconducting. Due to the intimate proximity between the two materials, the heterostructure exhibits giant volatile and nonvolatile, photoinduced resistivity changes with substantial shifts in the MIT transition temperatures. This approach can be extended to other judicious combinations of strongly correlated materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Towards Automatic and Agile AI/ML Accelerator Design with End-to-End Synthesis

Domain-specific designs offer greater energy efficiency and performance gain than general-purpose processors. For this reason, modern system-on-chips have a significant portion of their silicon area with custom accelerators. However, designing hardware by hand is laborious and time-consuming, given the large design space and the performance, power, and area constraints that are not realized in the software. Moreover, domain-specific algorithms (e.g., machine learning models) are evolving quickly, challenging the accelerator design further. To address these issues, this paper presents SODA Synthesizer, an automated open-source high-level ML framework to Verilog modular compiler targeting AI/ML Application-Specific Integrated Circuits (ASICs) accelerators. SODA tightly couples the Multi- Level Intermediate Representation (MLIR) compiler infrastructure [24] and open-source HLS approaches. Thus, SODA can support various ML frameworks and algorithms and can perform optimizations that combine specialized architecture templates and conventional HLS to generate the hardware modules. In addition, SODA’s closed-loop design space exploration (DSE) engine allows developers to perform end-to-end design space explorations on different metrics and technology nodes.

Zhang, Jeff↗

Tausch: A halo exchange library for large heterogeneous computing systems using MPI, OpenCL, and CUDA

Exchanging halo data is a common task in modern scientific computing applications and efficient handling of this operation is critical for the performance of the overall simulation. Tausch is a novel header-only library that provides a simple API for efficiently handling these types of data movements. Tausch supports both simple CPU-only systems, but also more complex heterogeneous systems with both CPUs and GPUs. It currently supports both OpenCL and CUDA for communicating with GPGPU devices, and allows for communication between GPGPUs and CPUs. The API allows for drop-in replacement in existing codes and can be used for the communication layer in new codes. This paper provides an overview of the approach taken in Tausch, and a performance analysis that demonstrates expected and achieved performance. Here, we highlight the ease of use and performance with three applications: First Tausch is compared to the halo exchange framework from two Mantevo applications, HPCCG and miniFE, and then it is used to replace a legacy halo exchange library in the flexible multigrid solver framework Cedar.

97 MATHEMATICS AND COMPUTING↗

Assessing tension metrics with dark energy survey and Planck data

ABSTRACT Quantifying tensions – inconsistencies amongst measurements of cosmological parameters by different experiments – has emerged as a crucial part of modern cosmological data analysis. Statistically significant tensions between two experiments or cosmological probes may indicate new physics extending beyond the standard cosmological model and need to be promptly identified. We apply several tension estimators proposed in the literature to the dark energy survey (DES) large-scale structure measurement and Planck cosmic microwave background data. We first evaluate the responsiveness of these metrics to an input tension artificially introduced between the two, using synthetic DES data. We then apply the metrics to the comparison of Planck and actual DES Year 1 data. We find that the parameter differences, Eigentension, and Suspiciousness metrics all yield similar results on both simulated and real data, while the Bayes ratio is inconsistent with the rest due to its dependence on the prior volume. Using these metrics, we calculate the tension between DES Year 1 3 × 2pt and Planck, finding the surveys to be in ∼2.3σ tension under the ΛCDM paradigm. This suite of metrics provides a toolset for robustly testing tensions in the DES Year 3 data and beyond.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Dose Summation Strategies for External Beam Radiation Therapy and Brachytherapy in Gynecologic Malignancy: A Review from the NRG Oncology and NCTN Medical Physics Subcommittees

Definitive, nonsurgical management of gynecologic malignancies involves external beam radiation therapy (EBRT) and/or brachytherapy (BT). Summation of the cumulative dose is critical to assess the total biologic effective dose to targets and organs at risk. Cumulative dose calculation from EBRT and BT can be performed with or without image registration (IR) and biologic dose summation. Among these dose summation strategies, linear addition of dose-volume histogram (DVH) parameters without IR is the global standard for composite dose reporting. This approach stems from an era without image guidance and simple external beam and brachytherapy treatment approaches. With technological advances, EBRT and high-dose-rate BT have evolved to allow for volume-based treatment planning and delivery. Modern conformal therapeutic radiation involves volumetric or intensity modulated EBRT, capable of simultaneously treating multiple targets at different specified dose levels. Therefore, given the complexity of modern radiation treatment, the linear addition of DVH parameters from EBRT and high-dose-rate BT is challenging to represent the combined dose distribution. Deformable image registration (DIR) between EBRT and image guided brachytherapy (IGBT) data sets may provide a more nuanced calculation of multimodal dose accumulation. However, DIR is still nascent in this regard, and needs further development for accuracy and efficiency for clinical use. Biologic dose summation can combine physical dose maps from EBRT and each IGBT fraction, thereby generating a composite DVH from the biologic effective dose. However, accurate radiobiologic parameters are tissue-dependent and not well characterized. A combination of voxel-based DIR and biologic weighted dose maps may be the best approximation of dose accumulation but remains invalidated. The purpose of this report is to review dose summation strategies for EBRT and BT, including conventional equivalent dose in 2-Gy fractions dose summation without image registration, physical dose summation using 3-dimensional rigid IR and DIR, and biologic dose summation. We also provide general clinical workflows for IGBT with a focus on cervical cancer.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Star formation and AGN feedback in the local Universe: Combining LOFAR and MaNGA

The effect of active galactic nuclei (AGN) on their host galaxies – in particular their levels of star formation – remains one of the key outstanding questions of galaxy evolution. Successful cosmological models of galaxy evolution require a fraction of energy released by an AGN to be redistributed into the interstellar medium to reproduce the observed stellar mass and luminosity function and to prevent the formation of over-massive galaxies. Observations have confirmed that the radio-AGN population is energetically capable of heating and redistributing gas at all phases, however, direct evidence of AGN enhancing or quenching star formation remains rare. With modern, deep radio surveys and large integral field spectroscopy (IFS) surveys, we can detect fainter synchrotron emission from AGN jets and accurately probe the star-forming properties of galaxies, respectively. Here, in this paper, we combine data from the LOw Frequency ARray Two-meter Sky Survey (LoTSS) with data from one of the largest optical IFS surveys, Mapping Nearby Galaxies at Apache Point Observatory (MaNGA), to probe the star-forming properties of 307 local (z < 0.15) galaxies that host radio-detected AGN (RDAGN). We compare our results to a robust control sample of non-active galaxies that each match the stellar mass, redshift, visual morphology, and inclination of a RDAGN host. We find that RDAGN and control galaxies have broad star-formation rate (SFR) distributions, typically lie below the SFMS, and have negative stellar light-weighted age gradients. These results indicate that AGN selected based on their current activity are not responsible for suppressing their host galaxies’ star formation. Rather, our results support the maintenance mode role that radio AGN are expected to have in the local Universe.

79 ASTRONOMY AND ASTROPHYSICS↗

First identification of resonant states in the {sup 45}V(p, {gamma}){sup 46}Cr reaction and their influence on the production of {sup 44}Ti in Core Collapse Supernovae

The observation of active 44Ti in supernova remnants offers the potential to solve one of the most debated questions in modern astrophysics, the exact underlying explosion mechanism of core collapse supernovae (CCSNe). In particular, a comparison between the predicted synthesized yield of 44Ti and the ejected flux is expected to allow for a determination of the mass cut of the star. Unfortunately, such comparisons are currently severely hindered by large uncertainties in the nuclear reactions governing the production and destruction of 44Ti. On this note, the astrophysical 45Vðp; γÞ46Cr reaction is expected to play a decisive role. We present the first experimental information on resonant states in the 45Vðp; γÞ46Cr reaction. Excitation energies have been measured with high precision for ten previously unknown, low-spin, proton-unbound states in 46Cr, and spin-parity-assignment constraints have been made using observed γ-decay patterns and shell-model calculations. We find that an l ¼ 1 resonance at 714 keV dominates the 45Vðp; γÞ46Cr reaction for the entire temperature range of α-rich freezeout. Moreover, the nominal rate is found to be ∼25–50% of previous estimates for T > 1.2 GK, indicating an increase in 44Ti production in CCSNe of ∼20–35%.

Cousins, C↗

First observation of 28 O

Subjecting a physical system to extreme conditions is one of the means often used to obtain a better understanding and deeper insight into its organization and structure. In the case of the atomic nucleus, one such approach is to investigate isotopes that have very different neutron-to-proton ($N/Z$) ratios than in stable nuclei. Light, neutron-rich isotopes exhibit the most asymmetric $N/Z$ ratios and those lying beyond the limits of binding, which undergo spontaneous neutron emission and exist only as very short-lived resonances (about 10 -21 s), provide the most stringent tests of modern nuclear-structure theories. Here we report on the first observation of 28 O and 27 O through their decay into 24 O and four and three neutrons, respectively. The 28 O nucleus is of particular interest as, with the $Z$ = 8 and $N$ = 20 magic numbers, it is expected in the standard shell-model picture of nuclear structure to be one of a relatively small number of so-called ‘doubly magic’ nuclei. Furthermore, both 27 O and 28 O were found to exist as narrow, low-lying resonances and their decay energies are compared here to the results of sophisticated theoretical modelling, including a large-scale shell-model calculation and a newly developed statistical approach. In both cases, the underlying nuclear interactions were derived from effective field theories of quantum chromodynamics. Finally, it is shown that the cross-section for the production of 28 O from a 29 F beam is consistent with it not exhibiting a closed $N$ = 20 shell structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Using Systems Theoretic Process Analysis and Causal Analysis to Map and Manage Organizational Information to Enable Digitalization and Information Automation

The overarching goal of this Light Water Reactor Sustainability Program–supported research and development project is to provide planning tools and comprehensive guidance to utilities considering or undertaking full nuclear plant modernization. The results of this research will provide the nuclear industry with a comprehensive and usable solution, including guidance, lessons learned, methods, and planning tools. This research is currently working to provide guidance on digitalization and information automation to enable the evolution of data to information, insight, and action—thereby allowing utilities to operate safely and cost-competitively with all other electrical generation sources. Light Water Reactor Sustainability Program researchers have also recently started investigating how human and technology integration principles, information automation, and digitalization enable data evolution. These researchers are currently in the process of validating the use of System-Theoretic Process Analysis to define high-level safety constraints in the United States Nuclear Regulatory Commission’s problem identification and resolution process (i.e., a plant compliance information gathering activity). The next step in this research, which is described in the following sections of this report, is to map out data evolution in a use case to identify inefficiencies in another aspect of plant compliance information gathering and communication activities—event investigations and root cause analyses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Downscaling atmospheric chemistry simulations with physically consistent deep learning

Abstract. Recent advances in deep convolutional neural network (CNN)-based super resolution can be used to downscale atmospheric chemistry simulations with substantially higher accuracy than conventional downscaling methods. This work both demonstrates the downscaling capabilities of modern CNN-based single image super resolution and video super-resolution schemes and develops modifications to these schemes to ensure they are appropriate for use with physical science data. The CNN-based video super-resolution schemes in particular incur only 39 % to 54 % of the grid-cell-level error of interpolation schemes and generate outputs with extremely realistic small-scale variability based on multiple perceptual quality metrics while performing a large (8×10) increase in resolution in the spatial dimensions. Methods are introduced to strictly enforce physical conservation laws within CNNs, perform large and asymmetric resolution changes between common model grid resolutions, account for non-uniform grid-cell areas, super-resolve lognormally distributed datasets, and leverage additional inputs such as high-resolution climatologies and model state variables. High-resolution chemistry simulations are critical for modeling regional air quality and for understanding future climate, and CNN-based downscaling has the potential to generate these high-resolution simulations and ensembles at a fraction of the computational cost.

54 ENVIRONMENTAL SCIENCES↗

Achievements and Prospects of Global Broadband Seismographic Networks After 30 Years of Continuous Geophysical Observations

Global seismographic networks (GSNs) emerged during the late nineteenth and early twentieth centuries, facilitated by seminal international developments in theory, technology, instrumentation, and data exchange. The mid- to late-twentieth century saw the creation of the World-Wide Standardized Seismographic Network (1961) and International Deployment of Accelerometers (1976), which advanced global geographic coverage as seismometer bandwidth increased greatly allowing for the recording of the Earth's principal seismic spectrum. The modern era of global observations and rapid data access began during the 1980s, and notably included the inception of the GEOSCOPE initiative (1982) and GSN (1988). Through continual improvements, GEOSCOPE and the GSN have realized near-real time recording of ground motion with state-of-art data quality, dynamic range, and timing precision to encompass 180 seismic stations, many in very remote locations. Data from GSNs are increasingly integrated with other geophysical data (e.g., space geodesy, infrasound and Interferometric Synthetic Aperture Radar). Globally distributed seismic data are critical to resolving crust, mantle, and core structure; illuminating features of the plate tectonic and mantle convection system; rapid characterization of earthquakes; identification of potential tsunamis; global nuclear test verification; and provide sensitive proxies for environmental changes. As the global geosciences community continues to advance our understanding of Earth structure and processes controlling elastic wave propagation, GSN infrastructure offers a springboard to realize increasingly multi-instrument geophysical observatories. Here, we review the historical, scientific, and monitoring heritage of GSNs, summarize key discoveries, and discuss future associated opportunities for Earth Science.

58 GEOSCIENCES↗

Advected carbon younger than the sediment fuels microbial metabolism in a pumped deep aquifer

Ever deeper wells are drilled worldwide to pump potable groundwater. Recent studies argue that overpumping compresses clays and releases reactive dissolved organic carbon (DOC), which in turn drives a series of microbial reactions that affect groundwater potability including arsenic concentrations. Here, we use a novel method to measure the radiocarbon ages of microbial RNA to determine the source of reactive or metabolizable DOC and argue against the compression of clays as the sole source of carbon. We show that microbial RNA (5,230; 5,550; 6,250 yr; n = 3 wells), DOC (280-10,800 yr; n = 13), and methane (modern to 6,240 yr; n = 3), from an overpumped deep aquifer in Bangladesh are much younger than the overlying clay layers deposited during the Pleistocene over 12,000 years ago. Mass-balance indicates that at least half of the carbon incorporated into RNA has to come from reactive DOC or methane that is advected downward via vertical recharge. This metabolism of advected organic carbon could have implications for the quality of water pumped from deep aquifers.

Biological and medical sciences↗

Strictly Enforcing Invertibility and Conservation in CNN-Based Super Resolution for Scientific Datasets

Abstract Recently, deep convolutional neural networks (CNNs) have revolutionized image “super resolution” (SR), dramatically outperforming past methods for enhancing image resolution. They could be a boon for the many scientific fields that involve imaging or any regularly gridded datasets: satellite remote sensing, radar meteorology, medical imaging, numerical modeling, and so on. Unfortunately, while SR-CNNs produce visually compelling results, they do not necessarily conserve physical quantities between their low-resolution inputs and high-resolution outputs when applied to scientific datasets. Here, a method for “downsampling enforcement” in SR-CNNs is proposed. A differentiable operator is derived that, when applied as the final transfer function of a CNN, ensures the high-resolution outputs exactly reproduce the low-resolution inputs under 2D-average downsampling while improving performance of the SR schemes. The method is demonstrated across seven modern CNN-based SR schemes on several benchmark image datasets, and applications to weather radar, satellite imager, and climate model data are shown. The approach improves training time and performance while ensuring physical consistency between the super-resolved and low-resolution data. Significance Statement Recent advancements in using deep learning to increase the resolution of images have substantial potential across the many scientific fields that use images and image-like data. Most image super-resolution research has focused on the visual quality of outputs, however, and is not necessarily well suited for use with scientific data where known physics constraints may need to be enforced. Here, we introduce a method to modify existing deep neural network architectures so that they strictly conserve physical quantities in the input field when “super resolving” scientific data and find that the method can improve performance across a wide range of datasets and neural networks. Integration of known physics and adherence to established physical constraints into deep neural networks will be a critical step before their potential can be fully realized in the physical sciences.

54 ENVIRONMENTAL SCIENCES↗

A Functional All-Hazard Approach to Critical Infrastructure Dependency Analysis

The critical infrastructures protection landscape is a vast and varied pattern of independent, but interconnected infrastructure systems that are essential to the function of our modern society. The U.S. policy on critical infrastructure protection has been continually evolving since the “President’s Commission on Critical Infrastructure Protection” was published in 1997. In response to these policies, federal, state, and local governments, along with research institutions, have invested a substantial amount of time and effort into identifying and analyzing critical infrastructure, their functions, and dependencies/interdependencies to better understand their vulnerabilities. To date, the ability to assess vulnerabilities, resiliency, and priorities for protecting interdependent critical infrastructure systems from an all-hazards perspective remains a difficult problem. In this paper we introduce the All-Hazards Analysis (AHA) methodology, which provides an integrated functional basis across infrastructure systems, through the implementation of a common language and a scalable level of decomposition to effectively evaluate the resilience of interconnected infrastructure systems. AHA models infrastructure systems as directed multidimensional graphs, which enable the evaluation of cross-sector interdependencies prior to, during, and after disruptive events. Finally, and by design, AHA enables the cross linking of data taxonomies to enable more effective data sharing, such as the National Critical Functions (NCF) and Infrastructure Data Taxonomy (IDT).

02 PETROLEUM↗

Uniform Recalibration of Common Spectrophotometry Standard Stars onto the CALSPEC System Using the SuperNova Integral Field Spectrograph

Abstract We calibrate spectrophotometric optical spectra of 32 stars commonly used as standard stars, referenced to 14 stars already on the Hubble Space Telescope–based CALSPEC flux system. Observations of CALSPEC and non-CALSPEC stars were obtained with the SuperNova Integral Field Spectrograph over the wavelength range 3300–9400 Å as calibration for the Nearby Supernova Factory cosmology experiment. In total, this analysis used 4289 standard-star spectra taken on photometric nights. As a modern cosmology analysis, all presubmission methodological decisions were made with the flux scale and external comparison results blinded. The large number of spectra per star allows us to treat the wavelength-by-wavelength calibration for all nights simultaneously with a Bayesian hierarchical model, thereby enabling a consistent treatment of the Type Ia supernova cosmology analysis and the calibration on which it critically relies. We determine the typical per-observation repeatability (median 14 mmag for exposures ≳5 s), the Maunakea atmospheric transmission distribution (median dispersion of 7 mmag with uncertainty 1 mmag), and the scatter internal to our CALSPEC reference stars (median of 8 mmag). We also check our standards against literature filter photometry, finding generally good agreement over the full 12 mag range. Overall, the mean of our system is calibrated to the mean of CALSPEC at the level of ∼3 mmag. With our large number of observations, careful cross-checks, and 14 reference stars, our results are the best calibration yet achieved with an integral-field spectrograph, and among the best calibrated surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗