Numerical simulation and analysis of generalized difference method on triangular networks for electrical impedance tomography

Electrical impedance tomography is an inverse problem of elliptic differential equations. Numerical methods based on combining generalized difference method and Levenberg–Marquardt iteration on a planar domain are proposed. Positive semi-definiteness and existence of solution of the generalized difference scheme are proved. Element geometry matrix is introduced to shortcut calculation and standardize computer program. A series of numerical experiments verify the reliability of its mathematical model and the feasibility of the algorithm. A class of electrical current patterns is proposed to minimize the number of direct problems to be solved in each iteration. These methods have been applied successfully in practical simulation of electrical impedance tomography.     

Identification of an inclusion in multifrequency electric impedance tomography

The multifrequency electrical impedance tomography is considered to image a conductivity inclusion inside a homogeneous background medium by injecting one current. An original spectral decomposition of the solution of the forward conductivity problem is used to retrieve the Cauchy data corresponding to the extreme case of perfect conductor. Using results based on the unique continuation, we then prove the uniqueness of multifrequency electrical impedance tomography and obtain rigorous stability estimates. Our results in this paper are quite surprising in inverse conductivity problem since in general infinitely many input currents are needed to obtain the uniqueness in the determination of the conductivity.     

Mixed impedance transmission problems for vibrating layered elastic bodies

Direct scattering problems for partially coated piecewise homogenous and inhomogeneous layered obstacles in linear elasticity lead to mixed impedance transmission problems for the steady‐state elastic oscillation equations. For a piecewise homogenous isotropic composite body, we employ the potential method and reduce the mixed impedance transmission problem to an equivalent system of boundary pseudodifferential equations. We give a detailed analysis of the corresponding pseudodifferential operators, which live on the interface between the layers and on a proper submanifold of the boundary of the composite elastic body, and establish uniqueness and existence results for the original mixed impedance transmission problem for arbitrary values of the oscillation frequency parameter; this is crucial in the study of inverse elastic scattering problems for partially coated layered obstacles. We also investigate regularity properties of solutions near the collision curves, where the different boundary conditions collide, and establish almost best Hölder smoothness results. Further, we analyze the asymptotic behavior of the stress vector near the collision curve and derive explicit formulas for the stress singularity exponents. The case of Lipschitz surfaces is briefly treated separately. In the case of a composite body containing homogeneous or inhomogeneous finite anisotropic inclusions, we develop an alternative hybrid method based on the so‐called nonlocal approach and reduce the mixed transmission problem to an equivalent functional‐variational equation with a sesquilinear form that ‘lives’ on a bounded part of the layered composite body and its boundary. We show that this sesquilinear form is coercive and that the corresponding variational equation is uniquely solvable. Copyright © 2015 John Wiley & Sons, Ltd.     

利用远场模式的完全与不完全数据反演声波阻尼区域

提出一种方法,利用远场模式的完全数据与不完全数据反演声波阻尼区域,证明了方法的收敛性,并给出若干数值例子.     

Electromagnetic scattering by a smooth convex impedance cone

The problem of the diffraction of an electromagnetic planewave by a convex cone of arbitrary smooth cross-section withimpedance (Leontovich) boundary conditions is studied. The vectorproblem is reduced to that for the Debye potentials. By meansof Kontorovich–Lebedev integrals, two spectral functionsare introduced and the corresponding boundary value problemis formulated. The spectral functions for the potentials arefound to satisfy the Helmholtz equations on the unit sphereand to be coupled through non-traditional boundary conditionsof the impedance type with shifts on the spectral variable.The use of the Green theorem permits us to establish an integralformulation of the boundary value problem for the spectral functions.The formal asymptotic solution of the problem is then givenfor the case of a narrow cone. For this, two different methodsare given: a method of perturbation applied to the spectralintegral equations and an adaptation of the method of matchingthe asymptotic series in spectral domain. Both methods leadto the same closed-form result for the leading term of the scatteringdiagram asymptotics.     

反演声波阻尼系数的另一种方法

提出了一种方法,利用正则化方法和积分方程,由散射波的近场数据反演时间调和声波阻尼系数.给出了该方法收敛性的证明及数值例子,算法与数值例子表明这种方法不仅简单而且很有效.     

Stability of the determination of the surface impedance of an obstacle from the scattering amplitude

We prove a stability estimate of logarithmic type for the inverse problem consisting in the determination of the surface impedance of an obstacle from the scattering amplitude. We present a simple and direct proof, which is essentially based on an elliptic Carleman inequality. Copyright © 2013 John Wiley & Sons, Ltd.     

Pulsatile flow of Herschel–Bulkley fluid through catheterized arteries – A mathematical model

The pulsatile flow of blood through catheterized artery has been studied in this paper by modeling blood as Herschel–Bulkley fluid and the catheter and artery as rigid coaxial circular cylinders. The Herschel–Bulkley fluid has two parameters, the yield stress θ and the power index n. Perturbation method is used to solve the resulting quasi-steady nonlinear coupled implicit system of differential equations. The effects of catheterization and non-Newtonian nature of blood on yield plane locations, velocity, flow rate, wall shear stress and longitudinal impedance of the artery are discussed. The existence of two yield plane locations is investigated and their dependence on yield stress θ, amplitude A, and time t are analyzed. The width of the plug core region increases with increasing value of yield stress at any time. The velocity and flow rate decrease, whereas wall shear stress and longitudinal impedance increase for increasing value of yield stress with other parameters held fixed. On the other hand, the velocity, flow rate and wall shear stress decrease but resistance to flow increases as the catheter radius ratio (ratio of catheter radius to vessel radius) increases with other parameters fixed. The results for power law fluid, Newtonian fluid and Bingham fluid are obtained as special cases from this model.     

The direct electromagnetic scattering problem by a mixed impedance screen in chiral media

In this article solvability results for the direct electromagnetic scattering problem for a mixed perfectly conducting-impedance screen in a chiral environment is studied. In particular, incident time-harmonic electromagnetic waves in a chiral medium upon a partially coated open surface Γ (the ‘screen’), that satisfies an impedance boundary condition on one side and a perfectly conducting boundary condition on the other side, are considered. We introduce the Beltrami fields, appropriate boundary integral relations for these fields are proved and via them a uniqueness result is established. A variational method in a suitable functional space setting is considered and using a Calderon type operator for the chiral case, existence for the scattering problem is established.     

利用远场模式的不完全数据反演声波阻尼系数

从阻尼边界条件声波散射问题的散射场远场模式的部分数据信息出发给出了反演声波阻尼系数的一种新方法,该问题既是非线性的又是不适定的,这里利用Tikhonov正则化方法将问题转化为一个最优化问题,成功地处理了第一类算子方程的不适定性及该问题的非线性性,给出了具体的数值方法并对其收敛性进行了严格地证明,数值结果表明该方法是非常准确且简单易行的．     

Cloaking via impedance boundary condition for the 2-D Helmholtz equation

We consider control problems for the 2-D Helmholtz equation in an unbounded domain with partially coated boundary. Dirichlet boundary condition is given on one part of the boundary and the impedance boundary condition is imposed on another its part. The role of control in control problem under study is played by boundary impedance. Quadratic tracking–type functionals for the field play the role of cost functionals. Solvability of control problems is proved. The uniqueness and stability of optimal solutions with respect to certain perturbations of both cost functional and incident field are established.     

Generalized impedance boundary condition at high frequency for a domain with thin layer: the circular case

Consider a conducting disk surrounded by a thin dielectric layer submitted to an electric field at the pulsation ω. The conductivity of the layer grows like ω^1?γ, γ∈[0,1], when the pulsation ω?tends to infinity. Using a pseudodifferential approach on the torus, we build an equivalent boundary condition with the help of an appropriate factorization of Helmholtz operator in the layer. This generalized impedance condition approximates the thin membrane in the high frequency limit for small thickness of the layer. L ²-error estimates are given and we illustrate our results with numerical simulations. This work extends, in the circular geometry, previous works of Lafitte and Lebeau (Lafitte O. Lebeau G. 1993, Équations de Maxwell et opérateur d’impédance sur le bord d’un obstacle convexe absorbant. Comptes Rendus de l ' Académic dis Science, Paris, Série I, Mathématiques, 316(11), 1177–1182); (Lafitte O.D., 1999, Diffraction in the high frequency regime by a thin layer of dielectric material. I. The equivalent impedance boundary condition. SIAM Journal on Applied Mathematics, 59(3), 1028–1052 (electronic)) in which γ?identically equals zero.     

Dynamic response of cable-stayed bridges to moving vehicles using the structural impedance method

The importance of dynamic interactions between cable-stayed bridges and heavy moving vehicles, such as trucks and locomotives, has been recognised by bridge engineers for a long time. A structural impedance algorithm is developed for analysing the dynamic response of cable-stayed bridges subjected to traversing vehicles. The bridge deck is modelled as an elastic plate, and the cables are idealised as springs for simplicity. The vehicles are modelled as a series of masses with suspension systems moving with different speeds and accelerations. A comprehensive computer program, CABLESIM, is developed for the static and dynamic analyses of a cable-stayed bridge. The accuracy of the numerical procedure and its computer implementation is verified with the available analytical and experimental results. A parametric study is conducted to investigate the effects of vehicle velocity, girder depth, different types of cable arrangements, and traffic load on the dynamic response of the deck. The numerical results are expected to be important in assessing the dynamics of cable-stayed bridge components and in determining the safety and allowable traffic conditions.     

Reconstruction and uniqueness of an inverse scattering problem with impedance boundary

Consider an inverse scattering problem for a scatterer D ⊂ R³ with impedance type boundary condition. By defining the scatterer shape in a completely new way, we give a constructive method to recover the scatterer shape with unknown impedance coefficient. The uniqueness for this inverse problem is also obtained.     

The asymptotic behaviour of weak solutions to the forward problem of electrical impedance tomography on unbounded three‐dimensional domains

The forward problem of electrical impedance tomography on unbounded domains can be studied by introducing appropriate function spaces for this setting. In this paper we derive the point‐wise asymptotic behaviour of weak solutions to this problem in the three‐dimensional case. Copyright © 2008 John Wiley & Sons, Ltd.     

Passive impedance systems with losses of scattering channels

We consider a new model of passive impedance system with minimum losses of scattering channels and bilaterally stable evolutionary semigroup. In the case of discrete time, the passive linear stationary bilaterally stable impedance system Σ is regarded as a part of a certain minimal lossless transmission system with a -unitary system operator and a bilaterally (J ₁, J ₂)-inner (in a certain weak sense) transfer function in the unit disk whose 22-block coincides with the impedance matrix of the system Σ, belongs to the Carathéodory class, and possesses a pseudoextension. If the external space of the system Σ is infinite-dimensional, then, instead of the last property, we consider more complicated necessary and sufficient conditions for the impedance matrix of the system Σ. We study passive bilaterally stable impedance realizations with minimum losses of scattering channels (minimum, optimal, *-optimal, minimum and optimal, and minimum and *-optimal). __________ Translated from Ukrains’kyi Matematychnyi Zhurnal, Vol. 59, No. 5, pp. 618–649, May, 2007.     

The impedance boundary value problem for the time harmonic Maxwell equations

The existence of a unique solution to Maxwell's equations defined in an exterior domain with impedance boundary condition is established for all frequencies. This is accomplished by reducing this problem to that of solving a system of singular integral equations and then regularizing this system such that the Riesz theory is applicable. We also consider the inverse problem in which it is desired to determine the impedance from a knowledge of the far field pattern. By restricting the impedance to lie a priori in a compact set results are obtained on the existence, uniqueness, and stability of the solution to this inverse scattering problem.     

Solution of the inverse spectral problem for an impedance with integrable derivative part II

The inverse eigenvalue problem for a Sturm-Liouville equation in impedance form with Dirichlet boundary conditions on a unit interval is solved. The solution of this nonlinear problem requires the investigation of a combined map from the logarithmic derivative of the impedance, assumed to be square integrable, to two sequences of spectral data. The first is the sequence of Dirichlet eigenvalues, shown to be locally bounded with square roots differing from the sequence of integral multiples of π by a square summable sequence. The second sequence has as first term the mean of the logarithmic derivative of the impedance. Each term in the remainder of the sequence is the logarithm of the product of the impedance and the derivative of an eigenfunction evaluated at an endpoint. It is shown that this second sequence is a locally bounded square summable sequence. The combined map is a real analytic isomorphism. A uniqueness result follows from the asymptotics and analyticity shown in the first paper and properties of Riesz bases. The existence is shown by giving a closed form of the solution. It is also shown that the same map solves the inverse Sturm-Liouville problem in impedance form with Neumann boundary conditions. © 1993 John Wiley & Sons, Inc.     

Surface wave scattering by a rectangular impedance cylinder located on a reactive plane

The electromagnetic scattering of the surface wave by a rectangular impedance cylinder located on an infinite reactive plane is considered for the case that the impedances of the horizontal and vertical sides of the cylinder can have different values. Firstly, the diffraction problem is reduced into a modified Wiener–Hopf equation of the third kind and then solved approximately. The solution contains branch‐cut integrals and two infinite sets of constants satisfying two infinite systems of linear algebraic equations. The approximate analytical or numerical evaluations of corresponding integrals and numerical solution of the linear algebraic equation systems are obtained for various values of parameters such as the surface reactance of the plane, the vertical and horizontal wall impedances, the width and the height of the cylinder. Copyright © 2004 John Wiley & Sons, Ltd.     

利用散射波远场模第P个傅立叶系数重建声波阻尼系数

阐述了利用声波散射远场模Fourier展开的第P个傅立叶系数（声散射远场模的不完全信息）,重建声阻抗系数的一种非线性最优化方法．并给出了该方法收敛性的证明,其数值例子说明这种方法的有效性和可行性．