Finite difference on grids with nearly uniform cell area and line spacing for the wave equation on complex domains

A finite difference time-dependent numerical method for the wave equation, supported by recently derived novel elliptic grids, is analyzed. The method is successfully applied to single and multiple two-dimensional acoustic scattering problems including soft and hard obstacles with complexly shaped boundaries. The new grids have nearly uniform cell area (J-grids) and nearly uniform grid line spacing (αγ-grids). Numerical experiments reveal the positive impact of these two grid properties on the scattered field convergence to its harmonic steady state. The restriction imposed by stability conditions on the time step size is relaxed due to the near uniformity cell areas and grid line spacing. As a consequence, moderately large time steps can be used for relatively fine spatial grids resulting in greater accuracy at a lower computational cost. Also, numerical solutions for wave problems inside annular regions of complex shapes are obtained. The use of the new grids results in late time stability in contrast with other classical finite difference time-dependent methods.     

Finite difference reaction diffusion equations with nonlinear boundary conditions

This article is concerned with numerical solutions of finite difference systems of reaction diffusion equations with nonlinear internal and boundary reaction functions. The nonlinear reaction functions are of general form and the finite difference systems are for both time-dependent and steady-state problems. For each problem a unified system of nonlinear equations is treated by the method of upper and lower solutions and its associated monotone iterations. This method leads to a monotone iterative scheme for the computation of numerical solutions as well as an existence-comparison theorem for the corresponding finite difference system. Special attention is given to the dynamical property of the time-dependent solution in relation to the steady-state solutions. Application is given to a heat-conduction problem where a nonlinear radiation boundary condition obeying the Boltzmann law of cooling is considered. This application demonstrates a bifurcation property of two steady-state solutions, and determines the dynamic behavior of the time-dependent solution. Numerical results for the heat-conduction problem, including a test problem with known analytical solution, are presented to illustrate the various theoretical conclusions. © 1995 John Wiley & Sons, Inc.     

Time domain scattering and inverse scattering problems in a locally perturbed half-plane

This paper is concerned with efficient numerical methods for solving the time-dependent scattering and inverse scattering problems of acoustic waves in a locally perturbed half-plane. By symmetric continuation, the scattering problem is reformulated as an equivalent symmetric problem defined in the whole plane. The retarded potential boundary integral equation method is modified to solve the forward problem. Then we consider the inverse scattering problem of determinating the local perturbation from the measured scattered data. The time domain linear sampling method is employed to deal with the inverse problem. The computation schemes proposed in this paper are relatively simple and easy to implement. Several numerical examples are presented to show the effectiveness of the proposed methods.     

相邻多个浅圆弧凹陷地形对平面SH波散射的级数解

本文采用解析法研究相邻多个线圆弧凹陷地形对平面SH波散射问题。文中由分离变量法把相邻多个圆弧凹陷对平面SH波的多重散射表示为各局部坐标中的级数之和,再利用Graf加法公式的内域和外域表达式进一步表示某个局部坐标中的双重无穷级数形式。问题最后可归结为求解一组无穷型的线性代数方程。文末给出了相邻两个等直径浅圆弧凹陷地形的多种深宽比对地面运动的影响的计算结果,并讨论了浅圆凹陷地形对波的屏蔽及它们之间的相互作用。     

Furtivity and Masking Problems in Time-Dependent Electromagnetic Obstacle Scattering

In this paper, we consider furtivity and masking problems in time-dependent three-dimensional electromagnetic obstacle scattering. That is, we propose a criterion based on a merit function to minimize or to mask the electromagnetic field scattered by a bounded obstacle when hit by an incoming electromagnetic field and, with respect to this criterion, we drive the optimal strategy. These problems are natural generalizations to the context of electromagnetic scattering of the furtivity problem in time-dependent acoustic obstacle scattering presented in Ref. 1. We propose mathematical models of the furtivity and masking time-dependent three-dimensional electromagnetic scattering problems that consist in optimal control problems for systems of partial differential equations derived from the Maxwell equations. These control problems are approached using the Pontryagin maximum principle. We formulate the first-order optimality conditions for the control problems considered as exterior problems defined outside the obstacle for systems of partial differential equations. Moreover, the first-order optimality conditions derived are solved numerically with a highly parallelizable numerical method based on a perturbative series of the type considered in Refs. 2–3. Finally, we assess and validate the mathematical models and the numerical method proposed analyzing the numerical results obtained with a parallel implementation of the numerical method in several experiments on test problems. Impressive speedup factors are obtained executing the algorithms on a parallel machine when the number of processors used in the computation ranges between 1 and 100. Some virtual reality applications and some animations relative to the numerical experiments can be found in the website http://www.econ.unian.it/recchioni/w10/.     

