LAYER STRIPPING METHOD FOR POTENTIAL INVERSION

1.IntroductionThepotentialinversionproblemofthefollowingPlasmawaveequationisdiscussedinthispaper:Thatis,givinganimpulseatthesurfacez~0,todeterminethewavefieldpandpotentialvfromtheimpulseresponseh.TherearethreekindsofinverseproblemsofthisPlasmawaveequation:(1)TOdeterminethedifferentialequationfromitsspectralfunction[1];(2)Todeterminethepotentialfromthewavefunctionformatlargedistance.Itistheso--calledinversescatteringproblem[2,3];(3)TOdeterminethepotentialfromtheresponseontheboundarytoaunitim…     

Moving collocation method for a reaction‐diffusion equation with a traveling heat source

We consider a reaction‐diffusion equation with a traveling heat source on an unbounded domain. The numerical simulation of the problem is difficult because of the moving singularity, the blow‐up phenomenon, and the delta function in the equation. Because we are only interested in the solution behavior near the heat source, we choose a bounded moving domain which contains the heat source and has the same speed as the source. Local absorbing boundary conditions are constructed on the boundaries of the moving domain. Then, we transform the moving domain to a fixed one. At last, a special moving collocation method is adopted. The new method is much simpler than the existing moving finite difference methods. Moreover, numerical experiments illustrate the accuracy and efficiency of our moving collocation method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Error analysis for solving the Korteweg‐de Vries equation by a Legendre pseudo‐spectral method

A Legendre pseudo‐spectral method is proposed for the Korteweg‐de Vries equation with nonperiodic boundary conditions. Appropriate base functions are chosen to get an efficient algorithm. Error analysis is given for both semi‐discrete and fully discrete schemes. The numerical results confirm to the theoretical analysis. © (2000) John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 513–534, (2000)     

Spectral method for vorticity‐stream function form of Navier–Stokes equations with slip boundary conditions

In this paper, we propose a spectral method for the vorticity‐stream function form of the Navier–Stokes equations with slip boundary conditions. The numerical solutions fulfill the incompressibility and the physical boundary conditions automatically. The stability and convergence of the proposed methods are proven. Numeric results demonstrate the efficiency of suggested algorithm. Copyright © 2011 John Wiley & Sons, Ltd.     

Spectral smoothed boundary methods: The role of external boundary conditions

The spectral smoothed boundary method (SSBM) is a recently proposed numerical method to approximate the solution of partial differential equations in irregular domains with no‐flux boundary conditions by means of Fourier spectral methods. In this article we explore the robustness and accuracy of the scheme under variations of the artificial boundary conditions that must be imposed on the boundary of the enlarged domain in which the problem is solved. As a test model, we present quantitative numerical results based on a problem of propagation of waves of electrical activity in cardiac tissue for which the method is relevant. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Taylor‐Galerkin B‐spline finite element method for the one‐dimensional advection‐diffusion equation

The advection‐diffusion equation has a long history as a benchmark for numerical methods. Taylor‐Galerkin methods are used together with the type of splines known as B‐splines to construct the approximation functions over the finite elements for the solution of time‐dependent advection‐diffusion problems. If advection dominates over diffusion, the numerical solution is difficult especially if boundary layers are to be resolved. Known test problems have been studied to demonstrate the accuracy of the method. Numerical results show the behavior of the method with emphasis on treatment of boundary conditions. Taylor‐Galerkin methods have been constructed by using both linear and quadratic B‐spline shape functions. Results shown by the method are found to be in good agreement with the exact solution. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Singular boundary method based on time‐dependent fundamental solutions for active noise control

Active noise control is an efficient strategy of noise control. A numerical wave shielding model to inhibit wave propagation, which can be considered as an extension of traditional active noise control, is established using the singular boundary method using time‐dependent fundamental solutions in this study. Two empirical formulas to evaluate the origin intensity factors with Dirichlet and Neumann boundary conditions are derived respectively. In comparison with other similar numerical methods, the method can obtain highly accurate results using very few boundary nodes and small CPU time. These meet the major technical requirements of simulation of active noise control. The subsequent numerical experiments show that the proposed model can shield efficiently from the wave propagation for both inner and exterior problems. By applying the newly derived empirical formulas, the CPU time of the singular boundary method is further reduced significantly, which makes the method a competitive new and efficient meshless method. In addition, the singular boundary method makes active noise control in an online manner via time‐dependent fundamental solutions as its basis functions.     

