Modeling of Laser Pulse Heating of Solid Dielectric with a Fractal Structure

The nonlinear partial differential equation which describes heating of dielectric materials with fractal structure under the influence of pulse laser beam is received. The known equations from literature can be considered as special cases of the obtained equation. Fractality is considered in calculations of thermal conductivity coefficient, isochoric heat capacity and intensity of substance thermal radiation. The numerical decision of the equation is received. It is shown that the fractal properties of dielectric solid can essentially influence in heating process.


Introduction
Absorption of laser radiation occurs at exposure to laser radiation on matter. This process is directly related to the dielectric properties of matter and causes the heating and phase transitions. It is believed that the destruction of an optically transparent dielectric material by laser radiation is often precisely thermal nature [1]. Process of a laser heating of various solid bodies up to them destruction has the thermal nature is more often [2].
Dielectric materials with fractal structure have a great fundamental and applied interest. It is enough to mention that many polymers, ceramic compounds and carbon nanocomposites fall into such materials. Investigation of thermal and optical properties of these materials is important in technology of thin film coating for energy-saving technologies, electronics and medicine. Various questions of thermophysics of fractals were discussed in works [3][4][5].
Heating of non-polar dielectrics is caused, mainly, by phonons transfer and is described by partial differential equations of parabolic type [6]. These equations are sometimes modified using the fractional integro-differentiation [7] to simulate the heating of fractal structures which, however, involves considerable mathematical difficulties.
In this paper we propose simple model of heating of solid dielectric solid state. The physical approach based on nonlinear partial differential equation. The fractal properties of material are considered in the calculation of radiative energy losses and thermal conductivity coefficient depending on temperature which leads to a nonlinear heat equation.

Physical Characteristics of Solid Dielectric Fractals
Consider the important physical characteristics of the sample depending on the fractal structure. In the work [8] was obtained expression for the isochoric heat capacity which we here write in the form where D is the spectral fractal dimensionality, B k and A N are the Boltzmann and Avogadro constants, is the molar mass, T and θ are the absolute and characteristic temperatures. The thermal conductivity coefficient in the quantum statistical approach can be written [9]: where  is the Planck's constant, ω is the phonons frequency, is the characteristic frequency, const p τ = is the mean free time of phonons, V is the molar volume, ( ) g ω is the phonon states density, s c is the speed of sound and N is the number of phonons. Using the Bose-Einstein distribution and the expression for the phonon states density [8] from equation (2) we find 2 ( / ) 3 is the volume concentration of atoms.
During the heating of sample by laser beam the quantity of absorbed energy should be spent on thermal radiation. The substance will continuously cool by radiative energy losses even in the absence of convective heat transfer. We calculate the integrated intensity of thermal radiation to take into account this mechanism. For this we use the familiar expression: where c is the speed of light in a vacuum, ω is the frequency of radiation, σ is the photon absorption cross sections, E is the spectral density of blackbody radiation per unit volume and unit solid angle interval. The absorption cross section of some dielectric body having volume V is given by [10]: where 0 ε is the permittivity of free space, E is the electric field vector, A is the amplitude of the flux vector of the incident wave power density, Imε is the imaginary part of dielectric function which characterizes the absorption of radiation.
If the size of a single fractal cluster smaller than the wavelength of the incident radiation we can assume that the boundary conditions for electromagnetic wave coincides with the boundary conditions for solid particles. That is some averaging of the polarizability and the cluster behaves as isotropic with a density and dielectric properties depending on the fractal dimensionality. If we accept and make allowance for that integrand in (5) does not depend on the volume we get Im V c ω ε σ = For Planck distribution we have ( ) Substituting (6), (7) in (4), we find the total intensity of radiation We note that the integral in (8) diverges if at low frequencies the imaginary component of the dielectric function has a power asymptotic and tends to infinity faster than 4 1/ ω .
The specific properties dealing with the frequency dispersion of polarization and dielectric function are characteristic of polymers, composite materials and glasses having fractal structure. For the description of complex dielectric susceptibility and, accordingly, complex permittivity of these materials is commonly used Cole-Cole, Cole-Davidson and Havriliak-Negami formulas [11]. The adequacy of these equations has been confirmed repeatedly by experiments on the dielectric and capacitance spectroscopy, nuclear magnetic resonance, neutron scattering. The Havriliak-Negami complex dielectric function has a most general view: where a and b are the semi-empirical parameters taken into account the fractal structure implicitly, d τ is the dielectric relaxation time, 1 ε and 2 ε are the high-frequency and low-frequency dielectric permittivity ( 2 1 ε ε > ). Expression (9) describes of Debye dielectric relaxation processes at where it has been taken into account that Im 0 ε > . After substituting (10) into (8) and subsequent integration with the condition where ( ) This function is characteristic of many physical processes in the presence of self-similarity and fractality. It is pertinent to note that the power-temperature dependence with fractional exponents appear in fractal Debay theory of solid state heat capacity [8].

