Detection of small inhomogeneities via direct sampling method in transverse electric polarization

Various studies have confirmed the possibility of identifying the location of a set of small inhomogeneities via a direct sampling method; however, when their permeability differs from that of the background, their location cannot be satisfactorily identified. However, no theoretical explanation for this phenomenon has been verified. In this study, we demonstrate that the indicator function of the direct sampling method can be expressed by the Bessel function of order one of the first kind and explain why the exact locations of inhomogeneities cannot be identified. Numerical results with noisy data are exhibited to support our examination.     

Improvement of direct sampling method in transverse electric polarization

Recently, the mathematical structure of the direct sampling method (DSM) was investigated in transverse-electric (TE) polarization, and the reason why the exact locations of small inhomogeneities cannot be retrieved through traditional DSM has been revealed. In this paper, we present an improved DSM for identifying the exact locations of small inhomogeneities in TE polarization. Furthermore, we investigate a multi-frequency indicator function to obtain a better result. Corresponding mathematical analysis and simulations are performed to show the feasibilities of the proposed improvement techniques.     

Two effective post-filtering strategies for improving direct sampling methods

In this paper, we introduce two novel strategies to reduce artifacts in the direct sampling type methods (DSM). The newly proposed techniques can essentially reduce the artifacts and provide more accurate and reliable physical profiles of the scatterers compared with the original DSM. The techniques can find wide applications in the inverse scattering problems. Moreover, the novel techniques exhibit several strengths: direct, stable, robust, and ease of implementation.     

Spectral analysis of the finite Hankel transform and circular prolate spheroidal wave functions

