SEARCH · Engineering Papers

Results for “qualitative methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Benchmarking the multiconfigurational Hartree method by the exact wavefunction of two harmonically trapped bosons with contact interaction

Highlights: • Comparison of exact and numerical single-particle density matrices. • Benchmark for multiconfigurational Hartree. • Convergence for large number of orbitals. We consider two bosons in a one-dimensional harmonic trap, interacting by a contact potential, and compare the exact solution of this problem to a self-consistent numerical solution by using the multiconfigurational time-dependent Hartree (MCTDH) method. We thereby benchmark the predictions of the MCTDH method with a few-body problem that has an analytical solution for the most commonly experimentally realized interaction potential in ultracold quantum gases. It is found that exact ground state energy and first order correlations are accurately reproduced by MCTDH up to the intermediate dimensionless coupling strengths corresponding to typical background scattering lengths of magnetically trapped ultracold dilute Bose gases. For larger couplings, established for example by (a combination of) Feshbach resonances and optical trapping, the MCTDH approach overestimates the depth of the trap-induced correlation dip of first order correlations in position space, as well as overestimates the fragmentation, defined as the average relative occupation of orbitals other than the energetically lowest one. We anticipate that qualitatively similar features in the correlation function may arise for larger particle numbers, paving the way for a quantitative assessment of the accuracy of MCTDH by experiments with ultracold atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Quality Ranking System for Molten Salt Thermal Property Data

A simple ranking system for use with data and property values compiled in the Molten Salt Thermal Database-Thermophysical Properties (MSTDB-TP) is presented. The approach includes assessments of five features of data collection: the method used to measure the property value, calibration of devices used to make the measurements, the quantified uncertainty of the results, control of the environmental conditions, and determination of the salt composition. The verifiability of how the property value was determined is also assessed based on information provided or cited within the data set. The quality of each facet is ranked as High, Moderate, or Incomplete based on objective metrics. The overall quality of the dataset is ranked based on rankings of the five data features and the reported property value, with rankings of A and B indicating the data and property value are deemed suitable for quantitative use and U indicating uncertainties that restrict the recommended use to qualitative purposes. The system is suitable for assessing data already in the MSTDB-TP and screening new entries.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

High-Energy X-ray Diffraction Microscopy for Nuclear Forensics FY2022 Project Report

Morphological information on nuclear material has been identified using visible light and scanning electron microscopy. These identify qualitative differences in particle morphology. Three-dimensional imaging of materials through alternating scanning electron microscopy imaging and focused ion beam milling has also been used. Unfortunately, these techniques are time- and labor-intensive, with significant sample preparation required and lengthy analysis times. Further, the resulting 3D images are qualitative, require manual identification, and do not capture statistically-representative populations. High energy X-ray 3D imaging using a direct-beam or diffracted-beam (High-Energy Diffraction Microscopy) have been developed at the Advanced Photon Source and can produce quantitative information on grains (phase, location, etc.) and pores (size distribution, sphericity) in a material. These techniques require only minutes to characterize a sample volume and are non-destructive, thus suitable for a wide range of existing samples and for confirmatory analyses to be carried out using conventional microscopy techniques. In this first year of the project, all uranium oxide samples were synthesized and characterized using conventional analyses by the analytical chemistry laboratory. Conventional analysis methods included powder x-ray diffraction, scanning electron microscopy, impurity analysis via inductively coupled plasma mass spectrometry, and infrared spectroscopy. Impurity analysis shows a drop in boron content from UO 3 to the lowest U 3 O 8 calcination temperature, but otherwise no appreciable difference in any sample. Analysis of diffraction data shows a flip of peaks from UO 3 dominated for the 600 °C calcined sample to U 3 O 8 dominated at 700 °C and 800 °C. Analysis of scanning electron microscopy images shows that with increased calcination temperature the size distribution of particles seems to increase and broaden. Both of these last findings are in line with previously published data, though this work used significantly fewer particles to simply show similar trends instead of getting truly quantitative particle analysis. Infrared analysis similarly shows ingrowth of U 3 O 8 as calcination temperature is increased, along with depression of peaks associated with UO 3 and water. Samples were prepared for analysis at the Advanced Photon Source at beamline 1-ID. It is anticipated that analysis will occur in November of 2022. AI/ML techniques to de-noise data coming out of 1-ID during the analyses was also developed during this time using previously gathered data. Preliminary results using a self-supervision technique called Noise2Selfshow good de-noising of data. Once the uranium oxide samples are analyzed, real data will be used to test the de-noising and other AI/ML techniques that may be developed in the second year of the project.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1996

