Stability properties of the ENO method

We review the currently available stability properties of the ENO reconstruction procedure, such as its monotonicity and non-oscillatory properties, the sign property, upper bounds on cell interface jumps and a total variation-type bound. We also outline how these properties can be applied to derive stability and convergence of high-order accurate schemes for conservation laws.Comment: To appear in Handbook of Numerical Methods for Hyperbolic Problem

Rhinoceros Horn Libation Cup

On display in the “Wonders of Nature and Artifice” exhibit at Gettysburg College is an exquisitely carved Chinese rhinoceros horn cup decorated with many images of animals, from dragons to tortoises.The rhinoceros horn has been noted by the Chinese as early as the T’ang dynasty (618-907) to have magical properties, and it was believed that when a poisonous liquid was poured into a rhino horn, the horn would change colors to alert to the presence of poison.Due to these magical properties, rhinoceros horns have been regarded as especially valuable. [excerpt

Recent developments in shock-capturing schemes

The development of the shock capturing methodology is reviewed, paying special attention to the increasing nonlinearity in its design and its relation to interpolation. It is well-known that higher-order approximations to a discontinuous function generate spurious oscillations near the discontinuity (Gibbs phenomenon). Unlike standard finite-difference methods which use a fixed stencil, modern shock capturing schemes use an adaptive stencil which is selected according to the local smoothness of the solution. Near discontinuities this technique automatically switches to one-sided approximations, thus avoiding the use of discontinuous data which brings about spurious oscillations

Multi-resolution analysis for ENO schemes

Given an function, u(x), which is represented by its cell-averages in cells which are formed by some unstructured grid, we show how to decompose the function into various scales of variation. This is done by considering a set of nested grids in which the given grid is the finest, and identifying in each locality the coarsest grid in the set from which u(x) can be recovered to a prescribed accuracy. This multi-resolution analysis was applied to essentially non-oscillatory (ENO) schemes in order to advance the solution by one time-step. This is accomplished by decomposing the numerical solution at the beginning of each time-step into levels of resolution, and performing the computation in each locality at the appropriate coarser grid. An efficient algorithm for implementing this program in the 1-D case is presented; this algorithm can be extended to the multi-dimensional case with Cartesian grids

Uniformly high-order accurate non-oscillatory schemes, 1

The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws was begun. These schemes share many desirable properties with total variation diminishing schemes (TVD), but TVD schemes have at most first order accuracy, in the sense of truncation error, at extreme of the solution. A uniformly second order approximation was constucted, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell

Some physical implications of recent solar wind measurements

The physical implications of the existence at about 1 AU of a quiet solar wind particle flux about 90 percent larger than that suggested in the past is investigated within the framework of the two-fluid solar wind model equations. During the spherically symmetric radial expansion of the quiet solar wind, the particle flux is conserved quantity. It is found that a pure collisional two-fluid model provides good particle density and streaming velocity at 1 AU, but predicts too large an electron temperature and too small a proton temperature. When noncollisional contributions to the transport coefficients are incorporated in the model equations, a complete satisfactory agreement with the available observations is obtained. Upper limits to the effective coupling between electrons and protons, as well as to the effective proton thermal conductivity, and both upper and lower limits to the effective electron thermal conductivity in the quiet solar wind, required to provide agreement with observations, are given

Optimal clustering of frequency-constrained maintenance jobs with shared set-ups

Since maintenance jobs often require one or more set-up activities, joint execution or clustering of maintenance jobs is a powerful instrument to reduce shut-down costs. We consider a clustering problem for frequency-constrained maintenance jobs, i.e. maintenance jobs that must be carried out with a prescribed (or higher) frequency. For the clustering of maintenance jobs with identical, so-called common set-ups, several strong dominance rules are provided. These dominance rules are used in an efficient dynamic programming algorithm which solves the problem in polynomial time. For the clustering of maintenance jobs with partially identical, so-called shared set-ups, similar but less strong dominance rules are available. Nevertheless, a surprisingly well-performing greedy heuristic and a branch and bound procedure have been developed to solve this problem. For randomly generated test problems with 10 set-ups and 30 maintenance jobs, the heuristic was optimal in 47 out of 100 test problems, with an average deviation of 0.24% from the optimal solution. In addition, the branch and bound method found an optimal solution in only a few seconds computation time on average

A study of the XY model by the Monte Carlo method

The massively parallel processor is used to perform Monte Carlo simulations for the two dimensional XY model on lattices of sizes up to 128 x 128. A parallel random number generator was constructed, finite size effects were studied, and run times were compared with those on a CRAY X-MP supercomputer

Multi-Dimensional ENO Schemes for General Geometries

A class of ENO schemes is presented for the numerical solution of multidimensional hyperbolic systems of conservation laws in structured and unstructured grids. This is a class of shock-capturing schemes which are designed to compute cell-averages to high order accuracy. The ENO scheme is composed of a piecewise-polynomial reconstruction of the solution form its given cell-averages, approximate evolution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is based on an adaptive selection of stencil for each cell so as to avoid spurious oscillations near discontinuities while achieving high order of accuracy away from them