pankin | Integrated Modeling of Tokamak Plasmas

Full name:: Alexei Pankin
Nickname:: pankin
Website:
Description:

Posts by pankin

Two low-beta NSTX discharges that represent the collisionality scan are being analyzed

I have started the analysis of two low-beta H-mode NSTX discharges suggested by Stan Kaye. These discharges are

141031A10 – high collisionality
141040A05 – low collisionality

I have reproduced the old analysis runs. Also, I started to do predictive runs with the Multi-Mode Model MMM7.1:

141031P00/141040P00 – analysis runs. The results are the same as in the original TRANSP runs
141031P01/141040P01 – FIrst predictive runs with the default MMM7.1 settings. The effect of ExB flow shear in the DRIBM model is disabled
141031P02/141040P02 – Number of NZONES is increased from 40 to 80
141031P03/141040P03 – The ExB flow shear for DRIBM is enambed. The ExB factors for DRIBM and Weiland are set to the same level:
FACEXB          = 0.5     ! factor for ExB flow shear
LMMM07(4)    = 0 ! Switch for disabling ExB shear eects in DRIBM model
CMMM07(3)    = 0.5 ! ExB shear Multiplier in DRIBM model
141031P04/141040P04 – The ExB factors for DRIBM and Weiland are set to the default level that is typically used for other tokamaks:
FACEXB          = 1.0     ! factor for ExB flow shear
LMMM07(4)    = 0 ! Switch for disabling ExB shear eects in DRIBM model
CMMM07(3)    = 1.0 ! ExB shear Multiplier in DRIBM model

Some observations, conclusions, and plans

I have several observations:

The default settings for $k_yrho$ are suppose to correspond the most unstable modes. However, these setting were found in scans for conventional tomakak parameters (primary for H-mode parameters at DIII-D). These default settings do not necessarily correspond to the most unstable modes at NSTX.
The modes with different $k_yrho$ might be more unstable in different plasma regions. The DRIBM model with default parameter $k_yrho=0.2$ does not show any anomalous transport in the plasma core for $rho< 0.3$ . However, we found that there might be anomalous transport driven by smaller scale instabilities that have $k_yrho=0.6$ .

In order to address these observations, several improvements in the MMM8.1 model can be considered. These improvements include:

Loops over spectrum of $k_yrho$ needs to be implemented for DRIBM and Weiland models. Individual values of $k_yrho$ that correspond to maximum growth rates need to be selected for every flux surface for every numerical iteration. Depending on the $k_yrho$ ranges and granularity, these change will make the Multi-Mode model order of magnitude slower, but it will increase the robustness of the results.
The ExB flow shear effect should depend on $k_yrho$ value. There is no theory-based model at this time that can be implemented in MMM8.1. We can develop such a model by running gyro-kinetic simulations (for example, using the GYRO code). We can perform several scans with respect to $k_yrho$ values. Using these scan, we can develop a parametrized dependencies of the ExB flow shear factor with respect to $k_yrho$ .

These developments introduce significant changes to the model and are not currently funded. They might be included as new tasks to the NSTX proposal.

Next steps:

Compile the scan with respect to $k_yrho$ and ExB flow shear factor
Run the PTRANSP simulations with adjusted $k_yrho$ and ExB flows shear factors from t=0.364 sec to t=0.482 sec and to compare the predicted temperature profiles with the experimental observations.
Compare the diffusivities from the Multi-Mode model with the diffusivities from other models (TGLF).
Verify the ETG model results with other models (GLF23, TGLF and possibly GYRO).

Several new cases that complete the scan with respect to (k_y rho)

I have completed two additional cases with $k_yrho=0.6$ (case 141716H29/30) and with $k_yrho=6$ (case 141716H27/28) for DRIBM. Probably, the later case stretches the applicability of DRIBM model somewhat. I used the ExB flow shear factor 0.3 for the case 141716H29/30 and 0.1 for the case 141716H27/28. There very a small electron thermal transport predicted in the case 141716H27/28. The results for the case 141716H29/30 are given below:
There is large a level of anomalous transport for modes with $k_yrho=0.6$ . The level of anomalous transport exceeds predictions from the model that use the default settings $k_yrho=0.2$ and $k_yrho=0.316$ for Weiland.

Summary of recent TRANSP runs

Here is the summary of recent PTRANSP runs

