Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tunnel Testing”

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 199 records · Page 11

Gravitational wave spectrum of chain inflation

Chain inflation is an alternative to slow-roll inflation in which the inflaton tunnels along a large number of consecutive minima in its potential. In this work we perform the first comprehensive calculation of the gravitational wave (GW) spectrum of chain inflation. In contrast to slow-roll inflation the latter does not stem from quantum fluctuations of the gravitational field during inflation, but rather from the bubble collisions during the first-order phase transitions associated with vacuum tunneling. Our calculation is performed within an effective theory of chain inflation which builds on an expansion of the tunneling rate capturing most of the available model space. The effective theory can be seen as chain inflation’s analog of the slow-roll expansion in rolling models of inflation. The near scale-invariance of the scalar power spectrum translates to a quasiperiodic shape of the inflaton potential in chain inflation, with the tunneling rate changing very slowly during the e-folds leading to cosmic microwave background observables. We show that chain inflation produces a very characteristic double-peak GW spectrum: a faint high-frequency peak associated with the gravitational radiation emitted during inflation, and a strong low-frequency peak associated with the graceful exit from chain inflation (marking the transition to the radiation-dominated epoch). There exist very exciting prospects to test the gravitational wave signal from chain inflation at the aLIGO-aVIRGO-KAGRA network, at LISA and /or at pulsar timing array experiments. A particularly intriguing possibility we point out is that chain inflation could be the source of the stochastic gravitational wave background recently detected by NANOGrav, PPTA, EPTA, and CPTA. We also show that the gravitational wave signal of chain inflation is often accompanied by running/ higher running of the scalar spectral index to be tested at future cosmic microwave background experiments. Published by the American Physical Society 2024

Freese, Katherine↗

Towards replacing physical testing of granular materials with a Topology-based Model

In the study of packed granular materials, the performance of a sample (e.g., the detonation of a high-energy explosive) often correlates to measurements of a fluid flowing through it. The “effective surface area,” the surface area accessible to the airflow, is typically measured using a permeametry apparatus that relates the flow conductance to the permeable surface area via the Carman-Kozeny equation. This equation allows calculating the flow rate of a fluid flowing through the granules packed in the sample for a given pressure drop. However, Carman-Kozeny makes inherent assumptions about tunnel shapes and flow paths that may not accurately hold in situations where the particles possess a wide distribution in shapes, sizes, and aspect ratios, as is true with many powdered systems of technological and commercial interest. To address this challenge, we replicate these measurements virtually on micro-CT images of the powdered material, introducing a new Pore Network Model based on the skeleton of the Morse-Smale complex. Pores are identified as basins of the complex, their incidence encodes adjacency, and the conductivity of the capillary between them is computed from the cross-section at their interface. We build and solve a resistive network to compute an approximate laminar fluid flow through the pore structure. Here, we provide two means of estimating flow-permeable surface area: (i) by direct computation of conductivity, and (ii) by identifying dead-ends in the flow coupled with isosurface extraction and the application of the Carman-Kozeny equation, with the aim of establishing consistency over a range of particle shapes, sizes, porosity levels, and void distribution patterns.

36 MATERIALS SCIENCE↗

Perovskite Photovoltaic Devices with Carbon-Based Electrodes Withstanding Reverse-Bias Voltages up to –9 V and Surpassing IEC 61215:2016 International Standard

One of the key challenges of perovskite photovoltaics (PV) is the long-term stability. Although efforts are made to improve the lifetime of perovskite PV devices, their degradation under reverse-bias conditions is barely addressed. Herein, perovskite solar cells with carbon-based electrodes are presented which demonstrate superior resilience against reverse-bias-induced degradation. Although their breakdown voltage is identified to be at approximately -3.6 V, cells do not degrade until the applied reverse-bias exceeds -9 V. Two main degradation mechanisms are identified: 1) iodine loss due to hole tunneling into perovskite, which takes place even at low reverse-bias but decomposes the perovskite only after long time durations; and 2) rapid heating at large reverse-bias leading to formation of PbI2, which starts at shunts and then follows the path of the least resistance for the cell current, which is primarily influenced by the electrode sheet resistances. Finally, perovskite solar modules with carbon-based electrodes are demonstrated, which are subjected to a "hotspot" test described in the IEC 61215:2016 international standard at an accredited module testing laboratory. Passing this accelerated test for the first time confirms the superior stability of perovskite PV devices with carbon-based electrodes and highlights their large industrialization potential.

14 SOLAR ENERGY↗

A multi-Physics Experiment for Low-Yield Nuclear Explosion Monitoring

A series of multi-physics experiments, referred to as Physics Experiment 1 (PE1) is underway at the United States’ Nevada National Security Site (NNSS). The PE1 series includes detonations of three underground chemical explosions in P-tunnel, with fully coupled (PE1 A), partially decoupled (PE1 D L ), and fully decoupled (PE1 B) emplacements. Canisters with gas tracers are imbedded in the explosives, and the tracers are released when the canister is destroyed by the detonation. A dedicated electromagnetic (EM) experiment (EMX) generates well-characterized EM signals at an underground location near the chemical explosive experiments. A series of atmospheric experiments (METEX, REACT, and METREX) release smoke and radioactive tracers around Aqueduct Mesa to test gas transport in complex topography. Each of the chemical explosive experiments includes a network of sensors to record seismic, acoustic, and electromagnetic waves, measurement of atmospheric conditions, and air sample collection for measurement of tracer concentration. EMX records EM signals underground and on the surface of Aqueduct Mesa. METEX, REACT, and METREX include measurement of atmospheric condition, as well as tracking smoke releases. REACT and METREX add low-level radioactive gas tracers to the atmospheric releases.

58 GEOSCIENCES↗

Study the existing form of copper ( p-type oxide/segregation ) and its release mechanism from the passive film of Ti-7Cu alloy

For this work, the existing form of copper and its release mechanism from the passive film of Ti-7Cu alloy was studied by Quantitative analysis of depth profile and series of electrochemical tests. Results were interpreted within the Point Defect Model-II. Copper does not form any p-type (cuprous or cupric) oxide but preferentially oxidises at the metal/film interface and segregated as Cu interstitials ($Cu^+_i$) and become substituted as $Cu^{x'}_{Ti}$ in the oxygen-deficient TiO 2-x barrier layer having n-type semiconducting character. Eventually, the segregated species of Cu preferentially dissolute through the film as Cu ions into the solution at film/solution interface via the Esaki-tunnelling effect.

Siddiqui, Muhammad Ali↗

Next Generation Ta-STJ Sensor Arrays for BSM Physics Searches

The Beryllium Electron capture in Superconducting Tunnel junctions (BeEST) experiment uses superconducting tunnel junction (STJ) sensors to search for physics beyond the standard model (BSM) with recoil spectroscopy of the 7 Be EC decay into 7 Li. A pulsed UV laser is used to calibrate the STJs throughout the experiment with ∼20 meV precision. Phase-III of the BeEST experiment revealed a systematic calibration discrepancy between STJs. We found these artifacts to be caused by resistive crosstalk and by intensity variations of the calibration laser. For phase-IV of the BeEST experiment, we have removed the crosstalk by designing the STJ array so that each pixel has its own ground wire. We now also use a more stable UV laser for calibration. The new STJ arrays were fabricated at STAR Cryoelectronics and tested at LLNL and FRIB. They have the same high energy resolution of ∼1-2~eV in the energy range of interest below 100~eV as before, and they no longer exhibit the earlier calibration artifacts. We discuss the design changes and the STJ array performance for the next phase of the BeEST experiment.

Laser Calibration↗

Assessing Options for Particulate Measurement Tools to Advance Wood Heaters (CRADA Final Report)

For approximately 40 years, particulate matter (PM) from residential wood heaters has been certified in North America using an EPA dilution tunnel with a filter-based gravimetric measurement. While reasonably repeatable, this approach reports only a single value for the entire burn sequence and does not provide information on PM at different points within the sequence, which would support improved wood heater designs. PM monitors that report with high time resolution offer a potential solution, but their performance for various types of heaters and across heater operating conditions has not been extensively evaluated. In this 18-month collaborative project, the Hearth, Patio & Barbecue Association (HPBA) and Lawrence Berkeley National Laboratory (LBNL) tested five continuous PM instruments spanning three measurement principles—light-scattering, inertial mass (TEOM), and triboelectric—against gravimetric and size-resolved references on three residential heater types (cordwood, catalytic cordwood, and pellet) across startup, high-burn, and low-burn operating phases.

42 ENGINEERING↗

Analytical nonadiabatic coupling and state-specific energy gradient for the crystal field Hamiltonian describing lanthanide single-ion magnets

Paramagnetic molecules with a metal ion as an electron spin center are promising building blocks for molecular qubits and high-density memory arrays. However, fast spin relaxation and decoherence in these molecules lead to a rapid loss of magnetization and quantum information. Nonadiabatic coupling (NAC), closely related to spin-vibrational coupling, is the main source of spin relaxation and decoherence in paramagnetic molecules at higher temperatures. Predicting these couplings using numerical differentiation requires a large number of computationally intensive ab initio or crystal field electronic structure calculations. To reduce computational cost and improve accuracy, we derive and implement analytical NAC and state-specific energy gradient for the ab initio parametrized crystal field Hamiltonian describing single-ion molecular magnets. Our implementation requires only a single crystal field calculation. In addition, the accurate NACs and state-specific energy gradients can be used to model spin relaxation using sophisticated nonadiabatic molecular dynamics, which avoids the harmonic approximation for molecular vibrations. To test our implementation, we calculate the NAC values for three lanthanide complexes. Finally, the predicted values support the relaxation mechanisms reported in previous studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Performance Comparisons for Artificially Propagated and Wild Pacific Lamprey Juveniles and Larvae

ABSTRACT Artificially propagated Pacific lamprey ( Entosphenus tridentatus ) are produced for restoration and for use in dam passage studies to reduce the demand for wild fish. Such uses require that animals are representative of their wild counterparts. Previous work indicated that this is true for Pacific lamprey larvae and juveniles reared in the hatchery with respect to the length of sustained swimming. However, more subtle differences in behaviour and performance that lamprey need to survive have not been assessed. In this study, artificially propagated and wild fish were compared in laboratory tests under no‐flow conditions to examine light avoidance, burrowing speed, burst swim speed, volitional routine swim speed and time to come to rest. Most larvae burrowed in less than a minute, and we found highly significant differences ( p < 0.001) between artificially propagated and wild larvae burrowing times, a critical escape behaviour. This could have implications for studies of larval entrainment at irrigation diversion canals or in turbine boils at dams. Interestingly, all of the wild juveniles tested came to rest quickly after introduction to the chamber (1.5 min), while artificially propagated lamprey swam robotically near the surface and 48% did not come to rest in the first 10 min (median time to rest = 9.5 min). In contrast, wild juveniles quickly (median = 1.47 min) sought areas near the tank bottom and attached strongly with their oral disc. Such behavioural differences could have important survival consequences for artificially propagated lamprey as they approach turbine intakes, bypass screens and irrigation diversion headgates. This study highlights the need to conduct behavioural assays that examine subtleties of fish behaviour that can be missed with traditional swim tunnel comparisons.

Frick, Kinsey [Fish Ecology Division, Northwest Fi↗

Electric current paths in a Si:P delta-doped device imaged by nitrogen-vacancy diamond magnetic microscopy

The recently-developed ability to control phosphorous-doping of silicon at an atomic level using scanning tunneling microscopy, a technique known as atomic precision advanced manufacturing (APAM), has allowed us to tailor electronic devices with atomic precision, and thus has emerged as a way to explore new possibilities in Si electronics. In these applications, critical questions include where current flow is actually occurring in or near APAM structures as well as whether leakage currents are present. In general, detection and mapping of current flow in APAM structures are valuable diagnostic tools to obtain reliable devices in digital-enhanced applications. In this report, we used nitrogen-vacancy (NV) centers in diamond for wide-field magnetic imaging (with a few-mm field of view and micron-scale resolution) of magnetic fields from surface currents flowing in an APAM test device made of a P delta-doped layer on a Si substrate, a standard APAM witness material. We integrated a diamond having a surface NV ensemble with the device (patterned in two parallel mm-sized ribbons), then mapped the magnetic field from the DC current injected in the APAM device in a home-built NV wide-field microscope. The 2D magnetic field maps were used to reconstruct the surface current densities, allowing us to obtain information on current paths, device failures such as choke points where current flow is impeded, and current leakages outside the APAM-defined P-doped regions. Analysis on the current density reconstructed map showed a projected sensitivity of ~0.03 A m -1 , corresponding to a smallest-detectable current in the 200 μm wide APAM ribbon of ~6 μA. These results demonstrate the failure analysis capability of NV wide-field magnetometry for APAM materials, opening the possibility to investigate other cutting-edge microelectronic devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Understanding the Surface Structure and Catalytic Activity of SnO x /Au(111) Inverse Catalysts for CO 2 and H 2 Activation

Carbon dioxide hydrogenation is a promising approach for the reduction of greenhouse gas pollution via the production of fuels and high-value chemicals utilizing C1 chemistry. In this process, the activation of nonpolar molecules, CO 2 and H 2 , at mild conditions is challenging. In this study, we report a well-defined inverse SnO x /Au(111) catalyst that shows the ability to activate both CO 2 and H 2 at room temperature. Scanning tunneling microscopy (STM) and ambient pressure X-ray photoemission spectroscopy (AP-XPS) are combined to understand the surface structure, growth mode, chemical state, and activity of SnO x /Au(111) surfaces. Nanostructures of SnO x at the sub-monolayer level were prepared by depositing Sn on Au(111) followed by O 2 oxidation. For the as-prepared SnO x /Au(111), two-dimensionally formed SnO x thin films on a Au(111) substrate were observed with STM of two different moieties, discernible based on their height: clusters (~0.4 Å) and nanoparticles (NPs, 1–2.5 Å), which are assigned to Sn–Au alloys and SnO x , respectively, in corroboration with XPS analysis. Furthermore, SnO x /Au(111) was annealed under UHV to test its thermal stability. Upon annealing at 400–600 K, a disappearance of SnO x NPs and reappearance of highly dispersed Sn clusters were clearly noticeable from the STM and XPS results, identifying the thermal decomposition of SnO x and subsequent formation of Sn–Au alloys on the surface due to the recombination of Sn clusters with Au. We investigated the reactivity of the SnO x /Au(111) surfaces toward CH 4 , CO 2 , and H 2 . The SnO x /Au(111) surfaces have excellent CO 2 and H 2 activation abilities even at room temperature with negligible reactivity for methane activation. Our AP-XPS results show that H 2 can be activated on the SnO x NPs by the reduction to Sn. For CO 2 , the activation and further dissociation are identified by a reoxidation of Sn with newly formed Sn–O bonds and the formation of surface carbon. Therefore, we propose that SnO x is a potential catalyst or additive to achieve CO 2 hydrogenation under mild conditions.

36 MATERIALS SCIENCE↗

Strongly-overdoped $\mathrm{La}$ 2- $_{x}$$\mathrm{Sr}$$_{x}$$\mathrm{CuO}$ 4 : Evidence for Josephson-coupled grains of strongly-correlated superconductor

The interpretation of how superconductivity disappears in cuprates at large hole doping has been controversial. To address this issue, here we present an experimental study of single-crystal and thin film samples of La 2-x Sr x CuO 4 (LSCO) with x ≥0.25. In particular, measurements of bulk susceptibility on LSCO crystals with x = 0.25 indicate an onset of diamagnetism at T c1 =38.5 K, with a sharp transition to a phase with full bulk shielding at T c2 =18 K, independent of field direction. Strikingly, the in-plane resistivity only goes to zero at T c2 . Inelastic neutron scattering on 0.25 crystals confirms the presence of low-energy incommensurate magnetic excitations with reduced strength compared to lower doping levels. The ratio of the spin gap to T c2 is anomalously large. Our results are consistent with a theoretical prediction for strongly overdoped cuprates by Spivak, Oreto, and Kivelson, in which superconductivity initially develops within disconnected self-organized grains characterized by a reduced hole concentration, with bulk superconductivity occurring only after superconductivity is induced by proximity effect in the surrounding medium of higher hole concentration. Beyond the superconducting-to-metal transition, local differential conductance measurements on an LSCO thin film suggest that regions with pairing correlations survive, but are too dilute to support superconducting order. Future experiments will be needed to test the degree to which these results apply to overdoped cuprates in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A unified thermionic and thermionic-field emission (TE–TFE) model for ideal Schottky reverse-bias leakage current

We present a unified thermionic emission (TE) and thermionic-field emission (TFE) model for the ideal reverse-bias leakage current in Schottky junctions. The unified TE–TFE analytical model advances upon previous analytical TFE models by Murphy–Good and Padovani–Stratton, which are the two most widely adopted models by the community, in two major aspects: (i) the applicability of the TFE expression therein is extended to near-zero surface electric fields by an error-function correction, allowing for the calculation of the total current by a nontrivial sum of TE (over-the-barrier current) and TFE (below-the-barrier current) contributions; therefore, an accurate description of the TE-to-TFE transition region is captured analytically for the first time; (ii) image-force lowering is considered with much-simpler correction terms. Comparisons with the reference numerical model show that the unified TE–TFE model has excellent accuracy, as well as a 10 000× reduction in computation time. The unified model is further tested against experimental data from Schottky barrier diodes based on Si, 4H-SiC, GaN, and β-Ga 2 O 3 , revealing accurate extractions of barrier heights and correct descriptions of the ideal reverse leakage characteristics. With the extended applicable range, improved accuracy, and computational efficiency, the unified TE–TFE model is highly valuable for the design and analysis of devices based on Schottky junctions, as well as for potential integration in technology computer-aided design (TCAD) tools.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Materials response under hypersonic flow conditions probed by multi-modal diagnostics in a benchtop wind tunnel

In this LDRD “Materials response under hypersonic flow conditions probed by multi-modal diagnostics in a benchtop wind tunnel” a peer-reviewed lab-scale hypersonic wind tunnel (WT) was planned to be assembled in-house, along with time and space resolved and related flow and materials diagnostics to derive high quality measurements and the science that is able to inform the development of our modeling and simulation codes for the hypersonic regime. The baseline capability to achieve Mach 3 flow was scheduled to be completed in Q4 2020, in the first year of the LDRD. However due to the COVID-19 delays and travel restrictions starting in March, 2020 this task could not be completed to carry out the materials response studies in-house, or those planned with external partners that were engaged during this LDRD (NASA Langley, NASA Ames, Texas A&M university, and Virginia Tech). Still, using our hypersonic computational fluid simulation codes, many of the hypersonic components have been designed (for Mach 3, 5, and 7) and procured (or in procurement phase), and a new dedicated hypersonic laboratory was facilitized, 2 WT test sections have been manufactured, and a Mach 3 nozzle was 3D printed for prototyping. In parallel, unique materials diagnostics have been tested under static conditions and fully specified and awaiting validation in our partner facilities. Due to the new hypersonic initiative announced in January, 2020 we proposed a related scope increase in hypersonics research requiring energy interactions and new diagnostics well beyond what could be supported by an LDRD ER. A decision was made to end the current LDRD ER and formally propose a new LDRD SI that leverages the current effort, which will be starting in FY21. As a result of the uniqueness of this LDRD work, our early efforts have been encouraged and well received in the hypersonic research community, and is now part of multiple unsolicited internal and external project proposals expected to be funded in the coming years, which is encouraging to make LLNL a leader in experimentally probing extreme physics of local energetics critical to hypersonics.

36 MATERIALS SCIENCE↗

Validation of Calibrated k–ε Model Parameters for Jet-in-Crossflow

Previous efforts determined a set of calibrated, optimal model parameter values for Reynolds-averaged Navier–Stokes (RANS) simulations of a compressible jet in crossflow (JIC) using a $k–ε$ turbulence model. These parameters were derived by comparing simulation results to particle image velocimetry (PIV) data of a complementary JIC experiment under a limited set of flow conditions. Here, a $k–ε$ model using both nominal and calibrated parameters is validated against PIV data acquired from a much wider variety of JIC cases, including a realistic flight vehicle. The results from the simulations using the calibrated model parameters showed considerable improvements over those using the nominal values, even for cases that were not used in the calibration procedure that defined the optimal parameters. This improvement is demonstrated using a number of quality metrics that test the spatial alignment of the jet core, the magnitudes of multiple flow variables, and the location and strengths of vortices in the counter-rotating vortex cores on the PIV planes. These results suggest that the calibrated parameters have applicability well outside the specific flow case used in defining them and that with the right model parameters, RANS solutions for the JIC can be improved significantly over those obtained from the nominal model.

42 ENGINEERING↗

Contributions to the BeEST sterile neutrino search

Sterile neutrinos are hypothetical right-handed leptons that have been postulated as extensions to the Standard Model to help explain the origin of neutrino masses and the matter-antimatter asymmetry in the Universe. Sterile neutrinos with keV-scale masses are also a promising candidate for dark matter. LLNL and the Colorado School of Mines are developing an experiment to search for sterile neutrinos by measuring the \Beryllium-7 Electron capture in Superconducting Tunnel junctions (STJs)" (the \BeEST" Experiment). When 7 Be implanted into STJ radiation detectors decays by electron capture, the neutrino escapes and the energy of the recoiling 7 Li nucleus can be measured with high accuracy. Using conservation of energy and momentum, the mass of the emitted neutrino can then be found. Sterile keV neutrinos would generate additional peaks in the spectrum, and the experiment uses high-resolution STJs to observe such a potential signal. During my internship, I contributed to different aspects of the BeEST experiment. I designed a UHV setup for X-ray spectroscopy analysis with STJs. It will allow high-resolution X-ray spectroscopy of different materials to examine the broadening in the BeEST spectra, the origins of which are unknown but may be related to material effects inside the STJ detectors. The design accounts for the specific demands of the BeEST project, but provides a degree of flexibility that makes it useful for other potential STJ projects. Additionally, to test the accuracy of the electronics and data analysis routine, an arbitrary waveform generator (AWG) will be used to generate simulated spectra that contain sterile neutrino signals of a specified mass and admixture. For this I programmed an AWG to output random STJ signals whose amplitudes follow the appropriate user-de ned probability distribution. The AWG signals can be input into the 32-channel preamplifier to simulate any spectrum of interest and test the uniformity of the amplifier chain, as well as the sensitivity of the data analysis.

42 ENGINEERING↗

Integrated Circuit Design for Radiation-Hardened Charge-Sensitive Amplifier Survived up to 2 Mrad

According to the continuous development of metal-oxide semiconductor (MOS) fabrication technology, transistors have naturally become more radiation-tolerant through steadily decreasing gate-oxide thickness, increasing the tunneling probability between gate-oxide and channel. Unfortunately, despite this radiation-hardened property of developed transistors, the field of nuclear power plants (NPPs) requires even higher radiation hardness levels. Particularly, total ionizing dose (TID) of approximately 1 Mrad could be required for readout circuitry under severe accident conditions with 100 Mrad around a reactor in-core required. In harsh radiating environments such as NPPs, sensors such as micro-pocket-fission detectors (MPFD) would be a promising technology to be operated for detecting neutrons in reactor cores. For those sensors, readout circuits should be fundamentally placed close to sensing devices for minimizing signal interferences and white noise. Therefore, radiation hardening ability is necessary for the circuits under high radiation environments. This paper presents various integrated circuit designs for a radiation hardened charge-sensitive amplifier (CSA) by using SiGe 130 nm and Si 180 nm fabrication processes with different channel widths and transistor types of complementary metal-oxide-semiconductor (CMOS) and bipolar CMOS (BiCMOS). These circuits were tested under γ–ray environment with Cobalt-60 of high level activity: 490 kCi. The experiment results indicate amplitude degradation of 2.85%–34.3%, fall time increase of 201–1730 ns, as well as a signal-to-noise ratio (SNR) of 0.07–11.6 dB decrease with irradiation dose increase. These results can provide design guidelines for radiation hardening operational amplifiers in terms of transistor sizes and structures.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Performance evaluation of cosmic ray muon trajectory estimation algorithms

Muons, being elementary particles with minimal interaction with nuclear materials and abundant at sea level, have sparked interest in utilizing them for imaging various applications, such as mining [Borselli et al., Sci. Rep. 12, 22329 (2022)], volcano imaging [Nagamine et al., Nucl. Instrum. Meth. A, 356, 585(1995)], and underground tunnel detection [Guardincerri et al., Pure Appl. Geophys. 174, 2133 (2017)]. Recently, their use in nuclear nonproliferation and safeguard verification has gained attention, particularly in cargo screening for nuclear waste smuggling [Baesso et al., J. Instrum. 9, C10041 (2014)], source localization [L. J. Schultz et al., Nucl. Instrum. Meth. A 519, 687 (2004)], and locating nuclear fuel debris in reactors [Borozdin et al., Phys. Rev. Let. 109, 152501 (2012)]. However, the resolution of muon image reconstruction techniques is limited due to multiple Coulomb scattering (MCS) within the target object. To achieve robust muon tomography, it is crucial to develop efficient and flexible physics-based algorithms that can model the MCS process accurately and estimate the most probable trajectory of muons as they pass through the target object. To address this limitation, in this study, a novel algorithmic approach utilizing the Bayesian probability theory and Gaussian approximation of MCS is chosen. Different energy levels, materials, and target sizes were considered in the evaluations. The results demonstrate that the Generalized Muon Trajectory Estimation (GMTE) algorithm offers significant improvements over currently used algorithms. Across all test scenarios, the GMTE algorithm demonstrated ~50% and 38% increase in precision compared to Straight Line Path (SLP) and Point of Closest Approach (PoCA) algorithms, respectively. Furthermore, it exhibited 10%–35% and 10%–15% increases in muon flux utilization for high and medium Z materials, respectively, compared to the PoCA algorithm. In conclusion, the extensive simulations confirm the enhanced performance and efficiency of the GMTE algorithm, offering improved resolution and reduced measurement time for cosmic ray muon imaging compared to the current SLP and PoCA algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