Thermal Energy and Luminosity Characterization of an Advanced Ignition System Using a Non-Intrusive Methodology in an Optically Accessible Calorimeter

: To restrain the environmental impact of modern SI engines, igniters must guarantee stable combustions with low cycle-to-cycle variability in extreme operating conditions (high EGR, ultra-lean), via high energy release in the combustion chamber. The direct measurement of this energy is not trivial and requires a controlled environment. Luminosity detection is a non-intrusive diagnostic technique to indirectly measure the thermal energy released by the discharge on optically accessible apparatus. This work compares energy and luminosity produced by a plasma igniter in a constant volume vessel at realistic working conditions (ignition at 8 bar and air as a medium). A calibration factor can be deﬁned to describe the thermal energy behavior as a function of the discharge luminosity and to give an assessment of such approach for its use in optically accessible engine. This study shows that thermal energy and luminosity are inﬂuenced by the gas type and related by a linear relationship for both air and nitrogen. The presence of oxygen resulted in discharges with reduced energy delivery to the medium and a lower discharge luminosity compared to nitrogen. This work outcome could improve the use of a non-intrusive methodology, based on luminosity detection, to characterize the igniter performance, exploitable for 3D-CFD.


Introduction
The role of internal combustion engines (ICE) has been predominant for ground transportation in recent decades.In the foreseeable future, ICE will still hold a key role in transportation, industry and power generation, justifying the manufacturers endeavor to design more efficient and cleaner engines [1,2].The path toward a more sustainable transport sector determined the enforcement of ever-stringent pollutant emissions and greenhouse gas regulations in order to limit the impact of fossil fuels on a local and global scale.Such targets require Spark-Ignition (SI) engines to adopt combustion strategies involving increased boost level, water injection [3,4], lean burn and/or high exhaust gas recirculation (EGR) diluted mixture [5,6].In these challenging operating conditions, the reduced mixture reactivity demands high ignition energies to be delivered by the igniter to guarantee an effective combustion initiation [7][8][9].Hence, the determination of the amount of energy deposited into the medium is of pivotal importance to characterize the ignition system capability of providing a robust and stable ignition process [10,11].In fact, the thermal energy released by the discharge accentuates the chemical kinetics of the fuel oxidation chain reactions by increasing the local gas temperature, thus mightily influencing flame kernel formation [12,13].
Energies 2023, 16, 520 2 of 22 Generally speaking, the igniter thermal energy characterization requires the adoption of dedicated controlled environments, such as the constant volume vessel [14 -16].In this context, the discharge produces a pressure rise inside the chamber, which can be considered proportional to the released energy, in the approximation of adiabaticity, and measured by employing extremely sensible pressure sensors [14].This approach cannot be exploited in real engine applications because of the inherent dynamic pressure variation occurring during the compression stroke.Therefore, the small pressure oscillation caused by the discharge event could not be appropriately quantified.Moreover, the approximation of adiabaticity is hardly acceptable and the imperfect combustion chamber sealing, resulting in the leakage of a portion of the compressed charge, renders the determination of the working gas volume quite challenging.
Ricci et al. [17], in their previous work, studied the correlation among the luminous emission and the thermal energy deposited into nitrogen by a plasma-assisted igniter (PAI), in an optically accessible pressure-based calorimeter.A linear correlation between the two main outputs of the discharge has been found.The main purpose was to furnish an innovative optical method to estimate the thermal energy released by the igniter through the estimation of the discharge luminosity.
In ICE research field, luminosity detection is a well-known approach to providing quantitative information about combustion inception and flame kernel formation where the indicated analysis encounters its limit [18][19][20].Optical access engines are widely adopted to study the flame propagation, the combustion stability, the heat release rate and as diagnostic tool for the combustion process and pollutant emissions [21][22][23].Predominantly, these analyses are performed exploiting high-speed cameras, whose output frames permit not only luminosity emission to be assessed but also 2D-spatial information to be maintained [24,25].Another technique to collect discharge or flame luminosity consists of placing an optical sensor-a photodetector-inside the combustion chamber, to provide luminous emission magnitude of certain area but losing the spatial information [26,27].
Within this context, the work herein wants to mirror the methodology in [17] in order to demonstrate if a direct relationship between luminosity and the amount of deposited thermal energy could be extended from a quasi-inert atmosphere (nitrogen) to a real oxidant one (air), and how the latter could affect the igniter discharge.As a test medium, synthetic air has been employed with the aim of investigating near-like engine oxidant at engine relevant pressure.A PAI was chosen to carry out the experimental campaign, and specifically, a barrier discharge igniter (BDI) prototype has been tested in an optical access calorimeter by simultaneously acquiring the pressure variation inside the chamber and the event luminosity due to the discharge.
Differently from traditional igniters, PAIs are innovative ignition systems based on the generation of non-equilibrium plasma, in the form of streamers (ionization waves), and characterized by the absence of a breakdown phase as well as a relatively larger ignition volume [28,29].On one hand, the physical features of the discharge promote ignition not only through thermal effect but also by chemical kinetic effect, i.e., the production of active radicals and excited species [30].On the other hand, a wider discharge volume renders such igniters well-suited to collect luminosity data and then to correlate those with thermal energy measurements.Several concepts have been studied to date: nanosecond pulse discharge [15], microwave assisted [31], radio-frequency corona ignition [32][33][34] and barrier discharge ignition [35,36].
Among those, the BDI proved its capability of promoting flame kernel formation and enhance combustion stability, which are attributes of major interest for combustion in lean, ultra-lean or highly EGR diluted mixture conditions [37,38].
The selection of a BDI relates to its discharge modulation capability, in terms of peak electrode voltage and time duration, and its dielectric layer (barrier).The first permits several luminosity and thermal energy levels to be investigated, while the second prevents streamer-to-arc transition [38].In addition to that, the igniter choice was driven from the Energies 2023, 16, 520 3 of 22 fact that numerical results suggested how the radical production effect is less predominant than the thermal deposition one [39].
Thermal effect is the dominant component for a traditional igniter such as a spark plug, since they are associated with the production of thermal or equilibrium plasma (because the plasma is in thermal equilibrium with the free electrons and therefore it has their high temperature).The latter is not capable of producing a substantial ionization (i.e., radicals production), that proved to be a pivotal factor in enhancing flame kernel formation in hard-to-ignite mixtures when non-equilibrium plasma igniters were adopted.However, the thermal effect is still produced by the PAIs and still plays a role in the overall ignition process.On the other hand, in the absence of spectral analysis, the impact of each ignition mechanism is not distinguishable and cannot be quantified separately: the overall igniter ability to start the fuel mixture combustion allows to determine the main ignition mechanism only, since it heavily affects the flame onset speed (fast speed for the low temperature non-equilibrium plasma related effects, slow speed for the thermal effect).Additionally, the overall luminosity is a sum of contributions, due both to the thermal effect and to the electronic transitions, distinctive of the PAI physics, in the (not only) visible light spectrum.
The luminosity and the ability to ignite a fuel mixture are then functions of these effects (thermal effect as well as ionization effect) and therefore can be related each other.Even if this relationship is precisely unknown, is fixed once the particular igniter kind is chosen, allowing for the luminosity calibration.
Consequently, the present work has two major intents: the first one is to fully analyze the luminosity-thermal energy correlation in the presence of air, an in-cylinder engine mixture component (fresh air, fuel vapor, exhaust residual) present in optical and metal engines, and assess the impact of oxygen on the discharge behavior, since it is by far the most active species in the combustion chamber; the second is to evaluate if luminosity can be exploited as a non-intrusive diagnostic parameter, in applications like optical access engine, to quantify the thermal energy released into the mixture.Moreover, since the oxygen presence affects dramatically the discharge morphology, in terms of shape and branching, the present paper shows a novel and quantitative approach to characterize the physical development of the streamers on the igniter surface to highlight and evaluate the oxygen contribution.Finally, a comparison with pure N2 [17] is also proposed, not only in terms of luminous emission, but also considering the main discharge features.

