Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Source term”

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 1,423 records · Page 79

Effects of Low-Carbon Energy Adoption on Airborne Particulate Matter Concentrations With Feedbacks to Future Climate Over California

California plans to reduce emissions of long-lived greenhouse gases (GHGs) through adoption of new energy systems that will also lower concentrations of short-lived absorbing soot contained in airborne particulate matter (PM). Here we examine the direct and indirect effects of reduced PM concentrations under a low-carbon energy (GHG-Step) scenario on radiative forcing in California. Simulations were carried out using the source-oriented WRF/Chem (SOWC) model with 12 km spatial resolution for the year 2054. Furthermore, the avoided aerosol emissions due to technology advances in the GHG-Step scenario reduce ground level PM concentrations by ~8.85% over land compared to the Business-as-Usual (BAU) scenario, but changes to meteorological parameters are more modest. Top-of-atmospheric forcing predicted by the SOWC model increased by 0.15 W m -2 , surface temperature warmed by 0.001 K, and planetary boundary layer height (PBLH) increased by 2.20 cm in the GHG-Step scenario compared to the BAU scenario. PM climate feedbacks are small because the significant changes in ground level PM concentrations associated with the GHG-Step scenario are limited to the first few hundred meters of the atmosphere, with little change for the majority of the vertical column above that level. As an order-of-magnitude comparison, the long-term effects of global reductions in GHG emissions (RCP8.5–RCP4.5) lowered average surface temperature over the California study domain by approximately 0.76 K. The effects of long-lived climate pollutants such as CO 2 are much stronger than the effects of short-lived climate pollutants such as PM soot over California in the year 2054.

54 ENVIRONMENTAL SCIENCES↗

Investigation of potential aerosol transmission and infectivity of SARS-CoV-2 through central ventilation systems

Here we evaluate the concentrations and probabilities of infection for both building interior and exterior exposure sources using a well-mixed model in a connected multizone building. As a central hub of many community and economic activities, buildings provide social connectivity, but the COVID-19 pandemic has reduced social connectivity due to concerns of viral spread within buildings. Although single zone models of infectious spread are well studied, the impact of aerosolized spread of SARS-CoV-2 via air handling systems in multizone buildings remains unexplored. Here we evaluate the influence of filtration, air exchange rates, and the fraction of outdoor air on the probability of infection using the well-known well-mixed modeling approach for a multizone. We find filtration lowers the concentration and probability of infection in both source and connected rooms provided at least some air is recirculated, but that probability is not zero. Filtration has no influence without recirculation or unless the outdoor air contains virus. We find that increasing the air exchange rate removes virus from the source room faster but also increases the rate of exposure to connected rooms. Therefore, slower air exchange rates reduce infectivity in connected rooms at shorter durations, but higher air exchange rates reduce infectivity at longer durations. We further find that when outdoor air is virus free, increasing the fraction of outdoor air is helpful, but, when outdoor air is infective, pathogen exposure inside can persist for hours after a short-term release.

60 APPLIED LIFE SCIENCES↗

Spray and combustion characteristics of pure hydrothermal liquefaction biofuel and mixture blends with diesel

Hydrothermal liquefaction (HTL) is the process in which biological compounds can be converted into biocrude diesel; various feedstocks, such as algae and wet biowaste, can be utilized. Almost 20 million dry tons of wet biowaste is deposited in U.S. municipal wastewaters annually. If recovered and converted into biodiesel, it can act as a massive source of renewable energy. Preliminary experiments on a diesel engine show that HTL fuel can achieve similar performance and emission levels to diesel. In order to further understand the influence of the HTL fuel properties, the HTL fuel and blends with diesel are observed in a constant volume chamber. Previously, 10% and 20% of the HTL fuel were blended with diesel by volume and demonstrated properties that roughly matched diesel in terms of combustion. In this study, 50% and 100% HTL fuel blends are tested on the same apparatus. It was discovered that the workability of HTL is almost the same as that of diesel and is less sensitive to ambient temperatures. When the ambient temperature is reduced from 1200 K to 800 K, the change ratio of ignition delay is only 62% of diesel. By analyzing the combustion process, it can be seen that there is no need to make major modifications to the engine control strategy if HTL fuel is applied to diesel engines. Though, under low temperature and low oxygen conditions, the soot luminosity peak of HTL100 is 3.12 times than that of diesel.

09 BIOMASS FUELS↗

Multiphase Processing of the Water-Soluble and Insoluble Phases of Biomass Burning Organic Aerosol

Biomass burning is one of the most significant sources of organic aerosol in the atmosphere. Biomass burning organic aerosol (BBOA) has been observed to undergo liquid– liquid phase separation (LLPS) to give core–shell morphology with the hydrophobic phase encapsulating the hydrophilic phase, potentially impacting the evolution of light-absorbing components, i.e., brown carbon (BrC), through multiphase processes. Here, we demonstrate how multiphase processing differs between the watersoluble (i.e., hydrophilic) and insoluble (i.e., hydrophobic) phases of BBOA in terms of reactive uptake of ozone in a coated-wall flow tube. Effects of relative humidity (RH) and ultraviolet (UV) irradiation were investigated. Experimental timeseries were used to inform simulations using multilayer kinetic modeling. Among non-irradiated thin films, the uptake coefficient was greatest for the water-soluble phase at 75% RH (3 × 10 –5 , corresponding to a diffusion coefficient of BrC, D BrC , of 3 × 10 –9 cm 2 s –1 ) and least for the same phase at 0% RH (1 × 10 –5 , corresponding to D BrC of 1 × 10 –10 cm 2 s –1 ). The uptake coefficient for the water-insoluble phase fell between these two (about 1.5 × 10 –5 ), regardless of RH, and the corresponding D BrC increased only slightly (8 × 10 –10 cm 2 s –1 at 0% RH to 9 × 10 –10 cm 2 s –1 at 75% RH). The uptake coefficients of both phases at 0% RH decreased significantly after UV irradiation, consistent with a transition from viscous liquid to solid and supported by qualitative microscopy observations. Modeling multiphase ozone oxidation of primary BrC components in the atmosphere demonstrated, first, that LLPS may extend the lifetime of water-soluble BBOA encapsulated by water-insoluble species by a factor of 1.5 at moderate to high RH and, also, that UV irradiation may extend the lifetime of both phases by more than a factor of 2.5.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Open-source simulation platform for air source heat pump integrated with thermal energy storage

Here, this article introduces a modular simulation platform for assessing thermal energy storage (TES) integrated with air source heat pumps (ASHP). The Python platform is an open-source library that includes classes for modeling air-air and air-water heat pumps, TES devices, and the heating load of residential buildings. To validate the ASHP model, the study utilized experimental data obtained from a commercial heat pump evaluated at the National Renewable Energy Laboratory (NREL). The results indicate a mean deviation of 0.7% for COP across the operating range, with a maximum relative deviation of 12.6%. In terms of system heating capacity, the model had an average deviation of 4.3% compared to experimental results, with a maximum deviation of 8.2%. Three classes were implemented for modeling distinct types of TES devices: a generic TES based on energy balances, a sensible isothermal water tank, and a stratified water tank. Details of the mathematical models are provided, along with their respective strengths and limitations. An example is provided showcasing the integration of a residential 10 kWh thermal storage unit with an ASHP operating in Denver, CO. The comparison of two different discharge criteria for the TES unit highlights the importance of control strategies in the system performance.

25 ENERGY STORAGE↗

Comparative scintillation performance of EJ-309, EJ-276, and a novel organic glass

An organic glass scintillator developed by Sandia National Laboratories was characterized in terms of its light output and pulse shape discrimination (PSD) properties and compared to commercial liquid (EJ-309) and plastic (EJ-276) organic scintillators. The electron light output was determined through relative comparison of the 137 Cs Compton edge location. The proton light yield was measured using a double time-of-flight technique at the 88-Inch Cyclotron at Lawrence Berkeley National Laboratory. Using a tunable broad-spectrum neutron source and an array of pulse-shape-discriminating observation scintillators, a continuous measurement of the proton light yield was performed for EJ-309 (200 keV–3.2 MeV), EJ-276 (170 keV–4.9 MeV), and the organic glass (50 keV–20 MeV) . Finally, the PSD properties of the organic glass, EJ-309, and EJ-276 were evaluated using an AmBe source and compared via a figure-of-merit metric. The organic glass exhibited a higher electron light output than both EJ-309 and EJ-276. Its proton light yield and PSD performance were comparable to EJ-309 and superior to that of EJ-276. Finally, with these performance characteristics, the organic glass scintillator is well poised to replace current state-of-the-art PSD-capable scintillators in a range of fast neutron detection applications.

