Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Time-integration”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Measurement of the time-integrated 𝐶⁢𝑃 asymmetry in 𝐷 0 → 𝐾$^{0}_{S}$𝐾$^{0}_{S}$ decays using opposite-side flavor tagging at Belle and Belle II

We measure the time-integrated 𝐶⁢𝑃 asymmetry in 𝐷 0 → 𝐾$^{0}_{S}$𝐾$^{0}_{S}$ decays reconstructed in 𝑒 + ⁢𝑒 − → $c\bar{c}$ events collected by the Belle and Belle II experiments. The corresponding data samples have integrated luminosities of 980 and 428 fb −1 , respectively. To infer the flavor of the 𝐷 0 meson, we exploit the correlation between the flavor of the reconstructed decay and the electric charges of particles reconstructed in the rest of the 𝑒 + ⁢𝑒 − → $c\bar{c}$ event. This results in a sample which is independent from any other previously used at Belle or Belle II. The result, 𝐴 𝐶⁢𝑃 ⁡(𝐷 0 → 𝐾$^{0}_{S}$𝐾$^{0}_{S}$)=(1.3±2.0±0.2)%, where the first uncertainty is statistical and the second systematic, is consistent with previous determinations and with 𝐶⁢𝑃 symmetry.

CP violation↗
DOE OSTI2025

Measurement of the time-integrated C P asymmetry in D 0 → K S 0 K S 0 decays using Belle and Belle II data

We measure the time-integrated C P asymmetry in D 0 → K S 0 K S 0 decays reconstructed in e + e − → c c ¯ events collected by the Belle and Belle II experiments. The corresponding data samples have integrated luminosities of 980 and 428 fb − 1 , respectively. The D 0 decays are required to originate from the D * + → D 0 π + decay, which determines the charm flavor at production time. A control sample of D 0 → K + K − decays is used to correct for production and detection asymmetries. The result, ( − 1.4 ± 1.3 ( stat ) ± 0.1 ( syst ) ) % , is consistent with previous determinations and with C P symmetry. Published by the American Physical Society 2025

Adachi, I. (ORCID:0000000322870173)↗
DOE OSTI2025

Measurement of the time-integrated 𝐶⁢𝑃 asymmetry in 𝐷 0 → 𝜋 0 ⁢𝜋 0 decays at Belle II

We measure the time-integrated 𝐶⁢𝑃 asymmetry, 𝐴 𝐶⁢𝑃 , in 𝐷 0 → 𝜋 0 ⁢𝜋 0 decays reconstructed in 𝑒 + ⁢𝑒 − → $𝑐⁢\bar{𝑐}$ events collected by Belle II during 2019–2022. The data corresponds to an integrated luminosity of 428 fb −1 . The 𝐷 0 decays are required to originate from the flavor-conserving 𝐷 *+ → 𝐷 0 ⁢𝜋 + decay to determine the charm flavor at production time. Control samples of 𝐷 0 → 𝐾 − ⁢𝜋 + decays, with or without an associated pion from a 𝐷 *+ decay, are used to correct for detection asymmetries. The result, 𝐴 𝐶⁢𝑃 ⁡(𝐷 0 → 𝜋 0 ⁢𝜋 0 ) = (0.30 ± 0.72 ± 0.20)%, where the first uncertainty is statistical and the second systematic, is consistent with 𝐶⁢𝑃 symmetry.

CP violation↗
DOE OSTI2026

Time-integrated Southern-sky Neutrino Source Searches with 10 yr of IceCube Starting-track Events at Energies Down to 1 TeV

In the IceCube Neutrino Observatory, a signal of astrophysical neutrinos is obscured by backgrounds from atmospheric neutrinos and muons produced in cosmic-ray interactions. IceCube event selections used to isolate the astrophysical neutrino signal often focus on the morphology of the light patterns recorded by the detector. The analyses presented here use the new IceCube Enhanced Starting Track Event Selection (ESTES), which identifies events likely generated by muon–neutrino interactions within the detector geometry, focusing on neutrino energies of 1–500 TeV with a median angular resolution of 1.4°. Selecting for starting-track events filters out not only the atmospheric-muon background but also the atmospheric-neutrino background in the southern sky. This improves IceCube’s muon–neutrino sensitivity to southern-sky neutrino sources, especially for Galactic sources that are not expected to produce a substantial flux of neutrinos above 100 TeV. In this work, the ESTES sample was applied for the first time to search for astrophysical sources of neutrinos, including a search for diffuse neutrino emission from the Galactic plane. No significant excesses were identified from any of the analyses; however, constraining limits are set on the hadronic emission from TeV gamma-ray Galactic plane objects and models of the diffuse Galactic plane neutrino flux.

