Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Accelerator 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 649 records · Page 36

SU-E-T-66: Characterization of Radiation Dose Associated with Dark Currents During Beam Hold for Respiratory-Gated Electron Therapy

Purpose: The main objective of this study was to estimate the radiation dose contributed by dark currents associated with the respiratory-gated electron therapy during beam hold. The secondary aim was to determine clinical benefits of using respiratory-gated electron therapy for left-sided breast cancer patients with positive internal mammary nodes (IMN). Methods: Measurements of the dark current-induced dose in all electron modes were performed on multiple Siemens and Varian linear accelerators by manually simulating beam-hold during respiratory gating. Dose was quantified at the machine isocenter by comparing the collected charge to the known output for all energies ranging from 6 to 18 MeV for a 10cm × 10cm field at 100 SSD with appropriate solid-water buildup. Using the Eclipse treatment planning system, we compared the additional dose associated with dark current using gated electron fields to the dose uncertainties associated with matching gated photon fields and ungated electron fields. Dose uncertainties were seen as hot and cold spots along the match line of the fields. Results: The magnitude of the dose associated with dark current is highly correlated to the energy of the beam and the amount of time the beam is on hold. For lower energies (6–12 MeV), there was minimal dark current dose (0.1–1.3 cGy/min). Higher energies (15–18 MeV) showed measurable amount of doses. The dark current associated with the electron beam-hold varied between linear accelerator vendors and depended on dark current suppression and the age of the linear accelerator. Conclusion: For energies up to 12 MeV, the dose associated with the dark current for respiratorygated electron therapy was shown to be negligible, and therefore should be considered an option for treating IMN positive left-sided breast cancer patients. However, at higher energies the benefit of respiratory gating may be outweighed by dose due to the dark current.

60 APPLIED LIFE SCIENCES↗

Enhanced ion acceleration in transition from opaque to transparent plasmas

Using particle-in-cell simulations, we investigate ion acceleration in the interaction of high intensity lasers with plasmas which transition from opaque to transparent during the interaction process. We show that the highest ion energies are achieved when the laser traverses the target around the peak intensity and re-heats the electron population responsible for the plasma expansion, enhancing the corresponding sheath electric field. This process can lead to an increase of up to 2x in ion energy when compared with the standard Target Normal Sheath Acceleration in opaque targets under the same laser conditions. A theoretical model is developed to predict the optimal target areal density as a function of laser intensity and pulse duration. A systematic parametric scan for a wide range of target densities and thicknesses is performed in 1D, 2D and 3D and shown consistent with the theory and with recent experimental results. Thus, these results open the way for a better optimization of the ion energy in future laser–solid experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MuSIKAL: Multiphysics Simulations and Knowledge Discovery through AI/ML Technologies

Under the MuSiKAL project, we developed a framework for a coastal digital twin (DT) platform capable of integrating diverse data resources, configuring multiscale model simulations, performing SciML‐accelerated predictions, with applications primarily driven by storm surge and heavily rainfall events impacting the Gulf Coast of the U.S.

54 ENVIRONMENTAL SCIENCES↗

Pulsed-focusing recirculating linacs for muon acceleration

Since the muon has a short lifetime, fast acceleration is essential for high-energy applications such as muon colliders, Higgs factories, or neutrino factories. The best one can do is to make a linear accelerator with the highest possible accelerating gradient to make the accelerating time as short as possible. However, the cost of such a single linear accelerator is prohibitively large due to expensive power sources, cavities, tunnels, and related infrastructure. As was demonstrated in the Thomas Jefferson Accelerator Facility (Jefferson Lab) Continuous Electron Beam Accelerator Facility (CEBAF), an elegant solution to reduce cost is to use magnetic return arcs to recirculate the beam through the accelerating RF cavities many times, where they gain energy on each pass. In such a Recirculating Linear Accelerator (RLA), the magnetic focusing strength diminishes as the beam energy increases in a conventional linac that has constant strength quadrupoles. After some number of passes the focusing strength is insufficient to keep the beam from going unstable and being lost. In this project, the use of fast pulsed quadrupoles in the linac sections was considered for stronger focusing as a function of time to allow more successive passes of a muon beam in a recirculating linear accelerator. In one simulation, it was shown that the number of passes could be increased from 8 to 12 using pulsed magnet designs that have been developed and tested. This could reduce the cost of linac sections of a muon RLA by 8/12, where more improvement is still possible. The expense of a greater number of passes and corresponding number of return arcs was also addressed in this project by exploring the use of ramped or FFAG-style magnets in the return arcs. A better solution, invented in this project, is to use combined-function dipole-quadrupole magnets to simultaneously transport two beams of different energies through one magnet string to reduce costs of return arcs by almost a factor of two. A patent application was filed for this invention and a detailed report published in Physical Review Special Topics. A scaled model using an electron beam was developed and proposed to test the concept of a dog bone RLA with combined-function return arcs. The efforts supported by this grant were reported in a series of contributions to particle accelerator conferences that are reproduced in the appendices and summarized in the body of this report.

