The Half-plane with an Edge Crack in Plane Elastostatics

The Wiener-Hopf technique is applied to the problem of a half-planewith a crack perpendicular to its edge. A closed solution isobtained giving the stress and displacement fields in the formof eigenfunction expansions, valid over the whole region, exceptthe neighbourhood of the tip of the crack. The behaviour inthis region is obtained by an asymptotic technique. The coefficientsof the expansions involve a generalized factorial function,and a method for the computation of this function is given,thus enabling numerical values of the stresses and displacementsto be obtained for the given cases of stretching and bendingthe half-plane.     

Self-adjoint Extensions for the Neumann Laplacian and Applications

A new technique is proposed for the analysis of shape optimization problems. The technique uses the asymptotic analysis of boundary value problems in singularly perturbed geometrical domains. The asymptotics of solutions are derived in the framework of compound and matched asymptotics expansions. The analysis involves the so-called interior topology variations. The asymptotic expansions are derived for a model problem, however the technique applies to general elliptic boundary value problems. The self-adjoint extensions of elliptic operators and the weighted spaces with detached asymptotics are exploited for the modelling of problems with small defects in geometrical domains, The error estimates for proposed approximations of shape functionals are provided.     

Matched Asymptotic Expansions in the Two-body Problem with Quick Loss of Mass

In studying models for the two-body problem with quick lossof mass a boundary layer problem arises for a third-order systemof non-linear ordinary differential equations. The models areidentified by a real parameter n with n ? 1. It turns out thatfor n = 1 asymptotic approximations of the solutions can beobtained by applying the method of matched asymptotic expansionsaccouonding to Vasil'eva or a multiple time scales method developedby O'Malley. For n> 1 these methods break down and it isshown that this is due to the occurrence of "unexpected" orderfunctions in the asymptotic expansions. The expansions for n> 1 are obtained by constructing an inner and outer expansionof the solution and matching these by the process of takingintermediate limits. The asymptotic validity of the matched expansions is provedby using an iteration technique; the proof is constructive sothat it provides us at the same time with an alternative wayof constructing approximations without using a matching technique.     

On canonicity of poset expansions

In this paper, we discuss the canonical extension of poset expansions. To obtain canonicity results on poset expansions, we study Ghilardi and Meloni’s canonicity methodology for Heyting algebras with unary modalities, raise the problem of extending the technique to poset expansions, and give a possible solution for the problem. Finally, we obtain a syntactic account of canonical inequalities on poset expansions consisting of constants, e_^{\epsilon_{\bot}} -additive operations, e^T{\epsilon^{\top}} -multiplicative operations, diamond, box, and strict adjoint pairs, and bounded poset expansions consisting of constants, e{\epsilon} -join preserving operations, e{\epsilon} -meet preserving operations, e{\epsilon} -additive operations, e{\epsilon} -multiplicative operations and adjoint pairs, which are more restricted than the case of lattice expansions, but can still account for Sahlqvist-like canonicity results.     

Second order probabilistic parametrix method for unbiased simulation of stochastic differential equations

In this article, following the paradigm of bias–variance trade-off philosophy, we derive parametrix expansions of order two, based on the Euler–Maruyama scheme with random partitions, for the purpose of constructing an unbiased simulation method for multidimensional stochastic differential equations. These formulas lead to Monte Carlo simulation methods which can be easily parallelized. The second order method proposed here requires further regularity of coefficients in comparison with the first order method but achieves finite moments even when Poisson sampling is used for the partitions, in contrast to Andersson and Kohatsu-Higa (2017). Moreover, using an exponential scaling technique one achieves an unbiased simulation method which resembles a space importance sampling technique which significantly improves the efficiency of the proposed method. A hint of how to derive higher order expansions is also presented.     

A Chebyshev expansion method for solving Nonlinear circuit equations

A new method using Chebyshev expansions is developed for solving nonlinear circuit equations. This technique is more efficient and accurate than existing simulators. Two particular case studies illustrate the methodology, and further applications are outlined.     

The Asymptotic Expansion and Numerical Verification of Van der Pol's Equation

A technique is presented in this paper to verify the order of accuracy of asymptotic expansion of Van der Pol's equation. The technique is focused on using numerical solutions as an independent means of verifying the validity of asymptotic expansions.     

Generalized contour models

The paper is devoted to a survey and systematic exposition of the technique of obtaining cluster expansions for lattice Gibbs fields in the low-temperature region in the case of a finite or countable number of basic states.     

Asymptotic error expansions for numerical solutions of one-dimensional problems with singularities

A systematic analysis is given on asymptotic error expansions for numerical solutions of one-dimensional problems whose solutions are singular. Numerical examples show a great improvement on the accuracy of numerical solutions by using the Richardson extrapolation technique.     

Perturbation analysis of singular subspaces and deflating subspaces

Summary. Perturbation expansions for singular subspaces of a matrix and for deflating subspaces of a regular matrix pair are derived by using a technique previously described by the author. The perturbation expansions are then used to derive Fr\'echet derivatives, condition numbers, and th-order perturbation bounds for the subspaces. Vaccaro's result on second-order perturbation expansions for a special class of singular subspaces can be obtained from a general result of this paper. Besides, new perturbation bounds for singular subspaces and deflating subspaces are derived by applying a general theorem on solution of a system of nonlinear equations. The results of this paper reveal an important fact: Each singular subspace and each deflating subspace have individual perturbation bounds and individual condition numbers. Received July 26, 1994     

Determination of periodic orbits of nonlinear oscillators by means of piecewise-linearization methods

An approximate method based on piecewise linearization is developed for the determination of periodic orbits of nonlinear oscillators. The method is based on Taylor series expansions, provides piecewise analytical solutions in three-point intervals which are continuous everywhere and explicit three-point difference equations which are P-stable and have an infinite interval of periodicity. It is shown that the method presented here reduces to the well-known Störmer technique, is second-order accurate, and yields, upon applying Taylor series expansion and a Padé approximation, another P-stable technique whenever the Jacobian is different from zero. The method is generalized for single degree-of-freedom problems that contain the velocity, and (approximate) analytical solutions are presented. Finally, by introducing the inverse of a vector and the vector product and quotient, and using Taylor series expansions and a Padé approximation, the method has been generalized to multiple degree-of-freedom problems and results in explicit three-point finite difference equations which only involve vector multiplications.     

Asymptotic Solutions of a Fourth Order Differential Equation

In this paper, we derive uniform asymptotic expansions of solutions to the fourth order differential equation where x is a real variable and λ is a large positive parameter. The solutions of this differential equation can be expressed in the form of contour integrals, and uniform asymptotic expansions are derived by using the cubic transformation introduced by Chester, Friedman, and Ursell in 1957 and the integration-by-part technique suggested by Bleistein in 1966. There are two advantages to this approach: (i) the coefficients in the expansion are defined recursively, and (ii) the remainder is given explicitly. Moreover, by using a recent method of Olde Daalhuis and Temme, we extend the validity of the uniform asymptotic expansions to include all real values of x .     

A numerical approach to cyclic-service queueing models

An iterative numerical technique for the evaluation of queue length distributions is applied to multi-queue systems with one server and cyclic service discipline with Bernoulli schedules. The technique is based on power-series expansions of the state probabilities as functions of the load of the system. The convergence of the series is accelerated by applying a modified form of the epsilon algorithm. Attention is paid to economic use of memory space.     

多孔介质中非Fick流的校正方法

本文首先给出H~1-模意义下多孔介质中非Fick流的矩形双线性元的渐进误差展开,进而通过插值后处理方法得到一种插值校正格式来提高有限元近似解的精度.     

Asymptotic expansions using blow-up

The method of matched asymptotic expansions and geometric singular perturbation theory are the most common and successful approaches to singular perturbation problems. In this work we establish a connection between the two approaches in the context of the simple fold problem. Using the blow-up technique [5], [12] and the tools of geometric singular perturbation theory we derive asymptotic expansions of slow manifolds continued beyond the fold point. Our analysis explains the structure of the expansion and gives an algorithm for computing its coefficients.     

Asymptotic expansions of integrals of two Bessel functions via the generalized hypergeometric and Meijer functions

Asymptotic expansions of certain finite and infinite integrals involving products of two Bessel functions of the first kind are obtained by using the generalized hypergeometric and Meijer functions. The Bessel functions involved are of arbitrary (generally different) orders, but of the same argument containing a parameter which tends to infinity. These types of integrals arise in various contexts, including wave scattering and crystallography, and are of general mathematical interest being related to the Riemann—Liouville and Hankel integrals. The results complete the asymptotic expansions derived previously by two different methods — a straightforward approach and the Mellin-transform technique. These asymptotic expansions supply practical algorithms for computing the integrals. The leading terms explicitly provide valuable analytical insight into the high-frequency behavior of the solutions to the wave-scattering problems.     

A Note on Simple Non-Zero Singular Values

The technique described in [4] is used to investigate the analyticity and to obtain second order perturbation expansions of simple non-zero singular values of a matrix analytically dependent on several parameters.     

Perturbation method for solving the nonlinear eigenvalue problem arising from fatigue crack growth problem in a damaged medium

An analytical solution of the nonlinear eigenvalue problem arising from the fatigue crack growth problem in a damaged medium in coupled formulation is obtained. The perturbation technique for solving the nonlinear eigenvalue problem is used. The method allows to find the analytical formula expressing the eigenvalue as the function of parameters of the damage evolution law. It is shown that the eigenvalues of the nonlinear eigenvalue problem are fully determined by the exponents of the damage evolution law. In the paper the third-order (four-term) asymptotic expansions of the angular functions determining the stress and continuity fields in the neighborhood of the crack tip are given. The asymptotic expansions of the angular functions permit to find the closed-form solution for the problem considered.     

On g-functions for Hermite function expansions

Summary We prove weighted L^p-inequalities for the gradient square function associated with the Poisson semigroup in the multi-dimensional Hermite function expansions setting. In the proof a technique of vector valued Calderón-Zygmund operators is used.     

Greedy expansions in convex optimization

This paper is a follow-up to the author’s previous paper on convex optimization. In that paper we began the process of adjusting greedy-type algorithms from nonlinear approximation for finding sparse solutions of convex optimization problems. We modified there the three most popular greedy algorithms in nonlinear approximation in Banach spaces-Weak Chebyshev Greedy Algorithm, Weak Greedy Algorithm with Free Relaxation, and Weak Relaxed Greedy Algorithm-for solving convex optimization problems. We continue to study sparse approximate solutions to convex optimization problems. It is known that in many engineering applications researchers are interested in an approximate solution of an optimization problem as a linear combination of elements from a given system of elements. There is an increasing interest in building such sparse approximate solutions using different greedy-type algorithms. In this paper we concentrate on greedy algorithms that provide expansions, which means that the approximant at the mth iteration is equal to the sum of the approximant from the previous, (m ? 1)th, iteration and one element from the dictionary with an appropriate coefficient. The problem of greedy expansions of elements of a Banach space is well studied in nonlinear approximation theory. At first glance the setting of a problem of expansion of a given element and the setting of the problem of expansion in an optimization problem are very different. However, it turns out that the same technique can be used for solving both problems. We show how the technique developed in nonlinear approximation theory, in particular, the greedy expansions technique, can be adjusted for finding a sparse solution of an optimization problem given by an expansion with respect to a given dictionary.