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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Electron Modulational Instability in the Strong Turbulent Regime for an Electron Beam Propagating in a Background Plasma

We study collective processes for an electron beam propagating through a background plasma using simulations and analytical theory. A new regime where the instability of a Langmuir wave packet can grow locally much faster than ion frequency is clearly identified. The key feature of this new regime is an electron modulational instability that rapidly creates a local Langmuir wave packet, which in its turn produces local charge separation and strong ion density perturbations because of the action of the ponderomotive force, such that the beam-plasma wave interaction stops being resonant. Three evolution stages of the process and observed periodic burst features are discussed. In conclusion, different physical regimes in the plasma and beam parameter space are demonstrated for the first time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Angular distribution of neutron production by proton and carbon-ion therapeutic beams

Carbon-ion beams are increasingly used in the clinical practice for external radiotherapy treatments of deep-seated tumors. At therapeutic energies, carbon ions yield significant secondary products, including neutrons, which may be of concern for the radiation protection of the patient and personnel. We simulated the neutron yield produced by proton and carbon-ion pencil beams impinging on a clinical phantom at three different angles: 15°, 45° and 90°, with respect to the beam axis. We validated the simulated results using the measured response of organic scintillation detectors. We compared the results obtained with FLUKA 2011.2 and MCNPX 2.7.0 based on three different physics models: Bertini, Isabel, and CEM. Over the different ions, energies, and angles, the FLUKA simulation results agree better with the measured data, compared to the MCNPX results. Simulations of carbon ions at low angles exhibit both the highest deviation from measured data and inter-model discrepancy, which is probably due to the different treatment of the pre-equilibrium stage. The reported neutron yield results could help in the comparison of carbon-ion and proton treatments in terms of secondary neutron production for radiation protection applications.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Laser-to-proton conversion efficiency studies for proton fast ignition

We assess the conversion efficiency from intense picosecond laser pulses to multi-MeV ion beams for a wide range of laser and target parameters, using 2D kinetic particle-in-cell simulations. Scalings are addressed in a quasi-one-dimensional geometry, leaving out beam divergence. Then, we study the conversion efficiency into a narrow spatial band along the laser axis for flat and hemispherical targets in large-scale 2D simulations. Combining these findings allows us to calculate the energy requirements for ignition of a compressed ICF target with an intense proton beam in a fast-ignition scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Sample spinning to mitigate polarization artifact and interstitial-vacancy imbalance in ion-beam irradiation

Accelerator-based ion-beam irradiation has been widely used to mimic the effects of neutron radiation damage in nuclear reactors. However, ion radiation is most often monodisperse in the incoming ions’ momentum direction, leading to excessive polarization in defect distribution, while the scattering under neutron irradiation is often more isotropic and has less radiation-induced polarization. Mitigation of the excess-polarization as well as the damage non-uniformity artifact might be crucial for making the simulation of neutron radiation by ion-beam radiation more realistic. In this work, a general radiation polarization theory in treating radiation as external polar stimuli is established to understand the natural material responses in different contexts, and the possibility to correct the defect polarization artifact in ion-beam irradiation. Inspired by Magic Angle Spinning in Nuclear Magnetic Resonance, we present a precise sample spinning strategy to reduce the point-defect imbalance effect in ion-beam irradiation. It can be seen that with optimized surface inclination angle and the axis of sample rotation, the vacancy-interstitial population imbalance, as well as the damage profile non-uniformity in a designated region in the target are both reduced. It is estimated that sample spinning frequency on the order of kHz should be sufficient to scramble the ion momentum monodispersity for commonly taken ion fluxes and dose rates, which is experimentally feasible.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Transport and losses of fusion-born alpha particles in the presence of tearing modes using the new Toroidal Accelerated Particle Simulator (TAPaS)

Here, the transport and losses of fusion-born alpha particles are studied in the presence of a single-helicity tearing mode, characterized by (m = 2, n = 1). The analysis is performed by means of the recently developed Toroidal Accelerated Particle Simulator (TAPaS). Although such modes have usually been believed to result only in a local flattening of the radial profiles, it is shown that the density profile can exhibit a global modification leading to significant losses of alpha particles. This is due to the fact that, although the magnetic field does not exhibit any chaotic behaviour, the trajectories of alpha particles do, as revealed by their Poincaré maps. Such results are in qualitative agreement with past observations and simulations of energetic particles generated by neutral beam injection in TFTR, DIII-D and AUG tokamaks. In-depth analysis is carried out to characterize the impact of the tearing mode on the transport and losses of fusion-born alpha-particles with a realistic density profile. The impact of the amplitude is evidenced. Moreover, the effect of the island rotation frequency is assessed based on a detailed analysis of the linear resonances in phase-space, in agreement with the simulation results. Finally, the probability density function of the exit time has been computed and the transport of alpha particles has been found to be anomalous.

anomalous transport↗

DOE OSTI2024

Numerical modeling of a proof-of-principle experiment on optical stochastic cooling at an electron storage ring

Cooling of beams circulating in storage rings is critical for many applications including particle colliders and synchrotron light sources. A method enabling unprecedented beam-cooling rates, optical stochastic cooling (OSC), was recently demonstrated in the Integrable Optics Test Accelerator (IOTA) electron storage ring at Fermilab [J. Jarvis , ]. This paper describes the numerical implementation of the OSC process in the particle-tracking program and discusses the validation of the developed model with available experimental data. The model is also employed to highlight some features associated with different modes of operation of OSC. The developed simulation tool should be valuable in guiding future configurations of optical stochastic cooling and, more broadly, modeling self-field-based beam manipulations. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

DOE OSTI2023

FAST Low-Energy Beamline Studies: Toward High-Peak 5D Brightness Beams for FAST-GREENS

The FAST beamline is the injector for the planned Gamma-Ray Electron ENhanced Source (GREENS) program, which aims to achieve the demonstration and first application of a high-efficiency, high-average-power free-electron laser at 515 nm. FAST-GREENS requires high 5D peak brightness; transverse normalized projected emittances of 3 mm-mrad and a peak current of 600 A are the minimum beam requirements for the FEL to reach the 10% efficiency goal. In this work, studies of the low-energy section of the FAST beamline are presented toward these ends, including preliminary measurements of beam compression and beam emittance. An effort toward developing a high-fidelity simulation model that could be later optimized for FAST-GREENS is presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

A Simulation of the Photoionization of H- Together With the Subsequent Tracking of the Liberated Electrons

The Proton Improvement Plan - II (PIP-II) is a new linear accelerator (LINAC) complex being built at Fermilab. It is based on superconducting radiofrequency cavities and will accelerate H- ions to 800 MeV kinetic energy before injection into the existing Booster ring. Measurements of the profile of the beam along the LINAC must be done by non-intercepting methods due to the superconducting cavities. The method chosen is photoionization of a small number of H- by a focused infrared laser, aka laserwire. The number of ionized electrons is measured as a function of laser position within the H- beam. To aid in the design of the collection mechanism, a simulation was written in MATLAB with input from the commercial electromagnetic simulation, CST. This simulation calculates the number and positions of the liberated electrons and tracks them through the magnetic collection and H- beam fields to the collection point. Results from this simulation for various points along the LINAC will be shown.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

A Simulation of the Photoionization of H- Together with the Subsequent Tracking of the Liberated Electrons

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

RaDIATE Collaboration Thermal shock studies

