Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “object reconstruction”

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

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

At least 271 records · Page 15

A comparative analysis of YOLOv8 and U-Net image segmentation approaches for transmission electron micrographs of polycrystalline thin films

Metallic thin films offer a platform to experimentally study the dynamics of microstructural evolution, but the required transmission electron microscopy (TEM)-based imaging generates complex images that are challenging to segment and quantify. This work provides a comparative analysis of a new YOLOv8 model and an established U-Net model for bright-field TEM images of polycrystals, employing a framework leveraging physical observables to evaluate performance against two hand-traced benchmark datasets. This methodology obviates the comparison of large, diversely structured, and manually labeled datasets that are required to assess performance on a per-image/per-pixel basis. It is found that the YOLOv8 model, adapted for real-time instance segmentation, has up to 43× faster inferencing (NVIDIA GeForce RTX 4090) compared to U-Net and reconstructs hand-traced grain size distributions (GSDs) with excellent fidelity, finding mean diameter within 3% for grains near an optimal magnification; for grains that deviate from the optimal pixel-diameter, the size of small- (large)-diameter grains is systematically over- (under)-estimated. This is partially mitigated by including scale-aware augmentations during training. Moreover, when the bias is corrected post-inference by a rigid shift in distribution, the YOLOv8 model reproduces ground truth GSDs with exceptional fidelity, with statistical tests indicating <5% probability that the distributions are distinct. Based on ground truth data, calibration curves pertaining to this shift can be constructed for a given model. This issue is not present in the U-Net model’s results, indicating that for quantitative measurements where the true size of objects is of interest, special procedures must be implemented for YOLO-based models.

36 MATERIALS SCIENCE↗

Reconstruction and identification of boosted di-τ systems in a search for Higgs boson pairs using 13 TeV proton-proton collision data in ATLAS

In this paper, a new technique for reconstructing and identifying hadronically decaying τ+τ- pairs with a large Lorentz boost, referred to as the di-τ tagger, is developed and used for the first time in the ATLAS experiment at the Large Hadron Collider. A benchmark di-τ tagging selection is employed in the search for resonant Higgs boson pair production, where one Higgs boson decays into a boosted $b\overline b$ pair and the other into a boosted τ+τ- pair, with two hadronically decaying τ-leptons in the final state. Using 139 fb -1 of proton-proton collision data recorded at a centre-of-mass energy of 13 TeV, the efficiency of the di-τ tagger is determined and the background with quark- or gluon-initiated jets misidentified as di-τ objects is estimated. The search for a heavy, narrow, scalar resonance produced via gluon-gluon fusion and decaying into two Higgs bosons is carried out in the mass range 1–3 TeV using the same dataset. No deviations from the Standard Model predictions are observed, and 95% confidence-level exclusion limits are set on this model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Feasibility study of a compact neutron resonance transmission analysis instrument

Neutron Resonance Transmission Analysis (NRTA) uses resonant absorption of neutrons to infer the absolute isotopic composition of a target object, enabling applications in a broad range of fields such as archaeology, enrichment analysis of nuclear fuel, and arms control treaty verification. In the past, NRTA involved large user facilities and complex detector systems. However, recent advances in the intensity of compact neutron sources have made compact neutron imaging designs increasingly feasible. This work describes the Monte Carlo (MC) based design of a compact epithermal NRTA radiographic instrument which uses a moderated, compact deuterium-tritium (DT) neutron source and an epithermal neutron detector. Such an instrument would have a wide range of applications, and would be especially impactful for such scenarios as nuclear inspection and arms control verification exercises, where system complexity and mobility may be of critical importance. The MC simulations presented in this work demonstrate accurate time-of-flight (TOF) reconstructions for transmitted neutron energies, capable of differentiating isotopic compositions of nuclear material with high levels of accuracy. A new generation of miniaturized and increasingly more intense neutron sources will allow this technique to achieve measurements with greater precision and speed, with significant impact on a variety of engineering and societal problems.

