SEARCH · Engineering Papers

Results for “time integrator”

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.

SourceResearch connectionFrom yearThrough year

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

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Search for Astrophysical Neutrinos from 1FLE Blazars with IceCube

The majority of astrophysical neutrinos have undetermined origins. The IceCube Neutrino Observatory has observed astrophysical neutrinos but has not yet identified their sources. Blazars are promising source candidates, but previous searches for neutrino emission from populations of blazars detected in ≳GeV gamma rays have not observed any significant neutrino excess. Recent findings in multimessenger astronomy indicate that high-energy photons, coproduced with high-energy neutrinos, are likely to be absorbed and reemitted at lower energies. Thus, lower-energy photons may be better indicators of TeV–PeV neutrino production. This paper presents the first time-integrated stacking search for astrophysical neutrino emission from MeV-detected blazars in the first Fermi Large Area Telescope low energy (1FLE) catalog using ten years of IceCube muon–neutrino data. The results of this analysis are found to be consistent with a background-only hypothesis. Assuming an E -2 neutrino spectrum and proportionality between the blazars MeV gamma-ray fluxes and TeV–PeV neutrino flux, the upper limit on the 1FLE blazar energy-scaled neutrino flux is determined to be 1.64 x 10 -12 TeV cm -2 s -1 at 90% confidence level. This upper limit is approximately 1% of IceCube’s diffuse muon–neutrino flux measurement.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

IceCube Search for Neutrinos Coincident with Gravitational Wave Events from LIGO/Virgo Run O3

Using data from the IceCube Neutrino Observatory, we searched for high-energy neutrino emission from the gravitational-wave events detected by the advanced LIGO and Virgo detectors during their third observing run. We did a low-latency follow-up on the public candidate events released during the detectors' third observing run and an archival search on the 80 confident events reported in the GWTC-2.1 and GWTC-3 catalogs. An extended search was also conducted for neutrino emission on longer timescales from neutron star containing mergers. Follow-up searches on the candidate optical counterpart of GW190521 were also conducted. We used two methods; an unbinned maximum likelihood analysis and a Bayesian analysis using astrophysical priors, both of which were previously used to search for high-energy neutrino emission from gravitational-wave events. No significant neutrino emission was observed by any analysis, and upper limits were placed on the time-integrated neutrino flux as well as the total isotropic equivalent energy emitted in high-energy neutrinos.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

The Atacama Cosmology Telescope: Systematic Transient Search of Single Observation Maps

We conduct a systematic search for astrophysical transients using data from the Atacama Cosmology Telescope. The data were taken from 2017 to 2022 in three frequency bands spanning 77 to 277 GHz. In this paper, we present a pipeline for transient detection using single-observation maps where each pixel of a map contains one observation with an integration time of approximately 4 minutes. We detect 34 transient events at 27 unique locations. All but two of the transients are associated with Galactic stars and exhibit a wide range of properties. We also detect an event coincident with the classical nova YZ Ret and one event consistent with a flaring active galactic nucleus. We notably do not detect any reverse shock emission from gamma-ray bursts, a nondetection that may be in tension with current models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Search for High-energy Neutrino Emission from Galactic X-Ray Binaries with IceCube

We present the first comprehensive search for high-energy neutrino emission from high- and low-mass X-ray binaries conducted by IceCube. Galactic X-ray binaries are long-standing candidates for the source of Galactic hadronic cosmic rays and neutrinos. The compact object in these systems can be the site of cosmic-ray acceleration, and neutrinos can be produced by interactions of cosmic rays with radiation or gas, in the jet of a microquasar, in the stellar wind, or in the atmosphere of the companion star. We study X-ray binaries using 7.5 yr of IceCube data with three separate analyses. In the first, we search for periodic neutrino emission from 55 binaries in the Northern Sky with known orbital periods. In the second, the X-ray light curves of 102 binaries across the entire sky are used as templates to search for time-dependent neutrino emission. Finally, we search for time-integrated emission of neutrinos for a list of 4 notable binaries identified as microquasars. In the absence of a significant excess, we place upper limits on the neutrino flux for each hypothesis and compare our results with theoretical predictions for several binaries. In addition, we evaluate the sensitivity of the next generation neutrino telescope at the South Pole, IceCube-Gen2, and demonstrate its power to identify potential neutrino emission from these binary sources in the Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Kinematic flow from the flow of cuts

