Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ORBITAL SIMULATOR”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Full-orbit simulations of fast-ion charge-exchange losses induced by neutral particles outside the plasma last closed flux surface

The interaction between neutrals outside tokamak plasmas and fast ions that are generated by neutral beam injection and pass outside the plasma is studied. It is shown that at beam energies (up to ~100~keV) fast-ion losses are significant, up to 20% in the studied case, because of charge exchange reactions between the neutrals and the fast ions. About half of the fast ions that neutralize outside the plasma get lost to the wall while the other half re-enter the plasma and ionize again, often on better confined orbits which is visible as an enhanced beam power deposition in the plasma core (rho ~< 0.7). As a result, the the plasma performance as measured from the beam-thermal neutrons is hardly changed despite the significant fast-ion losses toward the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

ORBIT simulations of fast ion power loads on the wall of the Divertor Tokamak Test

Abstract Neutral beam injection in tokamaks produces a population of fast ions, which interact with 3D magnetic fields in a variety of ways, often resulting in energetic particle losses in very short times to the wall. Careful design of neutral beams and active control of error fields helps to keep these losses to a minimum. Nevertheless, past experience in tokamaks in the 1980s and detailed simulations for future machines, such as International Thermonuclear Experimental Reactor, suggest that very localized fast ion losses (‘hot spots’) can be present, even if the overall losses are low. In this paper, we discuss this issue in the Divertor Tokamak Test (DTT) project, and we show that in the standard single-null full-power scenario of the DTT, fast ions produce two hot spots, corresponding approximately to the beam injection and exit toroidal angles: the former being mainly due to prompt losses/passing particles, while the latter is due to trapped ions. However, the maximum power load in these spots is of the order ∼ 63 kW m −2 , below the tolerance of plasma-facing components of the machine.

Physics↗
DOE OSTI2024

Orbit-following simulations of fast-ion transport and losses due to the Alfvén eigenmode burst in the Large Helical Device

Orbit-following simulations of fast-ion transport and losses with time-dependent electromagnetic perturbations are performed to clarify the roles of Alfvén eigenmodes (AEs) and the low-frequency magnetohydrodynamic (MHD) mode observed in the kinetic-MHD hybrid simulation of AE bursts in the Large Helical Device. Fast-ion pressure profile flattening in the kinetic-MHD hybrid simulation can be reproduced by an orbit-following simulation with only the primary single AE of the time-dependent amplitude following the kinetic-MHD hybrid simulation result, while orbit-following simulations with constant AE amplitude of average level during AE burst cannot reproduce the fast-ion pressure profile flattening observed. The effects of other modes are negligible on the fast-ion pressure profile flattening. The fast-ion losses in kinetic-MHD hybrid simulation can be reproduced by an orbit-following simulation with time-dependent amplitude when the low-frequency MHD mode is considered in addition to multiple AEs. This indicates the synergetic effect of multiple AEs and the low-frequency MHD mode on fast-ion losses.

Seki, Ryohsuke↗
DOE OSTI2021

Fast ion transport by sawtooth instability in the presence of ICRF-NBI synergy in JET plasmas

JET experiments have shown that the three-ion scenarios using waves in the ion cyclotron range of frequencies (ICRF) is an efficient way to build fast ion population through beam ion acceleration by radio frequency (RF) waves. Here, such a heating scheme is applied to plasmas with at least two thermal ion species. Analysis of mixed discharges with complex heating schemes requires a workflow that allows to model thermal and fast ion transport consistently. This paper is dedicated to modelling of a mixed plasma discharge with significant fraction of fast ions and contributes to development of fast ion transport models. For interpretive analysis with the TRANSP code a JET hydrogen-deuterium (HD) plasma discharge with neutral beam injection (NBI) and ICRF heating has been chosen. The task is complicated by NBI-ICRF synergy and plasma magnetohydrodynamic activity, like sawtooth crashes. D beam ions accelerated by RF waves form a high energy tail in fast ion distribution. Significant difference between the neutron rate computed by TRANSP and measured one is observed if the same diffusivity for electrons and ions is assumed. Sensitivity studies show that uncertainties in input plasma parameters and thermal ion transport models are crucial for modelling mixed plasma discharges and increased D transport is required to reach the plasma composition consistent with diagnostic measurements at the plasma edge. Fast ion redistribution by a sawtooth instability is characterised by non-resonant transport due to reconnection of magnetic field lines and resonant transport caused by resonance interaction between the instability and fast ions. With ORBIT simulations it has been shown that resonant interaction strongly affects fast ions of high energies, like beam ions accelerated by RF waves and fusion products. For the considered case, fast ion profiles simulated by ORBIT remain peaked after the sawtooth crashes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Towards accurate orbital-free simulations: A generalized gradient approximation for the noninteracting free energy density functional

