Home A generalised thermal LED-model and its applications
Article
Licensed
Unlicensed Requires Authentication

A generalised thermal LED-model and its applications

  • Ruben Stahlbaum EMAIL logo , Lars Röhe , Martin Kleimeyer and Bert Günther
Published/Copyright: June 30, 2022
Become an author with De Gruyter Brill

Abstract

Within the last 10 years the illuminants for automotive exterior lighting shifted nearly completely to LEDs. Due to being semiconductor devices, LEDs behave differently compared to incandescent lamps and xenon bulbs. The paper derives a generalized thermal and geometric LED model. This gains advantage because the data provided in data sheets is different from manufacturer to manufacturer and even from one manufacturer the data is not standardized. So the data is not prepared to be included easily in any development process. In this context “model” mainly refers to a calculation procedure. The data provided in data sheets often has to be digitized. Outgoing from this digitized data a model, based on a smart data combination and polynomial regression, is built up. This model is described in detail and an application to simulations by means of computational fluid dynamics (CFD) is discussed. In doing so a geometric simplification is suggested. This simplification is done in a manner to keep the thermal characteristic of the original LED. The model may be used in different applications such as simulations and design. It allows predicting the thermal status and light output during a virtual development phase, because it inherently calculates the thermal power and light output. This may lead to a more precise estimation of temperatures in lighting systems as well as a prediction of hot lumens.


Corresponding author: Ruben Stahlbaum, Volkswagen AG, Berliner Ring 2, 38440 Wolfsburg, Germany, E-mail:

  1. Author contributions: All the authors have accepted responsibility for the entire content of this submitted manuscript and approved submission.

  2. Research funding: None declared.

  3. Conflict of interest statement: The authors declare no conflicts of interest regarding this article.

Appendix A – LED model

Calculation of regression coefficients based on T j and I F

I F and T j are given in a definition interval according to data sheet. We will span a two dimensional surface of at least two perpendicular curves according to Figures 7 and 8. Consider the following according to Section 2.3.2:

U F(T j, I F) add two graphs

ϕ rel(T j, I F) multiply two graphs

Δx(T j, I F) add two graphs

Δy(T j, I F) add two graphs.

We perform the regression using a least squares approximation. Let us assume we have a set of points in two variables (x i , y i , z i ) for i = 1,…,n. Let r be the polynomial order in x-direction and s the order in y-direction. By setting

r i = z i v = 0 r w = 0 s α v w x i v y i w

for every i = 1,…,n we want to minimize

Q ( α ) = i = 1 n r i 2 = i = 1 n ( z i v = 0 r w = 0 s α v w x i v y i w ) 2 .

The expression Q can be formulated with the help of the n × (r + 1)(s + 1)-dimensional matrix X and the two vectors

X = ( 1 x 1 y 1 x 1 y 1 x 1 r y 1 s 1 x 2 y 2 x 2 y 2 x 2 r y 2 s 1 x n y n x n y n x n r y n s ) , α = ( α 00 α 10 α 01 α 11 α r s ) T , z = ( z 1 z n ) T

via

Q ( α ) = ( z X α ) T ( z X α ) .

We aim to derive α from the minimization problem Q(α)=!min. Hence, we solve /∂α ij Q(α) = 0 and obtain

Q ( α ) = 2 X T ( z X α ) = ! 0 .

The solution α satisfies

X T X α = X T z .

This enables us to get

α = ( X T X ) 1 ( X T z ) ,

if the ((r + 1)(s + 1))2-dimensional matrix X T X is invertible.

In practice, the LED suppliers give only two measured graphs. This is not enough to derive these full representations in Figures 7 and 8. By using only the two graphs, the regression polynomial oscillates heavily towards the boundaries. To avoid this physically unlikely behaviour, we use superposition to expand the given point set from the supplier linearly throughout the “black area” in the referenced figures.

Physical basis of characteristics

Electric characteristics

The IU-characteristic of an idealized diode was found by William Shockley of Bell Telephone Laboratories in 1949, [5, 10]. It is an exponential dependency of forward current I F on forward voltage U F,

(19) I F = I S ( exp ( U F n U T ) 1 ) .

Herein, n is the emission coefficient and U T the thermal voltage. The details are depending on the certain compound of the semiconductor materials forming the junction and not of importance for the further discussion because the characteristic is usually provided by the LED manufacturers, see upper left of Figure 4 for a red low power LED. These characteristics are usually obtained by measurements done by LED manufacturers.