Case #	Description	Comments
141716H09/10	ExB flow shear is enabled for W19 and disabled for DRIBM
141716H11/12	ExB flow shear is disabled for W19 and DRIBM
141716H13/14	MMM95 is used instead of MMM7.1	relatively small changes in diffusivities
141716H17/18	MMM7.1 with ExB flow shear effect for DRIBM	ITG/TEM diffusivities goes almost to zero with current levels ExB flow shear effects
141716H19/20	ExB flow shear factor is reduced from 1 to 0.5 for W19 and DRIBM	see comment below
141716H21/22	flow shear factor is increased from 0.5 to 0.75
141716H23/24	$k_yrho$ in DRIBM is increased from 0.2 to 2	see comments below
141716H25/26	flow shear factor is reduced from 0.75 to 0.5 for W19 and to 0.1 for DRIBM	see comments below

There are various definitions of $omega_{rm Etimes B}$ including the Hahm-Burrell and Waltz-Miller definitions. The flow shear rates computed using different definitions are often varied by a factor of 2 or more. The effect of ExB flow shear on particular modes in tokamak even is less certain and requires detailed studies. I plan to vary the ExB flow shear factor from 0.1 to 1 in my simulations depending on the scales of particular modes. In particular, I reduce the ExB flow shear factor to 0.1, when I change $k_yrho$ from 0.2 to 2, but I keet the ExB flow shear factor 0.5 for W19 that has $k_yrho=0.316$ in the PTRANSP simulations 141716H25 and 141716H26.

Some results from my PTRANSP simulations are below.
Case 141716H19/20

The effect of ExB flow shear is clearly visible if these results are compared with the results from the case 141761H11/12 that does not take into account the ExB flow shear effect for W19:

Without the ExB flow shear, the ITG/TEM modes become more unstable at later time. Account for the ExB flow shear makes these modes more stable at the later time.
Case 141716H23/24

Reduction of the characteristic scales of modes with $k_yrho$ increased from 0.2 to 2.0 resulted in a significant level of anomalous transport in the plasma core that has not predicted with the W19 or DRIBM model with the default $k_yrho$ parameters. The ExB flow shear factor should be probably less in these simulations because such short-scale modes should not be affected by the ExB flow shear at the same level as modes with $k_yrho=0.2$ . I reduced the ExB flow shear factor for these modes from 0.5 to 0.1 in the next case 141716H25/26.
Case 141716H25/26

Testing of the Weiland, ETG and DRIBM models for the NSTX discharge 141716

Several short PTRANSP simulations are performed to test the effect of the ExB flow shear on Weiland (ITG/TEM) and DRIBM diffusivities. The PTRANSP code has been run for a very short time to ensure that the plasma profiles do not change significantly from the experimental profiles. This choice of the PTRANSP parameters will allow direct comparison of the effective diffusivities, computed using the predictive transport models, with the diffusivities obtained in an interpretive TRANSP runs. [The TRANSP code does not support running predictive transport models when it runs in the interpretive mode.]

I used the following PTRANSP parameters in my simulation:

TINIT=0.482 !**** START TIME
FTIME=0.48201 !**** STOP TIME
DTMINT=4.e-6
DTMAXB=1.e-5
DTMAXT=1.e-5
DTINIT=1.e-5
TBONMIN=0.80 ! !**** MINIMUM BEAM TURN ON TIME

For the first set of simulations the start time was TINIT=0.364 and for the second set the start time was TINIT=0.482. For the both sets, the simulations run for 1e-5 sec. There were no sources in these simulations enabled. The objective was to estimate the predictive diffusivities for the experimental profiles and to validate some trends in the diffusivity profiles between the two time slices.

I found that the that the ETG diffusivities change very little between the two time slices. The result is somewhat unexpected because the plasma profiles are rather different:

	$Electron temperature profile$

While the electron density gradients change very little, there is a large increase of electron temperature gradient in the region from $rho=0.3$ to $rho=0.5$ . We assume that there is almost no flow shear effect in the ETG diffusivities and the ETG model does not include the ExB flow shear effect.

TRANSP modeling of NSTX discharge 141716

Two sets of simulations are proposed for the NSTX discharge 141716. The first set corresponds to the plasma parameters at t=0.364 sec. The second set corresponds to the plasma parameters at t=0.482 sec. The values of ExB flow shear was significantly increased between these two time slices in the regions from $rho=0.35$ to $rho=0.70$ and around $rho=0.80$ .

The objective is to test the applicability of the MMM 8.1 anomalous transport model under different values of the ExB flow shear. It is expected that the growth rates and effective diffusivities are much smaller at the later time, when the ExB flow shear rates are significantly larger. The NSTX experiment suggests that this statement is true for short scale instabilities as well.

Summary of anomalous diffusivity parameters that are used in the XGC0 modeling of four DIII-D discharges

DIII-D 132014

&tbl_param             ! (anomalous) turbulence diffusion
tbl_diffusion_on=1     ! 0: no turbulence diffusion, 1: simple turbulence diffusion
tbl_diffusion_mode=3   ! 1 = random walk model, 2 = advection-diffusion model
! mode=2 enables separate control of particle and heat diffusion
tbl_stop_time=500000   ! Stop the turbulence diffusion
tbl_period = 1       ! frequency of calling turbulence diffusion and fmcfm_call

