Limiting subhessians, limiting subjets and their calculus

We study calculus rules for limiting subjets of order two. These subjets are obtained as limits of sequences of subjets, a subjet of a function at some point being the Taylor expansion of a twice differentiable function which minorizes and coincides with at . These calculus rules are deduced from approximate (or fuzzy) calculus rules for subjets of order two. In turn, these rules are consequences of delicate results of Crandall-Ishii-Lions. We point out the similarities and the differences with the case of first order limiting subdifferentials.

     

On the Calculus of Limiting Subhessians

We discuss various qualification assumptions that allow calculus rules for limiting subhessians to be derived. Such qualification assumptions are based on a singular limiting subjet derived from a sequence of efficient subsets of symmetric matrices. We introduce a new efficiency notion that results in a weaker qualification assumption than that introduced in Ioffe and Penot (Trans Amer Math Soc 249: 789–807, 1997) and prove some calculus rules that are valid under this weaker qualification assumption. The work of A. Eberhard was supported by ARC research grant DP0664423.     

Mean field propagation of Wigner measures and BBGKY hierarchies for general bosonic states

Contrary to the finite dimensional case, Weyl and Wick quantizations are no more asymptotically equivalent in the infinite dimensional bosonic second quantization. Moreover neither the Weyl calculus defined for cylindrical symbols nor the Wick calculus defined for polynomials are preserved by the action of a nonlinear flow. Nevertheless taking advantage carefully of the information brought by these two calculuses in the mean field asymptotics, the propagation of Wigner measures for general states can be proved, extending to the infinite dimensional case a standard result of semiclassical analysis.     

Basic Cohomology for Singular Riemannian Foliations

 In this short note we prove that the basic cohomology of a singular Riemannian foliation of a compact manifold is a topological invariant and is finite dimensional. (Received 24 February 1998; in revised form 30 November 1998)     

On geometric vector fields of Minkowski spaces and their applications

As it is well-known, a Minkowski space is a finite dimensional real vector space equipped with a Minkowski functional F. By the help of its second order partial derivatives we can introduce a Riemannian metric on the vector space and the indicatrix hypersurface S:=F⁻¹(1) can be investigated as a Riemannian submanifold in the usual sense.Our aim is to study affine vector fields on the vector space which are, at the same time, affine with respect to the Funk metric associated with the indicatrix hypersurface. We give an upper bound for the dimension of their (real) Lie algebra and it is proved that equality holds if and only if the Minkowski space is Euclidean. Criteria of the existence is also given in lower dimensional cases. Note that in case of a Euclidean vector space the Funk metric reduces to the standard Cayley-Klein metric perturbed with a nonzero 1-form.As an application of our results we present the general solution of Matsumoto's problem on conformal equivalent Berwald and locally Minkowski manifolds. The reasoning is based on the theory of harmonic vector fields on the tangent spaces as Riemannian manifolds or, in an equivalent way, as Minkowski spaces. Our main result states that the conformal equivalence between two Berwald manifolds must be trivial unless the manifolds are Riemannian.     

Second order Lagrangian Twist systems: simple closed characteristics

We consider a special class of Lagrangians that play a fundamental role in the theory of second order Lagrangian systems: Twist systems. This subclass of Lagrangian systems is defined via a convenient monotonicity property that such systems share. This monotonicity property (Twist property) allows a finite dimensional reduction of the variational principle for finding closed characteristics in fixed energy levels. This reduction has some similarities with the method of broken geodesics for the geodesic variational problem on Riemannian manifolds. On the other hand, the monotonicity property can be related to the existence of local Twist maps in the associated Hamiltonian flow.

The finite dimensional reduction gives rise to a second order monotone recurrence relation. We study these recurrence relations to find simple closed characteristics for the Lagrangian system. More complicated closed characteristics will be dealt with in future work. Furthermore, we give conditions on the Lagrangian that guarantee the Twist property.

     

Nodal domains and growth of harmonic functions on noncompact manifolds

Harmonic functions are studied on complete Riemannian manifolds. A decay estimate is given for bounded harmonic functions of variable sign. For unbounded harmonic functions of variable sign, relations are derived between growth properties and nodal domains. On Riemannian manifolds of nonnegative Ricci curvature, it has been conjectured that harmonic functions, having at most a given order of polynomial growth, must form a finite dimensional vector space. This conjecture is established in certain special cases.     

The parameterization method in optimal control problems and differential–algebraic equations

The paper is devoted to the explanation of the numerical parameterization method (PM) for optimal control (OC) problems with intermediate phase constraint and to its circumstantiation for classical calculus of variation (CV) problems that arise in connection with singular ODEs or DAEs, especially in cases of their essential degeneracy. The PM is based on a finite parameterization of control functions and on derivation of the problem functional with respect to control parameters. The first and the second derivatives are calculated with the help of adjoint vector and matrix impulses. Results of the solution to one phase constrained OC and two degenerate CV problems, connected with singular DAEs nonreducible to the normal form, are presented.     

Two Dimensional Interface Problems for Elliptic Equations

We study the structure of solutions to the interface problems for second order quasi-linear elliptic partial differential equations in two dimensional space. We prove that each weak solution can be decomposed into two parts near singular points, a finite sum of functions in the form of cr^α log^m rφ(θ) and a regular one w. The coefficients c and the C^{1,α} norm of w depend on the H¹-norm and the C^{º,α}-norm of the solution, and the equation only.     

On the numerical solution of diffusion–reaction equations with singular source terms

A numerical study is presented of reaction–diffusion problems having singular reaction source terms, singular in the sense that within the spatial domain the source is defined by a Dirac delta function expression on a lower dimensional surface. A consequence is that solutions will be continuous, but not continuously differentiable. This lack of smoothness and the lower dimensional surface form an obstacle for numerical discretization, including amongst others order reduction. In this paper the standard finite volume approach is studied for which reduction from order two to order one occurs. A local grid refinement technique is discussed which overcomes the reduction.     

格点上的非交换微分运算及其应用

By introducing the noncommutative differential calculus on the function space of the infinite/finite set and construct a homotopy operator, one prove the analogue of the Poincare lemma for the difference complex. As an application of the differential calculus, a two dimensional integral model can be derived from the noncommutative differential calculus.     

Boundary integral method for thermoelastic screen scattering problem in ℝ3

We investigate a three‐dimensional mathematical thermoelastic scattering problem from an open surface which will be referred to as a screen. Under the assumption of the local finite energy of the unified thermoelastic scattered field, we give a weak model on the appropriate Sobolev spaces and derive equivalent integral equations of the first kind for the jump of some trace operators on the open surface. Uniqueness and existence theorems are proved, the regularity and the singular behaviour of the solution near the edge are established with the help of the Wiener–Hopf method in the halfspace, the calculus of pseudodifferential operators on the basis of the strong ellipticity property and Gårding's inequality. An improved Galerkin scheme is provided by simulating the singular behaviour of the exact solution at the edge of the screen. Copyright © 2000 John Wiley & Sons, Ltd.     

You can hear the local orientability of an orbifold

A Riemannian orbifold is a mildly singular generalization of a Riemannian manifold which is locally modeled on the quotient of a connected, open manifold under a finite group of isometries. If all of the isometries used to define the local structures of an entire orbifold are orientation preserving, we call the orbifold locally orientable. We use heat invariants to show that a Riemannian orbifold which is locally orientable cannot be Laplace isospectral to a Riemannian orbifold which is not locally orientable. As a corollary we observe that a Riemannian orbifold that is not locally orientable cannot be Laplace isospectral to a Riemannian manifold.     

The upwind finite difference fractional steps method for nonlinear coupled systems

For nonlinear coupled system of multilayer dynamics of fluids in porous media, a second‐order upwind finite‐difference fractional‐steps scheme applicable to parallel arithmetic are put forward, and two‐ and three‐dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, multiplicative commutation rule of difference operators, decomposition of high‐order difference operators, and prior estimates are adopted. Optimal order estimates in l² norm are derived to determine the error in the approximate solution. This method has already been applied to the numerical simulation of migration‐accumulation of oil resources. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Partition properties of the dense local order and a colored version of Milliken’s theorem

We study finite dimensional partition properties of the countable homogeneous dense local order (a directed graph closely related to the order structure of the rationals). Some of our results use ideas borrowed from the partition calculus of the rationals and are obtained thanks to a strengthening of Milliken’s theorem on trees.     

Variational splines for the numerical solution of singular differential equations

The numerical parametrization method (PM), originally created for optimal control problems, is specificated for classical calculus of variation problems that arise in connection with singular implicit (IDEs) and differential-algebraic equations (DAEs). The PM for IDEs is based on representation of the required solution as a spline with moving knots and on minimization of the discrepancy functional with respect to the spline parameters. Such splines are named variational splines. For DAEs only finite entering functions can be represented by splines, and the functional under minimization is the discrepancy of the algebraic subsystem. The first and the second derivatives of the functionals are calculated in two ways – for DAEs with the help of adjoint variables, and for IDE directly. The PM does not use the notion of differentiation index, and it is applicable to any singular equation having a solution. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

SOME EXTENSIONS OF THE MEAN CURVATURE FLOW IN RIEMANNIAN MANIFOLDS

Given a family of smooth immersions of closed hypersurfaces in a locally symmetric Riemannian manifold with bounded geometry,moving by mean curvature flow,we show that at the first finite singular time of mean curvature flow,certain subcritical quantities concerning the second fundamental form blow up.This result not only generalizes a result of Le-Sesum and Xu-Ye-Zhao,but also extends the latest work of Le in the Euclidean case.     

Simplified immersed interface methods for elliptic interface problems with straight interfaces

In this article, we propose simplified immersed interface methods for elliptic partial/ordinary differential equations with discontinuous coefficients across interfaces that are few isolated points in 1D, and straight lines in 2D. For one‐dimensional problems or two‐dimensional problems with circular interfaces, we propose a conservative second‐order finite difference scheme whose coefficient matrix is symmetric and definite. For two‐dimensional problems with straight interfaces, we first propose a conservative first‐order finite difference scheme, then use the Richardson extrapolation technique to get a second‐order method. In both cases, the finite difference coefficients are almost the same as those for regular problems. Error analysis is given along with numerical example. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 188–203, 2012     

A parametrization method for the numerical solution of singular differential equations

The paper explains the numerical parametrization method (PM), originally created for optimal control problems, for classical calculus of variation problems that arise in connection with singular implicit (IDEs) and differential-algebraic equations (DAEs) in frame of their regularization. The PM for IDEs is based on representation of the required solution as a spline with moving knots and on minimization of the discrepancy functional with respect to the spline parameters. Such splines are named variational splines. For DAEs only finite entering functions can be represented by splines, and the functional under minimization is the discrepancy of the algebraic subsystem. The first and the second derivatives of the functionals are calculated in two ways – for DAEs with the help of adjoint variables, and for IDE directly. The PM does not use the notion of differentiation index, and it is applicable to any singular equation having a solution.     

Line energies for gradient vector fields in the plane

In this paper we study the singular perturbation of by . This problem, which could be thought as the natural second order version of the classical singular perturbation of the potential energy by , leads, as in the first order case, to energy concentration effects on hypersurfaces. In the two dimensional case we study the natural domain for the limiting energy and prove a compactness theorem in this class. Received January 19, 1999 / Accepted February 26, 1999