|Message: Fano cavity simulations in a magnetic field?||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
I am performing fano cavity simulations (similar to extended example fano cavity 2) for my work to investigate charged particle transport in the presence of a magnetic field. I am following the methods of Malkov and Rogers (MedPhys 43 4447-4458 (2016)) for this work. Essentially, three geometries are considered for these simulations, two of which consist of cylindrical slabs with the thickness of the gas volume equal to 0.2 cm or 2 cm and the third geometry is a simplified ion chamber. The particle source for all geometries are monoenergetic electrons of various starting energies (0.01, 0.1, 1, and 10 MeV). For the slab geometries, a line source in the center of the slabs are considered whereas the source for the ion chamber geometry is uniform per unit mass throughout the geometry. The electron emission is isotropic for all geometries. I originally started using Geant4 10.02, but switched to geant4 v10.04.p01 as it was clear that significant improvements were made to the GoudsmitSaunderson MSC model and the UseSafetyPlus step limit type. The parameters used are shown in the attached table. Agreement within 0.16% was observed for all energies and geometries considered (see attached plot). With these parameters I started performing simulations with a 1.5 tesla magnetic field present. I followed extended example field01 and wrote my own magnetic field class to control all of the parameters related to magnetic field transport. I've run into an interesting problem and I was hoping someone could give me some ideas to help with data interpretation. Attached are the simulated dose/energy fluence in the two slab geometries for all the electron energies considered. The transport parameters related to the magnetic field are also attached. As you can see, there is a significant drop in the dose/energy fluence at particular energy-geometry combinations, specifically 0.1 MeV in the slab 0.2 cm geometry and 1 MeV in the slab 2 cm geometry. Furthermore, when the dose/energy fluence is simulated as a function of magnetic field strength for these two energy-geometry combinations (attached), the results are very puzzling. I am having trouble understanding why the simulated dose/energy fluence decreases with a magnetic field present and why the dose/energy fluence from some energy-geometry combinations are more sensitive to the presence of the magnetic field.
Any ideas to help explain this trend would be greatly appreciated.
Attachment: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2018/04/05/13.44-74485-fano_cavity_0T.PNG http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2018/04/05/13.44-50221-_transport_parameters.PNG http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2018/04/05/13.44-23224-ial_1,5T_slab_results.PNG http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2018/04/05/13.44-30325-_transport_parameters.PNG http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2018/04/05/13.44-30137-gnetic_field_strength.PNG