Solutions to Heat-Conduction Equation with Time-Dependent Boundary Conditions

1961 ◽  
Vol 28 (3) ◽  
pp. 251-252 ◽  
Author(s):  
Arthur E. Bergles ◽  
Joseph Kaye
Author(s):  
Siddharth Saurav ◽  
Sandip Mazumder

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.


1964 ◽  
Vol 60 (4) ◽  
pp. 897-907 ◽  
Author(s):  
M. Wadsworth ◽  
A. Wragg

AbstractThe replacement of the second space derivative by finite differences reduces the simplest form of heat conduction equation to a set of first-order ordinary differential equations. These equations can be solved analytically by utilizing the spectral resolution of the matrix formed by their coefficients. For explicit boundary conditions the solution provides a direct numerical method of solving the original partial differential equation and also gives, as limiting forms, analytical solutions which are equivalent to those obtainable by using the Laplace transform. For linear implicit boundary conditions the solution again provides a direct numerical method of solving the original partial differential equation. The procedure can also be used to give an iterative method of solving non-linear equations. Numerical examples of both the direct and iterative methods are given.


2016 ◽  
Vol 33 (1) ◽  
pp. 65-75 ◽  
Author(s):  
M. R. Talaee ◽  
V. Sarafrazi

AbstractThis paper is devoted to the analytical solution of three-dimensional hyperbolic heat conduction equation in a finite solid medium with rectangular cross-section under time dependent and non-uniform internal heat source. The closed form solution of both Fourier and non-Fourier profiles are introduced with Eigen function expansion method. The solution is applied for simple simulation of absorption of a continues laser in biological tissue. The results show the depth of laser absorption in tissue and considerable difference between the Fourier and Non-Fourier temperature profiles. In addition the solution can be applied as a verification branch for other numerical solutions.


If the solution, of the heat conduction equation θ τ ( 0 ) = θ ξ ξ ( 0 ) , ξ > 0 , τ > 0 of a chemically ‘inert’ material is known, then an approximate formula for the explosion time, ד expl. , of an explosive satisfying the heat conduction equation with zero order reaction, θ ד = θ ξξ +exp(-1/θ), ξ > 0, ד 0, and the same initial and boundary conditions as the ‘inert’, is given by the root of the equation, − ∂ θ ( 0 ) ( ξ , τ expt . ) / ∂ ξ | ξ − 0 = ∫ 0 ∞ exp ⁡ [ − 1 / θ ( 0 ) ( ξ , τ expl . ) ] d ξ provided 1/θ (0) (ξ, ד) is suitably expanded about the surface ξ = 0 such that the integrand vanishes as ξ→∞. Similar results hold for one-dimensional cylindrically and spherically symmetric problems. The derivation of the explosion criterion is based on observation of existing numerical solutions where it is seen that (i) almost to the onset of explosion, the solution θ(ξ, ד )does not differ appreciably from θ (0) (ξ, ד ) (ii) the onset of explosion is indicated by the appearance of a temperature maximum at the surface. Simple formulas for ד expl. readily obtainable for a wide variety of boundary conditions, are given for seven sample problems. Among these are included a semi-infinite explosive with constant surface flux, convective surface heat transfer, and constant surface temperature with and without subsurface melting. The derived values of ד expl. are in satisfactory agreement with those obtained from finite-difference solutions for the problems that can be compared.


Sign in / Sign up

Export Citation Format

Share Document