Investigation of CFRP Concrete Reinforcement Ageing Using Microwave Infrared Thermography Method (MIRT)

1. Introduction

The utilization of carbon-fiber-reinforced polymer (CFRP) in reinforcing concrete is increasingly favored for constructing resilient structures [1,2,3,4,5]. The process of CFRP reinforcement involves the application of a layer of carbon matting bound with epoxy resin onto the cement base. The primary focus in this technique is the quality of adhesion between cement and CFRP [6,7]. Insufficient adhesive bonding or delamination weakens the reinforcement in specific areas, making constructed structures susceptible to damage [8]. Hence, it is necessary to regulate the aforementioned structures in order to evaluate the bond between concrete and composite materials. Different non-destructive testing techniques can be used to assess the quality of the cement-composite connection. Possible methods include: radiography, ultrasonic imaging, and thermo-vision and electromagnetic techniques [9,10,11,12,13,14,15]. Given the substantial size of the buildings under scrutiny, it is preferable to employ techniques that enable comprehensive assessment of extensive areas. Active thermography is a rapid and non-contact method for analyzing large structures, which makes it valuable for field study. In this study, the active thermography with microwave excitation technique for quality assessment of the CFRP concrete reinforcement was employed. Active thermography utilizes an external energy source to create a temperature contrast in the specimen [16]. Numerous global research institutions investigate thermal imaging techniques for the investigation of large-scale objects, such as complex structures made of concrete. Passive thermography utilizes ambient solar radiation instead of an external heat source. This approach entails extended periods of monitoring and is influenced by daily changes and meteorological factors [16,17]. A variety of active thermography approaches employ halogen lamps, laser, flash lights, infrared radiators, and hot air as heat sources [18,19,20,21,22,23,24]. In this case, the outcomes are no longer contingent on the external condition, but the primary concern is the extended duration required for heating because of the size of the structure. It is crucial to recall that conduction is the process by which heat is transferred from a heated surface to the inside of an object. These stimulation methods are highly specialized, resulting in longer examination durations. Vibrothermography [25,26,27] employs mechanical energy and dissipation mechanisms to identify flaws. The acoustic wave stimulation propagates through the tested structure and attenuates in areas of damage and defects. While acoustic approaches have proven to be effective, they need direct physical contact between the exciter and the test object, which can pose challenges in some circumstances. These methods are most effective for identifying significant defects such as fractures. The utilization of microwave irradiation is proposed as the method of choice for heating in this study. The literature has utilized microwave heating to expedite the process of concrete curing, disinfect cement, and facilitate drilling or melting of concrete [28,29,30,31,32,33,34,35,36,37]. Microwaves have the potential to provide the necessary energy for infrared thermography, as indicated by several studies [20,38,39,40,41,42,43,44]. The primary advantage of microwave heating lies in its volumetric nature. This approach is efficient as it rapidly heats up a precise volume of the specimen at a given moment. Conversely, the heating ratio of this alternative energy source is influenced by both thermal and electrical properties, specifically the dielectric constant [38,45]. In case of conductive materials such as CFRP, microwave heating is limited by a shallower depth of penetration. Following the application of heat, a thermovision camera is utilized to visually assess the specimen. Delamination and adhesive loss can result in uneven heating of the sample, which is visible in the registered temperature distribution [20,46,47,48,49,50].

In this article the findings of our research on the study of a concrete samples reinforced with a carbon composite material are presented. The experiments involved submitting the material to microwave irradiation at a constant power level. The objective of the analysis was to validate the feasibility of utilizing the given methodology for identifying flaws in a carbon reinforcement in the form of lack of adhesive between the reinforcement and the concrete. First results of this analysis were presented in our previous articles [37,48]. An essential element of the current study involves utilizing a sample that underwent an extended aging period of 10 years. This prolonged duration facilitated the development of natural deterioration in the reinforcement, characterized by the fraying of the material at the edges and the chipping of the epoxy. During the initial phase of the study, we created numerical models with the aim of replicating material deterioration associated with aging. The models were utilized for comparison research and to ascertain the influence of this particular type of flaws on the efficacy of the approach. The numerical analysis presented in the study was validated through experimental verification, and the numerical and experimental findings were subsequently compared.

