Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Electron Beam Modeling and Simulation”

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

Turbulence and wave particle interactions in solar-terrestrial plasmas

Activities in the following study areas are reported: (1) particle and wave processes in solar flares; (2) solar convection zone turbulence; and (3) solar radiation emission. To investigate the amplification of cyclotron maser radiation in solar flares, a radio frequency. (RF) heating model was developed for the corona surrounding the energy release site. Then nonlinear simulations of compressible convection display prominent penetration by plumes into regions of stable stratification at the base of the solar convection zone, leading to the excitation of internal gravity waves there. Lastly, linear saturation of electron-beam-driven Langmuir waves by ambient density fluctuations, nonlinear saturation by strong turbulence processes, and radiation emission mechanisms are examined. An additional section discusses solar magnetic fields and hydromagnetic waves in inhomogeneous media, and the effect of magnetic fields on stellar oscillation.

Dulk, G. A.↗

Transport of kJ-laser-driven relativistic electron beams in cold and shock-heated vitreous carbon and diamond

We report experimental results on relativistic electron beam (REB) transport in a set of cold and shock-heated carbon samples using the high-intensity kilojoule-class OMEGA EP laser. The REB energy distribution and transport were diagnosed using an electron spectrometer and x-ray fluorescence measurements from aCu tracer buried at the rear side of the samples. The measured rear REB density shows brighter and narrower signals when the targets were shock-heated. Hybrid PIC simulations using advanced resistivity models in the target warm-dense-matter (WDM) conditions confirm this observation. We show that the resistivity response of the media, which governs the self generated resistive fields, is of paramount importance to understand and correctly predict the REB transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Secondary Beams at Jefferson Laboratory

This thesis explores the possibility of producing secondary beams at the Thomas Jefferson National Accelerator Facility (TJNAF). Their possible applications in future experiments are studied using advanced Monte Carlo simulations with state-of-the-art software. In Jefferson Lab’s Hall-A, the interaction of the high-current O(100) µA, medium-energy O(11) GeV electron beam produced by the Continuous Electron Beams Accelerator Facility (CEBAF) with the Beam-Dump (BD) generates a substantial shower of standard model particles. This interaction creates a high-intensity secondary muon and neutrino beam and hypothetical Light Dark Matter (LDM) particles. Additionally, with the proposed Compact Photon Source (CPS), Hall-D could offer the possibility of producing an intense muon beam alongside a photon beam.

Fulci, Antonino↗

Simulation of DIII-D disruption with argon pellet injection and runaway electron beam

Abstract The next generation of large tokamaks, including ITER, will be equipped with a disruption mitigation system (DMS) that can be activated if a disruption is deemed to be imminent. Introducing impurities by pellet (large or shattered) or massive gas injection has been shown to be an effective mitigation mechanism on many tokamaks. The goal of the mitigation is to lessen the thermal and electromagnetic loads from the disruption without generating enough high-energy (runaway) electrons to damage the device. Variations of this mitigation process with impurity injection are presently being tested on many experiments. We have modeled one such impurity injection experiment on DIII-D using the M3D-C1 nonlinear 3D extended MHD code (Jardin et al 2012 Comput. Sci. Discovery 6 014002), The model includes an argon large pellet injection and ablation model, impurity ionization, recombination, and radiation, and runaway electron formation and subsequent evolution, including both Dreicer and avalanche sources. We obtain reasonable agreement with the experimental results for the timescale of the thermal and current quench and for the magnitude of the runaway electron plateau formed during the mitigation. This is the first 3D full MHD simulation with pellets and REs to simulate the disruption process and it also provides a partial validation of the M3D-C1 DMS model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energy transport and dynamics

We report findings concerning energy transport and dynamics in flares during the impulsive and gradual phases based on new ground-based and space observations (notably from Yohkoh). A preheating sometimes occurs during the impulsive phase. Ca XIX line shifts are confirmed to be good tracers of bulk plasma motions, although strong blue shifts are not as frequent as previously claimed. They often appear correlated with hard X-rays but, for some events, the concept that electron beams provide the whole energy input to the thermal component seems not to apply. Theory now yields: new diagnostics of low-energy proton and electric beams; accurate hydrodynamical modeling of pulse beam heating of the atmosphere; possible diagnostics of microflares (based on X-ray line ratio or on loop variability); and simulated images of chromospheric evaporation fronts. For the gradual phase, the continual reorganization of magnetic field lines over active regions determines where and when magnetic reconnection, the mechanism favored for energy release, will occur. Spatial and temporal fragmentation of the energy release, observed at different wavelengths, is considered to be a factor as well in energy transport and plasma dynamics.

