Visualization for high-order discontinuous Galerkin CFD results

Thesis (S.M.)--Massachusetts Institute of Technology, Dept. of Aeronautics and Astronautics, 2007. === Includes bibliographical references (p. 69-71). === This thesis demonstrates a technique that reduces the perceived error introduced in the visualization of higher-order DG solutions with tradition...

Full description

Bibliographic Details
Main Author: Walfisch, David
Other Authors: Robert Haimes.
Format: Others
Language:English
Published: Massachusetts Institute of Technology 2008
Subjects:
Online Access:http://hdl.handle.net/1721.1/42196
Description
Summary:Thesis (S.M.)--Massachusetts Institute of Technology, Dept. of Aeronautics and Astronautics, 2007. === Includes bibliographical references (p. 69-71). === This thesis demonstrates a technique that reduces the perceived error introduced in the visualization of higher-order DG solutions with traditional tools. Most visualization techniques work by using the solution at the mesh vertices and linearly interpolating to color the interior of the element. For higher-order solutions (where the visualization mesh is the same as the solution mesh) this method renders the high-order solution linear. By inserting new nodes at element quadrature points where the difference between the actual solution and the linear rendering is above a user-defined tolerance, additional elements are created and used for the visualization process. In order to keep the counts low for this new mesh, after each insertion a local rearrangement is performed to readapt the parent element so that the total visualization error is reduced. The method introduced here has many advantages over isotropic adaptation used by some higher-order visualization techniques. Isotropic adaptation adapts all the elements regardless of error, thus creating a higher total element count and therefore requiring more memory and rendering time. In some cases isotropic elements are not ideal in representing the solution (ie: boundary layers, shocks, wakes, etc.). Lastly, by providing an option to define the maximum visualization error allows the user to specify how close the visualized solution is to the actual calculated one (at the expense of a denser visualization mesh). Second, this work introduces a new method to apply an accuracy maintaining post-processor on DG vector fields to improve on the standard streamlining algorithms. Time integration methods do not work accurately and may even fail on discontinuos fields. The post-processor smoothens the field and eliminates the discontinuity between elements thus resulting in more accurate streamlines. To keep down the computational cost of the method, the post-processing is done in a one dimensional manner along the streamline. === by David Walfisch. === S.M.