Igniter
The dielectric barrier discharge, studied in this work (Figure 1b), is produced by the so-called barrier discharge igniter.The device belongs to the advanced corona ignition system (acis) family [37] and it was provided by Federal Mogul Powertrain-a Tenneco group company.
A dedicated electronic system (ACIS Box) powers the igniter assembly with an input radiofrequency of about 1.04 MHz, corresponding to the resonance frequency of the equivalent RLC circuit [38] of Figure 1b.
The main control parameters are the driving voltage V d and the activation time t on [37], which are directly correlated with peak electrode voltage (V e ) and discharge duration [16], respectively.Once V d is set, the electronic system magnifies the voltage up to a proportional value (supplied voltage, V s ) and provides it to the coil.This latter amplifies the voltage to the firing end up to V e in order to produce the discharge.
According to the equivalent model (Figure 1a), the electrode voltage V e is proportional to the driving voltage V d .Since measuring V e is not allowed without causing the dielectric breakdown [35,36], the driving voltage V d is used as reference.
For an established gas mixture, density and pressure, the discharge arises in a precise voltage range, from corona inception voltage (V d CI ) to the maximum value allowed by the manufacturer (V d MAX ).Above the V d CI , the applied electric field surpasses a criti- cal threshold, promoting the discharge onset.Then, streamers (ionization waves) start propagating annularly from the zero-potential BDI thread on top of its dielectric globe (Figure 1c).The hemispheric alumina insulation covers the BDI inner electrode and hinders streamer-to-arc transition phenomenon near breakdown voltage values [36,40] in order to ensure the dielectric barrier discharge [41,42].A dedicated electronic system (ACIS Box) powers the igniter assembly with an input radiofrequency of about 1.04 MHz, corresponding to the resonance frequency of the equivalent RLC circuit [38] of Figure 1b.
The main control parameters are the driving voltage Vd and the activation time ton [37], which are directly correlated with peak electrode voltage (Ve) and discharge duration [16], respectively.Once Vd is set, the electronic system magnifies the voltage up to a proportional value (supplied voltage, Vs) and provides it to the coil.This latter amplifies the voltage to the firing end up to Ve in order to produce the discharge.
According to the equivalent model (Figure 1a), the electrode voltage Ve is proportional to the driving voltage Vd.Since measuring Ve is not allowed without causing the dielectric breakdown [35,36], the driving voltage Vd is used as reference.
For an established gas mixture, density and pressure, the discharge arises in a precise voltage range, from corona inception voltage (Vd CI ) to the maximum value allowed by the manufacturer (Vd MAX ).Above the Vd CI , the applied electric field surpasses a critical threshold, promoting the discharge onset.Then, streamers (ionization waves) start propagating annularly from the zero-potential BDI thread on top of its dielectric globe (Figure 1c).The hemispheric alumina insulation covers the BDI inner electrode and hinders streamer-toarc transition phenomenon near breakdown voltage values [36,40] in order to ensure the dielectric barrier discharge [41,42].

Experimental Setup
Tests were performed making use of an experimental set-up where a constant volume optical calorimeter was installed (Figure 2).The latter featured a volume of ν ch = 65.2 cm 3 and it was made of plexiglass to ensure low thermal conductivity (0.187 W•m −1 •K −1 ).The vessel was filled with synthetic air at room temperature (293 ± 1 • K), pressure regulation and flow exchange were managed by a pneumatic system equipped with a precision pressure reducer Festo LRP and one-way flow valves.It must be noted that a gas replacement after each discharge event (i.e., similarly to what occurs inside an internal combustion engine) was not feasible due to the time required for an effective medium replacement.Each discharge event was triggered by a TTL signal produced by an arbitrary wave generator (HP 33120A), at 10 Hz.With such frequency, the pressure level before a subsequent discharge is steady, so that each discharge is not affected by the previous one.Moreover, 10 Hz is coherent with the discharge frequency of the BDI at 1000-2000 rpm in a 4-stroke engine, a typical speed range in which the optical engine of our research group operated [23].The signal coming from an arbitrary wave generator was also utilized to trigger the acquisition of a Phantom V710 high-speed camera, placed in front of the calorimeter, which records Energies 2023, 16, 520 5 of 22 the igniter discharge natural luminosity.The voltage and the current at the igniter primary circuit were gathered by a Teledyne LeCroy PP020 passive probe measure and Teledyne LeCroy CP030 current probe (accuracy ±1.5% full scale), respectively.The pressure fluctuations due to the discharge event were recorded by a Kistler Type 7261 piezoelectric pressure transducer with 2200 pC•bar −1 of sensitivity and ≈10 −5 bar of resolution).The latter signal was amplified and converted by a Kistler Type 5011 charge amplifier into a proportional voltage signal.All the aforementioned signals were acquired by a fast oscilloscope Lecroy Wavesurfer 3000 with a sampling frequency of 10 MHz, about ten times larger than the BDI working frequency.For each test point, 25 consecutive events were recorded and stored.This experimental apparatus allowed us to generate BDI discharges and to match the subsequent pressure rise with the streamers evolution images.
replacement.Each discharge event was triggered by a TTL signal produced by an arbitrary wave generator (HP 33120A), at 10 Hz.With such frequency, the pressure level before a subsequent discharge is steady, so that each discharge is not affected by the previous one.Moreover, 10 Hz is coherent with the discharge frequency of the BDI at 1000-2000 rpm in a 4-stroke engine, a typical speed range in which the optical engine of our research group operated [23].The signal coming from an arbitrary wave generator was also utilized to trigger the acquisition of a Phantom V710 high-speed camera, placed in front of the calorimeter, which records the igniter discharge natural luminosity.The voltage and the current at the igniter primary circuit were gathered by a Teledyne LeCroy PP020 passive probe measure and Teledyne LeCroy CP030 current probe (accuracy ±1.5% full scale), respectively.The pressure fluctuations due to the discharge event were recorded by a Kistler Type 7261 piezoelectric pressure transducer with 2200 pC•bar −1 of sensitivity and ≈10 −5 bar of resolution).The latter signal was amplified and converted by a Kistler Type 5011 charge amplifier into a proportional voltage signal.All the aforementioned signals were acquired by a fast oscilloscope Lecroy Wavesurfer 3000 with a sampling frequency of 10 MHz, about ten times larger than the BDI working frequency.For each test point, 25 consecutive events were recorded and stored.This experimental apparatus allowed us to generate BDI discharges and to match the subsequent pressure rise with the streamers evolution images.In this work, synthetic air Linde 5.0 with a composition of 80% N 2 and 20% O 2 in volume was used [43].In the visible range, the discharge emission was mainly included in the 390-435 nm range for N 2 [37,38] and to the 525-845 nm range for oxygen species [44].
The high-speed camera (Table 1) was equipped with a 6-diopter close-up lens placed downstream of a 50 mm f/1.4 lens.Each frame comprised 256 × 256 pixels and sampled at the maximum framerate allowed at the chosen image resolution (herein 79,000 fps, around one every 12.66 ms).The minimum allowed distance between camera and igniter was used in order to improve the spatial resolution by reducing the minimum focal distance (54 µm/pixel of spatial resolution was obtained).
The experimental optical apparatus (imaging apparatus [45] and plexiglass calorimeter) is sensitive in all the above-mentioned ranges.Moreover, no special filters were used in this test campaign.Once the pressure signal (Figure 3a) was acquired, the deposited thermal energy (E r ) was computed starting from the internal pressure variation (∆P = P max − P min ) due to the discharge event, whose estimation is performed by in-house post-processing algorithms in MATLAB environment.[44].
The high-speed camera (Table 1) was equipped with a 6-diopter close-up lens placed downstream of a 50 mm f/1.4 lens.Each frame comprised 256 × 256 pixels and sampled at the maximum framerate allowed at the chosen image resolution (herein 79,000 fps, around one every 12.66 ms).The minimum allowed distance between camera and igniter was used in order to improve the spatial resolution by reducing the minimum focal distance (54 μm/pixel of spatial resolution was obtained).
The experimental optical apparatus (imaging apparatus [45] and plexiglass calorimeter) is sensitive in all the above-mentioned ranges.Moreover, no special filters were used in this test campaign.Once the pressure signal (Figure 3a) was acquired, the deposited thermal energy (Er) was computed starting from the internal pressure variation (ΔP =  −  ) due to the discharge event, whose estimation is performed by in-house post-processing algorithms in MATLAB environment.To correctly determine it, the recorded pressure signal (Figure 3a, grey line) was filtered by a 2 kHz low-pass filter [10] (Figure 3b, red line).
The determination of the cut-off frequency was performed by considering the FFT of the pressure signal with the aim of removing noise contribution [10].A value of 2 kHz was found to be adequate for this purpose (Figure 3b).After the low-pass filtering, the To correctly determine it, the recorded pressure signal (Figure 3a, grey line) was filtered by a 2 kHz low-pass filter [10] (Figure 3b, red line).
The determination of the cut-off frequency was performed by considering the FFT of the pressure signal with the aim of removing noise contribution [10].A value of 2 kHz was found to be adequate for this purpose (Figure 3b).After the low-pass filtering, the algorithm identifies the steady values before (P min ) and after (P max ) the discharge, to compute ∆P.The steady value determination would be much more unstable if performed on a raw signal, especially for low energy discharges (low voltage and/or low duration), where the signal-to-noise ratio is low.
Assuming the hypothesis of adiabaticity of the chamber (effective for some milliseconds), the first law of thermodynamics can be used to compute the deposited thermal energy (1): where ν ch is the chamber volume and γ is the specific heat ratio of the gas.To compute the average value E r and the standard deviation, 25 measurements were carried out for each tested point.

