Impact of tent shade on heat exposures and simulated heat strain for people experiencing homelessness.

Concurrent increases in homelessness and heat intensity, duration, and frequency translate to an urban heat risk trap for the unsheltered population. Homelessness is both a driver and consequence of poor health, co-creating distinct geographies with various risk factors that exacerbate heat vulnerability. We tested the efficacy of different tent shadings over identical tents often observed in the Phoenix area (white bedsheet, mylar, tarp, and aluminum foil) and compared them to a control tent (uncovered) and ambient conditions. We monitored all meteorological variables at all six locations, notably Mean Radiant Temperature (MRT). The in-tent microclimate variability was applied to complete statistical and physiological modeling including substance use on heat strain. Findings indicate that tent shadings resulted in significantly lower in-tent MRT during the day (p < 0.05), but exacerbated in-tent thermal risk during the night compared to the control tent and ambient conditions. Furthermore, we found evidence that the temperature metric matters, and using only either MRT or air temperature (Tair) to assess "heat" could lead to inconsistent conclusions about in-tent microclimate. Interactions between shade types and time significantly amplified in-tent thermal risk. Physiological modeling indicates a higher risk of heat strain (core temperature beyond 40˚C) for people using substances. Decision makers should promote testing different heat intervening strategies toward realizing effective means of protecting human life and preventing heat illnesses. This study illuminates the need for an interdisciplinary approach to studying tents as shelters that considers the total heat load with heat strain modeling.

A simple three-cylinder radiometer and low-speed anemometer to characterize human extreme heat exposure.

As populations and temperatures of urban areas swell, more people face extreme heat and are at increasing risk of adverse health outcomes. Radiation accounts for much of human heat exposure but is rarely used as heat metric due to a lack of cost-effective and accurate sensors. To this end, we fuse the concepts of a three-globe radiometer-anemometer with a cylindrical human body shape representation, which is more realistic than a spherical representation. Using cost-effective and readily available materials, we fabricated two combinations of three cylinders with varying surface properties. These simple devices measure the convection coefficient and the shortwave and longwave radiative fluxes. We tested the devices in a wind tunnel and at fourteen outdoor sites during July 2023's record-setting heat wave in Tempe, Arizona. The average difference between pedestrian-level mean radiant temperature (MRT) measured using research-grade 3-way net radiometers and the three-cylinder setup was 0.4 ± 3.0 °C ( ±  1 SD). At most, we observed a 10 °C MRT difference on a white roof site with extreme MRT values (70 °C to 80 °C), which will be addressed through discussed design changes to the system. The measured heat transfer coefficient can be used to calculate wind speed below 2 m·s-1; thus, the three cylinders combined also serve as a low-speed anemometer. The novel setup could be used in affordable biometeorological stations and deployed across urban landscapes to build human-relevant heat sensing networks.

Impact of human body shape on forced convection heat transfer.

Predicting human thermal comfort and safety requires quantitative knowledge of the convective heat transfer between the body and its surrounding. So far, convective heat transfer coefficient correlations have been based only upon measurements or simulations of the average body shape of an adult. To address this knowledge gap, here we quantify the impact of adult human body shape on forced convection. To do this, we generated fifty three-dimensional human body meshes covering 1st to 99th percentile variation in height and body mass index (BMI) of the USA adult population. We developed a coupled turbulent flow and convective heat transfer simulation and benchmarked it in the 0.5 to 2.5 m·s-1 air speed range against prior literature. We computed the overall heat transfer coefficients, hoverall, for the manikins for representative airflow with 2 m·s-1 uniform speed and 5% turbulence intensity. We found that hoverall varied only between 19.9 and 23.2 W·m-2 K-1. Within this small range, the height of the manikins had negligible impact while an increase in the BMI led to a nearly linear decrease of the hoverall. Evaluation of the local coefficients revealed that those also nearly linearly decreased with BMI, which correlated to an inversely proportional local area (i.e., cross-sectional dimension) increase. Since even the most considerable difference that exists between 1st and 99th percentile BMI manikins is less than 15% of hoverall of the average manikin, it can be concluded that the impact of the human body shape on the convective heat transfer is minor.

Mechanism of Oil-in-Liquid Metal Emulsion Formation.

