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 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Fermilab PIP-II machine protection system digitized data noise elimination scheme and its FPGA implementation

In Fermilab's PIP-II machine protection system, beam loss signals from various detectors are digitized at 125 MS/s. Noise from both high-frequency sources and low-frequency 60 Hz AC power equipment can contaminate the data. To suppress noise across these ranges—especially 60 Hz and its harmonics, which overlap with beam loss signal frequencies—advanced digital processing beyond standard filtering is required. Several real-time functional blocks were simulated and tested on an FPGA: (1) a dual time-constant discharging integrator filter, (2) a de-ripple baseline extraction and storage block, and (3) a fast-recovery discharging integrator. The nonlinear IIR integrator filter removes high-frequency noise and feeds into the baseline extractor. Upon detecting abrupt beam loss, it switches to a longer time constant to prevent baseline distortion. The de-ripple block calculates a valid baseline by averaging over multiple 60 Hz periods, storing results in a 4096-word FPGA RAM. This baseline is subtracted from raw data before integration by the fast-recovery block, which resets quickly after use. All blocks achieved expected performance and were successfully implemented on a low-cost FPGA.

Wu, Jinyuan [Fermilab]↗
DOE OSTI2022

Benchmarking and Exploring Parameter Space of the 2-Phase Bubble Tracking Model for Liquid Mercury Target Simulation

High intensity proton pulses strike the Spallation Neutron Source (SNS)’s mercury target to provide bright neutron beams. These strikes deposit extensive energy into the mercury and its steel vessel. Prediction of the resultant loading on the target is difficult when helium gas is intentionally injected into the mercury to reduce the loading and to mitigate the pitting damage on the vessel. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this complex multi-physics dynamics problem by including the bubble dynamics in the liquid mercury. We present a study comparing the measured target strains in the SNS target station with the simulation results of the solid mechanics simulation framework. We investigate a wide range of various physical model parameters, including the number of bubble families, bubble size distribution, viscosity, surface tension, etc. to understand their impact on simulation accuracy. Our initial findings reveal that using 8-10 bubble families in the model renders a simulation strain envelope that covers the experimental ones. Further optimization studies are planned to predict the strain response more accurately.

Lin, Lianshan↗
NASA NTRS1982

Numerical simulation of the flow and fuel-air mixing in an axisymmetric piston-cylinder arrangement

The implicit factored method of Beam and Warming was employed to describe the flow and the fuel-air mixing in an axisymmetric piston-cylinder configuration during the intake and compression strokes. The governing equations were established on the basis of laminar flow. The increased mixing due to turbulence was simulated by appropriately chosen effective transport properties. Calculations were performed for single-component gases and for two-component gases and for two-component gas mixtures. The flow field was calculated as functions of time and position for different geometries, piston speeds, intake-charge-to-residual-gas-pressure ratios, and species mass fractions of the intake charge. Results are presented in graphical form which show the formation, growth, and break-up of those vortices which form during the intake stroke and the mixing of fuel and air throughout the intake and compression strokes. It is shown that at bore-to-stroke ratio of less than unity, the vortices may break-up during the intake stroke. It is also shown that vortices which do not break-up during the intake stroke coalesce during the compression stroke. The results generated were compared to existing numerical solutions and to available experimental data.

Shih, T. I. P.↗
NASA NTRS1993

Simulation of TunneLadder traveling-wave tube cold-test characteristics: Implementation of the three-dimensional, electromagnetic circuit analysis code micro-SOS

The three-dimensional, electromagnetic circuit analysis code, Micro-SOS, can be used to reduce expensive time-consuming experimental 'cold-testing' of traveling-wave tube (TWT) circuits. The frequency-phase dispersion characteristics and beam interaction impedance of a TunneLadder traveling-wave tube slow-wave structure were simulated using the code. When reasonable dimensional adjustments are made, computer results agree closely with experimental data. Modifications to the circuit geometry that would make the TunneLadder TWT easier to fabricate for higher frequency operation are explored.

Kory, Carol L.↗
DOE OSTI2024

Simulation of a scintillator-based fast ion loss detector for steady-state operation in Wendelstein 7-X (invited)

