************************************************************************************************************** ************************************** MAGNETIC RIPPLE PARAMETERS ******************************************** ************************************************************************************************************** Create the local file RIPPLE_id_equil.mat by Y.Peysson CEA-DRFC <yves.peysson@cea.fr> and Joan Decker MIT-RLE (jodecker@mit.edu)
0001 %************************************************************************************************************** 0002 %************************************** MAGNETIC RIPPLE PARAMETERS ******************************************** 0003 %************************************************************************************************************** 0004 % Create the local file RIPPLE_id_equil.mat 0005 % 0006 %by Y.Peysson CEA-DRFC <yves.peysson@cea.fr> and Joan Decker MIT-RLE (jodecker@mit.edu) 0007 % 0008 clear all 0009 % 0010 id_equil = 'JETliketest'; 0011 id_ripple = 'rippletest'; 0012 % 0013 ripple.id = id_ripple; 0014 ripple.equil_id = id_equil; 0015 % 0016 %For the magnetic ripple losses, the loss frequency must be very large compared to the reference frequency. Therefore it is always enforced in the code to 1000*nhu_ref 0017 %Consequently, the ripple parameter is just a flag which indicates if magnetic ripple calculations are considered or not. Calculations are only valid for Tore Supra tokamak 0018 %so the first two letters must be TS to be recognized in the calculations. See main_dke_yp.m 0019 % 0020 ripple.mode = 1;%No magnetic ripple if empty 0021 % 0022 ripple.N = 18;%Number of toroidal magnetic field coils 0023 % 0024 % For the Yushmanov model in Reviews of Plasma Physics, p.117-241, "Diffusive Transport Processes Caused by Ripple in Tokamaks" by P.N. Yushmanov, 1990. 0025 % 0026 ripple.B0 = 3.2;%Vacuum toroidal magnetic field on the toroidal axis, (signed value as BPHI) (T) 0027 ripple.R0 = 3.05;%Toroidal major radius of the toroidal magnetic field windings(m) 0028 ripple.rn = 1.0;%Inner radius of the toroidal magnetic field windings (m) 0029 ripple.rx = 1.5;%Outer radius of the toroidal magnetic field windings (m) 0030 % 0031 % 0032 %************************************************************************************************************* 0033 % 0034 eval(['save RIPPLE_',id_equil,'_',id_ripple,'.mat ripple']);