47 OTHER INSTRUMENTATION↗

The Energetics of the Central Engine in the Powerful Quasar 3C 298

The compact steep-spectrum radio source 3C 298 (redshift of 1.44) has the largest 178 MHz luminosity in the Third Cambridge Revised Catalogue (3CR); its radio lobes are among the most luminous in the universe. The plasma state of the radio lobes is modeled with the aid of interferometric radio observations (in particular, the new Low Frequency Array observation and archival MERLIN data) and archival single-station data. It is estimated that the long-term time-averaged jet power required to fill these lobes with leptonic plasma is $\overline{Q}\,\approx 1.28\pm 0.51\times {10}^{47}\,\mathrm{erg}\,{{\rm{s}}}^{-1}$, rivaling the largest time-averaged jet powers from any quasar. Supporting this notion of extraordinary jet power is a 0.5–10 keV luminosity of ≈ 5.2 × 10 46 erg s –1 , comparable to luminous blazars, yet there is no other indication of strong relativistic beaming. We combine two new high signal-to-noise ratio optical spectroscopic observations from the Hobby-Eberly Telescope with archival Hubble Space Telescope, Two Micron All Sky Survey, and Galaxy Evolutionary Explorer data to compute a bolometric luminosity from the accretion flow of L bol ≈ 1.55 ± 0.15 × 10 47 erg s –1 . The ratio, $\overline{Q}/{L}_{\mathrm{bol}}\approx 1$, is the approximate upper limit for quasars. Characteristic of a large $\overline{Q}/{L}_{\mathrm{bol}}$, we find an extreme-ultraviolet (EUV) spectrum that is very steep (the "EUV deficit" of powerful radio quasars relative to radio-quiet quasars), and this weak ionizing continuum is likely a contributing factor to the relatively small equivalent widths of the broad emission lines in this quasar.

79 ASTRONOMY AND ASTROPHYSICS↗

Power at Sea: Use of Wave Energy Converters for Powering Offshore Oil and Gas Operations

New innovative services and goods based on artificial intelligence, next-generation robotics, autonomy, digitization, and electrification are taking root in the ocean economy. More typical competitive forces, such as value enhancement and safety improvement, are also present, as suppliers constantly strive to improve their competitive positioning. Concurrently, the Oil and Gas sector is experiencing mounting pressure to decrease its operational carbon emissions, forcing efforts to reduce or eliminate liquid fuel consumption offshore. Companies in the ocean economy are having to react to this combination of macro- and micro-level drivers and innovate at an ever more rapid pace. However, before many of these new capabilities can be implemented, new, clean, reliable power sources are needed, especially as liquid fuel grows increasingly out of favor. For many of these loads—both old and new and ranging from watts to megawatts—wave energy systems provide the best solution for remote power generation. As such, research, development and demonstration activities are underway to prove wave energy's ability to provide reliable, consistent energy supplies for certain offshore Oil and Gas activities, paving the way for near term commercialization. Furthermore, this paper outlines these drivers, competitive reaction to them, and why and how wave energy can be a preferable choice for power at sea.

16 TIDAL AND WAVE POWER↗

Measurement report: A comparison of ground-level ice-nucleating-particle abundance and aerosol properties during autumn at contrasting marine and terrestrial locations