43 PARTICLE ACCELERATORS↗

Influence of Longer Range Transfer of Vapor Interface Modified Caging Constraints on the Spatially Heterogeneous Dynamics of Glass-Forming Liquids

Based on the elastically collective nonlinear Langevin equation (ECNLE) theory of bulk relaxation in glass-forming liquids, and our recent ideas of how interface-nucleated modification of caging constraints are spatially transferred into the interior of a thick film, we present a force-based theory for dynamical gradients in thick films with one vapor interface that includes collective elasticity effects. Quantitative applications to the foundational hard sphere fluid and polymer melts of diverse fragilities are presented. In this work, we predict a roughly exponential spatial variation of the total activation barrier which is nonperturbatively modified to similar to 10 particle diameters into the film. This leads, to leading order, to the prediction of a reduced alpha relaxation time gradient of a double-exponential form characterized by a nearly constant spatial decay length, in qualitative accord with simulations. The relaxation acceleration at the surface grows exponentially with increasing packing fraction or decreasing temperature. The rate of increase with cooling of the alpha time strongly weakens upon approaching the vapor interface, with the top two layers remaining liquid-like down to similar to ~80-85% of the bulk glass transition temperature. These spatially heterogeneous changes of the temperature variation of the alpha time result in a large and long-range gradient of the local glass transition temperature. An average interfacial layer thickness relevant to ensemble-averaged dielectric and other experiments is also computed. It is predicted to be rather large at the bulk glass transition temperature, decreasing roughly linearly with heating. Remarkably, to leading order the ratio of this layer-averaged interfacial relaxation time to its bulk analogue is invariant to chemistry, volume fraction, and temperature and of modest absolute magnitude. Spatially inhomogeneous power law decoupling of the alpha relaxation time from its bulk value is predicted, with an effective exponent that decays to zero with distance from the free surface in a nearly exponential manner, trends which are in qualitative accord with recent simulations. Finally, this behavior and the double-exponential alpha time gradient are related and can be viewed as consequences of an effective quasi-universal factorization of the total barrier in films into the product of its bulk-temperature-dependent value times a function solely of location in the film.

36 MATERIALS SCIENCE↗

Multigroup Radiation-Hydrodynamics with a High-Order, Low-Order Method

Recent efforts at Los Alamos National Laboratory to develop a moment-based, scale-bridging [or high-order (HO)–low-order (LO)] algorithm for solving large varieties of the transport (kinetic) systems have shown promising results. A part of our ongoing effort is incorporating this methodology into the framework of the Eulerian Applications Project to achieve algorithmic acceleration of radiationhydrodynamics simulations in production software. By starting from the thermal radiative transfer equations with a simple material-motion correction, we derive a discretely consistent energy balance equation (LO equation). We demonstrate that the corresponding LO system for the Monte Carlo HO solver is closely related to the original LO system without material-motion corrections. We test the implementation on a radiative shock problem and show consistency between the energy densities and temperatures in the HO and LO solutions as well as agreement with the semianalytic solution. We also test the approach on a more challenging two-dimensional problem and demonstrate accuracy enhancements and algorithmic speedups. This paper extends a recent conference paper by including multigroup effects.

42 ENGINEERING↗

The Roles of Fluid Compression and Shear in Electron Energization during Magnetic Reconnection

