Search NASA⌕ Search

SEARCH · Search NASA

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 415 records · Page 23

First Estimation of Model Parameters for Neutrino-Induced Nucleon Knockout Using Simulation-Based Inference

To enable an accurate determination of oscillation parameters, accelerator-based neutrino experiments require detailed simulations of nuclear interaction physics in the GeV regime. While substantial effort from both theory and experiment is currently being invested to improve the fidelity of these simulations, their present deficiencies typically oblige experimental collaborations to resort to empirical tuning of simulation model parameters. As the precision requirements of the field continue to become more stringent, machine learning techniques may provide a powerful means of handling corresponding growth in the complexity of future neutrino interaction model tuning exercises. To study the suitability of simulation-based inference (SBI) for this physics application, in this paper we revisit a tuned configuration of the GENIE neutrino event generator that was originally developed by the MicroBooNE collaboration. Despite closely reproducing the adopted values of four physics parameters when confronted with the tuned cross-section predictions as input, we find that our trained SBI algorithm prefers modestly different values (within MicroBooNE's assigned uncertainties) and achieves slightly better goodness-of-fit when inference is run on the experimental data set originally used by MicroBooNE. We also find that our trained algorithm can create a fair approximation of an alternative neutrino scattering simulation, NuWro, that shares only a subset of its physics model parameters with GENIE.

Tame-Narvaez, Karla [Fermilab] (ORCID:000000022249↗

Distributions of Particles Accelerated by Strong Alfvénic Turbulence

This work presents a model for generating nonthermal power-law tails of particles’ energy probability density functions in turbulent collisionless plasmas, applicable to both nonrelativistic and relativistic scenarios. We propose that strong Alfvénic turbulence energizes plasma particles through curvature acceleration, particularly for particles with Larmor radii comparable to the scales of turbulence. When the energy density of the energized particles increases, the efficiency of the energy exchange process diminishes. As a result, the acceleration process saturates, leading to power-law distributions of particle momentum and energy. In the nonrelativistic case, the momentum probability density function scales as f(p)dp ∝ p −3 dp, while in the ultrarelativistic case, the energy probability density function scales as f(γ)dγ ∝ γ −3 dγ, where γ is the Lorentz factor. This model provides a unified framework for understanding particle acceleration in both energy regimes, complementing existing analytical approaches. The predicted scalings are consistent with available observations of energetic ion distributions in the heliosphere and with the findings from numerical simulations of ultrarelativistic particle acceleration in magnetically dominated plasma turbulence.

Alfven waves↗

Simulation study of high-current 7 Li 3+ beam acceleration with laser ion source and RFQ linac using direct plasma injection

Accelerator-based neutron sources (ABNS) utilizing the 7 Li(p,n) 7 Be reaction in inverse kinematics offer strong forward-directed neutron emission, making them attractive for compact and clean neutron source with low background neutrons and unwanted radiation. However, practical use of such systems requires lithium ion beam currents that exceed the capability of conventional ion accelerators by two orders of magnitude. In this study, we numerically designed and evaluated a high-current acceleration system based on a direct plasma injection scheme (DPIS), combining a laser ion source (LIS) and a radio-frequency quadrupole (RFQ) linac. The ion extraction optics and RFQ entrance were optimized using IGUN, OPERA, and GPT, demonstrating that over 1.2 A-class beam current can be injected into the RFQ. The RFQ structure was designed under realistic engineering constraints including surface field strength (Kilpatrick factor) and RF power. GPT simulations showed that a 370 mA 7 Li 3+ beam can be successfully accelerated within a ±10 % energy spread. To validate downstream compatibility, we also designed a simplified medium-energy beam transport (MEBT) section and an interdigital H-type (IH) linac, confirming successful acceleration of a 320 mA beam to the total energy of 14 MeV. These results support the feasibility of the DPIS and RFQ approach as a promising solution for compact neutron sources and other high-intensity ion beam applications.

43 PARTICLE ACCELERATORS↗

Particle-acceleration by turbulent magnetohydrodynamic reconnection

Test particles in a two dimensional, turbulent MHD simulation are found to undergo significant acceleration. The magnetic field configuration is a periodic sheet pinch which undergoes reconnection. The test particles are trapped in the reconnection region for times of order an Alfven transit time in the large electric fields that characterize the turbulent reconnection process at the relatively large magnetic Reynolds number used in the simulation. The maximum speed attained by these particles is consistent with an analytic estimate which depends on the reconnection electric field, the Alfven speed, and the ratio of Larmor period to the Alfven transit time.

Matthaeus, W. H.↗

The production of ion conics by oblique double layers

Magnetized test ions are subjected to acceleration through a numerically simulated oblique double layer in order to determine whether they emerge with velocity vectors aligned with or oblique to the ambient magnetic field. A criterion for oblique alignment, depending on the double-layer parameters and on the external magnetization, is obtained. When it is applied to observed and theoretical auroral double layers, this criterion predicts that accelerated heavy ions will be substantially less magnetic field aligned than will accelerated hydrogen ions, thus suggesting auroral double layers as a source of high-energy ion conics. Test particle simulations are also used to investigate the perpendicular heating of ions at low altitudes by the electric fields associated with moving auroral arcs. The rapid motion of small-scale structures in the arcs is suggested as a source of low-energy conical ion distributions, and the slow drifts of the entire arc forms are inferred to heat ionospheric ions.

Borovsky, J. E.↗

MBX V1.2: Accelerating Data-Driven Many-Body Molecular Dynamics Simulations

The MBX software provides an advanced platform for molecular dynamics simulations, leveraging state-of-the-art MB-pol and MB-nrg data-driven many-body potential energy functions. Developed over the past decade, these potential energy functions integrate physics-based and machine-learned many-body terms trained on electronic structure data calculated at the "gold standard" coupled-cluster level of theory. Recent advancements in MBX have focused on optimizing its performance, resulting in the release of MBX v1.2. While the inherently many-body nature of MB-pol and MB-nrg ensures high accuracy, it poses computational challenges. MBX v1.2 addresses these challenges with significant performance improvements, including enhanced parallelism that fully harnesses the power of modern multicore CPUs. In conclusion, these advancements enable simulations on nanosecond time scales for condensed-phase systems, significantly expanding the scope of high-accuracy, predictive simulations of complex molecular systems powered by data-driven many-body potential energy functions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Zero-RK Acceleration of Low-Life-cycle Carbon Fuel (LLFC)Simulations

Designing modern combustion systems now relies on computer models that predict how changes in design will affect performance. These models have replaced older methods that relied on the designer’s intuition or costly and time-consuming physical testing. By using improved models, design cycles can be shortened and cleaner and more efficient combustion devices can be created. This project aims to improve computer simulations of low-life-cycle carbon fuels (LLCFs) with the goal of making these simulations faster and more accurate for predicting combustion in vehicles.

02 PETROLEUM↗

Navier-Stokes simulation of the supersonic combustion flowfield in a ram accelerator

A computational study of the ram accelerator, a ramjet-in-tube device for accelerating projectiles to ultrahigh velocities, is presented. The analysis is performed using a fully implicit TVD scheme that efficiently solves the Reynolds-averaged Navier-Stokes equations and the species continuity equations associated with a finite rate combustion model. Previous analyses of this concept were based on inviscid assumptions. The present results indicate that viscous effects are of primary importance; in all the cases studied, shock-induced combustion always started in the boundary layer. The effects of Mach number, mixture composition, pressure, and turbulence are investigated for various configurations. Two types of combustion processes, one stable and the other unstable, were observed depending on the inflow conditions. In the unstable case, a detonation wave is formed, which propagates upstream and unstarts the ram accelerator. In the stable case, a solution that converges to steady-state is obtained, in which the combustion wave remains stationary with respect to the ram accelerator projectile. The possibility of stabilizing the detonation wave by means of a backward facing step is also investigated. In addition to these studies, two numerical techniques were tested. These two techniques are vector extrapolation to accelerate convergence, and a diagonal formulation that eliminates the expense of inverting large block matrices that arise in chemically reacting flows.

Yungster, Shaye↗

Navier-Stokes simulation of the supersonic combustion flowfield in a ram accelerator

A computational study of the ram accelerator, a ramjet-in-tube device for accelerating projectiles to ultrahigh velocities, is presented. The analysis is performed using a fully implicit TVD scheme that efficiently solves the Reynolds-averaged Navier-Stokes equations and the species continuity equations associated with a finite rate combustion model. The present results indicate that viscous effects are of primary importance in all the cases studied, shock-induced combustion always started in the boundary layer. The effects of Mach number, mixture composition, pressure and tubulence are investigated for various configurations. Two types of combustion processes, one stable and the other unstable, were observed depending on the inflow conditions. The possibility of stabilizing the detonation wave by means of a backward facing step is also investigated. Two numerical techniques were tested: vector extrapolation, to accelerate convergence, and a diagonal formulation that eliminates the expense of inverting large block matrices which arise in chemically reacting flows.

Yungster, Shaye↗

A First-Principles Model Based on Saturation of the Electron Cyclotron Drift Instability for Electron Transport in Hydrodynamics Simulations of Hall Thruster Plasmas

We present a computational model based on the hypothesis that the higher-than-expected electron mobility across magnetic field lines observed in Hall thrusters is due to the growth and eventual saturation of the electron cyclotron drift instability. The key improvement of this model with respect to previous work is that it decouples the saturation of the instability from the correlation that exists between mobility and magnitude of the wave perturbations in linear theory, thereby enabling first-principles simulations. The model has been incorporated in the 2-D (r-z) multi-fluids code Hall2De and simulations of the H6 thruster at 300 V, 20 A are numerically stable and achieve a steady-state solution at a computational cost that is not significantly higher than previous Hall2De simulations. The computed location of the acceleration region is found to be within 10% of the length of the acceleration channel compared to that inferred by experiments. The simulations also capture well the plasma gradients along the channel centerline of this thruster. To further establish the validity of this model, we plan to conduct simulations of other thrusters and/or operating conditions.

Lopez Ortega, Alejandro↗

Implementation of Acceleration Source Term in FUN3D

Implementation of the acceleration source term in FUN3D is described. The method is evaluated using analytical and idealized test cases in two and three dimensions. It is also compared with results obtained from the moving mesh and six-degree-of-freedom simulations. The method shows promise in simulating the flight of accelerating vehicles as they make their ascent through varying reference conditions.

Nash’at N Ahmad↗

Hemodynamic responses to seated and supine lower body negative pressure - Comparison with +Gz acceleration

The hemodynamic responses to LBNP in seated subjects and in subjects in supine body positions were compared and were correlated with hemodynamic changes which occurred during a simulated (by centrifugation) Shuttle reentry acceleration with a slow onset rate of 0.002 G/s and during gradual onset exposures to +3 Gz and +4 Gz. Results demonstrate that seated LBNP at a level of -40 mm Hg can serve as a static simulator for changes in the heart rate and in mean blood pressure induced by gradual onset acceleration stress occurring during Shuttle reentry. The findings also provide a rationale for using LBNP during weightlessness as a means of imposing G-loading on the circulation prior to reentry.

Polese, Alvese↗

Software Tool for Computing Maximum Von Mises Stress

The maximum Van Mises stress and stress direction are of interest far analyzing launch accelerations such as with the Mass Acceleration Curves developed by JPL. Maximum launch stresses can be combined with appropriate load cases at consistent locations with resulting stress tensors. Maximum Van Mises stress is also of interest for understanding maximum operational loading such as traverse events. - For example, planetary traversing simulations may prescribe bounding acceleration values during traverse for a rover such as Mars Science Lab (MSL) in (X,Y,Z) of the rover. - Such accelerations can be really in any directions for many parts such as a mast or head mounted components which can be in numerous configurations and orientations when traversing a planet surface.

software↗

A Unified Computational Model for Solar and Stellar Flares

We present a unified computational framework that can be used to describe impulsive flares on the Sun and on dMe stars. The models assume that the flare impulsive phase is caused by a beam of charged particles that is accelerated in the corona and propagates downward depositing energy and momentum along the way. This rapidly heats the lower stellar atmosphere causing it to explosively expand and dramatically brighten. Our models consist of flux tubes that extend from the sub-photosphere into the corona. We simulate how flare-accelerated charged particles propagate down one-dimensional flux tubes and heat the stellar atmosphere using the Fokker-Planck kinetic theory. Detailed radiative transfer is included so that model predictions can be directly compared with observations. The flux of flare-accelerated particles drives return currents which additionally heat the stellar atmosphere. These effects are also included in our models. We examine the impact of the flare-accelerated particle beams on model solar and dMe stellar atmospheres and perform parameter studies varying the injected particle energy spectra. We find the atmospheric response is strongly dependent on the accelerated particle cutoff energy and spectral index.

radiative transfer↗

A GPU ‐Accelerated 3D Unstructured Mesh Based Particle Tracking Code for Multi‐Species Impurity Transport Simulation in Fusion Tokamaks

ABSTRACT This paper presents the multi‐species global impurity transport capability developed in a GPU‐accelerated fully 3D unstructured mesh‐based code, GITRm, to simultaneously track multiple impurity species and handle interactions of these impurities with mixed‐material surfaces. Different computational approaches to model particle‐surface interaction or surface response have been developed and compared. Sheath electric field is taken into account by employing a fast distance‐to‐boundary calculation, which is carried out in parallel on distributed or partitioned meshes on multiple GPUs without the need for any inter‐process communication during the simulation. Several example cases, including two for the DIII‐D tokamak, that is, one with the SAS‐V divertor and the other with the collector probes, are used to demonstrate the utility of the current multi‐species capability. For the DIII‐D probe case, the capability of GITRm to resolve the spatial distribution of particles in localized regions, such as diagnostic probes, within non‐axisymmetric tokamak geometries is demonstrated. These simulations involve up to 320 million particles and utilize up to 48 GPUs.

Nath, Dhyanjyoti D. [Scientific Computation Resear↗