Next Article in Journal
Optimization of Process Parameters and Analysis of Microstructure and Properties of 18Ni300 by Selective Laser Melting
Next Article in Special Issue
Morphology and Anti-Corrosive Performance of Cr(III) Passivated Zn–Fe Alloy Coating on NdFeB Substrate
Previous Article in Journal
Study on the Influence of Defects on Fracture Mechanical Behavior of Cu/SAC305/Cu Solder Joint
Previous Article in Special Issue
Antimicrobial TiN-Ag Coatings in Leather Insole for Diabetic Foot
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Surface Topology on the Apparent Thermal Diffusivity of Thin Samples at LFA Measurements

by
Robert Szczepaniak
Faculty of Aviation, Polish Air Force University, Dywizjonu 303 Street No 35, 08-521 Deblin, Poland
Materials 2022, 15(14), 4755; https://doi.org/10.3390/ma15144755
Submission received: 7 May 2022 / Revised: 23 June 2022 / Accepted: 4 July 2022 / Published: 7 July 2022

Abstract

:
This paper deals with the problem of the influence of surface topography on the results of thermal diffusivity measurements when determined using the instantaneous surface heat source method, also called the pulse method. The analysis was based on numerical tests carried out using Comsol Multiphysics software. The results of experimental investigations on the actual material structure using an electron microscope, an optical microscope and a profilometer were used to develop a numerical model. The influence of the non-uniformity of the surface of the tested sample on the determined values of half-time of the thermal response of the sample’s rough surface to the impulse forcing on the opposing flat surface was determined by developing the data for simulated measurements. The effect of the position of the response data reading area on the obtained simulation results was also analyzed. The obtained results can be used to improve the accuracy of experimental heat transfer studies performed on thin-film engineering structures depending on the uniformity and parallelism of the material applied to engineering structures. The difference in half-life determination error results for various analyzed models can be as high as 16.7%, depending on the surface from which the responses of the heating impulse are read. With an equivalent model in which 10% of the material volume corresponds to the rough part as a single inclusion, hemisphere, the error in determining thermal diffusivity was equal to 3.8%. An increase in the number of inclusions with smaller weight reduces an error in the determination of thermal diffusivity, as presented in the paper.

1. Introduction