2. Materials and Methods

In this section, we present the numerical models and the utilized experimental setups. The numerical analysis was performed using Finite Element Method (FEM) in a COMSOL computing environment. In the subsequent section all the elements of the designed numerical models are presented.

2.1. Numerical Models

The numerical models were prepared in the commercial computing environment - COMSOL, in particular using Heat Transfer (HT) and Electromagnetic Waves (EMW) modules to fully analyze the phenomena of microwave heating. The propagation of an electromagnetic wave is defined as follows [52]:

$${\nabla }^2\mathbf{E}+\left(\mu ϵ{\omega }^2-i\mu \sigma \omega \right)\mathbf{E}=0,$$

(1)

where $\mu$ denotes the magnetic permeability, $ϵ$is the electric permittivity, $\sigma$ represents the electric conductivity, $\omega$ is the angular frequency, and $\mathbf{E}$ is the electric field vector.

The heat flow is characterized here by the equation [16]:

$$\rho C_p\left(\frac{\partial T}{\partial t}+\nabla T\right)+\nabla \mathit{q}=\mathrm{Q},$$

(2)

where $\rho$ is the material density, $C_p$denotes the material heat capacity, $q$ is the heat flux associated with the convection phenomenon, T is the temperature, and $\mathrm{Q}$ denotes the external heat source.

And the external heat source is equal to the resistive heat generated by the electromagnetic field, which is described as follows [52]:

$$\mathrm{Q}=1/2\mathrm{Re}\left(\right(\mathsf{\sigma }-\mathrm{j}\mathsf{\omega }\mathsf{ϵ}\left)\mathrm{E}\text{⋅}{\mathbf{E}}^{\mathbf{*}}\right),$$

(3)

where Re indicates the real part of the value. Figure 1 presents the geometry of the modelled horn antenna, with a rectangular waveguide. An electromagnetic wave was produced through the application of the port boundary condition on a specific wall within the waveguide. Here, the assumed mode was TE₁₀, the frequency was equal to 2.45 GHz and the total power was set to 1000 W. The energy was transmitted in the continuous manner.

To thoroughly examine the functioning of the tested system, the suggested model includes four distinct variations of the sample. The sample being modeled is a concrete block of 400 x 400 x 45 mm, which has been reinforced with a layer of CFRP measuring 220 x 400 x 1 mm. A 1 mm epoxy adhesive bond was simulated between the composite layer and the concrete. All material parameters important for the modeled objects have been collected in a table. Four variants of the sample were prepared (Figure 2):

• a sample without any defects (Figure 2 (a)) – (No Defect Sample - NDS);
• a sample with a defect in the form of lack of adhesive (air void in the epoxy glue layer with dimensions of 100 x 100 mm located centrally) (Figure 2 (b)) – (Central Defect Sample – CDS)
• a sample with defects simulating the aging of the material, located on the edges of the composite reinforcement layer (Figure 2 (c)) – (Ageing defects Sample – ADS)
• a sample with defects simulating aging and with a centrally located defect in the form of lack of adhesive (Figure 2 (d)) – (Central and Ageing Defects Sample – CADS)

The main challenge in designing the numerical model was to model defects resulting from the aging of the composite material. Such defects first appear on the edges of the reinforcement layer and are caused by crumbling of the glue and fraying of the carbon fibers. A unique approach to the problem was suggested in order to replicate a similar result in a numerical model. The aging changes were represented by a stochastic structure consisting of many cylindrical air spaces that spanned both the epoxy adhesive layer and the adjacent carbon composite layer. This method replicates both the deterioration of the adhesive and the weakening of the carbon fibers at the same time. Building such a structure in the COMSOL environment is possible thanks to the automation of the geometry creation process in Application Builder. The random structure of defects (Figure 3) was created with certain assumptions:

• the radius of each defect varies randomly between 0.1 and 1 cm. The height of the defect is either 0.01 cm (for 30% of defects) or 0.015 cm (for the remaining 70% of defects). In the first case, the entire defect is within the adhesive layer, while in the latter, it partially extends into the composite layer (Figure 3 (a));
• the defects were located exclusively at the periphery of the composite reinforcement, with their centers randomly distributed along a 1 cm margin encompassing the entire composite (Figure 3 (b));
• the total count of cylindrical defects is 100;
• overlapping defects are allowed, but their logical sum will be taken into account.

A whole structure is presented In Figure 3 (c) with the magnification of the chosen region.

Figure 4 (a) presents the whole geometry designed for one of the scenarios, specifically a sample exhibiting a defect in the form of non-adherence, together with defects presented on the borders of the composite layer (CADS). The system accurately replicates the physical laboratory setup, considering the precise location of the sample in relation to the antenna. An adequate angle was chosen to rotate the sample in relation to the antenna, allowing for the surface to be observed using a thermal imaging camera. The angular placement of the sample in relation to the antenna significantly impacts the final result of microwave heating. In our particular situation, the choice of this angle was arbitrarily, but it is worth investigating the selection of an angle that allows for both the viewing of the sample while it is being heated and the establishment of a potentially even temperature distribution throughout its surface. The accuracy of the solution in the finite element method is significantly affected by the quality of the mesh. The Figure (b) presents the mesh configuration of the system under analysis. The Figure (b) illustrates an instance where the total number of finite tetrahedral elements was 776 410;. The meshing process was carried out in a customized manner, with the maximum size of each element set to 1/10 of the wavelength. This criterion was applied to the regions encompassing the antenna and the surrounding sample, where the heat transfer problem was not handled. In the sample domain regions, where the heat transfer was resolved, especially in the thin layers of epoxy and carbon composite, as well as in the surrounding of the small aging defects, the elements number was noticeably increased, as depicted in the Figure (c).

2.2. Specimens and the Experimental Setups

An experimental campaign was carried out with a CFRP sample with defects, which corresponded to the case of the numerical simulation with both ageing and central defect (CADS). The central defect was represented by the absence of an epoxy adhesive layer with a thickness of 1 mm, over an area of 100 mm × 100 mm, as already used in the previous work [51]. The sample was exposed to laboratory ambiance conditions for 10 years. The experimental setup of the test is shown in Figure 5.

Figure 6 presents a comparison of the condition of the tested sample before the aging process with the condition of the same sample after 10 years of the natural ageing process. Some differences are noticeable, but the defects themselves are not visible.

The same microwave system 2.45 GHz as in the numerical model was used to heat the sample. The pyramidal horn antenna was placed at 550 mm from the sample, with a 35° direction. A heating power of 1000 W was used for heating (continuous heating mode) the sample for 2 minutes. An infrared camera was used to record thermograms with a frame rate of 1 image per second. In order to observe the whole surface of the specimen, the following position of the infrared camera was chosen for our experimental test: normal direction with a distance of 1 m from the sample. It is worthwhile to notice that the test was carried out inside a protection enclosure (as used in the previous works [36,38]). The operator must stay outside the enclosure during the test.

3. Numerical Modelling Results

In this section the results of numerical modeling for all the previously introduced scenarios will be presented. First, we will present the raw numerical results to assess their quality, and subsequently compare them with the data from the experiments.

Figure 7 illustrates the numerical result of subjecting the sample to microwave heating, specifically focusing on the case with a central defect and defects that imitate the effects of ageing (CADS). The presented thermogram was acquired after a simulation period of 100 seconds, which is equal to the heating period used in the experiments as well. It is evident that when the sample is positioned at a certain angle to the antenna, a higher temperature is observed at the surface closer to the antenna. The overall temperature rise equals to around 4 °C, which enables the initial identification of all the defects.

Further results depicted in the Figure 8 and Figure 9 provide a more comprehensive examination of the numerical modeling outcomes. To enhance the presentation of the data, it was determined that only a part of the sample, specifically the composite reinforcement, would be displayed. First, Figure 8 illustrates a comparison of the results obtained for two samples without the defects simulating ageing process – (NDS) and (CDS). As it can be noticed, in this case the proposed method has a high efficacy as the location of the defect may be identified even after 10 seconds of heating – nevertheless, it is important to acknowledge that during short heating, the difference in temperature between the background and the defects is small (approximately 0.3 °C for the given 10-second heating duration). While it is theoretically feasible to detect this with a highly sensitive thermal imaging camera, it could be exceedingly challenging to identify the defect under real measurement conditions. As the duration of heating increases, the difference between the background and the defect obviously grows. After 50 seconds, it reaches roughly 1.5 °C, and after 100 seconds, increases to around 2.1 °C. Figures 8(a) clearly demonstrate the significant imbalance in sample heating that is associated with the aforementioned angle at which the sample is positioned relative to the antenna. Furthermore, a distinct symmetrical temperature pattern can be observed on the surface being tested, which corresponds directly to the radiating characteristics of the antenna. Nevertheless, it is evident that these aspects do not completely eliminate the defect visibilty. Nevertheless, they have an impact on the overall image and can potentially hinder the identification of a flaw.

