Simulation with strongly reduced particle anomlous transport show expected slow evolution of plasma density profiles in the pedestal region (case XGC-220). There is still significant increase of the electron temperature in the SOL region.
Summary: There might be different ExB flow shear effects for the particle and thermal transport
Case XGC-219 has the same input as case XGC-216 but uses an improved expression for the ExB flow shear rate.
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 
The FMCFM interface is modified and particle diffusivity is updated for this case with tbl_d_coeff_in;:
fm_ion_diff_norm%particle(i) = tbl_d_coeff_in * fm_norm_diff
The particle diffusivity in the core region is selected to be 
The fm_wexb coefficient has been increased from 1e-6 to 1e-4. This factor determines the effect of ExB flow shear suppression on the anomalous transport:
There is no dramatical particle losses as in the case of MMM95 model implemented using the FMCFM interface.
The time step in the case XGC-213 has been increased from 2e-3 to 2e-2. The plasma profiles looks more noisy.
Summary: I should try to increase time step by a factor of 2 rather than 10.
Simulation results for the case XGC-209 that uses smoothing procedure introduced on January 12. In addition, effect of fm_wexb parameter on the particle transport is studied. The results show noticeable smoothing for the radial electric field profiles in the plasma core and very little change in the plasma density dynamics.
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 
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, 
where 
where R is a random number in the range from 0 to 1, 
The level of smoothness in the plasma core is controlled by 
The idea for smoothing of plasma profiles is based on routine implemented in the ASTRA code. The resulting function is based on finding of the minimum of the functional
with respect to U. The functional is design for the cylindrical case and the resulting function has zero derivative for x=0 and preserve boundary condition at x=1. The description of the method can be found at:
