Engineering Papers⌕ Search

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.

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

At least 235 records · Page 13

Emission phases of implosion sources for x-ray absorption fine structure spectroscopy

At the Laboratory for Laser Energetics' Omega Laser Facility, thin plastic shells were directly driven with ∼20 kJ resulting in a time-integrated x-ray yield of ∼1012 ph/eV/sr at 7 keV. Using temporally, spatially, and spectrally discriminating diagnostics, three x-ray emission phases were identified: corona emission produced by the laser ablation of the shell, core stagnation, and afterglow emission due to the expanding hot material after stagnation. The newly measured corona and afterglow emission phases account for ∼25% of the total x-ray signal and produce x-ray emission at a different time or larger radius than previously considered. The resulting implications of this additional emission for x-ray absorption fine structure spectroscopy are discussed. Finally, improvements to the laser drive intensity and uniformity produced a factor-of-2 increase in total x-ray emission while decreasing the size of the stagnated core.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Engineering spin coherence in core-shell diamond nanocrystals

Fluorescent diamond nanocrystals can host spin qubit sensors capable of probing the physical properties of biological systems with nanoscale spatial resolution. Sub-100 nm diamond nanosensors can readily be delivered into intact cells and even living organisms. However, applications beyond current proof-of-principle experiments require a substantial increase in sensitivity, which is limited by surface induced charge instability and electron-spin dephasing. In this work, we utilize engineered core-shell structures to achieve a drastic increase in qubit coherence times (T 2 ) from 1.1 to 35 μs in bare nanodiamonds to upward of 52 to 87 μs. We use electron-paramagnetic-resonance results to present a band bending model and connect silica encapsulation to the removal of deleterious mid-gap surface states that are negatively affecting the qubit’s spin properties. Combined with a 1.9-fold increase in particle luminescence these advances correspond to up to two-order-of-magnitude reduction in integration time. Probing qubit dynamics at a single particle level further reveals that the noise characteristics fundamentally change from a bath with spins that rearrange their spatial configuration during the course of an experiment to a more dilute static bath. The observed results shed light on the underlying mechanisms governing fluorescence and spin properties in diamond nanocrystals and offer an effective noise mitigation strategy based on engineered core-shell structures.

core-shell↗

Optical depth to reionization from perturbative 21 cm clustering

The optical depth τ is the least well determined parameter in the standard model of cosmology, and one whose precise value is important for both understanding reionization and for inferring fundamental physics from cosmological measurements. Here, we forecast how well future epoch of reionization experiments could constraint τ using a symmetries-based bias expansion that highlights the special role played by anisotropies in the power spectrum on large scales. Given a parametric model for the ionization evolution inspired by the physical behavior of more detailed reionization simulations, we find that future 21 cm experiments could place tight constraints on the timing and duration of reionization and hence constraints on τ that are competitive with proposed, space-based CMB missions provided they can measure k ≈ 0.1 h Mpc -1 with a clean foreground wedge across redshifts spanning the most active periods of reionization, corresponding to ionization fractions 0.2 ≲ x ≲ 0.8. Significantly improving upon existing CMB-based measurements with next-generation 21 cm surveys would require substantially longer observations (~ 5 years) than standard $\mathscr{O}$(1000 hour) integration times. Precise measurements of smaller scales will not improve constraints on τ until a better understanding of the astrophysics of reionization is achieved. In the presence of noise and foregrounds even future 21 cm experiments will struggle to constrain τ if the ionization evolution deviates significantly from simple parametric forms.

79 ASTRONOMY AND ASTROPHYSICS↗

IceCube Search for Neutrinos Coincident with Compact Binary Mergers from LIGO-Virgo’s First Gravitational-wave Transient Catalog

