Processing math: 100%
ARTICLE   Open Access    

Thermal hazard assessment and free radical inhibition of decomposition of tert-butyl perbenzoate

More Information
  • Tert-butyl perbenzoate (TBPB) is a common initiator widely used in polymerization processes, but the peroxide bond in its molecular structure is highly susceptible to breakage, leading to decomposition or even explosion. To explore the thermal behavior of TBPB and to inhibit the thermal hazard of free radicals generated during the reaction process, well-established calorimetric techniques were applied to measure the thermal stability of TBPB. The apparent activation energy of the TBPB decomposition reaction was also calculated using the Kissinger-Akahira-Sunose (KAS), Flynn-Wall-Ozawa (FWO), and Starink kinetic method. The thermal decomposition products of TBPB were determined by Fourier transform infrared spectroscopy (FTIR) experiment, and the qualitative analysis of free radicals generated during the reaction process was conducted by electron paramagnetic resonance spectroscopy (EPR) combined with free radical trapping technology. 2,2,6,6-tetramethylpiperidinooxy (TEMPO), a free radical trapping agent and inhibitor, was selected in this study as the thermal runaway inhibitor of the TBPB thermal decomposition reaction. Its inhibition effects on the corresponding free radicals and the thermal runaway of the TBPB decomposition reaction were verified. It is found that TEMPO can effectively reduce the potential thermal dangers and accident risks of TBPB, which provides a powerful reference for the prevention and management of thermal disasters during the production, storage, and transportation of TBPB.
  • In recent years, the widespread use of intelligent mobile terminals has been driving the process of societal intelligence. However, wired charging based on manual connections will severely limit this development[1]. Inductive power transfer (IPT) has emerged as a promising remedy for charging these intelligent terminals, enabling efficient power transmission from a transmitter (TX) to a receiver (RX). Presently, research focusing on single-TX single-RX IPT systems has achieved considerable maturity. Substantial advancements have been accomplished in coil design, compensation network configuration, output control, efficiency enhancement, and the detection of foreign objects[24]. These developments greatly support the expansion and application of this technology.

    Multi-RX charging has emerged as an attractive and challenging area of research. Varied application scenarios present distinct system characteristics, necessitating a focused approach towards system-level design adaptable to diverse power requirements and coupling conditions. In addressing these challenges, extensive discussions have centered around coupler designs, leading to the development of various structures aimed at enhancing and stabilizing the coupling amidst diverse misalignment scenarios[5,6]. Furthermore, the compensation networks have been recognized as crucial for attaining multiple objectives, including minimizing circulating energy, facilitating soft switching in active circuits, and ensuring output stability amidst load variations[7]. Notably, beyond passive circuitry, significant progress has been achieved in active circuit designs tailored for multi-RX scenarios[8,9]. These advancements in coupler and activating circuit designs serve as foundational requisites ensuring the steady-state performance of the system.

    A well-performed multi-RX system necessitates an intelligent control methodology to fully exploit its potential benefits. The key control objectives primarily encompass output regulation, power distribution among RXs, and efficiency optimization. For instance, in cases where each RX possesses a specific resonance frequency, employing multiple-frequency excitation allows for independent power control across multiple RXs[1012]. The overall efficiency is enhanced in the study by Kim et al.[13] through the computation of the optimal TX driving current, which heavily relies on the RX-side circuitry and real-time communication. Similarly, a communication-based control strategy has been devised for a system driven by a single Class E inverter, aiding in achieving simultaneous power and efficiency modulation[14]. Various input modes are formulated in the study of Zhao et al.[15] to optimize efficiency and cater to diverse power demands. However, these systems predominantly address specific scenarios, failing to comprehensively portray the inherent trade-offs and natural dependencies between achieved performance metrics (such as efficiency and power) and the requisite system complexity (e.g., real-time communication, information sensing, computational burden, and activation of circuits).

    This manuscript centers on upholding high efficiency within a single-TX multi-RX system, without knowing real coupling. Employing the single-TX scenario as a case study, the examination delves into the impact of terminal conditions—specifically, the driving current and load—on the resonant tank. This investigation then steers the discourse towards an in-depth analysis of existing approaches aimed at maximizing efficiency at different power. In the absence of real-time coupling knowledge, an optimal current is deduced based on pre-established coupling ranges. As the number of RX grows, the suggested driving logic exhibits considerable scalability and an elevation in system intricacy. The attained efficiency under complete knowledge of all couplings stands as the reference point, only marginally surpassing the efficiency attained through the proposed control strategy.

    The depicted single-TX single-RX IPT system, shown in Fig. 1, encompasses essential components such as the inverter, coupler, compensation unit, rectifier, and load, which collectively elucidate the fundamental control logic. The active circuits, including the inverter, rectifier, and DC/DC converters in both the front-end and post stages, enable input and output modulation for diverse control objectives. Within this system, Ltx and rtx denote the self-inductance and equivalent series resistance (ESR) of the TX coil, respectively, while Lrx and rrx represent the self-inductance and ESR of the RX coil. The symbol M signifies the mutual inductance between the coils, while the coefficient k is defined as the coil coupling coefficient: k=MLtxLrx.

    Figure 1.  Configuration of 1TX-1RX system.

    In Fig. 1, both the TX and RX sides should offer at least two control variables to adjust the output and optimize efficiency indirectly. Typically, these control variables encompass the duty cycle of the DC/DC converter, and the phase shift angle of the inverter, or rectifier. Different combinations of these variables yield various control configurations, as documented in prior studies[1618]. Under specific output conditions, such as maintaining a constant voltage (CV), the controllable circuits must satisfy the load's power requirements while ensuring high efficiency throughout the entire system. When dealing with a system incorporating an uncontrollable resonant tank, the diverse configurations and control strategies mentioned earlier can be simplified and clarified through a steady-state analysis of the central resonant tank. This analysis focuses specifically on the integration of compensation and coupler elements and hinges upon terminal conditions—specifically, the input from the TX side and the load output on the RX side.

    Despite the existence of various compensations, if assume that losses within the coupler exert a dominant influence on the overall losses, then the efficiency wouldn't be significantly impacted by compensation on the TX side. Consequently, this study is aimed at investigating the passive network depicted in Fig. 2, where TX compensation is disregarded, and a straightforward series compensation is implemented on the RX side. In this setup, Crx denotes the series compensating capacitor, fulfilling the condition 1/(ωCrx) = ωLrx, while itx represents the input current to the TX coil.

    Figure 2.  Equivalent circuit under different control strategies. (a) Optimal load. (b) Optimal excitation.

    In terms of the load viewpoint, controllable systems can be classified into two distinct types. The first type involves manipulating the RX-side circuitry to attain the optimal load while simultaneously adjusting the TX-side excitation to fulfill the output power requirements, as illustrated in Fig. 2a. In this arrangement, the rectifier and its subsequent circuit components are treated as a load resistance Rrx,opt when aiming for maximum coupler efficiency. The RX's output power during resonance is expressed as:

    Po=(ωMItxRrx,opt+rrx)2Rrx,opt (1)

    The optimal load resistance Rrx,opt for efficiency maximization has been well studied and derived as:

    Rrx,opt=rrx1+ω2M2rtxrrx (2)

    By solving Eqns (1) and (2) concurrently, it gives the root mean square (RMS) value of the necessary input current under the optimal condition, expressed as follows:

    Itx,ropt=PoωM2rrx+2rtxrrx+ω2M2rtxrrx+ω2M2rrxrtx (3)

    The efficiency of the coupler is:

    ηc,ropt=121+ω2M2rtxrrx+1 (4)

    which is loading and coupling dependent.

    The second type of system aims to provide optimal driving current to maximize efficiency, and its load model is illustrated in Fig. 2b. In this setup, the RX-side circuit maintains the desired output characteristics, such as a constant voltage (CV), while satisfying specific power requirements Po through RX-side control. Simultaneously, the TX-side circuit adjusts the input excitation to enhance efficiency. Despite the differing control methodologies between these two types, the eventual steady-state effects remain consistent. The primary difference lies in the approaches used to achieve this optimal state, involving distinctions in the activation of circuits, communication protocols, and the requisite information, particularly the coupling.

    In Fig. 2b, maximizing efficiency is equivalent to minimizing losses for a target power. The power loss is expressed as:

    Ploss=Ploss,tx+Ploss,rx=I2txrtx+I2rxrrx (5)

    The optimal excitation at this point can be determined by solving the equation dPlossdItx=0, yielding the coil's input current excitation Itx,iopt, which would further determine the corresponding optimum efficiency ηc,iopt. Based on derivation, it would find Itx,iopt = Itx,ropt and ηc,iopt = ηc,ropt. The fact that both types of systems share the same excitation in the steady state indicates that, despite different control methods and transient processes, the systems reach identical stable states, i.e., there exists only one global maximum efficiency point. If Rrx,opt rrx, the losses on the TX and RX coils are equal and can be expressed as:

    Ploss,txPloss,rx=Po1+ω2M2rtxrrx (6)

    It physically means the optimal efficiency indicates a loss balance condition.

    By employing optimal load control as illustrated in Fig. 2a, achieving Rrx,opt demands coupling information for the RX, while the TX relies on receiving feedback (usually via wireless communication) from the RX to regulate power[19]. Conversely, using the second approach depicted in Fig. 2b allows independent regulation of the RX through a linear controller. In the absence of communication, the TX must resort to either a non-linear tracking algorithm or a sensing-based coupling estimation approach to provide the necessary optimal driving current[20]. Consequently, a noticeable trade-off emerges involving the communication requirement, controller type, and sensing necessity.

    This paper presents a simplified excitation control approach based on the coupling range. The proposed method simplifies the control process by requiring determination of the coupling coefficient range during the coupler design phase. Following implementation, the system operates autonomously, obviating the necessity for real-time coupling estimation or nonlinear tracking algorithms. The core operational principle is demonstrated initially in a straightforward single-RX scenario, subsequently extending its applicability to a more intricate multi-RX scenarios.

    In a single-RX system, ensuring sufficient coupling generally requires the RX to be charged within a specified charging area. In such instances, the range of coupling coefficient variations is limited and becomes known upon the fabrication of the coupler. The controllable circuit on the RX side can uphold a constant voltage output (fulfilling the load power requirements), while the TX side employs a real-coupling-independent control. The excitation current is determined by the range of coupling variations.

    Figure 3 illustrates the efficiency curve of the coupler, demonstrating variations in input currents while regulating the RX side to achieve the target power. Within the graph, Mmin signifies the minimum mutual inductance, Mmax signifies the maximum mutual inductance, Po represents the output power, and the upper and lower boundaries are established by the coupling range. The solid black line represents the system response at maximum coupling, while the solid red line represents the system response at minimum coupling. Within the determined coupling range, the actual response curve of the system is denoted by the dashed blue line. Given offset conditions, the blue line remains constrained within the upper and lower boundaries.

    Figure 3.  Coupler efficiency at different k and Itx.

    In Fig. 3, L denotes the optimal efficiency point when the coupling is at its minimum. According to Eqn (3), this point is associated with an optimal input excitation, Itx,L, given by:

    Itx,L=PoωMmin2rrx+2rtxrrx+ω2M2minrtxrrx+ω2M2minrrxrtx (7)

    Here, H corresponds to the optimal point when the coupling is at its maximum, associated with an optimal excitation, Itx,H, given by:

    Itx,H=PoωMmax2rrx+2rtxrrx+ω2M2maxrtxrrx+ω2M2maxrrxrtx (8)

    If employing the non-linear perturbation observation method, the system requires continuous tracking of the optimal excitation between Itx,L and Itx,H[21]. Despite the absence of direct communication between the TX and RX sides, disturbances on the TX side can affect the stability of the RX side output. In scenarios without non-linear tracking, estimating the coupling based on TX side information becomes necessary to maintain the coupling-dependent optimal current[20].

    Within the specified coupling range, this paper proposes a strategy that utilizes the average of Itx,L and Itx,H as the driving current, denoted as:

    Itx,M = Itx,L+Itx,H2 (9)

    The median current value, Itx,M, is influenced by coil parameters, load power, and the coupling range. When the real coupling is Mmin, the operating point becomes L' using Itx,M as the driving current. The efficiency difference between L' and L is ΔηC,L. Conversely, at maximum coupling, the operating point is H', and the efficiency difference between H' and H is ΔηC,H. ΔηC,L and ΔηC,H depict the trade-off incurred by using the real-coupling-independent Itx,M.

    The proposed excitation control strategy eliminates the necessity for real-time tracking of the highest efficiency point and instead maintains operation at a compromise point. At this juncture, the TX side losses stabilize at I2tx,Mrtx, and different couplings determine varying currents on the RX-side coil. When the current on the RX-side coil reaches a certain value such that the losses at both ends of the coil become equal, i.e., I2tx,Mrtx=I2rxrrx, the system attains its actual maximum efficiency point.

    To provide a quantitative verification and study the impact of the coupling range, a simulation is employed for further elucidation of the qualitative analysis. The simulation utilizes the following parameters: f = 100 kHz, Ltx = 90 μH, Lrx = 75 μH, rtx = 0.192 Ω (Qtx = 250), rrx = 0.174 Ω (Qrx = 230), and Po = 20 W. As the coupling varies, the coupler efficiency changes with the input excitation, as depicted in Fig. 4. Region A delineates the operational area when k ranges within [0.1,0.2]. The simulation results indicate that with an increase in k, the corresponding optimal excitation current decreases.

    Figure 4.  Efficiency of the 1-RX system under different coupling conditions.

    If the driving current is controlled to be the optimal value, i.e., Itx = Itx,opt, the efficiency is depicted by the black dashed line in Fig. 5, which connects all peak points of Fig. 4. This reference curve serves as the basis for analyzing the influence when Itx changes to Itx,M. For instance, considering k [0.1, 0.2], the efficiency achieved by Itx,M—calculated using Eqns (7), (8), and (9)—is represented as case A (yellow solid line) in Fig. 5. It marginally falls below the reference curve, with a maximum drop of 0.5%. Expanding the coupling range to k [0.08, 0.2], corresponding to the region between the green and blue solid lines of Fig. 4, yields the efficiency achieved by Itx,M represented as case B in Fig. 5, exhibiting a maximum drop of 1%. When k [0.06, 0.2], as demonstrated in case C, it results in a maximum 2% efficiency drop. These simulation outcomes indicate that adopting Itx,M based on the coupling range demonstrates efficiency performance close to real-time tracking of the optimal excitation, showcasing favorable system characteristics.

    Figure 5.  Efficiency when the coupling range changes.

    Figure 6 illustrates a single-TX n-RX system. Lrx,j and rrx,j denote the self-inductance and ESR of the j-th RX coil. kj is the coupling coefficient between the TX and the j-th RX, satisfying: kj=Mj/LtxLrx,j. Crx,j is the series compensation capacitor, and Rrx,j is the equivalent AC load resistance. The resonance on the RX side satisfies the condition: ωLrx,j − 1/(ωCrx,j) = 0. In scenarios where charging occurs in a flat plane, having multiple RXs positioned over the TX charging area can induce cross-coupling between the RXs. Since the RX coil is enclosed within the device, the minimum coil clearance between adjacent devices becomes a variable during the design phase. Adequate clearance between coils is essential to minimize cross-coupling effects. Consequently, the analysis conducted above assumes a model where the RX cross-coupling is not considered.

    Figure 6.  Simplified circuit model of one-TX n-RX system.

    Similar to the single RX scenario depicted in Fig. 2, when employing different control logics in RXs, it becomes essential to individually examine the optimal conditions achieved. In scenarios where the system offers multiple control options, adjusting the load resistance at each RX becomes crucial to maximize coupler efficiency. As outlined in the study by Fu et al.[22] , the necessary load resistance values and their respective efficiencies are:

    Rrx,j,opt=rrx,j1+nj=1ω2M2jrtxrrx,j (10)
    ηc,ropt=121+nj=1ω2M2jrtxrrx,j+1 (11)

    In contrast to the single-RX system, the aforementioned multi-RX system possesses n power outputs and n optimal load resistance values. It offers n + 1 degrees of freedom for control, provided by the TX and RX sides (n on the RX side, 1 on the TX side). To achieve global optimal efficiency, i.e., satisfy Eqn (11), each load should adhere to the optimal load condition, i.e., Eqn (10), through its respective control circuit. In this scenario, the power distribution among the loads is constrained to follow the following proportions:

    P1:P2::Pn=M21rrx,1:M22rrx,2::M2nrrx,n (12)

    There are no degrees of freedom available for power regulation. Therefore, the control logic of Fig. 2a cannot be extended to a multi-RX scenario.

    When the RX-side circuit is dedicated to maintaining a constant voltage, the excitation current on the TX side becomes pivotal in enhancing the system's efficiency. This reflects the optimal excitation control employed in a scenario featuring a single TX and a single RX. In this context, maximizing the system's efficiency can be reformulated as a problem centered on minimizing losses while adhering to equality constraints, as articulated below:

    Minimize:   I2txrtx+nj=1(ωMjItxRrx,j+rrx,j)2rrx,jSubject to:   (ωMjItxRrx,j+rrx,j)2Rrx,j=Pj,    j=1,2,...,n (13)

    The method for solving such an optimization problem involves the Lagrange multiplier method. Equation (14) represents the Lagrangian, where L is the Lagrange function, and λi is the Lagrange multiplier.

    L=I2txrtx+nj=1(ωMjItxRrx,j+rrx,j)2rrx,j+nj=1λj((ωMjItxRrx,j+rrx,j)2Rrx,jPj) (14)

    If Rrx,optrrx, by solving LItx=0 and Lλj=0, the optimal input excitation current Itx, of the system under known coupling coefficients, which maximizes system efficiency, can be obtained as follows:

    Itx,opt=(1rtxnj=1P2jrrx,jω2M2j)14 (15)

    Substituting Eqn (15) into the efficiency expression yields:

    ηc,iopt=nj=1Pjnj=1Pj+2(nj=1P2jrrx,jω2M2j)12 (16)

    At this point, the coupler achieves an efficiency ηc,iopt while facing power constraints on the RX side, rather than the globally optimal efficiency ηc,ropt obtained without power constraints, expresses the system efficiency when operating at the global optimum, representing the theoretical efficiency limit of the system. Therefore, ηc,ioptηc,ropt, signifying that under the control strategy of optimal excitation, the efficiency is always lower than the efficiency under optimal load conditions. The losses on the TX and RX sides at this point are given by:

    Ploss,tx=Ploss,rx=(nj=1P2jrrx,jω2M2j)12 (17)

    This conclusion indicates that under optimal excitation, the coil conduction losses on the TX side and RX side are equal, consistent with the characteristics of a single-RX system [refer to Eqn (5)].

    From the above analysis, it can be concluded that for a general multi-RX system with different RX coils, various couplings, and distinct power requirements, when the coil parameters and output power are determined, the coupler efficiency depends solely on the coupling and the input excitation Itx. This characteristic is akin to the single-RX case. Despite the increase in the number of outputs, the control freedom of TX remains unchanged. The upper and lower boundaries of the efficiency curve are still determined by the coupling variation range. Therefore, the efficiency curve shown in Fig. 4 is still applicable to a multi-RX system. The red solid line represents the efficiency curve when the coupling between the TX and all RXs is minimized (i.e., Mj = Mmin). At this point, the current at the efficiency peak L is given by:

    Itx,L=(1rtxnj=1P2jrrx,jω2M2min)14 (18)

    The black solid line represents the efficiency curve when the coupling between the TX and all RXs is maximized (i.e., Mj = Mmax). The current at the corresponding efficiency peak is given by the following:

    Itx,H=(1rtxnj=1P2jrrx,jω2M2max)14 (19)

    It is possible to adopt the average current, which mirrors the approach used in the single-RX system, i.e.,

    Itx,M=Itx,L+Itx,H2 (20)

    which is also influenced by coil parameters, load power, and the coupling range. Given the known coupling range, the input excitation is determined without the need for real-time coupling information. When using the proposed excitation, the TX coil loss is Ploss,tx=I2tx,MRtx, while the RX-side losses vary with the coupling. When the coupling changes to the condition given by Eqn (17), where the sum of RX-side losses equals the TX side losses, the system reaches its optimal efficiency state, calculated as Eqn (16).

    In a single-RX system, the TX side must obtain specific coupling and load information through communication or detection means and then calculate the optimal excitation. Similarly, in a multi-RXs system, if coupling estimation or similar methods are employed, the difficulty of estimating information and the circuit cost will increase dramatically due to the addition of RX devices. Therefore, the proposed driving strategy would show great scalability and is particularly advantageous in multi-RX systems.

    Taking the example of a dual-RX system to study the influence of coupling, Table 1 shows the simulation parameters. The coupling range is k1 [0.1, 0.34] for RX1, and k2 [0.15, 0.45] for RX2. In Fig. 7, when meeting the power demand, the coupler efficiency changes with input excitation, and different colors represent different coupling conditions. The blue line represents the coil efficiency change curve when the coupling is strongest (k1 = 0.2, k2 = 0.25), the orange curve represents an intermediate state (k1 = 0.34, k2 = 0.45), and the black solid line is used to represent the weakest coupling (k1 = 0.1, k2 = 0.15). Similar to Fig. 4, stronger coupling leads to a smaller optimal excitation current and higher efficiency. Additionally, the power level and coil parameters also influence the shape and position of the efficiency curve.

    Table 1.  Parameters of 1TX-2RX IPT system.
    Symbol Value Symbol Value Symbol Value
    f 100 kHz Ltx 91.3 μH rtx & rrx 0.195 Ω
    Lrx,1 22.3 μH rrx,1 0.059 Ω k1 [0.1, 0.34]
    Lrx,2 26.4 μH rrx,2 0.074 Ω k2 [0.15, 0.45]
    P1 15 W P2 10 W
     | Show Table
    DownLoad: CSV
    Figure 7.  Coupler efficiency of a 2-RX system under different coupling conditions.

    Analogous to the single-RX system, for certain spatial offsets in practical applications, real-coupling-independent Itx,M is calculated by simultaneously solving Eqns (18), (19), and (20). Therefore, when Itx,M and Itx,opt are used, the efficiency difference can explain the high-efficiency operation.

    The coupling range denoted by region C of Fig. 7 means k1 [0.2, 0.34] and k2 [0.25, 0.45]. It is defined as Case C, and the corresponding simulation results are shown in Fig. 8a. The left graph illustrates the coupler efficiency ηc,M under excitation Itx,M, while the right graph displays the efficiency difference Δηc between the maximum coil efficiency ηc,iopt (using Itx,opt) and ηc,M (using Itx,M). Within the coupling range, the efficiency difference is less than 0.2%. As shown in Fig. 8b, with a further expansion of the coupling range in Case D (k1 [0.1, 0.34], k2 [0.15, 0.45]), the maximum efficiency drop is 0.8%. From these examples, it is evident that adopting the proposed current would maintain relatively high coupler efficiency without dramatically increasing the requirement of real-time coupling information.

    Figure 8.  Efficiency comparison. (a) Case C: ηs,M (left), Δηc (right); (b) Δηc (left), Δηc (right).

    An experimental platform was implemented using a single-TX dual-RX system to validate the theoretical conclusion. As illustrated in Fig. 9 the setup includes a DC source, a full-bridge inverter, a large square TX coil, two small RX coils (rectangular RX1 and square RX2), rectifiers, and loads. The TX side employs an LCC compensation structure, while the RX sides utilize series compensation. This compensation is convenient to use the DC input voltage for the TX coil current control due to the clamping effect. The coil parameters are provided in Table 1, and the compensation parameters are as follows: Lt = 32 μH, Ct = 110 nF, Ctx = 57 nF, Crx1 = 227 nF, Crx2 = 132.8 nF. The charging powers for the two loads are P1 = 10 W and P2 = 15 W.

    Figure 9.  Experimental setup.

    When two RXs are placed within a specific charging area in Fig. 9 (i.e., not moving outside the TX), the coupling variation ranges are measured as k1 [0.26, 0.34] and k2 [0.31, 0.45]. Using the driving current Eqn (20), with an effective value of Itx,M = 1.151 A, a corresponding DC voltage of Vdc = 23 V is provided. Note that the coil current is clamped by the input DC voltage in LCC compensation. Within the mentioned coupling variation range and power demands, the experimental waveforms of inverter output voltage vin, input current iin, and TX coil input current itx are depicted in Fig. 10. When employing Itx,M, vin, and itx remain constant, while iin exhibits slight fluctuations under coupling variation.

    Figure 10.  Experimental waveform.

    The efficiency analysis in the preceding section exclusively concentrates on the coupler itself, but the efficiency is also influenced by compensations, inverters, and rectifiers. And the characteristic traits of the efficiency curve primarily stem from the attributes of the coupler. Past studies have emphasized that a precise analysis of system efficiency necessitates the consideration of all loss models. However, when optimizing efficiency, the control mechanisms often resort to simplified analyses to deduce the terminal requisites for the resonant tank. Hence, in the ensuing analysis, a comprehensive system-level evaluation could be undertaken to yield results akin to those illustrated in Figs 7 and 8.

    Similar to the AC simulation in Fig. 7, the influence of the coupling could be studied in a complete DC to DC system. The overall efficiency ηs is measured by sweeping the input current as shown in Fig. 11. An optimal current exists, and it is still valid to have a conclusion that a stronger coupling leads to a lower optimal current.

    Figure 11.  Efficiency under different coupling conditions.

    When employing Itx,M, as illustrated in Fig. 10, and moving two RX coils to change their coupling, the input and output power were measured using a power analyzer. The system efficiency ηs,M is then obtained in Fig. 12a, which changes from 89.5% to 92.6%. Simultaneously, the optimal efficiency values ηc,sopt (obtained by sweeping the input current) under different couplings are measured and compared with ηs,M. The difference is illustrated in Fig. 12b. At higher coupling, ηs,M is closer to ηs,opt, differing by only 0.4%. As the coupling weakens, the efficiency difference gradually increases, reaching a maximum of 1.7%. The proposed control can maintain high-efficiency operations without the need for coupling estimation.

    Figure 12.  (a) System efficiency ηs,M. (b) Efficiency drop: Δηs = ηs,optηs,M.

    The similar trends in simulation and experiment indicate that the differences in efficiency are mainly attributed to the losses in the active circuits (inverter and rectifier) and conductor losses. The AC analysis in the previous sections would be used to guide the current control in a complete DC-to-DC system. The experimental results demonstrate that the proposed excitation based on the coupling range maintains a certain level of high efficiency, with stable input and simple control strategies.

    This paper proposes a simple driving logic to dramatically improve the scalability of a general single-TX multi-RX IPT system. The high-efficiency operation is achieved without knowing the real coupling and does not suffer from the coupling variation issue. Using the one-RX example, the coupler analysis unified the understanding of optimal conditions, which serves as the reference to justify the benefit of the proposed excitation. The same control logic is then extended to a general multi-RX case based on the derivation of optimal conditions. The claimed benefit is finally verified in a complete two-RX system.

  • The authors confirm contribution to the paper as follows: theoretical derivation, experimental validation, and manuscript writing: Wang X, He R; manuscript writing: Lin Z; conceptualization, funding acquisition, and supervision: Fu M. All authors reviewed the results and approved the final version of the manuscript.

  • The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

  • This work was supported by National Natural Science Foundation of China (Grant No. 52477013).

  • The authors declare that they have no conflict of interest.

  • [1]

    Zhou HL, Jiang JC, Huang AC. 2023. Thermal hazard assessment of tert-butyl perbenzoate using advanced calorimetric techniques and thermokinetic methods. Journal of Loss Prevention in the Process Industries 85:105166

    doi: 10.1016/j.jlp.2023.105166

    CrossRef   Google Scholar

    [2]

    Zhou HL, Jiang JC, Huang AC. 2024. Calorimetric and kinetic evaluation of thermal stability and process safety in tert-butyl peroxy-2-ethylhexanoate and tert-butyl peroxybenzoate. Process Safety and Environmental Protection 185:602−13

    doi: 10.1016/j.psep.2024.03.045

    CrossRef   Google Scholar

    [3]

    Huang AC, Huang CF, Xing ZX, Jiang JC, Shu CM. 2019. Thermal hazard assessment of the thermal stability of acne cosmeceutical therapy using advanced calorimetry technology. Process Safety and Environmental Protection 131:197−204

    doi: 10.1016/j.psep.2019.09.016

    CrossRef   Google Scholar

    [4]

    Milas NA, Surgenor DM. 1946. Studies in organic peroxides; t-amyl hydroperoxide and di-t-amyl peroxide. Journal of the American Chemical Society 68:643

    doi: 10.1021/ja01208a034

    CrossRef   Google Scholar

    [5]

    Wei W, Zhang C, Xu Y, Wan X. 2011. Synthesis of tert-butyl peresters from aldehydes by Bu4NI-catalyzed metal-free oxidation and its combination with the Kharasch–Sosnovsky reaction. Chemical Communications 47:10827−29

    doi: 10.1039/c1cc14602e

    CrossRef   Google Scholar

    [6]

    Zhang H, Dong DQ, Hao SH, Wang ZL. 2016. Bu4NI-catalyzed construction of tert-butyl peresters from alcohols. RSC Advances 6:8465−68

    doi: 10.1039/C5RA27500H

    CrossRef   Google Scholar

    [7]

    Chen X, Li Y, Wu M, Guo H, Jiang L, et al. 2016. An efficient method for the preparation of tert-butyl esters from benzyl cyanide and tert-butyl hydroperoxide under the metal free condition. RSC Advances 6:102023−27

    doi: 10.1039/C6RA20966A

    CrossRef   Google Scholar

    [8]

    Hashemi H, Saberi D, Poorsadeghi S, Niknam K. 2017. Temperature-controlled solvent-free selective synthesis of tert-butyl peresters or acids from benzyl cyanides in the presence of the TBHP/Cu(OAc)2 system. RSC Advances 7:7619−22

    doi: 10.1039/C6RA27921J

    CrossRef   Google Scholar

    [9]

    Singha R, Shit P. 2020. Sunlight assisted solvent free synthesis of tert-butylperesters. Synthetic Communications 50:2698−703

    doi: 10.1080/00397911.2020.1783560

    CrossRef   Google Scholar

    [10]

    Tseng JM, Lin YF. 2011. Evaluation of a tert-butyl peroxybenzoate runaway reaction by five kinetic models. Industrial & Engineering Chemistry Research 50:4783−87

    doi: 10.1021/ie100640t

    CrossRef   Google Scholar

    [11]

    Lv JY, Wei S, Chen WH, Chen GF, Chen LP, et al. 2012. Thermal kinetic analysis of tert-butyl peroxybenzoate under dynamic and adiabatic conditions. Advanced Materials Research 550−553:2782−85

    doi: 10.4028/www.scientific.net/AMR.550-553.2782

    CrossRef   Google Scholar

    [12]

    Wei TT, Qian XM, Yuan MQ. 2015. Thermal hazard analysis for tert-butyl peroxybenzoate contaminated by acid or alkali. CIESC Journal 66(10):3931−39

    doi: 10.11949/j.issn.0438-1157.20141395

    CrossRef   Google Scholar

    [13]

    Jiang JC, Li L, Jiang JJ, Wang Y, Lo SM, et al. 2019. Effect of ionic liquids on the thermal decomposition process of tert-butyl peroxybenzoate (TBPB) by DSC. Thermochimica Acta 671:127−33

    doi: 10.1016/j.tca.2018.11.017

    CrossRef   Google Scholar

    [14]

    Moane S, Raftery DP, Smyth MR, Leonard RG. 1999. Decomposition of peroxides by transition metal ions in anaerobic adhesive cure chemistry. International Journal of Adhesion and Adhesives 19:49−57

    doi: 10.1016/S0143-7496(98)00056-6

    CrossRef   Google Scholar

    [15]

    Lv J, Chen W, Chen L, Tian Y, Yan J. 2014. Thermal risk evaluation on decomposition processes for four organic peroxides. Thermochimica Acta 589:11−18

    doi: 10.1016/j.tca.2014.05.013

    CrossRef   Google Scholar

    [16]

    Dobbs AP, Jones P, Penny MJ, Rigby SE. 2009. Light-fluorous TEMPO: reagent, spin trap and stable free radical. Tetrahedron 65:5271−77

    doi: 10.1016/j.tet.2009.04.078

    CrossRef   Google Scholar

    [17]

    Barton DHR, Le Gloahec VN, Smith J. 1998. Study of a new reaction: trapping of peroxyl radicals by TEMPO. Tetrahedron Letters 39:7483−86

    doi: 10.1016/S0040-4039(98)01628-1

    CrossRef   Google Scholar

    [18]

    Makino K, Hagiwara T, Murakami A. 1991. A mini review: Fundamental aspects of spin trapping with DMPO. International Journal of Radiation Applications and Instrumentation. Part C: Radiation Physics and Chemistry 37:657−65

    doi: 10.1016/1359-0197(91)90164-W

    CrossRef   Google Scholar

    [19]

    Buettner GR. 1993. The spin trapping of superoxide and hydroxyl free radicals with DMPO (5,5-dimethylpyrroline-N-oxide): more about iron. Free Radical Research Communications 19:s79−s87

    doi: 10.3109/10715769309056s79

    CrossRef   Google Scholar

    [20]

    Villamena FA, Locigno EJ, Rockenbauer A, Hadad CM, Zweier JL. 2006. Theoretical and experimental studies of the spin trapping of inorganic radicals by 5,5-dimethyl-1-pyrroline N-oxide (DMPO). 1. carbon dioxide radical anion. The Journal of Physical Chemistry A 110:13253−58

    doi: 10.1021/jp064892m

    CrossRef   Google Scholar

    [21]

    Yao H, Jiang J, Li B, Ni L, Ni Y, et al. 2022. Investigation of pyrolysis kinetics, mechanism and thermal stability of tert-butyl peroxy-2-ethyl hexanoate. Process Safety and Environmental Protection 160:734−48

    doi: 10.1016/j.psep.2022.02.059

    CrossRef   Google Scholar

    [22]

    Tang Y, Li ZP, Zhou HL, Miao CF, Jiang JC, et al. 2023. Thermal stability assessment of nitrocellulose by using multiple calorimetric techniques and advanced thermokinetics. Journal of Thermal Analysis and Calorimetry 148:5029−38

    doi: 10.1007/s10973-022-11754-1

    CrossRef   Google Scholar

    [23]

    Huang AC, Li ZP, Liu YC, Tang Y, Huang CF, et al. 2021. Essential hazard and process safety assessment of para-toluene sulfonic acid through calorimetry and advanced thermokinetics. Journal of Loss Prevention in the Process Industries 72:104558

    doi: 10.1016/j.jlp.2021.104558

    CrossRef   Google Scholar

    [24]

    Liu YC, Zhou HL, Tang Y, Li Y, Zhai J, et al. 2023. Thermal hazard assessment by TGA, DSC, and ARC experimental and simulated thermokinetic approaches for trinitrophloroglucinol. Journal of Thermal Analysis and Calorimetry 148:5039−49

    doi: 10.1007/s10973-022-11649-1

    CrossRef   Google Scholar

    [25]

    Li ZP, Huang AC, Tang Y, Zhou HL, Liu YC, et al. 2022. Thermokinetic prediction and safety evaluation for toluene sulfonation process and product using calorimetric technology. Journal of Thermal Analysis and Calorimetry 147:12177−86

    doi: 10.1007/s10973-022-11384-7

    CrossRef   Google Scholar

    [26]

    Alavi SE, Ali Abdoli M, Khorasheh F, Bayandori Moghaddam A. 2020. Non-isothermal pyrolysis of used lubricating oil and the catalytic effect of carbon-based nanomaterials on the process performance. Journal of Thermal Analysis and Calorimetry 139:1025−36

    doi: 10.1007/s10973-019-08436-w

    CrossRef   Google Scholar

    [27]

    Wu ZH, Huang AC, Tang Y, Yang YP, Liu YC, et al. 2021. Thermal effect and mechanism analysis of flame-retardant modified polymer electrolyte for lithium-Ion battery. Polymers 13:1675

    doi: 10.3390/polym13111675

    CrossRef   Google Scholar

    [28]

    Zhou HL, Jiang JC, Huang AC, Tang Y, Zhang Y, et al. 2022. Calorimetric evaluation of thermal stability and runaway hazard based on thermokinetic parameters of O, O–dimethyl phosphoramidothioate. Journal of Loss Prevention in the Process Industries 75:104697

    doi: 10.1016/j.jlp.2021.104697

    CrossRef   Google Scholar

    [29]

    Li X, Yao H, Lu X, Chen C, Cao Y, et al. 2020. Effect of pyrogallol on the ring-opening polymerization and curing kinetics of a fully bio-based benzoxazine. Thermochimica Acta 694:178787

    doi: 10.1016/j.tca.2020.178787

    CrossRef   Google Scholar

    [30]

    Cao CR, Liu SH, Huang AC, Lee MH, Ho SP, et al. 2018. Application of thermal ignition theory of di(2, 4-dichlorobenzoyl) peroxide by kinetic-based curve fitting. Journal of Thermal Analysis and Calorimetry 133:753−61

    doi: 10.1007/s10973-018-7002-8

    CrossRef   Google Scholar

    [31]

    Wei R, Huang S, Wang Z, Wang C, Zhou T, et al. 2018. Effect of plasticizer dibutyl phthalate on the thermal decomposition of nitrocellulose. Journal of Thermal Analysis and Calorimetry 134:953−69

    doi: 10.1007/s10973-018-7521-3

    CrossRef   Google Scholar

    [32]

    Li ZP, Jiang JC, Huang AC, Tang Y, Miao CF, et al. 2021. Thermal hazard evaluation on spontaneous combustion characteristics of nitrocellulose solution under different atmospheric conditions. Scientific Reports 11:24053

    doi: 10.1038/s41598-021-03579-z

    CrossRef   Google Scholar

    [33]

    Huang AC, Huang CF, Tang Y, Xing ZX, Jiang JC. 2021. Evaluation of multiple reactions in dilute benzoyl peroxide concentrations with additives using calorimetric technology. Journal of Loss Prevention in the Process Industries 69:104373

    doi: 10.1016/j.jlp.2020.104373

    CrossRef   Google Scholar

    [34]

    Vyazovkin S, Burnham AK, Criado JM, Pérez-Maqueda LA, Popescu C, et al. 2011. ICTAC Kinetics Committee recommendations for performing kinetic computations on thermal analysis data. Thermochimica Acta 520:1−19

    doi: 10.1016/j.tca.2011.03.034

    CrossRef   Google Scholar

    [35]

    Shen S, Jiang J, Zhang W, Ni L, Shu CM. 2018. Process safety evaluation of the synthesis of tert-butyl peracetate. Journal of Loss Prevention in the Process Industries 54:153−62

    doi: 10.1016/j.jlp.2018.03.009

    CrossRef   Google Scholar

    [36]

    Zhang H, Jiang J, Fei M, Ni L, Hang Y. 2022. Thermal hazard characteristics and essential mechanism study of 1-hydroxybenzotriazole: Thermodynamic study combined DFT simulation. Process Safety and Environmental Protection 168:713−22

    doi: 10.1016/j.psep.2022.10.043

    CrossRef   Google Scholar

    [37]

    Zhang H, Jiang JC, Yan TY, Ni L, Liu SH. 2023. Thermal hazard risk and decomposition mechanism identification of 1-Hexyl-2,3-dimethylimidazolium nitrate: combined thermal analysis experiment and DFT emulation. Process Safety and Environmental Protection 172:38−47

    doi: 10.1016/j.psep.2023.01.065

    CrossRef   Google Scholar

  • Cite this article

    Zhang D, Li Z, Jiang J, Ni L, Chen Z. 2025. Thermal hazard assessment and free radical inhibition of decomposition of tert-butyl perbenzoate. Emergency Management Science and Technology 5: e001 doi: 10.48130/emst-0024-0029
    Zhang D, Li Z, Jiang J, Ni L, Chen Z. 2025. Thermal hazard assessment and free radical inhibition of decomposition of tert-butyl perbenzoate. Emergency Management Science and Technology 5: e001 doi: 10.48130/emst-0024-0029

Figures(12)  /  Tables(4)

Article Metrics

Article views(881) PDF downloads(299)

ARTICLE   Open Access    

Thermal hazard assessment and free radical inhibition of decomposition of tert-butyl perbenzoate

Emergency Management Science and Technology  5 Article number: e001  (2025)  |  Cite this article

Abstract: Tert-butyl perbenzoate (TBPB) is a common initiator widely used in polymerization processes, but the peroxide bond in its molecular structure is highly susceptible to breakage, leading to decomposition or even explosion. To explore the thermal behavior of TBPB and to inhibit the thermal hazard of free radicals generated during the reaction process, well-established calorimetric techniques were applied to measure the thermal stability of TBPB. The apparent activation energy of the TBPB decomposition reaction was also calculated using the Kissinger-Akahira-Sunose (KAS), Flynn-Wall-Ozawa (FWO), and Starink kinetic method. The thermal decomposition products of TBPB were determined by Fourier transform infrared spectroscopy (FTIR) experiment, and the qualitative analysis of free radicals generated during the reaction process was conducted by electron paramagnetic resonance spectroscopy (EPR) combined with free radical trapping technology. 2,2,6,6-tetramethylpiperidinooxy (TEMPO), a free radical trapping agent and inhibitor, was selected in this study as the thermal runaway inhibitor of the TBPB thermal decomposition reaction. Its inhibition effects on the corresponding free radicals and the thermal runaway of the TBPB decomposition reaction were verified. It is found that TEMPO can effectively reduce the potential thermal dangers and accident risks of TBPB, which provides a powerful reference for the prevention and management of thermal disasters during the production, storage, and transportation of TBPB.

    • Tert-butyl perbenzoate (TBPB) is a common organic peroxide, which is chemically active and has strong oxidation properties, so it is widely used as an initiator to participate in various polymerization reactions and is often used as an oxidizer to participate in peroxide and oxidation synthesis processes[1]. The molecular structure of TBPB (Fig. 1) contains a peroxide bond (−O−O−), which is highly sensitive to heat, impact, friction, and other causes. These can lead to the breakage of peroxy bond and thermal decomposition of TBPB, which in extreme circumstances can even result in combustion and explosion[2,3]. Therefore, the thermal instability of TBPB brings some difficulties and challenges in its preparation, and meanwhile, also brings some potential security risks in its storage, transportation, and application.

      Figure 1. 

      The chemical structure of TBPB.

      To solve the problems of complex processes, high thermal hazard, and low product conversion rate, many scholars were committed to developing different synthesis methods of TBPB. Previously, Milas & Surgenor reported that TBPB could be prepared by reacting carboxylic acid or benzoyl chloride with tert-butyl hydroperoxide (TBHP) in alkaline conditions[4]. Wei et al. developed a process route of preparing TBPB directly from aldehyde and TBHP in a water environment, which effectively improved the production rate of TBPB[5]. Subsequently, some studies were conducted to optimize the process conditions on this basis[6]. Then Chen et al. proposed a new efficient synthesis method of TBPB by phenylacetonitrile and TBHP in a metal-free and nitrogen environment[7], and later the solvent-free synthesis method of TBPB at room temperature was realized by Hashemi et al.[8]. Through a series of studies by scholars, the green, efficient, and safe preparation of TBPB is gradually developing[9].

      At the same time, many other scholars have explored the behavior and products of TBPB's thermal decomposition utilizing thermal analysis, calculated the kinetic parameters, and evaluated the thermal risk of the reaction[1,10,11]. Other research has studied the effects of common impurities in industrial processes such as acids, alkalis[12], ionic liquids[13], metal ions[14], and other peroxides[2,15] on the thermal stability and decomposition kinetics of TBPB.

      However, a lot of existing research mainly focuses on the apparent experimental perspective, and there are still deficiencies in analyzing the thermal decomposition mechanism of TBPB in different environments from the molecular microscopic level. There have been many reports on the thermal hazard and thermal runaway risk of TBPB, but there are few studies on the inhibition of thermal decomposition reaction at the microscopic level. Due to the lack of research in this part, this paper not only studied the thermal decomposition behavior and the apparent activation energy (Ea) of TBPB through traditional differential scanning calorimetry (DSC), accelerating rate calorimetry (ARC), and kinetic calculation but also used the Fourier transform infrared spectroscopy (FTIR), the electron paramagnetic resonance spectroscopy (EPR) combined with free radical trapping technology to study the generation of products and free radicals in the decomposition process of TBPB. Then, the free radical inhibitor was added to compare and analyze the inhibition effect of free radicals and thermal runaway.

      By further understanding the thermal decomposition behavior of TBPB and exploring the free radicals, functional groups, and decomposition products generated in the reaction process, this paper not only provides a strong reference for screening suitable free radical inhibitors of thermal decomposition but also has important significance for preventing thermal runaway of TBPB in practical application.

    • TBPB was purchased from Shanghai Macklin Biochemical Technology Co. Ltd. (Shanghai, China) at a concentration of 98% and was not required for purification when used. The free radical inhibitor selected was 2,2,6,6-tetramethylpiperidoxyl (TEMPO)[16,17], also purchased from the same company. TEMPO, illustrated in Fig. 2a, is a piperidine-type nitrogen oxide radical with the function of capturing free radicals and quenching singlet-oxygen. TEMPO was an orange crystal and had good solubility in TBPB. TBPB/TEMPO can be obtained by completely dissolving 3.0 mg TEMPO into 0.5 mL TBPB.

      Figure 2. 

      The chemical structures of (a) TEMPO, and (b) DMPO.

      In the EPR experiment, 5,5-dimethyl-1-pyrroline N-oxide (DMPO) (Fig. 2b) was used as a spin-trapping agent to combine with different free radicals generated during the thermal decomposition process of TBPB, forming different spin-trapping adducts[1820]. These adducts showed different absorption peaks after EPR detection, and the captured free radical information can be obtained by analyzing the EPR spectra.

      All materials were stored in the refrigerator, with TBPB and TEMPO stored at 2–8 °C and DMPO stored at −40 °C.

    • To study the exothermic behavior of TBPB's thermal decomposition and obtain the characteristic temperatures for calculating the Ea, the TA Q20 DSC manufactured in the United States was used to perform thermal decomposition experiments of TBPB under a nitrogen environment with a flow rate of 50 mL/min. Five different heating rates (β) were set at 2.0, 4.0, 6.0, 8.0, and 10.0 °C/min, and the experimental temperature was controlled within the range of 50–300 °C. All samples were placed in gold-plated crucibles, with an experimental sample mass of approximately 4.0 mg.

      To investigate the inhibitory effect of TEMPO on the exothermic reaction of TBPB, and to compare the differences in the exothermic behavior of TBPB's thermal decomposition reaction before and after adding the inhibitor, DSC was applied again to conduct calorimetric experiments on the TBPB/TEMPO. For the comparative experiments, the sample mass of about 5.0 mg was taken, and the gold-plated crucibles were also used under the same nitrogen environment. The sample was heated at β of 10.0 °C/min in the same range of 50–300 °C.

    • To collect and detect the products during the thermal decomposition process of TBPB in real-time, the FTIR and thermogravimetry from Netzsch company were combined to perform experiments[21]. In the test, about 5.0 mg TBPB was placed in an alumina crucible and the thermogravimetry was operated to heat the sample at β of 10.0 °C/min in the temperature range of 40–500 °C, with the same nitrogen conditions as in the DSC experiments. The thermal decomposition products of TBPB were collected through a tube connected to the thermogravimetry and detected by FTIR, then analyzed by the infrared absorption spectra.

    • ARC from the Young Instruments Company was adopted to test the thermal runaway behaviors of TBPB and TBPB/TEMPO under adiabatic conditions. Hastelloy sample balls were used in the experiments, and the sample mass was about 1.0 g. The testing method was selected as H-W-S mode (heat-wait-search)[2224], and the experimental temperature was also set as 50–300 °C with β of 5.0 °C/min.

    • The JES-X320 EPR from the JEOL company from Japan was employed for real-time monitoring of free radicals generated during the thermal decomposition process of TBPB and TBPB/TEMPO. The test sample was made by mixing TBPB or TBPB/TEMPO and the spin trapping agent DMPO in the volume ratio of 10:1. An appropriate volume of the sample was taken with a capillary and placed in a paramagnetic tube and put into the resonator of the EPR for microwave tuning. The experimental temperature was set at 110 °C by the variable temperature system of the EPR, then a suitable magnetic field range was selected, and finally, the microwave power was set at 1.0 mW to start the detection of free radicals generated during the decomposition of different samples. The electron paramagnetic resonance spectra produced by the adducts of free radicals and DMPO during the thermal decomposition reaction of TBPB and TBPB/TEMPO at 110 °C can be obtained by EPR experiments. Subsequently, by identifying EPR spectra morphology and analyzing the characteristic parameters such as the g-factor (giso) of the resonance absorption peaks, the hyperfine coupling constants of the hydrogen nuclei (AH), and that of the nitrogen nuclei (AN), the types of the radicals can be deduced.

    • In this study, three common isoconversional model-free methods, Kissinger-Akahira-Sunose (KAS) method, Flynn-Wall-Ozawa (FWO) method, and Starink method, were used to calculate Ea to determine the difficulty of TBPB's thermal decomposition reaction. The temperatures at different conversion rates (α = 0.10, 0.20, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80, and 0.90) were substituted into the calculation respectively, and then linear fitting was performed to obtain the slope of the straight line, which was the Ea of the reaction.

    • The KAS method is shown in Eqn (1)[2528]:

      ln(βT2)=lnAREaG(α)EaR1T (1)

      where, A is the pre-exponential factor, R is the universal gas constant (8.314 J/(mol K)), T is the absolute temperature (K) at the given α, and G(α) is the integral form of the kinetic mechanism function.

      Through the linear relationship between ln(β/T2) and 1/T, the corresponding slope and Ea can be derived.

    • The FWO method is a typical integral kinetic method and can be calculated as shown in Eqn (2)[2931].

      lgβ=lg(AEaRG(α))2.3150.4567EaRT (2)

      This method involves fewer variables and is more convenient to calculate, so it was widely used by many kinetic studies[32].

    • The Starink method is a new method obtained through continuous optimization on the basis of the Kissinger method, which is presented in Eqn (3)[3335].

      ln(βT1.8)=CS1.0037EaR1T (3)

      where Cs is a constant.

      The reaction Ea can be calculated by fitting the slope of the straight line from ln(β/T1.8) to 1/T.

    • As shown in Fig. 3, five heat flow curves at different heating rates were obtained by taking the reaction temperature as the X-axis and the heat flow, i.e., the heat release rate, as the Y-axis. From the curves, it can be seen that TBPB thermally decomposed and was exothermic in the temperature range of 100–210 °C. With the increase of β, the exothermic curves of TBPB shifted to high temperature on the right side of the horizontal coordinate. Comparing the maximum rate decomposition temperature (Tp) at each β, it can be seen that the Tp of TBPB showed an increasing trend with the increase of β, which gradually rose from 141.01 to 165.60 °C. When the temperature reached about 150 °C, the exothermic rate of the reaction reached its maximum and the decomposition reaction was at its most severe. By integrating the heat flow curve with respect to reaction time, the heat release (ΔH) throughout the whole decomposition at five different β can be obtained. The average ΔH during the whole decomposition could reach 924.59 J/g.

      Figure 3. 

      Exothermic DSC curves of TBPB at five β.

      From the above analysis, it can be concluded that TBPB has a high thermal hazard, once the thermal runaway occurs, it will release a large amount of heat and cause serious damage. Therefore, more attention should be paid to the strict control of ambient temperature during the production, storage, and transportation of TBPB.

    • As seen in Fig. 4, a linear fit was made using the KAS method based on the DSC experimental data. The slopes of the matching fitted straight lines were used to determine the Ea at various α. The average activation energy (¯Ea) of the TBPB thermal decomposition reaction was then obtained using the KAS method to be 96.80 kJ/mol by averaging all Ea values at different α. Figure 5 displays the results of linear fits using the FWO method with 1000/T as the horizontal coordinate and lgβ as the vertical coordinate. Following the aforementioned mathematical procedure, the ¯Ea value acquired by the FWO method was 98.72 kJ/mol. Similarly, the Starink method was also applied to calculate Ea values at each α of TBPB thermal decomposition, and the fitting results are shown in Fig. 6. The ¯Ea computed by this method was 97.14 kJ/mol.

      Figure 4. 

      The fitting results of TBPB's Ea by the KAS method.

      Figure 5. 

      The fitting results of TBPB's Ea by the FWO method.

      Figure 6. 

      The fitting results of TBPB's Ea by the Starink method.

      To observe the Ea change trend of TBPB thermal decomposition reaction at different α, the Ea change curves determined by three kinetic methods were plotted in Fig. 7. From the trends of the three curves in the figure, it can be seen that the Ea of TBPB decomposition reaction gradually decreases with the increasing α values. Therefore, it can be assumed that the energy required for the thermal decomposition reaction of TBPB decreases as the reaction continues. The Ea values calculated by the three methods (KAS, FWO, and Starink method) were very close to each other, especially the Ea curves of KAS and Starink methods were almost coincident. According to the calculation results listed in Table 1, the average fitting correlation coefficients (¯R2) of the three methods reached more than 0.96, proving that these three kinetic methods are well applicable to the Ea calculation of TBPB thermal decomposition reaction.

      Figure 7. 

      Changes of TBPB's Ea under different α.

      Table 1.  Ea values of TBPB obtained through three kinetic methods.

      Parameters Kinetic methods Average values
      KAS FWO Starink
      ¯Ea (kJ/mol) 96.80 98.72 97.14 97.55
      ¯R2 0.9636 0.9681 0.9641 0.9653

      Therefore, after averaging the values of ¯Ea obtained by the methods mentioned in Table 1, the activation energy for the thermal decomposition reaction of TBPB can be determined as 97.55 kJ/mol, which indicates that the reaction occurs relatively easily, and it is necessary to take effective measures to prevent the occurrence of TBPB thermal decomposition.

    • During the FTIR experiments, the gaseous products of TBPB were collected in real-time and scanned with several infrared spectra. Three typical temperatures were selected to analyze the products, and the infra-red spectra of the thermal decomposition products of TBPB under nitrogen conditions were plotted with wavenumber as the horizontal coordinate and absorbance as the vertical coordinate, as shown in Fig. 8. The three curves in the figure corresponded to the infra-red absorption spectra of the gaseous products of TBPB thermal decomposition at the peak temperature (160 °C), the end temperature of the exothermic phase (192 °C), and the end temperature of the experiment (500 °C), respectively. Through comparing the three curves, it can be found that at 160 °C, the spectrum showed the most obvious multiple absorption peaks, proving that the decomposition of TBPB was the most intense, and a large number of gaseous products were generated.

      Figure 8. 

      The infrared spectrogram of TBPB decomposition products.

      To determine the main chemical bonds and characteristic groups contained in the gas molecules, the characteristic absorption peaks of the gaseous products at 160 °C were analyzed. For alcohols (O−H bending in-plane: 1,310–1,410 cm−1; O−H stretching: 3,200–3,700 cm−1; C–O stretching: 1,000–1,250 cm−1), carboxylic acids, aldehydes, ketones (C=O stretching: 1,600–1850 cm−1), aromatic acids (C–O stretching: 1,205–1,290 cm−1), mono-substituted benzene ring (C=C skeleton vibration: 1,430–1,650 cm−1; C−H stretching in benzene ring: 3,030–3,125 cm−1; C−H bending in-plane: 1,000–1,250 cm−1; C−H bending out-of-plane: 690–770 cm−1), methyl group (C−H stretching: 29,62 ± 10 cm−1), and carbon dioxide (C=O=O antisymmetric stretching: 2,280–2,390 cm−1; C=O=O bending: 600–700 cm−1)[21,36,37]. It can be seen that the gaseous products of thermal decomposition of TBPB may mainly contain hydroxyl, carbonyl, carboxyl, benzene ring (mono-substituted), methyl, and other groups, and it is presumed that the decomposition products are mainly alkanes, aldehydes or ketones, aromatic acids and alcohols, which specifically may be tert-butanol, benzoic acid, etc. When the reaction proceeded to 192 °C or even higher temperatures, the gaseous products of the reaction became less and less, and only carbon dioxide remained at the end of the experiment.

    • The EPR experimental spectrum of the adducts formed by the free radicals generated during TBPB thermal decomposition captured by the spin-trapping agent DMPO are presented in Fig. 9a. The spectrum shown in the figure was formed by the superposition of the resonance absorption peaks of multiple radical adducts, which can be categorized into four groups of resonance absorption peaks, representing the four different radicals in the experimental sample. The characteristic parameters of these resonance absorption peaks are listed in Table 2. According to the data in the table, for absorption peaks 1 and 2, AH + 2AN ~ 35 Gauss, the giso values of these two absorption peaks were different, so it can be inferred that there were two different kinds of alkoxy radicals. For absorption peak 3, AH + 2AN ~ 50 Gauss, so it can be attributed to the alkyl radical. Absorption peak 4 was a radical produced by the oxide of the capture agent DMPO. As can be seen in Fig. 9a, the intensity of the absorption peak 3 of alkyl radical was very low, proving that few alkyl radicals may be produced in the thermal decomposition of TBPB at 110 °C.

      Figure 9. 

      EPR experimental and simulated spectra of TBPB decomposition radicals captured by DMPO.

      Table 2.  Characterization parameters of the resonance absorption peaks of TBPB decomposition radicals captured by DMPO.

      Absorption
      peaks
      Free radical attribution giso AN
      (Gauss)
      AH
      (Gauss)
      1 Alkoxy 1 2.0081 13.63 8.91
      2 Alkoxy 2 2.0072 13.63 9.12
      3 Alkyl 2.0075 14.06 20.49
      4 The oxide of DMPO 2.0073 13.95
      −: Not applicable.

      According to the above determination of the radical species, the spectrum was simulated and superimposed by combining the characteristic parameters, and compared with the experimental EPR spectrum, as plotted in Fig. 9b. The simulated curve overlapped well with the actual experimental curve, further verifying that two kinds of different alkoxy radicals and a small amount of alkyl radicals were produced when the TBPB thermal decomposition reaction occurred under the condition of 110 °C.

    • Contrasting the spectrum obtained from the EPR experiment of TBPB/TEMPO at 110 °C with the EPR experimental spectrum of pure TBPB without inhibitor, two curves can be acquired as displayed in Fig. 10. Different from the above curve, there were no corresponding alkoxy and alkyl resonance absorption peaks in the EPR spectrum of TBPB/TEMPO, and the obvious triple peaks on the curve represented the peak pattern of TEMPO itself in the solution, meaning that only one radical, TEMPO, remained in the reaction system under that condition. This phenomenon proves that the thermal decomposition of TBPB no longer produces alkoxy and alkyl radicals when TEMPO is added as the inhibitor. Therefore, TEMPO has an obvious inhibiting effect on the radicals produced by the thermal decomposition of the TBPB pure product.

      Figure 10. 

      Comparison of EPR spectra of TBPB before and after adding the inhibitor (TEMPO).

    • To investigate whether TEMPO has an inhibitory effect on the thermal runaway of TBPB, the exothermic behaviors of the thermal decomposition reaction of TBPB and TBPB/TEMPO were contrasted by a series of calorimetric measurements.

      The DSC exothermic curves of TBPB and TBPB/TEMPO at β of 10.0 °C/min are plotted in Fig. 11. The heat flow curve of TBPB/TEMPO was located below, and the heat flow value at Tp was significantly reduced. According to the specific exothermic parameters listed in Table 3, the initial heat release temperature (To) of TBPB/TEMPO increased, while Tp decreased, indicating that the addition of TEMPO shifted the exothermic interval backward and shortened the exothermic range of the reaction. Meanwhile, ΔH of TBPB/TEMPO was also greatly reduced compared with that of TBPB pure product, from 865.87 to 585.00 J/kg. Thus, it can be verified that TEMPO has a significant effect on reducing the thermal hazard of TBPB thermal decomposition reaction.

      Figure 11. 

      Comparison of thermal decomposition behaviors before and after adding the inhibitor (TEMPO).

      Table 3.  Comparison of the TBPB thermal decomposition parameters before and after adding inhibitor (TEMPO).

      Samples To (°C) Tp (°C) ΔH (J/g)
      TBPB 100.93 165.60 865.87
      TBPB/TEMPO 128.65 162.49 585.00

      The results of the thermal decomposition experiments of TBPB and TBPB/TEMPO under adiabatic conditions are presented in detail in Table 4, and the adiabatic temperature rise curves of these two samples are displayed in Fig. 12.

      Table 4.  Comparison of TBPB's ARC experiment results before and after adding inhibitor (TEMPO).

      Samples Toad
      (°C)
      Tfad
      (°C)
      ΔTad
      (°C)
      ΔTad*
      (°C)
      Poad
      (MPa)
      Pfad
      (MPa)
      ΔPad
      (MPa)
      ΔHad
      (J/g)
      TBPB 84.93 230.76 145.83 193.17 0.20 1.81 1.61 1989.65
      TBPB/TEMPO 92.03 204.73 112.70 149.25 0.28 1.46 1.18 1537.29

      Figure 12. 

      Comparison of the adiabatic temperature rises before and after adding the inhibitor (TEMPO).

      From the two curves, it can be seen that TBPB decomposed and released a large amount of heat from 748.3 to 1,016.0 min of the experiment, leading to a rapid temperature increase from the adiabatic decomposition onset temperature (Toad) 84.93 °C to the finishing temperature (Tfad) 230.76 °C, within a short period. Nevertheless, the decomposition reaction of TBPB/TEMPO occurred from 775.6 to 968.2 min, and the temperature accordingly increased from 92.03 to 204.73 °C. It demonstrated that TEMPO can also effectively delay and shorten the thermal decomposition of TBPB under adiabatic conditions.

      According to the data in Table 4, the measured adiabatic temperature rise (ΔTad) of TBPB/TEMPO was 33.13 °C lower than that of TBPB. Since the thermal inertia factor (φ) of the experiment was 1.325 if ignoring the heat transfer of the equipment and external environment, the adiabatic temperature rise after the modification (ΔTad*) of TBPB and TBPB/TEMPO was 193.17 and 149.25 °C respectively under ideal conditions. Through the comparative analysis from Table 4, when TEMPO acted as an inhibitor, ΔTad* of TBPB can be effectively reduced by 43.93 °C, and the heat release of adiabatic decomposition (ΔHad) was significantly reduced by nearly a quarter from the original 1,989.65 to 1,537.29 J/g. During the adiabatic decomposition of TBPB, a large number of gaseous products were also generated, and the gas would be affected by a sudden increase in temperature, so the pressure in the sample ball would also increase significantly. Compared with the TBPB pure product, the adiabatic pressure rise (ΔPad) of TBPB/TEMPO was reduced by 0.43 MPa. In summary, the addition of TEMPO can not only moderate the heat release of TBPB during the adiabatic decomposition, but also reduce the pressure release, thus greatly reducing the hazardous and destructive effects of the reaction.

    • A series of calorimetric tests and free radical trapping experiments were conducted to study the hazard of the thermal decomposition reaction of TBPB, the products, and the types of free radicals generated during the reaction process. At the same time, the inhibition ability of the free radical inhibitor TEMPO on the reaction radicals of TBPB and the inhibition effect of thermal runaway were measured and evaluated, and the major findings are as follows:

      Under normal conditions, TBPB thermal decomposition happened in the range of 100–210 °C and gave out a large amount of heat, the heat release could reach 924.59 J/g. The average activation energy of the whole reaction was 97.55 kJ/mol, and the energy required for the reaction decreased as the decomposition progressed. The thermal hazards of TBPB are enormous, and the consequences of a thermal runaway can be horrific.

      During the TBPB thermal decomposition reaction process, two kinds of alkoxy radicals and a kind of alkyl radical would be generated. The gaseous products of the thermal decomposition of TBPB tested and analyzed by infrared technology were mainly alkanes, aldehydes or ketones, aromatic acids, and alcohols.

      Under both conventional and adiabatic conditions, TEMPO had an obvious inhibition effect on the thermal decomposition of TBPB. The addition of trace amounts of TEMPO to TBPB could effectively inhibit alkoxy radicals and alkyl radicals generated during the decomposition process. At the same time, the heat release and the adiabatic temperature rise of the thermal decomposition reaction could be greatly reduced. TEMPO played a positive role in reducing the thermal danger and the risk of thermal runaway of TBPB.

      • This study was supported by the National Natural Science Foundation of China (No. 52274209, 21927815, 52334006, 51834007), Jiangsu Province '333' project(BRA2020001, Jiangsu Qing Lan Project, and Jiangsu Association for Science and Technology Youth Talent Support Program.

      • The authors confirm contribution to the paper as follows: study conception and design: Zhang D, Jiang J, Ni L; data collection: Zhang D, Li Z; analysis and interpretation of results: Li Z, Chen Z; draft manuscript preparation: Zhang D, Li Z. All authors reviewed the results and approved the final version of the manuscript.

      • All data generated or analyzed during this study are included in this published article.

      • The authors declare that they have no conflict of interest. Juncheng Jiang is the Editorial Board member of Emergency Management Science and Technology who was blinded from reviewing or making decisions on the manuscript. The article was subject to the journal's standard procedures, with peer-review handled independently of this Editorial Board member and the research groups.

      • Copyright: © 2025 by the author(s). Published by Maximum Academic Press on behalf of Nanjing Tech University. This article is an open access article distributed under Creative Commons Attribution License (CC BY 4.0), visit https://creativecommons.org/licenses/by/4.0/.
    Figure (12)  Table (4) References (37)
  • About this article
    Cite this article
    Zhang D, Li Z, Jiang J, Ni L, Chen Z. 2025. Thermal hazard assessment and free radical inhibition of decomposition of tert-butyl perbenzoate. Emergency Management Science and Technology 5: e001 doi: 10.48130/emst-0024-0029
    Zhang D, Li Z, Jiang J, Ni L, Chen Z. 2025. Thermal hazard assessment and free radical inhibition of decomposition of tert-butyl perbenzoate. Emergency Management Science and Technology 5: e001 doi: 10.48130/emst-0024-0029

Catalog

  • About this article

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return