Engineering Papers⌕ Search

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.

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

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Unified formalism for Thermal Quantum Field Theories: A geometric viewpoint

In this paper we study a unified formalism for Thermal Quantum Field Theories, i.e., for the Matsubara approach, Thermo Field Dynamics and the Path Ordered Method. To do so, we employ a mechanism akin to the Hawking effect which explores a relationship between the concept of temperature and spacetimes endowed with event-horizons. In particular, we consider an eight dimensional static spacetime, the so-called η–ξ spacetime, which we show to form an appropriate geometric background for generic Thermal Quantum Field Theories. Within this framework, the different formalisms of Thermal Field Theory are unified in a very natural way via various analytical continuations and the set of time-paths used in the Path Ordered Method is interpreted in geometric terms. We also explain reported inconsistencies inherent in the Thermo Field Dynamics through the appearance of horizons (and ensuing loss of information) in the η–ξ spacetime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

SU-C-201-06: Small Field Correction Factors for the MicroDiamond Detector in the Gamma Knife-Model C Derived Using Monte Carlo Methods

Purpose: To determine small field correction factors for PTW’s microDiamond detector in Elekta’s Gamma Knife Model-C unit. These factors allow the microDiamond to be used in QA measurements of output factors in the Gamma Knife Model-C; additionally, the results also contribute to the discussion on the water equivalence of the relatively-new microDiamond detector and its overall effectiveness in small field applications. Methods: The small field correction factors were calculated as k correction factors according to the Alfonso formalism. An MC model of the Gamma Knife and microDiamond was built with the EGSnrc code system, using BEAMnrc and DOSRZnrc user codes. Validation of the model was accomplished by simulating field output factors and measurement ratios for an available ABS plastic phantom and then comparing simulated results to film measurements, detector measurements, and treatment planning system (TPS) data. Once validated, the final k factors were determined by applying the model to a more waterlike solid water phantom. Results: During validation, all MC methods agreed with experiment within the stated uncertainties: MC determined field output factors agreed within 0.6% of the TPS and 1.4% of film; and MC simulated measurement ratios matched physically measured ratios within 1%. The final k correction factors for the PTW microDiamond in the solid water phantom approached unity to within 0.4%±1.7% for all the helmet sizes except the 4 mm; the 4 mm helmet size over-responded by 3.2%±1.7%, resulting in a k factor of 0.969. Conclusion: Similar to what has been found in the Gamma Knife Perfexion, the PTW microDiamond requires little to no corrections except for the smallest 4 mm field. The over-response can be corrected via the Alfonso formalism using the correction factors determined in this work. Using the MC calculated correction factors, the PTW microDiamond detector is an effective dosimeter in all available helmet sizes. The authors would like to thank PTW (Friedberg, Germany) for providing the PTW microDiamond detector for this research.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

Turning statistical physics models into materials design engines

Despite the success statistical physics has enjoyed at predicting the properties of materials for given parameters, the inverse problem, identifying which material parameters produce given, desired properties, is only beginning to be addressed. Recently, several methods have emerged across disciplines that draw upon optimization and simulation to create computer programs that tailor material responses to specified behaviors. However, so far the methods developed either involve black-box techniques, in which the optimizer operates without explicit knowledge of the material's configuration space, or require carefully tuned algorithms with applicability limited to a narrow subclass of materials. Here we introduce a formalism that can generate optimizers automatically by extending statistical mechanics into the realm of design. The strength of this approach lies in its capability to transform statistical models that describe materials into optimizers to tailor them. Here, by comparing against standard black-box optimization methods, we demonstrate how optimizers generated by this formalism can be faster and more effective, while remaining straightforward to implement. The scope of our approach includes possibilities for solving a variety of complex optimization and design problems concerning materials both in and out of equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fully self-consistent finite-temperature $GW$ in Gaussian Bloch orbitals for solids

In this work, we present algorithmic and implementation details for the fully self-consistent finite-temperature $GW$ method in Gaussian Bloch orbitals for solids. Our implementation is based on the finite-temperature Green's function formalism in which all equations are solved on the imaginary axis, without resorting to analytical continuation during the self-consistency. No quasiparticle approximation is employed and all matrix elements of the self-energy are explicitly evaluated. The method is tested by evaluating the band gaps of selected semiconductors and insulators. We show agreement with other, differently formulated, finite-temperature sc ⁢$GW$ implementations when finite-size corrections and basis-set errors are taken into account. By migrating computationally intensive calculations to graphics processing units, we obtain scalable results on large supercomputers with nearly optimal performance. Our work demonstrates the applicability of Gaussian orbital based sc⁢ $GW$ for ab initio correlated material simulations and provides a sound starting point for embedding methods built on top of $GW$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Stochastic fluctuations in relativistic fluids: Causality, stability, and the information current

We develop a general formalism for introducing stochastic fluctuations around thermodynamic equilibrium which takes into account, for the first time, recent developments in the causality and stability properties of relativistic hydrodynamic theories. The method is valid for any covariantly stable theory of relativistic viscous fluid dynamics derived from a covariant maximum entropy principle. We illustrate the formalism with some applications, showing how it could be used to consistently introduce fluctuations in a model of relativistic heat diffusion and in conformally invariant Israel-Stewart theory in a general hydrodynamic frame. Furthermore, the latter example is used to study the hydrodynamic frame dependence of the symmetric two-point function of fluctuations of the energy-momentum tensor.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Bayesian tomography and integrated data analysis in fusion diagnostics

In this article, a Bayesian tomography method using non-stationary Gaussian process for a prior has been introduced. The Bayesian formalism allows quantities which bear uncertainty to be expressed in the probabilistic form so that the uncertainty of a final solution can be fully resolved from the confidence interval of a posterior probability. Moreover, a consistency check of that solution can be performed by checking whether the misfits between predicted and measured data are reasonably within an assumed data error. In particular, the accuracy of reconstructions is significantly improved by using the non-stationary Gaussian process that can adapt to the varying smoothness of emission distribution. The implementation of this method to a soft X-ray diagnostics on HL-2A has been used to explore relevant physics in equilibrium and MHD instability modes. This project is carried out within a large size inference framework, aiming at an integrated analysis of heterogeneous diagnostics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Perspective on Kramers symmetry breaking and restoration in relativistic electronic structure methods for open-shell systems

Without rigorous symmetry constraints, solutions to approximate electronic structure methods may artificially break symmetry. In the case of the relativistic electronic structure, if time-reversal symmetry is not enforced in calculations of molecules not subject to a magnetic field, it is possible to artificially break Kramers degeneracy in open shell systems. This leads to a description of excited states that may be qualitatively incorrect. Despite this, different electronic structure methods to incorporate correlation and excited states can partially restore Kramers degeneracy from a broken symmetry solution. For single-reference techniques, the inclusion of double and possibly triple excitations in the ground state provides much of the needed correction. Formally, however, this imbalanced treatment of the Kramers-paired spaces is a multi-reference problem, and so methods such as complete-active-space methods perform much better at recovering much of the correct symmetry by state averaging. Using multi-reference configuration interaction, any additional corrections can be obtained as the solution approaches the full configuration interaction limit. A recently proposed “Kramers contamination” value is also used to assess the magnitude of symmetry breaking.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Bayesian Monte Carlo Evaluation of Imperfect (n, 233 U) Data and Model

Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Exact stochastic unraveling of an optical coherence dynamics by cumulant expansion

A numerically exact Monte Carlo scheme for calculation of open quantum system dynamics is proposed and implemented. The method consists of a Monte Carlo summation of a perturbation expansion in terms of trajectories in Liouville phase-space with respect to the coupling between the excited states of the molecule. The trajectories are weighted by a complex decoherence factor based on the second-order cumulant expansion of the environmental evolution. The method can be used with an arbitrary environment characterized by a general correlation function and arbitrary coupling strength. It is formally exact for harmonic environments, and it can be used with arbitrary temperature. Time evolution of an optically excited Frenkel exciton dimer representing a molecular exciton interacting with a charge transfer state is calculated by the proposed method. We calculate the evolution of the optical coherence elements of the density matrix and linear absorption spectrum, and compare them with the predictions of standard simulation methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Fermionic mean-field dynamics for spin systems beyond free fermions

We introduce the fermionized time-dependent Hartree–Fock (fTDHF), a real-time quantum dynamics method for spin-1/2 Hamiltonians following their mapping to fermions via the Jordan-Wigner transformation. fTDHF is formally equivalent to exact dynamics in the case of free fermions, and can efficiently handle non-local string operators arising from long-range interactions via transition matrix elements between non-orthogonal Slater determinants. We show that the fTDHF method can be implemented on a classical computer with a cost that scales polynomially with system size, and linearly with the time steps. We benchmark fTDHF against exact dynamics on three separate spin-1/2 models, representing adiabatic preparation of states with long-range correlations, disorder-driven observation of many-body localization, and particle production in the Schwinger model. For each of these systems, fTDHF is shown to reproduce the qualitative dynamics generated by the exact evolutions, while maintaining a simple physical picture due to its mean-field nature.

Dutta, Rishab↗
DOE OSTI2014

From the Weyl quantization of a particle on the circle to number–phase Wigner functions

A generalized Weyl quantization formalism for a particle on the circle is shown to supply an effective method for defining the number–phase Wigner function in quantum optics. A Wigner function for the state ϱ{sup ^} and the kernel K for a particle on the circle is defined and its properties are analysed. Then it is shown how this Wigner function can be easily modified to give the number–phase Wigner function in quantum optics. Some examples of such number–phase Wigner functions are considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Spherical coupled-cluster theory for open-shell nuclei

