Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shadows”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Quantitative investigation of surface structure and interatomic potential with impact-collision ion scattering spectroscopy

Helium ion beam interactions with materials have important implications for magnetic confinement fusion, material modification, and helium ion microscopy. These interactions depend on the precise physics of how helium ions channel into the materials, which can vary greatly based on the local crystalline orientation. Here, we performed a dedicated experiment to investigate helium ion channeling in a well-characterized tungsten single crystal. Time-of-flight impact-collision ion scattering spectroscopy was used to obtain multi-angle maps of the backscattering intensity for 3 keV He + → W(111). We found that the backscattering intensity profile arising from helium ion channeling could be well described by a shadow cone analysis. This analysis revealed that subsurface W atoms as deep as the ninth monolayer contributed to the backscattering intensity profile. Binary collision approximation simulations were performed with MARLOWE to model the experimental maps with sufficient accuracy to allow for quantitative comparisons using reliability factors. These quantitative comparisons were applied to investigate how the W lattice structure and He–W interatomic potential affect the multi-angle maps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Intensity modulated Ir-192 brachytherapy using high-Z 3D printed applicators

Gynecologic cancers are often asymmetric, yet current Ir-192 brachytherapy techniques provide only limited radial modulation of the dose. The shielded solutions investigated here solve this by providing the ability to modulate between highly asymmetric and radially symmetric dose distributions at a given location. To find applicator designs that can modulate between full dose and less than 50% dose, at the dimensions of the urethra, a 2D calculation algorithm was developed to narrow down the search space. Two shielding design types were then further investigated using Monte Carlo and Boltzmann-solver dose calculation algorithms. 3D printing techniques using ISO 10993 certified biocompatible plastics and 3D printable tungsten-loaded plastics were tested. It was also found that shadowing effects set by the shape of the shielding cannot be easily modulated out, hence careful design is required. The shielded applicator designs investigated here, allow for reduction of the dose by over 50% at 5 mm from the applicator surface in desired regions, while also allowing radially symmetric dose with isodose line deviations less than 0.5 mm from circular. The shielding designs were also chosen with treatment delivery time in mind. Treatment times for these shielded designs were found to be less than 1.4 times longer than a 6-channel unshielded cylinder for the equivalent fully symmetric dose distribution. The 2D calculation methods developed here provide a simple way to rapidly evaluate shielding designs, while the 3D printing techniques also allow for devices with novel shapes to be rapidly prototyped. Both TOPAS Monte Carlo and Acuros BV calculations show that significant dose shaping and organ at risk sparing can be achieved without significantly compromising the plan in regions that require the full dose.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Modeling of convective cells, turbulence, and transport induced by a radio-frequency antenna in the tokamak boundary plasma

The edge turbulence model Hermes (Dudson et al 2017 Plasma Phys. Control. Fusion 59 05401) is set up for plasma boundary simulations with an radiofrequency (RF) antenna, using parameters characteristic of a tokamak edge. Cartesian slab geometry is used with thin plate limiters representing the ion cyclotron range of frequency (ICRF) antenna side-wall limiters. Ad-hoc DC electric biasing of the limiters, motivated by calculations with VSim (Nieter et al 2004 J. Comput. Phys. 196 448), represents an induced RF sheath rectified potential in the plasma turbulence model. Flux-driven turbulence simulations demonstrate a realistic distribution of plasma profiles and fluctuations. There is a clear effect of the antenna sheath voltage leading to formation of convective cells; bias-induced convective transport flattens the scrape-off layer density profile and fluctuations penetrate into the shadow region of the limiters as the bias voltage increases. Turbulent transport for impurity ions is inferred by following ion trajectories in the simulated plasma turbulence fields, showing Bohm-like effective diffusion rates. All in all, the model elucidates the key physical phenomena governing the effects of ICRF-induced antenna biasing on the tokamak boundary plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Micro-trench measurements of the net deposition of carbon impurity ions in the DIII-D divertor and the resulting suppression of surface erosion

In this work, we report carbon impurity ion incident angles and deposition rates, along with silicon erosion rates, from measurements of micro-engineered trenches on a silicon surface exposed to L-mode deuterium plasmas at the DIII-D divertor. Post exposure ex-situ analysis determined elemental maps and concentrations, carbon deposition thicknesses, and erosion of silicon surfaces. Carbon deposition profiles on the trench floor showed carbon ion shadowing that was consistent with ERO calculations of average carbon ion angle distributions (IADs) for both polar and azimuthal angles. Measured silicon net erosion rates negatively correlated with the deposited carbon concentration at different locations. Differential erosion of surfaces on two different ion-downstream trench slope structures suggested that carbon deposition rate is affected by the carbon ion incident angle and significantly suppressed the surface erosion. The results suggest the C impurity ion incident angles, determined by the IADs and surface morphology, strongly affect erosion rates as well as the main ion (D, T, He) incident angles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Towards testing the theory of gravity with DESI: summary statistics, model predictions and future simulation requirements