Low-Dimensional Dynamical Models of Thermal Convection

A low-dimensional dynamic model for transitional buoyancy-driven flow in a differentially heated tall enclosure is presented. The full governing partial differential equations with the associated boundary conditions are solved by a spectral element method for a cavity of aspect ratio A=20. Proper orthogonal decomposition is applied to the oscillatory solution at Prandtl number Pr=P tau (omega) = 0.71 and Grashof number G tau (omega) = 3.2 x 10 (exp 4) to construct empirical eigenfunctions. Using the four most energetic empirical eigenfunctions for the velocity and temperature as basis functions and applying Galerkin's method, a reduced model consisting of eight nonlinear ordinary differential equations is obtained. Close to the 'design' conditions (P tau(omega) G tau(omega)), the low-order model (LOM) predictions are in excellent agreement with the predictions of the full model. In particular, the critical Grashof number at the onset of the first temporal flow instability (Hopf bifurcation) was well as the frequency and amplitude of oscillations at supercritical conditions are in excellent agreement with the predictions of the full model. Far from the 'design' conditions, the LOM predicts the existence of multiple stable steady solutions at large values of G tau, and a unique stable steady solution at small values of G tau, and exhibits hysteretic behavior that is qualitatively similar to that observed in direct numerical simulations based on the full model.

Liakopoulos, Anthony↗

NASA NTRS2005

Quantitative Analysis of Venus Radar Backscatter Data in ArcGIS

Ongoing mapping of the Ganiki Planitia (V14) quadrangle of Venus and definition of material units has involved an integrated but qualitative analysis of Magellan radar backscatter images and topography using standard geomorphological mapping techniques. However, such analyses do not take full advantage of the quantitative information contained within the images. Analysis of the backscatter coefficient allows a much more rigorous statistical comparison between mapped units, permitting first order selfsimilarity tests of geographically separated materials assigned identical geomorphological labels. Such analyses cannot be performed directly on pixel (DN) values from Magellan backscatter images, because the pixels are scaled to the Muhleman law for radar echoes on Venus and are not corrected for latitudinal variations in incidence angle. Therefore, DN values must be converted based on pixel latitude back to their backscatter coefficient values before accurate statistical analysis can occur. Here we present a method for performing the conversions and analysis of Magellan backscatter data using commonly available ArcGIS software and illustrate the advantages of the process for geological mapping.

Long, S. M.↗

DOE OSTI2020

Full circle mechanical dynamic characterization including experimental modal analysis and finite element analysis

During operation, it was observed that a specific mechanical system experienced undesirable vibration and it became necessary to understand and mitigate this phenomenon. This document investigates the tools, methodology, and results of the dynamic characterization of the system. The characterization makes use of the experimental modal analysis (EMA) methods of single input multiple output (SIMO) and single input single output (SISO). The validity of the theory of reciprocity is confirmed to minimize measurement error, cost, and time of repeat testing. Finite element analysis (FEA) is used in choosing transducer and modal impact locations to adequately characterize the system. Single degree of freedom (SDOF) and multiple degree of freedom (MDOF) curve fitting is used to fully characterize the system’s mode shapes and natural frequencies. The EMA characterization results are used to modify and validate the FEA model so that FEA can be used to model potential structural modifications to the system to mitigate the undesirable vibration. Structural modifications are chosen, implemented, and their effectiveness is quantified using EMA. Finally, a qualitative evaluation of the methodology of FEA validation by EMA and tuning of the model to match the experimental results is discussed.

42 ENGINEERING↗

DOE OSTI2023

Characterization of Fe-6Si Soft Magnetic Alloy Produced by Laser-Directed Energy Deposition Additive Manufacturing

