Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “attribute measurement”

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 37 records · Page 2

Wind 2H2 Project (CRADA Final Report)

The Xcel-NREL Wind Electrolysis R&D Project is designed to gain knowledge and insight about producing hydrogen from electricity for use as a vehicle fuel and as a “stored wind”. The project will analyze producing hydrogen directly from wind as well as from the electric grid. In addition, it will look at storing hydrogen and then using it to put electricity back onto the grid. Modifications were made to the original scope to also include electric vehicles integration with the electric grid. This report addresses Task 8 portion of the CRADA, focused on testing and measuring the attributes of plug-in electric vehicles.

33 ADVANCED PROPULSION SYSTEMS↗

New Capability at ORNL: High Precision Uranium Titration

Destructive analytical measurements establish a nuclear facility’s nuclear material inventory and inventory differences for Nuclear Material Accountancy and Control. A nuclear laboratory’s ability to perform accurate high-precision analytical measurements is key for tracking large inventories within a facility’s material balance areas and during production to track material movement through dynamic processes. For uranium, these measurements are made by using several established high-precision measurement protocols. These include isotope dilution mass spectrometry, gravimetry, and potentiometric titrations. Whatever measurement technique a nuclear lab chooses to use, reference materials (RMs) with certified attributes and accompanying uncertainties are used to calibrate measurement systems, and are the cornerstone for accurate results. In addition to calibration, RMs provide for metrological traceability, are used for method development and validation, and thus provide critical evaluations of the appropriateness and performance of analytical processes used. Evaluations may include validation that a method is fit-for-purpose, quantification of systematic and random biases, and the evaluation of long-term and short-term performance metrics. High-precision measurement techniques require that measured attributes be certified to a high degree of precision in the RMs used—ultimately, to a higher degree than that of the measurement technique itself. The US authority on the production of special nuclear material Certified Reference Materials is the NBL Program Office (NBLPO), formally known as New Brunswick Laboratory (NBL). The NBLPO is responsible for the sales and distribution of existing NBL certified reference materials (CRMs) and for the production of the next generation of nuclear RMs. To accomplish its mission, NBLPO is establishing key base capabilities within the DOE laboratory complex that formerly existed at the NBL laboratory. The Nuclear Analytical Chemistry (NAC) section within the Chemical Sciences Division (CSD) at Oak Ridge National Laboratory (ORNL) is currently working with NBLPO to set up laboratory and measurement capabilities to provide measurements and capabilities for production and/or recertification of existing and future CRMs for uranium assay. The NBL-developed high-precision titration (HPT) method is a critically-evaluated, extremely precise and accurate primary method utilized for the determination of uranium content in a variety of uranium materials. The HPT method, combined with detailed balance weighing protocols, provides for an analytical methodology that is unsurpassed in precision and one in which all sources of error have been evaluated, a requirement of CRM certification. HPT capability within the United States was lost with the closure of the labs at NBL. The NAC has been collaborating with NBLPO to stand up and demonstrate the capability to perform uranium assay via HPT at ORNL. HPT can produce results with an expanded uncertainty of approximately 0.01% for pure uranium compounds, with typical precisions of <0.006%. The major tasks required to stand up the method at ORNL were the refurbishment of a dedicated lab and equipment setup, procedure development, analyst training, establishing method-specific quality assurance, and qualification of an analyst to perform the method. This report summarizes these tasks, outlines the documents drafted, and gives the outcome of the qualification titrations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Geostatistical interpolation of streambed hydrologic attributes with addition of left censored data and anisotropy

Spatial geostatistical interpolation of point measurements of streambed attributes in the hyporheic zone may be constrained by the streambed anisotropy, and data density and spatial distribution may significantly impact the results. Spatial clustering and low spatial data density can be caused by bedrock outcropping at the streambed limiting installation of in-stream piezometers. This study examines parameter error variability of the geostatistical interpolation using anisotropic interpolation methods and increasing the data density by adding left censored values (i.e., data below measurement limit) to locations where measurements were limited by exposed bedrock lining the streambed. The reduction in relative standard error of the interpolation was determined for the spatial distributions of streambed attributes including hydraulic conductivity, seepage flux, and mercury solute flux measured in two different years along a study reach in East Fork Poplar Creek, Tennessee, USA. Here, two methods to impute the left censored values were compared including the conventional half the detection limit substitution method, and the Stochastic Approximation of Expectation-Maximization (SAEM) algorithm, which both had comparable results. Imputing left censored data increased the data density to recommended ranges, reduced data clustering, increased the spatial dependence for some attributes, and reduced the standard error for each of the three attributes. For the reach considered herein, addition of the left censored values resulted in a larger error reduction than the consideration of anisotropy within the interpolation, which confirms the benefit of data addition to increase data density within data-limited river corridors.

58 GEOSCIENCES↗

Pre-Transient Characterization of Historic EBR-II Pins for Transient Testing

Current interest in sodium-cooled fast reactor (SFR) designs, such as TerraPower’s Natrium Reactor, has highlighted the need for advanced reactor fuel technology development. Modern U-Zr and U- Pu-Zr pin designs are primary candidates to fuel SFRs and boast high fuel utilization capacity, increased fuel-cladding compatibility, and improved safety through inherent feedback mechanisms. Despite over 60 years of metallic fuel irradiation, uncertainties exist in the performance of the fuel system, particularly under transient overpower (TOP) and loss of flow (LOF) scenarios. Throughout historical testing within the Experimental Breeder Reactor II (EBR-II) and the Fast Flux Test Facility (FFTF), fuel behavior has demonstrated benign response to transient reactor conditions; however, accurate predictions of failure thresholds to inform operational limitations rely heavily on fuel composition, burnup, and irradiation history. In expanding TOP and LOF testing, the Transient Heat sink Overpower Response (THOR) Capsule will be used to test modern fuel technologies in a static sodium environment in the Transient Reactor Test (TREAT) Facility. The THOR capsule is highly instrumented and will provide time-dependent thermal behavior of SFR fuel pins subjected to accident conditions within TREAT. The THOR-Metallic (THOR- M) campaign aims to validate and expand historical TOP and LOF testing on high burnup U-Zr and U-Pu- Zr fuel alloys previously irradiated in EBR-II by running the rods to failure. This contribution focuses primarily on the pre-transient engineering-scale destructive and non- destructive characterization that has been conducted on both the test and sibling pins used for the TOP and LOF tests. All pins underwent visual examination, neutron radiography, element contact profilometry, and precise gamma scan. The sibling pins used for each test were further analyzed using gas assay, sampling, and recharge analysis (GASR), and optical microscopy. The results from each technique confirmed that the fuel pins were intact and devoid of any atypical developments when compared to historical data. Additionally, the analyzed measurements establish a baseline for comparison to post-transient analysis. Key fuel behaviors quanitifed include axial elongation of the fuel column, diametral strain of the pin, patterns in fluff structure geometry, changes in axial isotope distribution, evolution of constituent redistribution, porosity, and fission gas release. The pre-transient measurements and changes attributed to transient behavior from post-transient measurement will be compared to historical data to capture the behavioral dependence on composition, burnup, and irradiation history. Results from this work advance the initiatives of the THOR-M campaign, which aid in informing fuel performance models and establishing safety criteria for SFR operational limits. The novel combination of test environment, in-situ instrumentation, and comprehensive suite of characterization methods provides greater understanding of transient fuel behavior. Overall, information on the time and condition of pin failure for high burnup U-Pu-Zr will greatly expand the limited existing TOP and LOF test data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Pre-Transient Characterization of Historic EBR-II Pins for Transient Testing