Numerical methods for scattering problems from multi-layers with different periodicities

In this paper, we consider a numerical method to solve scattering problems with multi-periodic layers with different periodicities. The main tool applied in this paper is the Bloch transform. With this method, the problem is written into an equivalent coupled family of quasi-periodic problems. As the Bloch transform is only defined for one fixed period, the inhomogeneous layer with another period is simply treated as a non-periodic one. First, we approximate the refractive index by a periodic one where its period is an integer multiple of the fixed period, and it is decomposed by finite number of quasi-periodic functions. Then the coupled system is reduced into a simplified formulation. A convergent finite element method is proposed for the numerical solution, and the numerical method has been applied to several numerical experiments. At the end of this paper, relative errors of the numerical solutions will be shown to illustrate the convergence of the numerical algorithm.     

Positive solutions and dynamics of a finite difference reaction–diffusion system

We investigate the dynamics and methods of computation for some nonlinear finite difference systems that are the discretized equations of a time-dependent and a steady-state reaction–diffusion problem. The formulation of the discrete equations for the time-dependent problem is based on the implicit method for parabolic equations, and the computational algorithm is based on the method of monotone iterations using upper and lower solutions as the initial iterations. The monotone iterative method yields improved upper and lower bounds of the solution in each iteration, and the sequence of iterations converges monotonically to a solution for both the time-dependent and the steady-state problems. An important consequence of this method is that it leads to a bifurcation point that determines the dynamic behavior of the time-dependent problem in relation to the corresponding steady-state problem. This bifurcation point also determines whether the steady-state problem has one or two non-negative solutions, and is explicitly given in terms of the physical parameters of the system and the type of boundary conditions. Numerical results are presented for both the time-dependent and the steady-state problems under various boundary conditions, including a test problem with known analytical solution. These numerical results exhibit the predicted dynamic behavior of the time-dependent solution given by the theoretical analysis. Also discussed are the numerical stability of the computational algorithm and the convergence of the finite difference solution to the corresponding continuous solution of the reaction–diffusion problem. © 1993 John Wiley & Sons, Inc.     

夹层椭圆形板的1/3亚谐解

研究了夹层椭圆形板的非线性强迫振动问题。在以5个位移分量表示的夹层椭圆板的运动方程的基础上,导出了相应的非线性动力方程。提出一类强非线性动力系统的叠加-叠代谐波平衡法。将描述动力系统的二阶常微分方程,化为基本解为未知函数的基本微分方程和派生解为未知函数的增量微分方程。通过叠加-叠代谐波平衡法得出了椭圆板的1/3亚谐解。同时,对叠加-叠代谐波平衡法和数值积分法的精度进行了比较。并且讨论了1/3亚谐解的渐近稳定性。     

Dynamic analysis of anisotropic half space containing an elliptical inclusion under SH waves

Based on the methods of complex function, conformal mapping, and multipolar coordinate system, dynamic response of an elliptical inclusion embedded in an anisotropic half space is investigated. In order to find the solution of SH waves, the governing equation is transferred into its normalized form. Then, the scattering wave induced by the inclusion and the standing wave in the inclusion is deduced. Different incident wave angles and the corresponding anisotropy of the half space are considered to obtain the reflected waves. The elliptical inclusion is transferred into a unit circle by conformal mapping method, and then the undetermined coefficients in scattering wave and standing wave are solved by using the continuous condition at the boundary of the inclusion. Subsequently, the dynamic stress concentration factor (DSCF) around the inclusion is calculated and analyzed. Numerical results demonstrate that the distribution of the DSCF is mainly influenced by the incident wave angle and the incident wave number. It is affected by anisotropic parameters as well.     

