<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.2 20190208//EN" "https://jats.nlm.nih.gov/publishing/1.2/JATS-journalpublishing1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" article-type="research-article" dtd-version="1.2" xml:lang="en">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">IJTET</journal-id>
<journal-title-group><journal-title>International Journal of Thermal Engineering and Technology</journal-title></journal-title-group>
<publisher><publisher-name>Thomson &amp; Ryberg Publications</publisher-name></publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.62715/ijtet.2026.v1.i1.pp13-18</article-id>
<article-id pub-id-type="publisher-id">9f5f1b70ea</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Methodology Article</subject></subj-group></article-categories>
<title-group><article-title>Uncertainty Quantification in Transient Plane Source Measurements of Thermal Conductivity</article-title></title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes"><name><surname>Ngo</surname><given-names>Larissa</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1">*</xref></contrib>
<contrib contrib-type="author"><name><surname>Muthukumar</surname><given-names>Andre</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib>
<contrib contrib-type="author"><name><surname>Nielsen</surname><given-names>Meera</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib>
<aff id="aff1"><label>1</label>Department of Engineering, Chao Phraya School of Engineering, Japan</aff>
<aff id="aff2"><label>2</label>Department of Engineering, Balkan Institute of Engineering, Canada</aff>
<aff id="aff3"><label>3</label>Department of Engineering, Vistula Academy of Engineering, Qatar</aff>
</contrib-group>
<author-notes><corresp id="cor1">* Correspondence: <email>larissa.ngo@chao.edu</email></corresp></author-notes>
<pub-date pub-type="epub" iso-8601-date="2026-09-08"><day>08</day><month>09</month><year>2026</year></pub-date>
<volume>1</volume>
<issue>1</issue>
<fpage>13</fpage>
<lpage>18</lpage>
<history>
<date date-type="received" iso-8601-date="2026-05-11"><day>11</day><month>05</month><year>2026</year></date>
<date date-type="accepted" iso-8601-date="2026-07-15"><day>15</day><month>07</month><year>2026</year></date>
</history>
<permissions><copyright-statement>© 2026 The Author(s)</copyright-statement><copyright-year>2026</copyright-year><license license-type="open-access" xlink:href="https://creativecommons.org/licenses/by/4.0/"><license-p>This is an open access article distributed under the Creative Commons Attribution 4.0 International License (CC BY 4.0).</license-p></license></permissions>
<self-uri xlink:href="https://trpub.org/article.php?id=9f5f1b70ea"/>
<abstract><p>The transient plane source technique, widely known as the hot disk method, has become a routine way of measuring the thermal conductivity of solids, pastes and liquids because it is fast and asks little of the specimen. Its results are commonly quoted with an accuracy figure taken from the instrument documentation, which describes the best case and says nothing about a particular specimen, sensor and time window. The measurement is in fact an inverse problem: thermal conductivity and diffusivity are obtained by fitting an idealised conduction model to a temperature transient, and the fitted values inherit uncertainty from the inputs, from the fit and from any departure of the experiment from the model. This article proposes a stepwise procedure for constructing an uncertainty statement for a transient plane source result in line with the Guide to the Expression of Uncertainty in Measurement and its Monte Carlo supplement. The steps cover definition of the measurand and model, identification of input quantities, assignment of distributions, sensitivity analysis, propagation, checks on the validity of the model, and reporting. The procedure is assembled from published theory and standards. It is put forward as a proposal and has not been validated against measurements in this article.</p></abstract>
<kwd-group kwd-group-type="author">
<kwd>Transient plane source</kwd>
<kwd>Hot disk method</kwd>
<kwd>Thermal conductivity</kwd>
<kwd>Measurement uncertainty</kwd>
<kwd>Monte Carlo propagation</kwd>
<kwd>Sensitivity analysis</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="sec1"><title>1. Introduction</title>
<p>Thermal conductivity is among the harder material properties to measure well. Steady-state methods such as the guarded hot plate are capable of low uncertainty but need large, flat specimens and long settling times (Salmon, 2001). Transient methods trade some of that rigour for speed and convenience. The transient plane source technique introduced by Gustafsson (1991) uses a thin metal spiral as both heater and thermometer, clamped between two pieces of the specimen, and returns thermal conductivity and thermal diffusivity from a recording that lasts seconds to minutes. It is now the subject of an international standard (ISO, 2015) and is used across polymers, building materials, ceramics, battery components and insulation.</p>
<p>Ease of use has a cost. Papers that rely on the technique frequently report a conductivity with an uncertainty copied from the manufacturer&#x27;s general specification, or with the standard deviation of a few repeat shots. Neither is an uncertainty in the sense of the Guide to the Expression of Uncertainty in Measurement, the GUM (JCGM, 2008a). Repeatability captures only random scatter under unchanged conditions. A general specification cannot know whether the probing depth stayed inside the specimen, whether contact was good, or how well the specific heat was known.</p>
<p>This article sets out a procedure by which a laboratory can build an uncertainty statement for an individual transient plane source result. It is written for users of commercial instruments who have access to the recorded transient and to ordinary scripting tools. The procedure is proposed on the basis of the published theory of the method and of the metrological guidance documents; it has not been applied to experimental data here, and what would be needed to validate it is stated near the end.</p>
</sec>
<sec id="sec2"><title>2. Background</title>
<p>Two bodies of material underlie the procedure: the conduction model on which the technique rests, and the framework within which measurement uncertainty is evaluated.</p>
<sec id="sec2-2"><title>2.1. The measurement model</title>
<p>The sensor is a double spiral etched from nickel foil and sandwiched between thin insulating films, usually polyimide. A constant current heats it, and its resistance, which rises with temperature, is recorded through a bridge. The technique grew out of the transient hot-strip method (Gustafsson et al., 1979). In the ideal model the spiral is replaced by a set of concentric ring sources in an infinite homogeneous medium, and the solution of the conduction equation (Carslaw and Jaeger, 1959) gives the mean temperature rise of the sensor as ΔT(τ) = P0 D(τ) / (π3/2 r k). Here P0 is the heating power, r the sensor radius, k the thermal conductivity, D(τ) a tabulated dimensionless function and τ = (αt)1/2/r the dimensionless time formed with the thermal diffusivity α (He, 2005).</p>
<p>Because τ contains the unknown diffusivity, the evaluation is iterative. The software searches for the diffusivity, and for a small time correction, that make the recorded temperature rise a straight line against D(τ); the conductivity then follows from the slope. The earliest points are discarded, because they are dominated by the insulating film and the contact between sensor and specimen, and the recording must end before the heat reaches any outer boundary. Bohac et al. (2000) examined how the choice of time window affects the sensitivity of the fit to each parameter and showed that both properties are well determined only when the window spans a suitable range of dimensionless time.</p>
<p>Variants of the basic experiment exist for slabs, thin films and anisotropic materials (Gustavsson et al., 1994; Log and Gustafsson, 1995). Several of them require the volumetric heat capacity as an input, which adds a further measured quantity to the model.</p>
</sec>
<sec id="sec2-3"><title>2.2. The uncertainty framework</title>
<p>The GUM expresses a result as a function of input quantities, assigns each input a standard uncertainty from statistical analysis of repeated observations (Type A) or from other knowledge (Type B), and combines them through first-order sensitivity coefficients (JCGM, 2008a; Taylor and Kuyatt, 1994). The approach assumes that the model is nearly linear over the range of the uncertainties. Supplement 1 replaces the linear combination with numerical propagation of full probability distributions by Monte Carlo sampling and is recommended when linearity or normality is doubtful (JCGM, 2008b). A fitted, iterative model such as that of the transient plane source is a natural candidate for the second route. The engineering literature on experimental uncertainty offers the same ideas in different vocabulary (Moffat, 1988).</p>
</sec>
</sec>
<sec id="sec3"><title>3. The Proposed Procedure</title>
<p>The procedure has seven steps, shown in sequence in Figure 1. Table 1 lists the input quantities that the second and third steps are concerned with, since identifying them is where most of the judgement lies. The steps are then described in turn.</p>
<table-wrap id="tbl1"><label>Table 1</label><caption><p>Input quantities of a transient plane source measurement, how they act on the result and how their uncertainty may be assessed</p></caption><table frame="hsides" rules="groups"><thead><tr><th>Input quantity</th><th>How it enters the result</th><th>Basis for its uncertainty</th></tr></thead><tbody><tr><td>Heating power</td><td>Scales conductivity directly through the slope</td><td>Calibration of current and sensor resistance; lead losses</td></tr><tr><td>Temperature coefficient of resistance</td><td>Converts resistance change to temperature rise</td><td>Sensor calibration certificate; dependence on temperature</td></tr><tr><td>Sensor radius</td><td>Enters the slope and the dimensionless time</td><td>Manufacturing tolerance; effective versus nominal radius</td></tr><tr><td>Time window and time correction</td><td>Fixes which points are fitted and the fitted diffusivity</td><td>Repeat evaluation over alternative windows</td></tr><tr><td>Sensor heat capacity and insulation</td><td>Delays the early response; biases low-conductivity results</td><td>Model with and without correction; reference material</td></tr><tr><td>Contact between sensor and specimen</td><td>Adds an initial offset that the fit must absorb</td><td>Remounting the specimen between repeats</td></tr><tr><td>Volumetric heat capacity</td><td>Required input in slab and anisotropic modes</td><td>Uncertainty of the separate calorimetric measurement</td></tr><tr><td>Specimen temperature drift</td><td>Superimposes a baseline on the transient</td><td>Recorded drift before the heating pulse</td></tr></tbody></table></table-wrap>
<fig id="fig1"><label>Figure 1</label><caption><p>Sequence of the proposed procedure for building an uncertainty statement for a transient plane source result</p></caption><graphic xlink:href="IJTET_V1I1_2026_03_Figure1.png" mimetype="image" mime-subtype="png"/></fig>
<sec id="sec3-2"><title>3.1. Step 1: State the measurand and the model</title>
<p>The statement begins by saying what is being measured: for example, the thermal conductivity of a named material at a stated mean temperature, in a stated direction if the material is anisotropic, at a stated moisture condition if it is porous. The evaluation model is then written down in full, including the sensor geometry, the form of D(τ), the fitted parameters and the criterion by which the time window is chosen. Writing it out exposes which quantities are inputs with their own uncertainty and which are outputs of the fit.</p>
</sec>
<sec id="sec3-3"><title>3.2. Step 2: List the input quantities</title>
<p>Every quantity that enters the model, or that could make the experiment depart from it, is listed. Table 1 gives a starting list. Some inputs are instrumental and the same for every test with a given sensor. Others belong to the specimen and its mounting and change from test to test.</p>
</sec>
<sec id="sec3-4"><title>3.3. Step 3: Assign a distribution to each input</title>
<p>Each input receives a probability distribution. Where repeated observations exist, such as the noise of the temperature signal, a normal distribution with the observed standard deviation is appropriate. Where only limits are known, such as a tolerance on radius, the GUM recommends a rectangular distribution between the limits. Calibration certificates supply a standard uncertainty directly. The reasons for each assignment are recorded, because a reader can only judge the final figure by judging these choices. Correlations should be noted: the temperature coefficient and the power both depend on the same resistance measurement.</p>
</sec>
<sec id="sec3-5"><title>3.4. Step 4: Compute sensitivity coefficients</title>
<p>Before any propagation, each input is perturbed by a small amount and the evaluation is rerun on the recorded transient to find the change in conductivity. The resulting sensitivity coefficients show which inputs matter, and they can differ sharply from one material to another. They also reveal the condition of the fit. Parameter estimation theory shows that two parameters cannot be separated if their sensitivity curves are proportional over the time window (Beck and Arnold, 1977), and a near-proportionality between the sensitivities to conductivity and diffusivity is a warning that the window is poorly chosen. Where inputs interact strongly, a variance-based global analysis is more informative than one-at-a-time perturbation (Saltelli et al., 2008).</p>
</sec>
<sec id="sec3-6"><title>3.5. Step 5: Propagate the distributions</title>
<p>Values of all inputs are drawn from their distributions, the full evaluation, including the iteration and the fit, is repeated for each draw, and the resulting conductivities are collected. Supplement 1 gives guidance on the number of trials and on extracting a coverage interval from the output distribution (JCGM, 2008b). Noise in the transient is included by resampling the residuals of the fit or by adding synthetic noise of the observed variance. The linear GUM combination, using the coefficients from Step 4, should be computed alongside. Agreement between the two gives confidence in both, and disagreement shows that the model is too nonlinear for the simpler method.</p>
</sec>
<sec id="sec3-7"><title>3.6. Step 6: Check the validity of the model</title>
<p>Propagation only covers uncertainty within the model. Whether the model fits the experiment is a separate question, and three checks address it. The probing depth, conventionally taken as 2(αt)1/2 at the end of the window, must be smaller than the distance from the sensor edge to the nearest specimen boundary. The residuals of the fit should scatter randomly about zero; a systematic curve indicates boundary effects, contact resistance or drift. And the whole procedure should be run on a reference material of similar conductivity. If a check fails, the set-up is changed and the measurement repeated; a larger uncertainty is not a substitute.</p>
<p>The third check matters most for insulating materials. Zheng et al. (2020) showed that the standard evaluation can overestimate the conductivity of low-conductivity specimens because the sensor and its insulating film store and conduct a non-negligible share of the heat, which the ideal model ignores. A bias of this kind has to be corrected, or included as a separate component with its own justification.</p>
</sec>
<sec id="sec3-8"><title>3.7. Step 7: Report</title>
<p>The result is reported as a value with an expanded uncertainty and a stated coverage probability, together with the conditions of Step 1, the sensor, the power, the measurement time and the window used. The uncertainty budget, giving each input, its distribution, its sensitivity and its contribution, should accompany the result or be available on request, as laboratories working to recognised quality systems already do. Reporting the fitted diffusivity with its own uncertainty costs nothing extra and lets a reader test the result for consistency against an independent heat capacity.</p>
</sec>
</sec>
<sec id="sec4"><title>4. Fitting the Procedure to Practice</title>
<p>The effort required depends on access to the raw transient. Commercial software exports the recorded temperature rise, and the evaluation described in the standard can be reproduced in a short script, after which Steps 4 and 5 are automatic. A laboratory need not repeat the full analysis for every specimen. A reasonable practice is to carry it out once for each class of material and sensor, keep the resulting budget as a template, and rerun only the checks of Step 6 and the specimen-dependent inputs for routine tests. The published assessment of the transient hot-wire method by Hammerschmidt and Sabuga (2000) is a useful model of what a complete budget for a transient technique looks like.</p>
<p>Some materials need special attention. For anisotropic specimens the heat capacity input usually dominates the budget. For moist or porous materials the measurand of Step 1 must fix the moisture state, since heating can redistribute moisture during the test. For liquids the measurement time has to be short enough to avoid convection, and its onset is a validity check in the sense of Step 6.</p>
</sec>
<sec id="sec5"><title>5. Limitations and What Remains to Be Shown</title>
<p>The procedure is a proposal. It has not been tested here on measured transients, and three questions would need answers before it could be recommended as standard practice. The first is whether the intervals it produces are realistic: on certified reference materials, the reference value should fall within the stated interval about as often as the coverage probability implies. The second is whether different laboratories applying the procedure to the same material arrive at compatible budgets, which only an interlaboratory comparison can show. The third is how sensitive the outcome is to the Type B judgements of Step 3, in particular the treatment of contact and of the sensor&#x27;s own heat capacity.</p>
<p>There are also limits of principle. Monte Carlo propagation cannot reveal an effect that is absent from the model, so the checks of Step 6 carry real weight and depend on the availability of reference materials in the relevant conductivity range. Material inhomogeneity is a property of the specimen and not an uncertainty of the measurement, and it should be reported separately as the spread between specimens.</p>
</sec>
<sec id="sec6"><title>6. Conclusion</title>
<p>A transient plane source result is the output of a fit, and its uncertainty cannot be read from an instrument specification. The seven-step procedure proposed here treats the evaluation as a measurement model in the sense of the GUM, identifies the inputs that carry uncertainty, propagates them numerically through the complete fitting routine, and separates that calculation from checks on whether the ideal model applies. It asks for more work than quoting a nominal accuracy, mainly at the first use for a new class of material. Whether it yields intervals that hold up against reference materials and between laboratories remains to be demonstrated.</p>
</sec>
<sec id="sec7"><title>Conflict of interest</title>
<p>To be declared by the author(s) before publication.</p>
</sec>
<sec id="sec8"><title>Funding</title>
<p>To be declared by the author(s) before publication.</p>
</sec>
</body>
<back>
<ref-list><title>References</title>
<ref id="ref1"><label>1</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Beck, J.V. and Arnold, K.J.</string-name></person-group> (<year>1977</year>). <article-title>Parameter Estimation in Engineering and Science</article-title>. <source>Wiley, New York</source></mixed-citation></ref>
<ref id="ref2"><label>2</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Bohac, V., Gustavsson, M.K., Kubicar, L. and Gustafsson, S.E.</string-name></person-group> (<year>2000</year>). <article-title>Parameter estimations for measurements of thermal transport properties with the hot disk thermal constants analyzer</article-title>. <source>Review of Scientific Instruments</source>, <volume>71(6)</volume>, <fpage>2452-2455</fpage></mixed-citation></ref>
<ref id="ref3"><label>3</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Carslaw, H.S. and Jaeger, J.C.</string-name></person-group> (<year>1959</year>). <article-title>Conduction of Heat in Solids, 2nd edition</article-title>. <source>Oxford University Press, Oxford</source></mixed-citation></ref>
<ref id="ref4"><label>4</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Gustafsson, S.E.</string-name></person-group> (<year>1991</year>). <article-title>Transient plane source techniques for thermal conductivity and thermal diffusivity measurements of solid materials</article-title>. <source>Review of Scientific Instruments</source>, <volume>62(3)</volume>, <fpage>797-804</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1063/1.1142087">https://doi.org/10.1063/1.1142087</ext-link></mixed-citation></ref>
<ref id="ref5"><label>5</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Gustafsson, S.E., Karawacki, E. and Khan, M.N.</string-name></person-group> (<year>1979</year>). <article-title>Transient hot-strip method for simultaneously measuring thermal conductivity and thermal diffusivity of solids and fluids</article-title>. <source>Journal of Physics D: Applied Physics</source>, <volume>12(9)</volume>, <fpage>1411-1421</fpage></mixed-citation></ref>
<ref id="ref6"><label>6</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Gustavsson, M., Karawacki, E. and Gustafsson, S.E.</string-name></person-group> (<year>1994</year>). <article-title>Thermal conductivity, thermal diffusivity, and specific heat of thin samples from transient measurements with hot disk sensors</article-title>. <source>Review of Scientific Instruments</source>, <volume>65(12)</volume>, <fpage>3856-3859</fpage></mixed-citation></ref>
<ref id="ref7"><label>7</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Hammerschmidt, U. and Sabuga, W.</string-name></person-group> (<year>2000</year>). <article-title>Transient hot wire (THW) method: Uncertainty assessment</article-title>. <source>International Journal of Thermophysics</source>, <volume>21(6)</volume>, <fpage>1255-1278</fpage></mixed-citation></ref>
<ref id="ref8"><label>8</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>He, Y.</string-name></person-group> (<year>2005</year>). <article-title>Rapid thermal conductivity measurement with a hot disk sensor: Part 1. Theoretical considerations</article-title>. <source>Thermochimica Acta</source>, <volume>436(1-2)</volume>, <fpage>122-129</fpage></mixed-citation></ref>
<ref id="ref9"><label>9</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>ISO</string-name></person-group> (<year>2015</year>). <article-title>ISO 22007-2:2015, Plastics: Determination of thermal conductivity and thermal diffusivity, Part 2: Transient plane heat source (hot disc) method</article-title>. <source>International Organization for Standardization, Geneva</source></mixed-citation></ref>
<ref id="ref10"><label>10</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>JCGM</string-name></person-group> (<year>2008</year>). <article-title>Evaluation of measurement data: Guide to the expression of uncertainty in measurement (JCGM 100:2008)</article-title>. <source>Joint Committee for Guides in Metrology, Sèvres</source></mixed-citation></ref>
<ref id="ref11"><label>11</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>JCGM</string-name></person-group> (<year>2008</year>). <article-title>Evaluation of measurement data: Supplement 1 to the &quot;Guide to the expression of uncertainty in measurement&quot;, Propagation of distributions using a Monte Carlo method (JCGM 101:2008)</article-title>. <source>Joint Committee for Guides in Metrology, Sèvres</source></mixed-citation></ref>
<ref id="ref12"><label>12</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Log, T. and Gustafsson, S.E.</string-name></person-group> (<year>1995</year>). <article-title>Transient plane source (TPS) technique for measuring thermal transport properties of building materials</article-title>. <source>Fire and Materials</source>, <volume>19(1)</volume>, <fpage>43-49</fpage></mixed-citation></ref>
<ref id="ref13"><label>13</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Moffat, R.J.</string-name></person-group> (<year>1988</year>). <article-title>Describing the uncertainties in experimental results</article-title>. <source>Experimental Thermal and Fluid Science</source>, <volume>1(1)</volume>, <fpage>3-17</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/0894-1777(88)90043-X">https://doi.org/10.1016/0894-1777(88)90043-X</ext-link></mixed-citation></ref>
<ref id="ref14"><label>14</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Salmon, D.</string-name></person-group> (<year>2001</year>). <article-title>Thermal conductivity of insulations using guarded hot plates, including recent developments and sources of reference materials</article-title>. <source>Measurement Science and Technology</source>, <volume>12(12)</volume>, <fpage>R89-R98</fpage></mixed-citation></ref>
<ref id="ref15"><label>15</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Saltelli, A., Ratto, M., Andres, T., Campolongo, F., Cariboni, J., Gatelli, D., Saisana, M. and Tarantola, S.</string-name></person-group> (<year>2008</year>). <article-title>Global Sensitivity Analysis: The Primer</article-title>. <source>Wiley, Chichester</source></mixed-citation></ref>
<ref id="ref16"><label>16</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Taylor, B.N. and Kuyatt, C.E.</string-name></person-group> (<year>1994</year>). <article-title>Guidelines for Evaluating and Expressing the Uncertainty of NIST Measurement Results</article-title>. <source>NIST Technical Note 1297, National Institute of Standards and Technology, Gaithersburg</source></mixed-citation></ref>
<ref id="ref17"><label>17</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Zheng, Q., Kaur, S., Dames, C. and Prasher, R.S.</string-name></person-group> (<year>2020</year>). <article-title>Analysis and improvement of the hot disk transient plane source method for low thermal conductivity materials</article-title>. <source>International Journal of Heat and Mass Transfer</source>, <volume>151</volume>, <fpage>119331</fpage></mixed-citation></ref>
</ref-list>
</back>
</article>