Compact-Reconstruction Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws

A new class of non-linear compact interpolation schemes is introduced in this dissertation that have a high spectral resolution and are non-oscillatory across discontinuities. The Compact-Reconstruction Weighted Essentially Non-Oscillatory (CRWENO) schemes use a solution-dependent combination of lower-order compact schemes to yield a high-order accurate, non-oscillatory scheme. Fifth-order accurate CRWENO schemes are constructed and their numerical properties are analyzed. These schemes have lower absolute errors and higher spectral resolution than the WENO scheme of the same order. The schemes are applied to scalar conservation laws and the Euler equations of fluid dynamics. The order of convergence and the higher accuracy of the CRWENO schemes are verified for smooth solutions. Significant improvements are observed in the resolution of discontinuities and extrema as well as the preservation of flow features over large convection distances. The computational cost of the CRWENO schemes is assessed and the reduced error in the solution outweighs the additional expense of the implicit scheme, thus resulting in higher numerical efficiency. This conclusion extends to the reconstruction of conserved and primitive variables for the Euler equations, but not to the characteristic-based reconstruction. Further improvements are observed in the accuracy and resolution of the schemes with alternative formulations for the non-linear weights. The CRWENO schemes are integrated into a structured, finite-volume Navier-Stokes solver and applied to problems of practical relevance. Steady and unsteady flows around airfoils are solved to validate the scheme for curvi-linear grids, as well as overset grids with relative motion. The steady flow around a three-dimensional wing and the unsteady flow around a full-scale rotor are solved. It is observed that though lower-order schemes suffice for the accurate prediction of aerodynamic forces, the CRWENO scheme yields improved resolution of near-blade and wake flow features, including boundary and shear layers, and shed vortices. The high spectral resolution, coupled with the non-oscillatory behavior, indicate their suitability for the direct numerical simulation of compressible turbulent flows. Canonical flow problems -- the decay of isotropic turbulence and the shock-turbulence interaction -- are solved. The CRWENO schemes show an improved resolution of the higher wavenumbers and the small-length-scale flow features that are characteristic of turbulent flows. Overall, the CRWENO schemes show significant improvements in resolving and preserving flow features over a large range of length scales due to the higher spectral resolution and lower dissipation and dispersion errors, compared to the WENO schemes. Thus, these schemes are a viable alternative for the numerical simulation of compressible, turbulent flows

A Space-time Smooth Artificial Viscosity Method For Nonlinear Conservation Laws

We introduce a new methodology for adding localized, space-time smooth, artificial viscosity to nonlinear systems of conservation laws which propagate shock waves, rarefactions, and contact discontinuities, which we call the $C$-method. We shall focus our attention on the compressible Euler equations in one space dimension. The novel feature of our approach involves the coupling of a linear scalar reaction-diffusion equation to our system of conservation laws, whose solution $C(x,t)$ is the coefficient to an additional (and artificial) term added to the flux, which determines the location, localization, and strength of the artificial viscosity. Near shock discontinuities, $C(x,t)$ is large and localized, and transitions smoothly in space-time to zero away from discontinuities. Our approach is a provably convergent, spacetime-regularized variant of the original idea of Richtmeyer and Von Neumann, and is provided at the level of the PDE, thus allowing a host of numerical discretization schemes to be employed. We demonstrate the effectiveness of the $C$-method with three different numerical implementations and apply these to a collection of classical problems: the Sod shock-tube, the Osher-Shu shock-tube, the Woodward-Colella blast wave and the Leblanc shock-tube. First, we use a classical continuous finite-element implementation using second-order discretization in both space and time, FEM-C. Second, we use a simplified WENO scheme within our $C$-method framework, WENO-C. Third, we use WENO with the Lax-Friedrichs flux together with the $C$-equation, and call this WENO-LF-C. All three schemes yield higher-order discretization strategies, which provide sharp shock resolution with minimal overshoot and noise, and compare well with higher-order WENO schemes that employ approximate Riemann solvers, outperforming them for the difficult Leblanc shock tube experiment.Comment: 34 pages, 27 figure

ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement

We present the first high order one-step ADER-WENO finite volume scheme with Adaptive Mesh Refinement (AMR) in multiple space dimensions. High order spatial accuracy is obtained through a WENO reconstruction, while a high order one-step time discretization is achieved using a local space-time discontinuous Galerkin predictor method. Due to the one-step nature of the underlying scheme, the resulting algorithm is particularly well suited for an AMR strategy on space-time adaptive meshes, i.e.with time-accurate local time stepping. The AMR property has been implemented 'cell-by-cell', with a standard tree-type algorithm, while the scheme has been parallelized via the Message Passing Interface (MPI) paradigm. The new scheme has been tested over a wide range of examples for nonlinear systems of hyperbolic conservation laws, including the classical Euler equations of compressible gas dynamics and the equations of magnetohydrodynamics (MHD). High order in space and time have been confirmed via a numerical convergence study and a detailed analysis of the computational speed-up with respect to highly refined uniform meshes is also presented. We also show test problems where the presented high order AMR scheme behaves clearly better than traditional second order AMR methods. The proposed scheme that combines for the first time high order ADER methods with space--time adaptive grids in two and three space dimensions is likely to become a useful tool in several fields of computational physics, applied mathematics and mechanics.Comment: With updated bibliography informatio