The most interesting feature of this research is the analysis of several elements related to the occurrence of flaws in the sample associated with the process of aging. Our objective was to determine the degree to which the method of the active microwave thermography can be used to identify these flaws, as well as the assessment of the influence these flaws have on the ability to detect other defects in the sample.

The most interesting feature of this research is the analysis of several elements related to the occurrence of flaws in the sample associated with the process of aging. Our objective was to determine the degree to which the method of the active microwave thermography can be used to identify these flaws, as well as the assessment of the influence these flaws have on the ability to detect other defects in the sample.

Figure 9 presents a comparison of temperature distributions for ADS and CADS samples. Evidently, defects associated with the aging process are observable in the periphery of the composite layer. Defects that encompass both the epoxy layer and a portion of the composite layer heat up more compared to defects solely present in the epoxy layer. The provided information indicates that the proposed method enables the detection and quantification of the extent of aging changes. It is evident that the presence of aging changes does not significantly affect the detection of a defect located in the center. The defect remains clearly visible just like in the previous case shown in the Figure. The temperature differences between the background and the defect are approximately the same for subsequent time steps in both cases. It is important to mention that the temperature of aging flaws' thermal signatures is similar to the temperature of a centrally positioned defect's thermal signature (maximum 24.1 °C). If the sample is significantly damaged, it can pose challenges in identifying additional defects.

4. Experimental Results

In this section a selected set of experimental results will be presented. As it was shown before, the experiment consisted of heating a concrete sample that was reinforced with a carbon composite. The reinforcement had an artificially produced defect in its center, in the form of lack of adhesive. Furthermore, the sample underwent an aging process for a duration of 10 years, leading to the emergence of a number of defects along the edges of the composite layer. The Figure 10 depicts a representative outcome obtained by subjecting the sample to heating for a duration of 100 seconds. The presence of numerous defects caused by the aging process is evident on the periphery of the composite layer, resulting in higher temperature readings. Furthermore, within the central part of the reinforcement, one can notice the thermal signature of a defect that is positioned at the center.

Figure 11 illustrates a comparison of the temperature distribution at specific time intervals (10, 50, and 100 seconds) for the region including the composite reinforcement. The results are comparable to those obtained numerically which will be quantitatively evaluated in the next section, however it is evident that the experimental results show higher temperature values for defects associated to the sample aging process compared to the numerical results. This is influenced by various variables. In the numerical model, we encountered a random distribution of defects that may not necessarily correlate to this particular scenario. The depth of defects in the real object may be larger. The rise in temperature in the mentioned region could potentially be affected by wave reflections, which were not taken into account in the computer model.

5. Conclusions Results Analysis and Discussion

In this section, the examination of the acquired data and the attempt of the comparison the numerical and the experimental results will be presented. The quantitative influence of the presence of the defects caused by the object's ageing process on the capability to detect other sorts of defects will be evaluated. To conduct the analysis, the same region of interest (ROI) that covered the flaw located in the center for each numerical scenario was chosen. When analyzing the experimental results, not just the area that directly covered the defect, but also the surrounding area beneath the defect, which was presumed to be defect-free were considered. In the Figure 12, all the chosen regions are presented.

In order to assess the acquired data in a quantitative manner, an average value was computed for all designated regions of interest (Figure 12). Thus, the time - temperature characteristics were determined for the predicted mean temperature values in certain regions. First, the numerical results were compared for samples that did not have any aging defects (referred to as NDS and CDS samples), samples that had aging defects (referred to as ADS and CASD samples), and the sample that was physically tested in the experiment. In the comparison, both the regions containing the defect and the areas without the defect were taken into account. Figure 12 illustrates all the aforementioned characteristics. Figure 13 (a) presents a comparison of time - temperature characteristics for NDS and CDS samples. These characteristics were calculated by computing the average values from the areas marked in Figures 12 (a) and 12 (b). The largest disparity between the given features is roughly 1 °C. The comparison of characteristics between ADS and CADS samples is conducted by computing the averages for the regions shown in Figure 13 (b). Therefore, the existence of defects associated with the ageing of the material has little impact on the temperature distribution in the area of the defect situated at the center of the sample. Figure 13 (c) shows a comparison of analogous characteristics for a real sample that was evaluated via experimentation. In this case, the curves were derived by averaging the data from the specified areas shown in Figures 12 (e) and 12 (f). In this instance, the largest disparity between the characteristics is almost 2.3 °C which is a noteworthy difference in comparison to the numerical data. It is important to mention that the regions being compared in this context do not completely align with those being studied numerically.

The differences between the results shown in Figures 13(a) and 13(b) are minimal, therefore in order to quantitatively determine how the presence of aging defects affects the temperature distribution within the region encompassing centrally located defect, additionally the difference between the characteristics for the NDS and ADS samples has been plotted (which will show how the presence of aging defects affects the temperature distribution in the central area of the sample when there is no defect there) and for CDS and CADS samples (which will show how the presence of aging defects affects the temperature distribution in the central area of the sample when there is a defect there). Plots of the differences defined above are presented in the Figure 14. As it can be seen, the maximum values of $\mathsf{\Delta }\mathrm{T}$ are $2.5\cdot {10}^{-3}$ °C for the difference between NDS and ADS (Figure 13(a)) and $18\cdot {10}^{-3}$ °C for the difference in CDS and CADS (Figure 13(b)), respectively. This difference between samples without a central defect and with such a defect is significant, but not observable in real conditions.

Figure 15 illustrates a comparison between the characteristics acquired from the regions of interest for the matching numerical and experimental results. One graph displays the characteristics obtained for the CADS sample (Figure 12 (d)) and the real sample (Figure 11 (e)). The indicated outcome demonstrates a high level of consistency. The highest disparity between the characteristics is around 0.7 °C. The agreement observed confirms the accuracy of the assumptions made in the developed numerical model.

6. Conclusions

In this study, the feasibility of applying the microwave infrared thermography method (MIRT) for identifying flaws in a CFRP reinforcement in the form of lack of adhesive between the reinforcement and the concrete was examined. The most challenging element considered in the study involves using a sample that underwent an extended aging period of 10 years, which facilitated the development of natural deterioration in the reinforcement characterized by the fraying of the material at the edges and the chipping of the epoxy. Numerical models were used for comparison research and to ascertain the influence of the flaws, related to material deterioration associated with aging, on the efficacy of the MIRT, and were validated through experimental verification. For this purpose, four numerical models with four variants of the sample were prepared: a sample without any defects, a sample with a defect in the form of lack of adhesive (Central Defect Sample), a sample with defects simulating the aging of the material, located on the edges of the composite reinforcement layer (Ageing defects Sample), a sample with defects simulating aging and with a centrally located defect in the form of lack of adhesive (Central and Ageing Defects Sample). A stochastic structure consisting of many cylindrical air spaces that spanned both the epoxy adhesive layer and the adjacent carbon composite layer was used to represent the flaws associated with ageing process.

The results show that the proposed MIRT enables the detection and quantification of the extent of aging deterioration. The presence of aging changes does not significantly affect the detection of a defect located in the center of the sample. The comparison between the matching numerical and experimental results demonstrates a high level of consistency confirming the accuracy of the assumptions made in the developed numerical models.

It is worthwhile to notice that the experimental results show slightly higher temperature values for defects associated to the sample aging process compared to the numerical results, which may be caused by the random distribution of defects that may not necessarily correlate to the scenario, and the depth of defects in the real object may be larger. Moreover, the temperature increase in the mentioned region could potentially be affected by wave reflections, which were not considered in the numerical model. These remarks may be the keys for improving the accuracy of the proposed models in the future studies.

The present study offers potential perspectives in terms of applying the MIRT method to investigate the durability of different types of composites and other new construction materials under various environmental conditions.