| |
Records |
Links |
|
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 |
Figueroa, D.G.; Florio, A.; Torrenti, F.; Valkenburg, W. |
|
| |
Title |
CosmoLattice: A modern code for lattice simulations of scalar and gauge field dynamics in an expanding universe |
Type |
Journal Article |
| |
Year |
2023 |
Publication |
Computer Physics Communications |
Abbreviated Journal |
Comput. Phys. Commun. |
|
| |
Volume |
283 |
Issue |
|
Pages |
108586 - 13pp |
|
| |
Keywords |
Early universe; Real-time lattice simulations; Gauge -invariant lattice techniques |
|
| |
Abstract |
This paper describes CosmoGattice, a modern package for lattice simulations of the dynamics of interacting scalar and gauge fields in an expanding universe. CosmoGattice incorporates a series of features that makes it very versatile and powerful: i) it is written in C++ fully exploiting the object oriented programming paradigm, with a modular structure and a clear separation between the physics and the technical details, ii) it is MPI-based and uses a discrete Fourier transform parallelized in multiple spatial dimensions, which makes it specially appropriate for probing scenarios with well -separated scales, running very high resolution simulations, or simply very long ones, iii) it introduces its own symbolic language, defining field variables and operations over them, so that one can introduce differential equations and operators in a manner as close as possible to the continuum, iv) it includes a library of numerical algorithms, ranging from O(delta t(2)) to O(delta t(10)) methods, suitable for simulating global and gauge theories in an expanding grid, including the case of 'self-consistent' expansion sourced by the fields themselves. Relevant observables are provided for each algorithm (e.g. energy densities, field spectra, lattice snapshots) and we note that, remarkably, all our algorithms for gauge theories (Abelian or non-Abelian) always respect the Gauss constraint to machine precision. Program summary Program Title:: CosmoGattice CPC Library link to program files: https://doi .org /10 .17632 /44vr5xssc6 .1 Developer's repository link: http://github .com /cosmolattice /cosmolattice Licensing provisions: MIT Programming language: C++, MPI Nature of problem: The phenomenology of high energy physics in the early universe is typically characterized by non-linear dynamics, which cannot be captured accurately with analytical techniques. In order to fully understand the non-linearities developed in a given scenario, one needs to carry out lattice simulations. A number of public packages for lattice simulations have appeared over the years, but most of them are only capable of simulating scalar fields. However, realistic models of particle physics do contain other kind of field species, such as (Abelian or non-Abelian) gauge fields, whose non-linear dynamics can also play a relevant role in the early universe. Tensor modes representing gravitational waves are also naturally expected in many scenarios. Solution method: CosmoGattice represents a modern code for lattice simulations of scalar-gauge field theories in an expanding universe. It allows for the simulation of the evolution of interacting (singlet) scalar fields, charged scalar fields under U(1) and/or SU(2) gauge groups, and the corresponding associated Abelian and/or non-Abelian gauge fields. From version 1.1 onward, CosmoGattice also allows to simulate the production of gravitational waves. Simulations can be done either in a flat space-time background, or in a homogeneous and isotropic (spatially flat) expanding FLRW background. CosmoGattice provides symplectic integrators, with accuracy ranging from O (delta t(2)) up to O(delta t(10)), to simuate the non-linear dynamics of the appropriate fields in comoving three-dimensional lattices. The code is parallelized with MPI, and uses a discrete Fourier Transform parallelized in multiple spatial dimensions, which makes it a very powerful code for probing physical problems with well-separated scales. Moreover, the code has been designed as a `platform' to implement any system of dynamical equations suitable for discretization on a lattice. |
|
| |
Address |
[Figueroa, Daniel G.] CSIC, Inst Fis Corpuscular IFIC, Valencia, Spain, Email: f.torrenti@unibas.ch |
|
| |
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:000899506700008 |
Approved |
no |
|
| |
Is ISI |
yes |
International Collaboration |
yes |
|
| |
Call Number |
IFIC @ pastor @ |
Serial |
5451 |
|
| Permanent link to this record |
| |
|
| |
|
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 |
Valdes-Cortez, C.; Mansour, I.; Rivard, M.J.; Ballester, F.; Mainegra-Hing, E.; Thomson, R.M.; Vijande, J. |
|
| |
Title |
A study of Type B uncertainties associated with the photoelectric effect in low-energy Monte Carlo simulations |
Type |
Journal Article |
| |
Year |
2021 |
Publication |
Physics in Medicine and Biology |
Abbreviated Journal |
Phys. Med. Biol. |
|
| |
Volume |
66 |
Issue |
10 |
Pages |
105014 - 14pp |
|
| |
Keywords |
Monte Carlo simulations; brachytherapy; low energy physics; photoelectric effect |
|
| |
Abstract |
Purpose. To estimate Type B uncertainties in absorbed-dose calculations arising from the different implementations in current state-of-the-art Monte Carlo (MC) codes of low-energy photon cross-sections (<200 keV). Methods. MC simulations are carried out using three codes widely used in the low-energy domain: PENELOPE-2018, EGSnrc, and MCNP. Three dosimetry-relevant quantities are considered: mass energy-absorption coefficients for water, air, graphite, and their respective ratios; absorbed dose; and photon-fluence spectra. The absorbed dose and the photon-fluence spectra are scored in a spherical water phantom of 15 cm radius. Benchmark simulations using similar cross-sections have been performed. The differences observed between these quantities when different cross-sections are considered are taken to be a good estimator for the corresponding Type B uncertainties. Results. A conservative Type B uncertainty for the absorbed dose (k = 2) of 1.2%-1.7% (<50 keV), 0.6%-1.2% (50-100 keV), and 0.3% (100-200 keV) is estimated. The photon-fluence spectrum does not present clinically relevant differences that merit considering additional Type B uncertainties except for energies below 25 keV, where a Type B uncertainty of 0.5% is obtained. Below 30 keV, mass energy-absorption coefficients show Type B uncertainties (k = 2) of about 1.5% (water and air), and 2% (graphite), diminishing in all materials for larger energies and reaching values about 1% (40-50 keV) and 0.5% (50-75 keV). With respect to their ratios, the only significant Type B uncertainties are observed in the case of the water-to-graphite ratio for energies below 30 keV, being about 0.7% (k = 2). Conclusions. In contrast with the intermediate (about 500 keV) or high (about 1 MeV) energy domains, Type B uncertainties due to the different cross-sections implementation cannot be considered subdominant with respect to Type A uncertainties or even to other sources of Type B uncertainties (tally volume averaging, manufacturing tolerances, etc). Therefore, the values reported here should be accommodated within the uncertainty budget in low-energy photon dosimetry studies. |
|
| |
Address |
[Valdes-Cortez, Christian; Ballester, Facundo; Vijande, Javier] Univ Valencia UV, Dept Fis Atom Mol & Nucl, Burjassot, Spain, Email: javier.vijande@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:000655291500001 |
Approved |
no |
|
| |
Is ISI |
yes |
International Collaboration |
yes |
|
| |
Call Number |
IFIC @ pastor @ |
Serial |
4847 |
|
| Permanent link to this record |
| |
|
| |
|
Author |
Baran, J. et al; Brzezinski, K. |
|
| |
Title |
Feasibility of the J-PET to monitor the range of therapeutic proton beams |
Type |
Journal Article |
| |
Year |
2024 |
Publication |
Physica Medica |
Abbreviated Journal |
Phys. Medica |
|
| |
Volume |
118 |
Issue |
|
Pages |
103301 - 9pp |
|
| |
Keywords |
PET; Range monitoring; J-PET; Monte Carlo simulations; Proton radiotherapy |
|
| |
Abstract |
Purpose: The aim of this work is to investigate the feasibility of the Jagiellonian Positron Emission Tomography (J -PET) scanner for intra-treatment proton beam range monitoring. Methods: The Monte Carlo simulation studies with GATE and PET image reconstruction with CASToR were performed in order to compare six J -PET scanner geometries. We simulated proton irradiation of a PMMA phantom with a Single Pencil Beam (SPB) and Spread -Out Bragg Peak (SOBP) of various ranges. The sensitivity and precision of each scanner were calculated, and considering the setup's cost-effectiveness, we indicated potentially optimal geometries for the J -PET scanner prototype dedicated to the proton beam range assessment. Results: The investigations indicate that the double -layer cylindrical and triple -layer double -head configurations are the most promising for clinical application. We found that the scanner sensitivity is of the order of 10-5 coincidences per primary proton, while the precision of the range assessment for both SPB and SOBP irradiation plans was found below 1 mm. Among the scanners with the same number of detector modules, the best results are found for the triple -layer dual -head geometry. The results indicate that the double -layer cylindrical and triple -layer double -head configurations are the most promising for the clinical application, Conclusions: We performed simulation studies demonstrating that the feasibility of the J -PET detector for PET -based proton beam therapy range monitoring is possible with reasonable sensitivity and precision enabling its pre -clinical tests in the clinical proton therapy environment. Considering the sensitivity, precision and cost-effectiveness, the double -layer cylindrical and triple -layer dual -head J -PET geometry configurations seem promising for future clinical application. |
|
| |
Address |
[Baran, Jakub; Silarski, Michal; Chug, Neha; Coussat, Aurelien; Czerwinski, Eryk; Dadgar, Meysam; Dulski, Kamil; Eliyan, Kavya, V; Gajos, Aleksander; Kacprzak, Krzysztof; Kaplon, Lukasz; Korcyl, Grzegorz; Kozik, Tomasz; Kumar, Deepak; Niedzwiecki, Szymon; Panek, Dominik; Parzych, Szymon; del Rio, Elena Perez; Simbarashe, Moyo; Sharma, Sushil; Shivani; Skurzok, Magdalena; Stepien, Ewa L.; Tayefi, Keyvan; Tayefi, Faranak; Moskal, Pawel] Jagiellonian Univ, Fac Phys Astron & Appl Comp Sci, 11 Lojasiewicza St, PL-30348 Krakow, Poland, Email: jakubbaran92@gmail.com |
|
| |
Corporate Author |
|
Thesis |
|
|
| |
Publisher |
Elsevier Sci Ltd |
Place of Publication |
|
Editor |
|
|
| |
Language |
English |
Summary Language |
|
Original Title |
|
|
| |
Series Editor |
|
Series Title |
|
Abbreviated Series Title |
|
|
| |
Series Volume |
|
Series Issue |
|
Edition |
|
|
| |
ISSN |
1120-1797 |
ISBN |
|
Medium |
|
|
| |
Area |
|
Expedition |
|
Conference |
|
|
| |
Notes |
WOS:001178648400001 |
Approved |
no |
|
| |
Is ISI |
yes |
International Collaboration |
yes |
|
| |
Call Number |
IFIC @ pastor @ |
Serial |
5990 |
|
| Permanent link to this record |
