SEARCH · Engineering Papers

Results for “Technique Development”

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 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Light in the dark forest. Part I. An efficient optimal estimator for 3D Lyman-alpha forest power spectrum

The highly anisotropic nature of the Lyman-alpha (Lyα) forest data introduces a complex survey window function that complicates the measurement of the three-dimensional power spectrum ( P 3D ). In this paper, we present the first fully optimal estimator for P 3D , which exactly deconvolves the survey window function and marginalizes contaminated modes that distort the power spectrum. Our approach adapts optimal estimator techniques developed for the 2D cosmic microwave background data to the 3D case. To achieve computational feasibility, we employ the conjugate gradient method and implement the P 3 M formalism to handle large-scale and small-scale operations separately and efficiently. We validate our estimator using Monte Carlo mocks and Gaussian simulations, demonstrating its accuracy and computational efficiency. We confirm that mode marginalization eliminates distortions arising from quasar continuum errors and delivers robust power spectrum estimation, though it also inflates errors at large scales. This first implementation works in the flat-sky case; we discuss the remaining steps needed to generalize it to the curved-sky case. This formalism offers a foundation for the Lyα forest P 3D measurements and a new path toward cosmological constraints from the Lyα forest data.

Lyman alpha forest↗

DOE OSTI2023

Time-dependent SOLPS-ITER simulations of the tokamak plasma boundary for model predictive control using SINDy *

Abstract Time-dependent SOLPS-ITER simulations have been used to identify reduced models with the sparse identification of nonlinear dynamics (SINDy) method and develop model-predictive control of the boundary plasma state using main ion gas puff actuation. A series of gas actuation sequences are input into SOLPS-ITER to produce a dynamic response in upstream and divertor plasma quantities. The SINDy method is applied to identify reduced linear and nonlinear models for the electron density at the outboard midplane n e , s e p O M P and the electron temperature at the outer divertor T e , s e p d i v . Note that T e , s e p d i v is not necessarily the peak value of T e along the divertor. The identified reduced models are interpretable by construction (i.e. not black box), and have the form of coupled ordinary differential equations. Despite significant noise in T e , s e p d i v , the reduced models can be used to predict the response over a range of actuation levels to a maximum deviation of 0.5% in n e , s e p O M P and 5%–10% in T e , s e p d i v for the cases considered. Model retraining using time history data triggered by a preset error threshold is also demonstrated. A model predictive control strategy for nonlinear models is developed and used to perform feedback control of a SOLPS-ITER simulation to produce a setpoint trajectory in n e , s e p O M P using the integrated plasma simulator framework. The developed techniques are general and can be applied to time-dependent data from other boundary simulations or experimental data. Ongoing work is extending the approach to model identification and control for divertor detachment, which will present transient nonlinear behavior from impurity seeding, including realistic latency and synthetic diagnostic signals derived from the full SOLPS-ITER output.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Diagnosing inertial confinement fusion ignition

Abstract Fusion ignition by inertial confinement requires compression and heating of the fusion fuel to temperatures in excess of 5 keV and densities exceeding hundreds of g/cc. In August 2021, this scientific milestone was surpassed at the National Ignition Facility (NIF), when the Lawson criterion for ignition was exceeded generating 1.37MJ of fusion energy (Abu-Shawareb et al 2022 Phys. Rev. Lett. 129 075001), and then in December 2022 target gain >1 was realized with the production of 3.1MJ of fusion energy from a target driven by 2.0MJ of laser energy (Abu-Shawareb et al 2024 Phys. Rev. Lett. 132 065102). At the NIF, inertial confinement fusion research primarily uses a laser indirect drive in which the fusion capsule is surrounded by a high-Z enclosure (‘hohlraum’) used to convert the directed laser energy into a symmetric x-ray drive on the capsule. Precise measurements of the plasma conditions, x-rays, γ -rays and neutrons produced are key to understanding the pathway to higher performance. This paper discusses the diagnostics and measurement techniques developed to understand these experiments, focusing on three main topics: (1) key diagnostic developments for achieving igniting plasmas, (2) novel signatures related to thermonuclear burn and (3) advances to diagnostic capabilities in the igniting regime with a perspective toward developments for intertial fusion energy.

