Parallel computation of unsteady compressible flows with the EDICT

Abstract

Recently, the Enhanced-Discretization Interface-Capturing Technique (EDICT) was introduced for simulation of unsteady flow problems with interfaces such as two-fluid and free-surface flows. The EDICT yields increased accuracy in representing the interface. Here we extend the EDICT to simulation of unsteady viscous compressible flows with boundary/shear layers and shock/expansion waves. The purpose is to increase the accuracy in selected regions of the computational domain. An error indicator is used to identify these regions that need enhanced discretization. Stabilized finite-element formulations are employed to solve the Navier-Stokes equations in their conservation law form. The finite element functions corresponding to enhanced discretization are designed to have two components, with each component coming from a different level of mesh refinement over the same computational domain. The primary component comes from a base mesh. A subset of the elements in this base mesh are identified for enhanced discretization by utilizing the error indicator. A secondary, more refined, mesh is constructed by patching together the second-level meshes generated over this subset of elements, and the second component of the functions comes from this mesh. The subset of elements in the base mesh that form the secondary mesh may change from one time level to other depending on the distribution of the error in the computations. Using a parallel implementation of this EDICT-based method, we apply it to test problems with shocks and boundary layers, and demonstrate that this method can be used very effectively to increase the accuracy of the base finite element formulation.