SEARCH · Search NASA

Results for “Delaunay”

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 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

TetCrust WILD: Weighted Delaunay meshing in the wild

We introduce the first conforming Voronoi mesher capable of handling unclean input geometries encountered 'in the wild'. Our algorithm handles non-watertight, non-manifold, and self-intersecting inputs without requiring any manual cleanup, providing a standalone, end-to-end meshing framework. The framework enables simultaneous Voronoi and weighted Delaunay meshing, allowing users to generate polyhedral or tetrahedral elements based on application requirements, and produces high-quality elements

Ebeida, Mohamed Salah [Sandia National Laboratorie↗

DOE OSTI2024

Algorithm 1049: The Delaunay Density Diagnostic

Accurate approximation of a real-valued function depends on two aspects of the available data: the density of inputs within the domain of interest and the variation of the outputs over that domain. There are few methods for assessing whether the density of inputs is sufficient to identify the relevant variations in outputs—i.e., the “geometric scale” of the function—despite the fact that sampling density is closely tied to the success or failure of an approximation method. In this article, we introduce a general purpose, computational approach to detecting the geometric scale of real-valued functions over a fixed domain using a deterministic interpolation technique from computational geometry. The algorithm is intended to work on scalar data in moderate dimensions (2–10). Our algorithm is based on the observation that a sequence of piecewise linear interpolants will converge to a continuous function at a quadratic rate (in L 2 norm) if and only if the data are sampled densely enough to distinguish the feature from noise (assuming sufficiently regular sampling). We present numerical experiments demonstrating how our method can identify feature scale, estimate uncertainty in feature scale, and assess the sampling density for fixed (i.e., static) datasets of input–output pairs. Finally, we include analytical results in support of our numerical findings and have released lightweight code that can be adapted for use in a variety of data science settings.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Dark matter substructure or source model systematics? A case study of cluster lens Abell S1063

Mapping the small-scale structure of the universe through gravitational lensing is a promising tool for probing the particle nature of dark matter. Curved Arc Basis (CAB) has been proposed as a local lensing formalism in galaxy clusters, with the potential to detect low-mass dark matter substructure. In this work, we analyse the cluster lens Abell S1063 in search of dark matter substructure with the CAB formalism, using multiband imaging data from James Webb Space Telescope ( JWST ). We use two different source modelling methods: shapelets and pixel-based source reconstruction based on Delaunay triangulation. We find that source modelling systematics from shapelets result in a disagreement between CAB parameters measured from different filters. Source modelling with Delaunay significantly alleviates this systematic, as seen in the improvement in agreement across filters. We also find that inadequate complexity in source modelling can result in convincing spurious detections of dark matter substructure from strong gravitational lenses, as seen by our $\Delta \text{BIC} > 20$ measurement of a $M \sim 10^{10}$ ${\rm M}_{\odot }$ subhalo with shapelets, a spurious detection that is not reproduced with Delaunay source modelling. We demonstrate that multiband analysis with different JWST filters is key for disentangling source and lens model systematics from dark matter substructure detections.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

An invertible seven-dimensional Dirichlet cell characterization of lattices

Characterization of crystallographic lattices is an important tool in structure solution, crystallographic database searches and clustering of diffraction images in serial crystallography. Characterization of lattices by Niggli-reduced cells (based on the three shortest non-coplanar lattice vectors) or by Delaunay-reduced cells (based on four non-coplanar vectors summing to zero and all meeting at obtuse or right angles) is commonly performed. The Niggli cell derives from Minkowski reduction. The Delaunay cell derives from Selling reduction. All are related to the Wigner–Seitz (or Dirichlet, or Voronoi) cell of the lattice, which consists of the points at least as close to a chosen lattice point as they are to any other lattice point. The three non-coplanar lattice vectors chosen are here called the Niggli-reduced cell edges. Starting from a Niggli-reduced cell, the Dirichlet cell is characterized by the planes determined by 13 lattice half-edges: the midpoints of the three Niggli cell edges, the six Niggli cell face-diagonals and the four body-diagonals, but seven of the lengths are sufficient: three edge lengths, the three shorter of each pair of face-diagonal lengths, and the shortest body-diagonal length. These seven are sufficient to recover the Niggli-reduced cell.

Chemistry↗

DOE OSTI2024

T RI M E ++: Multi-threaded triangular meshing in two dimensions

We present T RI M E ++, a multi-threaded software library designed for generating two-dimensional meshes for intricate geometric shapes using the Delaunay triangulation. Multi-threaded parallel computing is implemented throughout the meshing procedure, making it suitable for fast generation of large-scale meshes. Three iterative meshing algorithms are implemented: the DistMesh algorithm, the centroidal Voronoi diagram meshing, and a hybrid of the two. We compare the performance of the three meshing methods in T RI M E ++, and show that the hybrid method retains the advantages of the other two. The software library achieves significant parallel speedup when generating large-scale meshes containing between 10 4 to 10 7 points. T RI M E ++ can handle complicated geometries and generates adaptive meshes of high quality.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Remark on Algorithm 1012: Computing Projections with Large Datasets

In ACM TOMS Algorithm 1012, the DELAUNAYSPARSE software is given for performing Delaunay interpolation in medium to high dimensions. When extrapolating outside the convex hull of the training set, DELAUNAYSPARSE calls the nonnegative least squares solver DWNNLS to compute projections onto the convex hull. However, DWNNLS and many other available sum-of-squares optimization solvers were not intended for usage with many variable problems, which result from the large training sets that are typical in machine learning applications. Thus, a new PROJECT subroutine is given, based on the highly customizable quadratic program solver BQPD. This solution is shown to be as robust as DELAUNAYSPARSE for projection onto both synthetic and real-world datasets, where other available solvers frequently fail. Although it is intended as an update for DELAUNAYSPARSE, due to the difficulty and prevalence of the problem, this solution is likely to be of external interest as well.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

On the statistical theory of self-gravitating collisionless dark matter flow: Scale and redshift variation of velocity and density distributions

The statistics of velocity and density fields are crucial for cosmic structure formation and evolution. Here, this paper extends our previous work on the two-point second-order statistics for the velocity field [Phys. Fluids 35, 077105 (2023)] to one-point probability distributions for both density and velocity fields. The scale and redshift variation of density and velocity distributions are studied by a halo-based non-projection approach. First, all particles are divided into halo and out-of-halo particles so that the redshift variation can be studied via generalized kurtosis of distributions for halo and out-of-halo particles, respectively. Second, without projecting particle fields onto a structured grid, the scale variation is analyzed by identifying all particle pairs on different scales $r$. We demonstrate that: (i) Delaunay tessellation can be used to reconstruct the density field. The density correlation, spectrum, and dispersion functions were obtained, modeled, and compared with the N-body simulation; (ii) the velocity distributions are symmetric on both small and large scales and are non-symmetric with a negative skewness on intermediate scales due to the inverse energy cascade on small scales with a constant rate $\varepsilon_u$; (iii) On small scales, the even order moments of pairwise velocity $\Delta u_L$ follow a two-thirds law $\propto{(-\varepsilon_ur)}^{2/3}$, while the odd order moments follow a linear scaling $\langle(\Delta u_L)^{2n+1}\rangle=(2n+1)\langle(\Delta u_L)^{2n}\rangle\langle\Delta u_L\rangle\propto{r}$; (iv) The scale variation of the velocity distributions was studied for longitudinal velocities $u_L$ or $u_L^{'}$, pairwise velocity (velocity difference) $\Delta u_L$=$u_L^{'}$-$u_L$ and velocity sum $\Sigma u_L$=$u^{'}_L$+$u_L$. Fully developed velocity fields are never Gaussian on any scale, despite that they can initially be Gaussian; (v) On small scales, $u_L$ and $\Sigma u_L$ can be modeled by a $X$ distribution to maximize the entropy of the system. The distribution of $\Delta u_L$ can be different; (vi) On large scales, $\Delta u_L$ and $\Sigma u_L$ can be modeled by a logistic or a $X$ distribution, while $u_L$ has a different distribution; (vii) the redshift variation of the velocity distributions follows the evolution of the $X$ distribution involving a shape parameter $\alpha(z)$ decreasing with time.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Filaments of the Slime Mold Cosmic Web and How They Affect Galaxy Evolution

Abstract We present a novel approach for identifying cosmic web filaments within theDisPerSEstructure identification framework, using cosmic density field estimates from the Monte Carlo Physarum Machine (MCPM), inspired by the slime mold organism. We apply our method to the IllustrisTNG (TNG100) cosmological simulations and investigate the impact of filaments on galaxies. The MCPM density field is superior to the Delaunay tessellation field estimator in tracing the true underlying matter distribution and allows filaments to be identified with higher fidelity, finding more low-prominence/diffuse filaments. Using our new filament catalogs, we find that ≳90% of galaxies are located within ∼1.5 Mpc of a filamentary spine, with little change in the median star formation activity with distance to the nearest filament. Instead, we uncover a differential effect of the local filament line density, Σ fil (MCPM)—the total MCPM overdensity per unit length along a filament segment—on galaxy formation: most galaxies are quenched and gas-poor near high-line density filaments atz≤ 1. At earlier times, the filamentary environment appears to have no effect on galactic gas supply and quenching. Atz= 0, quenching in log ( M * / M ⊙ ) ≳ 10.5 galaxies is mainly driven by mass, while lower-mass galaxies are significantly affected by the filament line density. Satellites are far more susceptible to filaments than centrals. The local environments of massive halos are not sufficient to account for the effect of filament line density on gas removal and quenching. Our new approach holds great promise for observationally identifying filaments from galaxy surveys such as SDSS and DESI.

Astronomy & Astrophysics↗

DOE OSTI2023

Intruder configurations in 29 Ne at the transition into the island of inversion: Detailed structure study of 28 Ne

Detailed γ-ray spectroscopy of the exotic neon isotope 28 Ne has been performed for the first time using the one-neutron removal reaction from 29 Ne on a liquid hydrogen target at 240 MeV/nucleon. Based on an analysis of parallel momentum distributions, a level scheme with spin-parity assignments has been constructed for 28 Ne and the negative-parity states are identified for the first time. The measured partial cross sections and momentum distributions reveal a significant intruder p-wave strength providing evidence of the breakdown of the N=20 and N=28 shell gaps. Only a weak, possible f-wave strength was observed to bound final states. Large-scale shell-model calculations with different effective interactions do not reproduce the large p-wave and small f-wave strength observed experimentally, indicating an ongoing challenge for a complete theoretical description of the transition into the island of inversion along the Ne isotopic chain.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Exploration of mass splitting and muon/tau mixing parameters for an eV-scale sterile neutrino with IceCube

We present the first three-parameter fit to a 3+1 sterile neutrino model using 7.634 years of data from the IceCube Neutrino Observatory on $ν$ μ + $\overline{ν}$ μ charged-current interactions in the energy range 500–9976 GeV. Our analysis is sensitive to the mass-squared splitting between the heaviest and lightest mass state (Δm$^{2}_{41}$), the mixing matrix element connecting muon flavor to the fourth mass state (|U μ4 | 2 ), and the element connecting tau flavor to the fourth mass state (|U $\tau$4 | 2 ). Predicted propagation effects in matter enhance the signature through a resonance as atmospheric neutrinos from the Northern Hemisphere traverse the Earth to the IceCube detector at the South Pole. The remaining sterile neutrino matrix elements are left fixed, with |U e4 | 2 = 0 and δ 14 = 0, as they have a negligible effect, and δ 24 = π is set to give the most conservative limits. The result is consistent with the no-sterile neutrino hypothesis with a probability of 4.3%. Profiling the likelihood of each parameter yields the 90% confidence levels: 2.4eV 2 < Δm$^{2}_{41}$) < 9.6eV 2 , 0.0081< |U μ4 | 2 < 0.10, and |U $\tau$4 | 2 < 0.035, which narrows the allowed parameter-space for |U $\tau$4 | 2 . However, the primary result of this analysis is the first map of the 3+1 parameter space exploring the interdependence of Δm$^{2}_{41}$, |U μ4 | 2 , and |U $\tau$4 | 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

First observation of 28 O

Subjecting a physical system to extreme conditions is one of the means often used to obtain a better understanding and deeper insight into its organization and structure. In the case of the atomic nucleus, one such approach is to investigate isotopes that have very different neutron-to-proton ($N/Z$) ratios than in stable nuclei. Light, neutron-rich isotopes exhibit the most asymmetric $N/Z$ ratios and those lying beyond the limits of binding, which undergo spontaneous neutron emission and exist only as very short-lived resonances (about 10 -21 s), provide the most stringent tests of modern nuclear-structure theories. Here we report on the first observation of 28 O and 27 O through their decay into 24 O and four and three neutrons, respectively. The 28 O nucleus is of particular interest as, with the $Z$ = 8 and $N$ = 20 magic numbers, it is expected in the standard shell-model picture of nuclear structure to be one of a relatively small number of so-called ‘doubly magic’ nuclei. Furthermore, both 27 O and 28 O were found to exist as narrow, low-lying resonances and their decay energies are compared here to the results of sophisticated theoretical modelling, including a large-scale shell-model calculation and a newly developed statistical approach. In both cases, the underlying nuclear interactions were derived from effective field theories of quantum chromodynamics. Finally, it is shown that the cross-section for the production of 28 O from a 29 F beam is consistent with it not exhibiting a closed $N$ = 20 shell structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

D-Egg: a dual PMT optical module for IceCube

The D-Egg, an acronym for "Dual optical sensors in an Ellipsoid Glass for Gen2," is one of the optical modules designed for future extensions of the IceCube experiment at the South Pole. The D-Egg has an elongated-sphere shape to maximize the photon-sensitive effective area while maintaining a narrow diameter to reduce the cost and the time needed for drilling of the deployment holes in the glacial ice for the optical modules at depths up to 2700 m. The D-Egg design is utilized for the IceCube Upgrade, the next stage of the IceCube project also known as IceCube-Gen2 Phase 1, where nearly half of the optical sensors to be deployed are D-Eggs. With two 8-inch high-quantum efficiency photomultiplier tubes (PMTs) per module, D-Eggs offer an increased effective area while retaining the successful design of the IceCube digital optical module (DOM). The convolution of the wavelength-dependent effective area and the Cherenkov emission spectrum provides an effective photodetection sensitivity that is 2.8 times larger than that of IceCube DOMs. The signal of each of the two PMTs is digitized using ultra-low-power 14-bit analog-to-digital converters with a sampling frequency of 240 MSPS, enabling a flexible event triggering, as well as seamless and lossless event recording of single-photon signals to multi-photons exceeding 200 photoelectrons within 10 ns. Mass production of D-Eggs has been completed, with 277 out of the 310 D-Eggs produced to be used in the IceCube Upgrade. In this paper, we report the design of the D-Eggs, as well as the sensitivity and the single to multi-photon detection performance of mass-produced D-Eggs measured in a laboratory using the built-in data acquisition system in each D-Egg optical sensor module.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Improved modeling of in-ice particle showers for IceCube event reconstruction

The IceCube Neutrino Observatory relies on an array of photomultiplier tubes to detect Cherenkov light produced by charged particles in the South Pole ice. IceCube data analyses depend on an in-depth characterization of the glacial ice, and on novel approaches in event reconstruction that utilize fast approximations of photoelectron yields. Here, a more accurate model is derived for event reconstruction that better captures our current knowledge of ice optical properties. When evaluated on a Monte Carlo simulation set, the median angular resolution for in-ice particle showers improves by over a factor of three compared to a reconstruction based on a simplified model of the ice. The most substantial improvement is obtained when including effects of birefringence due to the polycrystalline structure of the ice. When evaluated on data classified as particle showers in the high-energy starting events sample, a significantly improved description of the events is observed.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

The LED calibration systems for the mDOM and D-Egg sensor modules of the IceCube Upgrade: Design, production, testing and use in module calibration

The IceCube Neutrino Observatory, instrumenting about 1 km 3 of deep, glacial ice at the geographic South Pole, is due to be enhanced with the IceCube Upgrade. The IceCube Upgrade, to be deployed during the 2025/26 Antarctic summer season, will consist of seven new strings of photosensors, densely embedded near the bottom center of the existing array. Aside from a world-leading sensitivity to neutrino oscillations, a primary goal is the improvement of the calibration of the optical properties of the instrumented ice. This calibration will be applied to the entire archive of IceCube data, improving the angular and energy resolution of the detected neutrino events. For this purpose, the Upgrade strings include a host of new calibration devices. Aside from dedicated calibration modules, several thousand LED flashers have been incorporated into the photosensor modules. We describe the design, production, and testing of these LED flashers before their integration into the sensor modules as well as the use of the LED flashers during lab testing of assembled sensor modules.

Analogue electronic circuits↗

DOE OSTI2026

Evidence for a Spectral Break or Curvature in the Spectrum of Astrophysical Neutrinos from 5 TeV to 10 PeV

Here, we report improved measurements of the all flavor astrophysical neutrino spectrum with IceCube by combining complementary neutrino samples in two independent analyses. Both analyses show evidence of a harder spectrum at energies below ∼ 30 TeV compared to higher energies where the spectrum is well characterized by a power law. The spectrum is better described by a log parabola or a broken power law, the latter being the preferred model. Both, however, reject a single power law over an energy range 5 TeV–10 PeV with a significance > 4⁢𝜎, providing new constraints on properties of cosmic neutrino sources.

Abbasi, R. [Loyola University Chicago] (ORCID:0000↗

DOE OSTI2025

Probing the PeV region in the astrophysical neutrino spectrum using 𝜈 𝜇 from the Southern sky

IceCube has observed a diffuse astrophysical neutrino flux over the energy region from a few TeV to a few PeV. At PeV energies, the spectral shape is not yet well measured due to the low statistics of the data. This analysis probes the gap between 1 and 10 PeV by using high-energy downgoing muon neutrinos. Here, to reject the large atmospheric muon background, two complementary techniques are combined. The first technique selects events with high stochasticity to reject atmospheric muon bundles whose stochastic energy losses are smoothed due to high muon multiplicity. The second technique vetoes atmospheric muons with the IceTop surface array. Using 9 yrs of data, we found two neutrino candidate events in the signal region, consistent with expectation from background, each with relatively high signal probabilities. A joint maximum likelihood estimation is performed using this sample and an independent 9.5-yr sample of tracks to measure the neutrino spectrum. A likelihood ratio test is done to compare the single power-law (SPL) vs SPL+cutoff hypothesis; the SPL+cutoff model is not significantly better than the SPL. High-energy astrophysical objects from four source catalogs are also checked around the direction of the two events. No significant coincidence was found.

Abbasi, R. [Loyola University Chicago] (ORCID:0000↗

DOE OSTI2026

Improved measurements of the TeV-PeV extragalactic neutrino spectrum from joint analyses of IceCube tracks and cascades

The IceCube South Pole Neutrino Observatory has discovered the presence of a diffuse astrophysical neutrino flux at energies of TeV and beyond using neutrino induced muon tracks and cascade events from neutrino interactions. Here, we present two analyses sensitive to neutrino events in the energy range 1 TeV to 10 PeV, using more than 10 years of IceCube data. Both analyses consistently reject a neutrino spectrum following a single power-law with a significance of > 4⁢𝜎 in favor of a broken power law. We describe the methods implemented in the two analyses, the spectral constraints obtained, and the validation of the robustness of the results. Additionally, we report the detection of a muon neutrino in the medium energy starting events sample, or MESE, with an energy of 11.4$^{+2.46}_{−2.53}$ PeV, the highest energy neutrino observed by IceCube to date. The results presented here show insights into the spectral shape of astrophysical neutrinos, which has important implications for inferring their production processes in a multimessenger picture.

astroparticles↗

DOE OSTI2026

Probing neutrino emission at GeV energies from compact binary mergers with the IceCube Neutrino Observatory

The advent of multimessenger astronomy has allowed for new types of source searches by neutrino detectors. Here, we present the results of the search for 0.5–100 GeV astrophysical neutrinos detected with IceCube and emitted from compact binary mergers detected by the LIGO, Virgo, and KAGRA interferometers from their first run of observation (O1) to the end of the first part of the fourth (O4a). An innovative approach is used to lower the energy threshold to 0.5 GeV and to search for an excess of GeV neutrinos in time coincidence with astrophysical transient events. Furthermore, we use a statistical combination of all observations, a binomial test, to search for a subpopulation of neutrino emitters. No significant excess was found from the studied mergers, with a best post-trial 𝑝-value of 40%, and there is currently no hint of a population of GeV neutrino emitters found in the IceCube data (post-trial 𝑝−value=81%).

Abbasi, R. [Loyola University Chicago]↗