Moore, A. S. (ORCID:0000000221258882)↗

DOE OSTI2025

Global stellarator coil optimization with quadratic constraints and objectives

Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly $10$ 2 x faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Gauges, loops, and polynomials for partition functions of graphical models

Graphical models represent multivariate and generally not normalized probability distributions. Computing the normalization factor, called the partition function, is the main inference challenge relevant to multiple statistical and optimization applications. The problem is #P-hard that is of an exponential complexity with respect to the number of variables. Here, aimed at approximating the partition function, we consider multi-graph models where binary variables and multivariable factors are associated with edges and nodes, respectively, of an undirected multi-graph. We suggest a new methodology for analysis and computations that combines the Gauge function technique from Chertkov and Chernyak with the technique developed in Anari and Oveis Gharan 2017 arXiv:1702.02937; Gurvits 2011 arXiv:1106.2844; Straszak and Vishnoi 2017 55th Annual Allerton Conf. on Communication, Control, and Computing, based on the recent progress in the field of real stable polynomials. We show that the Gauge function, representing a single-out term in a finite sum expression for the partition function which achieves extremum at the so-called belief-propagation gauge, has a natural polynomial representation in terms of gauges/variables associated with edges of the multi-graph. Moreover, Gauge function can be used to recover the partition function through a sequence of transformations allowing appealing algebraic and graphical interpretations. Algebraically, one step in the sequence consists of the application of a differential operator over gauges associated with an edge. Graphically, the sequence is interpreted as a repetitive elimination/contraction of edges resulting in multi-graph models on decreasing in size (number of edges) graphs with the same partition function as in the original multi-graph model. Even though the complexity of computing factors in the sequence of the derived multi-graph models and respective Gauge functions grow exponentially with the number of eliminated edges, polynomials associated with the new factors remain bi-stable if the original factors have this property. Moreover, we show that BP estimations in the sequence do not decrease, each low-bounding the partition function.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Integrating Visualization Applications, such as ParaView, into HEP Software Frameworks for In-situ Event Displays

ParaView is a high performance visualization application not widely used in High Energy Physics (HEP). It is a long standing open source project led by Kitware and involves several Department of Energy (DOE) and Department of Defense (DOD) laboratories. Futhermore, it has been adopted by many DOE supercomputing centers and other sites. ParaView is unique in speed and efficiency by using state-of-the-art techniques developed by the academic visualization community that are often not found in applications written by the HEP community. In-situ visualization of events, where event details are visualized during processing/analysis, is a common task for experiment software frameworks. Kitware supplies Catalyst, a library that enables scientific software to serve visualization objects to client ParaView viewers yielding a real-time event display. Connecting ParaView to the Fermilab art framework will be described and the capabilities it brings discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Industrial deposition of wavelength-shifting films for liquid argon photon detection systems

The Deep Underground Neutrino Experiment (DUNE) Phase-II Far Detector is considering an approximately 2000 m 2 photon detection system to achieve a target mean light yield of 180 PE/MeV. Meeting this requirement demands scalable, cost-effective, and high-quality wavelength-shifter (WLS) coatings capable of converting 127 nm liquid-argon scintillation light into visible photons with controlled and reproducible optical performance. We report on the successful realization of an industrial physical vapor deposition (PVD) process forp-terphenyl (pTP) coatings, adapted from vacuum deposition techniques developed for OLED display manufacturing, to produce uniform WLS layers on large-area inorganic substrates, a task traditionally challenged by adhesion and uniformity issues at organic-inorganic interfaces. Surface characterization by profilometry and spectroscopic measurements demonstrates edge-region thickness variation below 10% and emission spectra consistent with high-quality pTP reference samples. The industrial process demonstrates reproducibility, scalability, and significantly reduced production time compared to laboratory-based methods, while maintaining optical characteristics consistent with established pTP reference samples. These results establish a viable pathway for mass production of high-performance pTP coatings for DUNE FD3 and future neutrino experiments, from a coating manufacturing and process standpoint. Detector-level performance validation, including quantitative VUV conversion efficiency measurements at 127 nm, is identified as future work.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Outlook for artificial intelligence and machine learning at the NSLS-II

