Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geant4”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

A pion-argon cross section measurement in the ProtoDUNE-SP experiment with cosmogenic muon

Neutrinos are tiny mysterious fundamental particles with small cross sections. Through neutrino physics, scientists across the world are trying to answer many intriguing questions about nature such as the dominance of matter over antimatter, CP violation in the lepton sector, number of supernovas in the early universe, etc. Detection of neutrinos requires massive particle detectors and intense neutrino beam owing to their small cross section. Deep Underground Neutrino Experiment (DUNE) is a next-generation neutrino experiment that is planned to start taking data beginning in 2026. DUNE will consist of 4 massive detectors, the first of which will be using single-phase liquid argon time projection chamber (LArTPC) technology. The ProtoDUNE-SP experiment is a prototype of the DUNE built at the CERN neutrino platform and uses the same detector technology that will be used in DUNE first module. The ProtoDUNE-SP experiment collected months of test beam and cosmic ray data beginning in September 2018. It was built to provide a testbed for the installation of detector parts for DUNE, showing long-term stability of the detector, understanding detector response for different test beam particles (including protons, pions, electrons, kaons, muons), and measurement of hadron-argon cross sections. When a particle passes through LArTPC electron-ion pairs are produced. To reconstruct the position and energy of a particle passing through the medium knowledge of ionization electron drift velocity is essential. The electron drift velocity is distorted by an excess positive charge built up in the detector, known as space charge. This study discusses a novel technique for measuring the ionization electron drift velocity using cosmic-ray muons. The technique uses tracks that travel the entire drift distance of the TPC for drift velocity determination. Secondly, the study discusses a method for converting the charge deposited into energy. The method is carried out in two step s. In th e first step detector response for energetic cosmic ray muons crossing the entire the TPC is used to make the charge deposition uniform throughout the TPC, and in the second step stopping cosmic-ray muons are used for determining the energy scale. Finally, the study discusses a pion-argon cross section measurement based on reweighting of Monte Carlo simulations using J. Calcutt's Geant4Reweight framework. Neutrinos cannot be directly detected; they are identified based on the interaction products. Pions are a common interaction product in a neutrino interaction. For precise modeling of neutrino event generators, it is essential to understand the pion-argon interaction. Pion-argon cross section measurement serves as an important input for neutrino interaction models. The results of the pion-argon total reaction cross section using the Geant4 reweighting technique are found to be in good agreement with Geant4 predictions. The many studies carried out in the ProtoDUNE-SP experi ment wil l be useful for current and future neutrino experiments using LArTPC technology including ICARUS, MicroBooNE, DUNE

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Energy Characterization of the UMD Electron Accelerator

Multiple different groups at Los Alamos National Laboratory have a cited interest in understanding the effects of ionizing radiation in materials. Of especial interest to the ISR division is understanding the effects in materials that are used in the space environment, where radiation is omnipresent and impossible to fully shield from. Development is under way to simulate the effects of radiation damage in detectors and electronics, primarily with particle transport codes such as Geant4 and RAMPART, a code that provides a user friendly interface to Geant4. However, these simulations need benchmarking to experimental tests, some of which have been performed at a variety of small electron accelerator facilities. One of these facilities is the Linear Electron Accelerator at the University of Maryland Radiation Facilities. This facility can provide electron beams on the order of 100 mA delivered in tight macropulse bunches, which allows for experiments where the behavior of electronics is monitored as a function of pulse number. While the beam rates per macropulse of this accelerator is monitored during an experiment via a current pick-off connected to a large graphite beamstop, the energy distribution of this electron beam is not measured during the experiment and has not been well characterized. The energy distribution of the beam is pertinent for these types of experiments, since dose is directly related to energy deposited in a medium. Further, these experiments are performed in open air and sample placement in the room can dictate what dose rate they receive per macropulse. We have created simulated dosemaps for the UMD vault room, but these simulations assume a given energy of the beam, both for calculating the dose and to simulate how much the beam diverges in the air. There have been two attempts now to characterize this accelerator beam energy with a plate type spectrometer, including the most recent experiment earlier this year to quantify the beam energy distribution. These results will also be compared to a dosimetry experiment performed by collaborators at the UMD facility, the details of which are described in Reference. Understanding the electron beam energy is important for future experiments, because the energy directly impacts the dose and radiation damage given to a sample and is important for properly simulating the electron beam when comparing simulation to experiment.

43 PARTICLE ACCELERATORS↗

Classification Algorithm for Dark Matter Search using Skipper-CCD

In this project, I used ROOT framework with C++ for data analysis and Geant4 for simulating different particles’ interaction with the detector to study the characteristic of their tracks. I prepared the tagged dataset for machine learning training both the real experiment data taken in MINOS and simulation data from Geant4.

Chen, Meng-Wei↗

In-memory Geometry Converter

The slides describe the Geant4-to-VecGeom in-memory geometry converter, including some details on its development and validation. (To be presented at the Geometry Session of the 2023 Geant4 Collaboration Meeting, in Hokkaido, Japan, Sept. 25-29, 2023)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Monte Carlo Simulation with CAD Interface for Calculation of 3D Maps of Residual Dose (CRADA)

Objective: To develop an easy-to-use software application to predict and mitigate radiation effects in research environment, space instruments, nuclear plants and medical facilities and help nonproliferation and national security efforts. Tech-X will develop standalone software libraries and command-line tools for ( 1) translating CAD into tessellated surfaces and tetrahedral meshes in GDML (for Geant4 and MARS 15), ROOT (for MARS 15) and HDF5 (for compact representation and for the visualization) formats, (2) healing CAD geometries to make them suitable for Monte Carlo simulations; (3) creating uniform and variable Cartesian and cylindrical meshes for detailed scoring; and ( 4) efficient Monte Carlo navigation in CAD geometries. JLAB will finish automation of simulations of residual dose in CAD geometries and integrate Tech-X software into Geant4 and MARS15. Finally, Tech-X will develop a Graphical User Interface to set up and heal CAD geometries, create input files, run and visualize simulations for residual dose. This application will run on local desktops, local and remote clusters and supercomputers and will be made available through public clouds, such as Amazon Web Services.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Current State of Neutron Source Term Calculations

SOURCES-4C has some significant deficiencies that are mostly due to problems with the default data libraries that it uses. Some research groups have modified these libraries to fix most of the problems, but not all of them. Adding modern data libraries to SOURCES-4C is one of the most promising approaches to resolving these issues. Other promising options include using either NEDIS or a GEANT4 tool. Both of these codes approximately match the accuracy of the modified SOURCES-4C code. The NEDIS code is Russian. The GEANT4 tool might be computationally slow and it might also no longer be actively maintained. There are a handful of other (α,n) codes, but they are not accurate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Heterogeneous Computing

To leverage the increasing heterogeneity in modern computing resources, Geant4 incorporates advanced software tools and a task-based framework (G4Tasking) that enables efficient parallelism at event, sub-event, and track levels. Ongoing R&D efforts focus on integrating GPUs into high-energy physics (HEP) simulations, including optical photon simulation with Opticks/NVIDIA OptiX, offloading electromagnetic particle transport using G4HepEM/AdePT and Celeritas, and employing advanced surface-based geometry models such as VecGeom2.0 and ORANGE. As Geant4 continues evolving toward high-performance computing (HPC) and heterogeneous architectures, it remains a key tool for large-scale simulations in HEP and beyond.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantification of LEU Holdup using gamma ray imaging and inverse transport solver

Holdup is the residual amount of special nuclear material (SNM) remaining in a processing facility after the bulk materials have been cleaned out. In commercial uranium processing facilities, quantification of holdup is a major challenge because of the highly variable shapes and sizes of the deposits. Any method that attempts to generalize and calibrate deposit shapes in order to quantify holdup will be prone to high uncertainties. Uncertainties on the order of ±50% are typical in holdup results. In international safeguards applications, a ±50% uncertainty can result in a large amount of material unaccounted for (MUF) thereby increasing the difficulty of detecting material diversion and facility misuse. An imaging-based methodology has been developed with the objective of significantly reducing this uncertainty by using the true deposit shape, instead of relying on oversimplified geometric assumptions. The project is a collaboration between ORNL, Y-12, and the University of Tennessee, Knoxville, TN. Uranium sources of known masses were measured using the Germanium Gamma-ray Imager (GeGI), a high-resolution imaging spectrometer, creating a pixelated map for each spectral bin. Two different gamma imaging methods are employed in this work: coded aperture imaging and Compton imaging. A validated MonteCarlo model of the detector has been developed using the GEANT4 code for determining the intrinsic response of the detector, its enclosure, and the coded aperture mask. An inverse transport solver based on the Markov Chain Monte-Carlo approach known as Differential Evolution Adaptive Metropolis (DREAM) is employed to use the measurement data from the image pixels (coded aperture or Compton) to solve for the mass of 235 U in the deposit. A reliable method based on the DREAM solver has been developed to flag the infinite thickness condition of a uranium deposit. The project team is working towards improving the image reconstruction for Compton imaging so that a better localization of the source can be achieved. Besides treating the coded aperture and Compton imaging methods independently, the project is also evaluating a combined method that uses the Compton scatter data from a coded aperture measurement. GEANT4 simulations are being performed to evaluate the combined approach. The impact on the DREAM optimization as the source thickness progressively approaches infinite thickness is being evaluated. A number of uranium sources available at ORNL have been measured, and the DREAM results have been tested and validated for the coded aperture imaging. A similar effort will be carried out to validate the Compton based method once the development of algorithms for better localization are complete. The imaging based quantification is very amenable to unattended monitoring of holdup accumulation at key measurement points. A proof of concept measurement has been completed to demonstrate this capability The current work used the high energy resolution imager GeGI. However, the approach and methodologies are applicable to other imagers such as the cadmium zin telluride (CZT) based imager manufactured by H3D, Inc.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of a Space Radiation Monte-Carlo Computer Simulation Based on the FLUKE and Root Codes

The radiation environment in space is a complex problem to model. Trying to extrapolate the projections of that environment into all areas of the internal spacecraft geometry is even more daunting. With the support of our CERN colleagues, our research group in Houston is embarking on a project to develop a radiation transport tool that is tailored to the problem of taking the external radiation flux incident on any particular spacecraft and simulating the evolution of that flux through a geometrically accurate model of the spacecraft material. The output will be a prediction of the detailed nature of the resulting internal radiation environment within the spacecraft as well as its secondary albedo. Beyond doing the physics transport of the incident flux, the software tool we are developing will provide a self-contained stand-alone object-oriented analysis and visualization infrastructure. It will also include a graphical user interface and a set of input tools to facilitate the simulation of space missions in terms of nominal radiation models and mission trajectory profiles. The goal of this project is to produce a code that is considerably more accurate and user-friendly than existing Monte-Carlo-based tools for the evaluation of the space radiation environment. Furthermore, the code will be an essential complement to the currently existing analytic codes in the BRYNTRN/HZETRN family for the evaluation of radiation shielding. The code will be directly applicable to the simulation of environments in low earth orbit, on the lunar surface, on planetary surfaces (including the Earth) and in the interplanetary medium such as on a transit to Mars (and even in the interstellar medium). The software will include modules whose underlying physics base can continue to be enhanced and updated for physics content, as future data become available beyond the timeframe of the initial development now foreseen. This future maintenance will be available from the authors of FLUKA as part of their continuing efforts to support the users of the FLUKA code within the particle physics community. In keeping with the spirit of developing an evolving physics code, we are planning as part of this project, to participate in the efforts to validate the core FLUKA physics in ground-based accelerator test runs. The emphasis of these test runs will be the physics of greatest interest in the simulation of the space radiation environment. Such a tool will be of great value to planners, designers and operators of future space missions, as well as for the design of the vehicles and habitats to be used on such missions. It will also be of aid to future experiments of various kinds that may be affected at some level by the ambient radiation environment, or in the analysis of hybrid experiment designs that have been discussed for space-based astronomy and astrophysics. The tool will be of value to the Life Sciences personnel involved in the prediction and measurement of radiation doses experienced by the crewmembers on such missions. In addition, the tool will be of great use to the planners of experiments to measure and evaluate the space radiation environment itself. It can likewise be useful in the analysis of safe havens, hazard migration plans, and NASA's call for new research in composites and to NASA engineers modeling the radiation exposure of electronic circuits. This code will provide an important complimentary check on the predictions of analytic codes such as BRYNTRN/HZETRN that are presently used for many similar applications, and which have shortcomings that are more easily overcome with Monte Carlo type simulations. Finally, it is acknowledged that there are similar efforts based around the use of the GEANT4 Monte-Carlo transport code currently under development at CERN. It is our intention to make our software modular and sufficiently flexible to allow the parallel use of either FLUKA or GEANT4 as the physics transport engine.

Pinsky, L. S.↗

Simulation of the ATIC-2 Silicon Matrix for Protons and Helium GCR Primaries at 0.3, 10, and 25 TeV/Nucleon

The energy deposition distribution for protons and helium galactic cosmic ray primaries at 0.3, 10, and 25 TeV/nucleon in the ATIC-2 silicon matrix detector are simulated with GEANT4. The GEANT3 geometrical model of ATIC developed by the University of Maryland was combined with a GEANT4 application developed for the Deep Space Test Bed (DSTB) detector package. The new code included relatively minor modifications to completely describe the ATIC materials and a more detailed model of the Silicon Matrix detector. For this analysis all particles were started as a unidirectional beam at a single point near the center of the Silicon Matrix front surface. The point was selected such that each primary passed through at least two of the overlapping silicon pixels.

Watts, J.↗

Depth Dose Distribution Study within a Phantom Torso after Irradiation with a Simulated Solar Particle Event at NSRL

The adequate knowledge of the radiation environment and the doses incurred during a space mission is essential for estimating an astronaut's health risk. The space radiation environment is complex and variable, and exposures inside the spacecraft and the astronaut's body are compounded by the interactions of the primary particles with the atoms of the structural materials and with the body itself Astronauts' radiation exposures are measured by means of personal dosimetry, but there remains substantial uncertainty associated with the computational extrapolation of skin dose to organ dose, which can lead to over- or underestimation of the health risk. Comparisons of models to data showed that the astronaut's Effective dose (E) can be predicted to within about a +10% accuracy using space radiation transport models for galactic cosmic rays (GCR) and trapped radiation behind shielding. However for solar particle event (SPE) with steep energy spectra and for extra-vehicular activities on the surface of the moon where only tissue shielding is present, transport models predict that there are large differences in model assumptions in projecting organ doses. Therefore experimental verification of SPE induced organ doses may be crucial for the design of lunar missions. In the research experiment "Depth dose distribution study within a phantom torso" at the NASA Space Radiation Laboratory (NSRL) at BNL, Brookhaven, USA the large 1972 SPE spectrum was simulated using seven different proton energies from 50 up to 450 MeV. A phantom torso constructed of natural bones and realistic distributions of human tissue equivalent materials, which is comparable to the torso of the MATROSHKA phantom currently on the ISS, was equipped with a comprehensive set of thermoluminescence detectors and human cells. The detectors are applied to assess the depth dose distribution and radiation transport codes (e.g. GEANT4) are used to assess the radiation field and interactions of the radiation field with the phantom torso. Lymphocyte cells are strategically embedded at selected locations at the skin and internal organs and are processed after irradiation to assess the effects of shielding on the yield of chromosome damage. The initial focus of the present experiment is to correlate biological results with physical dosimetry measurements in the phantom torso. Further on, the results of the passive dosimetry within the anthropomorphic phantoms represent the best tool to generate reliable data to benchmark computational radiation transport models in a radiation field of interest. The presentation will give first results of the physical dose distribution, the comparison with GEANT4 computer simulations based on a Voxel model of the phantom, and a comparison with the data from the chromosome aberration study.

Berger, Thomas↗