Luminosity Analysis
Depending on the discharge activation time of the device, the recorded events were composed by a certain number of frames (S), each of which comprised 256 × 256 pixels at 8 bits.In a similar manner to the thermal energy, the estimation of the discharge luminosity is performed by in-house MATLAB scripts.The luminosity level of each pixel was assumed proportional to the grey level intensity acquired by the camera [46] (Figure 4a).
sumed proportional to the grey level intensity acquired by the camera [46] (Figure 4a).
Since pixel saturation could influence the energy-luminosity correlation, such circumstance must be limited to properly interrelate the amount of thermal energy released from the discharge to its brightness.Because of this, an adequate selection of the camera diaphragm position was found to obtain a trade-off between restricting the number of saturating pixels (in case of high energy-brightness discharge), by reducing the luminosity absorption, and avoiding the loss of resolution for the low energy-brightness discharge events.Before image processing, the user circumscribes the discharge event by selecting, on the file acquired by the camera PC, the frames to be post-processed.In fact, the 20 frames before the start of discharge were needed to characterize the background, while 10 frames were required after the end of the discharge to make sure that the discharge was effectively over (Figure 4b).All the steps of the employed procedure for the luminosity estimation are described below: • Background noise characterization and removing.The average grey level value of each pixel in the 20 initial frames with no discharge was computed (Figure 5a).A map of the average noise in the 256 × 256 pixel matrix can be obtained by computing the mean raw image of these 20 frames.
The gray levels distribution of such image was found to be predominantly in the range 60-70 (Figure 5b).Once determined, noise was removed, for each frame of the specific series, by subtracting, pixel-by-pixel, the average background map (Figure 5c).Since pixel saturation could influence the energy-luminosity correlation, such circumstance must be limited to properly interrelate the amount of thermal energy released from the discharge to its brightness.Because of this, an adequate selection of the camera diaphragm position was found to obtain a trade-off between restricting the number of saturating pixels (in case of high energy-brightness discharge), by reducing the luminosity absorption, and avoiding the loss of resolution for the low energy-brightness discharge events.Before image processing, the user circumscribes the discharge event by selecting, on the file acquired by the camera PC, the frames to be post-processed.In fact, the 20 frames before the start of discharge were needed to characterize the background, while 10 frames were required after the end of the discharge to make sure that the discharge was effectively over (Figure 4b).
All the steps of the employed procedure for the luminosity estimation are described below:

•
Background noise characterization and removing.The average grey level value of each pixel in the 20 initial frames with no discharge was computed (Figure 5a).A map of the average noise in the 256 × 256 pixel matrix can be obtained by computing the mean raw image of these 20 frames.The gray levels distribution of such image was found to be predominantly in the range 60-70 (Figure 5b).Once determined, noise was removed, for each frame of the specific series, by subtracting, pixel-by-pixel, the average background map (Figure 5c).

•
Filtering: A 2D Gaussian filter with a filter order σ = 1 was applied to reduce the residual noise.A sensitivity campaign was preliminarily realized on the filter type (median, gaussian or no filter at all) and filter order to determine the best setting.Gaussian filter with σ value equal to 1 belongs to a typical range that ensures effective Gaussian filtering [47] and, for this work, it was found to be a good compromise between noise cancellation and boundary conservation.This step allows artifacts and occasional spare pixels to be excluded, which could wrongly influence the final results.

•
Thresholding: For all the recorded series, preliminary tests were carried out to estimate the fixed threshold (specifically found to be 10 out of 255) to determine the differences between background and the luminosity emission related to the corona streamer.Such threshold was chosen, taking into consideration the image grey level distribution (Figure 5c), in order to account for the effect of light diffusion, in regions around the streamers' body, on the computed discharge luminosity level.The value of 10 was selected as a suitable threshold to obtain a compromise in terms of luminosity over-or under-estimation.
• Average grey level l g estimation.Only for the frames in which the discharge is present, the scalar parameter "Average grey level" l g (2) is a measure of the luminous emission of the streamers in a period corresponding to the sensor exposure time, so it basically represents the average brightness of the frame.
where N is the number of pixels (256 × 256) and l j the luminosity of the jth pixel.• Filtering: A 2D Gaussian filter with a filter order σ = 1 was applied to reduce the residual noise.A sensitivity campaign was preliminarily realized on the filter type (median, gaussian or no filter at all) and filter order to determine the best setting.Gaussian filter with σ value equal to 1 belongs to a typical range that ensures effective Gaussian filtering [47] and, for this work, it was found to be a good compromise between noise cancellation and boundary conservation.This step allows artifacts and occasional spare pixels to be excluded, which could wrongly influence the final results.

•
Thresholding: For all the recorded series, preliminary tests were carried out to estimate the fixed threshold (specifically found to be 10 out of 255) to determine the differences between background and the luminosity emission related to the corona streamer.Such threshold was chosen, taking into consideration the image grey level distribution (Figure 5c), in order to account for the effect of light diffusion, in regions around the streamers' body, on the computed discharge luminosity level.The value of 10 was selected as a suitable threshold to obtain a compromise in terms of luminosity over-or under-estimation.

•
Average grey level lg estimation.Only for the frames in which the discharge is present, the scalar parameter "Average grey level" lg (2) is a measure of the luminous emission of the streamers in a period corresponding to the sensor exposure time, so it basically represents the average brightness of the frame.
where N is the number of pixels (256 × 256) and lj the luminosity of the j th pixel.Moreover, in order to estimate the surface of the streamers propagating on the BDI cupola, each frame was cropped around the igniter annular base.The aim was to not consider, as active streamer surface, the luminosity emission outside the cupola perimeter.A binarization processes is needed to estimate the streamer surface.The binarized area was computed by counting the white pixels for each black-and-white frame.This procedure allowed us to characterize the streamers branching around the BDI firing end with the purpose to investigate, from a different point of view, the discharge behavior when using a different medium.Figure 6a reports the procedure to obtain, at specific operating point, the average grey level lg of each frame (see up to Filtering step) and the surface streamers (Cropping + Binarization step).Moreover, in order to estimate the surface of the streamers propagating on the BDI cupola, each frame was cropped around the igniter annular base.The aim was to not consider, as active streamer surface, the luminosity emission outside the cupola perimeter.A binarization processes is needed to estimate the streamer surface.The binarized area was computed by counting the white pixels for each black-and-white frame.This procedure allowed us to characterize the streamers branching around the BDI firing end with the purpose to investigate, from a different point of view, the discharge behavior when using a different medium.Figure 6 reports the procedure to obtain, at specific operating point, the average grey level l g of each frame (see up to Filtering step) and the surface streamers (Cropping + Binarization step).The sequence of all consecutive lg values (black lines, Figure 7), characterizing a specific event, allow a trend to be obtained-Lg (3). = ( ,  , . . .,  ) (3) The sequence of all consecutive l g values (black lines, Figure 7), characterizing a specific event, allow a trend to be obtained-L g (3).

Test Campaign
The experimental campaign was conducted in synthetic air (80% N2-20% O2 by volume) at an internal pressure chamber of 8 bar absolute and at room temperature.The 8 bar pressure value was selected because it could be representative of charge conditions, within the combustion chamber of optical access and metal engines, when the discharge event occurred [38,39].
Furthermore, such value refers to a widely investigated pressure range which can be found in several calorimeter studies [14 -16].The test campaign was built by varying each operating parameter (Vd, ton), one at a time.The discharge duration plays an important role in ignition strategy, since it could act as a reducer for the high cycle-to-cycle variability, characteristic of the ultra-lean blends or high EGR levels [28,39].The driving voltage Vd is instead responsible for the corona development, which depends on the working conditions, namely, medium and internal pch [16,46].In the first part of the test campaign, a sweep of Vd maintaining a fixed ton was performed: the BDI activation time was set to 300 μs, while the driving voltage was varied from the inception voltage to the maximum allowed by the manufacturer (60 V), with 1 V steps.Concerning the activation time, a value of 300 μs resulted to be an effective choice to guarantee mixture ignition in the internal combustion engine [20,28].Figure 8   The average trend L g (4) was then obtained by averaging, frame-by-frame, the l g values of the 25 consecutive discharges, which represents the average luminosity trend for a specific test point (red line, Figure 7).
During the activation time of the igniter, the luminous emission depends not only on the steady state value but also on the transients (indicated by blue dashed lines in Figure 7), because of the discharge management carried out by the electronics.In order to account for the rising and falling luminosity transient and to characterize the average luminosity behavior of the discharge, the integral luminosity I Lg is defined.Such a parameter represents the frame-by-frame sum of L g levels (l 1 g , l 2 g , . . ., l S g ) within frame interval where the discharge is present (red transparency area, Figure 7).As a matter of fact, I Lg is directly proportional to the area under the L g curve, and it is used to compared different tested points.

Test Campaign
The experimental campaign was conducted in synthetic air (80% N 2 -20% O 2 by volume) at an internal pressure chamber of 8 bar absolute and at room temperature.The 8 bar pressure value was selected because it could be representative of charge conditions, within the combustion chamber of optical access and metal engines, when the discharge event occurred [38,39].Furthermore, such value refers to a widely investigated pressure range which can be found in several calorimeter studies [14 -16].The test campaign was built by varying each operating parameter (V d , t on ), one at a time.The discharge duration plays an important role in ignition strategy, since it could act as a reducer for the high cycle-to-cycle variability, characteristic of the ultra-lean blends or high EGR levels [28,39].The driving voltage V d is instead responsible for the corona development, which depends on the working conditions, namely, medium and internal p ch [16,46].In the first part of the test campaign, a sweep of V d maintaining a fixed t on was performed: the BDI activation time was set to 300 µs, while the driving voltage was varied from the inception voltage to the maximum allowed by the manufacturer (60 V), with 1 V steps.Concerning the activation time, a value of 300 µs resulted to be an effective choice to guarantee mixture ignition in the internal combustion engine [20,28].Figure 8  In the second part of the test campaign, a sweep of ton from 50 μs to 1500 μs was conducted.The minimum value was chosen as representative of a very small discharge duration, as found with other kinds of ACIS igniters [33].The maximum value selection prevents coil overheating with a fair safety margin.A time-step extension of 50 μs was adopted from minimum ton up to 300 μs to magnify what happens when a short-lasting discharge is extended, while a larger step of 300 μs was applied for the remaining tests.Three diverse driving voltages were chosen (35 V, 45 V, 55 V) to properly operate the BDI in a voltage range where an effective discharge was produced.
The test campaign performed in this work is summarized in Table 2.

Results and Discussions
The purpose of this experimental investigation was to examine, in a qualitative and quantitative manner, the behavior of a BDI igniter in air, as the control parameters vary.The luminosity and the energy release measurements, as well as the discharge morphology analysis, are herein proposed in order to highlight the effect of oxygen on the igniter performance.Once the igniter was characterized in air with varying operating conditions, the results were compared with those obtained using nitrogen as a reference medium in a previous work of our research group [17].In fact, differently from the weak interaction In the second part of the test campaign, a sweep of t on from 50 µs to 1500 µs was conducted.The minimum value was chosen as representative of a very small discharge duration, as found with other kinds of ACIS igniters [33].The maximum value selection prevents coil overheating with a fair safety margin.A time-step extension of 50 µs was adopted from minimum t on up to 300 µs to magnify what happens when a short-lasting discharge is extended, while a larger step of 300 µs was applied for the remaining tests.Three diverse driving voltages were chosen (35 V, 45 V, 55 V) to properly operate the BDI in a voltage range where an effective discharge was produced.
The test campaign performed in this work is summarized in Table 2.

Results and Discussions
The purpose of this experimental investigation was to examine, in a qualitative and quantitative manner, the behavior of a BDI igniter in air, as the control parameters vary.
The luminosity and the energy release measurements, as well as the discharge morphology analysis, are herein proposed in order to highlight the effect of oxygen on the igniter performance.Once the igniter was characterized in air with varying operating conditions, the results were compared with those obtained using nitrogen as a reference medium in a previous work of our research group [17].In fact, differently from the weak interaction properties of nitrogen [15], the presence of oxygen and the corresponding excited species, produced after discharge onset, could have a substantial influence not only on the global discharge behavior but also on the igniter electrode performance and durability [16].

Driving Voltage Sweep
The driving voltage V d is directly responsible for the entity of the streamer production and, consequently, for the thermal energy released into the surrounding medium as well as of the luminosity produced by the discharge: Figure 9 reports quantitative trends of energy released E r (a) and integral luminosity I Lg (b) for both air (black markers) and nitrogen (blue markers); whereas Figure 10 shows images of streamers propagating over the BDI cupola in presence of air (left column) and nitrogen (right column).The corresponding applied V d and the resulting E r and I Lg are also therein reported, all of them increasing from the top to the bottom of the figure.Each row reports a pair of images with approximately the same amount of released thermal energy.Using the BDI maximum driving voltage (established by the manufacturer, V d = 60 V), an opening step was to identify the inception voltage V CI in air at 8 bar absolute pressure, in order to determine the igniter operating range [40].All the events are considered capable of inducing a pressure variation inside the chamber if the uncertainty intervals of p max ch and p min ch do not overlap.Starting from this consideration, the corona inception voltage was identified as 32 V.It is worth mentioning that the high-speed camera records a non-zero luminous response even at lower voltages-up to 30 V (I Lg = 10).This occurrence was similarly reported with nitrogen (V d CI = 27 V), suggesting again that either the camera is more sensitive than the piezoelectric pressure sensor or different energy transfer patterns are present near inception conditions.In addition to that, such evidence confirms that a non-zero thermal energy release occurs up to 24 V (I Lg = 12) for nitrogen.However, air required a higher V d CI (Figure 10a) for streamers generations with respect to the one found in N 2 (Figure 10b).This behavior could be a result of the quenching effect of molecular oxygen on other excited species and of the ability of oxygen to easily capture the first emitted electrons, thus hindering corona inception at a low driving voltage [48].Regardless of the medium, in the whole analyzed range, the released energy and the luminosity raise as V d increases.Moving toward increasing V d values, a linear relationship (Figure 9) among the analyzed quantities, I Lg (a), E r (b), and V d can be recognized (see the linear data fits), with the exclusion of the low voltage points (filled markers).For the latter V d values, i.e. close to zero E r , the luminosity diverges from the linear trend because of the strong reduction in streamers number (average behavior).Close to V d CI , the energy released and the luminosity produced by the small streamers proportionally decrease as V d decreases, but, under a certain threshold, the single streamers vanish.Such a threshold, variable for each streamer, depends on the electrode surface local roughness that modulates the reduced electric field intensity, responsible, in turn, for the streamer generation and propagation [49,50].For that reason, the linear regressions (red dotted lines) in Figure 9 exclude those operating points presenting E r below 1 mJ, because of the inherent stochasticity of the phenomena.
Considering this, V d , energy data show that in nitrogen (Figure 10e) the thermal energy released to the medium is higher than the case of synthetic air (Figure 10f).Moving toward higher V d , the differences in released thermal energy ∆E r between the two gases become more pronounced.For instance, at V d = 34 V (first common not-filled point in Figure 9) a difference ∆E r of 4.3 mJ (nearly 40% with respect to air) is present, whereas at V d = 60 V (extreme operating point tested) an increment of ∆E r at about 50% was found, thus leading to a difference ∆E r = 6.5 mJ.The lower energy recorded in air could be justified by the fact that the presence of oxygen generates different energy transfer pathways to the medium matter, thus hampering the thermal energy deposition.As a matter of fact, the pulse input energy breakdown changes when oxygen is dispersed into the medium, i.e., energy is less dissipated as heat (thermalization process) but employed in production of plasma chemical species and stored within atoms of O and ozone [51].
Similarly, at the same V d , the luminous emission in nitrogen (Figure 10f) is stronger than in air (Figure 10c), and the higher the V d the larger the I Lg difference between the measurements (Figure 9a).The angular coefficient m of the linear regressions certifies this aspect by showing values equal to 54.97 V −1 in N 2 and about half less in air (m air = 28.06V −1 ).The difference in I Lg among the two gases is even more than the ones found for E r , for which m N2 = 0.58 mJ V −1 and m air = 0.49 mJ V −1 .This evidence is coherent with the results found by Ono et al. [52], in which images of streamers, propagating in mixtures with rising O 2 concentration, required a progressive increase in the luminous intensity amplification factor.Streamer's brightness levels could be related to the quenching effect of the excited species previously produced, which fall to a lower excitation or ground level emitting photons [53].As a consequence, this aspect could determine the presence or the absence of atomic or molecular lines in the emitted spectrum in the visible region of our camera [54,55].This reasoning could possibly justify not only the different luminosity level at the same driving voltage, but also synthetic air/N 2 cases characterized by same energy levels but different values of luminosity (for instance, see Figure 10c,d).From the results of Figure 9, a direct correlation between E r and I Lg can be inferred.Indeed, a linear relationship among thermal energy and discharge luminosity can be estimated for synthetic air (Figure 11a).As expected from Figure 9, the linearity is lost for low values of I Lg (a) and of E r , due to the different instruments' sensitivity and energy transfer modes.By excluding those points, the released thermal energy could be expressed as follows (Equation ( 5)): where m D equals to 0.011 mJ for N 2 and to 0.018 mJ for air (Figure 11b).Therefore, by determining the proper calibration factors, discharge luminosity measurement could be a viable means to compute the released thermal energy into the mixture.It is worth noticing that the lack of linearity in the low energy-luminosity region (points excluded from the fits) leads to an overall positive value of the intercept q.In fact, this evidence is counter-intuitive, since the intercepts should be negative because of the above-cited greater diagnostic sensitivity given by the luminosity emission with respect to the energy measurement.
Energies 2023, 16, x FOR PEER REVIEW 14 of 23 presence or the absence of atomic or molecular lines in the emitted spectrum in the visible region of our camera [54,55].This reasoning could possibly justify not only the different luminosity level at the same driving voltage, but also synthetic air/N2 cases characterized by same energy levels but different values of luminosity (for instance, see Figure 10c,d).
From the results of Figure 9, a direct correlation between  and ILg can be inferred.Indeed, a linear relationship among thermal energy and discharge luminosity can be estimated for synthetic air (Figure 11a).As expected from Figure 9, the linearity is lost for low values of ILg (a) and of  , due to the different instruments' sensitivity and energy transfer modes.By excluding those points, the released thermal energy could be expressed as follows (Equation ( 5)): where  equals to 0.011 mJ for N2 and to 0.018 mJ for air (Figure 11b).Therefore, by determining the proper calibration factors, discharge luminosity measurement could be a viable means to compute the released thermal energy into the mixture.It is worth noticing that the lack of linearity in the low energy-luminosity region (points excluded from the fits) leads to an overall positive value of the intercept q.In fact, this evidence is counterintuitive, since the intercepts should be negative because of the above-cited greater diagnostic sensitivity given by the luminosity emission with respect to the energy measurement.

Branching
In Figure 10e,f, streamers are fully propagated over the BDI hemisphere, thus giving the chance to discuss air and nitrogen test outcomes from a different perspective.More specifically, these images suggest a peculiar dissimilarity in terms of streamers structure.

Branching
In Figure 10e,f, streamers are fully propagated over the BDI hemisphere, thus giving the chance to discuss air and nitrogen test outcomes from a different perspective.More specifically, these images suggest a peculiar dissimilarity in terms of streamers structure.In fact, the discharge events in air seem to be characterized by a less-prominent branching than those in nitrogen, i.e., the streamers' contour is smoother.To quantitatively test the branching, the ratio between the streamers area and the corresponding perimeter is evaluated at a different voltage level, which will be referred to as aspect ratio.The latter is defined considering the ratio between the number of white pixels contained in the discharge area and those pixels belonging to the streamers' perimeter, i.e., this parameter is nondimensional.Figure 12a reports the cropping and binarization procedure adopted for the calculation, whose details were introduced in the Materials and Methods section.It is important to note that the circular image crop is necessary to correctly estimate the active discharge area which can otherwise be overestimated, due to light diffusion, at the base of the igniter cupola.Figure 12b depicts the ratio between the discharge binarized area and its corresponding perimeter as a function of the integral luminosity I Lg for three discharge events with equivalent luminosity (different colors).Each marker is representative of the average behavior of the discharge, i.e., for the corresponding tested point, the mean binarized area and perimeter are computed as the mean of the single frame, then along the single discharge event and then between the 25 recorded events, ending with the scalar numbers reported in Figure 12b. of the igniter cupola.Figure 12b depicts the ratio between the discharge binarized area and its corresponding perimeter as a function of the integral luminosity ILg for three discharge events with equivalent luminosity (different colors).Each marker is representative of the average behavior of the discharge, i.e., for the corresponding tested point, the mean binarized area and perimeter are computed as the mean of the single frame, then along the single discharge event and then between the 25 recorded events, ending with the scalar numbers reported in Figure 12b.As seen in Figure 12b, increasing values of ILg result in greater aspect ratios for either air (triangular markers) and nitrogen (square markers).Moreover, at different values of integral luminosity, there exists a relative difference among the two media, in terms of binarized-area-to-perimeter ratio, thus evidencing again the oxygen effect on the streamers number and propagation length [56,57].Specifically, for each of the three cases, the higher aspect ratio found in air highlights the presence of a larger perimeter of the streamers in nitrogen.In other words, streamers tend to branch more in pure N2 than in the case with 20% O2.This analysis supports the previously inferred discharge feature which can be qualitatively visualized in Figure 13, for the three mentioned luminosity levels (associated with the markers).As seen in Figure 12b, increasing values of I Lg result in greater aspect ratios for either air (triangular markers) and nitrogen (square markers).Moreover, at different values of integral luminosity, there exists a relative difference among the two media, in terms of binarized-area-to-perimeter ratio, thus evidencing again the oxygen effect on the streamers number and propagation length [56,57].Specifically, for each of the three cases, the higher aspect ratio found in air highlights the presence of a larger perimeter of the streamers in nitrogen.In other words, streamers tend to branch more in pure N 2 than in the case with 20% O 2 .This analysis supports the previously inferred discharge feature which can be qualitatively visualized in Figure 13, for the three mentioned luminosity levels (associated with the markers).Generally, streamers branching is an unpredictable phenomenon that leads to a deeply stochastic discharge morphology.The streamers' structure, both in time and in space, is connected to the fluctuations in the electron density in the barely ionized region in front of the streamer head [56,57].The above-mentioned behavior could be connected to the fact that streamers propagation is affected by the photo-ionization process, which, in turn, is strongly influenced by the presence of oxygen species [58].Upon being a lowelectron attacher (N2 is non-attaching), molecular oxygen substantially modifies the photoionization range of the medium, resulting in noticeable changes in the propagation features of streamers [59].Reducing the oxygen concentration causes the nonlocal photoionization decrease, and therefore local instabilities and more branching.According to Li et al. [60], the streamers produced by high-purity N2 showed the branching much earlier in time and a longer distance in space, at the same operating pressure than in air.The oxygen concentration smooths the electron number density and causes the development of an inception cloud-a quasi-spherical ionization structure close to the electrode tip [61]that avoids the streamers branching as long as they break up [60].Research in the literature on streamers morphology revealed coincidental observations, pointing out the effects on streamer structure of the addition of oxygen to pure nitrogen [62][63][64].

Activation Time Sweep
Activation time is the key control parameter affecting the total amount of thermal energy released by the discharge within the medium.In other words, discharge duration is directly correlated to the capability of the igniter to trigger combustion [33].In an opposite manner to conventional igniters, which are limited by the quantity of energy stored in the coil, the BDI can be continuously powered by the supply system; therefore, the released energy increases over time.Similarly, the cumulative value of luminosity, ILg, is expected to rise as the discharge proceeds, being proportional to the discharge duration.
In Figure 14a,c, thermal energy and luminosity measurements performed in air are plotted against the activation time for three different Vd.The same data representation is mirrored for the results found with nitrogen (Figure 14b,d).Regardless of the medium Generally, streamers branching is an unpredictable phenomenon that leads to a deeply stochastic discharge morphology.The streamers' structure, both in time and in space, is connected to the fluctuations in the electron density in the barely ionized region in front of the streamer head [56,57].The above-mentioned behavior could be connected to the fact that streamers propagation is affected by the photo-ionization process, which, in turn, is strongly influenced by the presence of oxygen species [58].Upon being a low-electron attacher (N2 is non-attaching), molecular oxygen substantially modifies the photoionization range of the medium, resulting in noticeable changes in the propagation features of streamers [59].Reducing the oxygen concentration causes the nonlocal photoionization decrease, and therefore local instabilities and more branching.According to Li et al. [60], the streamers produced by high-purity N2 showed the branching much earlier in time and a longer distance in space, at the same operating pressure than in air.The oxygen concentration smooths the electron number density and causes the development of an inception cloud-a quasi-spherical ionization structure close to the electrode tip [61]-that avoids the streamers branching as long as they break up [60].Research in the literature on streamers morphology revealed coincidental observations, pointing out the effects on streamer structure of the addition of oxygen to pure nitrogen [62][63][64].

