High-Order Methods for Computational Fluid Dynamics: A Brief Review of Compact Differential Formulations on Unstructured Grids

Popular high-order schemes with compact stencils for Computational Fluid Dynamics (CFD) include Discontinuous Galerkin (DG), Spectral Difference (SD), and Spectral Volume (SV) methods. The recently proposed Flux Reconstruction (FR) approach or Correction Procedure using Reconstruction (CPR) is based on a differential formulation and provides a unifying framework for these high-order schemes. Here we present a brief review of recent developments for the FR/CPR schemes as well as some pacing items

A Survey of the Isentropic Euler Vortex Problem Using High-Order Methods

The flux reconstruction (FR) method offers a simple, efficient, and easy to implement method, and it has been shown to equate to a differential approach to discontinuous Galerkin (DG) methods. The FR method is also accurate to an arbitrary order and the isentropic Euler vortex problem is used here to empirically verify this claim. This problem is widely used in computational fluid dynamics (CFD) to verify the accuracy of a given numerical method due to its simplicity and known exact solution at any given time. While verifying our FR solver, multiple obstacles emerged that prevented us from achieving the expected order of accuracy over short and long amounts of simulation time. It was found that these complications stemmed from a few overlooked details in the original problem definition combined with the FR and DG methods achieving high-accuracy with minimal dissipation. This paper is intended to consolidate the many versions of the vortex problem found in literature and to highlight some of the consequences if these overlooked details remain neglected

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

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:

De-Aliasing Through Over-Integration Applied to the Flux Reconstruction and Discontinuous Galerkin Methods

High-order methods are quickly becoming popular for turbulent flows as the amount of computer processing power increases. The flux reconstruction (FR) method presents a unifying framework for a wide class of high-order methods including discontinuous Galerkin (DG), Spectral Difference (SD), and Spectral Volume (SV). It offers a simple, efficient, and easy way to implement nodal-based methods that are derived via the differential form of the governing equations. Whereas high-order methods have enjoyed recent success, they have been known to introduce numerical instabilities due to polynomial aliasing when applied to under-resolved nonlinear problems. Aliasing errors have been extensively studied in reference to DG methods; however, their study regarding FR methods has mostly been limited to the selection of the nodal points used within each cell. Here, we extend some of the de-aliasing techniques used for DG methods, primarily over-integration, to the FR framework. Our results show that over-integration does remove aliasing errors but may not remove all instabilities caused by insufficient resolution (for FR as well as DG)

Higher-order Multi-dimensional Limiting Strategy for Correction Procedure via Reconstruction

AIAA SciTech, 52nd Aerospace Sciences Meeting 13-17 January 2014, National Harbor, MarylandHigher-order multi-dimensional limiting Process (MLP) [J. S. Park and C. Kim, Higher-order Multi-dimensional Limiting Strategy for Discontinuous Galerkin Methods in Compressible Inviscid and Viscous Flows, Comp. & Fluids, In press] is improved and applied to correction procedure via reconstruction (CPR) on unstructured grids. MLP, which has been originally developed in nite volume method (FVM), provides an accurate, robust and ecient oscillation-control mechanism in multiple dimensions for linear reconstruction. This limiting philosophy can be hierarchically extended into higher-order Pn reconstruction. The resulting algorithms, called the hierarchical MLP, facilitate the accurate capturing of detailed flow structures in both continuous and discontinuous regions. This algorithm has been developed in 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 nal accuracy is determined by the projection procedure and MLP limiting step. Through extensive numerical analyses and computations, it is demonstrated that the proposed limiting approach yields the desired accuracy and outstanding performances in resolving compressible inviscid and viscous flow features.This work is supported by NSL (National Space Lab.) program through the National Research Foundation of Korea funded by the Ministry of Education, Science and Technology (Grant 20090091724), by the third stage of the Brain Korea 21 Plus Project in 2013, and by Korea Ministry of Land, Transport and Maritime A airs as Haneul Project. This work is also supported by the KISTI Supercomputing Center (KSC-2012-C3-40) and by PLSI supercomputing resources of KISTI.OAIID:oai:osos.snu.ac.kr:snu2014-01/104/0000004648/1SEQ:1PERF_CD:SNU2014-01EVAL_ITEM_CD:104USER_ID:0000004648ADJUST_YN:NEMP_ID:A001138DEPT_CD:446CITE_RATE:0FILENAME:6%2e2014-0772.pdfDEPT_NM:기계항공공학부EMAIL:[email protected]: