Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface roughness”

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

Effects of micrometer-scale surface roughness on thermal infrared emittance spectra of silica glass

Surface roughness is known to decrease thermal infrared (TIR) absorption band intensity, but studies of the effect on geologically relevant samples are relatively limited. To determine the effect of surface roughness (with features smaller than ~2/3 of the wavelength) on TIR spectra, we investigated two glass compositions with prepared roughened surfaces: (1) high purity fused silica and (2) soda-lime glass (73 wt% SiO$_2$). We roughened the surfaces of the glasses by sandblasting and polishing with grit paper. The surfaces were characterized with scanning electron microscopy and stylus profilometry. We then analyzed the roughened glasses with TIR emittance spectroscopy. Micrometer-scale roughness causes a decrease in TIR absorption band intensity, relative to a specular surface. No significant changes in band shape or shifts in wavelength were detected. As roughness increases, empirical results show a logarithmic decrease in TIR absorption band intensity. The logarithmic trends of the two glass compositions are different; empirical roughness calibrations do not translate across different compositions. A linear, least-squares spectral deconvolution using two endmembers, specular and blackbody, predicts model spectra of roughened glass surfaces with relatively low error. This is of consequence to orbital TIR measurements of poorly constrained targets, such as the martian surface, because micrometer-scale roughness is adequately modeled by the addition of a blackbody spectrum to the deconvolution endmember matrix.

36 MATERIALS SCIENCE↗

Designing slippery rough surfaces to enhance dropwise condensation of low surface tension fluid

Enhancing low surface tension liquid condensation is critical for achieving high energy efficiency and reducing the size of thermal energy systems. Extensive research has focused on promoting dropwise condensation of these liquids using state-of-the-art coatings on plain surfaces. However, maintaining dropwise condensation with low surface tension fluids is challenging due to rivulet formation, resulting in wetted tails that transition to filmwise condensation at elevated heat fluxes. Here, to address this issue, we uncover the role of surface structures and surface chemistry in the dropwise condensation of low surface tension ethanol on slippery rough surfaces (SRS). High-performance dropwise condensation has been achieved on slippery microchannels grafted with perfluoropolyether. The SRS uniquely facilitates rapid lateral droplet removal, enabling faster directional droplet shedding without rivulet formation. The resulting higher droplet removal frequency on SRS leads to heat transfer coefficients 100 % and 500 % higher than conventional dropwise and filmwise condensation on plain surfaces, respectively. Our findings uncover the pivotal role of rapid droplet removal through slippery microchannels in sustaining dropwise condensation of low surface tension liquids. This study introduces a new paradigm for promoting dropwise condensation using engineered SRS that incorporates surface structure and surface chemistry. The work will provide fundamental design guidelines to design efficient and compact condensers that use costly or flammable low global warming potential refrigerants in future refrigeration systems.

Dropwise condensation↗

Surface Roughness Effects from Additive Manufacturing in High Efficiency Gas Turbine Combustion Systems

