Original Research ARTICLE
The Combined Thermal and Mechanical Effect of an Interstitial Gas on Thermal Rectification Between Periodically Grooved Surfaces
- Department of Mathematical Problems of Contact Mechanics, Pidstryhach Institute for Applied Problems of Mechanics and Mathematics, The National Academy of Sciences of Ukraine, Lviv, Ukraine
This paper presents a study on the combined thermal and mechanical effect of interstitial gas on thermal rectification between a periodically grooved surface and a flat one. To evaluate the interstitial medium influence, the analytico-numerical solution to the corresponding thermoelastic contact problem is constructed taking into account the effect of thermal strains on gap deformation. The results are provided for the Stainless Steel AISI 304—Aluminum Alloy A380 pair in the presence of air or krypton in the interface gaps. The effects of the gas thermal conductivity and pressure, the imposed pressure and heat flux, and the maximum groove height on the effective thermal contact resistance and the level of thermal rectification are analyzed. It is revealed that taking into account the mechanical effect of the gap filler leads to lower values of thermal rectification level. Also, the change of the gap filler thermal conductivity has more pronounced effect that the change of its pressure.
Today the tribological mechanical components of many modern devices are often regularly textured with grooves, dimples, protrusions, etc. of small height in order to improve their surface properties significantly (Etsion, 2004, 2005; Stepien, 2011). When textured surfaces are in contact, periodic arrays of gaps are present at the contact interface. The contacting textured surfaces are often subject to the action of not only normal and shear mechanical loads, but also thermal load. The thermal strains caused by the imposed heat flow can have a significant effect on the contact parameters of the bodies, especially dimensions of contact regions/interface gaps and thermal contact resistance.
The results of theoretical studies of thermoelastic contact between bodies of dissimilar materials in a vacuum environment (Barber, 1971; Dundurs and Panek, 1976; Stevenson et al., 1991; Kharitonov and Yakutin, 1997; Chumak and Martynyak, 2012) showed that the contact parameters and thermal contact resistance of the bodies depend not only on the heat flow density but also on the reversal of its direction. This phenomenon was repeatedly observed experimentally [a review of these experimental observations is available in Roberts and Walker (2011)] and is known as thermal rectification or the directional effect. The essence of thermal rectification is that the thermal contact resistance between two bodies is greater in one heat flow direction across the interface than in the reverse direction. Thermal rectification between textured surfaces is poorly investigated [the only study available was conducted by Chumak and Martynyak (2019)].
The interface gaps between textured surfaces are commonly filled with natural (gas, liquid) or operational (grease, coolant) substance. The interstitial medium exerts pressure on the contacting surfaces and offers additional thermal resistance, and thus can significantly affect the distribution of stress-strain state and temperature in the contacting bodies, and the level of thermal rectification.
The vast majority of investigations of the effect of the interface gap filler on the contact of bodies have been conducted taking into account either only filler pressure or only filler thermal conductivity.
The influence of the interstitial medium filling the interface gaps on the mechanical behavior of bodies has been studied by Machyshyn and Nagórko (2003), Slobodyan et al. (2014), Kozachok et al. (2017), Malanchuk et al. (2017), Serednytska et al. (2019) (gas-filled gaps), Evtushenko and Sulim (1981), Kuznetsov (1985), Kit et al. (2003), and Martynyak and Slobodyan (2009) (liquid-filled gaps). A number of works studied the influence of liquid bridges in interface gaps on contact interaction when the remainder of the gap is filled with a vacuum (Peng and Bhushan, 2002; Boer, 2007; Xue and Polycarpou, 2007; Persson, 2008; Wang et al., 2009) or with a gas (Martynyak and Slobodyan, 2008; Martynyak et al., 2009; Slobodyan, 2012; Kozachok et al., 2016; Kozachok and Martynyak, 2018). The normal contact approach and contact compliance of the bodies with periodically grooved or wavy surfaces and medium-filled interface gaps was investigated by Kozachok et al. (2018a,b).
The comprehensive review of the investigations of the thermal effect of the interface gap filler is provided in Chumak (2018) and Chumak and Martynyak (2019), where this effect and thermal rectification phenomena was studied for periodically grooved and wavy surfaces, respectively.
The works that take into account both gap filler pressure and thermal conductivity are limited to Yovanovich et al. (1982), Song et al. (1993), and Bahrami et al. (2004a,b). However, they did not take into account the influence of thermal strains that did not allow them to study thermal rectification.
To summarize, thermal rectification between textured surfaces has been poorly investigated. The previous studies of the contact between textured surfaces do not take into account the combined mechanical and thermal effect of the interstitial medium in the interface gaps and effect of thermal strains on the deformation of interface gaps at the same time. Considering only thermal or mechanical effect of the interface gap filler may greatly influence the outcome on the one hand, and neglecting the influence of thermal strains does not enable investigating thermal rectification phenomenon on the other hand.
The aims of the present study are:
(1) To construct the solution of the thermoelastic contact problem for a body with a periodically grooved surface and a body with a smooth surface in the presence of an interstitial gas in the interface gaps, which exerts pressure as well as offers thermal resistance.
(2) To investigate the combined pressure and thermal conductivity effect of the interstitial gas on thermal rectification between the contacting bodies.
Statement of Problem
Consider two bodies that are made from isotropic and dissimilar materials with Poisson's ratios ν1, ν2, shear moduli μ1, μ2, coefficients of linear thermal expansion α1, α2, and thermal conductivities λ1, λ2 (Figure 1). The surface of the lower body S1 is perfectly flat, and the surface of the upper body has regular surface texture in the form of periodically arranged grooves. The shape of the textured surface is described by the smooth periodic function r(x) of the coordinate x:
Here, b denotes the width of each groove, r0 is the maximum height of the grooves, d denotes the period of the surface pattern (d > b ).
The groove height is assumed to be small in comparison to the groove width, i.e., r0 ≪ b. This assumption will allow us to impose the boundary conditions not on the actual surfaces but write them on the plane y = 0 (Dundurs et al., 1973).
Since the problem is posed in the framework of linear thermoelasticity, assuming plane strain conditions, we consider the contact of two half-planes (Figure 2). They are compressed against each other by a remotely applied uniform pressure p, and a far-field uniform heat flux q is imposed in the direction normal to the interface.
Two opposite directions of heat flow are considered: heat flows from body S1 into body S2 when q > 0 (as in Figure 2), and heat flows from body S2 into body S1 when q < 0.
Due to the initial regular surface texture of the upper body, the interface between the solids consists of a periodic array of gaps and a periodic array of contacts.
Every gap is supposed to be filled with an equal amount of an ideal gas. The mechanical influence of the gap filler is simulated by pressure pg, which, according to the Pascal's law, is the same in each point of the gas. To simulate the thermal influence of the interstitial gas, we will use the following gap thermal resistance model (Chumak, 2018):
where h(x) denotes the gap height, and λg denotes the thermal conductivity of the gas. This model provides for gap deformation due to both mechanical and thermal loads. Also, it assumes that heat is transferred through the interstitial medium only by conduction in the y-direction.
The interface offers no thermal resistance outside the gaps.
The values of contact parameters lie within two limiting cases—frictionless contact (zero coefficient of friction) and complete stick (infinite coefficient of friction). As shown in Popov et al. (2019), the difference between the limiting values of the contact parameters is not large for normal contact of bodies. Based on this fact and since the allowance for friction forces will greatly increase complexity of the problem, we assume frictionless contact in the work.
If the magnitude of the imposed heat flow is to be the same for the contact with interface gaps and the contact without interface gaps, the additional far-field difference in temperatures ΔT must be supplied in the case of the periodic array of gaps. ΔT is clearly the cumulative effect of temperature disturbances arising out of each gap in the array, and equals to the average temperature jump across the interface (Das and Sadhal, 1997; Manners, 2000; Chumak, 2018):
where a denotes the interface gap width, γ(x) = T − (x, 0) − T+(x, 0) is the temperature jump across the interface, T(x, y) is the temperature, the superscripts + and − denote boundary values of a function on the x-axis in the upper and lower half-planes, respectively.
The effective thermal contact resistance, which is defined as Reff = ΔT/q (Madhusudana and Fletcher, 1986), is then
The boundary conditions at the interface are:
on the gap surfaces (|x − md| < 0.5a ):
in the contact regions (|x − md| ≥ 0.5a ):
The boundary conditions at the far-field are:
Here and further on, qx(x, y), qy(x, 0) are components of heat flow, σxx, σxy, σyy are stress components, uy is a displacement component, and the quantity is called the thermal distortivity of a material (Dundurs and Panek, 1976). The condition (11) describes application of linearly distributed stress to prevent the global warping of the contacting bodies and provide periodic contact along the whole interface (Panek and Dundurs, 1979).
Solution to Problem
Using the approach described in Chumak (2018), the considered thermoelastic contact problem has been reduced to the following system of singular integro-differential equations for the gap height h(ξ) and the temperature jump γ(ξ):
where , , , ,
The system (12, 13) is supplemented with the conditions (Chumak, 2018)
which ensure that the contact stresses and heat flows are bounded at the gap ends.
The system of singular integro-differential equations of a similar form was obtained by Chumak (2018) when solving the thermoelastic contact problem for wavy surfaces. To solve the system, the iterative procedure, developed by Martynyak and Chumak (2012), was adapted. For the sake of brevity, we omit the details, which can be found in Martynyak and Chumak (2012) and Chumak (2018), and provide only the final result:
Here, , α is supposed to be known value, Un(ξ) is the Chebyshev polynomial of the second kind, An, Bn, Cm are known constants. On each iteration step, these constants are determined from the corresponding systems of linear algebraic equations (for more details, see section 3.2 of Chumak, 2018).
The solution of the problem includes an unknown gas pressure pg. To calculate it, we can use the ideal gas law:
where mg denotes the mass of the gas, μg is the molar mass of the gas, Vg denotes the volume of the gas, Tg is the absolute temperature of the gas, and R denotes the universal gas constant (R = 8.3145 JK−1mol−1).
Expressing the gas volume in terms of the gap height h(ξ),
substituting (15) into (20) and performing the integration, we obtain
where l = 1 m.
We relate the absolute gas temperature to the temperatures of the gap surfaces as
where T0 = 273.15 K.
Since and (Chumak, 2018),
Substitution of (21) and (23) into (19) yields
Then, the applied pressure p can be calculated using (17).
Results and Discussion
Calculations are performed for the Stainless Steel AISI 304 (the material of the body S1)—Aluminum Alloy A380 (the material of the body S2) pair in the presence of air or krypton in the interface gaps. The properties of the materials and the interstitial gases are listed in Tables 1, 2, respectively. The other parameters are chosen to be: b = 1 mm, d = 4 mm. The notations that are used in the figures below are provided in Table 3.
Table 2. The properties of the interstitial gases (Dorf, 2003).
Since the thermal distortivity δ1 of Stainless Steel AISI 304 is greater than the thermal distortivity δ2 of Aluminum Alloy A380, q > 0 corresponds to the case of heat flowing into the material with the smaller thermal distortivity (blue lines in the figures) and q < 0 corresponds to the case of heat flowing into the material with the larger thermal distortivity (red lines in the figures).
Figure 3 shows the dependences of the effective thermal contact resistance Reff on the applied pressure p. The results are presented for the gas pressure pg = 5 MPa and two values of heat flux: q = 1 W/mm2 and q = −1 W/mm2. The effective thermal contact resistance changes upon reversal of the direction in which heat is flowing, i.e., thermal rectification between the contacting bodies takes place. At fixed values of input parameters, Reff is always greater when heat flows into the material with the larger thermal distortivity (red lines). For the two opposite directions of heat flow, Reff decreases as p increases and is always larger in the case of krypton. The greater is the maximum height r0 of the grooves, the greater is Reff.
It should be noted that red curves in Figure 3 do not intersect the axis p = 0. This is connected with the fact that the solution of the thermoelastic contact problem has been constructed under the assumption that the gap width cannot exceed the groove width (a ≤ b). However, the heat flow directed to the material with higher thermal distortivity has prevailing influence on the gap deformation at low applied pressures and causes spread of the gap outside the initial groove boundaries. In this case, the thermoelastic contact problem must be solved for a > b. This adds certain complexity to construction of the problem solution and goes beyond the scope of the investigations performed in the work.
The effective thermal contact resistance Reff vs. the gas pressure pg is shown in Figure 4 for p = 100 MPa and |q| = 1 W/mm2. An increase in pg leads to an increase in Reff.
Reff for pg = 10 MPa is greater Reff for pg = 0 MPa by
• 15.58% when r0 = 0.005 mm, the interstitial gas is air and q = 1 W/mm2;
• 14.24% when r0 = 0.005 mm, the interstitial gas is air and q = −1 W/mm2;
• 13.26% when r0 = 0.005 mm, the interstitial gas is krypton and q = 1 W/mm2;
• 10.96% when r0 = 0.005 mm, the interstitial gas is krypton and q = −1 W/mm2.
Reff for krypton is greater Reff for air by
• 35.62% when r0 = 0.005 mm, pg = 0 MPa, and q = 1 W/mm2;
• 37.24% when r0 = 0.005 mm, pg = 0 MPa, and q = −1 W/mm2;
• 33.85% when r0 = 0.005 mm, pg = 10 MPa, and q = 1 W/mm2;
• 34.85% when r0 = 0.005 mm, pg = 10 MPa, and q = −1 W/mm2.
As an absolute value of increases, Reff decreases when heat flows into the aluminum alloy and increases when heat flows into the stainless steel (Figure 5).
To assess a level of thermal rectification for a given |q|, we use the rectification index ε defined by the formula (Snaith et al., 1986; Madhusudana, 1993). Here, denotes the effective resistance in the case of heat flowing from the aluminum alloy to the stainless steel and denotes the effective resistance in the case of heat flowing from the stainless steel to the aluminum alloy.
ε increases with increasing p, besides more rapidly for lower pressures (Figure 6). When r0 = 0.005 mm, ε for air is greater than ε for krypton up to approx. p = 67 MPa, but ε for air is smaller than ε for krypton for larger values of p. When r0 = 0.01 mm, ε for air is greater than ε for krypton within the whole considered range of p.
ε for pg = 10 MPa is greater ε for pg = 0 MPa by
• 7.93% when r0 = 0.005 mm, the interstitial gas is air and |q| = 1 W/mm2;
• 11.44% when r0 = 0.005 mm, the interstitial gas is krypton and |q| = 1 W/mm2.
ε for krypton is greater ε for air by
• 11.16% when r0 = 0.005 mm, pg = 0 MPa and |q| = 1 W/mm2;
• 7.64% when r0 = 0.005 mm, pg = 10 MPa and |q| = 1 W/mm2.
The thermoelastic contact between a periodically grooved surface and a flat surface in the presence of an ideal gas in the interface gaps has been investigated. The combined thermal and mechanical effect of the interstitial medium was taken into account. The thermal influence of the gas was simulated by the gap thermal resistance that is directly proportional to the gap height and inversely proportional to the gas thermal conductivity. The mechanical influence of the gas was taken into account by its pressure.
The contact problem was reduced to a system of singular integro-differential equations for a temperature jump across the gap and a gap height. To solve the system, the iterative procedure, developed recently by the authors, was adapted. The relationship between the gas pressure and the applied pressure was obtained using the condition of boundedness of contact stresses at the gap ends. Expressing the gas volume in terms of the gap height and relating the absolute gas temperature to the temperatures of the gap surfaces, the ideal gas law was used to calculate the gas pressure.
To evaluate the effective thermal contact resistance and the level of thermal rectification, the Stainless Steel AISI 304—Aluminum Alloy A380 interface in the presence of air or krypton in the interface gaps was examined. It was revealed that the effective thermal contact resistance is greater when heat flows from the aluminum alloy to the stainless steel (into the material with the higher thermal distortivity) than for the reverse heat flow direction. The effective resistance for krypton-filled gaps is always larger than the effective resistance for air-filled gaps. The effective resistance increases as the applied pressure decreases or the gas pressure increases. An increase in the absolute value of heat flow leads to a decrease/an increase in the effective resistance when heat flows into the aluminum alloy/stainless steel. It is observed that the level of thermal rectification increases with increasing heat flow magnitude, decreasing gas pressure and increasing applied pressure. The rectification level changes variously with the thermal conductivity of the interstitial gas depending on the applied pressure.
Taking into account the mechanical effect of the gap filler leads to lower values of thermal rectification level. Also, the change of the gap filler thermal conductivity has more pronounced effect that the change of its pressure.
RM designed the model and the study, and supervised the findings of this work. KC wrote the manuscript with support from RM, performed the analytic calculations and the numerical simulations, and designed the figures. All authors discussed the results and contributed to the final manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Chumak, K., and Martynyak, R. (2012). Thermal rectification between two thermoelastic solids with a periodic array of rough zones at the interface. Int. J. Heat Mass. Transf. 55, 5603–5608. doi: 10.1016/j.ijheatmasstransfer.2012.05.049
Chumak, K. A., and Martynyak, R. M. (2019). Effective thermal contact resistance of regularly textured bodies in the presence of intercontact heat-conducting media and the phenomenon of thermal rectification. J. Math. Sci. 236, 160–171. doi: 10.1007/s10958-018-4103-7
Kozachok, O., and Martynyak, R. (2018). Contact problem for wavy surfaces in the presence of an incompressible liquid and a gas in interface gaps. Math. Mech. Solids 1–13. doi: 10.1177/1081286518781679
Kozachok, O. P., Slobodian, B. S., and Martynyak, R. M. (2017). Interaction of two elastic bodies in the presence of periodically located gaps filled with a real gas. J. Math. Sci. 222, 131–142. doi: 10.1007/s10958-017-3287-6
Kozachok, O. P., Slobodian, B. S., and Martynyak, R. M. (2018b). Contact of two elastic bodies with wavy topography in the presence of gases in gaps. Mater. Sci. 54, 194–201. doi: 10.1007/s11003-018-0173-4
Kozachok, O. P., Slobodyan, B. S., and Martynyak, R. M. (2016). Contact of elastic bodies in the presence of gas and incompressible liquid in periodic interface gaps. Mater. Sci. 51, 804–813. doi: 10.1007/s11003-016-9905-5
Machyshyn, I., and Nagórko, W. (2003). Interaction between a stratified elastic half-space and an irregular base allowing for the intercontact gas. J. Theor. Appl. Mech. 41, 271–288. Available online at: http://www.ptmts.org.pl/jtam/index.php/jtam/article/view/v41n2p271
Malanchuk, N., Slobodyan, B., and Martynyak, R. (2017). Frictional contact of two solids with a periodically grooved surface in the presence of an ideal gas in interface gaps. J. Theor. Appl. Mech. 55, 1181–1192. doi: 10.15632/jtam-pl.55.4.1181
Martynyak, R., and Chumak, K. (2012). Effect of heat-conductive filler on interface gap on thermoelastic contact of solids. Int. J. Heat Mass Transf. 55, 1170–1178. doi: 10.1016/j.ijheatmasstransfer.2011.09.053
Martynyak, R. M., and Slobodyan, B. S. (2008). Influence of liquid bridges in the interface gap on the contact of bodies made of compliant materials. Mater. Sci. 44, 147–155. doi: 10.1007/s11003-008-9068-0
Martynyak, R. M., Slobodyan, B. S., and Zelenyak, V. M. (2009). Pressure of an elastic half space on a rigid base with rectangular hole in the case of a liquid bridge between them. J. Math. Sci. 160, 470–477. doi: 10.1007/s10958-009-9511-2
Serednytska, K. R., Martynyak, R., and Chumak, K. (2019). The thermoelastic state of a bi-material with an open gas-filled interface crack. J. Theor. Appl. Mech. 57, 331–341. doi: 10.15632/jtam-pl/104587
Slobodyan, B. S. (2012). Pressure of an elastic body on a rigid base with a recess partially filled with a liquid that does not wet their surfaces. Mater. Sci. 47, 561–568. doi: 10.1007/s11003-012-9428-7
Slobodyan, B. S., Malanchuk, N. I., Martynyak, R. M., Lyashenko, B. A., and Marchuk, V. E. (2014). Local sliding of elastic bodies in the presence of gas in the intercontact gap. Mater. Sci. 50, 261–268. doi: 10.1007/s11003-014-9716-5
Keywords: thermal rectification, interstitial gas, periodically grooved surface, thermoelastic contact, mechanical effect, thermal effect
Citation: Chumak K and Martynyak R (2019) The Combined Thermal and Mechanical Effect of an Interstitial Gas on Thermal Rectification Between Periodically Grooved Surfaces. Front. Mech. Eng. 5:42. doi: 10.3389/fmech.2019.00042
Received: 06 February 2019; Accepted: 26 June 2019;
Published: 16 July 2019.
Edited by:Alessandro Ruggiero, University of Salerno, Italy
Reviewed by:Jiliang Mo, Southwest Jiaotong University, China
Yoshitaka Nakanishi, Kumamoto University, Japan
Valentin L. Popov, Technische Universität Berlin, Germany
Copyright © 2019 Chumak and Martynyak. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rostyslav Martynyak, email@example.com