1 / 9

Update on FIXFREE and CREATE-NL ITM General Meeting Lisbon, 13-15 September 2010

This update provides information on the achievements and open problems related to FIXFREE and CREATE-NL. It includes the modification of CREATE-NL to accept nonzero current density on the plasma boundary, as well as the adaptation of CREATE-NL to ITM needs. The update also discusses the problems with running CREATE-NL under Kepler and the coupling of CREATE-NL with the Transport CNL equilibrium code.

lorenzon
Download Presentation

Update on FIXFREE and CREATE-NL ITM General Meeting Lisbon, 13-15 September 2010

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. M. Mattei, G. Calabrò, E.Giovannozzi, C. Labate Update on FIXFREE and CREATE-NLITM General MeetingLisbon, 13-15 September 2010

  2. Achievements CREATE-NL has been modified so as to accept nonzero current density on the plasma boundary Open problems High current density gradient in profiles cannot be accepted (current density derivatives have to be compatible with the FEM mesh elements dimensions). Increasing number of elements may be a solution (high computational burden) completed (a more robust convergence) Smoothing can be a solution (approximated) completed (a more robust convergence) CREATE-NL adaptation to ITM needs

  3. Achievements CREATE-NL has been implemented on Gateway CREATE-NL can read and write CPOs (ITER test case) Open problems CREATE-NL cannot run under Kepler for two main reasons Matlab licenses Kepler-Matlab interfaces CREATE-NL as a Kepler Actor

  4. Open problems on Kepler (UAL version 4.07c) – Matlab: CREATE-NL as a Kepler Actor How to “read & use”? More help informations (ITM and ISIP) needed.

  5. Problems to be discussed in view of coupling with Transport CNL equilibrium code accepts plasma current density profiles in terms of ff’ and p’ functions of normalized psi; Transport codes need an evolution of q profile which is consistent with transport equations; As free boundary equilibrium code computes step by step the psi map, q profile is slightly modified at each step, causing mismatched and artificial transients in the transport code; Non-physical drift in q profile is induced; Intermediate iterations at each time step required to solve this problem. Simulations need a robust plasma shape control able to counteract variation in li and poloidal beta. Tools for designing such a controller are needed.

  6. Experience in coupling CREATE-NL to JETTO • The coupling between the codes occurs with a time step depending on the shape variation dynamics. Actually, to simulate ITER cases the time step is 0.1sec (without simulating eddy currents and vertical stability). JETTO works with a time step 1msec. Plasma current is assigned. • At each step (2 distinct problems): • JETTO “p” and “f” profiles to CREATE-NL • JETTO “f” profile is slightly modified to guarantee the “q” continuity • Shape (and eventually VS) control “corrects” the external currents to preserve the equilibrium • CREATE-NL equilibrium to JETTO

  7. FIXFREE (in poster session) Installed on the Gateway, on a private SVN archive. A makefile can build the full program and the actors too. As an independent program or as Kepler Actor. It needs in input the current in the coils and a profile of p’ and ff’ as function of psi. A few fixes are still needed in order to be compatible with ITM (i.e. the vectors are in reverse order from the outside of the plasma to the inside, etc.)

  8. FIXFREE example: FAST equilibrium • FAST input tree: • /afs/efda-itm.eu/imp5/user/egiovan/public/itmdb/itm_trees/test/4.07b/mdsplus/0 • euitm_30002.characteristics • euitm_30002.datafile • euitm_30002.tree • The output produced by FIXFREE is in: • euitm_30005.characteristics • euitm_30005.datafile • euitm_30005.tree • The kepler model is: • /afs/efda-itm.eu/imp5/user/egiovan/FIXFREE/src/modelloKepleroM.xml • Provisional manual how to use FIXFREE under Kepler is avalaible.

  9. FIXFREE example: FAST equilibrium

More Related