Core-Edge Coupled Transport Simulations

The goal of this work is to implement a direct scheme for coupled core-edge transport simulations. For the core the ETS, the core impurity code and the core neutrals code will be used. For the edge SOLPS-B2 will be used.

At the end of 2011, a manual coupling procedure had been implemented, as described in this document.

Activities 2012

TopPreviousNext

In 2012 work is ongoing to automate this procedure. The following steps have occurred:

  • Modified the main part of SOLPS-B2 to introduce

    • "b2mn_init", a routine to initialize SOLPS-B2
    • "b2mn_step", a routine that can be called repeatedly to advance the edge plasma simulation
    • "b2mn_fin", a routine to finalize SOLPS-B2

    These changes involved splitting "b2mndr" into "b2mndr_0", "b2mndr_1" and "b2mndr_2"

  • The introduction of a new routine "b2mn_ets" which is called with input CPOs "coreprof", "coreimpur" and "coreneutrals", and returns output CPOs "coreprof", "coreimpur", "coreneutrals" and "edge". The routine:

    • calls "b2mn_init" on the first call
    • based on the input core CPOs, alters the SOLPS-B2 boundary conditions before calling "b2mn_step", and then revises the boundary conditions in the output core CPOs
    • calls "b2mn_fin" at the end of the calculation

  • Changes to many of the input/output routines in SOLPS-B2 so that

    • input data is looked for in the present directory, it's parent, or the sibling "baserun" directory
    • output files that are already present no longer cause the code to stop with an error message

    The result of this is that SOLPS-B2 no longer needs to be run by the previously required Makefile which copied input files into a run directory.

  • The modification of the ETS Fortran workflow code "eq_ets_test" to include calls to "b2mn_ets".
  • The Fortran workflow runs coupled calculations, and verification runs are currently in progress.

    • A problem in the UAL which causes the message "No more slots available for arrays of structures" and the program then to hang has been identified
    • Fixes in the 4.09a ETS were back-ported from 4.10a, and problems with impurity transport coefficients identified and fixed. There may still be a problem with impurity energy losses.
    • D and D+C+Ar+Ne+He runs are both being done


last update: 2013-06-18 by denka