Engine Combustion Network – France
Open Access
Numéro
Oil Gas Sci. Technol. – Rev. IFP Energies nouvelles
Volume 75, 2020
Engine Combustion Network – France
Numéro d'article 35
Nombre de pages 24
DOI https://doi.org/10.2516/ogst/2020023
Publié en ligne 5 juin 2020

© M. Ben Houidi et al., published by IFP Energies nouvelles, 2020

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

The Engine Combustion Network (ECN) [1] has greatly contributed to accelerate the research efforts focused on spray atomization and combustion in both Diesel and Gasoline direct injection engines. In this context, experimental and simulation works are joined in a complementary way to improve our understanding of the complex physical and chemical phenomena taking place in fuel injection and combustion process. The efforts from the experimental side are intended to provide a robust and accurate database for model validation, and therefore the boundary conditions need to be accurately characterized. To this end, spray A reference condition has been defined to provide a detailed list of boundary conditions. The conformity of the experimental apparatus employed to these boundary conditions has been discussed in different publications for the first institutions participating to ECN [24]. The analysis of the results indicated that some discrepancies in the spray results are caused by different boundary conditions between these facilities. Other effects related to injector nozzle geometry [57] and to aging [8] have been extensively studied. The boundary conditions of Constant Pressure Flow (CPF) and Constant Volume Pre-burn (CVP) vessels used for spray A studies have been reviewed in [3]. The work reported a detailed comparison of the following facilities: the CPF of CMT [9], the CPF of Spray Combustion Laboratory (SCL) in Caterpillar, the CVP of IFPEN, Sandia National Laboratories (SNL) [2], Technical University of Eindhoven (TUE) [10], and Michigan Technological University (MTU) [11, 12]. The review in [3] concluded that the ambient gas was relatively homogeneous in CPF compared to CVP. Nevertheless, differences among the CVP were highlighted through the comparison of the bulk temperature Tbulk (defined as the average homogeneous temperature) with the actual measured temperature (called Tcore). Different ratios of Tcore/Tbulk were observed in these facilities and this demonstrated that the heat losses to the walls were different from one vessel to another. These ratios were easily correlated to the gas mixing intensity controlled by the CVP fan speed. The injector nozzle temperature, which is representative of the fuel temperature, was also pointed out in [3] as an important boundary parameter. It has been demonstrated that the liquid length and the soot formation are highly sensitive to this parameter [3, 9, 13]. Thus, many efforts were dedicated in the ECN community to standardize the measurement of the temperature in the injector nozzle. A dummy injector instrumented with a sensitive thermocouple located in the sac area was proposed for this purpose. This tool helped to optimize the cooling of the injectors in the CPF and CVP vessels. As a result, the fuel temperature in these facilities was on the target within ±15 K [3]. In CVP vessels, major and minor products are generated in the post-burn stage. Several studies [3, 13, 14] pointed out the potential effects of these species on the spray ignition, lift-off length and the soot oxidation.

Few studies [1517] have reported ambient velocity measurement in ECN facilities. In spray A target conditions, a near quiescent ambient gas is required and is fulfilled if the global velocity in the vessel is lower than 1 m/s. This threshold has been chosen for two reasons: (i) this value is significantly lower than the velocities of the spray entrained gases, (ii) lower values are hard to achieve in CVP and CPF, otherwise the temperature homogeneity is too deteriorated. Still, recent studies [13, 18, 19] have pointed the sensitivity of the spray to possible different turbulence levels.

ECN-France is an Agence Nationale de la Recherche (ANR) founded project with the goal of involving French research centers (PPRIME institute in Poitiers and PRISME institute in Orleans) to contribute to ECN database. To this end, an extensive campaign has been carried out to characterize their experimental facilities to meet the ECN spray A conditions. IFPEN led this project by providing assistance on the experimental methodologies and by characterizing all the injectors in their well-characterized CVP. Thus, the experimental data from IFPEN have been used as reference data for comparison.

This work has been conducted in two phases: (i) characterize in details the boundary conditions and, (ii) perform measurement at spray A conditions using the reference diagnostics that are standardized in ECN community [4]. The present paper (part 1) presents the results obtained in the first phase, while the second is discussed in a further publication (part 2) [20].

More in details, the current paper reviews the characterization of the boundary conditions in spray A experiments with focus on the ambient temperature and velocity measurements. Particularly, it introduces two new vessels that have been used to study spray A during the ECN France project. The Rapid Compression Machine (RCM) of PPRIME institute and the New One Shot Engine (NOSE) of PRISME institute are presented and compared to the CVP of IFPEN in which all the injectors used were first characterized. The target conditions tested are the reference spray A conditions for inert (O2 0%) and reactive conditions (O2 15% vol) [1]. The boundary conditions characterization will be the base for the discussion of the spray results that will be presented in the part 2 of this paper [20] in an approach similar to that of [4].

In the first part of this paper, we describe the new RCM vessels used for the spray A experimental studies and we compare them to the IFPEN CVP. In a second step, we present Particle Image Velocimetry (PIV) measurements focused on the Start Of Injection (SOI) and we compare the results with the data available in the literature. Later, we review the temperature measurements in ECN facilities with a quick overview on optical methods that can be of potential interest. We also review thin-wire thermocouple thermometry, as it was widely used for the measurement of ambient gas temperature in previous spray A studies. We present the temperature distributions of current vessels and compare them to those of CPF and CVP reported in literature. We also report time-resolved temperature measurement in the near nozzle area during the injection. In the remnant of this work, for the sake of simplicity, we may refer to the experimental facility of each institution by the name of the institution itself (e.g., instead of “in the rapid compression machine at PPRIME”, “in PPRIME”).

2 Experimental setups

In this section, the specific features of the experimental facilities employed by the different institutions are presented. Since the description of IFPEN facility can be found in the literature [2, 3, 8], only a general description is presented in this work. On the other hand, more attention will be given to the PPRIME and PRISME facilities, being their standardization the main objective of this paper.

2.1 IFPEN

The IFPEN CVP vessel has been previously presented and compared to other CPF and CVP vessels in [2, 3]. The layout of the chamber is presented in Figure 1. The vessel operation has already been extensively described in [3] and it is only briefly summarized here. In a pre-burn constant volume vessel, the fuel is injected in a chamber whose pressure, density, and temperature are representative of the thermo-dynamic conditions encountered in the combustion chamber of a real engine at the injection timing. This requires the ability to reach high temperatures and density, which are achieved by the combustion of a flammable mixture, generally referred as “pre-burn”. The experimental methodology can be summarized in four main phases:

  1. Sequential gas filling: the chamber is filled with the pre-burn gases in precise proportions.

  2. Pre-burn: a premixed combustion is triggered by a spark plug to increase the temperature and the pressure of the gas. The product composition (e.g., O2 content) is the experiment target value.

  3. Cool down: the heat flux towards the walls decreases progressively the gas temperature and pressure.

  4. Injection: when the gases reach the target temperature, the fuel is injected.

thumbnail Fig. 1

(a) 3D view and (b) cross section of IFPEN CVP.

It is important to underline that the time scale of the injection (~1.5 ms) is short compared to the cool down time scale (2–3 s). Therefore, the pressure and the temperature can be assumed constant during the injection event. The details of these phases are presented in [2, 8].

2.2. PPRIME

In PPRIME, a flat piston RCM is used for the study of ECN spray A. This device was previously used to study auto-ignition phenomena. More details about its standard configuration can be found in [2124]. The schematic of the RCM is shown in Figure 2. PPRIME RCM performs a single shot experiment in two steps:

  1. Rapid compression of approximately 35 ms (approximately 1050 cm3 displacement and a CR of 9:1).

  2. Constant volume stage of approximately 500 ms.

thumbnail Fig. 2

Schematic of PPRIME institute RCM. (a) General view: (1) hydraulic cylinder, (2) cam, (3) guiding wheel, (4) brake system, (5) return pneumatic cylinder, (6) piston/cylinder, (7) combustion chamber. (b) Zoomed section view on the piston cylinder assembly.

When controlling the initial gas temperature, pressure and composition, a wide range of conditions can be achieved at the end of compression while the volume is kept constant.

The RCM has optical access to the entire volume at TDC without any distortion, thus facilitating the use of optical diagnostics techniques. The ECN injector (#306.14) is centrally mounted on the upper side of the RCM. The same side is also equipped with: (i) a cartridge heater having an integrated thermocouple, (ii) an intake/exhaust valve and (iii) a Kistler 6125 piezoelectric pressure sensor. Three of the four lateral sides are equipped with Sapphire windows. The entire RCM cylinder body is covered with Watlow Silicon rubber heaters. All heaters are controlled to reach an initial gas temperature of 363 K. The injector is connected to a GM 97303659 rail through a 24 cm length high-pressure tube while the fuel return line is connected to a TESCOM 44-1700 backpressure regulator. A Kistler 4067E fuel pressure sensor is installed at 7 cm from the injector inlet. The fuel is pressurized using a Maximator GSF300 hydro-pneumatic pump. The RCM settings used to reach the ECN spray A test conditions are summarized later in Table 1. The flat piston used in the current RCM configuration is known to generate substantial temperature heterogeneities caused by the internal aerodynamics [23, 25]. The piston cylinder motion generates roll-up vortices that mixe the cold boundary layer with the warmer bulk gas. At the end of the compression, the piston stops at the TDC and the vortices get smaller by the effect of viscous dissipation, while making the gas temperature more homogeneous. Thus, the injection is placed later in the post-compression period (approximately 74 ms after TDC) to allow reaching a reasonably homogeneous and quiescent flow as required for the ECN spray A conditions. A typical RCM pressure profile is illustrated in Figure 3. At inert conditions (O2 mol% = 0), the comparison of cases with and without injections shows how the evaporated spray cools down the temperature inside the chamber as the in-cylinder pressure decreases after the injection. At reactive conditions (O2 mol% = 15), the spray combustion induces a significant pressure rise of approximately 10 bar. The relatively small volume of the RCM chamber (131 cm3 at TDC) explains this pressure increase. The ambient gas temperature was measured with thin wires (7.6 μm) type K thermocouples and using a pre-calibrated acquisition system (QuickDAQ DT9805). According to the hardware manufacturer, this device is capable of recording the junction temperature at a sampling rate of 10 kHz with an accuracy of 0.2 K. The thermocouple wires used have a high purity level that meet the type K standard, which has an error level of ±1 K. Further details about the measurement technique and the calculation of the gas temperature are reported later in Section 4.2.4.

thumbnail Fig. 3

Typical pressure profile of inert and reactive conditions at spray A conditions in PPRIME.

Table 1

Comparison of IFPEN, PPRIME, and PRISME vessels when set for spray A experiments.

2.3 PRISME

The experimental apparatus at PRISME is referred to as New One Shot Engine (NOSE). It was developed based on the Rapid Cycling Machine (RCYM) concept. This facility is designed to study the combustion at conditions relevant to recent Compression Ignition (CI) engines. A single-cylinder low-speed diesel engine (standard power: 13.4 kW at 750 rpm) is used as shown in Figure 4a. The engine displacement volume is 3354 cm3 (the bore and stroke are respectively 155 mm and 177.8 mm) and the CR is 15:1. The original cylinder head is replaced with a customized optical chamber equipped with four quartz windows (25 mm thick, 25 mm wide, and 80 mm long). Different cross-section views of this chamber are presented in Figures 4b4d. The chamber has a rectangular prism shape of 44 × 44 × 90 mm length (240 cm3 total volume). The ECN injector (#306.19 [8]) is centrally mounted on the upper side. The piston head is designed to achieve a near-quiescent gas condition at the top dead center, by optimizing the roof angle, which is set to reduce the squish effect. The engine coolant flow and temperature are set to control the initial temperature before compression. The optical cylinder head is equipped with four cartridge heaters placed in the four corners and a type K thermocouple is used to monitor the wall temperature. This system allows controlling the injector nozzle temperature. The NOSE combustion chamber has four gas connections: intake, air, vacuum and exhaust as illustrated in Figure 4. The intake connection is used to supply the vessel with working gases. Up to four gases (usually N2, O2, NO, and CO2) can be supplied while controlling the amount of each gas. Mass Flow Controllers MFC (Brooks 5860S combined with 5850S) are used for this purpose. A liquid MFC is also installed on the same line, which allows injecting liquids in the chamber. At the end of the experiment, the gases are released through the exhaust connection. Then, the chamber is purged with air through a dedicated connection for at least 10 s. Before starting a new experiment, the valves are closed and a vacuum pump is used to empty the chamber from all residual gases. The optical cylinder head is equipped with a piezoelectric pressure sensor (Kistler 7001). A National Instrument (NI) Compact RIO system is used to control the setup and to acquire the parameters of interest. Particularly, the in-cylinder pressure and temperature, the injector current and the encoder crank angle pulses are logged at 250 kHz. PRISME-NOSE is driven by a brushless (permanent magnet rotor) DC electric motor (Phase-Automation© U31340). The torque and engine speed are controlled with a dedicated electric driver (Phase-Automation© axM300-400) through the LabVIEW interface of the NOSE setup. First, a reverse mode is activated and the piston is moved from 180 CAD before TDC to 280 CAD before TDC. The system is then switched to velocity feedback mode based on a closed loop control. The high accuracy engine speed measured with the encoder is used to impose an optimized piston velocity profile (rapid compression < 40 ms and steady state at TDC 10 ms). This method reduced the piston static friction and smoothed the piston motion. It also provides a satisfying repeatability of the piston velocity profiles. A comparison of the target and measured piston velocity is presented in Figure 5 demonstrating the effectiveness of the closed loop feedback control-mode. To reach the ECN spray A target conditions, the PRISME-NOSE is operated with a coolant temperature set at 356 K and an initial pressure of 1.8 bar (for inert experiments with N2). The optical cylinder head is heated up to reach 363 K in the injector nozzle region. Figure 6a illustrates the pressure, temperature, and density profiles obtained with these initial conditions. Figure 6b demonstrates the good repeatability of experiments at the SOI. More details about the NOSE can be found in [26, 27]. The time-resolved gas temperature is measured with thin wires (12.7 μm) type K thermocouples and using an AD595CQ amplifier coupled to NI Analog Input module. Further details about the measurement technique are reported later in Section 4.2.5.

thumbnail Fig. 4

The NOSE chamber: (a) the overview of real image, (b) the cross-section view in middle plane of chamber, (c) the top-view without chamber head, and (d) the cross-section view A–A without chamber head.

thumbnail Fig. 5

NOSE operation: comparison of the feedback and command velocities as a function of time referenced to the initial position (280 CAD before TDC).

thumbnail Fig. 6

Thermodynamic condition of NOSE operation: (a) an example of pressure, temperature and density temporal evolution, (b) pressure and temperature evolutions from −20 ms to 20 ms after TDC for 10 repeated tests.

2.4 Comparison

The main features of the experimental facilities together with the settings for spray A experiments are compared in Table 1. Unlike CVP, PPRIME, and PRISME RCM reach the spray A target conditions through volumetric compression without the need for a pre-burn event. Therefore, the gas composition is relatively simpler in these devices. In PPRIME RCM, the lowest available CR is used to maximize the volume at the SOI (131 cm3). A controlled amount of Argon is added to N2 and O2 mixture to adjust the gas density and the heat capacity ratio, which enable reaching the target conditions. In PRISME NOSE, the compression ratio has been adapted to reach these conditions by using only the target ratio of O2 and N2. For reference, the heat capacity Cp at spray A conditions is reported in Table 1, since this parameter has a primary role in the resulting temperature of the fuel-gas mixture and therefore it might have a non-negligible impact on the spray liquid length, ignition delay and lift off length. Similarly to PPRIME RCM, the volume of PRISME NOSE chamber at the time of injection (240 cm3) is significantly smaller than the one used in IFPEN CVP. The RCM devices have an order of magnitude smaller volume and their wall/body temperatures are 100 K lower compared to IFPEN, while the injector nozzles are protruding similarly.

In the early ECN studies [2, 3], the spray experiments are always performed after a proper characterization of the boundary conditions. This step, often referred to as the “standardization”, assures the accuracy of boundary conditions by: (i) providing a detailed list of conditions, which are not often considered in other studies (e.g., fuel temperature, gas velocities, etc.), (ii) sharing and discussing the characterization methodologies employed, together with the related uncertainties. In this work, the attention is focused on the characterization of three main boundary conditions:

  • The gas velocity (this will be detailed in Sect. 3).

  • The injector nozzle temperature (this will be detailed in Sect. 4.1).

  • The gas temperature (this will be detailed in Sect. 4.2).

In the following, the results related to these characterizations will be presented and discussed.

3 Characterization of the gas velocity

3.1 PIV optical setup

3.1.1 IFPEN

The PIV setup used in IFPEN is detailed in [28], where inert experiments were performed. A 10 kHz dual cavity YLF laser was used to illuminate the seeded combustion chamber. A Photron SA1 CMOS camera equipped with a Nikkor 50 mm f/1.2 lens was used to acquire the images at 20 kHz. The laser pulses had a 2 μs time delay, thus this time constant was used to calculate high velocity levels. Lower velocities were calculated using the 100 μs time step. The ROI was centered on the spray axis and has 56 mm width and 74 mm long (image resolution 8 mm/pixel). A sketch of the layout is presented in Figure 7.

thumbnail Fig. 7

Schematic of the PIV setup in IFPEN.

3.1.2 PPRIME

In PPRIME, high-speed 2D PIV technique was performed to study the flow characteristics of the RCM at test conditions similar to those set for spray A experiments. A Mesa PIV Nd-Yag laser has been used to generate laser pulses at 10 kHz. The laser sheet was placed in the RCM center plane (coplanar to spray axis) and the chamber was seeded with Silicone oil droplets. A schematic of the PIV setup is presented in Figure 8 (more details are reported in [24]). The RCM settings are summarized in Table 2. The CR and the target density are similar to those used for spray A experiments (see Tab. 1). However, a lower target compression temperature is used to avoid the evaporation of the seeded oil droplets. This has a minor influence on the internal aerodynamics of the RCM. The characterization of the flow is focused on the injection period (tinj ± 10 ms). The velocity vectors are calculated with a time series pyramid sum of correlation algorithm [29] and using a decreasing window size of 128 × 128 – 32 × 32 with 50% overlap (vector density 1.28 per mm). The turbulence level was estimated based on a time averaging of 198 velocity fields and following the equations described in the next Section 3.2.

thumbnail Fig. 8

Schematic of the PIV setup in PPRIME.

Table 2

RCM test conditions for PIV experiments.

3.1.3 PRISME

In PRISME, similar high-speed PIV measurements were performed using a dual Hawk-HP Nd–Yag laser (9.4 mJ/pulse, 532 nm at 2.5 kHz) and a Phantom CMOS camera. The NOSE chamber was also seeded with Silicon oil droplets and the laser sheet was placed in the center plane (an illustration is presented in Fig. 9). The velocity vectors are calculated using again a time series pyramid sum of correlation algorithm but using a decreasing window size of 64 × 64 – 16 × 16 with 50% overlap.

thumbnail Fig. 9

Schematic of PIV setup in PRISME.

3.2 Velocity and turbulence at spray A conditions

In PPRIME RCM and PRISME NOSE, the 2D PIV measurements were used to estimate the turbulence levels at the SOI. The calculations are based on time averaging in the range SOI ± 10 ms (using 198 velocity fields in PPRIME and 50 in PRISME). In IFPEN CVP, the 2D PIV measurements were more focused on the aerodynamics of the spray [28]. Only a limited number of flow fields are available without injection. Thus, the turbulence was estimated based on the ensemble averaging and using 73 velocity fields measured at the SOI.

The following equations are used for the calculation of the average velocity, the RMS and the turbulent kinetic energy TKE:V¯(x,y)=1Nt=1NV(x,y,t),$$ \bar{V}\left(x,y\right)=\frac{1}{N}\sum_{t=1}^NV(x,y,t), $$(1)σ(x,y)=1N-1t=1N(V(x,y,t)-V¯(x,y))2,$$ \sigma \left(x,y\right)=\sqrt{\frac{1}{N-1}\sum_{t=1}^N{\left(V\left(x,y,t\right)-\bar{V}\left(x,y\right)\right)}^2}, $$(2)TKE=12(σ(x)2+σ(y)2+(σ(x)+σ(y)2)2).$$ \mathrm{TKE}=\frac{1}{2}\left({\sigma (x)}^2+{\sigma (y)}^2+{\left(\frac{\sigma (x)+\sigma (y)}{2}\right)}^2\right). $$(3)

The PIV results from PPRIME RCM, PRISME NOSE, and IFPEN CVP are presented in Figures 1014 and summarized in Table 3, where spatially averaged data are compared to results reported in [30].

thumbnail Fig. 10

Average and standard deviation of velocity magnitude at time of injection (IFPEN CVP using ensemble average).

thumbnail Fig. 11

Average and turbulent kinetic energy at time of injection (IFPEN CVP using ensemble average).

thumbnail Fig. 12

(a) Average and (b) standard deviation of velocity magnitude at time of injection ±10 ms (PPRIME).

thumbnail Fig. 13

(a) Average and (b) turbulent kinetic energy at time of injection ±10 ms (PPRIME).

thumbnail Fig. 14

(a) Average velocity magnitude field and (b) turbulent kinetic energy field at SOI (PRISME).

Table 3

Comparison of turbulence characteristics: Sandia CVP and Caterpillar CPV reported in [17, 30].

The data from all ECN facilities reported in Table 3 show that they comply with the spray A ambient gas near-quiescent condition. The standard condition implies that the maximum and average velocity at the SOI is lower than 1 m/s. However, different levels of velocities and turbulence are observed: the values in PPRIME RCM and Caterpillar CPF are an order of magnitude higher than in Sandia and IFPEN CVP. It is worth reminding that the turbulence level is dependent on the fan speed in CVP vessels. Higher fan speed would generate more temperature homogeneity but too high turbulence. As it was highlighted in [13], small differences in the turbulence level and the nozzle diameter between different ECN facilities are expected to cause a vapor distribution significantly different. Thus, the early stages of the fuel-air mixing might be affected and lead to different ignition behaviors. This will be further investigated in [20], where spray behavior in current facilities will be compared. Although all facilities comply with the near quiescent condition, the estimated TKE is significantly different from one vessel to another. This also encourages considering a more realistic velocity and turbulence level when simulating the spray A.

4 Characterization of temperature boundary conditions

4.1 Injector nozzle temperature

The injector nozzle temperature is an important parameter for spray characterization. It is considered highly representative of fuel temperature, which has a first order influence on spray evaporation and auto-ignition [9]. At spray A reference conditions, it is estimated that all the injected fuel, for injections up to 4 ms long, is originally contained in a distance lower than 4 mm from the nozzle sac [3, 9]. Therefore, considering the long time between two consecutive injections (>60 s) it can be considered that the fuel and the injector nozzle reach the same temperature. In current study, the injector nozzle temperature is measured with a dummy injector instrumented with a type K thermocouple. The same device was used in previous ECN characterizations [2, 3, 9]. At IFPEN, the CVP walls are heated up to ~473 K and thus, the injector is water cooled down to the reference 363 K. In PPRIME RCM and PRISME NOSE, the walls are heated respectively to 363 K and 373 K. This ensures that the injector nozzle temperature is on the target standard value. In Figure 15, the average nozzle temperatures measured with the dummy injector during a spray A experiment in IFPEN CVP and PPRIME RCM are compared. At the start of experiment, the temperature is well set on the target value. However, the pre-burn event in the CVP and the rapid compression of the RCM induce an increase of the nozzle temperature. At the instant of injection, the nozzle temperature on the CVP increases, up to 10 K higher than at the start of experiment. The increase is less pronounced on the RCM (approximately 1 K). Considering also the measurement performed in [2], it is reasonable to expect that the short time scale of the pre-combustion event does not affect significantly the fuel and nozzle temperature (<20 K). However, despite the small fuel volume interested by the injection, we cannot exclude some fuel temperature gradient during the injection. According to what is reported in [9], 10 K difference in the fuel temperature could cause a difference in liquid length of ~0.3 mm and 0.1 mm in lift-off length. The effects of boundary conditions on the spray characteristics will be further discussed in [20].

thumbnail Fig. 15

Comparison of the injector nozzle temperature in IFPEN CVP and PPRIME RCM. The time is referenced to the SOI.

4.2 Gas temperature measurements

4.2.1 General methodology

Many optical diagnostic techniques could be used for local temperature measurements at ECN test conditions. A brief review of these methods is presented in Table 4 [31]. In RCMs and CPFs, where oxygen free gases can be used, Laser Induced Fluorescence (LIF) with different tracers (acetone [32], toluene [33], or anisole [34]) could be considered. LIF methods can provide spatially resolved information, which would facilitate the characterization of the combustion vessels used for ECN spray studies. However, it could have limitation at relatively high temperature, as the tracer pyrolysis should be considered when calculating the temperature from the fluorescence signal. In addition, thermal quenching significantly reduces florescence signal at high temperatures (e.g., 900 K), increasing measurement uncertainties. Laser Rayleigh scattering [3537] is commonly used in gas phase thermometry. However, this technique can be limited by the background interference, which causes a too low signal to noise ratio. Laser absorption thermometry [38, 39] is relatively easy to install and can be time-resolved [40]. The measurement is performed in the line of sight, which makes the method not suitable for the assessment of temperature gradients. Recently, Laser Induced Grating Spectroscopy (LIGS) was successfully performed to measure temperature at ECN like conditions [41, 42]. This method measures local temperature with a relatively good accuracy and can be time resolved. Thin-wire thermocouples are easier to install and can perform similar measurements. However, they are more intrusive compared to LIGS.

Table 4

Summary of optical diagnostic techniques that could be used in ECN spray vessels [31].

4.2.2 Thin-wire thermocouple thermometry

Thermocouple thermometry is often used for the study of temperature fluctuations in turbulent flows [46, 47]. The proper design of these sensors and the use of extremely thin wires enable fast response with a relatively good accuracy. The study of their time response is usually performed based on the theory of frequency response [4851]. To have a fast response thermocouple, the sensor should be bare-bead with junction on wires as small as possible. It is recommended to ground isolate the junction to reduce the noise potentially caused by Electro-Magnetic Interference (EMI). The welding of the junction on extremely small wires can be achieved using capacitive discharge technique after an optimization of the discharge parameters. Recommended settings for this method are described in [52]. Thermocouples measure the temperature of the junction based on the Seebeck effect. To measure gas temperature, the heat exchange around the sensor must be modeled [47, 53]. Assuming energy conservation on the junction, the gas temperature Tg and the junction temperature Tj can be correlated. The main heat exchanges on the junction are through: (i) convection between gas and junction, (ii) radiation between walls and junction, (iii) conduction between junction and the prongs. The heat transfer by dissipation is often neglected at relatively low flow velocities. The radiation between the junction and the gas is usually negligible compared to the radiation to the wall. Catalytic reactions can occur between the sensor wires and the gas. In such cases, the heat transfer would depend on the main chemical reactions involved in this phenomenon. Catalytic reactions make the sensor more intrusive especially in unstable gas environments. This was the case when using type R thermocouples in some ECN CVP vessels [2]. The sensor wires can be coated to avoid their catalytic effect, as proposed in [47, 54, 55]. However, this would modify the convective and radiative heat transfers. Based on the previous assumptions, the gas temperature can be modeled as follows:Tg=Tj+(τh+τc)dTjdt+σεh(T14-Ts4),$$ {T}_{\mathrm{g}}={T}_{\mathrm{j}}+\left({\tau }_{\mathrm{h}}+{\tau }_{\mathrm{c}}\right)\frac{{\mathrm{d}T}_{\mathrm{j}}}{\mathrm{d}t}+\frac{{\sigma \epsilon }}{h}\left({T}_1^4-{T}_{\mathrm{s}}^4\right), $$(4)where: σ is the Stefan-Boltzmann constant, ε is the temperature dependent emissivity (typical value for thermocouple wire can be found in [56]), h is the convection heat-transfer coefficient, Ts is the temperature of wall surfaces surrounding the junction, and τh and τc are respectively the time constant of the convective heat-transfer contribution and the time constant of the conduction heat-transfer contribution.

In the early studies on cold wires thermometry [57, 58], a cooled length from conduction lc has been defined to describe the characteristic length of the conduction heat transfer in these applications. Applying this theory to thin-wire thermocouples would imply a design condition that will reduce substantially the errors from the conduction losses. It has been demonstrated in [47] that conduction losses from the junction to the prongs through thin-wires can be neglected if L/lc > 10. Where L is the length of the wires between the two prongs and lc is defined as follows:lc=dw2λwλg Nu,$$ {l}_{\mathrm{c}}=\frac{{d}_{\mathrm{w}}}{2}\sqrt{\frac{{\lambda }_{\mathrm{w}}}{{\lambda }_{\mathrm{g}}\enspace \mathrm{Nu}}}, $$(5)where dw is the wire diameter, λw and λg are respectively the wire and the gas thermal conductivity, Nu is the Nusselt number. For typical sensors and gas flow applications, this condition is equivalent to L/dw ≥ 400. With this design condition, the gas temperature can be calculated as:Tg=Tj+τhdTjdt+σεh(Tj4-Ts4),$$ {T}_{\mathrm{g}}={T}_{\mathrm{j}}+{\tau }_{\mathrm{h}}\frac{{\mathrm{d}T}_{\mathrm{j}}}{\mathrm{d}t}+\frac{{\sigma \epsilon }}{h}\left({T}_{\mathrm{j}}^4-{T}_{\mathrm{s}}^4\right), $$(6)where τh$ {\tau }_{\mathrm{h}}$ is:τh=ρj Cj dj4 h,$$ {\tau }_{\mathrm{h}}=\frac{{\rho }_{\mathrm{j}\enspace }{C}_{\mathrm{j}\enspace }{d}_{\mathrm{j}}}{4\enspace h}, $$(7)ρj is the thermocouple junction density, cj is the junction specific heat, dj is the junction diameter and h is the convective heat-transfer coefficient, which can be calculated from the Nusselt number:Nu=h djλg.$$ \mathrm{Nu}=\frac{h\enspace {d}_{\mathrm{j}}}{{\lambda }_{\mathrm{g}}}. $$(8)

Many Nusselt correlations can be found in the literature. Usually, the correlations depend on the flow configuration and the Reynolds (Re) and Prandtl (Pr) numbers. Collis and Williams [59] proposed for 0.02 < Re < 44:Nu=(0.24+0.51 Re0.45)(TfTg)0.17,$$ \mathrm{Nu}=\left(0.24+0.51\enspace {\mathrm{Re}}^{0.45}\right){\left(\frac{{T}_{\mathrm{f}}}{{T}_{\mathrm{g}}}\right)}^{0.17}, $$(9)where Tf is the film temperature defined by Collis and Williams [60] as (Tg + Tj)/2. Using the same definition of the reference temperature for the fluid physical properties, Kramers [61] proposed a generalized correlation:Nu=2+1.3 Pr0.15+0.66 Pr0.31Re0.5,$$ \mathrm{Nu}=2+1.3{\enspace \mathrm{Pr}}^{0.15}+0.66\enspace {\mathrm{Pr}}^{0.31}{\mathrm{Re}}^{0.5}, $$(10)and a specific one for the configurations where 0.01 < Re < 10 000 and 0.71 < Pr < 1000 [61]:Nu=0.42 Pr0.2+0.57 Pr0.33Re0.5.$$ \mathrm{Nu}=0.42{\enspace \mathrm{Pr}}^{0.2}+0.57\enspace {\mathrm{Pr}}^{0.33}{\mathrm{Re}}^{0.5}. $$(11)

Churchill and Bernstein [62] proposed a generalized correlation for forced convection heat transfer across a cylinder, which was validated on a wide range of conditions (Re Pr > 0.2):Nu=0.3+0.62 Pr1/3Re1/2[1+(0.4/Pr)2/3]1/4[1+(Re282 000)5/8]4/5.$$ \mathrm{Nu}=0.3+\frac{0.62{\enspace \mathrm{Pr}}^{1/3}{\mathrm{Re}}^{1/2}}{{\left[1+{(0.4/\mathrm{Pr})}^{2/3}\right]}^{1/4}}{\left[1+{\left(\frac{\mathrm{Re}}{282\mathrm{\enspace }000}\right)}^{5/8}\right]}^{4/5}. $$(12)

In the case of thin-wire thermocouples designed as described above, Re and Pr numbers are defined as follows:Re=U djνg,$$ \mathrm{Re}=\frac{U\enspace {d}_{\mathrm{j}}}{{\nu }_{\mathrm{g}}}, $$(13)Pr=μg Cpλg,$$ \mathrm{Pr}=\frac{{\mu }_{\mathrm{g}}\enspace {C}_{\mathrm{p}}}{{\lambda }_{\mathrm{g}}}, $$(14)where U is the gas velocity, νg is the kinematic viscosity, μg is the gas dynamic viscosity and Cp is the gas specific heat capacity. Thus, the calculation of the gas temperature from the measured junction temperature is highly affected by the wire diameter dw, the junction size dj/dw, the local gas velocity U and the Nu correlation.

In PPRIME, the thermocouples used are bare-bead and ground isolated. The junction is welded on a 7.6 μm type K wires, which are mounted in a strained position between similar wires of 200 μm diameter. An illustration of these sensors is presented in Figure 16. The distance between the two prongs is approximately 400 times greater than the wire diameter. In most cases, the junction is perfectly welded and its size is as small as the wires, as demonstrated in Figure 16b. In few cases, the junction has a spherical shape and is slightly bigger than the wires, as illustrated in Figure 16c. A test run has been performed in PPRIME at spray A target conditions, using simultaneously two thermocouples with wires diameters of 7.6 μm and 12.7 μm and mounted in the center of the RCM chamber. The comparison of the corrected temperature is presented in Figure 17 to investigate the main parameters contributing in the measurement error. The effects of dw, dj, U, and Nu correlations are highlighted here. In default error calculation, the velocity considered is time-resolved and estimated from PIV measurements and the Churchill correlation for Nu (Eq. (12)) is used. The comparison of the correction (TgTj) in Figure 17a is focused on the compression step where it is demonstrated that the temperature is homogeneous in both thermocouples locations [23, 25]. The correction is approximately twice higher when using the 12.7 μm thermocouple. In applications where temperature fluctuates significantly, it is expected to have higher corrections when using thicker wires. The use of smaller thermocouples substantially reduces the need for correction and improves the time response of the sensor. It is also important to consider the effective size of the thermocouple junction. For instance, this was the case for the welding presented in Figure 16c. In Figure 17b, it is demonstrated that an under-estimation of the junction size would result in an under-estimation of the correction and thus a systematic error on the estimation of the temperature. One of the difficulties of thermocouple thermometry is the proper estimation of the fluid velocity across the junction. In Figure 17c, three estimations of the gas velocity have been used for the calculation of temperature correction for a 7.6 μm thermocouple diameter: “U reference” a time resolved estimation based on PIV measurements in the location of the sensor; “U = 3.5 m/s” a constant velocity whose level is representative of the maximum velocity in the compression step; and “U = 0.2 m/s” a constant velocity whose level is representative of average velocity in the late post-compression step. It is demonstrated that an over-estimation of the gas velocity would lead to an under-estimation of the error correction. Conversely, an under-estimation of the velocity would result in a substantial increase of the error correction and thus an over-estimation of the temperature fluctuations. It is then very important to estimate properly the velocity across the sensor wires as this parameter highly influences the convective heat transfer coefficient h. The Nu correlation choice is also important for the estimation of this coefficient. The error corrections from the commonly used correlations in thermocouple thermometry are presented in Figure 17d. It is demonstrated that the use of the Kramers generalized correlation (Eq. (10)) would lead to the under-estimation of the temperature fluctuation in this case. All the other correlations have similar temperature correction. In configurations where the junction size is as small as the wires, it is more suitable to use the Churchill correlation (Eq. (13)).

thumbnail Fig. 16

Schematic of the 7.6 μm thermocouple sensor (PPRIME institute): (a) sensor prongs mounting, (b) zoom on the thermocouple junction welded successfully, (c) zoom on junction having a spherical shape.

thumbnail Fig. 17

Main parameters affecting the measured temperature correction: (a) effect of dw, (b) effect of dj/dw, (c) effect of the gas velocity U (dw = 7.6 μm), (d) effect of the Nu correlation (dw = 7.6 μm). PPRIME RCM test at spray A target conditions (time axis is referenced to the end of compression).

4.2.3 Previous gas temperature measurement in ECN facilities

In most ECN studies, thin-wire thermocouples were used for the measurement of ambient gas temperature. Mainly, two types of thermocouples were used (R and K). Type R thermocouples are easily available and are able to measure high temperatures (up to 1873 K), which is particularly needed in CVP vessels. Nevertheless, the catalytic effect of type R wires (87% Pt/13% Rh–Pt) have limited their use as uncontrolled pre-burn events may occur, particularly when the reactant gases are introduced sequentially in the CVP chamber [2]. Type K thermocouples, despite a lower maximum measured temperature (approximately 1520 K), are widely used in many applications and were suitable for measuring gas temperature in RCM and CPF. Different methods for “temperature correction” have been used in the ECN community and many recommendations have been highlighted to improve the measurement accuracy [2, 3]. Table 5 summarizes the different thermocouples and correction methods used in previous ECN studies. In CPF vessels, sheathed thermocouples were used since this type of facility provides steady thermodynamic conditions within the chamber. This configuration of the sensors is not suitable for measuring fast temperature fluctuation as the heat transfer by conduction has a relatively long time constant. Covering the thermocouple junction makes the sensor more robust, while it deteriorates its spatial resolution and time response.

Table 5

Summary of temperature measurement with thermocouples under spray A conditions [2, 3].

4.2.4 Gas temperature measurement in PPRIME

In PPRIME, the temperature is simultaneously measured at three different locations downstream to the injector nozzle as illustrated in Figures 18a and 18b (thermocouple 1–3). Complementary measurements have been performed in the near nozzle area using a further thermocouple (thermocouple 4, Fig. 18c). All thermocouples have a 7.6 μm wires and have been designed as recommended in Section 4.2.2. The temperature correction is calculated using equation (6) with a time-resolved velocity from PIV measurements and the Nu correlation of equation (13). An example of the corrected gas temperatures is presented in Figure 19. The measured temperature is compared to the adiabatic core temperature Tc, using the following equation:T0Tcγγ-1dTT=ln(PcP0),$$ {\int }_{{T}_0}^{{T}_{\mathrm{c}}}\frac{\gamma }{\gamma -1}\frac{\mathrm{d}T}{T}={ln}\left(\frac{{P}_{\mathrm{c}}}{{P}_0}\right), $$(15)where T0 is the initial temperature before compression, γ is the heat capacity ratio, P0 and Pc are respectively the initial and the compression pressures. This calculated temperature is representative of the maximum temperature in the RCM chamber [23, 25]. The local instantaneous measurements show excellent agreement with this temperature during the early stage of compression. In this time range, all sensors are located in the hot adiabatically compressed gas region. This demonstrates that the thermocouples have a relatively fast response and that the correction method is suitable for the current application. At later stages, the cold roll-up vortices generated by the piston motion move toward the upper side of the chamber leading to substantial temperature heterogeneities. As the piston is kept at TDC at the end of the compression, the vortices motion dissipates progressively while mixing the hot and colder gases. Thus, the injection timing has been placed at 74 ms after TDC, when the overall velocity level in the chamber is lower than 1 m/s and the temperature distribution is relatively homogeneous. The averaged gas temperature from five repeated tests measured in 10 different locations is 925 K. The standard deviation of the average temperature through the different positions is 18 K.

thumbnail Fig. 18

(a) Illustration of thermocouples positions, (b) picture of thermocouples mounted in the center plane of the PPRIME RCM, (c) position of the thermocouple in the near nozzle measurement tests.

thumbnail Fig. 19

Illustration of the measured temperature from a single test compared to adiabatic core temperature.

4.2.5 Gas temperature measurement in PRISME

In PRISME, four type K thin-wire thermocouples have been installed in the optical cylinder head of NOSE. The ambient gas temperature is measured at different locations as illustrated in Figures 20a and 20b, where the spatial coordinates are referenced to the injector tip position. The thermocouple wires have a nominal diameter of 13 μm approximately. The junction has a spherical shape of approximately 39 μm diameter, a picture of the bead is reported in Figure 20c. The temperature correction is calculated with equation (6) and the Nu correlation of Kramers equation (10) using an average constant velocity (U = 0.1 m/s, estimation based on PIV measurements presented in the Sect. 3.2). An example of ambient temperature measurement with thermocouple no. 2 is illustrated in Figure 21. The correction of measured temperature is significantly higher during the compression and later in the expansion stage. Most likely, this is explained by the important temperature rise/decay rate in these periods. The average and standard deviation of the measured temperature from 10 tests is presented in Figure 21b. Up to 5 ms after TDC, the temperature is stable. As the temperature start, to decrease, more heterogeneities are observed. In this period of interest (10 ms after TDC), the pressure and density decrease slightly, from 59.8 bar to 59.3 bar and from 22.8 kg/m3 to 21.7 kg/m3. The temperature fluctuation in this period is around 1.6% of the target value.

thumbnail Fig. 20

Location of the temperature measurement in PRISME. Four thin-wire thermocouples of 13 μm nominal diameter are used: (a) illustration of measurement positions, (b) picture of thermocouples mounted in the centerline of the NOSE chamber, and (c) picture of the thermocouple junction.

thumbnail Fig. 21

Illustration of typical temperature profile in PRISME (thermocouple 2): (a) example of density, pressure, and temperature profiles around TDC, (b) average temperature profile at 900 K target condition, shaded-error bar represents the standard deviation of 10 measurements.

4.2.6 Time resolved temperature measurement in the near-nozzle area

The temperature heterogeneities observed in CVP and RCM vessels are mainly generated by the heat transfer to the walls. In most cases, the injector nozzle protrudes few millimeters inside the vessels and this makes the temperature in the near wall area very important for the interpretation of the spray characterization results.

A thin-wire thermocouple, labeled thermocouple 4 in Figure 18, was placed at 0.3 mm from the injector nozzle in PPRIME. The temperature is measured during the injection under inert and reactive conditions and compared to temperature measurement without injection (compression test). In inert experiments, the exact position of thermocouple 4 is (0.3, 0.3) on the (x, y) coordinates illustrated in Figure 18a. In reactive experiments, the exact position of thermocouple 4 is (3.5, 0.3) as the sensor has been pulled away from the nozzle to avoid damaging it during the combustion at the EOI. The comparison of these measurements is presented in Figure 22, where the theoretical maximum temperature (obtained from the pressure measurement and the adiabatic core hypothesis) is also reported. Figure 22a shows that, in the inert case, the temperature significantly decreases immediately after the start of the injection (in average 50 K lower at the EOI). Considering that the thermocouple is not directly faced to the spray flow and that the temperature decrease is significantly larger than the one indicated by the pressure-derived calculation, it can be concluded that the thermocouple is measuring a local effect. This effect is most likely caused by the colder gases from the boundary layer, which move towards the spray region following the gas-entrainment up to the EOI. The increasing difference between the cases with and without injection observed 20–40 ms after SOI in inert cases is related to the evaporated fuel, which mixes with the compressed hot gases and generates further temperature heterogeneities that are transported and diffused to the entire RCM chamber. The temperature of the fuel-air mixture after the complete mixing (~40 ms after SOI) is approximately 200 K lower than the compression test case.

thumbnail Fig. 22

Temperature measurement in the near nozzle region: (a) at inert conditions, (b) at reactive conditions.

In reactive case, the fuel combustion induces a significant increase of the in-cylinder pressure and temperature. This is clearly demonstrated in Figure 22b where the measured temperature increases substantially after the spray ignition. A similar variation is also observed on the pressure-derived temperature. Due to the short ignition delay time and the position of thermocouple 4 in reactive cases, it is not possible to observe the cooling effect of the spray gas-entrainment. Conversely, the temperature of the entrained air is increasing progressively caused by the combustion event up to approximately the EOI. The spray combustion generates further temperature heterogeneities as the hot combustion products are mixed with air and are transported to the entire RCM chamber, which explains measuring higher temperature fluctuations at approximately 40 ms after SOI. These additional heterogeneities have been observed on Schlieren images (see Fig. A1), which consolidate the above interpretation and explain why the maximum measured temperature is higher than the pressure-derived adiabatic compression temperature at these instants.

ECN vessels have shown different levels of temperature heterogeneities. The standard deviation of the overall measured temperature is around 5% [30]. These non-homogeneities are often not considered in the spray modeling and simulation. So far, more efforts have been spent by the community on the development of adequate chemical kinetic models for the spray. The above temperature measurements during the injection demonstrate that the entrained gas temperature is highly dependent on these heterogeneities. This would influence the spray behavior, particularly the ignition delay and the lift-off length, which are sensitive to ambient gas temperature [63]. Thus, it would be more suitable to consider the specific temperature distribution of a vessel in the simulation works aiming to validate new models against experimental data.

4.2.7 Comparison of temperature homogeneity in ECN spray A vessels

The characterization of the boundary conditions in different ECN facilities has been reported in [2, 3]. The ambient gas temperature was pointed out as one of the key parameters that influences the spray behavior, particularly the ignition delay and the lift-off length. The reported temperature heterogeneities in [2, 3] is reviewed in Table 6 and compared to RCM results from the present study. In CPF and CVP, lower temperature levels are observed in the near nozzle region. This is illustrated in Figure 23a for IFPEN and Sandia CVP vessels; the temperature is approximately 6% lower at the nozzle location compared to downstream in the spray axis. RCMs data from PPRIME and PRISME show relatively more homogeneous levels through the spray axis. In PPRIME particular case, the combustion chamber is small (the piston is located at approximately 45 mm from the injector). Lower temperatures and higher heterogeneities (high test-to-test standard deviation) are observed downstream of the spray axis. In the transverse direction (Fig. 23b), significant temperature gradients are observed in CVP vessels and this is caused by the buoyancy effect [2]. This phenomenon is not observed in RCM vessels as the time constant needed to reach the spray A conditions is an order of magnitude shorter. In addition, the injectors are mounted in the upper side of the RCM chambers. Nevertheless, a small gradient is observed, as the temperature is lower near to the lateral walls, which are located at approximately 25 mm from the spray axis. Thus, in RCM devices, where the flow pattern is usually axisymmetric, the temperature is higher in the core region. The assessment of temperature heterogeneity through the standard deviation of measured temperature at SOI shows higher levels in the RCM of PPRIME. Most likely, the thermocouples used in PPRIME have a faster time response and thus they are more sensitive to temperature fluctuations. In addition, PPRIME RCM has the lowest wall/body temperature, compared to the other ECN facilities presented in Table 6.

thumbnail Fig. 23

Comparison of the averaged normalized temperature measured in different ECN facilities: (a) temperature heterogeneity in the spray axis, (b) temperature heterogeneity in the transverse axis at 40 mm from the injector nozzle. IFPEN and Sandia data from [3].

Table 6

Temperature heterogeneities levels at spray A conditions.

In most ECN facilities used for spray A experiments, the ambient temperature has been measured with thermocouples. The time response of these sensors is a key parameter that reflects how good the temperature fluctuations can be detected. The sensors should be as fast as possible to capture properly the temperature heterogeneity level in each vessel. Suitable thermocouples would be made of thin-wires (as small as possible) and should have a bare junction. The application of temperature correction is needed to estimate the gas temperature from the junction temperature and thus consider the time response of the sensor. In the current work, it has been demonstrated that a careful design and a proper calculation of the correction factors can provide a relatively accurate measurement. Sheathed thermocouples will not have a time response sufficiently short for temperature fluctuation measurements. The corrected temperature from such sensors should consider the conduction of heat through the sensor body and this reduces substantially the spatial resolution of the measurement. Despite the use of different sensors, all spray A vessels displayed a certain temperature heterogeneity; in region occupied by the spray axis, the temperature deviates of maximum 6% from the average temperature as a consequence of thermal stratification. Nevertheless, each facility has its own specific heterogeneity level and pattern. It can also be observed that RCM facilities have smaller heterogeneities compared to CVP vessels. This fact is probably related to the short time elapsed between the gas compression and the injection event and the lower maximum temperature. While in CVP vessels, approximately 700 ms separate the pre-burn event from the SOI. This delay is needed to cool the pre-burn product gases down to the target temperature. In RCM, the cool down is an order of magnitude shorter, thus the total wall heat-losses would be significantly lower.

The measurement of the entrained gas temperature during the injection demonstrates the importance of the near wall gas temperature in the overall spray mixing process and this highly encourages considering a more realistic temperature field when simulating spray A experiments.

The overall boundary conditions of IFPEN CVP, PPRIME RCM, and PRISME NOSE vessels are summarized in Table 7. The comparison of the measured temperature Tcore and the calculated average temperature Tbulk confirms that similar global temperature heterogeneities are expected in these devices. However, small differences in Tcore at SOI are observed, with PPRIME RCM having the highest value. These differences, caused by both spatial heterogeneity and test-to-test variation, make the estimation of a density representative of the experiments very challenging. The bulk and core densities reported in Table 7 show significantly different values because of the differences between calculated and measured average temperatures. It is worth reminding that considering a constant ambient pressure, the local gas density is expected to have a distribution directly correlated to the local temperature. Thus, it is more relevant to consider this local core density for the definition of a standard target value. This implies that the gas temperature should be characterized extensively in all the locations of the spray region. Particularly, the current work highlights the necessity to measure temperature in the near nozzle and wall regions. The gases in these areas may be entrained to the spray region during the injection and are expected to have a significant effect on the spray behavior.

Table 7

Summary of the boundary conditions in IFPEN, PPRIME and PRISME vessels.

In addition, due to the difference in gas composition in PPRIME, the target density was reached at pressures significantly lower than in IFPEN CVP and PRISME NOSE vessels. It is worth reminding that the spray A experiments have to be performed at target density, while the pressure is adjusted according to the gas composition used in the vessel. The effect of different pressures and densities on the spray behavior will be investigated in part 2 of the paper [20].

5 Summary and Conclusion

In the current paper, the characterization of boundary conditions in different ECN facilities at spray A conditions is reviewed. Particularly, two new vessels based on rapid compression are characterized in details and compared to other facilities employed in ECN community.

The PPRIME RCM and PRISME NOSE have an order of magnitude smaller volume and their body/wall temperatures are 100 K lower compared to IFPEN CVP. In these rapid compression devices, the injector nozzle temperatures at the time of injection are also lower (approximately 10 K lower than CVP). The pre-burn event induces a more pronounced increase of the nozzle temperature compared to the rapid compression event of RCM.

The review of gas temperature measurement demonstrated that thermocouple thermometry is competitive when compared to optical diagnostics, being comparatively easier to use and having a very good accuracy. The use of such sensors for the measurement of ambient gas temperature is also reviewed. It is recommended to use bare-bead thermocouples with wires and junction as small as possible to measure properly the temperature fluctuations. Thicker wires lead to a slower time response with the sensor, causing a larger difference between the gas temperature and the measured junction temperature. The methodology used to correct for this deviation has been presented and it has underlined the strong sensitivity of this correction to the characteristic value of the gas velocity through the junction. The correct estimation or measure of this value is therefore fundamental for the accuracy of the temperature measurement.

The comparison of the temperature distribution in the RCM and CVP showed different levels of heterogeneities with different patterns in each vessel. In all facilities, the standard deviation of the overall temperature distribution at SOI is approximately 5%. The gas temperature has been measured during the injection under inert and reactive conditions in PPRIME using a fast response thermocouple. In inert tests, the results demonstrated that colder gases are entrained to the spray-mixing region. The gases in the boundary (wall) layer are likely to influence the spray evaporation and mixing. In reactive tests, the spray combustion induced significant pressure rise and this increased the ambient gas temperature. All these results emphasize the necessity to consider the specific temperature distribution of each vessel when simulating the spray A.

The gas velocity measurements in IFPEN CVP, PPRIME RCM, and PRISME NOSE are presented. Different levels of velocity and turbulence are observed and this is expected to generate differences in the spray behavior between these vessels.

The second part of this work presented in [20], will investigate spray A formation and combustion in the facilities described here. The boundary conditions characterization presented in this work will be the basis for the analysis of the spray results.

Appendix

thumbnail Fig. A1

Raw Schlieren images showing thermocouple 4 position and the qualitative density gradients in the RCM chamber at different instants: test condition of Figure 22b.

Acknowledgments

The authors thank Herve Arlaud and Gildas Lalizel for providing the thin-wire thermocouples to measure the temperature and thank Alain Claverie for providing assistance during the PIV experiments in PPRIME institute. The authors thank Laurent Hermant for helping with the experiments at IFPEN. The authors greatly acknowledge the valuable discussions with Lyle Pickett, Koji Yasutomi, and Russ Fitzgerald when preparing for “ECN Standardization” topic during ECN 5 workshop. The authors acknowledge the National Research Agency (contract ANR-14-CE22-0015-01) for financial support to the ECN-France project. The authors acknowledge Region Centre Val de Loire (CPER 2007–2013 Energies du Futur) and FEDER for financial support to build the experimental set-up in PRISME institute.

References

  1. Engine Combustion Network (2019) [Online] Available from: https://ecn.sandia.gov/.
  2. Pickett L., Genzale C., Bruneaux G., Malbec L.M., Hermant L., Christiansen C., Schramm J. (2010) Comparison of diesel spray combustion in different high-temperature, high-pressure facilities, SAE Int. J. Engines 3, 2, 156–181.
  3. Meijer M., Somers B., Johnson J., Naber J., Lee S.Y., Malbec L.M., Bruneaux G., Pickett L., Bardi M., Payri R., Bazyn T. (2012) Engine Combustion Network (ECN): Characterization and comparison of boundary conditions for different combustion vessels, Atomization Sprays 22, 9, 777–806.
  4. Bardi M., Payri R., Malbec L.M., Bruneaux G., Pickett L.M., Manin J., Bazyn T., Genzale C. (2012) Engine Combustion Network: Comparison of spray development, vaporization and combustion in different combustion vessels, Atomization Sprays 22, 10, 807–842.
  5. Payri R., Molina S., Salvador F.J., Gimeno J. (2004) A study of the relation between nozzle geometry, internal flow and sprays characteristics in diesel fuel injection systems, KSME Int. J. 18, 1222.
  6. Payri R., Salvador F.J., Gimeno J., de-la-Morena J. (2009) Effects of nozzle geometry on direct injection diesel engine combustion process, Appl. Therm. Eng. 29, 10, 2051–2060.
  7. Kastengren A.L., Tilocco F.Z., Powell C.F., Manin J., Pickett L.M., Payri R., Bazyn T. (2012) Engine Combustion Network (ECN): Measurements of nozzle diameter and hydraulic behavior of diesel sprays, Atomization Sprays 22, 12, 1011–1052.
  8. Bardi M., Bruneaux G., Malbec L.M. (2016) Study of ECN injectors’ behavior repeatability with focus on, in: SAE Technical Paper 2016-01-0845.
  9. Payri R., Garcia-Oliver J.M., Bardi M., Manin J. (2012) Fuel temperature influence on diesel sprays in inert and reacting conditions, Appl. Therm. Eng. 35, 185–195.
  10. Meijer M., Galle J., Somers L.M.T., Griensven J.G.H., Verhelst S. (2013) High-speed characterization of ECN spray A using various diagnostic techniques, SAE Int. J. Engines 6, 2, 1238–1248.
  11. Nesbitt J.E. (2011) Diesel spray mixing limited vaporization with non-ideal and multi-component fuel thermophysical property effects, Michigan Technological University, Houghton, MI.
  12. Zhao L. (2018) An experimental and computational study of fuel spray interaction: Fundamentals and engine applications, Michigan Technological University, Houghton, MI.
  13. Pei Y., Davis M.J., Pickett L.M., Soma S. (2015) Engine Combustion Network (ECN): Global sensitivity analysis of Spray A for different combustion vessels, Combust. Flame 162, 2337–2347.
  14. Nesbitt J.E., Johnson S.E., Pickett L.M., Siebers D.L., Lee S.Y., Naber J.D. (2011) Minor species production from lean premixed combustion and their impact on auto-ignition of diesel surrogates, Energy Fuels 25, 3, 926–936.
  15. Payri R., Vieraa J.P., Wang H., Malbec L. (2016) Velocity field analysis of the high density, high pressure diesel spray, Int. J. Multiph. Flow 80, 69–78.
  16. Malbec L.M., Bruneaux G. (2010) Study of air entrainment of multi-hole diesel injection by particle image velocimetry – effect of neighboring jets interaction and transient behavior after end of injection, SAE Int. J. Engines 3, 1, 107–123.
  17. Sphicas P., Pickett L., Skeen S., Frank J. (2017) ECN experimental database [Online]. Available from: https://ecn.sandia.gov/cvdata/assets/datafiles/velocity/SprayAambient6.xlsx.
  18. Pei Y., Som S., Garcia-Oliver J.M. (2014) Bridging the gap: How spray details affect combustion [Online]. Available from: https://ecn.sandia.gov/wp-content/uploads/2014/12/ECN3_BridgeTheGap.pptx.
  19. Pastor J., Garcia-Oliver J.M., Garcia A., Zhong W., Mico C., Xuan T. (2017) An experimental study on diesel spray injection into a non-quiescent chamber, SAE Int. J. Fuels Lubricants 10, 2, 394–406.
  20. Hespel C., Ben Houidi M., Ajrouche H., Foucher F., Haidous Y., Moreau B., Nilaphai O., Rousselle C., Bellenoue M., Claverie A., Sotton J., Strozzi C., Bardi M., Bruneaux G., Malbec L.-M. (2020) Characterization of the ECN spray A in different facilities. Part 2: Spray vaporization and combustion, Oil Gas Sci. Technol. - Rev. IFP Energies nouvelles. Article submitted.
  21. Strozzi C., Sotton J., Mura A., Bellenoue M. (2008) Experimental and numerical study of the influence of temperature heterogeneities on self-ignition process of methane-air mixtures in a rapid compression machine, Combust. Sci. Technol. 180, 10–11, 1829–1857.
  22. Strozzi C., Mura A., Sotton J., Bellenoue M. (2012) Experimental analysis of propagation regimes during the autoignition of a fully premixed methane–air mixture in the presence of temperature inhomogeneities, Combust. Flame 159, 11, 3323–3341.
  23. Ben-Houidi M., Sotton J., Bellenoue M. (2016) Interpretation of auto-ignition delays from RCM in the presence of temperature heterogeneities: Impact on combustion regimes and negative temperature coefficient behavior, Fuel 186, 476–495.
  24. Ben-Houidi M., Sotton J., Claverie A., Strozzi C., Bellenoue M. (2016) Application of high-speed PIV and Toluene PLIF techniques to study aerodynamics and thermal stratification inside a flat piston rapid compression machine, in: 18th International Symposium on the Application of Laser and Imaging Techniques to Fluid Mechanics, Lisbon, Portugal.
  25. Ben-Houidi M., Sotton J., Bellenoue M., Strozzi C. (2019) Effects of low temperature heat release on the aerodynamics of a flat piston rapid compression machine: Impact on velocity and temperature fields, Proc. Combust. Inst. 37, 4, 4777–4785.
  26. Nilaphai O., Hespel C., Moreau B., Contino F., Bourgeois N., Chanchaona S., Foucher F., Mounaïm-Rousselle C. (2016) New high-pressure and high-temperature chamber for diesel spray characterization, in: 27th Annual Conference on Liquid Atomization and Spray Systems, Brighton, UK.
  27. Nilaphai O. (2018) Vaporization and combustion processes of alcohols and Acetone-Butanol-Ethanol (ABE) blended in n-dodecane for high pressure-high temperature conditions: Application to compression ignition engine, HAL Id: tel-02095626, Orleans, France.
  28. Eagle W.E., Musculus M.P.B., Malbec L.M., Bruneaux G. (2014) Measuring transient entrainment rates of a confined vaporizing diesel jet, in:ILASS Americas 26th Annual Conference on Liquid Atomization and Spray Systems, Portland, OR.
  29. Sciacchitano A., Scarano F., Wieneke B. (2012) Multi-frame pyramid correlation for time-resolved PIV, Exp. Fluids 53, 1087–1105.
  30. Ben-Houidi M., Yasutomi K., Fitzgerald R. (2017) ECN5 proceedings [Online]. Available from: https://ecn.sandia.gov/workshop/ECN5/ECN_5_Standardization_of_Spray_A_V_final.pdf.
  31. Childs P.R.N., Greenwood J.R., Long C.A. (2000) Review of temperature measurement, Rev. Sci. Instrum. 71, 8, 2959–2978.
  32. Griffiths J.F., MacNamara J.P., Sheppard C.G.W., Turton D.A., Whitaker B.J. (2002) The relationship of knock during controlled autoignition to temperature inhomogeneities and fuel reactivity, Fuel 81, 2219–2225.
  33. Strozzi C., Sotton J., Mura A., Bellenoue M. (2009) Characterization of a two dimensional temperature field within a rapid compression machine using a toluene planar laser induced fluorescence imaging technique, Meas. Sci. Technol. 20, 125403.
  34. Tran K.H., Guibert P., Morin C., Bonnety J., Pounkin S., Legros G. (2015) Temperature measurements in a rapid compression machine using anisole planar laser-induced fluorescence, Combust. Flame 162, 10, 3960–3970.
  35. Desgroux P., Gasnot L., Sochet L. (1995) Instantaneous temperature measurement in a rapid-compression machine using laser rayleigh scattering, Appl. Phys. B 61, 69–72.
  36. Clarkson J., Griffiths J.F., Macnamara J.P., Whitaker B. (2001) Temperature fields during the development of combustion in a rapid compression machine, Combust. Flame 125, 3, 1162–1175.
  37. Miles R.B., Lempert W.R., Forkey J.N. (2001) Laser Rayleigh scattering, Meas. Sci. Technol. 12, R33–R51.
  38. Uddi M., Das A.K., Sung C. (2012) Temperature measurements in a rapid compression machine using Mid-IR H2O absorption spectroscopy near 7.6 μm, Appl. Opt. 51, 22, 5464–5476.
  39. Das A.K., Uddi M., Sung C. (2012) Two-line thermometry and H2O measurement for reactive mixtures in rapid compression machine near 7.6 μm, Combust. Flame 159, 12, 3493–3501.
  40. Nasir E.F., Farooq A. (2017) Time-resolved temperature measurements in a rapid compression machine using quantum cascade laser absorption in the intrapulse mode, Proc. Combust. Inst. 36, 4453–4460.
  41. Förster F., Crua C., Davy M., Ewart P. (2019) Temperature measurements under diesel engine conditions using laser induced grating spectroscopy, Combust. Flame 199, 249–257.
  42. Förster F.J., Crua C., Davy M., Ewart P. (2017) Time-resolved gas thermometry by laser-induced grating spectroscopy with a high-repetition rate laser system, Exp. Fluids 58, 87.
  43. Patton R.A., Gabet K.N., Jiang N., Lempert W.R., Sutton J.A. (2012) Multi-kHz temperature imaging in turbulent non-premixed flames using planar Rayleigh scattering, Appl. Phys. B: Lasers Opt. 108, 377–392.
  44. Miller J.D., Slipchenko M.N., Meyer T.R., Stauffer H.U., Gord J.R. (2010) Hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering for high-speed gas-phase thermometry, Opt. Lett. 35, 14, 2430–2432.
  45. Sun B.N., Chang J., Lian J., Wang Z.L., Lv G.P., Liu X.Z., Wang W.J., Zhou S., Wei W., Jiang S., Liu Y.N., Luo S., Liu X.H., Liu Z., Zhang S.S. (2013) Accuracy improvement of Raman distributed temperature sensors based on eliminating Rayleigh noise impact, Opt. Commun. 306, 117–120.
  46. Heitor M.V., Moreira A.L.N. (1993) Thermocouples and sample probes for combustion studies, Prog. Energy Combust. Sci. 19, 259–278.
  47. Paranthoen P., Lecordier J. (1996) Mesures de temperature dans les ecoulements turbulents, Rev. Gen. Therm. 35, 283–308.
  48. Paranthoen P., Lecordier J., Petit C., Gajan P. (1983) Survey and Recent developments of frequency response studies of cold wires and fine wire thermocouples in turbulent heated flow, in: Eighth Austalasian Fluid Mechanics Conference, University of Newcastle, N.S.W.
  49. Tagawa M., Kato K., Ohta Y. (2005) Response compensation of fine wire temperature sensors, Rev. Sci. Instrum. 76, 094904.
  50. Arwatz G., Bahri C., Smits A.J., Hultmark M. (2013) Dynamic calibration and modeling of a cold wire for temperature measurement, Meas. Sci. Technol. 24, 125301.
  51. Dupont A., Paranthoen P., Lecordier J.C., Gajan P. (1984) Influence of temperature on the frequency response of fine-wire thermocouples over the range (300K–800K) in airflows, J. Phys. E: Sci. Instrum. 17, 808–812.
  52. Thierry L. (2006) Microsondes thermoelectriques, Techniques de l’ingenieur R 2 765, 1–24.
  53. Pitts W., Braun E., Peacock R., Mitler H., Johnsson E., Reneke P., Blevins L. (2003) Temperature uncertainties for bare-bead and aspirated thermocouple measurements in fire environments, in: Gritzo L., Alvares N. (eds), Thermal measurements: The foundation of fire standards, ASTM International, West Conshohocken, PA, pp. 3–15.
  54. Kent J.H. (1970) A noncatalytic coating for platinum-rhodium thermocouples, Combust. Flame 28, 301.
  55. Madson J.M., Theby E.A. (1983) SiO2 coated thermocouples, Combust. Sci. Technol. 36, 205–209.
  56. Caldwell F.R. (1962) Thermocouple materials, National Bureau of Standards Monograph No. 40, Washington, DC.
  57. Betchov R. (1948) L’influence de la conduction thermiquesur les anenometres a fils chauds, Proc. Koninklijke Nederlandse Akademie van Wetenschappen 51, 721.
  58. Corrsin S. (1963) Turbulence: experimental methods, Encycl. Phys. 3, 8, 524–590.
  59. Collis D.C., Williams M.J. (1959) Two dimensional convection from heated wires at low Reynolds numbers, J. Fluid Mech. 6, 4, 357–384.
  60. Collis D.C., Williams M.J. (1959) Two-dimensional convection from heated wires at low Reynolds Numbers, J. Fluid Mech. 6, 4, 357–384.
  61. Kramers H. (1946) Heat transfer from spheres to flowing media, Physica 12, 61–80.
  62. Churchill S.W., Bernstein M. (1977) A correlating equation for forced convection from gases and liquids to a circular cylinder in crossflow, J. Heat Trans. 99, 300–306.
  63. Pei Y., Som S., Garcia-Oliver J.M. (2014) ECN3 proceedings [Online]. Available from: https://ecn.sandia.gov/wp-content/uploads/2014/12/ECN3_BridgeTheGap.pptx.

All Tables

Table 1

Comparison of IFPEN, PPRIME, and PRISME vessels when set for spray A experiments.

Table 2

RCM test conditions for PIV experiments.

Table 3

Comparison of turbulence characteristics: Sandia CVP and Caterpillar CPV reported in [17, 30].

Table 4

Summary of optical diagnostic techniques that could be used in ECN spray vessels [31].

Table 5

Summary of temperature measurement with thermocouples under spray A conditions [2, 3].

Table 6

Temperature heterogeneities levels at spray A conditions.

Table 7

Summary of the boundary conditions in IFPEN, PPRIME and PRISME vessels.

All Figures

thumbnail Fig. 1

(a) 3D view and (b) cross section of IFPEN CVP.

In the text
thumbnail Fig. 2

Schematic of PPRIME institute RCM. (a) General view: (1) hydraulic cylinder, (2) cam, (3) guiding wheel, (4) brake system, (5) return pneumatic cylinder, (6) piston/cylinder, (7) combustion chamber. (b) Zoomed section view on the piston cylinder assembly.

In the text
thumbnail Fig. 3

Typical pressure profile of inert and reactive conditions at spray A conditions in PPRIME.

In the text
thumbnail Fig. 4

The NOSE chamber: (a) the overview of real image, (b) the cross-section view in middle plane of chamber, (c) the top-view without chamber head, and (d) the cross-section view A–A without chamber head.

In the text
thumbnail Fig. 5

NOSE operation: comparison of the feedback and command velocities as a function of time referenced to the initial position (280 CAD before TDC).

In the text
thumbnail Fig. 6

Thermodynamic condition of NOSE operation: (a) an example of pressure, temperature and density temporal evolution, (b) pressure and temperature evolutions from −20 ms to 20 ms after TDC for 10 repeated tests.

In the text
thumbnail Fig. 7

Schematic of the PIV setup in IFPEN.

In the text
thumbnail Fig. 8

Schematic of the PIV setup in PPRIME.

In the text
thumbnail Fig. 9

Schematic of PIV setup in PRISME.

In the text
thumbnail Fig. 10

Average and standard deviation of velocity magnitude at time of injection (IFPEN CVP using ensemble average).

In the text
thumbnail Fig. 11

Average and turbulent kinetic energy at time of injection (IFPEN CVP using ensemble average).

In the text
thumbnail Fig. 12

(a) Average and (b) standard deviation of velocity magnitude at time of injection ±10 ms (PPRIME).

In the text
thumbnail Fig. 13

(a) Average and (b) turbulent kinetic energy at time of injection ±10 ms (PPRIME).

In the text
thumbnail Fig. 14

(a) Average velocity magnitude field and (b) turbulent kinetic energy field at SOI (PRISME).

In the text
thumbnail Fig. 15

Comparison of the injector nozzle temperature in IFPEN CVP and PPRIME RCM. The time is referenced to the SOI.

In the text
thumbnail Fig. 16

Schematic of the 7.6 μm thermocouple sensor (PPRIME institute): (a) sensor prongs mounting, (b) zoom on the thermocouple junction welded successfully, (c) zoom on junction having a spherical shape.

In the text
thumbnail Fig. 17

Main parameters affecting the measured temperature correction: (a) effect of dw, (b) effect of dj/dw, (c) effect of the gas velocity U (dw = 7.6 μm), (d) effect of the Nu correlation (dw = 7.6 μm). PPRIME RCM test at spray A target conditions (time axis is referenced to the end of compression).

In the text
thumbnail Fig. 18

(a) Illustration of thermocouples positions, (b) picture of thermocouples mounted in the center plane of the PPRIME RCM, (c) position of the thermocouple in the near nozzle measurement tests.

In the text
thumbnail Fig. 19

Illustration of the measured temperature from a single test compared to adiabatic core temperature.

In the text
thumbnail Fig. 20

Location of the temperature measurement in PRISME. Four thin-wire thermocouples of 13 μm nominal diameter are used: (a) illustration of measurement positions, (b) picture of thermocouples mounted in the centerline of the NOSE chamber, and (c) picture of the thermocouple junction.

In the text
thumbnail Fig. 21

Illustration of typical temperature profile in PRISME (thermocouple 2): (a) example of density, pressure, and temperature profiles around TDC, (b) average temperature profile at 900 K target condition, shaded-error bar represents the standard deviation of 10 measurements.

In the text
thumbnail Fig. 22

Temperature measurement in the near nozzle region: (a) at inert conditions, (b) at reactive conditions.

In the text
thumbnail Fig. 23

Comparison of the averaged normalized temperature measured in different ECN facilities: (a) temperature heterogeneity in the spray axis, (b) temperature heterogeneity in the transverse axis at 40 mm from the injector nozzle. IFPEN and Sandia data from [3].

In the text
thumbnail Fig. A1

Raw Schlieren images showing thermocouple 4 position and the qualitative density gradients in the RCM chamber at different instants: test condition of Figure 22b.

In the text

Les statistiques affichées correspondent au cumul d'une part des vues des résumés de l'article et d'autre part des vues et téléchargements de l'article plein-texte (PDF, Full-HTML, ePub... selon les formats disponibles) sur la platefome Vision4Press.

Les statistiques sont disponibles avec un délai de 48 à 96 heures et sont mises à jour quotidiennement en semaine.

Le chargement des statistiques peut être long.