47 OTHER INSTRUMENTATION↗

Neural network-based single material beam-hardening correction for X-ray CT in Additive Manufacturing

Beam-hardening (BH) artifacts are ubiquitous in X-ray CT scans of dense metal additively manufactured (AM) parts. While linearization approaches are useful for correcting beam-hardened data from single material objects, they either require a calibration scan or detailed system and material composition information. In this paper, we introduce a neural network-based, material-agnostic method to correct beam-hardening artifacts. We train a neural network to map the acquired beam-hardened projection values and the corresponding estimated thickness of the part based on an initial segmentation to beam-hardening related parameters, which can be used to compute the coefficients of a linearizing correction polynomial. A key strength of our approach is that, once the network is trained, it can be used for correcting beam hardening from a variety of materials without any calibration scans or detailed system and material composition information. Furthermore, our method is robust to errors in the estimated thickness due to the typical challenge of obtaining an accurate initial segmentation from reconstructions impacted by BH artifacts. We demonstrate the utility of our method to obtain high-quality CT reconstructions from a collection of AM parts -- suppressing cupping and streaking artifacts

Rahman, Obaid↗

New techniques for jet calibration with the ATLAS detector

A determination of the jet energy scale is presented using proton–proton collision data with a centre-of-mass energy of $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 140 fb -1 collected using the ATLAS detector at the LHC. Jets are reconstructed using the ATLAS particle-flow method that combines charged-particle tracks and topo-clusters formed from energy deposits in the calorimeter cells. The anti-k t jet algorithm with radius parameter R = 0.4 is used to define the jet. Novel jet energy scale calibration strategies developed for the LHC Run 2 are reported that lay the foundation for the jet calibration in Run 3. Jets are calibrated with a series of simulation-based corrections, including state-of-the-art techniques in jet calibration such as machine learning methods and novel in situ calibrations to achieve better performance than the baseline calibration derived using up to 81 fb -1 of Run 2 data. The performance of these new techniques is then examined in the in situ measurements by exploiting the transverse momentum balance between a jet and a reference object. The b -quark jet energy scale using particle flow jets is measured for the first time with around 1% precision using $\gamma$+jet events.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Demonstration of object location, classification, and characterization by developed deep learning dust ablation trail analysis code package using plasma jets

Based on deep learning, a Dust Ablation Trail Analysis (DATA) code package is developed to detect dust ablation trails in tokamaks, which is intended to analyze a large amount data of tokamak dusts. To validate and benchmark the DATA code package, 2440 plasma jet images are exploited for the training and test of the deep learning DATA code package, since plasma jets resemble the shape and size of dust ablation clouds in tokamaks. After being trained by 1920 plasma jet images, the DATA code package is able to locate 100% plasma jets, classify plasma jets with the accuracy of >99.9%, and output image skeleton information for classified plasma jets. The DATA code package trained by the plasma jet images is also used to analyze the dust ablation trails captured in the Experimental Advanced Superconducting (EAST) tokamak with the satisfactory performance, further verifying its applicability in the fusion dust ablation investigation. Based on its excellent performance presented here, it is demonstrated that our DATA code package is able to automatically identify and analyze dust ablation trails in tokamaks, which can be used for further detailed investigations, such as the three-dimensional reconstruction of dusts and their ablation trails.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Background-free dual-mode optical and 13 C magnetic resonance imaging in diamond particles

Multimodal imaging—the ability to acquire images of an object through more than one imaging mode simultaneously—has opened additional perspectives in areas ranging from astronomy to medicine. In this paper, we report progress toward combining optical and magnetic resonance (MR) imaging in such a “dual” imaging mode. They are attractive in combination because they offer complementary advantages of resolution and speed, especially in the context of imaging in scattering environments. Our approach relies on a specific material platform, microdiamond particles hosting nitrogen vacancy (NV) defect centers that fluoresce brightly under optical excitation and simultaneously “hyperpolarize” lattice C 13 nuclei, making them bright under MR imaging. We highlight advantages of dual-mode optical and MR imaging in allowing background-free particle imaging and describe regimes in which either mode can enhance the other. Leveraging the fact that the two imaging modes proceed in Fourier-reciprocal domains (real and k-space), we propose a sampling protocol that accelerates image reconstruction in sparse-imaging scenarios. Our work suggests interesting possibilities for the simultaneous optical and low-field MR imaging of targeted diamond nanoparticles.

36 MATERIALS SCIENCE↗

Improvements to the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. This document summarizes the transport solver development efforts pursued during Fiscal Year 2022. We added the multiphysics transient capability for the coarse-mesh finite difference accelerated Richardson iteration for discontinuous finite element method discrete ordinates (DFEM-SN) scheme to support high-order heterogeneous transport simulations. HFEM (hybrid finite element method) - PN (spherical harmonics expansion) was completed and red-black iteration was added for solving the HFEM-PN system with both preconditioned Jacobian-free Newton Krylov and Richardson iteration solvers. The HFEM-PN scheme, as one of the low-order transport schemes, is expected for supporting routine design simulations. Pin power reconstruction capability was also designed and implemented with the Griffin ISOXML module to enhance all the low-order transport solvers for more accurate multiphysics simulations. Numerical results are presented for demonstrating the capabilities and verifying their performance, and future works are discussed.

97 MATHEMATICS AND COMPUTING↗

MatPhase: Material phase prediction for Li-ion Battery Reconstruction using Hierarchical Curriculum Learning

Li-ion Batteries (LIB), one of the most efficient energy storage devices, are used extensively in many industrial applications. These batteries consist of electrodes that are put together with heterogeneous material compositions. Imaging data of these battery electrodes obtained from X-ray tomography can explain the distribution of material constituents and allow reconstructions to study electron transport pathways. Such reconstructions of material constituents help quantify various associated properties of electrodes (e.g., volume-specific surface area, porosity) which determine the performance of batteries. These images often suffer from low image contrast between multiple material constituents, hence making it difficult for humans to distinguish and characterize these constituents through visual inspection. A minor error in detecting distributions of the material constituents can lead to magnified errors in the calculated parameters of material properties (e.g., porosity). We present MatPhase, a novel hierarchical curriculum learning technique to address the complex task of estimating material constituent distribution in battery electrodes. MatPhase comprises three modules: (i) an uncertainty-aware global model trained to yield inferences conditioned upon global knowledge of material distribution, (ii) a local model to capture relatively more fine-grained (local) distributional signals, (iii) an aggregator model to appropriately fuse the local and global effects towards obtaining the final distribution. On average, MatPhase improves prediction up to 8.5% relative to other sophisticated modeling pipelines and state-of-the-arts (SOTA) object detection models employed in the performance comparison.

Tabassum, Anika↗

Reconsidering tympanal-acoustic interactions leads to an improved model of auditory acuity in a parasitoid fly

Although most binaural organisms locate sound sources using neurological structures to amplify the sounds they hear, some animals use mechanically coupled hearing organs instead. One of these animals, the parasitoid fly Ormia ochracea (O. ochracea), has astoundingly accurate sound localization abilities. It can locate objects in the azimuthal plane with a precision of 2°, equal to that of humans, despite an intertympanal distance of only 0.5 mm, which is less than 1/100th of the wavelength of the sound emitted by the crickets that it parasitizes. O. ochracea accomplishes this feat via mechanically coupled tympana that interact with incoming acoustic pressure waves to amplify differences in the signals received at the two ears. In 1995, Miles et al developed a model of hearing mechanics in O. ochracea that represents the tympana as flat, front-facing prosternal membranes, though they lie on a convex surface at an angle from the flies' frontal and transverse planes. The model works well for incoming sound angles less than ±30° but suffers from reduced accuracy (up to 60% error) at higher angles compared to response data acquired from O. ochracea specimens. Despite this limitation, it has been the basis for bio-inspired microphone designs for decades. Here, we present critical improvements to this classic hearing model based on information from three-dimensional reconstructions of O. ochracea's tympanal organ. We identified the orientation of the tympana with respect to a frontal plane and the azimuthal angle segment between the tympana as morphological features essential to the flies' auditory acuity, and hypothesized a differentiated mechanical response to incoming sound on the ipsi- and contralateral sides that depend on these features. We incorporated spatially-varying model coefficients representing this asymmetric response, making a new quasi-two-dimensional (q2D) model. The q2D model has high accuracy (average errors of under 10%) for all incoming sound angles. This improved biomechanical model may inform the design of new microscale directional microphones and other small-scale acoustic sensor systems.

36 MATERIALS SCIENCE↗

Cannibals in the thick disk II. Radial-velocity monitoring of the young α -rich stars

Determining ages of stars for reconstructing the history of the Milky Way remains one of the most difficult tasks in astrophysics. This involves knowing when it is possible to relate the stellar mass with its age and when it is not. The young α-rich (YAR) stars present such a case in which we are still not sure about their ages because they are relatively massive, implying young ages, but their abundances are α-enhanced, which implies old ages. We report the results from new observations from a long-term radial-velocity-monitoring campaign complemented with high-resolution spectroscopy, as well as new astrometry and seismology of a sample of 41 red giants from the third version of APOKASC, which includes YAR stars. The aim is to better characterize the YAR stars in terms of binarity, mass, abundance trends, and kinematic properties. The radial velocities of HERMES, APOGEE, and Gaia were combined to determine the binary fraction among YAR stars. In combination with their mass estimate, evolutionary status, chemical composition, and kinematic properties, it allowed us to better constrain the nature of these objects. We found that stars with M < 1 M ⊙ were all single, whereas stars with M > 1 M ⊙ could be either single or binary. This is in agreement with theoretical predictions of population synthesis models. Studying their [C/N], [C/Fe], and [N/Fe], trends with mass, it became clear that many YAR stars do not follow the APOKASC stars, favoring the scenario that most of them are the product of mass transfer. Our sample further includes two likely undermassive stars, that is to say of such as low mass that they cannot have reached the red clump within the age of the Universe, unless their low mass is the signature of mass loss in previous evolutionary phases. These stars do not show signatures of currently being binaries. Both YAR and undermassive stars might show some anomalous APOGEE abundances for the elements N, Na, P, K, and Cr; although, higher-resolution optical spectroscopy might be needed to confirm these findings. Considering the significant fraction of stars that are formed in pairs and the variety of ways that makes mass transfer possible, the diversity in properties in terms of binarity, and chemistry of the YAR and undermassive stars studied here implies that most of these objects are likely not young.

79 ASTRONOMY AND ASTROPHYSICS↗

The incorporation of environmentally derived 87 Sr / 86 Sr and Sr/Ca in early otolith formation of Chinook salmon

Abstract Effective species management often requires understanding patterns of movement and habitat use. A common approach in identifying where individuals reside relies upon chemical tracers from the environment that are incorporated into an individual's tissues. For fish, isotopes in their otoliths, specifically the portion of their otolith formed during their larval stage, have been used to identify the natal origin. Complicating this work, however, is the fact that during this life stage, there is a shift in the source of isotopes deposited onto the growing otolith from maternally to environmentally derived. The objective of this study was to identify the portion of the otolith representing this transition to environmentally derived isotopes so as to accurately investigate questions of natal origin for a threatened population of fall Chinook salmon ( Oncorhynchus tshawytscha ). We exposed developing larvae to four treatments that differed in terms of their water strontium isotope ratio ( 87 Sr/ 86 Sr) and used change‐point analysis of otolith 87 Sr/ 86 Sr and strontium to calcium ratio (Sr/Ca) to identify the otolith radius corresponding to the transition to environmentally derived isotopes. Our results indicated this transition occurred, on average, at 132 μm ( 87 Sr/ 86 Sr; ±50 μm standard deviation) and 127 μm (Sr/Ca; ±29 μm) from the otolith core, which corresponded to the developmental time between hatching and exogenous feeding. A substantial proportion of our otoliths (i.e., 61%) did not show convergence between otolith and water 87 Sr/ 86 Sr by the end of the 113‐day experiment, which was likely due to the dietary contribution of marine‐based feed. Therefore, we were unable to recommend an otolith radius to target for the purposes of reconstructing natal origin apart from being beyond approximately 130 μm.

Chittaro, Paul↗

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

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

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Improved Cross Section Generation Capability of Griffin in FY22

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of deterministic steady-state transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication as well as transient solvers for spatial dynamics with the improved quasi-static method. Griffin uses cross section data in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. In recent years, the MC 2 -3 modules have been added to Griffin for fast reactor cross section generation, and the self-shielding application programming interface (SSAPI) was implemented in the ISOXML module for thermal reactor cross section generation. The on-the-fly slowing down method and double-heterogeneity treatment have been implemented to SSAPI and verified against particulate fuel-bearing graphite-moderated thermal reactor problems with high accuracy. This year, work has been focused on improving the cross section generation capability of ISOXML and streamlining the cross section generation procedures. In addition, the form function data were added to ISOXML in order to support the pin power reconstruction capability that was newly implemented in Griffin in this fiscal year. To facilitate the cross section generation using MC 2 -3 and SSAPI in Griffin, the cross section generation workflows have been set up for both fast and thermal spectrum reactors. The MOOSE action system tool was devised for fast spectrum problems, and the MOOSE stochastic tool was adopted to the branch calculation procedure for thermal spectrum problems. Meanwhile, to ensure the accuracy of group-constants, the thermal up-scattering kernel calculator accounting for resonance scattering was implemented in ISOXML, demonstrating the accurate computation of a Doppler-broadened scattering kernel of any Legendre order within a reasonable timescale. Other aspects of ISOXML, such as deletion solver and data, documentation, ISOXML file management, and the interface for Mixture, were improved as well.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Status of ICARUS NUMI Interaction Cross-Section Analysis

The ICARUS experiment, utilizing Liquid Argon Time Projection Chamber (LAr TPC) technology, has been installed at Fermilab in Chicago, Illinois, following its initial operation in Italy and subsequent refurbishment at CERN. ICARUS completed commissioning in June 2022. Currently, the experiment is in the phase of analyzing data from its two runs of physics data acquisition and gearing up for the third run. While its primary objective is to function as the far detector of the Short Baseline Neutrino program (SBN), seeking sterile neutrino signatures, ICARUS also offers diverse physics capabilities, including searches beyond the standard model and measurements of cross-sections. In addition to being exposed to the common Booster Neutrino (BNB) beamline, ICARUS also receives off-axis neutrinos from the Main Injector (NuMI) beam. Due to the off-axis angle between NuMI and ICARUS, coupled with contributions from both pion and kaon decays to neutrino fluxes, interactions of NuMI neutrinos within ICARUS can be detected over a range of several GeV in energy. These interactions present opportunities for crucial cross-section measurements and model tests within an energy range that overlaps both the SBN oscillation search and a portion of the DUNE spectrum. This poster presentation will delve into our efforts to conduct a muon-neutrino cross-section measurement, where the signal is defined by events with no pions produced in the final state of the interaction, along with some preliminary muon-neutrino inclusive measurements. Additionally, it will provide updates on the current status and future plans, including reconstruction, selection, and analysis procedures.

43 PARTICLE ACCELERATORS↗

Modernization efforts for the R -Matrix code SAMMY [Abstract]

