Transformation of n-heptane using an in-liquid submerged microwave plasma jet of argon

The reforming of hydrocarbons has gained much interest as a means to upgrade low-grade fuels and to produce value-added chemicals. Plasmas have been considered one of the potential ways to reform fuels to achieve more effective and cleaner combustion, particularly by producing various hydrocarbons, hydrogen carriers, and oxygenates as well as syngas. Here, we employed a submerged microwave plasma jet of argon to investigate its potential to transform n-heptane. We found that the product selectivities were mainly governed by the effective gas temperature, which we adjusted by changing the energy density of the argon stream. The transformation of n-heptane by this method mostly produced ethylene and acetylene, which is different than the products produced by pyrolysis or a chemical equilibrium composition. Such unique selectivities could be attributed to the rapid quenching of the microwave plasma jet upon direct contact with the colder liquid. The transformation of n-heptane was significantly affected by the interactions between the microwave plasma jet and the liquid n-heptane. To support our results, we include a detailed chemical analysis and discussion of the physical characterization of the microwave plasma jet using optical emission spectroscopy. © 2021 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/). https://doi.org/10.1063/5.0036041


I. INTRODUCTION
The reforming of fuel has been of special interest ever since the Fischer-Tropsch process, i.e., the liquefaction of syngas from coal, was first introduced in the early 20th century. Recently, the study of the reforming of fuel has diversified, partly in response to climate change issues, such as using CO 2 as feed or producing carbon-neutral (or carbon-free) fuels, and partly to cope with intermittent renewable electricity storing into chemical energy.
Since the early 2000s, nonthermal plasmas have been considered as a promising method for the on-board production of hydrogen-rich syngas, particularly because plasma devices are compact and respond quickly. 1 Dry reforming, partial oxidation, and steam reforming of methane have been popular chemical processes in the context of plasma-based fuel reforming. [2][3][4][5] Recently, much attention has been given to modeling the fundamental underlying plasma chemistry, [5][6][7][8] and recent experiments [9][10][11] have helped to identify the respective roles of energetic electrons and gas temperature in plasma chemistry; electrons initiate the chemical process by creating radicals, and the temperature determines the product distribution.
Recent advances in the field of in-liquid discharges have inspired the exploration of the reforming of liquid hydrocarbons. [12][13][14][15][16] These approaches practically do not require the evaporation of liquid fuels; thus, these can accommodate various kinds of feedstock and save the energy for the evaporation as compared to gas phase plasma methods. Studies have shown that the product composition of in-liquid plasma processes to transform liquid hydrocarbons can be tailored by adding chemical additives. [18][19][20] Particularly, methane and carbon dioxide were found to selectively increase the production of lighter hydrocarbons and oxygenates, respectively, in the reforming of n-dodecane 19 and n-heptane. 20 Thus, there exist at least two ways to control the product composition: by controlling the reaction temperature and the choice of chemical additives.
n-Heptane has been considered as a model liquid fuel; therefore, its reforming has been extensively studied. Pyrolysis-induced n-heptane transformation is one of the popular topics in this subject area. Bajus et al. 21 have proposed one of the well-known mechanisms, which is based on the Rice-Kossiakoff theory. The unimolecular scission of primary and secondary C-C bonds (C 7 H 16 → R + R 0 ; R and R 0 are alkyl radicals) and the cleavage of the C-H bonds to form heptyl radicals (C 7 H 16 → H + C 7 H 15 ) have been proposed as the two main initiation steps. Also, in the presence of small radicals, such as H, CH 3 , C 2 H 5 , and C 3 H 7 , hydrogen abstraction from C 7 H 16 can occur to produce C 7 H 15 and other light saturated molecules (H 2 , CH 4 , C 2 H 6 , and C 3 H 8 ). Then, further decomposition of C 7 H 15 radical occurs via the following four pathways: Pant et al. 22 have studied the kinetics and product distribution during n-heptane pyrolysis. They have reported that the n-heptane decomposition could be represented by a first-order reaction and proposed a comprehensive reaction scheme to address the various ways of product formation. Wang et al. 23 have studied n-heptane reforming using gliding arc discharges. Gaseous n-heptane mixed with air has been reformed into hydrogen-rich syngas. Particularly, Lebedev et al. 24 have performed a numerical study for microwave discharges in a small bubble of n-heptane vapor. They have found that the role of electrons in the dissociation of n-heptane could be important only during a very short initial stage of time (∼1 ms), indicating that a main role of microwave plasma is thermal. On the other hand, Reddy and Cha 20 have investigated the effect of additives (CH 4 , CO 2 , and H 2 O) on n-heptane reforming using a dielectric barrier discharge (DBD) reactor at room temperature, minimizing the role of plasma's thermal effect. They have shown that the additives could selectively increase lighter hydrocarbons or oxygenates in liquid phase. Recently, we have developed a submerged microwave plasma jet (MWPJ) 15 and a submerged arc plasma jet 16 to overcome the following technical challenges in general in-liquid pulsed discharges: (i) a significant dependence of discharges on the liquid's properties, such as permittivity and electrical conductivity, and (ii) a limited, relatively low level of electrical power to the liquid. We have decoupled plasma generation from the liquid's properties by injecting a remotely produced plasma jet into the liquid and achieved higher electrical power by using an arc torch or a microwave plasma jet. These submerged plasma jet systems have been verified in water, which has a relatively high relative permittivity (∼80) and a wide range of conductivity by decomposing organic compounds. 15,16 The submerged plasma jet systems are required to be further applied to transform liquid hydrocarbons, which have very low permittivity (∼2) and negligible conductivity.
In the present study, we employed a submerged microwave plasma jet (MWPJ) to transform liquid fuel. Main motivation to adopt the MWPJ was to investigate the effect of hotter plasma on product distribution, which is expected to be very different from the result in the DBD reactor with additives. 20 Particularly, because the submerged configuration was designed to quench the plasma by ambient liquid, we expected to find another way of controlling a product composition. Argon was used as the feed gas to the MWPJ, and n-heptane was used as the liquid fuel. To focus on the thermal effect of the MWPJ by maintaining the chemistry as simple as possible, we considered no additives other than Ar. We conducted a product analysis to obtain a conversion of n-heptane and product selectivities. In particular, we focused on the product selectivities to highlight contrasts between the submerged MWPJ and other results, such as pyrolysis and chemical equilibrium composition. To estimate gas temperature and electron density in the MWPJ, a spectroscopic analysis of light emission was also conducted. As a result, we presented a unique feature of the plasma chemistry in the submerged MWPJ system.

