SEARCH · Engineering Papers
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.
We use the Auxiliary Field Diffusion Monte Carlo (AFDMC) method to calculate ground state properties of nuclei up to A=16. We employ local chiral interactions up to next-to-next-to-leading order (N 2 LO), including consistent three-nucleon forces, fit to few-body observables.
An experimental and simulation study of warm dense carbon foams at ambient density (n e ~10 21 cm –3 ) is presented. This study of isochorically heated foams is motivated by their potential application in carbon-atmosphere white-dwarf envelopes, where there are modeling uncertainties due to the equation of state. The foams are heated on an approximately picosecond time scale with a laser-accelerated proton beam. The cooling and expansion of the heated foams can be modeled with appropriately initialized radiation-hydrodynamics codes; xRAGE code is used in this work. The primary experimental diagnostic is the streaked optical pyrometer, which images a narrow band of radiation from the rear surface of the heated material. Presented are xRAGE modeling results for both solid aluminum targets and carbonized resorcinol-formaldehyde foam targets, showing that the foam appears to cool slowly on the pyrometer because of partial transparency. So that simulations of cooling foam are processed properly, it is necessary to account for finite optical depth in the photosphere calculation, and the methods for performing that calculation are presented in depth.
The present work reports on the development of techniques for in-field fast-neutron radiography measurements using the associated-particle imaging (API) method. The API method employs alpha-neutron coincidences from the d+t→α+n reaction to enable fast neutron transmission imaging with excellent contrast using a wide cone beam. However, for field radiography applications, the API method is burdened by the need for the relative positions of the source and detector to be known. Fortunately, these relative positions can be inferred from transmission data. The inferred positions also enable accurate stitching of multiple images into a composite image even when using a low-resolution detector panel and acquiring images having few overlapping pixels. The developed techniques address analysis of measurements where (1) the source and detector panel are separately hand positioned rather than held in registration by a gantry, (2) multiple detector panel positions within the “coincident cone” of tagged neutrons are required to piece together an image of an item of interest, and (3) normalization measurements that have identical source-detector positioning but without the inspected object are not possible. The present work will describe the system calibrations (including timing calibrations and neutron direction calibrations) necessary for subsequent analyses, the method of locating the detector in the coincident cone of neutrons with millimeter precision using the timing and directions of coincident neutrons, the method of calculating a normalization image for a given detector panel position, and the method used to project multiple images into a common image.
Abstract The electronic transport behaviour of materials determines their suitability for technological applications. We develop a computationally efficient method for calculating carrier scattering rates of solid-state semiconductors and insulators from first principles inputs. The present method extends existing polar and non-polar electron-phonon coupling, ionized impurity, and piezoelectric scattering mechanisms formulated for isotropic band structures to support highly anisotropic materials. We test the formalism by calculating the electronic transport properties of 23 semiconductors, including the large 48 atom CH 3 NH 3 PbI 3 hybrid perovskite, and comparing the results against experimental measurements and more detailed scattering simulations. The Spearman rank coefficient of mobility against experiment ( r s = 0.93) improves significantly on results obtained using a constant relaxation time approximation ( r s = 0.52). We find our approach offers similar accuracy to state-of-the art methods at approximately 1/500th the computational cost, thus enabling its use in high-throughput computational workflows for the accurate screening of carrier mobilities, lifetimes, and thermoelectric power.
Theoretical ab initio studies of resonances in the unbound 7 He nucleus are presented. We perform no-core shell-model calculations with NN interactions Daejeon16 and JISP16 and utilize the SS-HORSE method to calculate the S matrix for two-body channels n + 6 He and n + 6 He* with 6He respectively in the ground and excited 2 + states as well as for the four-body democratic decay channel 4 He + n + n + n. The resonant energies and widths are obtained by numerical location of the S-matrix poles. Here, we describe all experimentally known 7 He resonances and suggest an interpretation of an observed wide resonance of unknown spin-parity.
The resistance of the electric machine’s windings is an important parameter when analyzing the machine’s performance, as conduction losses in the winding are often the largest loss component. However, this resistance changes with frequency due to the skin and proximity effects, making it challenging to accurately estimate these losses. In this article, a computationally efficient, general method of calculating the ac resistance of an electric machine winding with rectangular conductors is presented. The method uses a 1-D magnetoquasi-static model of the stator slot region to develop boundary transfer relations between individual turns. These relations are then used to create an impedance matrix from which the ac resistance can be calculated. This technique is more flexible than previous solutions since it can accommodate different winding configurations such as fractional-pitch windings. It is shown that a fractional-pitch winding has a lower ac resistance than a full-pitch winding. Finally, finite-element analysis simulations and experimental results are used to validate the proposed method.
We provide a method to calculate the rate of false vacuum decay induced by a black hole. The method uses complex tunneling solutions and consistently takes into account the structure of different quantum vacua in the black hole metric via boundary conditions. The latter are connected to the asymptotic behavior of the time-ordered Green’s function in the corresponding vacua. We illustrate the technique on a two-dimensional toy model of a scalar field with inverted Liouville potential in an external background of a dilaton black hole. We analytically derive the exponential suppression of tunneling from the Boulware, Hartle-Hawking and Unruh vacua and show that they are parametrically different. The Unruh vacuum decay rate is exponentially smaller than the decay rate of the Hartle-Hawking state, though both rates become unsuppressed at high enough black hole temperature. We interpret the vanishing suppression of the Unruh vacuum decay at high temperature as an artifact of the two-dimensional model and discuss why this result can be modified in the realistic case of black holes in four dimensions.
A new general method for calculating the alpha decay half-life is presented. The method predicts an a priori exact value for the beryllium-8 half-life. Beryllium-8 is an exception to the current alpha decay theory captured in the Geiger-Nuttal law. The new method predicts the beryllium-8 alpha half-life using only constants and measured isotopic mass. The method also reliably predicts all the heavier isotope alpha decay half-lives consistent with the Geiger-Nuttal law. With respect to current theory, the inability of the Geiger-Nuttal Law to predict the alpha-decay half-life in the case of beryllium-8 has led to consideration of other decay mechanisms for this isotope, such as fission for example. One result is that given the consistency of the new method presented here for all isotopes including an exact a priori result for beryllium-8, the evidence strongly suggests that the beryllium-8 decay is in fact an alpha decay. A second result is that the method definitively demonstrates that the entire rest mass of the two helium-4 electrons is converted to energy in the decay process and this energy becomes part of the emitted alpha particle kinetic energy.
Abstract Recent wildfires in the western United States have led to substantial economic losses and social stresses. There is a great concern that the new climatic state may further increase the intensity, duration, and frequency of wildfires. To examine temporal and spatial features of historical wildfire trends and future changes, a common fire index, the Keetch–Byram Drought Index (KBDI), is calculated over the contiguous United States (CONUS) and Alaska. We introduce an efficient initialization method and calculate 36 years (1982–2017) of historical KBDI at 4 km using a high-quality observational dataset. KBDIs are also calculated at 12 km using regional climate models and extended into the mid- and late-21st century. Based on the observational data, annual mean (and 95th percentile) KBDI over forested regions in the southwestern and northwestern CONUS has risen since 1982 at a rate of 5.2 (4.0) and 2.9 (3.3) per year, respectively, indicating a persistent drying trend fostering fire activities; the number of days exceeding the top 5% historical KBDI has increased by 16 and 25 d in the 36 years. Multiple regional climate model simulations project increasing wildfire potential and longer fire seasons over broader areas based on the estimated KBDI for the mid- and late-21st century. By the end of the century, most of the CONUS would experience about 90–189 more days per year exceeding the historical local maximum KBDI; areas with high KBDI (>600), whose occurrence correlates with large burned area, are projected to broaden by nearly 60 times in the southern CONUS. While both temperature and precipitation contribute to future KBDI changes, warming is the main driver of more frequent, intense and wide-spread extreme wildfires indicated by high KBDIs in future projections.
This article is part-I of a review of density-functional theory (DFT) that is the most widely used method for calculating electronic structure of materials. The accuracy and ease of numerical implementation of DFT methods has resulted in its extensive use for materials design and discovery and has thus ushered in the new field of computational material science. In this article, we start with an introduction to Schrödinger equation and methods of its solutions. After presenting exact results for some well-known systems, difficulties encountered in solving the equation for interacting electrons are described. How these difficulties are handled using the variational principle for the energy to obtain approximate solutions of the Schrödinger equation is discussed. The resulting Hartree and Hartree–Fock theories are presented along with results they give for atomic and solid-state systems. We then describe Thomas–Fermi theory and its extensions which were the initial attempts to formulate many-electron problem in terms of electronic density of a system. Having described these theories, we introduce modern DFT by discussing Hohenberg–Kohn theorems that form its foundations. We then go on to discuss Kohn–Sham (KS) formulation of DFT in its exact form. Next, local density approximation (LDA) is introduced and solutions of KS equation for some representative systems, obtained using the LDA, are presented. We end part-I of the review describing the contents of part-II.
Monte Carlo simulations of low-energy ( <50 keV) electron transport in matter are essential for a broad range of application fields. Several Monte Carlo codes have developed specialized treatments for this case, but a comprehensive validation of low-energy electron transport for general-purpose simulations remains lacking in the literature. One approach to accomplish this validation is calculation of stopping power using low-energy electron transport physics, as stopping power is a fundamental radiation transport quantity which must be simulated accurately for nearly any application. Here in this work, we use the Monte Carlo N-Particle (MCNP) radiation transport code with the single-event method for electron transport to calculate stopping powers of low-energy electrons (50 eV to 30 keV) in 41 elemental solids, 14 compound solids, and five rare gas solids, comparing simulation results to published semi-empirical stopping power calculations from optical measurements. In general, the simulations give good agreement (typically within ±10%) with semi-empirical stopping power values at higher energies: 300 eV and above for most elemental solids, 1 keV and above for compound solids, and 400 eV and above for rare gas solids. Agreement between MCNP and semi-empirical values is worse below these energies. The most significant source of error is the EPRDATA14 cross section data, which does not account for changes in electronic structure due to solid-state bonding, particularly in compound materials. The simplistic model of atomic excitation used to generate the EPRDATA14 cross sections is another key source of error. Additionally, the breakdown of the continuous slowing-down approximation introduces significant uncertainty at low energies, although this is a limitation of the calculation method and not of the simulation procedure. Accounting for these and other uncertainty sources, the single event method in MCNP is robust and able to give good accuracy for a variety of low-energy electron transport problems through diverse kinds of materials.
ISO/TS 6336-22 (Calculation of load capacity of spur and helical gears - Part 22: Calculation of micropitting load capacity) is the ISO technical specification containing a proposal for a method to calculate the risk of micropitting in gear sets. It assesses micropitting risk through a safety factor which is calculated as the minimum specific film thickness in the contact zone divided by a permissible specific film thickness. A sliding parameter is used to adjust the film thickness in the contact zone to the region of highest sliding, which is where micropitting has been observed. The permissible specific film thickness is best determined through experience or testing, but there is an option to estimate it based on the lubricant’s failure load stage in FZG testing. This paper is a case study of how the method behaves when applied to several gear sets and applications that have experienced micropitting in the field.
Abstract Accurate predictions of inclusive scattering cross sections in the linear response regime require efficient and controllable methods to calculate the spectral density in a strongly-correlated many-body system. In this work we reformulate the recently proposed Gaussian Integral Transform technique in terms of Fourier moments of the system Hamiltonian which can be computed efficiently on a quantum computer. One of the main advantages of this framework is that it allows for an important reduction of the computational cost by exploiting previous knowledge about the energy moments of the spectral density. For a simple model of medium mass nucleus like $$^{40}$$ 40 Ca and target energy resolution of 1 MeV we find an expected speed-up of $$\approx 125$$ ≈ 125 times for the calculation of the giant dipole response and of $$\approx 50$$ ≈ 50 times for the simulation of quasi-elastic electron scattering at typical momentum transfers.
This article discusses a unique agglomeration modeling methodology developed based on binary collisions to combine the effects of heterogeneity in ash chemistry and granular physics. A simple population balance is defined to find changes in the particle size distribution (PSD) of a fluidized bed. Thermodynamic equilibrium calculations and a computational fluid dynamics (CFD) code, are used to obtain hydrodynamic parameters. A method to calculate and use a distribution of collision frequencies in fluidized beds, corresponding to a poly-disperse particle size distribution, was developed in order to incorporate the particle level heterogeneities. A distribution of collision frequencies obtained for poly-dispersed particles showed a three orders of magnitude higher collision frequency amongst the smaller particles than the coarser ones, at the initiation of agglomeration. Ash agglomeration occurred when the slag amount (binder) was less than 10 wt% at temperatures <850 °C for Pittsburgh No. 8 coal.
The past decade has seen the emergence of ab initio computational methods for calculating phonon-limited carrier mobilities in semiconductors with predictive accuracy. More realistic calculations ought to take into account additional scattering mechanisms such as, for example, impurity and grain-boundary scattering. Here, in this paper, we investigate the effect of ionized-impurity scattering on the carrier mobility. We model the analytical impurity potential parameterized from first principles by a collection of randomly distributed Coulomb scattering centers, and we include this relaxation channel into the ab initio Boltzmann transport equation, as implemented in the EPW code. We demonstrate this methodology by considering silicon, silicon carbide, and gallium phosphide, for which detailed experimental data are available. Our calculations agree well with experiments over a broad range of temperatures and impurity concentrations. For each compound investigated here, we compare the relative importance of electron-phonon scattering and ionized-impurity scattering, and we critically assess the reliability of Matthiessen's rule. We also show that an accurate description of dielectric screening and carrier effective masses can improve quantitative agreement with experiments.
We introduce the transition-density formalism, an efficient and general method for calculating the interaction of external probes with light nuclei. One- and two-body transition densities that encode the nuclear structure of the target are evaluated once and stored. They are then convoluted with an interaction kernel to produce amplitudes, and hence observables. Here, by choosing different kernels, the same densities can be used for any reaction in which a probe interacts perturbatively with the target. The method therefore exploits the factorisation between nuclear structure and interaction kernel that occurs in such processes. We study in detail the convergence in the number of partial waves for matrix elements relevant in elastic Compton scattering on 3 He. The results are fully consistent with our previous calculations in Chiral Effective Field Theory. But the new approach is markedly more computationally efficient, which facilitates the inclusion of more partial-wave channels in the calculation. We also discuss the usefulness of the transition-density method for other nuclei and reactions. Calculations of elastic Compton scattering on heavier targets like 4 He are straightforward extensions of this study, since the same interaction kernels are used. And the generality of the formalism means that our 3 He densities can be used to evaluate any 3 He elastic-scattering observable with contributions from one- and two-body operators.
ABSTRACT Magnesium mono-sulphide (MgS) plays a crucial role in astrochemical processes, particularly in the interstellar medium (ISM), where metal-sulphur chemistry influences molecular formation and evolution. This work presents a detailed study of the rotational excitation and de-excitation of MgS through collisions with helium (He) atoms, the second most abundant species in the ISM after hydrogen, which includes both atomic (H) and molecular forms (H2). The focus on MgS–He collisions arises from He's high abundance, chemical inertness, and simpler electronic structure, which make it well suited for quantum scattering calculations. These characteristics establish He as an ideal candidate for initial studies, providing fundamental data for future investigations involving H2. The study uses quantum scattering methods to calculate the collisional rate coefficients over a broad temperature range. These rates are critical for interpreting observational data on MgS and predicting its abundance in space. The interaction potential between MgS and He is calculated using the rigid rotor approximation and the Jacobi coordinate system, employing the CCSD(T)-F12a/aug-cc-pVTZ method for accurate two-dimensional potential energy surface. The study explores the anisotropic nature of the MgS–He interaction, which favours odd Δj rotational transitions at low collision energies. The inelastic cross-sections for rotational transitions involving up to 16 rotational levels of MgS were computed up to 1000 cm−1, enabling the calculation of rate coefficients up to 150 K for Δj = ±1, ±2, and ± 3 rotational transitions. The results show that Δj = 1 transitions dominate at low temperatures, while Δj = 2 transitions become more significant at higher temperatures. This study provides valuable data for interpreting future astrophysical observations of MgS. The findings also propose new rotational transitions for MgS detection in space, enhancing our ability to track and study this molecule in various cosmic environments.
An efficient implementation for the relativistic exact two-component core–valence-separated equation-of-motion coupled-cluster singles and doubles (X2C-CVS-EOM-CCSD) method is reported. The explicit exclusion of pure valence excitations in the EOM-CCSD excited-state eigenvalue equations significantly improves the efficiency for calculations of core-excited states. Benchmark relativistic CVS-EOM-CC calculations with systematic inclusion of relativistic, correlation, and basis-set effects are shown to provide highly accurate results for core ionized and excited states involving heavy atoms.