Visualization for high-order discontinuous Galerkin CFD results

• Main Author: Walfisch, David
• Other Authors: Robert Haimes.
• Format: Others
• Language: English
• Published: Massachusetts Institute of Technology 2008
• Subjects: Aeronautics and Astronautics.
• Online Access: http://hdl.handle.net/1721.1/42196

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.