Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion length”

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 1,261 records · Page 70

Buoyancy Effects in Fully-Modulated, Turbulent Diffusion Flames

Pulsed combustion appears to have the potential to provide for rapid fuel/air mixing, compact and economical combustors, and reduced exhaust emissions. The objective of this experiment (PuFF, for Pulsed-Fully Flames) is to increase the fundamental understanding of the fuel/air mixing and combustion behavior of pulsed, turbulent diffusion flames by conducting experiments in microgravity. In this research the fuel jet is fully-modulated (i.e., completely shut off between pulses) by an externally controlled valve system. This gives rise to drastic modification of the combustion and flow characteristics of flames, leading to enhanced fuel/air mixing compared to acoustically excited or partially-modulated jets. Normal-gravity experiments suggest that the fully-modulated technique also has the potential for producing turbulent jet flames significantly more compact than steady flames with no increase in exhaust emissions. The technique also simplifies the combustion process by avoiding the acoustic forcing generally present in pulsed combustors. Fundamental issues addressed in this experiment include the impact of buoyancy on the structure and flame length, temperatures, radiation, and emissions of fully-modulated flames.

Hermanson, J. C.↗

Surface Processes of Faceted Growth

Atomic force microscopy and high precision optical interferometry allow to analyze the processes that are in the core of our present understanding of faceted crystal growth. Some of these processes will be reviewed. Recent experiments suggest that the Gibbs-Thomson Law (GTL) may not be applicable to the weakly fluctuating strongly polygonized steps, with low kink density. Propagation rate of a straight short step segment changes with the segment length much steeper than predicted by GTL. If the step fluctuations are still well developed, the typical transition length may be determined just by the radius of rounded comers of the polygonized step. If fluctuations do not have enough time to develop, the kinetic effects may be essential. Indeed, 'communications' between comers of a short straight step segment may occur and, thus, the GTL may be implemented only via kink exchange. E.g., dissolutions of this short segment at macroscopic equilibrium occurs as follows: The kinks are split from one comer as a result of fluctuations, and annihilate with kinks of the opposite sign generated by another comer. Similarly, an exchange via chain of kinks on a rough step is the mechanism to implement the GT Law. If there is a supersaturation with respect to a large crystal, the kinks not only diffuse along the step, but move back the generating comers. This 'wind' prevents communication between the comers and places upper limit supersaturation only below which GTL is applicable. Steps replace kinks in the 3D case, though, of course, another dimensionality brings about another physics. Another specific problem to think about is if fluctuations are always fast enough to provide enough kinks for linear dependence of step rate on supersaturation. Increasing azimuthal polygonization of growth hillocks is a sign of insufficient fluctuation rate. Morphological stability of vicinal faces is usually considered on the basis of mutually parallel steps. Azimuthal anisotropy and interlacing of step bunches inducing the well known bunch splitting is another challenge. Among other challenging problems are: existence of surface diffusion on the crystal/so1ution interface, step interaction in solutions, achievement of step flow modes from liquids, impurity effects. New phenomena essential in biomacromolecular crystallization is another interesting area.

Chernov, A.↗

Characterization of a Laminate Flat Plate Diffusion Flame in Microgravity using PIV, Visible and CH Emissions

Motivated by fire safety concerns and the advent of long-term micro-gravity facilities, a cooperative program has been developed to study the mechanisms and material properties that control flow assisted (co-current) flame spread. This program has used as a common fire scenario a reacting steady-state boundary layer. Preliminary studies explored the aerodynamics of a reacting boundary layer by simulating a condensed fuel by means of a gas burner. Stability curves for ethane air flames were obtained and different burning regimes were identified. An important feature of this study was the independent identification of the different mechanisms leading to the instability of the flow. It was observed that fuel injection velocity and thermal expansion independently contributed to the separation of the flow at the leading edge of the burner. The occurrence of separation resulted in complex three-dimensional flow patterns that have a dominant effect on critical fire safety parameters such as the stand-off distance and flame length. This work was extended to a solid fuel (PMMA) leading to a Sounding Rocket experiment (Mini-Texus-6). The solid phase showed similar flow patterns, mostly present at low flow velocities (<100 mm/s) but the results clearly demonstrated that the thermal balance at the pyrolyzing fuel surface is the dominant mechanism that controls both stand-off distance and flame length. This thermal balance could be described in a global manner by means of a total mass transfer or "B" number. This "B" number incorporates surface re-radiation, radiative feedback and in-depth heat conduction as first prescribed by Emmons. The mass transfer number becomes the single parameter that determines the evolution of these fire safety variables (flame length, stand-off distance) and therefore can be used as a ranking criterion to assess the flammability of materials. The particular configuration is representative of the NASA upward flame spread test (Test 1) therefore this approach can be used in the interpretation of the results obtained from this test. Nevertheless, complete validation of this approach has not been fully achieved due, mainly because all the measurements necessary to compare with the theoretical predictions have not been obtained. Following these studies two different directions have been taken. The first attempts to elucidate the details of the gas phase combustion reaction and the associated flow field by means of quantitative and qualitative measurements. The second approach, a more practical one, is to apply this methodology to the assessment of material flammability. The former is currently being conducted with a gas burner because it allows for easier control and longer experimentation time. The results obtained so far will be presented in more detail. The latter is a new program therefore only a brief summary of the objectives will be presented.

