List of Completed DIVIMP Development Tasks

The following entries describe development items that have been added to the DIVIMP code since the Development Task list was first adopted in the fall of 1999. Items are considered to be substantially complete by inclusion in this list, any portions of an item that may still need to be addressed are marked in red.

  1. DIVIMP code fixes
  2. Integrate the Reiser transport model from Peter Schwanke's code
  3. Integrate the code changes made at DIIID into the base code version
  4. Integrate Steve's latest changes and options into the base code version
  5. Web based case log and results access
  6. Documentation update
  7. Code development items related to getting the current version of DIVIMP/Eirene to work for Asdex Upgrade - Karl's visit from Feb 21st to March 3rd 2000
  8. Resolve issue of optional input items in the input file and decide on policy for this feature for future code development.
  9. Filterscope geometry details
  10. Add hydrocarbon break-up option.
  11. Continue re-organizing the output files splitting the impurity simulation output from the OSM output.
  12. Generate bolometry plots
  13. DIVIMP cases
  14. EIRENE clean up and completion of code integration
  15. Extract fluxes from EIRENE.
    The fluxes and average energies of D-atoms and D-molecules onto surfaces near/at the pressure gauges and pump openings need to be extracted so we can do hand checks with the p-gauge simulations directly produced by EIRENE.
  16. Change grid near pump opening so that the 'plasma tongue' goes into the pump opening. (This has been done, but the grid extensions have been hardcoded into DIVIMP for each of the DIII-D grids. That is, DIVIMP is not presently able to analyse the DIII-D grids and then see where the field lines need to be extended. So, for the general user, this feature is not really available.) (Steve)
  17. Turn on EIRENE n-n collisions in the voids.
  18. On Along Ring plots, use different symbols for the DTS points that are for the time closest to the time of interest, e.g. 1650 ms currently (the DTS points for other times involve greater uncertainties related to EFIT since EFIT has to be used twice for them.)
  19. Overlay SOL22 and SOL23 solutions in Along Ring plots to make direct comparisons.
  20. Utilize EIRENE puffing capabilities
    Something new: find out how to introduce extra puffing of deuterium into the PFZ (which will have to be pumped away at same rate as injected) as a way to fill in the Dalpha-valley. It's known that the walls are net sources of neutrals for certain periods in a discharge (also net sinks at other times) and it's been suggested that the PFZ wall may be a net source during the flat-tops of shots, perhaps being re-filled in ramp-ups and ramp-downs.
  21. Add the Psin values for each ring to the print-out.
  22. Report for Gary, Jon, etc updating the current status of DIVIMP modeling for shot 86575.
  23. Additional core plasma options based on experimental input.

    NOTE: The code already has core plasma condition options that allow the input data for each ring in the core to be specified and further allow this data to be used as the basis for a MARFE type core plasma. Is more development needed? NO

  24. Turn SOL23 into user-friendly, work-horse, main-line OSM.
  25. TN1517 supplement to Chiperp-Extractor:

    I'd like to add a new output related to the geometry of the separatrix flux surface, namely its effective area for cross-field diffusive transport, A_sep_eff. For both particle and power, the LOCAL fluxes across the separatrix surface involve delta_A_perp/lambda, where delta_A_perp is the element of surface area of the separatrix surface and lambda is the density or temperature decay length at the separatrix at that element. Now due to flux expansion, the local lambda will equal the reference lambda, lambda_0 , times the flux expansion ratio, Del/Del_0, where Del_0 is the cross-field separation at the reference location between the separatrix flux tube and the next flux tube outward, and Del is the same but for the particular surface element. Thus the effective surface area of the separatrix surface for diffusive transport can be defined to be the sum of delta_A_perp / (Del/Del_0), summed over all the elements of separatrix surface. For the cylindrical approximation of a tokamak one can just use the LENGTH of each segment of the separatrix flux line in the POLOIDAL PLANE and multiply by 2*pi*R_o to get delta_A_perp. For the toroidal geometry version one uses the local value of R. The reference location should be the outside midplane.

  26. Explore how 2D plots (to compare with Max's toroidal view camera) and filterscope profiles change as one changes the target mult factors, in order to see if this is going to be a way that we can confirm the target probe data.
  27. Compare the 4 cases of Gary's new EFITs done in the 2 different ways, and the 2 different ways of assigning the DTS points.
  28. Add to the standard print-outs the (R,Z) values where each Ring strikes the targets.
  29. To compare Max's toroidal view camera of CIII (and for new shots, of H-lines also), we'll have the basic comparison based on simply displaying the 2D plots from Max and from the code, side by side. However, as I mentioned, that's hard to turn into a quantitative comparison, and Gary has produced 'Along-R Slice Plots', at different values of Z. We can use his expt 'Slice Plots' and compare our code ones. (You've called these 'arbitrary slice plots').
  30. DIVIMP DEVELOPMENT ITEM: power deposition profile
    We have not been interested in producing plots of the deposited heat flux [W/m^2] because we're pretty sure we won't get very good agreement with the infrared, IR, camera results. That attitiude, however, is as bad as Gary just choosing to ignore the probe results. We'd better face up to the discrepancy with the IR results and be up front about it. Trying to sweep it under the carpet is, in the end, more discrediting than just owning up to it, and saying it's one of the things we don't match and don't know why. Fortunately, it's a long-standing problem on DIIID. About 10 years ago they first established the fact that they couldn't confirm the theoretical values of about 7 for the sheath heat transmission coeff, shtc: they combine the IR heat flux density with the probe Isat and Te, and spit out a shtc. They've refined that procedure over the years, but they still get values of about 1 for the shtc near the strike point. Further out the value rises to close to the theoretical value of about 7. A long-standing, stubborn mystery. Anyway, we better be up front about it. But to do that we have to produce code plots of deposited heat flux density to overlay on, or just put side by side with, the expt, IR one, e.g. Fig 2-4(a) in Gary's Memo. The code value of q_dep[W/m^2] consists of a number of parts:
    David Elder
    Last modified: Thu Aug 24 12:11:07 EDT 2000