SEARCH · Engineering Papers

Results for “Formal Methods”

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

SourceResearch connectionFrom yearThrough year

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

At least 811 records · Page 45

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Single-cone finite difference scheme for the (2+1)D Dirac von Neumann equation

An explicit finite difference scheme is presented for the von Neumann equation for (2+1)D Dirac fermions. It is founded upon a staggered space–time grid which ensures a single-cone energy dispersion and performs the time-derivative in one sweep using a three-step leap-frog procedure. It enables a space–time-resolved numerical treatment of the mixed-state dynamics of Dirac fermions within the effective single-particle density matrix formalism. Energy-momentum dispersion, stability and convergence properties are derived. Elementary numerical tests to demonstrate stability properties use parameters which pertain to topological insulator surface states. A method for the simulation of charge injection from an electric contact is presented and tested numerically. Potential extensions of the scheme to a Dirac–Lindblad equation, real-space–time Green's function formulations, and higher-order finite-difference schemes are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2013

Oscillatory features in the curvature power spectrum after a sudden turn of the inflationary trajectory

In the context of two-field inflation characterized by a light direction and a heavy direction, we revisit the question of the impact of the massive modes on the power spectrum produced after a turn in the inflationary trajectory. We consider in particular the resonant effect due to the background oscillations following a sharp turn. Working in the mass basis, i.e. in the basis spanned by the eigenvectors of the effective mass matrix for the perturbations, we provide an analytical estimate of the resonant effect, using the in-in formalism. In comparison with earlier estimates, we find the same the spectral dependence but a smaller amplitude. We also compute, again via the in-in formalism, the effect of the direct coupling between the light and heavy modes at the instant of the turn and confirm our previous results obtained via a different method.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

SU-E-T-169: Evaluation of Oncentra TPS for Nasopharynx Brachy Using Patient Specific Voxel Phantom and EGSnrc

Purpose: Nasopharnx carcinoma (NPC) treatment is being carried out using Ir-192 HDR seeds in Mehdieh Hospital in Hamadan, Iran. The Oncentra™ TPS is based on optimized TG-43 formalism which disregards heterogeneity in the treatment area. Due to abundant heterogeneity in head and neck, comparison of the Oncentra™ TPS dose evaluation and an accurate dose calculation method in NPC brachytherapy is the objective of this study. Methods: CT DICOMs of a patient with NPC obtained from Mehdieh Hospital used to create 3D voxel phantom with CTCREATE utility of EGSnrc code package. The voxel phantom together with Ir-192 HDR brachytherapy source were the input to DOSXYZnrc to calculate the 3D dose distribution. The sources were incorporate with type 6 source in DOSXYZnrc and their dwell times were taken into account in final dose calculations. Results: The direct comparison between isodoses as well as DVHs for the GTV, PTV and CTV obtained by Oncentra™ and EGSnrc Monte Carlo code are made. EGSnrc results are obtained using 5×10{sup 9} histories to reduce the statistical error below 1% in GTV and 5% in 5% dose areas. The standard ICRP700 cross section library is employed in DOSXYZnrc dose calculation. Conclusion: A direct relationship between increased dose differences and increased material density (hence heterogeneity) is observed when isodoses contours of the TPS and DOSXYZnrc are compared. Regarding the point dose calculations, the differences range from 1.2% in PTV to 5.6% for cavity region and 7.8% for bone regions. While Oncentra™ TPS overestimates the dose in cavities, it tends to underestimate dose depositions within bones.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-G-201-10: Experimental Determination of Modified TG-43 Dosimetry Parameters for the Xoft Axxent® Electronic Brachytherapy Source