For orbital-free ab initio molecular dynamics, especially on systems in extreme thermodynamic conditions, we provide the first pseudo-potential-adapted generalized gradient approximation (GGA) functional for the non-interacting free energy. This is achieved by systematic finite-temperature extension of our recent LKT ground state non-interacting kinetic energy GGA functional (Phys. Rev. B 98, 041111(R) (2018)). We test the performance of the new functional first via static lattice calculations on crystalline aluminum and silicon. Then we compare deuterium equation of state results against both path-integral Monte Carlo and conventional (orbital-dependent) Kohn-Sham results. The new functional, denoted LKTF, outperforms the previous best semi-local free energy functional, VT84F (Phys. Rev. B 88, 161108(R) (2013)), and provides modestly faster simulations. Finally, we discuss subtleties of identification of kinetic and entropic contributions to non-interacting free-energy functionals obtained by extension from ground state orbital-free kinetic energy functionals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Target and orbit feedback simulations of a muSR beamline at BNL

Well-polarized positive surface muons are a tool to measure the magnetic properties of materials since the precession rate of the spin can be determined from the observation of the positron directions when the muons decay. The use of the AGS complex at BNL has been explored for a muSR facility previously. Here we report simulations of a beamline with a target inside a solenoidal field, and of an orbit feed-back system with single muon beam positioning monitors based on technology available today

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Reinventing atomic magnetic simulations with spin-orbit coupling

We propose a powerful extension to the combined molecular and spin dynamics method that fully captures the coupling between the atomic and spin subsystems via spin-orbit interactions. Moreover, the foundation of this method lies in the inclusion of the local magnetic anisotropies that arise as a consequence of the lattice symmetry breaking due to phonons or crystallographic defects. By using canonical simulations of bcc iron with the system coupled to a phonon heat bath, we show that our extension enables the previously unachievable angular momentum exchange between the atomic and spin degrees of freedom.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Target and orbit feedback simulations of a muSR beam line at BNL

Well-polarized positive surface muons are a tool to measure the magnetic properties of materials since the precession rate of the spin can be determined from the observation of the positron directions when the muons decay. For a dc beam an ideal µSR flux for surface µ + should be about 40 kHz/mm 2 . In this report we show how this flux could be achieved in a beam line using the AGS complex at BNL for a source of protons. We also determined that an orbit feedback system with a pair of thin silicon position monitors and kickers would miss the desired flux by at least an order of magnitude, even with perfect time resolution and no multiple scattering.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Conservative magnetic moment of runaway electrons and collisionless pitch-angle scattering

Recently, the validity of the guiding-center approach to model relativistic runaway electrons in tokamaks has been challenged by full-orbit simulations that demonstrate the breakdown of the standard magnetic moment conservation. In this paper, we derive a new expression for the magnetic moment of relativistic runaway electrons with p∥>>⊥, which is conserved significantly better than the standard one. The new result includes one of the second-order corrections in the standard guiding-center theory which, in case of runaway electrons with p∥>>⊥, can peculiarly be of the same order as the lowest-order term. The better conservation of the new magnetic moment also explains the collisionless pitch-angle-scattering effect observed in full-orbit simulations since it allows momentum transfer between the perpendicular and parallel directions when the runaway electron is accelerated by an electric field. In conclusion, while the derivation of the second-order correction to the magnetic moment in general case would require the full extent of the relativistic second-order guiding-center theory, we exploit the Lie-perturbation method at the limit p∥Gtp⊥ which simplifies the computations significantly. Consequently, we present the corresponding guiding-center equations applicable to runaway electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

