The 3D edge transport code EMC3-EIRENE has been applied for the first time to the NSTX spherical tokamak. A new disconnected double null grid has been developed to allow the simulation of plasma where the radial separation of the inner and outer separatrix is less than characteristic widths (e.g. heat flux width) at the midplane. Modeling results are presented for both an axisymmetric case and a case where 3D magnetic field is applied in an n=3 configuration. In the vacuum approximation the perturbed field consists of a wide region of destroyed flux surfaces and helical lobes which are a mixture of long and short connection length field lines formed by the separatrix manifolds. This structure is reflected in coupled 3D plasma fluid (EMC3) and kinetic neutral particle (EIRENE) simulations. The helical lobes extending inside of the unperturbed separatrix are filled in by hot plasma from the core. The intersection of the lobes with the divertor results in a striated flux footprint pattern on the target plates. Profiles of divertor heat and particle fluxes are compared to experimental data, and possible sources of discrepancy are discussed.
IntroductionNon-axisymmetric fields affect both the transport and stability properties of tokamak plasmas. Even a small asymmetry in the ideal axisymmetric field, for example due to unavoidable field errors, can result in increased collisional transport [1], creation of magnetic islands [2], stochastic fields [3], and rotation damping [4]. 3D field structures are also generated intentionally by perturbing the plasma with control coils. Application of 3D fields can cause several effects that are of critical importance to future tokamaks such as ITER. ELMs are mitigated, suppressed, or triggered [5][6][7][8][9]. The pedestal profiles are affected, with some cases showing a flattening of the density and temperature profiles, while in others the pedestal pressure gradient increases. The heat and particle divertor loads are also affected, through strike point splitting [10][11][12] or changes in the divertor recycling properties, e.g., reattachment of the divertor plasma [13]. These effects can be seen as either detrimental behavior to be mitigated or as beneficial behavior to be exploited. With this viewpoint, the critical aspect is to have an understanding of the underlying physics so that reliable predictions can be made, and full advantage of potential benefits taken.Modeling of the above effects is complicated due to the 3D geometry and the coupling of parallel and cross-field transport. The opportunity exists, however, to apply to tokamaks advanced 3D codes originally developed for stellarator research. The EMC3-EIRENE code [14-16] is a coupled 3D Monte Carlo plasma fluid transport (EMC3) and kinetic neutral recycling and transport (EIRENE) package. The code includes transport of particles, electron and ion energy, and parallel momentum in stochastic fields with 3D plasma facing component geometries. EMC3-EIRENE has been applied to stellarators [15] and limiter tokamaks ...