Purpose: The establishment of an air kerma rate standard at NIST for the Xoft Axxent{sup ®} electronic brachytherapy source (Axxent{sup ®} source) motivated the establishment of a modified TG-43 dosimetry formalism. This work measures the modified dosimetry parameters for the Axxent{sup ®} source in the absence of a treatment applicator for implementation in Xoft’s treatment planning system. Methods: The dose-rate conversion coefficient (DRCC), radial dose function (RDF) values, and polar anisotropy (PA) were measured using TLD-100 microcubes with NIST-calibrated sources. The DRCC and RDF measurements were performed in liquid water using an annulus of Virtual Water™ designed to align the TLDs at the height of the anode at fixed radii from the source. The PA was measured at several distances from the source in a PMMA phantom. MCNP-determined absorbed dose energy dependence correction factors were used to convert from dose to TLD to dose to liquid water for the DRCC, RDF, and PA measurements. The intrinsic energy dependence correction factor from the work of Pike was used. The AKR was determined using a NIST-calibrated HDR1000 Plus well-type ionization chamber. Results: The DRCC was determined to be 8.6 (cGy/hr)/(µGy/min). The radial dose values were determined to be 1.00 (1cm), 0.60 (2cm), 0.42 (3cm), and 0.32 (4cm), with agreement ranging from (5.7% to 10.9%) from the work of Hiatt et al. 2015 and agreement from (2.8% to 6.8%) with internal MCNP simulations. Conclusion: This work presents a complete dataset of modified TG-43 dosimetry parameters for the Axxent{sup ®} source in the absence of an applicator. Prior to this study a DRCC had not been measured for the Axxent{sup ®} source. This data will be used for calculating dose distributions for patients receiving treatment with the Axxent{sup ®} source in Xoft’s breast balloon and vaginal applicators, and for intraoperative radiotherapy. Sources and partial funding for this work were provided by Xoft Inc. (a subsidiary of iCAD). This work was also supported by the Radiological Sciences T32 Training Grant through the University of Wisconsin-Madison Medical Physics department (5T32CA009206-37).

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

New charged-particle transport computational capability: the SIT code An L-4 milestone

We have developed a new high-fidelity code for direct transport of charged particles using the simple integral transport method. The code can be coupled to the outputs of any hydrodynamical simulation code in 1-D, 2-D or 3-D. In this report we summarize the formalism involved in treating complex transport problems. We present physical examples wherein we have used the code to calculate the transport of alpha particles. Future work is planned to study the sensitivity of hydrodynamical mix to charged-particle radiochemistry and reaction-in-flight neutrons for the complex inertial confinement fusion problems encountered at NIF and at the Z-machine.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2016

Numerical methods for the weakly compressible Generalized Langevin Model in Eulerian reference frame

A well established approach for the computation of turbulent flow without resolving all turbulent flow scales is to solve a filtered or averaged set of equations, and to model non-resolved scales by closures derived from transported probability density functions (PDF) for velocity fluctuations. Effective numerical methods for PDF transport employ the equivalence between the Fokker–Planck equation for the PDF and a Generalized Langevin Model (GLM), and compute the PDF by transporting a set of sampling particles by GLM (Pope (1985) [1]). The natural representation of GLM is a system of stochastic differential equations in a Lagrangian reference frame, typically solved by particle methods. A representation in a Eulerian reference frame, however, has the potential to significantly reduce computational effort and to allow for the seamless integration into a Eulerian-frame numerical flow solver. GLM in a Eulerian frame (GLMEF) formally corresponds to the nonlinear fluctuating hydrodynamic equations derived by Nakamura and Yoshimori (2009) [12]. Unlike the more common Landau–Lifshitz Navier–Stokes (LLNS) equations these equations are derived from the underdamped Langevin equation and are not based on a local equilibrium assumption. Similarly to LLNS equations the numerical solution of GLMEF requires special considerations. In this paper we investigate different numerical approaches to solving GLMEF with respect to the correct representation of stochastic properties of the solution. We find that a discretely conservative staggered finite-difference scheme, adapted from a scheme originally proposed for turbulent incompressible flow, in conjunction with a strongly stable (for non-stochastic PDE) Runge–Kutta method performs better for GLMEF than schemes adopted from those proposed previously for the LLNS. We show that equilibrium stochastic fluctuations are correctly reproduced.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

CAN-D: A Modular Four-Step Pipeline for Comprehensively Decoding Controller Area Network Data