Solving the non-isospectral Ablowitz–Ladik hierarchy via the inverse scattering transform and reductions

The non-isospectral Ablowitz–Ladik hierarchy is integrated by the inverse scattering transform. In contrast with the isospectral Ablowitz–Ladik hierarchy, the eigenvalues of the non-isospectral Ablowitz–Ladik equations in the scattering data are time-dependent. The multi-soliton solution for the hierarchy is presented. The reductions to the non-isospectral discrete NLS hierarchy and the non-isospectral discrete mKdV hierarchy and their solutions are considered.     

椭圆孔边裂纹对SH波的散射及其动应力强度因子

采用复变函数和Green函数方法求解具有任意有限长度的椭圆孔边上的径向裂纹对SH波的散射和裂纹尖端处的动应力强度因子．取含有半椭圆缺口的弹性半空间水平表面上任意一点承受时间谐和的出平面线源荷载作用时的位移解作为Green函数,采用裂纹“切割”方法,并根据连续条件建立起问题的定解积分方程,得到动应力强度因子的封闭解答．讨论了孔洞的存在对动应力强度因子的影响．     

Volume integral equations of the scattering problem of poroelasticity and their properties

Scattering of monochromatic waves on an isolated inhomogeneity (inclusion) in an infinite poroelastic medium is considered. Wave propagation in the medium and the inclusion are described by Biot's equations of poroelasticity. The problem is reduced to 3D‐integro‐differential equations for displacement and pressure fields in the region occupied by the inclusion. Properties of the integral operators in these equations are studied. Discontinuities of the fields on the inclusion boundary are indicated. The case of a thin inclusion with low permeability is considered. The corresponding scattering problem is reduced to a 2D integral equation on the middle surface of the inclusion. The unknown function in this equation is the pressure jump in the transverse direction to the inclusion middle surface. An inclusion with a thin layer of low permeability on its interface is considered. The appropriate boundary conditions on the inclusion interface are pointed out. Methods of numerical solution of the volume integral equations of the scattering problems of poroelasticity are discussed.     

The efficient solution of electromagnetic scattering for inhomogeneous media

We consider an electromagnetic scattering problem for inhomogeneous media. In particular, we focus on the numerical computation of the electromagnetic scattered wave generated by the interaction of an electromagnetic plane wave and an inhomogeneity in the corresponding propagation medium. This problem is studied in the VV polarization case, where some special symmetry requirements for the incident wave and for the inhomogeneity are assumed. This problem is reformulated as a Fredholm integral equation of second kind, which is discretized by a linear system having a special form. This allows to compute efficiently an approximate solution of the scattering problem by using iterative techniques for linear systems. Some numerical examples are reported.     

Strong Interactions between Solitary Waves Belonging to Different Wave Modes

Strong interactions between weakly nonlinear long waves are studied. Strong interactions occur when the linear long wave phase speeds are nearly equal although the waves belong to different modes. Specifically we study this situation in the context of internal wave modes propagating in a density stratified fluid. The interaction is described by two coupled Korteweg-deVries equations, which possess both dispersive and nonlinear coupling terms. It is shown that the coupled equations possess an exact analytical solution involving the characteristic “sech²” profile of the Korteweg-deVries equation. It is also shown that when the coefficients satisfy some special conditions, the coupled equations possess an n-solition solution analogous to the Korteweg-deVries n-solition solution. In general though the coupled equations are found not to be amenable to solution by the inverse scattering transform technique, and thus a numerical method has been employed in order to find solutions. This method is described in detail in Appendix A. Several numerical solutions of the coupled equations are presented.     

Models of superconducting–normal–superconducting junctions

A time-dependent Ginzburg–Landau-type model of a superconducting–normal–superconducting junction is presented. The existence and the uniqueness of the solutions are proved. When the data of the model are symmetric of some kinds, the solutions turns out to be symmetric of some kinds. In this symmetric case, an approximate model with the small thickness of the normal material in the middle of the junction as coefficients of a differential system is established for the sake of numerical computations. And also the existence and the uniqueness of the solution to this approximate model are set up. © 1998 by B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

Numerical methods for the linear Boltzmann transport equation in slab geometry

An iterative method for computing numerical solutions of a finite-difference system corresponding to the linear Boltzmann equation in slab geometry is presented. This iterative scheme gives a straightforward marching process starting from the given boundary and initial conditions. It is shown that with a suitable initial iteration the sequence of iterations converges monotonically to a unique solution of the finite-difference system. This monotone convergence leads to improved upper and lower bounds of the solution in each iteration, and to the well-posedness of the discrete system in the sense of Hadamard. It also leads to the convergence of the discrete system to the continuous system as the mesh size of the space–velocity–time variables approaches to zero. Under a mild restriction on the time-increment the discrete system is numerically stable, independent of the mesh-size of the space and velocity. An error estimate for the computed solution due to simultaneous initial and iteration error is obtained. Also given are some numerical results for the time-dependent and the steady-state solutions.     

Finite difference reaction-diffusion systems with coupled boundary conditions and time delays

This paper is concerned with finite difference solutions of a coupled system of reaction-diffusion equations with nonlinear boundary conditions and time delays. The system is coupled through the reaction functions as well as the boundary conditions, and the time delays may appear in both the reaction functions and the boundary functions. The reaction-diffusion system is discretized by the finite difference method, and the investigation is devoted to the finite difference equations for both the time-dependent problem and its corresponding steady-state problem. This investigation includes the existence and uniqueness of a finite difference solution for nonquasimonotone functions, monotone convergence of the time-dependent solution to a maximal or a minimal steady-state solution for quasimonotone functions, and local and global attractors of the time-dependent system, including the convergence of the time-dependent solution to a unique steady-state solution. Also discussed are some computational algorithms for numerical solutions of the steady-state problem when the reaction function and the boundary function are quasimonotone. All the results for the coupled reaction-diffusion equations are directly applicable to systems of parabolic-ordinary equations and to reaction-diffusion systems without time delays.     

Light scattering by a homogeneous sphere near a plane boundary

The scattering of electromagnetic waves by a homogeneous sphere near a plane boundary is presented in this paper. The vector wave equations derived from Maxwell’s equations are solved by means of the two orthogonal solutions to the scalar wave equation. Hankel transformation and Erdélyi’s formula are used to satisfy the planar boundary conditions and the determination of the unknown coefficients in the scattered field and internal fields is achieved by matching the electromagnetic boundary conditions on the surface of the sphere. Existence and uniqueness of the solution of the series involving these unknown coefficients are shown.     

The exact solutions to a Ragnisco-Tu hierarchy with self-consistent sources

The Ragnisco-Tu hierarchy with self-consistent sources is derived. The exact solutions of the hierarchy are obtained via the inverse scattering transform (IST). An explicit form for a solution of the Ragnisco-Tu equation is presented.     

Free vibration analysis of a finite-length isotropic solid elliptic cylinder using exact three dimensional elasticity theory

A three-dimensional analytical model based on Navier’s displacement equation of motion is developed to describe the free vibrations of a simply supported elastic isotropic solid elliptical cylinder of finite length. Helmholtz decomposition theorem is invoked and the classical technique of separation of variables in elliptical coordinates is employed to solve the subsequent uncoupled governing differential equations as products of ordinary and modified Mathieu functions. Inclusive numerical data are provided in a systematic manner for the first several symmetric/anti-symmetric eigen-frequencies as a function of cross sectional eccentricity for selected cylinder length ratios. The incidences of “mode crossing” and “frequency veering” among different modes of the same symmetry group as well as swapping of the corresponding mode shapes in the veering region are examined. Also, selected three dimensional deformation mode shapes are illustrated. The precision of calculations is checked by performing proper convergence studies, and the correctness of results is demonstrated using a commercial finite element software as well as by comparison with the information in the accessible literature. The collection of presented data is believed to be the first meticulous attempt to determine the free vibration frequencies of finite-length simply-supported elastic elliptical cylinders in an exact manner.