The R-Matrix code SAMMY is a widely used nuclear data evaluation code focused on the resolved range, which includes corrections for experimental effects. The code is still mostly written in Fortran 77, and uses a memory management system suitable for the time of its initial writing (1984). A modernization effort is under way to bring the code in-line with modern software development practices. A continuous-integration testing framework was added, automating the large existing set of test cases. It is run on every commit. The memory management was updated to current standard practices suitable for modern software analysis tools. The code can be obtained from https://code.ornl.gov/RNSD/SAMMY. The resonance parameters and covariance information are now stored in C++ objects shared by SAMMY and AMPX, the processing code that generates nuclear data libraries for SCALE. This allows for easier maintenance and access to the resonance parameters inside and outside of SAMMY. This feature is already used by accessing and changing parameters in memory in the Bayesian Monte Carlo Evaluation Framework for Cross Sections Nuclear Data and Integral Benchmark Experiments project, Further plans include the switch to the ENDF reading and writing routines in AMPX, as these routines are more robust, easier to maintain, and support more features. Of note here is support for the new GNDS format. Previously it wasn’t easy to share the full covariance matrix for evaluations containing more than one isotope due to limitations on the ENDF format; this is now supported in GNDS. The data are currently available in a binary SAMMY format and can be exported to GNDS to make them more widely available and sharable. The next step will be to use the same resonance processing code at 0K in AMPX and SAMMY as one of the available Reich-Moore R-Matrix formalism. The first step toward this goal is to isolate the reconstruction into a module that takes resonance parameters as its input and does not depend on SAMMY global parameters. This goal has been achieved and it should now be possible to more easily change the resonance formalism and add enhancements as the Phenomenological R-Matrix parameterization of direct, doorway, and compound nuclear reactions discussed elsewhere on this conference. This concerted modernization and enhancement effort provides multiple advantages to the nuclear data community. It will allow parameter optimization using enhanced formalisms, including experimental effects, that better match complex experimental data. Then those evaluated parameters can immediately be passed off to AMPX to be reconstructed with the exact same cross section model and be put into a data library for subsequent testing using SCALE and the Valid Benchmark suite or other suitable benchmark suites.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Accuracy Enhancement of Nuclear Power Plant Simulators Utilizing High Accuracy Simulation Predictions

More recently, reactor core simulators for core designs associated with commercial nuclear power plants that utilize what is believed to be higher fidelity models have been developed. Features such as neutronics models that utilize transport equation solvers with fine spatial meshes and many energy-groups, thermal-hydraulic models that utilize sub-channel solvers with fine spatial mesh and capable of treating a wide range of fluid conditions, and fuel-coolant chemistry interaction models capable of treating CRUD deposition are to be found in these higher fidelity core simulators. These reactor core simulators require access to higher performance computers, characterized by many processors, cores and large memory. So associated with utilization of these simulators is access to high performance computers and ability to accommodate in one’s workflow longer execution times. By contrast, currently used core simulators by the nuclear industry can execute on engineering workstations and have execution times of seconds to minutes. The desirability for having short execution times is not only desired for support of time critical tasks but supports the mental process of decision making by engineers. The goal of the work reported upon here has the objective of retaining the fidelity of higher fidelity models while retaining the ability to utilize engineering workstations. Beyond the core simulator goal, additional goals of this work include incorporating the just described core simulator capability into a Nuclear Steam Supply System (NSSS) simulator, and to incorporate the resulting capability into an environment supportive of design and operational decision making associated with nuclear power stations. The model selected for the core neutronics model is the NESTLE code, for the core thermal-hydraulic model is the CTF code utilizing coarse mesh, and for the NSSS model is the RELAP5-3D code. WSC’s proprietary 3KEYMASTERTM platform is being used to provide software coupling, user interface, visualization, and reporting. The NESTLE core neutronics simulator was first integrated with the CTF core thermal-hydraulic simulator using CTF developed communication commands which are also used for CTF to communicate with RELAP5-3D under WSC’s proprietary 3KEYMASTERTM platform. To assure NESTLE prediction consistency with higher fidelity core neutronic simulators, buffer codes have been created to automatically generate from output files written by the VERA core simulator the NESTLE nodal neutronic parameter’ library, geometry, and pin-power reconstruction input files, thereby avoiding a number of challenges associated with utilizing lattice physics codes and providing consistency with VERA predictions. To treat absorber rod effects a multi-set library is utilized, where a set refers to a specific absorber rod fully inserted pattern. A coarse spatial mesh CTF model was developed with features added that support using CTF as envisioned in the engineering quality simulator. A hybrid meshing approach was implemented to allow for automated construction of models with mixed levels of refinement. Specifically, a core model could resolve some assemblies at a nodal level (4 subchannels per assembly) and others at a pin-resolution (one subchannel per coolant subchannel in the assembly). The intention is that this will allow for better resolution of limiting conditions such as DNBR and PCT, which are based on local rod and subchannel conditions. Further development was done of features that enhance the capabilities for the envisioned engineering quality simulator that has been developed, but now for RELAP-3D. The RELAP5-3D code development includes ability to model more than 999 components and the addition of the cross-channels turbulence mixing model and the void drift model that are implemented in CTF, aiming to achieve closer prediction agreement of the two codes for transient simulations, specifically, more accurate matches of the overall mass, momentum, and energy exchanges of both the liquid and gas phases between the neighboring core assemblies. Graphics were also developed for the Instructor Station for this project under WSC’s proprietary 3KEYMASTERTM platform to facilitate design and operational decision making.

42 ENGINEERING↗

Toward Drilling the Perfect Geothermal Well: An International Research Coordination Network for Geothermal Drilling Optimization Supported by Deep Machine Learning and Cloud Based Data Aggregation

The EDGE project, supported by the U.S. Department of Energy Geothermal Technologies Office under award DE-EE0008793, established a data-driven framework for improving the efficiency, cost-effectiveness, and reliability of geothermal well drilling. The project focused on developing scalable data infrastructure, advanced machine learning and probabilistic models, and integrated analytics tools to support continuous drilling optimization. A central objective was to reduce geothermal drilling costs by up to seventy percent while minimizing the risk of well failure through predictive diagnostics and adaptive planning. Over the project period, a comprehensive data repository was designed and deployed, incorporating records from over one hundred geothermal wells across varied geological settings. This repository supported both structured and unstructured data and adhered to FAIR data principles, enabling provenance tracking, quality control, and standardized metadata. The project introduced automated ingestion pipelines and a cloud-hosted platform that facilitated access to raw, processed, and derived datasets. This infrastructure served as the foundation for model development and analysis. Machine learning workflows were developed to predict key drilling metrics including rate of penetration, non-productive time, and total drilling costs. Self-organizing maps and dimensionality reduction methods were used to uncover operational patterns and outliers, while supervised learning algorithms such as random forests and deep neural networks were applied to forecast performance outcomes. The models were validated on heterogeneous datasets from both U.S. and Icelandic fields, demonstrating variable but significant predictive accuracy. The results indicated that finer temporal resolution, inclusion of lithological data, and consistency in operational annotations could substantially improve model performance. The project also implemented process mining techniques to reconstruct state-transition models from drilling event logs. These models enabled the identification of deviations from optimal workflows and provided insights into recurring failure modes. Analysis of non-productive time highlighted the impact of equipment failures, geological challenges, and human factors, offering opportunities for targeted mitigation strategies. The EDGE Dashboard was developed as a web-based expert system integrating data visualization, model outputs, and user-driven queries. It provided an accessible interface for operators to explore historical data, evaluate predicted outcomes, and compare drilling scenarios. Initial feedback from project partners suggested that the dashboard could serve as a foundation for more advanced advisory and optimization tools. Overall, the EDGE project demonstrated the feasibility and value of applying modern data science techniques to geothermal drilling. It delivered a set of interoperable tools and models that can support more efficient, lower-risk well development. The findings point toward a viable path for transitioning from advisory analytics to semi-autonomous drilling systems, contingent on continued collaboration, expanded datasets, and field validation. The project results have immediate relevance for drilling operations, data management practices, and future geothermal R&D efforts aimed at achieving reliable, cost-competitive geothermal energy at scale.

15 GEOTHERMAL ENERGY↗