Joulain, P.↗

Using Strain in the Vadose Zone During Sinusoidal-Rate Well Tests to Characterize an Aquifer

Well tests with sinusoidal variation in pumping rate are appealing because of the ability to filter competing signals and they also can be set up with no net removal of water, which is an advantage at contaminated sites. Current methods of interpreting sinusoidal-rate well tests analyze the phase lag and amplitude of pressure signals in the aquifer. This approach is restricted by the availability of monitoring wells. An alternative approach could be to measure the stain in the vadose zone instead of pressure in the aquifer. The goal of this research is to evaluate the feasibility of using signals from in-situ sensors in the vadose zone to characterize sinusoidal-rate well tests. Well tests with a variety of periods and amplitudes were performed using a well completed in an unconfined saprolite aquifer with a water table at 9m depth at our test site near Clemson, South Carolina. Vertical strain was measured every meter between 1 to 7 m using a distributed optical fiber sensing technique called coherence length gated microwave photonics interferometry (CMPI). A novel packaging system was designed to couple the optical fiber to the saprolite. The results show that strain measured in the vadose zone during sinusoidal-rate pumping varies periodically with amplitudes from 0.1 to 3 microstrain, depending on the amplitude of the sinusoidal pumping rate. Hydraulic heads at similar radial distances were also periodic with changing amplitudes ranging up to ±0.2 meters. The strain in the vadose zone responds to pressure in the underlying aquifer—the strain is tensile when the pressure is falling, and compressive when the pressure is rising. The phase lag between the pumping rate and the strain was 3 min with a pumping period of 15 minutes, and this gives a hydraulic diffusivity of 0.06 m2/s. This is the same diffusivity obtained by analyzing a constant rate pumping test, and it is slightly larger than the diffusivity (0.044 m2/s) obtained by analyzing the phase lag of the pressure in a nearby monitoring well. Strain measurements in the vadose zone could be made with much less cost and effort than installing a deep monitoring well. These results suggest that strain measurements in the vadose zone could be used to supplement monitoring wells, thereby increasing the resolution, and reducing costs of analyzing well tests.

well test, sinusoidal, periodic pumping, strain↗

Signatures of Thermoreversible Associations in X-ray and Neutron Scattering from Dilute Polyzwitterion Solutions

In aqueous solutions of polyzwitterions (PZs), an interplay between dipole–dipole interactions and hydration of zwitterionic groups can lead to thermoreversible associations, which have been difficult to detect in experiments. Here, in this study, we investigated dilute aqueous solutions of poly(1-(3-sulphonatopropyl)-2-vinylpyridinium) (P2VPPS) using small-angle X-ray and neutron scattering (SAXS and SANS) to probe the structure and neutron spin-echo (NSE) spectroscopy to probe dynamics. The SAXS and SANS data show that the correlation length increases with an increase in the concentration of P2VPPS for three different molecular weights. The addition of 0.1 M NaCl to one of the solutions led to almost no dependence of the correlation length on the concentration. Such a concentration dependence of the correlation length suggests the formation of clusters driven by thermoreversible associations in the solutions. The NSE measurements show that the solutions with a larger correlation length display slower relaxation, reflecting the reduced mobility of larger clusters. These results should be considered as signatures of thermoreversible associations in the solutions of P2VPPS. To establish a quantitative link between local structure (clusters) and dynamics in dilute solutions of P2VPPS, we combined a thermoreversible gelation theory for the structure of PZ solutions (Li, S.-F.; Muthukumar, M. Theory of Thermoreversible Gelation and Anomalous Concentration Fluctuations in Polyzwitterion Solutions. J. Chem. Phys. 2024, 161, 024903) with a model for the dynamics of the clusters (generalized Zimm model), developed in this work. Using such a theoretical framework, we have predicted the distribution of clusters in the solutions probed with SANS and SAXS. With the distributions, the generalized Zimm model has been used to extract the diffusion constant of the clusters and their characteristic size from the NSE data, where the latter agrees with the values estimated from the SAXS/SANS data. These findings confirm the presence of thermoreversible associations and establish a quantitative link between local structure (clusters) and dynamics in dilute solutions of P2VPPS. Furthermore, with growing interest in technological applications, this work can provide useful insights into the structural and dynamical properties of other PZ solutions.

