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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Power Density Distributions in Electron Beams

Measuring and controlling the power density distribution of electron beams used for welding is critical for producing repeatable welds and for transferring welding parameters between different machines. On any given machine, the power density distribution is controlled by defocusing the beam relative to its sharpest focused condition. However, measuring the power density distribution can be difficult due to the intense nature of welding beams and is further complicated by imperfect electron optics that can distort the beam, making it difficult to quantify. The enhanced modified Faraday cup (EMFC) diagnostic method was used here for beam analysis that utilizes computed tomography to reconstruct the beam’s power density distribution. These results were compared to the International Standards Organization (ISO) method for characterizing laser beams using a second-moment D4σ calculation. For ideal Gaussian-shaped beams, both methods would give the same result. However, for imperfect beams, the calculated D4σ diameter was shown to be about 25% larger relative to the FWe2 diameter measured by the EMFC due to the heavier weighting of data in the tails of the beam by D4σ. Finally, although both methods produce repeatable welds, it is important to understand the differences in the reported beam diameters, divergence angles, and beam parameter products when transferring parameters between machines.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Accurate and efficient description of interacting carriers in quantum nanostructures by selected configuration interaction and perturbation theory

In this paper, we present a method to accurately and efficiently calculate many-body states of interacting carriers in quantum nanostructures based on a combination of iterative selection of configurations and perturbation theory. This method enables investigations of large excitonic complexes and multielectron systems with near full configuration interaction accuracy, even though only a small subspace of the full many-body Hilbert space is sampled, thus saving orders of magnitudes in computational resources. Important advantages of this method are that the convergence is controlled by a single parameter, the threshold, and that ground and excited states can be treated on an equal footing. On the example of InAsP nanowire quantum dots described using a million-atom tight-binding approach, we demonstrate the extreme efficiency of the method by numerical studies of large excitonic complexes filling up to the fourth electronic shell. We find that the method generally converges fast as a function of the threshold, profiting from a significant acceleration due to the perturbative corrections. The role of the choice of single-particle basis states is discussed. It is found that the algorithm converges faster in the Hartree-Fock basis only for highly charged systems, where Coulomb repulsion dominates. Finally, based on the observation that second-order perturbative energy corrections only depend on off-diagonal elements of the many-body Hamiltonian, we present a way to accurately calculate many-body states that requires only a relatively small number of Coulomb matrix elements.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Data from: Comparing traditional and Bayesian approaches to ecological meta-analysis

Despite the wide application of meta-analysis in ecology, some of the traditional methods used for meta-analysis may not perform well given the type of data characteristic of ecological meta-analyses. We reviewed published meta-analyses on the ecological impacts of global climate change, evaluating the number of replicates used in the primary studies (ni) and the number of studies or records (k) that were aggregated to calculate a mean effect size. We used the results of the review in a simulation experiment to assess the performance of conventional frequentist and Bayesian meta-analysis methods for estimating a mean effect size and its uncertainty interval. Our literature review showed that ni and k were highly variable, distributions were right-skewed, and were generally small (median ni =5, median k=44). Our simulations show that the choice of method for calculating uncertainty intervals was critical for obtaining appropriate coverage (close to the nominal value of 0.95). When k was low (<40), 95% coverage was achieved by a confidence interval based on the t-distribution that uses an adjusted standard error (the Hartung-Knapp-Sidik-Jonkman, HKSJ), or by a Bayesian credible interval, whereas bootstrap or z-distribution confidence intervals had lower coverage. Despite the importance of the method to calculate the uncertainty interval, 39% of the meta-analyses reviewed did not report the method used, and of the 61% that did, 94% used a potentially problematic method, which may be a consequence of software defaults. In general, for a simple random-effects meta-analysis, the performance of the best frequentist and Bayesian methods were similar for the same combinations of factors (k and mean replication), though the Bayesian approaches had higher than nominal (>95%) coverage for the mean effect when k was very low (k<15). Our literature review suggests that many meta-analyses that used z-distribution or bootstrapping confidence intervals may have over-estimated the statistical significance of their results when the number of studies was low; more appropriate methods need to be adopted in ecological meta-analyses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Methodology for Sensitivity Analysis of Homogenized Cross-Sections to Instantaneous and Historical Lattice Conditions with Application to AP1000® PWR Lattice

In the two-step method for nuclear reactor simulation, lattice physics calculations are performed to compute homogenized cross-sections for a variety of burnups and lattice configurations. A nodal code is then used to perform full-core analysis using the pre-calculated homogenized cross-sections. One source of uncertainty introduced in this method is that the lattice configuration or depletion conditions typically do not match a pre-calculated one from the lattice physics simulations. Therefore, some interpolation model must be used to estimate the homogenized cross-sections in the nodal code. This current study provides a methodology for sensitivity analysis to quantify the impact of state variables on the homogenized cross-sections. This methodology also allows for analyses of the historical effect that the state variables have on homogenized cross-sections. An application of this methodology on a lattice for the Westinghouse AP1000® reactor is presented where coolant density, fuel temperature, soluble boron concentration, and control rod insertion are the state variables of interest. The effects of considering the instantaneous values of the state variables, historical values of the state variables, and burnup-averaged values of the state variables are analyzed. Using these methods, it was found that a linear model that only considers the instantaneous and burnup-averaged values of state variables can fail to capture some variations in the homogenized cross-sections.

Price, Dean (ORCID:0000000309990111)↗
DOE OSTI2024

Quarkonium transport in weakly and strongly coupled plasmas

We report on progress in the nonperturbative understanding of quarkonium dynamics inside a thermal plasma. The time evolution of small-size quarkonium is governed by two-point correlation functions of chromoelectric fields dressed with an adjoint Wilson line, known in this context as generalized gluon distributions (GGDs). The GGDs have been calculated in both weakly and strongly coupled plasmas by using perturbative and holographic methods. Strikingly, the results of our calculations for a strongly coupled plasma indicate that the quarkonium dissociation and recombination rates vanish in the transport descriptions that assume quarkonium undergoes Markovian dynamics. However, this does not imply that the dynamics is trivial. As a starting point to explore the phenomenological consequences of the result at strong coupling, we show a calculation of the Y(1 S ) formation probability in time-dependent perturbation theory. This is a first step towards the development of a transport formalism that includes non-Markovian effects, which, depending on how close the as of yet undetermined nonperturbative QCD result of the GGDs is to the strongly coupled N = 4 SYM result, could very well dominate over the Markovian ones in quark-gluon plasma produced at RHIC and the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Simulated CO 2 storage efficiency factors for saline formations of various lithologies and depositional environments using new experimental relative permeability data

Saline formations are attractive geologic reservoirs for permanent carbon dioxide (CO 2 ) storage. Here, the U.S. Department of Energy's National Energy Technology Laboratory (DOE-NETL) has worked to develop and refine methods and tools for the calculation of CO 2 storage potential in subsurface reservoirs. DOE-NETL's CO 2 -SCREEN provides an online tool for executing these storage methods. CO 2 storage efficiency terms are input parameters in DOE-NETL's methods and equations embedded in the CO 2 -SCREEN, which assesses pore space available for CO 2 storage. In this work, a modeling workflow was initiated to refine two CO 2 storage efficiency terms - volumetric displacement (E V ) and microscopic displacement (E d ). The models are based on new experimental relative permeability data that are specific to homogenous lithology and depositional environments of key subsurface saline formations targeted for CO 2 storage. In future work, heterogenous features will be added to this initial modeling effort to update efficiency factors as described in DOE-NETL's methods and CO 2 -SCREEN tool. E V accounts for the volume utilized in the reservoir under the areal plume, while E d accounts for saturation values in the plume to assess efficiency of CO 2 storage at the pore scale. The results of this work are significant in that prior values were based on a limited geologically non-specific relative permeability data set that were collected prior to 2009. Specifically, we applied numerical simulations using TOUGH3 models to update CO 2 storage efficiency values for supercritical CO 2 injection into brine-saturated reservoirs for three lithologies (clastics, limestone, dolomite) and six depositional environments (Marginal Marine, Strand Plain, Deltaic Complex Fluvial, Aeolian, Shallow Marine, and Reef) that have a high potential for geologic CO 2 storage. Experimental relative permeability data in cores from these environments were utilized in the models with corresponding rock type/sedimentary environment. Results of this study showed that dolomite followed by limestone generated higher ranges of storage efficiency compared to clastics. The updated values provided a tighter efficiency range for clastics, lower P 10 but higher P 90 range for limestone, and higher P 10 and P 90 for dolomite. In general, tighter reservoirs with relatively low permeability and porosity were associated with higher E V and E d , showing efficient reservoir and pore utilization in these scenarios. High reservoir pressure and temperature associated with increasing depth increased the E V , and high CO 2 injection rates resulted in increases in E V and E d , while the impact of permeability anisotropy was minimal after the 30-year injection period.

03 NATURAL GAS↗
DOE OSTI2021

Second Linear Response Theory and the Analytic Calculation of Excited-State Properties

We present a method based on second linear response time-dependent density func- tional theory (TDDFT) to calculate permanent and transition multipoles of excited states, which are required to compute excited-state absorption/emission spectra, multi- photon optical processes, among others. In previous work, we examined computations based on second linear response theory in which linear response TDDFT was employed twice. In contrast, the present methodology only requires information from a single linear response calculation to compute the excited-state properties. These are eval- uated analytically through various algebraic operations involving electron repulsion integrals and excitation vectors. The present derivation focuses on full many-body wave functions instead of single orbitals, as in our previous approach. Here, we test the proposed method by applying it to several diatomic and triatomic molecules. This shows that the computed excited-state dipoles are consistent with respect to reference equation-of-motion coupled cluster calculations.

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

Competing magnetic and nonmagnetic states in monolayer VSe 2 with charge density wave

The field of two-dimensional ferromagnets has been reinvigorated by the discovery of VSe 2 monolayer grown on van der Waals substrates, which is reported to be ferromagnetic with a Curie point higher than 330 K. However, the ferromagnetic and nonmagnetic states of pristine monolayer VSe 2 are highly debated. Here, employing density functional theory, Wannier function calculations, and the band unfolding method, we explore the electronic structure of monolayer VSe 2 with a $\sqrt{3} × \sqrt{7}$ charge density wave (CDW). Certain qualitative aspects of the calculated unfolded band dispersion and unfolded Fermi surface of monolayer VSe 2 with the $\sqrt{3} × \sqrt{7}$ CDW in the nonmagnetic state agree well with previous angle-resolved photoemission spectroscopy results, albeit with uncertainty about whether these experiments probed single or multiple domains. Specifically, we find that an isolated CDW domain naturally induces a strong breaking of the threefold symmetry of the electronic structure. In addition we find that, relative to the undistorted structure, the CDW structure shows a strong competition between nonmagnetic and various magnetic states, with an energy difference less than 5 meV/f.u. For the CDW structure in the antiferromagnetic state, the band dispersions and Fermi surface are similar to those in the nonmagnetic state, while the unfolded bands of the ferromagnetic CDW state display a sizable exchange splitting. These results indicate the possibility of various antiferromagnetic fluctuations in VSe 2 to coexist and compete with ferromagnetic order and the experimentally reported CDW order. In this work, our calculations build insights for exploring the interplay between magnetism and CDW behaviors more generally in transition metal dichalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Charged-Particle Bound States in Periodic Boxes

We consider the binding energy of a two-body system with a repulsive Coulomb interaction in a finite periodic volume. We define the finite-volume Coulomb potential as the usual Coulomb potential, except that the distance is defined the shortest separation between the two bodies in the periodic volume. We investigate this problem in one and three-dimensional periodic boxes and derive the asymptotic behavior of the volume dependence for bound states with zero angular momentum in terms of Whittaker functions. We benchmark our results against numerical calculations and show how the method can be used to extract asymptotic normalization coefficients for charged-particle bound states. Furthermore, the results we derive here have immediate applications for calculations of atomic nuclei in finite periodic volumes for the case where the leading finite-volume correction is associated with two charged clusters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Calculation of permanent magnet arrangements for stellarators: a linear least-squares method

A problem arising in several engineering areas is to design magnets outside a volume that produce a desired magnetic field inside it. One instance of this problem is stellarator design, where it has recently been shown that permanent magnets can provide the required shaping of the magnetic field. Here we demonstrate a robust and efficient algorithm REGCOIL_PM to calculate the spatial distribution of these permanent magnets. The procedure involves a small number of fixed-point iterations, with a linear least-squares problem solved at each step. Furthermore, the method exploits the Biot–Savart Law's exact linearity in magnetization density and approximate linearity in magnet size, for magnets far from the target region. No constraint is placed on the direction of magnetization, so Halbach solutions are found naturally, and the magnitude of the magnetization can be made uniformly equal to a target value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Environmental Protection Agency Drinking Water Protective Action Guides Implementation Recommendations

This report documents the analysis conducted by Sandia National Laboratories on the effect of various alterations to published methodologies to calculate Derived Response Levels for Environmental Protection Agency Drinking Water Protective Action Guides. Specifically, this study sought to assess and provide recommendations on calculation of the Derived Response Level accounting for decay during the consumption period, assess the impact of decay on laboratory Minimal Detectable Concentration in water samples as compared to the Derived Response Level, make a recommendation on the calculation of Derived Response Level consistent with existing Public Protection Methods, and make a recommendation on the use of six age groups versus eight age groups based on available dose coefficients for calculation of the Derived Response Level. The authors analyzed these various factors using nominal radionuclide mixes from four scenarios and compared calculation of the Derived Response Level accounting for decay and no decay and then compared those results to the laboratory Minimal Detectable Concentrations. The authors concluded that decay should be included in the calculation of the Derived Response Level, existing Public Protection Methods should be employed to calculate the Derived Response Level, six age groups should be used versus eight, and the use of both decay and Public Protection Methods result in little to no concern for water samples meeting the Minimum Detectable Concentration requirements. The results of this study may be used in further developing and implementing a method for the Environmental Protection Agency Water Derived Response Level calculation in the Federal Radiological Monitoring and Assessment Center Assessment Manual.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Analysis of implant loss risk factors after simultaneous guided bone regeneration: A retrospective study of 5404 dental implants

Abstract Purpose The purpose was to analyze the risk factors for implant loss after simultaneous guided bone regeneration (GBR). Materials and Methods Patients who underwent implant placement with simultaneous GBR between January 2011 and December 2018 were screened for this study. The cumulative survival rate (CSR) was calculated using the life table method. Log‐rank test and Kaplan–Meier survival estimates were used to identify potential risk factors for implant loss. The association between the investigated variables and implant loss was determined using hazard ratios (HRs) obtained from a multivariate Cox regression analysis. Results A total of 3973 patients with 5404 implants were included in this study. The CSRs of the implants at 1, 5, and 10 years were 99.6%, 98.9%, and 98.7%, respectively. Male patient (HR = 2.94, 95% CI: 1.41–6.14), periodontitis (HR = 4.26, 95% CI: 2.05–9.86), tissue‐level implants (HR = 3.02, 95% CI: 1.30–6.98), narrow implants (HR = 2.71, 95% CI: 1.12–6.57), and implant length ≤10 mm (HR = 2.91, 95% CI: 1.41–6.02) significantly increased the risk of implant loss ( p < 0.05). The risk of implant loss was significantly higher in the maxillary posterior region (HR = 2.26, 95% CI: 1.04–4.90) than in the maxillary anterior region ( p < 0.05). Compared to Straumann, Nobel (HR = 4.07, 95% CI: 1.75–9.44) and other implant systems (HR = 14.23, 95% CI: 4.32–46.85) showed a significantly higher risk of implant loss ( p < 0.05). Conclusion Male patient, periodontitis, maxillary posterior region, Nobel implant system, other implant systems, tissue‐level implants, narrow implants, and implant length ≤10 mm were considered risk factors for implant loss after simultaneous GBR.

Shen, Xiaoting↗
DOE OSTI2025

Electron-Ion Dynamics with Time-Dependent Density Functional Theory: Towards Predictive Solar Cell Modeling

This project focussed on two aspects of computational modeling with a view toward application for photovoltaic design: (i) increased reliability of exchange-correlation functionals in time-dependent density functional theory (TDDFT) (a method of choice of the calculation of electronic spectra and dynamics), especially for time-resolved non-perturbative dynamics, and (ii) the development of a practical but rigorously-based method for coupling electronic and nuclear motion via the exact-factorization (EF) approach (a relatively new framework for developing approximations).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Predicting Thermoset Resin Cure Kinetics for Composite Processing

