Reconstructing small perturbations of scatterers from electric or acoustic far‐field measurements

In this paper, we consider the problem of determining the boundary perturbations of an object from far‐field electric or acoustic measurements. Assuming that the unknown scatterer boundary is a small perturbation of a circle, we develop a linearized relationship between the far‐field data and the shape of the object. This relationship is used to find the Fourier coefficients of the perturbation of the shape. Copyright © 2008 John Wiley & Sons, Ltd.     

Domain sensitivity analysis of the acoustic far‐field pattern

We consider acoustic scattering problems described by the mixed boundary value problem for the scalar Helmholtz equation in the exterior of a 2D bounded domain or in the exterior of a crack. The boundary of the domain is assumed to have a finite set of corner points where the scattered wave may have singular behaviour. The paper is concerned with the sensitivity of the far‐field pattern with respect to small perturbations of the shape of the scatterer. Using a modification of the method of adjoint problems, we obtain an integral representation for the Gâteaux derivative which contains only boundary values of functions easily computable by standard BEM and which depends explicitly on the perturbation of the boundary. In some cases, we show the direct influence of the singularities of the solution on the sensitivity of the far‐field pattern. In this way, we generalize the domain sensitivity analysis developed earlier for smooth domains by Hettlich, Kirsch, Kress, Potthast and others. Finally, we show that the same approach can be applied to scattering from 3D domains with smooth edges. Copyright © 2002 John Wiley & Sons, Ltd.     

On the Green function in visco‐elastic media obeying a frequency power‐law

In this work, we present an explicit expression for the Green function in a visco‐elastic medium. We choose Szabo and Wu's frequency power law model to describe the visco‐elastic properties and derive a generalized visco‐elastic wave equation. We express the ideal Green function (without any viscous effect) in terms of the viscous Green function using an attenuation operator. By means of an approximation of the ideal Green function, we address the problem of reconstructing a small anomaly in a visco‐elastic medium from wavefield measurements. Copyright © 2011 John Wiley & Sons, Ltd.     

Operator splitting for numerical solutions of the RLW Equation

In this study, the numerical behavior of the one-dimensional Regularized Long Wave (RLW) equation has been sought by the Strang splitting technique with respect to time. For this purpose, cubic B-spline functions are used with the finite element collocation method. Then, single solitary wave motion, the interaction of two solitary waves and undular bore problems have been studied and the effectiveness of the method has been investigated. The new results have been compared with those of some of the previous studies available in the literature. The stability analysis has also been taken into account by the von Neumann method.     

Finite Element Method and Discontinuous Galerkin Method for Stochastic Scattering Problem of Helmholtz Type in ℝ<Superscript><Emphasis Type="Italic">d</Emphasis></Superscript> (<Emphasis Type="Italic">d</Emphasis> = 2, 3)

In this paper, we address the finite element method and discontinuous Galerkin method for the stochastic scattering problem of Helmholtz type in ℝ ^d (d = 2, 3). Convergence analysis and error estimates are presented for the numerical solutions. The effects of the noises on the accuracy of the approximations are illustrated. Results of the numerical experiments are provided to examine our theoretical results. The first author is supported by NSF under grand number 0609918 and AFOSR under grant numbers FA9550-06-1-0234 and FA9550-07-1-0154, the second author is supported by NSFC(10671082, 10626026, 10471054), and the third author is supported by NNSF (No. 10701039 of China) and 985 program of Jilin University.     

Spatial behaviour of solutions of the dual‐phase‐lag heat equation

In this paper we study the spatial behaviour of solutions of some problems for the dual‐phase‐lag heat equation on a semi‐infinite cylinder. The theory of dual‐phase‐lag heat conduction leads to a hyperbolic partial differential equation with a third derivative with respect to time. First, we investigate the spatial evolution of solutions of an initial boundary‐value problem with zero boundary conditions on the lateral surface of the cylinder. Under a boundedness restriction on the initial data, an energy estimate is obtained. An upper bound for the amplitude term in this estimate in terms of the initial and boundary data is also established. For the case of zero initial conditions, a more explicit estimate is obtained which shows that solutions decay exponentially along certain spatial‐time lines. A class of non‐standard problems is also considered for which the temperature and its first two time derivatives at a fixed time T are assumed proportional to their initial values. Copyright © 2004 John Wiley & Sons, Ltd.     

Reconstruction of the shape of object with near field measurements in a half-plane

We consider a mathematical problem modelling some characteristics of near field optical microscope.We take a monofrequency line source to illuminate a sample with constant index of refraction and use the scattered field data measured near the sample to reconstruct the shape of it. Mixed reciprocity relation and factorization method are applied to solve our problem.Some numerical examples to show the feasibility of the method are presented.     

On the numerical solution of Korteweg-de Vries equation by the iterative splitting method

In this paper, we apply the method of iterative operator splitting on the Korteweg-de Vries (KdV) equation. The method is based on first, splitting the complex problem into simpler sub-problems. Then each sub-equation is combined with iterative schemes and solved with suitable integrators. Von Neumann analysis is performed to achieve stability criteria for the proposed method applied to the KdV equation. The numerical results obtained by iterative splitting method for various initial conditions are compared with the exact solutions. It is seen that they are in a good agreement with each other.     

Local existence of solutions of a three phase‐field model for solidification

In this article we discuss the local existence and uniqueness of solutions of a system of parabolic differential partial equations modeling the process of solidification/melting of a certain kind of alloy. This model governs the evolution of the temperature field, as well as the evolution of three phase‐field functions; the first two describe two different possible solid crystallization states and the last one describes the liquid state. Copyright © 2008 John Wiley & Sons, Ltd.     

Collocation discretization for an integral equation in ocean acoustics with depth‐dependent speed of sound

We analyze two collocation schemes for the Helmholtz equation with depth‐dependent sonic wave velocity, modeling time‐harmonic acoustic wave propagation in a three‐dimensional inhomogeneous ocean of finite height. Both discretization schemes are derived from a periodized version of the Lippmann‐Schwinger integral equation that equivalently describes the sound wave. The eigenfunctions of the corresponding periodized integral operator consist of trigonometric polynomials in the horizontal variables and eigenfunctions to some Sturm‐Liouville operator linked to the background profile of the sonic wave velocity in the vertical variable. Applying an interpolation projection onto a space spanned by finitely many of these eigenfunctions to either the unknown periodized wave field or the integral operator yields two different collocation schemes. A convergence estimate of Sloan [J. Approx. Theory, 39:97–117, 1983] on non‐polynomial interpolation allows to show converge of both schemes, together with algebraic convergence rates depending on the smoothness of the inhomogeneity and the source. Copyright © 2016 John Wiley & Sons, Ltd.     

双连通区域上的电磁波散射问题数值解

对于双连通区域上的电磁波散射问题,通过位势理论将其转化为边界积分方程组问题,然后采用Nystrom法和配置法对其离散求解,针对不同形状的障碍散射体,给出远场模式的数值解.     

Consistency of iterative operator‐splitting methods: Theory and applications

In this article, we describe a different operator‐splitting method for decoupling complex equations with multidimensional and multiphysical processes for applications for porous media and phase‐transitions. We introduce different operator‐splitting methods with respect to their usability and applicability in computer codes. The error‐analysis for the iterative operator‐splitting methods is discussed. Numerical examples are presented. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

New kink solutions for the van der Waals p‐system

The simple equation method and modified simple equation method are employed to seek exact traveling wave solutions to the (1 + 1)‐dimensional van der Waals gas system in the viscosity‐capillarity regularization form. Under the help of Mathematica, new classes of kink solutions are derived. Numerical simulations with special choices of the free parameters are displayed by three‐ and two‐dimensional plots. The two methods demonstrate simplicity, reliability, and efficiency.     

Time‐splitting methods with charge conservation for the nonlinear Dirac equation

In this work, four numerical time‐splitting methods are proposed for the (1 + 1)‐dimensional nonlinear Dirac equation. All of these methods (or schemes) are proved to satisfy the charge conservation in the discrete level. To enhance the computation efficiency, the block Thomas algorithm is adopted. Numerical experiments are given to test the accuracy order for these schemes, to simulate numerically the binary collision including two standing waves and two moving solitons, meanwhile, the dynamic properties for the nonlinear Dirac equation are discussed. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1582–1602, 2017     

Iterative operator‐splitting methods for nonlinear differential equations and applications

In this article, we consider iterative operator‐splitting methods for nonlinear differential equations with bounded and unbounded operators. The main feature of the proposed idea is the embedding of Newton's method for solving the split parts of the nonlinear equation at each step. The convergence properties of such a mixed method are studied and demonstrated. We confirm with numerical applications the effectiveness of the proposed scheme in comparison with the standard operator‐splitting methods by providing improved results and convergence rates. We apply our results to deposition processes. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1026–1054, 2011     

Application of the homotopy perturbation method to the modified regularized long‐wave equation

In this article, the homotopy perturbation method (HPM) is used to implement the modified regularized long‐wave (MRLW) equation with some initial conditions. The method is very efficient and convenient and can be applied to a large class of problems. In the last, three invariants of the motion are evaluated to determine the conservation properties of the system. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

On the convergence rate and optimization of a numerical method with splitting of boundary conditions for the stokes system in a spherical layer in the axisymmetric case: Modification for thick layers

The convergence rate of a fast-converging second-order accurate iterative method with splitting of boundary conditions constructed by the authors for solving an axisymmetric Dirichlet boundary value problem for the Stokes system in a spherical gap is studied numerically. For R/r exceeding about 30, where r and R are the radii of the inner and outer boundary spheres, it is established that the convergence rate of the method is lower (and considerably lower for large R/r) than the convergence rate of its differential version. For this reason, a really simpler, more slowly converging modification of the original method is constructed on the differential level and a finite-element implementation of this modification is built. Numerical experiments have revealed that this modification has the same convergence rate as its differential counterpart for R/r of up to 5 × 10³. When the multigrid method is used to solve the split and auxiliary boundary value problems arising at iterations, the modification is more efficient than the original method starting from R/r ～ 30 and is considerably more efficient for large values of R/r. It is also established that the convergence rates of both methods depend little on the stretching coefficient η of circularly rectangular mesh cells in a range of η that is well sufficient for effective use of the multigrid method for arbitrary values of R/r smaller than ～ 5 × 10³.     

Mixed boundary value problems of diffraction by a half‐plane with an obstacle perpendicular to the boundary

The paper is devoted to the analysis of wave diffraction problems modeled by classes of mixed boundary conditions and the Helmholtz equation, within a half‐plane with a crack. Potential theory together with Fredholm theory, and explicit operator relations, are conveniently implemented to perform the analysis of the problems. In particular, an interplay between Wiener–Hopf plus/minus Hankel operators and Wiener–Hopf operators assumes a relevant preponderance in the final results. As main conclusions, this study reveals conditions for the well‐posedness of the corresponding boundary value problems in certain Sobolev spaces and equivalent reduction to systems of Wiener–Hopf equations. Copyright © 2013 John Wiley & Sons, Ltd.     

Determination of the Robin coefficient in a nonlinear boundary condition for a steady‐state problem

In this paper, a constant heat transfer coefficient present in a nonlinear Robin‐type boundary condition associated with an elliptic equation is reconstructed uniquely from a single boundary energy measurement. Two types of such boundary energy measurement are considered, and solvability theorems for the solution of the resulting nonlinear inverse problems are provided. Further, one‐dimensional numerical results are presented and discussed. Copyright © 2008 John Wiley & Sons, Ltd.     

Analysis of the operator splitting scheme for the Cahn‐Hilliard equation with a viscosity term

In this paper, we consider a second‐order fast explicit operator splitting method for the viscous Cahn‐Hilliard equation, which includes a viscosity term αΔu_t (α ∈ (0, 1)) described the influences of internal micro‐forces. The choice α = 0 corresponds to the classical Cahn‐Hilliard equation whilst the choice α = 1 recovers the nonlocal Allen‐Cahn equation. The fundamental idea of our method is to split the original problem into linear and nonlinear parts. The linear subproblem is numerically solved using a pseudo‐spectral method, and thus an ordinary differential equation is obtained. The nonlinear one is solved via TVD‐RK method. The stability and convergence are discussed in L²‐norm. Numerical experiments are performed to validate the accuracy and efficiency of the proposed method. Besides, a detailed comparison is made for the dynamics and the coarsening process of the metastable pattern for various values of α. Moreover, energy degradation and mass conservation are also verified.