1 / 8

Strategy

ITER divertor and SOL have been modelled with SOLPS4.2 (A. Kukushkin, H. Pacher, G. Pacher, et al.) Large library of results Contains extended physics package for neutrals and pumping, including Eirene “face-lift” But no drifts, no currents, and little current physics development

borna
Download Presentation

Strategy

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. ITER divertor and SOL have been modelled with SOLPS4.2 (A. Kukushkin, H. Pacher, G. Pacher, et al.) Large library of results Contains extended physics package for neutrals and pumping, including Eirene “face-lift” But no drifts, no currents, and little current physics development The more recent SOLPS5 package includes the latest physics and is the currently supported version Treatment of X-point area much improved by use of indirect addressing SOLPS5.1 = SOLPS5.0 + Eirene “face-lift” So we are currently attempting to migrate the ITER library of SOLPS runs from 4.2 to 5.1 Paying attention to recover old results With benefits of new physics to advance in the design Documenting the changes from one version to another Transfer of ITER SOLPS4.2 simulations to SOLPS5.1X. Bonnin (CNRS-LIMHP), A. Kukushkin (ITER), D. Coster (IPP-Garching)

  2. Strategy • Continue a SOLPS4.2 run with SOLPS5.0 • Wrote a converter to take SOLPS4 input/output files into SOLPS5 format • Single time-step. Re-evaluate fluxes and examine output • Added new diagnsotics (radiation, bremsstrahlung) to SOLPS5.0 • Implemented the detailed particle balance treatment necessary to resolve the ITER pumping fluxes • Identified several differences/mistakes • Example: Luciani correction on closed field lines: • SOLPS4: D → D/(1 + 3e15 * T^2/(ne * Z^2)) (1) • SOLPS5: D → D/(1 + (3.75e16/L) * <T^2/(ne * Z^2)> ) (2) • Implemented a series of switches for „backward compatibility“ of SOLPS5

  3. Some of the identified issues

  4. Ion temperature equations SOLPS4: total energy SOLPS5: internal energy

  5. Flux comparisons (I): past the X-Point # fnisip fnisap feisip feisap feesip feesap 4.0_short -1.07243E+23 -5.33376E+22 -1.54480E+07 1.85345E+07 -2.42950E+07 3.50906E+07 5.0_with_4.0_short -1.13724E+23 -5.89033E+22 -9.21913E+06 -2.12669E+06 -4.90001E+07 6.80004E+07 5.0_with_4.0_area_fix0 -1.55764E+23 -3.39944E+22 -1.50248E+07 1.38652E+07 -5.53261E+07 7.99219E+07 5.0_with_4.0_luciani2 -1.55765E+23 -3.39939E+22 -1.50178E+07 1.38580E+07 -5.53236E+07 7.99277E+07 5.0_with_4.0_b2tqce -1.55780E+23 -3.39781E+22 -1.50314E+07 1.38742E+07 -4.80632E+07 7.39885E+07 5.0_with_4.0_b2tqca -1.55514E+23 -3.36086E+22 -1.40038E+07 1.23717E+07 -4.80571E+07 7.39860E+07 5.0_with_4.0_mode2 -1.55512E+23 -3.38129E+22 -1.61660E+07 1.85562E+07 -5.08657E+07 7.78907E+07 5.0_with_4.0_poteq2 -1.56334E+23 -3.31955E+22 -1.62630E+07 1.87585E+07 -5.14561E+07 7.88552E+07 5.0_with_4.0_22 -1.56308E+23 -3.30849E+22 -1.67913E+07 1.82750E+07 -5.08636E+07 7.62398E+07 5.0_with_4.0_23 -1.56273E+23 -3.32557E+22 -1.67585E+07 1.81844E+07 -2.59605E+07 3.65067E+07 5.0_with_4.0_24 -1.05084E+23 -4.92526E+22 -1.48725E+07 1.83610E+07 -2.40191E+07 3.55561E+07 5.0_with_4.0_25 -1.05084E+23 -4.92538E+22 -1.48775E+07 1.83843E+07 -2.40186E+07 3.55560E+07 5.0_with_4.0_no_solve -1.04748E+23 -4.92469E+22 -1.46850E+07 1.82429E+07 -2.36311E+07 3.53916E+07 5.0_with_4.0_solve_ua -1.05067E+23 -4.93021E+22 -1.47024E+07 1.82311E+07 -2.36289E+07 3.53778E+07 5.0_with_4.0_solve_co -1.04684E+23 -4.92552E+22 -1.46836E+07 1.82456E+07 -2.36302E+07 3.53885E+07 5.0_with_4.0_solve_ht -1.04748E+23 -4.92469E+22 -1.49470E+07 1.84051E+07 -2.41166E+07 3.55713E+07

  6. Flux comparisons (II): at the targets # fnixip fnixap feixip feixap feexip feexap 4.0_short -1.72205E+24 2.34399E+24 -2.16323E+06 7.64054E+06 -1.90619E+06 6.63314E+06 5.0_with_4.0_short -1.73176E+24 2.35863E+24 -1.71247E+06 6.65389E+06 -2.57813E+06 9.10311E+06 5.0_with_4.0_area_fix0 -3.60897E+24 6.67273E+24 -4.01878E+06 1.77044E+07 -3.99014E+06 1.41969E+07 5.0_with_4.0_luciani2 -3.60850E+24 6.67226E+24 -4.01433E+06 1.76911E+07 -3.98900E+06 1.41951E+07 5.0_with_4.0_b2tqce -3.58877E+24 6.64197E+24 -3.92338E+06 1.74966E+07 -3.71052E+06 1.34651E+07 5.0_with_4.0_b2tqca -3.58739E+24 6.63811E+24 -3.90465E+06 1.71188E+07 -3.70489E+06 1.34461E+07 5.0_with_4.0_mode2 -3.58745E+24 6.63466E+24 -3.90516E+06 1.70399E+07 -3.70624E+06 1.33874E+07 5.0_with_4.0_poteq2 -3.12039E+24 5.78374E+24 -3.59958E+06 1.64228E+07 -4.72017E+06 2.07636E+07 5.0_with_4.0_22 -3.09866E+24 5.76404E+24 -3.54744E+06 1.63860E+07 -4.58989E+06 2.04137E+07 5.0_with_4.0_23 -3.05312E+24 5.63017E+24 -3.33489E+06 1.55668E+07 -3.26749E+06 1.25012E+07 5.0_with_4.0_24 -1.10058E+24 1.34661E+24 -1.43528E+06 1.10517E+07 -1.72900E+06 6.18843E+06 5.0_with_4.0_25 -1.10043E+24 1.34665E+24 -1.43545E+06 1.10715E+07 -1.72847E+06 6.18636E+06 5.0_with_4.0_no_solve -1.72724E+24 2.35326E+24 -1.73297E+06 1.14022E+07 -1.91073E+06 6.65224E+06 5.0_with_4.0_solve_ua -1.73165E+24 2.35730E+24 -1.73532E+06 1.13927E+07 -1.91307E+06 6.64763E+06 5.0_with_4.0_solve_co -1.73789E+24 2.38036E+24 -1.74732E+06 1.15849E+07 -1.92438E+06 6.79571E+06 5.0_with_4.0_solve_ht -1.72724E+24 2.35326E+24 -1.90109E+06 1.13541E+07 -2.04373E+06 6.46947E+06 5.0_with_4.0_solve_uaco -1.74189E+24 2.37633E+24 -1.74749E+06 1.15346E+07 -1.92492E+06 6.73133E+06 5.0_with_4.0_solve_uaht -1.05917E+24 1.26341E+24 -1.46269E+06 1.13199E+07 -1.74057E+06 5.95530E+06 5.0_with_4.0_solve_coht -1.73815E+24 2.38033E+24 -1.90631E+06 1.14632E+07 -2.03978E+06 6.53116E+06

  7. Remaining issues • Ion temperature equation • Most terms in code can now flip between the two formulation according to a switch • Missing the viscous and kinetic energy terms in the solving step for the equation itself • SOLPS5.0 (Eirene 96 & 99) instead of SOLPS5.1 (Eirene „face-lift“) • Work currently underway to finalize the coupling (V. Kotov) • Not all differences between the two codes have been properly explained yet • Must complete the work with series of automated translation/resubmission/analysis scripts

  8. traces.1: Corrected src_ allocation and Zeff averaging bounds traces.2: As above, without interpolation of velocities in converter traces.3: As 1, with 'eirene_extrap'.eq.0, correction of sources entry argument list, new computation and scaling of residuals, fixed impurity radiation traces traces.4: With dtim=1e-29, as 3, with src00 arrays filled, src_r recoded, residuals and SSEEtot and SSEItot corrected traces.5: As 4, coded the volume calculation from SOLPS4.0 (but not used!) traces.6: As 5, with new (misread) bremsstrahlung treatment traces.7: Using new grid with SOLPS4.0 volumes, bremsstrahlung rate read properly, proper zeroing of 3rd instant of sources traces.8: Corrected bremsstrahlung traces, dt=1.e-9s. traces.9: Sum of absolute values of residuals. traces.10: Changed UA from UP to UU/PIT. traces.11: Corrected SX factor in the fluxes traces 12: Correction 'eirene_fixmap' in fort.31 at cuts traces.13: Correction 'eirene_fixuub'.eq.'2' in fort.31, upb=uub*rrb traces.14: Removed ion flux limiter traces.15: Running with same STRAHL data as SOLPS4.0 traces.16: Removed tailep and fix_cx settings from b2ar. New b2frates. traces.17: SOLPS4.0 Luciani formulation introduced: b2trcl_lluciani.eq.2 traces.18: Introduced SOLPS4.0 formulation of electron classical parallel transport traces.19: Added SOLPS4.0 formulation of ion classical parallel transport traces.20: Run with mode=2, no_solve.eq.1, area_fix=3, no hybrid schemes traces.21: Same as 20, but with more Eirene histories for better statistics and pot=3.1*te/qe traces.22: Same as 21, but with lesser (SOLPS4.0) constant accuracy and corrected no_current treatment, SOLPS4.0 upwind scheme for energy (i&e), SOLPS4.0 treatment for electron flux limits traces.23: Same as 22, but fixed fch_[3|5]2 no_current treatment and with SOLPS4.0-like treatment for particle fluxes computation traces.24: Same as 23, after clean-up of uses of vol,hx,hy,gs (only 1 Eirene call) traces.25: Code as 24, with ion flux limits off, and solving

More Related