Quantifying the efficacy of first aid treatments for burn injuries using mathematical modelling and in vivo porcine experiments

First aid treatment of burns reduces scarring and improves healing. We quantify the efficacy of first aid treatments using a mathematical model to describe data from a series of in vivo porcine experiments. We study burn injuries that are subject to various first aid treatments. The treatments vary in the temperature and duration. Calibrating the mathematical model to the experimental data provides estimates of the thermal diffusivity, the rate at which thermal energy is lost to the blood, and the heat transfer coefficient controlling the loss of thermal energy at the interface of the fat and muscle. A limitation of working with in vivo experiments is the difficulty of measuring variations in temperature across the tissue layers. This limitation motivates us to use a simple, single layer mathematical model. Using the solution of the calibrated mathematical model we visualise the temperature distribution across the thickness of the tissue. With this information we propose a novel measure of the potential for tissue damage. This measure quantifies two important factors: (i) the volume of tissue that rises above the threshold temperature associated with the accumulation of tissue damage; and (ii) the duration of time that the tissue remains above this threshold temperature.

First aid treatment of burns reduces scarring and improves healing. We quantify the efficacy of first aid treatments using a mathematical model to describe data from a series of in vivo porcine experiments. We study burn injuries that are subject to various first aid treatments. The treatments vary in the temperature and duration. Calibrating the mathematical model to the experimental data provides estimates of the thermal diffusivity, the rate at which thermal energy is lost to the blood, and the heat transfer coefficient controlling the loss of thermal energy at the interface of the fat and muscle. A limitation of working with in vivo experiments is the difficulty of measuring variations in temperature across the tissue layers. This limitation motivates us to use a simple, single layer mathematical model. Using the solution of the calibrated mathematical model we visualise the temperature distribution across the thickness of the tissue. With this information we propose a novel measure of the potential for tissue damage. This measure quantifies two important factors: (i) the volume of tissue that rises above the threshold temperature associated with the accumulation of tissue damage; and (ii) the duration of time that the tissue remains above this threshold temperature.
Accidental burn injuries occur most frequently in children, and are a major cause of injury worldwide 1 . Burn injuries worsen after the initial insult 2 , as thermal energy and associated tissue destruction spreads into surrounding tissues. However, if first aid treatment is administered immediately at the scene or prior to qualified medical treatment, burn patients have improved healing outcomes 3 . Optimal first aid treatment conditions have been previously determined using standardised, controlled porcine burn models [4][5][6][7] . Porcine models are used because pig skin is anatomically and physiologically similar to human skin 8,9 , and pig skin responds to therapeutic agents in a similar way to human skin 10 . Recommended first aid treatment is the immediate administration of cool, running water for 20 minutes. This decreases tissue injury, increases the rate of wound healing and reduces scarring 4,6,7,11 . Ethically and morally, there is a limit to the number of different burn and treatment conditions that can be tested using live animals. Therefore, animal experiments cannot be used to examine all potential treatment temperatures and durations.
In this work we use a porcine burn model, illustrated in Fig. 1, to study a series of standardised burns created by exposing the surface of the skin to a temperature of 92 °C for 15 seconds. The propagation of thermal energy through the skin is observed by measuring the temperature under the skin using a subdermal temperature probe. Using this approach we examine the effects of various first aid treatments by applying cooling water of various temperatures (0, 2 and 15 °C), for different durations (10, 20, 30 and 60 minutes). This data set provides a wealth of information about how thermal energy propagates through the skin, and clearly illustrates that different first aid treatments have a measurable effect. However, taken in isolation, even a high-quality reproducible experimental data set such as this does not provide sufficient information to quantify the efficacy of different first aid treatments as this requires more detailed spatial and temporal information about the propagation of thermal energy through the skin.
To address this gap in our knowledge we examine a suite of experimental data describing the effect of various first aid treatment strategies using a mathematical model that describes the spatial and temporal distribution of thermal energy within the skin [12][13][14][15][16] . Owing to the geometry of the experiments, we use a one-dimensional model 1. Calibrate a simplified one-dimensional mathematical model describing the spatial and temporal distribution of thermal energy in living porcine tissues to a suite of experimental data describing various first aid treatment designs. The model is calibrated to provide estimates of: (i) the thermal diffusivity, α; (ii) the loss rate of thermal energy to the blood, β; and (iii) the loss rate of thermal energy at the interface of the subdermal fat and muscle, γ; 2. Repeat the model calibration to obtain separate estimates of α, β and γ for different first aid treatment designs, allowing us to examine whether our estimates of α, β and γ are similar for the different first aid treatment conditions, or whether they might vary significantly between the different first aid treatment conditions; 3. Use the calibrated model to extend the experimental data set by visualising the spatial and temporal distribution of temperature across the tissues; 4. Use the calibrated model to propose a new, biologically motivated and quantifiable measurement of the potential damage associated with the various first aid treatment conditions; and 5. Finally, to use the new measure to quantitatively compare the effectiveness of various first aid treatment conditions.

