解洞穴逆散射问题的一种正则化Newton方法

We consider the inverse problem to determine the shape of a open cavity embedded in the infinite ground plane from knowledge of the far-field pattern of the scattering of TM polarization.For its approximate solution we propose a regularized Newton iteration scheme.For a foundation of Newton type methods we establish the Fréchet differentiability of solution to the scattering problem with respect to the boundary of the cavity.Some numerical examples of the feasibility of the method are presented.     

Regularized Newton iteration method for a penetrable cavity with internal measurements in inverse scattering problem

In this paper, we consider the inverse scattering problem of determining the shape of a cavity with a penetrable inhomogeneous medium of compact support from one source and a knowledge of measurements placed on a curve inside the cavity. First, the boundary value problem of the partial differential equations can be transformed into an equivalent system of nonlinear and ill-posed integral equations for the unknown boundary. Then, we apply the regularized Newton iterative method to reconstruct the boundary and prove the injectivity for the linearized system. Finally, we present some numerical examples to show the feasibility of our method.     

Huygens’ principle and iterative methods in inverse obstacle scattering

The inverse problem we consider in this paper is to determine the shape of an obstacle from the knowledge of the far field pattern for scattering of time-harmonic plane waves. In the case of scattering from a sound-soft obstacle, we will interpret Huygens’ principle as a system of two integral equations, named data and field equation, for the unknown boundary of the scatterer and the induced surface flux, i.e., the unknown normal derivative of the total field on the boundary. Reflecting the ill-posedness of the inverse obstacle scattering problem these integral equations are ill-posed. They are linear with respect to the unknown flux and nonlinear with respect to the unknown boundary and offer, in principle, three immediate possibilities for their iterative solution via linearization and regularization. In addition to presenting new results on injectivity and dense range for the linearized operators, the main purpose of this paper is to establish and illuminate relations between these three solution methods based on Huygens’ principle in inverse obstacle scattering. Furthermore, we will exhibit connections and differences to the traditional regularized Newton type iterations as applied to the boundary to far field map, including alternatives for the implementation of these Newton iterations.     

Inverse scattering from an open arc

A Newton method is presented for the approximate solution of the inverse problem to determine the shape of a sound-soft or perfectly conducting arc from a knowledge of the far-field pattern for the scattering of time-harmonic plane waves. Fréchet differentiability with respect to the boundary is shown for the far-field operator, which for a fixed incident wave maps the boundary arc onto the far-field pattern of the scattered wave. For the sake of completeness, the first part of the paper gives a short outline on the corresponding direct problem via an integral equation method including the numerical solution.     

An overfilled cavity problem for Maxwell's equations

This paper is concerned with the mathematical analysis of the scattering of a time‐harmonic electromagnetic plane wave by an open and overfilled cavity that is embedded in a perfect electrically conducting infinite ground plane, where the electromagnetic wave propagation is governed by the Maxwell equations. Above the flat ground surface and the open aperture of the cavity, the space is assumed to be filled with a homogeneous medium with a constant permittivity and permeability, whereas the interior of the cavity is filled with some inhomogeneous medium with a variable permittivity and permeability. The scattering problem is modeled as a boundary value problem over a bounded domain, with transparent boundary condition proposed on the hemisphere enclosing the inhomogeneity represented by the cavity. The existence and uniqueness of the weak solution for the model problem are established by using a variational approach. The perfectly matched layer (PML) method is investigated to truncate the unbounded electromagnetic cavity scattering problem. It is shown that the truncated PML problem attains a unique solution. An explicit error estimate is given between the solution of the original scattering problem and that of the truncated PML problem. The error estimate implies that the PML solution converges exponentially to the original cavity scattering problem by increasing either the PML medium parameter or the PML layer thickness. The convergence result is expected to be useful for determining the PML medium parameter in the computational electromagnetic scattering problem. Copyright © 2012 John Wiley & Sons, Ltd.     

A Newton iteration process for inverse eigenvalue problems

Summary We suppose an inverse eigenvalue problem which includes the classical additive and multiplicative inverse eigenvalue problems as special cases. For the numerical solution of this problem we propose a Newton iteration process and compare it with a known method. Finally we apply it to a numerical example.     

A posteriori parameter choice strategies for some Newton type methods for the regularization of nonlinearill-posed problems

This paper treats a class of Newton type methods for the approximate solution of nonlinear ill-posed operator equations, that use so-called filter functions for regularizing the linearized equation in each Newton step. For noisy data we derive an aposteriori stopping rule that yields convergence of the iterates to asolution, as the noise level goes to zero, under certain smoothness conditions on the nonlinear operator. Appropriate closeness and smoothness assumptions on the starting value and the solution are shown to lead to optimal convergence rates. Moreover we present an application of the Newton type methods under consideration to a parameter identification problem, together with some numerical results. Received November 29, 1996 / Revised version received April 25, 1997     

Shape reconstruction of an inverse Stokes problem

This paper deals with the shape reconstruction of a viscous incompressible fluid driven by the Stokes flow. For the approximate solution of the ill-posed and nonlinear problem we propose a regularized Newton method. A theoretical foundation for the Newton method is given by establishing the differentiability of the initial boundary value problem with respect to the interior boundary curve in the sense of a domain derivative. The numerical examples show that our theory is useful for practical purpose and the proposed algorithm is feasible.     

Newton Method for Solving the Multi-Variable Fuzzy Optimization Problem

In this article, we propose the Newton method to find a non-dominated solution of an unconstrained multi-variable fuzzy optimization problem. For this purpose, we use the Hukuhara differentiability of fuzzy-valued functions and partial order relation on set of fuzzy numbers.     

A hybrid method for sound-hard obstacle reconstruction

We are interested in solving the inverse problem of acoustic wave scattering to reconstruct the position and the shape of sound-hard obstacles from a given incident field and the corresponding far field pattern of the scattered field. The method we suggest is an extension of the hybrid method for the reconstruction of sound-soft cracks as presented in [R. Kress, P. Serranho, A hybrid method for two-dimensional crack reconstruction, Inverse Problems 21 (2005) 773–784] to the case of sound-hard obstacles. The designation of the method is justified by the fact that it can be interpreted as a hybrid between a regularized Newton method applied to a nonlinear operator equation with the operator that maps the unknown boundary onto the solution of the direct scattering problem and a decomposition method in the spirit of the potential method as described in [A. Kirsch, R. Kress, On an integral equation of the first kind in inverse acoustic scattering, in: Cannon, Hornung (Eds.), Inverse Problems, ISNM, vol. 77, 1986, pp. 93–102. Since the method does not require a forward solver for each Newton step its computational costs are reduced. By some numerical examples we illustrate the feasibility of the method.     

A hybrid method for inverse cavity scattering problem for shape

In this paper, we give a hybrid method to numerically solve the inverse open cavity scattering problem for cavity shape, given the scattered solution on the opening of the cavity. This method is a hybrid between an iterative method and an integral equations method for solving the Cauchy problem. The idea of this hybrid method is simple, the operation is easy, and the computation cost is small. Numerical experiments show the feasibility of this method, even for cases with noise.     

Shape reconstruction for the time‐dependent Navier‐Stokes flow

This article deals with the shape reconstruction of a bounded domain with a viscous incompressible fluid driven by the time‐dependent Navier‐Stokes equations. For the approximate solution of the ill‐posed and nonlinear problem we propose a regularized Newton method. A theoretical foundation for the Newton method is given by establishing the differentiability of the initial boundary value problem with respect to the interior boundary curve in the sense of the domain derivative. Numerical examples indicate the feasibility of our method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Evolution of phase transitions in methane hydrate

We consider a simplified model of methane hydrates which we cast as a nonlinear evolution problem. For its well-posedness we extend the existing theory to cover the case in which the problem involves a measurable family of graphs. We represent the nonlinearity as a subgradient and prove a useful comparison principle, thus optimal regularity results follow. For the numerical solution we apply a fully implicit scheme without regularization and use the semismooth Newton algorithm for a solver, and the graph is realized as a complementarity constraint (CC). The algorithm is very robust and we extend it to define an easy and superlinearly convergent fully implicit scheme for the Stefan problem and other multivalued examples.     

Aspects of electromagnetic scattering in chiral media

In this work, we study two operators that arise in electromagnetic scattering in chiral media. We first consider electromagnetic scattering by a chiral dielectric with a perfectly conducting core. We define a chiral Calderon‐type surface operator in order to solve the direct electromagnetic scattering problem. For this operator, we state coercivity and prove compactness properties. In order to prove existence and uniqueness of the problem, we define some other operators that are also related to the chiral Calderon‐type operator, and we state some of their properties that they and their linear combinations satisfy. Then we sketch how to use these operators in order to prove the existence of the solution of the direct scattering problem. Furthermore, we focus on the electromagnetic scattering problem by a perfect conductor in a chiral environment. For this problem, we study the chiral far‐field operator that is defined on a unit sphere and contains the far‐field data, and we state and prove some of its properties that are preliminaries properties for solving the inverse scattering problem. Copyright © 2016 John Wiley & Sons, Ltd.     

Shape inverse problem for the two‐dimensional unsteady Stokes flow

In this article, the shape inverse problem for the two‐dimensional unsteady Stokes flow has been presented. We employ Piola transformation to bypass the divergence free condition for the flow and prove the differentiability of the solution to the initial boundary value problem. For the approximate solution of the ill‐posed and nonlinear problem, we propose a regularized Gauss‐Newton method. The numerical examples show that our theory is useful for practical purpose and the proposed algorithm is feasible. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

The interior inverse scattering problem for cavities with an artificial obstacle

The interior inverse scattering by an impenetrable cavity is considered. Both the sources and the measurements are placed on a curve or surface inside the cavity. As a rule of thumb, both the direct and the inverse problems suffer from interior eigenvalues. The interior eigenvalues are removed by adding an artificial obstacle with impedance boundary condition to the underlying scattering system. For this new system, we prove a reciprocity relation for the scattered field and a uniqueness theorem for the inverse problem. Some new techniques are used in the arguments of the uniqueness proof because of the Lipschitz regularity of the boundary of the cavity. The linear sampling method is used for this new scattering system for reconstructing the shape of the cavity. Finally, some numerical experiments are presented to demonstrate the feasibility and effectiveness of the linear sampling method. In particular, the introduction of the artificial obstacle makes the linear sampling method robust to frequency.     

New open modified Newton Cotes type formulae as multilayer symplectic integrators

New modified open Newton Cotes integrators are introduced in this paper. For the new proposed integrators the connection between these new algorithms, differential methods and symplectic integrators is studied. Much research has been done on one step symplectic integrators and several of them have obtained based on symplectic geometry. However, the research on multistep symplectic integrators is very poor. Zhu et al. [1] studied the well known open Newton Cotes differential methods and they presented them as multilayer symplectic integrators. Chiou and Wu [2] studied the development of multistep symplectic integrators based on the open Newton Cotes integration methods. In this paper we introduce a new open modified numerical method of Newton Cotes type and we present it as symplectic multilayer structure. The new obtained symplectic schemes are applied for the solution of Hamilton’s equations of motion which are linear in position and momentum. An important remark is that the Hamiltonian energy of the system remains almost constant as integration proceeds. We have applied also efficiently the new proposed method to a nonlinear orbital problem and an almost periodic orbital problem.     

Numerical simulation for the shape reconstruction of a cavity

We consider the inverse problem of reconstructing the interior boundary curve of a cavity from the knowledge of the measurements on the exterior boundary. The domain derivative of the corresponding operator is presented, and this allows the investigation of the regularized Newton method for the solution of the ill‐posed and nonlinear problem. Numerical examples indicate the feasibility of our method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Numerical solution of a quasilinear parabolic problem

A combined approach of linearisation techniques and finite difference method is presented for obtaining the numerical solution of a quasilinear parabolic problem. The given problem is reduced to a sequence of linear problems by using the Picard or Newton methods. Each problem of this sequence is approximated by Crank-Nicolson difference scheme. The solutions of the resulting system of algebraic equations are obtained by using Block-Gaussian elimination method. Two numerical examples are solved by using both linearisation procedures to illustrate the method. For these examples, the Newton method is found to be more effective, especially when the given nonlinear problem has steep gradients.     

CONVERGENCE OF THE CRANK-NICOLSON/NEWTON SCHEME FOR NONLINEAR PARABOLIC PROBLEM

In this paper, the Crank-Nicolson/Newton scheme for solving numerically secondorder nonlinear parabolic problem is proposed. The standard Galerkin finite element method based on P2 conforming elements is used to the spatial discretization of the problem and the Crank-Nicolson/Newton scheme is applied to the time discretization of the resulted finite element equations. Moreover, assuming the appropriate regularity of the exact solution and the finite element solution, we obtain optimal error estimates of the fully discrete CrankNicolson/Newton scheme of nonlinear parabolic problem. Finally, numerical experiments are presented to show the efficient performance of the proposed scheme.