1 / 9

Enhancements To Full-Wave TORIC for HHFW Modeling

Enhancements To Full-Wave TORIC for HHFW Modeling. P. Bonoli, PSFC, Cambridge, MA M. Brambilla, IPP, Garching, Germany C.K. Phillips, PPPL, Princeton, NJ NSTX Research Forum HHFW / EBW Planning Session September 11-13, 2002. Applicability of TORIC to HHFW Regime.

kalkin
Download Presentation

Enhancements To Full-Wave TORIC for HHFW Modeling

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. Enhancements To Full-Wave TORIC for HHFW Modeling P. Bonoli, PSFC, Cambridge, MA M. Brambilla, IPP, Garching, Germany C.K. Phillips, PPPL, Princeton, NJ NSTX Research Forum HHFW / EBW Planning Session September 11-13, 2002

  2. Applicability of TORIC to HHFW Regime • TORIC uses SCK wave equation • Valid for (ki)2 < 1 • But in NSTX we can have (ki)2 ≈ 50 • k≈  / VA, Ti = 0.5 keV, ne = 4  1019 m-3, Bt = 0.25 T • Using the ion FLR wave equation in this regime gets wrong wave polarization and subsequently the wrong wave damping • ELD  |Ez|2 • TTMP  |Ey|2 • Cross - term  Ey Ez

  3. TORIC Modifications for HHFW Regime • Reformulated dielectric tensor elements in coefficients of SCK equations using full Bessel function expansion [Brambilla, 2002]: • Approach used by Belgian group in the SPRUCE code • Can then evaluate cyclotron damping on minority hydrogen, following Ono [Physics of Plasmas 2, 4075 (1995)]

  4. Applicability of TORIC to HHFW Regime - Tests Performed • Compare the following: • TORIC with ion FLR wave equation • METS 1D with full ion Bessel expansion (Phillips) • Ray tracing with full ion Bessel expansion (Rosenberg) • TORIC with modified full-wave equation

  5. Ti (0) = 0.5 keV Te (0) = 1.5 keV ne (0) = 4.1  1019 m-3 Bt = 0.32 T Ip = 800 kA Profiles: ne  [1 - (r/a)4 ] Te  [1 - (r/a)4 ] Ti  [1 - (r/a)2 ] f0 = 30 MHz n = 22 k(ANT) ≈ 14 m-1 Parameters for Comparison of Models • Plasma Composition • 90% Deuterium • 10% Hydrogen

  6. TORIC Results with ion FLR Wave Equation Results in substantial disagreement with 1-D integral wave code METS and toroidal ray tracing (HPRT)

  7. P (1/m3) 0.8 0.8 P (1/m3) 0.4 0.4 0 0 1.2 1.4 1.6 1.0 Major Radius (m) Models predict off-axis wave absorption on electrons 2D effect: k// evolution due to Bpol HPRT ray tracing code METS 1D full-wave code

  8. TORIC Results with Modified Wave Equation Similar results for electron damping obtained with ORA in TORIC , but algorithm uses B = 0

  9. Planned Work for FY2003 • Implement TORIC with modified wave equation within the TRANSP - FPPRF module [collaboration with D. McCune, M. Brambilla, and F. Meo]: • Can then perform routine transport analysis of NSTX discharges using HHFW heating. • Immediately begin HHFW current drive studies using the modified full-wave module: • Code already coupled to Ehst - Karney efficiency parameterization.

More Related