tbl_d_profile_on=-1 ! 1: spatially varying D_turb profile (tanh), 0: spatially constant profile in pedestal region
tbl_slope_on=0 ! D_turb turned on only within finite poloidal angles (+45 deg) around outside midplane
tbl_d_profile_psicen=0.93D0 ! D_turb profile center
tbl_d_profile_psiwid=0.02D0 ! D_turb profile full width
tbl_d_profile_psicen2=0.999D0 ! D_turb profile center
tbl_d_profile_psiwid2=-0.01D0 ! D_turb profile full width

! particle diffusivities (random walk diffusion coefficients)
tbl_D_coeff_in=0.2D0     ! m^2/sec, inside (pedestal top) value of D_turb
tbl_D_coeff=0.03D0      ! m^2/sec, outside (pedestal bottom) value of D_turb
tbl_D_coeff2=0.6D0     ! m^2/sec, D_turb for outside of separatrix (SOL+private flux region)

! ion thermal diffusivities
tbl_therDi_coeff_in=.070D0    ! m^2/sec, inside (pedestal top) value of thermal D_turb for ions
tbl_therDi_coeff=-0.050D0    ! m^2/sec, outside (pedestal bottom) value of thermal D_turb for ions
tbl_therDi_coeff2=0.9D0    ! m^2/sec, thermal D_turb for outside of separatrix (SOL+private flux region) for ions

! electron thermal diffusivities
tbl_therDe_coeff_in=0.8D0    ! m^2/sec, inside (pedestal top) value of thermal D_turb for electrons
tbl_therDe_coeff=-0.1D0    ! m^2/sec, outside (pedestal bottom) value of thermal D_turb for electrons
tbl_therDe_coeff2=0.9D0    ! m^2/sec, thermal D_turb for outside of separatrix (SOL+private flux region) for electrons

DIII-D 132016

tbl_d_profile_on=-1 ! 1: spatially varying D_turb profile (tanh), 0: spatially constant profile in pedestal region
tbl_slope_on=0 ! D_turb turned on only within finite poloidal angles (+45 deg) around outside midplane
tbl_d_profile_psicen=0.94D0 ! D_turb profile center
tbl_d_profile_psiwid=0.02D0 ! D_turb profile full width
tbl_d_profile_psicen2=0.999D0 ! D_turb profile center
tbl_d_profile_psiwid2=-0.01D0 ! D_turb profile full width

! particle diffusivities (random walk diffusion coefficients)
tbl_D_coeff_in=0.13D0     ! m^2/sec, inside (pedestal top) value of D_turb
tbl_D_coeff=0.04D0      ! m^2/sec, outside (pedestal bottom) value of D_turb
tbl_D_coeff2=0.3D0     ! m^2/sec, D_turb for outside of separatrix (SOL+private flux region)

! ion thermal diffusivities
tbl_therDi_coeff_in=0.1D0    ! m^2/sec, inside (pedestal top) value of thermal D_turb for ions
tbl_therDi_coeff=-0.05D0    ! m^2/sec, outside (pedestal bottom) value of thermal D_turb for ions
tbl_therDi_coeff2=0.7D0    ! m^2/sec, thermal D_turb for outside of separatrix (SOL+private flux region) for ions

! electron thermal diffusivities
tbl_therDe_coeff_in=0.25D0    ! m^2/sec, inside (pedestal top) value of thermal D_turb for electrons
tbl_therDe_coeff=-0.05D0    ! m^2/sec, outside (pedestal bottom) value of thermal D_turb for electrons
tbl_therDe_coeff2=0.7D0    ! m^2/sec, thermal D_turb for outside of separatrix (SOL+private flux region) for electrons
/

DIII-D 132017

tbl_d_profile_on=-1 ! 1: spatially varying D_turb profile (tanh), 0: spatially constant profile in pedestal region
tbl_slope_on=0 ! D_turb turned on only within finite poloidal angles (+45 deg) around outside midplane
tbl_d_profile_psicen=0.91D0 ! D_turb profile center
tbl_d_profile_psiwid=0.02D0 ! D_turb profile full width
tbl_d_profile_psicen2=0.995D0 ! D_turb profile center
tbl_d_profile_psiwid2=-0.01D0 ! D_turb profile full width

! particle diffusivities (random walk diffusion coefficients)
tbl_D_coeff_in=0.13D0     ! m^2/sec, inside (pedestal top) value of D_turb
tbl_D_coeff=0.015D0      ! m^2/sec, outside (pedestal bottom) value of D_turb
tbl_D_coeff2=0.3D0     ! m^2/sec, D_turb for outside of separatrix (SOL+private flux region)

