On surface radiation conditions for high-frequency wave scattering

A new approximation of the logarithmic derivative of the Hankel function is derived and applied to high-frequency wave scattering. We re-derive the on surface radiation condition (OSRC) approximations that are well suited for a Dirichlet boundary in acoustics. These correspond to the Engquist–Majda absorbing boundary conditions. Inverse OSRC approximations are derived and they are used for Neumann boundary conditions. We obtain an implicit OSRC condition, where we need to solve a tridiagonal system. The OSRC approximations are well suited for moderate wave numbers. The approximation of the logarithmic derivative is also used for deriving a generalized physical optics approximation, both for Dirichlet and Neumann boundary conditions. We have obtained similar approximations in electromagnetics, for a perfect electric conductor. Numerical computations are done for different objects in 2D and 3D and for different wave numbers. The improvement over the standard physical optics is verified.     

Classical solution of a one-dimensional parabolic boundary value problem with nonlinear boundary conditions and moving boundary

We consider a nonlinear parabolic boundary value problem of the Stefan type with one space variable, which generalizes the model of hydride formation under constant conditions. We suggest a grid method for constructing approximations to the unknown boundary and to the concentration distribution. We prove the uniform convergence of the interpolation approximations to a classical solution of the boundary value problem. (The boundary is smooth, and the concentration distribution has the necessary derivatives.) Thus, we prove the theorem on the existence of a solution, and the proof is given in constructive form: the suggested convergent grid method can be used for numerical experiments.     

Quartic spline method for solving fourth order obstacle boundary value problems

In this paper, we use uniform quartic polynomial splines to develop a new method, which is used for computing approximations to the solution and its first, second as well as third derivatives for a system of fourth order boundary value problems associated with obstacle, unilateral and contact problems. It is shown that the present method is of order two and gives approximations which are better than those produced by other collocation and finite difference methods. Numerical examples are presented to illustrate the applicability of the new method.     

Some mathematical problems in seriation

Some mathematical aspects of seriation are studied in this paper. Certain conditions on an abundance or an incidence matrix have been given in the past which imply that there exists a permutation of its rows so that the resulting matrix is a Q matrix (in which case the original matrix is said to be a pre-Q). These types of results have applications to chronologically ordering archaeological provenances under certain circumstances. Unfortunately these conditions are deficient both theoretically and practically, in that for much archaeological data the conditions are not necessarily true yet the corresponding provenances do have chronological orderings. Here we are able to generalize these results in two ways. First we are able to establish necessary and sufficient conditions on the rows of a matrix for it to be pre-Q. These conditions are local in that they concern only certain triples and quadruples of the rows. Secondly, we are able to interpret seriation in terms of a ternary relation R on a set A and prove the results in this general context. In this form the theorem says that if only certain of the triples and quadruples are R-strings, then the whole set A is an R-string, and so has a linear order consistent with the ternary relation R. This would appear to generalize a theorem of P. C. Fishburn. Both aspects of the generalization mean that the results stated herein have a wider applicability than those given heretofore. Possibly more importantly than this is that they lead to numerical invariants, called the fixing number and the related linear rigidity, of such an R-string on A. The archaeological interpretation of these is given in the paper and data supplied which illustrates this point. Finally various other conditions on products and representations of relations are stated which imply that A is an R-string. One of these generalizes and completes a theorem of D. G. Kendall.     

A parameter robust Petrov–Galerkin scheme for advection–diffusion–reaction equations

In this paper, a singularly perturbed convection diffusion boundary value problem, with discontinuous diffusion coefficient is examined. In addition to the presence of boundary layers, strong and weak interior layers can also be present due to the discontinuities in the diffusion coefficient. A priori layer adapted piecewise uniform meshes are used to resolve any layers present in the solution. Using a Petrov–Galerkin finite element formulation, a fitted finite difference operator is shown to produce numerical approximations on this fitted mesh, which are uniformly second order (up to logarithmic terms) globally convergent in the pointwise maximum norm.     

一类具非线性边界条件的高阶方程的奇摄动问题

通过引入伸展变量和非常规的渐近序列{∈}）,运用合成展开法,对一类具非线性边界条件的非线性高阶微分方程的奇摄动问题构造了形式渐近解,再运用微分不等式理论证明了原问题解的存在性及所得渐近近似式的一致有效性．     

Integration by Parts Formulae for Wiener Measures on a Path Space between two Curves

This paper is concerned with the integration by parts formulae for the pinned or the standard Wiener measures restricted on a space of paths staying between two curves. The boundary measures, concentrated on the set of paths touching one of the curves once, are specified. Our approach is based on the polygonal approximations. In particular, to establish the convergence of boundary terms, a uniform estimate is derived by means of comparison argument for a sequence of random walks conditioned to stay between two polygons. Applying the Brascamp–Lieb inequality, the stochastic integrals of Wiener type are constructed relative to the three-dimensional Bessel bridge or the Brownian meander. Supported in part by the JSPS Grant (B)(1)14340029     

An energy-momentum conserving scheme for Hamiltonian wave equation based on multiquadric trigonometric quasi-interpolation

Based on multiquadric trigonometric quasi-interpolation, the paper proposes a meshless symplectic scheme for Hamiltonian wave equation with periodic boundary conditions. The scheme first discretizes the equation in space using an iterated derivative approximation method based on multiquadric trigonometric quasi-interpolation and then in time with an appropriate symplectic scheme. This in turn yields a finite-dimensional semi-discrete Hamiltonian system whose energy and momentum (approximations of the continuous ones) are invariant with respect to time. The key feature of the scheme is that it conserves both the energy and momentum of the Hamiltonian system for both uniform and scattered centers, while classical energy-momentum conserving schemes are only for uniform centers. Numerical examples provided at the end of the paper show that the scheme is efficient and easy to implement.     

Polynomial Basis Multiplication over GF(2 m )

In this paper, we describe, analyze and compare various multipliers. Particularly, we investigate the standard modular multiplication, the Montgomery multiplication, and the matrix–vector multiplication techniques.     

Positive Sampling in Wavelet Subspaces

This paper is devoted to the discussion of a “hybrid” sampling series, a series of translates of a nonnegative summability function used in place of an orthogonal scaling function. The coefficients in the series are taken to be sampled values of the function to be approximated. This enables one to avoid the integration which arises in the other series. The approximations based on this hybrid series have certain desirable convergence properties: they are locally uniformly convergent for locally continuous functions, they have quadratic uniform convergence rate for functions in certain Sobolev spaces, they are locally bounded when the function is locally bounded and therefore, in particular, Gibbs' phenomenon is avoided. Numerical experiments are given to illustrate the theoretical results and to compare these approximations with the scaling function approximations.     

A collocation method for solving a class of singular nonlinear two-point boundary value problems

In this paper, an iterative algorithm for solving singular nonlinear two-point boundary value problems is formulated. This method is basically a collocation method for nonlinear second-order two-point boundary value problems with singularities at either one or both of the boundary points. It is proved that the iterative algorithm converges to a smooth approximate solution of the BVP provided the boundary value problem is well posed and the algorithm is applied appropriately. Error estimates for uniform partitions are also investigated. It has been shown that, for sufficiently smooth solutions, the method produces order h⁴ approximations. Numerical examples are provided to show the effectiveness of the algorithm.     

On the existence of mosaic-skeleton approximations for discrete analogues of integral operators

Exterior three-dimensional Dirichlet problems for the Laplace and Helmholtz equations are considered. By applying methods of potential theory, they are reduced to equivalent Fredholm boundary integral equations of the first kind, for which discrete analogues, i.e., systems of linear algebraic equations (SLAEs) are constructed. The existence of mosaic-skeleton approximations for the matrices of the indicated systems is proved. These approximations make it possible to reduce the computational complexity of an iterative solution of the SLAEs. Numerical experiments estimating the capabilities of the proposed approach are described.     

The design of linear algebra and geometry

Conventional formulations of linear algebra do not do justice to the fundamental concepts of meet, join, and duality in projective geometry. This defect is corrected by introducing Clifford algebra into the foundations of linear algebra. There is a natural extension of linear transformations on a vector space to the associated Clifford algebra with a simple projective interpretation. This opens up new possibilities for coordinate-free computations in linear algebra. For example, the Jordan form for a linear transformation is shown to be equivalent to a canonical factorization of the unit pseudoscalar. This approach also reveals deep relations between the structure of the linear geometries, from projective to metrical, and the structure of Clifford algebras. This is apparent in a new relation between additive and multiplicative forms for intervals in the cross-ratio. Also, various factorizations of Clifford algebras into Clifford algebras of lower dimension are shown to have projective interpretations.As an important application with many uses in physics as well as in mathematics, the various representations of the conformal group in Clifford algebra are worked out in great detail. A new primitive generator of the conformal group is identified.     

On certain order constrained Chebyshev rational approximations

Some rational approximations which share the properties of Padé and best uniform approximations are considered. The approximations are best in the Chebyshev sense, but the optimization is performed over subsets of the rational functions which have specified derivatives at one end point of the approximation interval. Explicit relationships between the Padé and uniform approximations are developed assuming the function being approximated satisfies easily verified constraints. The results are applied to the exponential function to determine the existence of best uniform A-acceptable approximations.     

On the uniqueness of minimizing harmonic maps to a closed hemisphere

We study the uniqueness of minimizing harmonic maps to a closed hemisphere. We are able to describe the boundary data for which there are more than one minimizer, and to describe in these cases the corresponding set of minimizers. This is a limiting case for a former result of W.Jäger and H.Kaul.This article was processed by the author using the Springer-Verlag T_EX PJourlg macro package 1991.     

Scattering of elastic waves by nonaxisymmetric three-dimensional dipping layer

Using a boundary method, we investigated the scattering of elastic plane harmonic SH, SV, P, and Rayleigh waves by three-dimensional nonaxisymmetric dipping layers embedded in an elastic half-space. The valley was subjected to incident Rayleigh wave and oblique incident SH, SV, and P waves. The method utilized spherical wave functions to express the unknown scattered field. These functions satisfy the equation of motion and radiation conditions at infinity but they do not satisfy the stress-free boundary conditions at the surface of the half-space. The boundary and continuity conditions are imposed locally in the least-square-sense at several points on the layer interface and on the surface of the half-space. A comparative study was done to examine the validity and limitations of the two-dimensional approximations (antiplane and plane strain models) of three-dimensional models. It is demonstrated that the two-dimensional approximations may be inadequate to represent actual displacement field for three-dimensional irregularities.     

Cubic spline method for solving two-point boundary-value problems

In this paper, we use uniform cubic spline polynomials to derive some new consistency relations. These relations are then used to develop a numerical method for computing smooth approximations to the solution and its first, second as well as third derivatives for a second order boundary value problem. The present method outperforms other collocations, finite-difference and splines methods of the same order. Numerical illustrations are provided to demonstrate the practical use of our method.     

Best approximation and moduli of smoothness: Computation and equivalence theorems

In this paper we investigate the trigonometric series with the β-general monotone coefficients. First, we study the uniform convergence criterion. The estimates of best approximations and moduli of smoothness of the series in uniform metrics are obtained in terms of coefficients. These results imply several important relations between moduli of smoothness of different orders (in particular, Marchaud-type inequality) and best approximations.     

Convergence of the numerical solution of the initial-boundary value problem for a heat equation with a corner singularity in derivatives of the solution

An estimate O(τ + h ²)ln(j + 1) of the convergence rate for the solution of a four-point implicit difference scheme used for approximations on a uniform grid of a one-dimensional heat equation is obtained, with the provision that the boundary and initial data are subject at corner points only to the continuity condition, with no other compatibility conditions being satisfied. A discrete Green’s function is used to obtain an a priori estimate of the grid solution in terms of the appropriate negative norm of the right-hand side.     

Bundle automorphisms for fiberG-structures of finite type

In this paper we study geometric settings where a Lie group preserving a measurable field of measurable Riemannian metrics on the fibers of a smooth fiber bundle must actually preserve a measurable field of smooth Riemannian metrics. For ergodic actions on bundles with compact fiber this will imply that the standard fiber is a homogeneous space for a compact Lie group. In particular we show this conclusion holds for a semisimple Lie group of higher real rank (or a lattice subgroup) preserving a finite measure and either a field of connections or pseudo-Riemannian metrics when the fiber is compact and of low dimension.Research completed while a member of the University of Chicago Mathematics Department.