Next Article in Journal
Bridge Inspection and Defect Recognition with Using Impact Echo Data, Probability, and Naive Bayes Classifiers
Next Article in Special Issue
Experimental Study on the Out-of-Plane Behavior of Brick Masonry Walls Strengthened with Mortar and Wire Mesh: A Pioneer Study
Previous Article in Journal
Model Bridge Span Traversed by a Heavy Mass: Analysis and Experimental Verification
Previous Article in Special Issue
BAT Algorithm-Based ANN to Predict the Compressive Strength of Concrete—A Comparative Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Defining the Thermal Features of Sub-Surface Reinforcing Fibres in Non-Polluting Thermo–Acoustic Insulating Panels: A Numerical–Thermographic–Segmentation Approach

1
Institute of Process Equipment and Control Engineering, Zhejiang University of Technology, Hangzhou 310023, China
2
Independent Researcher, 67100 Santa Rufina di Roio, Italy
3
Engineering School of Nondestructive Testing, Tomsk Polytechnic University, 634050 Tomsk, Russia
4
Department of Chemical Engineering, National Tsing Hua University, Hsinchu 300044, Taiwan
5
Department Industrial and Systems Engineering, Khalifa University, Abu Dhabi 127788, United Arab Emirates
6
Department of Industrial and Information Engineering and Economics (DIIIE), University of L’Aquila, 67100 L’Aquila, Italy
*
Authors to whom correspondence should be addressed.
Submission received: 4 August 2021 / Revised: 5 September 2021 / Accepted: 8 September 2021 / Published: 10 September 2021
(This article belongs to the Special Issue Sustainability of Building Materials and Structures)

Abstract

:
Natural fibres present ozone-friendly solutions in the field of construction. The attenuation of the sound and heat losses is an important feature in such type of materials above all, when used in non-woven fabrics and fibre-reinforced composites. Hemp fibres show robust insulation performance; this research work should be considered beneficial to the development of a non-destructive thermographic methodology, which can address the thermal barrier (typically applied on multi-layer panel) effects. The intent is to assess the integrity of the sub-surface reinforcing glass fibres; such integrity state will help confer the rigidity and the resistance to mechanical stresses. The testing proposed in this study can be further developed in a laboratory right after the manufacturing process of similar type of components. The testing needs preliminary numerical simulations to help guide the selection of the appropriate pre- and post-processing algorithms combined with or without segmentation operators. A set of numerical and experimental tests were performed through controlled thermal stimulation while recording the thermal responses. The study also highlights the advantages, disadvantages, and future development of the presented technique and methodologies.

1. Introduction

Non-destructive testing (NDT) of composites, including active infrared thermography (IRT), remains an important issue in civil engineering [1]. IRT needs a thermal stimulus to generate thermal imprints [1]. In many construction structures, the composite core is a reinforcing grid made up of basalt [2], carbon [3], or glass fibres [4]. Infrared NDT has proven to be efficient in inspecting glass fibres, but there is still a lack of knowledge in using IRT techniques to assess the composite integrity, when combining glass fibres with hemp fibres scattered in a matrix, to be used as protective coating.
In recent years, a variety of IRT techniques have attracted the attention of many researchers thanks to its capabilities in fatigue assessment, defect detection, and damage repair applications for fibre components or products. The application of IRT for inspecting wind turbine blades, made of GFRP and sandwich panels, was investigated in [5], where surface defects, such as air bubbles, edge bonding, etc., were efficiently detected in blades. In [6], the use of lock-in IRT method for detecting small temperature changes caused by energy dissipation and crack propagation during the tensile test was studied. The evaluation of the IRT effectiveness in identifying different types of defects in masonry structures was performed by Angiuli et al. [7]. Thermographic surveys were also repeated after repair works to identify an acceptance criteria that can be used on areas subjected to total or partial repair. Pitarresi et al. suggested combining techniques of a modified transverse crack tension test coupon and IRT in order to evaluate the static and fatigue behaviour of GFRP composites [8]. In 2020, Xu et al. explored the feasibility of eddy current pulsed thermography to detect sub-surface defects in GFRP composites, a typical non-conductive material. The method introduced a metal part as the heat-generating source to produce the induction heat [9]. The work of Rani and Mulaveesala highlighted advantages of pulse compression in modulated thermal wave imaging, which can be carried out with a moderate peak-power heat source during a certain time period with improved test resolution and sensitivity [10].
Beyond the use of a single IRT method, the combined application of IRT with other NDT techniques, such as acoustic emission (AE), ultrasound imaging, and three-dimensional (3D) computed tomography, demonstrated advantages in numerous scenarios and became an effective solution for NDT practitioners. Karger-Kocsis et al. studied two types of glass fibre composites by applying the AE and IRT [11,12]. Czigány et al. applied AE and IRT on glass fibre-reinforced poly (ethylene terephthalate) composite sheets, concluding that both techniques are suitable for damage detection and assessment of the fibre structure [12]. Ibarra-Castanedo et al. investigated the possibility of detecting defects and assessing the impact severity of delaminations in glass fibre-reinforced polymer (GFRP) by using pulsed IRT and vibro-thermography [13]. Katunin et al. tested various techniques, including piezoelectric ultrasound, IRT, and vibration-based NDT methods, to analyse their applicability for the inspection of aircraft components [14]. Palumbo et al. characterized GFRP adhesive joints by applying ultrasonic imaging and lock-in IRT analysis for the assessment of adhesion quality in composites subjected to static tensile mechanical tests and accelerated aging cycles [15]. Crupi et al. investigated the failure mode and the integrated damage of the sandwich panels by using three-dimensional (3D) computed tomography and IRT. The combined experimental and theoretical studies demonstrated that the use of glass–epoxy reinforcement in aluminium honeycomb sandwiches enhances both energy absorption and load carrying capacities [16]. In 2021, Quattrocchi et al. compared air-coupled ultrasonic testing and active IRT by performing experimental tests on two reference composite panels (GFRP and CFRP) with synthetic defects having different sizes and located at various depths. The pros and cons of the two techniques were highlighted in terms of defect size sensitivity [17].
The continuous development of detection technologies is important to ensure the accommodation of new applications with high precision requirements. Obviously, it is difficult to achieve the desired results via a single detection technique. At the same time, with the development of machine learning algorithms, more image processing methods are used as a post-processing of IRT scans. In order to enhance and expand IRT detection performance and capabilities, Ghorbel et al. investigated the evolution of temperature distributions in glass fibre-reinforced polyamides. In addition to IRT, digital image correlation (DIC) was used to quantify deformed zones and correlate them to heat dissipation sources [18]. Omer and Sultan proposed a new method for post-processing of IRT data by converting temperature read-outs into a pulse train and applying the discrete Fourier transform. Promising results were obtained even for small defects located at the depth of 10 mm below the surface [19]. To increase the detection capability of IRT, a novel IR image processing technique involving the complex Morlet wavelet phase analysis was developed in [20]. It was demonstrated, using the concept of signal-to-noise ratio, that the proposed technique was able to separate the low- and high-frequency noise components from the thermal image signal, thus reliably identifying defects at different depths. Yanjie et al. used the phase Fourier analysis to determine defect location and applied the logarithmic second derivative time to calculate defect depth by using the technique of long pulse thermography [21]. The experimental results were compared with the numerical simulations of a GFRP panel containing fabricated defects [21]. Vijaya Lakshmi et al. proposed a classification- and a regression-based quantitative post-processing technique to characterize sub-surface defects in carbon fibre-reinforced polymers (CFRP) and GFRP specimens by using quadratic frequency modulated thermal wave imaging [22]. In 2019, Sutthaweekul et al. determined the location of flat bottom holes in coated GFRP pipes by using synthetic aperture radar tomography in combination with principal component analysis (PCA), while the defect depth was characterized by analysing the time-of-flight characteristics of the background-subtracted responses [23].
Most of the aforementioned studies are based on using IRT technologies as the primary approach. In fact, numerical simulations and mathematical models are often used in conjunction with experiments to guide and calibrate the test process or assist in developing the pre- and post-processing algorithms. For example, Pracht and Swiderski described the results of numerical simulations in the IRT inspection of multi-layer composite structures, including GFRP [24]. Bernegger et al. presented an approach to quantify delaminations in semi-transparent GFRP. The experimental data was reconstructed by fitting experimental and theoretical profiles using a one-dimensional mathematical model. The delamination depths were estimated with a reasonable accuracy [25]. Hu et al. investigated the characterization and imaging of localized thickness loss (LTL) in GFRP using microwave NDT. A two-dimensional (2D) finite element model with the Ka-band open-ended waveguide and GFRP sample subject to LTL were set up and adopted for this analysis [26].
Herein, we describe specific numerical–experimental tests conducted on a multi-layer specimen. The results were integrated into a segmentation procedure to overcome existing shortcomings in traditional methodologies. Considering the current limited scientific knowledge about non-polluting thermo–acoustic materials [27,28,29] combined with GFRP (that falls into the field of “science and technology for the built environment”), the novelty of presented research will be described in terms of the numerical and experimental techniques and accompanying results. The simulation in Comsol Multiphysics® software was performed following a purposeful approach aimed at the minimization of computational cost, and the thermal tomography technique developed at Tomsk Polytechnic University (Russia) allowed the detection of detachments at different depths. Surprisingly, the use of local tangent space alignment (LTSA) technique, when applied to the thermographic data and then combined with the Canny segmentation operator, allowed a complete shape (morphology) retrieval of a grid pattern (made by glass fibres), which was located beneath an insulator material in hemp fibres.
The specimen has therefore been analysed by using two numerical models and two different thermal stresses. This helps the reader to understand pros and cons of each of them as the added value of the work. The useful considerations regarding the masking features of noise in respect to the thermal visualization of sub-surface features are presented in the text; also, a comparison between innovative post-processing algorithms applied to the analysis of thermal images (to our knowledge, this is a novel technique applied to composites in the construction sector) is also presented.
The remainder of the article is organized as follows: In Section 2, the description of the multi-layer specimen and main information about numerical–thermographic tests are supplied. Section 3 and Section 5 describe the modalities of numerical simulations performed by Comsol Multiphysics® and ThermoCalc-3D software; Section 4 and Section 6 describe the different sets of produced thermal images. Section 7 concludes the study by highlighting the outcome of testing and simulation and discussing the results.

2. The Multi-Layer Specimen, the Numerical–Thermographic Tests

The numerical model was developed to study heat conduction peculiarities in a specimen made of materials used in civil construction. Specifically, the specimen was mainly fabricated with a load-bearing base of an expanded Styrofoam (see point 5 both in Table 1 and in Figure 1a) and a cement milk layer (see point 4 both in Table 1 and in Figure 1a). A reinforcing grid of fiberglass was embedded in the cement milk layer (see point 3 both in Table 1 and in Figure 1a). This aimed to homogenize possible mechanical stresses that might appear in the base structure and propagate toward the more superficial finishing layers. Moreover, thermal stresses resulting from having different materials with different expansion effects can appear in multi-layer structures as the in-plane deformations of varying magnitude. The studying of such material is important due to its favourable mechanical (i.e., very flexible, etc.) and chemical (i.e., not subject to oxidative phenomena, etc.) nature, yet no established testing regimes are available in the literature.
Above the grid and the cement milk, a layer of cement mortar (see point 2 both in Table 1 and in Figure 1a) was applied while leaving an open space (see Figure 1). In this space, the grid is visually detectable (see point 3 both in Table 1 and in Figure 1a). Such sample design allowed evaluation of the thermal insulation effects conditioned by the cement mortar. Above the cement mortar, the layer of hemp fibres (see point 1 both in Table 1 and in Figure 1a) was present in a part of the inspected sample. The length of this layer was only 8 cm, but it covered the entire width of the specimen. This cover layer highlights the thermal insulation effect of the hemp fibres. For complete details, a numerical description of the stratigraphy in terms of size and thermophysical properties is presented in Table 1.
Figure 1 shows a couple of visual images (a,b), and the numerical models analysed in the first test (c,d).
The specimen was placed on a low-conductive support and illuminated by two projectors equipped with Osram Siccatherm® 250 (W) lamps and located at −45° and +45°, while the axis of each lamp was oriented towards the centroid of the specimen.
The thermal imaging camera was placed on the bisector of the total angle at a greater distance than that between the specimen and the projectors. In this setup, the camera viewing direction did not intercept the light beams from the projectors during the heating period. The testing scheme was simulated in Comsol Multiphysics®. In particular, two modelling cases were investigated, of which choice was dictated by the high computational cost required for analysing the model, which will be described below.
In the first modelling, the measurement environment was reproduced by placing the projectors and the specimen at appropriate distances by preserving the real projector-specimen angles.
Figure 2a shows the first test scheme.
In the second test (Figure 2b), the specimen was investigated by applying the classical one-sided thermal non-destructive testing (TNDT) procedure. The sample was heated with a 2.5 kW halogen lamp for 5 s. The temperature distributions were recorded using an Optris PI450 IR imager with an acquisition frequency of 15 Hz.

3. Modelling via Comsol Multiphysics® Software

