Commit 4f58e9cc authored by Caneses Marin, Juan Francisco's avatar Caneses Marin, Juan Francisco
Browse files

Fully developed Power constrained RF operator. This commit is prior to... 

Fully developed Power constrained RF operator. This commit is prior to proceeding with Fueling constrained Re-injection operator
parent cd34e1df

InputFiles/xp_IonSlowingDown.in

+0 −1
Original line number Diff line number Diff line
@@ -75,7 +75,6 @@ params%zRes1 = 2.4 
params%zRes2      = 3.0                                                         ! Defines interval where RF is present

params%kper       = 29920.                                                      ! Perpendicular wave number of EBW

params%kpar       = 5864.                                                       ! Parallel wave number of EBW

params%Ew         = 10000                                                       ! Wave electric field magnitude [V/m]

params%Prf        = 300000                                                      ! Mean RF power absorbed via cyclotron interaction

params%n_harmonic = 4                                                           ! Cyclotron harmonic

InputFiles/xp_TestCollisions.in

+0 −1
Original line number Diff line number Diff line
@@ -75,7 +75,6 @@ params%zRes1 = 2.4 
params%zRes2      = 3.0                                                         ! Defines interval where RF is present

params%kper       = 29920.                                                      ! Perpendicular wave number of EBW

params%kpar       = 5864.                                                       ! Parallel wave number of EBW

params%Ew         = 10000                                                       ! Wave electric field magnitude [V/m]

params%Prf        = 300000                                                      ! Mean RF power absorbed via cyclotron interaction

params%n_harmonic = 4                                                           ! Cyclotron harmonic

InputFiles/xp_TestDrag.in

+1 −2
Original line number Diff line number Diff line 
&params_nml

! Simulation name:

! ===============

params%fileDescriptor = '2021_02_22d',

params%fileDescriptor = '2021_02_22e',



! Magnetic field input data:

! =========================
@@ -75,7 +75,6 @@ params%zRes1 = 0. 
params%zRes2      = +1.5                                                        ! Defines interval where RF is present

params%kper       = 125.                                                        ! Perpendicular wave number of EBW

params%kpar       = 30.                                                         ! Parallel wave number of EBW

params%Ew         = 10000                                                       ! Wave electric field magnitude [V/m]

params%Prf        = 2000000                                                      ! Mean RF power absorbed via cyclotron interaction

params%n_harmonic = 2                                                           ! Cyclotron harmonic

src/Modules.f90

+1 −1
Original line number Diff line number Diff line
@@ -244,7 +244,7 @@ TYPE paramsTYP 
  INTEGER(i4) :: IC_Type

  REAL(r8)    :: IC_zp_mean, IC_Ep, IC_xip, IC_zp_std, IC_Tp

  ! RF heating operator conditions:

  REAL(r8)    :: f_RF, kpar, kper, Ew, zRes1, zRes2, Prf

  REAL(r8)    :: f_RF, kpar, kper, zRes1, zRes2, Prf

  INTEGER(i4) :: n_harmonic

  ! Electric potential conditions:

  REAL(r8) :: s1, s2, s3, phi1, phi2, phi3

src/MoveParticlePack.f90

+4 −3
Original line number Diff line number Diff line
@@ -337,12 +337,13 @@ mean_dkep_per = 0.5*(e_c/Ma)*(besselterm*tau_rf)**2. 
! Populate plasma structure:

plasma%Erf_hat(i) = mean_dkep_per

plasma%doppler(i) = params%kpar*(upar0)/Omega

plasma%E3_hat(i)  = plasma%Erf_hat(i)*(1. + plasma%doppler(i))



RETURN

END SUBROUTINE RFoperatorTerms



! =======================================================================================================

SUBROUTINE RFOperator(i,Prf_0,plasma,fieldspline,params)

SUBROUTINE RFOperator(i,plasma,fieldspline,params,q3_hat)

! =======================================================================================================

USE local

USE spline_fits
@@ -351,7 +352,7 @@ USE dataTYP 


IMPLICIT NONE

! Define interface variables:

REAL(r8)               , INTENT(IN)    :: Prf_0

REAL(r8)               , INTENT(IN)    :: q3_hat

INTEGER(i4)            , INTENT(IN)    :: i

TYPE(plasmaTYP)        , INTENT(INOUT) :: plasma

TYPE(paramsTYP)        , INTENT(IN)    :: params
@@ -378,7 +379,7 @@ kep_par0 = kep0*xip0**2. 
kep_per0 = kep0*(1. - xip0**2.)



! Calculate the mean RF kick:

Ew2 = (params%Prf/Prf_0)

Ew2 = (params%Prf/q3_hat)

mean_dkep_per = plasma%Erf_hat(i)*Ew2



! Calculate the change in perp, parallel and total energy: