Optimal finite difference grids and rational approximations of the square root I. Elliptic problems

The main objective of this paper is optimization of second‐order finite difference schemes for elliptic equations, in particular, for equations with singular solutions and exterior problems. A model problem corresponding to the Laplace equation on a semi‐infinite strip is considered. The boundary impedance (Neumann‐to‐Dirichlet map) is computed as the square root of an operator using the standard three‐point finite difference scheme with optimally chosen variable steps. The finite difference approximation of the boundary impedance for data of given smoothness is the problem of rational approximation of the square root on the operator's spectrum. We have implemented Zolotarev's optimal rational approx‐imant obtained in terms of elliptic functions. We have also found that a geometrical progression of the grid steps with optimally chosen parameters is almost as good as the optimal approximant. For bounded operators it increases from second to exponential the convergence order of the finite difference impedance with the convergence rate proportional to the inverse of the logarithm of the condition number. For the case of unbounded operators in Sobolev spaces associated with elliptic equations, the error decays as the exponential of the square root of the mesh dimension. As an example, we numerically compute the Green function on the boundary for the Laplace equation. Some features of the optimal grid obtained for the Laplace equation remain valid for more general elliptic problems with variable coefficients. © 2000 John Wiley & Sons, Inc.     

Integral equation methods for scattering by infinite rough surfaces

In this paper, we consider the Dirichlet and impedance boundary value problems for the Helmholtz equation in a non‐locally perturbed half‐plane. These boundary value problems arise in a study of time‐harmonic acoustic scattering of an incident field by a sound‐soft, infinite rough surface where the total field vanishes (the Dirichlet problem) or by an infinite, impedance rough surface where the total field satisfies a homogeneous impedance condition (the impedance problem). We propose a new boundary integral equation formulation for the Dirichlet problem, utilizing a combined double‐ and single‐layer potential and a Dirichlet half‐plane Green's function. For the impedance problem we propose two boundary integral equation formulations, both using a half‐plane impedance Green's function, the first derived from Green's representation theorem, and the second arising from seeking the solution as a single‐layer potential. We show that all the integral equations proposed are uniquely solvable in the space of bounded and continuous functions for all wavenumbers. As an important corollary we prove that, for a variety of incident fields including an incident plane wave, the impedance boundary value problem for the scattered field has a unique solution under certain constraints on the boundary impedance. Copyright © 2003 John Wiley & Sons, Ltd.     

An Inverse Problem of Calderón Type with Partial Data

A generalized variant of the Calderón problem from electrical impedance tomography with partial data for anisotropic Lipschitz conductivities is considered in an arbitrary space dimension n ≥ 2. The following two results are shown: (i) The selfadjoint Dirichlet operator associated with an elliptic differential expression on a bounded Lipschitz domain is determined uniquely up to unitary equivalence by the knowledge of the Dirichlet-to-Neumann map on an open subset of the boundary, and (ii) the Dirichlet operator can be reconstructed from the residuals of the Dirichlet-to-Neumann map on this subset.     

Detecting inclusions with a generalized impedance condition from electrostatic data via sampling

In this paper, we derive a sampling method to solve the inverse shape problem of recovering an inclusion with a generalized impedance condition from electrostatic Cauchy data. The generalized impedance condition is a second order differential operator applied to the boundary of the inclusion. We assume that the Dirichlet‐to‐Neumann mapping is given from measuring the current on the outer boundary from an imposed voltage. A simple numerical example is given to show the effectiveness of the proposed inversion method for recovering the inclusion. We also consider the inverse impedance problem of determining the impedance parameters for a known material from the Dirichlet‐to‐Neumann mapping assuming the inclusion has been reconstructed where uniqueness for the reconstruction of the coefficients is proven.     

Iterative method for solving a three-dimensional electrical impedance tomography problem in the case of piecewise constant conductivity and one measurement on the boundary

The problem of electrical impedance tomography in a bounded three-dimensional domain with a piecewise constant electrical conductivity is considered. The boundary of the inhomogeneity is assumed to be unknown. The inverse problem is to determine the surface that is the boundary of the inhomogeneity from given measurements of the potential and its normal derivative on the outer boundary of the domain. An iterative method for solving the inverse problem is proposed, and numerical results are presented.     

Asymptotic expansions for the Sturm-Liouville problem by homotopy perturbation method

The asymptotic formulae for the eigenvalues and eigenfunctions of Sturm-Liouville problem with the Dirichlet boundary conditions when the potential is square integrable on [0, 1] are obtained by using homotopy perturbation method.     

Uniqueness For An Inverse Boundary Value Problem For Dirac Operators

The uniqueness of both the inverse boundary value problem and inverse scattering problem for Dirac equation with a magnetic potential and an electrical potential are proved. Also, a relation between the Dirichlet to Dirichlet map for the inverse boundary value problem and the scattering amplitude for the inverse scattering problem is given     

The factorization method for a class of inverse elliptic problems

In this paper the factorization method from inverse scattering theory and impedance tomography is extended to a class of general elliptic differential equations in divergence form. The inverse problem is to determine the interface ?Ω of an interior change of the material parameters from the Neumann‐Dirichlet map. Since absorption is allowed a suitable combination of the real and imaginary part of the Neumann‐Dirichlet map is needed to explicitely characterize Ω by the data. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nagumo type existence results of Sturm-Liouville BVP for impulsive differential equations

In this paper, the method of upper and lower solutions and the Schauder degree theory are employed in the study of Sturm-Liouville boundary value problems for second order impulsive differential equations. We obtain the existence of at least three solutions to the problem under the assumption that the nonlinear term f satisfies a Nagumo condition with respect to the first order derivative.     

Necessary and sufficient conditions for the solvability of the inverse problem for the matrix Sturm-Liouville operator

The matrix Sturm-Liouville operator on a finite interval with Dirichlet boundary conditions is studied. Properties of its spectral characteristics and the inverse problem of recovering the operator from these characteristics are investigated. Necessary and sufficient conditions on the spectral data of the operator are obtained. Research is conducted in the general case, with no a priori restrictions on the spectrum. A constructive algorithm for solving the inverse problem is provided.     

About stability of difference equations with continuous time and fading stochastic perturbations

The long term behavior of solutions of difference equations with continuous time and fading stochastic perturbations is investigated. It is shown that if the level of stochastic perturbations fades on the infinity, for instance, if it is given by square summable sequence, then an asymptotically stable and a square summable solution of a deterministic difference equation remains to be an asymptotically mean square stable and a mean square summable solution by stochastic perturbations.     

The Dirichlet problem for the minimal hypersurface equation on arbitrary domains of a Riemannian manifold

We show that the Dirichlet problem for the minimal hypersurface equation defined on arbitrary C ² bounded domain Ω of an arbitrary complete Riemannian manifold M is solvable if the oscillation of the boundary data is bounded by a function \({\mathcal{C}}\) that is explicitely given and that depends only on the first and second derivatives of the boundary data as well as the second fundamental form of the boundary \({\partial\Omega}\) and the Ricci curvature of the ambient space M. This result extends Theorem 2 of Jenkins-Serrin (J Reine Angew Math 229:170–187,1968) about the solvability of the Dirichlet problem for the minimal hypersurface equation defined on bounded domains of the Euclidean space. We deduce that the Dirichlet problem for the minimal hypersurface equation is solvable for any continuous boundary data on a mean convex domain. We also show existence and uniqueness of the Dirichlet problem with boundary data at infinity—exterior Dirichlet problem—on Hadamard manifolds.     

A stabilized implicit fractional-step method for the time-dependent Navier–Stokes equations using equal-order pairs

A stabilized implicit fractional-step method for numerical solutions of the time-dependent Navier–Stokes equations is presented in this paper. The time advancement is decomposed into a sequence of two steps: the first step has the structure of the linear elliptic problem; the second step can be seen as the generalized Stokes problem. The two problems satisfy the full homogeneous Dirichlet boundary conditions on the velocity. On the other hand, a locally stabilized term is added in the second step of the schemes. It allows one to enhance the numerical stability and efficiency by using the equal-order pairs. Convergence analysis and error estimates for the velocity and pressure of the schemes are established via the energy method. Some numerical experiments are also used to demonstrate the efficiency of this new method.     

n维耦合振子格点系统的同步性

本文通过在相空间中引入新范数的方法研究了Dirichlet,Neumann,周期三个不同边界条件下带有周期外力的n维二阶耦合振子格点系统的同步性.在Dirichlet边界条件下,如果系统非线性项的一阶偏导数有界,则当耦合系数足够大时,系统是有界耗散的并且任意两个解之间是同步的.在Neumann与周期边界条件下,如果不同子系统的外力之间的差和不同子系统的非线性项之间的差都比较小,并且系统是有界耗散的,则当耦合系数足够大时,系统任意一个解的任意两个分量之间是渐近同步的.这两种情况下,当耦合系数c_1→+∞,c_2→+∞时,系统任意一个解的任意两个分量之间是同步的.     

Inverse boundary value problem for ocean acoustics

For an ocean with constant depth and rigid bottom which contains compactly supported inhomogeneity of the water sound velocity, we prove uniqueness for the identification of the inhomogeneity from the Dirichlet‐to‐Neumann (DtN) map on the surface of a bounded region containing the inhomogeneity. The DtN map is the map which maps the pressure applied on the boundary of this region to the corresponding flux (displacement). In an analogous geometric configuration and with similar boundary conditions, the uniqueness for the inverse electroconductivity problem from the DtN map (i.e. voltage‐to‐current map) can be proved in the same framework. Copyright © 2001 John Wiley & Sons, Ltd.     

非线性项含一阶导数广义Sturm-Liouville边值问题的多个正解

By using fixed point theorems, we consider multiplicity of positive solutions for second-order generalized Sturm-Liouville boundary value problem, where the first order derivative is involved in the nonlinear term explicitly. We show the existence of multiple positive solutions for the problems. Example is given to illustrate the main results of the article.     

Stability of inverse problems in an infinite slab with partial data

In this paper, we study the stability of two inverse boundary value problems in an infinite slab with partial data. These problems have been studied by Li and Uhlmann for the case of the Schrödinger equation and by Krupchyk, Lassas, and Uhlmann for the case of the magnetic Schrödinger equation. Here, we quantify the method of uniqueness proposed by Li and Uhlmann and prove a log–log stability estimate for the inverse problems associated to the Schrödinger equation. The boundary measurements considered in these problems are modeled by partial knowledge of the Dirichlet-to-Neumann map: in the first inverse problem, the corresponding Dirichlet and Neumann data are known on different boundary hyperplanes of the slab; in the second inverse problem, they are known on the same boundary hyperplane of the slab.     

On the rate of convergence for perforated plates with a small interior Dirichlet zone

The aim of the paper is to compare the asymptotic behavior of solutions of two boundary value problems for an elliptic equation posed in a thin periodically perforated plate. In the first problem, we impose homogeneous Dirichlet boundary condition only at the exterior lateral boundary of the plate, while at the remaining part of the boundary Neumann condition is assigned. In the second problem, Dirichlet condition is also imposed at the surface of one of the holes. Although in these two cases, the homogenized problem is the same, the asymptotic behavior of solutions is rather different. In particular, the presence of perturbation in the boundary condition in the second problem results in logarithmic rate of convergence, while for non-perturbed problem the rate of convergence is of power-law type.     

一类带裂缝散射问题的边界积分方程方法

声波的散射问题中,如果散射体由不可穿透障碍物和可穿透裂缝两部分组成,障碍物表面分别满足第一类和第三类边界条件,裂缝两边满足不同的第二类边界条件,通过位势理论,可以将此混合问题转化为边界积分方程,通过Fredholm算子理论可以得到这个边界积分方程解的存在性和唯一性,从而获得原问题解的存在和唯一性.     

On the rate of convergence for perforated plates with a small interior Dirichlet zone

The aim of the paper is to compare the asymptotic behavior of solutions of two boundary value problems for an elliptic equation posed in a thin periodically perforated plate. In the first problem, we impose homogeneous Dirichlet boundary condition only at the exterior lateral boundary of the plate, while at the remaining part of the boundary Neumann condition is assigned. In the second problem, Dirichlet condition is also imposed at the surface of one of the holes. Although in these two cases, the homogenized problem is the same, the asymptotic behavior of solutions is rather different. In particular, the presence of perturbation in the boundary condition in the second problem results in logarithmic rate of convergence, while for non-perturbed problem the rate of convergence is of power-law type.