Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Beam Simulations”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Impact of increased smoothing by spectral dispersion bandwidth on stimulated Brillouin scattering in laser driven Hohlraums

Experiments were conducted at the National Ignition Facility (NIF) to investigate the impact of increased smoothing by spectral dispersion (SSD) bandwidth on the production of stimulated Brillouin scattering (SBS) within an indirect-drive inertial confinement fusion (ICF) Hohlraum. This was done in a subscale gold Hohlraum driven by 192 laser beams depositing 1.1 MJ of energy. The laser bandwidth was increased from 45 to 118 GHz (before frequency tripling) on the 30° and 50° cones, where backscatter could be measured using the Full Aperture Backscatter Station (FABS). It was expected that this 2.6-fold increase in bandwidth would suppress SBS generated within the Hohlraum plasma and lower the backscattered SBS energy on the 50° cones by over a factor of four. Experimental results, however, show that this SSD change only reduced the 50° cone SBS during the main capsule drive by -18±31% and -4.5±7.8% over the entire pulse. This is small compared to expected shot-to-shot SBS reproducibility (∼30%), such that the result can be considered within normal performance fluctuations. New 3D parallel paraxial code (pF3D) simulations, accounting for beam refractive intensification reproduce this result, suggesting that closer to 300 GHz of bandwidth would have been required to mitigate SBS to the expected level. Delivering such a high bandwidth is not feasible at NIF when operating at high peak power and would potentially prevent NIF’s ability to use cross-beam energy transfer (CBET) for implosion symmetry tuning.

Physics↗
DOE OSTI2025

High efficiency RF sources developments

Calabazas Creek Research, Inc. (CCR) and its collaborators are developing high efficiency RF sources operating from a few hundred MHz to C-Band and power levels from tens to hundreds of kilowatts with the goal of providing MW-relevant sources. The efficiencies approach or exceed 80% with projected costs as low as $0.50/ Watt. Sources under development include magnetrons with phase and amplitude control, single and multi-beam klystrons, multi-beam power grid tubes, and multiple beam IOTs. A magnetron system achieved more than 80% efficiency with fast amplitude control using modulation of the phase locking signal. This would be a low cost, high efficiency RF source for superconducting accelerators. An L-Band, single beam klystron was built with simulated efficiency of 80%. The klystron has yet to be tested to confirm the simulation results. CCR is currently developing a multi-beam klystron to produce more than 200 kW CW at 80% efficiency. Also in development is a multiple beam triode to produce 200 kW CW from 300 MHz to approximately 1 GHz. Not only does the simulated efficiency exceed 75%, but it would be the lowest cost RF source in this frequency range. Finally, CCR recently concluded research for a multiple beam IOT at 700 MHz using third harmonic drive to boost efficiency toward 85%. Successful development and transition to production of these sources will significantly alter the cost/performance landscape for RF power generation.

43 PARTICLE ACCELERATORS↗
NASA NTRS1989

A current disruption mechanism in the neutral sheet for triggering substorm expansions

Two main areas were addressed in support of an effort to understand mechanism responsible for the broadband electrostatic noise (BEN) observed in the magnetotail. The first area concerns the generation of BEN in the boundary layer region of the magnetotail whereas the second area concerns the occassional presence of BEN in the neutral sheet region. For the generation of BEN in the boundary layer region, a hybrid simulation code was developed to perform reliable longtime, quiet, highly resolved simulations of field aligned electron and ion beam flow. The result of the simulation shows that broadband emissions cannot be generated by beam-plasma instability if realistic values of the ion beam parameters are used. The waves generated from beam-plasma instability are highly discrete and are of high frequencies. For the plasma sheet boundary layer condition, the wave frequencies are in the kHz range, which is incompatible with the observation that the peak power in BEN occur in the 10's of Hz range. It was found that the BEN characteristics are more consistent with lower hybrid drift instability. For the occasional presence of BEN in the neutral sheet region, a linear analysis of the kinetic cross-field streaming instability appropriate to the neutral sheet condition just prior to onset of substorm expansion was performed. By solving numerically the dispersion relation, it was found that the instability has a growth time comparable to the onset time scale of substorm onset. The excited waves have a mixed polarization in the lower hybrid frequency range. The imposed drift driving the instability corresponds to unmagnetized ions undergoing current sheet acceleration in the presence of a cross-tail electric field. The required electric field strength is in the 10 mV/m range which is well within the observed electric field values detected in the neutral sheet during substorms. This finding can potentially account for the disruption of cross-tail current and its diversion to the ionosphere to form the substorm current wedge. Furthermore, a number of features associated with substorm expansion onset can be understood based on this substorm onset scenario.

