Cooling of electric machines is a key to increasing power density and improving reliability. This paper focuses on the design of a machine using a cooling jacket wrapped around the stator. The thermal contact resistance (TCR) between the electric machine stator and cooling jacket is a significant factor in overall performance and is not well characterized. This interface is typically an interference fit subject to compressive pressure exceeding 5 MPa. An experimental investigation of this interface was carried out using a thermal transmittance setup using pressures between 5 and 10 MPa. The results were compared to currently available models for contact resistance, and one model was adapted for prediction of TCR in future motor designs.

Introduction

Cooling of electric machines is an important factor in increasing power density. Components within an electric machine such as the magnets and wire insulation experience reduced performance or reliability when exposed to temperatures beyond their thermal specifications [1]. Therefore, heat must be removed from these components to limit the temperature to which they are exposed. Typically, heat must flow through several components and material interfaces or contacts before the heat can be extracted from the motor through convective cooling. For this reason, the convective cooling approach for the machine can impact the path of heat transfer through the electric machine, resulting in certain materials or interfaces having a larger impact on the overall thermal management of the electric machine. Examples of some thermal management approaches that have been applied to electric machines are summarized in Ref. [2].

The choice of convective cooling technology depends on the electric machine's performance requirements, intended application, and available coolants. The research described in this paper focused on a new electric machine design cooled with a high-performance cooling jacket integrated into the machine case using water–ethylene glycol, as shown in Fig. 1.

Fig. 1
Electric machine cross section highlighting stator-to-case contact
Fig. 1
Electric machine cross section highlighting stator-to-case contact
Close modal

Because the machine of interest used a high-performance cooling jacket for heat removal, the heat flow path through the machine was predominantly radial. A parameter sensitivity study was performed on a preliminary design of the electric machine using a thermal finite element analysis model to determine which thermal resistances were most significant. To perform the study, the thermal properties of each material or interface were decreased by 20%, and the change in maximum temperature was reported. The thermal sensitivity analysis indicated that the dominant thermal resistances for heat conduction through the machine were the stator lamination stack and the thermal contact resistance (TCR) between the stator and case [3]. For this reason, there was a need to experimentally measure the in-plane lamination thermal conductivity and the TCR of the stator-to-case interface.

The coolant jacket for the electric machine was mounted to the stator through an interference fit between the coolant jacket and stator. For the machine of interest, the stator-to-case interference fit resulted in a compressive pressure on the joint ranging from 5.52 to 9.65 MPa. The laminated silicon steel used in the stator gives the outer diameter of the stator a ridged surface. The inner diameter of the case in contact with the stator is a machined aluminum surface. Examples from a similar electric machine are shown in Fig. 2, showing the stator's outer diameter surface and the cooling jacket's inner diameter surface.

Fig. 2
Top: machine stator surface. Bottom: case interior surface.
Fig. 2
Top: machine stator surface. Bottom: case interior surface.
Close modal

There are a number of models and experimental data for calculating TCR. The simplest use an effective air gap with suggestions for values based on surface roughness [4,5]. Using these techniques yields a pressure-independent TCR ranging from 840 to 1400 mm2 K/W. Theoretical models will take into account factors such as the pressure on the joint, the hardness of the two materials, the surface roughness, and interstitial gas [6,7]. Using such a model gives a pressure-dependent TCR ranging from 354 to 244 mm2 K/W for the 5.52–9.65 MPa pressure range. The disadvantage of such a model is the requirement to obtain upward of 15 parameters, some of which may be obscure or difficult to obtain or require specialized equipment to measure. Finally, extrapolating from previously published experimental data [8] gives a TCR of 157–149 mm2 K/W for the same pressure range.

The large range in the calculated TCR values (149–1400 mm2 K/W) causes a challenge in the electric machine thermal design. A TCR of 150 mm2 K/W corresponds to a temperature drop across the stator-to-jacket interface of 5 °C–10 °C. A TCR of 1400 mm2 K/W across the same interface results in a temperature drop about seven times higher. Therefore, obtaining an accurate value for the stator-to-case TCR is critical in the thermal design of case-cooled electric machines.

This paper presents an experimental method to measure thermal contact resistance at pressures representing interference fits. Pressures from 5.52 to 9.65 MPa and two types of commonly used machine laminations were studied. The experimental data were used to select a model that can be used with easily obtained parameters, while giving an accurate result that can be used in relevant electric machine thermal models.

Experimental Apparatus and Procedures

To measure the TCR of the case-to-stator interference fit, a high-pressure transmittance setup was constructed in accordance with ASTM International (ASTM) Method D5470 [9] using a high-capacity hydraulic clamp, as shown in Fig. 3 [10]. The setup includes a lamination coupon with the laminations stacked horizontally. The laminations are sandwiched between two 4-mm-thick aluminum contact plates that simulate the aluminum case. The aluminum contact plates on the top and bottom of the fixture are pressed against copper metering blocks with thermal grease applied at the interface. The setup was validated against both an existing thermal transmittance setup at the National Renewable Energy Laboratory for lower pressures described by Narumanchi et al. [11] and a xenon flash diffusivity measurement system using the 4-mm-thick aluminum contact plates to be used in the experiment. The setup uses resistance temperature detectors (RTDs) placed into the copper metering blocks located on the top and bottom of the test sample. Heat is provided to one side of the sample through electrical heaters inserted into an aluminum heater block. The test apparatus is thermally isolated from the large pneumatic press using a cold plate on the bottom. Both cold plates are cooled with a bath circulator using plain water.

Fig. 3
High-pressure thermal transmittance setup showing cutaway view of sample stabilizing rig
Fig. 3
High-pressure thermal transmittance setup showing cutaway view of sample stabilizing rig
Close modal

For the experiments, the laminations were manufactured with a tolerance of 0.1 mm. A stabilizing rig was built to hold the laminations together. The stabilizing rig provided a method to hold the laminations level without welding or adhesion such that the surface was “self-leveling.” The stabilizing rig consisted of an aluminum frame with thermally insulating Teflon bumpers held against the lamination stack with set screws.

The objective of the experiment was to measure the TCR between the steel laminations and the aluminum contact plate in addition to the lamination thermal conductivity parallel to the in-plane direction of the laminations. The primary test factors for the experiment included pressure and lamination material. The sample temperature ranged from 45 °C to 115 °C. A temperature gradient must be applied to the sample to obtain results using the thermal transmittance method, and the temperature was limited by the safe operating temperature of the heaters in the heater block shown in Fig. 3. The highest safe temperature was used to maximize the heat flux across the sample to improve the accuracy of the measurement results. Electric machines will operate across a wide range of temperatures, but the impact of temperature on the thermal contact resistance was not included as part of the design space of this experiment.

Test Pressures.

The pressures applied during the tests were intended to represent an interference fit between the case and stator where the stator diameter is larger than the inner diameter of the case. The pressure was calculated based on the geometry of the UQM prototype electric machine. The tested pressures covered a range of values that encompassed the pressure estimate of the UQM prototype machine, which included 5.52 MPa, 6.89 MPa, 8.27 MPa, and 9.65 MPa. Each data point at the listed pressure was repeated three times to aid in the measurement uncertainty analysis. Thermal contact resistance is affected by hysteresis or the loading history of the contact. Thermal cycling, load cycling, and extended time under load can reduce thermal contact resistance [12]. To minimize hysteresis effects, measurements were taken immediately after thermal steady-state was reached. Steady-state was defined by a maximum temperature change of less than 0.03 °C for a 10-min period. Pressures were only increased throughout each test (i.e., a higher-pressure test could not precede a test at lower pressure).

The high pressures involved necessitated investigation into whether the surface of the contact plate would become indented from the lamination coupon surface. This was done by using a laser profilometer to scan the surface and then using spatial two-dimensional fast Fourier transform to transform the surface data into the frequency domain. Analyzing the surface in the frequency domain allows differentiation between tool marks and indentation caused by the lamination coupon surface. If the experiment caused an indentation in the aluminum, it would necessitate replacing the aluminum for each experimental test to ensure repeatability. Whether or not the aluminum plates could be reused affects the number of experiments that can be performed due to the potential expense of having to cut or refinish plates for each repetition should they not be reusable. The laminations were also tested for possible work hardening by comparing thermal properties before and after application of pressure.

Lamination Materials.

Two types of laminations were used: 29-gauge M15 material and 0.2-mm high-frequency material from JFE Steel Corporation. Both laminations were laser cut to get the shapes needed. The different edges are the result of different feed rates and/or intensity used by the vendor in the cutting process to ensure samples were not damaged due to thinner laminations being more brittle. Figure 4 shows the edges of the two lamination materials included in the experimental measurements. The JFE material is thinner, and the resulting edge has a serrated appearance.

Fig. 4
Edge views of lamination materials. The microscope stage is labeled for clarity.
Fig. 4
Edge views of lamination materials. The microscope stage is labeled for clarity.
Close modal

Surface Properties of Contact Plates.

Two surfaces for the contact plates were tested and compared to determine the significance of the surface finish. The first was a factory-ground surface with no special treatment with an average surface roughness of 0.6 μm. In an actual machine case, the inner surface is lathed to a specified tolerance. Ideally, the second surface would exactly replicate a lathed machine pattern. However, it was not possible to exactly duplicate a lathed machine pattern on a flat surface without specialized equipment. Instead, the second surface approximated a lathed surface by machining the aluminum contact block surface with a fly cut. Fly cutting results in surface tooling (ridges) similar to what would be obtained with lathing; using a very large bit results in nearly straight ridges. The feed speed gives control over the surface roughness and was adjusted to match the surface roughness of the stator case, which was 1.6 μm. Figure 5 shows a surface topography map of the fly cut surface obtained using a laser profilometer. The data were postprocessed in matlab to obtain the area-average surface roughness.

Fig. 5
Ten-mm square sample area surface profile of the contact plate
Fig. 5
Ten-mm square sample area surface profile of the contact plate
Close modal

Data Analysis Procedure.

The test apparatus measures the thermal resistance of the total stack-up between the two metering blocks shown in Fig. 3. The stack-up can be represented by Eq. (1). The TCR is obtained by subtracting the thermal resistances due to the other layers.
(1)

RGL represents the thermal grease layer between the copper metering blocks and the contact plates shown in Fig. 3. The grease layer improves the thermal contact between the copper metering block and contact plate. RCP represents the thermal resistance of the aluminum contact plate shown in Fig. 3. The thermal resistance of the contact plate was characterized using a xenon flash transient technique. RLS represents the thermal resistance of the lamination stack coupon shown in Fig. 3.

The challenge with the thermal transmittance setup is that the measurement uncertainty is a fixed percentage of Rtot. Therefore, care must be taken to minimize the subtractions from Rtot to obtain the value of interest, in this case, TCR. Care must also be taken to minimize variations in Rtot due to variations of individual layers (e.g., grease layer thickness).

The grease layer between the metering block and contact plate was characterized at all four pressures using only a contact plate and obtained using Eq. (2). Using a contact plate instead of simply characterizing the grease layer independently simulates the experimental conditions of the grease layer exactly between the copper metering block and the aluminum contact plate
(2)
While estimates and measurements for the bulk lamination material thermal conductivity were available, the direction-dependent thermal conductivity of the lamination stack parallel with the laminations was not known. For this reason, it was not possible to subtract the thermal resistance of the lamination stack as shown in Eq. (1) with sufficient confidence. To overcome this challenge, three sizes of laminations were used to generate a thermal resistance versus sample height curve similar to what is shown in Fig. 6. By fitting a first-order curve, the lamination stack's thermal resistance can be extracted independently of other parameters using Eq. (3), and TCR can be extracted from the stack using Eq. (4)
(3)
(4)
Fig. 6
Thermal resistance as a function of lamination coupon thickness
Fig. 6
Thermal resistance as a function of lamination coupon thickness
Close modal

Results and Discussion

The experimental results are described in the following sections. To validate the experimental design and the test procedure, a series of preliminary tests was performed. The results of the preliminary testing provided confidence in the experimental approach. Once the experimental design was finalized, the experimental results were recorded as described in the following.

Preliminary Test Results.

Before the high-pressure tests were conducted, an initial test was performed on the existing thermal transmittance setup to check if the random uncertainty of the measurement would be acceptable before constructing the high-pressure setup. The test yielded an acceptable uncertainty and work proceeded to build the high-pressure experimental setup.

The thermal resistance of the grease layer in Eq. (4) was determined experimentally as described previously. The results of the measurements are shown in Table 1. The result is the average of three experiments. The values are small compared to the expected measurements. It is worth noting that the thermal resistance decreases with pressure, which follows the expected behavior.

Table 1

Meter block-to-contact plate thermal resistance with thermal grease interface

Test pressure (MPa)Grease thermal resistance (mm2 K/W)
5.528.4
6.898.1
8.277.8
9.657.6
Test pressure (MPa)Grease thermal resistance (mm2 K/W)
5.528.4
6.898.1
8.277.8
9.657.6

To test for indentations in the contact plates, the experiment was performed at 9.65 MPa for the same temperature and time conditions used in the experiment. The surface was then analyzed using a spatial two-dimensional fast Fourier transform technique. Transforming the surface data into spatial frequency space allows simple differentiation between a potential indentation due to a lamination stack and existing tool marks that both run in the same direction. There was no indication of any indentation or change in surface properties due to the experiment. Samples were also periodically analyzed throughout the experiment to ensure that no cumulative damage occurred to the plate surface that would affect results. In addition, the laminations were tested on the xenon flash before and after the preliminary test to ensure no thermal property changes due to potential work hardening had occurred.

To highlight the effect the lathed surface on the interior of the case had on the contact resistance, a series of tests was run using a smooth contact plate and the fly-cut contact plate. The results are shown in Fig. 7. For the TCR calculation in Fig. 7, a value for the lamination coupon was calculated based on the known bulk thermal properties of the laminations as a preliminary estimate. As seen in Fig. 7, the surface finish significantly impacts the measured contact resistance. For this reason, the tests proceeded with the fly-cut surface to more closely represent the lathed surface properties.

Fig. 7
Comparison of contact plate surface finish effect on TCR
Fig. 7
Comparison of contact plate surface finish effect on TCR
Close modal

Thermal Property Experimental Results.

Figure 8 shows a plot of the data for the M15 29-gauge lamination coupons at the lowest of the four pressures tested. Three data points were collected for each of the three coupon thicknesses. The data were fit with a first-order curve with the curve fit data points weighted according to the systematic error of the experiment. The thermal properties for both the lamination stack and the stator-to-case contact resistance were extracted using the slope and intercept of the fit, respectively, using Eqs. (3) and (4).

Fig. 8
Total thermal resistance measurements for M15 29-gauge coupons at 5.52 MPa with extrapolation to 0 lamination coupon thickness
Fig. 8
Total thermal resistance measurements for M15 29-gauge coupons at 5.52 MPa with extrapolation to 0 lamination coupon thickness
Close modal

As described previously, the thermal resistances of the grease and contact plate are subtracted from the intercept to obtain the contact resistance. To calculate the uncertainty of the measurements, the curve fit is weighted on systematic error, which includes RTD calibration error, estimated spatial temperature variation calculated using a thermal finite element model of the setup, and RTD location error. The RTDs were calibrated to a National Institute of Standards and Technology traceable reference probe using a static temperature calibration bath. The curve fit yields random error and the 95% confidence interval (U95) is calculated [12,13]. Figure 9 summarizes the results for both the stator-to-case contact and the bulk thermal conductivity for the laminations as calculated from the slope-intercept technique.

Fig. 9
Top: stator-to-case TCR results. Bottom: lamination effective thermal conductivity results. Error bars represent the 95% confidence interval.
Fig. 9
Top: stator-to-case TCR results. Bottom: lamination effective thermal conductivity results. Error bars represent the 95% confidence interval.
Close modal

The results shown in Fig. 9 are generally consistent with expectations, which give confidence in the technique. The thermal conductivity of the laminations does not vary with pressure and is near the bulk value measured with the xenon flash technique of 22 W/m K for both materials. The contact resistance decreases with increasing pressure. The contact resistance is significantly lower than what is predicted by effective air gap models, but agrees well with pressure-dependent models for TCR. The uncertainty in the results shown by the error bars in Fig. 9 precludes making further conclusions about the results. However, there is a clear trend of pressure dependence, and the different surface of the JFE as shown in Fig. 4 may impact TCR.

Theoretical Model.

The model described in Madhusudana [7] appeared to give the best agreement with the experimental results for TCR in the 5–10 MPa pressure range. The purpose here is not to propose a new model or a physically comprehensive model but to summarize the existing model as it pertains to a stator-to-case TCR. The model includes both solid [6] and fluid [14] components to calculate the inverse of TCR or thermal contact conductance (TCC)
(5)
The solid spot conductance (hs) is described by Eq. (6) as
(6)

where k is the harmonic mean of the thermal conductivities of the two materials; C and n are coefficients, which are 1.13 and 0.94, respectively [6]; σ is the centerline average roughness of the two contacting surfaces and is 1.25σRMS; σRMS is the root mean square of the average surface roughness of the two contacting surfaces; tan Θ is the mean slope of the surface asperities; P is the applied pressure; and H is the contact microhardness of the softer of the two materials. P and H have units of MPa.

The surface roughness of the case plate was 1.6 μm. The surface roughness of the laminations was determined using the same technique and yielded 10.9 μm and 11.8 μm for the M15 and JFE materials, respectively. Any roughness tester could be used to measure the surface roughness provided it has the range needed. Measuring tan Θ is nontrivial, and motor designers are unlikely to possess the necessary equipment. As an alternative, tan Θ can be calculated based on the surface roughness correlations [15]. However, the surface roughness is outside the valid range for the correlation to be valid. Furthermore, the correlation assumes a Gaussian, random surface, which is not necessarily the case for this contact. Because of this, tan Θ was determined empirically to be 0.12. Typical values for tan Θ range between 0.03 and 0.18 [7,15]. Note that bulk properties for k are used (22 W/m K for both types of laminations), not the effective properties reported in Fig. 9, which are a function of the stacking factor and include a small amount of interstitial air. The aluminum thermal conductivity was measured to be 195 W/m K. The aluminum microhardness was estimated to be 930 MPa, derived from typical Brinell hardness for Al 6061-T6. Surface microhardness is dependent on the history of the surface, including how it was finished, heat treatment, or work hardening. The actual surface microhardness may differ from the estimated value and requires specialized equipment to measure. The gaseous fluid conductance (hg) is described by Eq. (7) as
(7)
where kg is the thermal conductivity of the interstitial gas, in this case air at 80 °C. δ is calculated using Eq. (8) [14]
(8)

Figure 10 shows a comparison of the model to the data for both sets of laminations. The model shows good agreement with the data. Although not included in the plot, the model agrees with the preliminary test at 0.22 MPa. The range of 150–250 mm2 K/W is on the lower end of the estimates of 149–1400 mm2 K/W for the interface. More importantly, it reduces the uncertainty of the value by 90%, which leads to more accurate machine design.

Fig. 10
Top: Comparison of TCR model to M15 29-gauge data. Bottom: comparison of TCR model to JFE (0.2 mm) data. Root mean square surface roughness of the contacting surfaces used in the model is noted on the plot. Error bars indicate 95% confidence interval for the data.
Fig. 10
Top: Comparison of TCR model to M15 29-gauge data. Bottom: comparison of TCR model to JFE (0.2 mm) data. Root mean square surface roughness of the contacting surfaces used in the model is noted on the plot. Error bars indicate 95% confidence interval for the data.
Close modal

Thermal contact resistance calculated by the model agrees reasonably well (17% difference) with the experimental results reported by Kulkarni et al. [8] for 13.28 MPa, but the model results are lower than their results when extrapolating to higher pressures (21–39 MPa). Note that values were not given for the material and surface properties, so it is possible they differed from what is presented here.

Conclusions

The TCR between machine laminations and an aluminum cooling jacket under interference loading was measured using a simulated stack within a high-pressure ASTM D5470 setup. The results for two commonly used machine lamination types are presented and compared to existing contact resistance models. It was found that the TCR ranged from 150 to 250 mm2 K/W for interference fit pressures between 5.52 and 9.65 MPa. A model was presented with parameters that can be obtained using standard equipment and agrees closely with the experimental data.

Pressure-independent models (such as effective air gap models) are valid for pressures below 500 kPa. For pressures exceeding 500 kPa, the air gap model overestimates TCR proportional to the pressure. For the interference fit pressures investigated, the pressure-independent models overestimated TCR by a factor of four to seven compared to the pressure-dependent model presented. The model presented was validated to 10 MPa, but should be valid for any pressure well below the material yield point.

The stator-to-case TCR is a critical parameter in high-performance and compact machine designs, especially those that depend on case cooling. Without an accurate value for the stator-to-case TCR, the cooling system may be optimized incorrectly or additional expense invested in alternative cooling or oversizing the motor that may not be necessary. This work helps electric machine designers accurately estimate the stator-to-case TCR for a range of electric machines without having to resort to expensive experimental tests or overly conservative estimates.

Acknowledgment

We acknowledge financial support for the work provided by Susan Rogers and Steven Boyd, Technology Managers for the Electric Drive Technologies Program, Vehicle Technologies Office, U.S. Department of Energy, Office of Energy Efficiency and Renewable Energy (EERE). The significant contributions from Doug DeVoto (NREL) and Charlie King (NREL) to the project are acknowledged.

This work was supported by the U.S. Department of Energy under Contract No. DE-AC36-08GO28308 with the Alliance for Sustainable Energy, LLC, the Manager and Operator of the National Renewable Energy Laboratory.

Nomenclature

C =

constant

H =

surface microhardness

hg =

gaseous fluid conductance

hs =

solid spot conductance

k =

thermal conductivity

n =

exponential constant

P =

pressure

R =

thermal resistance

RTD =

resistance temperature detector

tanΘ =

mean asperity slope

TCC =

thermal contact conductance (1/TCR)

TCR =

thermal contact resistance

δ =

mean surface separation

σ =

surface roughness

Subscripts
CP =

contact plate

GL =

grease layer

LS =

lamination stack

RMS =

root-mean-square

Tot =

total

References

1.
Stone
,
G. C.
,
Boulter
,
E. A.
,
Culbert
,
I.
, and
Dhirani
,
H.
,
2004
,
Electrical Insulation for Rotating Machines: Design, Evaluation, Aging, Testing, and Repair
,
Wiley
, Hoboken, NJ.
2.
Popescu
,
M.
,
Staton
,
D. A.
,
Boglietti
,
A.
,
Cavagnino
,
A.
,
Hawkins
,
D.
, and
Goss
,
J.
,
2016
, “
Modern Heat Extraction Systems for Power Traction Machines—A Review
,”
IEEE Trans. Ind. Appl.
,
52
(
3
), pp.
2167
2175
.
3.
Lutz
,
J.
,
2014
, “Unique Lanthide-Free Motor Construction,” 2014 Annual Merit Review, Washington, DC, accessed Nov. 29, 2016, http://energy.gov/sites/prod/files/2014/07/f17/ape044_lutz_2014_o.pdf
4.
Lindström
,
J.
,
1999
,
Thermal Model of a Permanent-Magnet Motor for a Hybrid Electric Vehicle
,
Chalmers University of Technology
, Gothenburg, Sweden.
5.
Staton
,
D.
,
Boglietti
,
A.
, and
Cavagnino
,
A.
,
2005
, “
Solving the More Difficult Aspects of Electric Motor Thermal Analysis in Small and Medium Size Industrial Induction Motors
,”
IEEE Trans. Energy Convers.
,
20
(
3
), pp.
620
628
.
6.
Mikić
,
B. B.
,
1974
, “
Thermal Contact Conductance; Theoretical Considerations
,”
Int. J. Heat Mass Transfer
,
17
(
2
), pp.
205
214
.
7.
Madhusudana
,
C.
,
2014
,
Thermal Contact Conductance
, 2nd ed.,
Springer
, Cham,
Switzerland
.
8.
Kulkarni
,
D. P.
,
Rupertus
,
G.
, and
Chen
,
E.
,
2012
, “
Experimental Investigation of Contact Resistance for Water Cooled Jacket for Electric Motors and Generators
,”
IEEE Trans. Energy Convers.
,
27
(
1
), pp.
204
210
.
9.
ASTM D09 Committee
,
2012
, “Standard Test Method for Thermal Transmission Properties of Thermally Conductive Electrical Insulation Materials,” ASTM International, West Conshohocken, PA, Standard No.
ASTM D5470-06
.https://www.astm.org/DATABASE.CART/HISTORICAL/D5470-06.htm
10.
Ley
,
J.
,
2015
, “Unique Lanthide-Free Motor Construction,” 2015 Annual Merit Review, Washington, DC, accessed Nov. 29, 2016, http://energy.gov/sites/prod/files/2015/06/f24/edt044_gilbert_2015_o.pdf
11.
Narumanchi
,
S.
,
Mihalic
,
M.
,
Kelly
,
K.
, and
Eesley
,
G.
,
2008
, “
Thermal Interface Materials for Power Electronics Applications
,”
11th Intersociety Conference on Thermal and Thermomechanical Phenomena in Electronic Systems
(
ITHERM
), Orlando, FL, May 28–31, pp.
395
404
.
12.
Kirkup
,
L.
,
2002
,
Data Analysis With Excel®: an Introduction for Physical Scientists
,
Cambridge University Press
, Cambridge, UK.
13.
Dieck
,
R. H.
,
2007
,
Measurement Uncertainty: Methods and Applications
,
ISA
, Research Triangle Park, NC.
14.
Madhusudana
,
C. V.
, and
Fletcher
,
L. S.
,
1981
, “Gas Conductance Contribution to Contact Heat Transfer,”
AIAA
Paper No. 81-1163.
15.
Antonetti
,
V. W.
,
Whittle
,
T. D.
, and
Simons
,
R. E.
,
1993
, “
An Approximate Thermal Contact Conductance Correlation
,”
ASME J. Electron. Packag.
,
115
(
1
), pp.
131
134
.