Abstract We describe the current and future plans for using artificial intelligence and machine learning (AI/ML) methods at the National Synchrotron Light Source II (NSLS-II), a scientific user facility at the Brookhaven National Laboratory. We discuss the opportunity for using the AI/ML tools and techniques developed in the data and computational science areas to greatly improve the scientific output of large scale experimental user facilities. We describe our current and future plans in areas including from detecting and recovering from faults, optimizing the source and instrument configurations, streamlining the pipeline from measurement to insight, through data acquisition, processing, analysis. The overall strategy and direction of the NSLS-II facility in relation to AI/ML is presented.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Modelling the stellar halo with RR-Lyrae stars

ABSTRACT A seven-parameter distribution function (DF) is fitted to $20\, 000$ RR-Lyrae stars for which only astrometric data are available. The observational data are predicted by the DF in conjunction with the gravitational potential of a self-consistent model Galaxy defined by DFs for the dark halo, the bulge, and a four-component disc. Tests of the technique developed to deal with missing line-of-sight velocities show that adding such velocities tightens constraints on the DF only slightly. The recovered model of the RR-Lyrae population confirms that the population is flattened and has a strongly radially biased velocity distribution. At large radii, its density profile tends to ρ ∼ r−4.5 but no power law provides a good fit inside the solar sphere. The model is shown to provide an excellent fit to the data for stars brighter than r = 16.5 but at certain longitudes it predicts too few faint stars at Galactocentric radii $\sim 20\, \mathrm{kpc}$, possibly signalling that the halo is not axisymmetric. The DF is used to predict the velocity distribution of BHB stars for which space velocities are available. The z components are predicted successfully but too much anisotropy in the vRvϕ plane is expected.

Li, Chengdong↗

DOE OSTI2022

Prospects for future studies using deep imaging: analysis of individual Galactic cirrus filaments

The presence of Galactic cirrus is an obstacle for studying both faint objects in our Galaxy and low surface brightness extragalactic structures. With the aim of studying individual cirrus filaments in Sloan Digital Sky Survey (SDSS) Stripe 82 data, we develop techniques based on machine learning and neural networks that allow one to isolate filaments from foreground and background sources in the entirety of Stripe 82 with a precision similar to that of the human expert. Our photometric study of individual filaments indicates that only those brighter than 26 mag arcsec −2 in the SDSS r band are likely to be identified in SDSS Stripe 82 data by their distinctive colours in the optical bands. We also show a significant impact of data processing (e.g. flat-fielding, masking of bright stars, and sky subtraction) on colour estimation. Analysing the distribution of filaments’ colours with the help of mock simulations, we conclude that most filaments have colours in the following ranges: 0.55 ≤g − r ≤ 0.73 and 0.01 ≤ r − i ≤ 0.33. Our work provides a useful framework for an analysis of all types of low surface brightness features (cirri, tidal tails, stellar streams, etc.) in existing and future deep optical surveys. For practical purposes, we provide the catalogue of dust filaments.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

How negative feedback and the ambient environment limit the influence of recombination in common envelope evolution