Schmieder, Brigitte↗

Measurements of Beam Spin Asymmetries in p+p0 and p´p0 Dihadron Production at CLAS12

Semi-Inclusive Deep Inelastic Scattering (SIDIS) is a powerful experimental tool for studying the internal structure and dynamics of the proton, revealing how quarks and gluons are distributed and interact within it. SIDIS describes a process where an elec tron scatters off one of the constituent quarks within the proton, causing it to undergo hadronization, creating multiple hadrons in the final state. Through factorization, the full process can be split into probabilistic components: one which describes the internal structure of the proton using Parton Distribution Functions (PDFs), and another which describes the hadronization process using Fragmentation Functions (FFs). These functions are non-perturbative quantities of Quantum Chromodynamics (QCD), meaning they cannot be calculated directly from first principles and must instead be extracted from experimental measurements. Acommon approach for accessing PDFs and FFs using SIDIS is to measure asymmetries. In this context, asymmetries correspond to subtle differences in the angular distribution of outgoing particles that arise when the spin orientation of the incoming beam or target is reversed. Because many of these effects only appear when spin is involved, they isolate specific, nuanced properties of the proton’s spin-structure that are otherwise hidden in spin averaged measurements. In practice, they show up as specific azimuthal modulations (e.g., sin ¿R, sin(¿h ´ ¿R)), whose amplitudes isolate convolutions of PDFs and FFs at leading and subleading twist. Non-zero asymmetries of these angular distributions can be traced back to unique combinations of PDFs and FFs, offering a way to probe them directly. In this work, we measure SIDIS by analyzing high energy electron-proton scattering events using the CLAS12 detector at Jefferson Lab. This study focuses on subset of SIDIS referred to as dihadron SIDIS, where pairs of hadrons — here p+p0 and p´p0 — are observed. We analyzed these dihadrons using detector data collected during Fall 2018 and Spring 2019, where longitudinally polarized electrons from the CEBAF accelerator were incident on a liquid hydrogen target. A photon classifier using a Gradient Boosted Trees (GBTs) architecture was trained using Monte Carlo simulations to reduce the amount of iv false combinatorial background p0’s. When deployed on experimental data, the model in creases our dihadron statistics by up to five-fold compared to previous CLAS12 p0 analyses. This work reports the first measurements of beam spin asymmetries for p+p0 and p´p0 dihadron production in SIDIS. The measured asymmetries offer new insights to the spin-dependent structure and dynamics within the proton, as well as the spin-dependent properties of quark fragmentation. Non-zero twist-3 sin¿R amplitudes are observed, pro viding sensitivity to the subleading twist PDF e(x). The PDF e(x) encodes quark-gluon correlations within the proton — a property that is otherwise inaccessible at leading twist. Additionally, this work measured significant twist-2 modulations carried by sin(¿h ´ ¿R) and sin(2¿h ´2¿R), providing experimental access to the helicity dihadron fragmentation function (DiFF) GK 1 . Because there is no equivalent quark helicity-dependent FF in single pion SIDIS, the DiFF GK 1 offers a unique lens into novel spin-dependent fragmentation. For instance, the twist-2 modulations observed in this study are enhanced by vector mesons created during fragmentation — a behavior predicted by phenomenological models. This study broadens our understanding of dihadron fragmentation, revealing new details about the flavor and charge dependence of hadronization.

Matousek, Gregory [Duke Univ., Durham, NC (United ↗

Overview of interpretive modelling of fusion performance in JET DTE2 discharges with TRANSP

In the paper we present an overview of interpretive modelling of a database of JET-ILW 2021 D-T discharges using the TRANSP code. The main aim is to assess our capability of computationally reproducing the fusion performance of various D-T plasma scenarios using different external heating and D-T mixtures, and to understand the performance driving mechanisms. We find that interpretive simulations confirm a general power-law relationship between increasing external heating power and fusion output, which is supported by absolutely calibrated neutron yield measurements. A comparison of measured and computed D-T neutron rates shows that the calculations' discrepancy depends on the absolute neutron yield. The calculations are found to agree well with measurements for higher performing discharges with external heating power above ~20 MW, while low-neutron shots display an average discrepancy of around +40% compared to measured neutron yields. A similar trend is found for the ratio between thermal and beam-target fusion, where larger discrepancies are seen in shots with dominant beam-driven performance. We compare the observations to studies of JET-ILW D discharges, to find that on average the fusion performance is well modelled over a range of heating power, although an increased unsystematic deviation for lower-performing shots is observed. The ratio between thermal and beam-induced D-T fusion is found to be increasing weakly with growing external heating power, with a maximum value of ≳1 achieved in a baseline scenario experiment. An evaluation of the fusion power computational uncertainty shows a strong dependence on the plasma scenario type and fusion drive characteristics, varying between ±25% and 35%. D-T fusion alpha simulations show that the ratio between volume-integrated electron and ion heating from alphas is ≲10 for the majority of analysed discharges. Alphas are computed to contribute between ~15% and 40% to the total electron heating in the core of highest performing D-T discharges. An alternative workflow to TRANSP was employed to model JET D-T plasmas with the highest fusion yield and dominant non-thermal fusion component because of the use of fundamental radio-frequency heating of a large minority in the scenario, which is calculated to have provided ~10% to the total fusion power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sodium tracer measurements of an expanded dense aluminum plasma from e-beam isochoric heating

Spatially and temporally resolved visible absorption spectroscopy is performed on sodium D-lines present as surface contaminants on an expanded dense aluminum plasma plume. An 80-ns FWHM, intense, relativistic electron beam deposits 5.4 J into a 100- μ m-thick Al foil, which isochorically heats and subsequently hydrodynamically expands the material through the warm dense matter state and into a classical-like plasma state, with a coupling parameter of approximately 0.2 and a degeneracy parameter of approximately 270. The Na contamination, carried along with the expanding plume, shows saturated absorption features in the dense Al continuum for λ > 450 nm. X-ray photoelectron spectroscopy and laser-induced breakdown spectroscopy confirm Na is a surface contaminant with an atomic concentration of ~ 0.1 % when interrogating identical foil samples. Additioanlly, a spectroscopic-quality radiation transport model is used to post-process 2D hydrodynamic simulations to interpret the plasma conditions based on the measured Na 3p-3s doublet line profiles. A sodium number density of 3 × 10 15 cm -3 best matches the experimental spectra, which originate from a dense surface plasma with n e = 3.0 ± 0.8 × 10 18 cm -3 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of the imaging neutral particle analyzer in nearly MHD quiescent plasmas using injected beam ions on DIII-D

The imaging neutral particle analyzer (INPA) is a scintillator-based diagnostic that provides energy and radially resolved measurements of confined fast ions. To verify operation, understanding, and modeling of the diagnostic performance, distributions of deuterium beam ions in low density, nearly MHD-quiescent plasmas are measured on the DIII-D tokamak using the INPA. Here, the distribution is altered by changing the pitch-angle scattering rate with electron cyclotron heating (ECH). Experimental images are mostly reproduced with simulation of core fast ions charge exchanging with injected beam and edge cold neutrals in classical cases. Image comparisons between the experiment and simulation show agreement with less than 25% discrepancies near the beam injection energy of 50 keV. Other strong passive signals, observed in a localized region on the image, are suggested to come from charge exchange between cold neutrals and fast ions around the plasma boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Suppression of Hydrogen Emission in an X-Class White-Light Solar Flare

We present unique NUV observations of a well-observed X-class flare from NOAA 12087 obtained at the Ondrejov Observatory. The flare shows a strong white-light continuum but no detectable emission in the higher Balmer and Lyman lines. Reuven Ramaty High-Energy Solar Spectroscopic Imager and Fermi observations indicate an extremely hard X-ray spectrum and gamma-ray emission. We use the RADYN radiative hydrodynamic code to perform two types of simulations: one where an energy of 3 x 10(exp 11) erg/sq cm/s is deposited by an electron beam with a spectral index of approx. = 3, and a second where the same energy is applied directly to the photosphere. The combination of observations and simulations allows us to conclude that the white-light emission and the suppression or complete lack of hydrogen emission lines is best explained by a model where the dominant energy deposition layer is located in the lower layers of the solar atmosphere, rather than the chromosphere.

Prochazka, Ondrej↗

Random Forest Prediction of Crystal Structure from Electron Diffraction Patterns

Transmission electron microscopy (TEM) diffraction patterns are regularly used to determine the structure of crystalline materials. Electron diffraction is the most common method to solve for unknown or partially known crystal structures, as it provides direct and interpretable feedback on the orientation of crystal grains under the beam [1]. However, it remains a challenge to determine the crystal structure of a new material or even a new phase of an existing material. Analysis of such materials commonly requires manual exploration and comparison with simulated diffraction patterns. This is often a time consuming process with no obvious start point when many similar structures are possible, and this method cannot be used to determine crystal structure or orientation from structures not included in the diffraction libraries. Therefore, we have developed a machine learning model to determine the crystal structure of a material from its electron diffraction pattern.

36 MATERIALS SCIENCE↗

Fatigue life prediction of powder bed fused–laser beam AlSi10Mg: Incorporating critical defects via crystal plasticity modelling

The current study provides a microstructurally-based computational framework to predict the fatigue life of additive manufactured (AM), i.e., powder bed fused–laser beam (PBF-LB), AlSi10Mg specimens using the crystal plasticity finite element method (CPFEM). The fractography analysis, electron backscatter diffraction (EBSD), uniaxial and cyclic responses, and fatigue life of specimens were used to inform the computational framework. CPFE simulation was used to compute fatigue indicator parameters (FIPs) as fatigue driving forces. A new fatigue criterion is introduced based on FIPs, which was calibrated using experimental fatigue data. The proposed fatigue measure was evaluated versus the specimens with critical defects of various sizes and locations subjected to different stress amplitudes. The results show that the developed framework can capture the fatigue life of samples with different critical defect locations and sizes along with different stress amplitudes for both high-cycle fatigue (HCF) and very high-cycle fatigue (VHCF) regimes.

Additive manufacturing↗

DSMC and continuum analyses of low-density nozzle flow

Two different approaches, the direct-simulation Monte Carlo (DSMC) method based on molecular gas dynamics and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gas dynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental dam, in which density and rotational temperature variations along the centerline and at various locations inside a low density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chan-Hong↗

Low-Density Nozzle Flow by the Direct Simulation Monte Carlo and Continuum Methods

Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chang-Hong↗

Langmuir turbulence in the auroral ionosphere 2: Nonlinear theory and simulations

A theoretical interpretation of sounding rocket measurements of intense Langmuir wave fields (less than or equal to 500 mV/m) driven by a stream of 20 eV to 4 keV electrons in the lower auroral zone is developed. This interpretation is based on the ability of the 10 microseconds sampling rate of the wave detector to temporally resolve the structure of the Langmuir wave field envelope. A modified form of the Zakharov equations is used to numerically study beam-driven Langmuir turbulence in the presence of a moderate magnetic field (OMEGA (sub e) approximately equals Omega (sub pe). Strong Landau damping on observed nonthermal scattered electrons, which is treated in a companion paper (Newman et al., this issue), plays an important role by inhibiting backscatter cascade and the development of strong turbulence. A parameterized model of the linear electron stream-driven wave instability is introduced, which incorporates limited quasilinear plateau formation. A reasonable set of parameters is found that yields semiquantitative agreement between observed properties of the Langmuir fields and the results of Zakharov equation simulations, including the amplitude and characteristic frequency of the electric field envelope modulations.

Newman, D. L.↗

Cross-field Streaming Instabilities in the Lower Hybrid Frequency Near the Front Pole Regions of a Magnetically Shielded Hall Thruster

The mechanism that led to the erosion of the inner front pole cover in wear tests of the Hall-effect Rocket with Magnetic Shielding (HERMeS) at 300 V and 12.5 kW has been elusive. Several classical mechanisms have been considered yet none has explained conclusively the source of the discrepancy between wear measurements and simulations. In this article we consider a new mechanism. The electrostatic dispersion relation for an unbounded, homogeneous plasma in a fixed magnetic field is solved numerically in the near-plume region of this thruster. The plasma conditions have been obtained from Hall2De (r-z) multifluid-PIC simulations. We find that the modified two-stream instability can be excited in most regions of the domain, with the fastest growth approaching the fluid-limit value of ωLH/2 near the magnetic poles where ωpe/ωce~1. In this region the relative drift between magnetized electrons and beam ions perpendicular to the magnetic field can exceed 3uTi, while that between counter-streaming beam and cathode ions can exceed 4uTi. The latter suggests that the ion-ion cross-field lower hybrid instability may also be active here. Unlike the classical ion-acoustic instability, lower hybrid instabilities can heat ions to temperatures that are comparable to the electron temperature. Such anomalous heating would enhance the average energy of these ions and, in turn, the erosion of the nearby surfaces they sputter. An idealized model that scales the anomalous ion heating with ωLH has been incorporated in Hall2De and the simulations now yield significantly improved agreement with the measured wear rates and IVDFs, as reported in our companion paper.

Ortega, Alejandro Lopez↗

Initial Investigation of Coupled Electron-Photon Transport with SCEPTRE and RAPTURE (Rev.1.1)

The SCEPTRE1 code models coupled photon and electron transport in material media using discrete ordinates methods. The code can use different solvers in 1D and 3D in parallel, promising users the ability to run problems that are significant both in physics and computational scale. One capability that is perhaps unique at this time in the NNSA complex is SCEPTRE’s coupled photon– electron transport physics in a discrete ordinates approach. For example, incident x–rays can generate a significant electron population, whose transport can be simulated within the same code. In this context, fully–coupled means that photons can produce electrons and electrons can produce photons; uncoupled or no–coupling means photons do not produce electrons. Partial–coupling, an intermediate option in which photons can produce electrons but electrons do not produce photons, is not considered in the present work. The MCNP code can treat this physics with a Monte Carlo approach. MCNP should be viewed as complementary to SCEPTRE in both computational method (Monte Carlo versus discrete ordinates) as well as data libraries. SCEPTRE and MCNP should provide important checks on each other, especially from complementary modeling frameworks of a S n (discrete ordinates) mesh based method versus Monte Carlo particle based method. In order to learn more about SCEPTRE, and how it might inform engineering work at the Laboratory, we used it to carry out a set of initial studies. This was a series of 1D simulations of photon beams hitting Aluminum slabs. Our studies mainly used SCEPTRE directly, and made little use of the closely associated RAPTURE code (more on this below). The goals of these studies were first, to familiarize ourselves with operating the code, and second, to begin to assess the importance of coupled electron–photon transport. To support these studies, we began to develop a set of scripts to initialize inputs, run tools, process outputs, and plot results. An initial version of these scripts has been committed to the Lightning workflow repository. This memo relates the results of those studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Experimental and computational study of phase space dynamics in strongly coupled plasmas with steep density gradients

Understanding how plasmas thermalize when density gradients are steep remains a fundamental challenge in plasma physics, with direct implications for fusion experiments and astrophysical phenomena. Standard hydrodynamic models break down in these regimes, and kinetic theories make predictions that have never been directly tested. Here, we present the first detailed phase-space measurements of a strongly coupled plasma as it evolves from sharp density gradients to thermal equilibrium. Using laser-induced fluorescence imaging of an ultracold calcium plasma, we track the complete ion distribution function f(x,v,t). We discover that commonly used kinetic models (Bhatnagar–Gross–Krook and Lenard–Bernstein) overpredict thermalization rates, even while correctly capturing the initial counterstreaming plasma formation. Our measurements reveal that the initial ion acceleration response scales linearly with electron temperature, and that the simulations underpredict the initial ion response. In our geometry we demonstrate the formation of well-controlled counterpropagating plasma beams. This experimental platform enables precision tests of kinetic theories and opens new possibilities for studying plasma stopping power and flow-induced instabilities in strongly coupled systems.

Bergeson, Scott (ORCID:0000000231249226)↗