Using the IceCube Neutrino Observatory, we search for high-energy neutrino emission coincident with compact binary mergers observed by the LIGO and Virgo gravitational-wave (GW) detectors during their first and second observing runs. In this paper, we present results from two searches targeting emission coincident with the sky localization of each GW event within a 1000 s time window centered around the reported merger time. One search uses a model-independent unbinned maximum-likelihood analysis, which uses neutrino data from IceCube to search for pointlike neutrino sources consistent with the sky localization of GW events. The other uses the Low-Latency Algorithm for Multi-messenger Astrophysics, which incorporates astrophysical priors through a Bayesian framework and includes LIGO-Virgo detector characteristics to determine the association between the GW source and the neutrinos. No significant neutrino coincidence is seen by either search during the first two observing runs of the LIGO-Virgo detectors. We set upper limits on the time-integrated neutrino emission within the 1000 s window for each of the 11 GW events. These limits range from 0.02 to 0.7 GeV cm -2 . We also set limits on the total isotropic equivalent energy, E iso , emitted in high-energy neutrinos by each GW event. These limits range from 1.7 × 10 51 to 1.8 × 10 55 erg. We conclude with an outlook for LIGO-Virgo observing run O3, during which both analyses are running in real time.

79 ASTRONOMY AND ASTROPHYSICS↗

A Python Library for Radiance Matrix-based Simulation Control and EnergyPlus Integration

Radiance matrix-based methods enable efficient parametric simulations, allowing users to vary sky conditions, fenestration systems, and other model parameters at a minimal cost to computation. However, the steep learning curve and complex workflow hinder the widespread adoption of matrix-based methods. The frads Python library with a series of command-line tools was developed to automate the entire matrix-based simulation process, lowering entry barriers and reducing human error. Co-simulation between EnergyPlus and Radiance was also enabled using the Python library from EnergyPlus. Key Innovations • Command-line based automation of Radiance matrix-based simulation methods • Python library facilitates broader adoption of Radiance matrix-based simulation methods • Radiance EnergyPlus run-time integration enabling the modeling of advanced control systems Practical Implications The frads library, with associated command-line tools, provides practitioners with the capability to easily adopt and use Radiance matrix-based simulation methods for various daylighting, solar control, and energy-related evaluations. Frads' current form is designed for 1) users familiar with a command-line interface and 2) software developers to integrate the matrix-based methods into existing software packages.

Wang, Taoning↗

Deep Spectroscopy with DESI for Photometric Redshift Training and Calibration

Deep spectroscopic samples can improve photometric redshift (photo-z) estimates and reduce uncertainties on redshift distributions. Such improvements can increase the cosmological constraining power of large imaging-based experiments such as the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST) and mitigate what may be a limiting systematic effect. We present results from the “DESI-Deep pilot” program, which was designed to assess the capability of the Dark Energy Spectroscopic Instrument (DESI) on the 4m Mayall telescope to measure redshifts of galaxies as faint as expected lensing samples for early LSST data (m i ≤ 24.5). We find that DESI is remarkably efficient at this task, with redshift success rates comparable to the results of observations from 10 m class telescopes with only ∼2 × longer integration time (rather than ∼8 × longer as would be expected from aperture-area scaling), while simultaneously achieving ∼30 times larger multiplexing. We also find that the signal-to-noise ratio of the spectra scales as expected for background-limited observations even for the longest exposure times (∼7 hr) and faintest targets in the program. These results demonstrate that DESI could provide the definitive redshift sample for the early years of LSST with a modest investment of observing time. Based upon the results of this program, we provide updated predictions for the time required to collect benchmark samples for photo-z training and calibration using a variety of spectroscopic facilities. Finally, we describe a potential “DESI-Deep” survey designed to train and calibrate photo-z’s for imaging experiments, and provide forecasts of its impact on cosmological inference.