ABSTRACT We perform 3D hydrodynamical simulations to study recombination and ionization during the common envelope (CE) phase of binary evolution, and develop techniques to track the ionic transitions in time and space. We simulate the interaction of a $2\, \mathrm{M_\odot }$ red giant branch primary and a $1\, \mathrm{M_\odot }$ companion modelled as a particle. We compare a run employing a tabulated equation of state (EOS) that accounts for ionization and recombination, with a run employing an ideal gas EOS. During the first half of the simulations, ∼15 per cent more mass is unbound in the tabulated EOS run due to the release of recombination energy, but by simulation end the difference has become negligible. We explain this as being a consequence of (i) the tabulated EOS run experiences a shallower inspiral and hence smaller orbital energy release at late times because recombination energy release expands the envelope and reduces drag, and (ii) collision and mixing between expanding envelope gas, ejecta and circumstellar ambient gas assists in unbinding the envelope, but does so less efficiently in the tabulated EOS run where some of the energy transferred to bound envelope gas is used for ionization. The rate of mass unbinding is approximately constant in the last half of the simulations and the orbital separation steadily decreases at late times. A simple linear extrapolation predicts a CE phase duration of ${\sim}2\, {\rm yr}$, after which the envelope would be unbound.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Probabilistic inference of the structure and orbit of Milky Way satellites with semi-analytic modelling

Semi-analytic modelling furnishes an efficient avenue for characterizing dark matter haloes associated with satellites of Milky Way-like systems, as it easily accounts for uncertainties arising from halo-to-halo variance, the orbital disruption of satellites, baryonic feedback, and the stellar-to-halo mass (SMHM) relation. We use the SatGen semi-analytic satellite generator, which incorporates both empirical models of the galaxy–halo connection as well as analytic prescriptions for the orbital evolution of these satellites after accretion onto a host to create large samples of Milky Way-like systems and their satellites. By selecting satellites in the sample that match observed properties of a particular dwarf galaxy, we can infer arbitrary properties of the satellite galaxy within the cold dark matter paradigm. For the Milky Way’s classical dwarfs, we provide inferred values (with associated uncertainties) for the maximum circular velocity v max and the radius r max at which it occurs, varying over two choices of baryonic feedback model and two prescriptions for the SMHM relation. While simple empirical scaling relations can recover the median inferred value for v max and r max , this approach provides realistic correlated uncertainties and aids interpretability. We also demonstrate how the internal properties of a satellite’s dark matter profile correlate with its orbit, and we show that it is difficult to reproduce observations of the Fornax dwarf without strong baryonic feedback. Furthermore, the technique developed in this work is flexible in its application of observational data and can leverage arbitrary information about the satellite galaxies to make inferences about their dark matter haloes and population statistics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Ultrafast dynamics of spin and orbital correlations in quantum materials: an energy- and momentum-resolved perspective

Many remarkable properties of quantum materials emerge from states with intricate coupling between the charge, spin and orbital degrees of freedom. Ultrafast photo-excitations of these materials hold great promise for understanding and controlling the properties of these states. Here we introduce time-resolved resonant inelastic X-ray scattering (tr-RIXS) as a means of measuring the charge, spin and orbital excitations out of equilibrium. These excitations encode the correlations and interactions that determine the detailed properties of the states generated. After outlining the basic principles and instrumentation of tr-RIXS, we review our first observations of transient antiferromagnetic correlations in quasi two dimensions in a photo-excited Mott insulator and present possible future routes of this fast-developing technique. The increasing number of X-ray free electron laser facilities not only enables tackling long-standing fundamental scientific problems, but also promises to unleash novel inelastic X-ray scattering spectroscopies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Color glass condensate meets high twist expansion

We establish the correspondence between two well-known frameworks for quantum chromodynamics (QCD) multiple scattering in nuclear media: the color glass condensate (CGC) and the high-twist (HT) expansion formalism. We argue that a consistent matching between both frameworks, in their common domain of validity, is achieved by incorporating the subeikonal longitudinal momentum phase in the CGC formalism, which mediates the transition between coherent and incoherent scattering. We perform a detailed calculation and analysis of direct photon production in proton-nucleus scattering as a concrete example to establish the matching between HT and CGC up to twist-4, including initial- and final-state interactions, as well as their interferences. The techniques developed in this work can be adapted to other processes in electron-nucleus and proton-nucleus collisions, and they provide a potential avenue for a unified picture of dilute-dense dynamics in nuclear media.

