77 research outputs found

    Non-Oscillatory Hierarchical Reconstruction for Central and Finite Volume Schemes

    Get PDF
    This is the continuation of the paper "central discontinuous Galerkin methods on overlapping cells with a non-oscillatory hierarchical reconstruction" by the same authors. The hierarchical reconstruction introduced therein is applied to central schemes on overlapping cells and to nite volume schemes on non-staggered grids. This takes a new nite volume approach for approximating non-smooth solutions. A critical step for high order nite volume schemes is to reconstruct a nonoscillatory high degree polynomial approximation in each cell out of nearby cell averages. In the paper this procedure is accomplished in two steps: first to reconstruct a high degree polynomial in each cell by using e.g., a central reconstruction, which is easy to do despite the fact that the reconstructed polynomial could be oscillatory; then to apply the hierarchical reconstruction to remove the spurious oscillations while maintaining the high resolution. All numerical computations for systems of conservation laws are performed without characteristic decomposition. In particular, we demonstrate that this new approach can generate essentially non-oscillatory solutions even for 5th order schemes without characteristic decomposition.The research of Y. Liu was supported in part by NSF grant DMS-0511815. The research of C.-W. Shu was supported in part by the Chinese Academy of Sciences while this author was visiting the University of Science and Technology of China (grant 2004-1-8) and the Institute of Computational Mathematics and Scienti c/Engineering Computing. Additional support was provided by ARO grant W911NF-04-1-0291 and NSF grant DMS-0510345. The research of E. Tadmor was supported in part by NSF grant 04-07704 and ONR grant N00014-91-J-1076. The research of M. Zhang was supported in part by the Chinese Academy of Sciences grant 2004-1-8

    A reconstructed discontinuous Galerkin method based on a Hierarchical WENO reconstruction for compressible flows on tetrahedral grids

    Get PDF
    A reconstructed discontinuous Galerkin (RDG) method based on a hierarchical WENO reconstruction, termed HWENO (P1P2) in this paper, designed not only to enhance the accuracy of discontinuous Galerkin methods but also to ensure the nonlinear stability of the RDG method, is presented for solving the compressible Euler equations on tetrahedral grids. In this HWENO (P1P2) method, a quadratic polynomial solution (P-2) is first reconstructed using a Hermite WENO reconstruction from the underlying linear polynomial (P-1) discontinuous Galerkin solution to ensure the linear stability of the RDG method and to improve the efficiency of the underlying DG method. By taking advantage of handily available and yet invaluable information, namely the derivatives in the DG formulation, the stencils used in the reconstruction involve only von Neumann neighborhood (adjacent face-neighboring cells) and thus are compact. The first derivatives of the quadratic polynomial solution are then reconstructed using a WENO reconstruction in order to eliminate spurious oscillations in the vicinity of strong discontinuities, thus ensuring the nonlinear stability of the RDG method. The developed HWENO (P1P2) method is used to compute a variety of flow problems on tetrahedral meshes to demonstrate its accuracy, robustness, and non-oscillatory property. The numerical experiments indicate that the HWENO (P1P2) method is able to capture shock waves within one cell without any spurious oscillations, and achieve the designed third-order of accuracy: one order accuracy higher than the underlying DG method

    WENO schemes on unstructured meshes using a relaxed a posteriori MOOD limiting approach

    Get PDF
    In this paper a relaxed formulation of the a posteriori Multi-dimensional Optimal Order Detection (MOOD) limiting approach is introduced for weighted essentially non-oscillatory (WENO) finite volume schemes on unstructured meshes. The main goal is to minimise the computational footprint of the MOOD limiting approach by employing WENO schemes—by virtue of requiring a smaller number of cells to reduce their order of accuracy compared to an unlimited scheme. The key characteristic of the present relaxed MOOD formulation is that the Numerical Admissible Detector (NAD) is not uniquely defined for all orders of spatial accuracy, and it is relaxed when reaching a 2nd-order of accuracy. The augmented numerical schemes are applied to the 2D unsteady Euler equations for a multitude of test problems including the 2D vortex evolution, cylindrical explosion, double-Mach reflection, and an implosion. It is observed that in many events, the implemented MOOD paradigm manages to preserve symmetry of the forming structures in simulations, an improvement comparing to the non-MOOD limited counterparts which cannot be easily obtained due to the multi-dimensional reconstruction nature of the schemes

    Multi-dimensional Limiting Strategy for Higher-order CFD Methods - Progress and Issue (Invited)

    Get PDF
    The present paper deals with the progress of multi-dimensional limiting process (MLP) and discuss the issues for further improvements. MLP, which has been originally developed in finite volume method (FVM), provides an accurate, robust and efficient oscillationcontrol mechanism in multiple dimensions for linear reconstruction. This limiting philosophy can be hierarchically extended into higher-order Pn approximation or reconstruction. The resulting algorithm, called the hierarchical MLP, facilitates the capturing of detailed flow structures while maintaining the formal order-of-accuracy in smooth region and providing accurate non-oscillatory solutions across discontinuous region. This algorithm has been developed within the modal DG framework, but it also can be formulated into a nodal framework, most notably the CPR framework. Troubled-cells are detected by applying the MLP concept, and the final accuracy is determined by the projection procedure and the hierarchical MLP limiting step. Through extensive numerical analyses and computations ranging from scalar conservation laws to fluid systems, it is demonstrated that the proposed limiting approach yields the outstanding performances in capturing compressible inviscid and viscous flow features. Further issues are also mentioned to improve and extend the current approach for higher-order simulations of high-Reynolds number compressible flows.Authors appreciate the financial supports by the EDISON program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (NRF-2011-0020559) and by NSL (National Space Laboratory) program through the National Research Foundation of Korea funded by the Ministry of Education, Science and Technology (NRF-2014M1A3A3A02034856). This work is also partially supported by the RoK ST&R project of Lockheed Martin Corporation. Authors also appreciate the computing resources provided by the KISTI Supercomputing Center(KSC-2014-C3-054).OAIID:RECH_ACHV_DSTSH_NO:420150000004648007RECH_ACHV_FG:RR00200003ADJUST_YN:EMP_ID:A001138CITE_RATE:FILENAME:6.2015-3199.pdfDEPT_NM:기계항공공학부EMAIL:[email protected]_YN:FILEURL:https://srnd.snu.ac.kr/eXrepEIR/fws/file/a984d649-4b23-435b-adc9-df9aa0c8aa46/linkCONFIRM:

    A high-order finite-volume method for atmospheric flows on unstructured grids

    Get PDF
    This paper presents an extension of a Weighted Essentially Non-Oscillatory (WENO) type schemes for the compressible Euler equations on unstructured meshes for stratified atmospheric flows. The schemes could be extended for regional and global climate models dynamical cores. Their potential lies in their simplicity; accuracy; robustness; non-oscillatory properties; versatility in handling any type of grid topology; computational and parallel efficiency. Their defining characteristic is a non-linear combination of a series of high-order reconstruction polynomials arising from a series of reconstruction stencils. In the present study an explicit Strong Stability Preserving (SSP) Runge-Kutta 3rd-order method is employed for time advancement. The WENO schemes (up to 5th-order) are applied to the two dimensional and three dimensional test cases: a 2D rising thermal bubble; the 2D density current and the 3D Robert smooth bubble. The parallel performance of the schemes in terms of scalability and efficiency is also assessed

    A discontinuous Galerkin method based on a Taylor basis for the compressible flows on arbitrary grids

    Get PDF
    A discontinuous Galerkin method based on a Taylor basis is presented for the solution of the compressible Euler equations on arbitrary grids. Unlike the traditional discontinuous Galerkin methods, where either standard Lagrange finite element or hierarchical node-based basis functions are used to represent numerical polynomial solutions in each element, this DG method represents the numerical polynomial solutions using a Taylor series expansion at the centroid of the cell. Consequently, this formulation is able to provide a unified framework, where both cell-centered and vertex-centered finite volume schemes can be viewed as special cases of this discontinuous Galerkin method by choosing reconstruction schemes to compute the derivatives, offer the insight why the DG methods are a better approach than the finite volume methods based on either TVD/MUSCL reconstruction or essentially non-oscillatory (ENO)/weighted essentially non-oscillatory (WENO) reconstruction, and has a number of distinct, desirable, and attractive features, which can be effectively used to address some of shortcomings of the DG methods. The developed method is used to compute a variety of both steady-state and time-accurate flow problems on arbitrary grids. The numerical results demonstrated the superior accuracy of this discontinuous Galerkin method in comparison with a second order finite volume method and a third-order WENO method, indicating its promise and potential to become not just a competitive but simply a superior approach than its finite volume and ENO/WENO counterparts for solving flow problems of scientific and industrial interest
    corecore