the evolution and interaction of two shock accelerated
play

The evolution and interaction of two shock-accelerated, unstable gas - PowerPoint PPT Presentation

The evolution and interaction of two shock-accelerated, unstable gas cylinders Chris Tomkins, Kathy Prestridge, Paul Rightley, Mark Marr-Lyon, Robert Benjamin Summer students: James Doyle and Michael Schneider Dynamic Experimentation Division,


  1. The evolution and interaction of two shock-accelerated, unstable gas cylinders Chris Tomkins, Kathy Prestridge, Paul Rightley, Mark Marr-Lyon, Robert Benjamin Summer students: James Doyle and Michael Schneider Dynamic Experimentation Division, Los Alamos National Laboratory Peter Vorobieff Department of Mechanical Engineering, University of New Mexico Cindy Zoldi Applied Physics Division, Los Alamos National Laboratory

  2. PIV PIV Fog Fog generator generator IC Experimental setup: Gas cylinders shock tube SF 6 air air suction N Y D

  3. Overview • Examine interaction of planar shock with 2 gas cylinders, separated spanwise. • S = 1.2D to 2.0D. Shock S (D = cylinder diameter) D • Goal: Investigate the evolution of the interacting, RM-unstable cylinders. Issues of interest include: • What is the effect of the interaction on the resulting flow morphologies? On the initial vorticity deposition? On the post-shock vortex development? • How sensitive is the flow evolution to the initial separation S?

  4. Single shock-accelerated cylinder Double-cylinder “vortex blob” simulation (a) (b) (c) 8 8 8 7 7 7 Spanwise Spanwise Spanwise 2D 6 6 6 5 5 5 2 1 0 -1 -2 2 1 0 -1 -2 2 1 0 -1 -2 Streamwise Streamwise Streamwise

  5. Double-cylinder interaction: weak S � 2.0D Shock S � 1.8D

  6. Double-cylinder interaction: moderate S � 1.6D Shock S � 1.5D

  7. Double-cylinder interaction: strong S � 1.4D Shock S � 1.2D

  8. PIV images: double cylinder Two-frame cross-correlation, flow left to right, 6 th pulse • • S = 2.0D. Note non-uniform seeding. t o t 1 Shock

  9. Double-cylinder velocity field: PIV 12 10 • Double-cylinder data, 6 th pulse, 10 m/s S = 2.0D 8 • Two-frame cross- 2 + v' 2 ) 1/2 (u' correlation z (mm) (m/s) 6 40.0 (Christensen et al., 35.0 2000) 30.0 25.0 20.0 4 • Not smoothed 15.0 10.0 • Contours are 5.0 fluctuating velocity 2 magnitude 0 Shock 12 10 8 6 4 2 0 x (mm)

  10. Double-cylinder vorticity field 12 • Same realization 11 • Vorticity 10 contours � y (1/s) 9 • Not smoothed 60 8 48 36 7 24 z (mm) 12 Shock 6 6 4 -4 5 -6 -12 4 -24 • And the ratio of -36 3 -48 circulations is -60 2 � outer 3 1 � 0 12 10 8 6 4 2 0 � inner x (mm)

  11. Correlation-based ensemble averaging Match one image (template) to each individual realization. Desire optimum match between template 2 and image, i.e. minimize mean sq. error: e I ( x ) I ( x x ) dA o � t � � � x This requires maximizing I ( x ) I ( x x ) dA w.r.t. o o � t � � D Do for each realization, then extract and average (Soloff, 1997) Yields cond. avg.: I ( x x ) | x o o � � � This avg. becomes the new template. I t Properties: x D -Minimizes dependence on initial choice of template. x I o -Converges quickly.

  12. Correlation-based ensemble average S � 1.2D, Ensemble average Shock S � 1.2D, Individual realization

  13. Fluctuating intensity fields, S = 1.2D t 750 s t 470 s � � � � Total Fluctuating Total Fluctuating

  14. RMS of fluctuating intensity RMS Intensity vs. Time for several values of S/D. 40 RMS of Initial Conditions 35 1.5 1.4 30 1.3 RMS intensity, � I Normalized IC RMS Intensity 25 1.2 Mean 1.1 20 1 0.9 15 0.8 S/D = 1.2 S/D = 1.4 10 0.7 S/D = 1.5 S/D = 1.6 0.6 S/D = 1.8 5 S/D = 2.0 0.5 1.2 1.4 1.6 1.8 2 S/D 0 0 100 200 300 400 500 600 700 800 Time after shock passage ( � sec)

  15. Small-scale activity: single cylinder 800 � m A 12 10 A 8 Streamwise (mm) 6 1.0mm B 4 2 B 0 2 4 6 8 10 Spanwise (mm)

  16. Vorticity and swirling strength 12 12 11 11 10 10 9 9 � z 8 (1/s) 8 � ci Streamwise (mm) Streamwise (mm) 50 300 7 7 40 200 30 100 6 6 20 -100 10 -200 -10 5 5 -300 -20 -30 4 4 -40 -50 3 3 2 2 1 1 0 0 0 5 10 0 5 10 Spanwise (mm) Spanwise (mm)

  17. Conclusions • The degree of cylinder-cylinder interaction, and hence the resulting flow morphology, is highly sensitive to the initial cylinder separation. • Different separations may lead to weak, moderate, or strong interactions. • An idealized “vortex blob” simulation leads to very different flow morphologies than experiment, suggesting that the inner vortices are weakened by interaction. • Vorticity fields calculated from high-resolution PIV measurements confirm that the inner vortices are significantly weaker, even for S/D = 2.0: / 3 � � � outer inner

  18. Conclusions • A correlation-based ensemble averaging procedure effectively captures the large and intermediate scales of the flow, providing confirmation of the experimental repeatability, and permitting decomposition of the density field into mean and fluctuating components. • The RMS intensity fluctuations based on this decomposition are substantially greater for the case of “moderate” interaction than for the “strong” or “weak” interaction cases, despite comparable initial RMS values. • High-resolution PIV data resolves mm-scale vortices being convected around the vortex cores.

Recommend


More recommend