Gallium-based liquid metals (LMs) combine metallic properties with the deformability of a liquid, which makes them promising candidates for a variety of applications. To broaden the range of physical and chemical properties, a variety of solid additives have been incorporated into the LMs in the literature. In contrast, only a handful of secondary fluids have been incorporated into LMs to create foams (gas-in-LM) or emulsions (liquid-in-LM). LM foams readily form through mixing of LM in air, facilitated by the formation of a native oxide on the LM. In contrast, LM breaks up into microdroplets when mixed with a secondary liquid such as silicone oil. Stable silicone oil-in-LM emulsions form only during mixing of the oil with LM foam. In this work, we investigate the fundamental mechanism underlying this process. We describe two possible microscale mechanisms for emulsion formation: (1) oil replacing air in the foam or (2) oil creating additional features in the foam. The associated foam-to-emulsion density difference demonstrates that emulsions predominantly form through the addition of oxide-covered silicone oil capsules to the LM foam. We demonstrate this through density and surface wettability measurements and multiscale imaging of LM foam mixed with varied silicone oil contents in air or nitrogen environments. We also demonstrate the presence of a continuous silicone oil film on the emulsion surface and that this oil film prevents the embrittlement of contacting aluminum.

Human body radiation area factors for diverse adult population.

Radiation accounts for a significant fraction of the human body and environment heat exchange and strongly impacts thermal comfort and safety. The direct radiative exchange between an individual and a source or sink can be quantified using the effective (feff) and projected radiation area factors (fp). However, these factors have not been quantified for half of the population of the USA with an above-average body mass index (BMI). Here, we address this gap by developing thirty male and thirty female computational manikin models that cover the 1 to 99 percentile variation in height and BMI of adults in the USA. The radiative simulations reveal that the feff and the fp angular distributions are nearly independent of gender, height, and BMI. Appreciable relative differences from the average models only emerge for manikins with BMI above 80th percentile. However, these differences only occur at low zenith angles and, in absolute terms, are small as compared to variations induced by, for example, the zenith angle increase. We also use the manikin set to evaluate whether the body shape impacts the quality of human representation with several levels of geometrical simplification. We find that the "box/peg" body representation, which is based on the hemispherical fp average, is independent of the body shape. In turn, the fp distributions averaged over the azimuth angle range, representing the rotationally symmetric humans, are only impacted to the same degree as for the anatomical manikins. We also show that the anatomical manikins can be closely approximated by the multi-cylinder and sphere representation, at least from a radiation perspective. The developed anatomical manikin set is freely available and can be used to compute how body shape impacts a variety of external heat transport processes.

Analysis of thermocouple-based finger contact temperature measurements.

The understanding of heat conduction during finger contact with cooler or hotter objects is important for designing many electronic devices and for setting safety standards in a variety of occupational settings. In the most common experimental approach to study this process, a micro-thermocouple is placed at the finger-object interface. However, the interpretation of what this measurement corresponds to is not clear. To this end, we develop a three-dimensional thermal simulation of the finger-thermocouple-substrate configuration. The model predictions match finger cooling measurements in eight distinct cases available in prior literature (finger pressed with 1 N or 9.8 N against a steel or an aluminum substrate held at -2 °C or -10 °C). We demonstrate that the thermocouple can be represented accurately as a truncated sphere with emerging cylindrical wires while a multilayer block model of the finger provides similar results to an anatomically representative model. Our simulations show that in the eight previously studied cooling cases, the average surface temperature of skin that is in contact with the substrate follows nearly the same but offset cooling trend as the thermocouple tip temperature. The value of the offset is predominantly determined by the substrate material, with the thermocouple tip temperature being lower than the average skin surface temperature by 1-5 °C and 3-10 °C for steel and aluminum substrate cases, respectively. This temperature difference results in a moderate to an extreme thermocouple underprediction of the time necessary for the skin surface to reach the experimental safety threshold of 1 °C. Consequently, from the perspective of the safety related applications the thermocouple measurement provides a conservative limit on the contact duration and thereby is suitable for such purposes, but for applications requiring accurate skin temperature measurements alternative experimental approaches should be used.

Gallium oxide-stabilized oil in liquid metal emulsions.

