Next Article in Journal
A Dataset and a Comparison of Classification Methods for Valve Plate Fault Prediction of Piston Pump
Previous Article in Journal
Research on the Application of Pruning Algorithm Based on Local Linear Embedding Method in Traffic Sign Recognition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameter Extraction for a SPICE-like Delphi4LED Multi-Domain Chip-Level LED Model with an Improved Nelder–Mead Method

Department of Electron Devices, Budapest University of Technology and Economics, Magyar tudósok krt. 2, Bld. Q, H-1117 Budapest, Hungary
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2024, 14(16), 7186; https://doi.org/10.3390/app14167186
Submission received: 10 July 2024 / Revised: 12 August 2024 / Accepted: 13 August 2024 / Published: 15 August 2024
(This article belongs to the Section Electrical, Electronics and Communications Engineering)

Abstract

:
In this paper, a novel method is presented to estimate the parameters of the SPICE-like multi-domain model of light-emitting diode (LED) chips developed and proposed by the Delphi4LED project. The proposed estimation algorithm employs a modified Nelder-Mead method, as the gradient methods and the original version of Nelder-Mead fail to properly handle this problem. By using the new, modified Nelder-Mead method presented in this paper the parameters are estimated faster, compared to the previously used brute-force algorithm-based parameter extraction process, allowing the same precision of the SPICE-like multi-domain LED model. The modification of the parameter extraction procedure also allows speeding up and simplifying the multi-domain LED characterization method proposed earlier by the Delphi4LED project. The speed and robustness of the new model eliminate the need for time-consuming junction temperature control during measurements by employing a novel extraction strategy that seeks the global minimum, rather than relying on the composition of marginal minima.

1. Introduction

In recent decades, LED devices have taken over other light sources in almost every application, such as automotive headlights, interior lighting, streetlamps, as well as sports lighting. The main reason for this rapid change in the use of light sources is that the efficiency/efficacy of LEDs and LED-based luminaires have surpassed those of highly efficient gas discharge lamps, such as high-pressure sodium lamps. Besides the higher efficiency, LEDs offer higher reliability and longer lifetimes. Additionally, they allow for easy electronic control, including dimming, built-in diagnostics, and health-monitoring options. Unlike the classical, phased-out light sources, the operation of LEDs is strongly temperature dependent with both electrical and light-output characteristics influenced by junction temperature, as is typical in semiconductor technology. The thermal-aware design of LED-based lighting solutions is still a challenge that needs novel solutions [1,2,3]. The design of LED packages and luminaires must consider electrical [4,5], optical [6,7], and thermal [8,9] aspects.
Prior to the LED era, luminaire design did not require a close cooperation between the electrical, optical, and thermal development teams. These different design aspects were considered independently, as cross-domain effects were insignificant. However, when using LEDs as light sources, the temperature significantly influences electrical and optical behavior. As a result, thermal, optical, and electrical design teams must work closely together in the product development workflow [10].
For the thermal management of electronic equipment (also known as electronics cooling), computer model-based simulations became daily practice decades ago. This was facilitated by the availability of CFD (computational fluid dynamics)-based numerical thermal simulation tools and the semiconductor vendors’ ability to create so-called boundary condition, independent, compact thermal models of semiconductor device packages [11]. The achievements of the past European research project DELPHI [10,11,12,13,14,15] became widespread both among EDA (Electronic Design Automation) users and semiconductor vendors. A decade after the completion of the DELPHI project, the compact thermal modeling methodology became an industrial standard [16,17].
DELPHI compact models were primarily aimed at the package design of digital ICs for CFD-based thermal simulations of the steady-state or stationary state of digital electronics. These compact thermal models provide the advantage of accurately predicting the junction temperatures of ICs during system-level numerical simulations without requiring the sharing of proprietary information about package design details from semiconductor vendors. A natural extension of the DELPHI compact thermal models was toward thermal transient simulation capabilities, achieved by the European research project PROFIT [18]. This development became important for power semiconductor device packages, where the Cauer-type thermal RC models of the dominantly one-directional heat-flow path proved to be boundary condition independent as well [19]. Such compact thermal models are well suited for the representation of the thermal impedance of the packages of power LEDs as well and can also be easily combined with the electrical models of LEDs, as suggested in [19].
By reducing the detailed mechanical CAD model of an LED luminaire to a multi-port compact model and combining this with the LED chip + package model, complete LED luminaires can also be simulated. Thus, a natural extension of the compact modeling concept for LED packages was the combination of the chip-level multi-domain model of LEDs and the dynamic compact thermal model of the mechanical structure of an LED package in the form of a netlist suitable for application in an electrical circuit simulation program, such as different versions of SPICE. The main target of the Delphi4LED H2020 ECSEL European research project [20] was to develop appropriate multi-domain LED measurement and modeling methodologies, along with an LED product development workflow in which such compact models are used. A major result of the Delphi4LED project was the refinement and reformulation of a prior multi-domain LED model that was implemented and tested under different environments, with the major purpose of supporting the virtual prototyping of LED-based luminaires [21,22]. This multi-domain model has 27 parameters that should be fitted to a set of measured isothermal IVL characteristics of LEDs.
According to the original recommendations of the Delphi4LED project, the isothermal IVL (current–voltage–luminance) characteristics had to cover the typical operational domain of a power LED determined by 33 pairs of I F , T J forward currents and junction temperatures [23]. The huge amount of measurement time needed to obtain the recommended amount of test data and the problems associated with the brute-force fitting techniques used to extract the model parameters have impeded the widespread adoption and use of the proposed multi-domain LED model. Therefore, one of the major targets of the AI-TWILIGHT H2020 ECSEL project [1] (which followed Delphi4LED) was to improve both the measurement processes and the model parameter extraction, both in terms of speed and robustness. Important aspects of the parameter extraction are the number of LED package samples to be measured and a method of finding a single set of model parameters to describe the “LED-type” represented by the measured samples, which was already addressed by the Delphi4LED project [23]. The aim of this paper is to present our recent results that aim to properly address these issues.
As will be shown in Section 2, this model is based on nonlinear formulas with a wide range of possible values in the 27-dimensional parameter space that makes the parameter estimation a real challenge in computing. Our approach, described here, aims to provide a robust, fast, and precise method to extract the model parameters using the least possible sets of test data. Using the extracted parameters and the thermal model described above (the process of thermal model generation is well-established and thus not further detailed here), a digital twin of the LED chip can be developed in a SPICE-like format. This digital twin can be seamlessly integrated with the driver’s model at the electrical port and with the thermal model of the surrounding structure at the thermal port.

2. LED Modeling: The Governing Equations

Many research groups have dealt with the problem of modeling LED chips and luminaires. Two main approaches exist: empirical modeling, which creates models to fit data, and physics-based modeling, which builds models grounded in physical principles. The problem with the former is justifying its generalization to other structures beyond the one used for fitting. The latter cannot be used without the detailed knowledge of the physical structure of the LED chips that is not likely to be publicly shared by the LED manufacturers.
Theoretical models for the optical behavior are established for GaN-based LEDs in [24] and organic LEDs (OLEDs) in [25]. General behavioral models for both LEDs and OLEDs are presented in [26]. Empirical models for optical behavior are presented in [27] and in [28].
In the Delphi4LED multi-domain model, a hybrid approach is used. This is a quasi-black-box model, where the model is based on physical considerations, where the parameters are not calculated from the actual structure, but are identified from measurement results [29]. This hybridization allows for the generalization of the model, as the inter- and extrapolation are based on physical principles, and no information is provided on the actual structure of the LED. Simulation with this model is computationally effective [30] and so it can be used for design and control purposes as well [24].

2.1. The Delphi4LED Chip-Level Multi-Domain LED Model [31]