Lui, A. T. Y.↗
DOE OSTI2023

Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas

Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Implementation of compound refractive lenses for large field-of-view x-ray phase-contrast imaging during hypervelocity impact experiments

Synchrotron x-ray phase-contrast imaging (XPCI) offers time-resolved visualization of dynamic compression phenomena, but its intrinsically small field-of-view (FOV) limits the time that key features remain in frame. A novel approach to enlarge the FOV is achieved by positioning a two-dimensional parabolic compound refractive lens (CRL) upstream of the sample to deliberately defocus the white beam. Ray-tracing simulations and XPCI measurements show that this CRL configuration can expand the beam by ∼50% vertically and ∼15% horizontally based on the full width at half-maximum of the beam. Implementing the CRL, however, attenuates the photon flux and lowers signal-to-noise ratio (SNR). Task-based analysis using a calibration grid (30 μm dots) showed that both setups fail to consistently meet the Rose criterion (SNR ≥ 5) for features of this size in single-bunch imaging. Extrapolating the measured SNR Rose values suggests that the minimum consistently detectable feature lies closer to 30–40 μm for the standard XPCI setup and above 40 μm for CRL-XPCI. Despite this limitation, the CRL configuration nearly doubles the illuminated area, enabling simultaneous tracking of front and rear observations of boron carbide targets subjected to rod and sphere impacts at 1.0–2.6 km/s. Image tracking algorithms and photonic Doppler velocimetry were used to measure penetration and rear-surface velocity histories. Together, these measurements capture crack fronts, penetration, and material breakout, offering new benchmark data for validating high-strain-rate constitutive models of ceramic materials.

Ceramic materials↗
DOE OSTI2023

A simple, flexible technique for RF cavity wake-field calculations

It is typical in the accelerator field to model machine components, especially RF cavities, as parallel RLC resonators. In the interest of simulating and diagnosing beam instabilities, knowledge of the time-domain voltage waveform over an equivalent resonator by a bunch current often proves useful. This waveform may be found by convolving the bunch current with the RLC resonator impulse response. While analytical and quasi-analytical expressions are available in this regime for common distributions such as the Gaussian, analogous results for less standard distributions are difficult to obtain using direct methods, which opens the door for the development of a more generalized technique. In this paper, a formulation is created that allows for the simple computation of the time-domain voltage waveform of and RLC resonator. The formulation uses the Cauchy Residue Theorem to extract the convolution result from the Fourier Domain, and it only requires that the current distribution Fourier Transform be holomorphic and known at one specific evaluation point. This greatly simplifies the computation of the time domain voltage for a large amount of bunch distributions both common and uncommon.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Plasma-based particle sources

High-brightness particle beams generated by advanced accelerator concepts have the potential to become an essential part of future accelerator technology. In particular, high-gradient accelerators can generate and rapidly accelerate particle beams to relativistic energies. The rapid acceleration and strong confining fields can minimize irreversible detrimental effects to the beam brightness that occur at low beam energies, such as emittance growth or pulse elongation caused by space charge forces. Due to the high accelerating gradients, these novel accelerators are also significantly more compact than conventional technology. Advanced accelerators can be extremely variable and are capable of generating particle beams with vastly different properties using the same driver and setup with only modest changes to the interaction parameters. So far, efforts have mainly been focused on the generation of electron beams, but there are concepts to extend the sources to generate spin-polarized electron beams or positron beams. The beam parameters of these particle sources are largely determined by the injection and subsequent acceleration processes. Although, over the last decade there has been significant progress, the sources are still lacking a sufficiently high 6-dimensional (D) phase-space density that includes small transverse emittance, small energy spread and high charge, and operation at high repetition rate. This is required for future particle colliders with a sufficiently high luminosity or for more near-term applications, such as enabling the operation of free-electron lasers (FELs) in the X-ray regime. Major research and development efforts are required to address these limitations in order to realize these approaches for a front-end injector for a future collider or next-generation light sources. In particular, this includes methods to control and manipulate the phase-space and spin degrees-of-freedom of ultrashort plasma-based electron bunches with high accuracy, and methods that increase efficiency and repetition rate. These efforts also include the development of high-resolution diagnostics, such as full 6D phase-space measurements, beam polarimetry and high-fidelity simulation tools. A further increase in beam luminosity can be achieve through emittance damping. Emittance cooling via the emission of synchrotron radiation using current technology requires kilometer-scale damping rings. For future colliders, the damping rings might be replaced by a substantially more compact plasma-based approach. Here, plasma wigglers with significantly stronger magnetic fields are used instead of permanent-magnet based wigglers to achieve similar damping performance but over a two orders of magnitude reduced length.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Advancing Cross-Sectional Scanning Electron Microscopy of Perovskite Solar Cells

