A fully coupled opto-electro-thermal model to investigate silicon solar cells under real operating conditions

. In this work, a fully coupled opto-electro-thermal model for crystalline silicon solar cells is presented. Based on a detailed set of material properties, the developed model allows us to predict and analyse the solar cell behaviour under real operating conditions in a standalone framework. The results show the potential of our model to study the in ﬂ uence of the cell design on its real operating performance, thus giving a new opportunity for silicon solar cell optimisation. Speci ﬁ cally, the doping level is found to impact both the operating temperature and the temperature coef ﬁ cient, showing that two cells with the same power conversion ef ﬁ ciency in standard test conditions can have a very different ef ﬁ ciency under real operating conditions. We also demonstrate the model capability to assess in detail the in ﬂ uence of environmental conditions, such as the solar spectrum, which also impacts the temperature coef ﬁ cient. As the latter is not required by our material-based approach but is a simulation output, this work opens the way to more reliable outdoor prediction. Moreover, the various perspectives and challenges associated with the proposed detailed multiphysics simulation of solar cells are discussed, providing important guidelines for future studies. non-active regions. The temperature-dependent optical and electrical properties introduce an additional coupling between the different physics.


Introduction
Silicon solar cells are designed to efficiently absorb a large part of solar photons but for most of them convert only a limited proportion of sunlight into electricity. Under real operating conditions (ROC) solar cells therefore operate at much higher temperatures [1] than in standard test conditions (STC), i.e., with the reference AM 1.5 spectrum under 1000 W.m À2 and a cell temperature fixed to 25°C. However, this overheating reduces at first the power conversion efficiency as well as the overall photovoltaic (PV) module lifetime [2]. For commercial crystalline silicon (c-Si) devices, the relative efficiency loss is typically between 0.3 and 0.5% for every 1°C increase [3,4]. This combined loss of performance and lifetime significantly reduces the energy yield [5]. Thus, there is a need to develop solar cell models that allow to accurately address these thermal effects in order to both make reliable outdoor prediction and tackle the elevated temperature issues.
As illustrated in Figure 1, the performance of solar cells under ROC involves a complex interplay between their optical, electrical and thermal properties. Current modelling approaches are therefore often based on empirical correlations [6] or on global heat balance analysis [7][8][9]. Although the latter approach is constructed on a physical basis, it generally relies on empirical device parameters such as the temperature coefficients in [8]. Little work has been done to predict the thermal and electrical behaviour directly using material properties as input parameters for simulation. The first milestone for silicon cells was set in 2006 with the work of Vaillon et al. [10]. However, the heuristic model used for the heat sources did not include some important thermoelectric self-heating phenomena such as the Joule and Peltier effect. Using the rigorous heat treatment in semiconductors developed by Lindefelt [11], our previous work [12] has shown the importance of using a complete set of heat sources to consistently describe the electro-thermal coupling. However, the electrical model of reference [12] is based on a depletion zone approximation and the optical part is treated analytically. These two assumptions severely limit the model's ability to capture the behaviour of real silicon devices. Since then, Shang et al. [13] has proposed a model for GaAs (Gallium Arsenide) cells that include complete semiconductor carrier transport equations and the calculation of the optical properties by electromagnetic simulations. Similar studies have been used very recently to study perovskite and CZTS (Copper Zinc Tin Sulfide) structures [14,15]. However, the latter three works do not incorporate detailed material parameters as required for c-Si cells [16], such as the temperature-dependent Klaassen mobility model [17]. Furthermore, they do not achieve a full coupling between the optical and electro-thermal physics which is found to be of great importance.
In this paper, we present a fully coupled opto-electrothermal (OET) model of c-Si solar cells to calculate the temperature and the electrical output in given environmental conditions. The COMSOL Multiphysics software is used to solve the problem using the finite element method. The simulation solves for the complete Maxwell's and carrier transport equations, together with a detailed thermal model which includes all the main self-heating processes occurring in the solar cell (Joule heating, thermalisation, charge carrier recombination, Peltier effect, Thomson cooling) and solar absorption in nonactive layers. Detailed input sets are used for all physical modules, such as a temperature-dependent complex refractive index of c-Si which is found to be crucial for accurate simulation results. Radiative and convective cooling are also implemented. Notably, the former is treated in a free-standing way by computing the optical properties up to mid-infrared wavelengths.
The paper starts with a description of the physical model and the simulation procedure used to solve the coupled OET equations. Then, we investigate the importance of our fully coupled approach by comparing different coupling scenarios on a benchmark c-Si cell. Following that, we unveil the potential of our model to study the influence of the cell design on its real operating performance, thus offering a new route for silicon solar cell optimisation. This is illustrated by showing the influence of the doping level on the electrical and thermal behaviour. We also explore the influence of the solar spectrum on the real operating performance of the solar cell, showing that our model allows for an in-depth investigation of the impact of the environmental conditions. Although a simple solar cell is considered for the sake of illustration, the presented model could be straightforwardly extended to handle more complex designs by means of a higher computational cost. In this respect, the last section aims to provide a clear vision about the prospects and challenges associated with such detailed multiphysics simulation. This work paves the way for a standalone material-based simulation of the opto-electro-thermal behaviour of silicon solar cells under real operating conditions.

Physical models
The various governing equations are implemented using the Wave Optics, Semiconductor, and Heat Transfer in Solids modules of COMSOL Multiphysics. As mentioned in the introduction, these mainly originate from the work of Couderc et al. [12] and Shang et al. [13]. Here, we briefly review the equations and the underlying physics. We also present the various boundary conditions and environmental parameters used in the present paper.

Optical model
The optical absorption is calculated by solving the Helmholtz equation in the frequency domain where E is the electric field, k 0 the free space wavenumber, and n is the complex refractive index of the solar cell materials. The output is the local electromagnetic field, which is then used to compute the optical generation rate G of electrons and holes with E g being the band gap energy of silicon, and g E (E) the optical generation profile for every excitation energy. As for the boundary conditions, Floquet periodic conditions are

Electrical model
The electrical output of the solar cell is obtained by solving the Poisson equation together with the steady-state current continuity equations using the electron (J n ) and hole (J p ) drift-diffusion current densities Here, F denotes the electrostatic potential, q the electron charge, e the permittivity, n (p) the electron (hole) concentration, C the impurity charge concentration, m n (m p ) the electron (hole) mobility, E c (E v ) the conduction (valence) band energy, D n (D p ) the electron (hole) diffusion coefficient, N c (N v ) the effective conduction (valence) band density of states, D n,th (D p,th ) the thermal diffusion coefficient for electrons (holes), and T the temperature. It should be noted that these equations also take into account non-isothermal effects (last term of Eqs. (6) and (7)), which is not usually the case in the simulation of solar cell devices [18,19].
As for charge carrier generation and recombination, we consider optical generation G and non-radiative recombination R arising from Auger and Shockley-Read-Hall (SRH) mechanisms. Radiative recombination are not taken into account because of their minor role in siliconbased devices. For charge carrier collection we assume perfect ohmic contacts on both sides of the device.

Thermal model
The temperature distribution through the device is computed with the heat transfer in solids interface, which solves for the stationary heat equation where k is the thermal conductivity, and H the local heat source. In the non-active layers of the solar cell (metallic contacts, dielectric coatings) and for above band gap wavelengths the electromagnetic power dissipation is used as the heat source. In the active layer (i.e., the crystalline silicon layer), the heat source is the sum of four bulk heat processes (H = H th + H Joule + H NRR + H Tho ) with H th the thermalisation heat source, H Joule the Joule heating, H NRR the non-radiative recombination heating, H Tho the Thomson effect, J the total current density from electrons and holes, k b the Boltzmann constant, p n (p p ) the Peltier coefficients for electrons (holes). These expressions originate from the theoretical work of Lindefelt [11]. The latter has developed a rigorous framework for the treatment of heat in semiconductor devices by deriving the various source terms directly from the Boltzmann transport equation and the fundamental relations of linear irreversible thermodynamics. At the boundary of the active layer we also account for surface recombination [13] and Peltier effect at both the n-type (H Peltier,n ) and the p-type (H Peltier,p ) contacts [13,20] where E Fn is the quasi-Fermi level for electrons and E Fp the quasi-Fermi level for holes. Finally, as boundary conditions for the entire device, we include convective (H con ) and radiative (H rad ) heat transfer at both the front and rear sides where T a is the ambient temperature, s is the Stefan-Boltzmann constant, h c is the convective heat transfer coefficient, and e the total hemispherical emissivity. Heat conduction with the front side contact is not taken into account.

