|
|
Ferrer-Sanchez, A., Martin-Guerrero, J., Ruiz de Austri, R., Torres-Forne, A., & Font, J. A. (2024). Gradient-annihilated PINNs for solving Riemann problems: Application to relativistic hydrodynamics. Comput. Meth. Appl. Mech. Eng., 424, 116906–18pp.
Abstract: We present a novel methodology based on Physics-Informed Neural Networks (PINNs) for solving systems of partial differential equations admitting discontinuous solutions. Our method, called Gradient-Annihilated PINNs (GA-PINNs), introduces a modified loss function that forces the model to partially ignore high-gradients in the physical variables, achieved by introducing a suitable weighting function. The method relies on a set of hyperparameters that control how gradients are treated in the physical loss. The performance of our methodology is demonstrated by solving Riemann problems in special relativistic hydrodynamics, extending earlier studies with PINNs in the context of the classical Euler equations. The solutions obtained with the GA-PINN model correctly describe the propagation speeds of discontinuities and sharply capture the associated jumps. We use the relative l(2) error to compare our results with the exact solution of special relativistic Riemann problems, used as the reference ''ground truth'', and with the corresponding error obtained with a second-order, central, shock-capturing scheme. In all problems investigated, the accuracy reached by the GA-PINN model is comparable to that obtained with a shock-capturing scheme, achieving a performance superior to that of the baseline PINN algorithm in general. An additional benefit worth stressing is that our PINN-based approach sidesteps the costly recovery of the primitive variables from the state vector of conserved variables, a well-known drawback of grid-based solutions of the relativistic hydrodynamics equations. Due to its inherent generality and its ability to handle steep gradients, the GA-PINN methodology discussed in this paper could be a valuable tool to model relativistic flows in astrophysics and particle physics, characterized by the prevalence of discontinuous solutions.
|
|
|
|
Oliver, S., Rodriguez Bosca, S., & Gimenez-Alventosa, V. (2024). Enabling particle transport on CAD-based geometries for radiation simulations with penRed. Comput. Phys. Commun., 298, 109091–11pp.
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.
|
|
|
|
Angles-Castillo, A., Perez, A., & Roldan, E. (2024). Bright and dark solitons in a photonic nonlinear quantum walk: lessons from the continuum. New J. Phys., 26(2), 023004–16pp.
Abstract: We propose a nonlinear quantum walk model inspired in a photonic implementation in which the polarization state of the light field plays the role of the coin-qubit. In particular, we take profit of the nonlinear polarization rotation occurring in optical media with Kerr nonlinearity, which allows to implement a nonlinear coin operator, one that depends on the state of the coin-qubit. We consider the space-time continuum limit of the evolution equation, which takes the form of a nonlinear Dirac equation. The analysis of this continuum limit allows us to gain some insight into the existence of different solitonic structures, such as bright and dark solitons. We illustrate several properties of these solitons with numerical calculations, including the effect on them of an additional phase simulating an external electric field.
|
|
|
|
LHCb Collaboration(Aaij, R. et al), Jaimes Elles, S. J., Jashal, B. K., Martinez-Vidal, F., Oyanguren, A., Rebollo De Miguel, M., et al. (2024). Prompt and nonprompt ψ(2S) production in pPb collisions at √sNN = 8.16 TeV. J. High Energy Phys., 04(4), 111–52pp.
Abstract: The production of psi(2S) mesons in proton-lead collisions at a centre-of-mass energy per nucleon pair of root s(NN) = 8.16TeV is studied with the LHCb detector using data corresponding to an integrated luminosity of 34 nb(-1). The prompt and nonprompt psi(2S) production cross-sections and the ratio of the psi(2S) to J/psi cross-section are measured as a function of the meson transverse momentum and rapidity in the nucleon-nucleon centre-of-mass frame, together with forward-to-backward ratios and nuclear modification factors. The production of prompt psi(2S) is observed to be more suppressed compared to pp collisions than the prompt J/psi production, while the nonprompt productions have similar suppression factors.
|
|
|
|
Martin-Luna, P., Esperante, D., Prieto, A. F., Fuster-Martinez, N., Rivas, I. G., Gimeno, B., et al. (2024). Simulation of electron transport and secondary emission in a photomultiplier tube and validation. Sens. Actuator A-Phys., 365, 114859–10pp.
Abstract: The electron amplification and transport within a photomultiplier tube (PMT) has been investigated by developing an in-house Monte Carlo simulation code. The secondary electron emission in the dynodes is implemented via an effective electron model and the Modified Vaughan's model, whereas the transport is computed with the Boris leapfrog algorithm. The PMT gain, rise time and transit time have been studied as a function of supply voltage and external magnetostatic field. A good agreement with experimental measurements using a Hamamatsu R13408-100 PMT was obtained. The simulations have been conducted following different treatments of the underlying geometry: three-dimensional, two-dimensional and intermediate (2.5D). The validity of these approaches is compared. The developed framework will help in understanding the behavior of PMTs under highly intense and irregular illumination or varying external magnetic fields, as in the case of prompt gamma-ray measurements during pencil-beam proton therapy; and aid in optimizing the design of voltage dividers with behavioral circuit models.
|
|
