Inverse resistivity problem: Geoelectric uncertainty principle and numerical reconstruction method

Mathematical model of vertical electrical sounding by using resistivity method is studied. The model leads to an inverse problem of determination of the unknown leading coefficient (conductivity) of the elliptic equation in R²$R^2$ in a slab. The direct problem is obtained in the form of mixed BVP in axisymmetric cylindrical coordinates. The additional (available measured) data is given on the upper boundary of the slab, in the form of tangential derivative. Due to ill-conditionedness of the considered inverse problem the logarithmic transformation is applied to the unknown coefficient and the inverse problem is studied as a minimization problem for the cost functional, with respect to the reflection coefficient. The Conjugate Gradient method (CGM) is applied for the numerical solution of this problem. Computational experiments were performed with noise free and random noisy data.     

Adaptive error feedback regulation problem for 1D wave equation

By using the adaptive control approach, we solve the error feedback regulator problem for the one‐dimensional wave equation with a general harmonic disturbance anticollocated with control and with two types of disturbed measurements, ie, one collocated with control and the other anti‐collocated with control. Different from the classical error feedback regulator design, which is based on the internal mode principle, we give the adaptive servomechanism design for the system by making use of the measured tracking error (and its time derivative) and the estimation mechanism for the parameters of the disturbance and of the unknown reference. Constructing auxiliary systems and observer and applying the backstepping method for infinite‐dimensional system play important roles in the design. The control objective, which is to regulate the tracking error to zero and to keep the states bounded, is achieved.     

Seismic inverse problems: determining seismic wave speeds using arrival times

The speeds of propagation of seismic waves are computed by detecting the effect of the waves that propagate when a pulse of energy, in the form of a small explosion, is released near the surface of the ground. The pulse is detected and recorded by underground sensors. The times between the initiation of the pulse and the point when it reaches the various sensors are denoted the travel times. In this paper we consider two questions. Is it possible to compute the wave speed profile from accurately measured travel times? If so, how can this be done?     

Computing Green's function of the initial-boundary value problem for the wave equation in a layered cylinder

A new analytical method for the approximate computation of the time-dependent Green's function for the initial-boundary value problem of the three-dimensional wave equation on multi-layered bounded cylinder is suggested in this paper. The method is based on the derivation of eigenvalues and eigenfunctions for an ordinary differential equation with piecewise constant coefficients, and an approximate computation of Green's function in the form of the Fourier series with a finite number of terms relative to the orthogonal set of the derived eigenfunctions. The computational experiment confirms the robustness of the method.     

Reconstruction of the corrosion boundary for the Laplace equation by using a boundary collocation method

In this paper, we consider the identification of a corrosion boundary for the two-dimensional Laplace equation. A boundary collocation method is proposed for determining the unknown portion of the boundary from the Cauchy data on a part of the boundary. Since the resulting matrix equation is badly ill-conditioned, a regularized solution is obtained by employing the Tikhonov regularization technique, while the regularization parameter is provided by the generalized cross-validation criterion. Numerical examples show that the proposed method is reasonable and feasible.     

Distributed feedback design for systems governed by the wave equation

This article deals with the geometric control of a one-dimensional non-autonomous linear wave equation. The idea consists in reducing the wave equation to a set of first-order linear hyperbolic equations. Then, based on geometric control concepts, a distributed control law that enforces the exponential stability and output tracking in the closed-loop system is designed. The presented control approach is applied to obtain a distributed control law that brings a stretched uniform string, modelled by a wave equation with Dirichlet boundary conditions, to rest in infinite time by considering the displacement of the middle point of the string as the controlled output. The controller performances have been evaluated in simulation by considering both tracking and disturbance rejection problems. The robustness of the controller has also been studied when the string tension is subjected to sudden variations.     

Uniqueness and numerical scheme for the Robin coefficient identification of the time-fractional diffusion equation

