Understanding the Structure of the Turbulent Mixing Layer in Hydrodynamic Instabilities

When a heavy fluid is placed above a light fluid, tiny vertical perturbations in the interface create a characteristic structure of rising bubbles and falling spikes known as Rayleigh-Taylor instability. Rayleigh-Taylor instabilities have received much attention over the past half-century because of their importance in understanding many natural and man-made phenomena, ranging from the rate of formation of heavy elements in supernovae to the design of capsules for Inertial Confinement Fusion. We present a new approach to analyze Rayleigh-Taylor instabilities in which we extract a hierarchical segmentation of the mixing envelope surface to identify bubbles and analyze analogous segmentations of fields on the original interface plane. We compute meaningful statistical information that reveals the evolution of topological features and corroborates the observations made by scientists. We also use geometric tracking to follow the evolution of single bubbles and highlight merge/split events leading to the formation of the large and complex structures characteristic of the later stages. In particular we (i) Provide a formal definition of a bubble; (ii) Segment the envelope surface to identify bubbles; (iii) Provide a multi-scale analysis technique to produce statistical measures of bubble growth; (iv) Correlate bubble measurements with analysis of fields on the interface plane; (v) Track the evolution of individual bubbles over time. Our approach is based on the rigorous mathematical foundations of Morse theory and can be applied to a more general class of applications

[1]  Rephael Wenger,et al.  Volume Tracking Using Higher Dimensional Isocontouring , 2003, IEEE Visualization.

[2]  Han-Wei Shen,et al.  Efficient isosurface tracking using precomputed correspondence table , 2004, VISSYM'04.

[3]  Andrew W. Cook,et al.  The mixing transition in Rayleigh–Taylor instability , 2004, Journal of Fluid Mechanics.

[4]  Valerio Pascucci,et al.  Time-varying reeb graphs for continuous space-time data , 2004, SCG '04.

[5]  Valerio Pascucci,et al.  Morse-smale complexes for piecewise linear 3-manifolds , 2003, SCG '03.

[6]  Andrzej Szymczak Subdomain aware contour trees and contour evolution in time-dependent scalar fields , 2005, International Conference on Shape Modeling and Applications 2005 (SMI' 05).

[7]  H. Edelsbrunner,et al.  Foundations of Computational Mathematics: Minneapolis, 2002: Jacobi Sets , 2004 .

[8]  Jack Snoeyink,et al.  Simplifying flexible isosurfaces using local geometric measures , 2004, IEEE Visualization 2004.

[9]  Xavier Tricoche,et al.  Topological Methods for Flow Visualization , 2005, The Visualization Handbook.

[10]  Hans-Peter Seidel,et al.  Saddle connectors - an approach to visualizing the topological skeleton of complex 3D vector fields , 2003, IEEE Visualization, 2003. VIS 2003..

[11]  Hans-Peter Seidel,et al.  Topological methods for 2D time-dependent vector fields based on stream lines and path lines , 2005, IEEE Transactions on Visualization and Computer Graphics.

[12]  Timothy T. Clark,et al.  A numerical study of the statistics of a two-dimensional Rayleigh–Taylor mixing layer , 2003 .

[13]  Valerio Pascucci,et al.  A Practical Approach to Two-Dimensional Scalar Topology , 2007, Topology-based Methods in Visualization.

[14]  Hans Hagen,et al.  Topology tracking for the visualization of time-dependent two-dimensional flows , 2002, Comput. Graph..

[15]  Paul E. Dimotakis,et al.  Transition stages of Rayleigh–Taylor instability between miscible fluids , 2002, Journal of Fluid Mechanics.

[16]  Bernd Hamann,et al.  A topological hierarchy for functions on triangulated surfaces , 2004, IEEE Transactions on Visualization and Computer Graphics.

[17]  Hans Hagen,et al.  A topology simplification method for 2D vector fields , 2000 .

[18]  Herbert Edelsbrunner,et al.  Hierarchical Morse—Smale Complexes for Piecewise Linear 2-Manifolds , 2003, Discret. Comput. Geom..

[19]  Jesse Freeman,et al.  in Morse theory, , 1999 .

[20]  Xin Wang,et al.  Tracking scalar features in unstructured data sets , 1998 .

[21]  Hans-Peter Seidel,et al.  Extracting higher order critical points and topological simplification of 3D vector fields , 2005, VIS 05. IEEE Visualization, 2005..

[22]  Han-Wei Shen,et al.  Volume tracking using higher dimensional isosurfacing , 2003, IEEE Visualization, 2003. VIS 2003..

[23]  Xin Wang,et al.  Tracking and Visualizing Turbulent 3D Features , 1997, IEEE Trans. Vis. Comput. Graph..

[24]  Bernd Hamann,et al.  Topology-based simplification for feature extraction from 3D scalar fields , 2005, VIS 05. IEEE Visualization, 2005..

[25]  Dan Oron,et al.  Three-dimensional multimode Rayleigh–Taylor and Richtmyer–Meshkov instabilities at all density ratios , 2003 .

[26]  Hans J. W. Spoelder,et al.  Visualization of time-dependent data with feature tracking and event detection , 2001, The Visual Computer.

[27]  Martin Isenburg,et al.  Large mesh simplification using processing sequences , 2003, IEEE Visualization, 2003. VIS 2003..

[28]  Deborah Silver,et al.  Visualizing features and tracking their evolution , 1994, Computer.

[29]  Xavier Tricoche,et al.  Tracking of vector field singularities in unstructured 3D time-dependent datasets , 2004, IEEE Visualization 2004.

[30]  Chandrajit L. Bajaj,et al.  Time-varying contour topology , 2006, IEEE Transactions on Visualization and Computer Graphics.

[31]  Hans Hagen,et al.  Exploring scalar fields using critical isovalues , 2002, IEEE Visualization, 2002. VIS 2002..

[32]  Lambertus Hesselink,et al.  Visualizing vector field topology in fluid flows , 1991, IEEE Computer Graphics and Applications.

[33]  Mark Segal,et al.  Fast shadows and lighting effects using texture mapping , 1992, SIGGRAPH.

[34]  Robert van Liere,et al.  Collapsing flow topology using area metrics , 1999, Proceedings Visualization '99 (Cat. No.99CB37067).