Dey, Biprateep [University of Toronto; University ↗

Understanding latent timescales in neural ordinary differential equation models of advection-dominated dynamical systems

The neural ordinary differential equation (ODE) framework has shown considerable promise in recent years in developing highly accelerated surrogate models for complex physical systems characterized by partial differential equations (PDEs). For PDE-based systems, state-of-the-art neural ODE strategies leverage a two-step procedure to achieve this acceleration: a nonlinear dimensionality reduction step provided by an autoencoder, and a time integration step provided by a neural-network based model for the resultant latent space dynamics (the neural ODE). This work explores the applicability of such autoencoder-based neural ODE strategies for PDEs in which advection terms play a critical role. More specifically, alongside predictive demonstrations, physical insight into the sources of model acceleration (i.e., how the neural ODE achieves its acceleration) is the scope of the current study. Such investigations are performed by quantifying the effects of both autoencoder and neural ODE components on latent system time-scales using eigenvalue analysis of dynamical system Jacobians. To this end, the sensitivity of various critical training parameters – de-coupled versus end-to-end training, latent space dimensionality, and the role of training trajectory length, for example – to both model accuracy and the discovered latent system timescales is quantified. Furthermore, this work specifically uncovers the key role played by the training trajectory length (the number of rollout steps in the loss function during training) on the latent system timescales: larger trajectory lengths correlate with an increase in limiting neural ODE time-scales, and optimal neural ODEs are found to recover the largest time-scales of the full-order (ground-truth) system. Demonstrations are performed across fundamentally different unsteady fluid dynamics configurations influenced by advection: (1) the Kuramoto–Sivashinsky equations (2) Hydrogen-Air channel detonations (the compressible reacting Navier–Stokes equations with detailed chemistry), and (3) 2D Atmospheric flow.

Advection-dominated dynamical systems↗

Efficient lattice QCD computation of radiative-leptonic-decay form factors at multiple positive and negative photon virtualities

In previous work [D. Giusti, Methods for high-precision determinations of radiative-leptonic decay form factors using lattice QCD, Phys. Rev. D 107, 074507 (2023)], we showed that form factors for radiative leptonic decays of pseudoscalar mesons can be determined efficiently and with high precision from lattice QCD using the “three-dimensional (3D) method,” in which three-point functions are computed for all values of the current insertion time and the time integral is performed at the data-analysis stage. Here, we demonstrate another benefit of the 3D method: the form factors can be extracted for any number of nonzero photon virtualities from the same three-point functions at no extra cost. We present results for the $D_s → ℓνγ*$ vector form factor as a function of photon energy and photon virtuality, for both positive and negative virtuality, for a single ensemble with 340 MeV pion mass and 0.11 fm lattice spacing. In our analysis, we separately consider the two different time orderings and the different quark flavors in the electromagnetic current. We discuss in detail the behavior of the unwanted exponentials contributing to the three-point functions, as well as the choice of fit models and fit ranges used to remove them for various values of the virtuality. While positive photon virtuality is relevant for decays to multiple charged leptons, negative photon virtuality suppresses soft contributions and is of interest in QCD-factorization studies of the form factors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Spectrally separable photon-pair generation in dispersion engineered thin-film lithium niobate

Existing nonlinear-optic implementations of pure, unfiltered heralded single-photon sources do not offer the scalability required for densely integrated quantum networks. Additionally, lithium niobate has hitherto been unsuitable for such use due to its material dispersion. We engineer the dispersion and the quasi-phasematching conditions of a waveguide in the rapidly emerging thin-film lithium niobate platform to generate spectrally separable photon pairs in the telecommunications band. Such photon pairs can be used as spectrally pure heralded single-photon sources in quantum networks. We estimate a heralded-state spectral purity of >94% based on joint spectral intensity measurements. Further, a joint spectral phase-sensitive measurement of the unheralded time-integrated second-order correlation function yields a heralded-state purity of (86 ± 5)%.

36 MATERIALS SCIENCE↗

Lifting the Garage Door on Spawn, An Open-Source BEM-Controls Engine

Spawn is the latest whole-building energy simulation engine developed by the US Department of Energy, National Labs and industry. Whereas EnergyPlus was designed as a successor to DOE-2, Spawn is not a direct successor of–nor is it intended as an imminent replacement for– EnergyPlus. Instead, Spawn reuses parts of EnergyPlus while supporting new use cases in HVAC and controls. Spawn is intended to provide several capabilities that significantly advance beyond EnergyPlus. It is intended to support the evaluation of novel HVAC and district energy systems in a more physically realistic way. Critically, it can model control in a physically realistic way, using portable specifications that can be compiled for execution on control platforms. Spawn is also intended to support co-simulation in an intrinsic way to enable integration with third-party models. This paper describes the software architecture of Spawn from model authoring to compilation and simulation. It explains how Spawn reuses the envelope and daylighting modules of EnergyPlus and couples them to HVAC and control models from the Modelica Buildings Library using the Functional Mockup Interface (FMI) standard. It presents a number of examples that: i) validate Spawn’s coupled simulation approach by comparing its results to those of EnergyPlus, ii) illustrate the Spawn methodology for modeling and simulating HVAC systems, and iii) evaluate the performance of Spawn’s Quantized State System (QSS) time integration algorithms

Wetter, Michael↗

