Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Averages”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quantifying and simulating the weather forecast uncertainty for advanced building control

Weather forecast uncertainty is unavoidable despite technological advancements. Accurately quantifying and modelling this uncertainty is essential for developing and comparing advanced building controllers. In this study, we present a structured approach using a first-order autoregressive model (AR(1)) to model uncertainty in ambient temperature and global solar irradiation (GHI) forecasts. We analyzed weather data from four cities and employed Jensen–Shannon divergence (JSD) to evaluate the similarity between synthetic and actual forecast errors. The average JSD values for temperature are 0.027 (Berkeley), 0.021 (Leuven), 0.018 (Berlin), and 0.008 (Oslo), and for GHI, the average JSD values are 0.016 (Berkeley), 0.058 (Leuven), and 0.013 (Berlin). The low JSD values indicate a high similarity between the synthetic and real forecast error distributions. Further, our approach successfully generates synthetic weather forecasts that mirror the statistical properties of actual forecasts. The implementation of our method for uncertain forecast generation is being added to the BOPTEST framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mapping dynamical ejecta and disk masses from numerical relativity simulations of neutron star mergers

Abstract We present fitting formulae for the dynamical ejecta properties and remnant disk masses from the largest to date sample of numerical relativity simulations. The considered data include some of the latest simulations with microphysical nuclear equations of state (EOS) and neutrino transport as well as other results with polytropic EOS available in the literature. Our analysis indicates that the broad features of the dynamical ejecta and disk properties can be captured by fitting expressions, that depend on mass ratio and reduced tidal parameter. The comparative analysis of literature data shows that microphysics and neutrino absorption have a significant impact on the dynamical ejecta properties. Microphysical nuclear EOS lead to average velocities smaller than polytropic EOS, while including neutrino absorption results in larger average ejecta masses and electron fractions. Hence, microphysics and neutrino transport are necessary to obtain quantitative models of the ejecta in terms of the binary parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Midisuperspace foam and the cosmological constant

Abstract Wheeler’s conjectured ‘spacetime foam’—large quantum fluctuations of spacetime at the Planck scale—could have important implications for quantum gravity, perhaps even explaining why the cosmological constant seems so small. Here I explore this problem in a midisuperspace model consisting of metrics with local spherical symmetry. Classically, an infinite class of ‘foamy’ initial data can be constructed, in which cancellations between expanding and contracting regions lead to a small average expansion even if Λ is large. Quantum mechanically, the model admits corresponding stationary states, for which the probability current is also nearly zero. These states appear to describe a self-reproducing spacetime foam with very small average expansion, effectively hiding the cosmological constant.

Astronomy & Astrophysics↗
DOE OSTI2025

Microfluidic cells for the 1–10 2 MPa pressure range

Thin membrane-delimited fluid cells supporting up to 1 at (0.1 MPa) of pressure are well known and commercially available for use in vacuum chambers of electron, photon, or various particle beam microscopies or spectroscopies. Hereby, we report on the development of fluid cells capable of working at 1–10 MPa, extending the analysis domain for investigating chemical, biochemical, or physical processes at pressures of interest in chemical synthesis, underwater biochemistry studies or underground geology. We explored ways to optimize cell membranes to better resist pressure beyond simply increasing the thickness or decreasing the size of the membranes, using finite element analysis and experimental validation via membrane bulging experiments and failure statistics. Fluid cell prototypes were fabricated using ∼75 nm-thick SiN x membranes, engineered to withstand 4.7 MPa (average value), compared to regular (un-engineered) membranes withstanding only 3.4 MPa (average value). The fluid cell prototypes include eight microchannels for feeding/evacuating the fluids and applying pressure into micro-reaction chambers, two electrodes for electrochemical or conduction measurements in the sample, and a possible pressure or temperature sensor, customizable for specific experiments.

hi-pressure↗
DOE OSTI2023

Dosimetric delivery validation of dynamically collimated pencil beam scanning proton therapy

Pencil beam scanning (PBS) proton therapy target dose conformity can be improved with energy layer-specific collimation. One such collimator is the dynamic collimation system (DCS), which consists of four nickel trimmer blades that intercept the scanning beam as it approaches the lateral extent of the target. While the dosimetric benefits of the DCS have been demonstrated through computational treatment planning studies, there has yet to be experimental verification of these benefits for composite multi-energy layer fields. The objective of this work is to dosimetrically characterize and experimentally validate the delivery of dynamically collimated proton therapy with the DCS equipped to a clinical PBS system. Approach. Optimized single field, uniform dose treatment plans for 3 × 3 × 3 cm 3 target volumes were generated using Monte Carlo dose calculations with depths ranging from 5 to 15 cm, trimmer-to-surface distances ranging from 5 to 18.15 cm, with and without a 4 cm thick polyethylene range shifter. Treatment plans were then delivered to a water phantom using a prototype DCS and an IBA dedicated nozzle system and measured with a Zebra multilayer ionization chamber, a MatriXX PT ionization chamber array, and Gafchromic™ EBT3 film. Main results. For measurements made within the SOBPs, average 2D gamma pass rates exceeded 98.5% for the MatriXX PT and 96.5% for film at the 2%/2 mm criterion across all measured uncollimated and collimated plans, respectively. For verification of the penumbra width reduction with collimation, film agreed with Monte Carlo with differences within 0.3 mm on average compared to 0.9 mm for the MatriXX PT. Significance. We have experimentally verified the delivery of DCS-collimated fields using a clinical PBS system and commonly available dosimeters and have also identified potential weaknesses for dosimeters subject to steep dose gradients.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Numerical simulation of the hot-tail runaway electron production mechanism using CQL3D and comparison with Smith–Verwichte analytical model

Abstract The hot-tail mechanism of runaway electron (RE) production (Harvey et al 2000 Phys. Plasmas 7 4590) is the primary source of RE in the case of rapidly cooling tokamak plasma. Quantifying this mechanism is very important as it can provide most of the post-thermal-quench (TQ) current, or a seed current for the secondary source of RE through the avalanche mechanism. An analytic model which omits pitch-angle scattering is often used in literature for estimating the hot-tail RE density (Smith and Verwichte 2008 Phys. Plasmas 15 072502). In the present study, we use the CQL3D bounce-averaged Fokker–Planck code (Harvey and McCoy 1992 Proc. IAEA Technical Committee Meeting on Advances in Simulation and Modeling of Thermonuclear Plasmas p 527) to test the limits of validity of the model. In particular, we examine the cases of Z = 1 and Z = 18 ions, for sets of different initial temperature, density, electric field and the characteristic time of temperature decay. We show that for Z = 1 plasma, the ratio of RE density computed by CQL3D to that estimated from the model is within 0.6–6.0 in studied cases. For the Z = 18 case, this factor is systematically a much smaller number, typically 0.02–0.6. We suggest a simple correction to the model that narrows down the range of this ratio to 0.3–3.8 in all of the cases, including Z = 1 and Z = 18 plasmas.

Physics↗
DOE OSTI2021

High resolution density pedestal measurements during edge localized modes by short-pulse reflectometry in the TCV tokamak

Here, this publication presents high spatio-temporal resolution (mm µs -1 ) density profile measurements of the pedestal top during type I, III, and small edge localized mode (ELM) H-mode plasmas in the tokamak à configuration variable (TCV). These measurements were performed using a novel short-pulse reflectometer. Average inter-ELM density profiles are obtained via conditional averaging using the D α trace as ELM indicator. Changes to the pedestal density profile gradients prior to type-III ELMs reveal unique pedestal dynamics leading to the ELM crash which can provide important experimental data for validation of non-linear MHD ELM simulations. The small-ELM scenario is found to feature a ~25–35 kHz quasi-coherent density fluctuation near the separatrix ρ ψ 0.993–1.05 not observed during a similar type-I ELM discharge. This oscillation is also found in low-field-side magnetic pick-up probes displaying a ballooning character and n = +1 toroidal mode number. This oscillation could help explain the markedly different pedestal dynamics observed in the small-ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the relationship between the multi-region relaxed variational principle and resistive inner-layer theory

We show that the variational energy principle of the multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by 'ideal interfaces', and in each volume the magnetic field relaxes to a Taylor state, where the pressure gradient $\nabla p = 0$. The MRxMHD model is implemented in the Stepped-Pressure Equilibrium Code (SPEC) so that the equilibrium solution in each region is computed while preserving the force balance across the interfaces. As SPEC computes the Hessian matrix (a discretized stability matrix), the stability of an MRxMHD equilibrium can also be computed with SPEC. In this article, using SPEC, we investigate the effect of local pressure gradients and the $\nabla p = 0$ in the vicinity of the resonant surface of a tearing mode. For low-beta plasma, we have been able to illustrate a relationship between the resistive singular-layer theory (Coppi et al 1966 Nucl. Fusion6 101; Glasser et al 1975 Phys. Fluids18 875–88) and the MRxMHD model. Within the singular layer, the volume-averaged magnetic helicity and the flux-averaged toroidal flux are shown to be the invariants for the linear tearing modes in SPEC simulations. Our technique to compute MRxMHD stability is first tested numerically in a cylindrical tokamak and its application in toroidal geometry is demonstrated. We demonstrate an agreement between the stability boundary obtained with SPEC simulation and the resistive inner-layer theories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Plasma edge and scrape-off layer turbulence in gyrokinetic simulations of negative triangularity plasmas

