Thermal analysis of suspended single droplet evaporation measurements with a coupled lumped parameter model

The measurement data of single droplet evaporation experiments are often biased due to the extra heat input through the fiber suspension and the presence of thermal radiation in hot environments. This encumbers model validation for heat and mass transfer simulations of liquid droplets. In this paper, a thermal analysis of this measurement layout is presented with a coupled lumped parameter model, considering heat conduction through the suspension. The model was validated by experimental data from the literature and good agreements were found. The thermal analysis focused on fiber material and geometry, and thermal radiation properties. Calculations were performed on a broad range of ambient conditions for liquids with different volatility characteristics. Temporal squared droplet diameter- and temperature-profiles, furthermore, droplet stationary evaporation rate were used to characterize vaporization phenomena. The thermal balance of the droplet is dominated by the convective heat rate from the environment in the early stage of evaporation. The effect of heat conduction through the fiber becomes important at the end of the droplet lifetime when the droplet size is decreased. Temperature sensor suspension may seriously bias droplet temperature due to the larger thermal conductivity compared to quartz fiber. Large droplets in high-temperature environments show significant sensitivity to thermal radiation properties, which should be considered in measurements and model validation.


Introduction
Transportation is the largest source of greenhouse gas (GHG) emissions in the European Union.Concerning passenger or light-duty transportation, electric vehicles represent a maturing technology.The question here is the source of electricity production and the long-term environmental impact of battery technology.Heavy-duty transportation, including shipping, aviation, and road freight, has a 46% share in the GHG emissions of transportation [1].Consequently, alternative powertrains are necessary to achieve deep decarbonization.However, heavy-duty transportation has strict constraints on cargo space and payload, therefore, energy sources with high energy density are required to provide economic operation.Consequently, conventional and renewable liquid fuels will remain in the portfolio in the foreseeable future [2], which necessitates heat and mass transfer simulations in engine design.However, calculations must be validated by experimental results.
Comprehensive measurement data on single droplet evaporation is available in the literature.There are measurement layouts, where the droplet is deposited on a hot surface [3,4].However, most experimental setups can be categorized into two groups: stagnant suspended droplet and falling droplet, both aim to obtain the temporal evolution of droplet diameter.Several experimental datasets are published corresponding to the stagnant suspended droplet method for liquids with significantly different fuel properties.Among others, these include n-alkanes, like n-heptane [5][6][7], n-decane [8,9], n-hexadecane [10], vegetable oils [11], jet fuel [12], water and emulsions [13,14], and binary mixtures containing polar and non-polar components [15].The reason for the prevalence of this layout is that the effect of forced convection on droplet evaporation can be eliminated since the droplet is stagnant.Therefore, the influence of thermophysical and transport properties of the fuel can be evaluated independently of operating conditions.Moreover, the optical apparatus detecting the temporal evolution of droplet size can be arranged conveniently and spherical symmetry of the droplet can be provided at a reasonable level.Despite its advantages, the setup has several drawbacks.Thermal radiation may affect the rate of droplet evaporation due to the high-temperature wall of the measurement section and the hot ambient gas.Furthermore, the droplet is fixed to a fiber suspension or temperature sensor.In the latter case, the temporal variation of droplet temperature can be detected.However, both cases seriously affect the thermal balance of the droplet since the suspension has a significantly higher thermal conductivity than the vapor-gas mixture around the droplet surface.Therefore, it acts as an additional heat input, intensifying evaporation.
To overcome the bias in the thermal balance of suspended droplet, the evaporation process of falling droplets is also measured.However, systematic experimental data for a wide range of ambient pressure and temperature is rare in the literature due to the more problematic arrangement of optical apparatus [16].The velocity of the ambient gas is adjusted, resulting in a low relative velocity between the falling droplet and the gas flow.In this manner, the effect of forced convection enhancing evaporation can be decreased.However, thermal radiation can remain important due to the hot surfaces of the measurement chamber.
The present paper focuses on the suspended droplet measurement method since this layout is used more frequently.However, a systematic evaluation of its thermal biases is missing.Figure 1 summarizes the concept of a typical experimental setup based on ref. [11].A pressure regulator and a temperature controller operating a heating rod are used to adjust the ambient pressure and temperature in the measurement chamber.High temperature may lead to droplet burning, therefore, inert gas, which is usually nitrogen, is used to avoid oxidation and focus on heat and mass transfer phenomena.The droplet is placed to the suspension with a thin needle before insertion into the chamber, then it is moved inside with a droplet elevator operated by a stepper motor.In order to avoid preliminary droplet heat-up, the inlet of the chamber is cooled with water in a closed cycle.The optical access to the measurement position is provided by glass windows on the sidewalls.Backlight is applied from one side and the temporal variation of droplet size is recorded from the other side with a high-speed camera, thus binary images are obtained from the droplet by adjusting the threshold.Data is collected with a proper data acquisition (DAQ) device and transferred to a computer for post-processing.The pixel to distance conversion is done by calibration with a known object.After detecting the droplet boundary in the processed images, a high-order polynomial is fitted to it.Rotational symmetry is assumed and the volume of the droplet is determined.The diameter of the droplet is the diameter of a sphere equivalent in volume.Spherical symmetry is a cornerstone in most droplet evaporation models.To preserve this shape, an aircraft moving on a parabolic flight or a free-falling capsule, containing the suspended droplet can be used to carry out experiments.However, they are very expensive and therefore limited in the literature.Another difficulty is that the stable vapor boundary around the droplet obstacles evaporation in microgravity environment [10].
It is rather challenging to quantify the total uncertainty of the experiments.However, it is important to highlight the most important biasing factors, affecting the measured temporal variation of droplet size.The rotational symmetry may be violated due to the presence of the suspension.Curve fitting to the boundary of the droplet in the processed images is also a source of uncertainty due to the resolution of the images.Yang and Wong [17] reported an average ± 15 μm uncertainty, which is 1.5% of the initial diameter and 2.7% of the diameter at the end of the measurement.Despite the water-cooled inlet, an uncontrolled heat-up of the suspended droplet may occur at the beginning of the experiments.Moreover, temperature sensors have larger thermal conductivity than quartz fibers, which can result in the further enhancement of the extra heat input.
The additional heat to the droplet through the fiber suspension enhances vaporization, thus, it can seriously affect model validation.Consequently, this effect should be considered.Generally, two concepts are used for this correction.The first one introduces an empirical correction factor to the stationary evaporation rate [7].The stationary evaporation rate characterizes the temporal surface decrease of the droplet and it is discussed in the next section.Correction factors are determined in terms of fiber suspension diameter for measurements under identical conditions.The reference data corresponds to the smallest fiber diameter and cross-fiber suspension arrangements are usually used for this purpose [7].This method can provide sufficient corrections, however, the empirical factors are confined to the experimental setup and ambient conditions.In the other method, the additional heat is considered as a source term in the thermal balance of the droplet.Several modeling concepts exist for this purpose, like lumped parameter modeling [18], one-dimensional approaches [13,17,19], and multi-dimensional approaches [20][21][22][23].However, systematic analyses on a broad range of ambient conditions, fuel volatility, and suspension material and geometry are scarce in the literature.Therefore, the novelty and the aim of this paper is to fill this gap with a coupled lumped parameter modeling approach to provide a sufficient estimate of the general thermal biasing effects in single droplet evaporation measurements.The focus is on constructing a model with reasonable computational demand applicable for comprehensive parameter analysis.Experimental data is used to evaluate the model and quantify its limitations.Focusing on qualitative, rather than quantitative thermal analysis, general recommendations can be made for further model validation and measurement planning.The applied and tested evaluation methodology of this analysis can be the base for further advanced models.Moreover, droplet evaporation models are also used in non-combustionrelated fields corresponding to different ambient condition range [24], further necessitating the thermal evaluation of droplet vaporization measurements used for validation.The coupled model is presented in the next section.