Quantum simulation of real-space dynamics

Quantum simulation is a prominent application of quantum computers. While there is extensive previous work on simulating finite-dimensional systems, less is known about quantum algorithms for real-space dynamics. We conduct a systematic study of such algorithms. In particular, we show that the dynamics of a d-dimensional Schrödinger equation with η particles can be simulated with gate complexity O ~ (ηdFpoly(log(g'/ϵ))), where ϵ is the discretization error, g' controls the higher-order derivatives of the wave function, and F measures the time-integrated strength of the potential. Compared to the best previous results, this exponentially improves the dependence on ϵ and g' from poly(g'/ϵ) to poly(log(g'/ϵ)) and polynomially improves the dependence on T and d, while maintaining best known performance with respect to η. For the case of Coulomb interactions, we give an algorithm using η 3 (d + η)Tpoly(log(ηdTg'/(Δϵ)))/Δ one- and two-qubit gates, and another using η 3 (4d) d/2 Tpoly(log(ηdTg'/(Δϵ)))/Δ one- and two-qubit gates and QRAM operations, where T is the evolution time and the parameter Δ regulates the unbounded Coulomb interaction. We give applications to several computational problems, including faster real-space simulation of quantum chemistry, rigorous analysis of discretization error for simulation of a uniform electron gas, and a quadratic improvement to a quantum algorithm for escaping saddle points in nonconvex optimization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multiscale Assessment of Caprock Integrity for Geologic Carbon Storage in the Pennsylvanian Farnsworth Unit, Texas, USA

Leakage pathways through caprock lithologies for underground storage of CO2 and/or enhanced oil recovery (EOR) include intrusion into nano-pore mudstones, flow within fractures and faults, and larger-scale sedimentary heterogeneity (e.g., stacked channel deposits). To assess multiscale sealing integrity of the caprock system that overlies the Morrow B sandstone reservoir, Farnsworth Unit (FWU), Texas, USA, we combine pore-to-core observations, laboratory testing, well logging results, and noble gas analysis. A cluster analysis combining gamma ray, compressional slowness, and other logs was combined with caliper responses and triaxial rock mechanics testing to define eleven lithologic classes across the upper Morrow shale and Thirteen Finger limestone caprock units, with estimations of dynamic elastic moduli and fracture breakdown pressures (minimum horizontal stress gradients) for each class. Mercury porosimetry determinations of CO2 column heights in sealing formations yield values exceeding reservoir height. Noble gas profiles provide a “geologic time-integrated” assessment of fluid flow across the reservoir-caprock system, with Morrow B reservoir measurements consistent with decades-long EOR water-flooding, and upper Morrow shale and lower Thirteen Finger limestone values being consistent with long-term geohydrologic isolation. Together, these data suggest an excellent sealing capacity for the FWU and provide limits for injection pressure increases accompanying carbon storage activities.

58 GEOSCIENCES↗

Optical and x-ray characterization of the Daedalus ultrafast x-ray imager

The Daedalus ultrafast x-ray imager is the latest generation in Sandia’s hybrid CMOS detector family. With three frames along an identical line of sight, 1 ns minimum integration time, a higher full well than Icarus, and added features, Daedalus brings exciting new capabilities to diagnostic applications in inertial confinement fusion and high energy density science. In this work, we present measurements of time response, dynamic range, spatial uniformity, pixel cross-talk, and absolute x-ray sensitivity using pulsed optical and x-ray sources. Here, we report a measured 1.5 Me — full well, pixel sensitivity at 9.58 × 10 —7 V/e — , and an estimate of spatial uniformity at ~5% across the sensor array.

47 OTHER INSTRUMENTATION↗

Dose-Weighted Radiographic Source-Spot Size

The spatial resolution of a flash radiograph is determined by the time-integrated size of the radiation source. Since the radiation pulse includes the beam head and tail, which may be substantially larger than the beam body, the radiation source is larger than the focused flattop of the electron beam pulse. Moreover, the radiation dose rate is weighted by both beam current and electron energy, and it is the radiation dose that forms the radiograph image. Therefore, in order to estimate radiographic resolution from computer simulations of a focused electron beam, one should use a dose weighted moment analysis of the radiation source distribution. This approach is used to show that, in the absence of beam effects other than space-charge and emittance, radiographic resolution for Scorpius exceeds requirements, regardless of reasonable timing options for current and gap-voltage pulses.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

PETSc Users Manual (Rev. 3.13)

This manual describes the use of PETSc for the numerical solution of partial differential equations and related problems on high-performance computers. The Portable, Extensible Toolkit for Scientific Computation (PETSc) is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all message-passing communication. PETSc includes an expanding suite of parallel linear solvers, nonlinear solvers, and time integrators that may be used in application codes written in Fortran, C, C++, and Python. PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or dbx, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than \rolling them" yourself.

97 MATHEMATICS AND COMPUTING↗

Neural Networks for Nuclear Reactions in MAESTROeX

We demonstrate the use of neural networks to accelerate the reaction steps in the MAESTROeX stellar hydrodynamics code. A traditional MAESTROeX simulation uses a stiff ODE integrator for the reactions; here, we employ a ResNet architecture and describe details relating to the architecture, training, and validation of our networks. Our customized approach includes options for the form of the loss functions, a demonstration that the use of parallel neural networks leads to increased accuracy, and a description of a perturbational approach in the training step that robustifies the model. We test our approach on millimeter-scale flames using a single-step, 3-isotope network describing the first stages of carbon fusion occurring in Type Ia supernovae. We train the neural networks using simulation data from a standard MAESTROeX simulation, and show that the resulting model can be effectively applied to different flame configurations. This work lays the groundwork for more complex networks, and iterative time-integration strategies that can leverage the efficiency of the neural networks.

79 ASTRONOMY AND ASTROPHYSICS↗

Greenhouse production contributes to pesticide occurrences in Swedish streams

Greenhouse and other covered cultivation systems have increased globally over the past several decades, leading to considerably improved product quality and productivity per land area unit. However, there is a paucity in information regarding the environmental impacts of covered production systems, especially regarding pesticides entering the surrounding environment. Aiming to address this knowledge gap, we collected grab samples downstream of greenhouses from seven Swedish streams every 14 days during a 12 month period. In three of the streams, samples were also taken upstream of the greenhouses and in four of the streams time-integrated samples were collected by TIMFIE samplers in the period between grab sampling occasions. The samples were analyzed for 28 substances (27 that were permitted for use in greenhouse production systems in Sweden and one degradation product to a permitted substance). Pesticide use journals were collected from the greenhouse producers for the 12 month period. The results were examined for indications of greenhouse contributions to detection frequencies, maximum and average concentrations, and potential ecotoxicicity in several ways: (1) comparing locations downstream of greenhouses with registered use of a substance with those without registered use, (2) comparing results from this study with those from the Swedish environmental monitoring program of pesticides in surface water from catchments with no greenhouses from the same period and region, (3) comparing concentration trends with registered pesticide application times in the greenhouses, and (4) comparing up- and downstream concentrations. The results strongly suggest that greenhouse applications do contribute to pesticide occurrences, maximum and median concentrations for most of the pesticides included in this study, and to potential toxicity to aquatic organisms for several of them, most notably imidacloprid, acetamiprid, carbendazim, and pirimicarb.

54 ENVIRONMENTAL SCIENCES↗

Gamma-Ray Imaging to Monitor Accumulating Holdup

Over time, nuclear materials can accumulate in process equipment and keeping track of such material is of interest to many stakeholder communities including safeguards, process monitoring, criticality safety, and radiological safety. Passive gamma-ray imaging provides a visual and quantitative tool that can be applied to this problem, allowing one to determine the amount of material present from the images. In this work we emulate the accumulation of highly enriched uranium (HEU) in a large rectangular stainless-steel duct by stacking 5, 27 × 50 cm2 cards, each of which has approximately 11 g of 235U in oxide form distributed across its surface, behind each other. A series of images were collected using a coded-aperture gamma-ray imager, starting with a single card, adding a second card, etc. Long integrations were collected to allow evaluating the performance of the approach as a function of integration time. The images are analyzed using first principles to determine the amount of material in each image. The results are compared to the known amount of material in each configuration, and as a function of time. The analysis includes a first-principles error propagation and to further evaluate the uncertainty, a boot-strap sampling approach is also applied to the data. The work offers a proof of concept for unattended monitoring of uranium holdup with a quantitative gamma-ray imaging technique.

Ziock, Klaus-Peter↗