Abstract Commercial electrical steels, Fe-Si alloys with < 4 wt.% Si, are inexpensive and efficient materials for electrical power conversion. Further efficiency improvements require increasing the silicon concentration to 6 wt.%, at which point the material becomes brittle and difficult to form by conventional rolling and sheet fabrication methods. Additive manufacturing stands to overcome challenges with commercial manufacturing techniques by leveraging near-net-shape fabrication. The wide array of process conditions provides additive manufacturing with increased flexibility, enabling control over the microstructure and mechanical properties. This work explores the microstructures and magnetic properties of ring-shaped Fe-Si alloys produced using concentric and cross-hatch tool paths on a laser-directed energy deposition additive manufacturing system. Concentric-built samples exhibit elongated grain structures while cross-hatch-built samples comprise lower aspect ratio grain structures. Thermal finite element analysis simulations model the stress conditions produced by the different scan path geometries. Microhardness measurements probe the mechanical properties as a function of anneal temperature, providing a qualitative understanding of the intergranular defect density. Soft magnetic properties measured under quasistatic and AC conditions show frequency- and microstructure-dependent coercivity and permeability. Finally, analysis of the core loss quantifies how the build strategies and thermal treatments influence efficiency in electrical power conversion applications. Understanding the influences of scan path geometry and thermal treatment provides a pathway towards application of additively manufactured soft magnetic materials.

Adamczyk, Jesse M.↗

NASA NTRS2015

L'CO/LFIR Relations with CO Rotational Ladders of Galaxies Across the Herschel SPIRE Archive

We present a catalog of all CO (carbon monoxide) (J (total angular momentum quantum number) equals 1-0 through J equals 13-12), [CI], and [NII] lines available from extragalactic spectra from the Herschel SPIRE (Spectral and Photometric Imaging Receiver) Fourier Transform Spectrometer (FTS) archive combined with observations of the low-J CO lines from the literature and from the Arizona Radio Observatory. This work examines the relationships between L (sub FIR (Far Infra Red)), L prime (sub CO), and L (sub CO) / L (sub CO, 1-0). We also present a new method for estimating probability distribution functions (PDFs) from marginal signal-to-noise ratio Herschel FTS spectra, which takes into account the instrumental "ringing" and the resulting highly correlated nature of the spectra. The slopes of Log (L (sub FIR)) versus Log (L prime (sub CO)) are linear for all mid- to high-J CO lines. The mid- to high-J CO luminosity relative to CO J equals 1-0 increases with increasing L (sub FIR), indicating higher excitement of the molecular gas, though these ratios do not exceed approximately 100. The luminosities relative to CO J equals 1-0 remain relatively at from J equals 6-5 through J equals13-2, across many orders of magnitude of L (sub FIR). Qualitative comparisons to current theoretical models do not match these flat SLED (spectral line energy distributions) shapes, indicating the need for more comprehensive modeling of the excitation processes of warm molecular gas in nearby galaxies.

Submillimeter Data↗

DOE OSTI2023

Development of fracture diagnostic methods for fluid distribution based on quantitative interpretation of distributed acoustic sensing and distributed temperature sensing

Multistage hydraulic fracturing design on horizontal wells has significantly evolved with larger fluid volume, more fracturing stages, and tighter perforation cluster spacing to efficiently stimulate unconventional reservoirs. From the published field observations, the recent fracturing design results in complex fracture networks or swarm of fractures. Fracture treatment evaluation is extremely challenging in such a case because of the large amount of variables in well completion and stimulation design. Combined measurements from different technologies can help in fracture diagnosis. Fluid distribution, either during fracture injection or during production, directly relates to the stimulation efficiency at the cluster level and at the stage level. Because it is unlikely in the real world to distribute the injected fluid uniformly among all the clusters, we need diagnostic techniques to generate the flow profile along a lateral. Fiber-optic measurements, such as distributed acoustic sensing (DAS) and distributed temperature sensing (DTS), are currently used to diagnose downhole flow conditions. This technology allows us to qualitatively confirm the fluid flow profile and other issues occurring downhole during fracturing such as leakage through plugs. For optimizing a fracturing design, we also need to understand how the design parameters are correlated with the stimulation efficiency. In this study, we combine two sets of models of DAS and DTS data interpretation for injected fluid volume distribution. The DAS is interpreted based on an empirical correlation between fluid flow rates and frequency band energy from the acoustic signals. The DTS is interpreted by performing temperature history match-based thermal energy conservation. Because of the completely different physics behind the interpretations, the confirmation of two interpretations provides confidence in fluid distribution.