Gallium based liquid metals (LM) have prospective biomedical, stretchable electronics, soft robotics, and energy storage applications, and are being widely adopted as thermal interface materials. The danger of gallium corroding most metals used in microelectronics requires the cumbersome addition of "barrier" layers or LM break-up into droplets within an inert matrix such as silicone oil. Such LM-in-oil emulsions are stabilized by native oxide on the droplets but have decreased thermal performance. Here we show that mixing of the silicone oil into an LM-air foam yields emulsions with inverted phases. We investigate the stability of these oil-in-LM emulsions through a range of processing times and oil viscosities, and characterize the impact of these parameters on the materials' structure and thermal property relationships. We demonstrate that the emulsion with 40 vol% of 10 cSt silicone oil provides a unique thermal management material with a 10 W m-1 K-1 thermal conductivity and an exterior lubricant thin film that completely prevents corrosion of contacting aluminum.

Improved methods for estimating mean radiant temperature in hot and sunny outdoor settings.

Thermal comfort research has utilized various sensors and models to estimate the mean radiant temperature (MRT) experienced by a human, including the standard black globe thermometer (SGT), acrylic globe thermometers (AGT), and cylindrical radiation thermometers (CRT). Rather than directly measuring radiation, a temperature is measured in the center of these low-cost sensors that can be related to MRT after theoretically accounting for convection. However, these sensors have not been systematically tested under long-term hot and clear conditions. Further, under variable weather conditions, many issues can arise due to slow response times, shape, inaccuracies in material properties and assumptions, and color (albedo, emissivity) inconsistencies. Here, we assess the performance of MRT produced by various heat transfer models, with and without new average surface temperature ([Formula: see text]) correction factors, using five instruments-the SGT (15 cm, black), tan and black CRTs, gray and black 38 mm AGTs-compared to 3D integral radiation measurements. Measurements were taken on an unobscured roof throughout summer-to-early-fall months in Tempe, Arizona, examining 58 full-sun days. Deviations without correcting for asymmetrical surface heating-found to be the main cause of errors-reached ± 15-20 °C MRT. By accounting for asymmetric heating through [Formula: see text] calculations, new corrective algorithms were derived for the low-cost sensor models. Results show significant improvements in the estimated MRT error for each sensor (i.e., ∆MRTmodel - IRM) when applying the [Formula: see text] corrections. The tan MRTCRT improved from 1.9 ± 6.2 to -0.1 ± 4.4 °C, while the gray AGT and SGT showed improvements from -1.6 ± 7.2 to -0.4 ± 6.3 °C and - 6.6 ± 6.4 to - 0.03 ± 5.7 °C, respectively. The new corrections also eliminated dependence on other meteorological factors (zenith, wind speed). From these results, we provide three simple equations for CRT, AGT, and SGT correction for future research use under warm-hot and clear conditions. This study is the most comprehensive empirical assessment of various low-cost instruments with broad applicability in urban climate and biometeorological research.

Fundamentals of soft thermofluidic system design.

The soft composition of many natural thermofluidic systems allows them to effectively move heat and control its transfer rate by dynamically changing shape (e.g. dilation or constriction of capillaries near our skin). So far, making analogous deformable "soft thermofluidic systems" has been limited by the low thermal conductivity of materials with suitable mechanical properties. By remaining soft and stretchable despite the addition of filler, elastomer composites with thermal conductivity enhanced by liquid-metal micro-droplets provide an ideal material for this application. In this work, we use these materials to develop an elementary thermofluidic system consisting of a soft, heat generating pipe that is internally cooled with flow of water and explore its thermal behavior as it undergoes large shape change. The transient device shape change invalidates many conventional assumptions employed in thermal design making analysis of this devices' operation a non-trivial undertaking. To this end, using time scale analysis we demonstrate when the conventional assumptions break down and highlight conditions under which the quasi-static assumption is applicable. In this gradual shape modulation regime the actuated devices' thermal behavior at a given stretch approaches that of a static device with equivalent geometry. We validate this time scale analysis by experimentally characterizing thermo-fluidic behavior of our soft system as it undergoes axial periodic extension-retraction at varying frequencies during operation. By doing so we explore multiple shape modulation regimes and provide a theoretical foundation to be used in the design of soft thermofluidic systems undergoing transient deformation.

Oxide-mediated mechanisms of gallium foam generation and stabilization during shear mixing in air.

