veltest.nsl

add edmexperiment "Velocity Test Experiment" {
	set phase_steps 20 #phase average loops
#	set bounces 10000 #number of bounces to run the simulation over
	set lifetime 20	#The lifetime of a particle (seconds)

	#Set the distance to 'lift' the particle off of the bottles
	#surface after a collision has happened. This is in to fix
	#a few rather nasty floating-point rounding errors in the
	#raytracer.
	set collision_compensation_distance 1e-15 # (Meters)
	
	#Gravity!
#	set gravity on
	
	add midpointsolver "midsolver" {
		set steptime 3e-5 # (Seconds)
	}

	ignore add rungekutta_solver "rksolver" {
		set initial_steptime 1e-6
		set tolerance 0.001
		set minimum_step 1e-7
	}
	
	add container "bottle" {
		add volume_cylinder "mainbottle" {
			set type negative
			set radius 0.235 # (Meters)
			set height 0.12 # (Meters)
			set reflection diffuse
		}
	}

	add particle neutron {
		set gamma -29.1646870E+6  # (Hertz per tesla)

		#vary velocity from 1 to 200 in 50 steps
		#set velocity 3 # (Meters per second)
		#vary velocity from 100 to 500 in 3 steps
		set velocity 500
		#This will 'spread out' the velocity to a random point on a
		#maxwell-boltzmann distribution of velocities
		#set maxwelldistribution on
		# This sets the maximum height for cutting off the velocity
		#set maxwelliancutoff 0.79 # (Meters)
		
		# set up the velocity VECTOR (this will be scaled to velocity)
		set vy 0 # (Meters per second)
		set vx 1 # (Meters per second)
		set vz 0 # (Meters per second)
		
		# Set up the starting position
		#set startvolume mainbottle
		#set startposition edge
		set z 0.1 	# (Meters)
		set x 0		# (Meters)
		set y 0		# (Meters)
		
		# Set up initial spin angle
		set start_spin_polar_angle 0.5 # * pi (Radians)
		set spin_phase 0 # * pi (Radians)
		
		#set mass 1.6749e-27
	}
	
	add bfield "bfield" {
		set b0z 	1.0E-6 # (Tesla)
		add linear_zgradient "bfield gradient" {
			set D_DZ 1e-9 # (Tesla per meter)
		}
	}
	
	add efield "Efield" {
		set e0z		8.e5 # (Volts per meter)
	}
	
	add edmreporter "edmreporter" {
		set output_file "veltest.txt"
		set volaverage_dbdz true
	}
}
tree