The process of coating the surface of an object or substrate with a very thin layer is referred to as coating. The layer can be a type of thin polymer sheet, paint, or varnish used for protective/decorative purposes. Most industrial products go through the coating process not only to prevent corrosion, but also to make them more attractive. The coating consists of the development of a thin layer that can be varnished or polymeric on a fabric or substrate, etc. [1].
Novas et al. [2] described the characteristics of the contributions made by scientists around the world in the field of solar coatings between 1957 and 2019. Photovoltaic systems depend greatly on the physical and chemical properties of their materials, the wavelength of the captured light, its intensity and angle of incidence, surface or texture characteristics, and the presence or absence of superficial coatings. In addition to these factors, temperature, pressure, machinability, durability, price, and lifetime cost are important in material selection.
To prevent dangerous failures, Cui et al. [3] developed a thin-film thermocouple temperature sensor to measure the temperature of the rolling elements of the bearing in real time while the train is running, bearing in mind the characteristics of the rapid change in the temperature of the bearings.
He et al. [4] investigated the coating technology of the DLC film for the orbiting screw of scroll compressors in air conditioning systems in order to reduce the friction coefficient and improve film adhesion to the substrate, wear resistance, and compressor efficiency. The friction coefficient of the DLC film is much lower than that of an anodic oxide film, as a result of which the DLC film has a good wear-reducing effect and excellent self-lubricating properties.
Mohammed et al. [5] investigated ZnS thin layers prepared by deposition in a chemical bath. They investigated the structural properties of films deposited and annealed at a temperature of 300 °C. This research found that relatively small grains were randomly formed and the deposited film had a heterogeneous surface with some cracks. Additionally, it was observed that the annealed samples had higher content of Zn than S, and the white color of the deposited thin layers did not change after annealing.
Poddighe et al. [6] addressed the problem of producing hydrophobic thin films from the liquid phase. Interest in the production of hydrophobic surfaces is constantly growing due to their wide application in several industries. Thin liquid-phase layers can be deposited on a variety of surfaces using a wide variety of techniques, and the design of the precursor solution offers the possibility to fine-tune the hydrophobic properties of the coating layers. The general trend is to design multi-functional films that have different properties in addition to being hydrophobic. The control of surface wettability is a key technological issue in several areas, such as microelectronics, separation membranes, car windows, self-cleaning surfaces, motion indicators, and biotechnology [7,8,9,10,11,12]. The leading technologies used for hydrophobic/hydrophilic thin film deposition and surface modification include chemical vapor deposition [13], laser ablation [14,15], and plasma treatment [16]. Superhydrophobic surfaces have a variety of uses in self-cleaning, anti-icing and anti-stick applications. Haj Ibrahim et al. [17] investigated the effect of surface topography on hydrophobic coatings both numerically and experimentally. Profilometry was used to create a numerical representation of the surface.
In addition to physical deposition techniques, the wet chemistry route, in particular sol–gel processing [18], has also met with great interest in the production of highly controlled hydrophobic-hydrophilic systems. Li et al. [19] proposed and improved an easy and cost-effective method involving acid etching and stearic acid self-assembly to successfully produce Al-based superhydrophobic fin tube heat exchangers. The 3D topography analysis showed suitable micro-nanostructures, whereas XPS and FTIR spectra showed chain self-assembly, essential for the realization of a superhydrophobic surface.
For temperature measurement, Dong et al. [20] produced a weak thin AlN layer with numerous defects. It was found that the thin film diffraction peak (002) increased monotonously with an increasing annealing temperature and an annealing time. This phenomenon is attributed to the evolution of defects in the AlN film network. Therefore, the relationship between defects and annealing can be expressed by the shift (002) of the diffraction peak, which can be used for temperature measurement. In addition, an algorithm for temperature interpretation was established, and software for temperature interpretation using MATLAB was also built. The temperature interpretation is realized by the software with a relative error of less than 7%. This study is of great importance in promoting the accurate measurement of the temperature at the surface of the high-temperature component. An accurate measurement of the surface temperature of a turbine blade and the control of its temperature distribution is an important basis for diagnosing turbine blade failure. The methods of measuring the surface temperature of turbine blades mainly include thin-film thermocouples [21,22,23], temperature indicating paint [24], infrared radiation [25,26,27], and irradiation crystals [28,29,30]. Based on the advantages of irradiating crystals, the authors [20] proposed an easy and cheap method of temperature measurement using a thin-film crystal. It was noted that the crystal quality could be strengthened after annealing [31,32,33], indicating that the AlN thin layer is a promising material for temperature measurement.
On the other hand, the thickness of the coating dramatically affects the functionality of the coatings. Accordingly, techniques used to determine thickness are of great importance for research and technology related to coatings. The use of thin films has become ubiquitous in many areas of the science and industry sectors. Coatings are widely used to obtain a synergistic effect between the characteristics of the substrate and the coating material. They can improve the physical, chemical, and aesthetic properties and lower the cost of the final product. For all these reasons, thickness measurement in composite materials is mandatory, both for obtaining the desired properties in the final component, and in order to keep costs under control [34,35].
The aim of this work is to address the problem of heat transfer in unevenly shaped samples based on previously published experimental studies in a temperature-sensitive paint (TSP) material [36]. It also deals with ongoing research on pressure-sensitive paint (PSP), which is dedicated to wind tunnel testing [37,38]. Due to the unproblematic application of this paint on surfaces, it is also used for industrial and scientific measurements [39,40,41]. The use of paint, particularly in wind tunnels, involves measurements that are quick-changing but with not much variation in the temperature measurement range; therefore, the authors of article [36] presented the thermophysical properties of the paint used in experimental studies in this field. However, during the measurements, the authors observed a large influence of the roughness of the prepared test samples on the measurement of thermal diffusivity determined using the classical method [42,43,44], as well as when using the modified method [45,46,47], especially of the sputtered or applied surface opposite to the flat substrate, similarly to Refs. [36,48,49,50].
In the case of heat transfer studies based on the determination of the global surface temperature distribution, it is necessary to observe, in the analysis, the thermal properties of all layers [36], as well as the very topography of the surface opposite to the area of influence of the laser pulse. Therefore, the thermal properties of the layers should be known as accurately as possible, as should the surface topography itself. TSP, similar to PSP, consists of molecular sensors embedded in a binder, usually polymer-based, which can affect the production of samples with a small cross-sectional dimension that automatically implies the problem of obtaining perfectly flat sample surfaces for experimental testing.
The determination of surface topography is particularly important when studying thin-film structures and the physical properties of materials whose property depends on the thickness of the sample (in particular, if the thickness is a square function, as is the case of determining thermal diffusivity by the laser flash method (LFA—Laser Flash Analysis), where the thickness is to the second power). Therefore, the determination of the thermal diffusivity of materials aerosolized on surfaces is not straightforward and is subject to high methodological error. By using an optical profilometer, it is possible to determine the topography and, at a later stage, to include it in numerical simulations to determine the accuracy of the thermal diffusivity parameter.
The use of numerical calculations by means of commercial software is increasingly popular, especially in mechanical engineering [51,52]. Comsol Multiphysics software is very frequently used for calculations in the area of heat transfer, as well as for the multi-physics simulation platform [53,54]. Recently it has also been particularly popular for heat transfer calculations in thin-film structures [55,56,57,58,59,60,61,62]. The influence of waviness is of great importance in calculations of heat transfer phenomena in thin films and in complex flow models used in, among other things, aviation technology, e.g., rocket nozzles [63], and in heat conduction itself in microchannel flows [64], which is a frequently tackled problem of modern mechanical engineering.
This study attempts to investigate the thermophysical properties of the Temperature Sensitive Paint (TSP) coating applied by spraying with an airbrush using numerical tests of substitute models. The error in determining the half-life was determined, and directly translates into the error in determining the thermal diffusivity of the coating.

2. Subject of Numerical Research

In the research presented here, a numerical analysis of heat transfer was carried out in models corresponding approximately to the structure of paint sprayed on a flat surface using UniTemp temperature-sensitive paint, manufactured by Innovative Scientific Solutions Incorporated (ISSI, Dayton, OH, USA). The thickness of the prepared model and the surface topography were determined to be massive and 90:10 (flat part of the sample: upper part—irregularities on the sample), which was simulated as inclusions in the form of hemispheres. Because of the problem of surface accuracy, the surfaces of this layer were checked using an optical profilometer (Figure 1a,b) and an electron microscope (Figure 1e,f). For testing purposes (for numerical calculations), a unit layer and maximum thickness deviations of approximately 10% were modeled (one variant of up to 20% was also made for checking). To illustrate the maximum deviation from parallelism of the sample, the roughness of the sputtered material is shown (Figure 1c,d). After a visual analysis, several geometric models were prepared. They are presented later in this article. The samples were modeled in such a manner that on one side there was a flat surface (in line with real-world testing), while on the opposite side, the surface was non-parallel, which corresponds to the replacement model of structure non-uniformity in real-world testing [36].
For the empirical studies, several samples were made with very different surfaces. In one batch of test samples for real-world testing, the coating thickness was not uniform and varied from 27 to 32 micrometers (Figure 1e,f). The inclusions stood out significantly, as can be seen in the electron microscope images. The densities and dimensions of the geometric TSP carriers were determined from the representative coating thickness. In a further batch of samples, the surface roughness of the TSP layer was examined and characterized by 2D and 3D surface scanning, using an FRT MicroProf 100 (FRT GmbH, Bergisch-Gladbach, Germany) optical profilometer. The morphology of the 3D surface is shown in Figure 2a. Figure 2b shows the results of the surface mapping. The surface profile along the red line in Figure 1c is shown in Figure 1d. The spectrum average is 45 µm, while 80% of the data points are within ±15 µm of this value.
Because the original material was applied to a flat substrate, one surface of the material created for the simulation is flat (the effect of the laser pulse causing the temperature rise will be simulated on this surface). A numerical simulation of impulsive forcing, applied to the determination of thermal diffusivity, such as in [65], will make it possible to assess the applicability of the finite element method to the influence of the parallelism of the surface of the material under investigation.

3. Methodology for Determining Thermal Diffusivity

The thermal diffusivity parameter, or thermal diffusivity coefficient, characterizes the behavior of a physical object under conditions of undetermined heat transfer and is a parameter determined experimentally in solid-state physics using an indirect method to specify thermal conductivity. Its value, through thermomechanical coupling, determines the magnitude of the loads to which a structure subjected to rapidly changing temperature may be exposed [66]. This parameter is difficult to determine for thin films due to the limitations of standard measurement methods, particularly for nonparallel sample surfaces. In physical terms, diffusivity is the ratio between the heat-transport properties and the heat storage capacity of a medium.
One of the most common methods for the determination of this parameter is the pulsed surface forcing method, which is based on the solution of a second-order Fourier differential equation with appropriate boundary and initial conditions in the sample, on which a laser pulse is applied to the front surface and the thermal response on the surface opposite to the thermal forcing is examined [42].
Despite the increasing number of modern laser techniques and non-contact temperature measurements in high-accuracy diffusometers [65], it is still quite important to prepare samples with flat parallel surfaces for testing. In accordance with the test methodology, the samples have parallel surfaces (Figure 2), and the measurement of the temperature rise is measured surface-to-surface from the side opposite to the laser forcing. Currently, there is no methodology for developing test results in software if there are uneven inclusions. Hence, unless it is possible to prepare flat test specimens, an estimate of the measurement error of such materials should be introduced. Therefore, this paper specifies the percentage error in the determination of the thermal diffusivity value of a material with known thermophysical parameters based on the example of previous experimental studies using the commercial software Comsol Multiphysics by determining the half-life.
In Parker’s pulse method, there is a simple relationship between sample thickness and half-life defined by Relation (1). For the purpose of this work, the half-life will be determined in accordance with the methodology shown in Figure 2.
The thermal diffusivity in Parker’s method is determined from the following dependence [42]:
a = 1.38 · l 2 π 2 · t 0.5
where l is the sample thickness, t0.5 is the time after which half of the maximum excess temperature on the sample’s rear surface has been reached.

4. Numerical Model

Several geometrical models of the samples were prepared for the numerical tests. All models were designed in such a manner that the elements occurring above the parallel surfaces represent 10% of the total sample mass. Assuming that the material is homogeneous, and that the experimental tests were carried out for axially symmetrical samples, for the purposes of numerical calculations, the sample has a unit dimension in each direction. The model was assumed to have isothermal edges on the outer sides, and the laser pulse was simulated by releasing energy onto the bottom surface. As in the experimental study, the samples were covered with flake graphite to improve the absorption of the laser flash energy on the shot side as well as on the opposite side to improve the registration of the infrared thermal response on the readout side of the response in the numerical simulation. The reading took place under ideal conditions in which no additional well-conductive layers were required.
A few sample elements corresponding to 10% of the mass of the whole sample were prepared. The tested models are shown in Figure 3a–e. The geometric dimensions are given in Table 1.
Geometric models were created in the manner previously described, and with test parameter settings as below:
The subdomain settings were as follows:
  • − General equation of heat transfer used in the model in Comsol Multiphysics:
ρ c p T t + ( k T ) = 0
  • − Thermophysical properties of the tested material (Table 2):
Table 2. Thermophysical properties of the tested material.
Table 2. Thermophysical properties of the tested material.
No Thermophysical Parameter SymbolValueUnit
1Thermal conductivityk0.3W/m K
2Densityρ1500kg/m3
3Heat capacity at constant pressurec1100J/kg K
4TemperatureT(t0)273.15k1
5Inward heat flux0100,000 (<0.01)W/m2
  • − Element settings for temperature: Lagrange–Quadratic.
The boundary settings were as follows:
All planes were aligned symmetrically except for the bottom plane, where a laser pulse of less than 10 ms with an inward heat flux of 105 W/m2 was simulated (Figure 4).
Computational grid settings. A calculation grid with tetrahedral elements was used (Figure 5).

5. Findings of Numerical Research

Numerical tests were carried out for seven different models (Table 1) with pre-determined boundary and starting conditions. In the numerical test procedure, the bottom surface corresponds to a sample aerosolized onto a flat metal plate, so it is flat, and the laser pulse was simulated on it. The temperature waveform from the nonparallel top surface was read at 0.0001 s intervals in the first temperature rise range (up to 0.02 s) and then the reading frequency was reduced to 0.0005 s.
With the models prepared in this way, half-life was checked for each model and for the parallel surfaces of the model (without the sphere part), as well as for the sphere itself. In addition, results were also given for a model with 20 per cent volume in the hemisphere element (Figure 6a,b). The obtained half-life results are included against the figures for each model. A comparison of extreme temperature increments at two points, i.e., at the parallel surface and at the highest point on the hemisphere, is also presented, in combination with the temperature change at the surface affected by the thermal forcing (Figure 6c). An exemplary model of the temperature distribution in the sample is also presented (Figure 6d). Next, the values for all the models tested are presented in a graphical form (Figure 6e) as a function of the number of hemispheres depicted as black dots on the graph. The red dots correspond to the time value t0.5 for the sample with 20% inclusion volume for both the sphere and the base and similarly the green color to the time value t0.5 for the sample with 10% share volume.
The value of the half-life for the flat-parallel sample is 13.2 ms, while each deviation from parallelism introduces an error and, depending on the number of inclusions, this parameter changes. The values of these changes are shown in the table below (Table 3).

6. Discussion

The presented procedure complements the determination of thermal diffusivity error for rough or uneven samples in the broad procedure of pulsed thermal testing, which further develops Parker’s method. In Parker’s method, various additional real effects accompanying the study of the heat conduction phenomenon are taken into account, except the influence of surface non-uniformity. Therefore, the methodology presented in this paper makes it possible to estimate the error of the determination of the thermal diffusivity parameter. With the additional technique of imaging the topography with, e.g., an optical profilometer, it is possible to determine an uneven volume in the sample to determine the spread of the error in the determination of the correct value.
The combination of numerical modeling and the finite element method with experimental tests and the use of a profilometer will allow a more accurate determination of thermal diffusivity values especially for materials that are problematic to manufacture with regard to the recommendations of the test methodology. Numerical measurements including surface non-uniformity can contribute to increasing the accuracy of thermal diffusivity measurements by pulse forcing methods. With this approach, it will be possible to determine the thermal properties of a material not only for perfectly manufactured samples, but for a whole range of materials that are difficult to apply in such a way as to produce perfectly flat opposite surfaces, such as paints. In previous research [36], the authors presented results assuming the parallelism of the samples [36], providing information on roughness and, as a complement to this, the potential for error due to surface irregularity. The method presented here shows the effect of bumps in the form of hemispheres whose volume corresponds to 10% (or 20%) of the volume of the sample roughness made for empirical testing. Any form of surface roughness can be used during testing once the geometric model has been implemented in the numerical calculations.
The obtained results illustrate the difference between the thermal diffusivity of the investigated material and the effective thermal diffusivity (TD), i.e., on the basis of the TD obtained from processing the measurement data of the material structure study when the structure geometry differs from the model one, the obtained results have to be interpreted as effective.
This research also provides a quantitative result. The differences in the determined diffusivity value of the resultant structure can be interpreted as actual accuracy data. The differences in the determination methodology, which involves using surface topography test data to specify the volumetric contribution of surface structures and, on this basis, determining the deviation of the test result from the true value, can be used in real-world testing.
The methodology for determining differences, which involves using the data of surface topography to specify the volume contribution of surface structures and, on this basis, determining the deviations of the test outcome from the true value, can be exploited in real testing. In this way, it is possible to determine the value of the corrections to the measurement results. In scientific works conducted based on numerical calculations of thin layers, authors have tried to simplify numerical models due to the extremely large number of computational points, which does not always translate into improvement of the obtained results, as presented in the work of Yin S et al. [67]. The obtained results are correlated with the author’s own research [36], and represent a quantitative correction for empirical data. Therefore, this paper focuses on the introduction of simplified computational models to illustrate the phenomenon of heat transfer itself and the determination of calculation errors in determining the half-life resulting from surface roughness.
In addition, the author tried to prove that the error in determining thermal diffusivity (indirectly by determining the half-life) is less than 3% when the entire volume in the non-uniform sample is replaced with the average thickness. It should be borne in mind that it is not possible to produce perfectly flat surfaces, and the roughness does not have to be uniform in accordance with the Gaussian distribution during the production of thin layers, as presented, inter alia, in the work of Z Ebrahiminejad et al. [68]. The author hopes that the methodology presented by him will be used to make computational corrections to the thermophysical parameter, thermal diffusivity, during the implementation of experimental research on similar subjects.

7. Conclusions

This paper presents the results of a study on the thermophysical properties of TSP coating applied by airbrush spraying using numerical replacement model studies. The error in determining half-life was determined, which translates directly into the error in determining the thermal diffusivity of the material. The numerical studies of the TSP paint were complemented by the results for TSP structures obtained by the paint casting method, which initiated the analysis of the issue. This analysis has contributed to a better understanding and determination of the thermal diffusivity of a material with a non-uniform TSP shell geometry and, at the same time, can be a methodology for testing any other material with non-planar topography. The difference in half-life determination error results for the different analyzed models can be as high as 16.7%, depending on the surface from which the responses are read on the surface opposite the impulse forcing surface. When reading the temperature from the entire rear surface, the error with the analyzed models was equal to 3.8% for the replacement model, where 10% of the volume was in one hemisphere. However, the greater the number of inclusions and the closer the surface is to the ideal sample, the smaller is the error will be. Thus, it was confirmed that this method is particularly suitable for materials that are layered on top of other materials, with the simultaneous problem of maintaining the opposite surfaces in parallel, which is crucial in determining the exact heat transfer parameters. The presented methodology as well as the test results themselves may be helpful in improving the accuracy of heat transfer tests using TSP, PSP, or other methods.

Funding