II. EXPERIMENTAL SETUP
A submerged MWPJ system to transform liquid n-heptane was employed, which consisted of a microwave power system and a quartz reservoir with a MWPJ inlet (Fig. 1). A solid-state microwave generator (Sairem, GMS200W) delivered continuous microwaves at 2.45 GHz to a surfatron device (Sairem, Surfatron60) via a coaxial cable up to 200 W of adjustable microwave power, P dis , which was estimated by substracting the reflected power from the incident power. A quartz tube with an outer and an inner diameter of 6.0 and 4.0 mm, respectively, was inserted at the center of the surfatron to guide the plasma gas. The surfatron was water-cooled, and the outer part of the quartz tube was air-cooled. Argon (99.999% purity) was selected as the plasma gas, which was supplied through the quartz tube. The flow rate of argon, Q Ar , was varied in a range of 1-5 l/min using a mass flow controller. One end of the quartz tube was welded to the quartz reservoir (4.6 × 4.6 cm 2 square cylinder), which was filled with n-heptane (V hep = 200 ml, 9.45-cm height in the reservoir), such that the Ar MWPJ was ejected into the liquid n-heptane directly. A gas outlet was placed at the upper part of the reservoir. We conducted the experiment in atmospheric pressure.
We investigated the physical characteristics of the submerged MWPJ using optical emission spectroscopy (OES). Particularly, in order to estimate the gas temperature, the rotational temperature of C 2 was calculated by fitting the experimental spectrum of the C 2 Swan band (Δν = 0) using the method developed by Cardoso et al. 17 A spectrometer (Princeton Instruments, SP2750), which was equipped with a grating (900 grooves/mm, blazed at visible) and an intensified charge-coupled device (ICCD) camera (Princeton Instruments, PI-MAX3), was used to acquire the spectral data of the emissions. The exposure time for each acquired spectrum was 1 ms, and 50 spectra were averaged for each measured condition. A fiber optics was installed at the outside of the reservoir to remotely collect the light emission from the inlet of the MWPJ to the reservoir, and the measuring volume was 10 mm in diameter. Thus, we obtained temporally and spatially averaged spectra. A high-speed camera (LaVision, Image Pro HS) was used to investigate the dynamic behavior of the plasma-bubble interactions. A lightemitting diode lamp (Luminus, SST-90) illuminated the reservoir from the opposite side to the high-speed camera to clearly illuminate the interface between gas and liquid. The exposure time and the framing rate of the high-speed imaging were set to 0.47 ms and 2 kHz, respectively.
We evaluated the chemical effects of the submerged MWJP on the transformation of n-heptane based on the resulting chemical products. The gas products were analyzed using a gas chromatography, GC (Agilent, GC 7890A), and the liquid products were identified by gas chromatography-mass spectrometry (Agilent, GC-MS/FID 7890A-5975C). The total plasma processing time, t proc , with MWPJ was maintained at 5 min for all tested conditions. Note that, at P dis = 200 W and t proc ¼ 5 min, the increase in the temperature of liquid n-heptane was ∼40 K, which may have negligible influence on the chemical reaction. The chemical analysis was repeated five times for each experimental condition, and the results showed less than ±3% deviation from the averaged concentrations for each species. Thus, we plotted the data based on the averaged concentrations.

III. RESULTS AND DISCUSSION
A. Chemical assessment of the submerged MWPJ system We investigated the transformation characteristics of the submerged MWPJ to liquid n-heptane with particular emphasis on the conversion of n-heptane and the selectivities of the produced components. To assess the conversion and the selectivities from the analytical data obtained from the GC and the GC-MS analyses, we developed the following equations.
For a species x in the gas products, its outlet molar flow rate, _ n x, gas , can be expressed as _ n x,gas ¼ C x,gas Á ( _ n Ar /C Ar ), where _ n Ar denotes the inlet molar flow rate of Ar. C x,gas and C Ar denote the outlet concentrations of x and Ar measured by the GC, respectively. The total produced moles of x, n x,gas , can be calculated as n x,gas ¼ _ n x,gas Á t proc . Thus, the total amount of carbon atoms (N C,gas , mole-based) in the gas products can be obtained from where i is the number of carbon atoms in the species x.
For the remaining liquid in the reservoir, which included unconverted n-heptane, the total amount of carbon atoms (N C,res , mole-based) is calculated using a carbon balance as N C,res ¼ 7 Á n C7H16 À N C,gas , where n C7H16 is the initial molar amount of n-heptane in the reservoir. Thus, the produced moles of y species (n y,liq ) in the liquid products can be estimated by n y,liq ¼ (N C,res Á C y,liq )/j, where j is the number of carbon atoms in the produced liquid species y, and C y,liq is the carbon-based concentration of the product y obtained from the GC-MS/FID analysis. Thus, the total amount of carbon atoms (N C,liq , mole-based) only as the form of transformed liquid, except for unconverted n-heptane, can be calculated as From Eqs (1) and (2), the molar amount of converted n-heptane (n C7H16,conv ) and the rate of n-heptane conversion (X hep ) can be obtained as In the gaseous product, we identified H 2 , CH 4 , C 2 H 2 , C 2 H 4 , C 2 H 6 , C 3 H 6 , C 4 H 8 , C 3 H 8 , C 4 H 10 , and C 5 H 10 , whereas in the liquid product, we found C 4 H 8 , C 5 H 10 , and C 6 H 12 . Note that no solid particles were found. The selectivity of the produced H 2 (S H2 ) is given by S H2 ¼ 1 8 n H 2 ,gas n C7H16,conv : For a gaseous and liquid hydrocarbon product x and y having i and j carbon atoms, respectively, their selectivities based on carbon atom, S x and S y , are S y ¼ j 7 n y,liq n C7H16,conv : Based on these calculations, we found that the energy density (or specific input energy) is not a useful parameter to explain n-heptane conversion in the submerged MWPJ system, unlike most of the hydrocarbon reactions in gaseous discharges that have been described in previous studies. [9][10][11] Energy density is defined as a density of the delivered electrical energy into a control volume and, in this study, can be defined for two different control volumes: one is the liquid n-heptane and the other is the flowing argon. The energy density for the liquid n-heptane, ED hep = P dis ⋅t proc /V hep , may play an important role in characterizing X hep in a similar way to the decomposition of organic compounds that are dissolved in water using a submerged arc jet system. 16 The conversion we employed here seemed to be reasonably proportional to ED hep when P dis is varied at Q Ar = 2 l/min. Meanwhile, X hep was significantly affected by Q Ar , which does not affect ED hep [ Fig. 2(a)]. We found one order of difference as X hep = 20.2 and 2.77 μmol/s at Q Ar = 2 and 5 l/min, respectively, for fixed P dis = 200 W. This indicates that the effectiveness of P dis to X hep is influenced by the quality of the MWPJ and the interaction between the MWPJ and the liquid.
Although the energy density in the argon flow, ED Ar = P dis /Q Ar , was considered instead of ED hep , we found no generalized trend for the conversion of n-heptane in terms of ED Ar [ Fig. 2(b)]. A strong correlation of X hep with increased ED Ar for the varying P dis conditions at Q Ar = 2 l/min was still found. However, for the varying Q Ar cases at P dis = 200 W, X hep showed a maximum at ED Ar = 6 kJ/l (Q Ar = 2 l/min), with unchanged X hep in the range of 4 ≤ ED Ar ≤ 12 kJ/l (1 ≤ Q Ar ≤ 3 l/min).
Interestingly, we found the selectivities of the products to be well correlated with the ED Ar , as highlighted by the shaded area in Fig. 3. We identified H 2 , alkane (CH 4 , C 2 H 6 , C 3 H 8 , and C 4 H 10 ), alkene (C 2 H 4 , C 3 H 6 , C 4 H 8 , and C 5 H 10 ), and alkyne (C 2 H 2 ) present in the gaseous product, whereas in the liquid product, we only found alkene (C 4 H 8 , C 5 H 10 , and C 6 H 12 ). The submerged MWPJ was highly selective to produce C 2 H 4 followed by C 2 H 2 , C 3 H 6 , C 4 H 8 , H 2 , and CH 4 . Note that although S H2 (<7% for all tested conditions) was comparable with S C4H8 , we have not shown this information in Fig. 3 for the sake of visibility and clarity.
Thus, the well-correlated results in the shaded area raise an interesting point; the reaction temperature might strongly affect the product composition since ED Ar , by definition, may be related to the gas temperature in the MWPJ system. More specifically, our results indirectly confirm that most of the reactions occurred in the gas phase with the evaporated n-heptane in the MWPJ. In our previous studies on the basic reforming processes of methane-partial oxidation, dry reforming, and steam reforming-using a temperature-controlled dielectric barrier discharge (DBD) reactor, [9][10][11] we found that the gas temperature mainly controlled the product selectivity, and the reactant conversion was initiated by electron-induced chemistry, depending on electron temperature and density.
However, our selectivity results at two conditions, P dis = 200 W:Q Ar = 5 l/min and P dis = 50 W: Q Ar = 2 l/min, showed significant deviation from each other. The selectivities of C 2 H 4 , C 2 H 2 , and CH 4 abruptly dropped as Q Ar increased to 5 l/min from 4 l/min at P dis = 200 W (from ED Ar = 3 kJ/l to 2.4 kJ/l) resulting in increased S C3H6 and S C4H8 . These results indicate that the effect of the argon flow on the selectivities and conversion is complicated. Thus, we provide in-depth discussion in Secs. III B and III C by elaborating on the physical characteristics of the MWPJ and detailing the role of ED Ar on the selectivities, respectively.

