SEARCH · Engineering Papers

Results for “direct simulation Monte Carlo methods”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗

DOE OSTI2025

Isopotential Electron Titration of Ammonia Charge Transfer on Metal Catalysts

Electron transfer between adsorbates and surfaces determines the binding strength and reactivity of chemical moieties on materials designed for separations and catalysis. To quantify electron exchange, the extent of charge transfer resulting from ammonia adsorption on a Ru surface was measured by isopotential electron titration (IET) on a Ru catalytic condenser, where isopotential conditions were maintained between Ru and silicon separated by an insulating HfO 2 layer during gas phase ammonia adsorption. Charge transfer upon ammonia adsorption on a Ru catalytic condenser increased from 40 to 1200 nC/cm 2 at 75 and 225 °C, respectively. Charge transfer measurements provided a direct estimate of ammonia adsorption thermodynamics on Ru without knowing surface coverages a priori, revealing an adsorption enthalpy of −53 ± 10 kJ/mol and entropy of −61 ± 26 J/mol·K. Combining experimentally-measured charge transfer with kinetic Monte Carlo simulations informed adsorbate surface coverages determined that 0.058 electrons were transferred to the Ru surface for each molecular ammonia adsorption event (δ NH 3 = 0.058 ± 0.005 e – /NH 3 *), consistent with calculated Bader charges. The ability to measure the extent of charge transfer for adsorbed species provides a fundamental descriptor to understand existing and new chemically functional surfaces, providing a foundational method for the emerging field of thermochemical surface coulometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Role of hot electrons in shock ignition constrained by experiment at the National Ignition Facility

Shock ignition is a scheme for direct drive inertial confinement fusion that offers the potential for high gain with the current generation of laser facility; however, the benefits are thought to be dependent on the use of low adiabat implosions without laser–plasma instabilities reducing drive and generating hot electrons. A National Ignition Facility direct drive solid target experiment was used to calibrate a 3D Monte Carlo hot-electron model for 2D radiation-hydrodynamic simulations of a shock ignition implosion. The α = 2.5 adiabat implosion was calculated to suffer a 35% peak areal density decrease when the hot electron population with temperature T ℎ = 55 keV and energy E ℎ = 13 kJ was added to the simulation. Optimizing the pulse shape can recover ~1/3 of the peak areal density lost due to a change in shock timing. Here, despite the harmful impact of laser–plasma instabilities, the simulations indicate shock ignition as a viable method to improve performance and broaden the design space of near ignition high adiabat implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Analysis and Compensation of Tx and Rx IQ Imbalances in Uplink IM-OFDMA Systems

Index modulation-based orthogonal division multiple access (IM-OFDMA) has recently been proposed as a potential technique for future wireless communication systems due to its superior spectral efficiency and error performance over conventional multiple access schemes. However, its performance is still under investigation by researchers in light of a variety of scenarios and assumptions. Following this direction, in this Article, the individual and joint effects of the transmitter (Tx) and receiver (Rx) in-phase and quadrature imbalances (IQI) on the bit-error-rate (BER) performance of the uplink IM-OFDMA scheme are investigated by considering whether each user's RF front-end is identical or not. Moreover, to reduce the detrimental effect of IQI, a preamble-based estimation and compensation method is proposed for IM-OFDMA. Closed-form expressions for the average BER of IM-OFDMA are obtained by considering the physical effect of the Tx and Rx IQI. Also, using Monte Carlo simulations, the derived expressions are verified under different system configurations. Furthermore, analytical and simulation results reveal that Tx and Rx IQI cause an error floor in the BER performance of IM-OFDMA schemes. However, the proposed estimation and compensation method not only reduces the impact of IQI but also leads to a better BER performance compared to the case of no IQI, which mainly refers to the frequency diversity caused by the IQI.

42 ENGINEERING↗

DOE OSTI2020

Predicting the Phase Stability of Multicomponent High-Entropy Compounds

A generic method to estimate the relative feasibility of formation of high-entropy compounds in a single phase, directly from first principles, is developed. As a first step, the relative formation abilities of 56 multicomponent, AO, oxides were evaluated. These were constructed from five cation combinations chosen from A = {Ca, Co, Cu, Fe, Mg, Mn, Ni, Zn}. Candidates for multicomponent oxides are predicted from descriptors related to the enthalpy and configurational entropy obtained from the mixing enthalpies of two-component oxides. The utility of this approach is evaluated by comparing the predicted combinations with the experimentally realized entropy-stabilized oxide, (MgCoCuNiZn)O. In the second step, Monte Carlo simulations are utilized to investigate the phase composition and local ionic segregation as a function of temperature. This approach allows for the evaluation of potential secondary phases, thereby making realistic predictions of novel multicomponent compounds that can be synthesized.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Pairing Correlations across the Superfluid Phase Transition in the Unitary Fermi Gas

In the two-component Fermi gas with a contact interaction, a pseudogap regime can exist at temperatures between the superfluid critical temperature $T_c$ and a temperature $T^\ast > T_c$. This regime is characterized by pairing correlations without superfluidity. However, in the unitary limit of infinite scattering length, the existence of this regime is still debated. To help address this, we have applied finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods to study the thermodynamics of the superfluid phase transition and signatures of the pseudogap in the spin-balanced homogeneous unitary Fermi gas. We present results at finite filling factor $\textit{ν}$ ≃ 0.06 for the condensate fraction, an energy-staggering pairing gap, the spin susceptibility, and the heat capacity, and compare them to experimental data when available. In contrast to previous AFMC simulations, our model space consists of the complete first Brillouin zone of the lattice, and our calculations are performed in the canonical ensemble of fixed particle number. The canonical ensemble AFMC framework enables the calculation of a model-independent gap, providing direct information on pairing correlations without the need for numerical analytic continuation. We use finite-size scaling to estimate $T_c$ at the corresponding filling factor. In this work, we find that the energy-staggering pairing gap vanishes above $T_c$, showing no pseudogap effects, and that the spin susceptibility shows a substantially reduced signature of a spin gap compared to previously reported AFMC simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Examination of synthetic gas puff imaging diagnostic data from a gyrokinetic turbulence code

A synthetic gas puff imaging (GPI) diagnostic has been developed for the purpose of validating the three-dimensional gyrokinetic turbulence code XGC. The synthetic diagnostic is described and applied to XGC simulations of two Alcator C-Mod discharges. The turbulence characteristics deduced from the resulting simulated GPI frames, using analysis techniques similar to those applied to experimental data, are compared with locally derived characteristics extracted directly from the XGC output. The comparison of the two is shown to be potentially impacted by misalignment between the GPI view and the magnetic field, the dependence of the light emission on the electron density and temperature, and spatial and temporal variations in the neutral gas cloud induced by the turbulent plasma fluctuations. We conclude that quantitative and, in some cases even qualitative, validation of turbulence simulations need to account for these effects. While we cannot directly compare our results with experimental data due to the absence of high quality GPI data from the shots and times simulated by XGC, we do relate the overall characteristics obtained from the synthetic GPI analysis to published Alcator C-Mod GPI data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Robust verification of stochastic simulation codes

We introduce a robust verification tool for computational codes, which we call Stochastic Robust Extrapolation based Error Quantification (StREEQ). Unlike the prevalent Grid Convergence Index (GCI) [1] method, our approach is suitable for both stochastic and deterministic computational codes and is generalizable to any number of discretization variables. Building on ideas introduced in the Robust Verification [2] approach, we estimate the converged solution and orders of convergence with uncertainty using multiple fits of a discretization error model. In contrast to Robust Verification, we perform these fits to many bootstrap samples yielding a larger set of predictions with smoother statistics. Here, bootstrap resampling is performed on the lack-of-fit errors for deterministic code responses, and directly on the noisy data set for stochastic responses. This approach lends a degree of robustness to the overall results, capable of yielding precise verification results for sufficiently resolved data sets, and appropriately expanding the uncertainty when the data set does not support a precise result. For stochastic responses, a credibility assessment is also performed to give the analyst an indication of the trustworthiness of the results. Furthermore, this approach is suitable for both code and solution verification, and is particularly useful for solution verification of high-consequence simulations..

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Rotational symmetry relation for efficient response function generation in the coarse mesh transport method COMET

The coarse mesh transport code COMET is a continuous energy hybrid stochastic-deterministic neutronics solver with high fidelity and formidable computation speed in solving reactor core problems. Its method is based on the incident flux expansion theory. In this work, we take advantage of the local geometric symmetry in many reactor cores lattices (e.g., fuel lattices and reflector blocks) to develop relations among the flux response expansion coefficients for symmetric surfaces to further improve the computational efficiency of the COMET response function generation tool (method). This is done by a rigorous derivation of the transformation matrices for the angular and spatial expansion moments resulting from a rotation of a coarse mesh by an arbitrary angle. The relations for the response coefficients for the symmetric surfaces can be then written as the Kronecker product of those transformation matrices. The method is implemented into COMET and tested on two advanced high temperature reactor (AHTR) full-length single assembly benchmark problems. The COMET results using the response function library based on the symmetry relations were compared to those using the library directly generated by continuous energy Monte Carlo for all surfaces. It was found that the eigenvalues and stripe-wise fission densities using the two libraries are in statistical agreement as expected. This indicates that the new method maintains the high fidelity of the original COMET method while improving the computational efficiency in the response function generation by 270% to 400%, depending on the local geometric symmetry. This method also reduces the size of the response function library by the same magnitude (270% to 400%). (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

An intermediate model for fitting triplet–triplet annihilation in phosphorescent organic light emitting diode materials

Triplet–triplet annihilation (TTA) is one of the primary contributors to efficiency roll-off and permanent material degradation in phosphorescent organic light-emitting diodes. Here, the two limiting case models typically used to quantify this quenching mechanism are multi-step Dexter and single-step Förster, which, respectively, assume ideal Fickian diffusion or perfect trapping of triplet excitons. For device-relevant guest doping levels (typically 5–12 vol. %), both significant diffusion of excitons and trapping due to spatial and energetic disorder exist, so neither conventional model fits experimental data well. We develop and validate an intermediate TTA model, which is a weighted average of the limiting cases of pure radiative decay (no TTA) and multi-step Dexter based TTA that returns an effective TTA rate constant and a parameter quantifying the portion of well-isolated excitons. Kinetic Monte-Carlo simulations and time-resolved photoluminescence measurements of an archetype host–guest system demonstrate that our intermediate model provides significantly improved fits with more realistic physical values, is more robust to variations in experimental conditions, and provides an analysis framework for the effects of trapping on TTA.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2025

Deciphering the Scattering of Mechanically Driven Polymers Using Deep Learning

Here, we present a deep learning approach for analyzing two-dimensional scattering data of semiflexible polymers under external forces. In our framework, scattering functions are compressed into a three-dimensional latent space using a Variational Autoencoder (VAE), and two converter networks establish a bidirectional mapping between the polymer parameters (bending modulus, stretching force, and steady shear) and the scattering functions. The training data are generated using off-lattice Monte Carlo simulations to avoid the orientational bias inherent in lattice models, ensuring robust sampling of polymer conformations. The feasibility of this bidirectional mapping is demonstrated by the organized distribution of polymer parameters in the latent space. By integrating the converter networks with the VAE, we obtain a generator that produces scattering functions from given polymer parameters and an inferrer that directly extracts polymer parameters from scattering data. While the generator can be utilized in a traditional least-squares fitting procedure, the inferrer produces comparable results in a single pass and operates 3 orders of magnitude faster. This approach offers a scalable automated tool for polymer scattering analysis and provides a promising foundation for extending the method to other scattering models, experimental validation, and the study of time-dependent scattering data.

Ding, Lijie [Oak Ridge National Laboratory (ORNL),↗

DOE OSTI2022

Spectral Radiative Properties of Polydispersed SiO 2 Particle Beds

The focus of this work is on the measurement and analysis of the radiative properties of polycrystalline SiO 2 particle beds with various layer thicknesses. The particles are polydispersed with average diameters of 222, 150, and 40 μm . The spectral, directional–hemispherical reflectance and transmittance of the particle bed are measured at wavelengths from 0.4 to 1.8 μm using a monochromator, and the reflectance measurement is extended to 15 μm using a Fourier-transform infrared spectrometer. Particles are closely packed between two transparent windows for measuring the radiative properties. In the visible and near-infrared region up to 1.8 μm, the inverse adding–doubling method yields the effective absorption and scattering coefficients. The results suggest that short wavelength absorption needs to be included in modeling the behavior of particle beds due to multiple scattering. A discrete-scale Monte Carlo ray-tracing method is developed to model the radiative properties by assuming monodispersed spherical particles, and the simulated results compare well with measurements. The effective absorption and scattering coefficients of the particle beds obtained from the independent scattering theory are compared to those from the inverse method. As a result, the impact of dependent scattering on the packed beds is observed for smaller-sized particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Nonthermal electron and ion acceleration by magnetic reconnection in large laser-driven plasmas

Magnetic reconnection is a fundamental plasma process that is thought to play a key role in the production of nonthermal particles associated with explosive phenomena in space physics and astrophysics. Experiments at high-energy-density facilities are starting to probe the microphysics of reconnection at high Lundquist numbers and large system sizes. We have performed particle-in-cell (PIC) simulations to explore particle acceleration for parameters relevant to laser-driven reconnection experiments. We study particle acceleration in large system sizes that may be produced soon with the most energetic laser drivers available, such as at the National Ignition Facility. In these conditions, we show the possibility of reaching the multi-plasmoid regime, where plasmoid acceleration becomes dominant. Our results show the transition from X point to plasmoid-dominated acceleration associated with the merging and contraction of plasmoids that further extend the maximum energy of the power-law tail of the particle distribution for electrons. We also find for the first time a system-size-dependent emergence of nonthermal ion acceleration in driven reconnection, where the magnetization of ions at sufficiently large sizes allows them to be contained by the magnetic field and energized by direct X point acceleration. For feasible experimental conditions, electrons and ions can attain energies of ϵ max , e / k B T e > 100 and ϵ max , i / k B T i > 1000 . Using PIC simulations with binary Monte Carlo Coulomb collisions, we study the impact of collisionality on plasmoid formation and particle acceleration. The implications of these results for understanding the role reconnection plays in accelerating particles in space physics and astrophysics are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Lattice QED in an external magnetic field: Evidence for dynamical chiral symmetry breaking

We simulate quantum electrodynamics (QED) in a strong constant homogeneous external magnetic field on a Euclidean space-time lattice using the rational hybrid Monte Carlo method, developed for simulating lattice quantum chromodynamics (QCD). Our primary goal is to measure the chiral condensate in the limit when the input electron mass m is zero. We observe a nonzero value, indicating that the external magnetic field catalyzes chiral symmetry breaking as predicted by approximate truncated Schwinger-Dyson methods. Such behavior is associated with dominance by the lowest Landau level which causes the effective dimensional reduction from 3 + 1 dimensions to 1 + 1 dimensions for charged particles (electrons and positrons) where the attractive forces of QED can produce chiral symmetry breaking with a dynamical electron mass and associated chiral condensate. Since our lattice simulations use bare (lattice) parameters, while the Schwinger-Dyson analyses work with renormalized quantities, direct numerical comparison will require renormalization of our lattice results. Published by the American Physical Society 2024

Kogut, J. B. (ORCID:000000030151548X)↗

DOE OSTI2024

Determination of the mean energy of fast electron losses and anisotropies through thick-target emission on WEST

Abstract A new method to obtain the mean energy of fast electron losses in fusion plasmas using a versatile multi-energy hard x-ray (HXR) detector is presented. The method is based on measuring the thick-target emission of tungsten in the divertor region produced by fast electron losses interacting with the target and modeling the tungsten spectra by a Monte Carlo code which simulates the interaction between a beam of electrons and a solid target. The mean energy of the fast electron losses is determined through the comparison between the experimental and synthetic emission. The results show that fast electron losses during lower hybrid current drive discharges at WEST have a mean energy of 90–140 keV and represent only 2% of the total heat flux at the target. Additionally, anisotropic HXR emission has been detected for the first time at the WEST core and edge plasma, with opposite directions. It is due to the forward-peak emission of two distinctive populations of fast electrons: co-current fast electrons in the core and counter-current fast electron losses at the inner strike point. In view of future experiments like ITER where electron cyclotron current drive will generate a fast electron population, this technique could serve as a real-time monitor of fast electron losses and eventually feed an actuator on the current drive generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Machine Learning Assisted Safety Modeling and Analysis of Advanced Reactors

With the advances in computational power and numerical methods, analysts can now rely on first-principle simulations to predict ultra-fine details in a variety of applications. Advances in machine learning (ML) have produced algorithms that can now learn high-level abstractions via hierarchical models. This project aims to leverage advances in ML techniques and the available high-resolution simulation data to develop a novel modeling and simulation (M\&S) methodology for reactor safety analysis. While application-agnostic ML techniques are available, complex physics constraints need to be incorporated into ML techniques to build ML-based closures for computationally efficient predictive simulations. This project intends to develop a physics-guided data-driven multi-scale methodology for M\&S of advanced reactors. The project focuses on thermal fluid (T/F) phenomena, which play major roles in advanced reactor safety. Specifically, we propose a data-driven coarse-mesh turbulence model based on local flow features for the transient analysis of thermal mixing and stratification in a sodium-cooled fast reactor (SFR). The model has a coarse-mesh setup to ensure computational efficiency, while it is trained by fine-mesh computational fluid dynamics (CFD) data with Reynolds-averaged Navier-Stokes (RANS) turbulence model to ensure accuracy. Three different neural networks are developed and tested for loss-of-flow transients in the hot pool of SFR, i.e. the densely connected convolutional neural network (DCNN), long-short-term-memory network based on proper orthogonal decomposition (POD-LSTM), and the DCNN informed by LSTM (DCNN-LSTM). The performances of these three neural networks are evaluated based on baseline models. The DCNN-LSTM model has been chosen for further hyperparameter optimization. Furthermore, based on a simplified two-dimensional case, uncertainty quantification (UQ) of the developed ML-based closure are investigated with three methods, i.e. Monte Carlo dropout, deep ensemble, and Bayesian neural network. The developed ML-based turbulent viscosity closure relation based on deep ensemble is then integrated into the system analysis module SAM and serves as a term in the conservation equations. Such a SAM-ML based procedure guarantees that the obtained results are consistent with the physical constraints of the thermal-fluid system. The SAM-ML simulation on the same loss-of-flow transient showed comparable accuracy with the CFD simulation but with a much coarser mesh setup. Last but not least, the ML-based closure improvement with the support of higher-fidelity data from large eddy simulation (LES) is discussed. As a first step towards this direction, a baseline LES simulation is performed to obtain comparable data with RANS results. Based on the early results, future investigation on further improving the ML-based closure is discussed. We believe the developed approach that combines scientific machine learning with nuclear system analysis code can benefit the advanced reactor community as more accurate safety analyses will better characterize reactor safety margins and reduce licensing efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Newton trust-region methods with primary variable switching for simulating high temperature multiphase porous media flow

Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Monte Carlo simulation of the passage of $\gamma$-rays and $\alpha$-particles in CsI

In this work, theoretical and computational methods for simulating the creation of ionization tracks by fast ions in solids were applied to the passage of a-particles in CsI, an inorganic scintillator commonly used for radiation detection. The methods were implemented in a Monte Carlo program to simulate the interaction of $\alpha$-particles, with incident energies of up to 1 MeV, with CsI. The simulations followed the fate of individual electron-hole pairs and thus allowed for a detailed description of the microscopic structure of ionization tracks created by incident radiation. Simulations were also performed with $\gamma$-rays of the same energy to compare and contrast the ionization tracks obtained with both types of particle. Intrinsic properties such as the mean energy per electron-hole pair, Fano factor, maximum theoretical light yield, and spatial distributions of electron-hole pairs were computed for both $\alpha$-particles and $\gamma$-rays. $\alpha$-particles created cylindrical tracks that were initially aligned with the incident direction and with initial radii of a few nanometers, whereas $\gamma$-rays showed significant scattering, resulting in probability distributions with lower intensities and much greater radial extents.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