As next-generation accelerator target facilities (High Energy Physics, Spallation Sources, ... ) become increasingly more powerful and intense, high power target systems face key technical challenges, such as radiation damage and thermal shock. Those ultimately degrade the performance and lifetime of targets and have been identified as the leading cross-cutting challenges of high-power target facilities. In order to operate reliable beam-intercepting devices in the framework of energy and intensity increase for next generation accelerators, the RaDIATE Collaboration (Radiation Damage In Accelerator Target Environment), established in 2012 and managed by Fermilab, brings together existing expertise in nuclear material and accelerator targets from 20 international institutions, including CERN, to execute a coordinated strategy for high power targetry R&D. In this context, several thermal shock studies were performed at CERN's HiRadMat (High-Radiation to Materials) facility, that took a key step towards improving our knowledge on target damage tolerance. The first experiment HRMT-24, supported by EURCARD2 and completed in 2015, successfully validated the Johnson-Cook strength model developed at SwRI on beryllium S200FH (used for beam window material), providing a better confidence in simulating the thermal shock response of current and future S200FH beryllium components. HRMT-43, supported by ARIES and completed in 2018, tested various materials (Be, C, SiC, Si, Ti and ceramic nanofiber). It was a first and unique test which included pre-irradiated specimens from high energy proton beam irradiation to identify thermal shock response differences between non-irradiated and previously irradiated materials. Real-time measurement of dynamic thermomechanical response of graphite helped to benchmark numerical simulations.

43 PARTICLE ACCELERATORS↗

DOE OSTI2022

MicroBooNE BNB Inclusive Overlay Sample (With Wire Info)

MicroBooNE samples are provided for collaborative development in two different formats: HDF5, targeting the broadest audience, and artroot, targeting users that are familiar with the software infrastructure of Fermilab neutrino experiments and more in general of HEP experiments. The HDF5 files are stored on Zenodo, together with a list of artroot files accessible with xrootd.This sample includes simulated interactions of neutrinos from the Booster Neutrino Beam (BNB), overlaid on top of cosmic ray data. The sample is inclusive, i.e. it includes all types of neutrinos and interactions, with relative abundance matching our nominal flux and cross section models. Interactions are simulated in in the whole cryostat volume.The HDF5 files in this sample include the information at the wire waveform level (after deconvolution and finding of regions of interest). As this information significantly increases the file size, this sample contains about 17% of the events of the corresponding sample without wire information.More documentation, including detailed description of content, recipes, and example usage, at https://github.com/uboone/OpenSamples.Suggested text for acknowledgment is the following:We acknowledge the MicroBooNE Collaboration for making publicly available the data sets [data set DOIs] employed in this work. These data sets consist of simulated neutrino interactions from the Booster Neutrino Beamline overlaid on top of cosmic data collected with the MicroBooNE detector [2017 JINST 12 P02017].In addition, we request that software products resulting from the usage of the datasets are also made publicly available.

Abratenko, Polina↗

DOE OSTI2022

MicroBooNE BNB Inclusive Overlay Sample (No Wire Info)

MicroBooNE samples are provided for collaborative development in two different formats: HDF5, targeting the broadest audience, and artroot, targeting users that are familiar with the software infrastructure of Fermilab neutrino experiments and more in general of HEP experiments. The HDF5 files are stored on Zenodo, together with a list of artroot files accessible with xrootd.This sample includes simulated interactions of neutrinos from the Booster Neutrino Beam (BNB), overlaid on top of cosmic ray data. The sample is inclusive, i.e. it includes all types of neutrinos and interactions, with relative abundance matching our nominal flux and cross section models. Interactions are simulated in in the whole cryostat volume.The HDF5 files in this sample do not include the information at the wire waveform level ("NoWire" label), allowing for larger number of events to be included in the data set.More documentation, including detailed description of content, recipes, and example usage, at https://github.com/uboone/OpenSamples.Suggested text for acknowledgment is the following:We acknowledge the MicroBooNE Collaboration for making publicly available the data sets [data set DOIs] employed in this work. These data sets consist of simulated neutrino interactions from the Booster Neutrino Beamline overlaid on top of cosmic data collected with the MicroBooNE detector [2017 JINST 12 P02017].In addition, we request that software products resulting from the usage of the datasets are also made publicly available.

Abratenko, Polina↗

DOE OSTI2025

Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC

Plasma-based acceleration (PBA) has emerged as a promising candidate for the accelerator technology used to build a future linear collider and/or an advanced light source. In PBA, a trailing or witness particle beam is accelerated in the plasma wave wakefield (WF) created by a laser or particle beam driver. The WF is often nonlinear and involves the crossing of plasma particle trajectories in real space and thus particle-in-cell methods are used. The distance over which the drive beam evolves is several orders of magnitude larger than the wake wavelength. This large disparity in length scales is amenable to the quasi-static approach. Three-dimensional (3D), quasi-static (QS), particle-in-cell (PIC) codes, e.g., QuickPIC, have been shown to provide high fidelity simulation capability with 2-4 orders of magnitude speedup over 3D fully explicit PIC codes. In PBA, the witness beam needs to be matched to the focusing forces of the WF to reduce the emittance growth. In some linear collider designs, the matched spot size of the witness beam can be 2 to 3 orders of magnitude smaller than the spot size (and wavelength) of the wakefield. Such an additional disparity in length scales is ideal for mesh refinement where the WF within the witness beam is described on a finer mesh than the rest of the WF. A mesh refinement scheme is described that has been implemented into the 3D QS PIC code, QuickPIC. Very fine (high) resolution is used in a small spatial region that includes the witness beam and progressively coarser resolutions in the rest of the simulation domain. A fast multigrid Poisson solver has been implemented for the field solve on the refined meshes and a Fast Fourier Transform (FFT) based Poisson solver is used for the coarse mesh. The code has been parallelized with both MPI and OpenMP, and the parallel scalability has also been improved by using pipelining. A preliminary adaptive mesh refinement technique is described to optimize the computational time for simulations with an evolving witness beam size. Several test problems are used to verify that the mesh refinement algorithm provides accurate results. Additionally, the results are benchmarked against highly resolved simulations exhibiting near-azimuthal symmetry, performed using QPAD—a novel hybrid QS PIC code that uses a PIC description in the coordinates (r, ct – z) and a gridless description in the azimuthal angle, Φ.

Linear collider↗

DOE OSTI2021

As-Built Simulation of the High Flux Isotope Reactor

The Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) is an 85 MWt flux trap-type research reactor that supports key research missions, including isotope production, materials irradiation, and neutron scattering. The core consists of an inner and an outer fuel element containing 171 and 369 involute-shaped plates, respectively. The thin fuel plates consist of a U 3 O 8 -Al dispersion fuel (highly enriched), an aluminum-based filler, and aluminum cladding. The fuel meat thickness is varied across the width of the involute plate to reduce thermal flux peaks at the radial edges of the fuel elements. Some deviation from the designed fuel meat shaping is allowed during manufacturing. A homogeneity scan of each fuel plate checks for potential anomalies in the fuel distribution by scanning the surface of the plate and comparing the attenuation of the beam to calibration standards. While typical HFIR simulations use homogenized fuel regions, explicit models of the plates were developed under the Low-Enriched Uranium Conversion Program. These explicit models typically include one inner and one outer fuel plate with nominal fuel distributions, and then the plates are duplicated to fill the space of the corresponding fuel element. Therefore, data extracted from these simulations are limited to azimuthally averaged quantities. To determine the reactivity and physics impacts of an as-built outer fuel element and generate azimuthally dependent data in the element, 369 unique fuel plate models were generated and positioned. This model generates the three-dimensional (i.e., radial–axial–azimuthal) plate power profile, where the azimuthal profile is impacted by features within the adjacent control element region and beryllium reflector. For an as-built model of the outer fuel element, plate-specific homogeneity data, 235 U loading, enrichment, and channel thickness measurements were translated into the model, yielding a much more varied azimuthal power profile encompassed by uncertainty factors in analyses. These models were run with the ORNL-TN and Shift Monte Carlo tools, and they contained upwards of 500,000 cells and 100,000 unique tallies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

Modal Analysis and Testing of a Cantilever Beam: Results and Guide

This report discusses modal analysis and testing of a cantilever beam including an introduction to modal testing, simulation, analysis, and a comparison of results. An accompanying presentation can be found from the LLNL CASIS website. Modal testing and analysis are essential for characterizing the response of a structure during vibration environments. Determining natural frequencies and mode shapes helps determine whether resonances will be reached during operations and if so, how the structure will respond and how the response can be tuned. The goal of this work is to convey the basics of modal testing and analysis through a simple project. This can be applied later to more complex, applicable structures. This includes a variety of work including data collection, finite element analysis and data processing. For this project a cantilever beam was chosen because it is a simple, well-characterized structure.

42 ENGINEERING↗

DOE OSTI2023

Laser-direct-drive fusion target design with a high- Z gradient-density pusher shell

Laser-direct-drive fusion target designs with solid deuterium-tritium (DT) fuel, a high-Z gradient-density pusher shell (GDPS), and a Au-coated foam layer have been investigated through both 1D and 2D radiation-hydrodynamic simulations. Compared with conventional low-Z ablators and DT-push-on-DT targets, these GDPS targets possess certain advantages of being instability-resistant implosions that can be high adiabat (α ≥ 8) and low hot-spot and pusher-shell convergence (CR hs ≈ 22 and CR PS ≈ 17), and have a low implosion velocity (v imp < 3 × 10 7 cm/s). Using symmetric drive with laser energies of 1.9 to 2.5 MJ, 1D LILAC simulations of these GDPS implosions can result in neutron yields corresponding to >~50–MJ energy, even with reduced laser absorption due to the cross-beam energy transfer (CBET) effect. Two-dimensional DRACO simulations show that these GDPS targets can still ignite and deliver neutron yields from 4 to ~10 MJ even if CBET is present, while traditional DT-push-on-DT targets normally fail due to the CBET-induced reduction of ablation pressure. If CBET is mitigated, these GDPS targets are expected to produce neutron yields of >20 MJ at a driven laser energy of ~2 MJ. Here are the key factors behind the robust ignition and moderate energy gain of such GDPS implosions are as follows: (1) The high initial density of the high-Z pusher shell can be placed at a very high adiabat while the DT fuel is maintained at a relatively low-entropy state; therefore, such implosions can still provide enough compression ρR >1 g/cm 2 for sufficient confinement; (2) the high-Z layer significantly reduces heat-conduction loss from the hot spot since thermal conductivity scales as ~1/Z; and (3) possible radiation trapping may offer an additional advantage for reducing energy loss from such high-Z targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

JSPEC – An Open Source Program for IBS and Electron Cooling Simulation

The intrabeam scattering can affect the accumulation, the lifetime, and the property of a high-intensity beam. Electron cooling is a method to mitigate the intrabeam scattering effect. JSPEC (JLab Simulation Package on Electron Cooling) is an open-source program developed at Jefferson Lab, which includes various numerical models and friction force formulas for intrabeam scattering and electron cooling simulations. JSPEC has been benchmarked with BETACOOL and experimental data. In this report, we will introduce what features JSPEC provides to the users and how it carries out the computations. Numerical examples are presented to demonstrate the performance and the validity of JSPEC.

Zhang, He↗

DOE OSTI2025

Steerable terahertz beams using surface waves on an active metasurface

The development of dynamic components for controlling wave fronts in the sub-terahertz region of the electromagnetic spectrum has emerged as a frontier research topic for many applications in sensing and communications. One approach which has attracted much attention involves the use of active metasurfaces, tiled arrays of sub-wavelength elements with properties that can be reconfigured via external actuation. In nearly all cases, these metasurfaces are employed as either transmissive or reflective elements, taking advantage of their strong and tunable interaction with free-space electromagnetic waves. These interactions can be significantly enhanced through the use of surface waves propagating parallel to the metasurface array, although very few studies have exploited this option. Here, we integrate a metasurface into the interior of a parallel-plate waveguide in a configuration explicitly designed to exploit this surface-wave geometry. We show that varying the electrical properties of the active metasurface changes the wave vector of the guided mode, and thereby alters the emission angle of radiation out-coupled through a leaky-wave slot aperture. These results, which are consistent with numerical simulations, represent a new approach to broadband beam steering suitable for the sub-terahertz spectral range.

77 NANOSCIENCE AND NANOTECHNOLOGY↗