Laser Heating Modeling
Let us consider the following processes in the modeling of heating: 1) the absorption of laser radiation, 2) reflection of the laser radiation from the surface, 3) thermal radiation from the heated sample. Convective heat exchange with the surrounding space is not considered. This is true if the irradiation of the material is carried out in vacuum conditions. Convective heat exchange at the experiments in air with laser pulses of short duration is also negligible due to low heat. Under the condition of (1), (3), (11), as well as Beer-Lambert-Bouguer law the one-dimensional heat conduction equation takes the form Writing equation (12) we was also assumed that the fractal matter volume varies slightly at heating, so We assume that the same temperature is maintained at the ends of the irradiated sample which is equal to the ambient temperature:  (14) can not be solved analytically so it is interesting to consider one special case. We assume that the thin surface layer heat ( where s u is the maximum temperature, which is set in a time 0 t . Carrying out the integration in (16) we obtain If radiation heat exchange can be neglected then the solution (17) is simplified to If we also neglect the ambient temperature 0 u and make transition to dimensional quantities then equation (18) leads to where W=I L t 0 Ω is the laser energy, is the isobaric heat capacity corresponding to the generalized law of Dulong-Petit [5]. Formula (19) allows us to estimate approximately the temperature of the sample surface which is installed as a result of pulsed laser irradiation.
Computing experiment was carried out by means of MathCad program. Built-in operator Pdesolve was used to solve the boundary problem (14) for equation (12). Figures  1-6 shows some dependencies which were calculated to demonstrate the influence of physical parameters of a fractal solid body on the process of laser pulse heating. Each pair of figures relates to a change in a one selected parameter (it is indicated above) at fixed values of other parameters (they are indicated below.). Estimates show that the limiting value of the relative temperature / m m u T θ = at which there is melting is located in a sufficient narrow range from 2 to 2.5 for some dielectrics (quartz, sapphire, ruby, diamond). Therefore, the intensity of laser radiation in the simulation of heating should be selected in accordance with the condition  Fig. 1 shows that increasing of laser radiation intensity leads to increase of temperature and to reducing transition time of heating. Laser radiation absorption at the initial moments of time causes a nonuniform temperature profile and increases the time of onset of the stationary mode of heating as shown in Fig. 2. The results of modeling presented in Fig. 3 show that the sample can cool down effectively due to radiation heat exchange. In the transition from frequency independent dielectric to the complex permittivity described by Havriliak-Negami equation (9) the exponent m in the temperature dependence of the intensity of thermal radiation varies from 5 to 4. Similar change leads to a slight increase in heating temperature (Fig. 4). In low-dimensional systems atoms or molecules are connected to each other weakly and therefore more mobile. For this reason reduction of the fractal dimensionality leads to an increase in temperature as illustrated in Fig. 5. In addition, taking into account the dependence ( ) G u in equation (12) considerably reduces the steady-state mode of heating. Fig. 6 shows that the activity of laser radiation during a finite period of time 0 t leads to a heat pulse in the sample. At the same time rising and falling edges of this pulse is blurred due to unsteady heating and cooling processes.

Conclusions
Materials with fractal structure are varied physical and chemical properties that make them promising in real structures. We have implemented a comprehensive approach to the analysis of thermal kinetics processes in dielectric fractal solids. The model is valid in the absence of gas-dynamic effects and effects associated with the ionization and plasma formation. In present paper the following conclusions were formulated: 1) analytical expressions for the intensity of thermal radiation and thermal conductivity coefficient of dielectric fractal materials are obtained, 2) numerical solution of the heat equation shows that the fractal properties of the sample can significantly affect the heating process.