Estimation of thermal transmittance based on temperature measurements with the application of perturbation numbers

Fast estimation of thermal transmittance based on temperature measurements is uncertain, and the obtained results can be burdened with a large error. Nevertheless, such attempts should be undertaken merely due to the fact that a precise measurement by means of heat flux measurements is not always possible in field conditions (resentment of the residents during the measurements carried out inside their living quarters), and the calculation methods do not allow for the nonlinearity of thermal insulation, heat bridges or other fragments of building envelope of diversified thermal conductivity. The present paper offers the estimation of thermal transmittance and internal surface resistance with the use of temperature measurements (in particular with the use of thermovision). The proposed method has been verified through tests carried out on a laboratory test stand built in the open space, subjected to the influence of real meteorological conditions. The present elaboration involves the estimation of thermal transmittance by means of temperature measurements. Basing on the mentioned estimation, the authors present correction coefficients which have impact on the estimation accuracy. Furthermore, in the final part of the paper, various types of disturbance were allowed for using perturbation numbers, and the introduced by the authors “credibility area of thermal transmittance estimation” was determined.


Introduction
The estimated determination of thermal transmittance based on the measurements is burdened with substantial uncertainty when defining heat loss from a building envelope and consequently the energy efficiency in buildings. On the other hand, the calculation methods which consist in mathematical modeling resort to certain assumptions or simplifications, which also leads to some sort of uncertainty when estimating heat loss. The calculations do not allow for thermal diversity of walls, or, in other words, they do not take into consideration the fact that there are places of different thermal conductivity. When the thermal transmittance U is in reality higher than estimated, the conservation of energy is in fact lower than estimated, which in the macro scale (e.g. the whole country) can be significant. Many energy performance certificates can be incorrectly compiled due to the incorrect estimation of thermal transmittance U. Such a situation was observed by Lowe who described it in his paper [1].
Over the recent years, many research works have been published in which in-situ measurements of thermal transmittance have been attempted, as for example in the work [2]. Another work covering the said issue was presented by Stevens [3], in which the consumption of energy measured in 93 apartments with respect to the measured value of thermal transmittance was presented. In the paper presented by Rye [4], the author declares that traditional buildings are energyinefficient, in particular with respect to the heat loss through walls. Also Evangelsiti [5] declares that an accurate assessment of the thermal transmittance U of walls is indispensable to calculate the annual consumption of energy. In his work Evangelisti analyzed three walls of the building having different thermal properties and he compared the obtained results with the calculation methods.
Basing on the quoted works as well as on other works [6][7][8][9], we know that the thermal properties of a building have direct impact on the annual consumption of energy. We also know that the thermal conductivity of a wall depends on the thermal conductivity of each layer of that wall, and therefore, the thermal resistance of the whole wall can be modeled in the way analogous to the resistance of resistors in electrical Jerzy Skrzypczyk died before publication of this work was completed. current flow models. Such a thermal-electrical analogy has been extensively described in literature and applied for the calculation of heat exchange by the envelopes whereof thermal properties are known. Nevertheless, the real (measured) value of thermal transmittance U is frequently different from the calculated value, and it is definitely different when we take into consideration e.g. the nonlinearity of insulation layers or heat bridges. Therefore, in order to estimate the real heat loss through the outer envelope of the building, the measurements are essential [10,11]. The generally acknowledged nondestructive measurement method of thermal resistance is the direct measurement of heat flux [12,13]. The measurements of thermal resistance can fail to reflect the resistance of the whole envelope if the measurement points are selected at places having thermal conductivity different from the conductivity of the whole envelope. Therefore, it is worth while complementing this measurement with the thermographic analysis [14]. Measurement methods can be also applied for the analysis of historical buildings. Walker and Pavia [15] have been carrying out in-situ research investigation studies involving thermal insulation of seven envelopes with internal insulation on a historical wall. We must emphasize here that the analyses involving the measured values of thermal transmittance have not been carried out for the last two years. In 2009 Jimenez [16] was investigating the said issue, and Dudek [17] demonstrated in a series of in-situ measurements of masonry walls that the discrepancy between the measurements and the calculation methods reached 25%. Naveros and Jiménez [18] investigated also the potentials and limitations of the regression method based on the average values applied for the thermal analysis of the performance of real-size building elements tested in the external dynamic atmospheric conditions.
Other researchers, e.g. Cesaratto [19], using FEM, have been simulating the impact of various external interference, or variable input data of temperature on the in-situ measurement results of thermal conductivity. They presented a general assessment of deviations involving the obtained results for various cases.
According to the authors [20], energy audits and the control of energy performance are very important in terms of energy conservation in buildings. Therefore, it is essential to be able to estimate thermal transmittance of building envelopes. Ficco as well as other researches mentioned in this Introduction base their measurements principally on the measurement of heat flux through the wall. Although the latter is investigating the impact of the meteorological conditions of the environment on the measurement results involving the estimation of thermal transmittance U and on the uncertainty of that estimation.
The latest research studies [21] present the design of a more straightforward and compact version of the traditional Hot-Box apparatus (it measures U-value) which instead determines the thermal resistance of the samples of building envelope.
The work [22] et al. presents the description of in-situ monitoring measurements, using the measurements of temperature and heat flux flows as a more precise measurement of the thermal performance of walls. However, the authors indicated some limitations, as for example the fact that with the use of standard methods, the measurements in winter must be carried out over the time period of two weeks. The work Krause [23] also presents thermovisual and numerical studies, but they are focused on the identification of thermal bridges. There are empirical models for the prediction of thermal conductivity estimation. To provide an example, we can mention the model for cement-based porous sensible heat storage materials and naturally occurring crystalline rock formations as a function of temperature [24].
And in the work [25], they applied the measurement data and Autoregressive models with exogenous (ARX) for the identification of physical parameters of the tested walls in real meteorological conditions.
One of the measurement methods of temperature is thermovision. In fact, we can state that the thermovisual method is burdened with the risk of incorrect interpretation of results. As we can read in the work of Fox et al. [26], the thermographic image of the building is reflecting the thermal condition at a given moment of time, and it does not reflect changes in time. They also noticed that the building was hardly ever in the state of thermal equilibrium, which might result in the incorrect interpretation of building defects in effect of the application of the standard thermographic method. For the same reason Taylor [27] is supporting the thermographic research studies with computer simulations, trying to find out how these two techniques can be used together. Numerical modeling and impact analysis of various factors on the final results necessitates that a lot of data is Bmemorized^. Is it particularly noticeable when creating the models of bigger, more complex layered systems. In such cases a considerable part of the operating memory must be engaged to ensure proper operation of respective programs. When we have to investigate a larger number of uncertain parameters, the task is becoming even more complicated. Hence, new algorithms are being created all the time to ensure better remembering capacity by the program of possibly the greatest amount of data with the minimum employment of memory. The analyses involving the calculation of temperature fields are most commonly carried out with the use of the finite elements method and boundary element method, as well as with the application of finite difference method and elementary balance method [28]. We can more and more frequently witness the application of neural networks, fractals theory or perturbation methods in different branches of engineering. Artificial neural networks have been applied among others to solve the issue of inverse heat conduction, with the assumption of functional dependence of heat conduction on temperature, and to predict the effective heat conduction in porous materials [29,30]. Perturbation methods are applied almost in all branches of science, including technology. For example, in the work [31] the material properties of thermoplastic materials were investigated. Perturbation methods are also used to analyze the stability of dynamic systems. Such investigation studies have been carried out for example by Nowoświat [32], and they involved the control and stability of Duffing oscillator. Perturbation methods have been applied to investigate the buckling of springy beams [33]. Very interesting modifications of this method and the explanation of perturbation numbers have been proposed by Skrzypczyk [34]. Perturbation methods are applied among others to obtain an approximate analytical solution involving the problems of one-dimensional thermal transmittance through porous materials. Such materials are characterized by a very low design thermal conductivity, which can play the role of a small parameter [28,35].
All the methods offered by the researchers presented in this Introduction necessitate long-term measurements and very often interfere with the privacy of people of the involved apartments. Consequently, such measurements in real conditions in many cases might be impossible to implement. Therefore, the authors of the present work are putting forward a proposition enabling fast estimation of thermal transmittance. Such a fast diagnosis can be applied to estimate heat loss through building envelopes in real conditions when the measurement with a heat meter is not feasible. Furthermore, with the use of the proposed fast estimation, we can determine an equivalent thermal transmittance for envelopes having diversified thermal activity.
The application of approximate analytical methods, which encompass perturbation techniques, can be found useful to solve the problems of heat transport. The basic assumptions of the traditional perturbation theory necessitate that the model should be transformed into the dimensionless form in such a way so that the parameters and variables involving the behavior of the system could be determined. As the next step, we select a parameter (referred to as small parameter or perturbation value) which is small as compared to the other ones and labeled as ε. The present work is attempting to apply an alternative perturbation method which consists in the application of n-perturbation numbers defined by Skrzypczyk [34] and used by Winkler-Skalna [36]. The overview of methods instigated the authors to undertake works involving the estimation capability of thermal transmittance by means of temperature measurement. Furthermore, using the algebra of perturbation numbers, we determined in the work a so called Bcredibility area^which almost certainly contains the measurement result. The application of the algebra of perturbation numbers for the estimation of thermal transmittance U is an innovative approach in such applications.