In the past, clean combustion systems were characterized by high system development costs due to growing complexity and escalating manufacturing costs from conventional manufacturing processes. As a result, high efficiency concepts are difficult to design and more difficult to be cost-effectively manufactured. In recent years, with the introduction of additive manufacturing (AM) technologies, the rapid prototyping and mass production processes of clean combustion systems are promising to be significantly simplified with significant reduction in terms of engine manufacturing cost and engine energy cost. Compared with conventional manufacturing built parts, the AM process enabled a simpler design to be adopted for the nozzle, reducing the number of required braze and weld joints from twenty-five to just five. The resulting nozzle was 25% lighter and five times more durable and contributed to a 15% reduction in fuel burn in comparison with the previous model produced. Applying these improvements to a fleet of 5,000 turbofan engines at 150 kN thrust would result in fuel cost savings of $ \$ $6B annually (at $ \$ $5/gallon Jet-A fuel price, 700 gallon/hour fuel consumption rate, and 2,300 operational hours per year), and reduce CO2 emissions by more than 25 million tons per year. These advances support core Department of Energy (DOE) missions in energy efficiency, improving productivity, and environmental sustainability. Maximizing the benefit of these new capabilities will require high prediction capability of high speed turbulent flow with wall modeled Large Eddy Simulation (LES). GE Aviation maintains that advanced simulation technology and supercomputing is required in order to provide the appropriate boundary conditions to realize the maximum potential of AM and Ceramic Matrix Composite to reduce cooling flow, a first order penalty on the Brayton cycle. The algorithm developed using the ANSYS/Fluent software would provide the foundation for entirely new avenues of research and development with potential multi-billion dollar impact to the US economy, and significant reduction in carbon based emissions across the aerospace and power generation industries. The formulation of this iWLES (integral wall model for LES) is generic and allows to capture the changes in flow dynamics that have a significant impact on the wall bounded flow characteristics, such as swirler effective area, bulk swirl number, local pressure distribution, exit velocity profile, and turbulence kinetic energy profile. A periodic channel flow with rough flat plate is simulated using LES (Wall-Adapting Local Eddy-viscosity) model to verify the implementation of iWLES in the Fluent User Defined Function (UDF). As the flow fields are highly sensitive due to surface roughness of the wall bounded flows in the combustion systems, there is significant potential to conduct further LES studies focusing on the turbulence boundary layer interaction. Accurately capturing near wall physics will elucidate the impact of rough surfaces on combustor flow and aero-thermal interactions.

42 ENGINEERING↗

Surface Roughness Effects from Additive Manufacturing in High Efficiency Gas Turbine Combustion Systems

In the past, clean combustion systems were characterized by high system development costs due to growing complexity and escalating manufacturing costs from conventional manufacturing processes. As a result, high efficiency concepts are difficult to design and more difficult to be cost-effectively manufactured. In recent years, with the introduction of additive manufacturing (AM) technologies, the rapid prototyping and mass production processes of clean combustion systems are promising to be significantly simplified with significant reduction in terms of engine manufacturing cost and engine energy cost. Compared with conventional manufacturing built parts, the AM process enabled a simpler design to be adopted for the nozzle, reducing the number of required braze and weld joints from twenty-five to just five. The resulting nozzle was 25% lighter and five times more durable and contributed to a 15% reduction in fuel burn in comparison with the previous model produced. Applying these improvements to a fleet of 5,000 turbofan engines at 150 kN thrust would result in fuel cost savings of $6B annually (at $5/gallon Jet-A fuel price, 700 gallon/hour fuel consumption rate, and 2,300 operational hours per year),and reduce CO2 emissions by more than 25 million tons per year. These advances support core Department of Energy (DOE)missions in energy efficiency, improving productivity, and environmental sustainability. Maximizing the benefit of these new capabilities will require high prediction capability of high speed turbulent flow with wall modeled Large Eddy Simulation (LES). GE Aviation maintains that advanced simulation technology and supercomputing is required in order to provide the appropriate boundary conditions to realize the maximum potential of AM and Ceramic Matrix Composite to reduce cooling flow, a first order penalty on the Brayton cycle. The algorithm developed using the ANSYS/Fluent software would provide the foundation for entirely new avenues of research and development with potential multi-billion dollar impact to the US economy, and significant reduction in carbon based emissions across the aerospace and power generation industries. The formulation of this iWLES (integral wall model for LES) is generic and allows to capture the changes in flow dynamics that have a significant impact on the wall bounded flow characteristics, such as swirler effective area, bulk swirl number, local pressure distribution, exit velocity profile, and turbulence kinetic energy profile. A periodic channel flow with rough flat plate is simulated using LES (Wall-Adapting Local Eddy-viscosity) model to verify the implementation of iWLES in the Fluent User Defined Function (UDF). As the flow fields are highly sensitive due to surface roughness of the wall bounded flows in the combustion systems, there is significant potential to conduct further LES studies focusing on the turbulence boundary layer interaction. Accurately capturing near wall physics will elucidate the impact of rough surfaces on combustor flow and aero-thermal interactions.