Reference solar cell structure and environmental parameters
The simulated reference solar cell ( Fig. 2) is made of a 10 mm thick p-doped crystalline silicon substrate. Charge carrier separation is realised with an n + front side emitter and a p + back surface field (BSF). A silicon nitride coating and a full-area aluminium layer are also added for enhanced solar absorption. The various material parameters and models are listed in Table 1. For the sake of simplicity, front contacts are evaluated using a numerical shading factor in the optical model and the electro-thermal model is solved assuming 1D charge carrier transport. The solar cell has a simulated efficiency h = 14.68% (open circuit voltage V oc = 671 mV, short-circuit current J sc = 26.03 mA.cm À2 ) under standard test conditions. The low short-circuit current compared to the maximum value of about 44 mA.cm À2 is due to the lack of light trapping through surface texturing and so large reflectance.
As the source of light, we use the standard AM 1.5 solar spectrum, except for Section 3.3, where the effect of the solar spectrum is investigated. Note that we considered normal incidence only, so that for a planar structure the optical response is independent of the polarisation state of the incident plane waves. The convective heat transfer coefficient is taken equal to 5 W. m À2 .K À1 for the front side and 2 W.m À2 .K À1 for the rear side. The total hemispherical emissivity e is estimated according to Kirchhoff's law, by computing the infrared absorptivity through wave optics simulation [34]. Any angular and temperature dependence is neglected in this work for simplicity. The radiative heat transfer on the rear side is also not considered due to the low emissivity of aluminium [35].

Coupling procedure
The equations governing the OET behaviour are coupled together through an iterative procedure: 1) The optical part is calculated by summing the contribution of each wavelength contained in the incident solar spectrum. 2) The electrical and thermal parts are solved selfconsistently using the output of step 1).
These two steps are repeated until the absolute temperature difference between two iterations is smaller than a target value (0.05°C in this work). In the first iteration, a constant temperature distribution (equal to the ambient temperature T a ) is used for the optical calculation. A flow chart is provided in the Appendix, representing the iterative procedure in more detail. The COMSOL direct solver called MUMPS is used for the three physical modules. Figure 3 shows the two main outputs of the fully coupled model, namely the current-voltage characteristic and the solar cell temperature under ROC. As compared to STC, the operating temperature under ROC is not arbitrarily set to 25°C but depends on the thermal equilibrium state reached by the solar cell, which itself depends on the set of environmental condition (solar spectrum, ambient temperature T a , convective heat transfer coefficient h), the solar cell  design, and the bias voltage V [28]. Thus, the operating temperature in ROC is much higher than 25°C. As a consequence, the current-voltage characteristic shows a decrease in V oc and a slight increase in J sc , which are two well-known consequences of increasing cell temperature [36,37]. Moreover, the effect of the bias voltage on the cell temperature can be observed. Indeed, in agreement with previous simulation results [28], the lowest cell temperature is reached near the maximum power point while the temperature is highest at the short-circuit and the opencircuit voltage. This illustrates that the model is able to capture the essential phenomena observed on the device scale from material-based inputs.

Results and discussions
To account for all these, attention must be paid to the coupling between the three physical modules. In the following, we first compare different coupling scenarios to highlight the importance of our fully coupled approach. We then show how this fully coupled model could allow us to study the influence of the solar cell design and the environmental conditions on its performance under ROC. For example, we will see how these two aspects can change the temperature sensitivity coefficient. Thanks to our material-based approach, the latter can be calculated and analysed in a standalone way. Finally, we discuss the numerous perspectives that emerge from this work and the associated challenges. Figure 4 illustrates the impact of different coupling scenarios on the cell temperature and short-circuit current, using the environmental conditions given in Methods section. In fact, the optical, thermal and electrical material properties of solar cells are strongly coupled together in the active layer. In this section, we therefore propose to study how the level of coupling quantitatively impacts the model output (cf., flow chart of the different coupling scenarios in the Appendix). The first scenario of Figure 4 ("w/o coupling") corresponds to the case where there is no coupling, meaning that the operating conditions are those of the standard test conditions and therefore the temperature is equal to 25°C.

The importance of a fully coupled simulation
In the second scenario ("w/ET coupling only"), only the electrical and thermal modules are coupled meaning that there is no feedback of these latter to the optical response (i.e., the complex refractive index at ambient temperature is used) as done in previous multiphysics simulation work [ [13][14][15]. Due to the many heat sources present in the cell, this leads to a much higher temperature than the ambient temperature (+63.2°C). As depicted in Table 2, the most important heat source is the Joule effect, which occurs mainly near the p-n junction. The second heat source is thermalisation, which occurs near the front surface due to the absorption profile of silicon. Thirdly, solar absorption in the non-active materials also plays an important role in the heating of the cell. It mostly arises from the backside metallic contact because many solar photons are not fully absorbed after a single pass through 10 mm of silicon. The Peltier effect at electrical contacts also accounts for a small fraction of the total heating.
The third scenario ("w/full OET coupling") corresponds to the complete coupling achieved in this work. This coupling appears to be necessary to account for the increase in J sc (+1.16 mA.cm À2 ) with temperature mentioned earlier. In the second scenario, a slight decrease in J sc (-0.02 mA.cm À2 ) is even observed due to the lower mobility of the charge carriers at higher temperature. This effect is outbalanced by the temperature dependence of the complex refractive index of c-Si enabled in the third scenario. In fact, the band gap energy decreases as the temperature rises, resulting in more band-to-band solar absorption. The absorption profile of the cell at different temperatures given in the Appendix shows this effect in more detail. This causes the temperature to increase even more (+1.1°C) and slightly changes the distribution of heat sources and sinks. Especially Joule and thermalisation heating are more pronounced than in the second scenario due to a higher electric current and more photogenerated carriers, respectively. The small decrease in non-active . Current density and cell temperature under short-circuit conditions in different coupling scenarios. In the first scenario ("w/o coupling"), the solar cell is operating under STC (its temperature is set to 25°C). In the second scenario ("w/ET coupling only"), the thermal module is coupled to the electrical module, but the optical response is calculated assuming a complex refractive index equal to 25°C. In the third scenario ("w/full OET coupling"), all the three physical modules are coupled together. absorption is due to the band gap shift, which leads to less absorption in the rear aluminium. As we can see, this example highlights the importance of a full OET coupling to accurately study the electrical and thermal output of the solar cell.

Influence of the solar cell design
The developed model allows an in-depth analysis of the influence of the solar cell design on its behaviour under real operating conditions. Here, we focus on the influence of the doping level. For this purpose, we consider a cell with a higher doping level (2.6 Â 10 18 cm À3 for the base, 9 Â 10 19 cm À3 peak doping for the emitter and the BSF) than the reference cell, referred to the "heavily doped cell" in the following. As can be seen in Figure 5a, both solar cells exhibit the same power conversion efficiency (14.68% for the reference cell, 14.67% for the heavily doped cell) under STC. However, under ROC, the heavily doped cell has a better efficiency (12.48%) than the reference cell (12.03%). One reason behind this difference is the lower operating temperature induced by the heavily doped absorber, as depicted in Figure 5d. Indeed, a higher doping level leads to more absorption by free charge carriers in the mid-infrared spectral range [33,25]. Therefore, as stated by Kirchhoff's law, the heavily doped cell exhibits a higher spectral emissivity (Fig. 5b). The result is a hemispherical e emissivity equal to 0.45 for the heavily doped cell, against 0.33 for the reference cell. From Figure 5e, it should be noted that the higher doping level also leads to increased non-radiative recombination (+41 W.m À2 at short-circuit) and more non-active solar absorption (+14 W.m À2 at short-circuit). This mitigates the impact of thermal emissivity on temperature. Nevertheless, the heavily doped cell operates about 2°C lower than the reference cell mostly due to the high sensitivity of the operating temperature to the radiative cooling heat sink. The other reason for the difference under ROC is the better temperature sensitivity of the heavily doped cell. The temperature sensitivity is generally described through the temperature coefficient [38] h where h stc and T stc are respectively the efficiency and the cell temperature in standard test conditions (T stc = 25°C). The value of b can be calculated either from the I-V curves at different temperatures (Fig. 5c), or directly by comparing the solar cell efficiency under real conditions with that under standard conditions. For the two designs considered here, both methods lead to a temperature coefficient equal to 0.34%/°C for the reference cell and 0.30%/°C for the heavily doped cell. Indeed, the doping level is known to increase the open-circuit voltage [39] which in turn decreases (in absolute value) the temperature coefficient [36]. Finally, it is interesting to note how the cell design may change the energy distribution between the different heat sources as seen on Figure 5e. The most notable is the Peltier effect, which changes from a heat source to a sink for the heavily doped cell. In this case, near the contacts, the electron (hole) quasi-Fermi level lies higher in the conduction (valence) band due to the elevated surface doping. From equations (13) and (14), in this configuration, the Peltier effect leads to heat extraction because the absolute difference E c À E Fn at the n-type contact (and E v À E Fp at the p-type contact) is greater than 1.5k b T (cf., band diagram in the Appendix). Figure 5e also confirms that the radiative heat sink of the heavily doped cell is larger than that of the reference cell, due to its higher thermal emissivity.

Influence of environmental conditions
In addition to the influence of the cell design, our model allows us to study in detail the impact of the environmental conditions on the solar cell performance. Here, we propose to investigate the effect of the air mass which is a major variability factor of the terrestrial solar spectra. In order to do so, we compare the model output for two different solar cell spectra (the standard AM 1.5G spectrum and a AM 3G spectrum). The AM 3G spectrum was simulated using the SMARTS 2.9.5 software [40], with the same set of parameters given in reference [12]. Figure 6 shows the current-voltage characteristic of the reference solar cell under both spectra. When the cell temperature is set to 25°C, the J sc and therefore the V oc are lower for the AM 3G spectrum than for the AM 1.5G spectrum. This is no longer the case under ROC, where the V oc under the AM 3G spectrum is higher than under the AM 1.5G spectrum. This is due to the lower cell temperature at V oc (68.7°C under the AM 3G spectrum, 89.3°C under the AM 1.5G spectrum). This lower temperature is caused, on the one hand, by the lower irradiance and, on the other hand, by the better spectral match between the cell and the AM 3G spectrum, which has a lower fraction of high energy photons. Under ROC, this results in a better efficiency under the AM 3G spectrum (+0.71%).
It is also noteworthy that the efficiency under the AM 3G spectrum is different from that estimated using the temperature coefficient calculated in Section 3.2. With the temperature coefficient b = 0.34 %/°C and a temperature at MPP equal to 59.9°C, the efficiency we expect with equation (17) is 12.94%. Nonetheless, the efficiency under ROC is 12.74%. In fact, the temperature sensitivity depends on the incident spectrum [38,41], even though this property is often not taken into account to predict the solar cell efficiency under ROC. The resulting errors might propagate and have a significant impact on, for example, the annual energy yield analysis. Our model avoids this difficulty because it does not require a temperature coefficient as an input. Nevertheless, it allows us to predict the latter and, if needed, could therefore also be used to feed energy yield simulation approaches with a spectrum-dependent b.

Future prospects and challenges
This paper provides a baseline model to study the behaviour of c-Si solar cells under real operating conditions. This type of material-based approach can be used to address a broad scope of scientific topics covering: (1) The influence of the solar cell design (cf., Sect. 3.2); (2) The role of environmental conditions (cf., Sect. 3.3). Put together, this paves the way for solar cell optimisation with thermal criteria as well as for accurate prediction of the OET behaviour in given environmental conditions. We note further that the proposed model could be set up to simulate other devices than c-Si solar cells, only requiring a change in the cell design and the material parameters used as inputs. Moreover, the model offers the opportunity to study thermal mitigation strategies such as radiative cooling [42] in a standalone simulation. By adding fluid dynamics simulations, it would also be possible to predict convective heat exchange from wind speed and installation conditions [43]. Finally, we believe that this model can serve as a core model to evaluate the behaviour of an entire solar module or to perform energy yield estimates [44]. These perspectives are accompanied by major challenges, both in terms of physical modelling and numerical simulation.
As far as multiphysics modelling is concerned, the biggest challenge is to choose the right assumptions for the different physics. This will constrain the complexity and the accuracy of the simulation. One should note that this difficulty is twofold, because it concerns both the modelling of the solar cell and its surrounding environment. For example, in this work, we have treated thermal radiation as a boundary condition but it has been suggested that a participating media approach may be more accurate [1]. As shown in Section 3.1, special attention should also be paid to the degree of coupling. Another concern is about material properties, which are not always readily available. Notably, there is a need for temperature-dependent refractive index values up to the mid-infrared. Even for crystalline silicon only very few models are available [25,33]. Simulation results also need to be compared with experiments. This might be challenging, because environmental conditions can be difficult to manage and also because of the great model complexity. Therefore, our study strongly supports experimental investigation under controlled environmental conditions and a thorough validation of the various building blocks of the model.
The numerical challenges mainly relate to the computational effort associated with the model. The first obstacle originates from the relatively large device thickness compared to the solar spectrum wavelengths. As the thickness increases, so does the number of mesh elements required for the optical simulation. The associated computational cost would therefore be quite significant if our model is used to evaluate, for example, the OET behaviour of a solar module. This is also the case for solar cell architectures that require two-or three-dimensional processing. This issue might be solved by implementing computational cost reduction strategies, such as adopting a multi-scale approach [45,46]. As for using our model to feed simulations for energy yield estimations, this requires quick and easy access to the simulation outputs. This is usually done by using the outputs of the material-based model to create a surrogate model, suitable for the device scale (e.g. a diode model for the electrical part). The possible drawback of this approach is a loss of accuracy and generality. A promising way to overcome this problem might be to employ some recently developed machine learning algorithms [47].

Conclusion
This study has presented a fully coupled opto-electrothermal model for crystalline silicon solar cells operating in real environmental conditions. Based on detailed governing equations and material parameters for the three physical modules, it enables us to predict the electrical output and the temperature of the solar cell by using material-based parameters as the only input. The thermal physics includes all the thermoelectric self-heating processes as well as solar absorption in non-active layers, the latter having been found to be significant for the solar cell considered in this work. A detailed comparison between different coupling scenarios revealed the importance of our fully coupled approach to properly account for the temperature dependence of the optical properties of c-Si.
The results unveiled the potential of our model to study the influence of the solar cell design and the environmental conditions on the solar cell performance under ROC. First, depending on the doping level, we have shown that two solar cells with the same power conversion efficiency under STC can exhibit different efficiencies (about 4% relative difference for the considered environmental conditions) under ROC. The underlying reason was found to be twofold: (1) Different cell temperatures due to different thermal emissivities; (2) Different temperature coefficients attributed to the increase in open-circuit voltage with the doping level. Secondly, our study shed light on the dependence of the temperature coefficient on the solar spectrum. This shows that our model can be used to study in detail the influence of environmental conditions, thus enabling more reliable predictions under ROC than models based on a single temperature coefficient calculated under the AM 1.5G spectrum.
For future work, the prospects and challenges of this study have also been discussed. The challenges relate to both the physical modelling and the numerical simulation. For the former, the challenges mainly concern the finetuning of the model to match the desired level of accuracy. As for the numerical challenges, they are related to the high computational cost when considering thick silicon devices or with two-or three-dimensional structures. Some promising pathways to overcome these limitations have been presented.

Appendix A: Multiphysics coupling scheme
The coupling procedure that has been used to solve for the fully coupled OET behaviour ("w/full OET coupling") is represented as a flow chart in the upper part of Figure 6. The graph shows the two-step procedure for calculating the solar cell operating temperature T and the current density J, for a given bias voltage V. To obtain the current-voltage or temperature-voltage characteristic, this procedure has to be repeated between the short-circuit and open-circuit conditions. In the fully coupled scenario, the electrical and thermal models are coupled together in step 2, and the resulting temperature is then fed back to the optical model, until self-consistent convergence is reach. Figure A.1 also illustrates the other two coupling scenarios that were studied in Section 3.1, namely the "w/o coupling" and the "w/ ET coupling only" scenarios. For the "w/o coupling" scenario, we can see that the thermal module is not used, thus assuming an operating temperature equal to the ambient temperature for both the optical and the electrical parts. For the "w/o coupling" scenario, the thermal module is used, but the operating temperature is still equal to the ambient temperature for the optical part. The values of E c À E Fn + 1.5k b T and E v À E Fp + 1.5k b T involved in the Peltier effect are also depicted, showing that for the heavily doped cell they are negative at the surface, leading to heat extraction. Indeed, in this case, the absolute difference E c À E Fn at the n-type contact (and E v À E Fp at the p-type contact) is greater than 1.5k b T.  Author contribution statement J. Dumoulin developped the COMSOL model and carried out the simulations. The supervision of the project was ensured by M. Amara and E. Drouard. As for the physical model, M. Amara largely contributed to set up the governing equations of the thermal model. The manuscript was written by J. Dumoulin. All the authors discussed the results and actively contributed to manuscript revisions.