The Delphi4LED multi-domain model is aimed at modeling the electrical, optical, and thermal properties of LEDs at the chip level. The model is an improved version of the original model proposed earlier [21]. It relies on Shockley’s model of the IV characteristics of pn-junctions. It is a SPICE-like model, meaning that it can be implemented using standard netlist components of an electrical-only, generic version of the SPICE nonlinear circuit simulation program. As mentioned earlier, the electrical operation depends on the temperature, while the generated heat depends on the electrical operating point and the amount of light emitted. The optical behavior depends primarily on the electrical conditions, but it is also influenced by the LED’s operating temperature. The optical losses, eventually, also result in heat generation within an LED package. As the temperature sensitivity of the electrical and optical behavior is high, the model must reflect the link between these domains and the thermal domain.
In the Delphi4LED model, the electrical, optical, and thermal domains are represented as branches (Figure 1). These branches are represented by electrical circuits that can be simulated by any version of the SPICE program. However, according to the modular modeling approach suggested and used in the Delphi4LED project, the model of the thermal impedance of the junction-to-ambient heat-flow path of the LED package (denoted as Z T h in Figure 1) is not part of the chip-level LED model.
The circuit model representation of the branches basically represents their governing equations that we briefly describe below. Each of the elements can be found in Figure 1.
The thermal branch (Figure 1a) is simple, consisting of a heat source terminated by an external model of the Z T h thermal impedance of the package. From the point of view of the steady-state value of the LED chip’s junction temperature, the steady-state value of this thermal impedance should be considered: R T h J A = Z T h t = . The junction temperature can be determined by solving the following energy conservation equation:
P H T J , I F , V F T J T a R T h J A = 0
where P H is the total heat generated in the LED package (also known as the heating power), which depends on the electrical operating point determined by the applied forward current   ( I F ) , the total forward voltage drop ( V F ), and the junction temperature ( T J ). T a is the ambient temperature. For modeling the thermal impedance, different approaches have emerged [9], also within the Delphi4LED project [31,32]. The heating power can be calculated from the applied electrical power ( I F V F , yellow route in Figure 1) subtracting the optical energy leaving the system ( Φ e , orange route in Figure 1).
The electrical branch (Figure 1b) consists of a current source, a series resistance, an ideal diode, and a voltage generator. This represents the LED as a semiconductor diode, seen barely as an electrical component. Since, in practice, LEDs are typically driven by a current source, our multi-domain LED model follows this and is based on a current-driven approach; thus, as indicated in Figure 1, it is driven by the I F current source, forcing the forward current into the electrical model. The next netlist item in the model is a built-in SPICE diode model, characterized by Shockley’s ideal diode equation. In the SPICE diode model used, the calculation of the secondary effects is switched off. The most important effect is represented by the inherent electrical series resistance, R S . In the current driven modeling approach, the ideal diode behavior and the effect of the electrical series resistance are jointly described by
V F = m e U T l n I F I 0 e + 1 + I F R S
where V F is the forward voltage, m e is the diode’s ideality factor, U T = k B T / q is the thermal voltage, I F is the forward current, I 0 e is the saturation current, and R S is the series resistance, the temperature dependence of which is merged into the overall model of the diode’s temperature sensitivity. It is a known fact that, under forced constant forward-current operating conditions, the forward voltage of the diode diminishes with the increasing junction temperature. This temperature-dependent forward voltage change is represented by the Δ V e voltage-controlled voltage generator connected in series with the temperature independent diode and R S models. The controlling voltage of this voltage generator is the voltage of the junction node,   T J , as indicated by the red route in Figure 1. The temperature dependence of the voltage of the Δ V e generator is approximated as follows:
Δ V e = a e I F 2 + b e I F + c e T J 2 T r e f 2 + d e I F 2 + e e I F + f e T J T r e f
where T r e f is the reference temperature for which the parameters of the ideal diode model describing the I F V F relationship and the constant value of the R S series resistance were identified. Note that this temperature dependence of the change in the diode’s overall forward voltage corresponds to the quadratic relationship of the so-called K-factor function proposed by the most recent thermal testing standard, JESD51-51A, published by JEDEC [33]. Symbols a e f e used in Equation (3) are model parameters to be identified.
The model of the optical branch (Figure 1c) is based on the idea that the electrons comprising the total forward current ( I F ) can be divided into two parts, according to the recombination processes to which they are subject to in the active area of the LED’s pn-junction. In this way, we separate the total forward current into a component associated with the radiative recombination processes (direct bandgap transitions of the charge carriers) giving rise to photon emission, and into another component that is associated with non-radiative recombination processes, eventually resulting in heat generation. The current–voltage characteristics of these forward current components can both be described by Shockley’s diode equation [31].
To improve model the optical losses, the ideal light-emitting diode through which the charge carriers pass to perform radiative recombination processes is complemented by a series resistance, denoted as R r in Figure 1. Thus, similarly to the model of the overall electrical behavior of the LED, the optical branch also comprises a temperature independent electrical series resistance ( R r ) and a temperature independent ideal diode model, supplemented by a voltage-controlled voltage generator ( Δ V r ). The controlling voltage of this generator, again, is the voltage representing the junction temperature of the LED. The model equation is similar to Equation (2), but the characteristic here is presented in a voltage-controlled manner, describing the dependence of this radiative diode on the V r forward voltage of the internal pn-junction of the LED (excluding the effect of the R S electrical series resistance):
I r = I 0 r e x p V r Δ V r I r R r m r U T 1
where I r is the radiative current component of the overall forward current of the entire LED, I 0 r is the saturation current, and m r is the ideality factor of the radiative diode. The Δ V r voltage of the voltage-controlled voltage source represents the temperature-dependent changes of the characteristic of this diode. The junction temperature dependence of the Δ V r generator is described with a similar bi-quadratic approximation, as seen earlier for the Δ V e generator:
Δ V r = a r I r 2 + b r I r + c r T J 2 T r e f 2 + d r I r 2 + e r I r + f r T J T r e f
The total emitted optical power of the LED, Φ e (also known as the total emitted radiant flux), is directly proportional to the radiative current:
Φ e = I r Δ V r + m r U T ln I r I 0 r + 1
With the emitted optical power, the value of the P H heating power can be calculated as:
P H = V F I F Φ e
With this last step, the model is closed and the law of the conservation of energy is satisfied.
With additional models, the spectral power distribution or another important light output property, the total emitted luminous flux can be also modeled. The Delphi4LED multi-domain LED model contains a simple, empirical approximation for the total luminous flux, Φ v , with a bi-quadratic formula as follows. The total luminous flux can be determined by:
Φ v = Φ e K I F , T J
where K is the so-called efficacy of the source of the radiation (for the LED) that can be well approximated by a bi-quadratic formula of the forward current and temperature for most LEDs (in which the luminous flux is a relevant light output property), as follows:
K I F , T J = ( a L T J 2 + b L T J + c L ) I F 2 + ( d L T J 2 + e L T J + f L ) I F + ( g L T J 2 + h L T J + l L ) )
Overall, 27 parameters (excluding the thermal impedance that is directly measured and modeled on the package level rather than being part of this model) should be fitted to the measurement data to establish a particular LED model (circuit schematic/netlist with parameters).

2.2. Measurement