99 GENERAL AND MISCELLANEOUS↗

Dendrite-free Al recycling via electrodeposition using ionic liquid electrolytes: The effects of deposition temperature and cathode surface roughness

In this report, the electrodeposition of Al from aluminum scrap alloys (A2020) on copper cathode substrates with varied surface roughness under different deposition temperatures was studied using low-temperature AlCl 3 -1-butyl-3-methyl-imidazolium chloride (BMIC) ionic liquid electrolytes. The bulk electrodeposition of Al was carried out under a voltage of 1.5 V at a stirring rate of 120 rpm using a fixed ionic liquid electrolyte concentration (molar ratio AlCl 3 : BMIC = 2:1). The effects of deposition temperature (range from 80 °C to 140 °C) and surface roughness of Cu cathode substrates (polished by 320, 600, 800, 1200 grits SiC sandpapers and mirror polishing process) on the morphology of deposited Al, current density, current efficiency and energy consumption, were investigated. The Al deposits were characterized using scanning electron microscopy (SEM) with energy-dispersive X-ray spectroscopy (EDS), X-ray diffraction (XRD), profilometer, and electrochemical measurements for current density, current efficiency, and energy consumption. It is demonstrated that the deposition temperature and surface roughness of Cu electrodes play a critical role in the nucleation and growth of Al deposits. Higher deposition temperature promotes the diffusion and/or migration of Al 2 Cl 7 - ions and then enhances the current density and efficiency during the electrodeposition of Al. Smoother surface of Cu electrodes is preferred for the formation of dendrite-free Al deposits. Typically, on the mirror polished Cu electrode, no Al dendrite structure was observed, and only plate-like Al deposits were formed at the deposition temperature of 100 °C. Pure metallic Al was successfully deposited on Cu electrodes in AlCl 3 + BMIC ionic liquid electrolytes for all experiments with a current efficiency range from 72% to 99% and energy consumption of 4.6–6.3 kW h/kg Al.

54 ENVIRONMENTAL SCIENCES↗

Investigating Future Arctic Sea Ice Loss and Near–Surface Wind Speed Changes Related to Surface Roughness Using the Community Earth System Model

The Arctic is undergoing a pronounced and rapid transformation in response to changing greenhouse gasses, including reduction in sea ice extent and thickness. There are also projected increases in near-surface Arctic wind. This study addresses how the winds trends may be driven by changing surface roughness and/or stability in different Arctic regions and seasons, something that has not yet been thoroughly investigated. We analyze 50 experiments from the Community Earth System Model Version 2 (CESM2) Large Ensemble and five experiments using CESM2 with an artificially decreased sea ice roughness to match that of the open ocean. We find that with a smoother surface there are higher mean wind speeds and slower mean ice speeds in the autumn, winter, and spring. The artificially reduced surface roughness also strongly impacts the wind speed trends in autumn and winter, and we find that atmospheric stability changes are also important contributors to driving wind trends in both experiments. In contrast to the clear impacts on winds, the sea ice mean state and trends are statistically indistinguishable, suggesting that near-surface winds are not major drivers of Arctic sea ice loss. Two major results of this work are: (a) the near-surface wind trends are driven by changes in both surface roughness and near-surface atmospheric stability that are themselves changing from sea ice loss, and (b) the sea ice mean state and trends are driven by the overall warming trend due to increasing greenhouse gas emissions and not significantly impacted by coupled feedbacks with the surface winds.

54 ENVIRONMENTAL SCIENCES↗

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗

On the Counter‐intuitive Heterogeneous Electron Transfer Barrier Properties of Alkanethiolate Monolayers on Gold: Smooth versus Rough Surfaces