field theory↗

One-dimensional numerical study of charged particle trajectories in turbulent electrostatic wave fields

The paper describes a numerical simulation of electron trajectories in weak random electric fields under conditions that are approximately true for Langmuir waves whose wavelength is much longer than the Debye length. Two types of trajectory calculations were made: (1) the initial particle velocity was made equal to the mean phase velocity of the waves, or (2) it was equal to 0.7419 times the mean velocity of the waves, so that the initial velocity differed substantially from all phase velocities of the wave spectrum. When the autocorrelation time is much greater than the trapping time, the particle motion can change virtually instantaneously from one of three states - high-velocity, low-velocity, or trapped state - to another. The probability of instantaneous transition from a high- or low-velocity state becomes small when the difference between the particle velocity and the mean phase velocity of the waves becomes high in comparison to the trapping velocity. Diffusive motion becomes negligible under these conditions also.

Graham, K. N.↗

Absence of magnetic trapping on closed field lines at Neptune

The large offset of Neptune's magnetic field in an offset tilted dipole (OTD) approximation allows some fraction of magnetospheric field lines to have monotonically increasing field magnitude over their entire length in the magnetosphere. The usual magnetic trapping is impossible on such field lines, meaning that all charged particles on these field lines will precipitate after mirroring at most once. These field lines can be thought of as the part of the global magnetic anomaly where no particles on the corresponding drift shells remain trapped for an entire rotation of the planet. In this paper we examine both the morphology and size of the affected surface area in the OTD2 and O8 magnetic field models. For regions where there are indications that part of the aurora could be due to monotonic field line precipitation we analyze whether these field lines could be sufficiently populated to produce auroral signatures by considering photoelectron production rates and radial diffusion. We conclude that monotonic field line precipitation may contribute to the observed Neptune aurora but does not explain the global phenomenon.

Paranicas, C.↗

Vortex Generators in a Streamline-Traced, External-Compression Supersonic Inlet

Vortex generators within a streamline-traced, external-compression supersonic inlet for Mach 1.66 were investigated to determine their ability to increase total pressure recovery and reduce total pressure distortion. The vortex generators studied were rectangular vanes arranged in counter-rotating and co-rotating arrays. The vane geometric factors of interest included height, length, spacing, angle-of-incidence, and positions upstream and downstream of the inlet terminal shock. The flow through the inlet was simulated numerically through the solution of the steady-state, Reynolds-averaged Navier-Stokes equations on multi-block, structured grids using the Wind-US flow solver. The vanes were simulated using a vortex generator model. The inlet performance was characterized by the inlet total pressure recovery and the radial and circumferential total pressure distortion indices at the engine face. Design of experiments and statistical analysis methods were applied to quantify the effect of the geometric factors of the vanes and search for optimal vane arrays. Co-rotating vane arrays with negative angles-of-incidence positioned on the supersonic diffuser were effective in sweeping low-momentum flow from the top toward the sides of the subsonic diffuser. This distributed the low-momentum flow more evenly about the circumference of the subsonic diffuser and reduced distortion. Co-rotating vane arrays with negative angles-of-incidence or counter-rotating vane arrays positioned downstream of the terminal shock were effective in mixing higher-momentum flow with lower-momentum flow to increase recovery and decrease distortion. A strategy of combining a co-rotating vane array on the supersonic diffuser with a counter-rotating vane array on the subsonic diffuser was effective in increasing recovery and reducing distortion.

computational fluid dynamics↗

An Encounter with the Ion and Electron Diffusion Regions at a Flapping and Twisted Tail Current Sheet

We analyze data returned by the Magnetospheric Multiscale mission (MMS) constellation during a rapid (~1.5 s) traversal of a flapping and reconnecting current sheet (CS) in the near-Earth magnetotail (X ~ -20 RE). The CS was highly tilted, with its normal pointing strongly duskward. Its extreme thinness was confirmed by a curvature analysis of the magnetic field lines. The event was associated with a guide field of 8% of the reconnecting components. From the pitch angle distributions of low-energy electrons we infer that the crossing occurred earthward of the X-line. Traveling practically normal to the CS, MMS encountered an ion diffusion region (IDR) in which was embedded an electron diffusion region (EDR). IDR signatures included breaking of the ion frozen-in condition in the presence of Hall B and E fields. EDR signatures included a strong out of-plane current associated with a super-Alfvenic electron jet, evidence of positive energy transfer, and a temperature anisotropy (TeII > Te1) which disappeared at the field reversal. Derived scale sizes normal to the CS are: ~6.9 de (EDR) and ~0.4 di (IDR) (40 and 100 km). We estimate the average dimensionless reconnection rate as 0.077 + 0.050. The observations and inferences are supported by PIC numerical simulations. We find very good agreement in the reconnection rates. We also discuss the effects of asymmetries in the density, temperature and magnetic field strength on the Hall fields and length of the outflow jets. The event is associated with a substorm onset which began 7 min after the MMS observations.

C J Farrugia↗

Characterisation of nano-assemblies inside mesopores using neutron scattering

Adsorption of molecular and nanoscale matter in mesoporous materials is important in filtration and chromatography as well as membrane processes. However, the morphology and distribution of self-assembled structures of adsorbate formed within nanoconfined environments is largely unknown. Small Angle Neutron Scattering (SANS) with porous matrix matching the neutron scattering length density of the solvent has the potential to provide detailed information on the self-assemblies formed in pore spaces. However, analysis of such SANS profiles remains a challenge. In this study, we extend the SANS analysis method previously developed by Findenegg and co-workers to include interparticle correlations, providing a qualitative characterisation of the adsorption state of surfactants and nanoparticles in the cylindrical pores of silica nanomaterials. We find that the pore filling fraction and the self-assembled state of the adsorbate governs the scattering profile. We apply the model to two materials namely, triethyleneglycol monohexyl ether (C 6 E 3 ) surfactant and gold nanoparticles adsorbed in SBA-15 mesoporous silica. In contrast to the C 6 E 3 system, Bragg scattering dominates the diffuse scattering in gold nanoparticle system, which we attribute to the immobile internal structure of the nanoparticles. In conclusion, this new SANS analysis method has potential to elucidate structures of soft matter nanoassemblies inside mesopores.

74 ATOMIC AND MOLECULAR PHYSICS↗

Computer simulation of a wind tunnel test section with discrete finite-length wall slots

A computer simulation of a slotted wind tunnel test section which includes a discrete, finite-length wall slot representation with plenum chamber constraints and accounts for the nonlinear effects of the dynamic pressure of the slot outflow jet and of the low energy of slot inflow air was developed. The simulation features were selected to be those appropriate for the intended subsequent use of the simulation in a wall interference assessment procedure using sparsely located wall pressure measurements. It is demonstrated that accounting for slot discreteness is important in interpreting wall pressure measured between slots, and that accounting for nonlinear slot flow effects produces significant changes in tunnel-induced velocity distributions and, in particular, produces a longitudinal component of tunnel-induced velocity due to model lift. A characteristic mode of tunnel flow interaction with constraints imposed by the plenum chamber and diffuser entrance is apparent in simulation results and is derived analytically through a simplified analysis.

Kemp, W. B., Jr.↗

A Microwave Photonics Optical Fiber Method for Measuring Distributed Strain for Hydrologic Applications in the Vadose and Saturated Zones

Distributed strain measurements appear to hold significant potential for monitoring hydrologic processes. Coherence-length-gated Microwave Photonics Interferometry (CMPI) is a distributed sensing technique that measures strain in an optical fiber by reading optical interference phase changes in the microwave domain. The technique provides 10nε resolution when cm spatial resolution is applied. CMPI was used to measure the strain caused by small periodic variations in air pressure (25 Pa amplitude and 4 Hz) in a sand-filled laboratory column used to represent barometric loading in the vadose zone. The strain varied as a periodic function with an amplitude that decreased and a phase that increased with depth. The distribution of amplitude and phase of the strain depended on the water saturation, permeability, and presence of a barrier at the top of the sand. A theoretical analysis suggests that the air pressure diffusivity estimated from pressure data is similar to the diffusivity estimated from the distributed strain. These results indicated that the CMPI distributed strain measurement system could be used to improve monitoring and characterization of the vadose and underlying saturated zone.

Hua, Liwei↗

Evaluation of light collection from highly scattering media using wavelength-shifting fibers

Opaque scintillators are designed to have a short scattering length such that scintillation photons are probabilistically confined to a small region of space about their origin. The benefit of this feature is that information on the interaction event topology can be recorded with greater fidelity than in traditional highly transparent media with sensors at large distances from the light production region. Opaque scintillator detectors rely on wavelength-shifting fibers to extract the scintillation light; however, the efficiency of light collection has not yet been directly measured in experiment. We measured the efficiency of light collection as a function of the optical parameters of an opaque liquid and the distance from the origin of the light to the fiber. We use the experimental data to validate a Monte Carlo model of light transport and collection and discuss a simple diffusion model that reproduces the results of Monte Carlo simulation with high fidelity. Furthermore, this combination of validated models has the potential for use in predictions of performance in various designs of future opaque scintillator detectors such as LiquidO.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Convergence properties of finite-difference hydrodynamics schemes in the presence of shocks

We investigate the asymptotic convergence of finite-difference schemes for the Euler equations when the limiting solution contains shocks. The Lax-Wendroff theorem guarantees that certain conservative schemes converge to correct, physically valid solutions. We focus on two one-dimensional operator-split schemes with explicit artificial-viscosity terms. One, an internal-energy scheme, does not satisfy the assumptions of Lax-Wendroff; the other, a conservative total-energy scheme, does. With viscous lengths chosen proportional to the grid size, we find that both schemes converge to their zero-grid-size limits at the theoretically expected rate, but only the conversative scheme converges toward correct solutions of the inviscid fluid equations. We show that the difference in their behaviors results directly from the presence of shocks in the limiting solution. Empirically, we find that when the viscous lenghts tend toward zero more slowly than the grid size, however the nonconservative scheme also converges toward correct solutions. We characterize the asymptotic behavior of the total-energy scheme in a particular problem in which a shock forms. As the grid is refined, a Cauchy error approaches the expected rate of change slowly. We show that the changes in the artificial viscosity alter the diffusion of small-amplitude waves. The differences associated with such waves make the dominant contribution to the Cauchy error. We formulate an analytic model to relate the rate of approach to the effect of varying diffusion in waves and find quantitative agreement with our numerical results.

Kimoto, Paul A.↗

Carbon-Nanotube-Carpet Heat-Transfer Pads

Microscopic thermal-contact pads that include carpet-like arrays of carbon nanotubes have been invented for dissipating heat generated in integrated circuits and similarly sized single electronic components. The need for these or other innovative thermal-contact pads arises because the requisite high thermal conductances cannot be realized by scaling conventional macroscopic thermal-contact pads down to microscopic sizes. Overcoming limitations of conventional thermal-contact materials and components, the carbon-nanotube thermal-contact pads offer the high thermal conductivities needed to accommodate the high local thermal power densities of modern electronic circuits, without need for large clamping pressures, extreme smoothness of surfaces in contact, or gap-filling materials (e.g., thermally conductive greases) to ensure adequate thermal contact. Moreover, unlike some conventional thermal-contact components, these pads are reusable. The figure depicts a typical pad according to the invention, in contact with a rough surface on an electronic component that is to be cooled. Through reversible bending and buckling of carbon nanotubes at asperities on the rough surface, the pad yields sufficiently, under relatively low contact pressure, that thermal contact is distributed to many locations on the surface to be cooled, including valleys where contact would not ordinarily occur in conventional clamping of rigid surfaces. Hence, the effective thermal-contact area is greater than that achievable through scaling down of a macroscopic thermal-contact pad. The extremely high longitudinal thermal conductivities of the carbon nanotubes are utilized to conduct heat away from potential hot spots on the surface to be cooled. The fibers protrude from a layer of a filler material (Cu, Ag, Au, or metal-particle- filled gels), which provides both mechanical support to maintain the carbon nanotubes in alignment and thermal conductivity to enhance the diffusion of concentrated heat from the nanotubes into the larger adjacent volume of a heat sink. The array of carbon nanotubes, the filler material, and the heat sink are parts of a unitary composite structure that is fabricated as follows: 1. Using techniques that have been reported previously, the array of substantially perpendicularly oriented carbon nanotubes is grown on a metal, silicon, or other suitable thermally conductive substrate that is intended to become the heat sink. 2. By means of chemical vapor deposition, physical vapor deposition, plasma deposition, ion sputtering, electrochemical deposition, or casting from a liquid phase, some or all of the interstitial volume between carbon nanotubes is filled with the aforementioned layer of mechanically supporting, thermally conductive material. 3. To cause the carbon nanotubes to protrude the desired length from the filler material, an outer layer of filler is removed by mechanical polishing, chemical mechanical polishing, wet chemical etching, electrochemical etching, or dry plasma etching.