Foaming of gallium-based liquid metals improves their processability and-seemingly in contrast to processing of other metal foams-can be achieved through shear-mixing in air without addition of solid microparticles. Resolving this discrepancy, systematic processing-structure-property characterization demonstrates that many crumpled oxide particles are generated prior to air bubble accumulation.

Droplet-train induced spatiotemporal swelling regimes in elastomers.

In this work, we perform a combined experimental and numerical analysis of elastomer swelling dynamics upon impingement of a train of solvent droplets. We use time scale analysis to identify spatiotemporal regimes resulting in distinct boundary conditions that occur based on relative values of the absorption timescale and the droplet train period. We recognize that when either timescale is significantly larger than the other, two cases of quasi-uniform swelling occur. In contrast, when the two timescales are comparable, a variety of temporary geometrical features due to localized swelling are observed. We show that the swelling feature and its temporal evolution depends upon geometric scaling of polymer thickness and width relative to the droplet size. Based on this scaling, we identify six cases of localized swelling and experimentally demonstrate the swelling features for two cases representing limits of thickness and width. A finite element model of local swelling is developed and validated with experimental results for these two cases. The model is subsequently used to explore the swelling behavior in the rest of the identified cases. We show that depending upon the lateral dimension of the sample, swelling can locally exhibit mushroom, mesa, and cap like shapes. These deformations are magnified during the droplet-train impact but dissipate during post-train polymer equilibration. Our results also show that while swelling shape is a function of lateral dimensions of the sample, the extent of swelling increases with the elastomer sample thickness.

Dropwise Condensation on Soft Hydrophobic Coatings.

Promoting dropwise condensation (DWC) could improve the efficiency of many industrial systems. Consequently, a lot of effort has been dedicated to finding durable materials that could sustainably promote DWC as well as finding routes to enhance the heat transfer rate during this phase change process. Motivated by previous reports of substrate softening increasing droplet nucleation rate, here we investigated how mechanical properties of a substrate impact relevant droplet-surface interactions and DWC heat transfer rate. Specifically, we experimentally quantified the effect of hydrophobic elastomer's shear modulus on droplet nucleation density and shedding radius. To quantify the impact of substrate softening on heat transfer through individual droplets, we combined analytical solution of elastomer deformation induced by droplets with finite element modeling of the heat transfer process. By substituting these experimentally and theoretically derived values into DWC heat transfer model, we quantified the compounding effect of the substrate's mechanical properties on the overall heat transfer rate. Our results show that softening of the substrates below a shear modulus of 500 kPa results in a significant reduction in the condensation heat transfer rate. This trend is primarily driven by additional thermal resistance of the liquid posed by depression of the soft substrate.

Microscale Mechanism of Age Dependent Wetting Properties of Prickly Pear Cacti (Opuntia).

Cacti thrive in xeric environments through specialized water storage and collection tactics such as a shallow, widespread root system that maximizes rainwater absorption and spines adapted for fog droplet collection. However, in many cacti, the epidermis, not the spines, dominates the exterior surface area. Yet, little attention has been dedicated to studying interactions of the cactus epidermis with water drops. Surprisingly, the epidermis of plants in the genus Opuntia, also known as prickly pear cacti, has water-repelling characteristics. In this work, we report that surface properties of cladodes of 25 taxa of Opuntia grown in an arid Sonoran climate switch from water-repelling to superwetting under water impact over the span of a single season. We show that the old cladode surfaces are not superhydrophilic, but have nearly vanishing receding contact angle. We study water drop interactions with, as well as nano/microscale topology and chemistry of, the new and old cladodes of two Opuntia species and use this information to uncover the microscopic mechanism underlying this phenomenon. We demonstrate that composition of extracted wax and its contact angle do not change significantly with time. Instead, we show that the reported age dependent wetting behavior primarily stems from pinning of the receding contact line along multilayer surface microcracks in the epicuticular wax that expose the underlying highly hydrophilic layers.

Surface and Wetting Properties of Embiopteran (Webspinner) Nanofiber Silk.

