Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model checking”

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 289 records · Page 16

What does fault tolerant Deep Learning need from MPI?

Deep Learning (DL) algorithms have become the {\em de facto} Machine Learning (ML) algorithm for large scale data analysis. DL algorithms are computationally expensive -- even distributed DL implementations which use MPI require days of training (model learning) time on commonly studied datasets. Long running DL applications become susceptible to faults -- requiring development of a fault tolerant system infrastructure, in addition to fault tolerant DL algorithms. This raises an important question: {\em What is needed from MPI for designing fault tolerant DL implementations?} In this paper, we address this problem for permanent faults. We motivate the need for a fault tolerant MPI specification by an in-depth consideration of recent innovations in DL algorithms and their properties, which drive the need for specific fault tolerance features. We present an in-depth discussion on the suitability of different parallelism types (model, data and hybrid); a need (or lack thereof) for check-pointing of any critical data structures; and most importantly, consideration for several fault tolerance proposals (user-level fault mitigation (ULFM), Reinit) in MPI and their applicability to fault tolerant DL implementations. We leverage a distributed memory implementation of Caffe, currently available under the Machine Learning Toolkit for Extreme Scale (MaTEx). We implement our approaches by extending MaTEx-Caffe for using ULFM-based implementation. Our evaluation using the ImageNet dataset and AlexNet neural network topology demonstrates the effectiveness of the proposed fault tolerant DL implementation using OpenMPI based ULFM.

Amatya, Vinay C.↗

Collaborative Research: Explosive reconnection in relativistic magnetically-dominated plasmas (Final Report)

Magnetic fields can be the dominant component of astrophysical plasmas, so that the magnetic energy density might exceed even the rest-mass energy density of matter. In this extreme (and largely unexplored) regime the magnetic field controls the overall plasma evolution, dissipation, and acceleration of non-thermal particles. This plasma regime, applicable to a variety of astrophysical sources - magnetars, pulsars and pulsar wind nebulae (PWNe), jets of Active Galactic Nuclei (AGNs) and Gamma-Ray Bursters (GRBs) - is dramatically different from laboratory plasmas, the magnetospheres of planets, and the interplanetary plasma. Relativistic astrophysical sources then provide an unique opportunity to study the fundamental plasma physics of magnetically-dominated plasmas; a novel and fast-evolving field of theoretical research which, by investigating energy conversion and particle energization processes in plasmas, is of vital importance to the Fusion Energy Sciences DoE program. Data coming from astrophysical high-energy missions, especially the Crab Nebula flares recently observed by the Fermi and AGILE satellites, suggest that the acceleration of non-thermal particles to the highest energies occurs in magnetic reconnection events - a major change of paradigm in high-energy astrophysics. Most importantly, observations demand that particle acceleration should proceed extremely fast (with accelerating electric field of the order of the magnetic field) and on macroscopic scales (much larger, e.g., than the microscopic plasma skin depth). We are conducting studies of the microphysics of magnetically-dominated plasmas focussing in particular on the highly dynamic regime of explosive reconnection and associated particle acceleration in relativistic plasmas. We are studying the stability and explosive plasma dynamics, particle acceleration and radiation production in a number of idealized plasma configurations that approximate relevant astrophysical sources (like the magnetic ABC structures and interacting flux tubes, as well as generalizations of analytical models of X-point collapse to relativistic plasmas). The well-studied case of the Crab Nebula is taken as a prototypical example for the application of the model. We are combining analytical studies of explosive magnetic dynamics and dissipation in relativistic plasmas with particle-in-cell (PIC) simulations and fluid simulations. The theoretical model, fluid and particle-in-cell simulations are cross-checked for agreement and convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-E-T-51: Bayesian Network Models for Radiotherapy Error Detection