quantum chromodynamics↗

DOE OSTI2026

Exploring the keV-scale physics potential of CUORE

We present the analysis techniques developed to explore the keV-scale energy region of the Cryogenic Underground Observatory for Rare Events (CUORE) experiment, based on more than 2 metric ton yr of data collected over five years. By prioritizing a stricter selection over a larger exposure, we are able to optimize data selection for thresholds at 10 keV and 3 keV with 691 kg yr and 11 kg yr of data, respectively. We study how the performance varies among the 988-detector array with different detector characteristics and data-taking conditions. We achieve an average baseline resolution of 2.54 ±0.14 keV FWHM and 1.18 ±0.02 keV FWHM for the data selection at 10 keV and 3 keV, respectively. The analysis methods employed reduce the overall background by about an order of magnitude, reaching 2.06±0.05 counts/(keV kg days) and 16±2 counts/(keV kg days) at the thresholds of 10 keV and 3 keV. We evaluate for the first time the near-threshold reconstruction efficiencies of the CUORE experiment, and find these to be 50 ±2% and 26 ±4% at 10 keV and 3 keV, respectively. This analysis provides crucial insights into rare decay studies, new physics searches, and keV-scale background modeling with CUORE. We demonstrate that ton-scale cryogenic calorimeters can operate across a wide energy range, from keV to MeV, establishing their scalability as versatile detectors for rare event and dark matter physics. These findings also inform the optimization of future large mass cryogenic calorimeters to enhance the sensitivity to low-energy phenomena.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2019

Probing Stark-induced nonlinear phase variation with opto-optical modulation

We extend the recently developed technique of opto-optical modulation (OOM) to probe state-resolved ac-Stark-induced phase variations of a coherently excited ensemble of helium atoms. In a joint experimental and theoretical study, we find that the spatial redirection of the resonant emission from the OOM process is different for the low-lying 1s2p state as compared with the higher-lying Rydberg states, and that this redirection can be controlled through the spatial characteristics of the infrared (IR) probe beam. In particular, we observe that the intensity dependence of the IR-induced Stark phase on the 1s2p emission is nonlinear, and that the phase accumulation changes sign for moderate intensities. Our results suggest that OOM, combined with precise experimental shaping of the probe beam, could allow future measurements of Stark-induced phase shifts of excited states.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2019

Direct detection of multiple backward volume modes in yttrium iron garnet at micron scale wavelengths

Using a set of wave-vector-specific multielement antennas, we have characterized the dispersion of spin waves in an yttrium iron garnet film at submicron lengths and resolved the dispersion relations of multiple backward volume modes, particularly in the region of their minima The techniques developed now facilitate the characterization of spin waves at length scales limited only by available lithography and at a spectral resolution that generally exceeds that of Brillouin scattering. The data obtained are in excellent agreement with theoretical predictions based on a model Hamiltonian.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Connection between asymptotic normalization coefficients and resonance widths of mirror states

Resonance widths are important characteristics of unstable nuclear states playing a crucial role in nuclear structure, nuclear reactions, and nuclear astrophysics. However, often it is difficult or impossible to determine the resonance widths from experiments, especially when required nuclei are unstable or the resonance widths are too small. In this paper, I suggest an indirect method of determination of resonance widths using transfer reactions populating bound states. From this reaction, one can extract the amplitudes of the tails of the bound-state wave functions. These amplitudes are called the asymptotic normalization coefficients (ANCs). The ANC method became one of the benchmarks of the Cyclotron Institute. Once the ANC is experimentally determined, one can find the resonance width of the mirror nuclei using the technique developed in this paper. Mirror nuclei are the ones in which the total number of nucleons is the same but the number of the protons of one of the nuclei is equal to the number of the neutrons of another nucleus. In this paper, 5 different examples of the application of the developed method are presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