Current interest in sodium-cooled fast reactor (SFR) designs, such as TerraPower’s Natrium Reactor, has highlighted the need for advanced reactor fuel technology development. Modern U-Zr and U- Pu-Zr pin designs are primary candidates to fuel SFRs and boast high fuel utilization capacity, increased fuel-cladding compatibility, and improved safety through inherent feedback mechanisms. Despite over 60 years of metallic fuel irradiation, uncertainties exist in the performance of the fuel system, particularly under transient overpower (TOP) and loss of flow (LOF) scenarios. Throughout historical testing within the Experimental Breeder Reactor II (EBR-II) and the Fast Flux Test Facility (FFTF), fuel behavior has demonstrated benign response to transient reactor conditions; however, accurate predictions of failure thresholds to inform operational limitations rely heavily on fuel composition, burnup, and irradiation history. In expanding TOP and LOF testing, the Transient Heat sink Overpower Response (THOR) Capsule will be used to test modern fuel technologies in a static sodium environment in the Transient Reactor Test (TREAT) Facility. The THOR capsule is highly instrumented and will provide time-dependent thermal behavior of SFR fuel pins subjected to accident conditions within TREAT. The THOR-Metallic (THOR- M) campaign aims to validate and expand historical TOP and LOF testing on high burnup U-Zr and U-Pu- Zr fuel alloys previously irradiated in EBR-II by running the rods to failure. This contribution focuses primarily on the pre-transient engineering-scale destructive and non- destructive characterization that has been conducted on both the test and sibling pins used for the TOP and LOF tests. All pins underwent visual examination, neutron radiography, element contact profilometry, and precise gamma scan. The sibling pins used for each test were further analyzed using gas assay, sampling, and recharge analysis (GASR), and optical microscopy. The results from each technique confirmed that the fuel pins were intact and devoid of any atypical developments when compared to historical data. Additionally, the analyzed measurements establish a baseline for comparison to post-transient analysis. Key fuel behaviors quanitifed include axial elongation of the fuel column, diametral strain of the pin, patterns in fluff structure geometry, changes in axial isotope distribution, evolution of constituent redistribution, porosity, and fission gas release. The pre-transient measurements and changes attributed to transient behavior from post-transient measurement will be compared to historical data to capture the behavioral dependence on composition, burnup, and irradiation history. Results from this work advance the initiatives of the THOR-M campaign, which aid in informing fuel performance models and establishing safety criteria for SFR operational limits. The novel combination of test environment, in-situ instrumentation, and comprehensive suite of characterization methods provides greater understanding of transient fuel behavior. Overall, information on the time and condition of pin failure for high burnup U-Pu-Zr will greatly expand the limited existing TOP and LOF test data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Prediction of Probabilistic Shock Initiation Thresholds of Energetic Materials Through Evolution of Thermal-Mechanical Dissipation and Reactive Heating

The ignition threshold of an energetic material (EM) quantifies the macroscopic conditions for the onset of self-sustaining chemical reactions. The threshold is an important theoretical and practical measure of material attributes that relate to safety and reliability. Historically, the thresholds are measured experimentally. In this work, we present a new Lagrangian computational framework for establishing the probabilistic ignition thresholds of heterogeneous EM out of the evolutions of coupled mechanical-thermal-chemical processes using mesoscale simulations. Furthermore, the simulations explicitly account for microstructural heterogeneities, constituent properties, and interfacial processes and capture processes responsible for the development of material damage and the formation of hotspots in which chemical reactions initiate. The specific mechanisms tracked include viscoelasticity, viscoplasticity, fracture, post-fracture contact, frictional heating, heat conduction, reactive chemical heating, gaseous product generation, and convective heat transfer. To determine the ignition threshold, the minimum macroscopic loading required to achieve self-sustaining chemical reactions with a rate of reactive heat generation exceeding the rate of heat loss due to conduction and other dissipative mechanisms is determined. Probabilistic quantification of the processes and the thresholds are obtained via the use of statistically equivalent microstructure sample sets (SEMSS). The predictions are in agreement with available experimental data.

36 MATERIALS SCIENCE↗

Measurement and analysis of GPU-accelerated applications with HPCToolkit

To address the challenge of performance analysis on the US DOE’s forthcoming exascale supercomputers, Rice University has been extending its HPCToolkit performance tools to support measurement and analysis of GPU-accelerated applications. To help developers understand the performance of accelerated applications as a whole, HPCToolkit’s measurement and analysis tools attribute metrics to calling contexts that span both CPUs and GPUs. To measure GPU-accelerated applications efficiently, HPCToolkit employs a novel wait-free data structure to coordinate monitoring and attribution of GPU performance. To help developers understand the performance of complex GPU code generated from high-level programming models, HPCToolkit constructs sophisticated approximations of call path profiles for GPU computations. To support fine-grained analysis and tuning, HPCToolkit uses PC sampling and instrumentation to measure and attribute GPU performance metrics to source lines, loops, and inlined code. To supplement fine-grained measurements, HPCToolkit can measure GPU kernel executions using hardware performance counters. To provide a view of how an execution evolves over time, HPCToolkit can collect, analyze, and visualize call path traces within and across nodes. Finally, on NVIDIA GPUs, HPCToolkit can derive and attribute a collection of useful performance metrics based on measurements using GPU PC samples. Here, we illustrate HPCToolkit’s new capabilities for analyzing GPU-accelerated applications with several codes developed as part of the Exascale Computing Project.

97 MATHEMATICS AND COMPUTING↗

The effect of interfacial phenomena on gas solubility measurements in molten salts

The behavior of fission gases in molten fuel salt reactors governs activity transport from the reactor and can also affect the performance of the reactor itself. The gas solubility can be described thermodynamically by Henry’s law. However, the coupling of the condensed and gas phases depends on the interfacial area, which is difficult to measure or even to estimate. Surfaces of materials in the reactor will include disperse phases in the salt and porosity within the structural materials, covering a range of compositions and sizes. These attributes can affect measurements of fundamental properties such as gas solubility. Methods to obtain gas solubility, surface tension, interfacial energies, and bubble gas transport are reviewed. Recent data from manometric experiments are interpreted based on xenon sorption onto salt-wetted quartz.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Mechanical Characterization of the NIF Ignition Target TMPSA Bonding Flexure

The Thermo-Mechanical Package Sub-Assembly (TMPSA) provides a critical mechanical, thermal, and electrical interface between a silicon arm and a TMP aluminum can in a NIF ignition target. During assembly, sixteen silicon pads are bonded to the aluminum using a fixture that positions the components and applies a repeatable prescribed displacement through a compliant flexure. The flexure converts fixture interference into displacement and reaction force. Because bondline thickness variation must be maintained within +/-4 µm, consistent flexure behavior is important to the assembly process. With a recent string of TMPSAs exhibiting low bond strength, the flexures were inspected. Despite being manufactured to the same specifications, flexures were found to exhibit variation in measured stiffness. Additionally, the measured stiffness values did not always follow the presumed beam mechanics model. This work addresses two questions: whether the bonding fixture is working as intended, and whether the previously made stiffness measurements are accurate representations of the use case. The flexures are analyzed using Euler-Bernoulli beam theory, measured beam dimensions, finite element analysis, and tolerance stack-up calculations. The analysis shows that the fixed-guided beam mechanics model appropriately represents the flexure during TMPSA bonding, but the chisel-tip stiffness measurement method introduces a deformation to the inner ring of the flexure that is not represented during use. The measured stiffness values should therefore be interpreted as test-condition stiffness values rather than direct measurements of operational flexure stiffness. The discrepancy is therefore attributed primarily to the measurement boundary condition rather than to failure of the fixed-guided beam model. Recommendations are provided for GD&T, dimensional inspection, and a representative stiffness testing method to better control bondline variation.

42 ENGINEERING↗

Evaluating proxies for the drivers of natural gas productivity using machine-learning models

We report the extensive development of unconventional reservoirs using horizontal drilling and multistage hydraulic fracturing has generated large volumes of reservoir characterization and production data. The analysis of this abundant data using statistical methods and advanced machine-learning (ML) techniques can provide data-driven insights into well performance. Most predictive modeling studies have focused on the impact that different well completion and stimulation strategies have on well production but have not fully exploited the available in situ rock property data to determine its role in reservoir productivity. We have used machine-learning techniques to rank rock mechanical properties, microseismic attributes, and stimulation parameters in the order of their significance for predicting natural gas production from an unconventional reservoir. The data for this study came from a hydraulically fractured well in the Marcellus Shale in Monongalia County, West Virginia. The data classes included measurements aggregated by well completion stage that included (1) gas production, (2) well-log-derived measurements including bulk density, elastic moduli, shear impedance, compressional impedance, brittleness, and gamma measurements, (3) microseismic attributes, (4) long-period long-duration (LPLD) event counts, (5) fracture counts, and (6) stimulation parameters that included the fluid injection volume and average pumping pressure. To identify observable proxies for the drivers of gas production, we evaluated five commonly used ML approaches including multivariate adaptive regression spline, Gaussian mixture model, random forest, gradient boosting, and neural network. We selected five variables including LPLD event count, seismogenic b-value, hydraulic diffusivity, cumulative moment, and fluid volume as the features most likely to impact gas productivity at the stage level in the study area. The data-driven selection of these parameters for their importance in determining gas production can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs. Plain language summary: We use machine-learning methods and data-driven selection of reservoir parameters to rank and better understand their importance in determining gas production, which can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs.

