A unified approach to a posteriori error estimation using element residual methods

Summary This paper deals with the problem of obtaining numerical estimates of the accuracy of approximations to solutions of elliptic partial differential equations. It is shown that, by solving appropriate local residual type problems, one can obtain upper bounds on the error in the energy norm. Moreover, in the special case of adaptiveh-p finite element analysis, the estimator will also give a realistic estimate of the error. A key feature of this is the development of a systematic approach to the determination of boundary conditions for the local problems. The work extends and combines several existing methods to the case of fullh-p finite element approximation on possibly irregular meshes with, elements of non-uniform degree. As a special case, the analysis proves a conjecture made by Bank and Weiser [Some A Posteriori Error Estimators for Elliptic Partial Differential Equations, Math. Comput.44, 283–301 (1985)].     

Regularity of Solutions to Bilinear Stochastic Wave Equations

Abstract

The paper is concerned with strong solutions of bilinear stochastic wave equations in ? ^d , of which the coefficients contain semimartingale white noises with spatial parameters. For the Cauchy problems, the existence and spatial regularity of solutions in Sobolev spaces are proved under appropriate conditions. The dependence of solution regularity on the smoothness of the random coefficients is ascertained. The proofs are based on stochastic energy inequalities, the semigroup method and certain submartingale inequalities. Regularity results are also obtained for the special case of Wiener semimartingales.     

Second Order Scalar Functional Differential Equations with Piecewise Constant Arguments

The study of functional differential equations with piecewise constant arguments usually results in a study of certain related difference equations. In this paper we consider certain neutral functional differential equations of this type and the associated difference equations. We give conditions under which such equations with almost periodic time dependence will have unique almost periodic solutions, and for certain autonomous cases, we obtain certain stability results and also conditions for chaotic behavior of solutions. We are particularly concerned with such equations which are partially discretized versions of non-forced Duffing equations.     

*-RINGS AND THEIR RADICALS

ABSTRACT

This paper is to investigate the class * of all prima rings without proper non-zero prima factors and to characterize the ‘minimal special radical’ containing the class * as well as the ‘maximal special radical’ the semisimple class of which contains the class *. This let us answer certain open questions put in [1] and [4].     

Cohomology of the moduli of parabolic vector bundles

The purpose of this paper is to compute the Betti numbers of the moduli space ofparabolic vector bundles on a curve (see Seshadri [7], [8] and Mehta & Seshadri [4]), in the case where every semi-stable parabolic bundle is necessarily stable. We do this by generalizing the method of Atiyah and Bott [1] in the case of moduli of ordinary vector bundles. Recall that (see Seshadri [7]) the underlying topological space of the moduli of parabolic vector bundles is the space of equivalence classes of certain unitary representations of a discrete subgroup Γ which is a lattice in PSL (2,R). (The lattice Γ need not necessarily be co-compact). While the structure of the proof is essentially the same as that of Atiyah and Bott, there are some difficulties of a technical nature in the parabolic case. For instance the Harder-Narasimhan stratification has to be further refined in order to get the connected strata. These connected strata turn out to have different codimensions even when they are part of the same Harder-Narasimhan strata. If in addition to ‘stable = semistable’ the rank and degree are coprime, then the moduli space turns out to be torsion-free in its cohomology. The arrangement of the paper is as follows. In § 1 we prove the necessary basic results about algebraic families of parabolic bundles. These are generalizations of the corresponding results proved by Shatz [9]. Following this, in § 2 we generalize the analytical part of the argument of Atiyah and Bott (§ 14 of [1]). Finally in § 3 we show how to obtain an inductive formula for the Betti numbers of the moduli space. We illustrate our method by computing explicitly the Betti numbers in the special case of rank = 2, and one parabolic point.     

On a characterization of monotone likelihood ratio experiments

Summary Pfanzagl (1962,Zeit. Wahrscheinlichkeitsth.,1, 109–115) showed that a dominated family of probability measures has monotone likelihood ratios with respect to some real valued statistic if there exists a set of tests which has certain nice properties. A similar characterization was given by Dettweiler (1978,Metrika,25, 247–254), who did not assume domination. However, Pfanzagl's result is not a special case of the one proved by Dettweiler. We present a theorem which comprises the results of both authors. Our proof shows that not all conditions introduced by them are needed. Furthermore, we investigate the question concerning the generality we get if we do not assume domination.     

On Ω-limit Sets of Discrete-time Dynamical Systems

The paper investigates z -limit sets for discrete-time dynamical systems of the form x n +1 = f n +1 ( x n ), n S 0, with each f n mapping an interval I of R into itself. For autonomous systems, i.e. f n = f for all n , and f continuous on I =[ a , b ], the case that all z -limit sets consist of one point only is characterized by several equivalent conditions, one being that f has no 2-periodic points. The non-autonomous case assumes that the functions f n converge uniformly to a continuous function f X that has no 2-periodic points. It is shown that the z -limit sets are closed intervals consisting of fixed points of f X only. Under certain conditions these closed intervals contain exactly one point each. This allows a treatment of certain discrete-time dynamical systems in R n .     

ADJOINT PROBLEMS IN THE CALCULUS OF VARIATIONS OF MULTIPLE INTEGRALS

Abstract

A general theory of adjoint variational problems is formulated for essentially arbitrary Lagrangians involving m independent and n dependent variables, together with the first derivatives of the latter, This approach contains as a special case the theory of Haar [4], in which the Lagrangian may depend solely on the derivatives of a single dependent function of two arguments. Because of the eventual occurrence of possibly incompatible sets of integrability conditions, the basic theory is developed against the background of non-integrable m-dimensional subspaces, which is in sharp contrast to the traditional approach to the calculus of variations. Relatively self-adjoint Lagrangians are defined and completely characterized in terms of an arbitrary Riemannian metric. In the course of the general theory certain geometric object fields are encountered in a very natural manner, some of which had arisen previously in the canonical formalism proposed by Caratheodory [2]. Accordingly the analysis of the present paper may serve to shed some light on this conceptually extremely difficult formalism.     

Constructing convex solutions via Perron’s method

Abstract In this paper we construct convex solutions for certain elliptic boundary value problems via Perron’s method. The solutions constructed are weak solutions in the viscosity sense, and our construction follows work of Ishii (Duke Math. J., 55 (2) 369–384, 1987). The same general approach appears in work of Andrews and Feldman (J. Differential Equations, 182 (2) 298–343, 2002) in which they show existence for a weak nonlocal parabolic flow of convex curves. The time independent special case of their work leads to a one dimensional elliptic result which we extend to two dimensions. Similar results are required to extend their theory of nonlocal geometric flows to surfaces. The two dimensional case is essentially different from the one dimensional case and involves a regularity result (cf. Theorem 3.1), which has independent interest. Roughly speaking, given an arbitrary convex function (which is not smooth) supported at one point by a smooth function of prescribed Hessian (which is not convex), one must construct a third function that is both convex and smooth and appropriately approximates both of the given functions. Keywords: Viscosity solutions, Elliptic partial differential equations, Perron procedure, Convexity, Regularity, Fully nonlinear, Monge-Ampere Mathematics Subject Classification (2000:) 35J60, 53A05, 52A15, 26B05     

Calibrated Subbundles in Noncompact Manifolds of Special Holonomy

This paper is a continuation of Math. Res. Lett. 12 (2005), 493–512. We first construct special Lagrangian submanifolds of the Ricci-flat Stenzel metric (of holonomy SU(n)) on the cotangent bundle of Sⁿ by looking at the conormal bundle of appropriate submanifolds of Sⁿ. We find that the condition for the conormal bundle to be special Lagrangian is the same as that discovered by Harvey–Lawson for submanifolds in Rⁿ in their pioneering paper, Acta Math. 148 (1982), 47–157. We also construct calibrated submanifolds in complete metrics with special holonomy G₂ and Spin(7) discovered by Bryant and Salamon (Duke Math. J. 58 (1989), 829–850) on the total spaces of appropriate bundles over self-dual Einstein four manifolds. The submanifolds are constructed as certain subbundles over immersed surfaces. We show that this construction requires the surface to be minimal in the associative and Cayley cases, and to be (properly oriented) real isotropic in the coassociative case. We also make some remarks about using these constructions as a possible local model for the intersection of compact calibrated submanifolds in a compact manifold with special holonomy. Mathematics Subject Classification (2000): 53-XX, 58-XX.     

Coincidences of Nonself Operators by a Simpler Algorithm

Abstract

In this article, we introduce a new class of contractive mappings and study analytical and computational aspects of a special case of Jungck-Khan iterative algorithm generated by this class of mappings. In particular, we improve upon strong convergence, rate of convergence and data dependence results existing in the current literature. Analytical as well numerical illustrative examples are given to support the new results.     

ON THE MINIMUM INTRINSIC 1-VOLUME OF VORONOI CELLS IN LATTICE UNIT SPHERE PACKINGS

A typical 3-dimensional (in short '3D') Voronoi cell of a 3Dlattice has six families of parallel edges. We call any six representants of these six families the generating edges of the Voronoi cell. The sum s of lengths of generating edges of a Voronoi-cell of a lattice unit sphere packing in the 3-dimensional Euclidean space is a special case of intrinsic 1-volumes of 3Dzonotopes with inradius 1 which are investigated accurately in [B]. However, the minimum of this value is unknown even in this special case. As the regular rhombic dodecahedron shows optimal properties in many similar problems, it was reasonable to conjecture that it also has the minimal s value. In this note we present a construction of a lattice unit ball packing whose Voronoi cell possesses an intrinsic 1-volume strictly less than the one of the proper regular rhombic dodecahedron, hence providing a smaller upper bound for s than it was conjectured. A further issue of the note is a formula for edge-lengths of Voronoi cells of lattice unit ball packings that can be used efficiently in similar calculations. This revised version was published online in August 2006 with corrections to the Cover Date.     

The Kleisli Category is a Category Of Fractions

We consider a generalised notion of category of fractions associated with a class Φ of permitted replacement situations between zig-zags in a category C. The Kleisli category associated with a monad [T, η, μ] in C is shown to be a special case in which the arrows {ηX X ? C} become left-invertible after passage to fractions.     

Energy bounds for entire nodal solutions of autonomous superlinear equations

We prove new lower bounds for the energy of sign changing solutions of two different types of autonomous superlinear equations on . These bounds increase the energy range for which Palais-Smale sequences of related variational problems have a nontrivial weak limit. The results are obtained by transforming sign changing solutions of low energy into a convenient form and then applying different variants of the moving plane method. Mathematics Subject Classification (2000) 35J60, 35J20     

Visibility graphs of staircase polygons and the weak Bruhat order,I: From visibility graphs to maximal chains

The recognition problem for visibility graphs of simple polygons is not known to be in NP, nor is it known to be NP-hard. It is, however, known to be inPSPACE. Further, every such visibility graph can be dismantled as a sequence of visibility graphs of convex fans. Any nondegenerated configuration ofn points can be associated with amaximal chain in the weak Bruhat order of the symmetric groupS _n . The visibility graph ofany simple polygon defined on this configuration is completely determined by this maximal chain via a one-to-one correspondence between maximal chains andbalanced tableaux of a certain shape. In the case of staircase polygons (special convex fans), we define a class of graphs calledpersistent graphs and show that the visibility graph of a staircase polygon is persistent. We then describe a polynomial-time algorithm that recovers a representative maximal chain in the weak Bruhat order from a given persistent graph, thus characterizing the class of persistent graphs. The question of recovering a staircase polygon from a given persistent graph, via a maximal chain, is studied in the companion paper [4]. The overall goal of both papers is to offer a characterization of visibility graphs, of convex fans. The research of J. Abello was supported by NSF Grants Nos. DCR 8603722 and DCR 8896281. This research was done while K. Kumar was at the Department of Computer Science, Texas A & M University.     

RESTRICTING THE COMPARISON FUNCTOR OF AN ADJUNCTION TO PROJECTIVE OBJECTS

Abstract

This paper deals with projectives (in the sense of K.A.Hardie [5] relative to a right adjoint functor U: A → K. We answer the question, raised by R.-E. Hoffmann [6] p. 135, of knowing under what conditions there exists an equivalence between Proj u and Proj U^r, induced by the comparison functor Φ: A → K^T, where T denotes the monad induced by U. In the case, that U is an algebraic functor we also give necessary and sufficient conditions for the re gular projective objects to coincide with the U-projectives. Finally, we delineate how these results nay be applied in certain familiar situations.     

Combined time and energy optimal trajectory planning with quadratic drag for mixed discrete-continuous task planning

ABSTRACT

The problem of mixed discrete-continuous task planning for mechanical systems, such as aerial drones or other autonomous units, can be often treated as a sequence of point-to-point trajectories. In this work, the problem of optimal trajectory planning under a combined completion time and energy criterion, for a straight point to point path for a second-order system with quadratic under state (velocity) and control (acceleration) constraints is considered. The solution is obtained and proved to be optimal using the Pontryagin Maximum Principle. Simulation results for different cases are presented and compared with a customary numerical optimal control solver.     

Blowup for the compressible isothermal Euler equations with non-vacuum initial data

ABSTRACT

In this paper, we study the compressible isothermal Euler equations with non-vacuum initial data. First, we prove the property of finite propagation to this Cauchy problem by using local energy estimates. Second, we establish the blowup results of the multi-dimensional case in radial symmetry and the one-dimensional case in non-radial symmetry by making assumptions on the initial velocity. Third, we present the blowup results of the three-dimensional case in non-radial symmetry by making assumptions on the initial momentum.     

The rodrigues type representations for a certain class of special functions

Summary An attempt is made here to present a systematic introduction to and several applications of a certain method of obtaining Rodrigues type representations for a fairly wide variety of sequences of special functions. The main results, contained in Theorems2 and3 below, are shown to apply not only to the Bessel polynomials, the classical orthogonal polynomials including, for instance, Hermite, Jacobi (and, of course, Gegenbauer, Legendre, and Tchebycheff), and Laguerre polynomials, and to their various generalizations studied in recent years, but also to such other special functions as the Bessel function and a certain class of generalized hypergeometric functions. Entrata in Redazione il 25 giugno 1977. This work was partially supported by the National Research Council of Canada under grants A-7353 and A-4027. For a preliminary report of this paper see Notices Amer. Math. Soc.,24 (1977), p. A-238, Abstract no. 77T-B43.