B. Physical characterization of MWPJ
We characterized the submerged MWPJ using optical emission spectroscopy to investigate its composition, temperature, and electron density. Based on the spectra obtained at P dis = 200 W and Q Ar = 2 l/min (Fig. 4), we identified the spectral lines of Ar (in the range of 700-850 nm), C 2 (Swan: head bands at ∼474, 517, and 564 nm and C 2 -Deslandre-D 0 Azambuja bands bands: head bands at ∼361, 385, and 410 nm), CH (A-X: head band a ∼431 nm, B-X: head band at ∼387 nm, and C-X: head band at ∼315 nm), Hα (at 656.3 nm), and Hβ (at 486.1 nm). NH (A-X: head band at ∼336 nm) was also detected, which we attributed to impurities in Ar or to dissolved air in n-heptane. The observed emission spectra indicated that the initial n-heptane was successfully decomposed by the plasma to hydrogen radicals and CH radicals. Note that since n-heptane is colorless, there is no significant interference in the visible range.
As a first approximation, the estimated rotational temperature of C 2 might be stretched to represent the gas temperature (T g ) for a given condition. 25,26 Figure 5 shows a typical calculated synthetic spectrum, which fits for the experimentally measured C 2 Swan band (Δν = 0), to estimate the rotational temperature of C 2 . We found a good agreement between the simulated and the experimental spectra for all tested conditions; as a result, the estimated gas temperatures were ranged in 1915-2732 K and 1683-2467 K for fixed P dis = 200 W and fixed Q Ar = 2 l/min cases, respectively. Note that the synthetic spectra showed a margin of uncertainty at ±50 K, as we found the best fit with the experimental one.
The energy balance might provide a simple explanation for the increase in gas temperature. Although P dis is dissipated into chemical formation, photon emission and the sensible enthalpy of both gas and liquid, a proportionality between the input power and a resulted temperature change is still valid: P dis ∼ ρQ Ar c p ΔT g , where ρ is the gas density, c p is the specific heat, and ΔT g is the temperature difference. Therefore, ΔT g ∼ P dis /Q Ar can be obtained explaining inverse proportionality of the temperature to Q Ar and its linear proportionality to P dis . Recalling that the parameter P dis /Q Ar = ED Ar , we plotted the estimated gas temperatures in terms of ED Ar in Fig. 6.
As a result, we found that the estimated gas temperature adequately represented ED Ar . Thus, the discrepancies in the  Fig. 3 can be attributed to something other than the gas temperature in the MWPJ.

selectivities shown in
To determine the electron density, we considered Stark broadening of Hα and found that the influence of the electron density on the kinetics of n-heptane transformation was neutral. By fitting the experimentally obtained Hα line with a Voigt profile, we determined the Lorentzian broadening (Δλ FWHM ) induced by the Stark effect considering the broadenings induced by van der Waals, Doppler, and the spectrometer in the fit; thus, the electron density in the MWPJ could be estimated using Gigosos's table. 27 We found that the electron density was in the range of 2.5 × 10 15 -4 × 10 15 cm −3 for all tested conditions. This range of electron densities seemed to be insignificantly effected by a fourfold increase in P dis (50-200 W) and a fivefold variation of Q Ar (1-5 l/min) and ED Ar (2.4-12 kJ/l); thus, electron density's effect on the chemistry should be rather uniform throughout the tested cases. This indicates a less significant role of electron impact reactions in "hotter" plasmas, and a previous study 24 for n-heptane decomposition using an direct in-liquid microwave plasma with a metallic electrode has supported the importance of thermochemistry.
The large differences in the selectivities at Q Ar = 5 l/min and P dis = 200 W shown in Fig. 3 may be attributed to significantly altered plasma dynamics and interactions with the liquid. To explore the causes behind the exceptional selectivities at Q Ar = 5 l/min and P dis = 200 W (Fig. 3) as well as the n-heptane conversions for the cases with varying Q Ar at fixed P dis = 200 W (Fig. 2), we investigated the dynamic interaction between the MWPJ and the liquid. As shown in Fig. 7, for the cases with fixed P dis = 200 W, increased Q Ar resulted not only in a longer plasma channel, but also in a stronger interaction between the plasma and the liquid. At low Q Ar , the plasma existed within a limited region of the confined bubble, and there seemed to be no direct interaction between the plasma and the liquid. However, at Q Ar = 5 l/min, FIG. 6. Estimated gas temperatures (rotational temperature of C 2 ) were well represented by ED Ar for various P dis and Q Ar .

FIG. 7.
Typical instantaneous images for each tested condition; the upper row shows varying Q Ar at P dis = 200 W, while the lower row represents varying P dis at Q Ar = 2 l/min. the void region surrounding the plasma channel was significantly reduced, showing a pinched plasma zone in the liquid. Therefore, we could expect a rapid quenching of the MWPJ due to a direct contact with liquid n-heptane. This quenching process might drastically modify the overall reacting conditions for the n-heptane transformation. Note that the varying P dis conditions (bottom row in Fig. 7) did not show such a pinching phenomenon, demonstrating that the void gas region-a buffer between the plasma and the liquid-for all tested conditions mildly influenced plasma and bubble behavior due to increased P dis .

C. Gas temperature toward product selectivity
We would expect that the gas temperature in the MWPJ has a primary effect on product selectivity regardless of Q Ar and P dis , since previous methane reforming studies [9][10][11] have shown that the gas temperature strongly influences selectivities. We replotted the product selectivities as a function of T g in Fig. 8. Similar to the selectivity result with ED Ar (Fig. 3), an overlap of the selectivities in a T g range of 1950-2450 K was found between both experimental sets, i.e., with varying P dis at Q Ar = 2 l/min and with varying Q Ar at P dis = 200 W. These results were as we expected because the estimated gas temperature showed a strong positive correlation with ED Ar, as shown in Fig. 6.
Results from the chemical equilibrium calculation suggest that a much lower temperature range is required to produce the experimentally obtained product composition. We calculated the chemical equilibrium based on the Gibbs free energy minimization method in MATLAB with pure n-heptane as the initial composition for a range of the equilibrium temperature, T eq . With the obtained molar concentrations of the products at various temperatures, we calculated the selectivity of each product to compare with our experimental results. The chemical equilibrium dictated the product in a range of T eq = 1950-2450 K to have only H 2 and C 2 H 2 , as shown in Fig. 9, where we observed more than 40% C 2 H 4 and moderate levels of C 2 H 2 and C 3 H 8 in the submerged MWPJ system (Fig. 8). To match the experimentally observed product composition to the chemical equilibrium result, we found that 950 < T eq < 1100 K would be the most probable range of temperature. Similar to the experimental result in Fig. 8, we found a peak of S C2H4 with high selectivity (∼30%) together with increasing S C2H2 and decreasing S C3H6 and S C4H8 as the temperature increased. Because there is a steep temperature gradient in the MWPJ due to significant heat loss to the liquid n-heptane, the overall reaction temperature could be reduced compared to the estimated T g using light emission.
Furthermore, the unusual selectivities at P dis = 200 W and Q Ar = 5 l/min might also be due to severer quenching of the MWPJ compared to other Q Ar , as shown in Fig. 7. Due to enhanced heat loss, caused by the direct contact with the liquid and turbulence mixing, the effective reaction temperature of this case should be further decreased compared to other cases. When we projected the temperature below 950 K, to 800 K for example, we obtained a significant drop in S C2H4 and S C2H2 and an increase in S C3H6 and S C4H8 in the equilibrium calculation (Fig. 9), similar to the result with the MWPJ (Fig. 8).
We must highlight a significant contrasting point in terms of CH 4 and C 2 H 4 . According to the chemical equilibrium in 950 < T eq < 1100 K, methane is the most probable product; however, our MWPJ experiment showed less than 5% selectivity for CH 4 , producing C 2 H 4 with more than 40% selectivity. Previous chemical kinetic studies 28,29 have shown that, in a pyrolysis of n-heptane, the formation of CH 4 is much slower than that of C 2 H 4 . In this regard, when a reaction time is limited by a system, the formation of CH 4 is expected to be reduced. The present submerged MWPJ underwent spatial and temporal, rapid quenching of the reaction due to the cold liquid boundary. Thus, there might be not sufficient time for methane to form.  Note that studies for a thermal dissociation (pyrolysis) of n-heptane, both experimentally 31 and theoretically, 28 have found that C 2 H 4 , H 2 , CH 4 , C 3 H 6 , and C 2 H 6 , and often H 2 and C 2 H 4 were the most dominant products, with barely any report on the presence of C 2 H 2 . Meanwhile, C 2 H 2 is the second populous species in the present study. This might be attributed to that the effective reaction temperature of the present experiment happened to correspond to the favorable range of temperature for the formation of C 2 H 2 production as shown in Fig. 9. However, as one-dimensional modeling study 30 for microwave plasma bubble in n-heptane has shown a highly dynamic nature of product distributions in time and space, a detailed modeling study should be required to understand the underlying plasma-chemical effect.
Therefore, we concluded that the observed consistent selectivities, irrespective of various experimental conditions, indicate the dominance of related thermochemistry for the tested ranges of power and flow rate. In addition, due to the nature of the submerged MWPJ, high-temperature jet, and rapid quenching of the chemical reaction, we could obtain the different product composition than that obtained by the pyrolysis and the chemical equilibrium.

IV. CONCLUSIONS
We experimentally investigated a submerged MWPJ of argon in liquid n-heptane to study its transformation capability with liquid substance and to prove its unique character in controlling product selectivity due to the rapid quenching of chemical reactions caused by heat interaction with the liquid. We found that product selectivities are controlled by gas temperature, and thus, the energy density, as defined by the electrical power and the flow rate of argon, was the key parameter for this selectivity. Conversion of n-heptane was primarily proportional to the electrical power but was significantly coupled with the MWPJ and the liquid as the flow rate of argon increased. Comparing the experimental results with the chemical equilibrium compositions, the obtained product selectivity seemed to be governed by the related thermochemistry (pyrolysis). However, the primary products with high selectivities were ethylene and acetylene, which was in contrast to the methane or hydrogen/ethylene via the chemical equilibrium or the pyrolysis, respectively. We determined that this was most likely due to the nature of the submerged MWPJ-rapid quenching. Furthermore, the measured electron density also indirectly indicated that the role of electron-induced chemistry with the submerged MWPJ on the selectivity was insignificant compared to the related thermochemistry.

AUTHORS' CONTRIBUTIONS
A.H. performed the experiment for the physical analysis. J.-L.L. performed the experiment for the chemical analysis. M.S.C. conceived the experiment and constructed the main idea of the study. All authors contributed to the preparation of the manuscript.