Maximum norm stability of difference schemes for parabolic equations on overset nonmatching space-time grids

In this paper, theoretical results are described on the maximum norm stability and accuracy of finite difference discretizations of parabolic equations on overset nonmatching space-time grids. We consider parabolic equations containing a linear reaction term on a space-time domain which is decomposed into an overlapping collection of cylindrical subregions of the form , for . Each of the space-time domains are assumed to be independently grided (in parallel) according to the local geometry and space-time regularity of the solution, yielding space-time grids with mesh parameters and . In particular, the different space-time grids need not match on the regions of overlap, and the time steps can differ from one grid to the next. We discretize the parabolic equation on each local grid by employing an explicit or implicit -scheme in time and a finite difference scheme in space satisfying a discrete maximum principle. The local discretizations are coupled together, without the use of Lagrange multipliers, by requiring the boundary values on each space-time grid to match a suitable interpolation of the solution on adjacent grids. The resulting global discretization yields a large system of coupled equations which can be solved by a parallel Schwarz iterative procedure requiring some communication between adjacent subregions. Our analysis employs a contraction mapping argument.

Applications of the results are briefly indicated for reaction-diffusion equations with contractive terms and heterogeneous hyperbolic-parabolic approximations of parabolic equations.

     

Small Solutions of Quadratic Diophantine Equations

We consider quadratic diophantine equations of the shape for a polynomial Q(X₁, ..., X_s) Z[X₁, ..., X_s] of degree 2.Let H be an upper bound for the absolute values of the coefficientsof Q, and assume that the homogeneous quadratic part of Q isnon-singular. We prove, for all s 3, the existence of a polynomialbound _s(H) with the following property: if equation (1) hasa solution x Z^s at all, then it has one satisfying For s = 3 and s = 4 no polynomial bounds _s(H) were previouslyknown, and for s 5 we have been able to improve existing boundsquite significantly. 2000 Mathematics Subject Classification11D09, 11E20, 11H06, 11P55.     

On the Differential Properties of Temlyakov-Type and Temlyakov-Bavrin-Type Integrals

We establish a formula that expresses the value of the density of the Temlyakov-type integral with defining domain of type A in terms of the integral itself by means of a differential operator. We extend this operator connection to the Temlyakov–Bavrin-type integral of order k. Our study bases on the method of linear differential operators with variable coefficients which is developed by A. V. Nelaev.     

Dirac Operators on Hypersurfaces of Manifolds with Negative Scalar Curvature

We give a sharp extrinsic lower bound for the first eigenvaluesof the intrinsic Dirac operator of certain hypersurfaces boundinga compact domain in a spin manifold of negative scalar curvature.Limiting-cases are characterized by the existence, on the domain,of imaginary Killing spinors. Some geometrical applications, as anAlexandrov type theorem, are given.     

On Solvability of a Boundary Value Problem for Singular Integral Equations

This paper deals with the solvability of boundary value problems for singular integral equations of the form (i)-(ii).By an algebraic method we reduce the problem (i)-(ii) to a system of linear algebraic equations which gives all solutions in a closed form.AMS Subject Classification: 47G05, 45GO5, 45E05     

Consistent Initialization for Nonlinear Index-2 Differential–Algebraic Equation: Large Sparse Systems in MATLAB

An important component of any initial-value solver for higher-index differential–algebraic equations consists in the computation of consistent initial values. In a recent paper [5], an algorithm is proposed which is applicable to a very general class of index-2 systems. Unfortunately, the computational expense is rather high. We present a modification of this approach, which gives rise to a MATLAB implementation capable of handling systems of moderate dimension (several thousands of unknowns). The algorithm is illustrated by examples.     

The Cauchy Problem for the System of Maxwell Equations

We consider the problem of analytic continuation of a solution to the system of Maxwell equations in a bounded spatial domain from data on part of the boundary of the domain. We construct an approximate solution to the problem using the Carleman matrix method.     

IRK Methods for Index 2 and 3 DAEs: Starting Algorithms

When semi-explicit differential-algebraic equations are solved with implicit Runge-Kutta methods, the computational effort is dominated by the cost of solving the non-linear systems. That is why it is important to have good starting values to begin the iterations. In this paper we study a type of starting algorithms, without additional computational cost, in the case of index-2 and index-3 DAEs. The order of the starting values is defined, and by using DA-series and rooted trees we obtain their general order conditions. If the RK method satisfies some simplifying assumptions, then the maximum order can be obtained.     

On a Problem of Karpilovsky

Let G be a finite elementary group. Let ⁿ (G) denote the nth power of the augmentation ideal (G) of the integral group ring G. In this paper, we give an explicit basis of the quotient group Q_n(G) = ⁿ(G)/ⁿ⁺¹ (G) and compute the order of Qⁿ (G).2000 Mathematics Subject Classification: 16S34, 20C05     

振荡型积分对极大函数的估计

首先将2维KaKeya型的极大函数通过Fourier变换转化成振荡型积分的形式,对此振荡型积分作相应的估计。通过Littlewood—Paley方法将此振荡型积分分解成不同的算子,在估计极大算子时,通过Newton—Leibnig公式,将其转化为相应的积分形式,最后通过经典的振荡型积分的估计,得到最后的结果,本方法强调振荡型积分的方法和以前几何测度论的方法不同。