A Computationally Efficient Algorithm for Sampling the Rudd Differential Cross Section

Monte Carlo radiation transport codes such as RITRACKS or Geant4 are used to simulate the interaction of ions with matter. These codes rely on sampling algorithms to determine interactions and various physical properties of particles involved in the simulations. It is crucial to develop efficient sampling algorithms since Monte Carlo radiation transport simulations can be time consuming. This work presents an efficient sampling algorithm to determine the energy of secondary electrons following ion-water interactions. The applicability and intended use of the algorithm are discussed in detail, and it is shown that the new algorithm is up to 6X10 4 times faster than the method currently used in Geant4-DNA.

Floriane Poignant↗

Effect of Radiation from Low-Level Radioactive Waste Leachate on Antioxidant Depletion in HDPE Geo-membrane - 16087

Low-level radioactive waste (LLW) and mixed waste (MW) disposal facilities in the US are required to have a service life in excess of 1000 yr. Understanding the rate of degradation of high-density polyethylene (HDPE) geo-membranes (GMs) used in LLW and MW facilities is necessary to assess their service life when conducting a performance assessment for a disposal facility. In LLW and MW disposal facilities, GMs may be exposed to radionuclides that may accelerate antioxidant depletion and degradation. This paper evaluates how radiation dosage affects antioxidant depletion in HDPE GMs and makes inferences regarding the impact of radiation on GM service life. Specimens of 0.04-mm-, 0.1-mm-, and 2-mm-thick HDPE GMs were exposed to a sealed source of Am-241 to simulate alpha radiation from LLW leachate (e.g., uranium). The specimens were exposed to the sealed source for 1-50 h, followed by oxidative induction time (OIT) testing to determine depletion of antioxidants in response to exposure. Penetration of alpha particles in the HDPE GM exposed to the sealed source was evaluated with the program Geometry And Tracking (GEANT4). The GEANT4 modeling indicates that physical dose deposition for alpha particles largely occurs on the superficial layers and decreased monotonically with depth to a negligible dose (∼ 28 μm). For a 2-mm-thick HDPE GM, there was no change in OIT after exposure, whereas in much thinner GMs (e.g., 0.04 mm), significant antioxidant depletion occurred after exposure, reflecting the shallow penetration of alpha radiation at the surface of HDPE GM. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Neutron detector response modeling in NOvA

Neutrons can present a significant challenge for neutrino experiments in which energy reconstruction is critical. With the ability to escape detection completely and with a weak correlation between their kinetic energy and any eventual energy deposition, it is difficult to fully account for neutrons produced in neutrino interactions. This in turn leads to significant model dependence when evaluating neutron-related systematic uncertainties. The NOvA experiment is a long-baseline neutrino oscillation experiment with a high-statistics sample of antineutrino data collected by its near detector. We report an excess relative to data of simulated neutron candidates with low energy depositions when using standard Geant4 physics lists. The simulation excess is traced to an overabundance of secondary photons produced from interactions of neutrons with kinetic energy greater than \SI{20}{\mega\eV}. Improved agreement with data is obtained by applying the data-driven neutron-on-carbon \menate model for neutrons between \SI{20}{\mega\eV} and ${\sim}$\SI{100}{\mega\eV}. With \menate, the residual oversimulation is more uniform across the calorimetric neutron energy spectrum, suggesting possible overproduction of primary neutrons by the GENIE neutrino interaction generator. These results motivate the adoption of \menate-supplemented Geant4 simulation as the nominal simulation in the production of future \nova simulation.

Abubakar, S.↗

Evaluation of cosmogenic production of 39 Ar and 42 Ar for rare-event physics using underground argon

Underground argon (UAr) with lower cosmogenic activities of 39 Ar and 42 Ar has been planned as a detector in detecting scintillation light and charge collection using time projection chambers for dark matter searches and as a veto detector in suppressing backgrounds for neutrinoless double beta decay (0$vββ$ ) experiments. Long-lived radioactive isotopes, 39 Ar and 42 Ar, can also be produced on the surface when UAr is pumped out from a deep well. Understanding the production of long-lived isotopes in Ar is important for utilizing UAr for dark matter and 0$vββ$ experiments in terms of its production, transportation, and storage. Ar exposure to cosmic rays at sea-level is simulated using Geant4 for a given cosmic ray muon, neutron, and proton energy spectrum. Here we report the simulated cosmogenic production rates of 39 Ar, 42 Ar, and other long-lived isotopes at sea-level from fast neutrons, high energy muons, and high energy protons. Total production rates of 938.53/kg Ar day and 5.81 × 10 -3 /kg Ar day for 39 Ar and 42 Ar are found from our simulation. Utilizing these production rates, we set a time limit of 954 days constrained by the production of 39 Ar for UAr to be on the surface before it compromises the sensitivity for a dark matter experiment. Similarly, a time limit of 1702 days constrained by the production of 42 Ar is found for a 0$vββ$ experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comparison of modeled and measured performance of a GSO crystal as gamma detector

We have modeled, tested, and installed a large, cerium-activated Gd{sub 2}SiO{sub 5} crystal scintillator for use as a detector of gamma rays. We present the measured detector response to two types of incident photons: nearly monochromatic photons up to 40 MeV, and photons from a continuous Compton backscattering spectrum up to 200 MeV. Our GEANT4 simulations, developed to determine the analyzing power of the Compton polarimeter in Hall A of Jefferson Lab, reproduce the measured spectra well.

Geant4 Gamma detection↗

CsI(Tl) pulse shape discrimination with the Belle II electromagnetic calorimeter as a novel method to improve particle identification at electron–positron colliders

Here we describe the implementation and performance of CsI(Tl) pulse shape discrimination for the Belle II electromagnetic calorimeter, representing the first application of CsI(Tl) pulse shape discrimination for particle identification at an electron–positron collider. The pulse shape characterization algorithms applied by the Belle II calorimeter are described. Control samples of $γ, μ^+, π^±, K^±$ and $ρ/ \bar{ρ}$ are used to demonstrate the significant insight into the secondary particle composition of calorimeter clusters that is provided by CsI(Tl) pulse shape discrimination. Comparisons with simulation are presented and provide further validation for newly developed CsI(Tl) scintillation response simulation techniques, which when incorporated with GEANT4 simulations allow the particle dependent scintillation response of CsI(Tl) to be modelled. Comparisons between data and simulation also demonstrate that pulse shape discrimination can be a new tool to identify sources of improvement in the simulation of hadronic interactions in materials. The $K^0_L$ efficiency and photon-as-hadron fake-rate of a multivariate classifier that is trained to use pulse shape discrimination is presented and comparisons are made to a shower-shape based approach. CsI(Tl) pulse shape discrimination is shown to reduce the photon-as-hadron fake-rate by over a factor of 3 at photon energies of 0.2 GeV and over a factor 10 at photon energies of 1 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Benchmarking fast-neutron response of 7 Li-enriched C 7 LYC scintillators: Simulations and experiment

New insights into the fast neutron response of 7 Li-enriched C 7 LYC scintillators have been obtained through detailed Geant4 simulations of 35 Cl(n,p) and 35 Cl(n, α) reactions that are primarily responsible for the fast neutron response of C 7 LYC. The light output of a 1” x1” C 7 LYC detector was simulated using the Birks equation for light quenching in scintillators. The simulations were benchmarked against neutron scattering data from a 12 C target obtained with a white neutron beam at the Los Alamos LANSCE facility and fourteen 1”x1” C 7 LYC detectors. Incident neutron energies extracted from time-of-flight measurements, and appropriate recoil corrections for the different detector angles were used to map the response of C 7 LYC to mono-energetic neutrons. Resonances in the 35 Cl(n,p) reaction cross-section below 3 MeV incident neutron energy were directly observed in C 7 LYC. The experimental data and simulations are compared with ENDF databases, as well as recent direct cross-section measurements of neutrons scattered from a 35 Cl target. In conclusion, the results suggest a scrutiny of the 35 Cl(n,p) and 35 Cl(n, α) cross sections at low energy in the ENDF library.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