80 likes | 239 Views
Equilibrium and Profile Advance (EPA) Status & Future Plans. L. P. Ku, S. Jardin, D. McCune Princeton Plasma Physics Laboratory SWIM Project Meeting Oct. 15-17, 2007 Oak Ridge National Laboratory. Outline. Status EPA in IPS EPA Functionality
E N D
Equilibrium and Profile Advance (EPA)Status & Future Plans L. P. Ku, S. Jardin, D. McCune Princeton Plasma Physics Laboratory SWIM Project Meeting Oct. 15-17, 2007 Oak Ridge National Laboratory
Outline • Status • EPA in IPS • EPA Functionality • Built-in flexibility for various modes of operation • TSC modernization/upgrade • Applications • ITER runaways, ITER hybrid discharge, CMOD current ramp-up • Future Plans
Status – EPA in IPSEPA has been fully integrated into the IPS framework. • Sample driver scripts and simulation configuration files written, tested and posted. • epa_driver.py, epa_lmhd_driver.py • epa_tsc.py, epa_tsc_cmod.py, epa_tsc_iter.py • epa_sim.conf, epa_iter_sim.conf, epa_iter+lmhd_sim.conf … • Available in normal mode of IPS operation or in stand-alone mode to generate libraries of plasma states for IPS to run in re-play mode. • Useful for component development, testing, comparative studies for resolution requirement and porting other components to new computer platforms
Status – EPA in IPS (cont) • In addition to updating the plasma state, EPA writes • GEQDSK file (.geq) • for the plasma state software to reconstruct equilibrium and write to the plasma state • for use by the RF and Fokker Planck components • JSOLVER file (.jso) • for use by the LMHD component to refine equilibrium in stability analyses.
EPA – FunctionalityBuilt-in options in EPA (TSC) increase the versatility and flexibility • The present EPA uses TSC that • performs free-boundary self-consistent transport evolution in any part of a tokamak discharge. • models PF coils and 2-D passive structures with circuit equations. Arbitrary power supply models can be used. • Runs with feedback controls: e.g. radial, vertical positions, shape, plasma current, stored energy. • EPA (TSC) runs in a variety combination of modes for the heating source and current drive: • normal IPS mode • internal coupling to TRANSP • internal analytic descriptions • data file input • automatic switching to internal analytic mode if other modes fail • new option recently added for CQL3D coupling.
EPA-Functionality (cont) • EPA(TSC) has been upgraded and improved: • Conversion to F90 recently completed. Code structure modularized to increase the flexibility for future modifications. TSC is now under SVN source control. • Transport model – GLF23 hardened to alleviate problems of numerical instability, hence greatly increased the reliability of results.
EPA – ApplicationSuccessful coupling to Fokker-Planck, linear MHD and RF (TORIC, AORSA) components accomplished. • ITER runaways • Coupling of EPA and CQL-3D in IPS normal mode of operation for the first 10 s of a discharge (without correcting the electric fields due to the runaways) indicated the runaways may be significant. • a tighter coupling of CQL-3D to EPA is called for, either internal or external to EPA(TSC), for next round of study. • ITER burn • Simulation of a hybrid discharge for 250 s, using internal TRANSP coupling with 25 ms steps has been done to generate a sequence of plasma states. • These states are retrievable for re-play and have been used by RF or LMHD components for various testing. • CMOD current ramp-up • Simulation of a RF heated discharge from 0.75 s to 2 s has been developed. • Cases for both normal IPS mode and re-play mode are available.
EPA – Future Plans • EPA in IPS • Porting to Jaguar • Checking consistency with other drivers in the integrated opearation • Coupling to AORSA, TORIC • Coupling to LMHD • Coupling to CQL3D (GENRAY) • Implement interface with machine description files • EPA functionality • Add prediction of plasma rotation profile • Include fast ion pressure in equilibrium evolution • Improve density prediction • Implement new GLF23 (TGLF23, available in Dec) and standardized transport interface • Implement options for edge pedestal model (including ELMs)