Organic–inorganic perovskites are an emerging class of photovoltaic materials. Despite achieving power conversion efficiencies surpassing 26%, the challenge of perovskite stability including degradation during exposure to operational conditions such as light, heat, humidity, water, oxygen, and electric fields is well known. Related, perovskite instability has limited high-resolution electron imaging and characterization techniques that can be used for understanding degradation mechanisms. Furthermore, we demonstrate perovskite device cross-section preparation using mechanical polishing in a water-free environment with cryogenic Ar ion milling. Scanning electron microscopy was then used in both backscattered electron and secondary electron imaging modes to obtain information about layer structure, grain aggregate structure, and compositional heterogeneity. Monte Carlo CASINO simulations inform optimum beam conditions and image acquisition parameters and the effects of accelerating voltage, dwell times, and frame averaging for practical image acquisition are reported.

14 SOLAR ENERGY↗
DOE OSTI2021

Stopping-power enhancement from discrete particle-wake correlations in high-energy-density plasmas

Three-dimensional (3D) simulations of electron beams propagating in high-energy-density plasmas using the quasistatic Particle-in-Cell (PIC) code QuickPIC demonstrate a significant increase in stopping power when beam electrons mutually interact via their wakes. Each beam electron excites a plasma wave wake of wavelength ~ 2 π c / ω p e , where c is the speed of light and ω p e is the background plasma frequency. Here we show that a discrete collection of electrons undergoes a beam-plasma-like instability caused by mutual particle-wake interactions that causes electrons to bunch in the beam, even for beam densities n b for which fluid theory breaks down. This bunching enhances the beam's stopping power, which we call “correlated stopping,” and the effect increases with the “correlation number” N b ≡ n b ( c / ω p e ) 3 . For example, a beam of monoenergetic 9.7 MeV electrons with N b = 1 / 8 , in a cold background plasma with n e = 10 26 cm – 3 (450 g cm – 3 DT), has a stopping power of 2.28 ± 0.04 times the single-electron value, which increases to 1220 ± 5 for N b = 64 . The beam also experiences transverse filamentation, which eventually limits the stopping enhancement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Distributed Stochastic Optimization of a Neural Representation Network for Time-Space Tomography Reconstruction

4D time-space reconstruction of dynamic events or deforming objects using X-ray computed tomography (CT) is an important inverse problem in non-destructive evaluation. Conventional back-projection based reconstruction methods assume that the object remains static for the duration of several tens or hundreds of X-ray projection measurement images (reconstruction of consecutive limited-angle CT scans). However, this is an unrealistic assumption for many in-situ experiments that causes spurious artifacts and inaccurate morphological reconstructions of the object. To solve this problem, we propose to perform a 4D time-space reconstruction using a distributed implicit neural representation (DINR) network that is trained using a novel distributed stochastic training algorithm. Our DINR network learns to reconstruct the object at its output by iterative optimization of its network parameters such that the measured projection images best match the output of the CT forward measurement model. Here, we use a forward measurement model that is a function of the DINR outputs at a sparsely sampled set of continuous valued 4D object coordinates. Unlike previous neural representation architectures that forward and back propagate through dense voxel grids that sample the object's entire time-space coordinates, we only propagate through the DINR at a small subset of object coordinates in each iteration resulting in an order-of-magnitude reduction in memory and compute for training. DINR leverages distributed computation across several compute nodes and GPUs to produce high-fidelity 4D time-space reconstructions. We use both simulated parallel-beam and experimental cone-beam X-ray CT datasets to demonstrate the superior performance of our approach.