Results
We consider six different first aid treatments, and we refer to these as six different experimental designs because each first aid treatment involves altering a particular variable, such as the duration of time that the cooling water   , to account for variability between the different experimental replicates. For each experimental design we obtain least-squares estimates of α, β and γ, as summarised in Table 2. Since we have six experimental designs, the calibration procedure gives six values of each parameter, and we report the sample mean and sample standard deviation of each parameter in Table 2. Individual estimates of α, β and γ for each experimental design are reported in the Supplementary Information. Results in Table 2 reveal several insightful trends. For example, the estimate of the thermal diffusivity, α = 0.014 mm 2 /s, is consistent with previously determined estimates in several studies 17,18 . Furthermore, the variability in the estimates of α among the six different experimental designs is very low. For example, the sample standard deviation for α is an order of magnitude smaller than the sample mean, giving a relatively small coefficient of variation, approximately 28%. This degree of variability is extremely small compared to the reported variability in estimates of diffusivities in other biological contexts 19,20 . Another insightful result is that the estimate of β is very close to zero. This means that our model calibration procedure suggests that perfusion of thermal energy into the blood supply is negligible in these experiments. Furthermore, our estimates of β and γ suggest that loss of thermal energy from the system is dominated by heat transfer to the muscle layer at x = H rather than to perfusion into the blood. data. This includes both the initial relatively rapid rise in temperature as a result of the application of the burn, as well as the slower cooling over a longer period of time when the cooling treatment is applied at the surface. Results in Fig. 2c, for T c = 0 °C and t d = 1200 seconds, also shows that the calibrated solution of the model captures the main features of the measured H t ( , )  data. Comparing results in Fig. 2b,c we clearly see the impact of the different first aid treatments as T(H, t) reduces to a lower temperature, and at a faster rate when the 0 °C treatment is applied compared to the 15 °C treatment. Since the least-squares estimates of α, β and γ do not vary too much between the different experimental designs outlined Table 1, from this point forward we will use the average values of these parameters as reported in terms of the sample mean in Table 2.
An interesting feature of our experimental data set is that the temperature of the tissue is 92 °C at the surface of the skin, at x = 0 mm, during the first 15 seconds of the experiment. However, the temperature at the base of the tissue, at x = 4 mm, never rises above approximately 38 °C. This highlights one of the limitations of taking a purely experimental approach. Tissue damage due is thought to occur when the temperature rises above some threshold temperature, and this threshold temperature is reported to be approximately 43 °C 21,22 . One of the aspects of interest in this study is to examine, and visualise, the spatial and temporal distribution of temperature throughout   Table 2. Parameter estimates, showing the sample mean and sample standard deviation for α, β and γ using the six estimates of each parameter associated with the six different experimental designs.
the tissue for all 0 ≤ x ≤ 4 mm. Understanding how the temperature varies across the skin layers would allow us to quantitatively assess how different first aid treatments affect the temperature of the tissues relative to the damage threshold. Since our experimental data provides information at one depth only, at x = 4 mm, it is unclear what the temperature distribution across the entire layer of skin is, and it is also unclear how different treatments affect the distribution of temperatures within the skin layer based on the experimental data alone. Results in Fig. 2d,e show the solution of the calibrated mathematical model, T(x, t), to provide additional insight into the details of the spatial and temporal variation in temperature across the depth of the skin. These results correspond to the experimental designs associated with the results in Fig. 2b,c, respectively. The temperature profiles during the first part of the experiment when the burn is created, 0 < t < 15 seconds, is identical in both cases. However, the temperature profiles differ during the latter part of the experiment, t > 15 seconds, where the different first aid treatments are applied. A qualitative comparison of the T(x, t) profiles in Fig. 2d,e indicates that the main differences are in the upper layer of the tissue, close to the surface at x = 0 mm, whereas the measured differences are at the base of the tissue, at x = 4 mm. Furthermore, the differences between the 15 °C and 0 °C treatments becomes more pronounced at later times. This observation suggests that an improved understanding of the efficacy of the various first aid treatments can be obtained by focusing on the details of the temperature distribution across the entire layer of skin, rather than at just one location.
To provide more detail about how the various first aid treatments affect the spatial and temporal variation of temperature within the tissue we plot the solution of the calibrated mathematical model, T(x, t), using the average values of α, β and γ, given in Table 2, as a series of space-time diagrams 23 in Fig. 3a-c. These plots are constructed for a range of first aid treatment temperatures, T c , including some of the conditions explored here experimentally (T c = 0 and T c = 15 °C), as well as other putative first aid treatment temperatures that extrapolate to conditions that were not examined experimentally (T c = 25 °C). All results in Fig. 3 correspond to the same duration of treatment, t d = 1200 seconds, and results are presented by plotting T(x, t) only where the temperature exceeds the 43 °C 21,22 threshold. Exploring the effects of varying T c in Fig. 3a-c provides very practical insight. For example, results in Fig. 3a,b correspond to T c = 0 and 15 °C, respectively. These conditions could correspond to treating a burn injury with ice cold water and tap water in a cool climate, respectively. In contrast, results in Fig. 3c correspond to T c = 25 °C, which could correspond to treating a burn injury with tap water in a warm, tropical environment. In each plot we truncate the T(x, t) profile, showing only the spatial and temporal distribution for T(x, t) ≥ 43 °C as this is thought to be the critical temperature where tissue damage occurs 21,22 . Therefore, the grey regions in Fig. 3a-c correspond to regions where there ought to be no tissue damage, whereas the coloured regions denote the region where tissue damage might occur.
Each space-time diagram in Fig. 3a-c includes a horizontal line at t = 15 seconds. This line denotes the transition from the first part of the experiment where the burn is created, to the latter part of the experiment where the first aid is applied. Since all experiments are initiated with the same burn condition, the space-time diagrams in Fig. 3a-c are identical for t < 15 seconds. The differences associated with the various first aid treatments is visually discernible in the upper region of the space-time diagrams, where t > 15 seconds. Here we see that the spatial and temporal extent of the region where T(x, t) > 43 °C reduces as T c is reduced. Each space-time diagram  Table 2. Each space-time diagram shows T(x, t), for T x t ( , ) 43 > °C as indicated by the color bars, and grey where < T x t ( , ) 43 °C. (d) shows the potential damage, PD, given by Equation (7). Values of PD are given for T c increasing from 0 °C to 35 °C, in increments of 5 °C. To calculate PD we evaluate the double integral in Equation (6) using the rectangle rule by discretizing the function T * (x, t) with δx = 0.05 mm and δt = 0.0125 seconds.
in Fig. 3a-c also includes a vertical line at x = 4 mm to denote the depth at which experimental measurements are made in our experimental design. Interestingly, as we note previously for the results in Fig. 2d,e, the temperature at x = 4 mm never rises above the threshold of 43 °C 21,22 . However, comparing the spatial and temporal extent of the region where T(x, t) > 43 °C in Fig. 3a,c suggests that different first aid treatment options lead to clear differences in terms of which part of the tissue rises above the threshold of 43 °C 21,22 , and how long the tissue remains above this threshold. Overall, the main effect is that cooler first aid treatments lead to a reduction in the area of the space-time diagram where T(x, t) > 43 °C. While this result is intuitively obvious, the quantitative details of how rapidly this area decreases with decreasing T c is not obvious from intuition alone.
To provide quantitative insight into how the potential for tissue damage varies with the experimental design, we calculate PD, given by Equations (6 and 7), for a suite of first aid treatments, and results are presented in Fig. 3d. These results show how PD* increases with T c , and the details of this dependence reveals some practical information. For example, since PD is an increasing function of T c , there is always a benefit in treating a burn injury with cooler water. However, the slope of the curve in Fig. 3d reveals further insight. The results in Fig. 3d suggests that there is a relatively large reduction in PD if the temperature of the treatment water is reduced from 35 °C to 30 °C. In this case a 5 °C reduction in the temperature of the cooling water reduces PD from 1.57 to 1.35. However, there is a relatively smaller reduction in PD if the temperature of the treatment water is reduced from 5 °C to 0 °C. In this case, a 5 °C reduction in the temperature of the cooling water reduces PD from just 1.02 to 1.00. In summary, the benefit of reducing the temperature of the cooling water is most pronounced for modest to warm first aid treatment temperatures.

Discussion
In this work we present a combined experimental and mathematical modelling study to quantify the efficacy of various first aid treatments using a suite of data from in vivo porcine experiments. In particular, we examine a series of consistent burn injuries that are then subject to various designs of first aid treatments by varying both the temperature and duration of the first aid treatment. Using this extensive experimental data set, we calibrate a mathematical model to reveal estimates of the thermal diffusivity, the rate at which thermal energy is lost to the blood, and the heat transfer coefficient determining the rate at which thermal energy is lost at the lower boundary of the tissue, at the interface of the fat and muscle layers. This calibration exercise provides precise estimates of the in vivo thermal diffusivity, and suggests that the loss of thermal energy in the system is dominated by loss to the muscle rather than loss to the blood. Previous mathematical modelling studies have assumed that thermal energy loss to the blood is negligible 12 , and here we provide novel in vivo evidence to support such assumptions.
The key benefit of working with the mathematical modelling framework is that our extensive experimental data set shows the temperature at just one position at the bottom of the skin tissues, at x = 4 mm, which is at the interface of the fat and muscle. Although burns are created by exposing the surface of the skin, at x = 0 mm, to a temperature of 92 °C, the experimental data at the base of the skin tissues, x = 4 mm, does not record any increase in temperature above approximately 38 °C. This is important because it is thought that damage to tissues occurs when the temperature increases above a threshold temperature of approximately 43 °C 21,22 and while we know that the temperature of the surface of the skin rises well above this threshold, all experimental measurements of temperature at the base of the tissue remain below the threshold. Therefore, using the experimental data alone, it is unclear how much of the tissue rises above the threshold temperature, and it is also unclear how long the temperature of the tissue remains above the threshold temperature. In contrast to the experimental data, the solution of the calibrated mathematical model illustrates the temperature distribution right across the thickness of the tissue and shows precisely how different first aid treatments affect the volume of tissue that rises above 43 °C, and the duration of time that the tissue rises above this threshold. To quantify this we propose a measure of the potential tissue damage, PD, and we show how different first aid treatment designs influence this measure.
The focus of this work is to combine an in vivo experimental data set and a mathematical model to provide quantitative insight into how first aid treatments affect the spatial and temporal distribution of temperature within living skin. Developing a quantitative understanding of the heat transfer process is an essential first step before we can attempt to develop a quantitative understanding of how first aid treatments affect tissue damage. Therefore, this study focuses on the heat transfer processes only. As a consequence, the modelling presented here indicates that first aid treatment is likely to be more beneficial if the treatment water is cooler. One might be tempted to conclude from this study that applying ice at 0 °C is the optimal first aid treatment. However, previous investigations focusing on tissue damage and healing outcomes suggest that ice must never be applied to a burn injury as it is not beneficial for healing 5 , may lead to additional tissue damage 24 , and might also increase the risk of hypothermia. However, our study suggests that if tap water in a tropical climate is warmer than 25 °C, it may be beneficial to apply cooler refrigerated water rather than local tap water for maximum effect. We suggest that further studies that aim to quantify how first aid treatments affect spatial and temporal tissue damage are warranted, and we leave this as a topic for future investigation.
There are several ways that this work could be extended by exploring different mathematical modelling techniques and by invoking different assumptions in the mathematical model. For example, all parameter estimates presented in this work are obtained using a typical constrained nonlinear least-squares optimisation routine that is widely available in MATLAB 25 . This kind of approach is standard throughout the mathematical biology literature, and has been used by us 26 and many others 27,28 . An alternative approach would be to use a computational inference approach, such as approximate Bayesian computation (ABC) 29 . The benefit of working with ABC is that this approach provides a means of obtaining point estimates of the unknown parameters as well as a measurement of the spread about the point estimates 20 . However, the limitation of working with an ABC approach is that it can involve a major increase in computational effort 29 .
Another possible extension of the current work would be to re-analyse the experimental data using a more complicated multilayer reaction-diffusion model [30][31][32] . This approach would involve dealing with distinct values of the thermal capacity, the thermal conductivity and the loss rate to the blood in the three layers of the skin, as well as specifying additional interface conditions describing how the thermal energy is transferred between the epidermis and dermis, and between the dermis and the fat layers 30 . These additional interface conditions may also involve specifying additional parameters, such as contact transfer coefficients between adjacent layers 30 . Taking this kind of more complicated approach would more than triple the number of unknown parameters from three to at least ten, and it is unclear whether the temperature data that we have, at just one location, would be sufficient to provide reliable estimates for this increased number of parameters. However, we do not pursue this kind of extension in the present work and we are confident in our parameter estimates for two reasons. Firstly, the current modelling strategy leads to parameter estimates that are consistent with previous, separate studies 17,18 . Secondly, our parameter estimates do not vary significantly between the six different first aid treatment conditions, and this is precisely what we would expect since α, β and γ act like material properties that are often thought of as constants. Therefore, we leave the extension of re-analysing the data using a more complicated multilayer reaction-diffusion model for future consideration.

Methods
All data generated or analysed during this study are included in this article and the associated Supplementary Information files.

Animal experimental data. Animal studies were ethically approved by the University of Queensland
Animal Ethics Committee (UQAECP&CH/202/06, UQAECMED/RCH/376/08), and all animals were treated humanely and in accordance with the Australian code of practice for the care and use of animals for scientific purposes. The data presented was initially obtained as part of other studies 6,11 . Data from a total of 37 Large White juvenile pigs are included. The animals were approximately 15-20 kg, or 8 weeks old, when data was collected.
Animals were under a general anaesthetic for all procedures and received appropriate analgesia to minimise any suffering. Anaesthesia was inducted with an intramuscular dose of 13 mg/kg ketamine hydrochloride (Ketamine 100 mg/mL, Parnell Laboratories, Alexandria, Australia) and 1 mg/kg xylazine (Xylazil 20 mg/mL, Ilium, Troy Laboratories, Sydney, Australia) and maintained with halothane via a size 4 laryngeal mask airway. Hair was clipped from the back and flanks, and the skin rinsed with clean water before wounding. Buprenorphine hydrochloride at 0.01 mg/kg (Temgesic 0.3 mg/mL, Reckitt Benckiser, West Ryde, Australia) was administered as an analgesic on induction. The subdermal temperature was measured during burn creation and first aid treatment via a temperature probe inserted under the skin. A 14 G 2.1 × 45 mm cannula was obliquely inserted from outside the wound area and advanced through the fat layer under the dermis until the tip was in the centre of the burn area. The needle was removed from the cannula and a type K thermocouple (Radiospares Components Pty Ltd., Smithfield, Australia) was inserted and taped into position. This insertion method was used to minimise the possibility of direct heat transfer onto the probe from the burning device and to minimise damage to the tissue. A digital 54II Fluke thermometer (Fluke Australia Pty Ltd., North Melbourne, Australia) logged temperature measurements every 15 seconds once the burning device was applied and during the course of the first aid treatment. One wound was created in the middle of each flank using a technique described previously 33 . Briefly, a Pyrex laboratory Schott (Mainz, Germany) Duran 500 mL, 75 mm diameter bottle was used (Fig. 1a) which had the glass bottom removed and replaced with plastic wrap secured with tape. A volume of 300 mL sterile water was microwaved until it was approximately 95 °C. The temperature of the water inside the bottle was monitored with a digital thermometer (N19-Q1436 Dick Smith, Australia, range 50 to 100 °C). When the water was at 92 °C, the device was placed on the mid flank for 15 seconds.
Administration of first aid treatment. Treatment commenced as soon as possible after the burn was created. The following different temperature treatments were applied to 25 animals for 20 minutes: • cool running water at 15 °C.
• cool running water at 2 °C; and • granulated ice at 0 °C.
Cool running water was applied at a rate of 1.6 L/min to the wound (Fig. 1b) from a temperature controlled waterbath (Grant S26, Jencons Scientific Limited, Leighton Buzzard, UK) to maintain the treatment at a constant temperature. In addition to varying the temperature of the first aid treatments, we also applied 15 °C cool running water to 20 animals for various durations, including: • 10 minutes; • 20 minutes; • 30 minutes; and, • 60 minutes; During extended treatment periods, the anaesthetised animals were warmed with hot water bottles and core temperatures were monitored.

Measurement of probe depth.
To determine the depth of the temperature probe within the skin, ultrasound imaging (LOGIQe with a L10-22MHz-RS Ultra High Frequency Transducer, GE Healthcare) was used on a different set of animals of similar weight and age with the temperature probe inserted using the same technique. The probes were found to be sitting in the fat layer, 2.6 ± 0.6 mm below the dermis (n = 29 wounds). Additionally, a section of tissue with a probe in situ was dissected out to macroscopically visualise the location of the probe placement (Fig. 1c). A biopsy of normal skin which was fixed in formalin, paraffin embedded, sectioned and stained with Haemotoxylin and Eosin was also obtained as a comparison of the macroscopic and histological appearance of the skin (Fig. 1d). The total thickness of the skin (epidermis and dermis) was determined through histological assessment of biopsies collected on the day of treatment from similar weight animals 17 , and values adjusted for the distortion created by excision and histological processing to relate it to in vivo skin thickness where the skin is stretched 34 . The total skin thickness of the epidermis and dermis was estimated to be 1.4 mm. In summary, the depth of the probe under the surface of the skin is approximately 4.0 mm. Naturally, this depth would vary slightly from animal to animal, and between different locations on the same animal. We will deal with this variability by treating the depth of the skin as a constant value of 4.0 mm, but then averaging our experimental measurements over data from many identically prepared experimental burns.

Mathematical modelling methods.
To model the transfer of heat through skin we use a modified one-dimensional Pennes bioheat equation 14 . This model is a linear heat equation with a source term to account for heat loss to the blood supply (perfusion) 12 . The temporal and spatial distribution of temperature, T(x, t), is governed by, for 0 < x < H and t > 0, subject to the following initial and boundary conditions, , 0 15, , 15 1 5 , where H is the total depth of the tissue including the epidermis, dermis and fat, t d is the duration for which the treatment is applied, α is the thermal diffusivity, β is a heat transfer coefficient governing the rate at which thermal energy is lost to the blood supply, T 0 is a reference temperature which we take to be approximately equal to the temperature of the blood and muscle, T h is the temperature of the heat source used to create the burn, T c is the temperature of the first aid treatment, and γ is a heat transfer coefficient governing the rate at which thermal energy is lost to the adjacent muscle at x = H. When we apply the model to our experimental data we adopt units of millimetres, seconds and degrees Celsius for the length, time and temperature scales, respectively. As with any mathematical modelling study, the application of this model to our experimental data involves several simplifying assumptions. Here we list, and discuss, the key assumptions and provide some justification for invoking these assumptions in this study: Assumption 1: The adoption of a one-dimensional model. The diameter of the burning device, shown in Fig. 1a,b, is approximately 75 mm. In contrast, the depth of the tissue is just 4 mm, as shown in Fig. 1c,d. Therefore, since the lateral length scale of the thermal energy source is much greater than the depth of the tissue, the motion of thermal energy under the centreline of the circular heat source will be in the vertical direction only. This situation can be described, quite accurately, using a simplified one-dimensional model 35 . We acknowledge, however, that if the measurements were made under the skin towards the edge of the thermal energy source, then it would be more appropriate to use a two-or three-dimensional mathematical model.

Assumption 2:
The adoption of a constant initial condition. All results in this work make the simplifying assumption that the initial temperature is given by a constant that is independent of depth, x. An alternative way of specifying the initial condition would be to solve the steady state analogue of Equations (2)-(5), and to use this solution to specify the initial condition that may depend on position, x. However, we avoid taking this approach because measurements reported in our previous work show that the initial spatial variation in temperature across the depth of the tissue is relatively small. For example, the difference between the initial subdermal temperature at x = 4 mm and the initial surface temperature at x = 0 mm is less than 1 °C 17 . Therefore, we choose to neglect this initial variation, and this is reasonable since the variation in temperatures caused by applying the source of thermal energy is much larger, as we have T(0, t) = 92 °C during the first 15 seconds of the experiment. Therefore, we expect that our results are relatively insensitive to our neglect of the small initial spatial variation in T(x, 0). In summary, since our measured initial subdermal temperatures at x = 4 mm is always approximately 36 °C, we set T(x, 0) = 36 °C, for 0 ≤ x ≤ 4 mm.

Assumption 3: The adoption of a homogeneous material.
A key assumption in this work is that we treat α and β as constants. Although skin is known to be layered [36][37][38] , as we show explicitly in Fig. 1d,  meaning that PD is the potential damage function relative to the potential damage incurred by treating the burn with granulated ice, at 0 °C applied for a duration of t d seconds.