Gyrokinetic simulations in the long-wavelength or drift-kinetic limit are carried out of DIII-D inner-wall-limited (IWL) plasmas to investigate the effect of triangularity on edge and scrape-off layer (SOL) turbulence. The effect of neutral interactions and triangularity on plasma blobs is explored due to the impact blobs can have in setting the SOL width or introducing impurities through interactions with plasma-facing components. Seeded blob simulations with neutrals in shaped SOL scenarios demonstrate that increasing elongation, triangularity, or Shafranov shift decreases radial blob velocities, but neutral interactions have a minor effect. Fully turbulent simulations of DIII-D IWL plasmas include both open- and closed-field-line regions. The negative triangularity (NT) simulation has lower average core T e , lower normalized T e fluctuations, and lower fluxes, but a greater number of coherent structures (blobs) identified with increased size and velocity, on average. Density and electron temperature profiles are within a factor of 2 of experimental values. In conclusion, the increased trapped electron particle fraction in NT simulations is consistent with previous studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Real-time reconstruction and control of pedestal-top electron density using RMP and gas puff at KSTAR

We report the experimental results of controlling the pedestal-top electron density by applying resonant magnetic perturbation (RMP) with in-vessel control coils and the main gas puff in the 2024-2025 KSTAR experimental campaign. The density is reconstructed using a parameterized $ψ_N$ grid and five channels of line-averaged density measured by the two-colored interferometer (TCI). The reconstruction procedure is accelerated by deploying a multi-layer perceptron to run in approximately 120 µ s and is sufficiently fast for real-time control. A proportional-integral controller was adopted, with the controller gains estimated from the system identification procedure. The experimental results demonstrate that the developed controller can follow a dynamic target while exclusively using both actuators. The absolute percentage errors between the electron density at $ψ_N$ = 0.89 and the target were approximately 1.5% median and a 2.5% average, respectively. The developed controller can even lower the density by using the pump-out mechanism under RMP, and it can follow a more dynamic range of density targets than a single actuator controller. The developed controller will enable experimental scenario exploration within a shot by dynamically setting the density target or maintaining a constant electron density within a discharge.

EFIT↗
DOE OSTI2026

Multimachine analysis of intermittent fluctuations in the scrape-off layer of magnetically confined fusion plasmas

Far scrape-off layer fluctuations, also known as blobs/filaments, are compared for the first time across six tokamaks: Alcator C-Mod, DIII-D, TCV, KSTAR, MAST, and MAST-U, under the context of stochastic modeling in Ohmic and low confinement mode scenarios. In this model, which well describes these fluctuations, measurements from Langmuir probes are modeled as a superposition of uncorrelated, two-sided exponential temporal pulses. The key parameters of this model are the mean waiting times between consecutive pulses, the pulse duration times and the mean of the pulse amplitudes. The study examines the influence of line-averaged density, plasma current, magnetic field and aspect ratio. The mean amplitudes and pulse waiting times increase with the line-averaged density across devices, while the pulse duration times remain constant. Whereas the variation in waiting times is gradual across the density scan, the mean amplitudes increase abruptly as the density profile broadens and flattens for Alcator C-Mod. The blob waiting times decrease with increasing plasma current and the toroidal magnetic field, while they increase with aspect ratio for conventional tokamaks. The duration times follow the same trends as the waiting times in these parameters but vary more slowly, resulting in the trends in the waiting times dominating the intermittency level of the fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Erosion and cathodic arc plasma of Nb–Al cathodes: composite versus intermetallic

Many properties of cathodic arcs from single-element cathodes show a correlation to the cohesive energy of the cathode material. For example, the burning voltage, the erosion rate, or, to a lesser extent, plasma properties like electron temperatures or average ion energy and charge states. For multi-element cathodes, various phases with different cohesive energies can initially be present in the cathode, or form due to arc exposure, complicating the evaluation of such correlations. To test the influence of morphology and phase composition of multi-element cathodes on cathodic arc properties, a Nb–Al cathode model system was used that includes: pure Nb and Al cathodes; intermetallic Nb 3 Al, Nb 2 Al and NbAl 3 cathodes; and three composite Nb–Al cathodes with atomic ratios corresponding to the stoichiometric ratios of the intermetallic phases. Pulsed cathodic arc plasmas from these cathodes were examined using a mass-per-charge and energy-per-charge analyzer, showing that charge-state-resolved ion energy distributions of plasmas from the intermetallic and corresponding composite cathodes are nearly identical. An examination of converted layers of eroded cathodes using x-ray diffraction and scanning electron microscopy indicates the formation of a surface layer with similar phase composition for intermetallic and their corresponding composite cathode types. The average arc voltages do not follow the trend of cohesive energies of Nb, Al and intermetallic Nb–Al phases, which have been calculated using density functional theory. Possible reasons for this effect are discussed based on the current knowledge of multi-element arc cathodes and their arc plasma available in literature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Time-resolved electron energy distribution functions at the substrate during a HiPIMS discharge with cathode voltage reversal

Abstract The time-dependent plasma properties of a high-power impulse magnetron sputtering plasma are investigated which include a positive ‘kick’ pulse on the sputtering target 2 μ s after the main negative pulse, this reversing the voltage on the cathode. At a substrate 15 cm distant from the magnetron, the time-dependent electron energy distribution function (EEDF), plasma potential, potential commute time and plasma diffusion properties are measured using a single Langmuir probe. Results show that the positive pulse on the target expels plasma and raises the plasma potential across the chamber on the order of 1 to 2 μ s, which is the time scale of the electron diffusion. The EEDF at the substrate fits a Druyvesteyn distribution during the main negative pulse rising slightly in average energy over time. The distribution is still Druyvesteyn and at the very start of the positive pulse, but then loses the higher energy electrons and drops in average electron energy as the positive pulse progresses. A Boltzmann equation solver, BOLSIG+, was used to predict the EEDF at the substrate during the positive pulse and it agrees best with the measurements assuming a value of 0.2 Td for the E / N (electric field/gas number density).

Physics↗
DOE OSTI2020

Curvature perturbations from stochastic particle production during inflation

We calculate the curvature power spectrum sourced by spectator fields that are excited repeatedly and non-adiabatically during inflation. In the absence of detailed information of the nature of spectator field interactions, we consider an ensemble of models with intervals between the repeated interactions and interaction strengths drawn from simple probabilistic distributions. We show that the curvature power spectra of each member of the ensemble shows rich structure with many features, and there is a large variability between different realizations of the same ensemble. Such features can be probed by the cosmic microwave background (CMB) and large scale structure observations. They can also have implications for primordial black hole formation and CMB spectral distortions. The geometric random walk behavior of the spectator field allows us to calculate the ensembleaveraged power spectrum of curvature perturbations semi-analytically. For sufficiently large stochastic sourcing, the ensemble averaged power spectrum shows a scale dependence arising from the time spent by modes outside the horizon during the period of particle production, in spite of there being no preferred scale in the underlying model. We find that the magnitude of the ensemble-averaged power spectrum overestimates the typical power spectra in the ensemble because the ensemble distribution of the power spectra is highly non-Gaussian with fat tails.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Measurement of the small-scale 3D Lyman- α forest power spectrum

Small-scale correlations measured in the Lyman-α (Lyα) forest encode information about the intergalactic medium and the primordial matter power spectrum. In this article, we present and implement a simple method to measure the 3-dimensional power spectrum, P 3D , of the Lyα forest at wavenumbers k corresponding to small, ~ Mpc scales. In order to estimate P 3D from sparsely and unevenly distributed data samples, we rely on averaging 1-dimensional Fourier Transforms, as previously carried out to estimate the 1-dimensional power spectrum of the Lyα forest, P 1D . Further, this methodology exhibits a very low computational cost. We confirm the validity of this approach through its application to Nyx cosmological hydrodynamical simulations. Subsequently, we apply our method to the eBOSS DR16 Lyα forest sample, providing as a proof of principle, a first P 3D measurement averaged over two redshift bins z = 2.2 and z = 2.4. This work highlights the potential for forthcoming P 3D measurements, from upcoming large spectroscopic surveys, to untangle degeneracies in the cosmological interpretation of P 1D .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Beam modulation and bump-on-tail effects on Alfvén eigenmode stability in DIII-D