58 GEOSCIENCES↗

Spatiotemporal Pattern Recognition in the PMU Signals in the WECC system

Phasor measurement unit (PMU) data has been used by multiple power system applications, including state estimation, post event analysis, oscillation detection, model validation, and many others. Still, due to its big data nature and availability to general research institutions, comprehensive understanding of the spatiotemporal patterns and underlying mechanisms are incomplete. This study applies a set of signal processing and machine learning approaches aiming at deciphering the characteristic behaviors of multiple phasor measurement units (PMUs) attributes (e.g., voltage, frequency, rate of change of frequency, phase angle), including their auto-correlation, cross-dependence, similarities and discrepancies across units and temporal scales, and distributions of anomalies and their linkages to potential external factors such as weather events. Data analytics are applied to PMUs from the U.S. Western Electricity Coordinating Council (WECC) system. The PMU measurements, recorded events, outages, and weather extremes are all from real world datasets. The findings from the study and mechanistic understanding of the PMU dynamics help provide guidance on system control or preventing blackouts. The derived metrics can be directly used for adjusting or filtering simulated PMU data used for advanced algorithm development.

Hou, Zhangshuan↗

Non-resonant picosecond three-wave mixing in the gas phase

We report on the experimental observation of non-resonant, second-order optical sum-frequency generation (SFG) in five different atomic and molecular gases. The measured signal is attributed to a SFG process by characterizing its intensity scaling and its polarization behavior. We show that the electric quadrupole mechanism cannot explain the observed trends and suggest a mechanism based on symmetry breaking along the incident beam path arising from laser-induced species ground state number density gradients. Our results demonstrate that the SFG is about four orders of magnitude stronger than the third-harmonic generation (THG) and independent from any externally applied electric fields. In conclusion, these features make this method suitable for gas number density measurements at the picosecond time scale in reactive flows and plasmas.

LaCombe, Grayson↗

Applications of Measuring and Valuing Resilience in Energy Systems

The electricity sector is vulnerable to numerous hazards that are being exacerbated by climate change, which can cause an increase in the hazards' frequency and intensity. Consumers, system regulators, system operators, and communities are now preparing to mitigate the increased risks posed by climate change. New York State's energy infrastructure resilience can be increased with targeted investments including but not limited to installing emergency backup systems, integrating microgrid solutions, weatherizing buildings, increasing energy efficiency, adding redundancy, investing in restoration and recovery, and hardening critical components. Such investments can reduce the likelihood, impact, and consequences of disruptive events but can also increase capital and operating costs. A barrier to prioritizing investments in resilience is that there is no widely established method for quantifying and assigning the benefits of resilience investments across various stakeholders. Decision-makers need better information detailing the value of resilience improvements. Developing methods to quantify, value, and price resilience helps meet resilience needs in an effective manner that also supports broader societal welfare. This report lays out considerations for quantifying and valuing resilience, discusses the current state of resilience valuation tools, and provides case studies of resilience projects that demonstrate how resilience attributes could be measured while highlighting broader, project-specific challenges to increasing resilience. Further, we present insights into methods and challenges to measuring, valuing, and enacting resilience investments.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Measurement of energy resolution with the NEXT-White silicon photomultipliers

The NEXT-White detector, a high-pressure gaseous xenon time projection chamber, demonstrated the excellence of this technology for future neutrinoless double beta decay searches using photomultiplier tubes (PMTs) to measure energy and silicon photomultipliers (SiPMs) to extract topology information. This analysis uses $^{83m}$Kr data from the NEXT-White detector to measure and understand the energy resolution that can be obtained with the SiPMs, rather than with PMTs. The energy resolution obtained of (10.9 ± 0.6)%, full-width half-maximum, is slightly larger than predicted based on the photon statistics resulting from very low light detection coverage of the SiPM plane in the NEXT-White detector. The difference in the predicted and measured resolution is attributed to poor corrections, which are expected to be improved with larger statistics. Furthermore, the noise of the SiPMs is shown to not be a dominant factor in the energy resolution and may be negligible when noise subtraction is applied appropriately, for high-energy events or larger SiPM coverage detectors. These results, which are extrapolated to estimate the response of large coverage SiPM planes, are promising for the development of future, SiPM-only, readout planes that can offer imaging and achieve similar energy resolution to that previously demonstrated with PMTs.[graphic not available: see fulltext]

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Application of multi-method-multi-model inference to radiation related solid cancer excess risks models for astronaut risk assessment

The impact of including model-averaged excess radiation risks (ER) into a measure of radiation attributed decrease of survival (RADS) for the outcome all solid cancer incidence and the impact on the uncertainties is demonstrated. It is shown that RADS applying weighted model averaged ER based on AIC weights result in smaller risk estimates with narrower 95% CI than RADS using ER based on BIC weights. Further a multi-method-multi-model inference approach is introduced that allows calculating one general RADS estimate providing a weighted average risk estimate for a lunar and a Mars mission. For males the general RADS estimate is found to be 0.42% (95% CI: 0.38%; 0.45%) and for females 0.67% (95% CI: 0.59%; 0.75%) for a lunar mission and 2.45% (95% CI: 2.23%; 2.67%) for males and 3.91% (95% CI: 3.44%; 4.39%) for females for a Mars mission considering an age at exposure of 40 years and an attained age of 65 years. It is recommended to include these types of uncertainties and to include model-averaged excess risks in astronaut risk assessment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Structural, Electronic, and Thermal Properties of CdSnAs 2

Structural, electrical, and thermal properties of CdSnAs 2 , with analyses from temperature-dependent transport properties over a large temperature range, are reported. Phase-pure microcrystalline powders were synthesized that were subsequently densified to a high-density homogeneous polycrystalline specimen for this study. Temperature-dependent transport indicates n-type semiconducting behavior with a very high and nearly temperature independent mobility over the entire measured temperature range, attributed to the very small electron effective mass of this material. The Debye model was successfully applied to model the thermal conductivity and specific heat. Furthermore, this work contributes to the fundamental understanding of this material, providing further insight and allowing for investigations into altering this and related physical properties of these materials for technological applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Preface to special topic: The High Repetition Rate Frontier in High-Energy-Density Physics

High-repetition-rate (HRR) experiments can collect large datasets with high temporal, spatial, and/or parametric resolution or large numbers of repeat measurements for statistics. HRR experiments also enable new experimental designs, including active feedback control loops and novel diagnostics, that can improve the reproducibility as well as the quantity of measurements. Together, these attributes make HRR experiments ideal for performing high-quality repeatable science. Until recently, these techniques have not been applied to high-energy-density–physics (HEDP) experiments, which are typically restricted to repetition rates of a few per day. However, recent advancements in lasers, pulsed-power drivers, target fabrication, and diagnostics are starting to change this fact, opening an exciting new frontier of HRR HEDP experiments. A mini-conference on this subject at the 2021 meeting of the American Physical Society Division of Plasma Physics brought together members of this growing community. As a result, the “High Repetition Rate Frontier in High-Energy-Density Physics” special topic in Physics of Plasmas highlights current progress in this exciting area.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Phonon Olympics: Phonon property and lattice thermal conductivity benchmarking from open-source packages

Three widely used open-source packages for determining phonon properties and lattice thermal conductivities (ALAMODE, phono3py, and ShengBTE) are benchmarked by teams of expert users and the package developers. The phonons for Ge, RbBr, monolayer MoSe 2 , and AlN are modeled at zero temperature, and they scatter through three-phonon and phonon-isotope processes, with thermal conductivities obtained from the linearized Peierls–Boltzmann transport equation with input from density functional theory calculations. Over a wide range of temperatures, the thermal conductivities calculated by the teams fall within at most ±15% of their mean values for each of the four materials. The phonon frequencies, obtained from the harmonic force constants, do not show large differences between the calculations, indicating that the modal heat capacities and group velocities are not responsible for the thermal conductivity variations. It is the lifetimes associated with three-phonon scattering, obtained from the cubic force constants, that drive the variations. The many decisions required to calculate the cubic force constants (e.g., supercell size, atomic displacement, neighbor cutoff, and application of symmetries) make identification of the precise origin of the thermal conductivity variations challenging. The calculated thermal conductivities do not generally show agreement with experimental measurements, which is attributed to the limitations of the density functional theory calculations. Guidance for the development of best practices is provided, which will help to standardize protocols needed for building thermal conductivity databases. The results provide a baseline for future benchmarking of other packages and more advanced calculations.

McGaughey, Alan J. H. [Carnegie Mellon Univ., Pitt↗