Runge–Kutta convolution quadrature methods for well-posed equations with memory

Runge–Kutta based convolution quadrature methods for abstract, well-posed, linear, and homogeneous Volterra equations, non necessarily of sectorial type, are developed. A general representation of the numerical solution in terms of the continuous one is given. The error and stability analysis is based on this representation, which, for the particular case of the backward Euler method, also shows that the numerical solution inherits some interesting qualitative properties, such as positivity, of the exact solution. Numerical illustrations are provided.     

A fictitious domain approach to the numerical solution of PDEs in stochastic domains

We present an efficient method for the numerical realization of elliptic PDEs in domains depending on random variables. Domains are bounded, and have finite fluctuations. The key feature is the combination of a fictitious domain approach and a polynomial chaos expansion. The PDE is solved in a larger, fixed domain (the fictitious domain), with the original boundary condition enforced via a Lagrange multiplier acting on a random manifold inside the new domain. A (generalized) Wiener expansion is invoked to convert such a stochastic problem into a deterministic one, depending on an extra set of real variables (the stochastic variables). Discretization is accomplished by standard mixed finite elements in the physical variables and a Galerkin projection method with numerical integration (which coincides with a collocation scheme) in the stochastic variables. A stability and convergence analysis of the method, as well as numerical results, are provided. The convergence is “spectral” in the polynomial chaos order, in any subdomain which does not contain the random boundaries.     

Lagrangian skeletons,periodic geodesic flows and symplectic cuttings

In this paper, we investigate symplectic manifolds endowed with a Morse–Bott function with only two critical submanifolds, one of which is Lagrangian while the other one is symplectic.     

Asymptotic Distribution of Nodes for Near-Optimal Polynomial Interpolation on Certain Curves in R 2

Let E\subset \Bbb R ^s be compact and let d _n ^E denote the dimension of the space of polynomials of degree at most n in s variables restricted to E . We introduce the notion of an asymptotic interpolation measure (AIM). Such a measure, if it exists , describes the asymptotic behavior of any scheme τ _n ={ \bf x _k,n } _k=1 ^dnE , n=1,2,\ldots , of nodes for multivariate polynomial interpolation for which the norms of the corresponding interpolation operators do not grow geometrically large with n . We demonstrate the existence of AIMs for the finite union of compact subsets of certain algebraic curves in R ² . It turns out that the theory of logarithmic potentials with external fields plays a useful role in the investigation. Furthermore, for the sets mentioned above, we give a computationally simple construction for ``good' interpolation schemes. November 9, 2000. Date revised: August 4, 2001. Date accepted: September 14, 2001.     

超软X射线流气式正比计数管

介绍了一种用于测量183～933eV超软X射线的圆柱形、侧窗式、流气式正比计数管,工作气体是0.11MPa的P-10气体或氦气-丙烷混合气体。计数管内径为φ25mm,直径为φ0.3mm的入射窗是由厚度80～90μgcm~2聚乙烯甲醛制成的。该计数管的特点:(1)薄窗,对软X射线透过率高。(2)流气式,工作寿命长。(3)能量分辨率好。(4)计数率高(1×10~(14)个/s)。(5)可测能区宽(0.183～10keV)。(6)可以方便更换窗膜材料、厚度及窗口直径。近几年来该计数管已经为高强度低能X光源提供较好监测。     

Maximum likelihood estimate for discontinuous parameter in stochastic hyperbolic systems

The purpose of this paper is to study the identification problem of a spatially varying discontinuous parameter in stochastic hyperbolic equations. In previous works, the consistency property of the maximum likelihood estimate (MLE) was explored and the generating algorithm for MLE proposed under the condition that an unknown parameter is in a sufficiently regular space with respect to spatial variables.In order to show the consistency property of the MLE for a discontinuous coefficient, we use the method of sieves, i.e. the admissible class of unknown parameters is projected into a finite-dimensional space. For hyperbolic systems, we cannot obtain a regularity property of the solution with respect to a parameter. So in this paper, the parabolic regularization technique is used. The convergence of the derived finite-dimensional MLE to the infinite-dimensional MLE is justified under some conditions.     

Nonstationary curves in Hilbert spaces and their correlation functions II

Partially supportted by Grant MM-22/1991 of MESC     

Optical current transducer with bulk type Bi12SiO20 faraday sensor for power systems

Two types of optical current transducers (OCTs) have a bulk Faraday sensor inserted into the gap of an iron core and a porcelain insulator with optical fiber. The sensor consists of Bi₁₂SiO₂₀ (BSO) single crystal, a polarizer, and an analyzer. The OCTs satisfied the target performance requirement for fault location and metering and demonstrated maintained performance at some power utilities in Japan and the US. We have developed a fault location system that immediately detects the fault current with the OCTs, there by locating the fault section. The OCT can easily replace the existing support insulators for the disconnecting switch without any modifications to structure height or bus-bar. For metering requiring 0.3% class accuracy, use of a BSO with right optical rotatory power combined with BSO with left optical rotatory power results in a Faraday sensor with improved temperature characteristics. The OCT demonstrated 0.3% class accuracy for metering described in the current transformer Specifications of IEEE C57–13, 1993.     

Nonnegative quadratic estimation and quadratic sufficiency in general linear models

Notions of linear sufficiency and quadratic sufficiency are of interest to some authors. In this paper, the problem of nonnegative quadratic estimation for β^′Hβ+hσ² is discussed in a general linear model and its transformed model. The notion of quadratic sufficiency is considered in the sense of generality, and the corresponding necessary and sufficient conditions for the transformation to be quadratically sufficient are investigated. As a direct consequence, the result on (ordinary) quadratic sufficiency is obtained. In addition, we pose a practical problem and extend a special situation to the multivariate case. Moreover, a simulated example is conducted, and applications to a model with compound symmetric covariance matrix are given. Finally, we derive a remark which indicates that our main results could be extended further to the quasi-normal case.     

On the convergence of Newton's method when estimating higher dimensional parameters

In this paper, we consider the estimation of a parameter of interest where the estimator is one of the possibly several solutions of a set of nonlinear empirical equations. Since Newton's method is often used in such a setting to obtain a solution, it is important to know whether the so obtained iteration converges to the locally unique consistent root to the aforementioned parameter of interest. Under some conditions, we show that this is eventually the case when starting the iteration from within a ball about the true parameter whose size does not depend on n. Any preliminary almost surely consistent estimate will eventually lie in such a ball and therefore provides a suitable starting point for large enough n. As examples, we will apply our results in the context of M-estimates, kernel density estimates, as well as minimum distance estimates.