The Effect of Data Transformations on Scalar Field Topological Analysis of High-Order FEM Solutions
Pith reviewed 2026-05-24 20:23 UTC · model grok-4.3
The pith
The choice of transformation method for high-order FEM data changes the topological features extracted from scalar fields.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that transforming high-order FEM simulation data using different methodologies produces diverse behaviors in topological analysis. The empirical comparison of two commonly used transformation methodologies together with the L-SIAC filter shows that each produces distinct topological features from the same high-order scalar fields, and therefore pipelines built with current topological tools require explicit consideration of the transformation choice.
What carries the argument
Data transformation methodologies that convert discontinuous high-order polynomial data into continuous low-order representations for input to topological analysis tools.
If this is right
- Topological features extracted from HO-FEM scalar fields vary with the transformation method applied before analysis.
- The L-SIAC filter produces results that differ from the two standard transformation approaches in at least some test cases.
- Analysis pipelines must evaluate the effect of the chosen transformation on the specific features of interest.
- Inconsistent topological structures can appear from the same simulation data when different preparation steps are used.
Where Pith is reading between the lines
- Directly extending topological tools to accept high-order discontinuous inputs would remove the transformation step and its observed variability.
- Quantifying how each transformation preserves or distorts specific topological invariants could guide method selection for particular analysis goals.
- Applying the same comparison to vector or tensor fields from the same simulations would test whether the sensitivity is limited to scalar cases.
Load-bearing premise
The two commonly used transformation methodologies together with the L-SIAC filter are representative of possible transformations and the diverse behaviors will generalize to other high-order FEM datasets and topological tools.
What would settle it
Running the same topological extraction on a high-order FEM dataset whose analytic topology is known in advance and finding identical features after every transformation method would falsify the claim of diverse behaviors.
Figures
read the original abstract
High-order finite element methods (HO-FEM) are gaining popularity in the simulation community due to their success in solving complex flow dynamics. There is an increasing need to analyze the data produced as output by these simulations. Simultaneously, topological analysis tools are emerging as powerful methods for investigating simulation data. However, most of the current approaches to topological analysis have had limited application to HO-FEM simulation data for two reasons. First, the current topological tools are designed for linear data (polynomial degree one), but the polynomial degree of the data output by these simulations is typically higher (routinely up to polynomial degree six). Second, the simulation data and derived quantities of the simulation data have discontinuities at element boundaries, and these discontinuities do not match the input requirements for the topological tools. One solution to both issues is to transform the high-order data to achieve low-order, continuous inputs for topological analysis. Nevertheless, there has been little work evaluating the possible transformation choices and their downstream effect on the topological analysis. We perform an empirical study to evaluate two commonly used data transformation methodologies along with the recently introduced L-SIAC filter for processing high-order simulation data. Our results show diverse behaviors are possible. We offer some guidance about how best to consider a pipeline of topological analysis of HO-FEM simulations with the currently available implementations of topological analysis.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents an empirical study evaluating the effects of two commonly used data transformation methodologies and the L-SIAC filter when preparing high-order FEM simulation outputs (typically polynomial degree up to 6, with element-boundary discontinuities) for input to topological analysis tools designed for linear continuous data. The central claim is that these transformations produce visibly diverse behaviors in the resulting topological outputs, and the authors provide guidance on pipeline construction using currently available tools.
Significance. If the reported differences are reproducible and not artifacts of the specific test cases, the work is significant for the growing community using topological analysis on HO-FEM data. It supplies concrete evidence that preprocessing choices are not neutral and must be documented, which is a practical contribution even if the study does not claim the three methods exhaust the design space.
major comments (2)
- [Results] Results section: the claim of 'diverse behaviors' rests on visual comparison of topological features; without quantitative metrics (e.g., persistence diagram distances, Betti number differences, or statistical tests across multiple datasets), it is difficult to judge whether the observed differences are systematic or case-specific.
- [Methods] Methods: the description of the two 'commonly used' transformation methodologies lacks sufficient implementation detail (exact interpolation orders, handling of quadrature points, and boundary treatment) to permit independent reproduction or extension of the comparison.
minor comments (2)
- [Introduction] Abstract and introduction: the statement that 'there has been little work evaluating the possible transformation choices' would benefit from one or two additional citations to related preprocessing studies in visualization or CFD.
- [Figures] Figure captions: several figures showing topological outputs would be clearer if they explicitly labeled the transformation method, polynomial degree, and mesh resolution used for each panel.
Simulated Author's Rebuttal
We appreciate the referee's feedback on our manuscript. The comments highlight areas where we can enhance the presentation of our empirical study. We respond to each major comment below.
read point-by-point responses
-
Referee: [Results] Results section: the claim of 'diverse behaviors' rests on visual comparison of topological features; without quantitative metrics (e.g., persistence diagram distances, Betti number differences, or statistical tests across multiple datasets), it is difficult to judge whether the observed differences are systematic or case-specific.
Authors: We agree that quantitative metrics would strengthen the evaluation. In the revised manuscript we will add Betti number differences computed on the presented test cases to provide a quantitative complement to the visual comparisons. revision: yes
-
Referee: [Methods] Methods: the description of the two 'commonly used' transformation methodologies lacks sufficient implementation detail (exact interpolation orders, handling of quadrature points, and boundary treatment) to permit independent reproduction or extension of the comparison.
Authors: We thank the referee for noting this gap. The revised version will expand the Methods section with the requested implementation details on interpolation orders, quadrature point handling, and boundary treatment for both transformation methodologies. revision: yes
Circularity Check
No significant circularity: purely empirical comparison
full rationale
The paper conducts an empirical study comparing the effects of two common data transformation methodologies and the L-SIAC filter on topological analysis outputs for high-order FEM scalar fields. No mathematical derivations, parameter fitting, predictions from first principles, or uniqueness theorems are present. The central claim rests on direct observation of differing topological features (e.g., critical points, persistence diagrams) across the tested transformations on provided datasets. Self-citations, if any, are incidental and not load-bearing for the reported behaviors. The work is self-contained as an observational benchmark against external tools and data.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
A. Acharya and V . Natarajan. A parallel and memory efficient algorithm for constructing the contour tree. In 2015 IEEE Pacific Visualization Symposium (PacificVis), pp. 271–278. IEEE, 2015. 9
work page 2015
- [2]
-
[3]
F. Al Sulti and K. Ohkitani. V ortex merger and topological changes in two-dimensional turbulence. Physical Review E, 86(1):016309, 2012
work page 2012
-
[4]
T. F. Banchoff. Critical points and curvature for embedded polyhedral surfaces. The American Mathematical Monthly, 1970
work page 1970
-
[5]
S. Biasotti, D. Giorgio, M. Spagnuolo, and B. Falcidieno. Reeb graphs for shape analysis and applications. Theoretical Computer Science, 2008
work page 2008
-
[6]
J. H. Bramble and A. H. Schatz. Higher order local accuracy by averaging in the finite element method. Mathematics of Computation, 31(137):94– 111, 1977
work page 1977
- [7]
-
[8]
C. D. Cantwell, D. Moxey, A. Comerford, A. Bolis, G. Rocco, G. Men- galdo, D. De Grazia, S. Yakovlev, J.-E. Lombard, D. Ekelschot, B. Jordi, H. Xu, Y . Mohamied, C. Eskilsson, B. Nelson, P. V os, C. Biotto, R. M. Kirby, and S. J. Sherwin. Nektar++ : An open-source spectral/hp element framework. Computer Physics Communications, 192:205–219, Jul 2015. doi: ...
-
[9]
H. Carr and J. Snoeyink. Representing interpolant topology for contour tree computation. In Topology-Based Methods in Visualization II , pp. 59–73. Springer, 2009
work page 2009
-
[10]
H. Carr, J. Snoeyink, and U. Axen. Computing contour trees in all dimen- sions. Computational Geometry, 24(2):75–94, 2003
work page 2003
-
[11]
H. Carr, J. Snoeyink, and M. van de Panne. Flexible isosurfaces: Simplify- ing and displaying scalar topology using the contour tree. Computational Geometry, 43(1):42 – 58, 2010
work page 2010
- [12]
-
[13]
G. Chen, K. Mischaikow, R. S. Laramee, and E. Zhang. Efficient morse decompositions of vector fields. IEEE Trans. Visualization and Computer Graphics, 14(4), 2008
work page 2008
-
[14]
B. Cockburn, M. Luskin, C.-W. Shu, and E. S ¨uli. Post-processing of Galerkin methods for hyperbolic problems. In Discontinuous Galerkin Methods, pp. 291–300. Springer, 2000
work page 2000
-
[15]
B. Cockburn, M. Luskin, C.-W. Shu, and E. S ´uli. Enhanced accuracy by post-processing for finite element methods for hyperbolic equations. Mathematics of Computation, 72(242):577–606, 2003
work page 2003
-
[16]
D. Cohen-Steiner, H. Edelsbrunner, and J. Harer. Stability of persistence diagrams. In Symp. on Comp. Geom., 2005
work page 2005
-
[17]
K. Cole-McLaughlin, H. Edelsbrunner, J. Harer, V . Natarajan, and V . Pas- cucci. Loops in Reeb graphs of 2-manifolds. In Symp. on Comp. Geom., pp. 344–350, 2003
work page 2003
-
[18]
M. Day, J. Bell, P.-T. Bremer, V . Pascucci, V . Beckner, and M. Lijew- ski. Turbulence effects on cellular burning structures in lean premixed hydrogen flames. Combustion and Flame, 156:1035–1045, 2009
work page 2009
-
[19]
L. De Floriani, U. Fugacci, F. Iuricich, and P. Magillo. Morse complexes for shape segmentation and homological analysis: discrete models and algorithms. Comp. Grap. For., 2015
work page 2015
-
[20]
S. E. Dillard, V . Natarajan, G. H. Weber, V . Pascucci, and B. Hamann. Topology-guided tessellation of quadratic elements. International Journal of Computational Geometry & Applications, 19(02):195–211, 2009
work page 2009
-
[21]
J. Docampo-S´anchez, J. K. Ryan, M. Mirzargar, and R. M. Kirby. Multi- dimensional filtering: Reducing the dimension through rotation. SIAM Journal on Scientific Computing, 39(5):A2179–A2200, 2017
work page 2017
-
[22]
H. Doraiswamy and V . Natarajan. Output-sensitive construction of reeb graphs. IEEE Trans. on Vis. and Comp. Graph., 2012
work page 2012
-
[23]
H. Doraiswamy and V . Natarajan. Computing reeb graphs as a union of contour trees. IEEE Trans. on Vis. and Comp. Graph., 2013
work page 2013
-
[24]
H. Edelsbrunner and J. Harer. Computational Topology: An Introduction. American Mathematical Society, 2009
work page 2009
-
[25]
H. Edelsbrunner, D. Letscher, and A. Zomorodian. Topological persistence and simplification. Disc. Compu. Geom., 2002
work page 2002
-
[26]
T. Etiene, L. Nonato, C. Scheidegger, J. Tienry, T. Peters, V . Pascucci, R. Kirby, and C. Silva. Topology verification for isosurface extraction. IEEE Transactions on Visualization and Computer Graphics, 18(6):952– 965, June 2012. doi: 10.1109/TVCG.2011.109
-
[27]
N. Fehn, W. A. Wall, and M. Kronbichler. Robust and efficient discontinu- ous Galerkin methods for under-resolved turbulent incompressible flows. Journal of Computational Physics, 372:667–693, 2018
work page 2018
-
[28]
R. Forman. A user’s guide to discrete Morse theory. Adv. in Math., 1998
work page 1998
-
[29]
D. Guenther, R. Alvarez-Boto, J. Contreras-Garcia, J.-P. Piquemal, and J. Tierny. Characterizing molecular interactions in chemical systems.IEEE Trans. on Vis. and Comp. Graph., 2014
work page 2014
-
[30]
C. Gueunet, P. Fortin, J. Jomier, and J. Tierny. Contour forests: Fast multi-threaded augmented contour trees. In IEEE LDAV, 2016
work page 2016
-
[31]
A. Gyulassy, P. Bremer, R. Grout, H. Kolla, J. Chen, and V . Pascucci. Stability of dissipation elements: A case study in combustion. Comp. Graph. For., 2014
work page 2014
-
[32]
A. Gyulassy, P. T. Bremer, B. Hamann, and V . Pascucci. A practical approach to morse-smale complex computation: Scalability and generality. IEEE Trans. on Vis. and Comp. Graph., 2008
work page 2008
-
[33]
A. Gyulassy, A. Knoll, K. Lau, B. Wang, P. Bremer, M. Papka, L. A. Curtiss, and V . Pascucci. Interstitial and interlayer ion diffusion geometry extraction in graphitic nanosphere battery materials. IEEE Trans. on Vis. and Comp. Graph., 2015
work page 2015
-
[34]
A. Gyulassy, V . Natarajan, M. Duchaineau, V . Pascucci, E. Bringa, A. Hig- ginbotham, and B. Hamann. Topologically Clean Distance Fields. IEEE Trans. on Vis. and Comp. Graph., 2007
work page 2007
- [35]
-
[36]
Z. W. H.T. Huynh and P. Vincent. High-order methods for computational fluid dynamics: A brief review of compact differential formulations on unstructured grids. Computers & Fluids, 98:209–220, 2014
work page 2014
-
[37]
A. Jallepalli, J. Docampo-S ´anchez, J. K. Ryan, R. Haimes, and R. M. Kirby. On the treatment of field quantities and elemental continuity in FEM solutions. IEEE Transactions on Visualization and Computer Graphics, 24(1):903–912, 2017
work page 2017
-
[38]
A. Jallepalli, R. Haimes, and R. M. Kirby. Adaptive characteristic length for L-SIAC filtering of fem data. Journal of Scientific Computing , pp. 1–22, 2018
work page 2018
-
[39]
A. Jallepalli and R. M. Kirby. Efficient algorithms for the line-siac filter. Journal of Scientific Computing, 2019
work page 2019
-
[40]
G. E. Karniadakis and S. J. Sherwin. Spectral/hp element methods for Computational Fluid Dynamics (Second Edition). Oxford University Press, 2005
work page 2005
- [41]
-
[42]
D. E. Laney, P. Bremer, A. Mascarenhas, P. Miller, and V . Pascucci. Understanding the structure of the turbulent mixing layer in hydrodynamic instabilities. IEEE Trans. on Vis. and Comp. Graph., 2006
work page 2006
-
[43]
X. Li, J. K. Ryan, R. M. Kirby, and C. Vuik. Smoothness-Increasing Accuracy-Conserving (SIAC) filters for derivative approximations of dis- continuous Galerkin (DG) solutions over nonuniform meshes and near boundaries. Journal of Computational and Applied Mathematics, 294:275– 296, 2016
work page 2016
-
[44]
J.-E. W. Lombard, D. Moxey, S. J. Sherwin, J. F. A. Hoessler, S. Dhanda- pani, and M. J. Taylor. Implicit large-eddy simulation of a wingtip vortex. AIAA Journal, 54(2):506–518, 2016
work page 2016
-
[45]
H. Mirzaee, L. Ji, J. K. Ryan, and R. M. Kirby. Smoothness-increasing accuracy-conserving (SIAC) postprocessing for discontinuous Galerkin solutions over structured triangular meshes. SIAM Journal on Numerical Analysis, 49(5):1899–1920, 2011
work page 1920
-
[46]
H. Mirzaee, J. King, J. K. Ryan, and R. M. Kirby. Smoothness-increasing accuracy-conserving filters for discontinuous Galerkin solutions over un- structured triangular meshes. SIAM Journal on Scientific Computing , 35(1):A212–A230, 2013
work page 2013
-
[47]
H. Mirzaee, J. K. Ryan, and R. M. Kirby. Efficient implementation of smoothness-increasing accuracy-conserving (SIAC) filters for discontinu- ous Galerkin solutions. Journal of Scientific Computing, 52(1):85–112, 2012
work page 2012
-
[48]
H. Mirzaee, J. K. Ryan, and R. M. Kirby. Smoothness-increasing accuracy- conserving (SIAC) filters for discontinuous Galerkin solutions: Applica- tion to structured tetrahedral meshes. Journal of Scientific Computing , 58(3):690–704, 2014
work page 2014
-
[49]
M. Mirzargar, A. Jallepalli, J. K. Ryan, and R. M. Kirby. Hexagonal smoothness-increasing accuracy-conserving filtering. Journal of Scientific Computing, 73(2-3):1072–1093, 2017
work page 2017
- [50]
-
[51]
D.-M. Nguyen and J. Peters. Nonuniform discontinuous Galerkin filters via shift and scale. SIAM Journal on Numerical Analysis, 54(3):1401–1422, 2016
work page 2016
-
[52]
G. Norgard and P.-T. Bremer. Robust computation of morse–smale complexes of bilinear functions. Computer Aided Geometric Design , 30(6):577–587, 2013
work page 2013
-
[53]
G. Nucha, G.-P. Bonneau, S. Hahmann, and V . Natarajan. Computing con- tour trees for 2d piecewise polynomial functions. In Computer Graphics Forum, vol. 36, pp. 23–33. Wiley Online Library, 2017
work page 2017
-
[54]
S. Parsa. A deterministic o(m log m) time algorithm for the reeb graph. In Symp. on Comp. Geom., 2012
work page 2012
-
[55]
V . Pascucci and K. Cole-McLaughlin. Parallel computation of the topology of level sets. Algorithmica, 38(1):249–268, 2004
work page 2004
-
[56]
V . Pascucci, G. Scorzelli, P. T. Bremer, and A. Mascarenhas. Robust on-line computation of Reeb graphs: simplicity and speed. ACM Trans. on Graph., 2007
work page 2007
- [57]
-
[58]
J. K. Ryan, X. Li, R. M. Kirby, and K. Vuik. One-sided position-dependent smoothness-increasing accuracy-conserving (SIAC) filtering over uniform and non-uniform meshes. Journal of Scientific Computing, 64(3):773–817, 2015
work page 2015
-
[59]
N. Shivashankar, P. Pranav, V . Natarajan, R. van de Weygaert, E. P. Bos, and S. Rieder. Felix: A topology based framework for visual exploration of cosmic filaments. IEEE Trans. on Vis. and Comp. Graph., 2016. http: //vgl.serc.iisc.ernet.in/felix/index.html
work page 2016
-
[60]
A. Sidot. Spectral/hp Element Methods Applied To Vortex Structures In A Low Incidence Delta Wing Wake Compared To Experiments. M.S. Thesis in Aeronautics, Imperial College London, September 2017
work page 2017
-
[61]
B. S. Sohn and C. L. Bajaj. Time varying contour topology. IEEE Trans. on Vis. and Comp. Graph., 2006
work page 2006
-
[62]
T. Sousbie. The persistent cosmic web and its filamentary structure: Theory and implementations. Royal Astronomical Society, 2011. http: //www2.iap.fr/users/sousbie/web/html/indexd41d.html
work page 2011
-
[63]
M. Steffen, S. Curtis, R. M. Kirby, and J. K. Ryan. Investigation of smoothness-increasing accuracy-conserving filters for improving stream- line integration through discontinuous fields. IEEE Transactions on Visu- alization and Computer Graphics, 14(3):680–692, 2008
work page 2008
-
[64]
S. Tarasov and M. Vyali. Construction of contour trees in 3d in o(n log n) steps. In Symp. on Comp. Geom., 1998
work page 1998
-
[65]
D. M. Thomas and V . Natarajan. Multiscale symmetry detection in scalar fields by clustering contours. IEEE Trans. on Vis. and Comp. Graph. , 2014
work page 2014
- [66]
- [67]
-
[68]
J. Tierny and V . Pascucci. Generalized topological simplification of scalar fields on surfaces. IEEE Trans. on Vis. and Comp. Graph., 2012
work page 2012
-
[69]
X. Tricoche, C. Garth, G. Kindlmann, E. Deines, G. Scheuermann, M. Ruetten, and C. Hansen. Visualization of intricate flow structures for vortex breakdown analysis. In IEEE Visualization 2004, pp. 187–194. IEEE, 2004
work page 2004
-
[70]
M. van Kreveld, R. van Oostrum, C. Bajaj, V . Pasucci, and D. Schikore. Contour trees and small seed sets for isosurface traversal. In Symp. on Comp. Geom., 1997
work page 1997
-
[71]
P. van Slingerland, J. K. Ryan, and C. Vuik. Position-dependent smoothness-increasing accuracy-conserving (SIAC) filtering for improv- ing discontinuous Galerkin solutions. SIAM Journal on Scientific Comput- ing, 33(2):802–825, 2011
work page 2011
-
[72]
D. Walfisch, J. K. Ryan, R. M. Kirby, and R. Haimes. One-sided smoothness-increasing accuracy-conserving filtering for enhanced stream- line integration through discontinuous fields.Journal of Scientific Com- puting, 38(2):164–184, 2009
work page 2009
- [73]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.