Abstract. Ice-nucleating particles (INPs) are an essential class of aerosols found worldwide that have far-reaching but poorly quantified climate feedback mechanisms through interaction with clouds and impacts on precipitation. These particles can have highly variable physicochemical properties in the atmosphere, and it is crucial to continuously monitor their long-term concentration relative to total ambient aerosol populations at a wide variety of sites to comprehensively understand aerosol–cloud interactions in the atmosphere. Hence, our study applied an in situ forced expansion cooling device to measure ambient INP concentrations and test its automated continuous measurements at atmospheric observatories, where complementary aerosol instruments are heavily equipped. Using collocated aerosol size, number, and composition measurements from these sites, we analyzed the correlation between sources and abundance of INPs in different environments. Toward this aim, we have measured ground-level INP concentrations at two contrasting sites, one in the Southern Great Plains (SGP) region of the United States with a substantial terrestrially influenced aerosol population and one in the Eastern North Atlantic Ocean (ENA) region with a primarily marine-influenced aerosol population. These measurements examined INPs mainly formed through immersion freezing and were performed at a ≤ 12 min resolution and with a wide range of heterogeneous freezing temperatures (Ts above −31 °C) for at least 45 d at each site. The associated INP data analysis was conducted in a consistent manner. We also explored the additional offline characterization of ambient aerosol particle samples from both locations in comparison to in situ data. From our ENA data, on average, INP abundance ranges from ≈ 1 to ≈ 20 L−1 (−30 °C ≤ T ≤ −20 °C) during October–November 2020. Backward air mass trajectories reveal a strong marine influence at ENA with 75.7 % of air masses originating over the Atlantic Ocean and 96.6 % of air masses traveling over open water, but analysis of particle chemistry suggests an additional INP source besides maritime aerosols (e.g., sea spray aerosols) at ENA. In contrast, 90.8 % of air masses at the SGP location originated from the North American continent, and 96.1 % of the time, these air masses traveled over land. As a result, organic-rich SGP aerosols from terrestrial sources exhibited notably high INP abundance from ≈ 1 to ≈ 100 L−1 (−30 °C ≤ T ≤ −15 °C) during October–November 2019. The probability density function of aerosol surface area-scaled immersion freezing efficiency (ice nucleation active surface site density; ns) was assessed for selected freezing temperatures. While the INP concentrations measured at SGP are higher than those of ENA, the ns(T) values of SGP (≈ 105 to ≈ 107 m−2 for −30 °C ≤ T ≤ −15 °C) are reciprocally lower than ENA for approximately 2 orders of magnitude (≈ 107 to ≈ 109 m−2 for −30 °C ≤ T ≤ −15 °C). The observed difference in ns(T) mainly stems from varied available aerosol surface areas, Saer, from two sites (Saer,SGP > Saer,ENA). INP parameterizations were developed as a function of examined freezing temperatures from SGP and ENA for our study periods.

54 ENVIRONMENTAL SCIENCES↗

Measurement and calibration of non-linear shear terms in galaxy cluster fields

Introduction Galaxy cluster lensing is a powerful tool for measuring the mass of galaxy clusters, but accurate shear measurement and calibration are critical to obtaining reliable results. This study focuses on the measurement and calibration of weak lensing shears to improve mass estimates in cluster lensing. To deal with the problem, we first developed an image simulation pipeline, jedisim, which utilizes galaxy images extracted from the Hubble Space Telescope (HST) Ultra Deep Field (UDF) and the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Methods The simulations represent realistic galaxy distributions and morphologies as input sources. The foreground halo with a Navarro–Frenk–White (NFW) profile is constructed such that the lensing signals of background galaxies can be measured by the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST) Science Pipelines. By comparing the measured reduced shear g meas and the true reduced shear g true , we observe non-linearity up to g ≲ 0.6 . We fit polynomials to the data with quadratic correction adequate to g ≲ 0.4 . Meanwhile, we conduct mass estimates using the pzmassfitter code on four different clusters. Results The mass estimate results are significantly improved after applying the shear calibration derived from the present work—from 4.954 ± 0.504 × 1 0 14 M ⊙ to 10.507 ± 0.498 × 1 0 14 M ⊙ after calibration for a simulated cluster with the mass of 10 × 1 0 14 M ⊙ . In multiple cases of validation, the estimated results are all consistent with true cluster mass. Discussion This study yields the first relationship between reality and shape measurement of the LSST Science Pipelines and serves as the first step toward the overall goal of mass calibration in cluster lensing. By addressing the challenges in shear measurement and calibration, we aim to enhance the accuracy and reliability of mass estimates in galaxy cluster lensing studies.

Astronomy & Astrophysics↗

Lignin-derived electrochemical energy materials and systems