Activation Time Sweep
Activation time is the key control parameter affecting the total amount of thermal energy released by the discharge within the medium.In other words, discharge duration is directly correlated to the capability of the igniter to trigger combustion [33].In an opposite manner to conventional igniters, which are limited by the quantity of energy stored in the coil, the BDI can be continuously powered by the supply system; therefore, the released energy increases over time.Similarly, the cumulative value of luminosity, I Lg , is expected to rise as the discharge proceeds, being proportional to the discharge duration.
In Figure 14a,c, thermal energy and luminosity measurements performed in air are plotted against the activation time for three different V d .The same data representation is mirrored for the results found with nitrogen (Figure 14b,d).Regardless of the medium investigated, it can be readily observed that a linear trend for either E r or I Lg against t on is present, at all the driving voltages analyzed.The energy release trend highlights the capability of the power supply system to seamlessly support the igniter energy request throughout the operating range analyzed, whereas the second result confirms that I Lg is proportional to the discharge duration.As seen in [17], below 300 µs the sampling step was intensified to 50 µs in order to analyze the above-discussed trends for low activation time.The purpose is to discern if the overall released energy and luminosity are affected by the rising transient time from streamers inception to full stable propagation.Hence, Figure 14 contains subplots displaying the E r or I Lg values up to the first 300 µs.Focusing on the latter permits us to confirm that a linear relationship also exists for short-lasting discharges.Except for a few cases (i.e., see in Figure 14c the point at t on = 150 µs and V d = 55 V), the experimental points of both E r and I Lg lie on the corresponding linear regression fits, which shows a coefficient of determination of the regression line R 2 higher than 0.998, in any case.For both gases, at same t on , it is confirmed that a higher V d leads to an increased amount of released energy.At same V d and t on , air features a lower level of E r in relation to the nitrogen case, as discussed in the previous section.It is worth highlighting that, as the igniter activation time increases, the broadening between iso-voltage curves becomes gradually substantial among nitrogen and air, which is a direct consequence of the linearity found against V d and t on .
Energies 2023, 16, x FOR PEER REVIEW 17 of 23 investigated, it can be readily observed that a linear trend for either E or ILg against ton is present, at all the driving voltages analyzed.The energy release trend highlights the capability of the power supply system to seamlessly support the igniter energy request throughout the operating range analyzed, whereas the second result confirms that ILg is proportional to the discharge duration.As seen in [17], below 300 μs the sampling step was intensified to 50 μs in order to analyze the above-discussed trends for low activation time.The purpose is to discern if the overall released energy and luminosity are affected by the rising transient time from streamers inception to full stable propagation.Hence, Figure 14 contains subplots displaying the  or ILg values up to the first 300 μs.Focusing on the latter permits us to confirm that a linear relationship also exists for short-lasting discharges.Except for a few cases (i.e., see in Figure 14c the point at ton = 150 μs and Vd = 55 V), the experimental points of both  and ILg lie on the corresponding linear regression fits, which shows a coefficient of determination of the regression line R 2 higher than 0.998, in any case.For both gases, at same ton, it is confirmed that a higher Vd leads to an increased amount of released energy.At same Vd and ton, air features a lower level of  in relation to the nitrogen case, as discussed in the previous section.It is worth highlighting that, as the igniter activation time increases, the broadening between iso-voltage curves becomes gradually substantial among nitrogen and air, which is a direct consequence of the linearity found against Vd and ton.The above-discussed trends also suggest that the energy release mechanism and the photon emission process do not change their qualitative characteristic against time either in the presence of oxygen or not (constant rate of thermal energy and luminous intensity).This supports again, even in the case of synthetic air, the quest for a direct relationship between deposited energy and emitted luminosity.Figure 15 reports plots of released thermal energy against the integral luminosity value for the previous three analyzed voltage.Graphs on the left-hand side depict the whole data range of Er and ILg, whereas those on the right-hand side offer a corresponding close-up view to better appreciate the relation The above-discussed trends also suggest that the energy release mechanism and the photon emission process do not change their qualitative characteristic against time either in the presence of oxygen or not (constant rate of thermal energy and luminous intensity).This supports again, even in the case of synthetic air, the quest for a direct relationship between deposited energy and emitted luminosity.Figure 15 reports plots of released thermal energy against the integral luminosity value for the previous three analyzed voltage.Graphs on the left-hand side depict the whole data range of E r and I Lg , whereas those on the right-hand side offer a corresponding close-up view to better appreciate the relation among the two discussed quantities, especially for synthetic air.As it can be seen, all the tested points could be approximated on two linear fitting lines (one for each tested gas), attesting to the existence of an energy-luminosity correlation.
among the two discussed quantities, especially for synthetic air.As it can be seen, all the tested points could be approximated on two linear fitting lines (one for each tested gas), attesting to the existence of an energy-luminosity correlation.For both gases, a reduction in the angular coefficient m is found as Vd increases.This decrease results to be higher for the air case (see Table 3).o At each Vd, nitrogen showed lower m values than air, since, for the same ILg, a lower amount of thermal energy has been recorded.Vice versa, for the same  , brighter streamers are present in nitrogen with respect to those analyzed in synthetic air.o Positive values of the intercept q, either in air or in nitrogen, are due to the lack of linearity in the low energy-luminosity region (points excluded as already mentioned for fit calculation).For sake of completeness, the following considerations, concerning the results of Figure 15, are reported below: For both gases, a reduction in the angular coefficient m is found as V d increases.This decrease results to be higher for the air case (see Table 3).At each V d , nitrogen showed lower m values than air, since, for the same I Lg , a lower amount of thermal energy has been recorded.Vice versa, for the same E r , brighter streamers are present in nitrogen with respect to those analyzed in synthetic air.Positive values of the intercept q, either in air or in nitrogen, are due to the lack of linearity in the low energy-luminosity region (points excluded as already mentioned for fit calculation).Based on the above considerations, the different slopes found in air and nitrogen emphasize that it is possible to determine calibration factors among thermal energy and luminosity.Their correct determination and their dependence on the tested media are of paramount importance if a similar approach is to be replicated on an optically accessible engine.To conclude, it must be pointed out that, as mentioned in the previous section, streamers luminosity emission is strongly influenced by the chemical species produced during the discharge events.As a consequence, a spectroscopical and chemiluminescent analysis could be suitable to gain further insights into the discharge physics, so as to profoundly quantify the luminosity-thermal energy correlation.
Actually, the developed method of image post-processing is also quite strict in the parameter value selection, because it consists of the implementation and adaptation of different algorithms imported from other fields of the image analysis.Furthermore, changing the adopted parameters also has an impact on the theory behind these methods, so the statistical meaning of the results can lose its validity.
Currently, our group is involved in the development of a new kind of CNN-based (convolutional neural network, a class of artificial neural network) algorithm to remove algorithm dependence from the particular parameters chosen, so as to improve the universality of the method and to widely reduce the use complexity and the user arbitrariness.The recognition of the flame evolution is carried out by PA, regardless of any user's decision, which could make the proposed method potentially suitable for any type of application requiring a high degree of objectivity.The new algorithm, applied to simplified but realistic images, shows promising results [65,66].

Conclusions
The present study analyzed the discharge characteristics of a barrier discharge igniter by simultaneously recording pressure data and natural luminosity images employing an optically accessible calorimeter.Test are performed in synthetic air at 8 bar absolute, room temperature and by executing a parametric investigation of the main setting parameters of the igniter, namely, driving voltage and activation time.The data gathered are compared with the results obtained with nitrogen as a test medium in a previous work of the same group.The relevant outcomes of the study are:

•
For a fixed activation time in synthetic air test, thermal energy and discharge luminosity show a linear relationship with the driving voltage in a similar manner to nitrogen.At the lowest voltages (close-to-inception conditions), non-linear behavior among the previous quantities is commonly present in both media.

•
Considering the same driving voltage, in synthetic air, a simultaneous reduction in released thermal energy and luminosity is detected with respect to pure nitrogen.

•
In common with N 2 , in synthetic air a direct linear correlation between deposited thermal energy and luminosity can be found.It must be noted that different calibration factors between the two media are identified.

•
The presence of oxygen modifies the structure of the streamers, resulting in a thicker and less branched morphology in the case of synthetic air with respect to N 2 .
In this experimental investigation, air was chosen since it approaches the in-cylinder mixture composition (fresh air, fuel vapor, exhaust residual) present in optical and metal engines.In principle, it could be feasible to exploit the outcome of such non-intrusive diagnostic technique to characterize the igniter's energetic behavior in an optical access engine.Future works are aimed at proving and supporting the above assumption, by evaluating the influence of different media (mixture of air and carbon dioxide for EGR simulation) and igniters on the energy-luminosity correlation.

Energies 2023 , 23 Figure 1 .
Figure 1.(a) 1: inductor, 2: connection, 3: firing end of the igniter and corresponding detail of the firing end.(b) BDI lumped-parameter circuit.The coil is represented by an inductive component L1, a resistance R1 which takes into account skin and proximity effects losses and ohmic losses, and a parasitic capacitance C1 between the turns.The firing end mostly represents the capacitance CL of the igniter assembly.A resistance RL, in parallel, is introduced to account for the ohmic losses of air and bulk insulation.(c) Streamers development around the BDI cupola as driving voltage increases from the inception up to the maximum allowed voltage.

Figure 1 .
Figure 1.(a) 1: inductor, 2: connection, 3: firing end of the igniter and corresponding detail of the firing end.(b) BDI lumped-parameter circuit.The coil is represented by an inductive component L 1 , a resistance R 1 which takes into account skin and proximity effects losses and ohmic losses, and a parasitic capacitance C 1 between the turns.The firing end mostly represents the capacitance C L of the igniter assembly.A resistance R L , in parallel, is introduced to account for the ohmic losses of air and bulk insulation.(c) Streamers development around the BDI cupola as driving voltage increases from the inception up to the maximum allowed voltage.