A microscopic description of nuclei is important to understand the nuclear shell model from fundamental principles. This is difficult to achieve for more than the lightest nuclei without an effective approximation scheme. Define and evaluate an approximation scheme that can be used to study nuclei that are described as two particles attached to a closed (sub-)shell nucleus. The equation-of-motion coupled-cluster formalism has been used to obtain ground- and excited-state energies. This method is based on the diagonalization of a non-Hermitian matrix obtained from a similarity transformation of the many-body nuclear Hamiltonian. A chiral interaction at the next-to-next-to-next-to leading order (N 3 LO) using a cutoff at 500 MeV was used. The ground-state energies of 6 Li and 6 He were in good agreement with a no-core shell-model calculation using the same interaction. Several excited states were also produced with overall good agreement. Only the $J^π = 3^+$ excited state in 6 Li showed a sizable deviation. The ground-state energies of 18 O, 18 F, and 18 Ne were converged but underbound compared to experiment. Moreover, the calculated spectra were converged and comparable to both experiment and shell-model studies in this region. Some excited states in 18 O were high or missing in the spectrum. It was also shown that the wave function for both ground and excited states separates into an intrinsic part and a Gaussian for the center-of-mass coordinate. Spurious center-of-mass excitations are clearly identified. Results are converged with respect to the size of the model space and the method can be used to describe nuclear states with simple structure. Especially the ground-state energies were very close to what has been achieved by exact diagonalization. To obtain a closer match with experimental data, effects of three-nucleon forces, the scattering continuum, as well as additional configurations in the coupled-cluster approximations are necessary.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

On the question of area weighting in cell-centered hydrodynamics

This work presents the theory behind a cell-centered Lagrangian formulation (CCH) in axial (RZ) geometry suitable for solid materials on unstructured polyhedral cells. There are two common discretizations of the 2D RZ equations. It is widely believed the finite volume equations in RZ are conservative, but do not preserve rotational symmetry. Conversely, it is also widely believed that the area weighted equations are flawed because they preserve symmetry at the cost of rigorous conservation. In this work, we show that both conservation and symmetry preservation are possible. We present a first principles derivation of the discrete solid dynamics equations. We employ an adaptation of the FV method to show that this formulation not only preserves symmetry, but is indeed formally conservative. We further show that these equations are appropriate for an infinitesimal control volume at the centroid of the cell. The resulting equations are closely related to the canonical area-weighted ones. This analysis lends theoretical support to the fluid formulation of Maire. Although the principal focus of the paper is on the Lagrange formulation, symmetry properties of advective terms appropriate to Eulerian and arbitrary Lagrange-Euler (ALE) formulations are also considered.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Symposium: Quantum Theory, Quantum Materials, Quantum Computing (Final Report, DE-SC0022551)

Quantum mechanics is fundamental to understanding, predicting, and controlling many physical processes in a vast range of scientific and technological areas. In molecules and materials, the quantum mechanics of Coulombic many-body systems determines far more than the ground state and bulk properties. Quantum mechanics governs magnetism, electron and spin transport, chemical reactions and catalysis, superconductivity, photosynthesis, etc. Each plays a vital role in modern society. As digital technologies push both temporal and length-scale limits, fundamental quantum principles become the current focus for a new class of technologies: quantum computing, quantum sensing, quantum encryption, etc. This oncoming quantum information sciences wave is manifest recently in frontier foci of all six Department of Energy Office of Science programs. The Sanibel Symposium offers a proven, influential platform for facilitating and advancing cross-fertilization among theoretical and computational chemists and physicists, computer scientists, and materials engineers. The Symposia have a long, successful record of driving progress on theory and computation of materials, nanostructures, and their molecular constituents. The Symposia thus are well-suited for advancing the quantum information science thrust in the Office of Basic Energy Sciences (BES). With that QIS thrust as the focus, we propose five sessions of the 61st Sanibel Symposium on the theme Quantum Theory, Quantum Materials, Quantum Computing. From a longer perspective, this thematic program for 2022 is the second year of a five year plan of program segments addressed at challenges in quantum information sciences and quantum materials. The rest of the Symposium program will be of considerable relevance, with, for example, sessions on quantum molecular and material dynamics, and new directions in coupled cluster theory. New ideas for methods and applications can be expected to emerge from the stimulus of the formal presentations and informal discussions. The proposed award is targeted at students, post-docs, and junior scientists. Conscious effort, characteristic of the Sanibel Symposia, to stimulate the participation of women and under-represented minorities will continue.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Security Self-Assessment Toolkit for Nuclear Materials Facilities: Focus on Insider Threat Mitigation