The thermal voltage U T = k B T j/q relates the kinetic energy of the charges to the potential energy of the total charge q of the P–N-junction. This temperature dependence of the forward voltage is also provided by the manufacturers; see upper right of Figure 4 for the same red LED. The voltage difference, ΔU F = U FU F, typ, is based on a typical voltage, U F, typ, which is explained Section 3.

Luminous characteristics

Prerequesite for efficient light emission is a high probability of radiating recombinations within the P–N-junction, [11]. This is usually the case in semiconductor materials with a direct band gap, e.g. GaAs. The wave length λ of the emitted radiation depends on the energy difference, E g, between valence band and conduction band

(20) λ = h c E g .

It forms a distribution, as can be seen as red curve in Figure 1 for the red LED. The half width and the peak wave length of the spectrum depend, among others, on the temperature. The band gap and thus the energy of emitted photons sink with rising temperature. According to [4] this behaviour was found to be approximately linear for junction temperatures above 290 K according to Varshni formula

(21) E g ( T j ) E g , ref α ( T j T ref ) .

This leads to larger wave lengths and thus to a kind of redshift of the colour of the direct radiation as seen in Figure 5.

Extrapolations

Regression can only be done within the interval of data. If temperature or current outside this interval are used or reached there still shall exist a physically appropriate solution. Thus an extrapolation will be necessary. Let Ω be the domain of definition for the functions of interest according to equations (7) to (10).

Around the rectangular area Ω eight regions appear as can be seen in Figure 9. First we take the vertical and horizontal parts Ω T −, Ω T +, Ω I − and Ω I +.

Extrapolation of voltage

We take equation (7) and calculate the derivative with respect to T j and I F

(22) U F ( T j , I F ) I F = i = 0 n k = 0 m k c U , i k T j i I F k 1

(23) U F ( T j , I F ) T j = i = 0 n k = 0 m i c U , i k T j i 1 I F k .

We use the following physically motivated functions.

  1. In Ω T we perform a linear ansatz mx + a and for every (T, I) ∈ Ω T we find

U F ( T j , min , I F ) T j ( T j T j , min ) + U ( T j , min , I F ) .

  1. In Ω T+ we use an exponential ansatz a ⋅ e b/x and find for every (T j, I F) ∈ Ω T+

U F ( T j , max , I F ) e f e f

with

f = T j U ( T j , max , I F ) U ( T j , max , I F ) T j , max .

  1. In Ω I we do exponential extrapolation a ⋅ e b/x until zero I F = 0 and find for all (T j, I F) ∈ Ω I

U F ( T j , I F , min ) e f 1 e f 2

with

f 1 = I F U F ( T j , I F , min ) U F ( T j , I F , min ) I F , min

and

f 2 = I F U F ( T j , I F , min ) U F ( T j , I F , min ) ( I F , min ) 2 I F .

  1. In Ω I+ we do exponential extrapolation a ⋅ e b/x and find for all (T j, I F) ∈ Ω I+

U F ( T j , I F , max ) e f 1 e f 2

with

f 1 = I F U F ( T j , I F , max ) U F ( T j , I F , max ) I F , max

and

f 2 = I F U F ( T j , I F , max ) U F ( T j , I F , max ) I .

The remaining areas Ω−−/−+/+−/++ comprise two bounding areas which can be used for extrapolation directly.

As an example, we derive a value U at (T, I) ∈ Ω++. From the extrapolated areas above we can observe that U maxmax = U(T j,max, I F,max), U I max = U(T, I F,max), and U T max = U(T j,max, I) from Ω I + resp. Ω T + and the desired value via

U ( T , I ) = U maxmax + [ U I  max U maxmax ] + [ U T max U maxmax ] .

Extrapolation light flux

We take equation (8) and calculate the derivative with respect to T j and I F

(24) Φ rel ( T j , I F ) I F = i = 0 n k = 0 m k c Φ , i k T j i I F k 1

(25) Φ rel ( T j , I F ) T j = i = 0 n k = 0 m i c Φ , i k T j i 1 I F k .

We use the following extrapolations motivated by physics.

  1. In Ω T we do linear extrapolation mx + b and find for every (T j, I F) ∈ Ω T

T j Φ rel ( T j , min , I F ) ( T j T j , min ) + Φ rel ( T j , min , I F ) .

  1. In Ω T+ we apply an exponential extrapolation a ⋅ e b/x and get for every (T j, I F)∈ Ω T+

Φ rel ( T j , max , I F ) e f 1 e f 2

with

f 1 = T j Φ rel ( T j , max , I F ) Φ rel ( T jmax , I F ) T j , max