Beam modulation effects on Alfvén eigenmode stability have been investigated in a recent DIII-D experiment and show that variations in neutral beam modulation period can have an impact on the beam driven Alfvén eigenmode spectrum and resultant fast ion transport despite similar time-averaged input power. The experiment was carried out during the current ramp phase of L-mode discharges heated with sub-Alfvénic 50–80 kV deuterium neutral beams that drive a variety of Alfvén eigenmodes unstable. The modulation period of two interleaved beams with different tangency radii was varied from shot to shot in order to modify the relative time dependent mix of the beam pitch angle distribution as well as the persistence of a bump-on-tail feature near the injection energy (a feature confirmed by imaging neutral particle analyzer measurements). As the beam modulation period is varied from 7 ms to 30 ms on/off (typical full energy slowing down time of τ slow ≈ 50 ms at mid-radius), toroidicity-induced Alfvén eigenmodes (TAEs) located in the outer periphery of the plasma become intermittent and coincident with the more tangential beam. Core mode activity changes from reversed shear Alfvén eigenmodes (RSAEs) to a mix of RSAE and beta-induced Alfvén eigenmodes. Discharges with 30 ms on/off period do not have a persistent bump-on-tail feature, have the lowest average mode amplitude and least fast ion transport. Furthermore, detailed analysis of an individual TAE using TRANSP kick modeling (Monte Carlo evolution of the distribution function with probabilistic 'kicks' by the AEs) and the resistive MHD code with kinetic fast ions, MEGA, find no strong role of energy gradient drive due to bump-on-tail features. Instead, the observed TAE modulation with interleaved beams is likely a pitch angle dependent result combined with slowing down of the tangential beam between pulses. For the conditions investigated, bump-on-tail contributions to TAE drive were found to be 5% or less of the total drive at any given time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

2D measurements of parallel counter-streaming flows in the W7-X scrape-off layer for attached and detached plasmas

Abstract Investigations of particle parallel flow velocities have been carried out for the scrape-off layer (SOL) of the Wendelstein 7-X (W7-X) stellarator, in order to gain insights on the SOL transport properties during attached and detached plasma scenarios. The experimental evidence is based on the coherence imaging spectroscopy (CIS) diagnostic, able to measure 2D impurity emission intensity and flow velocity. The impurity monitored by CIS is C 2+ , characterized by a line-emission intensity observed to be linearly proportional to the total plasma radiated power in both attached and detached plasmas. The related C 2+ velocity shows a strong dependence on the line-averaged electron density while remaining insensitive to the input power. During attached plasmas, the velocity increases with increasing line-averaged density. The tendency reverses in the transition to and during detachment, in which the velocity decreases by at least a factor of 2. The sharp drop in velocity, together with a rise in line-emission intensity, is reliably correlated to the detachment transition and can therefore be used as one of its signatures. The impurity flow velocity appears to be well coupled with the main ions’ one, thus implying the dominant role of impurity-main ion friction in the parallelimpurity transport dynamics. In view of this SOL impurity transport regime, the CIS measurement results are here interpreted with the help of EMC3-Eirene simulations, and their major trends are already explainable with a simple 1D fluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Disruption thermal load mitigation with shattered pellet injection on the Joint European Torus (JET)

Disruption mitigation remains a critical, unresolved challenge for ITER. To aid in addressing this challenge, a shattered pellet injection (SPI) system was installed on JET and experiments conducted at a range of thermal energy fractions and stored energies in excess of 7 MJ. The primary goals of these experiments were to investigate the efficacy of the SPI on JET and the ability of the plasma to assimilate multiple pellets. Single pellet injections produced a saturation in total radiated energy (Wrad) with increasing injected neon content, suggesting total radiation of stored thermal energy. Further increases in injected neon quantities resulted in reduced cooling times and current quench (CQ) durations, indicating higher impurity assimilation. No significant variation in CQ duration or W rad was observed when varying the deuterium content at fixed neon quantities. Additionally, higher assimilation, inferred by shorter CQ durations, was measured when a mechanical punch was used to launch the pellets and this was attributed to a lower pellet velocity leading to higher solid content in the pellet plume and larger fragments penetrating deeper into the plasma. Radiation asymmetries averaged over the cooling time were inferred from Emis3D and ranged from 1.6 to 1.9. Asymmetries averaged over the entire disruption sequence were found to increase at higher thermal energy fractions. The radiated energy fractions decreased with increasing thermal energy fractions but this trend was eliminated when toroidal asymmetries were accounted for with Emis3D. Pure deuterium pellets were able to produce cooling times of up to 75 ms with a gradual loss in thermal stored energy of up to 80%. Experiments with multiple pellet injection indicated W rad can be increased through pellet superposition and density can be increased with an additional D2 injection without a reduction in W rad . KPRAD modelling accurately reproduced the cooling times and the CQ duration at high thermal energies. Assimilation estimates from KPRAD indicated CQ rates scale strongly whilst W rad scales weakly and saturates with assimilated neon content. Comparable W rad can be achieved with lower assimilated neon quantities as longer cooling times are attained. Thus reduced neon content can be preferential in a thermal load mitigation scheme as it may reduce radiation asymmetries and prevent flash melting.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