Added an MPEX.in input file (40940c7b) · Commits · Kumar, Atul / LinearFokkerPlanck Axisymmetric

BfieldData/Bfield_c.tif

0 → 100644

+75.9 KiB

File added.

No diff preview for this file type.

View file

BfieldData/MPEX_Bfield_scenario_14.tiff

0 → 100644

+330 KiB

Loading image diff...

InputFiles/xp_MPEX.in

0 → 100644

+92 −0

Original line number	Diff line number	Diff line
		&params_nml
		! Simulation name:
		! ===============
		params%fileDescriptor = 'MPEX_2021_06_22l',

		! Magnetic field input data:
		! =========================
		params%repoDir ="/home/nfc/myRepos/LinearFokkerPlanck_Axisymmetric",
		params%BFieldFileDir = "/BfieldData",
		params%BFieldFile = "/MPEX_Bfield_scenario_14.txt",
		params%nz = 501, ! Number of points in BFieldFile

		! Simulation conditions:
		! =====================
		params%NC = 500000, ! Total number of particles
		params%NS = 40000, ! Total number of time steps
		params%dt = 2E-8, ! Time step in [s]
		params%G = 1E+22, ! Real particle injection rate

		! Time steps to record:
		! ====================
		params%jstart = 1, ! start frame
		params%jend = 40000, ! end frame
		params%jincr = 1000, ! increment

		! Domain geometry:
		! ====================
		params%dtheta = 6.2832 ! Radians
		params%r1 = 0.00 ! Inner radius at reference location [m]
		params%r2 = 0.05 ! Outer radius at reference location [m]
		params%zmax = +7.75, ! Boundary 2 (right)
		params%zmin = -1.95, ! Boundary 1 (left)
		params%NZmesh = 400, ! Mesh number of elements

		! Physics:
		! ======================
		params%iSave = .true.,
		params%iPush = .true., ! Enable RK4 particle pusher
		params%iColl = .true., ! Collisions or no collisions
		params%iHeat = .true., ! Enable RF heating operator
		params%iPotential = .false., ! Enable electric potential
		params%iDrag = .true., ! Enable plasma drag

		! Collision operator conditions:
		! ==============================
		params%Te0 = 25., ! eV
		params%Ti0 = 25., ! eV
		params%ne0 = 0.1E+19, ! m**-3
		params%Aion = 1., ! Ion mass
		params%Zion = 1., ! Charge number for main ion species
		params%species_a = 2, ! =1 for electrons, =2 for ions
		params%elevel = 10, ! Energy at which slowing down power is calculated E = elevel*Te0
		params%CollOperType = 2, ! 1 for Boozer-Only, 2 for Boozer-Kim

		! Particle Boundary Conditions:
		! ============================
		params%BC_Type = 1, ! 1: Isotropic plasma source, 2: NBI, 3: Periodic
		params%BC_zp_mean = 0.0, ! Mean particle injection
		params%BC_zp_std = 0.4, ! STD of particle injection
		params%BC_Ep = 2000., ! Drift kinetic energy
		params%BC_Tp = 25, ! Thermal kinetic energy
		params%BC_xip = 0.707, ! Mean pitch angle where xip = cos(theta) = vpar/v

		! Particle Initial Conditions:
		!===================
		params%IC_Type = 1, ! 1: uniform load, 2: gaussian load
		params%IC_zp_mean = 0.0, ! Mean particle injection
		params%IC_zp_std = 0.4, ! STD of particle injection
		params%IC_Ep = 0., ! Drift kinetic energy
		params%IC_Tp = 25, ! Thermal kinetic energy
		params%IC_xip = 0.707, ! Mean pitch angle where xip = cos(theta) = vpar/v

		! RF heating operator conditions:
		! ===============================
		params%f_RF = 17.073E+6 ! RF frequency
		params%zRes1 = +4.5 ! Defines interval where RF is present
		params%zRes2 = +5.6 ! Defines interval where RF is present
		params%kper = 100. ! Perpendicular wave number of EBW
		params%kpar = 20. ! Parallel wave number of EBW
		params%Prf = 50000. ! Mean RF power absorbed via cyclotron interaction
		params%n_harmonic = 1 ! Cyclotron harmonic

		! Electric potential conditions:
		! =============================
		params%s1 = 0,
		params%s2 = 1.3,
		params%s3 = 4.0,
		params%phi1 = 0.,
		params%phi2 = 0,
		params%phi3 = 0.,

		/

InputFiles/xp_TestDrag.in

+11 −11