The parameters of the Delphi4LED model were fitted to a sufficiently detailed set of so-called isothermal IVL characteristics. These characteristics were obtained using compliant procedures outlined in JEDEC standards JESD51-51A and JESD51-52A [33,34], which involved combined thermal and radiometric/photometric measurements. During these tests, the LEDs under test (DUT) had their junction temperatures set and kept at a pre-specified value while the LEDs were driven by a pre-specified value of the forward current. Thermal impedance, light output characteristics (total emitted radiant flux and spectral power distribution), and forward voltage were measured. Forward current values covered almost three decades, from the mA range up to the 1000 mA range, and junction temperatures varying from room temperature up to 100–120 °C.
The measurement setup utilized SIEMENS SIMCENTER T3Ster, Budapest, Hungary and SIMCENTER, TeraLED equipment (compliant with JEDEC JESD51-51A and JESD51-52A standards [33,34]), completed with a spectroradiometer (Instrumentsystems CAS-140CT, Munich, Germany).
According to the standards, the measurements were conducted in the hot steady-state of the DUT LED, when the junction temperature stabilized after powering. Waiting for the thermal steady state was time-consuming (minutes for each measurement point), but it was necessary to characterize the diode under stable conditions as required by the CIE’s recommendations for the optical testing of LEDs [35].
The measurement control software of the joint T3Ster-TeraLED setup enabled achieving the LEDs’ hot steady state at a pre-specified value of their T J junction temperature and capturing the IVL characteristics at such pre-defined, constant-junction temperature values. The problem with such a temperature control was that, once the forward current of the DUT LED changed, its heating power also changed; thus, the temperature of the cold plate holding the DUT LED had to be adjusted in order to maintain the pre-defined junction temperature. The settling time of this temperature control mechanism took several minutes, representing a significant portion of the overall measurement time to characterize a single LED package. One of the goals of this research is to find a way to establish the LED chip-level model without the requirement of the pre-defined junction temperatures during the measurement. A sample of a measurement result is shown in Figure 2. In subsequent sections, the forward voltage ( V F ) , the radiative flux ( Φ e ) , and the luminous flux ( Φ V ) will be analyzed at each forward current ( I F ) and junction temperature ( T j u n c t i o n ) . This measurement was aligned with the aforementioned standards. Therefore, the junction temperature goal was set according to the reference (first column). Not all reference temperatures can be achieved by the measurement apparatus, though. For instance, a junction temperature of 110 °C cannot be achieved with a 100 mA current due to the LED’s limited heating power. In the original extraction method, such data points were excluded, as discussed in Section 3.1. One other goal of this research is to be able to use these points as well for model extraction purposes.
As detailed in Section 3.2, the new parameter extraction method allows using isothermal characteristics in terms of the ambient temperature control, ( T a ), instead of the characteristics measured at a constant junction temperature, ( T J ). The benefit is that the temperature control loop of the measurement system is much simpler, resulting in shorter overall measurement times, while the level of detail and accuracy of the measurements are maintained, and all the measured data can be used in the extraction phase. Since the chip-level model lacks a thermal component, the junction temperature measurement is still essential.

3. LED Model Parameter Extraction Methods

In the Delphi4LED project, a simple parameter extraction method was implemented in an Excel spreadsheet. This macro was used to find the 27 parameters of the multi-domain model [23]. In this section, the original method is presented first as a baseline for comparing the results obtained with our new, more advanced parameter extraction.
In parallel with the development of the novel extraction algorithm, the extraction strategy was also reformed, enabling us to realize our goals mentioned in the Measurement Section.

3.1. Original Extraction Strategy

The first parameter extraction approach was presented in detail in [22] and briefly described in this section.
As a first step, the reference junction temperature, T r e f , is defined, at which the parameters of the temperature-independent diodes included in the model of Figure 1b,c are valid. At T J = T r e f , the d V F _ e l and d V F _ r a d temperature-controlled voltage generators, which describe the temperature dependence of the two electrical-only diode models in Figure 1, according to Equations (3) and (5), present 0 as the output voltage.
The extraction process starts with the electrical branch. At a high forward current, the series resistance dominates in the overall voltage drop across the LED: I F R S m U T ln ( I F / I 0 + 1 ) ; thus, I F R S   V F . Therefore, the slope of the high current range of the measured I F V F characteristic is a good approximation for the value of R S . Once the serial resistance is estimated, the measured I F V F curve can be stripped of its effect, resulting in an almost ideal diode characteristic. The ideality factor and the saturation current can be estimated by curve fitting Equation (2) to the series of V F values modified this way, yielding the estimated values of m e and I 0 e (values of I F are pre-defined ones, set during the measurements). This initial estimation of the parameters for Equation (2) is quite rough, so the initially obtained value of R S has to be fine-tuned in order to achieve a better overall fit of Equation (2) to the measured dataset. This approach requires an exhaustive search in the parameter spaces of R S , I 0 e , and m e . A good fit at this stage of the parameter extraction procedure is critical since the V r voltage drop on the entire “radiative diode” denoted as branch (c) in Figure 1 is calculated as V r = V F R S I F . The I r current across this “radiative diode” can be calculated by solving Equation (4), with Δ V r = 0 , corresponding to the T J = T r e f assumption.
The optical branch is very similar to the electrical one, but for the I r V r relationship, Equation (4) provides an implicit formula. The relationship between the “experimental values” of I r and the measured Φ e total radiant flux values provided by the optical measurements during the physical LED characterization process is provided by Equation (6).
The six parameters of the Δ V e ( T J ) relationship are determined by quadratic fitting. However, for any T J different from the T r e f reference temperature, a fitting procedure for the of R S , I 0 e , and m e parameters also needs to be performed, requiring substantial computation resources and long execution times, without assuring the best global optimum. The same applies to the parameters of the Δ V r ( T J ) relationship and the R r , I 0 r , and m r parameters of the “radiative diode” denoted as branch (c), Figure 1. The optimal parameters for the entire multi-domain LED model were obtained using a brute-force method [22]. The parameters of the simple model of the luminous flux—Equations (8) and (9)—are determined by quadratic fitting. This process is straightforward, but takes a long time in its 9-dimensional parameter space.

3.2. A Novel Parameter Extraction Strategy and Its Benefits

The main problems with the original extraction method were:
  • T J had to be set during the measurement, requiring a long settling time to meet the target value of T J within an error margin less than 0.5 °C by adjusting the temperature of the DUT (device under test) LED holder of the measurement system. Since the measurement system determines T J indirectly, relying also on an experimentally determined V F T J relationship obtained from a prior calibration procedure, the measured value of T J is loaded by a complex measurement uncertainty budget.
  • Including the IVL characteristics measured at T J = T r e f temperature into the global optimization would assure a better overall fit of the model, as the information contained in these characteristics would also contribute to the modeling of temperature dependence.
  • The sequential extraction of the parameters leads to a set of marginal optima, but not to the global optimum for all the nine parameters for the electrical and optical branches. An algorithm is required to estimate the global optimum in the 9-dimensional parameter space.
Since we aim at using a fully automatic parameter extraction method, we should avoid relying on an a priori distinction between the low- and high-current regions of the characteristics. To achieve such an algorithm, a new optimization method is needed. First trials with usual gradient methods failed, because the gradient in the direction of the saturation current direction is at least 15 orders of magnitude (the typical value of the saturation current is 10 32 10 15 ) higher than the gradients in any other dimension of the global parameter space. To handle this, semi-heuristic approaches are needed. We selected the Nelder–Mead optimization algorithm as our method [36].

3.3. The Nelder–Mead Method

The Nelder–Mead method is a gradient-free, simplex-based direct search method. Let us use the following notations:
  • Let Ν denote the number of parameters to be extracted;
  • let f ( x ) denote the cost function to be minimized during the search for the best fitting parameters in the Ν -dimensional parameter space, and;
  • let x denote a point in the Ν -dimensional parameter space.
During the search for the global minimum of the cost function f ( x ), an initial guess of Ν + 1 point (vectors) was made, each being a candidate for the final set of parameters, where f ( x ) reaches its minimum. This collection of Ν + 1 parameter sets is referred to as a simplex: S = x 0 , x 1 , x 2 x Ν R Ν . The evaluation and modification of S —as the optimization process involves the following steps:
  • Identify the data points for which the f ( x ) cost function provides the largest, the second largest, and the smallest values. Denote these points as follows: ( x max | max i f x i ) , ( x sec | max i max f x i ) , and ( x min | min i max , sec f x i ) .
  • Calculate the coordinates of the centroid of the N best vectors: x c = 1 Ν i m a x x i .
  • Replace x m a x with a new point in the vertex by creating a line between x c and x m a x , and try out different points on this line. The points are called the reflection, the expansion, and the contraction points of x m a x with respect to x c , calculated as follows:
    • Reflection: x r = x c + α x c x m a x . The reflection point replaces the maximum point if f x m i n f x r < f x s e c .
    • Expansion: If f x r < f x m i n , there is a chance to find a better point if we move elsewhere on the line, so x e = x c + γ x m a x x c . If x e is better than x r , the maximum point is replaced with that value; otherwise, we keep the reflection point.
    • Contraction: If f x m a x > f x r > f x s e c , the contraction point will be x c o = x c + β x r x c ; otherwise, if f x r > f x m a x , then x c o = x c + β x m a x x c . If f ( x c o ) < f ( x m a x ) , the condition holds.
    • If no replacement point is found, shuffle all the vectors except the best, x i = x m i n + σ x i x m i n .
The Nelder–Mead method can be terminated under one of these conditions:
  • A lack of improvement after numerous iterations.
  • If the minimum of f ( x ) is close enough to 0.
  • The number of function evaluations or iteration steps reaches its preset maximum number.
The function to be minimized is the following for the electrical part:
E e = i V c a l c i V m e a s i 2 V m e a s i 2
where the summation is applied to every data point of the measured I–V characteristics, with V m e a s i denoting the i -th measured forward voltage and V c a l c i denoting the i -th forward voltage value calculated for the same operating point, where V m e a s i is measured. The division by the measured voltage square is not necessary, but it is needed to avoid low voltage points, especially for the optical part of the model, for which the cost function to be minimized is as follows (with notations similar to the ones applied in Equation (10), with Φ e representing radiant flux values):
E r = i Φ e , c a l c i Φ e , m e a s i 2 Φ e , m e a s i  
Unfortunately, the original Nelder–Mead algorithm fails for the cost functions described in Equations (10) and (11), i.e., it is stuck at a local minimum. With a few modifications, however, we succeeded to overcome this problem.

3.4. Modified Nelder–Mead Method

In a 9-dimensional space, 10 points can easily approach a false minimum (i.e., a local minimum), or the vertex can contract too much into a set of points that are too close to each other, resulting in the premature termination of the optimization process. A possible approach to overcome this problem is to increase the number of data points in the parameter space from Ν + 1 to a higher value. As found in the literature [36,37,38], a potential solution is to combine the Nelder–Mead method with particle swarm optimization (PSO) or with other methods, as described in other papers [39,40]. These approaches certainly provide better results, but the computational effort for this problem is very high, as the number of iterations, which involves solving implicit equations, increases linearly with the number of trial data points used within the Nelder–Mead method.
In our present approach, the Nelder–Mead algorithm is modified to use more points in the simplex. We found that this increased the success rate of the optimization, i.e., the termination of the algorithm with convergence to a true optimum. The initialization (step 1 described in Section 3.3) is modified as follows: The vectors of measured data included in the simplex are chosen randomly from a growing hyper-volume around the initial vector x 0 :
x i k = 1.001 + i 10 N r a n d 1 2 , 1 2 x 0 [ k ]
where r a n d ( ) is a random value from a uniform distribution. The initialization ensures that the points are not too close to each other. The resulting initial simplex is shown in Figure 3.
Increasing the number of vectors in the simplex will not greatly increase the number of cost function evaluations, as only one vector changes over an iteration, except for contraction (see step 3b in Section 3.3), where all vectors change, and the cost function must be evaluated for all of them. Therefore, this modification does not significantly slow down the original Nelder–Mead method. Although the number of cost function evaluations per iteration does not change significantly, the overall optimization’s convergence becomes slower since the central point cannot move as quickly as it could with only Ν + 1 elements in the simplex. The number of iterations is assumed to increase linearly with the number of vectors in the simplex.
The convergence probability of the method is defined as the ratio of successful evaluations to all evaluations. The extraction was performed 100 times with different random simplexes on 10 different LED samples of various kinds for testing purposes. The success rate is defined as the ratio of successful extractions (finding the global minimum, obtained by an exhaustive search) to the total number of attempts. The best heuristic modifications and initial parameter settings of the Nelder–Mead algorithm were identified through experiments on more than 100 sets of measured LED characteristics obtained from the archived data of the Delphi4LED project.
First, we investigated how the Nelder–Mead algorithm behaved with the randomly selected initial Ν + 1 points (for the parameter extraction of the electrical branch Ν = 9 ). The success rate of the algorithm was approximately 50%, even if we followed the original extraction strategy. Increasing the number of points in the simplex resulted in an increased success rate, reaching 96% when the number of points was 4 Ν + 1 = 40 . Further increasing the number of points provided even better success rates; we found that 5 Ν + 1 = 50 yielded a success rate of 98%. The results are summarized in Table 1 and shown in Figure 4, where we can see that the error with the parameters from an unsuccessful extraction is also reduced as the number of vectors increases.
As an illustration, the excess error with respect to the best case of 5 Ν + 1 = 50 points is shown for the parameters of the optical branch of the multi-domain model. For the extraction of the parameters of the optical branch, the number of objective function evaluations for the original Nelder–Mead simplex with Ν + 1 = 10 data points in the parameter space was cca. 14,000 and for the 5-fold increased simplex size of 5 Ν + 1 = 50 , the number of evaluations of the objective function increased to cca. 52,000, slightly less than the expected linear increase.
The success rate of the Nelder–Mead algorithm depends on the actual initial simplex. Therefore, it is essential to investigate how the initial simplex ( x 0 ) influences the success rate. This dependence was investigated by increasing and decreasing the values of the initial parameters by an order of magnitude, except the values of the series resistances and the ideality factors in order to avoid non-physical solutions or no solution for Equation (4). The initial values of the series resistances can be reduced to 0.1 and 0.05, respectively, for electric and optic branches. Ideality factors must be greater than 1.
Decreasing the initial simplex with the above limitations has no adverse effect on the extraction process, other than a slight increase in computational resources, due to 20–30% more iterations being needed to achieve a global optimum. Increasing the parameters, however, greatly increases the probability of the successful execution of the new parameter extraction procedure. These statements are based on our initial numerical experiments using 100+ sets of archived LED measurement data detailed in Table 2.
Reducing the values of the parameter values forming the initial central vector in the parameter space had an adverse effect on the success rate if the parameter extraction process had a 10-fold decrease moderated to a 5-fold reduction. These results are also summarized in Table 2. These experimental results suggest that there are many local minima of the objective function at parameter sets with reduced values. Decreasing the parameter values defining the central vector also narrows the possible range within which additional vectors in the parameter space can be randomly selected. Additionally, overly aggressive upscaling the parameter values results in a decreased success rate. This indicates that using typical parameter sets of the chip-level multi-domain LED model as initial values ensures a high likelihood of success in the extraction process. However, the tested heuristics are not sufficiently robust against higher errors in the measurement data or against unusual LED characteristics. Further steps need to be taken to make our solver robust enough for unsupervised working conditions.
The conclusions drawn from our numerical experiments regarding the heuristic modifications of the original Nelder–Mead method are as follows:
  • Including more random vectors in the simplex increases the success rate of the parameter extraction process, but involves an increased computational cost.
  • Using points from the parameter space with values much less than those forming a valid parameter set for a typical LED chip decreases the success rate of the parameter extraction process.
  • An initial parameter set with parameter values not too far above the values from the expected solution of the parameter extraction process increases the success rate, involving only a slight increase in the computational cost.
  • The above heuristic modifications result in high success rates for the parameter extraction, but do not assure high robustness.
