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 1,657 records · Page 92

Fluid simulation of the bias effect in inductive/capacitive discharges

Computer simulations are performed for an argon inductively coupled plasma (ICP) with a capacitive radio-frequency bias power, to investigate the bias effect on the discharge mode transition and on the plasma characteristics at various ICP currents, bias voltages, and bias frequencies. When the bias frequency is fixed at 13.56 MHz and the ICP current is low, e.g., 6 A, the spatiotemporal averaged plasma density increases monotonically with bias voltage, and the bias effect is already prominent at a bias voltage of 90 V. The maximum of the ionization rate moves toward the bottom electrode, which indicates clearly the discharge mode transition in inductive/capacitive discharges. At higher ICP currents, i.e., 11 and 13 A, the plasma density decreases first and then increases with bias voltage, due to the competing mechanisms between the ion acceleration power dissipation and the capacitive power deposition. At 11 A, the bias effect is still important, but it is noticeable only at higher bias voltages. At 13 A, the ionization rate is characterized by a maximum at the reactor center near the dielectric window at all selected bias voltages, which indicates that the ICP power, instead of the bias power, plays a dominant role under this condition, and no mode transition is observed. Indeed, the ratio of the bias power to the total power is lower than 0.4 over a wide range of bias voltages, i.e., 0–300 V. Besides the effect of ICP current, also the effect of various bias frequencies is investigated. It is found that the modulation of the bias power to the spatiotemporal distributions of the ionization rate at 2 MHz is strikingly different from the behavior observed at higher bias frequencies. Furthermore, the minimum of the plasma density appears at different bias voltages, i.e., 120 V at 2 MHz and 90 V at 27.12 MHz.

42 ENGINEERING↗

Time-dependent modeling of gamma-ray flares in blazar PKS1510–089

Here we present a new approach for constraining luminous blazars, incorporating fully time-dependent and self-consistent modeling of bright γ-ray flares of PKS 1510–089, resolved with Fermi-LAT, in the framework of the internal shock scenario. The results of our modeling imply the location of the γ-ray flaring zone to be outside of the broad-line region, namely around ≃0.3 pc from the core for a free-expanding jet with the opening angle Γ θ{sub jet}≃1 (where Γ is the jet bulk Lorentz factor), up to ≃3 pc for a collimated outflow with Γ θ{sub jet}≃0.1. Moreover, under the Γ θ{sub jet}≃1 condition, our modeling indicates the maximum efficiency of the jet production during the flares, with the total jet energy flux strongly dominated by protons and exceeding the available accretion power in the source. This is in contrast to the quiescence states of the blazar, characterized by lower jet kinetic power and an approximate energy equipartition between different plasma constituents. We demostrate how strictly simultaneous observations of flaring PKS 1510–089 at optical, X-ray, and GeV photon energies, on hourly timescales, augmented by extensive simulations as presented in this paper, may help to impose further precise constraints on the magnetization and opening angle of the emitting region. In addition, our detailed modeling implies that a non-uniformity of the Doppler factor across the jet, caused by the radial expansion of the outflow, may lead to a pronounced time distortion in the observed γ-ray light curves, resulting, in particular, in asymmetric flux profiles with substantially extended decay phases.

79 ASTRONOMY AND ASTROPHYSICS↗

Monte Carlo computer simulation of sedimentation of charged hard spherocylinders

