-
Notifications
You must be signed in to change notification settings - Fork 59
RPT PARAMETERS
Radioactive particle tracking (RPT) is a non-intrusive velocimetry method which is used to study the hydrodynamics of single and multiphase systems. Launching a RPT simulation for photon count calculation in Lethe requires a solver which is rpt_3d, a parameter file, two files including detector and particle positions inside the vessel. The particle positions file includes either the experimental calibration positions or a set of generated points inside the vessel by the user. Detector positions file contains the position of detector face canter and the position of a point inside the detector on its axis. In Lethe-RPT code the middle bottom of the cylinder is considered as the origin. An example of these two files can be find here This section aims at describing the various parameters available within Lethe-RPT.
In the parameter file format, the parameters are established one by one using the following syntax for instance: 'set reactor radius= 0.1' would fix the reactor radius to 0.1 m. The arguments can be either doubles, integers, or a choice between a predefined condition. In the parameter files, comments are preceded by the sharp symbol (e.g., '#comment').
The parameter file is composed of different subsections. In the following, the principal subsections of a RPT parameter template file are explained.
This subsection contains the general information for photon count calculation which is described in detail in this part. “Particle position file” defines the file that Lethe uses to read the particle positions inside the vessel. Setting the option “export counts” to true allows exporting the results in a .csv file. The entry “counts file” specifies the mentioned .csv file name. The number of Monte-Carlo iteration (the number of traced gamma-ray from each particle position to the detector within the defined solid angle) is defined by “monte carlo iteration”. Setting “random number seed” to auto allows us to generate random numbers for the Monte-Carlo calculation. Reactor radius defines the radius of the cylindrical tank in which simulations are carried out. “Peak-to-total” ratio is the relation between the area of photopeak and the total area in the spectrum. The amount of time for which the RPT hardware records the photon count at each position is specified by “sampling time”. A sample of this subsection is shown below:
subsection rpt parameters
set particle positions file = positions1.particle
set export counts = true
set monte carlo iteration = 1000000
set random number seed = auto
set reactor radius = 0.1
set peak-to-total ratio = 0.4
set sampling time = 1
end
This subsection contains information regarding the photon count calculation which has also been used in the detector’s parameter tunning process. Setting “tuning” to false allows Lethe-RPT to calculates the photon counts at each particle position with respect to the detector. In order to tune the detector’s parameter this option musted be set to true. “Gamma-rays emitted” indicates the number gamma-rays emitted by each disintegration. RPT technique works based on the attenuation of the light which travels from the tracer particle to the detector through the reactor. Therefore, the attenuation coefficient of the detector’s crystal and that of the materials inside the reactor must be specified. These parameters respectively are defined as “attenuation coefficient detector” and “attenuation coefficient reactor”. The “dead time” specifies detector dead time parameter and “activity” declares the radioactive source activity.
subsection Parameter tuning
set tuning = false
set dead time = 2.7e-6
set activity = 9.5e6
set gamma-rays emitted = 2
set attenuation coefficient reactor = 7.1
set attenuation coefficient detector = 21.477
end
This subsection contains the specific information of the detector. “Detector position file” defines the files that Lethe uses to read the detector positions. This file includes the position of detector face center and the position of a point inside the detector on its axis. The options “radius” and “length” respectively defines the detector’s crystal radius and length.
subsection Detector parameters
set detector positions file = positions1.detector
set radius = 0.0381
set length = 0.0762
end