The wavefunction coefficients of conformally coupled scalars in power-law FRW cosmologies satisfy differential equations governed by a set of simple combinatorial rules known as the kinematic flow. In this paper we derive the kinematic flow, expressed using a set of differential forms referred to as the cut basis, from a geometric perspective, relying solely on the cosmological hyperplane arrangement and without invoking bulk physics. Each element of the cut basis corresponds to the positive geometry associated to an independent cut of the physical FRW-form and can be labeled by decorating (minors of) the truncated Feynman graph with an acyclic orientation. We provide a straightforward prescription to associate a logarithmic differential form to each element of the cut basis by considering its corresponding decorated graph. Moreover, we show that the residues of the physical FRW-form are canonical forms of certain graphical zonotopes labeled by the same set of decorated graphs. These zonotopes control the cut combinatorics -- flow of cuts -- of the physical FRW-form and the cut basis (by construction). Using the theory of relative twisted cohomology and intersection theory, we derive a closed form formula for the differential equations of the cut basis. We also introduce combinatorial rules that compute the kinematic differential of any basis element without explicit calculation. The combinatorics of our differential equations is a natural consequence of the flow of cuts and is equivalent (up to rescaling) to the kinematic flow for the recently studied time integral basis. In particular, our differential equations decouple into exponentially many sectors, one for each way of cutting a subset of edges of the graph.

General Relativity and Quantum Cosmology↗

DOE OSTI2026

Chemical ionization mass spectrometry utilizing benzene cations for measurements of volatile organic compounds and nitric oxide

We evaluate the capability of chemical ionization mass spectrometry (CIMS) using benzene cations as reagent ions (benzene CIMS) for detecting atmospheric trace gases. We characterize the ionization pathways and product ion distributions for 27 analytes spanning diverse chemical classes. To interpret the complex ion chemistry involving two reagent ions (C 6 H$^{+}_{6}$ and (C 6 H 6 )$^{+}_{2}$) and multiple ionization pathways (charge transfer, proton transfer, adduct formation, and hydride abstraction), we introduce a thermodynamics-based framework that classifies analytes into three categories based on their ionization energy (IE), relative to those of benzene monomer (9.24 eV) and dimer (8.69 eV). Each class exhibits distinct ionization mechanisms and product ions. Analytes with IE smaller than 8.69 eV (low IE) undergo charge transfer with both reagent ions; analytes with IE between 8.69 and 9.24 eV (mid IE) undergo charge transfer with C 6 H$^{+}_{6}$ and potential adduct formation with (C 6 H 6 )$^{+}_{2}$; analytes with IE larger than 9.24 eV (high IE) could undergo adduct formation, proton transfer, or hydride abstraction. Analytes within each class also show similar sensitivity, enabling sensitivity estimation for compounds lacking calibration standards. In addition to volatile organic compounds (VOCs), benzene CIMS detects nitric oxide (NO) with a detection limit of 5 pptv for 1 min integration time, exceeding the performance of most commercial NOx analyzers. Field deployments in Chicago and St. Louis demonstrate good agreement with reference NO measurements. Isoprene measurements show good agreement with a co-located gas chromatography–photoionization detector (GC-PID) in St. Louis, but exhibit substantial positive bias in Chicago, likely due to interferences from anthropogenic VOCs in the polluted urban environment. These results highlight the potential of benzene CIMS for concurrent measurements of NO, VOCs, and their oxidation products using a single instrument, while also underscoring challenges in complex atmospheric conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Exploring the impacts of unprecedented climate extremes on forest ecosystems: hypotheses to guide modeling and experimental studies

Abstract. Climatic extreme events are expected to occur more frequently in the future, increasing the likelihood of unprecedented climate extremes (UCEs) or record-breaking events. UCEs, such as extreme heatwaves and droughts, substantially affect ecosystem stability and carbon cycling by increasing plant mortality and delaying ecosystem recovery. Quantitative knowledge of such effects is limited due to the paucity of experiments focusing on extreme climatic events beyond the range of historical experience. Here, we present a road map of how dynamic vegetation demographic models (VDMs) can be used to investigate hypotheses surrounding ecosystem responses to one type of UCE: unprecedented droughts. As a result of nonlinear ecosystem responses to UCEs that are qualitatively different from responses to milder extremes, we consider both biomass loss and recovery rates over time by reporting a time-integrated carbon loss as a result of UCE, relative to the absence of drought. Additionally, we explore how unprecedented droughts in combination with increasing atmospheric CO2 and/or temperature may affect ecosystem stability and carbon cycling. We explored these questions using simulations of pre-drought and post-drought conditions at well-studied forest sites using well-tested models (ED2 and LPJ-GUESS). The severity and patterns of biomass losses differed substantially between models. For example, biomass loss could be sensitive to either drought duration or drought intensity depending on the model approach. This is due to the models having different, but also plausible, representations of processes and interactions, highlighting the complicated variability of UCE impacts that still need to be narrowed down in models. Elevated atmospheric CO2 concentrations (eCO2) alone did not completely buffer the ecosystems from carbon losses during UCEs in the majority of our simulations. Our findings highlight the consequences of differences in process formulations and uncertainties in models, most notably related to availability in plant carbohydrate storage and the diversity of plant hydraulic schemes, in projecting potential ecosystem responses to UCEs. We provide a summary of the current state and role of many model processes that give way to different underlying hypotheses of plant responses to UCEs, reflecting knowledge gaps which in future studies could be tested with targeted field experiments and an iterative modeling–experimental conceptual framework.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Further improvement and evaluation of nudging in the E3SM Atmosphere Model version 1 (EAMv1): simulations of the mean climate, weather events, and anthropogenic aerosol effects

Abstract. A previous study on the use of nudging in E3SM Atmosphere Model version 1 (EAMv1) had an unresolved issue; i.e., a simulation nudged to EAMv1's own meteorology showed non-negligible deviations from the free-running baseline simulation over some of the subtropical marine stratocumulus and trade cumulus regions. Here, we demonstrate that the deviations can be substantially reduced by (1) changing where the nudging tendency is calculated in the time integration loop of a nudged EAM simulation so as to improve consistency with the free-running baseline and (2) increasing the frequency of the constraining data so as to better capture strong sub-diurnal variations. The fact that modification (2) improves the climate representativeness of the nudged simulations has motivated us to investigate whether the use of newer reanalysis products with higher data frequency can help improve nudged hindcast simulations by better capturing the observed weather events. To answer this question, we present simulations conducted at EAMv1's standard horizontal resolution (approximately 1∘) with nudging towards 6-hourly ERA-Interim reanalysis or 6-hourly, 3-hourly, or hourly ERA5 reanalysis. These simulations are evaluated against the climatology of free-running EAMv1 simulations as well as reanalyses, satellite retrievals, and in situ measurements from the Atmospheric Radiation Measurement user facility. For the 1∘ EAMv1 simulations, we recommend using the relocated nudging tendency calculation and the ERA5 reanalysis at 3-hourly or higher frequency. Simulations used for estimating the anthropogenic aerosol effects often use nudging to help discern signal from noise. The sensitivity of such estimates to the configuration of nudging is investigated in EAMv1, again using the standard 1∘ horizontal resolution. We find that, when estimating the global mean effects, the frequency of constraining data has relatively small impacts, while the choice of nudged variables can change the results substantially. The nudging of air temperature (in addition to horizontal winds) has two non-negligible effects. First, when the constraining data come from reanalysis, the nudging-induced mean bias correction can cause significant changes in the simulated clouds and hence substantially different estimates of the aerosol effects. The impact of the mean bias correction on ice cloud formation has been noted in previous studies and is also seen in EAMv1. For applications like ours, where the preferred configurations of nudging are those capable of providing results consistent with the multi-year free-running simulations, the consequence of the mean bias correction is undesirable. The second important impact of temperature nudging is a significant suppression of adjustments to aerosol forcing, which also causes changes in the estimated aerosol effects. This effect can be seen in simulations nudged to either reanalysis or EAM's own meteorology. These results suggest that nudging horizontal winds but not temperature is a better choice for estimating the anthropogenic aerosol effects.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Numerical coupling of aerosol emissions, dry removal, and turbulent mixing in the E3SM Atmosphere Model version 1 (EAMv1) – Part 1: Dust budget analyses and the impacts of a revised coupling scheme

Abstract. An earlier study evaluating dust life cycle in the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1) has revealed that the simulated global mean dust lifetime is substantially shorter when higher vertical resolution is used, primarily due to significant strengthening of dust dry removal in source regions. This paper demonstrates that the sequential splitting of aerosol emissions, dry removal, and turbulent mixing in the model's time integration loop, especially the calculation of dry removal after surface emissions and before turbulent mixing, is the primary reason for the vertical resolution sensitivity reported in that earlier study. Based on this reasoning, we propose a revised numerical process coupling scheme that requires the least amount of code changes, in which the surface emissions are applied before turbulent mixing instead of before dry removal. The revised scheme allows newly emitted particles to be transported aloft by turbulence before being removed from the atmosphere, and hence better resembles the dust life cycle in the real world. Sensitivity experiments show that the revised process coupling substantially weakens dry removal and strengthens vertical mixing in dust source regions. It also strengthens the large-scale transport from source to non-source regions, strengthens dry removal outside the source regions, and strengthens wet removal and activation globally. In transient simulations of the years 2000–2009 conducted using 1∘ horizontal grid spacing, 72 vertical layers, and unchanged tuning parameters of emission strength, the revised process coupling leads to a 40 % increase in the global total dust burden and an increase of dust lifetime from 1.8 to 2.5 d in terms of 10-year averages. Weakened dry removal and increased mixing ratios are also seen for other aerosol species that have substantial surface emissions, although the changes in mixing ratio are considerably smaller for the submicron species than for dust and sea salt. Numerical experiments confirm that the revised coupling scheme significantly reduces the strong and non-physical sensitivities of model results to vertical resolution in the original EAMv1. This provides a motivation for adopting the revised scheme in EAM as well as for further improvements on the simple revision presented in this paper.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

MAESTROeX

MAESTROeX is a massively parallel, finite-volume C++/F90 solver for low Mach number astrophysical flows. The code utilizes a low Mach number equation set allowing for more efficient, long-time integration of highly subsonic flows compared to compressible approaches. The recommended range of applicability is for flows where the Mach number does not exceed~0.1.

Harpole, Alice↗

DOE OSTI2022

CabanaPD

CabanaPD is a meshfree peridynamics application built with Cabana and Kokkos. Kokkos enables performance portability across hardware architectures and Cabana provides particle capabilities including multi-node MPI support. The main components of CabanaPD are particle initialization, neighbor list generation, force and energy computation, time integration, and multi-node particle communication. In addition, options for creating pre-cracked regions and particle boundary conditions are available. CabanaPD currently enables two common force models: prototype microelastic brittle (PMB) and linear peridynamic solid (LPS). For both of these models, versions with and without fracture as well as linearized model options are available. CabanaPD is designed to be extensible for addition of more complex force models, boundary conditions, etc.

Reeve, Sam↗

DOE OSTI2021

Searching for Dark Matter with Superconducting Qubits

Dark matter is an enduring mystery in our quest to understand the fundamental constituents of our universe. Low mass bosons, such the axion or hidden photon, are compelling dark matter candidates. We leverage their potential interactions with electromagnetic fields, whereby the dark matter (of unknown mass) on rare occasion converts into a single photon, to devise a method of detecting these candidates. Current dark matter searches operating at microwave frequencies use a resonant cavity to coherently accumulate the field sourced by the dark matter and a near standard quantum limited (SQL) linear amplifier to read out the cavity signal. To further increase sensitivity to the dark matter signal and enable future searches, sub-SQL detection techniques are required.In this thesis, I report the development of a novel microwave photon counting technique and a new exclusion limit on hidden photon dark matter. We operate a superconducting qubit to make repeated quantum non-demolition measurements of cavity photons and apply a hidden Markov model analysis to reduce the noise to 15.7 dB below the quantum limit, with overall detector performance limited by a residual background of real photons. With the present device, we perform a hidden photon search and constrain the kinetic mixing angle to less than 1.68 x 10-15 in a band around 6.011 GHz (24.86 micro-eV) with an integration time of 8.33 s. This demonstrated noise reduction technique enables future dark matter searches to be sped up by a factor of 1300. By coupling a qubit to an arbitrary quantum sensor, more general sub-SQL metrology is possible with the techniques presented in this work.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Simulation of Reverse Osmosis Membrane Compaction Using Material Point Method (MPM)

Access to fresh drinking water in the future requires effective management of industrial wastewater and water purification from available resources. In this study, we present the simulation methodology and the analysis of a Reverse Osmosis (RO) membrane under various pressure conditions using the Material Point Method (MPM). In contrast to other numerical methods, MPM solves the continuum governing equations on material points in a Lagrangian framework. The method does not require a grid connecting the material points hence making it suitable to simulate large deformations during membrane compaction. The time integration is carried out using the explicit Euler method, while the spatial discretization is performed using linear or cubic-spline shape functions. A series of planar images containing detailed pore structures obtained from X-ray tomography experiments is converted to a three-dimensional collection of material points to simulate the membrane. Compressive loads are applied to the top layer of the membrane to simulate the experiments. The membrane deformation and pore size distribution before and after load application are reported and compared with the experimental measurements. The presentation discusses the numerical methods used, the performance of the solver on high-performance computing machines, and the results of membrane compaction in detail.

