Showing papers by "Mark H. Carpenter published in 1995"

The theoretical accuracy of Runge-Kutta time discretizations for the initial boundary value problem: a study of the boundary error

Abstract: The conventional method of imposing time dependent boundary conditions for Runge-Kutta (RK) time advancement reduces the formal accuracy of the space-time method to first order locally, and second order globally, independently of the spatial operator. This counter intuitive result is analyzed in this paper. Two methods of eliminating this problem are proposed for the linear constant coefficient case: 1) impose the exact boundary condition only at the end of the complete RK cycle, 2) impose consistent intermediate boundary conditions derived from the physical boundary condition and its derivatives. The first method, while retaining the RK accuracy in all cases , results in a scheme with much reduced CFL condition, rendering the RK scheme less attractive. The second method retains the same allowable time step as the periodic problem. However it is a general remedy only for the linear case. For non-linear hyperbolic equations the second method is effective only for for RK schemes of third order accuracy or less. Numerical studies are presented to verify the efficacy of each approach.

Spectral methods on arbitrary grids

Abstract: Stable and spectrally accurate numerical methods are constructed on arbitrary grids for partial differential equations. These new methods are equivalent to conventional spectral methods but do not rely on specific grid distributions. Specifically, we show how to implement Legendre Galerkin, Legendre collocation, and Laguerre Galerkin methodology on arbitrary grids.

Spatial direct numerical simulation of high-speed boundary-layer flows part I: Algorithmic considerations and validation

Abstract: A highly accurate algorithm for the direct numerical simulation (DNS) of spatially evolving high-speed boundary-layer flows is described in detail and is carefully validated. To represent the evolution of instability waves faithfully, the fully explicit scheme relies on non-dissipative high-order compact-difference and spectral collocation methods. Several physical, mathematical, and practical issues relevant to the simulation of high-speed transitional flows are discussed. In particular, careful attention is paid to the implementation of inflow, outflow, and far-field boundary conditions. Four validation cases are presented, in which comparisons are made between DNS results and results obtained from either compressible linear stability theory or from the parabolized stability equation (PSE) method, the latter of which is valid for nonparallel flows and moderately nonlinear disturbance amplitudes. The first three test cases consider the propagation of two-dimensional second-mode disturbances in Mach 4.5 flat-plate boundary-layer flows. The final test case considers the evolution of a pair of oblique second-mode disturbances in a Mach 6.8 flow along a sharp cone. The agreement between the fundamentally different PSE and DNS approaches is remarkable for the test cases presented.