Controller area networks (CANs) are a broadcast protocol for real-time communication of critical vehicle subsystems. Original equipment manufacturers of passenger vehicles hold secret their mappings of CAN data to vehicle signals, and these definitions vary according to make, model, and year. Without these mappings, the wealth of real-time vehicle information hidden in the CAN packets is uninterpretable, severely impeding vehicle-related research, including CAN cybersecurity and privacy studies, aftermarket tuning, efficiency and performance monitoring, and fault diagnosis to name a few. Guided by the four-part CAN signal definition, we present CAN-D (CAN-Decoder), a modular, four-step pipeline for identifying each signal's boundaries (start bit and length), endianness (byte ordering), signedness (bit-to-integer encoding), and by leveraging diagnostic standards, augmenting a subset of the extracted signals with meaningful, physical interpretation. En route to CAN-D, we provide a comprehensive review of the CAN signal reverse engineering research. All previous methods ignore endianness and signedness, rendering them incapable of decoding many standard CAN signal definitions. Incorporating endianness grows the search space from 128 to 4.72E21 signal tokenizations and introduces a web of changing dependencies. In response, we formulate, formally analyze, and provide an efficient solution to an optimization problem, allowing identification of the optimal set of signal boundaries and byte orderings. In addition, we provide two novel, state-of-the-art signal boundary classifiers—both of which are superior to previous approaches in precision and recall in three different test scenarios—and the first signedness classification algorithm, which exhibits a $>$ 97% F-score. Altogether, CAN-D is the only solution with the potential to extract any CAN signal that is also the state of the art. In evaluation on 10 vehicles of different makes, CAN-D's average $\ell ^1$ error is five times better (81% less) than all previous methods and exhibits lower average error, even when considering only signals that meet prior methods’ assumptions. Finally, CAN-D is implemented in lightweight hardware, allowing for an on-board diagnostic (OBD-II) plugin for real-time in-vehicle CAN decoding.

42 ENGINEERING↗

DOE OSTI2018

Light and unbound nuclei: An introduction to ab initio methods in nuclear structure and reaction theory

Light nuclei display a wide variety of phenomena: from stable bound states, to rare unstable isotopes and transient resonances, to mutual dynamical interactions during a reaction. Here, after an overview of these phenomena, basic formalism and concepts in ab initio nuclear structure and reaction theory are introduced, followed by selected examples of how powerful techniques for the description of the nuclear force and the solution of the nuclear quantum many-body problem combined with cutting-edge high-performance computing are enabling a comprehensive understanding of light exotic nuclei and thermonuclear reactions that play an important role in the Cosmos and on Earth.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2014

Line coupling effects in the isotropic Raman spectra of N{sub 2}: A quantum calculation at room temperature

We present quantum calculations of the relaxation matrix for the Q branch of N{sub 2} at room temperature using a recently proposed N{sub 2}-N{sub 2} rigid rotor potential. Close coupling calculations were complemented by coupled states studies at high energies and provide about 10 200 two-body state-to state cross sections from which the needed one-body cross-sections may be obtained. For such temperatures, convergence has to be thoroughly analyzed since such conditions are close to the limit of current computational feasibility. This has been done using complementary calculations based on the energy corrected sudden formalism. Agreement of these quantum predictions with experimental data is good, but the main goal of this work is to provide a benchmark relaxation matrix for testing more approximate methods which remain of a great utility for complex molecular systems at room (and higher) temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Density-functional theory for systems with noncollinear spin: Orbital-dependent exchange-correlation functionals and their application to the Hubbard dimer

Here, a class of orbital-dependent exchange-correlation potentials for applications in noncollinear spin-density-functional theory is developed. Starting from the optimized effective potential formalism for the exact exchange potential—generalized to the noncollinear case—correlation effects are added via a self-consistent procedure inspired by the Singwi-Tosi-Land-Sjölander method. The orbital-dependent exchange-correlation potentials are applied to the Hubbard dimer in uniform and noncollinear magnetic fields and compared to exact diagonalization and to the Bethe-ansatz local spin-density approximation. The Singwi-Tosi-Land-Sjölander method gives the overall best performance for total energies, densities, and magnetizations, particularly in the weakly to moderately correlated regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Learning new physics from data: A symmetrized approach

Thousands of person years have been invested in searches for new physics (NP), the majority of them motivated by theoretical considerations. Yet, no evidence of beyond the Standard Model physics has been found. This suggests that model-agnostic searches might be an important key to explore NP, and help discover unexpected phenomena which can inspire future theoretical developments. A possible strategy for such searches is identifying asymmetries between data samples that are expected to be symmetric within the Standard Model. We propose exploiting neural networks (NNs) to quickly fit and statistically test the differences between two samples. Our method is based on an earlier work, originally designed for inferring the deviations of an observed dataset from that of a much larger reference dataset. We present a symmetric formalism, generalizing the original one, avoiding fine-tuning of the NN parameters and any constraints on the relative sizes of the samples. Our formalism could be used to detect small symmetry violations, extending the discovery potential of current and future particle physics experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

ARCS: Agentic Retrieval-Augmented Code Synthesis with Iterative Refinement

Agentic Retrieval-Augmented Code Synthesis with Iterative RefinementIn supercomputing, efficient and optimized code generation is essential to leverage high-performance systems effectively. We have developed Agentic Retrieval-Augmented Code Synthesis (ARCS), an advanced framework for accurate, robust, and efficient code generation, completion, and translation. ARCS integrates Retrieval-Augmented Generation (RAG) with Chain-of-Thought (CoT) reasoning to systematically break down and iteratively refine complex programming tasks. An agent-based RAG mechanism retrieves relevant code snippets, while real-time execution feedback drives the synthesis of candidate solutions. This process is formalized as a state-action search tree optimization, balancing code correctness with editing efficiency. Evaluations on the Geeks4Geeks and HumanEval benchmarks demonstrate that ARCS significantly outperforms traditional prompting methods in translation and generation quality. By enabling scalable and precise code synthesis, ARCS offers transformative potential for automating and optimizing code development in supercomputing applications, enhancing computational resource utilization

Bhattarai, Manish [Los Alamos National Labs]↗

DOE OSTI2021

Anomaly Attribution with Likelihood Compensation

This paper addresses the task of explaining anomalous predictions of a black-box regression model. When using a black-box model, such as one to predict building energy consumption from many sensor measurements, we often have a situation where some observed samples may significantly deviate from their prediction. It may be due to a sub-optimal black-box model, or simply because those samples are outliers. In either case, one would ideally want to compute a responsibility score indicative of the extent to which an input variable is responsible for the anomalous output. In this work, we formalize this task as a statistical inverse problem: Given model deviation from the expected value, infer the responsibility score of each of the input variables. We propose a new method called likelihood compensation (LC), which is founded on the likelihood principle and computes a correction to each input variable. To the best of our knowledge, this is the first principled framework that computes a responsibility score for real valued anomalous model deviations. We apply our approach to a real-world building energy prediction task and confirm its utility based on expert feedback.

Idé, Tsuyoshi↗

DOE OSTI2013

Cobalt-iron red-ox behavior in nanostructured La{sub 0.4}Sr{sub 0.6}Co{sub 0.8}Fe{sub 0.2}O{sub 3-{delta}} cathodes

Nano-sized La{sub 0.4}Sr{sub 0.6}Co{sub 0.8}Fe{sub 0.2}O{sub 3-{delta}} (LSCF) perovskite samples (prepared by a conventional acetate route and a novel acetate synthesis with HMTA additives), were tested simulating a red-ox cycle. The crystallography was studied by X-ray Powder Diffraction (XPD) and the changes in the oxidation state of the perovskite B-site were evaluated by synchrotron X-ray Absorption Near Edge Spectroscopy (XANES). After a reducing treatment, LSFC particles show the appearance of a new phase that coexists with the original one. The structural change is accompanied by a Co and Fe formal oxidation states decrease, although Fe remains always closer to 4+ and Co closer to 3+. The treatment produces a B-site valence average reduction from 3.52+ to 3.26+ and the formation of oxygen vacancies. A re-oxidation treatment under O{sub 2} rich atmosphere at 800 Degree-Sign C for 10 h shows that the change is reversible and independent of the two chemical methods used to synthesize the LSCF nano-particles. - Graphical abstract: XANES and XPD measurements in nanostructured LSCF before (black) and after (red/green) a red/ox cycle. Highlights: Black-Right-Pointing-Pointer Red-ox treatments in LSCF nano-particles cause a reversible reaction. Black-Right-Pointing-Pointer XPD analyses show that a new 'reduced' phase coexist with the oxidize one. Black-Right-Pointing-Pointer The B-site formal oxidation state decreases and the {delta} increases upon reduction. Black-Right-Pointing-Pointer Fe remains in a higher valence (closer to 4+) than Co (close to 3+). Black-Right-Pointing-Pointer The behavior seems to be independent of the synthesis method used.

36 MATERIALS SCIENCE↗

DOE OSTI2014

Method for discovering relationships in data by dynamic quantum clustering