Purpose: To develop a probabilistic model of radiotherapy plans using Bayesian networks that will detect potential errors in radiation delivery. Methods: Semi-structured interviews with medical physicists and other domain experts were employed to generate a set of layered nodes and arcs forming a Bayesian Network (BN) which encapsulates relevant radiotherapy concepts and their associated interdependencies. Concepts in the final network were limited to those whose parameters are represented in the institutional database at a level significant enough to develop mathematical distributions. The concept-relation knowledge base was constructed using the Web Ontology Language (OWL) and translated into Hugin Expert Bayes Network files via the the RHugin package in the R statistical programming language. A subset of de-identified data derived from a Mosaiq relational database representing 1937 unique prescription cases was processed and pre-screened for errors and then used by the Hugin implementation of the Estimation-Maximization (EM) algorithm for machine learning all parameter distributions. Individual networks were generated for each of several commonly treated anatomic regions identified by ICD-9 neoplasm categories including lung, brain, lymphoma, and female breast. Results: The resulting Bayesian networks represent a large part of the probabilistic knowledge inherent in treatment planning. By populating the networks entirely with data captured from a clinical oncology information management system over the course of several years of normal practice, we were able to create accurate probability tables with no additional time spent by experts or clinicians. These probabilistic descriptions of the treatment planning allow one to check if a treatment plan is within the normal scope of practice, given some initial set of clinical evidence and thereby detect for potential outliers to be flagged for further investigation. Conclusion: The networks developed here support the use of probabilistic models into clinical chart checking for improved detection of potential errors in RT plans.

60 APPLIED LIFE SCIENCES↗

LeWRON: Agentic Analysis of Electroweak Phase Transitions

The electroweak phase transition (EWPT) is a central topic in particle physics and cosmology, connecting collider phenomenology, baryogenesis, and gravitational-wave observatories. Its analysis requires a technically demanding, convention-sensitive, and model-dependent pipeline, from constructing the finite-temperature effective potential to tracking thermal histories, computing bubble nucleation rates, and predicting gravitational-wave spectra. We present LeWRON (Learning ElectroWeak phase tRansitiON), an agentic framework that orchestrates this pipeline starting from an input Lagrangian. LeWRON combines audited toolbox construction with an Explorer module that uses the generated model-specific code for further analysis, including scans and plots. Intermediate analytic outputs are checked by auditor agents and stored as structured artifacts, enabling reproducible human inspection and downstream use through both a command-line interface and a public Python API. The framework supports a reproduction mode, which infers conventions from the literature and reproduces published results, and a discovery mode, which guides users through structured checkpoints for new models. We demonstrate LeWRON across representative beyond-the-Standard-Model scenarios and release the code on GitHub.

Wang, Isaac R. [Fermilab] (ORCID:000000030789218X)↗

Numerical Modelling of Radionuclide Migration in the Context of a Near-Surface LILW Disposal Facility - 18618

The reference solution for the management of low and intermediate level short-lived waste (LILW) in France consists in their emplacement in concrete vaults built at the near-surface. The corresponding disposal facility under operation by the French waste management organization (Andra), is located in North-eastern France and planned to be closed in 2060. The safety case developed by Andra has to show the evidence that, amongst other aspects, the radionuclide transfer through the confinement barriers will occur in such a way that the human and environmental exposure will be as low as reasonably achievable (ALARA) and in accordance with the dose constraints. For this purpose, modelling is used to integrate uncertainties about the characteristics of concrete components, waste and water table and their evolution in time and space. At last, modelling is finally used to investigate the radiological impact of specific scenarios in which the failure of safety barriers occurs. In the frame of its technical review of the safety assessment, the French Institute for Radiological Protection and Nuclear Safety (IRSN) also performs such kind of numerical modelling of radionuclide fate from the waste to the hydrogeological outlets, in order to check whether the assessment performed by Andra takes into account sufficiently conservative assumptions. The present work aims to evaluate the role of components of the vaults on radionuclide transport during the post-closure period over a few thousand years. For this purpose, IRSN carries out numerical calculations with MELODIE software, which simulates water flow and transport of species, in order to simulate radionuclide migration from waste packages emplaced in a concrete vault to the surrounding groundwater. The developed transport model focuses on 3-H radionuclide and integrates Andra's main assumptions (release model from the waste packages, chemical and hydraulic properties of the vault components, shrinkage and cracks in concrete components). An adaptive mesh method, based on an a posteriori error estimator, was used to model the cracks and the consideration of a strong heterogeneity in the hydraulic and transport parameters values in the different components of the model. This method shows an improvement of results of about 40 times, whereas the mesh size was increased by only 2%. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Development of a point-kinetic verification scheme for nuclear reactor applications