36 MATERIALS SCIENCE↗
DOE OSTI2022

X-ray pulse generation with ultra-fast flipping of its orbital angular momentum

A method to temporally tailor the properties of X-ray radiation carrying Orbital Angular Momentum (OAM) is presented. In simulations, an electron beam is prepared with a temporally modulated micro-bunching structure which, when radiating at the second harmonic in a helical undulator, generates OAM light with a corresponding temporally modulated intensity. This method is shown to generate attosecond pulse trains of OAM light without the need for any additional external optics, making the wavelength range tunable. In addition to the OAM pulse train, the method can be adapted to generate radiation where the handedness of the OAM mode may also be temporally modulated (flipped).

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

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↗
DOE OSTI2025

Current-Regulated Arc Modulator Optimization

n this study, the performance of a current-regulated arc modulator was investigated with a focus on its role in initiating and sustaining plasma discharge within the Magnetron Body of the LINAC system. The analysis centered on how switching components, circuit topology, and feedback loop architecture influence critical factors such as energy efficiency, discharge stability, and long-term plasma containment. Particular attention was given to variations in pulse termination behavior, as observed through oscilloscope traces, which revealed inconsistencies affecting the duty factor and cathode temperature. These fluctuations have downstream effects on the cesium-coated cathode surface, thereby impacting H⁻ ion production and beam reliability. Simulation-based testing in LTspice was used to evaluate noise suppression techniques and arc current regulation schemes, revealing how optimized snubber networks, improved pulse shaping, and feedback stability can mitigate modulator-induced noise. The results ide

Campos, Nathan↗
NASA NTRS1979

Gyrodampers for large space structures

The problem of controlling the vibrations of a large space structures by the use of actively augmented damping devices distributed throughout the structure is addressed. The gyrodamper which consists of a set of single gimbal control moment gyros which are actively controlled to extract the structural vibratory energy through the local rotational deformations of the structure, is described and analyzed. Various linear and nonlinear dynamic simulations of gyrodamped beams are shown, including results on self-induced vibrations due to sensor noise and rotor imbalance. The complete nonlinear dynamic equations are included. The problem of designing and sizing a system of gyrodampers for a given structure, or extrapolating results for one gyrodamped structure to another is solved in terms of scaling laws. Novel scaling laws for gyro systems are derived, based upon fundamental physical principles, and various examples are given.

Aubrun, J. N.↗
NASA NTRS1981

Space assembly methodology

Large space structure assembly analysis techniques are defined and simulation activities are described. The simulations included are: an extravehicular activity assembly simulation; a fabricated beam assembly series using a beam generating machine end caps, and cross beam brackets; deployment of a deployable truss, using the neutral buoyancy remote manipulator system with crewman assistance; and a series aboard the KC-135 zero g aircraft.

Stokes, J. W.↗
NASA NTRS1987

Electron stimulated desorption of atomic oxygen from silver

The electron stimulated desorption (ESD) of neutral oxygen atoms from polycrystalline silver and of oxygen ions from Ag(110) has been studied. Polycrystalline Ag charged with (16)O2 and (18)O2 and bombarded by low-energy electrons (approx 100 eV) under ultrahigh vacuum (UHV) conditions emitted O atom flux levels of 1 x 10 to the 12th power/sq cm/s at a Ag temperature of 300 C. The flux was detected with a quadrupole mass spectrometer operating in the appearance potential mode. The neutral cross section at about 100 C was determined to be 7 x 10 to the -19 sq cm. Ancillary experiments conducted in a UHV chamber equipped with a cylindrical mirror analyzer and rigged for ion energy distribution and ion angular distribution were used to study O ions desorbed from Ag(110). Two primary O(+) energies of 2.4 and 5.4 eV were detected from the Ag(110) after having been dosed with 2500 L of (16)O2. It also appears that in both experiments there was strong evidence for directionality of the emitted flux. The results of this study serve as a proof of concept for the development of a laboratory atomic oxygen beam generator that simulates the gas flux environment experienced by orbiting vehicles.

Outlaw, R. A.↗