Electrochemical energy storage systems such as supercapacitors, rechargeable batteries and fuel cells have been proven the most effective technologies for energy conversion, storage, and management at different scales. Although a large number of electrochemical energy technologies have been developed in the past and they will continue to be optimized in terms of cost, lifetime, and performance, there is a substantial growing demand for advanced electrochemical energy systems. To deploy these advanced systems, the electrode and electrolyte materials with higher performance, longer life, and lower cost, must be developed. Lignin is the second most abundant natural polymer after cellulose, a byproduct from emerging cellulosic biorefineries, and a waste product from pulp and paper industries. Numerous researches have successfully demonstrated that lignin from different sources can be used as precursors or feedstocks for preparing high-performance electrochemical energy materials and components such as electrodes, electrolytes, membrane separators, and additives. Moreover, techno-economic analyses indicate that it is possible to prepare cost-effective carbons from lignin at engineering scales, compared to current carbon products. These facts suggest that scalable conversion of lignin into high-value energy materials will offer a promising pathway to not only promote the utilization and valorization of lignin but also boost the development of the advanced electrochemical energy systems. This review presents state of the arts of renewable energy materials derived from various lignin and their applications in electrochemical energy systems with emphasis on supercapacitors, rechargeable batteries, and fuel cells. Meanwhile, this article also aims to carve out the critical barriers for lignin-derived high-performance materials for energy applications, intending to identify viable approaches for synthesis of sustainable new energy materials.

09 BIOMASS FUELS↗

A long-term pilot-scale algal cultivation on power plant flue gas – Cultivation stability and biomass accumulation

A long-term culture of a selected freshwater microalgal assemblage on power plant flue gas was realized using a pilot-scale photobioreactor. Different culture conditions of light cycle, harvesting amount, pH stabilization, and sulfur addition were tested during a 12-month continuous operation at a natural gas fired power plant. The results concluded that the daily harvesting amount of 60% with no light-cycle and NH 4 NO 3 as the nitrogen source and pH stabilization agent, resulted in optimal cultivation performance. Biomass concentration and biomass productivity reached 1.07 ± 0.23 g L -1 and 0.43 ± 0.14 g L -1 d -1 , respectively. During the year-long cultivation period, the algal assemblage was relatively stable with average abundances of 88 ± 4% and 12 ± 4% for algae and bacteria, respectively, and no significant differences between individual case scenarios. The stable culture fostered a stable biomass composition. Protein, carbohydrate, and lipid contents of the algal biomass were 58.60 ± 2.41%, 19.5 ± 2.44%, and 9.70 ± 1.68%, respectively.

09 BIOMASS FUELS↗

Comprehensive assessment of deep reinforcement learning approaches for economic dispatch in nuclear-driven microgrids

As the electrical grid integrates more variable renewable energy sources such as wind and solar, the demand for distributed and flexible systems to address this increased variability becomes critical. Nuclear-driven microgrids provide a promising solution by offering stable generation to complement intermittent renewables, ensuring grid reliability and operating efficiency. This paper proposes a recurrent deep reinforcement learning framework for optimal economic dispatch in a nuclear-powered microgrid integrating renewable energy sources, small modular reactors, battery storage systems, and balance-of-plant dynamics. A three-agent control architecture is developed, where demand and renewable energy agents act as forecasters, and a reinforcement learning-based dispatch agent performs real-time energy allocation. A nonlinear programming formulation is first used to generate an optimal baseline for benchmarking. The proposed dispatch controller, based on Proximal Policy Optimization enhanced with Long Short-Term Memory networks, exploits temporal correlations in system dynamics by taking advantage of the time series used as inputs to improve policy robustness under uncertainty. Comparative analysis against established deep reinforcement learning methods, including Proximal Policy Optimization with a feedforward architecture, Soft Actor-Critic, and Twin Delayed Deep Deterministic Policy Gradient, demonstrates superior performance. Numerical results indicate that the proposed controller achieves a 0.39% cost reduction relative to the nonlinear programming benchmark and outperforms other learning-based methods by generating additional revenue of up to 0.35%. All reinforcement learning controllers compute dispatch actions in less than 0.3 s, resulting in a computational speedup of more than three orders of magnitude over the nonlinear programming baseline. The findings of this paper highlight their applicability for real-time operation and control in nuclear-integrated microgrids under volatile operating conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Controls on the 36 Cl/Cl input ratio of paleo-groundwater in arid environments: New evidence from 81 Kr/Kr data

In this study, measurements of the long-lived 81 Kr and 36 Cl radioisotopes in groundwater from the Negev Desert (Israel) were used to assess the 36 Cl/Cl input ratios and Cl – contents for paleorecharge into the Nubian Sandstone Aquifer (NSA). The reconstructed Cl – content of the recharge flux was on the order of 300–400 mg/L. An initial 36 Cl/Cl ratio of 50 × 10 –15 was assessed for the groundwater replenishment in the Negev Desert since the late Pleistocene, in agreement with the 36 Cl/Cl ratios in recent local rainwater. This is despite possible changes in the climatic conditions and the 36 Cl production rates in the atmosphere over this timeframe. This similarity in values is explained by the major role played by the erosion and weathering of near-surface materials in the desert environment that dominate the hydrochemistry of rains, floods, and the consequent groundwater recharge. Spatial variation in the reconstructed initial 36 Cl/Cl ratio is accounted for by the differences in the mineral aerosol sources for specific recharge areas of the NSA. Accordingly, regional variations in the 36 Cl/Cl input in groundwater reservoirs surrounding the Mediterranean Sea indicate various processes that govern the 36 Cl/Cl system. Finally, the results of this study highlight the great advantage of integrating 81 Kr age information in evaluating the initial 36 Cl/Cl and Cl – input, which is essential for the calibration of 36 Cl radioisotope as an available long-term dating tool for a given basin.

54 ENVIRONMENTAL SCIENCES↗

Immunophenotypic characterization and therapeutics effects of human bone marrow- and umbilical cord-derived mesenchymal stromal cells in an experimental model of sepsis

Highlights: • Mesenchymal stromal cells (MSCs) may be a therapeutic option for sepsis. • Bone marrow (BM) vs. umbilical cord (UC) -MSCs induce different M1/M2 gene profiles. • In a CLP model of sepsis, only BM-MSCs improved survival and bacterial clearance. Sepsis is a complicated multi-system disorder characterized by a dysregulated host response to infection. Despite substantial progress in the understanding of mechanisms of sepsis, translation of these advances into clinically effective therapies remains challenging. Mesenchymal Stromal Cells (MSCs) possess immunomodulatory properties that have shown therapeutic promise in preclinical models of sepsis. The therapeutic effects of MSCs may vary depending on the source and type of these cells. In this comparative study, the gene expression pattern and surface markers of bone marrow-derived MSCs (BM-MSCs) and umbilical cord-derived MSCs (UC-MSCs) as well as their therapeutic effects in a clinically relevant mouse model of polymicrobial sepsis, cecal ligation and puncture (CLP), were investigated. The results showed remarkable differences in gene expression profile, surface markers and therapeutic potency in terms of enhancing survival and pro/anti-inflammatory responses between the two MSC types. BM-MSCs improved survival concomitant with an enhanced systemic bacterial clearance and improved inflammatory profile post CLP surgery. Despite some improvement in the inflammatory profile of the septic animals, treatment with UC-MSCs did not enhance survival or bacterial clearance. Overall, the beneficial therapeutic effects of BM-MSCs over UC-MSCs may likely be attributed to their pro-inflammatory function, and to some extent anti-inflammatory features, reflected in their gene expression pattern enhancing macrophage polarization to M1/M2 phenotypes resulting in a balanced pro- and anti-inflammatory response against polymicrobial sepsis.

60 APPLIED LIFE SCIENCES↗

Modeling the Dynamics of Energetic Protons in Earth's Inner Magnetosphere

Recently, Van Allen Probes data have revealed the interesting energy-dependent dynamics of energetic protons in Earth's inner magnetosphere. To quantify the role of radial diffusion and charge exchange to the observed energy-dependent dynamics of ring current protons, a 1D radial diffusion model with charge exchange loss is implemented. The observed proton flux over the long-term period of November 2012 to September 2013 is first converted to phase space densities, which are then simulated by our model driven by an outer boundary condition at L* = 5.5 derived from the data. The simulation results show that our model generally captures the transport and acceleration of energetic protons at μ = 30, 50, and 80 MeV/G and K = 0.11 G 1/2 R E , suggesting that radial diffusion is the dominant source mechanism for >75 keV protons at L* = 3.5–5.5. Additionally, the observed fast decay of protons at lower μ and slow decay at higher μ are well captured by the model, demonstrating the dominant role of charge exchange in explaining the observed energy-dependent proton decay. For higher μ protons, prompt losses of protons on the time scale of hours are observed over a wide range of L*, which is too fast to be explained by charge exchange. Some of the prompt losses at high L* regions are reproduced by the model with outward radial diffusion to the outer boundary. However, many prompt losses observed at lower L* regions are not captured by the model, which could be due to other loss mechanisms including electromagnetic ion cyclotron wave scattering and field line curvature scattering.

79 ASTRONOMY AND ASTROPHYSICS↗

QCOR; A Language Extension Specification for the Heterogeneous Quantum-Classical Model of Computation

Quantum computing (QC) is an emerging computational paradigm that leverages the laws of quantum mechanics to perform elementary logic operations. Existing programming models for QC were designed with fault-tolerant hardware in mind, envisioning stand-alone applications. However, the susceptibility of near-term quantum computers to noise limits their stand-alone utility. To better leverage limited computational strengths of noisy quantum devices, hybrid algorithms have been suggested whereby quantum computers are used in tandem with their classical counterparts in a heterogeneous fashion. This modus operandi calls out for a programming model and a high-level programming language that natively and seamlessly supports heterogeneous quantum-classical hardware architectures in a single-source-code paradigm. Motivated by the lack of such a model, we introduce a language extension specification, called QCOR, which enables single-source quantum-classical programming. Programs written using the QCOR library–based language extensions can be compiled to produce functional hybrid binary executables. After defining QCOR’s programming model, memory model, and execution model, we discuss how QCOR enables variational, iterative, and feed-forward QC. Additionally, QCOR approaches quantum-classical computation in a hardware-agnostic heterogeneous fashion and strives to build on best practices of high-performance computing. The high level of abstraction in the language extension is intended to accelerate the adoption of QC by researchers familiar with classical high-performance computing.

97 MATHEMATICS AND COMPUTING↗

Scalable Pattern Matching in Metadata Graphs via Constraint Checking

Pattern matching is a fundamental tool for answering complex graph queries. Unfortunately, existing solutions have limited capabilities: They do not scale to process large graphs and/or support only a restricted set of search templates or usage scenarios. Moreover, the algorithms at the core of the existing techniques are not suitable for today’s graph processing infrastructures relying on horizontal scalability and shared-nothing clusters, as most of these algorithms are inherently sequential and difficult to parallelize. In this article we present an algorithmic pipeline that bases pattern matching on constraint checking. The key intuition is that each vertex and edge participating in a match has to meet a set of constraints implicitly specified by the search template. These constraints can be verified independently and typically are less expensive to compute than searching the full template. The pipeline we propose generates these constraints and iterates over them to eliminate all the vertices and edges that do not participate in any match, thus reducing the background graph to a subgraph that is the union of all template matches—the complete set of all vertices and edges that participate in at least one match. Additional analysis can be performed on this annotated, reduced graph, such as full match enumeration, match counting, or computing vertex/edge centrality. Furthermore, a vertex-centric formulation for constraint checking algorithms exists, and this makes it possible to harness existing high-performance, vertex-centric graph processing frameworks. This technique (i) enables highly scalable pattern matching in metadata (labeled) graphs; (ii) supports arbitrary patterns with 100% precision; (iii) enables tradeoffs between precision and time-to-solution, while always selects all vertices and edges that participate in matches, thus offering 100% recall; and (iv) supports a set of popular data analytics scenarios. We implement our approach on top of HavoqGT, an open-source asynchronous graph processing framework, and demonstrate its advantages through strong and weak scaling experiments on massive scale real-world (up to 257 billion edges) and synthetic (up to 4.4 trillion edges) labeled graphs, respectively, and at scales (1,024 nodes / 36,864 cores), orders of magnitude larger than used in the past for similar problems. This article serves two purposes: First, it synthesises the knowledge accumulated during a long-term project. Second, it presents new system features, usage scenarios, optimizations, and comparisons with related work that strengthen the confidence that pattern matching based on iterative pruning via constraint checking is an effective and scalable approach in practice. The new contributions include the following: (i) We demonstrate the ability of the constraint checking approach to efficiently support two additional search scenarios that often emerge in practice, interactive incremental search and exploratory search. (ii) We empirically compare our solution with two additional state-of-the-art systems, Arabsque and TriAD. (iii) We show the ability of our solution to accommodate a more diverse range of datasets with varying properties, e.g., scale, skewness, label distribution, and match frequency. (iv) We introduce or extend a number of system features (e.g., work aggregation, load balancing, and the ability to cap the generated traffic) and design optimizations and demonstrate their advantages with respect to improving performance and scalability. (v) We present bottleneck analysis and insights into artifacts that influence performance. (vi) We present a theoretical complexity argument that motivates the performance gains we observe.

97 MATHEMATICS AND COMPUTING↗