In this paper, we develop two practical methods for the computation of the eigenvalues as well as the eigenfunctions of the finite Hankel transform operator. These different eigenfunctions are called circular prolate spheroidal wave functions (CPSWFs). This work is motivated by the potential applications of the CPSWFs as well as the development of practical methods for computing their values. Also, in this work, we should prove that the CPSWFs form an orthonormal basis of the space of Hankel band-limited functions, an orthogonal basis of L²([0,1]) and an orthonormal system of L²([0,+∞[). Our computation of the CPSWFs and their associated eigenvalues is done by the use of two different methods. The first method is based on a suitable matrix representation of the finite Hankel transform operator. The second method is based on the use of an efficient quadrature method based on a special family of orthogonal polynomials. Also, we give two Maple programs that implement the previous two methods. Finally, we present some numerical results that illustrate the results of this work.     

Bounds on differences of adjacent zeros of Bessel functions and iterative relations between consecutive zeros

Bounds for the distance between adjacent zeros of cylinder functions are given; and are such that ; stands for the th positive zero of the cylinder (Bessel) function , , .

These bounds, together with the application of modified (global) Newton methods based on the monotonic functions and , give rise to forward ( ) and backward ( ) iterative relations between consecutive zeros of cylinder functions.

The problem of finding all the positive real zeros of Bessel functions for any real and inside an interval , 0$">, is solved in a simple way.

     

A global Newton method for the zeros of cylinder functions

The zeros of cylinder functions C _u (x)=cos α, J _u (x) - sin α, Y _u(x) coincide with those of the ratios H _u (x)=C _u (x)/C _u-1 (x) except, perhaps, at x = 0. We show monotonicity properties of H u(x) and f _u (x) = x ^2v-1 H u(x) and their derivatives for x > 0. We then build a Newton-Raphson iterative method based on the monotonic function f u(x) which is shown to be convergent, for any real values of u and α and any starting value x ₀ > 0, to an sth positive root c ,s of C _u (x) = 0, s being such that c _,s and x0 belong to the same interval (c _u-1 ,s', c _{u -1} ,s'+1]. We also show applications of the method. In particular, taking advantage of the fact that the ratio H _u (x) for first kind Bessel functions J u(x) can be evaluated by using a continued fraction, a very simple algorithm is built; it becomes especially efficient for low values of u and s and it allows the evaluation of the real zeros for arbitrary orders u, positive or negative. This revised version was published online in August 2006 with corrections to the Cover Date.     

Application of lower bound direct method to engineering structures

Direct methods provide elegant and efficient approaches for the prediction of the long-term behaviour of engineering structures under arbitrary complex loading independent of the number of loading cycles. The lower bound direct method leads to a constrained non-linear convex problem in conjunction with finite element methods, which necessitates a very large number of optimization variables and a large amount of computer memory. To solve this large-scale optimization problem, we first reformulate it in a simpler equivalent convex program with easily exploitable sparsity structure. The interior point with DC regularization algorithm (IPDCA) using quasi definite matrix techniques is then used for its solution. The numerical results obtained by this algorithm will be compared with those obtained by general standard code Lancelot. They show the robustness, the efficiency of IPDCA and in particular its great superiority with respect to Lancelot.     

On sampling theory associated with the resolvents of singular Sturm-Liouville problems

This paper is concerned with the sampling theory associated with resolvents of eigenvalue problems. We introduce sampling representations for integral transforms whose kernels are Green's functions of singular Sturm-Liouville problems provided that the singular points are in the limit-circle situation, extending the results obtained in the regular problems.

     

Integrals involving products of Airy functions, their derivatives and Bessel functions

A new integral representation of the Hankel transform type is deduced for the function Fn(x,Z)=Zn−1Ai(x−Z)Ai(x+Z) with x∈R, Z>0 and n∈N. This formula involves the product of Airy functions, their derivatives and Bessel functions. The presence of the latter allows one to perform various transformations with respect to Z and obtain new integral formulae of the type of the Mellin transform, K-transform, Laplace and Fourier transform. Some integrals containing Airy functions, their derivatives and Chebyshev polynomials of the first and second kind are computed explicitly. A new representation is given for the function ²|Ai(z)| with z∈C.     

A Bessel collocation method for numerical solution of generalized pantograph equations

This article is concerned with a generalization of a functional differential equation known as the pantograph equation which contains a linear functional argument. In this article, we introduce a collocation method based on the Bessel polynomials for the approximate solution of the pantograph equations. The method is illustrated by studying the initial value problems. The results obtained are compared by the known results. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Mathematical basis of the linear sampling method for partially coated obstacles

We consider the inverse scattering problem of determining the shape of a partially coated obstacle D. To this end, we solve a scattering problem for the Helmholtz equation where the scattered field satisfies mixed Dirichlet–Neumann-impedance boundary conditions on the Lipschitz boundary of the scatterer D. Based on the analysis of the boundary integral system to the direct scattering problem, we propose how to reconstruct the shape of the obstacle D by using the linear sampling method.     

An improved time domain linear sampling method for Robin and Neumann obstacles

We consider inverse obstacle scattering problems for the wave equation with Robin or Neumann boundary conditions. The problem of reconstructing the geometry of such obstacles from measurements of scattered waves in the time domain is tackled using a time domain linear sampling method. This imaging technique yields a picture of the scatterer by solving a linear operator equation involving the measured data for many right-hand sides given by singular solutions to the wave equation. We analyse this algorithm for causal and smooth impulse shapes, we discuss the effect of different choices of the singular solutions used in the algorithm, and finally we propose a fast FFT-based implementation.     

Numerical analysis of a fast integration method for highly oscillatory functions

The integration of systems containing Bessel functions is a central point in many practical problems in physics, chemistry and engineering. This paper presents a new numerical analysis for the collocation method presented by Levin for and gives more accurate error analysis about the integration of systems containing Bessel functions. The effectiveness and accuracy of the quadrature is tested for Bessel functions with large arguments. AMS subject classification (2000)  65D32, 65D30     

Computation of the eigenvalues of Sturm-Liouville problems with parameter dependent boundary conditions using the regularized sampling method

The purpose in this paper is to compute the eigenvalues of Sturm-Liouville problems with quite general separated boundary conditions nonlinear in the eigenvalue parameter using the regularized sampling method, an improvement on the method based on Shannon sampling theory, which does not involve any multiple integration and provides higher order estimates of the eigenvalues at a very low cost. A few examples shall be presented to illustrate the power of the method and a comparison made with the the exact eigenvalues obtained as squares of the zeros of the exact characteristic functions.

     

Interval optimization based line sampling method for fuzzy and random reliability analysis

For structural system with fuzzy variables as well as random variables, a novel algorithm for obtaining membership function of fuzzy reliability is presented on interval optimization based Line Sampling (LS) method. In the presented algorithm, the value domain of the fuzzy variables under the given membership level is firstly obtained according to their membership functions. Then, in the value domain of the fuzzy variables, bounds of reliability of the structure are obtained by the nesting analysis of the interval optimization, which is performed by modern heuristic methods, and reliability analysis, which is achieved by the LS method in the reduced space of the random variables. In this way the uncertainties of the input variables are propagated to the safety measurement of the structure, and the membership function of the fuzzy reliability is obtained. The presented algorithm not only inherits the advantage of the direct Monte Carlo method in propagating and distinguishing the fuzzy and random uncertainties, but also can improve the computational efficiency tremendously in case of acceptable precision. Several examples are used to illustrate the advantages of the presented algorithm.     

On the optimal control method in quaternionic analysis

Gaveau?s optimal control method for real and complex Monge–Ampere operators is generalized to that for quaternionic Monge–Ampere operator. It is also applied to investigate quaternionic regular functions: the characterization of the Silov boundary of a smooth quaternionic pseudoconvex domain.     

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 systematization of the saddle point method. Application to the Airy and Hankel functions

The standard saddle point method of asymptotic expansions of integrals requires to show the existence of the steepest descent paths of the phase function and the computation of the coefficients of the expansion from a function implicitly defined by solving an inversion problem. This means that the method is not systematic because the steepest descent paths depend on the phase function on hand and there is not a general and explicit formula for the coefficients of the expansion (like in Watson's Lemma for example). We propose a more systematic variant of the method in which the computation of the steepest descent paths is trivial and almost universal: it only depends on the location and the order of the saddle points of the phase function. Moreover, this variant of the method generates an asymptotic expansion given in terms of a generalized (and universal) asymptotic sequence that avoids the computation of the standard coefficients, giving an explicit and systematic formula for the expansion that may be easily implemented on a symbolic manipulation program. As an illustrative example, the well-known asymptotic expansion of the Airy function is rederived almost trivially using this method. New asymptotic expansions of the Hankel function Hn(z) for large n and z are given as non-trivial examples.