Another approach to the heat flow for Plateau problem

For the minimal surfaces in Rⁿ with Plateau boundary condition and establish the global existence and uniqueness of the flow as well as the continuous dependence of the initial datum.     

Approximation of a laminated microstructure for a rotationally invariant,double well energy density

Summary. We give error estimates for the approximation of a laminated microstructure which minimizes the energy for a rotationally invariant, double well energy density . We present error estimates for the convergence of the deformation in the convergence of directional derivatives of the deformation in the “twin planes,” the weak convergence of the deformation gradient, the convergence of the microstructure (or Young measure) of the deformation gradients, and the convergence of nonlinear integrals of the deformation gradient. Received July 25, 1995 / Revised version received November 20, 1995     

Discontinuous Galerkin method with a posterioriL_p(0,t_i) error estimate for second-kind Volterra problems

Summary. We consider a discontinuous Galerkin finite element method applied in time to a model Volterra equation of the second kind. A residual-based computable Galerkin-error estimate is derived for . This estimate does not explicitly contain the time step and therefore the time step control must be based on a heuristic criterion, the estimate can then be used to demonstrate the integrity, or otherwise, of the finite element solution. After performing some numerical experiments we conclude that this approach is at least competetive with classical discretizations since it is computationally simple to implement, but has the added advantage of reliable error feedback. Received June 25, 1995     

Lower bounds for the discrepancy of triples of inversive congruential pseudorandom numbers with power of two modulus

This paper deals with the inversive congruential method with power of two modulusm for generating uniform pseudorandom numbers. Statistical independence properties of the generated sequences are studied based on the distribution of triples of successive pseudorandom numbers. It is shown that there exist parameters in the inversive congruential method such that the discrepancy of the corresponding point sets in the unit cube is of an order of magnitude at leastm ^–1/3. The method of proof relies on a detailed analysis of certain rational exponential sums.     

Numerical integration of 2‐D integrals based on local bivariate C 1 quasi‐interpolating splines

In this paper cubature formulas based on bivariate C ¹ local polynomial splines with a four directional mesh [4] are generated and studied. Some numerical results with comparison with other methods are given. Moreover the method proposed is applied to the numerical evaluation of 2‐D singular integrals defined in the Hadamard finite part sense. Computational features, convergence properties and error bounds are proved. This revised version was published online in August 2006 with corrections to the Cover Date.     

Third order nonoscillatory central scheme for hyperbolic conservation laws

Summary. A third-order accurate Godunov-type scheme for the approximate solution of hyperbolic systems of conservation laws is presented. Its two main ingredients include: 1. A non-oscillatory piecewise-quadratic reconstruction of pointvalues from their given cell averages; and 2. A central differencing based on staggered evolution of the reconstructed cell averages. This results in a third-order central scheme, an extension along the lines of the second-order central scheme of Nessyahu and Tadmor \cite{NT}. The scalar scheme is non-oscillatory (and hence – convergent), in the sense that it does not increase the number of initial extrema (– as does the exact entropy solution operator). Extension to systems is carried out by componentwise application of the scalar framework. In particular, we have the advantage that, unlike upwind schemes, no (approximate) Riemann solvers, field-by-field characteristic decompositions, etc., are required. Numerical experiments confirm the high-resolution content of the proposed scheme. Thus, a considerable amount of simplicity and robustness is gained while retaining the expected third-order resolution. Received April 10, 1996 / Revised version received January 20, 1997     

Smoothing minimally supported frequency wavelets: Part II

The main purpose of this paper is to give a procedure to “mollify” the low-pass filters of a large number ofMinimally Supported Frequency (MSF) wavelets so that the smoother functions obtained in this way are also low-pass filters for an MRA. Hence, we are able to approximate (in the L ² -norm) MSF wavelets by wavelets with any desired degree of smoothness on the Fourier transform side. Although the MSF wavelets we consider are bandlimited, this may not be true for their smooth approximations. This phenomena is related to the invariant cycles under the transformation x ↦2x (mod2π). We also give a characterization of all low-pass filters for MSF wavelets. Throughout the paper new and interesting examples of wavelets are described.     

Higher order symplectic RK and RKN methods using perturbed collocation

In this paper, we derive a one-parameter family of symplectic Runge-Kutta-Nyström methods of order 2s-1 and a two-parameter family of symplectic Runge-Kutta methods of order 2s-2.