Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TOKAMAK”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Study of the L-mode tokamak plasma “shortfall” with local and global nonlinear gyrokinetic δf particle-in-cell simulation

The δ f particle-in-cell code GEM is used to study the transport “shortfall” problem of gyrokinetic simulations. In local simulations, the GEM results confirm the previously reported simulation results of DIII-D [Holland et al., Phys. Plasmas 16, 052301 (2009)] and Alcator C-Mod [Howard et al., Nucl. Fusion 53, 123011 (2013)] tokamaks with the continuum code GYRO. Namely, for DIII-D the simulations closely predict the ion heat flux at the core, while substantially underpredict transport towards the edge; while for Alcator C-Mod, the simulations show agreement with the experimental values of ion heat flux, at least within the range of experimental error. Global simulations are carried out for DIII-D L-mode plasmas to study the effect of edge turbulence on the outer core ion heat transport. The edge turbulence enhances the outer core ion heat transport through turbulence spreading. However, this edge turbulence spreading effect is not enough to explain the transport underprediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Quasistationary Plasma Predator-Prey System of Coupled Turbulence, Drive, and Sheared E × B Flow During High Performance DIII-D Tokamak Discharges [A New, Quasi-stationary Plasma Predator-Prey System of Coupled Turbulence, Drive, and Sheared E × B Flow During High Performance DIII-D Tokamak Discharges]

A new, long-lived limit cycle oscillation (LCO) regime has been observed in the edge of near zero torque high-performance DIII-D tokamak plasma discharges. These LCOs are localized and comprised of density turbulence, gradient drives, and E X B velocity shear damping ( E and B are the local radial electric and total magnetic fields). Density turbulence sequentially acts as a predator (via turbulence transport) of profile gradients and a prey (via shear suppression) to the E X B velocity shear. Reported here for the first time, a unique spatiotemporal variation of the local E X B velocity which is found to be essential for the existence of this system. The LCO system is quasi-stationary, existing from 3 to 12 plasma energy confinement times (~30 to 900 LCO cycles) limited by hardware constraints. In conclusion, this plasma system appears to contribute strongly to the edge transport in these high-performance and transient-free plasmas as evident from oscillations in transport relevant edge parameters at LCO timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

New regime of low ion collisionality in the neoclassical equilibrium of tokamak plasmas

Here, the neoclassical description of an axisymmetric toroidal plasma equilibrium is formulated for an unconventionally low ordering of the collisionality that suits realistic thermonuclear fusion conditions. This requires a drift-kinetic analysis to the second order of the ion Larmor radius, which yields a new contribution to the leading solution for the non-Maxwellian part of the ion distribution function if the equilibrium geometry is not up-down symmetric. An explicit geometrical factor weighs this second Larmor-radius order, low-collisionality effect that modifies the neoclassical ion parallel flow, and the ion contribution to the bootstrap current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The locking and unlocking thresholds for tearing modes in a cylindrical tokamak

The locking and unlocking thresholds for tearing modes are in general different. In this work, the physics origin for this difference is illustrated from theory analysis, and a numerical procedure is developed to find both locking and unlocking thresholds. In particular, a new scaling law for the unlocking threshold that is valid in both weak and strong rotation regimes has been derived from the lowest amplitude of the RMP (resonant magnetic perturbation) allowed for the locked-mode solution. Above the unlocking threshold, the criterion for the phase-flip instability is extended to identify the entire locked-mode states. Two different regimes of the RMP amplitude in terms of the accessibility of the locked-mode states have been found. In the first regime, the locked-mode state may or may not be accessible depending on the initial conditions of an evolving island. In the second regime, the locked-mode state can always be reached regardless of the initial conditions of the tearing mode. The lowest RMP amplitude for the second regime is determined to be the mode-locking threshold. The different characteristics of the two regimes above the unlocking threshold reveal the underlying physics for the gap between the locking and unlocking thresholds and provide an explanation for the closely related and widely observed hysteresis phenomena in island evolution during the sweeping process of the RMP amplitude up and down across that threshold gap.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Discontinuous Galerkin methods for plasma physics in the scrape-off layer of tokamaks

A new parallel discontinuous Galerkin solver, called ArcOn, is developed to describe the intermittent turbulent transport of filamentary blobs in the scrape-off layer (SOL) of fusion plasma. The model is comprised of an elliptic subsystem coupled to two convection-dominated reaction–diffusion–convection equations. Upwinding is used for a class of numerical fluxes developed to accommodate cross product driven convection, and the elliptic solver uses SIPG, NIPG, IIPG, Brezzi, and Bassi–Rebay fluxes to formulate the stiffness matrix. A novel entropy sensor is developed for this system, designed for a space–time varying artificial diffusion/viscosity regularization algorithm. Some numerical experiments are performed to show convergence order on manufactured solutions, regularization of blob/streamer dynamics in the SOL given unstable parameterizations, long-time stability of modon (or dipole drift vortex) solutions arising in simulations of drift-wave turbulence, and finally the formation of edge mode turbulence in the scrape-off layer under turbulent saturation conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Verification of nonlinear particle simulation of radio frequency waves in tokamak

Nonlinear simulation model for radio frequency waves in fusion plasmas has been developed and verified using fully kinetic ion and drift kinetic electron. Ion cyclotron motion in the toroidal geometry is implemented using Boris push in the Boozer coordinates. Linear dispersion relation and nonlinear particle trapping are verified for the lower hybrid wave and ion Bernstein wave (IBW). Parametric decay instability is observed where a large amplitude pump wave decays into an IBW sideband and an ion cyclotron quasimode (ICQM). The ICQM induces an ion perpendicular heating, with a heating rate proportional to the pump wave intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Non-perturbative measurement of cross-field thermal diffusivity reduction at the O-point of 2/1 neoclassical tearing mode islands in the DIII-D tokamak

Neoclassical tearing modes (NTMs) often lead to the decrease of plasma performance and can lead to disruptions, which makes them a major impediment in the development of operating scenarios in present toroidal fusion devices. Recent gyrokinetic simulations predict a decrease of plasma turbulence and cross-field transport at the O-point of the islands, which in turn affects the NTM dynamics. In this paper, a heat transport model of magnetic islands employing spatially non-uniform cross-field thermal diffusivity (χ{sub ⊥}) is presented. This model is used to derive χ{sub ⊥} at the O-point from electron temperature data measured across 2/1 NTM islands in DIII-D. It was found that χ{sub ⊥} at the O-point is 1 to 2 orders of magnitude smaller than the background plasma transport, in qualitative agreement with gyrokinetic predictions. As the anomalously large values of χ{sub ⊥} are often attributed to turbulence driven transport, the reduction of the O-point χ{sub ⊥} is consistent with turbulence reduction found in recent experiments. Finally, the implication of reduced χ{sub ⊥} at the O-point on NTM dynamics was investigated using the modified Rutherford equation that predicts a significant effect of reduced χ{sub ⊥} at the O-point on NTM saturation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Integrated modeling applications for tokamak experiments with OMFIT

We developed one modeling framework for integrated tasks (OMFIT) is a comprehensive integrated modeling framework which will enable physics codes to interact in complicated workflows, and support scientists at all stages of the modeling cycle. The OMFIT development follows a unique bottom-up approach, where the framework design and capabilities organically evolve to support progressive integration of the components that are required to accomplish physics goals of increasing complexity. OMFIT provides a workflow for easily generating full kinetic equilibrium reconstructions that are constrained by magnetic and motional Stark effect measurements, and kinetic profile information that includes fast-ion pressure modeled by a transport code. It was found that magnetic measurements can be used to quantify the amount of anomalous fast-ion diffusion that is present in DIII-D discharges, and provide an estimate that is consistent with what would be needed for transport simulations to match the measured neutron rates. OMFIT was used to streamline edge-stability analyses, and evaluate the effect of resonant magnetic perturbation (RMP) on the pedestal stability, which have been found to be consistent with the experimental observations. Moreover, the development of a five-dimensional numerical fluid model for estimating the effects of the interaction between magnetohydrodynamic (MHD) and microturbulence, and its systematic verification against analytic models was also supported by the framework. OMFIT was used for optimizing an innovative high-harmonic fast wave system proposed for DIII-D. For a parallel refractive index n(parallel to) > 3, the conditions for strong electron-Landau damping were found to be independent of launched n(parallel to) and poloidal angle. OMFIT has been the platform of choice for developing a neural-network based approach to efficiently perform a non-linear multivariate regression of local transport fluxes as a function of local dimensionless parameters. Transport predictions for thousands of DIII-D discharges showed excellent agreement with the power balance calculations across the whole plasma radius and over a broad range of operating regimes. Concerning predictive transport simulations, the framework made possible the design and automation of a workflow that enables self-consistent predictions of kinetic profiles and the plasma equilibrium. It is found that the feedback between the transport fluxes and plasma equilibrium can significantly affect the kinetic profiles predictions. Finally, such a rich set of results provide tangible evidence of how bottom-up approaches can potentially provide a fast track to integrated modeling solutions that are functional, cost-effective, and in sync with the research effort of the community.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

2D analysis of tokamak divertor-plasma detachment-bifurcation with operational parameters and geometries

UEDGE simulations with density scans for various input power, transport coefficients and outer poloidal leg length are performed to study the conditions for the existence of a bifurcation-like drop of T e at the outer strike point, commonly referred to as a detachment cliff, when transitioning to a detached plasma from an attached plasma in the outer divertor as the upstream density increases (McLean et al., 2015). The simulation results show that a detachment cliff tends to occur with a higher power input regardless of diffusivities and leg length. Further analysis of change of plasma profiles at a cliff indicate that, in addition to the sharp reduction of the E x B drift fluxes in the outer divertor studied in Jaervinen et al., (2018), the substantial change of the Mach number in the outer divertor and the decrease of the outer mid-plane T e due to the radiation front moving across the separatrix into the confinement region above the X-point consistently occur for all UEDGE density scans that have a detachment cliff. UEDGE time-dependent simulation of the evolution of a detachment cliff shows that the rapid increase of radiation above the X-point occurs in a time scale of ~0.3–0.5, which could possibly be the trigger for the formation of a detachment cliff, quicker than the Mach number change in a time scale of ~1 ms and the drop of T e in a time scale of ~2–3 ms in the outer divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Modeling of transport phenomena in tokamak plasmas with neural networks

A new transport model that uses neural networks (NNs) to yield electron and ion heat ux pro les has been developed. Given a set of local dimensionless plasma parameters similar to the ones that the highest delity models use, the NN model is able to efficiently and accurately predict the ion and electron heat transport pro les. As a benchmark, a NN was built, trained, and tested on data from the 2012 and 2013 DIII-D experimental campaigns. It is found that NN can capture the experimental behavior over the majority of the plasma radius and across a broad range of plasma regimes. Although each radial location is calculated independently from the others, the heat ux pro les are smooth, suggesting that the solution found by the NN is a smooth function of the local input parameters. This result supports the evidence of a well-de ned, non-stochastic relationship between the input parameters and the experimentally measured transport uxes. Finally, the numerical efficiency of this method, requiring only a few CPU-μs per data point, makes it ideal for scenario development simulations and real-time plasma control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

In situ measurement of low-Z material coating thickness on high Z substrate for tokamaks

Rutherford backscattering of energetic particles can be used to determine the thickness of a coating of a low-Z material over a heavier substrate. Simulations indicate that 5 MeV alpha particles from an {sup 241}Am source can be used to measure the thickness of a Li coating on Mo tiles between 0.5 and 15 μm thick. Using a 0.1 mCi source, a thickness measurement can be accomplished in 2 h of counting. This technique could be used to measure any thin, low-Z material coating (up to 1 mg/cm{sup 2} thick) on a high-Z substrate, such as Be on W, B on Mo, or Li on Mo. By inserting a source and detector on a moveable probe, this technique could be used to provide an in situ measurement of the thickness of Li coating on NSTX-U Mo tiles. A test stand with an alpha source and an annular solid-state detector was used to investigate the measurable range of low-Z material thicknesses on Mo tiles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