# Numerical modelling of the transient heat transport in a thin gold film using the fuzzy lattice Boltzmann method with α-cuts

### Alicja Piasecka Belkhayat

,### Anna Korczak

Journal of Applied Mathematics and Computational Mechanics |
Download Full Text |

NUMERICAL MODELLING OF THE TRANSIENT
HEAT TRANSPORT IN A THIN GOLD FILM USING
THE FUZZY LATTICE BOLTZMANN METHOD WITH *a*-CUTS

Alicja Piasecka Belkhayat, Anna Korczak

Institute of
Computational Mechanics and Engineering

Silesian University of Technology, Poland

alicja.piasecka@polsl.pl, anna.korczak@polsl.pl

**Abstract.** In
this paper a description of heat transfer in one-dimensional crystalline solids
is presented. The fuzzy lattice Boltzmann method based on the Boltzmann
transport equation is used to simulate the nanoscale heat transport in thin
metal films. The fuzzy coupled lattice Boltzmann equations for electrons and
phonons are applied to analyze the heating process of thin metal films via a
laser pulse. Such an approach in which the parameters
appearing in the problem analyzed are treated as constant values is widely
used. Here,
the model with fuzzy values of relaxation times and an electron-phonon coupling
factor
is taken into account. The problem formulated has been solved by means of the
fuzzy
lattice Boltzmann method using the *a*-cuts and the rules of
directed interval arithmetic.
The application of *a*-cuts allows one to avoid
complicated arithmetical operations in the fuzzy numbers set. In the final part
of the paper the results of numerical computations
are shown.

*Keywords:** fuzzy lattice
Boltzmann method, **a**-cuts, directed interval
arithmetic*

1. Introduction

The estimation of heat transfer processes
proceeding in metal micro-domains (e.g. thin metal films) subjected to an
ultrafast laser pulse is of vital importance
in microtechnology applications. It should be pointed out that taking into
account the extreme temperature gradients, extremely short duration and the domain
dimen-sions expressed in nanometers, the macroscopic heat conduction equation
basing on the Fourier law cannot be used [1-3]. The generalization of this law
resulting from the introduction of the delay time between the heat flux and
temperature gradient (relaxation time τ* _{q}*) leads to the well known Cattaneo-Vernotte equation [4] belonging to the group of hyperbolic PDE. The
alternative, more general, mathemati-
cal model called the dual phase lag equation (DPLE) takes into account
two delay times concerning both the heat
flux and the temperature gradient (the thermalization time τ

*). The DPLE contains a second order time derivative and higher order mixed derivative in both time and space (e.g. [5-8]). The other approach to the micro-scale heat transfer analysis involves the use of two-temperature models. The two-tempe-rature hyperbolic (parabolic) model consists of equations describing the temporal and spatial evolution of the lattice and electrons temperatures, the lattice and electron heat fluxes [9-11].*

_{T}In this paper the thermal processes
proceeding in the micro-domains are analyzed using the fuzzy lattice Bolzmann
method [12, 13]. The heat transfer problems are usually solved assuming
that the equations appearing in the mathematical model and all parameters in
these equations are deterministic. Such an assumption does not give an exact
image of the thermal processes met in the engineering practice. It seems more
natural to take into account uncertainties related, for example, to the
material parameters. Here, the fuzzy lattice Bolzmann method is presented using the *a*-cuts of fuzzy numbers with the approach of the directed interval arithmetic [14]. The interval values of relaxation times and boundary conditions are
taken into account. The application of *a*-cuts
allows one to avoid very complicated arithmetical operations in the fuzzy
numbers set because the *a*-cuts are closed
intervals. In this case, the mathematical operations are defined according to
the rules of the directed interval arithmetic performed for every *a*-cut, which is treated as an interval number. In this arithmetic, a set of proper
intervals is extended by improper intervals, and all arithmetic operations and
functions are also extended. The main advantage of the directed interval
arithmetic upon the usual interval arithmetic is that the obtained temperature
intervals are much narrower and their width does not increase in time.

2. Boltzmann transport equation

During the heating of thin metal films via
laser pulse the electrons are energized and they subsequently transfer the
energy to phonons via coupling between them. The Boltzmann transport equations
for the coupled model (1D problem) with two kinds of carriers: electrons (*e*)
and phonons (*ph*) take the following form [15]:

(1) |

(2) |

where are the carrier distribution functions, are the equilibrium distribution functions given by the Bose-Einstein statistic for phonons [16] and Fermi-Dirac statistic for electrons [17], are the frequency-dependent carrier propagation speeds, are the frequency-dependent carrier relaxation times and are the carrier generation rates due to external sources such as laser heating.

The BTEs (1) and (2) can be transformed into equivalent carrier energy density equations [15]

(3) |

(4) |

where are the carrier energy densities, * *are the
equilibrium carrier energy densities and are the carrier energy sources related to
a unit of
volume. The equations (3) and (4)
must be supplemented by the adequate boundary-
-initial conditions.

The electron and phonon energy densities at their equivalent nonequilibrium temperatures are given by the formulas

(5) |

(6) |

where are the carrier temperatures, *k** _{b}* is the Boltzmann constant, is the Fermi energy, is the electron density, Θ

*is the Debye temperature of the solid, is the number density of oscillators [15].*

_{D}The electron and phonon energy sources are calculated using the following expressions [15]:

(7) |

(8) |

where is the power density deposited by the external source function
associated with the laser irradiation [18] and *G* is the electron-phonon coupling
factor which characterizes the energy exchange between electrons and phonons.

3. Fuzzy numbers

The ground of the mathematical rules used in this paper is given by the fuzzy set theory. This approach is not common in solving heat transfer problems and that is why some of the definitions used in this concept must be explained [19].

First of all, the definition of a fuzzy set
will be introduced. The fuzzy set in
a non empty universal set ** **() can be expressed by a set of pairs
consisting of the elements and a certain
degree of pre-assumed membership
of the form

(9) |

where function is defined as

(10) |

In fuzzy sets, each element is mapped to [0,1] by membership function , where [0,1] means real numbers between 0 and 1 (including 0,1). Consequently, a fuzzy set is a ‘vague boundary set’ compared with a crisp set.

For every can be considered three types of membership to the fuzzy set :

· - full membership to the fuzzy set, ,

· - lack of membership to the fuzzy set, ,

· - partial membership to the fuzzy set.

The α-cut set in
universal set _{ }is made up of members whose membership is not
less than α for every [20]

(11) |

The value α is arbitrary and this *a*-cut set is
a crisp set. This set is determined by the following characteristic function:

(12) |

Every
fuzzy set can be defined as a sum of all its *a*-cuts

(13) |

where is a fuzzy set in the universe whose membership function is the following:

(14) |

Arithmetical operations are generally very complicated. Among the infinite quantity of possible fuzzy sets that can be qualified as fuzzy numbers, some types of membership functions are of particular importance.

Due to its rather simple membership function of a linear type, triangular fuzzy numbers are one of the most frequently used fuzzy numbers. In this paper the triangular fuzzy numbers will be used to solve the Boltzmann transport equations specified for various mathematical models.

A triangular fuzzy number is a set with the following membership function [21]:

(15) |

where is the core of the number, are the left and the right end of the number, respectively. A triangular fuzzy number can be written as .

One of the ways to avoid very complicated
arithmetic operations performed on fuzzy numbers is to apply *a*-cuts of
fuzzy numbers. In this case, the mathematical operations are defined according
to the rules of the directed interval arithmetic
performed for every *a*-cut.

The *a*-cut of a fuzzy number defined by a pair of functions and is
called a set of closed intervals [22]

(16) |

which satisfies the following conditions

(17) |

where () is a limited, monotonic function for every

Every fuzzy number can be presented as a sum
of all its own *a*-cuts

(18) |

The *a*-cut of a triangular fuzzy number is the set of all closed
intervals in the form

(19)

The decomposition of a fuzzy number allows
one to make the mathematical
operations on closed intervals which are *a*-cuts. In this
situation, complicated arithmetic operations can be omitted and it is possible
to apply the interval arithmetic for every *a*-cut.
The mathematical operations are simplified, because they are done only on the
ends of the intervals.

For the *a*-cuts
of two fuzzy numbers and the
following mathematical
operations can be defined () [22]:

– addition

(20) |

– subtraction

(21) |

– multiplication by a scalar ()

(22) |

– multiplication

(23) |

– inverse ()

(24) |

– division ()

(25) |

Applying *a*-cuts of the
fuzzy numbers allows one to use directed interval
arithmetic [13, 23].

4. Fuzzy lattice Boltzmann method

The lattice Boltzmann method is a discrete representation of the Boltzmann transport equations. The fuzzy form of the BTEs for 1D problems can be expressed as [24]

(26) |

(27) |

where are the fuzzy values
of carrier energy densities for electrons and phonons, respectively, * *are the fuzzy equilibrium
carrier energy densities and are the fuzzy
relaxation times.

The fuzzy values of the electron and phonon energy sources are calculated using the following formulas according to the rules of directed interval arithmetic [25]

(28) |

(29) |

where is the electron-phonon coupling factor.

For the one-dimensional model the discrete set of propagation velocities in two lattice directions (1 and 2) for electrons and phonons is defined as (see Fig. 1)

(30) |

(31) |

Fig. 1. One-dimensional 2-speed lattice Boltzmann model

The equations (26) and (27) should be supplemented by the boundary
conditions [12], for example of the 1^{st} type
on the left boundary and the 2^{nd} type on the right boundary as
follows:

(32) |

and the initial conditions

(33) |

where and are the fuzzy values of boundary temperatures of electrons and phonons respectively, and are the fuzzy temperatures of electrons in the first and second node, and are the fuzzy temperatures of phonons in the first and second node, is the electron’s initial temperature and is the phonon’s initial temperature.

Taking into account equations (30) and (31), the set of four fuzzy partial differential equations is obtained

(34) |

Introducing discretizing form, time and position derivatives may be written as follows [13]:

(35) |

The total energy density for electrons and phonons is defined as the sum of discrete energy densities in all the lattice directions and takes the form

(36) |

The equilibrium electron energy density and phonon energy density is the same in all lattice directions and can be calculated using the formula [26]

(37) |

(38) |

After subsequent computations the fuzzy lattice temperature is determined using the formula

(39) |

5. Numerical examples

As a numerical example the heat transport in a
gold thin film of the dimensions
200 nm has been analyzed. The following input data have been introduced:
the relaxation time the Debye temperature the coupling factor G = 2.3·10^{16} W/(m^{3}K),
the external heat source the fuzzy boundary conditions of the
1^{st} type on the left boundary and the 2^{nd} type on the
right bound-
ary the initial temperature The lattice step and the time step have been assumed.

Figures 2 and 3 illustrate the fuzzy
temperature distribution of electrons in the thin gold film for the chosen
times and chosen values of parameter *α*. Figure 4 presents the interval values of temperatures
for the chosen parameter α at the node corresponding
to the depth of 60 nm after 0.2 ps.

In the Figures 5 and 6 the heating curves at the internal nodes 1 (40 nm), 2 (60 nm) and 3 (80 nm) for chosen values of the parameter α are shown.

Fig. 2. The fuzzy temperature
distribution for *a** *= 0

Fig. 3. The
fuzzy temperature distribution for *a** *= 0.75

Fig. 4. The interval values of
temperatures for the chosen parameter *a*
at the node 60 nm after 0.2 ps

Fig. 5. Heating curves for *a** *= 0.75 at the chosen nodes

Fig. 6. Heating
curves for *a** *= 0.25 at the chosen nodes

6. Conclusions

In the paper the
Boltzmann transport equation with the fuzzy values of the relaxa-
tion times and the boundary conditions has been
considered. The fuzzy version
of the lattice Boltzmann method with *a*-cuts for solving 1D problems has been
presented.

Such an approach allows one to avoid
complicated fuzzy arithmetic and treat the considered fuzzy numbers as interval
numbers. For bigger values of *a*, the temperature interval is narrower.
For *a* = 1, the wideness of the temperature interval
is equal to 0.

The generalization of LBM allows one to find the numerical solution in the fuzzy form and such information may be important especially for the parameters which are estimated experimentally, for example the relaxation time.

Acknowledgement

*This work is supported by the project
BK-255/RMT4/2016.*

References

[1] Zhang Z.M., Nano/Microscale Heat Transfer, McGraw-Hill, New York 2007.

[2] Chen G., Borca-Tasciuc D., Yang R.G., Nanoscale heat transfer, [in:] Encycl. of Nanoscience and Nanotechnology, CA: American Scientific Publishers, Valencia, 7, 2004, 429-359.

[3] Smith A.N., Norris P.M., Microscale heat transfer, [in:] Heat Transfer Handbook, eds. A. Bejan, D. Kraus, John Wiley & Sons, 2003, Ch.18, 1309-1358.

[4] Cattaneo C., A form of heat conduction equation which eliminates the paradox of instantaneous propagation, Comp. Rend. 1958, 247, 431-433.

[5] Majchrzak E., Mochnacki B., Greer A.L., Suchy J.S., Numerical modeling of short pulse laser interactions with multi-layered thin metal films, CMES: Computer Modeling in Engineering and Sciences 2009, 41, 2, 131-146.

[6] Majchrzak E., Kałuża G., Heat flux formulation for 1D dual-phase lag equation, Journal of Applied Mathematics and Computational Mechanics 2015, 14(1), 71-78, DOI: 10.17512/ jamcm.2015.1.07.

[7] Majchrzak E., Mochnacki B., Sensitivity analysis of transient temperature field in microdomains with respect to the dual phase lag model parameters, International Journal for Multiscale Computational Engineering 2014, 12(1), 65-77. DOI: 10.1615 /IntJMultCompEng. 2014007815.

[8] Mochnacki B., Paruch M., Estimation of relaxation and thermalization times in microscale heat transfer model, Journal of Theoretical and Applied Mechanics 2013, 51, 4, 837-845.

[9] Majchrzak E., Dziatkiewicz J., Analysis of ultashort laser pulse interactions with metal films using a two-temperature model, Journal of Applied Mathematics and Computational Mechanics 2015, 14(2), 31-39. DOI: 10.17512/jamcm.2015.2.04.

[10] Dziatkiewicz J., Kuś W., Majchrzak E., Burczyński T., Turchan Ł., Bioinspired identification of parameters in microscale heat transfer, International Journal for Multiscale Computational Engineering 2014, 12(1), 79-89.

[11] Tzou D.Y., Macro- to Microscale Heat Transfer: The Lagging Behavior, Taylor and Francis, 1997.

[12] Piasecka-Belkhayat A., Korczak A., Application of the interval lattice Boltzmann method for a numerical modeling of thin metal films irradiation by ultra short laser pulses, IAPGOŚ 2014, 4, 85-88.

[13] Piasecka Belkhayat A., Korczak A., Modelling of transient heat transport in one-dimensional crystalline solids using the interval lattice Boltzmann method, [in:] Recent Advances in Computational Mechanics, eds. T. Łodygowski, J. Rakowski, P. Litewka, Taylor & Francis Group, A Balkema Book, London 2014, 363-368.

[14] Mochnacki B., Piasecka Belkhayat A., Numerical modeling of skin tissue heating using the interval finite difference method, MCB: Molecular & Cellular Biomechanics 2013, 10, 3, 233-244.

[15] Escobar R.A., Ghai S.S., Jhon M.S., Amon C.H., Multi-length and time scale thermal transport using the lattice Boltzmann method with application to electronics cooling, Journal of Heat and Mass Transfer 2006, 49, 97-107.

[16] Pisipati S., Geer J., Sammakia B., Murray B.T., A novel alternate approach for multiscale thermal transport using diffusion in the Boltzmann Transport Equation, International Journal of Heat and Mass Transfer 2011, 54, 3406-3419.

[17] Ashcroft N.W., Mermin N.D., Solid State Physics, Harcourt College Publishers, New York 1976.

[18] Venkatakrishnan K., Tan B., Ngoi B.K.A., Femtosecond pulsed laser ablation of thin gold film, Optics & Laser Technology 2002, 34, 199-202.

[19] Hanss M., Applied Fuzzy Arithmetic, Springer-Verlag, Berlin-Heidelberg 2005.

[20] Otto K., Lewis A.D., Antonsson E., Approximation alpha-cuts with the vertex method, Fuzzy Sets and Systems 1993, 55, 43-50.

[21] Giachetti R.E., Young R.E., A parametric representation of fuzzy numbers and their arithmetic operators, Fuzzy Sets and Systems 1997, 91, 185-202.

[22] Guerra M.L., Stefanini L., Approximate fuzzy arithmetic operations using monotonic interpolations, Fuzzy Sets and Systems 2005, 150, 5-33.

[23] Piasecka Belkhayat A., Interval boundary element method for imprecisely defined unsteady heat transfer problems, D.Sc. Dissertation, Gliwice 2011.

[24] Ghai S.S., Kim W.T., Escobar R.A., et al., A novel heat transfer model and its application to information storage systems, Journal of Applied Physics 2005, 97, 10P703.

[25] Lee J.B., Kang K., Lee S.H., Comparison of theoretical models of electron-phonon coupling in thin gold films irradiated by femtosecond pulse lasers, Materials Transactions 2011, 52, 3, 547-553.

[26] Escobar R., Smith B., Amon C., Lattice Boltzmann modeling of subcontinuum energy transport in crystalline and amorphous microelectronic devices, Journal of Electronic Packaging 2006, 128 (2), 115-124.