Tokamak: JET

    Pulse number: 46664

    Contact person: Pam Stubberfield, Vassili Parail

    Institution: EFDA-JET

    Date of shot: 981216

    Analysis codes: TRANSP, EFIT for plasma boundary and initial q profile

    Run number: 46664C03

    Date of analysis: 990709 TRANSP run by Yuri Baranov

    Assumptions made for analysis:

    Ti profile from CX measurements is used from 45.00 to 46.35 secs, before
    this time period a Ti profile is provided by assuming Xi=Xneo.
    Te profile vs. radius is from KK3 measurements. 
    Ne profile from LIDAR measurements.1D bulk radiation is used in the analysis.
    Rotation velocity profile from CX measurements is provided.
    1D Zeff from visible bremsstrahlung to 45.00 secs then 2D Zeff profile
    B impurity. Edge ptcl confinement time is set to 0.3s.
    q profile provided by EFIT used from start of run at 42.0 secs
    until 43.0 secs. TRANSP field diffusion then provides the q profile.
    TRANSP run starts at 42.0 secs so no LHCD modelled.
 
    Shot description,purpose of the experiment: Optimised shear discharge

    Publication: Tala T.J.J et al, Plasma Phys. Contr. Fusion 43(2001) 507
    Additional information:

    Pulse #46664 is one of the best Optimised Shear JET discharges. 
    The current rises to Ip=3.2MA at a toroidal magnetic field Bt=3.45T. 
    Its characteristic feature is that plasma breakdown was assisted by
    LHCD with the power Plh=1.5MW applied from t=40.3sec until t=41.2 sec.
    The plasma current was ramping up throughout the pulse so that it has
    reached Ip=3.2MA at t=46.5 sec (shortly before disruption). Current
    ramp up has been accompanied by 1MW ICRH preheating to keep the
    q-profile flat or slightly reversed. The main heating power was
    applied at t=44.7 sec and reached its full capacity (16MW of NBI
    and 6MW ICRH) at t=45 sec. Plasma energy was rising throughout the
    main heating phase, reaching the level of 10MJ and confinement time 
    enhancement factor H89 exceeded 2.2 shortly before disruption. 
    The plasma developed a strong ITB during the main heating, which
    gradually expanded its footpoint to R=3.55m. 

    There are 24 1-D Ufiles and 37 2-D Ufiles.  All beam related quantities
    (energy, particle sources, currents etc) have been smoothed to
    suppress Monte-Carlo noise. In the 2d file the deuterium density is
    stored in NM1. 
    In the file JET_46664_0d.dat data is given at 45.10 secs. 

    Values for kappa and delta in the 0d, 1d and 2d files are from the TRANSP
    output derived from the boundary input definition provided by EFIT.
    kappa = (Zmax-Zmin)/(Rmax-Rmin), delta = (0.5(R(Zmax)+R(Zmin))-Rgeo)/a
    XPLIM is from XLOC ; SEPLIM is given by EFIT.
    PLTH is defined by WTH/TAUTH coming directly from TRANSP output and not 
    as defined in the standard list of variables.  
  
    The TRANSP run gives a simulated W_dia and neutron yield about
    10% lower than the measured values.