A quantitative theoretical framework has been created to model neutral beam injection and fast ion losses in the Wendelstein 7-X (W7-X) stellarator, including a novel method to develop synthetic diagnostics for fast ion loss detectors (FILDs) of many types, such as scintillating and Faraday Cup FILDs. This is the first time that this has been done in stellarator geometry with this level of fidelity, providing a way for fast ion losses to be predicted more precisely in future stellarator experiments and in W7-X. Simulations of the signal seen by a Faraday Cup FILD have been completed for multiple W7-X plasmas and show close agreement with the measured signals. This method is now applied to an actively water-cooled, scintillator-based FILD, which is currently in development to measure the fast ion loss distribution in W7-X in greater detail. The design makes use of a double slit to measure energy-and-pitch-angle-resolved losses of both co-going and counter-going fast ions. The diagnostic, which can be inserted to different radial positions, has been designed to withstand steady-state heat fluxes of up to 120 kW/m 2 along with additional transient heat loads of 100 kW/m 2 lasting for up to 20 s at a time. Simulations of W7-X standard magnetic configuration show up to 8 × 10 13 (s –1 cm –2 ) ion fluxes onto the sensor from each neutral beam source and no signal from the counter-going slit. Finally, these simulations will help inform experimental proposals for future W7-X campaigns after installation of this diagnostic.

3D printing↗
DOE OSTI2024

Performance prediction applying different reduced turbulence models to the SMART tokamak

The SMall Aspect Ratio Tokamak (SMART) is currently being commissioned at the University of Seville and will be able to compare the performance of positive and negative triangularity plasmas at low aspect ratio. Predictive simulations have been performed for different machine scenarios and heating schemes using the TRANSP code. The objectives of these simulations are to predict the parameters expected in positive triangularity plasmas, to guide diagnostic development, and to validate transport models. Several reduced turbulence models have been used to predict electron and ion temperatures for the operational phase 2. All models provide similar results from approximately mid-radius to the separatrix but important discrepancies are found in the core region. These positive triangularity results are compared with experiments from a similar size machine like GLOBUS-M2. The multi-mode model (MMM) shows the best agreement. Simulations with different boundary conditions have been performed and no strong differences have been observed between them. The impact of neutral beam injection (NBI) on the predicted profiles has also been addressed. Rotation reduces turbulence levels so higher temperatures are achieved when included in the simulations. Studying the different contributions to the thermal diffusivities, it is observed that electron temperature gradient (ETG) turbulence dominates at the plasma core while micro-tearing modes (MTM) dominate at the edge in the electron channel. In the ion channel, the neoclassical contribution is dominant at the core and at the very edge while the Weiland component, which includes ion temperature gradient mode (ITG), trapped electron mode (TEM), kinetic ballooning mode (KBM), peeling mode (PM) and collisionless and collision dominated magnetohydrodynamic (MHD) modes governs the mid-radius region. For phase 3, two plasmas with different electron densities have been studied. The case with lower density matches well a specific discharge of GLOBUS-M2. The higher density plasma shows high performance with β N ≈ 3.8.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Spinning disk calibration method and apparatus for laser Doppler velocimeter

A method and apparatus for calibrating laser Doppler velocimeters having one or more intersecting beam pairs are described. These velocimeters measure fluid velocity by observing the light scattered by particles in the fluid stream. Moving fluid particulates are simulated by fine taut wires that are radially mounted on a disk that is rotated at a known velocity. The laser beam intersection locus is first aimed at the very center of the disk and then the disk is translated so that the locus is swept by the rotating wires. The radial distance traversed is precisely measured so that the velocity of the wires (pseudo particles) may be calculated.

Snyder, P. K.↗
NASA NTRS2012

A New Electron Source for Laboratory Simulation of the Space Environment

We have developed a new collimated electron source called the Photoelectron Beam Generator (PEBG) for laboratory and spaceflight applications. This technology is needed to replace traditional cathodes because of serious fundamental weaknesses with the present state of the art. Filament cathodes suffer from numerous practical problems, even if expertly designed, including the dependence of electron emission on filament temperature, short lifetimes (approx 100 hours), and relatively high power (approx 10s of W). Other types of cathodes have solved some of these problems, but they are plagued with other difficult problems, such as the Spindt cathode's extreme sensitivity to molecular oxygen. None to date have been able to meet the demand of long lifetime, robust packaging, and precision energy and flux control. This new cathode design avoids many common pitfalls of traditional cathodes. Specifically, there are no fragile parts, no sensitivity to oxygen, no intrinsic emission dependencies on device temperature, and no vacuum requirements for protecting the source from contamination or damage. Recent advances in high-brightness Light Emitting Diodes (LEDs) have provided the key enabling technology for this new electron source. The LEDs are used to photoeject electrons off a target material of a low work-function, and these photoelectrons are subsequently focused into a laminar beam using electrostatic lenses. The PEBG works by illuminating a target material and steering photoelectrons into a laminar beam using electrostatic lenses

