skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Phonon Boltzmann Transport Equation Based Simulation of Frequency Domain Thermoreflectance Experiments
Abstract The Boltzmann Transport equation (BTE) was solved numerically in cylindrical coordinates and in time domain to simulate a Frequency Domain Thermo-Reflectance (FDTR) experiment. First, a parallel phonon BTE solver that accounts for all phonon modes, frequencies, and polarizations was developed and tested. The solver employs the finite-volume method (FVM) for discretization of physical space, and the finite-angle method (FAM) for discretization of angular space. The solution was advanced in time using explicit time marching. The simulations were carried out in time domain and band-based parallelization of the BTE solver was implemented. The phase lag between the temperature averaged over the probed region of the transducer and the modulated laser pump signal was extracted for a pump laser modulation frequency ranging from 20–200 MHz. It was found that with the relaxation time scales used in the present study, the computed phase lag is underpredicted when compared to experimental data, especially at smaller modulation frequencies. The challenges in solving the BTE for such applications are highlighted.  more » « less
Award ID(s):
2003747
PAR ID:
10442310
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Proceedings of the International Mechanical Engineering Congress and Exposition, IMECE2022
Volume:
V008T11A062
Page Range / eLocation ID:
95630
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, International Journal of Heat and Mass Transfer)
    The anisotropic Fourier Heat Conduction Equation (FHCE) and the multidimensional phonon Boltzmann Transport Equation (BTE) were solved numerically in cylindrical coordinates and in time domain to simulate a Time Domain Thermo-Reflectance (TDTR) experimental silicon/aluminum substrate/transducer setup. The out-of-phase response of the probe laser was predicted at various beam offset distances for a pump laser pulse frequency of 80 MHz and modulation frequency of 10 MHz and compared against experimental measurements for a silicon substrate. The isotropic FHCE was also solved for comparison. Results show that the isotropic FHCE with bulk thermal conductivity of 145 W/m/K significantly underpredicts the out-of-phase temperature difference, particularly at smaller beam offsets. With an isotropic thermal conductivity of 105 W/m/K, the computed results match experimental data at smaller beam offsets well, but overpredicts the experimental data at larger beam offsets. An almost-perfect match is obtained by using an anisotropic thermal conductivity wherein the radial (in-plane) thermal conductivity is set to 85 W/m/K and the axial (through-plane) conductivity is set to 130 W/m/K. The multidimensional frequency and polarization dependent phonon BTE is next solved. The BTE results for the out-of-phase temperature difference match experimental observations well at small and intermediate beam offsets, but overpredicts the experimental data at larger beam offsets. FHCE results are fitted to the BTE predictions, and the extracted (best fit) thermal conductivity is found to be 110 W/m/K. 
    more » « less
  2. ; (, Journal of Heat Transfer)
    Abstract The Fourier and the hyperbolic heat conduction equations were solved numerically to simulate a frequency-domain thermoreflectance (FDTR) experiment. Numerical solutions enable isolation of pump and probe laser spot size effects and use of realistic boundary conditions. The equations were solved in time domain and the phase lag between the temperature of the transducer (averaged over the probe laser spot) and the modulated pump laser signal was computed for a modulation frequency range of 200 kHz–200 MHz. Numerical calculations showed that extracted values of the thermal conductivity are sensitive to both the pump and probe laser spot sizes, while analytical solutions (based on Hankel transform) cannot isolate the two effects. However, for the same effective (combined) spot size, the two solutions are found to be in excellent agreement. If the substrate (computational domain) is sufficiently large, the far-field boundary conditions were found to have no effect on the computed phase lag. The interface conductance between the transducer and the substrate was found to have some effect on the extracted thermal conductivity. The hyperbolic heat conduction equation yielded almost the same results as the Fourier heat conduction equation for the particular case studied. The numerically extracted thermal conductivity value (best fit) for the silicon substrate considered in this study was found to be about 82–108 W/m/K, depending on the pump and probe laser spot sizes used. 
    more » « less
  3. ; (, Proceedings of the ASME Summer Heat Transfer Conference, SHTC2020)
    Abstract The Fourier heat conduction and the hyperbolic heat conduction equations were solved numerically to simulate a frequency-domain thermoreflectance (FDTR) experimental setup. Numerical solutions enable use of realistic boundary conditions, such as convective cooling from the various surfaces of the substrate and transducer. The equations were solved in time domain and the phase lag between the temperature at the center of the transducer and the modulated pump laser signal were computed for a modulation frequency range of 200 kHz to 200 MHz. It was found that the numerical predictions fit the experimentally measured phase lag better than analytical frequency-domain solutions of the Fourier heat equation based on Hankel transforms. The effects of boundary conditions were investigated and it was found that if the substrate (computational domain) is sufficiently large, the far-field boundary conditions have no effect on the computed phase lag. The interface conductance between the transducer and the substrate was also treated as a parameter, and was found to have some effect on the predicted thermal conductivity, but only in certain regimes. The hyperbolic heat conduction equation yielded identical results as the Fourier heat conduction equation for the particular case studied. The thermal conductivity value (best fit) for the silicon substrate considered in this study was found to be 108 W/m/K, which is slightly different from previously reported values for the same experimental data. 
    more » « less
  4. ; ; (, Journal of Computational Electronics)
    We calculate the electrical conductivity of suspended and supported monolayer MoS2 at terahertz (THz) frequencies by means of EMC–FDTD, a multiphysics simulation tool combining an ensemble Monte Carlo (EMC) solver for electron transport and a finite-difference time-domain (FDTD) solver for full-wave electrodynamics. We investigate the role of carrier and impurity densities, as well as substrate choice (SiO2 or hexagonal boron nitride, hBN), in frequency-dependent electronic transport. Owing to the dominance of surface-optical-phonon scattering, MoS2 on SiO2 has the lowest static conductivity, but also the weakest overall frequency dependence of the conductivity. In fact, at high THz frequencies, the conductivity of MoS2 on SiO2 exceeds that of either suspended or hBN-supported MoS2. We extract the parameters for Drude-model fits to the conductivity versus frequency curves obtained from microscopic simulation, which may aid in the experimental efforts toward MoS2 THz applications. 
    more » « less
  5. ; ; (, Optics Letters)
    We investigate the impact of the finite rise time of a spatiotemporal boundary inside a dispersive medium used for reflection and refraction of optical pulses. We develop a matrix approach in the frequency domain for analyzing such spatiotemporal boundaries and use it to show that the frequency range over which reflection can occur is reduced as the rise time increases. We also show that total internal reflection can occur even for boundaries with long rise times. This feature suggests that spatiotemporal waveguides can be realized through cross-phase modulation even when pump pulses have relatively long rise and fall times. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email