Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Material Point Method”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Monte Carlo Simulations of Crystal Defects in Open Ensembles

Zero- and two-dimensional crystal defects form in open statistical ensembles, such as the grand canonical, that are usually inaccessible with conventional simulation techniques. This longstanding challenge is overcome with a new Hamiltonian Monte Carlo method that samples energy-biased gradual transformations. In conclusion, the method enables free energy calculations for nonideal point defects and the direct prediction of finite-temperature interface structures.

Grain boundaries↗
DOE OSTI2021

Evaluation of Attila and MCNP computational methods for dose and exposure estimation

Radiation transport calculations are often used to estimate dose or exposure to components and personnel surrounding a radiation source. The sources for these calculations are decaying radionuclides within various nuclear materials. Historically, dose calculations use MCNP (Monte Carlo N-Particle) transport code as the primary particle transport tool without a secondary computational tool to validate the results from the MCNP simulations [1]. The goal of this study is to make an independent check of the Monte Carlo solution from MCNP6 Version 6.2.1 with the discrete ordinates solution from Attila 10.2.0 Beta 3. As an example problem for this study, water-filled, stainless-steel vessels, modeled with an unstructured mesh (UM) with both MCNP and Attila [2], are exposed to 252Cf and 60Co point sources. This report also includes a discussion of the limitations of unstructured mesh in a MCNP calculation.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Using molecular dynamics to predict the solidus and liquidus of mixed oxides ( Th , U ) O 2 , ( Th , Pu ) O 2 and ( Pu , U ) O 2

Molecular dynamics (MD) was used to establish a mechanistic basis for the experimentally observed reduction in liquidus and solidus temperatures below the melting point of the end-members for the mixed oxides (Th,U)O 2 , (Th, Pu) O 2 and (Pu, U) O 2 . This dip is found at additions of the oxide with higher melting point to the oxide with the lower melting point. There are many causes suggested for the dip; here the distribution of the cation Frenkel energy for the mixed oxides caused by the local environment is proposed as a contributor. Furthermore, a variant of the moving interface method which yields information on the position of the solidus and liquidus boundaries, is used to predict the phase diagrams of these systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Evidence for nitrogen binding to surface defects for topological insulator Bi 2 Se 3

Using scanning tunneling spectroscopy and theoretical simulations we have studied the effects of nitrogen gas exposure on the electronic density of states of Bi 2 Se 3 , a well-studied topological insulator. In carefully controlled measurements, Bi 2 Se 3 crystals were initially cleaved in a helium gas environment and then exposed to a 22 SCFH flow of ultra-high purity N 2 gas. We observe a resulting change in the spectral curves, with the exposure effect saturating after approximately 50 min, ultimately bringing the Dirac point about 50 meV closer to the Fermi level. Furthermore, these results are compared to density functional theory calculations, which support a picture of N 2 molecules physisorbing near Se vacancies. Furthermore, ab initio molecular dynamics simulations aided by a Blue Moon ensemble method reveal the dissociative adsorption of N 2 molecules which then bind strongly to Se vacancies at the surface. In this scenario, the binding of the N atom to a Se vacancy site removes the surface defect state created by the vacancy and changes the position of the Fermi energy with respect to the Dirac point.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Cool carriers: triplet diffusion dominates upconversion yield

Perovskites have gained popularity both as the active material in photovoltaics and as bulk triplet sensitizers for solid-state triplet–triplet annihilation upconversion (TTA-UC). Prior to widespread implementation into commercial photovoltaics, an in-depth understanding of the environmental influences on device performance is required. To this point, the temperature-dependent structure–function properties of TTA-UC within methylammonium formamidinium lead triiodide (MAFA)/rubrene UC devices are explored. A strong temperature dependence of the underlying UC dynamics is observed, where the maximum UC efficiency is achieved at 170 K, reflecting the competition between triplet diffusion length, diffusion rate, and triplet–triplet encounter events. Furthermore, a combination of spectroscopic and structural methods and theoretical modelling illustrates that despite the significantly increased carrier lifetime of the perovskite at low temperatures, the TTA-UC dynamics are not governed by the underlying sensitizer properties but rather limited by the underlying triplet diffusion.

Sullivan, Colette M.↗
DOE OSTI2025

Machine-learning-enabled on-the-fly analysis of RHEED patterns during thin film deposition by molecular beam epitaxy

Thin film deposition is a fundamental technology for the discovery, optimization, and manufacturing of functional materials. Deposition by molecular beam epitaxy (MBE) typically employs reflection high-energy electron diffraction (RHEED) as a real-time in situ probe of the growing film. However, the state-of-the-art for RHEED analysis during deposition requires human observation. Here, we present an approach using machine learning (ML) methods to monitor, analyze, and interpret RHEED images on-the-fly during thin film deposition. In the analysis workflow, RHEED pattern images are collected at one frame per second and featurized using a pretrained deep convolutional neural network. The feature vectors are then statistically analyzed to identify changepoints; these changepoints can be related to changes in the deposition mode from initial film nucleation to a transition regime, smooth film deposition, and in some cases, an additional transition to a rough, islanded deposition regime. The feature vectors are additionally analyzed via graph analysis and community classification. The graph is quantified as a stabilization plot, and we show that inflection points in the stabilization plot correspond to changes in the growth regime. The full RHEED analysis workflow is termed RHAAPsody and includes data transfer and output to a visual dashboard. We demonstrate the functionality of RHAAPsody by analyzing the precaptured RHEED images from epitaxial depositions of anatase TiO2 on SrTiO3(001) and show that the analysis workflow can be executed in less than 1 s. Our approach shows promise as one component of ML-enabled real-time feedback control of the MBE deposition process.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Impact of Integration Scheme on Performance of Anisotropic Plasticity Models

Given the prevalent role of metals in a variety of industries, schemes to integrate corresponding constitutive models in finite element applications have long been studied. A number of formulations have been developed to accomplish this task; each with their own advantages and costs. Often the focus has been on ensuring the accuracy and numerical stability of these algorithms to enable robust integration. While important, emphasis on these performance metrics may often come at the cost of computational expense potentially neglecting the needs of individual problems. In the current work, the performance of two of the most common integration methods for anisotropic plasticity -- the convex cutting plane (CCP) and closest point projection (CPP) -- across a variety of metrics is assessed; including accuracy and cost. A variety of problems are considered ranging from single elements to large representative simulations including both implicit quasistatic and explicit transient dynamic type responses. The relative performance of each scheme in the different instances is presented with an eye towards guidance on when the different algorithms may be beneficial.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Palladium at high pressure and high temperature: A combined experimental and theoretical study

Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid–solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Palladium at high pressure and high temperature: A combined experimental and theoretical study

Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid– solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Biobased chemical recycling: aminolysis of PET using renewable reagents and monomers to synthesize new semi-aromatic polyamides

Chemical recycling of PET is a method of depolymerizing polymer chains to monomeric components enabling the synthesis of second-generation materials with virgin-like quality. Commercial chemical recycling techniques rely upon high pressure methanolysis to create precursors capable of synthesizing a second-generation PET resin. However, despite the circular approach of methanolysis, a product with a very short lifespan and similar value is created. The approach of the current study is to utilize aminolysis as an ambient pressure technique to create precursors for higher value materials with longer lifespans to address the current crisis in plastic waste. Semi-aromatic polyamides (SAP) are desired in this circumstance because of their high melting point and heat resistance combined with good melt-processability similar to aliphatic polyamides. In this study SAPs were synthesized using precursors recovered from the aminolysis of PET employing biobased diamines and dicarboxylic acids. While aminolysis has been explored in previous studies, this work investigated the use of biobased components from castor oil: decamethylene diamine during recycling and sebacic acid during polymerization. Polymer synthesis resulted in the formation of SAPs similar to polyphthalamides (PPA) with novel structures given the aromatic portion from terephthalic acid (TPA) and aliphatic portion from the diamines and diacids. The synthesized materials exhibited excellent thermal stability with high glass transition temperatures. Novel polymers were created with varying aliphatic chain length to understand fundamental parameters needed to produce a valuable polymer from post-consumer waste.

High performance Applications↗
DOE OSTI2021

(U) Using the MCNP6 Perturbation Capability for Source Nuclide Density Sensitivities

The Taylor series (differential operator) perturbation method, as implemented in the PERT capability in MCNP6 (Ref. 1), can be used for first- and second-order sensitivity analyses in fixed-source problems. Unfortunately, the PERT capability does not yet work with point-detector (F5) or pulse-height (F8) tallies. Also, the second-order Taylor series term lacks the cross terms that would allow the calculation of mixed sensitivities. Finally, the PERT capability does not compute any sensitivities involving the source in a fixed-source problem. When neutrons or gamma rays are emitted due to radioactive decay, the source rate density depends intrinsically on the nuclide densities in the source material. When the PERT capability in MCNP6 is used to compute sensitivities of a response to nuclide and material densities, those sensitivities do not include the effect of the source.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Efficient and universal characterization of atomic structures through a topological graph order parameter

Abstract A graph-based order parameter, based on the topology of the graph itself, is introduced for the characterization of atomistic structures. The order parameter is universal to any material/chemical system and is transferable to all structural geometries. Four sets of data are used to validate both the generalizability and accuracy of the algorithm: (1) liquid lithium configurations spanning up to 300 GPa, (2) condensed phases of carbon along with nanotubes and buckyballs at ambient and high temperature, (3) a diverse set of aluminum configurations including surfaces, compressed and expanded lattices, point defects, grain boundaries, liquids, nanoparticles, all at nonzero temperatures, and (4) eleven niobium oxide crystal phases generated with ab initio molecular dynamics. We compare our proposed method to existing, state-of-the-art methods for the cases of aluminum and niobium oxide. Our order parameter uniquely classifies every configuration and outperforms all studied existing methods, opening the door for its use in a multitude of complex application spaces that can require fine structure-level characterization of atomistic graphs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Final Technical Report for "Cloud-based Low-Scaling Quantum Chemistry Simulations for Materials"

The objective of phase I was to develop the basic infrastructure for performing efficient hybrid DFT calculations on solid materials with a turnaround time amenable to the use in large-scale data-driven approaches. To fulfill the goal we have developed a novel algorithm that significantly reduces the cost of exchange matrix formation. The algorithm does so by making use of three ingredients (a) interpolative decomposition of the electron integrals (b) robust pseudospectral method and (c) occ-RI approach. Our published work demonstrates that the algorithm is orders of magnitude faster than any other hybrid-DFT method and is on the order of only 3-4 times slower than pure DFT. All steps of the algorithm developed during phase I can be accelerated to the point that, for large enough systems, the diagonalization of the Fock matrix will become the most expensive step. To treat such large systems we have created an interface to the ASCR-funded PEXSI library. The PEXSI method is used to compute the density matrix directly from the Fock matrix in a manner that preserves the sparsity of the local representation.

Shiozaki, Toru↗
DOE OSTI2021

In-situ synchrotron X-ray micro-diffraction investigation of ultra-low-strain deformation microstructure in laminated Ti-Al composites

An ultra-low-strain deformation microstructure was revealed for the first time non-destructively in the bulk interior of an annealed laminated Ti-Al composite-in the "fully recrystallized" Al layer-by a synchrotron-based micro-diffraction technique, namely differential aperture X-ray microscopy (DAXM), through real space mapping with a very high angular resolution (0.01 degrees). This ultra-low-strain deformation microstructure was found to result from the thermal stress, induced during cooling after annealing, due to the different coefficients of thermal expansion for the Ti and Al layers. The annealed sample was further tensile deformed to a strain of 1.66% and followed by in situ DAXM and analyzed by various misorientation mapping methods. The results pointed to the important effects of the initial microstructure and the interface constraint, as well as the grain size and crystal orientation, on the plastic deformation. A gradient in dislocation density from the layer interface to the center was found in the Al layer of the annealed sample, and this gradient increased slightly during tensile deformation. The variation of the dislocation density was further discussed based on the activation and interaction of dislocations in grains of different sizes and orientations during plastic deformation. The findings of this study provided valuable insights in understanding the constraint effect of the laminated metal composite and the design of novel composite materials. (c) 2020 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Crystallographic orientation↗
DOE OSTI2023

Spherical powders: Control over the size and morphology of powders for additive manufacturing and enriched stable isotope nuclear targets

Metal powders are a fundamental starting point for fabricating many types of nuclear targets. Elemental powder properties can differ drastically between batches, even when using the same method. Therefore, the variation in morphology and the size of metal powders can cause variable quality and produce inconsistent results with what are otherwise proven target manufacturing techniques. Additive manufacturing has additional requirements for higher quality and more uniform feedstock. The production of spheroidized powders with uniform, reproducible properties and a narrow size distribution represents unexplored opportunities for experiments. These opportunities include experimenting with solid metals that can now flow like liquids, new options for powder handling and dispensing, and new target fabrication methods using additive manufacturing. The Stable Isotope Materials and Chemistry Group at Oak Ridge National Laboratory obtained an AMAZEMET rePowder ultrasonic metal atomization tool for creating limited batches of fully dense, free flowing, spherical powders with a narrow size distribution of extremely rare materials. Early results are presented with materials that were produced. The team explores the anticipated limits of this instrument with extremely rare materials (e.g., enriched stable isotopes) and highlights research into new fabrication techniques that provide additional options benefitting the international nuclear target community.

Zach, Mike↗
DOE OSTI2020

Opportunities to Implement Solutions to Achieve Remedial Action Objectives in Consideration of Stakeholder Interests - 20365

A case study in the implementation of a soils excavation and removal project conducted at a Formerly Utilized Sites Remedial Action Program (FUSRAP) Maywood Superfund Site (FMSS) vicinity property; an active commercial business, in a densely populated area, with significant operational, technical and, logistical constraints that were expected to limit areas that could be remediated. The United States Army Corps of Engineers-led (USACE) team successfully navigated complex overlapping stakeholder interests to achieve a more efficient and complete removal of contaminated soils and debris while minimizing unnecessary excess costs to the Government. This paper centers on FUSRAP activities at the FMSS vicinity property located at 149-151 Maywood Avenue, Maywood, Bergen County, New Jersey. For most of its FUSRAP history, this ∼109,000 square meter (27-acre) property housed a now-demolished ∼26,000 square meter (6.5-acre) warehouse operated by Sears Logistics Services, a unit of the Sears Roebuck Company (Figure 2). Sears ended its property lease and vacated the warehouse in December 2016. Given that long history and for the purposes of this paper, the property will be referred to as the 'Sears property' or simply 'the property.' The Sears property is approximately 109,000 square meters (27-acres) and is currently zoned for commercial use. The property is bound to the north and northwest by 100 West Hunter Avenue (the Stepan Company), to the northeast by 205 Maywood Avenue, to the east by Maywood Avenue, to the south by 23 West Howcroft Road, and to the west by businesses on NJ Route 17 and the NJ Route 17 roadway. Until December 2016, the property housed a warehouse and distribution center operated by Sears Logistical Services. On-site structures were demolished by the property owner in 2017. The warehouse covered the north section of the property along the Stepan Company property line. A railroad spur from the adjacent property now known as the Maywood Interim Storage Site (MISS) ended at the northeast corner of the warehouse. Wetlands were located east of the warehouse; the rest of the property was covered with paved parking lots and grassed areas. During remediation at the property, opportunities frequently arose where the project team was able to coordinate effectively with stakeholders to sequence remediation activities to facilitate removal of otherwise inaccessible soils and identify opportunities for materials reuse when supported by residual radiological and chemical levels. Some specific opportunities include: - remediating individual truck loading dock bays to maintain tenant operations; - working along a busy highway and in a utility corridor containing a 30-inch high pressure gas main; - remediating Lodi Brook and associated wetlands requiring bypass pumping and other diversionary structures to maintain local community stormwater drainage; - protecting during remediation and supporting access to facilitate tenant/owner maintenance of critical fire protection and water supply system service lines buried in contaminated soils and necessary for safe warehouse operations; and - coordinating with stakeholders to ensure non-radiological contaminants of concern for the site (unaffiliated with FUSRAP) were addressed by the responsible party in a manner that maximized benefit to all parties while supporting an efficient overall site remediation program. Once the property was vacant and the warehouse structure and radiologically non-impacted above-grade structures were demolished, the project identified an opportunity to significantly reduce the volume of waste associated with the foundation of the former warehouse, a foundation suspected of being partially constructed in radiologically contaminated soils. The project developed and implemented supplemental radiological verification survey and sampling strategies based on radiological cross-contamination risk potential with process flow-charts and screening-level based decision points; a program that classified saw-cut sections of concrete based on observed residual surface radioactivity conditions using a combination of gamma sensitive sodium-iodide scans and beta sensitive Geiger-Muller direct measurements. The proposed approach was reviewed with stakeholders with feedback, including consideration of applicable State of New Jersey Site Remediation Program criteria and guidance before implementation. Once classified, additional sampling was performed at frequencies driven by screening. The sampling methods, stockpile sampling frequencies, criteria and, approach to data evaluation were developed in consideration of existing site cleanup standards, regional background ranges, State of New Jersey radiological remediation program guidance, and specific stakeholder input With the property remediation and survey efforts nearly completed, it is relevant to examine retrospectively the objectives, assumptions and limitations in the remedial design (i.e., what was planned) versus what was able to be accomplished through effective teaming. Benefits to other environmental remediation site programs include better understanding of how to work effectively with stakeholders to achieve win-win outcomes, and approaches to site remediation, waste minimization and materials reuse that were successful in their overall outcomes. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

In-situ Detection of Glovebox Glove Degradation Prior to Glove Failure; Effects on Safety, Waste, and Operational Costs - 20402

Glovebox glove pressure decay leak testing and analysis of the data can greatly improve glovebox safety, minimize Transuranic (TRU) glove waste, and minimize operational cost. It is well known, the weakest point of containment on a glovebox is the glove. Mitigating unplanned openings in gloves is critical to minimizing operational and safety costs in glovebox operations. Mitigating unplanned glove openings due to glove failures will be discussed in this paper. The lack of an engineered solution to determine a gloves operational life in nuclear operations, has lead to use of theoretical and manual subjective methods to determine the life and safety of a glovebox glove. Whether the glove is inspected visually prior to use, or gloves are changed regularly to theoretically avoid failure, both the operators PPE and the room are being exposed to contamination, or un-necessary amounts of TRU glove waste are being generated. Millions of dollars yearly in operational costs are spent on glove failures due to loss of production, incident analysis, regulatory audits, clean up, and paperwork. In a nuclear application, a glove failure could contaminate an operator and laboratory, shut down the operation for weeks, potentially create a regulatory audit, and potentially create a media frenzy; the damages to the organization can be astronomical. Elimination of glove failures is possible by tracking glove material degradation and setting limits to allow glove change when it becomes necessary, prior to glove failure. Analytical pressure decay leak testing allows the detection of material degradation, and in turn the ability to limit operational glove failures. A German based company called MK Versuchsanlagen, e.K., has developed an advanced Glove Integrity Testing System. The system, with its use of proprietary RFID and special software technology, is unique in its capabilities. The system is capable to perform highly accurate regular glove testing on any size glovebox line in minimal time. The system records and can analyze a tremendous amount of system data that in turn can determine, for example, a safe glove change time prior to a glove material failure. This paper will show how an advanced Glove Integrity Testing System can detect glove material degradation and help determine a safe change point prior to glove failure. Pressure decay curves of new gloves and recorded effects of accelerated aging on the glove over time will be shown. The use of this technology in nuclear facilities can greatly improve glovebox operational safety, minimize glovebox glove TRU waste, and save organizations tremendous costs in downtime, clean up, documentation, and potential regulatory review. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

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↗