Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALCULATION 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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Variational Self-Consistent Theory for Beam-Loaded Cavities

A new variational theory is presented for beam loading in microwave cavities. The beam{field interaction is formulated as a dynamical interaction whose stationarity according to Hamilton’s principle will naturally lead to steady-state solutions that indicate how a cavity’s resonant frequency, Q and optimal coupling coefficient will detune as a result of the beam loading. A driven cavity Lagrangian is derived from first principles, including the effects of cavity wall losses, input power and beam interaction. The general formulation is applied to a typical klystron input cavity to predict the appropriate detuning parameters required to maximize the gain (or modulation depth) in the average Lorentz factor boost, (Δ γ). Numerical examples are presented, showing agreement with the general detuning trends previously observed in the literature. The developed formulation carries several advantages for beam-loaded cavity structures. It provides a self-consistent model for the dynamical (nonlinear) beam{field interaction, a procedure for maximizing gain under beam-loading conditions, and a useful set of parameters to guide cavity-shape optimization during the design of beam-loaded systems. Enhanced clarity of the physical picture underlying the problem seems to be gained using this approach, allowing straightforward inclusion or exclusion of different field configurations in the calculation and expressing the final results in terms of measurable quantities. Two field configurations are discussed for the klystron input cavity, using finite magnetic confinement or no confinement at all. Formulating the problem in a language that is directly accessible to the powerful techniques found in Hamiltonian dynamics and canonical transformations may potentially carry an additional advantage in terms of analytical computational gains, under suitable conditions

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Electronic band gaps from quantum Monte Carlo methods

We develop a method for calculating the fundamental electronic gap of semiconductors and insulators using grand canonical quantum Monte Carlo simulations. Further, we discuss the origin of the bias introduced by supercell calculations of finite size and show how to correct the leading and subleading finite size errors either based on observables accessible in the finite-sized simulations or from density-functional theory calculations. Our procedure is applied to carbon, silicon, and molecular hydrogen crystals, and compared to experiment for carbon and silicon.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Landau level of fragile topology

Here, we study the Hofstadter butterfly and Landau levels of the twisted bilayer graphene (TBG). We show that the nontrivial fragile topology of the lowest two bands near the charge neutral point makes their Hofstadter butterfly generically connected with higher bands, closing the gap between the first and second conduction (valence) bands at a certain magnetic flux per unit cell. We also develop a momentum space method for calculating the TBG Hofstadter butterfly, from which we identify three phases where the Hofstadter butterflies of the lowest two bands and the higher bands are connected in different ways. We show this leads to a crossing between the $\textit{ν}$ = 4 Landau fan from the charge neutral point and the zero-field band gap at one flux per moiré unit cell, which corresponds to a magnetic field 25$θ^2$ T (twist angle $\textit{θ}$ in degrees). This provides an experimentally testable feature of the fragile topology. In general, we expect it to be a generic feature that the Hofstadter butterfly of the topological bands are connected with the Hofstadter spectra of other bands. We further show the TBG band theory with Zeeman splitting being the most sizable splitting could result in Landau fans at the charge neutral point and half fillings near the magic angle, and we predict their variations under an in-plane magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Nonlocal and nonadiabatic Pauli potential for time-dependent orbital-free density functional theory

Time-dependent orbital-free density functional theory is an efficient ab initio method for calculating the electronic dynamics of large systems. In comparison to standard time-dependent density functional theory, it computes only a single electronic state regardless of system size, but it requires an additional time-dependent Pauli potential term. Herein we propose a nonadiabatic and nonlocal Pauli potential whose main ingredients are the time-dependent particle and current densities. Our calculations of the optical spectra of metallic and semiconductor clusters indicate that nonlocal and nonadiabatic time-dependent orbital-free density functional theory performs accurately for metallic systems and semiquantitatively for semiconductors. This paper opens the door to wide applicability of time-dependent orbital-free density functional theory for nonequilibrium electron and electron-nuclear dynamics of complex materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Efficient time-dependent orbital-free density functional theory: Semilocal adiabatic response

Orbital-free density functional theory and its time-dependent extension are efficient ab initio methods for calculating the electronic structure and dynamics of large systems. Through the calculation of the optical spectra of selected clusters, we reach three important conclusions: (1) The quality of the spectra is strongly affected by the quality of the corresponding ground-state electron density; (2) the adiabatic part of the electronic response to external perturbations can be safely evaluated at the semilocal level; and (3) the nonadiabatic, current-dependent part of the time-dependent Pauli potential is key to recover correct spectral envelopes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hydrogen and deuterium tunneling in niobium

We use density functional methods to identify the atomic configurations of H and D atoms trapped by O impurities embedded in bulk Nb. The O atoms are located at the octahedral position in the Nb body-centered cubic (BCC) lattice, and H (D) atoms tunnel between two degenerate tetrahedral sites separated by a mirror plane. Using nudged elastic band (NEB) methods, we calculate the double-well potential for O-H and O-D impurities and the wave functions and tunnel splittings for H and D atoms. Our results agree with those obtained from analysis of heat capacity and neutron scattering measurements on Nb with low concentrations of O-H and O-D.

density functional theory↗
DOE OSTI2020

K-shell internal conversion coefficient for M4 decay of the 30.8 keV isomer in 93 Nb

The 30.8 keV (1/2) – isomeric state in $^{93}_{41}$Nb decays by an M4 transition to the (9/2) + ground state. Here, we have measured the K-shell internal conversion coefficient α K for this transition in order to test the validity of current methods of calculation for a case with relatively low atomic number, low transition energy, and a high value of α K . Taking the fluorescence yield of niobium to be 0.751(4), we obtain α K = 2.56(9) × 10 4 , a result that agrees with the Dirac-Fock calculations that take into account the presence of a K vacancy in the final state, and disagrees by 1.7 standard deviations with calculations that ignore the vacancy. We also determine the energy of the isomeric state to be 30.760(5) keV, a fourfold improvement in precision over the currently accepted value.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Improvements to the macroscopic-microscopic approach of nuclear fission

The well-established macroscopic-microscopic (mac-mic) description of nuclear fission enables the prediction of fission-fragment yields for a broad range of fissioning systems. In this work, we present several key enhancements to this approach. We improve upon the microscopic sector of nuclear potential-energy surfaces by magnifying the resolution of the Lipkin-Nogami equations and strengthening the Strutinsky procedure, thus reducing spurious effects from the continuum. We further present a novel deterministic method for calculating fission dynamics under the assumption of strongly damped nucleonic motion. Our technique directly determines the evolution of the scissioned shape distribution according to the number of random-walk steps rather than the statistical accumulation of fission events. We show that our new technique is equivalent to the Metropolis random walk pioneered over the past decade by Randrup and colleagues. It further improves upon it because we remove the need for altering the nuclear landscape via a biased potential. With our final improvement, we calculate fission fragments mass and charge distributions using particle number projection, which affords the simultaneous calculation of both mass- and charge-yield distributions. Fission fragments are thus calculated from the quantum-mechanical A -body states of the potential-energy surface rather than from the collective mass asymmetry variable α g of the finite-range liquid-drop model used in past work. We highlight the success of our enhancements by predicting the odd-even staggering and the charge polarization for the neutron-induced fission of 233 U and 235 U .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Quantum computation of nuclear observables involving linear combinations of unitary operators

Here we present the quantum computation of nuclear observables where the operators of interest are first decomposed in terms of the linear combination of unitaries. Then we utilize the Hadamard test and the linear combination of unitaries based methods to compute the expectation values. We apply these methods to calculate the electric quadrupole moment of the deuteron. The results are compared for the Jordan-Wigner transformation and Gray code encoding. We discuss the versatility of our approach that can be utilized in general to calculate several observables on a quantum computer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Scale setting the Möbius domain wall fermion on gradient-flowed HISQ action using the omega baryon mass and the gradient-flow scales 𝑡 0 and 𝑤 0

We report on a subpercent scale determination using the omega baryon mass and gradient-flow methods. The calculations are performed on 22 ensembles of 𝑁 𝑓 =2 +1 +1 highly improved, rooted staggered sea-quark configurations generated by the MILC and CalLat Collaborations. The valence quark action used is Möbius domain wall fermions solved on these configurations after a gradient-flow smearing is applied with a flowtime of 𝑡 gf = 1 in lattice units. The ensembles span four lattice spacings in the range 0.06 ≲ 𝑎 ≲0.15 fm, six pion masses in the range 130 ≲ 𝑚 𝜋 ≲ 400 MeV and multiple lattice volumes. On each ensemble, the gradient-flow scales 𝑡 0 /𝑎 2 and 𝑤 0 /𝑎 and the omega baryon mass 𝑎⁢𝑚 Ω are computed. The dimensionless product of these quantities is then extrapolated to the continuum and infinite volume limits and interpolated to the physical light, strange and charm quark mass point in the isospin limit, resulting in the determination of $\sqrt{t}_0$ = 0.1422⁢(14) fm and 𝑤 0 = 0.1709⁢(11) fm with all sources of statistical and systematic uncertainty accounted for. The dominant uncertainty in both results is the stochastic uncertainty, though for $\sqrt{t}_0$ there are comparable continuum extrapolation uncertainties. For 𝑤 0 , there is a clear path for a few-per-mille uncertainty just through improved stochastic precision, as recently obtained by the Budapest-Marseille-Wuppertal Collaboration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Dark matter isocurvature from curvature

Isocurvature fluctuations, where the relative number density of particle species spatially varies, can be generated from initially adiabatic, or curvature, fluctuations if the various species fall out of or were never in thermal equilibrium. The freezing of the thermal relic dark matter abundance is one such case, but for modes that are still outside the horizon the amplitude is highly suppressed and originates from the small change in the local expansion rate due to the local space curvature produced by the curvature fluctuation. We establish a simple separate-universe method for calculating this generation that applies to both freeze-in and freeze-out models, identify three critical epochs for this process, and give general scaling behaviors for the amplitude in each case: the freezing epoch, the kinetic decoupling epoch and matter-radiation equality. Freeze-out models are typically dominated by spatially modulated annihilation from the latter epochs and can generate much larger isocurvature fluctuations compared with typical freeze-in models, albeit still very small and observationally allowed by cosmic microwave background measurements. Here, we illustrate these results with concrete models where the dark matter interactions are vector or scalar mediated.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Absorption of charged particles in perfectly matched layers by optimal damping of the deposited current

Perfectly matched layers (PMLs) are widely used in particle-in-cell simulations, in order to absorb electromagnetic waves that propagate out of the simulation domain. However, when charged particles cross the interface between the simulation domain and the PMLs, a number of numerical artifacts can arise. Here, in order to mitigate these artifacts, we introduce a PML algorithm whereby the current deposited by the macroparticles in the PML is damped by an analytically derived optimal coefficient. The benefits of this algorithm are illustrated in practical simulations. In particular, it is shown that this algorithm is well suited for particles exiting the box in near-normal incidence, in the sense that the fields behave as if the exiting particle is propagating in an infinite vacuum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Resolving local structural motifs across the phase evolution of zinc titanates with computational x-ray absorption spectroscopy

Resolving the local structure motifs that characterize phase evolution as a function of composition is a key challenge in structure characterization of complex materials. Here, in this study, we combine first-principles simulations and x-ray absorption near-edge structures (XANES) analysis to gain insights into the structure evolution revealed by measurements across a combinatorial zinc titanate thin film, which was grown with smoothly varying composition over a wide range of the Ti:Zn ratio. Specifically, we propose a cluster blind-signal-separation (cBSS) method for XANES spectral analysis based on a library of the structures and spectra of representative local motifs. In addition to motifs from zinc titanate crystals, two types of Ti-defect models constructed in this study are key to the understanding of the structure characteristics in the Zn-rich region. The cBSS method makes use of both spectral clustering of the simulated site-XANES spectra library and the BSS procedure to construct high-fidelity spectral basis functions from an experimental spectral sequence. The method provides a rigorous measure of the spectral sensitivity and basis completeness. The results of the XANES analysis are corroborated with other experimental modalities, including x-ray diffraction and spectroscopic ellipsometry, to validate the cBSS method. The calculated motif weights resulting from fitting the XANES spectra with the cBSS basis probe the atomic structure characteristics of both crystalline and amorphous phases as a function of the Ti/Zn composition. The insights of the local structure motif evolution are pivotal to the understanding of the nonmonotonic trend in the optical gap, which may lead to potential applications through tuning the optical properties of zinc titanate. The workflow of the XANES spectral analysis developed in this work can be generalized to construct the structure-property relationship in a broad material space.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Matrix-Model Simulations Using Quantum Computing, Deep Learning, and Lattice Monte Carlo

Matrix quantum mechanics plays various important roles in theoretical physics, such as a holographic description of quantum black holes, and it underpins the only practical numerical approach to the study of complex high-dimensional supergravity theories. Understanding quantum black holes and the role of entanglement in a holographic setup is of paramount importance for the realization of a quantum theory of gravity. Moreover, a complete numerical understanding of the holographic duality and the emergence of geometric space-time features from microscopic degrees of freedom could pave the way for new discoveries in quantum information science. Euclidean lattice Monte Carlo simulations are the de facto numerical tool for understanding the spectrum of large matrix models and have been used to test the holographic duality. However, they are not tailored to extract dynamical properties or even the quantum wave function of the ground state of matrix models. Quantum computing and deep learning provide potentially useful approaches to study the dynamics of matrix quantum mechanics. If successful in the context of matrix models, these rapidly improving numerical techniques could become the new Swiss army knife of quantum gravity practitioners. In this paper, we perform the first systematic survey for quantum computing and deep-learning approaches to matrix quantum mechanics, comparing them to lattice Monte Carlo simulations. These provide baseline benchmarks before addressing more complicated problems. In particular, we test the performance of each method by calculating the low-energy spectrum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Protonation of histidine rings using quantum-mechanical methods

Histidine can be protonated on either or both of the two N atoms of the imidazole moiety. Each of the three possible forms occurs as a result of the stereochemical environment of the histidine side chain. In an atomic model, comparing the possible protonation states in situ, looking at possible hydrogen bonding and metal coordination, it is possible to predict which is most likely to be correct. A more direct method is described that uses quantum-mechanical methods to calculate, also in situ, the minimum geometry and energy for comparison, and therefore to more accurately identify the most likely proton­ation state.

Histidine↗
DOE OSTI2020

Solar/Photovoltaic DC Systems: Information for Electrical Workers and Firefighters

Solar photovoltaic (PV) systems are common and continue to gain popularity. This article helps electrical workers and firefighters understand basic operating principles and hazards of PV dc arrays. The effects of solar radiation and temperature on power output are discussed including; PV-module testing standards; methods for calculating the number of modules in a PV string; some National Electrical Code (NEC) guidelines that are aimed toward improving array safety; dc-array commissioning, periodic maintenance, and testing; arc-flash-hazard potential ways firefighters can work safely around PV arrays; potential risks posed by moonlight and other light sources; some unique dc-system hazards. We also present statistics on PV dc-array electrical incidents and injuries.

14 SOLAR ENERGY↗
DOE OSTI2024

Geomagnetically Induced Current Field Test on Large Grid-Connected Power Transformers: Analysis, Model Development, and Simulations

Geomagnetic-induced current (GIC) flow in power grids can cause undesirable effects such as transformer overheating, harmonics, higher reactive power demand, etc. Many simulation models have been developed to study these effects, but real-world verification on modern transformer designs is rare. Here, this paper presents the first long-duration GIC field test in the U.S. performed on high-voltage, grid-connected transformers featuring winding clamps and tie rods instead of conventional tie bars. Field measurements were taken to evaluate GIC effects. These measurements also aided in developing and validating thermal and electromagnetic transient (EMT) models of the transformers. During the test, significant current and voltage distortions were observed along with considerable transformer reactive power losses. Analysis of the field measurements showed that the transformers’ hottest spot was at the inner windings, and their k-factors were close to factory test and software default values. Thermal simulations indicated that the transformers would not violate their thermal limits even for a GIC waveform that peaks at about 200 A/phase. EMT simulations revealed that increased transformer loading may reduce GIC-induced reactive power demand and harmonics in certain scenarios. The study also highlighted potential inaccuracies in using the k-factor method to calculate transformer reactive power losses.

EMTDC↗