Insects of the order Embioptera, known as embiopterans, embiids, or webspinners, weave silk fibers together into sheets to make shelters called galleries. In this study, we show that silk galleries produced by the embiopteran Antipaluria urichi exhibit a highly hydrophobic wetting state with high water adhesion macroscopically equivalent to the rose petal effect. Specifically, the silk sheets have advancing contact angles above 150°, but receding contact angle approaching 0°. The silk sheets consist of layered fiber bundles with single strands spaced by microscale gaps. Scanning and transmission electron microscopy (SEM, TEM) images of silk treated with organic solvent and gas chromatography mass spectrometry (GC-MS) of the organic extract support the presence of a lipid outer layer on the silk fibers. We use cryogenic SEM to demonstrate that water drops reside on only the first layer of the silk fibers. The area fraction of this sparse outer silk layers is 0.1 to 0.3, which according to the Cassie-Baxter equation yields an effective static contact angle of ∼130° even for a mildly hydrophobic lipid coating. Using high magnification optical imaging of the three phase contact line of a water droplet receding from the silk sheet, we show that the high adhesion of the drop stems from water pinning along bundles of multiple silk fibers. The bundles likely form when the drop contact line is pinned on individual fibers and pulls them together as it recedes. The dynamic reorganization of the silk sheets during the droplet movement leads to formation of "super-pinning sites" that give embiopteran silk one of the strongest adhesions to water of any natural hydrophobic surface.

Inhibition of Condensation Frosting by Arrays of Hygroscopic Antifreeze Drops.

The formation of frost and ice can have negative impacts on travel and a variety of industrial processes and is typically addressed by dispensing antifreeze substances such as salts and glycols. Despite the popularity of this anti-icing approach, some of the intricate underlying physical mechanisms are just being unraveled. For example, recent studies have shown that in addition to suppressing ice formation within its own volume, an individual salt saturated water microdroplet forms a region of inhibited condensation and condensation frosting (RIC) in its surrounding area. This occurs because salt saturated water, like most antifreeze substances, is hygroscopic and has water vapor pressure at its surface lower than water saturation pressure at the substrate. Here, we demonstrate that for macroscopic drops of propylene glycol and salt saturated water, the absolute RIC size can remain essentially unchanged for several hours. Utilizing this observation, we demonstrate that frost formation can be completely inhibited in-between microscopic and macroscopic arrays of propylene glycol and salt saturated water drops with spacing (S) smaller than twice the radius of the RIC (δ). Furthermore, by characterizing condensation frosting dynamics around various hygroscopic drop arrays, we demonstrate that they can delay complete frosting over of the samples 1.6 to 10 times longer than films of the liquids with equivalent volume. The significant delay in onset of ice nucleation achieved by dispensing propylene glycol in drops rather than in films is likely due to uniform dilution of the drops driven by thermocapillary flow. This transport mode is absent in the films, leading to faster dilution, and with that facilitated homogeneous nucleation, near the liquid-air interface.

Inverted Leidenfrost-like Effect during Condensation.

Water droplets condensing on solidified phase change materials such as benzene and cyclohexane near their melting point show in-plane jumping and continuous "crawling" motion. The jumping drop motion has been tentatively explained as an outcome of melting and refreezing of the materials surface beneath the droplets and can be thus considered as an inverted Leidenfrost-like effect (in the classical case vapor is generated from a droplet on a hot substrate). We present here a detailed investigation of jumping movements using high-speed imaging and static cross-sectional cryogenic focused ion beam scanning electron microscope imaging. Our results show that drop motion is induced by a thermocapillary (Marangoni) effect. The in-plane jumping motion can be delineated to occur in two stages. The first stage occurs on a millisecond time scale and comprises melting the substrate due to drop condensation. This results in droplet depinning, partial spreading, and thermocapillary movement until freezing of the cyclohexane film. The second stage occurs on a second time scale and comprises relaxation motion of the drop contact line (change in drop contact radius and contact angle) after substrate freezing. When the cyclohexane film cannot freeze, the droplet continuously glides on the surface, resulting in the crawling motion.

How droplets nucleate and grow on liquids and liquid impregnated surfaces.

Condensation on liquids has been studied extensively in context of breath figure templating, materials synthesis and enhancing heat transfer using liquid impregnated surfaces. However, the mechanics of nucleation and growth on liquids remains unclear, especially on liquids that spread on the condensate. By examining the energy barriers of nucleation, we provide a framework to choose liquids that can lead to enhanced nucleation. We show that due to limits of vapor sorption within a liquid, nucleation is most favoured at the liquid-air interface and demonstrate that on spreading liquids, droplet submergence within the liquid occurs thereafter. We provide a direct visualization of the thin liquid profile that cloaks the condensed droplet on a liquid impregnated surface and elucidate the vapour transport mechanism in the liquid films. Finally, we show that although the viscosity of the liquid does not affect droplet nucleation, it plays a crucial role in droplet growth.

