Spectral Dynamics of the Velocity Gradient Field in Restricted Flows

We study the velocity gradients of the fundamental Eulerian equation, $\partial_t u +u\cdot \nabla u=F$, which shows up in different contexts dictated by the different modeling of $F$'s. To this end we utilize a basic description for the spectral dynamics of $\nabla u$, expressed in terms of the (possibly complex) eigenvalues, $\lambda=\lambda(\nabla u)$, which are shown to be governed by the Ricatti-like equation $\lambda_t+u\cdot \nabla\lambda+\lambda^2=$. We address the question of the time regularity of four prototype models associated with different forcing $F$. Using the spectral dynamics as our essential tool in these investigations, we obtain a simple form of a critical threshold for the linear damping model and we identify the 2D vanishing viscosity limit for the viscous irrotational dusty medium model. Moreover, for the $n$-dimensional restricted Euler equations we obtain $[n/2]+1$ global invariants, interesting for their own sake, which enable us to precisely characterize the local topology at breakdown time, extending previous studies in the $n=3$-dimensional case. Finally, as a forth model we introduce the $n$-dimensional restricted Euler-Poisson (REP)system, identifying a set of $[n/2]$ global invariants, which in turn yield (i) sufficient conditions for finite time breakdown, and (ii) characterization of a large class of 2-dimensional initial configurations leading to global smooth solutions. Consequently, the 2D restricted Euler-Poisson equations are shown to admit a critical threshold

The numerical viscosity of entropy stable schemes for systems of conservation laws

Discrete approximations to hyperbolic systems of conservation laws are studied. The amount of numerical viscosity present in such schemes, is quantified and related to their entropy stability by means of comparison. To this end, conservative schemes which are also entropy conservative are constructed. These entropy conservative schemes enjoy second-order accuracy; moreover, they admit a particular interpretation within the finite-element frameworks, and hence can be formulated on various mesh configurations. It is then shown that conservative schemes are entropy stable if and only if they contain more viscosity than the mentioned above entropy conservative ones

Entropy functions for symmetric systems of conservation laws

It is shown that symmetric systems of conservation laws are equipped with a one-parameter family of entropy functions. A simple symmetrizability criterion is used

Convenient stability criteria for difference approximations of hyperbolic initial-boundary value problems

The purpose of this paper is to achieve more versatile, convenient stability criteria for a wide class of finite-difference approximations to initial boundary value problems associated with the hyperbolic system u sub t = au sub x + Bu + f in the quarter-plane x greater than or equal to 0, t greater than or equal to 0. With these criteria, stability is easily established for a large number of examples, thus incorporating and generalizing many of the cases studied in recent literature

Recovering pointwise values of discontinuous data within spectral accuracy

The pointwise values of a function, f(x), can be accurately recovered either from its spectral or pseudospectral approximations, so that the accuracy solely depends on the local smoothness of f in the neighborhood of the point x. Most notably, given the equidistant function grid values, its intermediate point values are recovered within spectral accuracy, despite the possible presence of discontinuities scattered in the domain. (Recall that the usual spectral convergence rate decelerates otherwise to first order, throughout). To this end, a highly oscillatory smoothing kernel is employed in contrast to the more standard positive unit-mass mollifiers. In particular, post-processing of a stable Fourier method applied to hyperbolic equations with discontinuous data, recovers the exact solution modulo a spectrally small error. Numerical examples are presented

On the convergence of difference approximations to scalar conservation laws

A unified treatment of explicit in time, two level, second order resolution, total variation diminishing, approximations to scalar conservation laws are presented. The schemes are assumed only to have conservation form and incremental form. A modified flux and a viscosity coefficient are introduced and results in terms of the latter are obtained. The existence of a cell entropy inequality is discussed and such an equality for all entropies is shown to imply that the scheme is an E scheme on monotone (actually more general) data, hence at most only first order accurate in general. Convergence for total variation diminishing-second order resolution schemes approximating convex or concave conservation laws is shown by enforcing a single discrete entropy inequality

Convenient stability criteria for difference approximations of hyperbolic initial-boundary value problems

New convenient stability criteria are provided in this paper for a large class of finite difference approximations to initial-boundary value problems associated with the hyperbolic system u sub t = au sub x + Bu + f in the quarter plane x or = 0, t or = 0. Using the new criteria, stability is easily established for numerous combinations of well known basic schemes and boundary conditions, thus generalizing many special cases studied in recent literature

Stability analysis of spectral methods for hyperbolic initial-boundary value systems

A constant coefficient hyperbolic system in one space variable, with zero initial data is discussed. Dissipative boundary conditions are imposed at the two points x = + or - 1. This problem is discretized by a spectral approximation in space. Sufficient conditions under which the spectral numerical solution is stable are demonstrated - moreover, these conditions have to be checked only for scalar equations. The stability theorems take the form of explicit bounds for the norm of the solution in terms of the boundary data. The dependence of these bounds on N, the number of points in the domain (or equivalently the degree of the polynomials involved), is investigated for a class of standard spectral methods, including Chebyshev and Legendre collocations

Convergence of spectral methods for hyperbolic initial-boundary value systems

A convergence proof for spectral approximations is presented for hyperbolic systems with initial and boundary conditions. The Chebyshev collocation is treated in detail, but the final result is readily applicable to other spectral methods, such as Legendre collocation or tau-methods