ENGINEERING,MATHEMATICS AND COMPUTING↗

DOE OSTI2023

PI Loop Resonance Control for Dark Photon Experiment at 2 K Using a 2.6 GHz SRF Cavity

Two 2.6 GHz SRF cavities are being used for a dark photon search at the vertical test stand (VTS) in FNAL, for the second phase of the Dark SRF experiment. During testing at 2 K the cavities experience frequency detuning caused by microphonics and slow frequency drifts. The experiment requires that the two cavities have the same frequency within the cavity's bandwidth. These two cavities are equipped with frequency tuners consisting of three piezo actuators. The piezo actuators are used for fine-fast frequency tuning. A proportional-integral (PI) loop utilizing the three piezos on the emitter was used to stabilize the cavity frequency and match the receiver cavity frequency. The results from this implementation will be discussed. The integration time was also calculated via simulation.

43 PARTICLE ACCELERATORS↗

DOE OSTI2022

High-order accurate solutions of the point kinetics equations with the spectral deferred correction method

Solving initial value problems with higher-order methods can improve the accuracy of the simulation results or the efficiency of the calculation. In this paper, we apply the spectral deferred correction (SDC) method to solve the initial value problem of the point kinetics equations (PKE). SDC is a stable, robust, and efficient high-order time-integration method capable of an arbitrary order of accuracy. For our implementation we show that it is A-stable for orders up to 8 and the order of accuracy is verified for PKE problems with a range of different reactivities. A 5.-order SDC method was then implemented to solve the exact PKE (EPKE) in the Transient Multilevel (TML) method of MPACT. The error from solutions of the EPKE is shown to be negligible. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Quantum Enhanced Dark-Matter Search with Entangled Fock States in High-Quality Cavities

We present a quantum-enhanced protocol for detecting wave-like dark matter using an array of $N$ entangled superconducting cavities initialized in an $m$-photon Fock state. By distributing and recollecting the quantum state with an entanglement-distribution operation, the scan rate scales as $N^2(m+1)$ while thermal excitation is the dominant background, significantly outperforming classical single-cavity methods under matched conditions. We evaluate the robustness of our scheme against additional noise sources, including decoherence and beamsplitter infidelity, through theoretical analysis and numerical simulations. In practice, the key requirements, namely high-Q superconducting radio-frequency cavities that support long integration times, high-fidelity microwave beamsplitters, and universal cavity control, are already available on current experimental platforms, making the protocol experimentally feasible.

Freiman, Benjamin [Fermilab]↗

DOE OSTI2023

Time-lapse VSP integration and calibration of subsurface stress field utilizing machine learning approaches: A case study of the morrow B formation, FWU

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse vertical seismic profiling (VSP) integration. The selected study site is a region around the injector well located within Farnsworth field unit (FWU), where there is an ongoing CO 2 -enhanced oil recovery (EOR) operation. In our study, a site-specific rock physics model was created from extensive geological, geophysical, and geomechanical characterization through 3D seismic data, well logs, and core assessed as part of the 1D MEM conducted on the characterization well within the study area. Here, the Biot-Gassmann workflow was utilized to combine the rock physics and reservoir simulation outputs to determine the seismic velocity change due to fluid substitution. Modeled seismic velocities attributed to mean effective stress were determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from ultrasonic seismic velocity measurements. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer (PSO) was utilized to minimize a penalty function created between the modeled seismic velocities and the observed time-lapse VSP dataset. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.

58 GEOSCIENCES↗

DOE OSTI2025

Quantum chaos, integrability, and late times in the Krylov basis

Quantum chaotic systems are conjectured to display a spectrum whose fine-grained features (gaps and correlations) are well described by random matrix theory (RMT). We propose and develop a complementary version of this conjecture: quantum chaotic systems display a Lanczos spectrum whose local means and covariances are well described by RMT. To support this proposal, we first demonstrate its validity in examples of chaotic and integrable systems. We then show that for Haar-random initial states in RMTs the mean and covariance of the Lanczos spectrum suffice to produce the full long-time behavior of general survival probabilities including the spectral form factor, as well as the spread complexity. In addition, for initial states with continuous overlap with energy eigenstates, we analytically find the long-time averages of the probabilities of Krylov basis elements in terms of the mean Lanczos spectrum. This analysis suggests a notion of eigenstate complexity, the statistics of which differentiate integrable systems and classes of quantum chaos. Lastly, we clarify the relation between spread complexity and the universality classes of RMT by exploring various values of the Dyson index and Poisson distributed spectra.

Combinatorics↗