The scientific work was carried out in the framework of the project “Experimental and numerical studies of the basic physical parameters of composite materials used in aviation technology” No. 451/RN/2020, founded by the Polish Air Force University for 2021.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ellahi, R. Recent Trends in Coatings and Thin Film: Modeling and Application. Coatings 2020, 10, 777. [Google Scholar] [CrossRef]
  2. Novas, N.; Alcayde, A.; El Khaled, D.; Manzano-Agugliaro, F. Coatings in Photovoltaic Solar Energy Worldwide Research. Coatings 2019, 9, 797. [Google Scholar] [CrossRef] [Green Version]
  3. Cui, Y.; Gao, P.; Tang, W.; Mo, G.; Yin, J. Adaptive Thin Film Temperature Sensor for Bearing’s Rolling Elements Temperature Measurement. Sensors 2022, 22, 2838. [Google Scholar] [CrossRef] [PubMed]
  4. He, Z.; Ji, L.; Xing, Z. Experimental Investigation on the DLC Film Coating Technology in Scroll Compressors of Automobile Air Conditioning. Energies 2020, 13, 5103. [Google Scholar] [CrossRef]
  5. Mohammed, R.Y. Annealing Effect on the Structure and Optical Properties of CBD-ZnS Thin Films for Windscreen Coating. Materials 2021, 14, 6748. [Google Scholar] [CrossRef] [PubMed]
  6. Poddighe, M.; Innocenzi, P. Hydrophobic Thin Films from Sol–Gel Processing: A Critical Review. Materials 2021, 14, 6799. [Google Scholar] [CrossRef]
  7. Men, X.; Zhang, Z.; Yang, J.; Zhu, X.; Wanga, K.; Jiang, W. Spray-coated superhydrophobic coatings with regenerability. New J. Chem. 2011, 35, 881–886. [Google Scholar] [CrossRef]
  8. Manoudis, P.N.; Karapanagiotis, I. Modification of the wettability of polymer surfaces using nanoparticles. Prog. Org. Coat. 2014, 77, 331–338. [Google Scholar] [CrossRef]
  9. Dhere, S.L.; Latthe, S.S.; Kappenstein, C.; Pajonk, G.M.; Ganesan, V.; Rao, A.V.; Waghe, P.B.; Gupta, S.C. Transparent water repellent silica films by sol–gel process. Appl. Surf. Sci. 2010, 256, 3624–3629. [Google Scholar] [CrossRef]
  10. Prabhu, K.N.; Fernades, P.; Kumar, G. Effect of substrate surface roughness on wetting behaviour of vegetable oils. Mater. Des. 2009, 30, 297–305. [Google Scholar] [CrossRef]
  11. Shang, H.M.; Wang, Y.; Limmer, S.J.; Chou, T.P.; Takahashi, K.; Cao, G.Z. Optically transparent superhydrophobic silica-based films. Thin Solid Film. 2005, 472, 37–43. [Google Scholar] [CrossRef]
  12. Gu, G.; Dang, H.; Zhang, Z.; Wu, Z. Fabrication and characterization of transparent superhydrophobic thin films based on silica nanoparticles. Appl. Phys. 2006, 83, 131–132. [Google Scholar] [CrossRef]
  13. Liu, H.; Feng, L.; Zhai, J.; Jiang, L.; Zhu, D. Reversible Wettability of a Chemical Vapor Deposition Prepared ZnO Film between Superhydrophobicity and Superhydrophilicity. Langmuir 2004, 20, 5659–5661. [Google Scholar] [CrossRef] [PubMed]
  14. Luo, B.H.; Shum, P.W.; Zhou, Z.F.; Li, K.Y. Preparation of hydrophobic surface on steel by patterning using laser ablation process. Surf. Coat. Technol. 2018, 204, 1180–1185. [Google Scholar] [CrossRef]
  15. Thorhallsson, A.I.; Fanicchia, F.; Davison, E.; Paul, S.; Davidsdottir, S.; Olafsson, D.I. Erosion and Corrosion Resistance Performance of Laser Metal Deposited High-Entropy Alloy Coatings at Hellisheidi Geothermal Site. Materials 2021, 14, 3071. [Google Scholar] [CrossRef] [PubMed]
  16. Woodward, I.; Schofield, W.C.E.; Roucoules, V.; Badyal, J.P.S. Super-hydrophobic Surfaces Produced by Plasma Fluorination of Polybutadiene Films. Langmuir 2003, 19, 3432–3438. [Google Scholar] [CrossRef]
  17. Haj Ibrahim, S.; Wejrzanowski, T.; Przybyszewski, B.; Kozera, R.; García-Casas, X.; Barranco, A. Role of Surface Topography in the Superhydrophobic Effect—Experimental and Numerical Studies. Materials 2022, 15, 3112. [Google Scholar] [CrossRef]
  18. Zhou, T.; Cheng, X.; Pan, Y.; Li, C.; Gong, L.; Zhang, H. Mechanical performance and thermal stability of glass fiber reinforced silica aerogel composites based on co-precursor method by freeze drying. Appl. Surf. Sci. 2018, 437, 321–328. [Google Scholar] [CrossRef]
  19. Li, R.; Wang, Z.; Chen, M.; Li, Z.; Luo, X.; Lu, W.; Gu, Z. Fabrication and Characterization of Superhydrophobic Al-Based Surface Used for Finned-Tube Heat Exchangers. Materials 2022, 15, 3060. [Google Scholar] [CrossRef]
  20. Dong, L.; Li, Y.; Lv, J.; Jiang, H.; Zhang, W. Fabrication of Weak C-Axis Preferred AlN Thin Film for Temperature Measurement. Sensors 2021, 21, 5345. [Google Scholar] [CrossRef]
  21. Tillmann, W.; Kokalj, D.; Stangier, D.; Schoppner, V.; Malatyali, H. Effects of AlN and BCN Thin Film Multilayer Design on the Reaction Time of Ni/Ni-20Cr Thin Film Thermocouples on Thermally Sprayed Al2O3. Sensors 2019, 19, 3414. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Shi, P.; Wang, W.; Liu, D.; Zhang, J.; Ren, W.; Tian, B.; Zhang, J. Structural and Electrical Properties of Flexible ITO/In2O3 Thermocouples on PI Substrates under Tensile Stretching. ACS Appl. Electron. Mater. 2019, 1, 1105–1111. [Google Scholar] [CrossRef]
  23. Zhao, X.; Liang, X.; Jiang, S.; Zhang, W.; Jiang, H. Microstructure Evolution and Thermoelectric Property of Pt-PtRh Thin Film Thermocouples. Crystals 2017, 7, 96. [Google Scholar] [CrossRef] [Green Version]
  24. Li, Y.; Li, Z. The Research of Temperature Indicating Paints and its Application in Aero-engine Temperature Measurement. Procedia Eng. 2015, 99, 1152–1157. [Google Scholar]
  25. Usamentiaga, R.; Venegas, P.; Guerediaga, J.; Vega, L.; Molleda, J.; Bulnes, F.G. Infrared thermography for temperature measurement and non-destructive testing. Sensors 2014, 14, 12305–12348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Castro, P.; Lecuna, R.; Manana, M.; Martin, M.J.; Campo, D.D. Infrared Temperature Measurement Sensors of Overhead Power Conductors. Sensors 2020, 20, 7126. [Google Scholar] [CrossRef]
  27. Piccinini, F.; Martinelli, G.; Carbonaro, A. Reliability of Body Temperature Measurements Obtained with Contactless Infrared Point Thermometers Commonly Used during the COVID-19 Pandemic. Sensors 2021, 21, 3794. [Google Scholar] [CrossRef]
  28. Viswanathan, E.; Kanjilal, D.; Sivaji, K.; Ganapathy, S. Identification of sublattice damages in swift heavy ion irradiated N-doped 6H-SiC polytype studied by solid state NMR. J. Phys. Chem. B 2011, 115, 7766–7772. [Google Scholar] [CrossRef]
  29. Jin, E. Grain boundary effects on defect production and mechanical properties of irradiated nanocrystalline SiC. J. Appl. Phys. 2012, 111, 2303. [Google Scholar] [CrossRef]
  30. Ruan, Y.F.; Wang, P.F.; Huang, L.; Zhu, W. High-Temperature Sensor Based on Neutron-Irradiated 6H-SiC. Key Eng. Mater. 2011, 495, 335–338. [Google Scholar] [CrossRef]
  31. Kai, C.; Zang, H.; Ben, J.; Jiang, K.; Shi, Z.; Jia, Y.; Cao, X.; Lü, W.; Sun, X.; Li, D. Origination and evolution of point defects in AlN film annealed at high temperature. J. Lumin. 2021, 235, 118032. [Google Scholar] [CrossRef]
  32. Xiao, S.; Suzuki, R.; Miyake, H.; Harada, S.; Ujihara, T. Improvement mechanism of sputtered AlN films by high-temperature annealing. J. Cryst. Growth 2018, 502, 41–44. [Google Scholar] [CrossRef]
  33. Hagedorn, S.; Walde, S.; Mogilatenko, A.; Weyers, M.; Cancellara, L.; Albrecht, M.; Jaeger, D. Stabilization of sputtered AlN/sapphire templates during high temperature annealing. J. Cryst. Growth 2019, 512, 142–146. [Google Scholar] [CrossRef]
  34. Giurlani, W.; Berretti, E.; Innocenti, M.; Lavacchi, A. Measuring the Thickness of Metal Coatings: A Review of the Methods. Coatings 2020, 10, 1211. [Google Scholar] [CrossRef]
  35. Van Zele, M.; Watté, J.; Hasselmeyer, J.; Rijckaert, H.; Vercammen, Y.; Verstuyft, S.; Deduytsche, D.; Debecker, D.P.; Poleunis, C.; Van Driessche, I.; et al. Thickness Characterization Toolbox for Transparent Protective Coatings on Polymer Substrates. Materials 2018, 11, 1101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Panas, A.J.; Szczepaniak, R.; Stryczniewicz, W.; Omen, L. Thermophysical Properties of Temperature-Sensitive Paint. Materials 2021, 14, 2035. [Google Scholar] [CrossRef]
  37. Burner, A.W.; Liu, T.; Garg, S.; Bell, J.; Morgan, D. Unified model deformation and flow transition measurements. J. Aircr. 1999, 35, 898–901. [Google Scholar] [CrossRef]
  38. Liu, T. Pressure- and Temperature-Sensitive Paints, Encyclopedia of Aerospace Engineering; John Wiley and Sons, Ltd.: Hoboken, NJ, USA, 2011. [Google Scholar] [CrossRef]
  39. Zhou, W.; Peng, D.; Liu, Y.; Hu, H. Assessment of film cooling’s surface quantities using pressure- and temperature-sensitive paint: Comparisons between shaped and sand-dune inspired holes. Exp. Therm. Fluid Sci. 2019, 101, 16–26. [Google Scholar] [CrossRef]
  40. Watkins, A.N.; Buck, G.M.; Leighty, B.D.; Lipford, W.E.; Oglesby, D.M. Using Pressure- and Temperature-Sensitive Paint for Global Surface Pressure and Temperature Measurements on the Aftbody of a Capsule Vehicle. AIAA J. 2009, 47, 821–829. [Google Scholar] [CrossRef]
  41. Ghorbani-Tari, Z.; Chen, Y.; Liu, Y. End-wall heat transfer of a rectangular bluff body at different heights: Temperature-sensitive paint measurement and computational fluid dynamics. Appl. Therm. Eng. 2017, 122, 697–705. [Google Scholar] [CrossRef]
  42. Parker, W.J.; Jenkins, R.J.; Butler, C.P.; Abbott, G.L. Flash method of determining thermal diffusivity, heat capacity, and thermal conductivity. J. Appl. Phys. 1961, 32, 1679–1684. [Google Scholar] [CrossRef]
  43. Vozár, L.; Hohenauer, W. Flash method of measuring the thermal diffusivity, a review. High Temp. High Press. 2004, 36, 253–264. [Google Scholar] [CrossRef]
  44. Massard, H.; Pinto, C.S.C.; Couto, P.; Orlande, H.R.B.; Cotta, R.M. Analysis of flash method physical models for the measurement of thermal diffusivity of solid materials. In Proceedings of the 10th Brazilian Congress of Thermal Sciences and Engineering—ENCIT 2004, ABCM, Rio de Janeiro, Brazil, 29 November–3 December 2004; p. CIT04-0537. [Google Scholar]
  45. Terpiłowski, J.; Rudzki, R.; Szczepaniak, R.; Woroniak, W. Investigation of thermal diffusivity of Fe80Ni20 alloy by means of modified pulse method. J. Alloys Compd. 2018, 735, 560–568. [Google Scholar] [CrossRef]
  46. Terpiłowski, J.; Jóźwiak, S.; Rudzki, R.; Szczepaniak, R.; Woroniak, G. Investigation of Phase Transformation of Fe65Ni35 Alloy by the Modified Pulse Method. Materials 2020, 13, 3425. [Google Scholar] [CrossRef]
  47. Hemberger, F.; Ebert, H.P.; Fricke, J. Determination of the Local Thermal Diffusivity of Inhomogeneous Samples by a Modified Laser-Flash Method. Int. J. Thermophys. 2007, 28, 1509–1521. Available online: http://thermophysics.ru/pdf_doc/AutoPlay/Docs/CollectionOfManuscripts/ECTP2005paper96.pdf (accessed on 22 April 2022). [CrossRef]
  48. Panas, A.J.; Stryczniewicz, W.; Szczepaniak, R. Thermophysical property investigation of a thin layered paint. Thermochim. Acta 2018, 622, 100–107. [Google Scholar] [CrossRef]
  49. Stryczniewicz, W.; Panas, A.J. Numerical data processing from a laser flash experiment on thin graphite layer. Comput. Assist. Methods Eng. Sci. 2015, 22, 279–287. [Google Scholar]
  50. Grochalski, K.; Wieczorowski, M.; Pawlus, P.; H’Roura, J. Thermal Sources of Errors in Surface Texture Imaging. Materials 2020, 13, 2337. [Google Scholar] [CrossRef]
  51. Singh, S.; Saha, A.K. Numerical study of flow and heat transfer during a high-speed micro-drop impact on thin liquid films. Int. J. Heat Fluid Flow 2021, 89, 108808. [Google Scholar] [CrossRef]
  52. Su, Y. Modeling and Characteristic Study of Thin Film Based Biosensor Based on COMSOL. Hindawi Publ. Corp. Math. Probl. Eng. 2014, 11, 581063. [Google Scholar] [CrossRef]
  53. Zhou, J.; Huang, J.; Liao, J.; Guo, Y.; Zhao, Z.; Liang, H. Multi-field simulation and optimization of SiNx:H thin-film deposition by large-size tubular LF-PECVD. Sol. Energy 2021, 228, 575–585. [Google Scholar] [CrossRef]
  54. Yadav, H.N.S.N.; Kumar, M.; Kumar, A.; Das, M. COMSOL simulation of microwave plasma polishing on different surfaces. Mater. Today Proc. 2021, 45, 4803–4809. [Google Scholar] [CrossRef]
  55. Zandi, S.; Saxena, P.; Razaghi, M.; Gorji, N.E. Simulation of CZTSSe Thin-Film Solar Cells in COMSOL: Three-Dimensional Optical, Electrical, and Thermal Models. IEEE J. Photovolt. 2020, 10, 1503–1507. [Google Scholar] [CrossRef]
  56. Qin, Z.; Zeng, Y.; Norton, P.R.; Shoesmith, D.W. Modeling Hydrogen Permeation through a Thin TiO2 Film Deposited on Pd Using COMSOL Multiphysics. In Proceedings of the COMSOL Conference, Boston, MA, USA, 8–10 October 2009. [Google Scholar]
  57. Nan, C. Numerical Modeling of Flexible ZnO Thin-Film Transistors Using COMSOL Multiphysics. Master’s Thesis, University of Waterloo, Waterloo, ON, Canada, 2013. [Google Scholar]
  58. Tesar, J.; Semmar, N. Thermal effusivity of metallic thin films: Comparison between 1D multilayer analytical model and 2D numerical model using COMSOL. In Proceedings of the 5th European Thermal-Sciences Conference, Eindhoven, The Netherlands, 18–22 May 2008. [Google Scholar]
  59. Saxena, P.; Gorji, N.E. COMSOL Simulation of Heat Distribution in Perovskite Solar Cells: Coupled Optical–Electrical–Thermal 3-D Analysis. IEEE J. Photovolt. 2019, 9, 1693–1698. [Google Scholar] [CrossRef]
  60. Yantchev, V.; Turner, P.; Plessky, V. COMSOL modeling of SAW Resonators. In Proceedings of the IEEE International Ultrasonics Symposium Proceedings, Tours, France, 18–21 September 2016. [Google Scholar]
  61. Du, S.-Y.; Zhao, Y.-H. Numerical study of conjugated heat transfer in evaporating thin-films near the contact line. Int. J. Heat Mass Transf. 2012, 55, 61–68. [Google Scholar] [CrossRef]
  62. Al-Shurman, K.M.; Naseem, H. CVD Graphene Growth Mechanism on Nickel Thin Films. In Proceedings of the 2014 COMSOL Conference in Boston, Boston, MA, USA, 8–10 October 2014. [Google Scholar]
  63. Budakli, M.; Gambaryan-Roisman, T.; Stephan, P. Gas-driven thin liquid films: Effect of interfacial shear on the film waviness and convective heat transfer. Int. J. Therm. Sci. 2019, 146, 106077. [Google Scholar] [CrossRef]
  64. Zhou, L.; Zhou, S.; Du, X.; Yang, Y. Heat transfer characteristics of a binary thin liquid film in a microchannel with constant heat flux boundary condition. Int. J. Therm. Sci. 2018, 134, 612–621. [Google Scholar] [CrossRef]
  65. Min, S.; Blumm, J.; Lindemann, A. A new laser flash system for measurement of the thermophysical properties. Thermochim. Acta 2007, 45, 46–49. [Google Scholar] [CrossRef]
  66. Orłoś, Z. Naprężenia cieplne; PWN: Warszawa, Poland, 1991. [Google Scholar]
  67. Yin, S.; Guo, W.; Huang, D.; Zhang, H.; Wu, H.; Zhang, C. Simulation Study on Heat Transfer Characteristics and Cooling Water Design for the Casting Film. Huanan Ligong Daxue Xuebao/J. South China Univ. Technol. Nat. Sci. 2021, 49, 23–34. [Google Scholar] [CrossRef]
  68. Ebrahiminejad, Z.; Asgary, S.; Esmaili, P. Surface characterization of Cu-doped indium sulfide thin films. Indian J. Phys. 2022, 96, 1315–1319. [Google Scholar] [CrossRef]
Figure 1. Results of the TSP coating surface morphology investigation: (a,b) three-dimensional surface image, (c) typical profile (red line cross-section), (d) diagram of a normalized dimensional spectral analysis result, (e,f) SEM images of TSP layer fragment separation for microscopic structure inspection.
Figure 1. Results of the TSP coating surface morphology investigation: (a,b) three-dimensional surface image, (c) typical profile (red line cross-section), (d) diagram of a normalized dimensional spectral analysis result, (e,f) SEM images of TSP layer fragment separation for microscopic structure inspection.
Materials 15 04755 g001
Figure 2. Heat transfer model for laser radiation coordinate system and initial conditions: (a) semi-transparent sample, (b) temperature rise on back surface and determination of half-life, where: g—effective sample thickness, x—Cartesian coordinate of a one-dimensional sample.
Figure 2. Heat transfer model for laser radiation coordinate system and initial conditions: (a) semi-transparent sample, (b) temperature rise on back surface and determination of half-life, where: g—effective sample thickness, x—Cartesian coordinate of a one-dimensional sample.
Materials 15 04755 g002
Figure 3. Exemplary geometrical models in the calculation program Comsol: (a) 1 hemisphere, (b) 3 hemispheres, (c) 5 hemispheres, (d) cube, (e) part of hemisphere.
Figure 3. Exemplary geometrical models in the calculation program Comsol: (a) 1 hemisphere, (b) 3 hemispheres, (c) 5 hemispheres, (d) cube, (e) part of hemisphere.
Materials 15 04755 g003
Figure 4. Setting the boundary conditions for an exemplary geometric model with 4 hemispheres.
Figure 4. Setting the boundary conditions for an exemplary geometric model with 4 hemispheres.
Materials 15 04755 g004
Figure 5. An example of a computational mesh for a geometric model with 4 hemispheres..
Figure 5. An example of a computational mesh for a geometric model with 4 hemispheres..
Materials 15 04755 g005
Figure 6. Distribution of temperature rise on the rear surface of the samples: (a) flat surface (without sphere), (b) spherical surface, (c) temperature distribution in selected characteristic points in the model, (d) example of temperature distribution in the sample after a laser shot, (e) half-life for different variants tested numerically.
Figure 6. Distribution of temperature rise on the rear surface of the samples: (a) flat surface (without sphere), (b) spherical surface, (c) temperature distribution in selected characteristic points in the model, (d) example of temperature distribution in the sample after a laser shot, (e) half-life for different variants tested numerically.
Materials 15 04755 g006aMaterials 15 04755 g006b
Table 1. Geometric dimensions of numerical models in relation to unit value.
Table 1. Geometric dimensions of numerical models in relation to unit value.
No Type of InclusionNumber of Hemispheres Radius of Hemisphere [mm]Hemisphere Share in Total Mass [%]
1spherical000
2spherical10.365510
3spherical20.290110
4spherical30.253510
5spherical40.230010
6spherical50.213710
7spherical10.460520
8spherical part10.4810
9cube10.4642—cube edge10
Table 3. Half-life and error of its determination for different configurations of the tested material.
Table 3. Half-life and error of its determination for different configurations of the tested material.
Number of HemispheresPlace of Measurementt0.5 [ms]Error [%]
0Whole back surface13.20
1Whole back surface 13.73.8
2Whole back surface 13.63.0
3Whole back surface 13.41.5
4Whole back surface13.30.8
5Whole back surface13.260.5
1From the surface of the hemisphere 10%14.912.9
1From the surface of the hemisphere 20%14.59.8
1From the surface around a parallel sphere 10%12.18.3
1From the surface around a parallel sphere 20%11.016.7
1Whole rear surface 20%13.73.8
1Whole back surface (inclusion cube)12.83.0
1Whole back surface (part of sphere—10%)13.52.3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Szczepaniak, R. Effect of Surface Topology on the Apparent Thermal Diffusivity of Thin Samples at LFA Measurements. Materials 2022, 15, 4755. https://doi.org/10.3390/ma15144755

AMA Style

Szczepaniak R. Effect of Surface Topology on the Apparent Thermal Diffusivity of Thin Samples at LFA Measurements. Materials. 2022; 15(14):4755. https://doi.org/10.3390/ma15144755

Chicago/Turabian Style

Szczepaniak, Robert. 2022. "Effect of Surface Topology on the Apparent Thermal Diffusivity of Thin Samples at LFA Measurements" Materials 15, no. 14: 4755. https://doi.org/10.3390/ma15144755

APA Style

Szczepaniak, R. (2022). Effect of Surface Topology on the Apparent Thermal Diffusivity of Thin Samples at LFA Measurements. Materials, 15(14), 4755. https://doi.org/10.3390/ma15144755

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop