770 likes | 1.02k Views
2D Asymmetric Tensor Field Analysis and Visualization. Eugene Zhang, Darrel Palke, Harry Yeh, Zhongzang Lin, Guoning Chen, Robert S. Laramee. Introduction. Flow visualization has a wide range of applications in areo- and hydro-dynamics: Weather prediction Aircraft and missile design
E N D
2D Asymmetric Tensor Field Analysis and Visualization Eugene Zhang, Darrel Palke, Harry Yeh, Zhongzang Lin, Guoning Chen, Robert S. Laramee
Introduction • Flow visualization has a wide range of applications in areo- and hydro-dynamics: • Weather prediction • Aircraft and missile design • Natural disaster modeling: • Tsunami, hurricane, and tornado • Water quality study
Introduction • Existing techniques often focus on velocity [Laramee et al. 2004, Laramee et al. 2007] • Good for visualizing particle movement Trajectories Vector magnitudes Topology
Introduction • Basic types of non-translational motions Rotation (+/-) Expansion = Positive Isotropic Scaling Contraction = Negative Isotropic Scaling Anisotropic stretching
Introduction • Given a vector field , the local linearization at is: translation Velocity gradient: rotation, isotropic scaling, anisotropic stretching
Introduction • Rotation • Isotropic scaling • Anisotropic scaling
Introduction • Flow motions and physical meanings: [Batchelor 1967, Fischer et al. 1979, Ottino 1989, Sherman 1990] • Rotation: • Vorticity • Isotropic scaling: • volume change and/or stretching in the third dimension • Anisotropic stretching: • rate of angular deformation, related to energy dissipation and rate of fluid mixing
Introduction • Velocity gradient tensor has been used in vector field visualization • Singularity classification (source, sink, saddle, etc) [Helman and Hesselink 1991] • Periodic orbit extraction [Chen et al. 2007] • Attachment and separation detection [Kenwright 1998] • Vortex core identification [Sujudi and Haimes 1995, Jeong and Hussain 1995, Peikert and Roth 1999, Sadarjoen and Post 2000]
Introduction • However, tensor field structures were not investigated and applied to flow visualization • Velocity gradient is asymmetric • Past work in tensor field visualization focus on symmetric tensors [Delmarcelle and Hesselink 1994, Hesselink et al. 1997, Tricoche et al. 2001, Tricoche et al. 2003, Hotz et al. 2004, Zheng and Pang 2004, Zheng et al. 2005, Zhang et al. 2007]
Introduction • Symmetric tensors • have two families of mutually perpendicular real-valued eigenvectors • Asymmetric tensors • do not always have real-valued eigenvectors • when they do, the eigenvectors are not always mutually perpendicular
Introduction • Zheng and Pang provide analysis on 2D asymmetric tensor fields Circular points Complex domains (major dual-eigenvectors) Degenerate curves Real domains (major eigenvectors) Image courtesy: Xiaoqiang Zheng and Alex Pang
Introduction • Interesting questions • Orientation of flow rotation • Circular point • Dual-eigenvector computation • What roles do eigenvalues play? • What can these concepts tell us about the flow? • And how to present such information?
Introduction • Our contributions • Present geometric construction of dual-eigenvectors • Enable circular point classification • Incorporate flow orientation • Provide two techniques to better illustrate tensor structures inside complex domains • Pseudo-eigenvectors • Glyph packing • Describe analysis based on eigenvalues • Provide physical interpretation of our analysis
Outline • Review tensors • Eigenvector analysis • Eigenvalue analysis • Applications
Tensor Review • A 2x2 tensor is matrix • Typically asymmetric • Symmetric tensor • Anti-symmetric tensor • A tensor field is a continuous tensor-valued function
Tensor Review • Symmetric tensor • Eigenvalues • Eigenvectors: major and minor • Hyperstreamlines • Degenerate points • Tensor index
Tensor Review • Asymmetric tensors • Real domains and complex domains • Degenerate curves • Dual-eigenvectors (bisector) • Circular points
Our Analysis • Tensor reparameterization • Eigenvalue and eigenvector manifolds
Decomposition • Isotropic scaling: • Rotation: • Anisotropic stretching:
Decomposition • The set of 2x2 tensors can be parameterized by , , , and , such that • , and • This is a four-dimensional space • Can we focus on configuration spaces with lower-dimensions?
Decomposition • Eigenvalues only depend on , , and • Eigenvectors are dependent on , , and • Define eigenvector and eigenvalue manifolds
Eigenvector Manifold • Eigenvectors of is the same as is the same as which can rewritten as
Eigenvector Manifold Image credit: http://math.etsu.edu/MultiCalc/Chap3/Chap3-4/sphere1.gif
Eigenvector Manifold • Eigenvalues are constant along any latitude
Eigenvector Manifold • Along any latitude • Angular component of the eigenvectors and dual-eigenvectors linearly depend on
We can focus on any longitude and understand how eigenvectors/dual-eigenvectors change Eigenvector Manifold
Eigenvector Manifold • We can focus on any longitude and understand how eigenvectors/dual-eigenvectors change
Eigenvector Manifold • Dual-eigenvectors never change along any longitude • Major and minor dual-eigenvectors swap roles when , i.e., symmetric tensors • Major and minor dual-eigenvectors are the major and minor eigenvectors of the following symmetric matrix:
Eigenvector Manifold • Dual-eigenvectors of
Eigenvector Manifold • Degenerate (circular) points of • Number, location, index, orientation Major Eigenvectors of Symmetric Component Major Dual-Eigenvectors
Eigenvector Manifold • Degenerate (circular) points of • It matters which side of the Equator the point is in Major Eigenvectors of Symmetric Component Major Dual-Eigenvectors
Eigenvector Manifold • Poincaré-Hopf theorem (asymmetric tensors): • Given a continuous asymmetric tensor field defined on a closed surface S such that has only isolated degenerate points , then
Eigenvector Manifold • Which side of the Equator matters • Incorporate the Equator into asymmetric tensor topology
Visualization Black curves: White curves: Blue curves: Degenerate points The Equator Degenerate curves Eigenvector Manifold
Complex domains? Show ellipses Eigenvector Manifold
Eigenvector Manifold • How can we visualize the ellipses? • It is not enough to show or even and • We introduce pseudo-eigenvectors
Eigenvector Manifold • Visualization Blue = Dual-Eigenvectors Blue = Pseudo-Eigenvectors
Eigenvector Manifold Visualization Blue = Dual-Eigenvectors Blue = Pseudo-Eigenvectors
Eigenvector Manifold Visualization Blue = Dual-Eigenvectors Blue = Pseudo-Eigenvectors
Eigenvector Manifold Visualization Vector Field Tensor Field
Decomposition Isotropic scaling: Rotation: Anisotropic stretching:
Eigenvalue Manifold • Eigenvalue depend on , , and • We are interested in relatively strengths among the three components
Eigenvalue Manifold Dominant component