The ExB flow shear rate has been implemented as simple derivative of the radial electric field with respect to minor radius:
New implementation is according to Ref. [T.S. Hahm and K.H. Burrell Phys. Plasmas 2, 1648 (1995)]:
where and are the poloidal and toroidal components of the magnetic field.
Smoothing procedure adopted from the ASTRA code on Friday (see previous post) is tested for set of 200 random numbers. The plots are generated for the smoothness parameter (left figure), (center figure), and (left figure):
In order to introduce more smoothing in the plasma core which is needed to compensate additional numerical noise in the core introduced by volume effect, dependence on is introduced:
where is the coefficient that controls region where smoothing is applied. The following function is selected for the next test:
where R is a random number in the range from 0 to 1, and are the coefficient that are set to 0.95 and 0.001 correspondingly. This function reproduce XGC-0 results for the radial electric field with more noise in the plasma core and large potential well in the plasma edge. The goal of smoothing is to remove noise in the plasma core and to preserve the details of potential well in the plasma edge. The results on the figure below are obtained for smoothness parameters and .
The level of smoothness in the plasma core is controlled by and the region where soothing is applied is controlled by the coefficient . Smaller values of should result in more extended region where the smoothing is applied. Test below show results for the same , but set to 2.
Two new cases are submitted on FRANKLIN in order to investigate different settings for smoothing.
- Case XGC-208 is based on the case XGC-206, except of fm_nsmooth_chi that is increased from 1 to 2.
- Case XGC-209 is submitted to test new procedure to smooth the radial electric field profiles. In addition to the parameter fm_nsmooth that is kept 1 (as in the previous cases), the following settings are used:
fm_nsmooth_er = 4 ! number of smoothing for radial electric field in the core refion
fm_smooth_er_rhob = 0.95 ! rho where transition between different level of smoothing is applied
fm_smooth_er_coef = 75. ! coefficient that controls sharpness of transition between different level of smoothness for Er
fm_nsmooth_chi = 2 ! number of applications of diffusivity smoothing in fmcfm_call
New procedure for smoothing of radial electric field profile is implemented in the fmcfm_call routine. Different levels of smoothing is applied in the plasma core and in the plasma edge. The radial electric field profile in the plasma edge is smoothed fm_nsmooth times together other plasma profiles such as electron and ion temperatures and plasma density profiles. The radial electric field in the plasma core is smoothed fm_nsmooth_er times. The new smoothed profiles are merge together using the following formula:
and used in the computation of the flow shear. Here, is the coefficient that describes sharpness of the transition between two smoothed profiles The fm_smooth_er_coef is used to set this coefficient in the XGC-0 code through the namelist. The parameter controls the location of the transition. The fm_smooth_er_rhob variable is used to set this parameter in the code. The default values for these parameters are: and .