The first test was modelled using Catia® V5 processor reproducing both the heating sources and the specimen without reconstructing the specimen stratigraphy. The external control volume—open boundary—was also defined, and the corresponding governing equations will be explained below.
The modelling defined the electrical power ensured by the lamp filaments. In particular, having selected tungsten as the filament material, the electrical load curve was modelled in Comsol environment to provide Joule heating of the lamps in a time span of 60 s. In order to analyse thermal transients during the cooling phase, the software was developed to evaluate the thermal phenomena when the power supply ceased for 180 s, i.e., simulating the cooling stage. Noting that natural convection was also taken into account.
The model did also accommodate the evaluation of heat distribution and associated fluxes on both the projectors and the sample by taking into account radiation view factors.
All the information obtained from the specimen surface temperature mapping (seen as a temperature footprint) was transferred in the form of input data into the second numerical model. In the framework of the second model, it was possible to analyse the behaviour of the sample in stratigraphic terms without considering the physical environment. Such procedure seemed to be efficient and most appropriate because a hypothetic general model combining the optical features of heating and the entire sample stratigraphy would have resulted in a high computational cost in terms of the accepted numerical mesh (by requiring a RAM of at least 64 Gb [30]).
Figure 3 shows both the image of the mesh (Figure 3a) and the results from the numerical modelling (Figure 3b,c) inherent to the thermal distribution induced by the projectors on the specimen at certain computation times.
Figure 3a shows the tetrahedral mesh, which affects the heating sources with a notable thickening in terms of the number of nodes. Building a fairly dense mesh was necessary to analyse the behaviour of the heating effects in terms of the volumetric analysis of heat flux distribution. Similar mesh specifications were applied to the specimen. In particular, it included a triangular mesh, thus allowing the simulation on the sample sides without calculating sample internal structure. A volumetric tetrahedral mesh was set for the outer cube representing the control volume. For this mesh, a thinning factor of 1.5 was accepted, thus allowing the dense mesh both in the proximity of the agitation source and on the surfaces of the specimen and becoming sparse far from the sample. Such approach has considerably reduced the computational cost without a detailed simulation of the control volume external surfaces.
From Figure 3b, it is possible to evaluate the presence of the control volume represented by the external cube that had the unique function to limit the volume in the analysis of thermal fields. In order to simplify the display, it was decided to use the slice function with the section plane passing through the resistances of the heating sources.
The aforementioned modelling strategy allowed us to study the temperature trends of the tungsten filament. In addition, the isocontour function was used to evaluate the increase in air temperature and its distribution. It is evident that the temperature was in a wide range starting from the maximum value (detected at the tungsten filaments of the heating sources) down to the ambient air temperature of the control volume.
Figure 3c shows the thermal field as generated by the heating sources at 10 s. Both the high temperature from the tungsten filament and the effect of the environment are well demonstrated.
The concept to build two separate models was prompted by the complexity of the specimen shape and its unique internal structure. The apparent regularity of the external dimensions (parallelepiped) and the thickness of the layers constituting the model, however, has prevented the use of a sweep mesh. Therefore, the first model, i.e., the specimen with a swept envelope, was discretized by means of a triangular mesh including only 428 nodes, while the second model was used to meet some problems described below.
The main constraint that complicated the use of a regular mesh along the specimen thickness was the presence of the reinforcing grid. The latter, for a finite areal extension, has no layer of material on it, thus challenging implementation of a swept mesh for all layers of the model. In fact, there was no correspondence between the upper and lower surfaces in terms of the number of nodes.
A tetrahedral approach was used for producing the mesh, which would fit the reinforcing grid. Given its geometrical dimensions, both a nodal thickening equal to 0.0001 m and a resolution of narrow regions equal to 0.25 were necessary. These values were obtained through a series of iterations aimed to evaluate the most efficient thickening with reasonable computational efforts and associated processing cost.
Discretization of the cement milk layer also imposed some difficulties. The addition of the warp and weft of the fiberglass inside the cement milk resulted in the formation of a series of square-based parallelepipeds. The mesh of the cement milk layer required the same degree of accuracy as the fiberglass. Given the regularity of the bottom of the layer (a single parallelepiped), it was possible to reduce the number of nodes through controlling the maximum element grow rate parameter in the direction of the depth. Concerning the rest of the layers, the mesh was built following the automatic generation process. The model simulating the whole specimen had a number of nodal elements equal to 372 346. Because of the mesh complexity involving a high degree of freedom (DOF), it was difficult to ensure a single processing step, which would involve the physics of radiation, including both the projectors and the effects of global convection. In fact, the quantities calculated in the first model were added to the second model, as both the convection and the mutual irradiation of the specimen towards the outside ambient were validated. Furthermore, in the second model, the control volume was no longer present. The latter was not necessary in the first modelling effort, as only the mutual radiation towards infinity was involved. The simplification of the specimen in terms of the first model mesh and the thickening of the specimen DOF in the complete second model prevented the use of the multi-component technique typically provided by Comsol Multiphysics®. If the specimen remained unchanged, the sets of equations could have been solved separately within the same model at different times by avoiding the physical passage of the variables.
The first model consisted of the heat transfer physics with surface-to-surface radiation set of equations as follows:
ρ C p T t + ρ C p u · T 1 = · k T + Q + Q r ,
where ρ indicates the density in (kg/m3), Cp is the specific heat at constant pressure (J/kg·K), T the represents temperature (K), u represents the velocity field tensor (m/s) (in our case, term 1 of Equation (1) is null because the components are mutually immobile), Q indicates the generic heat source/sink (J), and Qr is the radiative heat source/sink (J). On the other hand, Equation (2) represents the basis of the radiative phenomena conditioned by the projectors:
S i · I i 1 = k I b T 2 β I i 3 + σ s 4 π j = 1 N ω j I j S i , S j 4 .
The term 4 of Equation (2) shows the presence of a discretization for the angular space; the term 4 is usually presented in the form of an extended integral from 0 → 4π to analyse all spatial directions. In this case, as there were only two projectors providing the thermal load, the term 4 degenerates into numerical quadratures of discrete directions. Within the summation, σs indicates the scattering coefficient divided by 4π, as the function is spatial. The remaining part, reported as a member of the summation, indicates the phase function. The latter evaluated the probability that a generic ray from the direction Si is projected in the direction Sj. The definition of the phase function depends on the material constituting a single layer, while ωj indicate the j-directions.
The first model does not include a detailed view of the specimen (i.e., internal structure), yet its main purpose is to compute the radiation effects of the projectors and environment appearing on the surfaces of the specimen. To fully understand this concept, it is necessary to refer to the corresponding problem by following Rosseland [31].
In Equation (2), the summation has N as its upper term, but in the considered case, it is equal to 2. In particular, the term 1 of Equation (2) represents the gradient of the radiative intensity I as the i-th component projected along the S direction. Note that, in the term 2 of Equation (2), k represents the absorbed fraction of the radiative intensity evaluated for the black body Ib(T). In the term 3 of Equation (2), the fraction β indicates the radiant intensity Ii with respect to the generic i-th direction.
Of particular interest are the boundary conditions of the control volume in the first model, which represent open boundaries being defined by Equation (3):
T = T 0 i f n · u = < 0 T · n = 0 i f n · u = 0
where n is the normal outgoing to the control volume.
Both the lighting bodies and the surfaces of the simplified specimen have boundary conditions of the type:
n · k T = ε G e b T ,
where G indicates the incoming radiative heat flux or irradiation (W/m2), eb(T) is the blackbody hemispherical total emissive power (W/(m2·T3)), and ε provides the surface emissivity, i.e., a dimensionless number in the range 0 ≤ ε ≤ 1.
In particular, as in the first model, the specimen is simplified in terms of geometry. The aim is to set the maximum value of G for each specimen surface; therefore, ε is equal to 1, i.e., being independent of the material. In Equation (4), the k parameter represents the thermal conductivity. Instead, the second model, where the sample is reproduced in all details involves the physics named Heat Transfer with Radiation in Participating Media. In the framework of this physics, the values of G produced by the calculation on the surfaces of the simplified specimen (see the first model) were considered as the input values for the entire model yet to be analysed. However, in the latter model, both the real emissivities of each material and the actual specimen shape (e.g., the geometry of the reinforcement grid) were taken into account to calculate the real thermal load. Then, the governing equation of G on the surfaces becomes:
q = G ε   e b T .
Equation (5) allows the calculation of the thermal energy on each surface exposed to the thermal load through considering the actual emissivity and other radiative effects of the media involved. In addition, the set of equations intended for the analysis of the porous materials calculated by means of the physics for porous average in scattered fibres, cement mortar, cement milk, and Styrofoam is given below:
ρ C p e f f T t + ρ C p u · T = · k e f f T + Q + Q v d + Q p ,
where,
ρ C p e f f 1 = θ p ρ p C p , p 2 + 1 θ p 3   ρ C p
and k is
k e f f = θ p k p + 1 θ p   k ,
where Qvd in Equation (6) indicates the energy of the viscous dissipations in the form of heat (J), while Qp is the point heat source (J). In Equation (7), the term 1 indicates the effective volumetric heat capacity at constant pressure, while the term 2 shows θp, i.e., the dimensionless value indicating the fraction of porosity present in the material, ρp is the density of the porous material (kg/m3), and Cp,p represents the specific heat capacity (J/(kg·K)) for the porous material. In Equation (7), the term 3 indicates the non-porous percentage of the material expressed as a dimensionless number. Finally, in Equation (8), kp indicates the thermal conductivity of the porous material.
Figure 4 shows the mesh linked to the second numerical model. This figure is a visual display demonstrating the meshing strategy deployed for the sample.

Numerical Results

The analysis was carried out on the second model, where both the heating and cooling phases can be studied (Figure 5 shows the specimen heating phase).
Figure 5a shows the temperature range calculated in the framework of the second model by using the first model input data. The entire volume of the specimen was subjected to thermal analysis. The slicing technique was not used to take advantage of the volume perspective. It is seen that after only 10 (s) of heating, the sample front surface was thermally loaded in a prevalent way. The calculated thermograms illustrate that the lowest temperature increase appeared on the bottom of the specimen. The insulating effect of hemp fibres was evident and prevailing in relation to the remaining layers. The layers of cement mortar (both in lower and upper positions) were more affected by the temperature increase than the cement milk layer. It is evident that the reinforcement grid, which was installed along the entire specimen, reached a higher local temperature in regard to the cement milk layer.
In Figure 5b, it is possible to notice a uniformity of the temperature field on the front surface, with a reduced thermal gradient among all points up to 1.1 (K). On the one hand, the grid below the cement milk layer became less evident due to the homogenization of the temperature differences between the surface layer and the remaining subsequent layers. On the other hand, it is possible to note the effect of temperature gradient increase along the specimen thickness due to the thermal load, which is mainly caused by the front surface. In Figure 5c, the sub-surface grid below the cement mortar layers (lower and upper position), is less evident. This depends on homogenization of the temperature related to the surface layers.
The diffusion of the thermal front can be also observed along the specimen depth. Figure 5d displays a further homogenization in the temperature map specific to the glass fibre grid; such trend was not evident below the cement milk layer, where the grid contours were clear and well defined. Throughout the remaining part of the specimen, the sub-surface reinforcement grid became less observable.
In Figure 5e, the diffusion process involved both layers of the cement mortar; therefore, it was difficult to define the grid geometry underneath them. On the other hand, the mesh under the cement milk layer was still visible. This can also be attributed to the thickness of this layer, which provided a lining to the glass fibre.
Figure 5f shows the temperature distribution at the end of the heating period. However, the maximum specimen temperature was observed later because of the thermal inertia of the projectors, which continued to illuminate the sample surface. This phenomenon is illustrated by Figure 6a showing the temperature distribution at 10 s after the projectors were switched off.
The temperature distribution in Figure 6b is fairly uniform except the cement milk layer. The hemp fibre area is characterized by a small temperature signal of −0.37 K that can be explained by natural convection. The similar temperature pattern with the temperature gradient on the surface of −0.07 K is shown in Figure 6c; also, the in-depth penetration of the temperature front is clearly seen. Respectively, in Figure 6d, the footprint of the reinforcement grid beneath the cement milk layer becomes less pronounced, and this phenomenon evolves in time (Figure 6e) to reach the temperature gradient between the glass fibres and cement milk of about 0.04 K at 120 s or 180 s after heating started (Figure 6f). In both Figure 6g,h the specimen displays a uniform temperature distribution for all the layers of interest. The reinforcement grid is barely visible below the cement milk layer, and the respective temperature difference is below 0.03 K. To further analyse the in-depth temperature, four different directions were identified. The corresponding segments shown in Figure 7 passed through the intersections of the grid weft and warp in order to calculate influence of glass fibres. The corresponding temperature profiles are presented in Figure 8.
As follows from Figure 8a, there are certain differences between the profiles calculated for the different materials at 10 s. Of particular interest is the comparison between the two layers of cement mortar where the temperature profiles are conditioned by the material thickness. Lower surface temperatures are observed on the surface of cement milk and hemp layers. Respectively, some temperature variations take place at the location of the reinforcing grid.
In general, the penetration depth of the temperature front can be evaluated to be around 8 mm; see the temperature profiles in Figure 8b–f. For example, the temperature increase at the depth of the reinforcing grid is negligible while comparing the two layers of cement mortar. However, between the cement milk layer and the hemp fibres, there is a temperature difference of 0.75 K. A noticeable increase of surface temperature (0.1 K) appears over hemp fibres to compare to the cement milk area (Figure 8c). The abovementioned features of heat conduction in the specimen remain the same for all times (Figure 8d,e).
In Figure 8e, the trends of the temperature field shows a very slight increase. On the one hand, the behaviour of the directrix that affected the hemp fibres can be noticed. On the other hand, no particular differences concerning the abovementioned trends are visible. This did not happen between the cement milk and the hemp fibre layers, as the latter increased its temperature until it exceeded the surface temperature of the former. The same effect was also evident in regard to the temperature of the reinforcement grid.
The profiles in Figure 8f correspond to the end of the heating period, demonstrating a general temperature increase on the sample surface, in particular, in the area over the hemp fibres where the temperature exceeded that of the cement mortar by about 0.5 K.
Figure 9 illustrates changes in the temperature profiles at the cooling phase along the same directions as shown in Figure 7. As mentioned above, for a certain period of cooling, the surface temperature was slightly increasing due to thermal inertia of the projectors (Figure 9a–d) with a slightly higher temperature over the hemp fibres. This can be explained by the higher thermal insulation properties of the hemp fibres. In general, all profiles were similar to those observed within the heating stage, but the temperature penetration depth apparently increased up to about 19 mm at 120 s. In their turn, the profiles in Figure 9e–f show that the residual heating from the projectors ceased at about 180 s; therefore, one may expect the appearance of the temperature inversion phenomenon conditioned by the layers’ different heat capacities. Note also that, at ~180 s, the temperature front reached the depth of 29 mm, which is about the specimen thickness.

4. Processing Experimental Thermograms: The First Session

Experimental IR image sequences (spatial format of 126 × 240), which were provided by active thermal NDT, represented 3D “data cubes” to be further transformed into 2D matrices (spatial format of 30,240 × 240); see Figure 10.
In Figure 10, X is the unfolded 2D matrix. At the next processing step, the data matrix was normalized by subtracting the average of each column (column centralization), and the normalized matrix was still denoted as X. In X, the rows and columns represented samples (thermal images) and variables (pixels), respectively. Then, three manifold learning methods, namely isometric feature mapping (ISOMAP), local linear embedding (LLE), and LTSA, were used to reduce the dimension of X, respectively. It should be noted that it was not necessary to analyse all thermal imagery by applying the aforementioned codes; only the first 60 thermal images were taken as the input for further data processing.
The matrices X were reduced to d-dimension by using the ISOMAP, LLE, and LTSA algorithms. In order to visualize the dimension reduction results, the partial least squares regression (PLSR) algorithm was chosen to approximate the mapping process from X to Y.

4.1. Experimental Results

4.1.1. ISOMAP

While using ISOMAP (Figure 11), the following parameters were accepted: d = 6 and neighbours = 4. After ISOMAP data processing, all six indicator vector (IV) images contained varying degrees of the inhomogeneous background, as well as noise, which obscured information on reinforcing grid. Some grid features can be observed in the IV2 image, but it is difficult to clearly detect the meridians and latitudes of the grid because of the noise. The IV3–IV5 images contained both the severe inhomogeneous background and noise, thus making demarcation at particular specimen sections difficult.

4.1.2. LLE

The parameters of the LLE algorithm were (Figure 12): d = 6 and neighbours = 6. Observing the six IV images of the LLE results, it can be said that the inhomogeneous background (caused by the inhomogeneity of the specimens) is attenuated in the images when compared to the ISOMAP results. In particular, the grid feature is more pronounced in the images of IV4–IV6 as compared to ISOMAP. This may be attributed to the fact that LLE is more sensitive to inhomogeneous information, as it retains the local intrinsic structure of the data unlike the ISOMAP method, which learns the global manifold of the data.

4.1.3. LTSA

The LTSA parameters were set as follows (Figure 13): d = 6 and neighbours = 9. Similarly to the previous two approaches, the 6 IVs of LTSA still had difficulty in identifying the grid lines. Compared with the ISOMAP and LLE results, the LTSA results had more severe inhomogeneous background and associated noise, while the grid features were clearer. In particular, the upper part of the IV2 image contained the weak inhomogeneous background, and the grid lines were more visible in regard to other images.
In summary, none of the results can clearly display the grid shown in area 1 (Figure 1a) because of noise. Therefore, it was necessary to use edge detection algorithms to further analyse the thermographic results.
Here, the authors select indication vector (IV) 2 in LTSA as a case, and use the Canny, Sobel and Laplacian operators to extract signal gradients (Figure 14).
Applying the edge detection algorithm resulted in a better display of the grid outline in the area 1 (Figure 1a). However, the Laplacian operator could only extract a small amount of grid information in the area 2 (Figure 1a), and the Sobel operator also showed a poor denoising ability. In contrast, the Canny operator was able not only to extract the grid information in the areas 1 and 2, but also to effectively remove the noise (Figure 14).
When implementing the Canny operator, some pre-treatments should be conducted to achieve better results. One method is to save the IV2 image (Figure 13) with a resolution of 300 dpi, and another method is to enlarge the IV2 image (Figure 13) directly where each pixel in the original image corresponds to a region of 8×8 in the enlarged image. The respective results are shown in Figure 15.
It is visually seen that the second method produced better results. It was compared with the ground truth image manually obtained by means of the Photoshop® software through the replicating area 3 above the areas 2 and 1 (Figure 1a). The ground truth image is shown on the left side of Figure 16. However, it is still difficult to give a quantitative comparison between the Canny and ground truth results, since the numbers and positions of grids did not perfectly fit. This means that the ground truth should be obtained in a different way in a future work, perhaps by removing the layers above the grid in the areas 1 and 2 in order to acquire a photograph covering the entire area.

5. Modelling via ThermoCalc-3D® Software

As shown in Figure 17a, the plane Section 1-1 (Figure 1b) was modelled in order to analyse the possibility to detect the fiberglass grid located beneath scattered hemp fibres and embedded into the layers of cement at the depth of ~1.5 mm from the sample surface. At such depth, air-filled defects can easily be detected due to their high thermal contrast in regard to the host material. The thermal properties of fiberglass were comparable to those of hemp fibres, and the detection of the grid seemed to be possible due its thermal contrast in respect to cement mortar. In the 3D model of Figure 17a (material thermal properties are given in the caption of the figure), the specimen was heated for 5 s using a uniform heat flux, and the purpose of the modelling was to evaluate the dimensionless temperature contrast C = DT/Tnd, where DT = TdTnd, with the subscripts “d” and “nd” specifying temperatures in defect and non-defect areas, as shown in Figure 17a [32].
The images of the contrast C are displayed in Figure 17b, and its time evolution, in Figure 17c. The contrast reaches the maximum within the heating period; then, it decays and finally experiences one more maxima at the end of heating that is a typical behaviour of the contrast in active TNDT (note that the differential signal DT is characterized by a single maxima). At longer times, the contrast becomes negative because of the higher heat capacity of fiberglass with respect to that of the cement mortar, which is the abovementioned signal inversion, but its amplitude is low; the reader may compare the first and last images in Figure 17b.
The maximum contrast is about C = 0.012, or 1.2%. This is a fairly low value to ensure reliable detection of the grid; for example, if the sample is heated up to 10 °C above the ambient, the temperature signal DT produced by the grid will be only 0.12 °C. In the ideal case of no noise, the grid pattern can be underlined, for example, by the thermal tomogram shown in Figure 17d; it is worth noting that the grid pattern in the tomogram is distorted because of heat diffusion. In general, any material can be characterized by the surface noise contrast C n = σ n / T n d = , which includes both the additive and multiplicative noise components; here, σ n is the standard deviation of noise within a chosen reference area, and T n d = is the mean temperature in the reference area [33]. Noise contrast varies in time and depends on an area chosen for calculating σ n . For hemp fibers, the value of C n is presented in Figure 17c with the dotted line. The noise is higher within the heating stage because of reflected radiation with a fixed value of about 1.2% within the cooling stage. By comparing the contrasts produced by the grid and the surface clutter (Figure 17c), one may conclude that the fiberglass grid cannot be detected on the noisy background (in some cases, a regular structure of a defect pattern may be helpful for identifying the corresponding structures).

6. Processing Experimental Thermograms: The Second Session

By applying the Principal Component Analysis (PCA) to the experimental results, one can clearly observe the non-uniform structure of the mortar layers (Figure 18), which appears due to varying thickness of the mortar, as well as some inter-layer detachments occurring on the boundaries of the fibre glass grid.
The analysis of the area 1 (Figure 1a and Figure 18c) shows a vague pattern of the grid structure underneath the hemp fibres, which act as a thermal insulator. As mentioned above, in one-sided TNDT, some depth-resolved results can be obtained by using the technique of dynamic thermal tomography [34]. The basic concept of this technique involves the analysis of time delays of differential temperature signals between a defect and sound areas. Each class of structural defects can be described by a particular calibration curve allowing the evaluation of a defect depth by the corresponding time delays. Material diffusivity and thickness also affect calibration results. Some thermal tomograms of the tested specimen are shown in Figure 19 to reveal no grid pattern, as predicted by the noise analysis. However, thermal tomography illustrates variations of mortar thickness and some peculiarities of inter-layer boundaries.

7. Discussion and Conclusions

New building enclosure materials, its assemblies, and systems used for minimizing and/or regulating space heating and cooling modes are today at the focus of several research projects and gaining more interest from many sectors.
In this work, a specific multi-layer composite (useful to limit heat losses from buildings) was analysed numerically and experimentally for integrity. In particular, Comsol Multiphysics® and ThermoCalc-3D® software and a couple of active IRT experiments were implemented to detect a grid in glass fibres positioned under an external layer made of hemp fibres. The numerical part revealed some matches, such as the evaluation of the thermal/signal inversion and exposed challenges in detecting the grid on the noisy background. In general, the ThermoCalc-3D software® can be considered more user-friendly in regard to the modelling performed on the Comsol platform, even if the latter approach was optimized to reduce computational cost.
In order to detect the reinforcing grid, some new data processing approaches were proposed and demonstrated. Namely, the IV2 and LTSA techniques were complemented by the use of the Canny, Sobel, and Laplacian operators. The best results were provided by the Canny algorithm for the presented case.
The one-sided TNDT (usually used in building and material inspections) [35,36] was also applied in this study, specifically in the second experimental session. Data processing was performed using the techniques of PCA and thermal tomography. No grid pattern, as predicted by the noise analysis, was visible by thermal tomography, although the varying thickness of the mortar and some peculiarities of inter-layer boundaries were revealed. The second component of PCA slightly revealed the presence of the sub-surface grid.
No quantitative analysis was possible at this stage of research as the extraction of an accurate ground truth image inherent to the grid pattern still requires an objective strategy to allow quantitative assessment.
In conclusion, we have demonstrated that: (a) the use of advanced software, such as Comsol Multiphysics, is useful to precisely simulate a multi-layered specimen within an accepted heat conduction model; (b) a thin external layer of hemp fibres ensures a significant shift of thermal phase in construction materials; (c) the use of ThermoCalc-3D® software allowed us to link the maximum thermal contrast with the noisy background and understand why a sub-surface fiberglass grid cannot be detected by using only optical heating; (d) among some advanced data processing approaches used, such as ISOMAP, LLE, and LTSA, the LTSA algorithm demonstrated better results; (e) the thermal tomography technique was able to reveal detachments between specimen layers better than other processing techniques; and (f) the LTSA algorithm provided the best input image for the subsequent segmentation procedure by means of the Canny operator.

Author Contributions

Conceptualization, S.P., Y.Y., V.V. and S.S.; methodology, K.L., A.O.C., M.O. and Y.L.; software, K.L., S.P., V.V. and S.S.; validation, Y.Y., M.O., V.V., Y.L. and S.S.; formal analysis, K.L., A.O.C., M.O. and Y.L.; investigation, S.P., A.O.C., M.O. and S.S.; resources, Y.Y., M.O., V.V., Y.L. and S.S.; data curation, K.L., S.P., A.O.C., V.V., Y.L. and S.S.; writing—original draft preparation, K.L., S.P., A.O.C., Y.Y., V.V., Y.L. and S.S.; writing—review and editing, K.L., S.P., A.O.C., Y.Y., M.O., V.V., Y.L. and S.S.; visualization, K.L., S.P., A.O.C., Y.Y., M.O., V.V., Y.L. and S.S.; supervision, Y.Y., M.O., V.V., Y.L. and S.S.; project administration, Y.Y., M.O., V.V., Y.L. and S.S.; funding acquisition, Y.Y., Y.L. and S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China (Grant No. 62022073). Yao was supported in part by Ministry of Science and Technology, ROC under Grant MOST 108-2221-E-007-068-MY3.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

2Dtwo-dimensional
3Dthree-dimensional
AEacoustic emission
CFRPcarbon fibre-reinforced polymers
DICdigital image correlation
FBHflat-bottom holes
GFRPglass fibre-reinforced polymer
IRTinfrared (IR) thermography
ISOMAPisometric feature mapping
IVindication vector
LCSlaminated composite structures
LLElocal linear embedding
LTLlocalized thickness loss
LTSAlocal tangent space alignment
NDTnon-destructive testing
PCAprincipal component analysis
PLSRpartial least squares regression
SARsynthetic aperture radar
TNDTthermal non-destructive testing

References

  1. Bisegna, F.; Ambrosini, D.; Paoletti, D.; Sfarra, S.; Gugliermetti, F. A qualitative method for combining thermal imprints to emerging weak points of ancient wall structures by passive infrared thermography—A case study. J. Cult. Herit. 2014, 15, 199–202. [Google Scholar] [CrossRef]
  2. Kalyanavalli, V.; Abilasha Ramadhas, T.K.; Sastikumar, D. Long pulse thermography investigations of basalt fiber reinforced composite. NDT E Int. 2018, 100, 84–91. [Google Scholar] [CrossRef]
  3. Fernandes, H.; Ibarra-Castanedo, C.; Zhang, H.; Maldague, X. Thermographic non-destructive evaluation of carbon fiber-reinforced polymer plates after tensile testing. J. Nondestruct. Eval. 2015, 34, 35. [Google Scholar] [CrossRef] [Green Version]
  4. Montanini, R.; Freni, F. Non-destructive evaluation of thick glass fiber-reinforced composites by means of optically excited lock-in thermography. Compos. Pt. A Appl. Sci. Manuf. 2012, 43, 2075–2082. [Google Scholar] [CrossRef]
  5. Zeng, Z.; Tao, N.; Feng, L.; Li, Y.; Zhang, C. Relative thermal contrast analysis in the inspection of wind turbine blades using pulsed thermography. Adv. Mater. Res. 2011, 301–303, 591–596. [Google Scholar] [CrossRef]
  6. Kim, J.W.; Kim, H.S.; Lee, D.G. Measurement of tensile strength and infrared thermography on unidirectional GFRP. Mater. Res. Bull. 2014, 58, 35–38. [Google Scholar] [CrossRef]
  7. Angiuli, R.; Corvaglia, P.; Largo, A.; Coricciati, A. Defect identification and acceptance of FRP and FRCM masonry reinforcement by infrared thermography survey. Key Eng. Mater. 2015, 624, 80–87. [Google Scholar] [CrossRef]
  8. Pitarresi, G.; Scalici, T.; Catalanotti, G. Infrared thermography assisted evaluation of static and fatigue Mode II fracture toughness in FRP composites. Compos. Struct. 2019, 226, 11220. [Google Scholar] [CrossRef]
  9. Xu, C.; Zhang, W.; Wu, C.; Xie, J.; Yin, X.; Chen, G. An improved method of eddy current pulsed thermography to detect subsurface defects in glass fiber reinforced polymer composites. Compos. Struct. 2020, 242, 112145. [Google Scholar] [CrossRef]
  10. Rani, A.; Mulaveesala, R. Depth resolved pulse compression favourable frequency modulated thermal wave imaging for quantitative characterization of glass fibre reinforced polymer. Infrared Phys. Technol. 2020, 110, 103441. [Google Scholar] [CrossRef]
  11. Karger-Kocsis, J.; Harmia, T.; Czigány, T. Comparison of the fracture and failure behavior of polypropylene composites reinforced by long glass fibers and by glass mats. Compos. Sci. Technol. 1995, 54, 287–298. [Google Scholar] [CrossRef]
  12. Karger-Kocsis, J.; Czigány, T. Effects of interphase on the fracture and failure behavior of knitted fabric reinforced composites produced from commingled GF/PP yarn. Compos. Pt. A Appl. Sci. Manuf. 1998, 29, 1319–1330. [Google Scholar] [CrossRef]
  13. Ibarra-Castanedo, C.; Avdelidis, N.P.; Grinzato, E.G.; Bison, P.G.; Marinetti, S.; Plescanu, C.C.; Bendada, A.; Maldague, X.P. Delamination detection and impact damage assessment of GLARE by active thermography. Int. J. Mater. Prod. Technol. 2011, 41, 5–16. [Google Scholar] [CrossRef]
  14. Katunin, A.; Dragan, K.; Dziendzikowski, M. Damage identification in aircraft composite structures: A case study using various non-destructive testing techniques. Compos. Struct. 2015, 127, 1–9. [Google Scholar] [CrossRef]
  15. Palumbo, D.; Tamborrino, R.; Galietti, U.; Aversa, P.; Tatì, A.; Luprano, V.A.M. Ultrasonic analysis and lock-in thermography for debonding evaluation of composite adhesive joints. NDT E Int. 2016, 78, 1–9. [Google Scholar] [CrossRef]
  16. Crupi, V.; Kara, E.; Epasto, G.; Guglielmino, E.; Aykul, H. Theoretical and experimental analysis for the impact response of glass fibre reinforced aluminium honeycomb sandwiches. J. Sandw. Struct. Mater. 2018, 20, 42–69. [Google Scholar] [CrossRef]
  17. Quattrocchi, A.; Freni, F.; Montanini, R. Comparison between air-coupled ultrasonic testing and active thermography for defect identification in composite materials. Nondestruct. Test. Eval. 2021, 36, 97–112. [Google Scholar] [CrossRef]
  18. Ghorbel, A.; Saitier, N.; Dhiab, A. Investigation of damage evolution in short glass fibers reinforced polyamide 6,6 under tensile loading using infrared thermography. Proc. Eng. 2011, 10, 2123–2128. [Google Scholar] [CrossRef] [Green Version]
  19. Omer, M.; Sultan, Q. Post processing of thermographic NDT data using Dirac comb. Exp. Tech. 2017, 41, 179–190. [Google Scholar] [CrossRef]
  20. Zhao, S.; Guo, X.; Liu, Y. Infrared modulated thermography based on complex Morlet wavelet phase. Infrared Laser Eng. 2018, 47, 0904002. [Google Scholar] [CrossRef]
  21. Yanjie, W.; Han, Y.; Zhang, S.; Zhang, D. Defects identification in fiber reinforced plastics using long pulse thermography. In Proceedings of the SPIE—The International Society for Optical Engineering—7th International Conference on Optical and Photonic Engineering, Baltimore, Thailand, 16 October 2019; Volume 11205, p. 1120513. [Google Scholar] [CrossRef]
  22. Vijaya Lakshmi, A.; Gopitilak, V.; Muzammil Parvez, M.; Subhani, S.K.; Ghali, V.S. Artificial neural networks based quantitative evaluation of subsurface anomalies in quadratic frequency modulated thermal wave imaging. Infrared Phys. Technol. 2019, 97, 108–115. [Google Scholar] [CrossRef]
  23. Sutthaweekul, R.; Tian, G.; Wang, Z.; Ciampa, F. Microwave open-ended waveguide for detection and characterisation of FBHs in coated GFRP pipes. Compos. Struct. 2019, 225, 111080. [Google Scholar] [CrossRef]
  24. Pracht, M.; Swiderski, W. Analysis of the possibility of non-destructive testing to detect in multi-layered composites reinforced fibers by optical IR thermography. Compos. Struct. 2019, 213, 204–208. [Google Scholar] [CrossRef]
  25. Bernegger, R.; Altenburg, S.J.; Maierhofer, C. Quantification of delaminations in demitransparent solids using pulsed thermography and mathematical 1D models. Int. J. Thermophys. 2020, 41, 67. [Google Scholar] [CrossRef] [Green Version]
  26. Hu, J.; Li, Y.; Tan, J.; Li, W.; Chen, Z. Characterization and imaging of localized thickness loss in GFRP with Ka-band microwave open-ended waveguides. Stud. Appl. Electromagn. Mech. 2020, 45, 218–223. [Google Scholar] [CrossRef]
  27. Ridolfi, S.; Crescenzi, S.; Zeli, F.; Perilli, S.; Sfarra, S. Evaluation of the heat changes in an ancient church because of restoration works: A microclimatic study supported by thermal images. Indoor Built Environ. 2021. Advance online publication. [Google Scholar] [CrossRef]
  28. Sfarra, S.; Perilli, S.; Guerrini, M.; Bisegna, F.; Chen, T.; Ambrosini, D. On the use of phase change materials applied on cork-coconut-cork panels: A thermophysical point of view concerning the beneficial effect in terms of insulation properties. J. Therm. Anal. Calorim. 2019, 138, 4061–4090. [Google Scholar] [CrossRef]
  29. Sfarra, S.; Perilli, S.; Ambrosini, D.; Paoletti, D.; Nardi, I.; de Rubeis, T.; Santulli, C. A proposal of a new material for greenhouses on the basis of numerical, optical, thermal and mechanical approaches. Constr. Build. Mater. 2017, 155, 332–347. [Google Scholar] [CrossRef]
  30. Comsol. Performing a mesh refinement study. Available online: https://www.comsol.ch/support/knowledgebase/1261 (accessed on 16 March 2020).
  31. Crnjac, P.; Škerget, L.; Ravnik, J.; Hriberšek, M. Implementation of the Rosseland and the P1 radiation models in the system of Navier-stokes equations with the boundary element method. Int. J. Comput. Methods Exp. Meas. 2017, 5, 348–358. [Google Scholar] [CrossRef] [Green Version]
  32. Vavilov, V.; Marinetti, S.; Pan, Y.; Chulkov, A. Detecting water ingress in aviation honeycomb panels: Qualitative and quantitative aspects. Polym. Test. 2016, 54, 270–280. [Google Scholar] [CrossRef]
  33. Vavilov, V.P. Noise-limited thermal/infrared nondestructive testing. NDT E Int. 2014, 61, 16–23. [Google Scholar] [CrossRef]
  34. Vavilov, V.P.; Shirayev, V.V.; Kuimova, M.V. Time- and phase-domain thermal tomography of composites. Photonics 2018, 5, 31. [Google Scholar] [CrossRef] [Green Version]
  35. Vavilov, V.P.; Chulkov, A.O.; Derusova, D.A.; Pan, Y. Thermal NDT research at Tomsk Polytechnic University. Quant. InfraRed Thermogr. J. 2016, 13, 128–143. [Google Scholar] [CrossRef]
  36. Vavilov, V. Thermal NDT: Historical milestones, stae-of-the-art and trends. Quant. InfraRed Thermogr J. 2014, 11, 66–83. [Google Scholar] [CrossRef]
Figure 1. (a) Test specimen (front view); (b) test specimen (lateral view) along with the plane section linked to Figure 17a; (c) numerical model involving surface thermal loading; (d) the specimen built in Comsol Multiphysics®.
Figure 1. (a) Test specimen (front view); (b) test specimen (lateral view) along with the plane section linked to Figure 17a; (c) numerical model involving surface thermal loading; (d) the specimen built in Comsol Multiphysics®.
Infrastructures 06 00131 g001aInfrastructures 06 00131 g001b
Figure 2. Thermographic testing: (a) first scheme (1st test), and (b) second scheme (2nd test).
Figure 2. Thermographic testing: (a) first scheme (1st test), and (b) second scheme (2nd test).
Infrastructures 06 00131 g002
Figure 3. The first model allows us to evaluate the thermal distribution induced on the specimen: (a) image of the mesh; (b) image of the thermal distribution at the beginning of heating; (c) thermal distribution at 10 s.
Figure 3. The first model allows us to evaluate the thermal distribution induced on the specimen: (a) image of the mesh; (b) image of the thermal distribution at the beginning of heating; (c) thermal distribution at 10 s.
Infrastructures 06 00131 g003
Figure 4. The second model mesh built in Comsol Multiphysics®.
Figure 4. The second model mesh built in Comsol Multiphysics®.
Infrastructures 06 00131 g004
Figure 5. Variations of specimen temperature distribution during the heating phase: (a) 10 s of heating, (b) 20s, (c) 30 s, (d) 40 s, (e) 50 s, and (f) 60 s.
Figure 5. Variations of specimen temperature distribution during the heating phase: (a) 10 s of heating, (b) 20s, (c) 30 s, (d) 40 s, (e) 50 s, and (f) 60 s.
Infrastructures 06 00131 g005
Figure 6. Variations of specimen temperature distribution during the cooling phase (total time in brackets): (a) 10 s of cooling (70 s), (b) 20 s (80 s), (c) 30 s (90 s), (d) 60 s (120 s), (e) 90 s (150 s), (f) 120 s (180 s), (g) 150 s (210 s), (h) 180 s (240 s).
Figure 6. Variations of specimen temperature distribution during the cooling phase (total time in brackets): (a) 10 s of cooling (70 s), (b) 20 s (80 s), (c) 30 s (90 s), (d) 60 s (120 s), (e) 90 s (150 s), (f) 120 s (180 s), (g) 150 s (210 s), (h) 180 s (240 s).
Infrastructures 06 00131 g006aInfrastructures 06 00131 g006bInfrastructures 06 00131 g006cInfrastructures 06 00131 g006d
Figure 7. Determining in-depth temperature distributions (see red segments): (a) cement mortar (upper position), (b) cement milk, (c) cement mortar (lower position), (d) hemp fibres.
Figure 7. Determining in-depth temperature distributions (see red segments): (a) cement mortar (upper position), (b) cement milk, (c) cement mortar (lower position), (d) hemp fibres.
Infrastructures 06 00131 g007
Figure 8. In-depth temperature profiles in the heating phase: (a) 10 s of heating, (b) 20 s, (c) 30 s, (d) 40 s, (e) 50 s, (f) 60 s.
Figure 8. In-depth temperature profiles in the heating phase: (a) 10 s of heating, (b) 20 s, (c) 30 s, (d) 40 s, (e) 50 s, (f) 60 s.
Infrastructures 06 00131 g008
Figure 9. In-depth temperature profiles in the cooling phase (total time in brackets): (a) 10 s of cooling (70 s), (b) 30 s (90 s), (c) 60 s (120 s), (d) 90 s (150 s), (e) 120 s (180 s), (f) 180 s (240 s).
Figure 9. In-depth temperature profiles in the cooling phase (total time in brackets): (a) 10 s of cooling (70 s), (b) 30 s (90 s), (c) 60 s (120 s), (d) 90 s (150 s), (e) 120 s (180 s), (f) 180 s (240 s).
Infrastructures 06 00131 g009
Figure 10. Unfolding the 3D matrix into 2D.
Figure 10. Unfolding the 3D matrix into 2D.
Infrastructures 06 00131 g010
Figure 11. ISOMAP analysis results.
Figure 11. ISOMAP analysis results.
Infrastructures 06 00131 g011
Figure 12. LLE analysis results.
Figure 12. LLE analysis results.
Infrastructures 06 00131 g012
Figure 13. LTSA analysis results.
Figure 13. LTSA analysis results.
Infrastructures 06 00131 g013
Figure 14. Edge detection results: (a) original image, (b) Canny result, (c) Sobel result, (d) Laplacian result.
Figure 14. Edge detection results: (a) original image, (b) Canny result, (c) Sobel result, (d) Laplacian result.
Infrastructures 06 00131 g014
Figure 15. Canny results based on some pre-treatments of IV2 image: (a) saving IV2 image with 300 dpi resolution, (b) saving enlarged (8 × 8 area) IV2 image.
Figure 15. Canny results based on some pre-treatments of IV2 image: (a) saving IV2 image with 300 dpi resolution, (b) saving enlarged (8 × 8 area) IV2 image.
Infrastructures 06 00131 g015
Figure 16. Comparison between the ground truth (left) and the Canny result (right).
Figure 16. Comparison between the ground truth (left) and the Canny result (right).
Infrastructures 06 00131 g016
Figure 17. Modelling active TNDT of the specimen shown in Figure 1a,b (material thermal properties: Scattered hemp fibres—K = 0.038 W m−1 K−1, C = 1700 J kg−1 K−1, ρ = 25 kgm−3; Cement mortar—K = 1.73 W m−1 K−1, C = 210 J kg−1 K−1, ρ = 900 kgm−3; Fiberglass—K = 0.035 W m−1 K−1, C=1030 J kg−1 K−1, ρ = 21 kgm−3; Cement milk—K = 1.4 W m−1 K−1, C=0.87 J kg−1 K−1, ρ = 1540 kgm−3; Expanded Styrofoam—K = 0.03 W m−1 K−1, C = 1450 J kg−1 K−1, ρ = 30 kgm−3). (a) 3D numerical model (ThermoCalc-3D software); (b) images of temperature contrast at three different times; (c) contrast evolution; (d) thermal tomogram of the grid layer.
Figure 17. Modelling active TNDT of the specimen shown in Figure 1a,b (material thermal properties: Scattered hemp fibres—K = 0.038 W m−1 K−1, C = 1700 J kg−1 K−1, ρ = 25 kgm−3; Cement mortar—K = 1.73 W m−1 K−1, C = 210 J kg−1 K−1, ρ = 900 kgm−3; Fiberglass—K = 0.035 W m−1 K−1, C=1030 J kg−1 K−1, ρ = 21 kgm−3; Cement milk—K = 1.4 W m−1 K−1, C=0.87 J kg−1 K−1, ρ = 1540 kgm−3; Expanded Styrofoam—K = 0.03 W m−1 K−1, C = 1450 J kg−1 K−1, ρ = 30 kgm−3). (a) 3D numerical model (ThermoCalc-3D software); (b) images of temperature contrast at three different times; (c) contrast evolution; (d) thermal tomogram of the grid layer.
Infrastructures 06 00131 g017aInfrastructures 06 00131 g017b
Figure 18. PCA-processing results. (a) Raw thermogram at 7 s; (b) 1st principal component; (c) area 1 (Figure 1a), 2nd component; (d) area 2 (Figure 1a), 3rd component; (e) area 3 (Figure 1a), 3rd component; (f) area 4 (Figure 1a), 3rd component.
Figure 18. PCA-processing results. (a) Raw thermogram at 7 s; (b) 1st principal component; (c) area 1 (Figure 1a), 2nd component; (d) area 2 (Figure 1a), 3rd component; (e) area 3 (Figure 1a), 3rd component; (f) area 4 (Figure 1a), 3rd component.
Infrastructures 06 00131 g018
Figure 19. Thermal tomograms of specimen “slices”: (a) layer coordinates: 0–2 mm, (b) 0–0.3 mm, (c) 0.3–0.6 mm, (d) 0.6–0.9 mm, (e) 0.9–1.2 mm.
Figure 19. Thermal tomograms of specimen “slices”: (a) layer coordinates: 0–2 mm, (b) 0–0.3 mm, (c) 0.3–0.6 mm, (d) 0.6–0.9 mm, (e) 0.9–1.2 mm.
Infrastructures 06 00131 g019
Table 1. Dimensional and thermophysical properties of the materials under analysis.
Table 1. Dimensional and thermophysical properties of the materials under analysis.
n.MaterialThickness [m]Conductivity [W/mK]Density [kg/m3]Specific Heat [J/kgK]Emissivity
ε
Porosity
%Vol.
1Scattered Hemp Fibres0.0010.0382517000.90.120
2Cement Mortar0.00051.739000.210.540.22
3Fiberglass0.00020.0352110300.750.44
4Cement Milk0.00011.415400.870.920.147
5Expanded Styrofoam0.02860.033014500.60.497
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, K.; Perilli, S.; Chulkov, A.O.; Yao, Y.; Omar, M.; Vavilov, V.; Liu, Y.; Sfarra, S. Defining the Thermal Features of Sub-Surface Reinforcing Fibres in Non-Polluting Thermo–Acoustic Insulating Panels: A Numerical–Thermographic–Segmentation Approach. Infrastructures 2021, 6, 131. https://0-doi-org.brum.beds.ac.uk/10.3390/infrastructures6090131

AMA Style

Liu K, Perilli S, Chulkov AO, Yao Y, Omar M, Vavilov V, Liu Y, Sfarra S. Defining the Thermal Features of Sub-Surface Reinforcing Fibres in Non-Polluting Thermo–Acoustic Insulating Panels: A Numerical–Thermographic–Segmentation Approach. Infrastructures. 2021; 6(9):131. https://0-doi-org.brum.beds.ac.uk/10.3390/infrastructures6090131

Chicago/Turabian Style

Liu, Kaixin, Stefano Perilli, Arsenii O. Chulkov, Yuan Yao, Mohammed Omar, Vladimir Vavilov, Yi Liu, and Stefano Sfarra. 2021. "Defining the Thermal Features of Sub-Surface Reinforcing Fibres in Non-Polluting Thermo–Acoustic Insulating Panels: A Numerical–Thermographic–Segmentation Approach" Infrastructures 6, no. 9: 131. https://0-doi-org.brum.beds.ac.uk/10.3390/infrastructures6090131

Article Metrics

Back to TopTop