Geochemistry & Geophysics↗

DOE OSTI2021

Application of molecular convergent close-coupling cross sections in a collisional radiative model for the triplet system of molecular hydrogen

Collisional radiative (CR) models for molecular hydrogen are of high relevance for performing qualitative and quantitative analysis of excited-state population densities measured in plasmas or predicting the dependence of plasma emission on parameter variations. Although the development of such models for H 2 started decades ago, major uncertainties still exist regarding the most important set of input parameters, namely the cross sections for electron-impact excitation. The deviations between cross sections from different datasets are particularly pronounced in the energy region close to the threshold energy, strongly increasing the uncertainty of CR models applied to low-temperature plasmas. This paper presents experimental validation of a set of newly calculated non ro-vibrationally resolved electron-impact cross sections calculated for the triplet system of H 2 using the molecular convergent close-coupling method in the adiabatic-nuclei formulation. These cross sections are implemented into a CR model based on the flexible solver Yacora. A first comparison of CR calculations with the different datasets to experimentally-determined population densities is performed at a planar ICP discharge for varying pressure (between 1 and 10 Pa) and RF power (between 700 and 1100 W). For the experimentally-accessible electron temperature and density range (2.5–10 eV and 1.8–3.3 × 10 16 m –3 , respectively), very good agreement between the model and experiment is obtained using the new data set, in contrast to previously used cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Examination of synthetic gas puff imaging diagnostic data from a gyrokinetic turbulence code

A synthetic gas puff imaging (GPI) diagnostic has been developed for the purpose of validating the three-dimensional gyrokinetic turbulence code XGC. The synthetic diagnostic is described and applied to XGC simulations of two Alcator C-Mod discharges. The turbulence characteristics deduced from the resulting simulated GPI frames, using analysis techniques similar to those applied to experimental data, are compared with locally derived characteristics extracted directly from the XGC output. The comparison of the two is shown to be potentially impacted by misalignment between the GPI view and the magnetic field, the dependence of the light emission on the electron density and temperature, and spatial and temporal variations in the neutral gas cloud induced by the turbulent plasma fluctuations. We conclude that quantitative and, in some cases even qualitative, validation of turbulence simulations need to account for these effects. While we cannot directly compare our results with experimental data due to the absence of high quality GPI data from the shots and times simulated by XGC, we do relate the overall characteristics obtained from the synthetic GPI analysis to published Alcator C-Mod GPI data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Orbital dependent complications for close vs well-separated electrons in diradicals

We investigate two limits in open-shell diradical systems: O3, in which the interesting orbitals are in close proximity to one another, and (C21H13)2, where there is a significant spatial separation between the two orbitals. In accord with earlier calculations, we find that standard density-functional approximations do not predict the open-shell character for the former case but uniformly predict the open-shell character for the latter case. We trace the qualitatively incorrect behavior in O3 predicted by these standard density functional approximations to self-interaction error and use the Fermi–Löwdin-orbital-self-interaction-corrected formalism to determine accurate triplet, closed-shell singlet, and open-shell broken-spin-symmetry electronic configurations. Analysis of the resulting many-electron overlap matrices allows us to unambiguously show that the broken-spin-symmetry configurations do not participate in the representation of the Ms = 0 triplet states and allows us to reliably extract the singlet–triplet splitting in O3 by analyzing the energy as a function of Fermi-orbital-descriptor permutations. The results of these analyses predict the percentage of open-shell character in O3, which agrees well with conventional wavefunction-based methods. While these techniques are expected to be required in cases near the Coulson–Fischer point, we find that they will be less necessary in diradical systems with well-separated electrons, such as (C21H13)2. Results based on energies from self-interaction-corrected generalized gradient, local density, and Hartree–Fock approximations and experimental results are in generally good agreement for O3. These results help form the basis for deriving extended Heisenberg-like Hamiltonians that are needed for descriptions of molecular magnets when there are competing low-energy electronic configurations.

Chemistry↗

NASA NTRSDate not supplied

An Eddy Current Case Study using NASA’s Transfer Function and Limited-Sample Probability of Detection Guidebooks

Probability of detection (POD) requirements for NASA fracture-critical human-rated systems can be met by (1) claiming ‘similarity’ to NASA Standard Nondestructive Evaluation (NDE) flaw sizes or (2) by inspector-specific POD demonstration, known as NASA Special NDE. Standard NDE flaw sizes are intended to represent the detection capability of most qualified inspectors, and were primarily based on POD studies in the development of the Space Shuttle Program Orbiter fracture control plan. However, these Standard NDE flaw sizes were based on fatigue cracks in flat panels of one material alloy over a limited range of panel thicknesses. While the POD study was comprehensive, the inspection conditions are relatively narrow. An evaluation of similarity between the Standard NDE POD study and flight component conditions is required to apply Standard NDE to specific flight component inspections. A similarity evaluation may include differences in materials, surface finish, component geometry (e.g., curvature, corners, welds), and inspection access. Similarity is typically based on a qualitative engineering evaluation. Recently, NASA published a quantitative methodology to assess similarity. If similarity is not supported by test and analysis, then the methodology provides a transfer function of flaw size for Special NDE POD demonstration. NASA recommends that Special POD demonstrations use flaws that are flight-like or more difficult to detect than flight component flaws, which are commonly induced fatigue cracks. Often, it is not feasible to induce fatigue cracks in flight component geometry, and therefore, a simpler representative geometry is used, commonly a flat panel, for POD demonstration. Inducing cracks in simple geometry specimens can be time-consuming and expensive, and thus, minimizing the number of flawed specimens is desired. Traditionally, NASA Special NDE demonstration was limited to MIL-HDBK-1823A methods or the binomial point-estimate method (PEM), commonly known as 29/29. Recently, NASA published a limited-sample POD (LS-POD) methodology for signal-response NDE techniques (e.g., eddy current) that provides POD demonstration with fewer specimens than previous approaches. LS-POD suggests a minimum of 10 flawed specimens, compared to 29 in the PEM and the suggested minimum of 40 from MIL-HDBK-1823A. The methodology also includes new guidance on evaluating the probability of false calls. In this presentation, these new transfer function and LS-POD methodologies are described and illustrated through an eddy current detection capability of far-side flaws on thin welded panels. A geometry-based transfer function is developed using electro-discharge machined (EDM) notches in flight component and simple geometry to estimate the Special NDE flaw size that is representative of the critical initial flaw size required by the fracture analysis. Then, fatigue cracks are induced in simple flat panels of the transferred flaw size, and the LS-POD methodology is used for Special NDE demonstration.

P. A. Parker↗

DOE OSTI2025

De Nive Quinquangula: Pentagonal Snowflake Metal Nanocrystals as an Example of Emergence Phenomenon in Nanotechnology Complex Systems

The formation of snowflakes (SFL) is one of the most captivating phenomena in nature. The fabrication of multimetallic gold (Au) nanoparticles containing copper (Cu) and iron (Fe) is reported, which adopt snowflake-like morphology after nanoparticle aggregation. The synthesis of these snowflake-like Au microcrystals occurs mainly by kinetic control. In this process, a surfactant-assisted salt reduction method is employed on a heated substrate. Here, the metal salts react within seconds to form nanoparticles. This process results in a hierarchical structure with pseudo-pentagonal symmetry and chiroptical activity. We discuss how large superstructures emerge from nanoparticle aggregation, resembling natural snowflakes. However, unlike ice snowflakes, which exhibit hexagonal symmetry, the Au snowflakes display pentagonal symmetry. It is proposed that the formation of these snowflake-like microcrystals is the result of the complexity involved in the crystal growth process and represents an example of the phenomenon of emergence, which has become very significant in several fields of modern science, including physics, biology, chemistry, economics, philosophy, and poetry. The term emergent is used to evoke the collective behavior of a large number of microscopic constituents that is qualitatively different than the behaviors of the individual constituents.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Andreev reflection at the altermagnet-superconductor interface

Altermagnets are a new class of magnetic materials, which exhibit large spin splitting, but due to the combined spin and real space group symmetry protection maintain zero net macroscopic magnetization. Such a characteristic may prove them to be superior in applications in superconducting heterostructures and thus here we investigate the Andreev reflection at the altermagnet/superconductor interface. We compare and contrast altermagnets to other magnetic materials, revealing qualitative differences in the behavior of altermagnetic junction depending on the Fermi surface orientation. We study the resonant states arising in setups with strong tunneling barriers and show that sensitivity to nonmagnetic disorder is also dependent on the orientation. In conclusion, our results provide a building block for altermagnetic superconducting heterostructures such as Josephson π junctions with superior properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1988

Liquid sprays and flow studies in the direct-injection diesel engine under motored conditions

A two dimensional, implicit finite difference method of the control volume variety, a two equation model of turbulence, and a discrete droplet model were used to study the flow field, turbulence levels, fuel penetration, vaporization, and mixing in diesel engine environments. The model was also used to study the effects of engine speed, injection angle, spray cone angle, droplet distribution, and intake swirl angle on the flow field, spray penetration and vaporization, and turbulence in motored two-stroke diesel engines. It is shown that there are optimum conditions for injection, which depend on droplet distribution, swirl, spray cone angle, and injection angle. The optimum conditions result in good spray penetration and vaporization and in good fuel mixing. The calculation presented clearly indicates that internal combustion engine models can be used to assess, at least qualitatively, the effects of injection characteristics and engine operating conditions on the flow field and on the spray penetration and vaporization in diesel engines.

Nguyen, Hung Lee↗

NASA NTRS1997

Ultrasonic Nondestructive Characterization of Adhesive Bonds

Qualitative measurements of adhesion or binding forces can be accomplished, for example, by using the reflection coefficient of an ultrasound or by using thermal waves (Light and Kwun, 1989, Achenbach and Parikh, 1991, and Bostrom and wickham, 1991). However, a quantitative determination of binding forces is rather difficult. It has been observed that higher harmonics of the fundamental frequency are generated when an ultrasound passes through a nonlinear material. It seems that such non-linearity can be effectively used to characterize the bond strength. Several theories have been developed to model this nonlinear effect (Adler and Nagy, 1991; Achenbach and Parikh, 1991; Parikh and Achenbach, 1992; and Hirose and Kitahara, 1992; Anastasi and Roberts, 1992). Based on a microscopic description of the nonlinear interface binding force, a quantitative method was presented by Pangraz and Arnold (1994). Recently, Tang, Cheng and Achenbach (1997) made a comparison between the experimental and simulated results based on this theoretical model. A water immersion mode-converted shear wave through-transmission setup was used by Berndt and Green (1997) to analyze the nonlinear acoustic behavior of the adhesive bond. In this project, the nonlinear responses of an adhesive joint was investigated through transmission tests of ultrasonic wave and analyzed by the finite element simulations. The higher order harmonics were obtained in the tests. It is found that the amplitude of higher harmonics increases as the aging increases, especially the 3dorder harmonics. Results from the numerical simulation show that the material nonlinearity does indeed generate higher order harmonics. In particular, the elastic-perfect plastic behavior generates significant 3rd and 5th order harmonics.

Qu, Jianmin↗

DOE OSTI2020

Microscopic Mechanisms of Glasslike Lattice Thermal Transport in Cubic Cu 12 Sb 4 S 13 Tetrahedrites

Materials based on cubic tetrahedrites (Cu 12 Sb 4 S 13 ) are useful thermoelectrics with unusual thermal and electrical transport properties, such as very low and nearly temperature-independent lattice thermal conductivity (κ L ). We explain the microscopic origin of the glasslike κ L in Cu 12 Sb 4 S 13 by explicitly treating anharmonicity up to quartic terms for both phonon energies and phonon scattering rates. We show that the strongly unstable phonon modes associated with trigonally coordinated Cu atoms are anharmonically stabilized above approximately 100 K and continue hardening with increasing temperature in accord with experimental data. This temperature-induced hardening effect reduces scattering of heat carrying acoustic modes by reducing the available phase space for three-phonon processes, thereby balancing the conventional ∝ T increase in scattering due to phonon population and yielding nearly temperature independent κ L . Furthermore, we find that very strong phonon broadening leads to a qualitative breakdown of the conventional phonon-gas model and modify the dominant heat transport mechanism from the particlelike phonon wave packet propagation to incoherent contributions described by the off-diagonal terms in the heat-flux operator, which are typically prevailing in glasses and disordered crystals. Furthermore, our work paves the way to a deeper understanding of glasslike thermal conductivity in complex crystals with strong anharmonicity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