Estimation of thermal transmittance
Solving the Fourier's heat conduction equation and using the definition of thermal resistance as the quotient of heat flow rate to the difference of the internal and external environment temperatures we obtain: where q heat flow rate, W/m 2 λ design thermal conductivity, W/mK θ i internal surrounding temperature, K θ e external surrounding temperature, K α i internal surface heat transfer coefficient, W/ m 2 K R si internal surface resistance, with R si ¼ 1 α i , m 2 K/W θ si internal surface temperature, K When using the Eq. (3) for the estimation of thermal transmittance on the basis of temperature, we might find it difficult to determine, or estimate the internal surface resistance R Si . Thermal transmittance can be also determined with the use of external surface resistance (from the external side of building envelope) and external surface temperature, but it is burdened with a much higher error. The error in question is effected by high changeability of the surrounding temperature on the external side of the envelope, by the changing speed of wind influencing the external heat transfer coefficient, by surface temperature interference from the outside (sunlight, shading, wind, etc.).
The Eqs. (2) and (3) are used with the assumption of stationary heat flow. However, due to the inability to stabilize the external atmospheric conditions, it is not possible to stabilize the temperature of envelope surface. Therefore, the changes were defined as dynamic ones.
Taking into account the above factors, the estimation of thermal transmittance was carried out with the use of Eq. (3).
The determination of thermal transmittance is carried out on the basis of the criteria equation where: Nu is Nusselt number, Re -Reynolds number, Pr -Prandtl number, Gr -Grashof number, K ggeometric similarity.
In the conditions of free convection, the speed of air circulation in rooms is from several to a dozen or so centimeters per second, and the difference between air temperature and envelope surface is most commonly within the range of 2 ÷ 8 K. In such conditions, the surface heat transfer coefficient through/ by convection on the surface of walls can be determined from the approximate Equation [37]: The surface heat transfer coefficient can be also approximated with other relations. For example, with the temperature difference between the air and the surface not higher than 5 K the above equation can have the following form: For the differences higher than 5 K, but not only, instead of the Eq. (6) we can use When we want to apply the calculation method of thermal transmittance by means of the Eq. (3), we take into account places in the envelope which have considerably different surface temperatures, as in Fig. 1.
then the equivalent thermal transmittance is:

Perturbation numbers
Basing on the works of Skrzypczyk [34], we can write that the n-perturbation number is defined as an ordered set of real numbers (x 0 , x 1 , x 2 , …, x n ) ∈ ℝ n . The set of n-perturbation numbers is labeled as R nε . If two numbers x, y ∈ R nε are nperturbation numbers written as x = (x 0 , x 1 , x 2 , …, x n ), y = (y 0 , y 1 , y 2 , …, y n ), x i , y i ∈ ℝ, i = 0, …, n, we can say that x ≡ y if and only if x i = y i for each i = 0, …, n. In the set of n-perturbation numbers R nε we introduced the operations of adding (+ ε ) and multiplication (• ε ), as well as neutral elements of adding and multiplication were defined. In the set R nε defined in such a way, which was labeled as the body of n-perturbation numbers, we retain the commutativity and connectivity of adding, and also the commutativity and connectivity of multiplication, as well as the separation principle. The presented above body R nε does not contain the body of real numbers R. However, we can demonstrate that real numbers can be treated as certain elements of the body R nε , with the algebraic operations and neutral elements of adding and multiplication being maintained [34].
Since the method has not been described in the world literature for the general audience, we will present here shortly the algebra of perturbation numbers defined by professor Skrzypczyk.
For three n-perturbation numbers x = (x 0 , x 1 , x 2 , …, x n ), y = (y 0 , y 1 , y 2 , …, y n ), z = (z 0 , z 1 , z 2 , …, z n ) and for the real number A we can write: the sum of two n-perturbation numbers: the sum of the real number and the n-perturbation number: the product of two n-perturbation numbers: the inverse of n-perturbation number: the division of two n-perturbation numbers: Fig. 1 Temperature distribution on the surface of envelope with the marked places of different surface temperatures Let D ⊂ R nε be an arbitrary subset. We can say that we have a definite function f nε , if to each number z ∈ D is assigned precisely one element from the set R nε . Then, we can say that f nε is the extension of the function defined on the subset D with the values contained in R nε . The n-perturbation function can be written as follows: f nε : D ⊂ R nε or, let f nε (z) = u(z) + ε 1 v(z) + ε 2 g(z) + … + ε n h(z), with u(.), v(.), g(.), …, h(.) standing for the real functions of the nperturbation variable z = z 0 + z 1 ε 1 + z 2 ε 2 + … + z n ε n or, in other words, n + 1 real variables z 0 , z 1 , z 2 , …, z n , which can be written as Hence u(.), v(.), g(.), …, h(.) will be labeled as ordinary real functions without the index ε. The functions u(.), v(.), g(.), …, h(.) will be labeled respectively as main part, first perturbation, second perturbation and n-th perturbation of the function f nε (.).

Experiment
The test stand was made as a small test building, located at the place having the following geographical coordinates 50°17′ 26.3″N + 18°40′34.7″E. The external walls were made from light curtain walls filled with mineral wool. The wall on the southern side was divided into three parts. Each part separated by wooden elements was filled with a different type of material making up the envelopes. The discussed test stand is presented in Fig. 2.
The stand was furnished with three test envelopes marked in Fig. 3 as 1,2,3. All the test envelopes had the same width of 25 cm. The envelopes were selected, taking into account the materials most commonly applied in the general construction of multifamily residential buildings.
The test envelopes (described in Table 1) were selected in the way ensuring that they have different thermal conductivity and that they can reflect their common application in building practice.
Thermal performance of the envelopes was tested by means of the equipment presented below: Main equipmentstationary: -Multi-channel recorder type MA56902M09TG3, -Heat flux transducer 250x250mm type FQA019C -3 items (one item for each of the tested envelopes), -Temperature sensor type Pt-100 -20 items (6 items for each of the tested envelope) + spare ones. -Stationary meteorological stationa meteorological multisensor with the heating type FMA510H.
The measurement of heat flux (Fig. 4) was carried out with the use of a heat meter FQ90119 (Ahlborn) of the surface of 0.25 × 0.25 m. The sensor was placed on the internal surface of the wall in its geometrical centre in order to minimize the impact of possible thermal interference (e.g. heat bridges, air infiltration, etc.). The surface of the wall was appropriately prepared to ensure optimal adhesion between the tested wall and the heat meter. The temperature was measured in two ways: internal temperature in the tested room as well as on the internal and external surfaces of the envelope, using the sensors PT100. The temperature of the outside air was measured by means of the meteorological station described above.
The applied testing (Fig. 5) methodology consisted in observing the measurement procedure for stationary tests comprised in the Standard [38].

Results of laboratory tests
Assuming the internal surface resistance R si ¼ 1 α i and using the Eqs. (5), (6) and (7), we will define the best, out of the proposed ones, estimation method of the said resistance, verified by the measurement (Fig. 6). To define the measured value of surface resistance we will use the transformation of the Eq. (3), which yields Knowing the value of surface resistance, we can estimate thermal transmittance from the Eq. (3). For that purpose, we apply the Eqs. (5)- (7). The thermal resistance estimated in that way was verified by the applied measurement method.
In order to define the measured value of thermal transmittance U, two approaches were applied. The first one (referred to as Measurement 1) consisted in the measurement of heat flow rate q and the temperature difference of the inside and outside air. The advantage of this method is that we can directly determine the U-value (with the internal and external heat transfer coefficients taken into consideration) from the obtained measurement results. The disadvantage is that the results are sensitive to fast changes of environmental conditions.
Depending on temperature fluctuations, the value of thermal transmittance U will be different, which leads to longlasting measurements and it hampers the possibility to compare the results between the stands located in different places.
The second approach (referred to as Measurement 2), which employs the measurement procedure described here, requires the measurement of heat flow rate q and the difference in temperatures on the internal and external surfaces of the envelope. Additionally, when determining thermal resistance of the envelope, the values of surface resistance should The studies were carried out at the stand as presented in The uncertainty of the measurements 1 and 2, as the measurements continuous in time, was calculated on the basis of the known error propagation methods obtained for the temperature measurements of: external environment, internal environment, external and internal surfaces of the wall as well as for heat flow rate.
Error estimation for the Eqs. 5, 6, 7 was based on the standard error determined for the particular temperatures and on the calculation of the standard uncertainty of composite number. The calculations were carried out in accordance with  Tables 2, 3 and 4. The standard deviation of the averaged internal temperature of the surrounding is σ θ i ¼ 1:93°C and the surface temperature σ θ si ¼ 2:12°C. The average temperature of the surrounding was θ i ¼ 19:24°C and the average temperature of the envelope surface was θ si ¼ 17:06°C.
Using the error propagation rules, we can estimate the measurement error for the relation (3). The results of such an estimation are presented in Table 2.
Similar results are presented for the wall 2.
The standard deviation of the averaged internal temperature of the surrounding is σ θ i ¼ 1:96°C and the surface temperatures σ θ si ¼ 2:3°C; σ θ si ¼ 3:1°C. The average temperature of the surrounding was θ i ¼ 19:25°C and the average temperature of the envelope surface was θ si ¼ 12:10°C.
Using the error propagation rules we can estimate the measurement error for the relation (3). The results of such an estimation are presented in Table 3.
The last measurement and estimation involved the wall 3. The standard deviation of the averaged internal temperature of the surrounding is σ θ i ¼ 2:6°C and the surface temperature σ θ si ¼ 3:1°C. The average temperature of the surrounding was θ i ¼ 22:15°C and the average temperature of the envelope surface was θ si ¼ 15:73°C .
Using the error propagation rules we can estimate the measurement error for the relation (3). The results of such an estimation are presented in Table 4.
We can see from the analysis of errors (Tables 2, 3 and 4) and from the presented graphs (Figs. 5, 6 and 7) that with the rise of thermal resistance of the investigated envelope, the Eqs. (5)-(7) yield the results closer to the measurements 1, 2. And together with the deterioration of thermal performance of the envelope, the results show lower agreement.

Analysis of results and discussion
As we already know, and which has also been confirmed in the present work, the thermal transmittance U determined with the use of Eq. (3) depends on the surface resistance. Furthermore, the estimation inaccuracy of thermal transmittance by means of the temperature method results from the sensitivity of the relation (3) to the slightest changes of this resistance. Therefore, the estimation of U on the basis of the , this estimation is satisfying. As we can see in Fig. 6 for the wall 1, the estimation of surface resistance with the Eqs. (5) or also (6) is closer to the tabular surface resistance [39] than the measurement. Admittedly, it does not bespeak of a wrong measurement method but rather of an inaccurate tabular data [39]. And looking at it the other way round, we can state that the estimation by means of Eqs. (5) or also (6) is closer to the measured value than the tabular value of resistance [39]. Therefore, although the estimation based on the Eqs. (5), (6) is worse than the measurement of heat flow rate, yet taking into account the measurement speed, the mentioned equations may be sufficient. It looks slightly worse on the wall 3 and the worst on the wall 2. We can observe that with the rise of temperature difference of the surrounding and the envelope, the estimation is becoming less precise. Let us survey now the estimation results of the thermal transmittance U.
First of all, we can observe that for each of the investigated walls we have a different temperature difference of the surrounding and surface, which is demonstrated in the figure below (Fig. 10).
We can observe that the thermal transmittance determined with the Eq. (3) with the use of Eqs. (5)- (7) in all the investigated cases is different. And independently of the thermal performance and also independently of the temperature difference of the surrounding and wall surface, the thermal transmittance U (3) with the application of Eq. (5) is closest to the measurement, and with the use of the Eq. (7) it is the furthest. And the difference between the estimation (3) and the measurement is progressively higher with the rise of temperature difference of the surrounding and surface.
Another interesting observation involves the fact that the graphs (3) in the function of time with the application of Eqs. (5)-(7) respectively are almost parallel. It is particularly visible on graphs 7b, 8b and 9b at the average temperatures of the surrounding and surface. In the same way the measurement graph in the function of time is almost parallel to the remaining ones. And for the rooms having a small temperature difference of the surrounding and wall surface (walls with thermal insulation), the graphs are almost overlapping.
We can therefore state that for well insulated walls, the measurement and estimation based on (3), and in particular with the application of the Eq. (5), do not differ significantly. Whereas for the case of poorly insulated walls, where the differences of surrounding and wall surface are large, we can find a coefficient correcting the estimation with the Eq. (3).
Such an assumption is significant particularly when it involves the application of Eq. (8) which can assume now the following form: where k 1 , k 2 , …, k n are correcting coefficients. And especially for places of the wall where the temperature difference of the surface and surrounding is small (good insulation, no discontinuities of insulation or heat bridge), we can assume k i = 1.
For example, for the walls investigated in the present work. The results of such a correction are presented in Fig. 11. The determined correction coefficients quite well reflect the following cases:  Error δU 0.03 0.03 0.35 0.29 0.22 k 1 for non-insulated envelopes U ≥ 2 W m 2 K , and k 2 for poorly insulated envelopes 1:0 W m 2 K ≤ U ≤ 2:0 W m 2 K ; We are tempted to claim here that when carrying out further research, it is possible to determine correction coefficients k for other cases, or to narrow them down and make them more specific for the cases that we suggested.
The results presented on the graphs 11 can be placed in the so called Bcredibility area^of results defined by the authors. For that purpose, perturbation methods will be applied.
Let us assume that in the Eq. (3) all parameters can be disrupted by various factors, with the assumption that each of the parameters will be dependent only on a part of the factors which may have impact on the final result. We assumed that the number of all factors may be n where of 11 were identified as follows: 1measurement errors, 2humidity of internal air, 3humidity of external air, 4temperature of internal air, 5temperature of external air, 6humidity content of the envelope, 7thermal conductivity index of the envelope, 8thickness of the envelope, 9surface of the envelope, 10impact of wind, 11insolation.
The measurements described in the present paper were devoid of the described by Li error sources, and yet the estimation results could be accepted as satisfying. We can accept as true that with respect to an inhabited dwelling place the estimation uncertainty with (3) using the Eqs. (5)-(7) and in particular (5) will not be higher than the measurement of heat flux, or at least not higher than the uncertainty described by Li. The latest investigation studies [41] also show a considerable impact of the thermal conductivity envelope and temperature amplitudes on thermal resistance U. The research studies described in the present paper also demonstrate the influence of temperature changes on the estimation of thermal transmittance. However, the proposed estimation method is different.

Conclusions
The present paper is attempting to estimate thermal transmittance U by means of temperature measurements. The research studies were carried out for three different wall constructions. The walls were selected in the way ensuring that they have a considerably different thermal resistance. Furthermore, the obtained estimation results were compared with the results based on the measurement of heat flow rate.
The results demonstrate that: & the thermal transmittance U can be estimated by means of temperature measurement. Furthermore, with respect to the envelopes of low thermal resistance, we can find a correction coefficient k which makes it possible to estimate the result in a quite precise way, & making use of the measurement of temperature and applying the Eq. (11), we can determine the equivalent thermal transmittance, & using the algebra of perturbation numbers we can determine a so called Bcredibility area^which almost certainly comprises the measurement result.  Fig. 11 Graphs of thermal transmittance in the function of time; U_measurement stands for the measured value, U_5 stands for the estimation with (3) using (5), U_k*(5) stands for the estimation with (3) using (5), taking into account the correction coefficient k 1 Fig. 12 Graphs of thermal transmittance as a function of time. Measurement stands for the measured value, Eq. 5 means the estimation by means of (3) using (5), credibility area was determined using the relation (13) We can also observe that the presented method can serve as an inspiration to use the thermovisual measurement for the estimation of thermal transmittance. Such an approach can replace the measurement by means of heat flow rate, which considerably facilitates the measurement without any substantial loss involving the results.
The elaboration of assumptions to apply thermovisual measurement as temperature measurement in the estimation of thermal transmittance U might be the subject of further research studies.