We study an inverse problem of determining the Robin coefficient of fractional diffusion equation from a nonlocal boundary condition. Based on the property of Caputo fractional derivative, the uniqueness is proved. The numerical schemes for the direct problem and the inverse problem are developed. Three examples are given to show the effectiveness of the presented methods.     

Exact and fast numerical algorithms for the stochastic wave equation

On the basis of integral representations we propose fast numerical methods to solve the Cauchy problem for the stochastic wave equation without boundaries and with Dirichlet boundary conditions. The algorithms are exact in a probabilistic sense.     

Numerical procedure for the determination of an unknown parameter in a parabolic equation

A numerical procedure for an inverse problem of determining unknown source parameter of one-dimensional parabolic equation subject to the specification of the solution at internal point along with the usual initial boundary conditions is considered. By using some transformation the problem is reformulated to a nonlocal parabolic problem. Some numerical examples using the proposed numerical procedure are presented.     

基于变系数Gardner方程的海洋内波研究

本文主要是基于同时含有二阶和三阶非线性项的变系数Gardner方程对海洋内孤立波的传播特性开展研究。在吕宋海峡海域,展示了下降型海洋内波的传播特性及其在SAR图像上的信号特征,并着重分析讨论了耗散项和微扰项对海洋内波所引起的表层流速变化的影响。     

Control of the wave equation at a rational point

The wave equation is controlled by a force at one point, which is a rational point of the interval [−1, 1]: geometric interpretations are given of the controllable subclasses.     

On the regional control problem of the bilinear wave equation

The present research deals with regional optimal control problem of the bilinear wave equation evolving on a spatial domain $\mathrm{\Omega }\subset {\mathrm{ℝ}}^n,n\ge 1$. Such an equation is excited by bounded controls that act on the velocity term. It addresses the tracking of a desired state all over the time interval $[0,T]$ only on a subregion $\omega$ of $\mathrm{\Omega }$ with minimum energy. Then, we prove that an optimal control exists and is characterized as a solution to an optimality system. Algorithm for the computation of such a control is given and successfully illustrated through simulations.     

Dirichlet feedback control for the stabilization of the wave equation: a numerical approach

In (J. Differential Equations 66 (1987) 340) a uniform stabilization method of the wave equation by boundary control à la Dirichlet has been discussed. In this article, we investigate the numerical implementation of the above stabilization process by a numerical scheme which mimics the energy decay properties of its continuous counterpart. The practical implementation of that scheme leads to a biharmonic problem of a new type which is solved by a method directly inspired by some related work of Glowinski and Pironneau on the solution of the Dirichlet problem for the biharmonic operator (SIAM Rev. 21(2) (1979) 167). Numerical experiments show that the decay properties of the energy are well-preserved by our numerical methodology.     

A linearized implicit finite-difference method for solving the equal width wave equation

A linearized implicit finite-difference method is presented to find numerical solutions of the equal width wave equation. The method has been used successfully to investigate the motion of a single solitary wave, the development of the interaction of two solitary waves and an undular bore. The obtained results are compared with other numerical results in the literature. A stability analysis of the scheme is also investigated.     

PETOOL: MATLAB-based one-way and two-way split-step parabolic equation tool for radiowave propagation over variable terrain

A MATLAB-based one-way and two-way split-step parabolic equation software tool (PETOOL) has been developed with a user-friendly graphical user interface (GUI) for the analysis and visualization of radio-wave propagation over variable terrain and through homogeneous and inhomogeneous atmosphere. The tool has a unique feature over existing one-way parabolic equation (PE)-based codes, because it utilizes the two-way split-step parabolic equation (SSPE) approach with wide-angle propagator, which is a recursive forward–backward algorithm to incorporate both forward and backward waves into the solution in the presence of variable terrain. First, the formulation of the classical one-way SSPE and the relatively-novel two-way SSPE is presented, with particular emphasis on their capabilities and the limitations. Next, the structure and the GUI capabilities of the PETOOL software tool are discussed in detail. The calibration of PETOOL is performed and demonstrated via analytical comparisons and/or representative canonical tests performed against the Geometric Optic (GO) + Uniform Theory of Diffraction (UTD). The tool can be used for research and/or educational purposes to investigate the effects of a variety of user-defined terrain and range-dependent refractivity profiles in electromagnetic wave propagation.

Program summary

Program title: PETOOL (Parabolic Equation Toolbox)Catalogue identifier: AEJS_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEJS_v1_0.htmlProgram obtainable from: CPC Program Library, Queen?s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 143 349No. of bytes in distributed program, including test data, etc.: 23 280 251Distribution format: tar.gzProgramming language: MATLAB (MathWorks Inc.) 2010a. Partial Differential Toolbox and Curve Fitting Toolbox requiredComputer: PCOperating system: Windows XP and VistaClassification: 10Nature of problem: Simulation of radio-wave propagation over variable terrain on the Earth?s surface, and through homogeneous and inhomogeneous atmosphere.Solution method: The program implements one-way and two-way Split-Step Parabolic Equation (SSPE) algorithm, with wide-angle propagator. The SSPE is, in general, an initial-value problem starting from a reference range (typically from an antenna), and marching out in range by obtaining the field along the vertical direction at each range step, through the use of step-by-step Fourier transformations. The two-way algorithm incorporates the backward-propagating waves into the standard one-way SSPE by utilizing an iterative forward–backward scheme for modeling multipath effects over a staircase-approximated terrain.Unusual features: This is the first software package implementing a recursive forward–backward SSPE algorithm to account for the multipath effects during radio-wave propagation, and enabling the user to easily analyze and visualize the results of the two-way propagation with GUI capabilities.Running time: Problem dependent. Typically, it is about 1.5 ms (for conducting ground) and 4 ms (for lossy ground) per range step for a vertical field profile of vector length 1500, on Intel Core 2 Duo 1.6 GHz with 2 GB RAM under Windows Vista.     

A matrix formulation to the wave equation with non-local boundary condition

Recently, various sequential numerical schemes have been proposed for the solution of non-classical hyperbolic initial value problems which involve non-local integral terms over the spatial domain. In this paper, we focus on the wave equation with the non-local boundary condition. Two matrix formulation techniques based on the shifted standard and shifted Chebyshev bases are proposed for the numerical solution. Several numerical examples and also some comparisons with another methods will be investigated to confirm the efficiency of this procedure.     

Non-polynomial splines method for numerical solutions of the regularized long wave equation

In this paper, we construct a numerical method based on cubic splines in tension for solving regularized long wave equation. The truncation error is analysed and the method shows that by choosing suitably parameters we can obtain various accuracy schemes. Numerical stability of the method has been studied by using a linearized stability analysis. Test problems are dealt with. The numerical simulations can validate and demonstrate the advantages of the method.     

Boundary stabilisation of the wave equation in the presence of singularities

We study the boundary stabilisation of the wave equation by a nonlinear feedback active on a part of the boundary in geometric situations for which the solutions have singularities. These singularities appear at the interfaces at which the mixed Neumann–Dirichlet boundary conditions meet. Under a simple geometrical condition concerning the orientation of the boundary, we obtain sharp energy decay rates under a general growth assumption on the feedback. We show that the singularities do not affect the energy decay rates and give examples.     

Determination of inner boundaries in modified Helmholtz inverse geometric problems using the method of fundamental solutions

In this paper, an inverse geometric problem for the modified Helmholtz equation arising in heat conduction in a fin, which consists of determining an unknown inner boundary (rigid inclusion or cavity) of an annular domain from a single pair of boundary Cauchy data is solved numerically using the method of fundamental solutions (MFS). A nonlinear minimisation of the objective function is regularised when noise is added into the input boundary data. The stability of numerical results is investigated for several test examples.