ON LOCAL STRUCTURAL STABILITY OF ONE-DIMENSIONAL SHOCKS IN RADIATION HYDRODYNAMICS

In this paper, we are concerned with the local structural stability of one-dimensional shock waves in radiation hydrodynamics described by the isentropic Euler-Boltzmann equations. Even though in this radiation hydrodynamics model, the radiative effects can be understood as source terms to the isentropic Euler equations of hydrodynamics, in general the radiation field has singularities propagated in an angular domain issuing from the initial point across which the density is discontinuous. This is the major difficulty in the stability analysis of shocks. Under certain assumptions on the radiation parameters, we show there exists a local weak solution to the initial value problem of the one dimensional Euler-Boltzmann equations, in which the radiation intensity is continuous, while the density and velocity are piecewise Lipschitz continuous with a strong discontinuity representing the shock-front. The existence of such a solution indicates that shock waves are structurally stable, at least local in time, in radiation hydrodynamics.

The Clustering Coefficient and the Diameter of Small-world Networks

The small-world network, proposed by Watts and Strogatz, has been extensively studied for the past over ten years. In this paper, a generalized smMl-world network is proposed, which extends severM small-world network models. Furthermore, some properties of a special type of generalized small-world network with given expectation of edge numbers have been investigated, such as the degree distribution and the isoperimetric number. These results are used to present a lower and an upper bounds for the clustering coefficient and the diameter of the given edge number expectation generalized small-world network, respectively. In other words, we prove mathematically that the given edge number expectation generalized small-world network possesses large clustering coefficient and small diameter.

An adaptive time stepping method with efcient error control for second-order evolution problems

This work is concerned with time stepping finite element methods for abstract second order evolution problems. We derive optimal order a posteriori error estimates and a posteriori nodal superconvergence error estimates using the energy approach and the duality argument. With the help of the a posteriori error estimator developed in this work, we will further propose an adaptive time stepping strategy. A number of numerical experiments are performed to illustrate the reliability and efficiency of the a posteriori error estimates and to assess the effectiveness of the proposed adaptive time stepping method.

Partially positive matrices

A real n × n symmetric matrix P is partially positive(PP) for a given index set I ? {1,..., n} if there exists a matrix V such that V(I, :) 0 and P = V VT. We give a characterization of PP-matrices. A semidefinite algorithm is presented for checking whether a matrix is partially positive or not. Its properties are studied. A PP-decomposition of a matrix can also be obtained if it is partially positive.