Gauss‐Lobatto‐Legendre‐Birkhoff pseudospectral approximations for the multi‐term time fractional diffusion‐wave equation with Neumann boundaryconditions

A second‐order finite difference/pseudospectral scheme is proposed for numerical approximation of multi‐term time fractional diffusion‐wave equation with Neumann boundary conditions. The scheme is based upon the weighted and shifted Grünwald difference operators approximation of the time fractional calculus and Gauss‐Lobatto‐Legendre‐Birkhoff (GLLB) pseudospectral method for spatial discretization. The unconditionally stability and convergence of the scheme are rigorously proved. Numerical examples are carried out to verify theoretical results.     

Numerical stability study and error estimation for two implicit schemes in a moving boundary problem

Two algorithms are described [Ferris D. H. (fixed time‐step method) and Gupta and Kumar (variable time‐step method)] that solve a mathematical model for the study of the one‐dimensional moving boundary problem with implicit boundary conditions. Landau's transformation is used, in order to work with a fixed number of nodes at each time‐step. The p.d.e. is discretized using an implicit finite difference scheme. The mathematical model describes the oxygen diffusion in absorbing tissues. An important application is the estimation of time‐variant radiation treatments of cancerous tumors. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 42–61, 2000     

STABILITY ANALYSIS OF FINITE ELEMENT METHODS FOR THE ACOUSTIC WAVE EQUATION WITH ABSORBINGBOUNDARY CONDITIONS (PART I)

1.IntroductionInthenumericalsimulationofwavepropagationinunboundedorsemi-unboundedmediumitisnecessarytointroduceartificialboundariestoobtainfinitecomputational'regions.Thensomeboundaryconditionshavetobeimposedontheseboundaries,whichshouldeliminatethe...     

Solving 2D time‐fractional diffusion equations by a pseudospectral method and Mittag‐Leffler function evaluation

Two‐dimensional time‐fractional diffusion equations with given initial condition and homogeneous Dirichlet boundary conditions in a bounded domain are considered. A semidiscrete approximation scheme based on the pseudospectral method to the time‐fractional diffusion equation leads to a system of ordinary fractional differential equations. To preserve the high accuracy of the spectral approximation, an approach based on the evaluation of the Mittag‐Leffler function on matrix arguments is used for the integration along the time variable. Some examples along with numerical experiments illustrate the effectiveness of the proposed approach. Copyright © 2016 John Wiley & Sons, Ltd.     

A method based on meshless approach for the numerical solution of the two‐space dimensional hyperbolic telegraph equation

In the current article, we investigate the RBF solution of second‐order two‐space dimensional linear hyperbolic telegraph equation. For this purpose, we use a combination of boundary knot method (BKM) and analog equation method (AEM). The BKM is a meshfree, boundary‐only and integration‐free technique. The BKM is an alternative to the method of fundamental solution to avoid the fictitious boundary and to deal with low accuracy, singular integration and mesh generation. Also, on the basis of the AEM, the governing operator is substituted by an equivalent nonhomogeneous linear one with known fundamental solution under the same boundary conditions. Finally, several numerical results and discussions are demonstrated to show the accuracy and efficiency of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

A numerical method based on integro‐differential formulation for solving a one‐dimensional Stefan problem

A numerical method based on an integro‐differential formulation is proposed for solving a one‐dimensional moving boundary Stefan problem involving heat conduction in a solid with phase change. Some specific test problems are solved using the proposed method. The numerical results obtained indicate that it can give accurate solutions and may offer an interesting and viable alternative to existing numerical methods for solving the Stefan problem. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Least‐square collocation and Lagrange multipliers forTaylor meshless method

A recently proposed meshless method is discussed in this article. It relies on Taylor series, the shape functions being high degree polynomials deduced from the Partial Differential Equation (PDE). In this framework, an efficient technique to couple several polynomial approximations has been presented in (Tampango, Potier‐Ferry, Koutsawa, Tiem, Int. J. Numer. Meth. Eng. vol. 95 (2013) pp. 1094–1112): the boundary conditions were applied using the least‐square collocation and the interface was coupled by a bridging technique based on Lagrange multipliers. In this article, least‐square collocation and Lagrange multipliers are applied for boundary conditions, respectively, and least‐square collocation is revisited to account for the interface conditions in piecewise resolutions. Various combinations of these two techniques have been investigated and the numerical results prove their effectiveness to obtain very accurate solutions, even for large scale problems.     

Absorbing boundary conditions for the wave equation and parallel computing

Absorbing boundary conditions have been developed for various types of problems to truncate infinite domains in order to perform computations. But absorbing boundary conditions have a second, recent and important application: parallel computing. We show that absorbing boundary conditions are essential for a good performance of the Schwarz waveform relaxation algorithm applied to the wave equation. In turn this application gives the idea of introducing a layer close to the truncation boundary which leads to a new way of optimizing absorbing boundary conditions for truncating domains. We optimize the conditions in the case of straight boundaries and illustrate our analysis with numerical experiments both for truncating domains and the Schwarz waveform relaxation algorithm.

     

The use of cubic B‐spline scaling functions for solving the one‐dimensional hyperbolic equation with a nonlocal conservation condition

Problems for parabolic partial differential equations with nonlocal boundary conditions have been studied in many articles, but boundary value problems for hyperbolic partial differential equations have so far remained nearly uninvestigated. In this article a numerical technique is presented for the solution of a nonclassical problem for the one‐dimensional wave equation. This method uses the cubic B‐spline scaling functions. Some numerical results are reported to support our study. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

An explicit spectral collocation method using nonpolynomial basis functions for the time‐dependent Schrödinger equation

We propose a spectral collocation method for the numerical solution of the time‐dependent Schrödinger equation, where the newly developed nonpolynomial functions in a previous study are used as basis functions. Equipped with the new basis functions, various boundary conditions can be imposed exactly. The preferable semi‐implicit time marching schemes are employed for temporal discretization. Moreover, the new basis functions build in a free parameter λ intrinsically, which can be chosen properly so that the semi‐implicit scheme collapses to an explicit scheme. The method is further applied to linear Schrödinger equation set in unbounded domain. The transparent boundary conditions are constructed for time semidiscrete scheme of the linear Schrödinger equation. We employ spectral collocation method using the new basis functions for the spatial discretization, which allows for the exact imposition of the transparent boundary conditions. Comprehensive numerical tests both in bounded and unbounded domain are performed to demonstrate the attractive features of the proposed method.     

A nonlinear approach to absorbing boundary conditions for the semilinear wave equation

We construct a family of absorbing boundary conditions for the semilinear wave equation. Our principal tool is the paradifferential calculus which enables us to deal with nonlinear terms. We show that the corresponding initial boundary value problems are well posed. We finally present numerical experiments illustrating the efficiency of the method.

     

The fully Sinc‐Galerkin method for time‐dependent boundary conditions

The fully Sinc‐Galerkin method is developed for a family of complex‐valued partial differential equations with time‐dependent boundary conditions. The Sinc‐Galerkin discrete system is formulated and represented by a Kronecker product form of those equations. The numerical solution is efficiently calculated and the method exhibits an exponential convergence rate. Several examples, some with a real‐valued solution and some with a complex‐valued solution, are used to demonstrate the performance of this method. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

Blow‐up phenomena in the model of a space charge stratification in semiconductors: analytical and numerical analysis

The initial‐boundary value problems for a Sobolev equation with exponential nonlinearities, classical, and nonclassical boundary conditions are considered. For this model, which describes processes in crystalline semiconductors, the blow‐up phenomena are studied. The sufficient blow‐up conditions and the blow‐up time are analyzed by the method of the test functions. This analytical a priori information is used in the numerical experiments, which are able to determine the process of the solution's blow‐up more accurately. The model derivation and some questions of local solvability and uniqueness are also discussed. Copyright © 2016 John Wiley & Sons, Ltd.