Original line number	Diff line number	Diff line
		&params_nml
		! Simulation name:
		! ===============
		params%fileDescriptor = '2021_03_28d',
		params%fileDescriptor = '2021_06_21b',

		! Magnetic field input data:
		! =========================
		params%repoDir ="/home/nfc/myRepos/LinearFokkerPlanck_Axisymmetric",
		params%BFieldFileDir = "/BfieldData",
		params%BFieldFile = "/Bfield_b.txt",
		params%BFieldFile = "/Bfield_c.txt",
		params%nz = 501, ! Number of points in BFieldFile

		! Simulation conditions:
		! =====================
		params%NC = 500000, ! Total number of particles
		params%NS = 20000, ! Total number of time steps
		params%dt = 0.5E-7, ! Time step in [s]
		params%dt = 1E-8, ! Time step in [s]
		params%G = 2.5E+20, ! Real particle injection rate

		! Time steps to record:
		@@ -28,8 +28,8 @@ params%jincr = 500,
		params%dtheta = 6.2832 ! Radians
		params%r1 = 0.03 ! Inner radius at reference location [m]
		params%r2 = 0.05 ! Outer radius at reference location [m]
		params%zmax = +3.0, ! Boundary 2 (right)
		params%zmin = -3.0, ! Boundary 1 (left)
		params%zmax = +1.45, ! Boundary 2 (right)
		params%zmin = -1.45, ! Boundary 1 (left)
		params%NZmesh = 200, ! Mesh number of elements

		! Physics:
		@@ -43,8 +43,8 @@ params%iDrag = .true.,

		! Collision operator conditions:
		! ==============================
		params%Te0 = 600., ! eV
		params%Ti0 = 600., ! eV
		params%Te0 = 25., ! eV
		params%Ti0 = 25., ! eV
		params%ne0 = 0.5E+19, ! m**-3
		params%Aion = 1., ! Ion mass
		params%Zion = 1., ! Charge number for main ion species
		@@ -56,9 +56,9 @@ params%CollOperType = 2,
		! ============================
		params%BC_Type = 1, ! 1: Isotropic plasma source, 2: NBI, 3: Periodic
		params%BC_zp_mean = 0.0, ! Mean particle injection
		params%BC_zp_std = 0.3, ! STD of particle injection
		params%BC_zp_std = 0.1, ! STD of particle injection
		params%BC_Ep = 2000., ! Drift kinetic energy
		params%BC_Tp = 600, ! Thermal kinetic energy
		params%BC_Tp = 25, ! Thermal kinetic energy
		params%BC_xip = 0.707, ! Mean pitch angle where xip = cos(theta) = vpar/v

		! Particle Initial Conditions:
		@@ -67,7 +67,7 @@ params%IC_Type = 1,
		params%IC_zp_mean = 0.0, ! Mean particle injection
		params%IC_zp_std = 0.3, ! STD of particle injection
		params%IC_Ep = 0., ! Drift kinetic energy
		params%IC_Tp = 10, ! Thermal kinetic energy
		params%IC_Tp = 25, ! Thermal kinetic energy
		params%IC_xip = 0.707, ! Mean pitch angle where xip = cos(theta) = vpar/v

		! RF heating operator conditions:

src/Modules.f90

+80 −71

Original line number	Diff line number	Diff line
		@@ -715,6 +715,8 @@ SUBROUTINE SaveData(output,dir1)
		OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		WRITE(8) output%tp
		CLOSE(unit=8)

		IF (.false.) THEN
		fileName = trim(trim(dir1)//'/'//'m.out')
		OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		WRITE(8) output%m
		@@ -747,6 +749,7 @@ SUBROUTINE SaveData(output,dir1)
		OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		WRITE(8) output%ddBp
		CLOSE(unit=8)
		END IF

		! Saving mesh-defined quantities:
		! -------------------------------------------------------------------------
		@@ -762,22 +765,26 @@ SUBROUTINE SaveData(output,dir1)
		OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		WRITE(8) output%nU
		CLOSE(unit=8)
		fileName = trim(trim(dir1)//'/'//'unU_mesh.out')
		OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		WRITE(8) output%unU
		CLOSE(unit=8)

		!fileName = trim(trim(dir1)//'/'//'unU_mesh.out')
		!OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		!WRITE(8) output%unU
		!CLOSE(unit=8)

		fileName = trim(trim(dir1)//'/'//'nUE_mesh.out')
		OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		WRITE(8) output%nUE
		CLOSE(unit=8)
		fileName = trim(trim(dir1)//'/'//'Ppar_mesh.out')
		OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		WRITE(8) output%Ppar
		CLOSE(unit=8)
		fileName = trim(trim(dir1)//'/'//'Pper_mesh.out')
		OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		WRITE(8) output%Pper
		CLOSE(unit=8)

		!fileName = trim(trim(dir1)//'/'//'Ppar_mesh.out')
		!OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		!WRITE(8) output%Ppar
		!CLOSE(unit=8)
		!fileName = trim(trim(dir1)//'/'//'Pper_mesh.out')
		!OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		!WRITE(8) output%Pper
		!CLOSE(unit=8)

		fileName = trim(trim(dir1)//'/'//'Tpar_mesh.out')
		OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		WRITE(8) output%Tpar
		@@ -790,10 +797,12 @@ SUBROUTINE SaveData(output,dir1)
		OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		WRITE(8) output%U
		CLOSE(unit=8)
		fileName = trim(trim(dir1)//'/'//'ddB_mesh.out')
		OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		WRITE(8) output%ddB
		CLOSE(unit=8)

		!fileName = trim(trim(dir1)//'/'//'ddB_mesh.out')
		!OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		!WRITE(8) output%ddB
		!CLOSE(unit=8)

		fileName = trim(trim(dir1)//'/'//'E_mesh.out')
		OPEN(unit=8,file=fileName,form="unformatted",status="unknown")
		WRITE(8) output%E