Pii: S0168-9274(99)00067-7

نویسندگان

  • Anne Gelb
  • Eitan Tadmor
چکیده

In this paper we construct, analyze and implement a new procedure for the spectral approximations of nonlinear conservation laws. It is well known that using spectral methods for nonlinear conservation laws will result in the formation of the Gibbs phenomenon once spontaneous shock discontinuities appear in the solution. These spurious oscillations will in turn lead to loss of resolution and render the standard spectral approximations unstable. The Spectral Viscosity (SV-) method (Tadmor, 1989) was developed to stabilize the spectral method by adding a spectrally small amount of high-frequencies diffusion carried out in the dual space. The resulting SV-approximation is stable without sacrificing spectral accuracy. The SV-method recovers a spectrally accurate approximation to the projection of the entropy solution; the exact projection, however, is at best a first order approximation to the exact solution as a result of the formation of the shock discontinuities. The issue of spectral resolution is addressed by post-processing the SV-solution to remove the spurious oscillations at the discontinuities, as well as increase the first-order—O(1/N) accuracy away from the shock discontinuities. Successful post-processing methods have been developed to eliminate the Gibbs phenomenon and recover spectral accuracy for the SV-approximation. However, such reconstruction methods require a priori knowledge of the locations of the shock discontinuities. Therefore, the detection of these discontinuities is essential to obtain an overall spectrally accurate solution. To this end, we employ the recently constructed enhanced edge detectors based on appropriate concentration factors (Gelb and Tadmor, 1999). Once the edges of these discontinuities are identified, we can utilize a post-processing reconstruction method, and show that the post-processed SV-solution recovers the exact entropy solution with remarkably high-resolution. We apply our new numerical method, the Enhanced SV-method, to two numerical examples, the scalar periodic Burgers’ equation and the one-dimensional system of Euler equations of gas dynamics. Both approximations exhibit high accuracy and resolution to the exact entropy solution.  2000 IMACS. Published by Elsevier Science B.V. All rights reserved.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Pii: S0168-9274(99)00020-3

We compare several methods for sensitivity analysis of differential–algebraic equations (DAEs). Computational complexity, efficiency and numerical conditioning issues are discussed. Numerical results for a chemical kinetics problem arising in model reduction are presented.  2000 IMACS. Published by Elsevier Science B.V. All rights reserved.

متن کامل

Pii: S0168-9274(99)00131-2

We study numerical integrators that contract phase space volume even when the ODE does so at an arbitrarily small rate. This is done by a splitting into two-dimensional contractive systems. We prove a sufficient condition for Runge–Kutta methods to have the appropriate contraction property for these two-dimensional systems; the midpoint rule is an example.  2000 IMACS. Published by Elsevier Sc...

متن کامل

Pii: S0168-9274(99)00082-3

A fast Chebyshev–Fourier algorithm for Poisson-type equations in polar geometries is presented in this paper. The new algorithm improves upon the algorithm of Jie Shen (1997), by taking advantage of the odd–even parity of the Fourier expansion in the azimuthal direction, and it is shown to be more efficient in terms of CPU and memory.  2000 IMACS. Published by Elsevier Science B.V. All rights ...

متن کامل

Pii: S0168-9274(99)00114-2

This paper concerns the iterative solution of the linear system arising from the Chebyshev–collocation approximation of second-order elliptic equations and presents an optimal multigrid preconditioner based on alternating line Gauss–Seidel smoothers for the corresponding stiffness matrix of bilinear finite elements on the Chebyshev–Gauss–Lobatto grid.  2000 IMACS. Published by Elsevier Science...

متن کامل

Pii: S0168-9274(00)00050-7

The aim of this work is to present a nonstandard linear finite element method for a planar elasticity problem. The error for the solution computed with this method is estimated with respect to H 1 ×H 1-norm and second-order convergence is shown.  2001 IMACS. Published by Elsevier Science B.V. All rights reserved.

متن کامل

Pii: S0168-9274(99)00148-8

We study the error propagation of time integrators of solitary wave solutions for the regularized long wave equation, ut +ux+ 2 (u)x −uxxt = 0, by using a geometric interpretation of these waves as relative equilibria. We show that the error growth is linear for schemes that preserve invariant quantities of the problem and quadratic for ‘nonconservative’ methods. Numerical experiments are prese...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2000