Development of the neutral model in the nonlinear MHD code JOREK: Application to e × B drifts in ITER PFPO-1 plasmas

JOREK Team

doi:10.1063/5.0135318

Development of the neutral model in the nonlinear MHD code JOREK: Application to e × B drifts in ITER PFPO-1 plasmas

JOREK Team

Research output: Contribution to journal › Article › peer-review

11 Scopus citations

Abstract

The prediction of power fluxes and plasma-wall interactions impacted by MHD processes during ITER operation [disruption, Edge Localized Modes (ELMs), 3D magnetic fields applied for ELM control, etc.] requires models that include an accurate description of the MHD processes themselves, as well as of the edge plasma and plasma-wall interaction processes. In this paper, we report progress on improving the edge plasma physics models in the nonlinear extended MHD code JOREK, which has capabilities to simulate the MHD response of the plasma to the applied external 3D fields, disruptions and ELMs. The extended MHD model includes E × B drifts, diamagnetic drifts, and neoclassical flows. These drifts can have large influences, on e.g., divertor asymmetries. Realistic divertor conditions are important for impurity sputtering, transport, and their effect on the plasma. In this work, we implemented kinetic and fluid neutral physics modules, investigated the influence of poloidal flows under divertor conditions in the ITER PFPO-1 (1.8T/5MA) H-mode plasma scenario, and compared the divertor plasma conditions and heat flux to the wall for both the fluid and kinetic neutral model (in JOREK) to the well-established 2D boundary plasma simulation code suite SOLPS-ITER. As an application of the newly developed model, we investigated time-dependent divertor solutions and the transition from attached to partially detached plasmas. We present the formation of a high-field-side high-density-region and how it is driven by poloidal E × B drifts.

Original language	English
Article number	042509
Journal	Physics of Plasmas
Volume	30
Issue number	4
DOIs	https://doi.org/10.1063/5.0135318
State	Published - Apr 1 2023
Externally published	Yes

Funding

This work was part of the research programme (No. 19KE01) of the Foundation for Nederlandse Wetenschappelijk Onderzoek Instituten (NWO-I), which was part of the NWO and has been carried out within the framework of the EUROfusion Consortium, funded by the European Union via the Euratom Research and Training Programme (Grant Agreement No. 101052200 – EUROfusion) and by the ITER Organization under Implementing Agreement No. 43-2174 with the Eindhoven University of Technology (TU/e) and with TU/e work order No. 107640. ITER is the Nuclear Facility, INB No. 174. Computations were performed using computer clusters Marconi-Fusion and the ITER Scientific Data and Computing Center (SDCC). The views and opinions expressed herein do not necessarily reflect those of the European Commission nor the ITER Organization nor NWO-I. Discussions with members of the ITER Organization (R. A. Pitts, S. D. Pinches, and X. Bonnin) are acknowledged.

Funders	Funder number
Foundation for Nederlandse Wetenschappelijk Onderzoek Instituten
NWO-I
Instituto Tecnológico y de Energías Renovables	43-2174
European Commission	101052200 – EUROfusion
Technische Universiteit Eindhoven	107640
Nederlandse Organisatie voor Wetenschappelijk Onderzoek

UN SDGs

This output contributes to the following UN Sustainable Development Goals (SDGs)

Access to Document

10.1063/5.0135318

Cite this

@article{ab4f99c29bf84c7881a4c0fe8f24b9f7,

title = "Development of the neutral model in the nonlinear MHD code JOREK: Application to e × B drifts in ITER PFPO-1 plasmas",

abstract = "The prediction of power fluxes and plasma-wall interactions impacted by MHD processes during ITER operation [disruption, Edge Localized Modes (ELMs), 3D magnetic fields applied for ELM control, etc.] requires models that include an accurate description of the MHD processes themselves, as well as of the edge plasma and plasma-wall interaction processes. In this paper, we report progress on improving the edge plasma physics models in the nonlinear extended MHD code JOREK, which has capabilities to simulate the MHD response of the plasma to the applied external 3D fields, disruptions and ELMs. The extended MHD model includes E × B drifts, diamagnetic drifts, and neoclassical flows. These drifts can have large influences, on e.g., divertor asymmetries. Realistic divertor conditions are important for impurity sputtering, transport, and their effect on the plasma. In this work, we implemented kinetic and fluid neutral physics modules, investigated the influence of poloidal flows under divertor conditions in the ITER PFPO-1 (1.8T/5MA) H-mode plasma scenario, and compared the divertor plasma conditions and heat flux to the wall for both the fluid and kinetic neutral model (in JOREK) to the well-established 2D boundary plasma simulation code suite SOLPS-ITER. As an application of the newly developed model, we investigated time-dependent divertor solutions and the transition from attached to partially detached plasmas. We present the formation of a high-field-side high-density-region and how it is driven by poloidal E × B drifts.",

author = "{JOREK Team} and Korving, {S. Q.} and Huijsmans, {G. T.A.} and Park, {J. S.} and A. Loarte",

note = "Publisher Copyright: {\textcopyright} 2023 Author(s).",

year = "2023",

month = apr,

day = "1",

doi = "10.1063/5.0135318",

language = "English",

volume = "30",

journal = "Physics of Plasmas",

issn = "1070-664X",

publisher = "American Institute of Physics",

number = "4",

}

TY - JOUR

T1 - Development of the neutral model in the nonlinear MHD code JOREK

T2 - Application to e × B drifts in ITER PFPO-1 plasmas

AU - JOREK Team

AU - Korving, S. Q.

AU - Huijsmans, G. T.A.

AU - Park, J. S.

AU - Loarte, A.

PY - 2023/4/1

Y1 - 2023/4/1

N2 - The prediction of power fluxes and plasma-wall interactions impacted by MHD processes during ITER operation [disruption, Edge Localized Modes (ELMs), 3D magnetic fields applied for ELM control, etc.] requires models that include an accurate description of the MHD processes themselves, as well as of the edge plasma and plasma-wall interaction processes. In this paper, we report progress on improving the edge plasma physics models in the nonlinear extended MHD code JOREK, which has capabilities to simulate the MHD response of the plasma to the applied external 3D fields, disruptions and ELMs. The extended MHD model includes E × B drifts, diamagnetic drifts, and neoclassical flows. These drifts can have large influences, on e.g., divertor asymmetries. Realistic divertor conditions are important for impurity sputtering, transport, and their effect on the plasma. In this work, we implemented kinetic and fluid neutral physics modules, investigated the influence of poloidal flows under divertor conditions in the ITER PFPO-1 (1.8T/5MA) H-mode plasma scenario, and compared the divertor plasma conditions and heat flux to the wall for both the fluid and kinetic neutral model (in JOREK) to the well-established 2D boundary plasma simulation code suite SOLPS-ITER. As an application of the newly developed model, we investigated time-dependent divertor solutions and the transition from attached to partially detached plasmas. We present the formation of a high-field-side high-density-region and how it is driven by poloidal E × B drifts.

AB - The prediction of power fluxes and plasma-wall interactions impacted by MHD processes during ITER operation [disruption, Edge Localized Modes (ELMs), 3D magnetic fields applied for ELM control, etc.] requires models that include an accurate description of the MHD processes themselves, as well as of the edge plasma and plasma-wall interaction processes. In this paper, we report progress on improving the edge plasma physics models in the nonlinear extended MHD code JOREK, which has capabilities to simulate the MHD response of the plasma to the applied external 3D fields, disruptions and ELMs. The extended MHD model includes E × B drifts, diamagnetic drifts, and neoclassical flows. These drifts can have large influences, on e.g., divertor asymmetries. Realistic divertor conditions are important for impurity sputtering, transport, and their effect on the plasma. In this work, we implemented kinetic and fluid neutral physics modules, investigated the influence of poloidal flows under divertor conditions in the ITER PFPO-1 (1.8T/5MA) H-mode plasma scenario, and compared the divertor plasma conditions and heat flux to the wall for both the fluid and kinetic neutral model (in JOREK) to the well-established 2D boundary plasma simulation code suite SOLPS-ITER. As an application of the newly developed model, we investigated time-dependent divertor solutions and the transition from attached to partially detached plasmas. We present the formation of a high-field-side high-density-region and how it is driven by poloidal E × B drifts.

UR - http://www.scopus.com/inward/record.url?scp=85146963659&partnerID=8YFLogxK

U2 - 10.1063/5.0135318

DO - 10.1063/5.0135318

M3 - Article

AN - SCOPUS:85146963659

SN - 1070-664X

VL - 30

JO - Physics of Plasmas

JF - Physics of Plasmas

IS - 4

M1 - 042509

ER -

Development of the neutral model in the nonlinear MHD code JOREK: Application to e × B drifts in ITER PFPO-1 plasmas

Abstract

Funding

UN SDGs

Access to Document

Other files and links

Fingerprint

Cite this