Abbasi, R. [Loyola University, Chicago, IL (United↗
DOE OSTI2025

Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers

Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.

Flood Simulation↗
DOE OSTI2026

A low-latitude, high-elevation 10Be production-rate calibration from the Rwenzori Mountains, Uganda

Cosmogenic nuclide surface-exposure dating has become a key tool in geomorphology. Robust application of the technique relies on determination of local nuclide production rates using independently dated calibration sites. However, there are few existing calibration sites in certain regions, particularly in low-latitude, high-elevation regions where cosmogenic nuclide production is theoretically more sensitive to changes in Earth's magnetic field and atmospheric conditions. We present a new cosmogenic beryllium-10 (10Be) production-rate calibration from the equatorial Rwenzori Mountains in Uganda. The calibration is based on radiocarbon dating of basal lake sediments (∼21.3 kyr BP) from Lake Mahoma (2900 m asl) and measured 10Be concentrations from boulders on the crest of a moraine which dams the lake. Our results indicate a time-integrated 10Be production rate of 15.83 ± 0.81 at/g/yr at Lake Mahoma. We then use these calibration data and two public online calculators to calculate ages for moraines in a separate Rwenzori valley for which there exists a minimum-limiting radiocarbon age constraint (∼12 kyr BP) on surficial glacial sediments. This enables us to assess the performance of different parameters used to scale production rates spatially and temporally (e.g., scaling frameworks, geomagnetic field reconstructions, atmospheric models) for surface-exposure age calculation in the region. Our results highlight the sensitivity of low-latitude cosmogenic nuclide surface-exposure ages to discrete calculation parameters, such as scaling framework or reconstructions of Earth's magnetic field, and the utility of determining production rates at sites where limiting age data exists for other nearby landforms.

Anderson, N [Dartmouth College]↗
DOE OSTI2025

Sensitized Near-Infrared Emission of SiGe Nanocrystals via Heterostructuring with Quasi Two-Dimensional Perovskite

The near-infrared (NIR) emission of silicon–germanium alloy nanocrystals (SiGe NCs) was sensitized by heterostructuring with a quasi two-dimensional (Q-2D) perovskite (CsPbBr 3 blended with butylammonium bromide (BABr)). Colloidal SiGe NCs were synthesized by using a nonthermal plasma method with chloride precursors. As-synthesized SiGe NCs showed no detectable photoluminescence (PL). Embedding the SiGe NCs within the perovskite matrix via spin coating led to a heterostructure exhibiting dual emissions: green emission at 520 nm from the perovskite and NIR emission at 1035 nm from the SiGe NCs. Time-integrated and resolved PL measurements and transient absorption spectroscopy revealed energy transfer from the perovskite to the SiGe NCs, which sensitized the NIR emission from the SiGe NCs and quenched the green PL for the perovskite. These results demonstrate that Ge alloying can effectively tune the bandgap of Si NCs and highlight the feasibility of heterostructuring perovskites and photoluminescent NCs to enhance or activate their PL. In conclusion, this approach broadens their potential applications in the NIR region.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Immersion Freezing in Particle-Based Aerosol-Cloud Microphysics: A Probabilistic Perspective on Singular and Time-Dependent Models

Cloud droplets containing immersed ice-nucleating particles (INPs) may freeze at temperatures above the homogeneous freezing threshold temperature in a process referred to as immersion freezing. In modeling studies, immersion freezing is often described using either so-called “singular” or “time-dependent” parameterizations. Here, we compare both approaches and discuss them in the context of probabilistic particle-based (super-droplet) cloud microphysics modeling. First, using a box model, we contrast how both parameterizations respond to idealized ambient cooling rate profiles and quantify the impact of the polydispersity of the immersed surface spectrum on the frozen fraction evolution. Presented simulations highlight that the singular approach, constituting a time-integrated form of a more general time-dependent approach, is only accurate under a limited range of ambient cooling rates. The time-dependent approach is free from this limitation. Second, using a prescribed-flow two-dimensional cloud model, we illustrate the macroscopic differences in the evolution in time of ice particle concentrations in simulations with flow regimes relevant to ambient cloud conditions. The flow-coupled aerosol-budget-resolving simulations highlight the benefits and challenges of modeling cloud condensation nuclei activation and immersion freezing on insoluble ice nuclei with super-particle methods. The challenges stem, on the one hand, from heterogeneous ice nucleation being contingent on the presence of relatively sparse immersed INPs, and on the other hand, from the need to represent a vast population of particles with relatively few so-called super particles (each representing a multiplicity of real particles). We discuss the critical role of the sampling strategy for particle attributes, including the INP size, the freezing temperature (for singular scheme) and the multiplicity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Addendum: Unified framework for open quantum dynamics with memory

This Addendum presents a detailed analysis of the discretization error in time-integration and time-derivative that appear in the Nakajima-Zwanzig equation. This was brought to our attention by Makri et al. [arXiv:2410.08239]. Our analysis in the Addendum shows that the relationship derived in our earlier work [Nat. Commun. 15, 8087 (2024)] is valid within the choice of discretization and is not contaminated by the discretization error.

Science & Technology - Other Topics↗
DOE OSTI2026

Measured Indoor Nitrogen Oxides in Homes with a Child with Asthma and Gas or Induction Electric Cooking

The Cooking Energy and Ventilation Impacts on Children's Asthma (CEVICA) study is a randomized control trial investigating the effects of replacing gas with induction electric ranges in the homes of children with asthma in California's San Joaquin Valley. Indoor air quality parameters and respiratory health indicators were measured over 2-week intensive at Baseline and after consecutive 3-month study phases, with half getting electric cooking at the start of Phase 1 and others in Phase 2. Indoor air measurements included time-integrated NO 2 and NO X by passive sampler and time-resolved NO 2 by electrochemical sensors. In the first 40 homes, data were collected during 61 intensives for gas cooking and 55 for induction. Cooking was identified using temperature sensors above the cooktop. Pollutant events were identified from sharp rises in concentrations Time integrated NOX species were lower with electric cooking, with mean differences of 14.2 ppb for NO 2 (95% CI: 9.6, 18.7; p < 0.001), 55.6 ppb for NO X (95% CI: 30.4, 80.9; p < 0.001), and 41.5 ppb for derived NO (95% CI: 19.7, 63.2; p < 0.001). A paired, within home analysis showed larger reductions in Group 2 (gas to induction) than Group 1 that remained electric across Phase 1 and Phase 2 for NO 2 (mean difference in As of 10.5 ppb; 95% CI: 5.6, 15.4; p < 0.001), NO X (32.3 ppb; 95% CI: 10.0, 54.5; p=0.007), and derived NO (21.8 ppb; 95% CI: 3.6, 40.0; p=0.022). Compared to electric cooking, gas had much higher rates of associated NO 2 events, and larger above-baseline NO 2 peaks. These preliminary results are consistent with prior findings that shifting to induction cooking can substantially lower indoor NO 2 compared with gas cooking.

Fang, Yi↗
DOE OSTI2026

Measured indoor PM2.5, black carbon, and oxidative potential before and after replacing gas with induction cooking in asthmatic households

Cooking is a major source of fine particulate matter (PM) in homes and evidence to date is inconclusive about the impact of cooking fuel on measures of residential PM exposure. The Cooking Energy and Ventilation Impacts on Children's Asthma (CEVICA) study measured cooking frequency, range hood use, indoor air quality (IAQ) and respiratory health indicators of children with asthma living in homes with gas stoves in California's San Joaquin Valley. Intensive measurements occurred over three 2-week intensive periods: at baseline and at the end of two consecutive 3-month study phases. Participants were randomly assigned to have their gas stoves replaced with electric induction at the start of Phase 1 or Phase 2. As part of the IAQ assessment, we collected particulate matter on Teflon filters using ultrasonic personal air samplers (UPAS). The filters were analyzed to quantify time-integrated PM 2.5 mass concentration, black carbon (BC) oxidative potential (OP). OP was measured using the SLF-OH assay, an acellular chemical assay that quantifies production of OH radicals in simulated lung fluid (SLF). Across all comparisons, PM 2.5 mass tended to be higher during gas cooking than induction electric cooking, although only one transition (Baseline to Phase 1) was statistically discernible at p<0.05. BC and mass-normalized OP showed no measurable differences across stove types or across phases. These preliminary findings suggest that stove fuel alone may not be the dominant determinant of indoor PM 2.5 concentrations in this sample, as previous studies have reported that indoor PM 2.5 during cooking can also be influenced by factors such as cooking practices, food and oil type, kitchen ventilation, and outdoor infiltration.

Lin, Zhuoying↗
DOE OSTI2025

VERITAS : A density-functional theory-based multiband kinetic model for understanding x-ray spectroscopy of dense plasmas

X-ray spectroscopy has long been a powerful diagnostic tool for hot, dilute plasmas, providing insights into plasma conditions by measuring line shifts and broadenings of atomic transitions. The technique critically depends on the accuracy of atomic physics models used to interpret spectroscopic measurements for inferring plasma properties such as free-electron density and temperature. Over the past decades, the atomic and plasma physics communities have developed robust atomic physics models to account for various processes in hot, dilute classical plasmas. While these models have been successful in that regime, their applicability becomes uncertain when interpreting x-ray spectroscopy experiments of above-solid-density plasmas. Given that finite-temperature density-functional theory (DFT) offers a more accurate description of dense plasma environments, we present the development of a DFT-based multi-band kinetic model, VERITAS, designed to improve the interpretation of x-ray spectroscopic measurements in high-density plasmas produced by laser-driven spherical implosions. This work details the VERITAS model and its application to both time-integrated and time-resolved x-ray spectra from implosion experiments on OMEGA. The advantages and limitations of the VERITAS model will also be discussed, along with potential directions for advancing x-ray spectroscopy of dense and superdense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

The role of focusing geometry in MeV x-ray production from petawatt laser–solid interaction

Relativistic laser–plasma interactions provide a compact and flexible route to generating bright, ultrashort pulses of MeV x rays, with applications in high-energy density science, nuclear physics, and radiography. Despite extensive study of intensity scaling in laser-driven electron acceleration, the role of focusing geometry and focal-volume effects in MeV radiation production remains insufficiently understood. Here, we present experimental and kinetic simulation results from the Texas Petawatt Laser (120 J, 140 fs) in which the focusing geometry (f/3 or f/1.5) is varied, while the laser energy and pulse duration are fixed. Experimentally, the f/3 geometry produces approximately three times more MeV radiation than the f/1.5 geometry, despite its twofold lower nominal vacuum intensity. This result is based on the time-integrated radiation per unit solid angle, as measured along the diagnostic line of sight. Three-dimensional particle-in-cell simulations reproduce this trend when a modest (10s of μm) effective focal plane shift is introduced, demonstrating that relativistic laser–plasma coupling is highly sensitive to focal geometry in the presence of a preplasma. This behavior is consistent with differences in interaction length and effective intensity at the critical surface between the two configurations. The sensitivity of the tightly focused f/1.5 configuration ($z_R ≈$⁠ 5 μ m) reflects the combined influence of thermal lensing and an extended preplasma, while the f/3 geometry (⁠$z_R ≈$ 80 μ m) remains comparatively robust. These results demonstrate a breakdown of conventional intensity scaling and identify focusing geometry as a critical control parameter for MeV electron and x-ray generation at petawatt powers. Implications for laser–plasma-based radiographic facilities are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Search for Λ − Λ ¯ oscillation in J / ψ → Λ Λ ¯ decay

Using ( 10087 ± 44 ) × 10 6 J / ψ decays collected by the BESIII detector at the BEPCII collider, we search for baryon number violation via Λ − Λ ¯ oscillation in the decay J / ψ → Λ Λ ¯ . No evidence for Λ − Λ ¯ oscillation is observed. The upper limit on the time-integrated probability of Λ − Λ ¯ oscillation is estimated to be 1.4 × 10 − 6 , corresponding to an oscillation parameter less than 2.1 × 10 − 18 GeV at 90% confidence level. Published by the American Physical Society 2025

Ablikim, M.↗
DOE OSTI2025

Graph-based Reversible Evaluation and Tangents Library

GRETL is a C++ library for evaluation, re-evaluation and algorithmic differentiation of functional operations on an arbitrary computational graph with limited memory usage. Similar to popular machine learning frameworks in Python, like PyTorch and JAX, it tracks and stores both operations and output data as functions are evaluated. Once this composition of functions is built up, the entire chain of operations can be back propagated to compute sensitivities of the final result with respect to any number of inputs. In contrast to most machine learning applications, memory usage becomes the bottleneck for back propagation in many physics applications, especially for time-dependent PDEs. Dynamic check pointing becomes essential. An important distinguishing feature of GRETL is its ability to limit the maximum memory usage by automatically dynamic checkpointing the data output for each graph operation (see Wang, Moin, Iaccarino, 2009). During backpropagation, parts of the graph that are no longer in memory are automatically re-evaluated from upstream checkpointed states as needed for derivative sensitivity calculations (or more precisely, for vector-Jacobian products). GRETL is particularly beneficial for applications, such as coupled multi-physics, where deriving adjoint-based sensitivities and managing checkpoint memory across modules becomes onerous. Cases which can be readily handled by the GRETL library include: different time-integration algorithms per physics (e.g., coupled predictor-corrector algorithms, IMEX, etc.), sub-cycling, asynchronous integrators, state dependent timestep sizes, iterative solvers and coupling algorithms, controller algorithms, and more.

Tupek, MichaelR [Lawrence Livermore National Labor↗
DOE OSTI2026

Spectral Flux Measurements on SNS Direct Geometry Spectrometers at IRP-2 End-of-Life

The Inner Reflector Plug (IRP) of the Spallation Neutron Source (SNS) was replaced during the SNS shutdown period beginning 25 November 2025. This IRP, IRP-2, was first placed in service in May 2018, and was in use for 42.5 GW-hrs of operation, somewhat longer than its original design life of 30 GW-hrs, or even a refined estimate of 33 GW-hrs for the water moderator poison plate and 39 GW-hrs on the decoupled-poisoned hydrogen moderator decoupler (Gallmeier, Lu, and Iverson 2018). The time-integrated power history for IRP-2 is tracked in Figure 1.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

All-sky Neutrino Point-source Search with IceCube Combined Track and Cascade Data

Despite extensive efforts, discovery of high-energy astrophysical neutrino sources remains elusive. We present an event-level simultaneous maximum likelihood analysis of tracks and cascades using IceCube data collected from 2008 April 6 to 2022 May 23 to search the whole sky for neutrino sources, and using a source catalog, for coincidence of neutrino emission with gamma-ray emission. This is the first time a simultaneous fit of different detection channels is used to conduct a time-integrated all-sky scan with IceCube. Combining all-sky tracks, with superior pointing power and sensitivity in the northern sky, with all-sky cascades, with good energy resolution and sensitivity in the southern sky, we have developed the most sensitive point-source search to date by IceCube that targets the entire sky. The most significant point in the northern sky aligns with NGC 1068, a Seyfert II galaxy, which, from the catalog search, shows a 3.5σ excess over background after accounting for trials. The most significant point in the southern sky does not align with any source in the catalog and is not significant after accounting for trials. A search for the single most significant Gaussian flare at the locations of NGC 1068, PKS 1424+240, and the southern highest-significance point shows results consistent with expectations for steady emission. Notably, this is the first time that a flare shorter than four years has been excluded as being responsible for NGC 1068’s emergence as a neutrino source. Our results show that combining tracks and cascades when conducting neutrino source searches improves sensitivity and can lead to new discoveries.

Abbasi, R. [Loyola University, Chicago, IL (United↗