## Efficient Morse Decompositions of Vector Fields

### Cached

### Download Links

- [eecs.oregonstate.edu]
- [web.engr.orst.edu]
- [web.engr.oregonstate.edu]
- DBLP

### Other Repositories/Bibliography

Citations: | 18 - 7 self |

### BibTeX

@MISC{Chen_efficientmorse,

author = {Guoning Chen and Konstantin Mischaikow and Robert S. Laramee and Eugene Zhang},

title = {Efficient Morse Decompositions of Vector Fields},

year = {}

}

### OpenURL

### Abstract

Abstract — Existing topology-based vector field analysis techniques rely on the ability to extract the individual trajectories such as fixed points, periodic orbits and separatrices which are sensitive to noise and errors introduced by simulation and interpolation. This can make such vector field analysis unsuitable for rigorous interpretations. We advocate the use of Morse decompositions, which are robust with respect to perturbations, to encode the topological structures of a vector field in the form of a directed graph, called a Morse connection graph (MCG). While an MCG exists for every vector field, it need not be unique. Previous techniques for computing MCG’s, while fast, are overly conservative and usually results in MCG’s that are too coarse to be useful for the applications. To address this issue, we present a new technique for performing Morse decomposition based on the concept of τ-maps, which typically provides finer MCG’s than existing techniques. Furthermore, the choice of τ provides a natural tradeoff between the fineness of the MCG’s and the computational costs. We provide efficient implementations of Morse decomposition based on τ-maps, which include the use of forward and backward mapping techniques and an adaptive approach in constructing better approximations of the images of the triangles in the meshes used for simulation.Furthermore, we propose the use of spatial τ-maps in addition to the original temporal τ-maps. These techniques provide additional tradeoffs between the quality of the MCG’s and the speed of computation. We demonstrate the utility of our technique with various examples in plane and on surfaces including engine simulation datasets. Index Terms — Vector field topology, Morse decomposition, τmaps, Morse connection graph, flow combinatorialization.

### Citations

8530 | Introduction to Algorithms - Cormen, Leiserson, et al. - 1990 |

1869 |
Numerical Recipes in C: The Art of Scientific Computing
- Press, Teukolsky, et al.
- 1992
(Show Context)
Citation Context ...ulting topology using an adaptive fourth-order Runge-Kutta integration, while the image to the right illustrates the topology of the same vector field using a second-order Runge-Kutta integration [2] =-=[23]-=-. This clearly demonstrates that the ECG’s rely on the employed numerical scheme. (The ECG’s in all the example flows are computed using the algorithms proposed by Chen et al. [4].) These observations... |

255 |
Numerical recipes: The art of scientific computing”, Cambridge University Press, NY. 84 Sriram Venkatesan Born September 28
- Press, Flannery, et al.
- 1989
(Show Context)
Citation Context ...lting topology using an adaptive fourth-order Runge-Kutta integration, while the image to the right illustrates the topology of the same vector field using a second-order Runge-Kutta integration [2], =-=[23]-=-. This clearly demonstrates that the ECGs rely on the employed numerical scheme. (The ECGs in all the example flows are computed using the algorithms proposed by Chen et al. [4].) These observations m... |

214 |
Computer methods for ordinary differential equations and differential-algebraic equations, chapter 5
- Ascher, Petzold
- 1998
(Show Context)
Citation Context ... resulting topology using an adaptive fourth-order Runge-Kutta integration, while the image to the right illustrates the topology of the same vector field using a second-order Runge-Kutta integration =-=[2]-=- [23]. This clearly demonstrates that the ECG’s rely on the employed numerical scheme. (The ECG’s in all the example flows are computed using the algorithms proposed by Chen et al. [4].) These observa... |

192 |
Visualization of vector field topology in fluid flows
- Helman, Hesselink
- 1991
(Show Context)
Citation Context ... followed by a summary and discussion of future work in Section VII. A2 R2 S1 A1 R1 R3 II. RELATED WORK Helman and Hesselink introduced to the visualization community the notion of flow topology [10] =-=[11]-=-. Since then, much research has been done on topological analysis of vector fields in the past two decades. To review all of it is beyond the scope of this paper. Here, we briefly review the most clos... |

188 | Isolated invariant sets and the Morse index - Conley |

152 | Texture synthesis on surfaces
- Turk
(Show Context)
Citation Context ...ursue the individual trajectory-based vector field topology without being concerned with the fact that the obtained ECGs may not be topologically rigor, such as, the applications in texture synthesis =-=[31]-=-, [33] and fluid simulation [27]. For such applications, an ECG can still be extracted from Morse decomposition as an additional step [4]. D. Flow Combinatorialization Based on τ-maps We now turn to t... |

139 | Texture synthesis over arbitrary manifold surfaces
- Wei, Levoy
- 2001
(Show Context)
Citation Context ...the individual trajectory-based vector field topology without being concerned with the fact that the obtained ECGs may not be topologically rigor, such as, the applications in texture synthesis [31], =-=[33]-=- and fluid simulation [27]. For such applications, an ECG can still be extracted from Morse decomposition as an additional step [4]. D. Flow Combinatorialization Based on τ-maps We now turn to the iss... |

134 |
L.: Representation and display of vector field topology in fluid flow data sets
- HELMAN, HESSELINK
- 1989
(Show Context)
Citation Context ... [16] in engine design, depending on their locations. Past work defines the topology of two-dimensional vector fields as fixed points and periodic orbits as well as the separatrices that connect them =-=[10]-=-, [4]. However, analysis and visualization of vector field topology based on individual trajectories can raise G. Chen is with the School of Electrical Engineering and Computer Science, Oregon State U... |

109 | Image Based Flow Visualization - Wijk - 2002 |

107 |
Dynamics and Bifurcations
- Hale, Koçac
- 1991
(Show Context)
Citation Context ...s examples in which proper interpretation can be difficult when performing analysis based on individual trajectories. Figure 1(a) shows an analytical vector field which contains pitchfork bifurcation =-=[9]-=-. The results shown in the two columns of (a) are obtained by computing sample vector values using two different meshes: (left) a regular triangulated mesh with 6144 triangles, and (right) a triangula... |

71 | Flows on surfaces of arbitrary topology
- STAM
- 2003
(Show Context)
Citation Context ...based vector field topology without being concerned with the fact that the obtained ECGs may not be topologically rigor, such as, the applications in texture synthesis [31], [33] and fluid simulation =-=[27]-=-. For such applications, an ECG can still be extracted from Morse decomposition as an additional step [4]. D. Flow Combinatorialization Based on τ-maps We now turn to the issue of flow combinatorializ... |

69 | DESBRUN M.: Discrete multiscale vector field decomposition - TONG, LOMBEYDA, et al. |

62 | Visualizing nonlinear vector field topology
- Scheuermann, Krüger, et al.
(Show Context)
Citation Context ...et al. [30] and Polthier and Preuß [21] give efficient methods to locate fixed points in a vector field. Scheuermann et al. use Clifford algebra to study the non-linear fixed points of a vector field =-=[25]-=- and present the approaches of visualizing higher-order fixed points and non-linear topology of a given vector field [24], [25]. Wischgoll and Scheuermann [34] present an algorithm for detecting perio... |

55 |
Continuous Topology Simplification of Planar Vector Fields
- Tricoche, Scheuermann, et al.
- 2001
(Show Context)
Citation Context ...rested readers to a number of surveys [15], [22]. A. Vector Field Topology Much work has been done to address the extraction of vector field topology in two-dimensional vector fields. Tricoche et al. =-=[30]-=- and Polthier and Preuß [21] give efficient methods to locate fixed points in a vector field. Scheuermann et al. use Clifford algebra to study the non-linear fixed points of a vector field [25] and pr... |

53 |
G.: Detection and visualization of closed streamlines in planar flows
- WISCHGOLL, SCHEUERMANN
(Show Context)
Citation Context ...on-linear fixed points of a vector field [25] and present the approaches of visualizing higher-order fixed points and non-linear topology of a given vector field [24], [25]. Wischgoll and Scheuermann =-=[34]-=- present an algorithm for detecting periodic orbits in planar flows. They also extend this work to 3D vector fields [35] and time-dependent flows [36]. Theisel et al. [28] present a mesh independent a... |

52 | Identifying vector fields singularities using a discrete hodge decomposition
- Polthier, Preuss
- 2003
(Show Context)
Citation Context ...f surveys [15], [22]. A. Vector Field Topology Much work has been done to address the extraction of vector field topology in two-dimensional vector fields. Tricoche et al. [30] and Polthier and Preuß =-=[21]-=- give efficient methods to locate fixed points in a vector field. Scheuermann et al. use Clifford algebra to study the non-linear fixed points of a vector field [25] and present the approaches of visu... |

51 |
Top scientific visualization research problems
- Johnson
- 2004
(Show Context)
Citation Context ... extracting vector field topology than the existing techniques. We point out that addressing such uncertainty in visualization was identified as one of the most important future challenges by Johnson =-=[12]-=-. In order to address this important challenge we present a rather different approach to the representation, extraction and visualization of flow topology. The representation of the global dynamics is... |

50 |
DOLEISCH H.: The state of the art in flow visualization: feature extraction and tracking
- POST, VROLIJK, et al.
(Show Context)
Citation Context ... past two decades. To review all of it is beyond the scope of this paper. Here, we briefly review the most closely related work to this paper. We refer interested readers to a number of surveys [15], =-=[22]-=-. A. Vector Field Topology Much work has been done to address the extraction of vector field topology in two-dimensional vector fields. Tricoche et al. [30] and Polthier and Preuß [21] give efficient ... |

48 | The Topology of Symmetric, Second-Order Tensor Fields - Delmarcelle, Hesselink - 1994 |

45 | Vector Field Design on Surfaces
- Zhang, Mischaikow, et al.
(Show Context)
Citation Context ...s interpretation of the given data. B. Morse Decomposition and Conley Theory Conley index theory and Morse decomposition have been introduced to the scientific visualization community by Zhang et al. =-=[38]-=- and Chen et al. [4], respectively. Morse decomposition focuses on extraction and analysis of the invariant sets of a S1 A3 R4 A2 R1sIEEE TVCG, VOL. ?,NO. ?, AUGUST 200? 4 flow. Fixed points and perio... |

33 | H.: Visual analysis and exploration of fluid flow in a cooling jacket
- LARAMEE, GARTH, et al.
- 2005
(Show Context)
Citation Context ...r prediction, tsunami and hurricane modeling, and airplane design and control. In particular, the existence of fixed points, recirculation zones or periodic orbits can be desirable [4] or undesirable =-=[16]-=- in engine design, depending on their locations. Past work defines the topology of two-dimensional vector fields as fixed points and periodic orbits as well as the separatrices that connect them [10],... |

32 | H.: Investigating swirl and tumble flow with a comparison of visualization techniques
- LARAMEE, WEISKOPF, et al.
- 2004
(Show Context)
Citation Context ...e extrapolated boundary velocity fields that are obtained through simulation of in-cylinder flow. Engineers are interested in knowing whether or not the flows on the surface follow the ideal patterns =-=[17]-=-. The results of applying our framework are shown in Figures 14 and 15, respectively. Figure 14 shows the results of the gas engine simulation data. The first column shows the results using the geomet... |

27 | Topological techniques for efficient rigorous computation in dynamics, Acta Numerica 11 Cambridge Univ - Mischaikow - 2002 |

27 |
Conley index
- Mischaikow, Mrozek
- 2002
(Show Context)
Citation Context ...d analysis of the invariant sets of a S1 A3 R4 A2 R1sIEEE TVCG, VOL. ?,NO. ?, AUGUST 200? 4 flow. Fixed points and periodic orbits are examples of invariant sets. An index called the Conley index [4] =-=[19]-=-, identifies the type of the invariant set. Efficient algorithms also exist for identifying the neighborhoods of the invariant sets [14]. Chen et al. [4] implement a Morse decomposition algorithm usin... |

27 |
Visualization of higher order singularities in vector fields
- Scheuermann, Hagen, et al.
- 1997
(Show Context)
Citation Context ...l. use Clifford algebra to study the non-linear fixed points of a vector field [25] and present the approaches of visualizing higher-order fixed points and non-linear topology of a given vector field =-=[24]-=-, [25]. Wischgoll and Scheuermann [34] present an algorithm for detecting periodic orbits in planar flows. They also extend this work to 3D vector fields [35] and time-dependent flows [36]. Theisel et... |

25 | E.: Vector field editing and periodic orbit extraction using Morse decomposition
- CHEN, MISCHAIKOW, et al.
(Show Context)
Citation Context ... (CFD) [15], weather prediction, tsunami and hurricane modeling, and airplane design and control. In particular, the existence of fixed points, recirculation zones or periodic orbits can be desirable =-=[4]-=- or undesirable [16] in engine design, depending on their locations. Past work defines the topology of two-dimensional vector fields as fixed points and periodic orbits as well as the separatrices tha... |

22 | Exploring Global Dynamics: A Numerical Algorithm Based on the Conley Index Theory - Eidenschink - 1996 |

20 | T.: Grid-independent detection of closed stream lines in 2D vector fields
- THEISEL, WEINKAUF
(Show Context)
Citation Context ...]. Wischgoll and Scheuermann [34] present an algorithm for detecting periodic orbits in planar flows. They also extend this work to 3D vector fields [35] and time-dependent flows [36]. Theisel et al. =-=[28]-=- present a mesh independent approach to compute periodic orbits. In general, previous topology-based techniques are based on individual trajectories. As we have seen in Figure 1, this may lead to ambi... |

14 | Locating Closed Streamlines in 3D Vector Fields
- Wischgoll, Scheuermann
- 2002
(Show Context)
Citation Context ...n-linear topology of a given vector field [24], [25]. Wischgoll and Scheuermann [34] present an algorithm for detecting periodic orbits in planar flows. They also extend this work to 3D vector fields =-=[35]-=- and time-dependent flows [36]. Theisel et al. [28] present a mesh independent approach to compute periodic orbits. In general, previous topology-based techniques are based on individual trajectories.... |

13 | 1986] Theory of Limit Cycles - Ye |

12 | An Algorithmic Approach to Chain Recurrence
- Kalies, Mischaikow, et al.
- 2005
(Show Context)
Citation Context ... and ECG, an MCG is stable with respect to perturbations, i.e. given sufficient information on errors of the vector field it is possible to make rigorous interpretations about the underlying dynamics =-=[14]-=-. In other words, a well defined error, ε > 0, can be bounded and included into the map of the flow domain. We demonstrate the stability of MCGs in Figure 1 (the last two rows). To perform Morse decom... |

12 |
DOLEISCH H.: The State of the Art
- POST, VROLIJK, et al.
(Show Context)
Citation Context ... past two decades. To review all of it is beyond the scope of this paper. Here, we briefly review the most closely related work to this paper. We refer interested readers to a number of surveys [15], =-=[22]-=-. 2.1 Vector Field Topology Much work has been done to address the extraction of vector field topology in 2D vector fields. Tricoche et al. [30] and Polthier and Preuß [21] give efficient methods to l... |

9 | C1-Interpolation for Vector Field Topology Visualization - Sheuermann, Tricoche, et al. - 1999 |

8 | A computational approach to conley’s decomposition theorem
- Ban, Kalies
- 2006
(Show Context)
Citation Context ...n F such as k = m and l = n. Finally, we extract the MCG from F by removing intermediate nodes from F as illustrated in Figure 3 (the bottom graph). The algorithm for MCG construction can be found in =-=[13]-=-. To visualize the MCG, we classify the nodes of the MCG into three types: Source Morse sets, Ri, are nodes absent of incoming edges in the MCG; Sink Morse sets, Ai, are nodes without outgoing edges i... |

8 | Topology Based Flow Visualization: The State of the Art - Laramee, Hauser, et al. - 2007 |

7 | Tracking closed streamlines in time-dependent planar flows
- Wischgoll, Scheuermann, et al.
- 2001
(Show Context)
Citation Context ...ector field [24], [25]. Wischgoll and Scheuermann [34] present an algorithm for detecting periodic orbits in planar flows. They also extend this work to 3D vector fields [35] and time-dependent flows =-=[36]-=-. Theisel et al. [28] present a mesh independent approach to compute periodic orbits. In general, previous topology-based techniques are based on individual trajectories. As we have seen in Figure 1, ... |

7 |
Set arithmetic and the enclosing problem
- Mrozek, Zgliczyński
(Show Context)
Citation Context ...ying vector field is known. Given a triangle T, one covers it with squares of size ". For each square S, define ðSÞ minf ðxÞjx 2 Sg and ðSÞ maxf ðxÞjx 2 Sg. Using rigorous enclosure techniques [1], =-=[20]-=-, one obtains an outer enclosure I S of the true image of the square S integrated forward for all times ðSÞ t ðSÞ. Then, I T [IS, where the union is taken over all squares S is an outer approximation... |

6 | Polygonal Approximation of Flows - Boczko, Kalies, et al. - 2007 |

5 |
Grid-Independent Detection of Closed Stream
- Theisel, Weinkauf
- 2004
(Show Context)
Citation Context ...]. Wischgoll and Scheuermann [34] present an algorithm for detecting periodic orbits in planar flows. They also extend this work to 3D vector fields [35] and time-dependent flows [36]. Theisel et al. =-=[28]-=- present a mesh independent approach to compute periodic orbits. In general, previous topology-based techniques are based on individual trajectories. As we have seen in Fig. 1, this may lead to ambigu... |

4 |
Set arithmetic and the enclosing problem in dynamics
- Mrozek, Zgliczyński
- 2000
(Show Context)
Citation Context ...r field is known. Given a triangle T one covers it with squares of size ε. For each square S define τ∗(S) = min{τ(x) | x ∈ S} and τ ∗ (S) = max{τ(x) | x ∈ S}. Using rigorous enclosure techniques [1], =-=[20]-=- obtains an outer enclosure IS of the true image of the square S integrated forward for all times τ∗(S) ≤ t ≤ τ ∗ (S). Then IT = ∪IS, where the union is taken over all squares S, is an outer approxima... |

2 | Polygonal approximation of flows. 2006. (forthcoming), The online version available at: http://www.math.fau.edu/kalies/papers/paf.pdf - Boczko, Kalies, et al. |

2 |
Flows on Surfaces of Arbitrary Topology,” Proc
- Stam
- 2003
(Show Context)
Citation Context ...ased vector field topology without being concerned with the fact that the obtained ECGs may not be topologically rigorous such as the applications in texture synthesis [31], [33] and fluid simulation =-=[27]-=-. For such applications, an ECG can still be extracted from Morse decomposition as an additional step [4].CHEN ET AL.: EFFICIENT MORSE DECOMPOSITIONS OF VECTOR FIELDS 853 3.4 Flow Combinatorializatio... |

2 | Image Based Flow Visualization,” Proc - Wijk - 2002 |

1 |
Topological analysis and visualization of asymmetric tensors (under review
- Zhang, Yeh, et al.
- 2007
(Show Context)
Citation Context ...s). To see this, we show Fig. 11. This example shows the comparison of the visualizations of the sample rate in the field (right) with the tensor analysis visualization results by Zhang et al. (left) =-=[39]-=-. The color coding of the density map uses red for the region with a larger sampling rate and blue for a lower sampling rate. In the tensor analysis results, they use green for positive isotropic scal... |