An efficient transport solver for tokamak plasmas

Park, Jin Myung; Murakami, Masanori; St. John, H. E.; Lao, Lang L.; Chu, M. S.; Prater, Ronald

doi:10.1016/j.cpc.2016.12.018

Title: An efficient transport solver for tokamak plasmas

Abstract

A simple approach to efficiently solve a coupled set of 1-D diffusion-type transport equations with a stiff transport model for tokamak plasmas is presented based on the 4th order accurate Interpolated Differential Operator scheme along with a nonlinear iteration method derived from a root- finding algorithm. Numerical tests using the Trapped Gyro-Landau-Fluid model show that the presented high order method provides accurate transport solution using a small number of grids with robust nonlinear convergence.

Authors:

Park, Jin Myung ^[1]; Murakami, Masanori ^[1]; St. John, H. E. ^[2]; Lao, Lang L. ^[2]; Chu, M. S. ^[2]; Prater, Ronald ^[2]

Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
General Atomics, San Diego, CA (United States)

Publication Date:: Tue Jan 03 00:00:00 EST 2017

Research Org.:: Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States); General Atomics, San Diego, CA (United States)

Sponsoring Org.:: USDOE Office of Science (SC), Fusion Energy Sciences (FES)

OSTI Identifier:: 1344254

Alternate Identifier(s):: OSTI ID: 1375926; OSTI ID: 1396489

Grant/Contract Number:: AC05-00OR22725; FC02-04ER54698

Resource Type:: Accepted Manuscript

Journal Name:: Computer Physics Communications

Additional Journal Information:: Journal Volume: 214C; Journal ID: ISSN 0010-4655

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; 97 MATHEMATICS AND COMPUTING; plasma simulation; diffusion equation; transport; tokamak

Citation Formats


                    Park, Jin Myung, Murakami, Masanori, St. John, H. E., Lao, Lang L., Chu, M. S., and Prater, Ronald. An efficient transport solver for tokamak plasmas.  United States: N. p., 2017. 
Web.  doi:10.1016/j.cpc.2016.12.018.

Copy to clipboard


                    Park, Jin Myung, Murakami, Masanori, St. John, H. E., Lao, Lang L., Chu, M. S., & Prater, Ronald. An efficient transport solver for tokamak plasmas.  United States.  https://doi.org/10.1016/j.cpc.2016.12.018

Copy to clipboard


                    Park, Jin Myung, Murakami, Masanori, St. John, H. E., Lao, Lang L., Chu, M. S., and Prater, Ronald. Tue .  
"An efficient transport solver for tokamak plasmas".  United States.  https://doi.org/10.1016/j.cpc.2016.12.018.  https://www.osti.gov/servlets/purl/1344254.

Copy to clipboard


                    
@article{osti_1344254,

  title        = {An efficient transport solver for tokamak plasmas},

  author       = {Park, Jin Myung and Murakami, Masanori and St. John, H. E. and Lao, Lang L. and Chu, M. S. and Prater, Ronald},

  abstractNote = {A simple approach to efficiently solve a coupled set of 1-D diffusion-type transport equations with a stiff transport model for tokamak plasmas is presented based on the 4th order accurate Interpolated Differential Operator scheme along with a nonlinear iteration method derived from a root- finding algorithm. Numerical tests using the Trapped Gyro-Landau-Fluid model show that the presented high order method provides accurate transport solution using a small number of grids with robust nonlinear convergence.},

  doi          = {10.1016/j.cpc.2016.12.018},

  journal      = {Computer Physics Communications},

  number       = ,

  volume       = 214C,

  place        = {United States},

  year         = {Tue Jan 03 00:00:00 EST 2017},

  month        = {Tue Jan 03 00:00:00 EST 2017}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1016/j.cpc.2016.12.018

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 12 works

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Works referenced in this record:

A gyro-Landau-fluid transport model
journal, July 1997

Waltz, R. E.; Staebler, G. M.; Dorland, W.
Physics of Plasmas, Vol. 4, Issue 7
DOI: 10.1063/1.872228

Predicting temperature and density profiles in tokamaks
journal, May 1998

Bateman, Glenn; Kritz, Arnold H.; Kinsey, Jon E.
Physics of Plasmas, Vol. 5, Issue 5
DOI: 10.1063/1.872848

A theory-based transport model with comprehensive physics
journal, May 2007

Staebler, G. M.; Kinsey, J. E.; Waltz, R. E.
Physics of Plasmas, Vol. 14, Issue 5
DOI: 10.1063/1.2436852

Neoclassical transport of impurities in tokamak plasmas
journal, September 1981

Hirshman, S. P.; Sigmar, D. J.
Nuclear Fusion, Vol. 21, Issue 9
DOI: 10.1088/0029-5515/21/9/003

On 1D diffusion problems with a gradient-dependent diffusion coefficient
journal, October 2008

Jardin, S. C.; Bateman, G.; Hammett, G. W.
Journal of Computational Physics, Vol. 227, Issue 20
DOI: 10.1016/j.jcp.2008.06.032

Tokamak profile prediction using direct gyrokinetic and neoclassical simulation
journal, June 2009

Candy, J.; Holland, C.; Waltz, R. E.
Physics of Plasmas, Vol. 16, Issue 6
DOI: 10.1063/1.3167820

Interpolated differential operator (IDO) scheme for solving partial differential equations
journal, May 1997

Aoki, Takayuki
Computer Physics Communications, Vol. 102, Issue 1-3
DOI: 10.1016/S0010-4655(97)00020-9

Progress toward fully noninductive discharge operation in DIII-D using off-axis neutral beam injection
journal, September 2013

Ferron, J. R.; Holcomb, C. T.; Luce, T. C.
Physics of Plasmas, Vol. 20, Issue 9
DOI: 10.1063/1.4821072

Steady state scenario development with elevated minimum safety factor on DIII-D
journal, August 2014

Holcomb, C. T.; Ferron, J. R.; Luce, T. C.
Nuclear Fusion, Vol. 54, Issue 9
DOI: 10.1088/0029-5515/54/9/093009

Fusion Nuclear Science Facility (FNSF) before Upgrade to Component Test Facility (CTF)
journal, August 2011

Peng, Y. K. M.; Canik, J. M.; Diem, S. J.
Fusion Science and Technology, Vol. 60, Issue 2
DOI: 10.13182/FST60-441

Integrated modelling of steady-state scenarios and heating and current drive mixes for ITER
journal, August 2011

Murakami, M.; Park, J. M.; Giruzzi, G.
Nuclear Fusion, Vol. 51, Issue 10
DOI: 10.1088/0029-5515/51/10/103006

On the heating mix of ITER
journal, November 2010

Wagner, F.; Becoulet, A.; Budny, R.
Plasma Physics and Controlled Fusion, Vol. 52, Issue 12
DOI: 10.1088/0741-3335/52/12/124044

Effect of finite aspect ratio on the neoclassical ion thermal conductivity in the banana regime
journal, January 1982

Chang, C. S.
Physics of Fluids, Vol. 25, Issue 9
DOI: 10.1063/1.863934

The tokamak Monte Carlo fast ion module NUBEAM in the National Transport Code Collaboration library
journal, June 2004

Pankin, Alexei; McCune, Douglas; Andre, Robert
Computer Physics Communications, Vol. 159, Issue 3
DOI: 10.1016/j.cpc.2003.11.002

Electron cyclotron current drive efficiency in general tokamak geometry
journal, October 2003

Lin-Liu, Y. R.; Chan, V. S.; Prater, R.
Physics of Plasmas, Vol. 10, Issue 10
DOI: 10.1063/1.1610472

Equilibrium analysis of current profiles in tokamaks
journal, June 1990

Lao, L. L.; Ferron, J. R.; Groebner, R. J.
Nuclear Fusion, Vol. 30, Issue 6
DOI: 10.1088/0029-5515/30/6/006

Works referencing / citing this record:

On benchmarking of simulations of particle transport in ITER
journal, June 2019

Na, Yong-Su; Koechl, F.; Polevoi, A. R.
Nuclear Fusion, Vol. 59, Issue 7
DOI: 10.1088/1741-4326/ab15e0

Machine-learning assisted steady-state profile predictions using global optimization techniques
journal, October 2019

Honda, M.; Narita, E.
Physics of Plasmas, Vol. 26, Issue 10
DOI: 10.1063/1.5117846

Integrated modeling of high β _N steady state scenario on DIII-D
journal, January 2018

Park, J. M.; Ferron, J. R.; Holcomb, C. T.
Physics of Plasmas, Vol. 25, Issue 1
DOI: 10.1063/1.5013021

Similar Records in DOE PAGES and OSTI.GOV collections:

Alfvén eigenmode stability and critical gradient energetic particle transport using the Trapped-Gyro-Landau-Fluid model

Journal Article Sheng, He ; Waltz, R. E. ; Staebler, G. M. - Physics of Plasmas

The Trapped-Gyro-Landau-Fluid (TGLF) transport model is a physically realistic and comprehensive theory based on a local quasilinear transport model fitted to linear and nonlinear GYRO gyrokinetic simulations [Staebler et al., Phys. Plasmas 14, 55909 (2007)]. This work presents the first use of the TGLF model to treat low-n Alfvén eigenmode (AE) stability and energetic particle (EP) transport. TGLF accurately recovers the local GYRO toroidicity-induced AE (TAE) and energetic particle mode (EPM) linear growth and frequency rates for a fusion alpha case. With a very high grid resolution, TGLF can quickly find the critical EP pressure gradient profile for stiff EP transport based on an AE linear threshold given the background thermal plasma profiles in DIII-D. The TGLF critical gradient profile using the recipe γ_AE=0, that is the linear AE growth rate without additional driving rates from the background plasma gradients, matches the more expensive linear GYRO results with a single worst toroidal mode number n. TGLF can easily find the minimum critical gradient profile with testing multiple ns. From a database of runs using a newly developed TGLFEP code, a rough but insightful parametric “power law” scaling for critical EP beta is demonstrated. An important toroidal stabilization condition on the EP pressure gradient p_EP/Lmore »« less
Cited by 18
https://doi.org/10.1063/1.4989716

Full Text Available
Steady-State Gyrokinetics Transport Code (SSGKT), A Scientific Application Partnership with the Framework Application for Core-Edge Transport Simulations, Final Report

Technical Report Fahey, Mark R. ; Candy, Jeff

This project initiated the development of TGYRO - a steady-state Gyrokinetic transport code (SSGKT) that integrates micro-scale GYRO turbulence simulations into a framework for practical multi-scale simulation of conventional tokamaks as well as future reactors. Using a lightweight master transport code, multiple independent (each massively parallel) gyrokinetic simulations are coordinated. The capability to evolve profiles using the TGLF model was also added to TGYRO and represents a more typical use-case for TGYRO. The goal of the project was to develop a steady-state Gyrokinetic transport code (SSGKT) that integrates micro-scale gyrokinetic turbulence simulations into a framework for practical multi-scale simulation ofmore »« less
https://doi.org/10.2172/1104936

Full Text Available
Predictions of the near edge transport shortfall in DIII-D L-mode plasmas using the trapped gyro-Landau-fluid model [Predictions of the near edge transport shortfall in DIII-D L-mode plasmas using the TGLF model]

Journal Article Kinsey, Jon E. ; Staebler, Gary M. ; Candy, Jefferey M. ; ... - Physics of Plasmas

Previous studies of DIII-D L-mode plasmas have shown that a transport shortfall exists in that our current models of turbulent transport can significantly underestimate the energy transport in the near edge region. In this paper, the Trapped Gyro-Landau-Fluid (TGLF) drift wave transport model is used to simulate the near edge transport in a DIII-D L-mode experiment designed to explore the impact of varying the safety factor on the shortfall. We find that the shortfall systematically increases with increasing safety factor and is more pronounced for the electrons than for the ions. Within the shortfall dataset, a single high current casemore »« less
Cited by 22
https://doi.org/10.1063/1.4905630

Full Text Available
Propagation of input parameter uncertainties in transport models

Journal Article Vaezi, P. ; Holland, C. ; Grierson, B. A. ; ... - Physics of Plasmas

Here, the many sources of uncertainty in validation studies of plasma turbulence in magnetically confined fusion devices are well-known. In this paper, we investigate how to efficiently transform uncertainties in experimentally derived transport model inputs into model prediction uncertainties, using the quasilinear trapped-gyro-Landau-fluid (TGLF) turbulent transport model. We use the rapidly converging and computationally inexpensive non-intrusive probabilistic collocation method (PCM) to propagate input parameter uncertainty probability distribution functions (PDFs) through TGLF, yielding PDFs of predicted transport fluxes. We observe in many cases that the flux PDFs exhibit significant non-normal features such as strong skewness, even when the input distributions weremore »« less
https://doi.org/10.1063/1.5053906

Full Text Available
Reduced electron thermal transport in low collisionality H-mode plasmas in DIII-D and the importance of TEM/ETG-scale turbulence

Journal Article Schmitz, L. ; Holland, C. ; Rhodes, T. L. ; ... - Nuclear Fusion

The first systematic investigation of core electron thermal transport and the role of local ion temperature gradient/trapped electron mode/electron temperature gradient (ITG/TEM/ETG)-scale core turbulence is performed in high temperature, low collisionality H-mode plasmas in the DIII-D tokamak. Wavenumber spectra of L-mode and H-mode density turbulence are measured by Doppler backscattering. H-mode wavenumber spectra are directly contrasted for the first time with nonlinear gyrokinetic simulation results. Core ITG/TEM-scale turbulence is substantially reduced/suppressed by E × B shear promptly after the L–H transition, resulting in reduced electron thermal transport across the entire minor radius. For small kθρs, both experiment and nonlinear gyrokineticmore »« less
https://doi.org/10.1088/0029-5515/52/2/023003

Similar Records

Title: An efficient transport solver for tokamak plasmas

Abstract

Citation Formats

A gyro-Landau-fluid transport model journal, July 1997

Predicting temperature and density profiles in tokamaks journal, May 1998

A theory-based transport model with comprehensive physics journal, May 2007

Neoclassical transport of impurities in tokamak plasmas journal, September 1981

On 1D diffusion problems with a gradient-dependent diffusion coefficient journal, October 2008

Tokamak profile prediction using direct gyrokinetic and neoclassical simulation journal, June 2009

Interpolated differential operator (IDO) scheme for solving partial differential equations journal, May 1997

Progress toward fully noninductive discharge operation in DIII-D using off-axis neutral beam injection journal, September 2013

Steady state scenario development with elevated minimum safety factor on DIII-D journal, August 2014

Fusion Nuclear Science Facility (FNSF) before Upgrade to Component Test Facility (CTF) journal, August 2011

Integrated modelling of steady-state scenarios and heating and current drive mixes for ITER journal, August 2011

On the heating mix of ITER journal, November 2010

Effect of finite aspect ratio on the neoclassical ion thermal conductivity in the banana regime journal, January 1982

The tokamak Monte Carlo fast ion module NUBEAM in the National Transport Code Collaboration library journal, June 2004

Electron cyclotron current drive efficiency in general tokamak geometry journal, October 2003

Equilibrium analysis of current profiles in tokamaks journal, June 1990

On benchmarking of simulations of particle transport in ITER journal, June 2019

Machine-learning assisted steady-state profile predictions using global optimization techniques journal, October 2019

Integrated modeling of high β N steady state scenario on DIII-D journal, January 2018

A gyro-Landau-fluid transport model
journal, July 1997

Predicting temperature and density profiles in tokamaks
journal, May 1998

A theory-based transport model with comprehensive physics
journal, May 2007

Neoclassical transport of impurities in tokamak plasmas
journal, September 1981

On 1D diffusion problems with a gradient-dependent diffusion coefficient
journal, October 2008

Tokamak profile prediction using direct gyrokinetic and neoclassical simulation
journal, June 2009

Interpolated differential operator (IDO) scheme for solving partial differential equations
journal, May 1997

Progress toward fully noninductive discharge operation in DIII-D using off-axis neutral beam injection
journal, September 2013

Steady state scenario development with elevated minimum safety factor on DIII-D
journal, August 2014

Fusion Nuclear Science Facility (FNSF) before Upgrade to Component Test Facility (CTF)
journal, August 2011

Integrated modelling of steady-state scenarios and heating and current drive mixes for ITER
journal, August 2011

On the heating mix of ITER
journal, November 2010

Effect of finite aspect ratio on the neoclassical ion thermal conductivity in the banana regime
journal, January 1982

The tokamak Monte Carlo fast ion module NUBEAM in the National Transport Code Collaboration library
journal, June 2004

Electron cyclotron current drive efficiency in general tokamak geometry
journal, October 2003

Equilibrium analysis of current profiles in tokamaks
journal, June 1990

On benchmarking of simulations of particle transport in ITER
journal, June 2019

Machine-learning assisted steady-state profile predictions using global optimization techniques
journal, October 2019

Integrated modeling of high β _N steady state scenario on DIII-D
journal, January 2018