In this article we present a NVT Monte Carlo computer simulation study of sedimentation of an electroneutral mixture of oppositely charged hard spherocylinders (CHSC) with aspect ratio L/σ = 5, where L and σ are the length and diameter of the cylinder and hemispherical caps, respectively, for each particle. This system is an extension of the restricted primitive model for spherical particles, where L/σ = 0, and it is assumed that the ions are immersed in an structureless solvent, i.e., a continuum with dielectric constant D. The system consisted of N = 2000 particles and the Wolf method was implemented to handle the coulombic interactions of the inhomogeneous system. Results are presented for different values of the strength ratio between the gravitational and electrostatic interactions, Γ = (mgσ)/(e{sup 2}/Dσ), where m is the mass per particle, e is the electron's charge and g is the gravitational acceleration value. A semi-infinite simulation cell was used with dimensions L{sub x} ≈ L{sub y} and L{sub z} = 5L{sub x}, where L{sub x}, L{sub y}, and L{sub z} are the box dimensions in Cartesian coordinates, and the gravitational force acts along the z-direction. Sedimentation effects were studied by looking at every layer formed by the CHSC along the gravitational field. By increasing Γ, particles tend to get more packed at each layer and to arrange in local domains with an orientational ordering along two perpendicular axis, a feature not observed in the uncharged system with the same hard-body geometry. This type of arrangement, known as tetratic phase, has been observed in two-dimensional systems of hard-rectangles and rounded hard-squares. In this way, the coupling of gravitational and electric interactions in the CHSC system induces the arrangement of particles in layers, with the formation of quasi-two dimensional tetratic phases near the surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Beam dynamic simulation and optimization of the CLIC positron source and the capture linac

The CLIC Positron Source is based on the hybrid target composed of a crystal and an amorphous target. Simulations have been performed from the exit of the amorphous target to the end of pre-injector linac which captures and accelerates the positrons to an energy of 200 MeV. Simulations are performed by the particle tracking code PARMELA. The magnetic field of the AMD is represented in PARMELA by simple coils. Two modes are applied in this study. The first one is accelerating mode based on acceleration after the AMD. The second one is decelerating mode based on deceleration in the first accelerating structure. It is shown that the decelerating mode gives a higher yield for the e{sup +} beam in the end of the Pre-Injector Linac.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

EFFECT OF CORONAL TEMPERATURE ON THE SCALE OF SOLAR CHROMOSPHERIC JETS

We investigate the effect of coronal temperature on the formation process of solar chromospheric jets using two-dimensional magnetohydrodynamic simulations of the region from the upper convection zone to the lower corona. We develop a new radiative magnetohydrodynamic code for the dynamic modeling of the solar atmosphere, employing an LTE equation of state, optically thick radiative loss in the photosphere, optically thin radiative loss in the chromosphere and the corona, and thermal conduction along the magnetic field lines. Many chromospheric jets are produced in the simulations by shock waves passing through the transition region. We find that these jets are projected farther outward when the coronal temperature is lower (similar to that in coronal holes) and shorter when the coronal temperature is higher (similar to that in active regions). When the coronal temperature is high, the deceleration of the chromospheric jets is consistent with the model in which deceleration is determined by the periodic chromospheric shock waves. However, when the coronal temperature is low, the gravitational deceleration becomes more important and the chromospheric jets approach ballistic motion.

79 ASTRONOMY AND ASTROPHYSICS↗

An Analytical Comparison of the Fidelity of "Large Motion" Versus "Small Motion" Flight Simulators in a Rotorcraft Side-Step Task

This paper presents an analytical and experimental methodology for studying flight simulator fidelity. The task was a rotorcraft bob-up/down maneuver in which vertical acceleration constituted the motion cue. The task considered here is aside-step maneuver that differs from the bob-up one important way: both roll and lateral acceleration cues are available to the pilot. It has been communicated to the author that in some Verticle Motion Simulator (VMS) studies, the lateral acceleration cue has been found to be the most important. It is of some interest to hypothesize how this motion cue associated with "outer-loop" lateral translation fits into the modeling procedure where only "inner-loop " motion cues were considered. This Note is an attempt at formulating such an hypothesis and analytically comparing a large-motion simulator, e.g., the VMS, with a small-motion simulator, e.g., a hexapod.

Hess, Ronald A.↗

Superfluid condensate fraction and pairing wave function of the unitary Fermi gas

The unitary Fermi gas is a many-body system of two-component fermions with zero-range interactions tuned to infinite scattering length. Despite much activity and interest in unitary Fermi gas and its universal properties, there have been great difficulties in performing accurate calculations of the superfluid condensate fraction and pairing wave function. In this paper, we present auxiliary-field lattice Monte Carlo simulations using a lattice interaction which accelerates the approach to the continuum limit, thereby allowing for robust calculations of these difficult observables. As a benchmark test, we compute the ground-state energy of 33 spin-up and 33 spin-down particles. As a fraction of the free Fermi gas energy $E_{\text{FG}}$, we find $E_0/E_{\text{FG}}$ = 0.369(2), 0.372(2), using two different definitions of the finite-system energy ratio, in agreement with the latest theoretical and experimental results. We then determine the condensate fraction by measuring off-diagonal long-range order in the two-body density matrix. We find that the fraction of condensed pairs is α = 0.43(2). Further, we also extract the pairing wave function and find the pair correlation length to be $ζ_pk_F$ = 1.8(3)ℏ, where $k_F$ is the Fermi momentum. Provided that the simulations can be performed without severe sign oscillations, the methods we present here can be applied to superfluid neutron matter as well as more exotic $\textit{P}$-wave and $\textit{D}$-wave superfluids.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Particle acceleration and plasma energization in substorms: MHD and test particle studies

The author organizes his slide presentation under the following topics: background, MHD simulation, orbit integration, typical orbits, spatial and temporal features, acceleration mechanisms, source locations, and source energies. Field-­aligned energetic particle fluxes are shown for 45-keV electrons and 80-keV protons. It is concluded that the onset from local thin current sheet is electron tearing. Acceleration is mainly from field collapse, governed by E y = -v x XB z : importance of localization; betatron acceleration (similar if nonadiabatic); 1st order Fermi, type B (or A; current sheet acceleration). There are two source regions (of comparable importance in magnetotail): - flanks, inner tail - drift entry - early, higher energy - outer plasma sheet - reconnection entry - later, lower energy. Both thermal and suprathermal sources are important, with limited energy range for acceleration

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The impact of positron range on PET resolution, evaluated with phantoms and PHITS Monte Carlo simulations for conventional and non-conventional radionuclides

The increasing interest and availability of non-standard positron-emitting radionuclides has heightened the relevance of radionuclide choice in the development and optimization of new positron emission tomography (PET) imaging procedures, both in preclinical research and clinical practice. Differences in achievable resolution arising from positron range can largely influence application suitability of each radionuclide, especially in small-ring preclinical PET where system blurring factors due to annihilation photon acollinearity and detector geometry are less significant. Here, some resolution degradation can be mitigated with appropriate range corrections implemented during image reconstruction, the quality of which is contingent on an accurate characterization of positron range.

43 PARTICLE ACCELERATORS↗

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory↗

Stable long range proton acceleration driven by intense laser pulse with underdense plasmas

Proton acceleration is investigated by 2.5-dimensional particle-in-cell simulations in an interaction of an ultra intense laser with a near-critical-density plasma. It was found that multi acceleration mechanisms contribute together to a 1.67 GeV collimated proton beam generation. The W-BOA (breakout afterburner based on electrons accelerated by a wakefield) acceleration mechanism plays an important role for the proton energy enhancement in the area far from the target. The stable and continuous acceleration maintains for a long distance and period at least several pico-seconds. Furthermore, the energy scalings are also discussed about the target density and the laser intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fingering convection induced by atomic diffusion in stars: 3D numerical computations and applications to stellar models