Krause, Linda Habash↗
DOE OSTI2022

X-ray FEL linear accelerator design via start-to-end global optimization

An X-ray Free Electron Laser (FEL) prefers using an electron beam with low emittance, small energy spread, and a high core current to generate coherent radiation through an undulator. In order to attain such a high brightness beam, the linear accelerator beam dynamics design generally involves separate photoinjector optimization and linac optimization. In this paper, we propose a new beam dynamics design strategy based on global optimization with fast start-to-end simulations from the photocathode to the end of the accelerator. The new start-to-end model significantly reduces the simulation time and makes the global optimization practical. The global optimization method avoids the need to choose a single solution based on bunch length at the injector exit for the linac optimization and helps find the solution with unfavorable bunch length at the injector exit but better phase space distribution that can result in better final electron beam phase space distribution at the entrance of the undulator. Here, using the start-to-end global optimization, we showed in an application example, with a 100 pC beam that good transverse emittance and over kilo-Ampere final core current can be attained using a photoinjector that consists of a VHF gun and boosting RF cavities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Particle-in-cell techniques for the study of space charge effects in the Advanced Cryogenic Gas Stopper

The Facility for Rare Isotope Beams, FRIB, will significantly increase the rates and varieties of rare isotopes compared to those presently available at existing rare-isotope facilities. The production of stopped beams via the current generation of linear gas stopping devices available is not well matched to the beam intensities that will be available at FRIB. Many experiments that make use of reaccelerated beams, from the ReA facility, greatly benefit from having the highest-intensity, rare-isotope beams possible. Some experiments would also benefit greatly from stopped beams with improved purity. In this paper, the intensity bottlenecks of FRIB’s newest linear gas stopper are studied via experiment and simulation in support of a future detailed conceptual design of a next-generation linear gas stopper that can provide rare-isotope beams with rates of 10^8 pps, a factor of 50 higher than the current state-of-the-art system, with improved beam purity.

43 PARTICLE ACCELERATORS↗
NASA NTRS1971

The interaction of O(plus) ions with the interior surface of a copper chamber

Modulated beams of 0(+), Ar(+), and Kr(+) in the 100-300 eV range are directed into a copper box simulating the ante-chamber of an orbiting mass spectrometer. An RF quadrupole mass spectrometer and phase sensitive detection extract the component of the internal mass spectrum correlated with the beam. Intense Ar and Kr signals are observed; however, no O or O2 is detectable, indicating loss of the primary O(+) beam to surface interactions. All four primary ions stimulate sizeable signals at masses 26 and 28. The relevance of these experiments to the interpretation of mass spectra obtained by orbiting satellites is discussed.

Siegel, M. W.↗
DOE OSTI2024

Validation of the NOνA experiment 2023-tuning on simulated neutrino-matter interactions

NOνA is a long-baseline neutrino oscillation experiment that utilizes a two-detector design to study the oscillations of muon neutrinos into electron neutrinos over a baseline of 810 km. The Near Detector (ND) measures the neutrino beam spectrum and composition before oscillation, which is then compared to the oscillated neutrino energy spectrum observed in the Far Detector (FD). In the ND, the neutrinos are detected through their interactions with the heavy target nuclei within the detector. NOνA employs the GENIE neutrino event generator for simulating these neutrino-nucleus interactions. However, the default GENIE prediction does not adequately reproduce the ND data. To address this, NOνA developed a tune of the neutrino interaction models within GENIE version 3:0:6 to minimize discrepancies between the simulated predictions and the observed data in the ND. This dissertation tests the NOνA’s 2023 tune of the GENIE neutrino cross-section simulations by performing a data/simulations comparison for the ND. The analysis employed datasets comprising $2.55\times10^{21}$ protons-on-target (POT) in neutrino beam mode and $1.14\times10^{21}$ POT in antineutrino beam mode. The NOνA tuning of neutrino-matter interaction simulations matches with ND data within the $1\sigma$ error band, overestimating muon neutrino and antineutrino charged current interactions by approximately 6 % and 9 %, respectively. Discrepancies were observed in the energy region dominated by Quasi-Elastic-like interactions. Systematic uncertainties associated with the modeling of the neutrino cross-section, particularly those pertaining Quasi-Elastic like interactions, contributed considerably to the overall error in the simulations. Furthermore, the reconstruction algorithm used in NOνA for particle classification demonstrated significant misidentifications between charged pions and protons, as well as a tendency to overlook additional pions or protons in multi-particle simulated events.

Cortés Parra, Camilo Andrés↗
NASA NTRS1987

Time- and position-varying infrared scene simulation

Rapidly changing infrared scenes can be generated as a function of position using thermally produced irradiance. This can be done most conveniently by using a variable dwell-time, variable diameter, and scanning laser beam to deposit a predetermined amount of laser energy, as a function of position, on the target surface. An infrared target simulator has been designed and tested to perform within stringent temperature simulation requirements: heating rate of over 100 K/sec is achieved, and temperature gradients of 200 K/mm can be maintained. Problems with heat transfer between adjacent elements at different temperatures are solved by modeling the thermal layer on the target surface as a low-pass filter, which cuts off all the frequencies above the thermal cutoff frequency. When the system approach is taken to the infrared simulation design, the optical system limited performance results in the best realizable performance. In this case, the thermal layer on the target surface is designed so that the thermal cutoff frequency is above the optical cutoff frequency.

Scholl, Marija S.↗
DOE OSTI2022

Whistler Waves Associated With Electron Beams in Magnetopause Reconnection Diffusion Regions

Abstract Whistler waves are often observed in magnetopause reconnection associated with electron beams. We analyze seven MMS crossings surrounding the electron diffusion region (EDR) to study the role of electron beams in whistler excitation. Waves have two major types: (a) Narrow‐band waves with high ellipticities and (b) broad‐band waves that are more electrostatic with significant variations in ellipticities and wave normal angles. While both types of waves are associated with electron beams, the key difference is the anisotropy of the background population, with perpendicular and parallel anisotropies, respectively. The linear instability analysis suggests that the first type of wave is mainly due to the background anisotropy, with the beam contributing additional cyclotron resonance to enhance the wave growth. The second type of broadband waves are excited via Landau resonance, and as seen in one event, the beam anisotropy induces an additional cyclotron mode. The results are supported by particle‐in‐cell simulations. We infer that the first type occurs downstream of the central EDR, where background electrons experience Betatron acceleration to form the perpendicular anisotropy; the second type occurs in the central EDR of guide field reconnection. A parametric study is conducted with linear instability analysis. A beam anisotropy alone of above ∼3 likely excites the cyclotron mode waves. Large beam drifts cause Doppler shifts and may lead to left‐hand polarizations in the ion frame. Future studies are needed to determine whether the observation covers a broader parameter regime and to understand the competition between whistler and other instabilities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

DUNE – Simulation Validation of Fermilab Detector Reconstruction

DUNE (Deep Underground Neutrino Experiment) is Fermilab’s flagship international experiment designed to study neutrinos by sending an intense beam from Illinois to detectors located 1,300 kilometers away at the Sanford Underground Research Facility (SURF) in South Dakota. To prepare for such a large-scale experiment, physicists develop detailed simulations to produce mock data sets which are analyzed by the CAFAna framework. During my internship, I developed software using the CAFAna framework to analyze simulated detector data and generated plots to make data trends easier to interpret and identify patterns. My analysis has uncovered inconsistencies in reconstructed neutrino tracks, duplicated reconstructed tracks causing sporadic spikes in the data, and unnatural differences in energy levels between interaction types. These analyses help verify that the improvements to detector simulations do not introduce unintended resolution errors and ensure proper reconstruction performance, supporting DUNE’s goal of making precise neutrino measurements and advancing the Department of Energy’s mission of fundamental scientific discovery.

