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 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Electron impact ionization in plasma technologies; studies on atomic boron and BN molecule

Electron impact ionization plays important role in plasma technologies. Relevant cross sections on atomic boron are required to understand the erosion processes in fusion experiments. Boronization of plasma exposed surfaces of tokomaks has proved to be an effective way to produce very pure fusion plasmas. This paper reports comprehensive theoretical investigations on electron scattering with atomic Boron and Boron Nitride in solid phases. Presently we determine total ionization cross-section Q{sub ion} and the summed-electronic excitation cross section ΣQ{sub exc} in a standard quantum mechanical formalism called SCOP and CSP-ic methods. Our calculated cross sections are examined as functions of incident electron energy along with available comparisons.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Atomic orbital-based SOS-MP2 with tensor hypercontraction. I. GPU-based tensor construction and exploiting sparsity

We present a tensor hypercontracted (THC) scaled opposite spin second order Møller-Plesset perturbation theory (SOS-MP2) method. By using THC, we reduce the formal scaling of SOS-MP2 with respect to molecular size from quartic to cubic. We achieve further efficiency by exploiting sparsity in the atomic orbitals and using graphical processing units (GPUs) to accelerate integral construction and matrix multiplication. The practical scaling of GPU-accelerated atomic orbital-based THC-SOS-MP2 calculations is found to be N{sup 2.6} for reference data sets of water clusters and alanine polypeptides containing up to 1600 basis functions. The errors in correlation energy with respect to density-fitting-SOS-MP2 are less than 0.5 kcal/mol for all systems tested (up to 162 atoms).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

Is sodium a superconductor under high pressure?

Superconductivity has been predicted or measured for most alkali metals under high pressure, but the computed critical temperature (T c ) of sodium (Na) at the face-centered cubic (fcc) phase is vanishingly low. Here we report a thorough, first-principles investigation of superconductivity in Na under pressures up to 260 GPa, where the metal-to-insulator transition occurs. Linear-response calculations and density functional perturbation theory were employed to evaluate phonon distributions and the electron-phonon coupling for bcc, fcc, cI16, and tI19 Na. Our results indicate that the maximum electron-phonon coupling parameter, λ, is 0.5 for the cI16 phase, corresponding to a theoretical peak in the critical temperature at T c ≈ 1.2 K. When pressure decreases or increases from 130 GPa, T c drops quickly. This is mainly due to the lack of p-d hybridization in Na even at 260 GPa. Since current methods based on the Eliashberg and McMillian formalisms tend to overestimate the Tc (especially the peak values) of alkali metals, we conclude that under high pressure—before the metal-to-insulator transition at 260 GPa—superconductivity in Na is very weak, if it is measurable at all.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2016

BRST quantization of cosmological perturbations

BRST quantization is an elegant and powerful method to quantize theories with local symmetries. In this article we study the Hamiltonian BRST quantization of cosmological perturbations in a universe dominated by a scalar field, along with the closely related quantization method of Dirac. We describe how both formalisms apply to perturbations in a time-dependent background, and how expectation values of gauge-invariant operators can be calculated in the in-in formalism. Our analysis focuses mostly on the free theory. By appropriate canonical transformations we simplify and diagonalize the free Hamiltonian. BRST quantization in derivative gauges allows us to dramatically simplify the structure of the propagators, whereas Dirac quantization, which amounts to quantization in synchronous gauge, dispenses with the need to introduce ghosts and preserves the locality of the gauge-fixed action.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Light in the dark forest. Part I. An efficient optimal estimator for 3D Lyman-alpha forest power spectrum

The highly anisotropic nature of the Lyman-alpha (Lyα) forest data introduces a complex survey window function that complicates the measurement of the three-dimensional power spectrum ( P 3D ). In this paper, we present the first fully optimal estimator for P 3D , which exactly deconvolves the survey window function and marginalizes contaminated modes that distort the power spectrum. Our approach adapts optimal estimator techniques developed for the 2D cosmic microwave background data to the 3D case. To achieve computational feasibility, we employ the conjugate gradient method and implement the P 3 M formalism to handle large-scale and small-scale operations separately and efficiently. We validate our estimator using Monte Carlo mocks and Gaussian simulations, demonstrating its accuracy and computational efficiency. We confirm that mode marginalization eliminates distortions arising from quasar continuum errors and delivers robust power spectrum estimation, though it also inflates errors at large scales. This first implementation works in the flat-sky case; we discuss the remaining steps needed to generalize it to the curved-sky case. This formalism offers a foundation for the Lyα forest P 3D measurements and a new path toward cosmological constraints from the Lyα forest data.

Lyman alpha forest↗

DOE OSTI2022

Nonequilibrium DMFT + CPA for correlated disordered systems

We present a solution for the nonequilibrium dynamics of an interacting disordered system. The approach adapts the combination of the equilibrium dynamical mean-field theory and the equilibrium coherent potential approximation methods to the nonequilibrium many-body formalism, using the Kadanoff-Baym-Keldysh complex time contour, for the dynamics of interacting disordered systems away from equilibrium. We use our time domain solution to obtain the equilibrium density of states of the disordered interacting system described by the Anderson-Hubbard model, bypassing the necessity for the cumbersome analytical continuation process. We further apply the nonequilibrium solution to the interaction quench problem for an isolated disordered system. Here, the interaction is abruptly changed from zero (noninteracting system) to another constant (finite) value at which it is subsequently kept. We observe via the time dependence of the potential, kinetic, and total energies the effect of disorder on the relaxation of the system as a function of final interaction strength. Finally, the real-time approach has the potential to shed light on the fundamental role of disorder in the nonequilibrium dynamics of interacting quantum systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Charge self-consistent density functional theory plus ghost rotationally invariant slave-boson theory for correlated materials

We present a charge self-consistent density functional theory combined with the ghost rotationally invariant slave-boson (DFT+gRISB) formalism for studying correlated materials. Here, this method is applied to SrVO 3 and NiO, representing prototypical correlated metals and charge-transfer insulators. For SrVO 3 , we demonstrate that DFT+gRISB yields an accurate equilibrium volume and effective mass close to experimentally observed values. Regarding NiO, DFT+gRISB enables the simultaneous description of charge-transfer and Mott-Hubbard bands, significantly enhancing the accuracy of the original DFT+RISB approach. Furthermore, the calculated equilibrium volume and spectral function reasonably agree with experimental observations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Secondary γ-ray decays from the partial-wave $\textit{T}$ matrix with an $\textit{R}$-matrix application to 15 N($p, α_1γ$) 12 C

The secondary γ rays emitted following a nuclear reaction are frequently relatively straightforward to detect experimentally. Despite the large volume of such data, a practical formalism for describing these γ rays in terms of partial-wave $\textit{T}$-matrix elements has never been given. The partial-wave formalism is applicable when $\textit{R}$-matrix methods are used to describe the reaction in question. This paper supplies the needed framework, and it is demonstrated by the application to the 15 N($p, α_1γ$) 12 C reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Resilience Metrics for the Electric Power System: A Performance-Based Approach.

Grid resilience is a concept related to a power system's ability to continue operating and delivering power even in the event that low probability, high-consequence disruptions such as hurricanes, earthquakes, and cyber-attacks occur. Grid resilience objectives focus on managing and, ideally, minimizing potential consequences that occur as a result of these disruptions. Currently, no formal grid resilience definitions, metrics, or analysis methods have been universally accepted. This document describes an effort to develop and describe grid resilience metrics and analysis methods. The metrics and methods described herein extend upon the Resilience Analysis Process (RAP) developed by Watson et al. for the 2015 Quadrennial Energy Review. The extension allows for both outputs from system models and for historical data to serve as the basis for creating grid resilience metrics and informing grid resilience planning and response decision-making. This document describes the grid resilience metrics and analysis methods. Demonstration of the metrics and methods is shown through a set of illustrative use cases.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2018

Thermodynamic, Kinetic, and Physical Properties of Nb-U (Niobium-Uranium)

Based on an improved coupling between ab initio energetics output, phenomenological thermodynamics based on the CALPHAD (CALculation of Phase Diagrams) approach, and the use of the Thermo-Calc software, together with experimental data (whenever available), a thermodynamic assessment of the Nb-U alloy system was carried out. We made use of the full array of modern electronic-structure techniques based on density functional theory, including ab initio fully relativistic full-potential linear-muffin orbital method (FR-FP-LMTO), methods, the Scalar Relativistic Korringa-Kohn-Rostoker Atomic-Sphere Approximation (SR-KKR-ASA), and the fully relativistic or scalar relativistic Exact muffin-orbital (SR- or FR-EMTO) method, both based on the Green function formalism that has the advantage of being able to treat chemically random configuration (i.e., solid solution) with the so-called coherent potential approximation (CPA). All electronic structure calculations were performed within the Generalized Gradient Approximation (GGA) for the treatment of electron correlations. The results from this technology serve as input for CALPHAD assessment of phase diagrams

36 MATERIALS SCIENCE↗

DOE OSTI2020

A Review of the Serrated-Flow Phenomenon and Its Role in the Deformation Behavior of High-Entropy Alloys

High-entropy alloys (HEAs) are a novel class of alloys that have many desirable properties. The serrated flow that occurs in high-entropy alloys during mechanical deformation is an important phenomenon since it can lead to significant changes in the microstructure of the alloy. In this article, we review the recent findings on the serration behavior in a variety of high-entropy alloys. Relationships among the serrated flow behavior, composition, microstructure, and testing condition are explored. Importantly, the mechanical-testing type (compression/tension), testing temperature, applied strain rate, and serration type for certain high-entropy alloys are summarized. The literature reveals that the serrated flow can be affected by experimental conditions such as the strain rate and test temperature. Furthermore, this type of phenomenon has been successfully modeled and analyzed, using several different types of analytical methods, including the mean-field theory formalism and the complexity-analysis technique. Importantly, the results of the analyses show that the serrated flow in HEAs consists of complex dynamical behavior. It is anticipated that this review will provide some useful and clarifying information regarding the serrated-flow mechanisms in this material system. Finally, suggestions for future research directions in this field are proposed, such as the effects of irradiation, additives (such as C and Al), the presence of nanoparticles, and twinning on the serrated flow behavior in HEAs.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Transforming U.S. Particle Physics Education: A Snowmass 2021 Study

The pursuit of knowledge in particle physics requires constant learning. As new tools become available, new theories are developed, and physicists search for new answers with ever-evolving methods. However, it is the case that formal educational systems serve as the primary training grounds for particle physicists. Graduate school (and undergraduate school to a lesser extent) is where researchers learn most of the technical skills required for research, develop scientific problem-solving abilities, learn how to establish themselves in their field, and begin developing their career. It is unfortunate, then, that the skills gained by physicists during their formal education are often mismatched with the skills actually required for a successful career in physics. We performed a survey of the U.S. particle physics community to determine the missing elements of graduate and undergraduate education and to gauge how to bridge these gaps. In this contributed paper, part of the 2021-22 Snowmass Community Planning Exercise, we report the results of this survey. We also recommend several specific community actions to improve the quality of particle physics education; the "community" here refers to physics departments, national labs, professional societies, funding agencies, and individual physicists.

Bitter, O.↗

DOE OSTI2022

Error analysis of numerical methods for thick diffusive neutron transport problems on Shishkin mesh

A thin layer will develop at the boundary if the incoming angular flux is anisotropic in thick diffusive neutron transport problems. Solving such singularly perturbed problems, which have non-smooth solutions with singularity near the boundary, is computationally challenging. Standard finite difference schemes on a uniform mesh cannot yield ε-uniform convergence, where ε is a small parameter, while it can be achieved on a suitable piecewise-uniform Shishkin mesh. We present a formal error analysis of the diamond difference (DD) method and step difference (SD) method for solving the S{sub N} neutron transport equation. The analysis can be extended to other finite difference methods. Numerical results are presented to confirm the error estimates and the advantages of the Shishkin mesh. (author)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2017

SS-HORSE method for studying resonances

A new method for analyzing resonance states based on the Harmonic-Oscillator Representation of Scattering Equations (HORSE) formalism and analytic properties of partial-wave scattering amplitudes is proposed. The method is tested by applying it to the model problem of neutral-particle scattering and can be used to study resonance states on the basis of microscopic calculations performed within various versions of the shell model.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2016

Clinical Significance of Accounting for Tissue Heterogeneity in Permanent Breast Seed Implant Brachytherapy Planning

Purpose: The inhomogeneity correction factor (ICF) method provides heterogeneity correction for the fast calculation TG43 formalism in seed brachytherapy. This study compared ICF-corrected plans to their standard TG43 counterparts, looking at their capacity to assess inadequate coverage and/or risk of any skin toxicities for patients who received permanent breast seed implant (PBSI). Methods and Materials: Two-month postimplant computed tomography scans and plans of 140 PBSI patients were used to calculate dose distributions by using the TG43 and the ICF methods. Multiple dose-volume histogram (DVH) parameters of clinical target volume (CTV) and skin were extracted and compared for both ICF and TG43 dose distributions. Short-term (desquamation and erythema) and long-term (telangiectasia) skin toxicity data were available on 125 and 110 of the patients, respectively, at the time of the study. The predictive value of each DVH parameter of skin was evaluated using the area under the receiver operating characteristic (ROC) curve for each toxicity endpoint. Results: Dose-volume histogram parameters of CTV, calculated using the ICF method, showed an overall decrease compared to TG43, whereas those of skin showed an increase, confirming previously reported findings of the impact of heterogeneity with low-energy sources. The ICF methodology enabled us to distinguish patients for whom the CTV V{sub 100} and V{sub 90} are up to 19% lower compared to TG43, which could present a risk of recurrence not detected when heterogeneity are not accounted for. The ICF method also led to an increase in the prediction of desquamation, erythema, and telangiectasia for 91% of skin DVH parameters studied. Conclusions: The ICF methodology has the advantage of distinguishing any inadequate dose coverage of CTV due to breast heterogeneity, which can be missed by TG43. Use of ICF correction also led to an increase in prediction accuracy of skin toxicities in most cases.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Momentum distributions for H 2 ( e , e ' p )

[Background] A primary goal of deuteron electrodisintegration is the possibility of extracting the deuteron momentum distribution. This extraction is inherently fraught with difficulty, as the momentum distribution is not an observable and the extraction relies on theoretical models dependent on other models as input. [Purpose] We present a new method for extracting the momentum distribution which takes into account a wide variety of model inputs thus providing a theoretical uncertainty due to the various model constituents. [Method] The calculations presented here are using a Bethe-Salpeter like formalism with a wide variety of bound state wave functions, form factors, and final state interactions. We present a method to extract the momentum distributions from experimental cross sections, which takes into account the theoretical uncertainty from the various model constituents entering the calculation. [Results] In order to test the extraction pseudo-data was generated, and the extracted "experimental'' distribution, which has theoretical uncertainty from the various model inputs, was compared with the theoretical distribution used to generate the pseudo-data. [Conclusions] In the examples we compared the original distribution was typically within the error band of the extracted distribution. The input wave functions do contain some outliers which are discussed in the text, but at least this process can provide an upper bound on the deuteron momentum distribution. Due to the reliance on the theoretical calculation to obtain this quantity any extraction method should account for the theoretical error inherent in these calculations due to model inputs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2014

Conductance matrix of multiterminal semiconductor devices with edge channels

A method for determining the conductance matrix of multiterminal semiconductor structures with edge channels is proposed. The method is based on the solution of a system of linear algebraic equations based on Kirchhoff equations, made up of potential differences U{sub ij} measured at stabilized currents I{sub kl}, where i, j, k, l are terminal numbers. The matrix obtained by solving the system of equations completely describes the structure under study, reflecting its configuration and homogeneity. This method can find wide application when using the known Landauer-Buttiker formalism to analyze carrier transport in the quantum Hall effect and quantum spin Hall effect modes. Within the proposed method, the contribution of the contact area resistances R{sub c} to the formation of conductance matrix elements is taken into account. The possibilities of practical application of the results obtained in developing analog cryptographic devices are considered.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Comparison of integral equations used to study ${T}_{cc}^{+}$ for a stable D *

We perform a detailed comparison between three formalisms used in recent studies of DD* scattering at heavier-than-physical pion masses, which aim to understand the properties of the doubly-charmed tetraquark, ${T}_{cc}^{+}$ (3875). These methods are the three-particle relativistic field theory (RFT) formalism, the two-body Lippmann-Schwinger (LS) equation with chiral effective field theory potentials, and the two-particle relativistic framework proposed by Baião Raposo and Hansen (BRH approach). In a simplified single-channel setting, we derive the conditions under which the infinite-volume integral equations from the RFT and BRH approaches reduce to the LS form. We present numerical examples showing that differences between these methods can be largely removed by adjusting short-range couplings. We also address a number of technical issues in the RFT approach.

Hadronic Spectroscopy↗