Data clustering is provided according to a dynamical framework based on quantum mechanical time evolution of states corresponding to data points. To expedite computations, we can approximate the time-dependent Hamiltonian formalism by a truncated calculation within a set of Gaussian wave-functions (coherent states) centered around the original points. This allows for analytic evaluation of the time evolution of all such states, opening up the possibility of exploration of relationships among data-points through observation of varying dynamical-distances among points and convergence of points into clusters. This formalism may be further supplemented by preprocessing, such as dimensional reduction through singular value decomposition and/or feature filtering.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Method for discovering relationships in data by dynamic quantum clustering

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Electrochemical Formation of Li-M-(M')-Si Phases Using Multivalent Electrolyte Salt Additives

Lithium-rich silicides (Li 15 Si 4 ), formed during the electrochemical lithiation of silicon, show high reactivity with electrolyte components that contribute to capacity decay, formal lithium loss, and low coulombic efficiency. Recently, the reactivity of lithium silicides was found to be suppressed by substituting a multivalent cation (i.e. Mg, Ca) for lithium that results in the room temperature formation of a ternary Li-M-Si phase. In this study, we explored a range of multivalent electrolyte salt additives (M = Ni, Cu, La, Ce, Sr, Ba, and Ca-Mg mixed salt) in a lithium-ion cell configuration and identified a room temperature electrochemical route to the formation of new ternary and quaternary lithium silicides. Using this method, both nickel and copper salts were found to plate onto the silicon electrode surface upon lithiation. Further, based on refined synchrotron XRD data, multivalent cations with an ionic radius similar to Na (~1.03 Å) or smaller can be inserted electrochemically into a formally cation-deficient Li 15 Si 4 host lattice to form new ternary (or quartenary) phases. The electrochemical synthesis of a new quaternary Li-M-M’-Si phase represents a facile route to preparing and scaling materials isostructural to the Heusler phase and electron-precise Li 14 MgSi 4 phase that results in enhanced cycling and calendar life performance.

25 ENERGY STORAGE↗

DOE OSTI2014

{sup 103}Pd strings: Monte Carlo assessment of a new approach to brachytherapy source design

Purpose: A new type of{sup 103}Pd source (CivaString and CivaThin by CivaTech Oncology, Inc.) is examined. The source contains {sup 103}Pd and Au radio-opaque marker(s), all contained within low-Z{sub eff} organic polymers that permit source flexibility. The CivaString source is available in lengths L of 10, 20, 30, 40, 50, and 60 mm, and referred to in the current study as CS10–CS60, respectively. A thinner design, CivaThin, has sources designated as CT10–CT60, respectively. The CivaString and CivaThin sources are 0.85 and 0.60 mm in diameter, respectively. The source design is novel and offers an opportunity to examine its interesting dosimetric properties in comparison to conventional {sup 103}Pd seeds. Methods: The MCNP5 radiation transport code was used to estimate air-kerma rate and dose rate distributions with polar and cylindrical coordinate systems. Doses in water and prostate tissue phantoms were compared to determine differences between the TG-43 formalism and realistic clinical circumstances. The influence of Ti encapsulation and 2.7 keV photons was examined. The accuracy of superposition of dose distributions from shorter sources to create longer source dose distributions was also assessed. Results: The normalized air-kerma rate was not highly dependent onL or the polar angle θ, with results being nearly identical between the CivaString and CivaThin sources for common L. The air-kerma strength was also weakly dependent on L. The uncertainty analysis established a standard uncertainty of 1.3% for the dose-rate constant Λ, where the largest contributors were μ{sub en}/ρ and μ/ρ. The Λ values decreased with increasing L, which was largely explained by differences in solid angle. The radial dose function did not substantially vary among the CivaString and CivaThin sources for r ≥ 1 cm. However, behavior for r < 1 cm indicated that the Au marker(s) shielded radiation for the sources having L = 10, 30, and 50 mm. The 2D anisotropy function exhibited peaks and valleys that corresponded to positions adjacent to {sup 103}Pd wells and Au markers, respectively. Dose distributions of both source types had minimal anisotropy in comparison to conventional {sup 103}Pd seeds. Contributions by 2.7 keV photons comprised ≤0.1% of the dose from all photons at positions farther than 0.13 mm from the polymer source surface. Differences between absorbed dose to water and prostate became more substantial as distance from the sources increased, with prostate dose being about 13% lower for r = 5 cm. Using a cylindrical coordinate system, dose superposition of small length sources to replicate the dose distribution for a long length source proved to be a robust technique; a 2.0% tolerance compared with the reference dose distribution did not exceed 0.1 cm{sup 3} for any of the examined source combinations. Conclusions: By design, the CivaString and CivaThin sources have novel dosimetric characteristics in comparison to Ti-encapsulated{sup 103}Pd seeds. The dosimetric characterization has determined the reasons for these differences through analysis using Monte Carlo-based radiation transport simulations.

62 RADIOLOGY AND NUCLEAR MEDICINE↗