Abstract Alkanethiolate monolayers formed on rough gold surfaces can, somewhat surprisingly, act as stronger barriers to heterogeneous electron transfer than those on smooth gold surfaces. This paper presents a possible explanation for this observation by constructing simple geometric models of a “rough” and “smooth” gold surface to examine how microscopic roughness differences can affect the nucleation/growth of the adlayer and size/density of structural defects. Expectedly, the number of defects predicted for adlayers formed on smooth gold is lower than any of those for rough gold. The counter‐intuitive result is that the sizes of a small portion of the defects in the adlayer on the smooth surface are larger than any of those found on the rough surface. The potential implications of these results are discussed.

Chemistry↗

Fabric controls on fracture surface roughness of an architected rock material

Fluid flow through fractures is intimately linked to the fracture surfaces that define the void geometry through which fluids flow. Thus, an understanding of what controls fracture surface roughness is essential to the development of models for predicting fluid transport through fractured rock. The difficulty in predicting surface roughness arises from the complexity of rock which is inherently heterogeneous and nonuniform in composition, fabric, and structural components, even when samples are acquired from the same rock mass. Here, a benchmarked-simulation approach motivated from geo-architected 3D printed synthetic gypsum rocks is used to provide insight into the competing contributions from fabric and layering on fracture roughness formation. Simulation results from a discrete element model (Particle Flow Code, Itasca Consulting Group, Inc.) clearly indicate that the relative orientation between mineral layers and in-layer mineral fabric, and the variability in mineral bonding strengths determine whether anisotropic corrugated surfaces or isotropic surfaces are formed. Weak mineral layers oriented perpendicular to the applied load resulted in strong roughness anisotropy. Peak failure loads were found to vary up to 30% depending on the strength of the mineral fabric at the location of fracture initiation, which provides insight into the observed high variability in strength values of natural rock. The uniqueness of induced fracture roughness and peak failure load is intimately linked to layering, mineral fabric, and their distribution in the rock. These findings have important implications for any architected material fabricated through serial printing of layers with local compositional heterogeneity.

3D printed rock↗

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part II: The Effect of Roughness on the Performance of a High-Pressure Turbine Blade

Abstract Turbine blades generally present surface roughness introduced in the manufacturing process or caused by in-service degradation, which can have a significant impact on aero-thermal performance. A better understanding of the fundamental physical mechanisms arising from the interaction between the roughness and the turbine flow at engine-relevant conditions can provide insights for the design of blades with improved efficiency and longer operational life. To this end, a high-fidelity numerical framework combining a well-validated solver for direct numerical simulation and a second-order accurate immersed boundary method is employed to predict roughness-induced aero-thermal effects on an LS89 high-pressure turbine (HPT) blade at engine-relevant conditions. Different amplitudes and distributions of surface roughness are investigated and a reference smooth-blade simulation under the same flow conditions is conducted for comparison. Roughness of increasing amplitude progressively shifts the blade suction side boundary layer transition upstream, producing larger values of the turbulent kinetic energy and higher total wake losses. The on-surface data-capturing capabilities of the numerical framework provide direct measurements of the heat flux and the skin friction coefficient, hence offering quantitative information between the surface topology and engineering-relevant performance parameters. This work may provide a benchmark for future numerical studies of turbomachinery flows with roughness.

Engineering↗

Optimizing Aluminum Bonding: Exploring Surface Roughness and Contact Angle Effects through Plasma and Acid Etching Contrasts

The use of adhesives to replace traditional techniques like welding and riveting has been on the rise in various industries as increasingly complex geometries and situations for bonding is required. It’s well known in the scientific community that an adhesive bond is only as efficient as the surface preparation used beforehand. Plasma treatments have been shown to gently clean the surface by only removing the top-most layer from the adherend and “charging” the surface instead of berating it. Whereas chemical treatments “etch” the surface in hopes of increasing surface area by increasing surface roughness. This method can be difficult to control quantitatively and can generate waste leading to serious environmental concerns. This study aims to address the encountered differences between plasma treatments and acid etching on aluminum substrates for adhesive bonding. Treating the surface of adherends plays a crucial role in promoting a good bond between adhesive and adherend. In this study, the effects of plasma and acid treatments are categorized using surface roughness measurements as well as contact angle which have been shown to correlate to mechanical strength. When evaluating different surface treatments, it is notable to mention that not only surface roughness and contact angle play a role in determining which method is most effective but also the consequences of the treatment itself. It was found that although chemical etching was the most effective in increasing roughness and decreasing contact angle, it generates a large amount of waste and is not sustainable compared to plasma work which uses a renewable process. The results from this work provide a deeper understanding of the relationship between surface treatments and surface roughness/contact angle and how this relates to a stronger adhesive bond.

36 MATERIALS SCIENCE↗

Optimizing Aluminum Bonding: Exploring Surface Roughness and Contact Angle Effects through Plasma and Acid Etching Contrasts.

The use of adhesives to replace traditional techniques like welding and riveting has been on the rise in various industries as increasingly complex geometries and situations for bonding is required. It’s well known in the scientific community that an adhesive bond is only as efficient as the surface preparation used beforehand. Plasma treatments have been shown to gently clean the surface by only removing the top-most layer from the adherend and “charging” the surface instead of berating it. Whereas chemical treatments “etch” the surface in hopes of increasing surface area by increasing surface roughness. This method can be difficult to control quantitatively and can generate waste leading to serious environmental concerns. This study aims to address the encountered differences between plasma treatments and acid etching on aluminum substrates for adhesive bonding. Treating the surface of adherends plays a crucial role in promoting a good bond between adhesive and adherend. In this study, the effects of plasma and acid treatments are categorized using surface roughness measurements as well as contact angle which have been shown to correlate to mechanical strength. When evaluating different surface treatments, it is notable to mention that not only surface roughness and contact angle play a role in determining which method is most effective but also the consequences of the treatment itself. It was found that although chemical etching was the most effective in increasing roughness and decreasing contact angle, it generates a large amount of waste and is not sustainable compared to plasma work which uses a renewable process. The results from this work provide a deeper understanding of the relationship between surface treatments and surface roughness/contact angle and how this relates to a stronger adhesive bond

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SURFACE ROUGHNESS REDUCTION AND PERFORMANCE OF Nb3Sn FILM FOR SRF APPLICATIONS

Nb3Sn offers better RF performance (Q and Eacc) than Nb at any given temperature because of its superior superconducting properties. The potential deployment of Nb3Sn-coated SRF cavities at 4 K, delivering similar performance to that of Nb cavities at 2 K, will be a transformational technology to enable new classes of SRF accelerator applications. A clean and smooth surface can enhance the performance of the Nb3Sn-coated cavity, typically, the attainable acceleration gradient. The reduction of surface roughness is often linked with a correlative reduction in average coating thickness and grain size. A new approach to enhance the nucleation centers using the pre-deposition of Sn on Nb is introduced for surface roughness reduction. Besides Sn supply's careful tuning, the temperature profiles were varied to engineer the coating surface with reduced roughness as low as ~40 nm in 20 µm × 20 µm AFM scans, one-third that of the typical coating. Electropolishing of coated sample shows surface a clear surface smoothening. A few sets of coating parameters were used to coat 1.3 GHz sin-gle-cell cavities to understand the effects of roughness variation on the RF performance. This presentation discusses ways to reduce surface roughness with results from a systematic analysis of the samples and Nb3Sn-coated single-cell cavities.

Pudasaini, Uttar↗

SURFACE ROUGHNESS REDUCTION AND PERFORMANCE OF Nb3Sn FILM FOR SRF APPLICATIONS

Nb3Sn offers better RF performance (Q and Eacc) than Nb at any given temperature because of its superior superconducting properties. The potential deployment of Nb3Sn-coated SRF cavities at 4 K, delivering similar performance to that of Nb cavities at 2 K, will be a transformational technology to enable new classes of SRF accelerator applications. A clean and smooth surface can enhance the performance of the Nb3Sn-coated cavity, typically, the attainable acceleration gradient. The reduction of surface roughness is often linked with a correlative reduction in average coating thickness and grain size. A new approach to enhance the nucleation centers using the pre-deposition of Sn on Nb is introduced for surface roughness reduction. Besides Sn supply's careful tuning, the temperature profiles were varied to engineer the coating surface with reduced roughness as low as ~40 nm in 20 µm × 20 µm AFM scans, one-third that of the typical coating. Electropolishing of coated sample shows surface a clear surface smoothening. A few sets of coating parameters were used to coat 1.3 GHz sin-gle-cell cavities to understand the effects of roughness variation on the RF performance. This presentation discusses ways to reduce surface roughness with results from a systematic analysis of the samples and Nb3Sn-coated single-cell cavities.

Pudasaini, Uttar↗

Robustness of Optical Response for Self-Assembled Plasmonic Metamaterials with Morphological Disorder and Surface Roughness

Bottom-up fabrication of metallized biotemplated nanostructures to form specific plasmonic nanoresonators holds promise as a means of achieving large-scale optical metamaterials. However, in contrast to top-down methods, the stochastic growth of self-assembled nanoresonators is prone to significant disorder and surface roughness, which naturally raise an important question about the robustness of their resonant properties in terms of structural imperfections. An aggregated-random-sphere model is developed to mimic the nucleated growth of metallized DNA origami assembly, leading to meta-atoms with realistic, experimentally observed morphological disorder and surface roughness. Using the well-known split-ring-resonator (SRR) motif as an example, the resonant properties of meta-atoms under different levels of roughness are investigated and a strong tolerance of optical response against morphological disorder is revealed. It is found that in SRRs, even with dramatic roughness introduced, the expected resonances are still observed, despite broadening line shapes compared to ideal smooth structure. Only for extreme disorder, which causes drastic segmentation of SRRs, does the resonant response disappear. Finally, the demonstrations are very encouraging for the prospects of bottom-up fabrication toward versatile functional metamaterials and metadevices.

36 MATERIALS SCIENCE↗

Load-Displacement Relation and Gap Distribution Between Rough Surfaces: Partial Differential Equations Approach

We develop a theoretical model to predict the load-displacement relation and probability density function for gaps between contacting rough surfaces. We derive partial differential equations from the previous model (Joe et al., 2018), and extend them to the non-adhesive contact problem. The predictions of the present theory are compared with numerical results using the Green's function molecular dynamics algorithm, and a good agreement is obtained.

adhesion↗

Interfacial toughness: dependence on surface roughness and test temperature

Interfacial toughness quantifies resistance to crack growth along an interface and in this investigation the toughness of an aluminum/epoxy interface was measured as a function of surface roughness and test temperature. The large strain response of the relatively ductile epoxy adhesive used in this study was also characterized. This epoxy adhesive exhibits intrinsic strain-softening after initial compressive yield and then deforms plastically at a roughly constant flow stress until it rapidly hardens at large compressive strains. Here, we find that interface toughness scales as the product of the temperature dependent epoxy yield strength and a length scale that characterizes surface roughness. The proposed scaling is based upon dimensional considerations of a model problem that assumes that the characteristic length scale of both the roughness and the crack-tip yield zone is small relative to the region dominated by the linear elastic asymptotic crack-tip stress field. Furthermore, the model assumes that interfacial failure occurs only after the epoxy begins to harden at large strains. The proposed relationship is validated by our interfacial toughness measurements.

36 MATERIALS SCIENCE↗

Extreme sensitivity of higher-order interatomic force constants and thermal conductivity to the energy surface roughness of exchange-correlation functionals

In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. In conclusion, with these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.

42 ENGINEERING↗