Different shades of oxide: from nanoscale wetting mechanisms to contact printing of gallium-based liquid metals.

Gallium-based liquid metals are of interest for a variety of applications including flexible electronics, soft robotics, and biomedical devices. Still, nano- to microscale device fabrication with these materials is challenging because, despite having surface tension 10 times higher than water, they strongly adhere to a majority of substrates. This unusually high adhesion is attributed to the formation of a thin oxide shell; however, its role in the adhesion process has not yet been established. In this work, we demonstrate that, dependent on dynamics of formation and resulting morphology of the liquid metal-substrate interface, GaInSn adhesion can occur in two modes. The first mode occurs when the oxide shell is not ruptured as it makes contact with the substrate. Because of the nanoscale topology of the oxide surface, this mode results in minimal adhesion between the liquid metal and most solids, regardless of substrate's surface energy or texture. In the second mode, the formation of the GaInSn-substrate interface involves rupturing of the original oxide skin and formation of a composite interface that includes contact between the substrate and pieces of old oxide, bare liquid metal, and new oxide. We demonstrate that in this latter mode GaInSn adhesion is dominated by the intimate contact between new oxide and substrate. We also show that by varying the pinned contact line length using varied degrees of surface texturing, the adhesion of GaInSn in this mode can be either decreased or increased. Lastly, we demonstrate how these two adhesion modes limit microcontact printing of GaInSn patterns but can be exploited to repeatedly print individual sub-200 nm liquid metal drops.

Dynamics of silver nanoparticle release from wound dressings revealed via in situ nanoscale imaging.

The use of silver nanoparticles (AgNPs) in textiles for enhanced anti-microbial properties has led to concern about their release and impact on both human and environmental health. Here a novel method for in situ visualization of AgNP release from silver-impregnated wound dressings is introduced. By combining an environmental scanning electron microscope, a gaseous analytical detector and a peltier cooling stage, this technique provides near-instantaneous nanoscale characterization of interactions between individual water droplets and AgNPs. We show that dressings with different silver application methods have very distinct AgNP release dynamics. Specifically, water condensation on dressings with AgNP deposited directly on the fiber surface resulted in substantial and rapid AgNP release. By comparison, AgNP release from wound dressing with nanoparticles grown, not deposited, from the fiber surface was either much slower or negligible. Our methodology complements standard bulk techniques for studying of silver release from fabrics by providing dynamic nanoscale information about mechanisms governing AgNP release from individual fibers. Thus coupling these nano and macro-scale methods can provide insight into how the wound dressing fabrication could be engineered to optimize AgNP release for different applications.

Particle-assisted formation of oil-in-liquid metal emulsions.

Gallium-based liquid metals (LMs) have surface tension an order of magnitude higher than water and break up into micro-droplets when mixed with other liquids. In contrast, silicone oil readily mixes into LM foams to create oil-in-LM emulsions with oil inclusions. Previously, the LM was foamed through rapid mixing in air for an extended duration (over 2 h). This process first results in the internalization of oxide flakes that form at the air-liquid interface. Once a critical fraction of these randomly shaped solid flakes is reached, air bubbles internalize into the LM to create foams that can internalize secondary liquids. Here, we introduce an alternative oil-in-LM emulsion fabrication method that relies on the prior addition of SiO2micro-particles into the LM before mixing it with the silicone oil. This particle-assisted emulsion formation process provides a higher control over the composition of the LM-particle mixture before oil addition, which we employ to systematically study the impact of particle characteristics and content on the emulsions' composition and properties. We demonstrate that the solid particle size (0.8µm to 5µm) and volume fraction (1%-10%) have a negligible impact on the internalization of the oil inclusions. The inclusions are mostly spherical with diameters of 20-100µm diameter and are internalized by forming new, rather than filling old, geometrical features. We also study the impact of the particle characteristics on the two key properties related to the functional application of the LM emulsions in the thermal management of microelectronics. In particular, we measure the impact of particles and silicone oil on the emulsion's thermal conductivity and its ability to prevent deleterious gallium-induced corrosion and embrittlement of contacting metal substrates.