Suspended droplet evaporation model
The thermal balance of the suspended droplet with the various heat sources and the main concept of the evaporation model are presented in Fig. 2. The fiber is horizontal and Fig. 1 Concept of a typical experimental setup for stagnant suspended droplet evaporation measurement based on ref. [11] the immersed part of it is equal to the instantaneous droplet diameter.Both the thermal balance of the droplet and that of the fiber suspension are considered as lumped parameter models, indicated by the red contours in Fig. 2. No temperature distribution is considered inside the droplet and along the fiber.The droplet surface temperature, T s , is concentrated to the droplet center.The fiber temperature, T f , corresponds to the dry surface of the suspension and only that part is considered in its thermal balance since the dry part is significantly longer than the immersed part.Coupling between the droplet and the fiber thermal balances is performed via the conduction heat rate, Qf .Generally, the temperature of the fiber increases more rapidly due to its lower heat capacity (~10 4 J/K for ~1 cm long SiO 2 fiber [5,25]).This results in a temperature difference between the fiber and the droplet during the evaporation process.
The key equations of the numerical model are presented next.The model was implemented into Matlab R2022a environment.All the liquid-and vapor-phase thermophysical and transport properties of the evaporating droplet are pressure-and temperature-dependent, acquired from the National Institute of Standards and Technology (NIST) database [25].The following equations are solved in each time-step.Two components are considered in the gas phase: fuel vapor and ambient gas.Vapor-ambient gas mixture properties, like specific heat capacity at constant pressure, c pg , thermal conductivity, k g , dynamic viscosity, μ g , and density, ρ g , are calculated for T ref reference temperature and Y v,ref reference vapor mass fraction according to the considerations detailed in ref. [26], assuming ideal mixing when Dalton's law and Amagat's law are valid and the conditions are far from critical.Vapor-phase properties are calculated for T ref , while liquid-phase properties correspond to the droplet temperature.Equations ( 1) and ( 2), solved with the explicit second-order Adams-Bashforth method, describe the thermal balances of the droplet and fiber suspension, respectively: where m d is the droplet mass, c p,L is the liquid-phase specific heat capacity at constant pressure, and T s is considered uniform inside the droplet.Overdot means time derivative.Qconv,d , Qevap , and Qrad,d are the convective heat rate from the environment to the droplet, heat rate of vaporization, and radiative heat rate from the environment to the droplet, respectively.m f is the fiber mass, c f is the specific heat capacity of fiber material, while Qconv,f and Qrad,f are the convective heat rate from the environment to the fiber and radiative heat rate from the environment to the fiber.The (1) where h d is the heat transfer coefficient between the droplet and surrounding gas, d is the droplet diameter and T ∞ is the ambient gas temperature.The heat rate of vaporization is calculated as: where ṁd is the mass flow rate of evaporation and L is the latent heat of vaporization.The latter is determined with the Watson equation [27].The conduction heat rate from the fiber, considering circular cross-section can be written as: where k f , and d f are the fiber thermal conductivity and fiber diameter.The instantaneous distance between the center of the droplet and the free surface of the fiber is represented by d/2.The radiative heat rate from the environment to the droplet is considered as: where (εφ) d describes the emissivity of the droplet and the view factor for the measurement configuration, σ 0 is the Stefan-Boltzmann constant, and T rad is the temperature of the radiation heat source.The latter can represent hightemperature measurement chamber walls, cold walls of the room, or ambient gas temperature.The convective heat rate from the environment to the fiber is calculated as: where h f is the heat transfer coefficient between the fiber and the environment, while l f is the instantaneous dry length of the fiber exposed to heat transfer from the environment.The radiative heat rate between the environment and fiber can be calculated as: where (εφ) f describes the emissivity of fiber and the view factor for the measurement setup.In accordance with the measurements, the initial droplet diameter is considered as an initial condition.Consequently, the occupied volume of suspension is subtracted from the calculated droplet volume, and the droplet mass is determined as follows: where ρ L is the droplet density.The instantaneous dry length of the fiber suspension exposed to the heat transfer from the environment is calculated as: where l 0 and d 0 are the initial dry length of the fiber and initial droplet diameter.The fiber is considered as a horizontal cylinder with d f diameter and l f length.h f is determined from Nusselt number correlations for natural and forced convection.The characteristic length is d f .The Reynolds number, required in case of forced convection is: where u is the relative velocity between the stagnant fiber and the gas flow.For forced convection, the Nusselt number for the fiber is calculated as [28]: while for natural convection, it is calculated as [29]: where the Prandtl number is: and the fiber Rayleigh number for heat transfer is: where the Grashof number for the fiber is: where g and β are the gravitational acceleration and the thermal expansion coefficient of ambient gas.The validity range of Eq. ( 12) is 0.2 ≤ Re f ⋅ Pr, while that of Eq. ( 13) is 0 < Pr < ∞.No information is available for Ra f .The total incoming heat rate of the droplet is defined as: In order to calculate the convective heat rate from the environment to the droplet by Eq. ( 3), h d needs to be determined with Nusselt number correlations [28][29][30].The characteristic (10) length is the instantaneous droplet diameter.For a stagnant non-evaporating sphere, the Nusselt number is: where p ∞ is the ambient pressure, p v,s is the vapor pressure, acquired from the NIST database, corresponding to T s , M a is the ambient gas molecular mass, and M v is the fuel molecular mass.Accounting for the Stefan flow, the mass flow rate of evaporation is: where overdot means time derivative.In order to calculate the ratio of convective mass transfer rate to diffusion rate, the Sherwood number needs to be determined with the following correlations [28].The droplet Sherwood number for natural convection is: where the Schmidt number is: while the droplet Rayleigh number for mass transfer is: The droplet Sherwood number for forced convection is: The mutual diffusion coefficient of fuel vapor and ambient gas is calculated with the method of Fuller et al. [32,33]: where M v,a is the average molar mass of the vapor-ambient gas mixture, while Σ v and Σ a are the sums of atomic and structural diffusion volume increments of vapor and ambient gas.Note that Eq. ( 33) was evaluated with reference data in ref. [34].
The d 2 -profile, shown in Fig. 3, characterizing the temporal evolution of droplet size is acquired by solving the equations in each time step.The stationary evaporation rate, λ st , is determined by fitting a line to the range of linear decrease in the d 2 -profile.In this manner, λ st characterizes the droplet surface decrease over time.The upper and lower limits of this fitting range are often arbitrary.However, the range of 0.15 ≤ (d/d 0 ) 2 ≤ 0.5 is frequently used [10].
where the droplet Rayleigh number for heat transfer is: where the Grashof number for the droplet is: The Nusselt number for a non-evaporating sphere in case of forced convection is: where the Reynolds number for the droplet is: The validity range of Eq. ( 18) is Ra T,d ≤ 10 11 and 0.7 ≤ Pr [28], while Eq. ( 21) is valid for 0.7 ≤ Pr ≤ 400 and 3.5 ≤ Re ≤ 7.6 ⋅ 10 4 [31] (available in Hungarian).The Nusselt number for the droplet, accounting for evaporation, is [30]: The Spalding heat transfer number is: where c p,v is the vapor-phase specific heat capacity at constant pressure, and the Lewis number is: characterizing the relation of thermal boundary layer thickness to the concentration boundary layer thickness, where D v,a is the mutual diffusion coefficient of fuel vapor and ambient gas.The Spalding mass transfer number is: where Y v,∞ is the mass fraction of vapor in the far field, which is considered zero for the single droplet case.Assuming vapor-liquid equilibrium and ideal gas conditions, the mass fraction of fuel vapor on the droplet surface is: that the lower boundary is usually limited by the experimental setup and droplet deformation, therefore, higher values may be applied.The corresponding range of fitting is 0.3 ≤ (d/d 0 ) 2 ≤ 0.5 for the results presented in this paper.If d 2 data are free from bias, the decrease corresponding to the stationary evaporation phase is linear.Therefore, the upper limit can be higher, as well (e.g., (d/d 0 ) 2 = 0.7).However, if fiber suspension is present, the limits of the fitting should be always provided in details since they can affect the numerical value of λ st due to the possible non-linear trend of the d 2 -profile resulting from the extra heat input.Furthermore, when the Knudsen number, relating the mean characteristic free path of the molecules to the droplet radius, is larger than 0.01, the gas phase cannot be regarded as a continuum and kinetic effects should be taken into account.This can be important at the final stage of evaporation or in the case of μm scale droplets generated by modern atomizers.In the case of the currently analyzed droplet size and ambient condition regime, the gas phase can be approximated as a continuum.For further details on kinetic modeling of droplet heating and evaporation, please see ref. [30].

Results and discussion
Section 3.1 presents the validation of the numerical model against experimental data obtained from the literature.Next, the results of the thermal analysis, focusing on various features of the experimental setup are discussed Section 3.2.

Model validation
The presented model was validated against experimental data of Nomura et al. [5], Yang and Wong [17], and Harada et al. [20] since the initial and boundary conditions and several details of their measurements are accurately discussed.These are summarized in Table 1.However, initial fiber temperature values were not detailed directly.Therefore, it was assumed identical to the initial droplet temperature, T s,0 .Next, the features of each experimental setup are summarized and the comparison of measurement data and the results of the model are presented.The published experimental results of Nomura et al. [5] are confined to a measurement chamber with an 80 mm inner diameter and 260 mm height.The ambient gas was nitrogen to prevent droplet burning.Four windows with a 20 mm diameter each provided visual access to the droplet.The ambient gas was heated by an electric furnace inside the chamber.N-heptane droplet was placed on the tip of a silica fiber (SiO 2 , k f = 1.4 W/(m ⋅ K) [35]), which was moved to the desired position by a droplet elevator.The insertion process required 0.16 s, which may led to an uncontrolled slight preheating of the droplet.T ∞ was measured with a thermocouple 4 mm above the test position.A microgravity  environment was used to perform the measurements.The tests were carried out in towers with a height of 5 m and 110 m and parabolic flights were used, as well, to acquire the desired conditions.The whole apparatus with the suspended droplet was covered to eliminate drag force and was placed in the tower, which was evacuated to low pressure.After the setup started to fall, the droplet was introduced.
The experimental setup was fixed to the floor of an aircraft in the case of parabolic flights.After microgravity conditions were achieved, the droplet was introduced.The evaporation process was recorded with a CCD camera.Droplet diameter was determined according to the concept discussed earlier in Fig. 1, assuming spherical symmetry.
The comparison of the experimental data of Nomura and the results of the model is presented in Fig. 4. Note that the time scale is divided by d 0 2 in accordance with the original published data.Solid lines indicate (εφ) d = 1, while dashed lines indicate (εφ) d = 0 to account for the uncertainty of (εφ) d .The former assumes the black body behavior of the droplet and considers the unity view factor, meaning that all the heat radiation from the environment reaches the droplet.The second extreme situation neglects thermal radiation.Obviously, 0 ≤ (εφ) d ≤ 1 condition is valid.Consequently, these boundaries should contain the experimental data.Note that the results showed no significant sensitivity to (εφ) f from the thermal balance of the fiber.Therefore, the effect of possible droplet transparency on the absorbed thermal radiation of fiber suspension was neglected in further calculations.The uncertainty of droplet insertion is indicated by the horizontal error bars.Model results show good agreement with reference data.In the case of (εφ) d = 1, the average relative deviation values are 8%, 7%, 7%, and 9% for 471 K, 555 K, 648 K, and 741 K, respectively.The possible influence of thermal radiation is indicated by higher deviation values at higher temperatures, when this phenomenon is neglected since thermal radiation is proportional to T 4 .The values are 7%, 22%, 31%, and 25% for 471 K, 555 K, 648 K, and 741 K, respectively.
A hot laminar gas flow generated by an electric heater was used in the experiments of Yang and Wong [17].The experimental analysis aimed to investigate the influence of suspension diameter on vaporization.Quartz (silica) fibers with 50 μm, 150 μm, and 300 μm were used and placed in the uniform laminar flow provided by a convergent nozzle.The published uncertainty of temperature measurement was 2 K due to thermal radiation and it was neglected in further calculations.Before the measurements, the droplet was protected from the hot gas flow by a water-cooled shield.At the beginning of the experiment, the shield was withdrawn and the droplet was exposed to the flow.The transient temperature history of the hot flow was also measured and the data was published, making it available for model validation.Flow velocity was measured, as well, with Laser Doppler Anemometry.However, no data was available for model validation.A high-speed camera at a framing rate of 500 fps was used to record the vaporization process.The droplets were considered ellipsoids and the reported uncertainty of diameter values was within ± 15 μm.
Figure 5 presents the comparison of the experimental data of Yang and the results of model calculations.Due to the absence of hot surfaces around the droplet, thermal radiation was neglected, thus, only dashed lines are present.The model provides reasonable accuracy since the calculations slightly underpredict droplet lifetime.However, the transient velocity history of gas flow was not considered in the model, which led to higher h d and enhanced evaporation.For 50 μm, 150 μm, and 300 μm, the average relative deviation between measurement data and calculations are 13%, 21%, and 15%, respectively.
Harada et al. [20] used a Pt-13%Rh/Pt thermocouple (Pt, k f = 71 W/(m ⋅ K)) [36] for the suspension of n-dodecane droplets.The fiber diameter was 50 μm.An alumina protection tube covered a significant part of the sensor.The ambient temperature was adjusted with an electric furnace.The droplet was suspended and placed in a water-cooled probe.Then the whole setup was inserted into the test position and the probe was moved away before the measurement.1000 fps frame rate was used for the high-speed camera to record the images.The droplet diameter was calculated from the area of an equivalent circle.
Figure 6 presents the comparison of experimental data of Harada and model calculations.Blue color corresponds to the d 2 -profile, while the red color indicates droplet temperature.Solid lines represent (εφ) d = 1, while dashed lines indicate that thermal radiation is neglected.Therefore, radiative heat Fig. 4 Comparison of experimental data of Nomura et al. [5] and results of the model.Boundary and initial conditions are presented in Table 1.Solid lines correspond to (εφ) d = 1, while dashed lines correspond to (εφ) d = 0. Uncertainty due to droplet insertion is indicated by the horizontal error bars transfer was considered as a sensitivity parameter again.A reasonable agreement can be observed with measurement data, similar to the multidimensional model of Harada.However, in their model, Harada rightly points out that the heat transfer between the suspension and droplet through the contact surface is a rate-determining factor and the corresponding heat transfer coefficient is a crucial parameter.Unfortunately, the literature has very limited information on accurately determining this coefficient, therefore, it is a potential future work.Note that the stationary evaporation phase is significantly influenced by radiation, indicated by the solid red line.On the other hand, this effect is less obvious from the d 2 -profile.It is also important to highlight that the higher k f of the thermocouple leads to increasing droplet temperature in the stationary evaporation phase.l f was considered in accordance with the protection tube.The average relative deviation values of the d 2 -profile are 22% and 31% for neglecting thermal radiation and (εφ) d = 1, respectively, and 18 K and 10 K, for the temperature-profiles, respectively.Overall, sufficient accuracy is provided by the coupled lumped parameter model.

Thermal analysis
Various features of the experimental setup, affecting the droplet thermal balance, are evaluated next with the presented numerical model.N-alkanes from n-hexane to n-dodecane, except for n-nonane and n-undecane, were analyzed to cover a broad range of fuel volatility.Furthermore, these compounds are often considered in experiments.Highfidelity data for their pressure-and temperature-dependent thermophysical and transport properties are available in the NIST database [25], therefore, the uncertainty resulting from these properties can be minimized.Nitrogen was considered as ambient gas in accordance with the experiments focusing on mass transfer phenomena.No gas flow was considered, therefore, Eqs. ( 13), ( 18) and ( 29) were used for the corresponding calculations, accounting for natural convection.It was assumed that the droplet was inserted into the measurement chamber while it was already suspended, therefore, the initial droplet and fiber temperatures were uniformly 300 K.The fiber material was quartz (SiO 2 ).
The share of different heat sources in the total heat rate, defined by Eq. ( 17), is presented in Fig. 7a.Characteristics of n-hexane and n-dodecane are compared to present the effect of fuel volatility.Blue color corresponds to n-hexane, while red color refers to n-dodecane. Figure 7b indicates the boundary conditions.Thermal radiation is considered with (εφ) d = 0.5.The droplet lifetimes are significantly different, thus, the time scale is nondimensional.t d 2 30% is the time elapsed until d 2 reduces to 30% of d 0 2 .Due to the larger droplet size at the beginning of the evaporation process, Qconv,d dominates.As vapori- zation progresses, the surface area of the droplet reduces due to mass transfer, thus, the share of convective heat rate decreases.However, depending on the experimental layout, droplet size, and ambient temperature, the share of thermal radiation increases, then decreases due to the reduction of droplet size.This is in agreement with the findings of Harada et al. [20].N-hexane is more volatile than n-dodecane, therefore, its droplet size decrease is faster than that of n-dodecane under identical conditions, resulting in a steeper decrease of Qconv,d at the beginning of the process and an earlier maximum of Qrad,d .In the early stage of evaporation, the share of Qf is marginal.[17] and results of the model.Boundary and initial conditions are presented in Table 1.Dashed lines correspond to (εφ) d = 0 Fig. 6 Comparison of experimental data of Harada et al. [20] and results of the model.Boundary and initial conditions are presented in Table 1.Solid lines correspond to (εφ) d = 1, while dashed lines correspond to (εφ) d = 0 1 3 However, as d decreases due to mass transfer, its share increases significantly by the end of the process.Therefore, the stationary evaporation rate is seriously affected.Consequently, d f /d 0 is a crucial parameter in measurements.Qf for n-hexane possesses a slightly higher share than that of n-dodecane.The reason is the following.N-hexane is more volatile and the temperature difference between the fiber and the droplet in the stationary evaporation regime is higher, shown in Fig. 7b.However, this effect is not significant.Due to the notable temperature difference between the fiber and the droplet, a significant temperature gradient occurs along the fiber, which is the most important limitation of the applied modeling approach for the fiber and the main reason for the qualitative rather than quantitative analysis.
As it was highlighted, d f /d 0 is an important parameter of experimental layouts.Therefore, its effect on evaporation characteristics is discussed next.To focus on the effect of heat conduction through the suspension, thermal radiation is neglected this time.l 0 = 1 cm was considered in the analysis and no notable effect of l 0 was recognized since the effect of d f is more dominant.Temporal d 2 -and droplet temperatureprofiles are presented in Fig. 8a   The effect of d f /d 0 on the relative deviation of λ st , extended to a wider ambient condition range for the investigated n-alkanes, is presented in Fig. 9 for 500 K, 700 K, and 900 K gas temperature and for 1 bar and 5 bar ambient pressure.The 5% deviation is indicated with a dashed line again.The relative deviation notably increases with fuel volatility at T ∞ = 500 K since the difference between T f and T s is higher, as detailed in Fig. 7. 500 K is a frequent lower limit for the gas temperature in experiments.Increasing T ∞ diminishes the effect of fuel volatility and the different n-alkanes show practically matching trends.For d f /d 0 < 5%, the relative deviation stays below 5%.900 K is a typical upper limit for the gas temperature in experiments.The effect of p ∞ is the following.The boiling point of the droplet increases with pressure, therefore, the stationary evaporation phase can be characterized by a higher T s .However, T f is not influenced by p ∞ and their difference decreases.This leads to a decrease in Qf , resulting in a decrease in the relative deviation, shown in Fig. 9d-f, compared to Fig. 9a-c.Consequently, increasing pressure decreases the effect of thermal bias through the fiber.If thermal radiation is considered in the thermal balance of the droplet, λ st without fiber increases due to enhanced vaporization.Therefore, the sensitivity of λ st on d f /d 0 decreases.Quartz (SiO 2 ) is the typical suspension material in experiments.However, if droplet temperature is of interest, thermocouples or resistance temperature detectors (RTDs) with extremely small diameters are used the acquire temperature data and act as droplet suspension, as well.Platinum (Pt) is a typical material to solder type R thermocouples and Fig. 9 Effect of fiber diameter-to-initial droplet diameter ratio on the relative deviation of stationary evaporation rate at different ambient conditions (a, b, and c correspond to 1 bar and 500 K, 700 K, and 900 K, respectively, while d, e, and f correspond to 5 bar and 500 K, 700 K, and 900 K, respectively) for various n-alkanes.Thermal radiation is neglected manufacture RTDs.Note that Harada et al. [20] also used a type R thermocouple to acquire droplet temperature, as discussed in Fig. 6.However, k f of Pt is a magnitude higher than that of SiO 2 , seriously affecting Qf .Figure 10a shows the effect of material selection (blue and red colors) on the temporal d 2and T s -profiles (solid and dashed lines, respectively) of an n-dodecane droplet, while the share of Qf in the total heat rate is presented in Fig. 10b.The boundary conditions are indicated in Fig. 10b.Thermal radiation is neglected this time.d f /d 0 = 5% was considered in order to minimize the effect of Qf as much as possible.Even though the thermal conductivity of Pt is significantly higher than that of SiO 2 , their volumetric heat capacities are similar [35,36].Droplet lifetime is significantly shorter for Pt and in the stationary regime, T s is higher by more than 20 K than in the case of quartz fiber.This can seriously affect temperature measurements.Furthermore, no actual stationary state can be observed for the temperature-profile of the Pt case, shown by the increasing red dashed line in Fig. 10a.The thermal balance is dominated by Qf for Pt, shown in Fig. 10b.However, its share remains much lower for SiO 2 during the vaporization process.Consequently, quartz fiber suspension is more favorable and the measured temperature value by the thermocouple suspension can be highly biased.
It is often troublesome to determine (εφ) d accurately for the actual experimental layout.However, the effect of thermal radiation significantly depends on the features of the measurement setup, such as T ∞ and d 0 .Two typical but extreme conditions are discussed next for an n-dodecane droplet at p ∞ = 1 bar.Heat conduction through the fiber is neglected this time to focus on the effect of radiative heat transfer.A larger droplet in higher T ∞ and a smaller droplet in lower T ∞ are considered, indicated with red and blue colors in Fig. 11.The curve parameter is (εφ) d .Figure 11a shows the temporal d 2 -and T s -profiles, while Fig. 11b presents the share of Qrad,d in the total heat Lines styles distinguish the various emissivity and view factor cases, while red and blue colors represent the different initial droplet sizes and ambient temperature conditions rate, where the time scale is non-dimensional again.The share barely exceeds 10% even for the black body assumption for the lower temperature case.However, the sensitivity is significant to (εφ) d for higher temperature and larger droplet size.When black body behavior is assumed, an overshooting tendency can be observed for the larger droplet and higher temperature case, when T s reaches a maximum and then starts to decrease, as shown in Fig. 11a.This behavior was also reported by Sazhin et al. [37] and by Harada et al. [20], where they attributed this maximum to the contribution of thermal radiation.Figure 12 shows the relative deviation of λ st in the function of (εφ) d for the same conditions.The relative deviation remains below 10% for the lower temperature case, however, it exceeds 90% for the high temperature and large droplet case in the extreme (εφ) d = 1 value.Consequently, the uncertainty of (εφ) d in typical experimental layouts can notably affect model validation.

Conclusions
A detailed thermal analysis of single droplet evaporation measurements was performed with a coupled lumped parameter model by revising the thermal balance of the droplet, accounting for the heat conduction through the fiber suspension.The model was validated against experimental data from the literature, showing reasonable agreement.Besides the temporal squared droplet diameter-, and droplet temperature-profiles, the stationary evaporation rate, λ st , was used as an indicator of vaporization characteristics for evaluation.Characteristics of C 6 -C 12 n-alkanes were analyzed in order to cover a broad range of fuel volatility.Pressure-and temperature-dependent thermophysical and transport properties were obtained from the database of the National Institute of Standards and Technology.Based on the results, the following qualitative conclusions can be derived: • The thermal balance of the droplet is dominated by the convective heat rate from the hot gas in the early stage of vaporization.As droplet size decreases, the share of conductive heat rate through the quartz suspension increases, notably enhancing vaporization in the stationary evaporation regime.• The relationship between the quartz fiber diameter-to-initial droplet diameter ratio and the relative deviation of λ st with respect to the case without fiber is non-linear.After a slight increase, a significant rise occurs as the ratio increases.This deviation decreases with increasing ambient pressure due to the increment of droplet boiling temperature.• Using temperature sensors for suspension can lead to serious bias in droplet temperature due to the typically higher thermal conductivity of the sensor, compared to quartz fiber.• Concerning thermal radiation, large (mm-scale) droplets in high-temperature environments show high sensitivity to droplet emissivity and view factor of the experimental setup.Consequently, radiative heat transfer should be carefully considered during measurements and model validation.• The presented coupled lumped parameter model provides reasonable accuracy validated by experimental data from the literature.Consequently, the applied parameter analysis and evaluation method can be the basis of further detailed investigations with more advanced models, where the key factor is the proper definition of the heat transfer coefficient on the contact surface between the suspension and the droplet.permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http:// creat iveco mmons.org/ licen ses/ by/4.0/.

Fig. 2 (
Fig. 2 (Left) Thermal balance of the suspended droplet setup and (right) the main concept of the evaporation model

Fig. 3
Fig. 3 Obtaining the stationary evaporation rate from the temporal d 2 -profile of the droplet

Fig. 5
Fig.5Comparison of experimental data of Yang and Wong[17] and results of the model.Boundary and initial conditions are presented in Table1.Dashed lines correspond to (εφ) d = 0 with no suspension and for different d f /d 0 ratios for an n-dodecane droplet at typical experimental conditions.The time scale is divided by d 0 2 in accordance with several published experimental data, like in Fig. 4. As d f /d 0 is increased, droplet temperature increases too, and droplet lifetime decreases.The relative deviation of λ st with respect to the case without fiber is presented in Fig. 8b in the function of d f /d 0 for the investigated n-alkanes at the same ambient conditions, which are presented in the figures.The 5% deviation value is highlighted with the dashed line.Fuels with different volatility show practically

Fig. 7 aFig. 8
Fig. 7 a Share of the different sources in the total incoming heat rate and b temperature difference between fiber suspension and droplet for n-dodecane (red) and n-hexane (blue)

Fig. 10
Fig. 10 Effect of fiber suspension material on the a d 2 -profile and droplet temperature and b share of conduction heat rate through the suspension in the total incoming heat rate for an n-dodecane droplet.Thermal radiation is neglected

Fig. 12
Fig.12 Effect of emissivity and view factor on the relative deviation of stationary evaporation rate for n-dodecane droplets.Heat conduction through the fiber suspension is neglected