A local error estimate of the method of multi-scale asymptotic expansions for elliptic problems with rapidly oscillatory coefficients

In this paper, on basis of [O.A. Oleinik, A.S. Shamaev, G.A. Yosifian, Mathematical Problems in Elasticity and Homogenization, North-Holland, Amsterdam, 1992], we present a local error estimate of the method of multi-scale asymptotic expansions for second order elliptic problems with rapidly oscillatory coefficients.     

Hyperasymptotic evaluation of the Pearcey integral via Hadamard expansions

We describe high-precision computations of the Pearcey integral Pe(x,y) for real x and y by means of Hadamard expansions. Numerical results for (x,y) situated in different regions of the x,y-plane are given to illustrate the levels of precision that can be achieved. Particular emphasis is given to computation in the neighbourhood of the two cusped curves associated with Pe(x,y) across which there is either a coalescence of saddles or a Stokes phenomenon.     

CBS constants for multilevel splitting of graph-Laplacian and application to preconditioning of discontinuous Galerkin systems

The goal of this work is to derive and justify a multilevel preconditioner of optimal arithmetic complexity for symmetric interior penalty discontinuous Galerkin finite element approximations of second order elliptic problems. Our approach is based on the following simple idea given in [R.D. Lazarov, P.S. Vassilevski, L.T. Zikatanov, Multilevel preconditioning of second order elliptic discontinuous Galerkin problems, Preprint, 2005]. The finite element space of piece-wise polynomials, discontinuous on the partition , is projected onto the space of piece-wise constant functions on the same partition that constitutes the largest space in the multilevel method. The discontinuous Galerkin finite element system on this space is associated to the so-called “graph-Laplacian”. In 2-D this is a sparse M-matrix with -1 as off diagonal entries and nonnegative row sums. Under the assumption that the finest partition is a result of multilevel refinement of a given coarse mesh, we develop the concept of hierarchical splitting of the unknowns. Then using local analysis we derive estimates for the constants in the strengthened Cauchy–Bunyakowski–Schwarz (CBS) inequality, which are uniform with respect to the levels. This measure of the angle between the spaces of the splitting was used by Axelsson and Vassilevski in [Algebraic multilevel preconditioning methods II, SIAM J. Numer. Anal. 27 (1990) 1569–1590] to construct an algebraic multilevel iteration (AMLI) for finite element systems. The main contribution in this paper is a construction of a splitting that produces new estimates for the CBS constant for graph-Laplacian. As a result we have a preconditioner for the system of the discontinuous Galerkin finite element method of optimal arithmetic complexity.     

Free Vibration Analysis of Rotating Nonlinearly Elastic Structures with Symmetry: An Efficient Group-Equivariance Approach

A study is made of the dynamics of oscillating systems with a slowly varying parameter. A slowly varying forcing periodically crosses a critical value corresponding to a pitchfork bifurcation. The instantaneous phase portrait exhibits a centre when the forcing does not exceed the critical value, and a saddle and two centres with an associated double homoclinic loop separatrix beyond this value. The aim of this study is to construct a Poincaré map in order to describe the dynamics of the system as it repeatedly crosses the bifurcation point. For that purpose averaging methods and asymptotic matching techniques connecting local solutions are applied. Given the initial state and the values of the parameters the properties of the Poincaré map can be studied. Both sensitive dependence on initial conditions and (quasi) periodicity are observed. Moreover, Lyapunov exponents are computed. The asymptotic expressions for the Poincaré map are compared with numerical solutions of the full system that includes small damping.     

Derivation of the Forchheimer Law Via Matched Asymptotic Expansions

Matched asymptotics are used formally in the linearized Navier–Stokes equations, the so-called Oseen approximation, to derive a nonlinear law for the flow of a single-phase, incompressible Newtonian fluid through a rigid porous medium. A suitable choice of the linear term gives, to order a Forchheimer-type law.     

Uniform asymptotics for the linearized Boltzmann equation describing sound wave propagation

The multiple-scale expansionmethod is used for constructing a uniformly applicable asymptotic approximation of the solution of the linearized Boltzmann equation for small Knudsen numbers. The asymptotic expansion is constructed for the particular example of a sound wave generated by a plane oscillation source and dissipating in a half-space. The simplicity of the problem makes it possible clearly to demonstrate the appearance of secular terms in the expansion and the introduction of multiple scales opens the way to eliminating them.     

A Cross-Entropy Approach to the Estimation of Generalized Linear Multilevel Models

In this article, we use the cross-entropy method for noisy optimization for fitting generalized linear multilevel models through maximum likelihood. We propose specifications of the instrumental distributions for positive and bounded parameters that improve the computational performance. We also introduce a new stopping criterion, which has the advantage of being problem-independent. In a second step we find, by means of extensive Monte Carlo experiments, the most suitable values of the input parameters of the algorithm. Finally, we compare the method to the benchmark estimation technique based on numerical integration. The cross-entropy approach turns out to be preferable from both the statistical and the computational point of view. In the last part of the article, the method is used to model the probability of firm exits in the healthcare industry in Italy. Supplemental materials are available online.     

Structure dependent sampling in compressed sensing: Theoretical guarantees for tight frames

Many of the applications of compressed sensing have been based on variable density sampling, where certain sections of the sampling coefficients are sampled more densely. Furthermore, it has been observed that these sampling schemes are dependent not only on sparsity but also on the sparsity structure of the underlying signal. This paper extends the result of Adcock, Hansen, Poon and Roman (arXiv:1302.0561, 2013) [2] to the case where the sparsifying system forms a tight frame. By dividing the sampling coefficients into levels, our main result will describe how the amount of subsampling in each level is determined by the local coherences between the sampling and sparsifying operators and the localized level sparsities – the sparsity in each level under the sparsifying operator.     

One‐dimensional q‐Dirac equation

In this paper, we introduce a q‐analog of 1‐dimensional Dirac equation. We investigate the existence and uniqueness of the solution of this equation. Later, we discuss some spectral properties of the problem, such as formally self‐adjointness, the case that the eigenvalues are real, orthogonality of eigenfunctions, Green function, existence of a countable sequence of eigenvalues, and eigenfunctions forming an orthonormal basis of . Finally, we give some examples.     

Asymptotic Axially Symmetric Deformations for Perfectly Elastic Neo-Hookean and Mooney Materials

For axially symmetric deformations of the perfectly elastic neo-Hookean and Mooney materials, formal series solutions are determined in terms of expansions in appropriate powers of 1/R, where R is the cylindrical polar coordinate for the material coordinates. Remarkably, for both the neo-Hookean and Mooney materials, the first three terms of such expansions can be completely determined analytically in terms of elementary integrals. From the incompressibility condition and the equilibrium equations, the six unknown deformation functions, appearing in the first three terms can be reduced to five formal integrations involving in total seven arbitrary constants A, B, C, D, E, H and k ², and a further five integration constants, making a total of 12 integration constants for the deformation field. The solutions obtained for the neo-Hookean material are applied to the problem of the axial compression of a cylindrical rubber tube which has bonded metal end-plates. The solution so determined is approximate in two senses; namely as an approximate solution of the governing equations and for which the stress free and displacement boundary conditions are satisfied in an average manner only. The resulting load-deflection relation is shown graphically. The solution so determined, although approximate, attempts to solve a problem not previously tackled in the literature.