In this paper, a new method that can be used for checking the proper implementation of time- or frequency-dependent neutron transport models and for verifying their ability to recover some basic reactor physics properties is proposed. This method makes use of the application of a stationary perturbation to the system at a given frequency and extraction of the point-kinetic component of the system response. Even for strongly heterogeneous systems for which an analytical solution does not exist, the point-kinetic component follows, as a function of frequency, a simple analytical form. The comparison between the extracted point-kinetic component and its expected analytical form provides an opportunity to verify and validate neutron transport solvers. The proposed method is tested on two diffusion-based codes, one working in the time domain and the other working in the frequency domain. As long as the applied perturbation has a non-zero reactivity effect, it is demonstrated that the method can be successfully applied to verify and validate time- or frequency-dependent neutron transport solvers. Although the method is demonstrated in the present paper in a diffusion theory framework, higher order neutron transport methods could be verified based on the same principles.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An extensive analysis of the triple W UMa type binary FI BOO

We present a detailed analysis of the interesting W UMa binary FI Boo in view of the spectroscopic signature of a third body through photometry, period variation, and a thorough investigation of solution uniqueness. We obtained new BVR{sub c}I{sub c} photometric data that, when combined with spectroscopic data, enable us to analyze the system FI Boo and determine its basic orbital and physical properties through PHOEBE, as well as the period variation by studying the times of the minima. This combined approach allows us to study the long-term period changes in the system for the first time in order to investigate the presence of a third body and to check extensively the solution uniqueness and the uncertainties of derived parameters. Our modeling indicates that FI Boo is a W-type moderate (f = 50.15% ± 8.10%) overcontact binary with component masses of M {sub h} = 0.40 ± 0.05 M {sub ☉} and M {sub c} = 1.07 ± 0.05 M {sub ☉}, temperatures of T {sub h} = 5746 ± 33 K and T {sub c} = 5420 ± 56 K, and a third body, which may play an important role in the formation and evolution. The results were tested by heuristic scanning and parameter kicking to provide the consistent and reliable set of parameters that was used to obtain the initial masses of the progenitors (1.71 ± 0.10 M {sub ☉} and 0.63 ± 0.01 M {sub ☉}, respectively). We also investigated the evolutionary status of massive components with several sets of widely used isochrones.

79 ASTRONOMY AND ASTROPHYSICS↗

The Newton constant and gravitational waves in some vector field adjusting mechanisms

At the present, there exist some Lorentz breaking scenarios which explain the smallness of the cosmological constant at the present era [1]–[2]. An important aspect to analyze is the propagation of gravitational waves and the screening or enhancement of the Newton constant G {sub N} in these models. The problem is that the Lorentz symmetry breaking terms may induce an unacceptable value of the Newton constant G {sub N} or introduce longitudinal modes in the gravitational wave propagation. Furthermore this breaking may spoil the standard dispersion relation ω= ck . In [3] the authors have presented a model suggesting that the behavior of the gravitational constant is correct for asymptotic times. In the present work, an explicit checking is made and we finally agree with these claims. Furthermore, it is suggested that the gravitational waves are also well behaved for large times. In the process, some new models with the same behavior are obtained, thus enlarging the list of possible adjustment mechanisms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SU-E-T-64: CG-Based Radiation Therapy Simulator with Physical Modeling for Avoidance of Collisions Between Gantry and Couch Or Patient

Purpose: It is time-consuming and might cause re-planning to check couch-gantry and patient-gantry collisions on a radiotherapy machine when using couch rotations for non-coplanar beam angles. The aim of this study was to develop a computer-graphics (CG)-based radiation therapy simulator with physical modeling for avoidance of collisions between gantry and couch or patient on a radiotherapy machine. Methods: The radiation therapy simulator was three-dimensionally constructed including a radiotherapy machine (Clinac iX, Varian Medical Systems), couch, and radiation treatment room according to their designs by using a physical-modeling-based computer graphics software (Blender, free and open-source). Each patient was modeled by applying a surface rendering technique to their planning computed tomography (CT) images acquired from 16-slice CT scanner (BrightSpeed, GE Healthcare). Immobilization devices for patients were scanned by the CT equipment, and were rendered as the patient planning CT images. The errors in the collision angle of the gantry with the couch or patient between gold standards and the estimated values were obtained by fixing the gantry angle for the evaluation of the proposed simulator. Results: The average error of estimated collision angles to the couch head side was -8.5% for gantry angles of 60 to 135 degree, and -5.5% for gantry angles of 225 to 300 degree. Moreover, the average error of estimated collision angles to the couch foot side was -1.1% for gantry angles of 60 to 135 degree, and 1.4% for gantry angles of 225 to 300 degree. Conclusion: The CG-based radiation therapy simulator could make it possible to estimate the collision angle between gantry and couch or patient on the radiotherapy machine without verifying the collision angles in the radiation treatment room.

60 APPLIED LIFE SCIENCES↗

TU-H-CAMPUS-TeP2-05: Selective Protection of Normal Tissue by Cerium Oxide Nanoparticles During Radiation Therapy

Purpose: Cerium oxide nanoparticles (CONPs) have unique pH dependent properties such that they act as a radical modulator. These properties may be used in radiation therapy (RT) to protect normal tissue. This work investigates the selective radioprotection of CONPs in-vitro and potential for in-situ delivery of CONPs in prostate cancer RT. Methods: i) Normal human umbilical vein endothelial cells (HUVEC) and human prostate cancer cells (PC-3) were treated with 0 or 2 ng/mL CONPs (NP size: 5 nm). 2 Gy of 100 kVp radiation was delivered to the cells 4 hours after the CONP treatment. Cell viability was checked 48 hours later using MTS assays. ii) A prostate tumor was modeled as a 2-cm diameter sphere. CONPs were proposed to be loaded in a hollow radiotherapy fiducial marker. The concentration profile for the CONPs within the tumor was modeled with a previously validated diffusion equation employed in other studies for nanoparticles 10 nm or less. Results: i) Without radiation, cell viability was above 90% when treated with 2 ng/mL CONPs for both HUVEC and PC-3. After irradiation, a slightly higher viability was observed in HUVEC with CONPs than the ones without CONPs, and this effect was not observed in PC-3. ii) Based on the calculations, 2 ng/mL of CONPs could be delivered to normal cells by diffusion with a 1 µg/mL initial concentration within two weeks. Conclusion: We conclude that CONPs can provide selective radioprotection. The delivery of needed concentrations of CONPs is feasible via in-situ release from radiotherapy biomaterials (e.g. fiducials) loaded with the CONPs.

60 APPLIED LIFE SCIENCES↗

Kip, Version 1.0

This high-performance ray tracing library provides very fast rendering, compact code, type flexibility via C++ 'generic programming' techniques, and ease of use through an application programming interface (API) that functions independently of any GUI, on-screen display, or other enclosing application. Kip supports constructive solid geometry models based on any of several built-in shapes and operators and allows for new, user-defined shapes and operators. Additional features include basic texturing, input/output of models using a simple human-readable file format with full error checking and detailed diagnostics, and support for shared data parallelism. Kip is written in pure, ANSI Standard C++, is entirely platform independent, and is very easy to use, as it requires no build system, configuration or installation scripts, wizards, special preprocessing, makefiles, or external libraries.

Staley, Martin↗

Utah FORGE: GES Well 16A(78)-32 and Well 16B(78)-32 Stimulation Seismic Event Catalogs

This dataset contains seismic event catalogs from the hydraulic stimulation of wells 16A(78)-32 and 16B(78)-32 at the Utah FORGE site in April 2024. The data was collected by Geo Energy Suisse (GES) using a variety of seismic monitoring technologies, including 3-component (3C) geophones and distributed acoustic sensing (DAS) systems. These technologies were deployed across several locations, including wells 16A, 16B, and Delano-1, with sensor arrays at multiple depths to capture microseismic activity during the stimulations. The catalogs provide both real-time and manually checked seismic event locations, with detailed parameters such as trigger conditions, velocity models, and data acquisition settings. The dataset includes information on the stimulation stages, event rates, and hydraulic injection conditions for each well, with a report detailing the data acquisition configuration and seismic event location methodologies. Users will need to reference the included report for a complete understanding of the sensor network, data processing techniques, and accuracy considerations.

15 GEOTHERMAL ENERGY↗

SolarAPP+ Pilot Analysis: Performance and Impact of Instant, Online Solar Permitting

The National Renewable Energy Laboratory (NREL) led a collaborative effort to develop the Solar Automated Permit Processing Plus (SolarAPP+), a no-cost solar permitting software solution to address residential solar photovoltaic (PV) permitting resource constraints and streamline solar adoption processes among authorities having jurisdiction (AHJs). The SolarAPP+ is an online portal that automates permit plan review, enabling an instant permit approval process for code-compliant residential PV systems. Based on national model building, electrical, and fire codes, the SolarAPP+ automatically performs a compliance check of permit inputs against code requirements and produces an inspection checklist that can be used to verify installation practices, workmanship, and adherence to the approved design. NREL conducted a two-phase pilot with five participating AHJs and 16 solar contractors, spanning from November 2020 through December 2021, evaluating the ability of SolarAPP+ to deliver instant permits and its impact on four critical areas within each AHJ: permit review timelines; solar adoption timelines; AHJ time saved; and inspection results. The pilot confirmed that SolarAPP+ reduced permit review times to less than one day, reduced solar adoption timelines by an average of 12 days, saved an estimated 2,067 staff hours, and had comparable inspection results to those found in traditional permitting. The pilot further identified key lessons for improving SolarAPP+, including ways to streamline SolarAPP+ adoption, development of training resources, and expansion of supported products.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

POSTER: Automatic Differentiation of Parallel Loops with Formal Methods

The accompanying poster to this short paper presents a combination of reverse mode AD and formal methods to enable efficient differentiation of (or backpropagation through) shared-memory parallel code. Compared to the state of the art, our approach can more often avoid the need for atomic updates or private data copies during the parallel derivative computation, even in the presence of unstructured or data-dependent data access patterns. This is achieved by gathering information about the memory access patterns from the input program, which is assumed to be correctly parallelized. This information is then used to build a model of assertions in a theorem prover, which can be used to check the safety of shared memory accesses during the parallel derivative computation

Automatic Differentiation↗

Conversational Grid Storage: Bridging Rucio and LLMs with Model Context Protocol

Experiments at Fermilab use Rucio to handle datasets that can be up to exabyte scale. However, navigating through Rucio’s syntax-heavy Command Line Interface (CLI) is a major workflow obstruction for researchers who just want to check quotas, track data identifiers (DIDs), or locate data sets. This project introduces a natural language interface. By building a containerized Model Context Protocol (MCP) server, an AI agent is created that translates plain English queries into data operations.

Akella, Kashyap [Fermilab; Illinois U., Urbana (ma↗

Planck 2015 results: XI. CMB power spectra, likelihoods, and robustness of parameters

This study presents the Planck 2015 likelihoods, statistical descriptions of the 2-point correlationfunctions of the cosmic microwave background (CMB) temperature and polarization fluctuations that account for relevant uncertainties, both instrumental and astrophysical in nature. They are based on the same hybrid approach used for the previous release, i.e., a pixel-based likelihood at low multipoles (ℓ< 30) and a Gaussian approximation to the distribution of cross-power spectra at higher multipoles. The main improvements are the use of more and better processed data and of Planck polarization information, along with more detailed models of foregrounds and instrumental uncertainties. The increased redundancy brought by more than doubling the amount of data analysed enables further consistency checks and enhanced immunity to systematic effects. It also improves the constraining power of Planck, in particular with regard to small-scale foreground properties. Progress in the modelling of foreground emission enables the retention of a larger fraction of the sky to determine the properties of the CMB, which also contributes to the enhanced precision of the spectra. Improvements in data processing and instrumental modelling further reduce uncertainties. Extensive tests establish the robustness and accuracy of the likelihood results, from temperature alone, from polarization alone, and from their combination. For temperature, we also perform a full likelihood analysis of realistic end-to-end simulations of the instrumental response to the sky, which were fed into the actual data processing pipeline; this does not reveal biases from residual low-level instrumental systematics. Even with the increase in precision and robustness, the ΛCDM cosmological model continues to offer a very good fit to the Planck data. The slope of the primordial scalar fluctuations, n s , is confirmed smaller than unity at more than 5σ from Planck alone. We further validate the robustness of the likelihood results against specific extensions to the baseline cosmology, which are particularly sensitive to data at high multipoles. For instance, the effective number of neutrino species remains compatible with the canonical value of 3.046. For this first detailed analysis of Planck polarization spectra, we concentrate at high multipoles on the E modes, leaving the analysis of the weaker B modes to future work. At low multipoles we use temperature maps at all Planck frequencies along with a subset of polarization data. These data take advantage of Planck’s wide frequency coverage to improve the separation of CMB and foreground emission. Within the baseline ΛCDM cosmology this requires τ = 0.078 ± 0.019 for the reionization optical depth, which is significantly lower than estimates without the use of high-frequency data for explicit monitoring of dust emission. At high multipoles we detect residual systematic errors in E polarization, typically at the μK 2 level; we therefore choose to retain temperature information alone for high multipoles as the recommended baseline, in particular for testing non-minimal models. Finally and nevertheless, the high-multipole polarization spectra from Planck are already good enough to enable a separate high-precision determination of the parameters of the ΛCDM model, showing consistency with those established independently from temperature information alone.

79 ASTRONOMY AND ASTROPHYSICS↗

Towards a model of pion generalized parton distributions from Dyson-Schwinger equations

We compute the pion quark Generalized Parton Distribution H{sup q} and Double Distributions F{sup q} and G{sup q} in a coupled Bethe-Salpeter and Dyson-Schwinger approach. We use simple algebraic expressions inspired by the numerical resolution of Dyson-Schwinger and Bethe-Salpeter equations. We explicitly check the support and polynomiality properties, and the behavior under charge conjugation or time invariance of our model. We derive analytic expressions for the pion Double Distributions and Generalized Parton Distribution at vanishing pion momentum transfer at a low scale. Our model compares very well to experimental pion form factor or parton distribution function data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fracture Network and Fluid Flow Imaging for Enhanced Geothermal Systems Applications from Multi-Dimensional Electrical Resistivity Structure

We have developed an algorithm for the inversion of magnetotelluric (MT) data to a 3D earth resistivity model based upon the finite element method. Hexahedral edge finite elements are implemented to accommodate discontinuities in the electric field across resistivity boundaries, and to accurately simulate topographic variations. All matrices are reduced and solved using direct solution modules which avoids ill-conditioning endemic to iterative solvers such as conjugate gradients, principally PARDISO for the finite element system and PLASMA for the parameter step estimate. Large model parameterizations can be handled by transforming the Gauss-Newton estimator to data-space form. Accuracy of the forward problem and jacobians has been checked by comparison to integral equations results and by limiting asymptotes. Inverse accuracy and performance has been verified against the public Dublin Secret Test Model 2 and the well-known Mount St Helens 3D MT data set. This algorithm we believe is the most capable yet for forming 3D images of earth resistivity structure and their implications for geothermal fluids and pathways.

15 GEOTHERMAL ENERGY↗