Figure 3 .
Figure 3. (a) FTT amplitude of the recorded pressure signal (grey signal) due to the igniter discharge and corresponding filtered signal (red line).(b) Corresponding pressure signals over the time.

Figure 3 .
Figure 3. (a) FTT amplitude of the recorded pressure signal (grey signal) due to the igniter discharge and corresponding filtered signal (red line).(b) Corresponding pressure signals over the time.

Figure 4 .
Figure 4. (a) Representation in MATLAB environment, in grey scale (0-255), of a frame acquired by the camera.(b) Frame selection process performed to identify the beginning and termination of the discharge event.

Figure 4 .
Figure 4. (a) Representation in MATLAB environment, in grey scale (0-255), of a frame acquired by the camera.(b) Frame selection process performed to identify the beginning and termination of the discharge event.

Figure 5 .
Figure 5. (a) Average grey level for each of the first 20 frame (black markers) and corresponding mean value of them (red dashed line).(b) Grey level distribution of the average noise frame.(c) Example of an original frame before (grey bars) and after (blue bars) the noise removing.

Figure 5 .
Figure 5. (a) Average grey level for each of the first 20 frame (black markers) and corresponding mean value of them (red dashed line).(b) Grey level distribution of the average noise frame.(c) Example of an original frame before (grey bars) and after (blue bars) the noise removing.

Figure 6 .
Figure 6.Workflow to obtain, at specific operating point, the average grey level and the binarized area for estimating the streamer surface around the igniter cupola.

Figure 6 .
Figure 6.Workflow to obtain, at specific operating point, the average grey level and the binarized area for estimating the streamer surface around the igniter cupola.

) 23 Figure 7 .
Figure 7. Trends of Lg (black lines) composed by the sequence of all the grey level values lg of each discharge and the corresponding mean trend  g (red line) of all the discharges.Note that the y-axis physical quantity is the average grey level lg in Equation (2).The parameter  is indicated by the red transparency area below the  g curve.
depicts the typical pressure signals arising from discharge events and the corresponding BDI globe images at Vd CI , intermediate Vd and Vd MAX .

Figure 7 .
Figure 7. Trends of L g (black lines) composed by the sequence of all the grey level values l g of each discharge and the corresponding mean trend L g (red line) of all the discharges.Note that the y-axis physical quantity is the average grey level l g in Equation (2).The parameter I Lg is indicated by the red transparency area below the L g curve.

23 Figure 8 .
Figure 8. Chamber relative pressure signal over time at three different driving voltages, in a range from Vd CI to Vd MAX .

Figure 8 .
Figure 8. Chamber relative pressure signal over time at three different driving voltages, in a range from V d CI to V d MAX .

Energies 2023 ,
16, x FOR PEER REVIEW 12 of 23the BDI cupola in presence of air (left column) and nitrogen (right column).The corresponding applied Vd and the resulting  and ILg are also therein reported, all of them increasing from the top to the bottom of the figure.Each row reports a pair of images with approximately the same amount of released thermal energy.

Figure 9 .
Figure 9. ILg (a) and  (b) versus the driving voltage in air (black markers) and nitrogen (blue markers) at all the voltage level investigated.Dashed black and blue vertical lines depict the voltage threshold over which a measurable release of thermal energy occurs.The dashed red lines represent a linear fit of ILg for both analyzed media.

Figure 9 .
Figure 9.I Lg (a) and E r (b) versus the driving voltage in air (black markers) and nitrogen (blue markers) at all the voltage level investigated.Dashed black and blue vertical lines depict the voltage threshold over which a measurable release of thermal energy occurs.The dashed red lines represent a linear fit of I Lg for both analyzed media.

Figure 9 .
Figure 9. ILg (a) and  (b) versus the driving voltage in air (black markers) and nitrogen (blue markers) at all the voltage level investigated.Dashed black and blue vertical lines depict the voltage threshold over which a measurable release of thermal energy occurs.The dashed red lines represent a linear fit of ILg for both analyzed media.

Figure 10 .
Figure 10.Natural luminosity images of streamers evolution in air (a,c,e) and nitrogen (b,d,f) at different driving voltages.

Figure 11 .
Figure 11.Correlation between  and ILg for all tested Vd in air (a) and comparison with nitrogen results (b).The dashed red lines represent a linear fit between  and ILg for both analyzed media.

Figure 11 .
Figure 11.Correlation between E r and I Lg for all tested V d in air (a) and comparison with nitrogen results (b).The dashed red lines represent a linear fit between E r and I Lg for both analyzed media.

Figure 12 .
Figure 12.(a) Igniter discharge image cropping for streamers binarized area and perimeter estimation; (b) Discharge aspect ratio as a function of the integral luminosity in pure N2 (square markers) and synthetic air (triangular markers).Same color markers (green, orange, red) represent discharge events with equivalent luminosity among the two media.

Figure 12 .
Figure 12.(a) Igniter discharge image cropping for streamers binarized area and perimeter estimation; (b) Discharge aspect ratio as a function of the integral luminosity in pure N 2 (square markers) and synthetic air (triangular markers).Same color markers (green, orange, red) represent discharge events with equivalent luminosity among the two media.

Figure 13 .
Figure 13.False colored iso-luminosity discharge images for air (left column) and nitrogen (right column).

Figure 13 .
Figure 13.False colored iso-luminosity discharge images for air (left column) and nitrogen (right column).

Figure 14 .
Figure 14.Thermal energy released (a,b) and luminous emission (c,d) plotted against the ton at three Vd level in air (black markers) and nitrogen (blue markers).Each graph contains a boxplot focused on the interval from 50 μs to 300 μs.

Figure 14 .
Figure 14.Thermal energy released (a,b) and luminous emission (c,d) plotted against the ton at three V d level in air (black markers) and nitrogen (blue markers).Each graph contains a boxplot focused on the interval from 50 µs to 300 µs.

Figure 15 .
Figure 15.Correlation between  and ILg in N2 (blue markers) and synthetic air (black markers) plotted against each activation time at the three voltage levels investigated.Left-hand side graphs report energy and luminosity data concerning the whole analyzed range and their linear fits (black and blue dashed lines).Plots on the right-hand side represent a magnification around the origin of the trend seen on the left-hand side.For sake of completeness, the following considerations, concerning the results of Figure 15, are reported below: oFor both gases, a reduction in the angular coefficient m is found as Vd increases.This decrease results to be higher for the air case (see Table3). oAt each Vd, nitrogen showed lower m values than air, since, for the same ILg, a lower amount of thermal energy has been recorded.Vice versa, for the same  , brighter streamers are present in nitrogen with respect to those analyzed in synthetic air.oPositive values of the intercept q, either in air or in nitrogen, are due to the lack of linearity in the low energy-luminosity region (points excluded as already mentioned for fit calculation).

Figure 15 .
Figure 15.Correlation between E r and I Lg in N 2 (blue markers) and synthetic air (black markers) plotted against each activation time at the three voltage levels investigated.Left-hand side graphs report energy and luminosity data concerning the whole analyzed range and their linear fits (black and blue dashed lines).Plots on the right-hand side represent a magnification around the origin of the trend seen on the left-hand side.

Table 3 .
Angular coefficients and intercept values corresponding to the linear fits of   versus ILg in air and nitrogen.

Table 3 .
Angular coefficients and intercept values corresponding to the linear fits of E r versus I Lg in air and nitrogen.