Li, Jun↗

A mechanistic interpretation of Nelson curves for PVP failures under high temperature hydrogen attack

As an empirically established design criterion, Nelson curves that relate the service temperature and the allowable hydrogen partial pressure have been developed and utilized for more than sixty years in pressure vessels and piping (PVP) safety design. Despite a relatively clear thermodynamic understanding of the high-temperature-hydrogen-attack (HTHA) problem, the detailed fracture process on the microstructural length scales, however, remains elusive, and a quantitative assessment of the PVP lifetime under HTHA from the available creep fracture dataset is still not possible. This work develops a microstructure-informed and micromechanics-based model by incorporating a synergy between hydrogen transport and intergranular-cavity-based fracture process. Based on the available creep lifetime data of C-0.5Mo steels, we are able to calibrate material constitutive parameters, and then conduct nonlinear finite element simulations that reveal a real-time stress-induced hydrogen diffusional transport along grain boundaries, coupled with a microstructure-explicit failure process, from which Nelson curves can be computed. Such failure analyses allow us to delineate two distinct regimes on the Nelson curves, i.e., dislocation-creep-controlled or grain boundary diffusion-assisted cavity growth. More importantly, we found that a small change of the pipe thickness and applied stresses can significantly shift these lifetime curves. However, these two parameters are usually not provided in Nelson curves, thus limiting their usage in material selection and safety design. In conclusion, this discrepancy can clearly be mitigated by extensive parametric studies from our micromechanical modeling/simulation framework.

36 MATERIALS SCIENCE↗

Generation of Z mode radiation by diffuse auroral electron precipitation

The generation of Z mode waves by diffuse auroral electron precipitation is investigated assuming that a loss cone exists in the upgoing portion of the distribution due to electron interactions with the atmosphere. The waves are generated at frequencies above, but very near, the local electron cyclotron frequency omega(e) and at wave normal angles larger than 90 deg. In agreement with Hewitt et al. (1983), the group velocity is directed downward in regions where the ratio of the upper hybrid frequency omega(pe) to Omega(e) is less than 0.5, so that Z mode waves excited above a satellite propagate toward it and away from the upper hybrid resonance. Z mode waves are excited in a frequency band between Omega(e) and about 1.02 Omega(e), and with maximum growth rates of about 0.001 Omega(e). The amplification length is about 100 km, which allows Z mode waves to grow to the intensities observed by high-altitude satellites.

Dusenbery, P. B.↗

Molecular tail chemistry controls thermal transport in fullerene films

Herein we report on the thermal conductivities of alkyl- and indene-group functionalized fullerene derivative thin films as measured via time domain thermoreflectance and steady state thermoreflectance. The thermal conductivities vary from 0.064 ± 0.007 W m – 1 K – 1 for [6,6]-phenyl C 61 -butyric acid methyl ester (PCBM) to 0.15 ± 0.017 W m – 1 K – 1 for indene- C 60 bisadduct at room temperature and do not exhibit significant temperature dependence from 300 to 375 K. In comparison to the thermal conductivity of PCBM, increasing the length of the alkyl chain, as in the case of [6,6]-phenyl C 61 -butyric acid butyl ester, and [6,6]-phenyl C 61 -butyric acid octyl ester leads to higher thermal conductivities. Likewise, increasing the number of alkyl chains attached to the fullerenes as in the case of bisadduct PCBM leads to a higher thermal conductivity compared to that of PCBM. We present atomistic insights into the role of chemical functionalization on the overall heat transfer in these fullerene derivatives by conducting molecular dynamics simulations and lattice dynamics calculations. The thermal conductivities predicted via our atomistic simulations qualitatively agree with the experimental trends for our fullerene derivatives. Lattice dynamics calculations reveal that one of the main factors dictating the ultralow thermal conductivities in fullerene derivatives is the large reduction in modal diffusivities in the molecular crystals as calculated from the Allen-Feldman model, thus providing an explanation for their largely reduced thermal conductivities as compared to that of C 60 crystals. The low diffusivities result from high degrees of localization of Einstein-like vibrations in fullerene derivatives due to the molecular side chains, providing the ability to dial-in the properties of these low thermal conductivity solids via molecular engineering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