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.     

Analysis of a meshless method for the time fractional diffusion-wave equation

In this paper a numerical technique is proposed for solving the time fractional diffusion-wave equation. We obtain a time discrete scheme based on finite difference formula. Then, we prove that the time discrete scheme is unconditionally stable and convergent using the energy method and the convergence order of the time discrete scheme is \(\mathcal {O}(\tau ^{3-\alpha })\). Firstly, we change the main problem based on Dirichlet boundary condition to a new problem based on Robin boundary condition and then, we consider a semi-discrete scheme with Robin boundary condition and show when \(\beta \rightarrow +\infty \) solution of the main semi-discrete problem with Dirichlet boundary condition is convergent to the solution of the new semi-discrete problem with Robin boundary condition. We consider the new semi-discrete problem with Robin boundary condition and use the meshless Galerkin method to approximate the spatial derivatives. Finally, we obtain an error bound for the new problem. We prove that convergence order of the numerical scheme based on Galekin meshless is \(\mathcal {O}(h)\). In the considered method the appeared integrals are approximated using Gauss Legendre quadrature formula. The main aim of the current paper is to obtain an error estimate for the meshless Galerkin method based on the radial basis functions. Numerical examples confirm the efficiency and accuracy of the proposed scheme.     

四元数分析中的T算子与两类边值问题

本文研究四元数分析中的非齐次 Ｄｉｒａｃ方程．引入了这类方程的分布解即 Ｔ算子,证明了Ｔ算子的一些性质并考察了非齐次Ｄｉｒａｃ方程的Ｄｉｒｉｃｈｌｅｔ边值问题,并将结果推广到高阶非齐次Ｄｉｒａｃ方程及这种方程的一类边值问题的情况．     

On a Class of Overdetermined Eigenvalue Problems

In this paper we present some new results of symmetry for inhomogeneous Dirichlet eigenvalue problems overdetermined by a condition involving the gradient of the first eigenfunction on the boundary. One specificity of the problem studied is the dependence of the equation and the boundary condition on the distance to the origin. The method of investigation is based on the use of continuous Steiner symmetrization together with some domain derivative tools. An application is given to the study of an overdetermined eigenvalue problem for a wedge-like membrane. © 1997 B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

Singular limit of a competition–diffusion system with large interspecific interaction

We consider a competition–diffusion system for two competing species; the density of the first species satisfies a parabolic equation together with an inhomogeneous Dirichlet boundary condition whereas the second one either satisfies a parabolic equation with a homogeneous Neumann boundary condition, or an ordinary differential equation. Under the situation where the two species spatially segregate as the interspecific competition rate becomes large, we show that the resulting limit problem turns out to be a free boundary problem. We focus on the singular limit of the interspecific reaction term, which involves a measure located on the free boundary.     

Existence and uniqueness of solutions of hyperbolic equations in divergence form with various boundary conditions on various parts of the boundary

We prove the existence and uniqueness of an energy class solution of an initial–boundary value problem for a semilinear equation in divergence form. We consider the case in which an inhomogeneous third boundary condition is posed on one part of the lateral surface of the cylinder in which the equation is studied and the homogeneous Dirichlet boundary condition is posed on the other part of the lateral surface.     

To the theory of asymptotically stable second-order accurate two-stage scheme for an inhomogeneous parabolic initial-boundary value problem

An asymptotically stable two-stage difference scheme applied previously to a homogeneous parabolic equation with a homogeneous Dirichlet boundary condition and an inhomogeneous initial condition is extended to the case of an inhomogeneous parabolic equation with an inhomogeneous Dirichlet boundary condition. It is shown that, in the class of schemes with two stages (at every time step), this difference scheme is uniquely determined by ensuring that high-frequency spatial perturbations are fast damped with time and the scheme is second-order accurate and has a minimal error. Comparisons reveal that the two-stage scheme provides certain advantages over some widely used difference schemes. In the case of an inhomogeneous equation and a homogeneous boundary condition, it is shown that the extended scheme is second-order accurate in time (for individual harmonics). The possibility of achieving second-order accuracy in the case of an inhomogeneous Dirichlet condition is explored, specifically, by varying the boundary values at time grid nodes by O(τ ²), where τ is the time step. A somewhat worse error estimate is obtained for the one-dimensional heat equation with arbitrary sufficiently smooth boundary data, namely, $O\left( {\tau ^2 \ln \frac{T} {\tau }} \right) $ , where T is the length of the time interval.     

A.V. Bitsadze's observation on bianalytic functions and the Schwarz problem

According to an observation of A.V. Bitsadze from 1948 the Dirichlet problem for bianalytic functions is ill-posed. A natural boundary condition for the polyanalytic operator, however, is the Schwarz condition. An integral representation for the solutions in the unit disc to the inhomogeneous polyanalytic equation satisfying Schwarz boundary conditions is known. This representation is extended here to any simply connected plane domain having a harmonic Green function. Some other boundary value problems are investigated with some Dirichlet and Neumann conditions illuminating that just the Schwarz problem is a natural boundary condition for the Bitsadze operator.     

Some sufficient conditions for the non-negativity preservation property in the discrete heat conduction model

In the course of the numerical approximation of mathematical models there is often a need to solve a system of linear equations with a tridiagonal or a block-tridiagonal matrices. Usually it is efficient to solve these systems using a special algorithm (tridiagonal matrix algorithm or TDMA) which takes advantage of the structure. The main result of this work is to formulate a sufficient condition for the numerical method to preserve the non-negativity for the special algorithm for structured meshes. We show that a different condition can be obtained for such cases where there is no way to fulfill this condition. Moreover, as an example, the numerical solution of the two-dimensional heat conduction equation on a rectangular domain is investigated by applying Dirichlet boundary condition and Neumann boundary condition on different parts of the boundary of the domain. For space discretization, we apply the linear finite element method, and for time discretization, the well-known Θ-method. The theoretical results of the paper are verified by several numerical experiments.     

A Tikhonov regularization method for Cauchy problem based on a new relaxation model

In this paper, we consider a Cauchy problem of recovering both missing value and flux on inaccessible boundary from Dirichlet and Neumann data measured on the remaining accessible boundary. Associated with two mixed boundary value problems, a regularized Kohn-Vogelius formulation is proposed. With an introduction of a relaxation parameter, the Dirichlet boundary conditions are approximated by two Robin ones. Compared to the existing work, weaker regularity is required on the Dirichlet data. This makes the proposed model simpler and more efficient in computation. A series of theoretical results are established for the new reconstruction model. Several numerical examples are provided to show feasibility and effectiveness of the proposed method. For simplicity of the statements, we take Poisson equation as the governed equation. However, the proposed method can be applied directly to Cauchy problems governed by more general equations, even other linear or nonlinear inverse problems.     

Numerical Solution of the Positional Boundary Control Problem for the Wave Equation with Unknown Initial Data

The problem of one-sided boundary Neumann control is considered for the onedimensional wave equation. Information about the initial state of the process is absent. Instead, the values of Dirichlet observations are received in real time at the controlled boundary. The aim is to bring the process to a complete rest by means of positional boundary controls. To solve this problem, we propose an efficient numerical algorithm with an optimal guaranteed damping time. Some results of numerical experiments are presented.     

A fast interface solver for the biharmonic Dirichlet problem on polygonal domains

Summary. In this paper we propose and analyze an efficient discretization scheme for the boundary reduction of the biharmonic Dirichlet problem on convex polygonal domains. We show that the biharmonic Dirichlet problem can be reduced to the solution of a harmonic Dirichlet problem and of an equation with a Poincaré-Steklov operator acting between subspaces of the trace spaces. We then propose a mixed FE discretization (by linear elements) of this equation which admits efficient preconditioning and matrix compression resulting in the complexity . Here is the number of degrees of freedom on the underlying boundary, is an error reduction factor, or for rectangular or polygonal boundaries, respectively. As a consequence an asymptotically optimal iterative interface solver for boundary reductions of the biharmonic Dirichlet problem on convex polygonal domains is derived. A numerical example confirms the theory. Received September 1, 1995 / Revised version received February 12, 1996     

Construction of polynomial solutions to the Dirichlet problem for the polyharmonic equation in a ball

An algorithm is proposed for the analytical construction of a polynomial solution to Dirichlet problem for an inhomogeneous polyharmonic equation with a polynomial right-hand side and polynomial boundary data in the unit ball.     

To the Solution of the Hilbert Problem with Infinite Index

In this paper, we obtain a generalization of the method of regularizing multipliers for the solution of the Hilbert boundary-value problem with finite index in the theory of analytic functions to the case of an infinite power-behaved index. This method is used to obtain a general solution of the homogeneous Hilbert problem for the half-plane, a solution that depends on the existence and the number of entire functions possessing mirror symmetry with respect to the real axis and satisfying some additional constraints related to the singularity characteristic of the index. To solve of the inhomogeneous problem, we essentially use a specially constructed solution of the homogeneous problem whereby we reduce the boundary condition of the Hilbert problem to a Dirichlet problem.     

Semi-implicit finite volume level set method for advective motion of interfaces in normal direction

In this paper a semi-implicit finite volume method is proposed to solve the applications with moving interfaces using the approach of level set methods. The level set advection equation with a given speed in normal direction is solved by this method. Moreover, the scheme is used for the numerical solution of eikonal equation to compute the signed distance function and for the linear advection equation to compute the so-called extension speed [1]. In both equations an extrapolation near the interface is used in our method to treat Dirichlet boundary conditions on implicitly given interfaces. No restrictive CFL stability condition is required by the semi-implicit method that is very convenient especially when using the extrapolation approach. In summary, we can apply the method for the numerical solution of level set advection equation with the initial condition given by the signed distance function and with the advection velocity in normal direction given by the extension speed. Several advantages of the proposed approach can be shown for chosen examples and application. The advected numerical level set function approximates well the property of remaining the signed distance function during whole simulation time. Sufficiently accurate numerical results can be obtained even with the time steps violating the CFL stability condition.     

Numerical method for determining the inhomogeneity boundary in the Dirichlet problem for Laplace’s equation in a piecewise homogeneous medium

The Dirichlet problem for Laplace’s equation in a two-dimensional domain filled with a piecewise homogeneous medium is considered. The boundary of the inhomogeneity is assumed to be unknown. The inverse problem of determining the inhomogeneity boundary from additional information on the solution of the Dirichlet problem is considered. A numerical method based on the linearization of the nonlinear operator equation for the unknown boundary is proposed for solving the inverse problem. The results of numerical experiments are presented.     

The Mixed Problem for the Laplace Equation in an Exterior Domain with an Arbitrary Partition of the Boundary

In this paper we propose a method for solving the mixed boundary-value problem for the Laplace equation in a connected exterior domain with an arbitrary partition of the boundary. All simple closed curves making up the boundary are divided into three sets. On the elements of the first set the Dirichlet condition is given, on the elements of the second set the third boundary condition is prescribed, and the third set, in turn, is divided into two subsets of simple closed arcs, with the Dirichlet condition prescribed on the elements of one of these subsets and the third boundary condition on the elements of the other subset. The problem is reduced to a uniquely solvable Fredholm equation of the second kind in a Banach space. The third boundary-value problem and the mixed Dirichlet--Neumann problem are particular cases of the problem under study.     

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.     

L-stable Simpson’s 3/8 rule and Burgers’ equation

In this paper we develop an unconditionally stable third order time integration formula for the diffusion equation with Neumann boundary condition. We use a suitable arithmetic average approximation and explicit backward Euler formula and then develop a third order L-stable Simpson’s 3/8 type formula. We also observe that the arithmetic average approximation is not unique. Then L-stable Simpson’s 3/8 type formula and Hopf-Cole transformation is used to solve Burger’s equation with Dirichlet boundary condition. It is also observed that this numerical method deals efficiently in case of inconsistencies in initial and boundary conditions.     

Numerical methods for a fixed domain formulation of the glacier profile problem with alternative boundary conditions

In this paper we develop a set of numerical techniques for the simulation of the profile evolution of a valley glacier in the framework of isothermal shallow ice approximation models. The different mathematical formulations are given in terms of a highly nonlinear parabolic equation. A first nonlinearity comes from the free boundary problem associated with the unknown basal extension of the glacier region. This feature is treated using a fixed domain complementarity formulation which is solved numerically by a duality method. The nonlinear diffusive term is explicitly treated in the time marching scheme. A convection dominated problem arises, so a characteristic scheme is proposed for the time discretization, while piecewise linear finite elements are used for the spatial discretization. The presence of infinite slopes in polar regimes motivates an alternative formulation based on a prescribed flux boundary condition at the head of the glacier instead a homogeneous Dirichlet one. Finally, several numerical examples illustrate the performance of the proposed methods.