Predicting a material’s thermal properties is not only useful for designing parts but also for defining processing parameters. This is especially applicable to wet-filament winding operations and the production of composite parts using thermoset res ins. For these, the formation of cross-linking bonds is affected by time and temperature, but the degree of cross-linking can be modeled using phenomenological methods. In the current study, the Kissinger and Flynn/Wall/Ozawa methods were both used to calculate kinetic parameters and produced values that were in agreement with each other. The kinetic parameters activation energy ( E a ) and pre-exponential factor (𝐴) calculated from exothermic peak data increased between an epoxy and resole phenol resin system, with a cyanate ester having the highest. The successful application of both methods suggests applicability to future studies for developing cure recipes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Temperature-dependent phonon lifetimes and thermal conductivity of silicon by inelastic neutron scattering and ab initio calculations

In this paper, inelastic neutron scattering on a single crystal of silicon was performed at temperatures from 100 to 1500 K. These experimental data were reduced to obtain phonon spectral intensity at all wave vectors $\overset{→}{Q}$ and frequencies ω in the first Brillouin zone. Thermal broadenings of the phonon peaks were obtained by fitting and by calculating with an iterative ab initio method that uses thermal atom displacements on an ensemble of superlattices. Agreement between the calculated and experimental broadenings was good, with possible discrepancies at the highest temperatures. Distributions of phonon widths versus phonon energy had similar shapes for computation and experiment. These distributions grew with temperature but maintained similar shapes. Parameters from the ab initio calculations were used to obtain the thermal conductivity from the Boltzmann transport equation, which was in good agreement with experimental data. Despite the high group velocities of longitudinal acoustic phonons, their shorter lifetimes reduced their contribution to the thermal conductivity, which was dominated by transverse acoustic modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Convex methods for physics

Software for performing calculations of physical systems, both equilibrium and dynamical, using convex constraints.

Scott, Lawrence↗
DOE OSTI2022

Ab initio self-consistent many-body theory of polarons at all couplings

We present a theoretical framework to describe polarons from first principles within a many-body Green’s function formalism. Starting from a general electron-phonon Hamiltonian, we derive a self-consistent Dyson equation in which the phonon-mediated self-energy is composed by two distinct terms. One term is the Fan-Migdal self-energy and describes dynamic electron-phonon processes, the other term is a contribution to the self-energy originating from the static displacements of the atomic nuclei in the polaronic ground state. The lowest-order approximation to the present theory yields the standard many-body perturbation theory approach to electron-phonon interactions in the limit of large polarons, and the ab initio polaron equations introduced in the limit of small polarons. Here, a practical recipe to implement the present unifying formalism in first-principles calculations is outlined. We apply our method to the Fröhlich model, and obtain remarkably accurate polaron energies at all couplings, in line with Feynman’s polaron theory and diagrammatic Monte Carlo calculations. We also recover the well-known results of Fröhlich and Pekar at weak and strong coupling, respectively. The present approach enables predictive many-body calculations of polarons in real materials at all couplings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

GPU Acceleration of Large-Scale Full-Frequency GW Calculations

Many-body perturbation theory is a powerful method to simulate electronic excitations in molecules and materials starting from the output of density functional theory calculations. By implementing the theory efficiently so as to run at scale on the latest leadership high-performance computing systems it is possible to extend the scope of GW calculations. Here, we present a GPU acceleration study of the full-frequency GW method as implemented in the WEST code. Excellent performance is achieved through the use of (i) optimized GPU libraries, e.g., cuFFT and cuBLAS, (ii) a hierarchical parallelization strategy that minimizes CPU-CPU, CPU-GPU, and GPU-GPU data transfer operations, (iii) nonblocking MPI communications that overlap with GPU computations, and (iv) mixed precision in selected portions of the code. A series of performance benchmarks has been carried out on leadership high-performance computing systems, showing a substantial speedup of the GPU-accelerated version of WEST with respect to its CPU version. Good strong and weak scaling is demonstrated using up to 25 920 GPUs. Finally, we showcase the capability of the GPU version of WEST for large-scale, full-frequency GW calculations of realistic systems, e.g., a nanostructure, an interface, and a defect, comprising up to 10 368 valence electrons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