and

f 2 = T j Φ rel ( T j , max , I F ) Φ rel ( T j , max , I F ) T j .

  1. In Ω I we perform an exponential extrapolation a ⋅ e b/x up to I F = 0, and get for every (T j, I F)∈ Ω I

Φ rel ( T j , I F , min ) e f 1 e f 2

with

f 1 = I F Φ rel ( T j , I F , min ) Φ rel ( T j , I F , min ) I F , min

and

f 2 = I F Φ rel ( T j , I F , min ) Φ rel ( T j , I F , min ) ( I F , min ) 2 I F .

  1. In Ω I+ we perform an extrapolation with the approach A ( x ) + b and we find for every (T, I)∈Ω I+

2 I F Φ rel ( T j , I F , max ) I F , max ( I F I F , max ) + Φ rel ( T j , I F , max ) .

The other domains Ω−−/−+/+−/++ are treated the same as in the previous part on the voltage U.

The extrapolation of the colours like shown in Figure 8 is treated only with linear extrapolation. Otherwise, they are derived in the same way as the values above.

Calculation algorithm

References

[1] K. Baran, M. Leśko, H. Wachta, and A. Różowicz, “Thermal modeling and simulation of high power led module,” in AIP Conference Proceedings, 2019, p. 2078.10.1063/1.5092051Search in Google Scholar

[2] H. M. Choi, L. Wang, S. Kang, J. Lim, and J. Choi, “Precise measurement of junction temperature by thermal analysis of light-emitting diode operated at high environmental temperature,” Microelectron. Eng., vol. 235, p. 111451, 2021, https://doi.org/10.1016/j.mee.2020.111451.Search in Google Scholar

[3] A. M. Colaco, C. P. Kurian, S. G. Kini, S. G. Colaco, and C. Johny, “Thermal characterization of multicolour led luminaire,” Microelectron. Reliab., vol. 78, pp. 379–388, 2017, https://doi.org/10.1016/j.microrel.2017.04.026.Search in Google Scholar

[4] A. Keppens, “Modelling and Evaluation of High-Power Light-Emitting Diodes for General Lighting,” Ph.D. Thesis Light and Lighting Laboratory Gent, 2010.Search in Google Scholar

[5] C. Kittel, Introduction to Solid State Physics, 8th ed. New York, NY, John Wiley and Sons, Inc, 2004.Search in Google Scholar

[6] D. Lee, H. Choi, S. Jeong, et al., “A study on the measurement and prediction of LED junction temperature,” Int. J. Heat Mass Transfer, vol. 127, pp. 1243–1252, 2018, https://doi.org/10.1016/j.ijheatmasstransfer.2018.07.091.Search in Google Scholar

[7] Y. Ohno, “Cie fundamentals for colour measurements,” in Digital Printing Technologies; IS&T’s NIP16, International Conference, CA, No. 16, 2000.10.2352/ISSN.2169-4451.2000.16.1.art00033_2Search in Google Scholar

[8] A. Poppe, “Multi-domain compact modeling of leds: an overview of models and experimental data,” Microelectron. J., vol. 46, pp. 1138–1151, 2015, https://doi.org/10.1016/j.mejo.2015.09.013.Search in Google Scholar

[9] A. Poppe, “Simulation of led based luminaires by using multi-domain compact models of leds and compact thermal models of their thermal environment,” Microelectron. Reliab., vol. 72, pp. 65–74, 2017, https://doi.org/10.1016/j.microrel.2017.03.039.Search in Google Scholar

[10] W. Shockley, “The theory of p–n junctions in semiconductors and p-n junction transistors,” Bell Syst. Tech. J., vol. 28, no. 3, pp. 435–489, 1949, https://doi.org/10.1002/j.1538-7305.1949.tb03645.x.Search in Google Scholar

[11] F. Thuselt, Physik der Halbleiterbauelemente. 2. Auflage, Heidelberg, Springer, 2011.10.1007/978-3-642-20032-8Search in Google Scholar

[12] L. Yang, S. Jang, W. Hwang, and M. Shin, “Thermal analysis of high power gan-based leds with ceramic package,” Thermochim. Acta, vol. 455, pp. 95–99, 2007, https://doi.org/10.1016/j.tca.2006.11.019.Search in Google Scholar

Received: 2022-04-22
Accepted: 2022-05-31
Published Online: 2022-06-30
Published in Print: 2022-12-16

© 2022 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 7.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/aot-2022-0017/pdf
Scroll to top button