Vershaw, Andre [Unlisted, US, IL; Fermilab] (ORCID↗
DOE OSTI2024

Monte Carlo Simulation of Electron Interactions in an MeV-STEM for Thick Frozen Biological Sample Imaging

A variety of volume electron microscopy techniques have been developed to visualize thick biological samples. However, the resolution is limited by the sliced section thickness (>30–60 nm). To preserve biological samples in a hydrated state, cryo-focused ion beam scanning electron microscopy has been developed, providing nm resolutions. However, this method is time-consuming, requiring 15–20 h to image a 10 μm thick sample with an 8 nm slice thickness. There is a pressing need for a method that allows the rapid and efficient study of thick biological samples while maintaining nanoscale resolution. The remarkable ability of mega-electron-volt (MeV) electrons to penetrate thick biological samples, even exceeding 10 μm in thickness, while maintaining nanoscale resolution, positions MeV-STEM as a suitable microscopy tool for such applications. Our research delves into understanding the interactions between MeV electrons and frozen biological specimens through Monte Carlo simulations. Single elastic scattering, plural elastic scattering, single inelastic scattering, and plural inelastic scattering events have been simulated. The electron trajectories, the beam profile, and the intensity change of electrons in each category have been investigated. Additionally, the effects of the detector collection angle and the focal position of the electron beam were investigated. As electrons penetrated deeper into the specimen, single and plural elastic scattered electrons diminished, and plural inelastic scattered electrons became dominant, and the beam profile became wider. Even after 10 μm of the specimen, 42% of the MeV electrons were collected within 10 mrad. This confirms that MeV-STEM can be employed to study thick biological samples.

electron specimen interaction↗
DOE OSTI2025

Results of a Geant4 benchmarking study for bio‐medical applications, performed with the G4‐Med system

Geant4, a Monte Carlo Simulation Toolkit extensively used in bio-medical physics, is in continuous evolution to include newest research findings to improve its accuracy and to respond to the evolving needs of a very diverse user community. In 2014, the G4-Med benchmarking system was born from the effort of the Geant4 Medical Simulation Benchmarking Group, to benchmark and monitor the evolution of Geant4 for medical physics applications. The G4-Med system was first described in our Medical Physics Special Report published in 2021. Results of the tests were reported for Geant4 10.5. Purpose In this work, we describe the evolution of the G4-Med benchmarking system. Methods The G4-Med benchmarking suite currently includes 23 tests, which benchmark Geant4 from the calculation of basic physical quantities to the simulation of more clinically relevant set-ups. New tests concern the benchmarking of Geant4-DNA physics and chemistry components for regression testing purposes, dosimetry for brachytherapy with a 125 I source, dosimetry for external x-ray and electron FLASH radiotherapy, experimental microdosimetry for proton therapy, and in vivo PET for carbon and oxygen beams. Regression testing has been performed between Geant4 10.5 and 11.1. Finally, a simple Geant4 simulation has been developed and used to compare Geant4 EM physics constructors and physics lists in terms of execution times. Results In summary, our EM tests show that the parameters of the multiple scattering in the Geant4 EM constructor G4EmStandardPhysics_option3 in Geant4 11.1, while improving the modeling of the electron backscattering in high atomic number targets, are not adequate for dosimetry for clinical x-ray and electron beams. Therefore, these parameters have been reverted back to those of Geant4 10.5 in Geant4 11.2.1. The x-ray radiotherapy test shows significant differences in the modeling of the bremsstrahlung process, especially between G4EmPenelopePhysics and the other constructors under study (G4EmLivermorePhysics, G4EmStandardPhysics_option3, and G4EmStandardPhysics_option4). These differences will be studied in an in-depth investigation within our Group. Improvement in Geant4 11.1 has been observed for the modeling of the proton and carbon ion Bragg peak with energies of clinical interest, thanks to the adoption of ICRU90 to calculate the low energy proton stopping powers in water and of the Linhard–Sorensen ion model, available in Geant4 since version 11.0. Nuclear fragmentation tests of interest for carbon ion therapy show differences between Geant4 10.5 and 11.1 in terms of fragment yields. In particular, a higher production of boron fragments is observed with Geant4 11.1, leading to a better agreement with reference data for this fragment. Conclusions Based on the overall results of our tests, we recommend to use G4EmStandardPhysics_option4 as EM constructor and QGSP_BIC_HP with G4EmStandardPhysics_option4, for hadrontherapy applications. The Geant4-DNA physics lists report differences in modeling electron interactions in water, however, the tests have a pure regression testing purpose so no recommendation can be formulated.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

Thermal Analysis and Simulations of the Abort Kicker Magnets Systems for the Electron-Ion Collider

The Abort Kicker Magnets system used in the current Relativistic Heavy Ion Collider (RHIC) to steer the circulating beam into the dump will be subjected to higher heat loads in the Electron-Ion Collider (EIC). After analyzing the existing abort kicker magnets and running thermal simulations in ANSYS it was concluded that they may not be suitable for use as-is in the future EIC due to heat and impedance concerns in the magnets. Possible solutions include adding a round titanium-coated ceramic beam tube to help solve the heat and impedance concerns, but this will reduce the limiting aperture of the kicker magnet. Another possible solution to meet all the performance requirements for EIC would be to add titanium-coated ceramic plates with water-cooling and tapered transitions that significantly improve the impedance and lower the heat in the magnets with less reduction to the aperture.

43 PARTICLE ACCELERATORS↗