Shortly after its discovery, General Relativity (GR) was applied to predict the behavior of our Universe on the largest scales, and later became the foundation of modern cosmology. Its validity has been verified on a range of scales and environments from the Solar system to merging black holes. However, experimental confirmations of GR on cosmological scales have so far lacked the accuracy one would hope for — its applications on those scales being largely based on extrapolation and its validity there sometimes questioned in the shadow of the discovery of the unexpected cosmic acceleration. Future astronomical instruments surveying the distribution and evolution of galaxies over substantial portions of the observable Universe, such as the Dark Energy Spectroscopic Instrument (DESI), will be able to measure the fingerprints of gravity and their statistical power will allow strong constraints on alternatives to GR. In this paper, based on a set of N-body simulations and mock galaxy catalogs, we study the predictions of a number of traditional and novel summary statistics beyond linear redshift distortions in two well-studied modified gravity models — chameleon f(R) gravity and a braneworld model — and the potential of testing these deviations from GR using DESI. These summary statistics employ a wide array of statistical properties of the galaxy and the underlying dark matter field, including two-point and higher-order statistics, environmental dependence, redshift space distortions and weak lensing. We find that they hold promising power for testing GR to unprecedented precision. The major future challenge is to make realistic, simulation-based mock galaxy catalogs for both GR and alternative models to fully exploit the statistic power of the DESI survey (by matching the volumes and galaxy number densities of the mocks to those in the real survey) and to better understand the impact of key systematic effects. Using these, we identify future simulation and analysis needs for gravity tests using DESI.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Determination of the characteristic magnetic pre-sheath length at divertor surfaces using micro-engineered targets on DiMES at DIII-D

The magnetic pre-sheath (MPS) width, L MPS , is a critical parameter to define the sheath potential, which controls the ion trajectory of low-Z species (D, T, He, and C), as well as the prompt re-deposition of high-Z species. In this work, to determine L MPS , we fabricated micro-trenches (30×30×4 µm) via focused ion beam (FIB) milling on a silicon surface and exposed them to L-mode deuterium plasmas in DIII-D via the Divertor Material Evaluation System (DiMES) removable sample exposure probe. The areal distribution of impurity depositions, mainly consisting of carbon, was measured by energy-dispersive X-ray spectroscopy (EDS) to reveal the deuterium ion shadowing effect on the trench floors. The carbon deposition profiles showed that the erosion was maximized for the azimuthal direction of φ = -40° (referenced to the toroidal magnetic field direction) as well as the polar angle of θ = 80°. A Monte Carlo equation-of-motion model, based on a collisionless MPS, was used to calculate the azimuthal and polar deuterium ion angle distributions (IADs) for a range of L MPS = k × ρ i , where ρ i is the ion gyro radius and k = 0.5-4. Then, gross erosion profiles were calculated by a Monte Carlo micro-patterning and roughness (MPR) code for ion sputtering using as input the calculated azimuthal and polar IADs for each value of k . Good agreement with the experimental C deposition profiles was obtained for the case k = 2.5-3.5. This result is consistent with a previous kinetic modeling prediction of k ~ 3, as well as previous analytical investigations that predicted the L MPS to be several ion gyro radii. A validation of theoretical sheath models supports its applicability to ITER and pilot plant divertors to successfully predict plasma-materials interactions.

Erosion↗
DOE OSTI2022

3D ion gyro-orbit heat load predictions for NSTX-U

Abstract High power tokamaks operate with divertor heat loads capable of destroying the plasma facing components (PFCs). High fidelity heat load predictions are necessary to ascertain the PFC state for design and during operation. Typical heat flux calculations are 2D, time invariant, and assume that power flows directly along the magnetic field lines (the optical approximation). These assumptions neglect the complex 3D geometries employed to protect the PFCs, the time varying nature of the plasma and PFC thermal state, and the helical trajectories of ions with finite Larmor radii (the gyro-orbit approximation). An integrated software framework, the heat flux engineering analysis toolkit (HEAT), was developed to generate time varying optical heat loads applied to real engineering computer aided design (CAD) (Looby et al 2022 Fusion Sci. Technol. 78 10–27). Recently, an ion-gyro orbit module has been added to HEAT. This module calculates the helical trajectories of ions as they gyrate about the magnetic field lines using kinetic theory macro-particles to accelerate the calculation. First, the new gyro-orbit module will be presented. Next, a comparison to existing research is performed. Finally, an analysis of the gyro-orbit heat loads for NSTX-U is presented for diverted discharges using the engineering CAD models utilized for PFC fabrication. Including these gyro-orbit effects can enhance the PFC performance by ‘smearing’ out the magnetic shadows associated with the castellated fish-scaled geometry. Simultaneously, the helical trajectories can degrade performance when they load narrow regions on edges and corners with high heat fluxes. Analysis of the trade-offs between these competing effects is included, and regions for further investigation are identified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Study of the electrothermal and MHD instabilities in exploding cylindrical foil liner

Abstract An experimental and numerical study of the plasma instabilities in an electrical exploding cylindrical Al liner is reported. The Al liner 3 mm in diameter and 10 μ m in thickness is exploded at the Qin-1 facility (450 ns, 400 kA). Various diagnostics, such as multi-frame laser shadowgraphy, an x-ray framing camera, and an x-ray backlighting system are developed. The different imaging systems are sensitive to plasma of different areal densities based on the comparison between the experiments and simulation, which reveal the dynamics of the exploding liner in more detail. The laser shadow images show the low-density plasma (∼1–2 × 10 −4 g cm −2 ) at the edge of the liner, and both the amplitude and wavelength of the plasma instabilities increase over time, which are considered to be magnetohydrodynamic (MHD) instabilities rather than electrothermal instabilities. During the ablation of the liner, quasi-periodic azimuthally correlated striations are directly observed in extreme ultraviolet (EUV) self-emission images. Meanwhile, the vertical filaments, which are electrothermal instabilities for plasma under the condition of ∂ η /∂ T < 0, are also observed in EUV self-emission images. The x-ray backlighting images of the exploding liner are obtained by placing an X -pinch load on the current-return path to serve as an x-ray point source (∼1 ns, ∼10 μ m). The x-ray backlighting results show the behavior of the high-density plasma (∼1.89 × 10 −3 g cm −2 ), which includes the transition from electrothermal to MHD instabilities. Finally, we realized a 2D MHD simulation of the exploding liner under experimental conditions, which shows good agreement with the results of the experimental perturbation.

Physics↗
DOE OSTI2025

Multi-X-point radiation and its dynamics in up/down asymmetry in plasma detachment regimes in Wendelstein 7-X

In Wendelstein 7-X (W7-X), stable plasma detachment in ECRH hydrogen plasmas is routinely achieved with the standard magnetic configuration, containing five magnetic islands in the scrape-off layer (SOL). Plasma detachment induced by intrinsic carbon impurities is characterized by a high radiation fraction (ƒ rad ) with significantly reduced divertor heat load and particle flux. In this work, we present the radiation dynamics during the detached plasma (DP) phase, such as the shift of radiation zones in radial and poloidal directions towards the X-points as ƒ rad increases (∼0.6–0.9), as well as the penetration of the radiation layer into the confinement region in the deep DP phase with ƒ rad approaching unity. In particular, a structure of multi-X-point radiation (multi-XPR) with an up/down asymmetry in the DP phase is highlighted, which is revealed by bolometer tomography and is further confirmed by video diagnostics. The multi-XPR structure forms helical 3D bands aligned with W7-X field periodicity. Field reversal experiments show that the brightest XPR displaces between the upper and lower SOL regions, suggesting the E× B drift effects. This paper presents the multi-XPR structure observed in the W7-X plasma in detail for the first time. A simplified model considering the influence of the poloidal E× B drift ( V d ) on the impurity flow in the SOL shows that the downstream drift toward the target or target shadow region ( V d < 0) decreases the impurity content, while the upstream drift toward the LCFS ( V d > 0) increases the impurity content. On this basis, the poloidal drift potentially leads to an up/down asymmetry of impurity density in the SOL despite the symmetry magnetic topology. The dynamics of the up/down asymmetry in the multi-XPR structure is also related to the magnitude V d /D (normalized to the impurity diffusivity), with an additional effect owing to the radial inward shift of the emission zone. These results provide new insights into impurity-induced detachment dynamics and provide a basis for improving 3D modeling of impurity transport, which typically does not consider drifts.

E × B drift↗
DOE OSTI2025

SHADOW4: the popular ray tracing revived for evolving synchrotron sources in fourth-generation storage rings

We present SHADOW4, a new version of the popular ray tracing code. The SHADOW kernel has been completely rewritten in Python applying modern concepts of software engineering. A new user interface is available in the OASYS ecosystem. The new tool has been designed and implemented preparing the future needs both in computing (cloud computing, AI integration) and in the transit to fourth generation sources and beyond.

Sanchez del Rio, Manuel↗
DOE OSTI2024

Not even 6 dB: Gaussian quantum illumination in thermal background

Abstract In analyses of target detection with Gaussian state transmitters in a thermal background, the thermal occupation is taken to depend on the target reflectivity in a way which simplifies the analysis of the symmetric quantum hypothesis testing problem. However, this assumption precludes comparison of target detection performance between an arbitrary transmitter and a vacuum state transmitter, i.e. ‘detection without illumination’, which is relevant in a bright thermal background because a target can be detected by its optical shadow or some other perturbation of the background. Using a target-agnostic thermal environment leads to the result that the oft-claimed 6 dB possible reduction in the quantum Chernoff exponent for a two-mode squeezed vacuum transmitter over a coherent state transmitter in high-occupation thermal background is an unachievable limiting value, only occurring in a limit in which the target detection problem is ill-posed. Further analyzing quantum illumination in a target-agnostic thermal environment shows that a weak single-mode squeezed transmitter performs worse than ‘no illumination’, which is explained by the noise-increasing property of reflected low-intensity squeezed light.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Characterization and thermometry of dissipatively stabilized steady states

In this work we study the properties of dissipatively stabilized steady states of noisy quantum algorithms, exploring the extent to which they can be well approximated as thermal distributions, and proposing methods to extract the effective temperature T. We study an algorithm called the relaxational quantum eigensolver (RQE), which is one of a family of algorithms that attempt to find ground states and balance error in noisy quantum devices. In RQE, we weakly couple a second register of auxiliary ‘shadow’ qubits to the primary system in Trotterized evolution, thus engineering an approximate zero-temperature bath by periodically resetting the auxiliary qubits during the algorithm’s runtime. Balancing the infinite temperature bath of random gate error, RQE returns states with an average energy equal to a constant fraction of the ground state. We probe the steady states of this algorithm for a range of base error rates, using several methods for estimating both T and deviations from thermal behavior. In particular, we both confirm that the steady states of these systems are often well-approximated by thermal distributions, and show that the same resources used for cooling can be adopted for thermometry, yielding a fairly reliable measure of the temperature. These methods could be readily implemented in near-term quantum hardware, and for stabilizing and probing Hamiltonians where simulating approximate thermal states is hard for classical computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Learning to classify quantum phases of matter with a few measurements

We study the identification of quantum phases of matter, at zero temperature, when only part of the phase diagram is known in advance. Following a supervised learning approach, we show how to use our previous knowledge to construct an observable capable of classifying the phase even in the unknown region. By using a combination of classical and quantum techniques, such as tensor networks, kernel methods, generalization bounds, quantum algorithms, and shadow estimators, we show that, in some cases, the certification of new ground states can be obtained with a polynomial number of measurements. An important application of our findings is the classification of the phases of matter obtained in quantum simulators, e.g. cold atom experiments, capable of efficiently preparing ground states of complex many-particle systems and applying simple measurements, e.g. single qubit measurements, but unable to perform a universal set of gates.

quantum machine learning↗
DOE OSTI2023

Regional infrasonic observations from surface explosions – influence of atmospheric variations and realistic terrain

SUMMARY A pair of 1-ton, conventional surface explosions were conducted at the Nevada National Security Site in the fall of 2020 producing seismoacoustic signatures observable hundreds of kilometres from the source location. Regional infrasonic observations include tropospheric ducting at large distances to the south, a wide stratospheric waveguide with signals observed more than 700 km to the east, and anomalous arrivals in the stratospheric shadow zone. Notable differences in propagation between the events are identified despite the explosions being conducted just two days apart due to a sharp temporal shift in the tropospheric winds as well as structural changes in the stratospheric winds. Propagation simulations of the two events have been completed using a combination of ray tracing and parabolic equation (PE) methods. Simulations have been conducted to quantify the impact of the temporal variations in the atmosphere as well as the influence of terrain on propagation. Temporal variations in reflection locations are found to produce notable changes in downrange propagation due to spatially varying terrain features. Finite frequency effects modelled by the PE are found to predict ensonification not included in corresponding 2D ray tracing simulations. Notable variations in predicted signal amplitude are found due to focusing by along-path and cross-path terrain gradients; though, the later of these is only modelled using fully 3D ray tracing analysis.

58 GEOSCIENCES↗
DOE OSTI2021

Morphostatistical characterization of the spatial galaxy distribution through Gibbs point processes

ABSTRACT This paper proposes a morphostatistical characterization of the galaxy distribution through spatial statistical modelling based on inhomogeneous Gibbs point processes. The galaxy distribution is supposed to exhibit two components. The first one is related to the major geometrical features exhibited by the observed galaxy field, here, its corresponding filamentary pattern. The second one is related to the interactions exhibited by the galaxies. Gibbs point processes are statistical models able to integrate these two aspects in a probability density, controlled by some parameters. Several such models are fitted to real observational data via the ABC shadow algorithm. This algorithm provides simultaneous parameter estimation and posterior-based inference, hence allowing the derivation of the statistical significance of the obtained results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Observation of Coherent 𝜙⁡(1020) Meson Photoproduction in Ultraperipheral PbPb Collisions at $\sqrt{𝑠_{NN}}$ = 5.36 TeV

The first observation of coherent 𝜙⁡(1020) meson photoproduction off heavy nuclei is presented using ultraperipheral lead-lead collisions at a center-of-mass energy per nucleon pair of 5.36 TeV. The data were collected by the CMS experiment and correspond to an integrated luminosity of 1.62 μ⁢𝑏 −1 . The 𝜙⁡(1020) meson signals are reconstructed via the 𝐾 + ⁢𝐾 − decay channel. The production cross section is presented as a function of the 𝜙⁡(1020) meson rapidity in the range 0.3 < |𝑦| < 1.0, probing gluons that carry a fraction of the nucleon momentum (𝑥) around 10 −4 . The observed cross section exhibits little dependence on rapidity and is significantly suppressed, by a factor of ∼5, compared to a baseline model that treats a nucleus as a collection of free nucleons. Theoretical models that incorporate the nuclear shadowing effect generally provide a better description of the 𝜙⁡(1020) data than those incorporating gluon saturation. This study establishes a powerful new tool for exploring nuclear effects and nuclear gluonic structure in the small-𝑥 regime at a unique energy scale bridging the perturbative and nonperturbative quantum chromodynamics domains.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Correlated purification for restoring 𝑁-representability in quantum simulation

Experimentally measured reduced density matrices (RDMs) often violate constraints that ensure they represent N-electron states—known as N-representability conditions—because of statistical and hardware noise. In this work, we present a correlated purification framework based on semidefinite programming to restore the accuracy of a noisy, unphysical two-electron RDM (2-RDM). The method performs a bi-objective optimization that minimizes both the many-electron energy and the nuclear norm of the correction to the measured 2-RDM. The nuclear norm, often employed in matrix completion, promotes low-rank corrections, while the energy term acts as a regularization term that can improve the purity of the ground state. While the method is particularly effective for ground states, it can also be applied to excited and nonstationary states by decreasing the weight of the energy relative to the error norm. In an application to fermionic shadow tomography of large hydrogen chains, correlated purification yields substantial reductions in both energy and 2-RDM error, achieving chemical accuracy across dissociation curves. This framework provides a robust strategy for tomography in many-body quantum simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Exponentially Reduced Circuit Depths Using Trotter Error Mitigation

Product formulas are a popular class of digital quantum simulation algorithms due to their conceptual simplicity, low overhead, and performance, which often exceeds theoretical expectations. Recently, Richardson extrapolation and polynomial interpolation have been proposed to mitigate the Trotter error incurred by the use of these formulas. This work provides a rigorous, general analysis of these techniques for computing time-evolved observables, simplifying the interpolation algorithm in the process, and shows that extrapolation generically improves the performance of product formulas for this task. We demonstrate that, to achieve error 𝜖 in a simulation of time 𝑇 using a 𝑝 ⁢th-order product formula with extrapolation, circuit depths of 𝑂⁡(𝑇 1+1/𝑝 ⁢polylog (1/𝜖)) are sufficient—an exponential improvement in the precision over product formulas alone. Furthermore, we prove that these algorithms achieve commutator scaling, and improve the 𝑇 complexity for the interpolation algorithm. By relaxing the requirement of performing exact Chebyshev interpolation, our simplified algorithm eliminates the need for fractional implementations of Trotter steps, reducing computational overhead. Finally, we show these techniques can be combined with the classical shadows method to estimate many time-evolved local observables. Taken together, our findings provide the strongest evidence yet for the utility of Trotter error-mitigation techniques in algorithmic applications.

quantum algorithms & computation↗