Time-Resolved Raman Spectroscopy in a Stratified-Charge Engine

Sep 23, 1980 - Setchell demonstrated the feasibility of using time-averaged Raman scattering in a combusting homogeneous charge engine.1 A stratified ...
0 downloads 0 Views 616KB Size
22 Time-Resolved Raman Spectroscopy in a Stratified-Charge Engine J. RAY SMITH Laser Probes for Combustion Chemistry Downloaded from pubs.acs.org by YORK UNIV on 12/03/18. For personal use only.

Sandia Laboratories, Livermore, CA 94550

The objectives of this research were to develop techniques to measure both the mean and fluctuating nitrogen density and temperature in a combusting stratified charge engine. Such data is necessary for analytical engine model verification. The method chosen to achieve these measurements was spontaneous vibrational Raman scattering by a pulsed frequency doubled YAG laser to get a time resolution of 10 nsec. The nitrogen density was determined from the Stokes signal and the temperature was determined from the ratio of the anti-Stokes to Stokes signal. Setchell demonstrated the feasibility of using time-averaged Raman scattering in a combusting homogeneous charge engine. A stratified charge engine with good optical access was recently developed, and its precombustion fuel-air distributions were determined by time-averaged Raman spectroscopy. The latter engine's precombustion velocity and turbulence fields were measured by laser Doppler velocimetry and its performance and emissions were quantified by conventional methods. The same engine design was used in the present study. 1

2

3

The short duration of the l a s e r pulse at 5321 Angstroms precluded any movement of the spectrometer g r a t i n g t o allow s p e c t r a l d e t a i l s to be r e s o l v e d nor were there s u f f i c i e n t photons to use a multichannel d e t e c t o r . Therefore the spectrometer g r a t i n g was f i x e d and the e n t i r e nitrogen Stokes band was i n t e g r a t e d by a p h o t o m u l t i p l i e r tube (PMT). S i m i l a r l y , the anti-Stokes spectrum was i n t e g r a t e d by a second p h o t o m u l t i p l i e r . The major problem that must be solved i n making Stokes nitrogen d e n s i t y measurements i n a t u r b u l e n t flow was pointed out by S e t c h e l l i t o be the temperature dependence of the Raman s c a t ­ t e r e d Stokes i n t e n s i t y . Because the t r a n s i t i o n p r o b a b i l i t y i s p r o p o r t i o n a l t o (v + 1) where ν i s the i n i t i a l v i b r a t i o n a l s t a t e , the i n t e g r a t e d Raman s c a t t e r e d i n t e n s i t y i s not a unique f u n c t i o n of number d e n s i t y . However, a t h e o r e t i c a l a n a l y s i s of the

0-8412-0570-l/80/47-134-259$05.00/0 © 1980 American Chemical Society

260

LASER

PROBES

FOR COMBUSTION

CHEMISTRY

Stokes v i b r a t i o n a l Raman spectrum of n i t r o g e n has l e d t o a method o f making n i t r o g e n number d e n s i t y measurements t h a t are e s s e n t i a l l y independent of temperature. A v a r i e t y o f spectrometer s l i t c o n v o l u t i o n s and c e n t e r wavelength s e t t i n g s were s t u d i e d t o determine t h e i r i n f l u e n c e on the i n t e g r a t e d Stokes Raman i n t e n s i t y versus temperature. As p o i n t e d out by Leonardo the b a s i c approach i s t o s e l e c t spectrometer s e t t i n g s which balance the increased t r a n s i t i o n p r o b a b i l i t y o f higher v i b r a t i o n a l s t a t e s against the decrease i n p o p u l a t i o n o f the ground s t a t e as the temperature r i s e s . F i g u r e 1 i s an example o f the r e s u l t s f o r a t r a p e z o i d a l s l i t o f 10 by 50 Angstroms with the c e n t e r wavelength v a r i e d from 6070 t o 6073 Angstroms. By s e t t i n g the spectrometer at 6072 Angstroms, i t was p o s s i b l e t o have the Stokes i n t e n s i t y vary by only ±2% while the temperature v a r i e d from 300 t o 1970 K e l v i n s . The s i g n i f i c a n c e o f t h i s r e s u l t i s t h a t i t i s not necessary t o make simultaneous temperature measurements i n a t u r b u l e n t flow f i e l d i n order t o make d e n s i t y measurements. A s i m i l a r a n a l y s i s o f the a n t i - S t o k e s t o Stokes i n t e n s i t y r a t i o expected from n i t r o g e n as a f u n c t i o n o f bandwidth c e n t e r p o s i t i o n o f the a n t i - S t o k e s s p e c t r a i s shown i n F i g u r e 2. I t appears extremely d i f f i c u l t t o use spontaneous v i b r a t i o n a l Raman s c a t t e r i n g f o r determining temperatures o f l e s s than 1000 Κ due to the small a n t i - S t o k e s s i g n a l . Although temperatures i n an engine are well above t h i s l e v e l a f t e r the flame passes through the s c a t t e r i n g volume, the present d e t e c t i o n system does not have s u f f i c i e n t background r e j e c t i o n t o make temperature measurements with good s i g n a l t o noise r a t i o s unless the s t r a t i f i ­ c a t i o n i s small and the equivalence r a t i o l e s s than 0.8. The engine i n t h i s work was designed t o s i m p l i f y the f l u i d mechanics f o r modeling purposes. The i n t a k e and exhaust v a l v e s , f u e l i n j e c t o r , spark plug and l a s e r input/output windows were l o c a t e d i n the c y l i n d e r s i d e w a l l s i n t h e c l e a r a n c e volume above the p i s t o n . The head c o n t a i n e d a 70 mm diameter, c l e a r aperture window. The l a s e r beam was passed through the small windows and s c a t t e r e d l i g h t was c o l l e c t e d a t r i g h t angles t o the beam through the l a r g e window i n the top o f the engine. The i n t a k e valve was shrouded which caused the a i r flow t o s w i r l . Propane at 3.35 MP a (485 p s i ) and 375 Κ was i n j e c t e d r a d i a l l y toward the c e n t e r o f the c y l i n d e r . The o v e r a l l equivalence r a t i o o f the data presented was 0.7 and the engine speed was 900 rpm. The experimental arrangement i s shown i n F i g u r e 3. A s h a f t encoder was used t o synchronize the l a s e r pulse t o the chosen crank angle. The Raman s c a t t e r e d l i g h t was imaged onto the s l i t of a 3/4 meter, s i n g l e spectrometer and detected by the cooled PMT's. Only the f i r s t s i x stages o f the PMT were used i n order to maintain l i n e a r i t y over a wide dynamic range. The PMT output charge was i n t e g r a t e d by a p r e a m p l i f i e r , pulse shaped by a spectroscopy a m p l i f i e r , d i g i t i z e d , and c o l l e c t e d by a m i n i ­ computer. The Raman s i g n a l was normalized by the l a s e r pulse energy, and by accounting f o r the a m p l i f i e r g a i n , the number o f

22.

SMITH

Stratified-Charge

Engine

261

10A

200

400

600

800

1000

1200

1400

1600

1800 2000 2200 2400

2600

TEMPERATURE!Κ) Figure 1.

Relative Stokes vibrational Raman intensity for nitrogen for a trape­ zoidal slit function and various center positions

262

LASER

PROBES

FOR COMBUSTION

CHEMISTRY

Figure 2. Intensity ratio of anti-Stokes to Stokes vibrational Raman scattering for a trapezoidal slit function. Center position of Stokes bandpass at 6072 A.

22.

SMITH

Stratified-Charge Engine

263

264

LASER PROBES FOR COMBUSTION CHEMISTRY

photoelectrons was recorded. One thousand measurements were taken during s u c c e s s i v e engine c y c l e s at each s e l e c t e d crank angle. Checks f o r p o s s i b l e s e n s i t i v i t y changes (due t o window t r a n s m i s ­ s i o n ) were made d u r i n g the course of the data a c q u i s i t i o n . The maximum s e n s i t i v i t y v a r i a t i o n observed over a one hour p e r i o d of lean engine o p e r a t i o n was 2.4 percent change i n the mean value. Sub-microsecond time r e s o l u t i o n was achieved by using a Quanta-Ray DCR-1A l a s e r having a pulse width of 10 nsec. Although t h i s l a s e r was capable of producing i n excess of 250 m i l l i j o u l e s per p u l s e , only 50 m i l l i j o u l e s were used i n the experiment. Above t h i s energy l e v e l the chances of window damage are g r e a t l y i n c r e a s e d . A l l of the data presented were gathered with l e s s than 20 MW/cm^ power d e n s i t i e s on the input/output windows. Gas breakdown was avoided by t i l t i n g the f o c u s i n g lens r e l a t i v e t o the beam a x i s thus i n t r o d u c i n g a l a r g e degree of astigmatism. T h i s gave a s c a t t e r i n g volume of 0.5 mm diameter by 1.25 mm l e n g t h . The length was determined by the 4x m a g n i f i c a t i o n of the c o l l e c t i o n o p t i c s and the spectrometer entrance s l i t height of 5 mm. A check of the l i n e a r i t y of the Raman s i g n a l versus both n i t r o g e n d e n s i t y and l a s e r beam energy well beyond the ranges of the experiment was w i t h i n two percent. The Raman s c a t t e r i n g process i s very weak and the number of photoelectrons produced per pulse w i l l obey Poisson s t a t i s t i c s . I f more than 100 p h o t o e l e c t r o n s are produced i n each event, the u n c e r t a i n t y (one standard d e v i a t i o n ) i n the actual number of p h o t o e l e c t r o n s , N, may be approximated by /N and the f r a c t i o n a l u n c e r t a i n t y i s ση = /Ν. Since n i t r o g e n d e n s i t y f l u c t u a t i o n s , σ-f, are not r e l a t e d t o the p h o t o e l e c t r o n f l u c t u a t i o n s , σ ρ , they w i l l combine randomly t o give a s i g n a l f l u c t u a t i o n , 2 2 1/2 σ = (σ.ρ + Op) ' . T h e r e f o r e the f r a c t i o n a l f l u c t u a t i o n s may be assessed by t a k i n g a s u f f i c i e n t l y l a r g e sample of measurements t o compute o and u s i n g the mean Ν value of the s i g n a l t o compute σ ρ . T h i s technique works well provided the actual f l u c t u a t i o n s are at l e a s t h a l f the s i z e of the Poisson s t a t i s t i c a l f l u c t u a t i o n s . Below t h i s l e v e l the u n c e r t a i n t y i n σρ begins t o dominate the computed of value. In t h i s experiment 200 t o 800 photoelectrons were c o l l e c t e d from each l a s e r pulse which gave s t a t i s t i c a l u n c e r t a i n t i e s from 3.4 t o 7 percent. The mean values of the r e l a t i v e n i t r o g e n d e n s i t y versus crank angle at the c e n t e r of the combustion chamber are shown i n F i g u r e 4. A l s o shown (dashed l i n e ) f o r comparison i s the n i t r o g e n d e n s i t y expected from the p i s t o n motion. The d e n s i t y i s r e l a t i v e t o the n i t r o g e n d e n s i t y i n a i r at atmospheric pressure. I g n i t i o n occurred at 358 crank angle degrees and the flame a r r i v e d at the s c a t t e r i n g volume at about 382 degrees. One observes the compres­ s i o n of the unburned gases ahead o f the flame f r o n t a f t e r i g n i t i o n . The r e l a t i v e f l u c t u a t i o n s i n n i t r o g e n d e n s i t y versus crank angle are i n d i c a t e d by the s o l i d curve shown i n F i g u r e 5. P r i o r t o i g n i t i o n , the f l u c t u a t i o n l e v e l was q u i t e low. The f l u c t u a 5

s

22.

SMITH

Stratified-Charge

265

Engine

X/R = 0 0 = 0.7 (PROPANE) 900 RPM INJECTION 3 1 0 - 3 3 5 IGNITION 358

λ\

300

320

340

360

380

-CYCLIC VARIATION SUBTRACTED

400

420

440

CRANK ANGLE, DEGREES

Figure 5.

One standard deviation of fluctuations in nitrogen number density

266

LASER PROBES FOR COMBUSTION CHEMISTRY

100 CRANK ANGLE = 384° 50

/

r- LATE FLAME ARRIVAL

0 50 - CRANK A N G L E 0 50 - CRANK A N G L E = 380° ^ E A R L Y FLAME ARRIVAL

100 - CRANK A N G L E = 372°

ΓΓΉ (1

BEFORE FLAME ARRIVAL

50

Π

"Κ,

0

100

200

300

ilk.

400

500

600

.

700

800

900

1000

PHOTOELECTRONS, (RELATIVE NITROGEN NUMBER DENSITY) Figure 6.

Histograms of nitrogen number density near time of flame arrival

22.

SMITH

Stratified-Charge

Engine

267

t i o n s reached a maximum near 382 degrees due t o a r r i v a l of the flame. The f l u c t u a t i o n s f e l l s l o w l y i n the post flame gases. The l a r g e peak value i s due t o the random v a r i a t i o n s i n the a r r i v a l time o f the flame f r o n t . Near 382 degrees some measure­ ments were made j u s t ahead o f the flame f r o n t where high d e n s i t i e s p r e v a i l and o t h e r measurements were made j u s t behind the flame f r o n t where low d e n s i t i e s p r e v a i l due t o high temperatures. C l e a r l y these l a r g e f l u c t u a t i o n s are due t o c y c l i c v a r i a t i o n s not t u r b u l e n t f l u c t u a t i o n s . The dashed curve i s an attempt t o remove t h i s c y c l i c v a r i a t i o n e f f e c t by using the most probable d e n s i t y value as the mean value o f a normal d i s t r i b u t i o n . The standard d e v i a t i o n of the d i s t r i b u t i o n i s determined from f i t t i n g the data to the s i d e o f the new mean t h a t has not been d i s t o r t e d by flame a r r i v a l . The r e d u c t i o n o f the apparent f l u c t u a t i o n s near the flame a r r i v a l crank angle i s dramatic. Both curves o f F i g u r e 5 have had the Poisson s t a t i s t i c a l f l u c t u a t i o n s s u b t r a c t e d . The histograms o f F i g u r e 6 represent the number o f measure­ ments versus the number o f photoelectrons at crank angles near flame a r r i v a l i n the s c a t t e r i n g volume. The b i n width i s ten photoelectrons and the t o t a l number o f events i s 1,000 f o r each histogram. The normal shape o f the histogram at 372 degrees i s t y p i c a l o f those from 300 t o 372 and 390 t o 420 degrees. At 380 degrees e a r l y flame a r r i v a l caused the long t a i l i n the d i s t r i b u ­ t i o n below the most probable d e n s i t y . S i m i l a r l y , the d i s t o r t i o n of the d i s t r i b u t i o n t o h i g h e r d e n s i t y values at 384 degrees was due t o l a t e flame a r r i v a l s . T h i s e x p l a n a t i o n o f the e f f e c t s o f c y c l i c v a r i a t i o n on the d i s t r i b u t i o n j u s t i f i e s the attempt t o separate them from the real d e n s i t y f l u c t u a t i o n s . The Stokes s i g n a l - t o - n o i s e r a t i o s were o f the order o f t h i r t y t o one even when the flame was i n the s c a t t e r i n g volume. It i s l i k e l y t h a t the i n c r e a s e i n f l u c t u a t i o n s immediately behind the flame f r o n t was flame induced t u r b u l e n c e . Improved f l u c t u a ­ t i o n measurements are expected by using the temperature d e r i v e d from the a n t i - S t o k e s channel f o r c o n d i t i o n a l sampling o f the d e n s i t y data. However, t h i s method cannot be used u n t i l the 5 microsecond i n t e g r a t i o n time o f present d e t e c t i o n e l e c t r o n i c s i s shortened t o reduce the background l u m i n o s i t y s i g n a l on the anti-Stokes channel. Acknowledgment T h i s w o r k s u p p o r t e d b y DOE a n d M o t o r V e h i c l e M a n u f a c t u r e r s Association.

References 1.

Setchell, R. Ε., 18th Annual Rocky Mountain Spectroscopy Conference, University of Denver, Aug. 2-3, 1976. 2. Johnston, S. C., SAE paper 790433, FEB. 1979.

268

3.

LASER PROBES FOR COMBUSTION CHEMISTRY

Johnston, S. C., Robinson, C. W., Rorke, W. S., Smith, J. R., and Witze, P. O., SAE paper 790092, Feb. 1979. 4. Setchell, R. E., 17th Aerospace Sciences Meeting, New Orleans, LA, Jan. 1979. 5. Leonard, D. Α., Project SQUID, Tech. Rep. AVCO-1-PU, 1972.

Received February 1, 1980.