! ion thermal diffusivities
tbl_therDi_coeff_in=.1D0    ! m^2/sec, inside (pedestal top) value of thermal D_turb for ions
tbl_therDi_coeff=0.001D0    ! m^2/sec, outside (pedestal bottom) value of thermal D_turb for ions
tbl_therDi_coeff2=0.5D0    ! m^2/sec, thermal D_turb for outside of separatrix (SOL+private flux region) for ions

! electron thermal diffusivities
tbl_therDe_coeff_in=0.35D0    ! m^2/sec, inside (pedestal top) value of thermal D_turb for electrons
tbl_therDe_coeff=-0.1D0    ! m^2/sec, outside (pedestal bottom) value of thermal D_turb for electrons
tbl_therDe_coeff2=0.9D0    ! m^2/sec, thermal D_turb for outside of separatrix (SOL+private flux region) for electrons
/

DIII-D 132018

&tbl_param ! (anomalous) turbulence diffusion
tbl_diffusion_on=1 ! 0: no turbulence diffusion, 1: simple turbulence diffusion

tbl_diffusion_mode=3   ! 1 = random walk model, 2 = advection-diffusion model
! mode=2 enables separate control of particle and heat diffusion
tbl_stop_time=500000   ! Stop the turbulence diffusion
tbl_period = 1       ! frequency of calling turbulence diffusion and fmcfm_call

tbl_d_profile_on=-1 ! 1: spatially varying D_turb profile (tanh), 0: spatially constant profile in pedestal region
tbl_slope_on=0 ! D_turb turned on only within finite poloidal angles (+45 deg) around outside midplane
tbl_d_profile_psicen=0.89D0 ! D_turb profile center
tbl_d_profile_psiwid=0.02D0 ! D_turb profile full width
tbl_d_profile_psicen2=0.995D0 ! D_turb profile center
tbl_d_profile_psiwid2=-0.01D0 ! D_turb profile full width

! particle diffusivities (random walk diffusion coefficients)
tbl_D_coeff_in=0.05D0     ! m^2/sec, inside (pedestal top) value of D_turb
tbl_D_coeff=0.002D0      ! m^2/sec, outside (pedestal bottom) value of D_turb
tbl_D_coeff2=1D0     ! m^2/sec, D_turb for outside of separatrix (SOL+private flux region)

! ion thermal diffusivities
tbl_therDi_coeff_in=0.5D0    ! m^2/sec, inside (pedestal top) value of thermal D_turb for ions
tbl_therDi_coeff=0.001D0    ! m^2/sec, outside (pedestal bottom) value of thermal D_turb for ions
tbl_therDi_coeff2=2.0D0    ! m^2/sec, thermal D_turb for outside of separatrix (SOL+private flux region) for ions

! electron thermal diffusivities
tbl_therDe_coeff_in=.8D0    ! m^2/sec, inside (pedestal top) value of thermal D_turb for electrons
tbl_therDe_coeff=-0.07D0    ! m^2/sec, outside (pedestal bottom) value of thermal D_turb for electrons
tbl_therDe_coeff2=0.2D0    ! m^2/sec, thermal D_turb for outside of separatrix (SOL+private flux region) for electrons
/

Effect of neutrals is studied for the discharge DIII-D 132016

The effect of neutrals is studied for the discharge DIII-D 132016. The heat load width is relatively weekly affected by the neutral collisions.

Summary, the effect of neutrals are rather small comparing to the anomalous effects.

Summary plots of divertor heat flux widths for DIII-D and Alcator C-Mod discharges

These summary plots DIII-D and Alcator C-Mod discharges are based on Julian and my simulations. NSTX discharges are simulated by Gunyoung. I have to redo all simulations for the DIII-D discharge 132018 due to insufficient resolution in lim_store_mz_in/out (the resulution has increased to 500/1500. It looks like that lambda_div is consistently decreasing with the plasma current.

The DIII-D plot has been updated on August 27. Parameters associated with the electron dynamics modeling has been corrected.

Poloidal magnetic field scaling for DIII-D discharge 132017

$B_{pol}$ has been scaled with the factor that was varied from 0.3 to 1.5. The neutral collisions, anomalous diffusion, heating, and radiation losses have been disable for this scan. Also, the electron background profile has been fixed in these simulations (col_varying_bg_elec=0). The divertor heat flux width, $lambda_{div}$ , is shown as a function of the $beta_p$ multiplier, $C_{beta_p}$ on the figure below.

The radial electric field profiles for three different values of C_bp are shown below:

C_bp=0.75

C_bp=1.00

C_bp=1.50