Theft or sabotage of weapons-usable nuclear materials is a global concern. To minimize this threat, establishing and maintaining an effective nuclear security regime is required to protect against criminal or other negligent acts. Use of a formalized insider threat mitigation program is one such security measure. Individuals who have or held authorized access to an organization's critical assets, such as nuclear materials, are considered "insiders." Insider threats, or insider adversaries, are motivated individuals who possess access, authority, and knowledge to conduct a malicious act or facilitate that of an external party. To thwart insider threats (both intentional and unintentional), organizations can formalize an enterprise-wide approach to identify and mitigate the unique risks presented by insiders. This report provides an approach to evaluate an insider threat mitigation program at facilities with nuclear materials. Formal program evaluations serve many purposes and can be designed using several different methods and techniques. This report presents a self-assessment approach to program evaluation whereby an organization can assess its strengths, identify key gaps, and set priorities for ongoing improvement efforts to mitigate insider threats. Results of the self-assessment can provide critical information to contribute to the continuous improvement of an organization’s insider threat mitigation program within eight specific domain areas.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2026

Comparative Assessment of Battery Carbon Footprint Calculation Frameworks to Support U.S. Battery Manufacturing

A “battery passport” is a digital record that provides comprehensive information about an individual battery across its life cycle. This concept has been introduced by Battery Regulation (EU) 2023/1542 [1], which applies to batteries for electric vehicles (EVs), industrial batteries with a capacities greater than 2 kWh, and batteries for light means of transport (LMT) greater than 2kWh sold in the European Union (EU) market – impacting manufacturers and exporters across multiple jurisdictions globally. Among its reporting mandates, a key feature of battery passports is the requirement for a carbon footprint (CF) calculation methodology supported by enhanced data granularity to enable traceable and verifiable CF results. While no other jurisdictions have yet adopted formal battery passport requirements like the EU’s, some are developing CF calculation methods to comply with the EU Battery Regulation or are creating CF-related regulations that could evolve in a similar direction, reflecting the growing importance of battery CF guidelines for manufacturers seeking to remain competitive in global markets.

Zhang, Jingyi [Argonne National Laboratory (ANL), ↗
DOE OSTI2015

Application of a New Ensemble Conserving Quantum Dynamics Simulation Algorithm to Liquid para-Hydrogen and ortho-Deuterium

Here, we apply the Feynman-Kleinert Quasi-Classical Wigner (FK-QCW) method developed in our previous work [Smith et al., J. Chem. Phys. 142, 244112 (2015)] for the determination of the dynamic structure factor of liquid para-hydrogen and ortho-deuterium at state points of (T = 20.0 K, n = 21.24 nm -3 ) and (T = 23.0 K, n = 24.61 nm -3 ), respectively. When applied to this challenging system, it is shown that this new FK-QCW method consistently reproduces the experimental dynamic structure factor reported by Smith et al. [J. Chem. Phys. 140, 034501 (2014)] for all momentum transfers considered. Moreover, this shows that FK-QCW provides a substantial improvement over the Feynman-Kleinert linearized path-integral method, in which purely classical dynamics are used. Furthermore, for small momentum transfers, it is shown that FK-QCW provides nearly the same results as ring-polymer molecular dynamics (RPMD), thus suggesting that FK-QCW provides a potentially more appealing algorithm than RPMD since it is not formally limited to correlation functions involving linear operators.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Projection-operator methods for classical transport in magnetized plasmas. Part 2. Nonlinear response and the Burnett equations [Projection-operator methods for classical transport in magnetized plasmas. Part 1. Nonlinear response and the Burnett equations]

The time-independent projection-operator formalism for the derivation of Burnett equations is extended and considered in the context of multispecies and magnetized plasmas. The procedure provides specific formulas for transport coefficients in terms of two-time correlation functions involving both two and three phase-space points. It is shown how to calculate those correlation functions in the limit of weak coupling. The results are used to demonstrate, with the aid of a particular non-trivial example, that the Chapman–Enskog methodology employed by Catto & Simakov (CS) (Phys. Plasmas, vol. 11, 2004, pp. 90–102) to calculate the contributions to the parallel viscosity driven by temperature gradients is consistent with formulas previously derived from the two-time formalism by Brey (J. Chem. Phys., vol. 79, 1983, pp. 4585–4598). Here, the work serves to unify previous work on plasma kinetic theory with formalism usually applied to turbulence. Additional contributions include discussions of (i) Braginskii-order interspecies momentum exchange from the point of view of two-time correlations; and (ii) a simple stochastic model, unrelated to many-body theory, that exhibits Burnett effects. Insights from that model emphasize the role of non-Gaussian statistics in the evaluation of Burnett transport coefficients, including the effects calculated by CS that stem from the nonlinear collision operator. Together, Parts 1 and 2 of this series provide an introduction to projection-operator methods that should be broadly useful in theoretical plasma physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