Particle acceleration in space and astrophysical reconnection sites is an important unsolved problem in studies of magnetic reconnection. Earlier kinetic simulations have identified several acceleration mechanisms that are associated with particle drift motions. Here, we show that, for sufficiently large systems, the energization processes due to particle drift motions can be described as fluid compression and shear, and that the shear energization is proportional to the pressure anisotropy of energetic particles. By analyzing results from fully kinetic simulations, we show that the compression energization dominates the acceleration of high-energy particles in reconnection with a weak guide field, and the compression and shear effects are comparable when the guide field is 50% of the reconnecting component. Spatial distributions of those energization effects reveal that reconnection exhausts, contracting islands, and island-merging regions are the three most important regions for compression and shear acceleration. This study connects particle energization by particle guiding-center drift motions with that due to background fluid motions, as in the energetic particle transport theory. It provides foundations for building particle transport models for large-scale reconnection acceleration such as those in solar flares.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Start-to-End Simulation of the Drive-Beam Longitudinal Dynamics for Beam-Driven Wakefield Acceleration

Collinear beam-driven wakefield acceleration (WFA) relies on shaped driver beam to provide higher accelerating gradient at a smaller cost and physical footprint. This acceleration scheme is currently envisioned to accelerate electron beams capable of driving free-electron laser *. Start-to-end simulation of drive-bunch beam dynamics is crucial for the evaluation of the design of accelerators built upon WFA. We report the start-to-end longitudinal beam dynamics simulations of an accelerator beamline capable of producing high charge drive beam. The generated wakefield when it passes through a corrugated waveguide results in a transformer ratio of 5. This paper especially discusses the challenges and criteria associated with the generation of temporally-shaped driver beam, including the beam formation in the photoinjector, and the influence of energy chirp control on beam transport stability.

43 PARTICLE ACCELERATORS↗

Particle Acceleration in Two Converging Shocks

Observations by spacecraft such as ACE , STEREO , and others show that there are proton spectral “breaks” with energy E {sub br} at 1–10 MeV in some large CME-driven shocks. Generally, a single shock with the diffusive acceleration mechanism would not predict the “broken” energy spectrum. The present paper focuses on two converging shocks to identify this energy spectral feature. In this case, the converging shocks comprise one forward CME-driven shock on 2006 December 13 and another backward Earth bow shock. We simulate the detailed particle acceleration processes in the region of the converging shocks using the Monte Carlo method. As a result, we not only obtain an extended energy spectrum with an energy “tail” up to a few 10 MeV higher than that in previous single shock model, but also we find an energy spectral “break” occurring on ∼5.5 MeV. The predicted energy spectral shape is consistent with observations from multiple spacecraft. The spectral “break,” then, in this case is caused by the interaction between the CME shock and Earth’s bow shock, and otherwise would not be present if Earth were not in the path of the CME.

79 ASTRONOMY AND ASTROPHYSICS↗

Orchestrating Spontaneous Emission With Metasurfaces: Recent Advances in Engineering Thermal, Luminescent, and Quantum Emissions

Metasurfaces have emerged as powerful tools for controlling spontaneous emission, offering unprecedented control over light-matter interactions at sub-wavelength scales. While metasurfaces are traditionally utilized for shaping coherent electromagnetic waves, they have recently extended their capabilities to control incoherent or spontaneous emission. This examines review how metasurfaces can enhance and precisely control properties of thermal, luminescent, and quantum emission. In thermal emission, metasurfaces enable control over spatial, temporal, and spin coherence, offering new possibilities for applications such as energy harvesting, radiative cooling and heat assisted ranging and detection. For luminescent emission, metasurfaces significantly improve emission rates, quantum efficiency, and directionality, driving innovations in lighting and display technologies. For controlling quantized spontaneous emission, metasurfaces are instrumental in enhancing single-photon sources and enabling novel functionalities in quantum states through photon-pair generation, which is vital for quantum communication, meteorology, and computing. Here, despite these advancements several challenges to increase the operational bandwidths, accelerate and develop simulation strategies, and fabrication complexities persist. Emerging trends are also discussed, such as dynamic metasurfaces and their integration with nanophotonic platforms, which could further expand the capabilities of light-emitting metasurfaces.

Luminescence↗

Thermodynamic Modeling of the Al-Ce-Cu-Mg-Si System and Its Application to Aluminum-Cerium Alloy Design

Recently discovered AlCe alloys have shown promise in a number of applications, but the propensity of Ce to react with Al and other alloying elements can complicate the phase equilibria and design approach. To solve this, the CALPHAD method is used to explore an alloy within the quinary Al-Ce-Cu-Mg-Si system by developing a thermodynamic database with self-consistent parameters. The database includes a description of all 10 binary systems and 8 ternary systems consisting of: (i) 6 Al-containing ternaries (Al-Ce-Cu, Al-Ce-Mg, Al-Ce-Si, Al-Cu-Mg, Al-Cu-Si and Al-Mg-Si); and (ii) 2 additional ternaries that include Mg and Si (i.e., Ce-Mg-Si and Cu-Mg-Si). The thermodynamic description for the Al-Ce-Mg and Al-Mg-Si systems were reassessed to ensure consistency with the binary systems and the Ce-Mg-Si system is presented for the first time and compared to theoretical data from DFT (Density Functional Theory). In addition to the ternary interactions, the quaternary compound Al3Cu2Mg9Si7 and solid solution extending from the ternary Al2CuMg phase (Al,Si)2CuMg are incorporated. The CALPHAD method is employed and leveraged through the use of a Materials Design Simulator (MDS) to accelerate the design of novel aluminum-cerium-based alloys. The combination of a CALPHAD-based framework with experimental efforts and industrial insight permits the development of three new Al-Ce alloys: Al-3.5Ce-0.4Mg-7Si (Ce-modified A356), Al-5Ce-1Cu-0.5Mg-10Si and Al-19Ce-0.9Mg-1.1Si.

36 MATERIALS SCIENCE↗

Next-generation anion exchange membrane water electrolyzers operating for commercially relevant lifetimes

Alkaline anion exchange membrane (AEM) water electrolysis has gained increasing attention due to its potential to achieve low-cost, high performance hydrogen production. However, most existing membranes are not durable in industrial settings. In this work, we demonstrate good performance relative to industrial parameters using Sustainion® anion exchange membranes. Long-duration tests showed stable performance of 1 A/cm 2 at 1.85 V with a degradation rate of less than 1 μV/h over 10,000 h. The projected lifetime is thus over 20 years. Daily on/off cycling performance over the course of 30 years was simulated experimentally through accelerated voltage shock tests, resulting in a performance loss of only 0.15 μV/cycle over 11,000 cycles. As shown through impact and crossover testing, an improvement in performance is achieved by the addition of zirconia to the polymer matrix and mechanically reinforcing the membrane.

08 HYDROGEN↗

Meso-scale modeling of irradiated concrete in test reactor

In this paper, we detail a numerical model accounting for the effects of neutron irradiation on concrete at the mesoscale. Irradiation experiments in test reactor (Elleuch et al.,1972), i.e., in accelerated conditions, are simulated. Concrete is considered as a two-phase material made of elastic inclusions (aggregate) subjected to thermal and irradiation-induced swelling and embedded in a cementitious matrix subjected to shrinkage and thermal expansion. The role of the hardened cement paste in the post-peak regime (brittle-ductile transition with decreasing loading rate), and creep effects are investigated. Radiation-induced volumetric expansion (RIVE) of the aggregate cause the development and propagation of damage around the aggregate which further develops in bridging cracks across the hardened cement paste between the individual aggregate particles. The development of damage is aggravated when shrinkage occurs simultaneously with RIVE during the irradiation experiment. The post-irradiation expansion derived from the simulation is well correlated with the experimental data and, the obtained damage levels are fully consistent with previous estimations based on a micromechanical interpretation of the experimental post-irradiation elastic properties (Le Pape et al.,2015). In conclusion, the proposed modeling opens new perspectives for the interpretation of test reactor experiments in regards to the actual operation of light water reactors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A pathway to synthesizing single-crystal Fe and FeCr films

Nuclear reactor environments provide a unique scientific and engineering challenge wherein materials must tolerate prolonged exposure to concurrent irradiation, elevated temperatures, and corrosive media. However, uncontrolled variability in material composition and structure often prohibits truly single-variable experiments that can reveal basic aspects of environmental damage. Here, Magnetron sputtering is used to provide a more controlled model system for these fundamental studies, yielding reproducible single-crystal Fe and FeCr thin films containing 8 and 18 at.% Cr. Electron microscopy is used to determine the systematic correlations between growth conditions and the resulting film microstructure and surface morphology. It is found that the substrate temperature and applied radio frequency (RF) bias can be tuned to obtain consistent homogeneous and single crystal films with a minimal amount of Ar impurities from the RF bias process. Epitaxial, single-crystal Fe films are obtained on MgO substrates at 500 °C with 10 Watt (W) RF bias deposition. However, when Cr is alloyed with Fe, higher substrate temperatures (600 °C) and applied RF biases (15 W) are required to achieve a similar epitaxial single-crystal FeCr film. Accelerated molecular dynamics simulations reveal that Cr impedes surface transport, explaining the need for higher temperature and bias during the growth of the Cr-bearing films.

36 MATERIALS SCIENCE↗

Dynamics of rapidly depressurized multiphase shock tubes

Rapid depressurization is a fluid phenomenon that occurs in many industrial and natural applications. Its behaviour is often complicated by the formation, propagation and interaction of waves. In this work, we perform computer simulations of the rapid depressurization of a gas–solid mixture in a shock tube. Our problem set-up mimics previously performed experiments, which have been historically used as a laboratory surrogate for volcanic eruptions. The simulations are carried out with a discontinuous Galerkin compressible fluid solver with four-way coupled Lagrangian particle tracking capabilities. The results give an unprecedented look into the complex multiphase physics at work in this problem. Different regimes have been characterized in a regime map that highlights the key observations. While the mean flow behaviour is in good agreement with experiments, the simulations show unexpected accelerations of the particle front as it expands. Additionally, a new lifting mechanism for gas bubble (void) growth inside the gas–solid mixture is detailed.

Mechanics↗

CoarsenConf: Equivariant Coarsening with Aggregated Attention for Molecular Conformer Generation

Molecular conformer generation (MCG) is an important task in cheminformatics and drug discovery. The ability to efficiently generate low-energy 3D structures can avoid expensive quantum mechanical simulations, leading to accelerated virtual screenings and enhanced structural exploration. Several generative models have been developed for MCG, but many struggle to consistently produce high-quality conformers for meaningful downstream applications. To address these issues, we introduce CoarsenConf, which coarse-grains molecular graphs based on torsional angles and integrates them into an SE(3)-equivariant hierarchical variational autoencoder. Through equivariant coarse-graining, we aggregate the fine-grained atomic coordinates of subgraphs connected via rotatable bonds, creating a variable-length coarse-grained latent representation. Our model uses a novel aggregated attention mechanism to restore fine-grained coordinates from the coarse-grained latent representation, enabling efficient generation of accurate conformers. Furthermore, we evaluate the chemical and biochemical quality of our generated conformers on multiple downstream applications, including property prediction and large-scale oracle-based protein docking. Overall, CoarsenConf generates more accurate conformer ensembles compared to prior generative models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Semi-Empirical Shadow Molecular Dynamics: A PyTorch Implementation

Here, extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) in its most recent shadow potential energy version has been implemented in the semiempirical PyTorch-based software PySeQM. The implementation includes finite electronic temperatures, canonical density matrix perturbation theory, and an adaptive Krylov subspace approximation for the integration of the electronic equations of motion within the XL-BOMB approach (KSA-XL-BOMD). The PyTorch implementation leverages the use of GPU and machine learning hardware accelerators for the simulations. The new XL-BOMD formulation allows studying more challenging chemical systems with charge instabilities and low electronic energy gaps. The current public release of PySeQM continues our development of modular architecture for large-scale simulations employing semi-empirical quantum-mechanical treatment. Applied to molecular dynamics, simulation of 840 carbon atoms, one integration time step executes in 4 s on a single Nvidia RTX A6000 GPU.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