Iron-rich layers are known to form in the stellar subsurface through a combination of gravitational settling and radiative levitation. Their presence, nature, and detailed structure can affect the excitation process of various stellar pulsation modes and must therefore be modeled carefully in order to better interpret Kepler asteroseismic data. In this paper, we study the interplay between atomic diffusion and fingering convection in A-type stars, as well as its role in the establishment and evolution of iron accumulation layers. To do so, we use a combination of three-dimensional idealized numerical simulations of fingering convection (which neglect radiative transfer and complex opacity effects) and one-dimensional realistic stellar models. Using the three-dimensional simulations, we first validate the mixing prescription for fingering convection recently proposed by Brown et al. (within the scope of the aforementioned approximation) and identify what system parameters (total mass of iron, iron diffusivity, thermal diffusivity, etc.) play a role in the overall evolution of the layer. We then implement the Brown et al. prescription in the Toulouse-Geneva Evolution Code to study the evolution of the iron abundance profile beneath the stellar surface. We find, as first discussed by Théado et al., that when the concurrent settling of helium is ignored, this accumulation rapidly causes an inversion in the mean molecular weight profile, which then drives fingering convection. The latter mixes iron with the surrounding material very efficiently, and the resulting iron layer is very weak. However, taking helium settling into account partially stabilizes the iron profile against fingering convection, and a large iron overabundance can accumulate. The opacity also increases significantly as a result, and in some cases it ultimately triggers dynamical convection. The direct effects of radiative acceleration on the dynamics of fingering convection (especially in the nonlinear regime) remain to be added in the future to improve the quantitative predictions of the model.

79 ASTRONOMY AND ASTROPHYSICS↗

Coupling of chemical kinetics and thermodynamics for simulations of leaching of cement paste in ammonium nitrate solution

A novel numerical method is proposed in this paper to model the leaching of various types of cement pastes in the ammonium nitrate solution by coupling chemical kinetics with chemical thermodynamics. Arrhenius equation and Eyring equation are applied to predict the rate constant of fast dissolution of portlandite in the ammonium nitrate, and a modified rate law is developed by incorporating the saturation ratio of the reaction. The estimations of leached depths from the proposed method comply with the well-known square root-time relationship and they also agree well with the reported experimental results. Moreover, the evolutions of porosity and major hydrate compositions during leaching are obtained. It is found that a simple linear relation exists for the acceleration factor of leaching in 6 M of ammonium nitrate with respect to leaching in pure water regardless of the type of the cement clinker. - Highlights: • Chemical kinetics and thermodynamics are coupled to simulate the accelerated leaching. • The proposed method generates good estimations for leached depths and changes of porosity. • The proposed method can predict the detailed distributions of all the major cement hydrates. • A simple linear relation exists for predicting the widely scattered acceleration factors.

36 MATERIALS SCIENCE↗

RF Gap Transformation in PARMILA

The PARMILA (Phase And Radial Motion in Ion Linear Accelerators) code is used to design ion linacs and simulate the beam dynamics therein. The original code was developed in the 1960’s and later the Los Alamos Accelerator Code Group developed a PC version that has been used worldwide. Recently, in an effort to modernize the code and use it as the basis for a high-performance beam dynamics simulator with the LANSCE linac, the details of the code have come under additional scrutiny. PARMILA is a “z code” that uses various transformations to approximately represent the effect of the beam line element on the beam. In the case of the RF gap, it treats the particle dynamics inside the gap as a combination of drift + kick at the mid-plane + drift. Recently, we uncovered an undocumented keyword “off-exact” in PARMILA that enables the code to use the full modified-Bessel functions instead of low-order approximations in the RF gap transformation. This change improved the accuracy of the transformation for off-axis particles. However, a discrepancy between the analytical formula and the code was found for the transformation of the transverse coordinates. The purpose of this tech note is to document the derivation of the gap transformation and the related DTL cell length calculation used in the PARMILA code.

43 PARTICLE ACCELERATORS↗

Reconstruction of implanted marker trajectories from cone-beam CT projection images using interdimensional correlation modeling

Purpose: Cone-beam CT (CBCT) is a widely used imaging modality for image-guided radiotherapy. Most vendors provide CBCT systems that are mounted on a linac gantry. Thus, CBCT can be used to estimate the actual 3-dimensional (3D) position of moving respiratory targets in the thoracic/abdominal region using 2D projection images. The authors have developed a method for estimating the 3D trajectory of respiratory-induced target motion from CBCT projection images using interdimensional correlation modeling. Methods: Because the superior–inferior (SI) motion of a target can be easily analyzed on projection images of a gantry-mounted CBCT system, the authors investigated the interdimensional correlation of the SI motion with left–right and anterior–posterior (AP) movements while the gantry is rotating. A simple linear model and a state-augmented model were implemented and applied to the interdimensional correlation analysis, and their performance was compared. The parameters of the interdimensional correlation models were determined by least-square estimation of the 2D error between the actual and estimated projected target position. The method was validated using 160 3D tumor trajectories from 46 thoracic/abdominal cancer patients obtained during CyberKnife treatment. The authors’ simulations assumed two application scenarios: (1) retrospective estimation for the purpose of moving tumor setup used just after volumetric matching with CBCT; and (2) on-the-fly estimation for the purpose of real-time target position estimation during gating or tracking delivery, either for full-rotation volumetric-modulated arc therapy (VMAT) in 60 s or a stationary six-field intensity-modulated radiation therapy (IMRT) with a beam delivery time of 20 s. Results: For the retrospective CBCT simulations, the mean 3D root-mean-square error (RMSE) for all 4893 trajectory segments was 0.41 mm (simple linear model) and 0.35 mm (state-augmented model). In the on-the-fly simulations, prior projections over more than 60° appear to be necessary for reliable estimations. The mean 3D RMSE during beam delivery after the simple linear model had established with a prior 90° projection data was 0.42 mm for VMAT and 0.45 mm for IMRT. Conclusions: The proposed method does not require any internal/external correlation or statistical modeling to estimate the target trajectory and can be used for both retrospective image-guided radiotherapy with CBCT projection images and real-time target position monitoring for respiratory gating or tracking.

60 APPLIED LIFE SCIENCES↗

Nonthermally Dominated Electron Acceleration during Magnetic Reconnection in a Low-beta Plasma

This work was motivated by electron acceleration during solar flares. After some introductory remarks on proposed particle acceleration mechanisms and questions needing answers, dynamic simulations and simulation results are presented including energy spectra and the formation of the power law distribution. In summary, magnetic reconnection is highly efficient at converting the free magnetic energy stored in a magnetic shear and accelerating electrons to nonthermal energies in low-β regime. The nonthermal electrons have a dominant fraction and form power-law energy spectra with spectral index p ~ 1 in low-β regime. Electrons are preferentially accelerated along the curvature drift direction along the electric field induced by the reconnection outflow. The results can be applied to explain the observations of electron acceleration during solar flares.

79 ASTRONOMY AND ASTROPHYSICS↗

A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

electric propulsion↗

Interstellar clouds in high-speed, supersonic flows: Two-dimensional simulations

We present a series of gasdynamical simulations of the interaction of a dense, cool interstellar cloud with a high-speed, supersonic wind that confines and accelerates the embedded cloud. Our goal is to attempt to determine if such clouds can survive various potentially disruptive instabilities, that occur at their peripheries, long enough to be accelerated to speeds which are comparable to the wind velocity. These simulations are performed using two-dimensional, Eulerian gas dynamics on both an axisymmetric (about the cloud axis) and 'slab' geometric grid. The spatial and temporal resolutions of the simulations are varied over a wide range to investigate the effects of small-scale instabilities on the overall acceleration of clouds and the development of large-scale, disruptive instabilities. Also, we study the effects of wind/cloud Mach number variations by changing the wind speed constant at about 12 km/s (which corresponds to a cloud temperature of 10,000 K). The current simulations track the evolution of clouds as they are accelerated to speeds approximately 4-5 times greater than their internal sound speeds. Furthermore, the models with the highest resolution were extended far beyond quasi-linear Rayleigh-Taylor growth times reaching 6-7 Rayleigh-Taylor growth times for the largest scale instabilities before being terminated because of the accumulation of errors at the rear grid boundary.

Schiano, A. V. R.↗