Search NASA⌕ Search

SEARCH · Search NASA

Results for “direct counts”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Radiation Mapping for an Unmanned Aerial Vehicle: Development and Simulated Testing of Algorithms for Source Mapping and Navigation Path Generation

Image reconstruction algorithms were developed for radiation source mapping and used for generating the search path of a moving radiation detector, such as one onboard an unmanned aerial vehicle. Simulations consisted of first assuming radioactive sources of varying complexity and estimating the radiation fields that would then be produced by that source distribution. Next, the "measurements" that would result from a pair of adjacent spatial locations were computed. A crude estimate of the source distribution likely to have produced such "measurements" was reconstructed based upon the limited measurements. Location of the next "measurement" was then determined as halfway between the location of the estimated source and the current "measurement." With each additional sample, improved source distribution reconstructions were made and used to inform the immediate direction of detector motion. Source reconstruction or mapping was formulated as an inverse problem solved with either maximum a posteriori or least squares (LS) regression deconvolution methods. Different amounts of noise were added to the simulated "measurements," allowing evaluation of the methods' performances as functions of signal-to-noise ratio of the measured map. As expected, methods that promote sparsity were better suited in reconstructing point sources. Reliable prior information of the source distribution also improved the reconstruction results, especially with distributed sources. With a non-negative least square algorithm and the suggested paths it generated, location of sources was successfully estimated to an accuracy of 0.014 m within nine iterations in a single-source scenario and 12 iterations in a two-source scenario, given a 10% error on the integrated counts and a Poisson distribution of the noise associated with the measured counts.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

Benchmarking an Embedded Adaptive Sampling Configuration Interaction Method for Surface Reactions: H 2 Desorption from and CH 4 Dissociation on Cu(111)

Embedded (emb-) correlated wavefunction (CW) theory enables accurate assessments of both ground- and excited-state reaction mechanisms involved in heterogeneous catalysis. Embedded multireference second-order perturbation theory (emb-MRPT2) based on reference wavefunctions generated via embedded complete active space self-consistent field (emb-CASSCF) theory is currently state-of-the-art. However, the factorial scaling of CASSCF limits the size of active space and the complexity of systems that can be studied. In this work, we assess the efficacy of an alternative CW method, adaptive sampling configuration interaction (ASCI)–which enables large active spaces to be used–for studying surface reactions. We couple ASCI with density functional embedding theory (DFET) and benchmark its performance for two reactions: H 2 desorption from and CH 4 dissociation on the Cu(111) surface. Unlike embedded complete active space second-order perturbation theory (emb-CASPT2) that accurately reproduces a measured H 2 desorption barrier, embedded ASCI, using a very large active space (though one that still comprises a small portion of the full set of orbitals) fails to do so. Adding an extra correlation term from embedded Møller–Plesset second-order perturbation theory (emb-MP2) improves the desorption barrier and endothermicity predictions. Thus, the inaccuracy of embedded ASCI comes from the missing dynamic correlation from the many other electrons and orbitals not included in the active space. For CH 4 dissociation, again embedded ASCI overestimates the dissociation barrier compared to emb-CASPT2 predictions. Adding dynamic correlation from emb-MP2 helps correct the barrier. However, this composite approach suffers from double counting of correlation within embedded ASCI followed by emb-MP2 calculations. We therefore conclude that the state-of-the-art emb-MRPT2 based on reference wavefunctions generated via emb-CASSCF remains the method of choice for studying surface reactions. emb-ASCI is useful when large active spaces beyond the limit of emb-CASSCF are essential, such as to study complex surface reactions with significant multiconfigurational character (static correlation) but weak dynamic correlation.

08 HYDROGEN↗
DOE OSTI2022

A fully Lagrangian, non-parametric bias model for dark matter halos

Here, we present a non-parametric Lagrangian biasing model and fit the ratio of the halo and mass densities at the field level using the mass-weighted halo field in the AbacusSummit simulations at z=0.5. Unlike the perturbative halo bias model that has been widely used in interpreting the observed large-scale structure traced by galaxies, we find a non-negative halo-to-mass ratio that increases monotonically with the linear overdensity δ 1 in the initial Lagrangian space. The bias expansion, however, does not guarantee non-negativity of the halo counts, and may lead to rising halo number counts at negative overdensities. The shape of the halo-to-mass ratio is unlikely to be described by a polynomial function of δ 1 and other quantities. Especially for massive halos with 6×10 12 h -1 M ⊙ , the halo-to-mass ratio starts soaring up at δ 1 >0, substantially different from the predictions of the bias expansion. We show that for the halo masses we consider (M>3×10 11 h -1 M ⊙ ) a non-parametric halo-to-mass ratio as a function of δ 1 and its local derivative ∇^2δ 1 can recover the halo power spectra to sub-percent level accuracy for wavenumbers k=0.01-0.1 h Mpc -1 given a proper smoothing scale to filter the initial density field, even though we do not fit the power spectrum directly. However, there is mild dependence of the recovery of the halo power spectrum on the smoothing scale and other input parameters. At k<0.01 h Mpc -1 and for massive halos with M>6x10 12 h -1 M ⊙ , our non-parametric model leads to a few percent overestimation of the halo power spectrum, indicating the need for larger or multiple smoothing scales. The halo-to-mass ratios obtained qualitatively agree with intuitions from extended Press-Schechter theory. We compare our framework to the bias expansion and discuss possible extensions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Stabilization of the 81-channel coherent beam combination using machine learning

We develop a rapidly converging algorithm for stabilizing a large channel-count diffractive optical coherent beam combination. An 81-beam combiner is controlled by a novel, machine-learning based, iterative method to correct the optical phases, operating on an experimentally calibrated numerical model. A neural-network is trained to detect phase errors based on interference pattern recognition of uncombined beams adjacent to the combined one. Due to the non-uniqueness of solutions in the full space of possible phases, the network is trained within a limited phase perturbation/error range. This also reduces the number of samples needed for training. Simulations have proven that the network can converge in one step for small phase perturbations. When the trained neural-network is applied to a realistic case of 360 degree full range, an iterative scheme exploits random walking at the beginning, with the accuracy of prediction on phase feedback direction, to allow the neural-network to step into the training range for fast convergence. This neural-network-based iterative method of phase detection works tens of times faster than the commonly used stochastic parallel gradient descent approach (SPGD) using a single-detector and random dither when both are tested with random phase perturbations.

Wang, Dan↗
DOE OSTI2026

Topology-Dependent Performance of Free-Space Photonic Quantum Networks Under Noise

Photonic quantum communication enables secure and high-fidelity information transfer beyond classical limits, with direct relevance to emerging quantum networks operating in free-space environments. While physical-layer models of depolarizing noise, Gamma–Gamma turbulence statistics, entanglement swapping, and decoy-state QKD security bounds are individually well established, prior work typically treats these components in isolation or under fixed network assumptions. In this work, we develop a unified topology-aware analytical framework that simultaneously integrates free-space optical link budgets, turbulence-induced visibility degradation, depolarizing qubit noise, multi-hop entanglement cascade dynamics, teleportation fidelity thresholds, CHSH nonlocality certification, and asymptotic decoy-state secret key rate bounds across star, mesh, and ring graph structures. Rather than introducing new physical channel models, we demonstrate that identical physical links exhibit fundamentally different end-to-end performance once embedded within different network topologies. Mesh architectures minimize visibility cascade through hop-count reduction but incur quadratic hardware scaling. Star topologies minimize link count but concentrate noise and synchronization overhead at the hub. Ring configurations offer linear hardware scaling with multiplicative fidelity degradation. The results establish topology as a first-order design parameter in near-term free-space quantum networks operating without full quantum repeater infrastructures. While motivated by distributed multi-agent architectures, the framework applies broadly to terrestrial, airborne, and satellite-assisted photonic quantum communication systems.

QKD↗
DOE OSTI2025

Global Methane Budget 2000–2020

Abstract. Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. CH4 is the second most important human-influenced greenhouse gas in terms of climate forcing after carbon dioxide (CO2), and both emissions and atmospheric concentrations of CH4 have continued to increase since 2007 after a temporary pause. The relative importance of CH4 emissions compared to those of CO2 for temperature change is related to its shorter atmospheric lifetime, stronger radiative effect, and acceleration in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in quantifying the factors responsible for the observed atmospheric growth rate arise from diverse, geographically overlapping CH4 sources and from the uncertain magnitude and temporal change in the destruction of CH4 by short-lived and highly variable hydroxyl radicals (OH). To address these challenges, we have established a consortium of multidisciplinary scientists under the umbrella of the Global Carbon Project to improve, synthesise, and update the global CH4 budget regularly and to stimulate new research on the methane cycle. Following Saunois et al. (2016, 2020), we present here the third version of the living review paper dedicated to the decadal CH4 budget, integrating results of top-down CH4 emission estimates (based on in situ and Greenhouse Gases Observing SATellite (GOSAT) atmospheric observations and an ensemble of atmospheric inverse-model results) and bottom-up estimates (based on process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations). We present a budget for the most recent 2010–2019 calendar decade (the latest period for which full data sets are available), for the previous decade of 2000–2009 and for the year 2020. The revision of the bottom-up budget in this 2025 edition benefits from important progress in estimating inland freshwater emissions, with better counting of emissions from lakes and ponds, reservoirs, and streams and rivers. This budget also reduces double counting across freshwater and wetland emissions and, for the first time, includes an estimate of the potential double counting that may exist (average of 23 Tg CH4 yr−1). Bottom-up approaches show that the combined wetland and inland freshwater emissions average 248 [159–369] Tg CH4 yr−1 for the 2010–2019 decade. Natural fluxes are perturbed by human activities through climate, eutrophication, and land use. In this budget, we also estimate, for the first time, this anthropogenic component contributing to wetland and inland freshwater emissions. Newly available gridded products also allowed us to derive an almost complete latitudinal and regional budget based on bottom-up approaches. For the 2010–2019 decade, global CH4 emissions are estimated by atmospheric inversions (top-down) to be 575 Tg CH4 yr−1 (range 553–586, corresponding to the minimum and maximum estimates of the model ensemble). Of this amount, 369 Tg CH4 yr−1 or ∼ 65 % is attributed to direct anthropogenic sources in the fossil, agriculture, and waste and anthropogenic biomass burning (range 350–391 Tg CH4 yr−1 or 63 %–68 %). For the 2000–2009 period, the atmospheric inversions give a slightly lower total emission than for 2010–2019, by 32 Tg CH4 yr−1 (range 9–40). The 2020 emission rate is the highest of the period and reaches 608 Tg CH4 yr−1 (range 581–627), which is 12 % higher than the average emissions in the 2000s. Since 2012, global direct anthropogenic CH4 emission trends have been tracking scenarios that assume no or minimal climate mitigation policies proposed by the Intergovernmental Panel on Climate Change (shared socio-economic pathways SSP5 and SSP3). Bottom-up methods suggest 16 % (94 Tg CH4 yr−1) larger global emissions (669 Tg CH4 yr−1, range 512–849) than top-down inversion methods for the 2010–2019 period. The discrepancy between the bottom-up and the top-down budgets has been greatly reduced compared to the previous differences (167 and 156 Tg CH4 yr−1 in Saunois et al. (2016, 2020) respectively), and for the first time uncertainties in bottom-up and top-down budgets overlap. Although differences have been reduced between inversions and bottom-up, the most important source of uncertainty in the global CH4 budget is still attributable to natural emissions, especially those from wetlands and inland freshwaters. The tropospheric loss of methane, as the main contributor to methane lifetime, has been estimated at 563 [510–663] Tg CH4 yr−1 based on chemistry–climate models. These values are slightly larger than for 2000–2009 due to the impact of the rise in atmospheric methane and remaining large uncertainty (∼ 25 %). The total sink of CH4 is estimated at 633 [507–796] Tg CH4 yr−1 by the bottom-up approaches and at 554 [550–567] Tg CH4 yr−1 by top-down approaches. However, most of the top-down models use the same OH distribution, which introduces less uncertainty to the global budget than is likely justified. For 2010–2019, agriculture and waste contributed an estimated 228 [213–242] Tg CH4 yr−1 in the top-down budget and 211 [195–231] Tg CH4 yr−1 in the bottom-up budget. Fossil fuel emissions contributed 115 [100–124] Tg CH4 yr−1 in the top-down budget and 120 [117–125] Tg CH4 yr−1 in the bottom-up budget. Biomass and biofuel burning contributed 27 [26–27] Tg CH4 yr−1 in the top-down budget and 28 [21–39] Tg CH4 yr−1 in the bottom-up budget. We identify five major priorities for improving the CH4 budget: (i) producing a global, high-resolution map of water-saturated soils and inundated areas emitting CH4 based on a robust classification of different types of emitting ecosystems; (ii) further development of process-based models for inland-water emissions; (iii) intensification of CH4 observations at local (e.g. FLUXNET-CH4 measurements, urban-scale monitoring, satellite imagery with pointing capabilities) to regional scales (surface networks and global remote sensing measurements from satellites) to constrain both bottom-up models and atmospheric inversions; (iv) improvements of transport models and the representation of photochemical sinks in top-down inversions; and (v) integration of 3D variational inversion systems using isotopic and/or co-emitted species such as ethane as well as information in the bottom-up inventories on anthropogenic super-emitters detected by remote sensing (mainly oil and gas sector but also coal, agriculture, and landfills) to improve source partitioning. The data presented here can be downloaded from https://doi.org/10.18160/GKQ9-2RHT (Martinez et al., 2024).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Energy Exascale Earth System Model v2.0.2

Second patch release of v2.0.0 Changes since v2.0.1 [Important change] Add and update SSP370 and SSP585 cases, add tests, fix use-case files, Change ocean and sea-ice IC for ARRM60to10. [EAM] allow thetaxx as a CAM_TARGET, enable northamericax4v1pg2_WC14to60E2r3 for AMIP, fix ndrop initialization, allow up to 15 history files. [EAMxx] allow use of readnl [EAM-MMF] allow transient SST case for C++ MMF, remove specific task/thread count for ESMT test, modify the variance transport diagnostic, reorg tests, [HOMME] Allow optional sponge layer. [ELM] add 2 land-atm compsets, update FATES to API 17.0.0, allow FATES sp mode, allow ELM harvest to drive FATES harvest, fix memleaks, reduce test build times, modified parameters for miscanthus and switchgrass based on calibration [MPAS-Ocean] turn on ocean BGC in BGC cases, fix interface locations for 60L PHC grid, add mode spec to conservation check streams, allow oRRS18to6v3 grid to run with JRA, update ocean and sea ice ICs for ARRM60to10 (needed for v2), GPU port of thickness tendency, add CMPASO-JRA1p4, [MPAS framework] Add new reproducible global sum module for MPAS components, Add mostRecentAccessTime attribute to streams [MPAS-landice] Update MALI version and Greenland mesh [MPAS-seaice] add single-cell test case, fix BGC restart, Adds omp critical directives for ice warnings seen on cori. [CIME] Fix component namelist creation bugs when NINST>1, add OpenACC, OpenMP and CXX GPU tests, fix nonBFB tests, stop using config_compilers, reduce ELM test build times with shared executable, update OpenMPI on Chrysalis. fix baseline handling, provenance handling with update to cime6.0.33, also use component-specific config_pes files and add the ones from master [Machine updates] cori modules after maint, Chrysalis to OpenMPI-4.1.3 [run_e3sm] replace default case name and group. [Externals] update SCORPIO to 1.3.2

ECP↗
DOE OSTI2020

Computationally Directed Discovery of MoBi 2

Incorporating bismuth, the heaviest element stable to radioactive decay, into new materials enables the creation of emergent properties such as permanent magnetism, superconductivity, and nontrivial topology. Understanding the factors that drive Bi reactivity is critical for the realization of these properties. Using pressure as a tunable synthetic vector, we can access unexplored regions of phase space to foster reactivity between elements that do not react under ambient conditions. Furthermore, combining computational and experimental methods for materials discovery at high-pressures provides broader insight into the thermodynamic landscape than can be achieved through experiment alone, informing our understanding of the dominant chemical factors governing structure formation. Herein, we report our combined computational and experimental exploration of the Mo–Bi system, for which no binary intermetallic structures were previously known. Using the ab initio random structure searching (AIRSS) approach, we identified multiple synthetic targets between 0–50 GPa. Highpressure in situ powder X-ray diffraction experiments performed in diamond anvil cells confirmed that Mo–Bi mixtures exhibit rich chemistry upon the application of pressure, including experimental realization of the computationally predicted CuAl 2 -type MoBi 2 structure at 35.8(5) GPa. Electronic structure and phonon dispersion calculations on MoBi 2 revealed a correlation between valence electron count and bonding in high-pressure transition metal–Bi structures as well as identified two dynamically stable ambient pressure polymorphs. Our study demonstrates the power of the combined computational–experimental approach in capturing highpressure reactivity for efficient materials discovery.

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

Low-dose high-resolution $\mathrm{TOF}$-$\mathrm{PET}$ using ionization-activated multi-state low-$\mathrm{Z}$ detector media

We propose PET scanners using low atomic number media that undergo a persistent local change of state along the paths of the Compton recoil electrons. Measurement of the individual scattering locations and angles, deposited energies, and recoil electron directions allows using the kinematical constraints of the 2-body Compton scattering process to perform a statistical time-ordering of the scatterings, with a high probability of precisely identifying where the gamma first interacted in the detector. In these cases the Line-of-Response is measured with high resolution, determined by the underlying physics processes and not the detector segmentation. There are multiple such media that act through different mechanisms. As an example in which the change of state is quantum-mechanical through a change in molecular configuration, rather than thermodynamic, as in a bubble chamber, we present simulations of a two-state photoswitchable organic dye, a ‘Switchillator’, that is activated to a fluorescent-capable state by the ionization of the recoil electrons. The activated state is persistent, and can be optically excited multiple times to image individual activated molecules. Energy resolution is provided by counting the activated molecules. Location along the LOR is implemented by large-area time-of-flight MCP-PMT photodetectors with single photon time resolution in the tens of ps and sub-mm spatial resolution. Simulations indicate a large reduction of dose.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Preferential Direction of Electron Transfers at a Dye–Metal Oxide Interface with an Insulating Fluorinated Self-Assembled Monolayer and MgO

Slowing non-desirable electron transfer reactions at metal oxide-dye interfaces is important for many technologies. Specifically, after an interfacial photoinduced charge separation event at a metal oxide-dye interface, it is critically important to limit the rate of electron transfer reactions back to sensitizers and to limit electron transfer reactions between the electrolyte and the metal oxide. Ruthenium-based dyes at metal oxide interfaces are widely used in many fields; however, these dyes often have poor surface insulation resulting in fast recombination kinetics with transition metal-based redox shuttles in an electrolyte. This work explores two semiconductor surface modification strategies designed to minimize recombination events of electrons in TiO 2 with oxidized redox shuttles using a fluorinated siloxane insulator (PFTS) and a metal oxide insulator (MgO) with a well-known Ru dye, B11. Additionally, the influence of these treatments on the rate and duration of photoinduced interfacial charge separation at the TiO 2 -dye interface was examined. The TiO 2 -dye-redox shuttle systems were studied via dye-sensitized solar cell current-voltage curve, incident photon-to-current conversion efficiency, small-modulated photovoltage transient, time-correlated single photon counting, and transient absorption spectroscopy measurements. MgO was found to decrease the rate of the electron transfer reaction from the metal oxide to the electrolyte, decrease the rate of the electron transfer reaction to the oxidized dye from TiO 2 , increase the electron transfer reaction rate from a reduced redox shuttle to an oxidized dye, and decrease the electron injection rate from the photoexcited dye to TiO2. Interestingly, PFTS was found to desirably improve on these rates relative to MgO or untreated TiO 2 . Here, a model based on electrostatic interactions is presented to explain the exceptional behavior of PFTS with DFT computational analysis of PFTS supporting this model.

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

Role of chemical disorder on radiation-induced defect production and damage evolution in NiFeCoCr

Understanding chemical disorder in many concentrated solid solution alloys (CSAs) at the levels of electrons and atoms has attracted increasing attention as a path forward to reveal and identify underlying mechanisms for extraordinary mechanical properties and improved radiation tolerance. Single-phase NiFeCoCr CSA is a common base for many high-entropy alloys (HEAs) that have shown improved mechanical strength and radiation tolerance. In this study, defect production and damage evolution in NiFeCoCr under ion irradiation at room temperature to dose over 20 dpa are determined using ion channeling technique along both <100> and <110> directions utilizing multiple probing beam energies. The results obtained from the multi-axial and multi-energy channeling analysis are compared with those previously obtained for Ni crystals irradiated under similar conditions. The influence of chemical complexity on defect production and clustering at early-stage under room temperature irradiation up to dose of 1 dpa is discussed based on positron annihilation spectroscopy results. Defect structure evaluation in Ni and NiFeCoCr is also discussed based on transmission electron microscopy results over a prolonged irradiation at both room and elevated temperatures. Compared with chemically complex NiFeCoCr, larger dislocation loops thus less lattice strain are expected to form in pure Ni. Moreover, the role of chemical disorder in this CSA is also investigated based on ab initio calculations using large supercells. Finally, to understand the impact of chemical complexity effect on defect structure evolution, this integrated research effort attempts to link the relatively large charge redistribution due to difference in valence electron counts resulting from alloying different 3d transition metal elements, moderate lattice distortion arising from similar adaptable atomic size, and notable suppressed or delayed damage evolution in NiFeCoCr.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phonon spectrum of underdoped HgBa 2 CuO 4+δ investigated by neutron scattering

The cuprates exhibit a prominent charge-density-wave (CDW) instability with a wave vector along [100], i.e., the Cu-O bond direction. Whereas CDW order is most prominent at moderate doping and low temperature, there exists increasing evidence for dynamic charge correlations throughout a large portion of the temperature-doping phase diagram. In particular, the signatures of incipient charge order have been observed as phonon softening and/or broadening near the CDW wave vector approximately halfway through the Brillouin zone. Most of this work is focused on moderately doped cuprates, for which the CDW order is robust, or on optimally doped samples, for which the superconducting transition temperature (T c ) attains its maximum. Here we present a time-of-flight neutron scattering study of phonons in simple-tetragonal HgBa 2 CuO 4+δ (T c = 55 K) at a low doping level where prior work showed the CDW order to be weak. We employ and showcase a new software-based technique that mines a large number of measured Brillouin zones for useful data in order to improve accuracy and counting statistics. Density-functional theory has not provided an accurate description of phonons in HgBa2CuO4+δ, yet we find the right set of parameters to qualitatively reproduce the data. The notable exception is a dispersion minimum in the longitudinal Cu-O bond-stretching branch along [100]. Here, this discrepancy suggests that, while CDW order is weak, there exist significant dynamic charge correlations in the optic phonon range at low doping, near the edge of the superconducting dome.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Space-Time Block Preconditioning for Incompressible Flow

Parallel-in-time methods have become increasingly popular in the simulation of time-dependent numerical PDEs, allowing for the efficient use of additional message passing interface processes when spatial parallelism saturates. Most methods treat the solution and parallelism in space and time separately. In contrast, all-at-once methods solve the full space-time system directly, largely treating time as simply another spatial dimension. All-at-once methods offer a number of benefits over separate treatment of space and time, most notably significantly increased parallelism and faster time to solution (when applicable). However, the development of fast, scalable all-at-once methods has largely been limited to time-dependent (advection-)diffusion problems. This paper introduces the concept of space-time block preconditioning for the all-at-once solution of incompressible flow. By extending well-known concepts of spatial block preconditioning to the space-time setting, we develop a block preconditioner whose application requires the solution of a space-time (advection-)diffusion equation in the velocity block, coupled with a pressure Schur complement approximation consisting of independent spatial solves at each time-step, and a space-time matrix-vector multiplication. The new method is tested on four classical models in incompressible flow. Finally, the results indicate perfect scalability in refinement of spatial and temporal mesh spacing, perfect scalability in nonlinear Picard iteration count when applied to a nonlinear Navier--Stokes problem, and minimal overhead in terms of number of preconditioner applications compared with sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Topical Report – Findings on Subtask 2.7 – Wet ESP and Aerosol Testing at Coal Creek Station

Growing concerns over the impact of CO 2 emissions from combustion sources on global climate change have prompted numerous research and development projects aimed at developing cost-effective technologies for CO 2 capture. One family of technologies being demonstrated at pilot and full scale globally is postcombustion carbon capture (PCCC) systems that employ amine-based solvents. The captured CO 2 can be compressed and permanently stored underground or used for enhanced oil recovery. The proximity of North Dakota’s lignite-fired fleet of power plants to potential CO 2 storage options creates a unique atmosphere for PCCC within the state. However, the unique components present in lignite flue gas present a challenge for large-scale PCCC at North Dakota power plants by contributing to aerosol formation. Aerosols can negatively impact the long-term performance of amine-based solvents for CO 2 capture. Amine-based solvents are volatile, and flue gas particulate provides nucleation sites where amine vapors can condense as aerosols. Because aerosols cannot be easily captured at the column outlet using conventional technologies, the amine-laden aerosols escape the system and lead to accelerated solvent losses. Moreover, particulate components can chemically react with amines to form degradation products that can permanently deactivate the amine, cause fouling, and lead to hazardous emissions. Many of the elements that have been shown to catalyze solvent degradation are present in lignite coals and can exacerbate solvent replacement economics. Understanding this issue is critical to the implementation of solvent-based CO 2 capture systems as applied to lignite-fired generation systems. The Energy & Environmental Research Center (EERC) designed and carried out this project to fully characterize aerosol behavior with various control technologies installed to better optimize aerosol mitigation technology for CO 2 capture. To meet the goal of this project, the following objectives were identified: Determine the effectiveness of a wet electrostatic precipitator (WESP) on mitigating formation of problematic aerosols at Great River Energy’s Coal Creek Station, upstream of the PCCC system. Determine the effectiveness of the Mitsubishi Heavy Industries (MHI) proprietary amine emission reduction unit (AERU) as a postcapture solvent recovery system for reducing aerosol emissions and extending solvent life downstream of the PCCC system. Determine the impact of aerosols on the efficiency and degradation products of both commercial and advanced solvents within the PCCC system. Work was conducted at Coal Creek Station Unit 1 using a slipstream of flue gas from the outlet of the plant’s flue gas desulfurization (FGD) unit. Flue gas was routed through a pilot-scale FGD unit to remove SO 2 to very low levels (~1 ppm) and then through a direct contact cooler (DCC) to further cool the gas and to remove moisture. The gas exiting the DCC was then optionally routed through a WESP before passing to the CO 2 absorber columns. The MHI solvent was used to scrub CO 2 from the slipstream through a set of two absorber columns. The rich solvent was regenerated in a stripper column by heating to drive off captured CO 2 . Flue gas exiting the absorber column was routed to MHI’s proprietary AERU to recover entrained solvent. The system operated using a catch-and-release method where the CO 2 was separated to provide data on the process, but the captured CO 2 was released back into the host site stack. Particulate was measured, collected, and analyzed from multiple locations throughout the pilot-scale system. Unlike the performance observed in prior work, the inlet FGD and DCC did not remove significant particulate matter from the flue gas. This appears to be due to a difference in the nature of the particulate. The DCC seemed to increase particulate size and count, most likely owing to water condensing onto the surfaces of fly ash particles. When the WESP was operated, it achieved >95% particulate capture. Very little particulate matter or indications of solvent were detected at the AERU outlet. When operating with advanced KS 21 solvent, the particulate material at the AERU outlet was even further decreased. Solvent analysis showed that some species derived from flue gas and ash were slowly concentrating in the solvent over the duration of the test. The levels observed were reported to be within expected ranges and were not of concern to MHI. A high-level techno-economic assessment of installing CO 2 capture at Coal Creek Station suggested that, when using a standard monoethanolamine (MEA)-based solution with simple heat integration, the energy penalty to net generation would be 34%. The bulk of this was due to steam losses for regenerating solvent, followed by parasitic electrical demand for CO 2 compression and then by increased parasitic load for pushing flue gas through the absorber column. These demands could be decreased with a more advanced solvent that exhibits lower heat of regeneration and lower pressure drop than does a simple MEA solution. Further energy could be saved with more thorough heat integration to recover useful energy from the steam used for solvent regeneration. Installing a WESP was predicted to increase the cost of electricity by nearly $5/MWh. This would become cost-effective if solvent losses were roughly 10 times the baseline estimate when not using a WESP but could be returned to baseline by installing the WESP. Piping CO 2 for storage in more favorable geology could help with carbon capture and storage economics. Although storing off-site would necessitate construction of a CO 2 transport pipeline, the cost of this pipeline might be more than offset by reducing the number of wells required, the depths of the wells required, and the electrical demand for the CO 2 compressor. Additional factors that favor off-site storage costs include smaller expected CO 2 plume sizes, which translates to less monitoring and fewer landowner agreements. More detailed assessment of the specific geology in the region under and around Coal Creek Station would be needed to accurately assess the costs and benefits of different storage site options. This subtask was cofunded through the Energy & Environmental Research Center–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission.

Strege, Joshua↗
DOE OSTI2024

Daily Arctic Lightning Strokes

In recent decades, lightning activity at high latitudes has increased. Tall thunderstorm clouds affect radiation balance directly, as well as indirectly through lightning-initiated fires and the resulting smoke. One can remotely sense lightning strokes over the globe through their VLF radio emission, and, with multiple receivers, it is possible to precisely locate lightning strokes. This technique makes it possible to continuously monitor arctic lightning--a capability not possible by other means. In this research effort, a World-Wide Lightning Location Network station, consisting of a VLF receiver, signal processing hardware, and analysis software, were installed at the North Slope of Alaska (NSA) facility and planned to operate for several years. This far north location is expected to improve the network's high-latitude detection efficiency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Metal Bioavailability and Ecotoxicity of Bioremediated Oils and Tailings by BioTiger{sup TM}, a Microbial Consortium

Oil Sands and Mature Fine Tailings: Oil sand reserves are a major source of oil for the United States. Oil sands are a mixture of sand, clay, water, and bitumen. The refining process requires large volumes of water and generates hazardous Mature Fine Tailings (MFTs) that are stored in engineered settling ponds. They are of major environmental concern due to their persistence and difficulty to naturally biodegrade. MFT Contaminants of Concern: Naphthenic Acids (NAs), Polycyclic Aromatic Hydrocarbons (PAHs), Benzene, Toluene, Ethylbenzene and Xylene (BTEX), Metals, Residual Bitumen. BioTiger{sup TM} (BT{sup TM}), the SRNL patented 12 component microbial consortia was found to cometabolically degrade some NAs and PAHs. Previous studies performing short term exposure (48 hours and 7 days) to BT{sup TM} have resulted in increased toxicity due to the partial degradation of PAHs forming toxic intermediates. Objective: Evaluate the ecotoxicity of BT{sup TM} remediated MFTs exposed to a 2 week period w/ yeast (Y) extract. This work will determine BT{sup TM}'s remediation of MFTs from Fort McMurray, Alberta. Acute toxicity tests will be performed under section 9 of EPA's Method for Measuring Acute Toxicity using the freshwater organism, Ceriodaphnia dubia. Monitor BT{sup TM} growth through Most Probable Number (MPN) counts, pH, and metal bioavailability will also be evaluated. Preparing tailing solutions for biodegradation: BT{sup TM} components were grown in R2A media and combined. Mother BT{sup TM} was centrifuged at 7000 rpm for 20 minutes and resuspended in Bushnell Haus. 33 g of tailings were added to 1100 mL Bushnell Haus along with yeast extract (11 g) to create a 3% tailing solution. Treatments were performed in triplicate. Preparing solutions for toxicity tests using C. dubia: Supernatant was collected immediately after centrifuging at 7000 rpm for 20 min. Supernatant was filtered via a 0.22 μm sterile system. All treatments except MFT, Y peaked at T=4. This decrease in microbial growth could be associated with toxicity from intermediate byproducts. Treatments containing only MFTs decreased in the initial 11 days before growing significantly. Final solution for MFT, Y treatment had the highest pH. This could potentially be due to increased microbial activity linked with yeast consumption. All treatments except MFT, Y peaked at T=4. This decrease in microbial growth could be associated with toxicity from intermediate byproducts. Treatments containing only MFTs decreased in the initial 11 days before growing significantly. Future Direction: Await results from Acute Toxicity Tests and Metal Bioavailability Quantification of biosurfactant production. Evaluate hydrocarbon degradation byproducts using an analytical chemistry approach. More toxicity experiments with variations in time and conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

An Update on Cardinal: Toward Full Core Pebble Simulations

The new DOE-NE Center of Excellence for Thermal-fluids applications in Nuclear Energy inaugurated in April 2018 considers and researches novel new solution strategies for historically challenging flow issues that still plague the current fleet of deployed Light Water nuclear Reactors (LWRs) as well as predicting various fluid flow and fluid related issues with advanced reactor technologies. Our advanced thermal-fluids research and development approach synergistically combines three natural, though overlapping, length and time scales in a hierarchal multi-scale approach to avoid the temptation and pitfalls of attempting to develop a single solve-all algorithm for physical fluid flow problems that will span nine orders of magnitude in spatial and temporal scales. To demonstrate the multi-scale philosophy of the center we focus on Fluoride Cooled High Temperature Reactors (FHRs), and in particular on the Berkley’s PB-FHR Mark-I design. The Fluoride salt cooled High temperature Reactor (FHR) is a class of advanced nuclear reactors that combine the robust coated particle fuel form from high temperature gas cooled reactors, direct reactor auxiliary cooling system (DRACS) passive decay removal of liquid metal fast reactors, and the transparent, high volumetric heat capacitance liquid Fluoride salt working fluids - such as Flibe - from molten salt reactors. This combination of fuel and coolant enables FHRs to operate in a high-temperature low-pressure design space that has beneficial safety and economic implications. The PB-FHR reactor relies on a pebble bed approach and pebble bed reactors are in a sense the poster child for multi-scale analysis. The lower length-scale simulator for pebble reactor cores comprises three physics: neutronics (OpenMC), thermal-fluids (Nek5000/NekRS) and fuel performance (BISON). As part of Center of Excellence ongoing research efforts we have developed Cardinal, a new tool platform for lower length-scale simulation. Cardinal tightly couples all three physics and leverages advances in MOOSE such as the MultiApp system and the concept of MOOSE-wrapped Apps. The present report provides an update on the development of Cardinal with an extension of Cardinal to GPUs. We also perform a first-of-a-kind demonstration simulation on Summit representing a 10x capability increase in terms of pebble count for Cardinal.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

The role of quantum computing in advancing scientific high-performance computing: A perspective from the ADAC institute

Quantum computing (QC) has gained significant attention over the past two decades due to its potential for speeding up classically demanding tasks. This transition from an academic focus to a thriving commercial sector is reflected in substantial global investments. While advancements in qubit counts and functionalities continue at a rapid pace, current quantum systems still lack the scalability for practical applications, facing challenges such as too high error rates and limited coherence times. Here, this perspective paper examines the relationship between QC and high-performance computing (HPC), highlighting their complementary roles in enhancing computational efficiency. It is widely acknowledged that even fully error-corrected QC will not be suited for all computational tasks. Rather, future compute infrastructures are anticipated to employ quantum acceleration within hybrid systems that integrate HPC and QC. While QC can enhance classical computing, traditional HPC remains essential for maximizing quantum acceleration. This integration is a priority for supercomputing centers and companies, sparking innovation to address the challenges of merging these technologies. The novelty of this work lies in its unique perspective, reflecting the collective insights of the Accelerated Data Analytics and Computing (ADAC) Institute, a global consortium of over 20 leading HPC centers. Recognizing the growing importance of QC, ADAC established a Quantum Computing Working Group in 2023 to foster collaboration and knowledge-sharing among its members. This paper synthesizes insights from the group’s collaborative efforts and incorporates findings from a member survey that captures shared experiences, ongoing projects, and strategic directions. By outlining the current landscape and challenges of QC integration into HPC ecosystems, this work offers HPC specialists practical and forward-looking guidance on the opportunities and implications of QC in computationally intensive endeavors.

Accelerated Data Analytics and↗