We found that the success rate of the parameter extraction method using these modifications strongly depends on the initial central vector. However, even without increasing the values in the initial set, a high success rate may still be achieved, as presented in Table 2. A thorough investigation of the failed extraction processes revealed that the process often becomes stuck in a local minimum, resulting in much higher cost function values than the global minimum.
To ensure robustness, an additional heuristic approach was added to the initialization step of the parameter extraction process: the number of the initial vectors was increased to 100, covering a significantly extended range, as follows:
x i k = 1.001 + r a n d r a n d i 100 , i 1 100 , r a n d i 1 100 , i 100 x 0 k ,
where i is the index between 1 and includes 100. The cost function for these points is evaluated, and the 30 points with the smallest cost function values are kept, for which we used the Nelder–Mead algorithm modified with the above-described heuristics. With this addition, the success rate of the parameter extraction process reached 100% in the case of 5200 extractions performed with measurement data from 52 different LED samples of various types, including color LEDs and phosphor-converted white LEDs. The choice of the central point has a much smaller effect than before. The success rate for increased parameter values remains 2600/2600 (50 extraction per LED sample); with 10-fold reduced values, the success rate was 2595/2600, nearly 100% This suggests that the Nelder–Mead algorithm, modified with these heuristics, is now robust enough and likely to work well for other LED types beyond those for which we had archived measurement data from the Delphi4LED project.

4. Results Obtained with the Modified Nelder–Mead Method

We investigated the behavior of our modified Nelder–Mead parameter extraction method with two kinds of benchmarks.
The first benchmark compares the overall error of the parameter fitting performed with the new method with the overall error of the prior brute-force approach used in the Delphi4LED project. This comparison was chosen because the brute-force-based method—though very slow and computationally intensive—provides very good fits (less than a 1% error for V F and less than a 5% error for Φ e ) for all the LED samples characterized in the Delphi4LED project [41]. Therefore, this comparison offers a valuable insight into the fit errors of the new parameter extraction procedure. However, caution is needed for this comparison, as there are several significant differences between the two methods.
  • The old method requires predefined junction temperatures, so the measurement points that do not meet this requirement are not taken into account.
  • Slight differences in junction temperatures (±0.5 °C) are neglected in the old method.
  • The old and the new parameter extraction methods also differ in the objective functions that they minimize during the fitting process. The old method uses the sum of square errors as the objective function to be minimized, while in the new method, the sum of the relative square error is used. Therefore, for the used figure of merit for the comparison of the two methods, the maximal relative error is much higher at low forward currents when applied to evaluate the fit error of the old method.
The new model uses the sum of the relative square error, because many end-users of the SPICE-like multi-domain LED model may require accurate simulations across a wide range of forward currents, spanning two decades (for our sample LEDs: from 20 mA up to 1000 mA). As a result, the LED model should be equipped with a parameter set that ensures accurate results, even below the greatest decade in the forward current ( I F < 100   m A ) . To address this need, we compared the performance of the two parameter extraction methods, both with and without including low-forward-current operating points. In the comparisons for the high current domain, we completely discarded the low-current operating points of the LEDs during the fitting procedures.
The parameter extraction procedures were applied to 52 high-power LED samples measured in the Delphi4LED project: 11 phosphor-converted white LED samples of the same type, 5 red LED samples of the same type, 5 amber LED samples of the same type, and 31 blue LED samples with 11 samples from the same LED type and 20 other samples of other types. The maximum relative error values obtained with the old and new parameter extraction procedures are shown in Figure 5 for the electrical branch, Figure 6 for the optical branch, and Figure 7 for luminous flux, applied to the high-forward-current operating points only (IF > 100 mA) and for the full range of operating points (involving the low-forward-current ones as well).
In the high-forward-current region (above 100 mA), where the “linearization effect” of the series resistance starts to appear in the LED characteristics, the strong dependence of the accuracy of the parameters of the pn-junction-related parts of the multi-domain LED model is increasingly suppressed as the forward current increases.
Consequently, the simulation results obtained using the parameter sets identified by the previous method and those identified by the new, modified Nelder–Mead-based parameter extraction process agree very well in this current range.
Considering the entire set of LED operating points, including forward currents below 100 mA, the simulation results obtained from the two different kinds of parameter sets show greater differences. Besides the stronger dependence of the accuracy of the model calculations on the accuracy of the parameter sets at lower forward currents, it is worth mentioning that the previous and the new parameter extraction processes use different objective functions.
The new method uses the relative square difference as the cost function, ensuring acceptable accuracy across the entire range of operating forward currents. The statistics for the maximum relative errors obtained by both parameter extraction processes are summarized in Table 3, where the statistics includes the relative errors of the model and the measurements for all operating points that are included in the extraction.
Consequently, the simulation results obtained using parameter sets identified by the previous method and those identified by the new, modified Nelder–Mead-based parameter extraction process agree very well in this current range.
The statistics show that the multi-domain LED model with parameter sets obtained through the new extraction methodology accurately describes the operation of various high-power LEDs, including color and phosphor-converted white types, with minimal errors. Specifically:
  • For the forward voltage, errors are below 0.5%.
  • For simulated radiant fluxes, errors are below 5%.
  • For the luminous flux, errors are below 8%.
These error levels are comparable to the industrial average uncertainties associated with measuring these LED properties.
The second benchmark evaluates how closely the measured LED characteristics are approximated by the multi-domain LED model equipped with the extracted parameter set. In these evaluations, we have to account for the effects of measurement uncertainties and noise, as their impact is absorbed into the overall error of the parameter fit. Since these measurement-related issues are consistent across any extraction process, and the theoretical accuracy of the multi-domain model is independent of the extraction method, the calculated error of the fit can be considered indicative solely of the quality of the parameter extraction process itself.
Figure 8 and Figure 9 illustrate the accuracy of parameter extraction by comparing the measurement and simulation results for two phosphor-converted white power LEDs—one with the best fit (Figure 8) and one with the worst fit (Figure 9). The parameter extraction for these LEDs was performed for the full sets of forward currents applied during measurements. However, for visibility, only the measurement values are plotted that have a junction temperature equal to the target. This allows us to display the isothermal lines and better present model errors. The extraction process, however, utilizes all measured values, which is a significant advantage of the new extraction methodology enabled by the new method. The difference of the results obtained by the model compared with the measurements in Figure 8 and Figure 9 shows only slight differences; however, the first has a maximum relative error of 1.5% and the second has a 4.7% radiant flux, but these differences appear in the low-current area. This confirms the hypothesis that the model adequately describes both the low- and high-current regions. There is only one original goal left that must be investigated: The new parameter extraction method should also effectively address sample-to-sample variations when the goal is to provide a single model parameter set representative of a given LED type or binning class. Because the extraction process is unknown for the specific LED samples, it can incorporate characteristics from all samples into the extended vertex used in the modified Nelder–Mead algorithm-based parameter extraction process. Figure 10 shows the extraction for one bin (white Nos. 1, 3, 4, 5, and 6). The resulting model provides a good fit, even with significant variance between the chips.
Note that the parameter sets obtained from the aggregate data of a larger number of LED samples of the same type differ from the parameter sets obtained by averaging individual samples. This highlights another advantage of our new approach: a significant reduction in required computational resources, particularly when extracting model parameter sets representative of an LED-type or binning class.
It is conspicuous but not surprising that the parameter sets obtained for multiple LED samples from the same type (and same binning class) provide a fit with more or less the same accuracy, as one LED can replace another, i.e., the model-to-measurement error is in the same range as the measurement-to-measurement error. This holds both for the prior brute-force method-based and for the new modified Nelder–Mead algorithm-based parameter extraction methods. Not only are the accuracies of the simulations similar for the LED samples of the same type, but the actual parameter values are also close to one another, e.g., the different I 0 values are within ± 10 % in a logarithmic scale. In the case of a larger population of phosphor-converted white LEDs of the same type (and presumed binning class), much larger sample-to-sample variations in parameter values were observed. This finding aligns with the previous results of the Delphi4LED project, which suggest that this population is a mix of two LED types exhibiting similar light output properties at the applied binning forward current and junction temperature [23].
The new parameter extractor was coded in Python 3.9. Extracting parameters for an LED sample measured at 33 operating points takes 10–18 s on a standard desktop computer (i5-12400, 3.2 GHz, 32 GB RAM) using a single core.

