| |
Records |
Links |
|
Author |
Gimenez-Alventosa, V.; Gimenez, V.; Oliver, S. |
|
| |
Title |
PenRed: An extensible and parallel Monte-Carlo framework for radiation transport based on PENELOPE |
Type |
Journal Article |
| |
Year |
2021 |
Publication |
Computer Physics Communications |
Abbreviated Journal |
Comput. Phys. Commun. |
|
| |
Volume |
267 |
Issue |
|
Pages |
108065 - 12pp |
|
| |
Keywords |
Radiation transport; Monte Carlo simulation; Electron-photon showers; Parallel computing; MPI; Medical physics |
|
| |
Abstract |
Monte Carlo methods provide detailed and accurate results for radiation transport simulations. Unfortunately, the high computational cost of these methods limits its usage in real-time applications. Moreover, existing computer codes do not provide a methodology for adapting these kinds of simulations to specific problems without advanced knowledge of the corresponding code system, and this restricts their applicability. To help solve these current limitations, we present PenRed, a general-purpose, standalone, extensible and modular framework code based on PENELOPE for parallel Monte Carlo simulations of electron-photon transport through matter. It has been implemented in C++ programming language and takes advantage of modern object-oriented technologies. In addition, PenRed offers the capability to read and process DICOM images as well as to construct and simulate image-based voxelized geometries, so as to facilitate its usage in medical applications. Our framework has been successfully verified against the original PENELOPE Fortran code. Furthermore, the implemented parallelism has been tested showing a significant improvement in the simulation time without any loss in precision of results. Program summary Program title: PenRed: Parallel Engine for Radiation Energy Deposition. CPC Library link to program files: https://doi .org /10 .17632/rkw6tvtngy.1 Licensing provision: GNU Affero General Public License (AGPL). Programming language: C++ standard 2011. Nature of problem: Monte Carlo simulations usually require a huge amount of computation time to achieve low statistical uncertainties. In addition, many applications necessitate particular characteristics or the extraction of specific quantities from the simulation. However, most available Monte Carlo codes do not provide an efficient parallel and truly modular structure which allows users to easily customise their code to suit their needs without an in-depth knowledge of the code system. Solution method: PenRed is a fully parallel, modular and customizable framework for Monte Carlo simulations of the passage of radiation through matter. It is based on the PENELOPE [1] code system, from which inherits its unique physics models and tracking algorithms for charged particles. PenRed has been coded in C++ following an object-oriented programming paradigm restricted to the C++11 standard. Our engine implements parallelism via a double approach: on the one hand, by using standard C++ threads for shared memory, improving the access and usage of the memory, and, on the other hand, via the MPI standard for distributed memory infrastructures. Notice that both kinds of parallelism can be combined together in the same simulation. Moreover, both threads and MPI processes, can be balanced using the builtin load balance system (RUPER-LB [30]) to maximise the performance on heterogeneous infrastructures. In addition, PenRed provides a modular structure with methods designed to easily extend its functionality. Thus, users can create their own independent modules to adapt our engine to their needs without changing the original modules. Furthermore, user extensions will take advantage of the builtin parallelism without any extra effort or knowledge of parallel programming. Additional comments including restrictions and unusual features: PenRed has been compiled in linux systems withg++ of GCC versions 4.8.5, 7.3.1, 8.3.1 and 9; clang version 3.4.2 and intel C++ compiler (icc) version 19.0.5.281. Since it is a C++11-standard compliant code, PenRed should be able to compile with any compiler with C++11 support. In addition, if the code is compiled without MPI support, it does not require any non standard library. To enable MPI capabilities, the user needs to install whatever available MPI implementation, such as openMPI [24] or mpich [25], which can be found in the repositories of any linux distribution. Finally, to provide DICOM processing support, PenRed can be optionally compiled using the dicom toolkit (dcmtk) [32] library. Thus, PenRed has only two optional dependencies, an MPI implementation and the dcmtk library. |
|
| |
Address |
[Gimenez-Alventosa, V] Univ Politecn Valencia, Inst Instrumentac Imagen Mol I3M, Ctr Mixto CSIC, Cami Vera S-N, Valencia 46022, Spain, Email: vicent.gimenez@i3m.upv.es; |
|
| |
Corporate Author |
|
Thesis |
|
|
| |
Publisher |
Elsevier |
Place of Publication |
|
Editor |
|
|
| |
Language |
English |
Summary Language |
|
Original Title |
|
|
| |
Series Editor |
|
Series Title |
|
Abbreviated Series Title |
|
|
| |
Series Volume |
|
Series Issue |
|
Edition |
|
|
| |
ISSN  |
0010-4655 |
ISBN |
|
Medium |
|
|
| |
Area |
|
Expedition |
|
Conference |
|
|
| |
Notes |
WOS:000678508900001 |
Approved |
no |
|
| |
Is ISI |
yes |
International Collaboration |
no |
|
| |
Call Number |
IFIC @ pastor @ |
Serial |
4907 |
|
| Permanent link to this record |
| |
|
| |
|
Author |
Oliver, S.; Rodriguez Bosca, S.; Gimenez-Alventosa, V. |
|
| |
Title |
Enabling particle transport on CAD-based geometries for radiation simulations with penRed |
Type |
Journal Article |
| |
Year |
2024 |
Publication |
Computer Physics Communications |
Abbreviated Journal |
Comput. Phys. Commun. |
|
| |
Volume |
298 |
Issue |
|
Pages |
109091 - 11pp |
|
| |
Keywords |
Radiation transport; PENELOPE physics; Monte Carlo simulation; PenRed; CAD; Triangular surface mesh |
|
| |
Abstract |
Geometry construction is a fundamental aspect of any radiation transport simulation, regardless of the Monte Carlo code being used. Typically, this process is tedious, time-consuming, and error-prone. The conventional approach involves defining geometries using mathematical objects or surfaces. However, this method comes with several limitations, especially when dealing with complex models, particularly those with organic shapes. Furthermore, since each code employs its own format and methodology for defining geometries, sharing and reproducing simulations among researchers becomes a challenging task. Consequently, many codes have implemented support for simulating over geometries constructed via Computer-Aided Design (CAD) tools. Unfortunately, this feature is lacking in penRed and other PENELOPE physics-based codes. Therefore, the objective of this work is to implement such support within the penRed framework. New version program summary Program Title: Parallel Engine for Radiation Energy Deposition (penRed) CPC Library link to program files: https://doi.org/10.17632/rkw6tvtngy.2 Developer's repository link: https://github.com/PenRed/PenRed Code Ocean capsule: https://codeocean.com/capsule/1041417/tree Licensing provisions: GNU Affero General Public License v3 Programming language: C++ standard 2011. Journal reference of previous version: V. Gimenez-Alventosa, V. Gimenez Gomez, S. Oliver, PenRed: An extensible and parallel Monte-Carlo framework for radiation transport based on PENELOPE, Computer Physics Communications 267 (2021) 108065. doi:https://doi.org/10.1016/j.cpc.2021.108065. Does the new version supersede the previous version?: Yes Reasons for the new version: Implements the capability to simulate on CAD constructed geometries, among many other features and fixes. Summary of revisions: All changes applied through the code versions are summarized in the file CHANGELOG.md in the repository package. Nature of problem: While Monte Carlo codes have proven valuable in simulating complex radiation scenarios, they rely heavily on accurate geometrical representations. In the same way as many other Monte Carlo codes, penRed employs simple geometric quadric surfaces like planes, spheres and cylinders to define geometries. However, since these geometric models offer a certain level of flexibility, these representations have limitations when it comes to simulating highly intricate and irregular shapes. Anatomic structures, for example, require detailed representations of organs, tissues and bones, which are difficult to achieve using basic geometric objects. Similarly, complex devices or intricate mechanical systems may have designs that cannot be accurately represented within the constraints of such geometric models. Moreover, when the complexity of the model increases, geometry construction process becomes more difficult, tedious, time-consuming and error-prone [2]. Also, as each Monte Carlo geometry library uses its own format and construction method, reproducing the same geometry among different codes is a challenging task. Solution method: To face the problems stated above, the objective of this work is to implement the capability to simulate using irregular and adaptable meshed geometries in the penRed framework. This kind of meshes can be constructed using Computer-Aided Design (CAD) tools, the use of which is very widespread and streamline the design process. This feature has been implemented in a new geometry module named “MESH_BODY” specific for this kind of geometries. This one is freely available and usable within the official penRed package1. It can be used since penRed version 1.9.3b and above. |
|
| |
Address |
[Oliver, S.] Univ Politecn Valencia, Inst Seguridad Ind Radiofis & Medioambiental ISIRY, Cami Vera S-N, Valencia 46022, Spain |
|
| |
Corporate Author |
|
Thesis |
|
|
| |
Publisher |
Elsevier |
Place of Publication |
|
Editor |
|
|
| |
Language |
English |
Summary Language |
|
Original Title |
|
|
| |
Series Editor |
|
Series Title |
|
Abbreviated Series Title |
|
|
| |
Series Volume |
|
Series Issue |
|
Edition |
|
|
| |
ISSN |
0010-4655 |
ISBN |
|
Medium |
|
|
| |
Area |
|
Expedition |
|
Conference |
|
|
| |
Notes |
WOS:001172840800001 |
Approved |
no |
|
| |
Is ISI |
yes |
International Collaboration |
yes |
|
| |
Call Number |
IFIC @ pastor @ |
Serial |
6077 |
|
| Permanent link to this record |
| |
|
| |
|
Author |
Zhang, X.; Xiao, Y.T.; Gimeno, B. |
|
| |
Title |
Multipactor Suppression by a Resonant Static Magnetic Field on a Dielectric Surface |
Type |
Journal Article |
| |
Year |
2020 |
Publication |
IEEE Transactions on Electron Devices |
Abbreviated Journal |
IEEE Trans. Electron Devices |
|
| |
Volume |
67 |
Issue |
12 |
Pages |
5723-5728 |
|
| |
Keywords |
Radio frequency; Dielectrics; Magnetic resonance; Discharges (electric); Surface discharges; Surface waves; Electrostatics; Monte Carlo simulation; multipactor discharge; orthogonal waves; resonant static magnetic field; secondary electron yield |
|
| |
Abstract |
In this article, we study the suppression of the multipactor phenomenon on a dielectric surface by a resonant static magnetic field. A homemade Monte Carlo algorithm is developed for multipactor simulations on a dielectric surface driven by two orthogonal radio frequency (RF) electric field components. When the static magnetic field is perpendicular to the tangential and normal RF electric fields, it is shown that if the normal electric field lags the tangential electric field by pi/2, the superposition of the normal and tangential electric fields will trigger a gyro-acceleration of the electron cloud and restrain the multipactor discharge effectively. By contrast, when the normal electric field is in advance of the tangential electric field by pi/2, the difference between the normal and tangential electric fields drives gyro-motion of the electron cloud. Consequently, two enhanced discharge zones are inevitable. The suppression effects of the resonant static magnetic field that is parallel to the tangential RF electric field or to the normal RF electric field are also presented. |
|
| |
Address |
[Zhang, Xue; Xiao, Yuting] Xiangtan Univ, Sch Automat & Elect Informat, Xiangtan 411105, Hunan, Peoples R China, Email: zhangxue.iecas@yahoo.com; |
|
| |
Corporate Author |
|
Thesis |
|
|
| |
Publisher |
Ieee-Inst Electrical Electronics Engineers Inc |
Place of Publication |
|
Editor |
|
|
| |
Language |
English |
Summary Language |
|
Original Title |
|
|
| |
Series Editor |
|
Series Title |
|
Abbreviated Series Title |
|
|
| |
Series Volume |
|
Series Issue |
|
Edition |
|
|
| |
ISSN |
0018-9383 |
ISBN |
|
Medium |
|
|
| |
Area |
|
Expedition |
|
Conference |
|
|
| |
Notes |
WOS:000594337700064 |
Approved |
no |
|
| |
Is ISI |
yes |
International Collaboration |
yes |
|
| |
Call Number |
IFIC @ pastor @ |
Serial |
4627 |
|
| Permanent link to this record |
| |
|
| |
|
Author |
Muñoz, E.; Barrio, J.; Bernabeu, J.; Etxebeste, A.; Lacasta, C.; Llosa, G.; Ros, A.; Roser, J.; Oliver, J.F. |
|
| |
Title |
Study and comparison of different sensitivity models for a two-plane Compton camera |
Type |
Journal Article |
| |
Year |
2018 |
Publication |
Physics in Medicine and Biology |
Abbreviated Journal |
Phys. Med. Biol. |
|
| |
Volume |
63 |
Issue |
13 |
Pages |
135004 - 19pp |
|
| |
Keywords |
Compton camera imaging; MLEM; Monte Carlo simulations; image quality |
|
| |
Abstract |
Given the strong variations in the sensitivity of Compton cameras for the detection of events originating from different points in the field of view (FoV), sensitivity correction is often necessary in Compton image reconstruction. Several approaches for the calculation of the sensitivity matrix have been proposed in the literature. While most of these models are easily implemented and can be useful in many cases, they usually assume high angular coverage over the scattered photon, which is not the case for our prototype. In this work, we have derived an analytical model that allows us to calculate a detailed sensitivity matrix, which has been compared to other sensitivity models in the literature. Specifically, the proposed model describes the probability of measuring a useful event in a two-plane Compton camera, including the most relevant physical processes involved. The model has been used to obtain an expression for the system and sensitivity matrices for iterative image reconstruction. These matrices have been validated taking Monte Carlo simulations as a reference. In order to study the impact of the sensitivity, images reconstructed with our sensitivity model and with other models have been compared. Images have been reconstructed from several simulated sources, including point-like sources and extended distributions of activity, and also from experimental data measured with Na-22 sources. Results show that our sensitivity model is the best suited for our prototype. Although other models in the literature perform successfully in many scenarios, they are not applicable in all the geometrical configurations of interest for our system. In general, our model allows to effectively recover the intensity of point-like sources at different positions in the FoV and to reconstruct regions of homogeneous activity with minimal variance. Moreover, it can be employed for all Compton camera configurations, including those with low angular coverage over the scatterer. |
|
| |
Address |
[Munoz, Enrique; Barrio, John; Bernabeu, Jose; Etxebeste, Ane; Lacasta, Carlos; Llosa, Gabriela; Ros, Ana; Roser, Jorge; Oliver, Josep F.] Univ Valencia, CSIC, Inst Fis Corpuscular, IFIC, Valencia, Spain, Email: Enrique.Munoz@ific.uv.es |
|
| |
Corporate Author |
|
Thesis |
|
|
| |
Publisher |
Iop Publishing Ltd |
Place of Publication |
|
Editor |
|
|
| |
Language |
English |
Summary Language |
|
Original Title |
|
|
| |
Series Editor |
|
Series Title |
|
Abbreviated Series Title |
|
|
| |
Series Volume |
|
Series Issue |
|
Edition |
|
|
| |
ISSN |
0031-9155 |
ISBN |
|
Medium |
|
|
| |
Area |
|
Expedition |
|
Conference |
|
|
| |
Notes |
WOS:000436390800004 |
Approved |
no |
|
| |
Is ISI |
yes |
International Collaboration |
no |
|
| |
Call Number |
IFIC @ pastor @ |
Serial |
3639 |
|
| Permanent link to this record |
| |
|
| |
|
Author |
Etxebeste, A.; Barrio, J.; Bernabeu, J.; Lacasta, C.; Llosa, G.; Muñoz, E.; Ros, A.; Oliver, J.F. |
|
| |
Title |
Study of sensitivity and resolution for full ring PET prototypes based on continuous crystals and analytical modeling of the light distribution |
Type |
Journal Article |
| |
Year |
2019 |
Publication |
Physics in Medicine and Biology |
Abbreviated Journal |
Phys. Med. Biol. |
|
| |
Volume |
64 |
Issue |
3 |
Pages |
035015 - 17pp |
|
| |
Keywords |
continuous crystals; NEMA NU 4-2008; positron emission tomography (PET); Monte Carlo simulations; image reconstruction; depth of interaction |
|
| |
Abstract |
Sensitivity and spatial resolution are the main parameters to maximize in the performance of a PET scanner. For this purpose, detectors consisting of a combination of continuous crystals optically coupled to segmented photodetectors have been employed. With the use of continuous crystals the sensitivity is increased with respect to the pixelated crystals. In addition, spatial resolution is no longer limited to the crystal size. The main drawback is the difficulty in determining the interaction position. In this work, we present the characterization of the performance of a full ring based on cuboid continuous crystals coupled to SiPMs. To this end, we have employed the simulations developed in a previous work for our experimental detector head. Sensitivity could be further enhanced by using tapered crystals. This enhancement is obtained by increasing the solid angle coverage, reducing the wedge-shaped gaps between contiguous detectors. The performance of the scanners based on both crystal geometries was characterized following NEMA NU 4-2008 standardized protocol in order to compare them. An average sensitivity gain over the entire axial field of view of 13.63% has been obtained with tapered geometry while similar performance of the spatial resolution has been proven with both scanners. The activity at which NECR and true peak occur is smaller and the peak value is greater for tapered crystals than for cuboid crystals. Moreover, a higher degree of homogeneity was obtained in the sensitivity map due to the tighter packing of the crystals, which reduces the gaps and results in a better recovery of homogeneous regions than for the cuboid configuration. Some of the results obtained, such as spatial resolution, depend on the interaction position estimation and may vary if other method is employed. |
|
| |
Address |
[Etxebeste, Ane; Barrio, John; Bernabeu, Jose; Lacasta, Carlos; Llosa, Gabriela; Munoz, Enrique; Ros, Ana; Oliver, Josef F.] Univ Valencia, CSIC, Inst Fis Corpuscular, Valencia, Spain, Email: ane.etxebeste@ific.uv.es |
|
| |
Corporate Author |
|
Thesis |
|
|
| |
Publisher |
Iop Publishing Ltd |
Place of Publication |
|
Editor |
|
|
| |
Language |
English |
Summary Language |
|
Original Title |
|
|
| |
Series Editor |
|
Series Title |
|
Abbreviated Series Title |
|
|
| |
Series Volume |
|
Series Issue |
|
Edition |
|
|
| |
ISSN |
0031-9155 |
ISBN |
|
Medium |
|
|
| |
Area |
|
Expedition |
|
Conference |
|
|
| |
Notes |
WOS:000457182500004 |
Approved |
no |
|
| |
Is ISI |
yes |
International Collaboration |
yes |
|
| |
Call Number |
IFIC @ pastor @ |
Serial |
3897 |
|
| Permanent link to this record |