SciDAC Center for Simulation of Wave-Plasma Interactions - Iterated Finite-Orbit Monte Carlo Simulations with Full-Wave Fields for Modeling Tokamak ICRF Wave Heating Experiments - Final Report

This final report describes the work performed under U.S. Department of Energy Cooperative Agreement DE-FC02-08ER54954 for the period April 1, 2011 through March 31, 2013. The goal of this project was to perform iterated finite-orbit Monte Carlo simulations with full-wall fields for modeling tokamak ICRF wave heating experiments. In year 1, the finite-orbit Monte-Carlo code ORBIT-RF and its iteration algorithms with the full-wave code AORSA were improved to enable systematical study of the factors responsible for the discrepancy in the simulated and the measured fast-ion FIDA signals in the DIII-D and NSTX ICRF fast-wave (FW) experiments. In year 2, ORBIT-RF was coupled to the TORIC full-wave code for a comparative study of ORBIT-RF/TORIC and ORBIT-RF/AORSA results in FW experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Method to integrate full particle orbit in toroidal plasmas

It is vital to integrate full particle orbit accurately when studying charged particle dynamics in electromagnetic waves with frequency higher than cyclotron frequency. We have derived a form of the Boris scheme using magnetic coordinates, which can be used effectively to integrate the cyclotron orbit in toroidal geometry over a long period of time. The proposed method has been verified by a full particle orbit simulation in toroidal geometry without high frequency waves. The full particle orbit calculation recovers guiding center banana orbit. This approach has better numeric properties than the conventional Runge-Kutta method for conserving particle energy and magnetic moment. The toroidal precession frequency is discovered to match that from guiding center simulation. Many other important phenomena in the presence of an electric field, such as E×B drift, Ware pinch effect and neoclassical polarization drift are also verified by the full orbit simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

A novel adaptive time stepping variant of the Boris–Buneman integrator for the simulation of particle accelerators with space charge

We show that adaptive time stepping in particle accelerator simulation is an enhancement for certain problems. The new algorithm has been implemented in the OPAL (Object Oriented Parallel Accelerator Library) framework. The idea is to adjust the frequency of costly self-field calculations, which are needed to model Coulomb interaction (space charge) effects. In analogy to a Kepler orbit simulation that requires a higher time step resolution at the close encounter, we propose to choose the time step based on the magnitude of the space charge forces. Inspired by geometric integration techniques, our algorithm chooses the time step proportional to a function of the current phase space state instead of calculating a local error estimate like a conventional adaptive procedure. Building on recent work, a more profound argument is given on how exactly the time step should be chosen. An intermediate algorithm, initially built to allow a clearer analysis by introducing separate time steps for external field and self-field integration, turned out to be useful by its own, for a large class of problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Constraining the range of Yukawa gravity interaction from S2 star orbits

We consider possible signatures for Yukawa gravity within the Galactic Central Parsec, based on our analysis of the S2 star orbital precession around the massive compact dark object at the Galactic Centre, and on the comparisons between the simulated orbits in Yukawa gravity and two independent sets of observations. Our simulations resulted in strong constraints on the range of Yukawa interaction Λ and showed that its most probable value in the case of S2 star is around 5000 - 7000 AU. At the same time, we were not able to obtain reliable constrains on the universal constant δ of Yukawa gravity, because the current observations of S2 star indicated that it may be highly correlated with parameter Λ in the range (0 < δ < 1). For δ > 2 they are not correlated. However, the same universal constant which was successfully applied to clusters of galaxies and rotation curves of spiral galaxies (δ = 1/3) also gives a satisfactory agreement with the observed orbital precession of the S2 star, and in that case the most probable value for the scale parameter is Λ ≈ 3000±1500 AU. Also, the Yukawa gravity potential induces precession of S2 star orbit in the same direction as General Relativity for δ > 0 and for δ < −1, and in the opposite direction for −1 < δ < 0. The future observations with advanced facilities, such as GRAVITY or/and European Extremely Large Telescope, are needed in order to verify these claims.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Fast ion transport in the quasi-single helical reversed-field pinch

The reversed-field pinch (RFP) can spontaneously transition from an axisymmetric magnetic topology to a 3D-helical geometry. Investigations on fast ion transport associated with energetic particle driven Alfvén instabilities, tearing mode induced stochasticity, and neoclassical effects have been performed on the Madison Symmetric Torus. STELLGAP produced shear-Alfvén continua seeded with V3FIT 3D-equilibrium reconstructions describe the response of Alfvénic bursting activity as a direct consequence of the equilibrium change on the fast ion resonance. Furthermore, far infrared interferometry resolved electron density perturbations associated with the bursts provide a spatial measurement of the mode structure and support the reconstructions. The bursts produce no global resonant fast ion transport; however, their disappearance at a high core-resonant amplitude implies other transport mechanisms at play. Additionally, neutral particle analysis and neutron signals suggest fast ion losses at sufficient core tearing mode strength, supporting the lack of Alfvénic activity. The guiding-center code ORBIT corroborates rapid fast ion loss times in the helical state largely as a consequence of remnant tearing modes. Additionally, ORBIT simulations demonstrate little neoclassical enhancement of particle transport. While superbanana orbits may exist, the growth in the core-resonant fast ion island and the associated secondary mode overlap govern the largest transport process, leading to robust fast ion losses in the 3D-RFP.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

The 3D flow field around an embedded planet

3D modifications to the well-studied 2D flow topology around an embedded planet have the potential to resolve long-standing problems in planet formation theory. We present a detailed analysis of the 3D isothermal flow field around a 5 Earth-mass planet on a fixed circular orbit, simulated using our graphics processing unit hydrodynamics code PEnGUIn. We find that, overall, the horseshoe region has a columnar structure extending vertically much beyond the Hill sphere of the planet. This columnar structure is only broken for some of the widest horseshoe streamlines, along which high altitude fluid descends rapidly into the planet’s Bondi sphere, performs one horseshoe turn, and exits the Bondi sphere radially in the midplane. A portion of this flow exits the horseshoe region altogether, which we refer to as the “transient” horseshoe flow. The flow continues as it rolls up into a pair of up–down symmetric horizontal vortex lines shed into the wake of the planet. This flow, unique to 3D, affects both planet accretion and migration. It prevents the planet from sustaining a hydrostatic atmosphere due to its intrusion into the Bondi sphere, and leads to a significant corotation torque on the planet, unanticipated by 2D analysis. In the reported simulation, starting with a Σ∼r{sup −3/2} radial surface density profile, this torque is positive and partially cancels with the negative differential Lindblad torque, resulting in a factor of three slower planet migration rate. Finally, we report 3D effects can be suppressed by a sufficiently large disk viscosity, leading to results similar to 2D.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Conceptual design of the ITER fast-ion loss detector

A conceptual design of a reciprocating fast-ion loss detector for ITER has been developed and is presented here. Fast-ion orbit simulations in a 3D magnetic equilibrium and up-to-date first wall have been carried out to revise the measurement requirements for the lost alpha monitor in ITER. In agreement with recent observations, the simulations presented here suggest that a pitch-angle resolution of ∼5° might be necessary to identify the loss mechanisms. Synthetic measurements including realistic lost alpha-particle as well as neutron and gamma fluxes predict scintillator signal-to-noise levels measurable with standard light acquisition systems with the detector aperture at ∼11 cm outside of the diagnostic first wall. At measurement position, heat load on detector head is comparable to that in present devices.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