5. Conclusions

In this paper, a novel parameter extraction procedure was presented that is based on a modified Nelder–Mead optimization algorithm. The modifications greatly improved the robustness of the optimization process, as, due to the applied heuristics, the success rate increased from 50% to 100%.
The reproducibility of the new method was tested with different initializations. Over 10,000 extraction runs were performed on the measured IVL characteristics of 52 high-power LEDs from various types.
Compared to the previous parameter extraction process used in the previous Delphi4LED project, the new parameter extraction procedure provides multiple advantages, where, at least by an order of magnitude faster (reduced computational resource need), the original method lasts for 2–6 h (and needs manual operations). The new method finds the required data within a minute (typically 10–18 s), while the accuracy does not decrease or fall into the range of the supposed model and measurement accuracy for both optical and electrical domains, that is, 1% and 5% for maximum relative errors for the forward voltage and emitted flux, respectively.
The new method demonstrates a high level of robustness, achieving a 100% success rate in a benchmark test. This allows for an autonomous, unsupervised operation, making it suitable for automated web-based model parameter extraction services, as it can be found on https://extractor.ai-twilight.org/ (accessed on 10 August 2024), which has been publicly available and in use since March 2023 without any robustness issues.
Additionally, the new parameter extraction process simplifies the use of measurement data from less complex and faster LED IVL tests. It eliminates the need for standard isothermal characterization processes that require precise junction temperature control. Instead, it is sufficient to wait for any equilibrium state, rather than ensuring the junction temperature aligns with a reference point. This approach allows for the use of all measured data points, even if the reference junction temperature is not achieved, thereby reducing the number of measurements needed to create the model. Moreover, parameter extraction for a specific LED type can be completed in a single session using the test data from multiple samples of that type.

Author Contributions

Conceptualization, M.N. and A.P.; methodology, M.N., software, M.N.; validation, J.H., G.H.; formal analysis, M.N.; investigation, M.N.; resources, J.H. and G.H.; data curation, J.H. and G.H.; writing—original draft preparation, A.P. and M.N; writing—review and editing, M.N. and A.K.A.; visualization, M.N.; supervision, A.P.; project administration, A.P.; funding acquisition, A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by European Union’s Horizon 2020 research and innovation program, grant number: 101007319; Hungarian government, grant number: 2019-2.1.3-NEMZ_ECSEL-2021-00008.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The work presented received funding from the European Union’s Horizon 2020 research and innovation program through the H2020 ECSEL project AI-TWILIGHT (grant agreement number: 101007319). Co-financing of the AI-TWILIGHT project by the Hungarian government through the 2019-2.1.3-NEMZ_ECSEL-2021-00008 grant from the National Research, Development and Innovation Fund is also acknowledged. We are also thankful for the additional support from project no. TKP2021-NVA-02, which has been implemented with the support provided by the Ministry of Culture and Innovation of Hungary from the National Research, Development and Innovation Fund, financed under the TKP2021-NVA funding scheme.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. AI-TWILIGHT—AI Powered Digital Twin for Lighting Infrastructure in the Context of Front-End Industry 4.0 Website of the AI-TWILIGHT H2020 ECSEL Project of the EU. Available online: https://ai-twilight.eu (accessed on 15 January 2024).
  2. Sharma, M.; Pandey, D.; Khosla, D.; Goyal, S.; Pandey, B.K.; Gupta, A.K. Design of a GaN-based Flip Chip Light Emitting Diode (FC-LED) with au Bumps & Thermal Analysis with different sizes and adhesive materials for performance considerations. Silicon 2022, 14, 7109–7120. [Google Scholar] [CrossRef]
  3. Ye, J.; Mo, S.; Jia, L.; Chen, Y. Experimental performance of a LED thermal management system with suspended microencapsulated phase change material. Appl. Therm. Eng. 2022, 207, 118155. [Google Scholar] [CrossRef]
  4. Wang, Y.; Alonso, J.M.; Ruan, X. A review of LED drivers and related technologies. IEEE Trans. Ind. Electron. 2017, 64, 5754–5765. [Google Scholar] [CrossRef]
  5. Zhang, X.; Yin, L.; Ren, K.; Zhang, J. Research on Simulation Design of MOS Driver for Micro-LED. Electronics 2022, 11, 2044. [Google Scholar] [CrossRef]
  6. Zhu, Y.; Tang, T.; Zhao, S.; Joralmon, D.; Poit, Z.; Ahire, B.; Li, X. Recent advancements and applications in 3D printing of functional optics. Addit. Manuf. 2022, 52, 102682. [Google Scholar] [CrossRef]
  7. Sun, C.C.; Lee, T.X. Optical Design for LED Solid-State Lighting: A Guide; IOP Publishing: Bristol, UK, 2022; pp. 1–34. [Google Scholar] [CrossRef]
  8. Thermal Management for LED Applications; Springer: New York, NY, USA; Volume 2. Lasance, C.J.; Poppe, A. (Eds.) Springer: New York, NY, USA; Volume 2. [CrossRef]
  9. Lasance, C.J. Ten years of boundary-condition-independent compact thermal modeling of electronic parts: A review. Heat Transf. Eng. 2008, 29, 149–168. [Google Scholar] [CrossRef]
  10. Martin, G.; Marty, C.; Bornoff, R.; Poppe, A.; Onushkin, G.; Rencz, M.; Yu, J. Luminaire Digital Design Flow with Multi-Domain Digital Twins of LEDs. Energies 2019, 12, 2389. [Google Scholar] [CrossRef]
  11. Rosten, H.I.; Parry, J.D.; Lasance, C.J.M.; Vinke, H.; Temmerman, W.; Nelemans, W.; Assouad, Y.; Gautier, T.; Slattery, O.; Cahill, C.; et al. Final report to SEMITHERM XIII on the European-funded project DELPHI-the development of libraries and physical models for an integrated design environment. In Proceedings of the Thirteenth Annual IEEE. Semiconductor Thermal Measurement and Management Symposium, Austin, TX, USA, 28–30 January 1997; pp. 73–91. [Google Scholar] [CrossRef]
  12. Rosten, H.I.; Lasance, C.J.M.; Parry, J.D. The world of thermal characterization according to DELPHI-Part I: Background to DELPHI. IEEE Trans. Compon. Packag. Manuf. Technol. Part A 1997, 20, 384–391. [Google Scholar] [CrossRef]
  13. Lasance, C.J.M.; Rosten, H.I.; Parry, J.D. The world of thermal characterization according to DELPHI-Part II: Experimental and numerical methods. IEEE Trans. Compon. Packag. Manuf. Technol. Part A 1997, 20, 392–398. [Google Scholar] [CrossRef]
  14. Pape, H.; Noebauer, G. Generation and verification of boundary independent compact thermal models for active components according to the DELPHI/SEED methods. In Proceedings of the 15th IEEE Semiconductor Thermal Measurement and Management Symposium (SEMI-THERM’99, San Diego, CA, USA, 9–11 March 1999; pp. 201–211. [Google Scholar] [CrossRef]
  15. Rosten, H.I.; Lasance, C.J. DELPHI: The development of libraries of physical models of electronic components for an integrated design environment. In Model Generation in Electronic Design; Springer: Boston, MA, USA, 1995; pp. 63–90. [Google Scholar] [CrossRef]
  16. JEDEC Standard JESD15-1.01; Compact Thermal Model Overview. JEDEC: Arlington County, VI, USA, 2008. Available online: https://www.jedec.org/system/files/docs/JESD15-1.01.pdf (accessed on 15 January 2024).
  17. JEDEC Standard JESD15-4; DELPHI Compact Thermal Model Guideline. JEDEC: Arlington County, VI, USA, 2008. Available online: https://www.jedec.org/system/files/docs/JESD15-4.pdf (accessed on 15 January 2024).
  18. Lasance, C.J.M. Final report on the EC-funded thermal project PROFIT. In Proceedings of the 9th International Workshop on THERMal INvestigations of ICs and Systems, Aix-en-Provence, France, 24–26 September 2003; pp. 283–289. [Google Scholar] [CrossRef]
  19. Poppe, A.; Farkas, G.; Szekely, V.; Horvath, G.; Rencz, M. Multi-domain simulation and measurement of power LED-s and power LED assemblies. In Proceedings of the 22nd IEEE SEMI-THERM Symposium, Dallas, TX, USA, 14–16 March 2006; pp. 191–198. [Google Scholar] [CrossRef]
  20. Website of the Delpi4LED H2020 ECSEL Project. Available online: https://delphi4led.org (accessed on 15 January 2024).
  21. Public Deliverable Report D3.2 of the Delph4LED Project, “Report on SPICE-Like Multi-Domain Modelling of LEDs on Chip Level”. Available online: https://delphi4led.org/pydio/public/c4a8bb (accessed on 15 January 2024).
  22. Poppe, A. Simulation of LED Based Luminaires by Using Multi-Domain Compact Models of LEDs and Compact Thermal Models of their Thermal Environment. Microelectron. Reliab. 2017, 72, 65–74. [Google Scholar] [CrossRef]
  23. Merelle, T.; Sari, J.; Di Bucchianico, A.; Onushkin, G.; Bornoff, R.; Farkas, G.; Gaal, L.; Hantos, G.; Hegedus, J.; Poppe, A. Does a single LED bin really represent a single LED type? In Proceedings of the 29th Session of the CIE, Washington, DC, USA, 14–22 June 2019; pp. 1204–1214. [Google Scholar] [CrossRef]
  24. Piprek, J. Efficiency models for GaN-based light-emitting diodes: Status and challenges. Materials 2020, 13, 5174. [Google Scholar] [CrossRef] [PubMed]
  25. Pohl, L.; Kollár, E.; Poppe, A.; Kohári, Z. Nonlinear electro-thermal modeling and field-simulation of OLEDs for lighting applications I: Algorithmic fundamentals. Microelectron. J. 2012, 43, 624–632. [Google Scholar] [CrossRef]
  26. Kordt, P.; van der Holst, J.J.M.; Al Helwi, M.; Kowalsky, W.; May, F.; Badinski, A.; Lennartz, C.; Andrienko, D. Modeling of organic light emitting diodes: From molecular to device properties. Adv. Funct. Mater. 2015, 25, 1955–1971. [Google Scholar] [CrossRef]
  27. Bender, V.C.; Marchesan, T.B.; Alonso, J.M. Solid-state lighting: A concise review of the state of the art on LED and OLED modeling. IEEE Ind. Electron. Mag. 2015, 9, 6–16. [Google Scholar] [CrossRef]
  28. Kamalakis, T.; Walewski, J.W.; Ntogari, G.; Mileounis, G. Empirical Volterra-series modeling of commercial light-emitting diodes. J. Light. Technol. 2011, 29, 2146–2155. [Google Scholar] [CrossRef]
  29. Padmasali, A.; Kini, S. Experimental investigation and empirical modelling of thermal and drive current effect on optical performance of commercial LEDs. Light. Res. Technol. 2020, 53, 581–593. [Google Scholar] [CrossRef]
  30. Van der Schans, M.; Yu, J.; Martin, G. Digital Luminaire Design Using LED Digital Twins—Accuracy and Reduced Computation Time: A Delphi4LED Methodology. Energies 2020, 13, 4979. [Google Scholar] [CrossRef]
  31. Poppe, A.; Farkas, G.; Gaál, L.; Hantos, G.; Hegedüs, J.; Rencz, M. Multi-Domain Modelling of LEDs for Supporting Virtual Prototyping of Luminaires. Energies 2019, 12, 1909. [Google Scholar] [CrossRef]
  32. Bornoff, R. Extraction of boundary condition independent dynamic compact thermal models of LEDs—A Delphi4LED methodology. Energies 2019, 12, 1628. [Google Scholar] [CrossRef]
  33. JEDEC Standard JESD 51-51A; Implementation of the Electrical Test Method for the Measurement of Real Thermal Resistance and Impedance of Light-Emitting Diodes with Exposed Cooling Surface. JEDEC: Arlington County, VI, USA, 2022. Available online: https://www.jedec.org/system/files/docs/JESD51-51A.pdf (accessed on 15 January 2024).
  34. JEDEC Standard JESD 51-52A; Guidelines for Combining CIE 127:-2007/CIE 225:2017 Total Flux Measurements with Thermal Measurements of LEDs with Exposed Cooling Surface Guidelines For Combining CIE 127-2007 total Flux Measurements with Thermal Measurements of Leds with Exposed Cooling Surface. JEDEC: Arlington County, VI, USA, 2022. Available online: https://www.jedec.org/system/files/docs/JESD51-52A.pdf (accessed on 15 January 2024).
  35. CIE 225:2017 Technical Report: “Optical Measurement of High-Power LEDs”. 2017. Available online: https://cie.co.at/publications/optical-measurement-high-power-leds (accessed on 2 July 2024). [CrossRef]
  36. Zahara, E.; Kao, Y.-T. Hybrid Nelder–Mead simplex search and particle swarm optimization for constrained engineering design problems. Expert Syst. Appl. 2009, 36, 3880–3886. [Google Scholar] [CrossRef]
  37. Mesbahi, T.; Khenfri, F.; Rizoug, N.; Bartholomeus, P.; Moigne, P.L. Combined Optimal Sizing and Control of Li-Ion Battery/Supercapacitor Embedded Power Supply Using Hybrid Particle Swarm–Nelder–Mead Algorithm. IEEE Trans. Sustain. Energy 2017, 8, 59–73. [Google Scholar] [CrossRef]
  38. Fakhouri, H.N.; Hudaib, A.; Sleit, A. Hybrid Particle Swarm Optimization with Sine Cosine Algorithm and Nelder–Mead Simplex for Solving Engineering Design Problems. Arab. J. Sci. Eng. 2020, 45, 3091–3109. [Google Scholar] [CrossRef]
  39. Chelouah, R.; Siarry, P. A hybrid method combining continuous tabu search and Nelder–Mead simplex algorithms for the global optimization of multiminima functions. Eur. J. Oper. Res. 2005, 161, 636–654. [Google Scholar] [CrossRef]
  40. Chelouah, R.; Siarry, P. Genetic and Nelder–Mead algorithms hybridized for a more accurate global optimization of continuous multiminima functions. Eur. J. Oper. Res. 2003, 148, 335–348. [Google Scholar] [CrossRef]
  41. Public Deliverable Report 3.3 of the Delph4LED Project. “Report on Parameter Extraction for SPICE-like Chip Level Multi-Domain LED Models”. Available online: https://delphi4led.org/pydio/public/abef7a (accessed on 15 January 2024).
Figure 1. The SPICE circuit schematic representation of the Delphi4LED multi-domain chip-level LED model: (a) the branch representing the heat generation of the LED, (b) the branch describing the electrical behavior of the whole LED as a diode, and (c) the branch modeling the light emission of the LED. Lines indicate how one branch influences the other branches. The circuit symbols used are those typically found in SPICE circuit schematics.
Figure 1. The SPICE circuit schematic representation of the Delphi4LED multi-domain chip-level LED model: (a) the branch representing the heat generation of the LED, (b) the branch describing the electrical behavior of the whole LED as a diode, and (c) the branch modeling the light emission of the LED. Lines indicate how one branch influences the other branches. The circuit symbols used are those typically found in SPICE circuit schematics.
Applsci 14 07186 g001
Figure 2. Part of an IVL measurement result for a PC-white LED (No. 1). Junction temperatures that do not meet the reference temperature are indicated with the rectangles.
Figure 2. Part of an IVL measurement result for a PC-white LED (No. 1). Junction temperatures that do not meet the reference temperature are indicated with the rectangles.
Applsci 14 07186 g002
Figure 3. Selection of vectors for the initial simplex for the modified Nelder–Mead optimization method.
Figure 3. Selection of vectors for the initial simplex for the modified Nelder–Mead optimization method.
Applsci 14 07186 g003
Figure 4. The difference in the value of the objective function during the extraction of the parameters of the optical branch of the multi-domain LED model (involving 9 parameters) for different simplex sizes, compared to the best case with 50 points. Cases using 40, 30, 20, and 10 data points (original Nelder–Mead algorithm) are shown. (Left) side: Full scale. (Right) side: A maximum limit of 10% is applied. The baselines are shifted by −1% each for visual clarity.
Figure 4. The difference in the value of the objective function during the extraction of the parameters of the optical branch of the multi-domain LED model (involving 9 parameters) for different simplex sizes, compared to the best case with 50 points. Cases using 40, 30, 20, and 10 data points (original Nelder–Mead algorithm) are shown. (Left) side: Full scale. (Right) side: A maximum limit of 10% is applied. The baselines are shifted by −1% each for visual clarity.
Applsci 14 07186 g004
Figure 5. Maximum relative errors for the electrical branch of the multi-domain LED model. The errors are shown for parameter sets estimated using the prior brute-force method from the Delphi4LED project and the novel modified Nelder–Mead method. Results are presented for different LED types, for forward currents above 100 mA, and for the full range of forward currents, including the low current regime.
Figure 5. Maximum relative errors for the electrical branch of the multi-domain LED model. The errors are shown for parameter sets estimated using the prior brute-force method from the Delphi4LED project and the novel modified Nelder–Mead method. Results are presented for different LED types, for forward currents above 100 mA, and for the full range of forward currents, including the low current regime.
Applsci 14 07186 g005
Figure 6. Maximum relative errors for the optical branch of the multi-domain LED model with parameter sets estimated using the prior brute-force method, compared to the modified Nelder–Mead-based parameter extraction process. The comparison is made for different LED types, including forward currents above 100 mA and the full range of forward currents, including those below 100 mA.
Figure 6. Maximum relative errors for the optical branch of the multi-domain LED model with parameter sets estimated using the prior brute-force method, compared to the modified Nelder–Mead-based parameter extraction process. The comparison is made for different LED types, including forward currents above 100 mA and the full range of forward currents, including those below 100 mA.
Applsci 14 07186 g006
Figure 7. Maximum relative errors for the luminous flux of the multi-domain LED model with parameter sets estimated using the prior brute-force method, compared to the modified Nelder–Mead-based parameter extraction process. The comparison is made for different LED types, including forward currents above 100 mA and the full range of forward currents, including those below 100 mA.
Figure 7. Maximum relative errors for the luminous flux of the multi-domain LED model with parameter sets estimated using the prior brute-force method, compared to the modified Nelder–Mead-based parameter extraction process. The comparison is made for different LED types, including forward currents above 100 mA and the full range of forward currents, including those below 100 mA.
Applsci 14 07186 g007
Figure 8. Measurement results for a phosphor-converted white power LED chip (No. 1.) (*) and the simulation with parameters estimated by novel the Nelder–Mead method (solid line); electrical (left), radiant flux (middle), and luminous flux (right) characteristics.
Figure 8. Measurement results for a phosphor-converted white power LED chip (No. 1.) (*) and the simulation with parameters estimated by novel the Nelder–Mead method (solid line); electrical (left), radiant flux (middle), and luminous flux (right) characteristics.
Applsci 14 07186 g008
Figure 9. Measurement results for a phosphor-converted power LED chip (No. 8.) (*) and the simulation with parameters estimated by the novel Nelder–Mead method (solid line); electrical (left), radiant flux (middle), and luminous flux (right) characteristics.
Figure 9. Measurement results for a phosphor-converted power LED chip (No. 8.) (*) and the simulation with parameters estimated by the novel Nelder–Mead method (solid line); electrical (left), radiant flux (middle), and luminous flux (right) characteristics.
Applsci 14 07186 g009
Figure 10. Measurement results for aggregated 5 white (*) and the simulation with parameters estimated by the novel Nelder–Mead method (solid line); electrical (left), radiant flux (middle), and luminous flux (right) characteristics.
Figure 10. Measurement results for aggregated 5 white (*) and the simulation with parameters estimated by the novel Nelder–Mead method (solid line); electrical (left), radiant flux (middle), and luminous flux (right) characteristics.
Applsci 14 07186 g010
Table 1. Dependence of the success rate of the modified Nelder–Mead algorithm on the size of the simplex.
Table 1. Dependence of the success rate of the modified Nelder–Mead algorithm on the size of the simplex.
Number of VectorsSuccess Rate [%]
N = 10
(original method)
51.4
N = 2086.8
N = 3092.4
N = 4095.8
N = 5098.5
Table 2. Dependence of the success rate of the modified Nelder–Mead algorithm on the choice of the initial central points.
Table 2. Dependence of the success rate of the modified Nelder–Mead algorithm on the choice of the initial central points.
Simplex SizeSuccess Rate with an
Increased Central Point (10×) [%]
Success Rate with a
Decreased Central Point (10×) [%]
Success Rate with a
Decreased Central Point (5×) [%]
N = 1058.1015.7
N = 2089.315.237.5
N = 3094.541.459.9
N = 4095.255.872.3
N = 5098.878.295.1
Table 3. Relative errors (REs) statistics for 52 different sets of LED samples.
Table 3. Relative errors (REs) statistics for 52 different sets of LED samples.
MethodAvg. of REs,
Electric
[%]
Variance of REs,
Electric
[10−4]
Avg. of REs,
Optic
[%]
Variance of REs,
Optic
[10−4]
Avg. of REs,
Luminous
[%]
Variance of REs Luminous
[10−4]
Brute   force   I F 0.1 A0.110.0020.880.4275.1912.48
Nelder Mead   I F 0.1 A0.110.0010.430.0290.620.053
Brute force, all currents1.550.28912.4818.8410.6814.60
Nelder–Mead, all currents0.250.0052.691.403.891.93
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Németh, M.; Hegedüs, J.; Hantos, G.; Abdulrazzaq, A.K.; Poppe, A. Parameter Extraction for a SPICE-like Delphi4LED Multi-Domain Chip-Level LED Model with an Improved Nelder–Mead Method. Appl. Sci. 2024, 14, 7186. https://doi.org/10.3390/app14167186

AMA Style

Németh M, Hegedüs J, Hantos G, Abdulrazzaq AK, Poppe A. Parameter Extraction for a SPICE-like Delphi4LED Multi-Domain Chip-Level LED Model with an Improved Nelder–Mead Method. Applied Sciences. 2024; 14(16):7186. https://doi.org/10.3390/app14167186

Chicago/Turabian Style

Németh, Márton, János Hegedüs, Gusztáv Hantos, Ali Kareem Abdulrazzaq, and András Poppe. 2024. "Parameter Extraction for a SPICE-like Delphi4LED Multi-Domain Chip-Level LED Model with an Improved Nelder–Mead Method" Applied Sciences 14, no. 16: 7186. https://doi.org/10.3390/app14167186

APA Style

Németh, M., Hegedüs, J., Hantos, G., Abdulrazzaq, A. K., & Poppe, A. (2024). Parameter Extraction for a SPICE-like Delphi4LED Multi-Domain Chip-Level LED Model with an Improved Nelder–Mead Method. Applied Sciences, 14(16), 7186. https://doi.org/10.3390/app14167186

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop