Computation of paraxial wave fields using transparent boundary conditions

The parabolic wave equation is solved numerically by applying transparent conditions used to confine the computational domain. A numerical implementation of the boundary conditions is proposed based on representing the incident wave as a superposition of Gaussian beams. A modification of the transparent conditions in the case of dielectric objects extending beyond the computational domain is described. Numerical examples are presented.     

Spectral Method for the Black-Scholes Model of American Options Valuation

In this paper, we devote ourselves to the research of numerical methods for American option pricing problems under the Black-Scholes model. The optimal exercise boundary which satisfies a nonlinear Volterra integral equation is resolved by a high-order collocation method based on graded meshes. For the other spatial domain boundary, an artificial boundary condition is applied to the pricing problem for the effective truncation of the semi-infinite domain. Then, the front-fixing and stretching transformations are employed to change the truncated problem in an irregular domain into a one-dimensional parabolic problem in [−1,1]. The Chebyshev spectral method coupled with fourth-order Runge-Kutta method is proposed for the resulting parabolic problem related to the options. The stability of the semi-discrete numerical method is established for the parabolic problem transformed from the original model. Numerical experiments are conducted to verify the performance of the proposed methods and compare them with some existing methods.     

Wavelet-Galerkin method for solving parabolic equations in finite domains

A novel wavelet-Galerkin method tailored to solve parabolic equations in finite domains is presented. The emphasis of the paper is on the development of the discretization formulations that are specific to finite domain parabolic equations with arbitrary boundary conditions based on weak form functionals. The proposed method also deals with the development of algorithms for computing the associated connection coefficients at arbitrary points. Here the Lagrange multiplier method is used to enforce the essential boundary conditions. The numerical results on a two-dimensional transient heat conducting problem are used to validate the proposed wavelet-Galerkin algorithm as an effective numerical method to solve finite domain parabolic equations.     

抛物型初边值问题的有限元与边界积分耦合的离散化及其误差分析

1．引言边界元方法是近二十几年来迅速发展起来的一类新的偏微分方程的数值方法．它的独特之处是将空间的维数降低一维,从而倍受工程技术人员的青睐,并在工程技术与计算数学领域得到越来越广泛的重视和研究．对椭圆型问题,边界元方法的理论与应用研究已取得丰硕成果;对发展型问题,近年来在理论方面的研究也已取得重要进展［6－11］．但边界元方法难以处理非均质问题,而有限元对各类问题及各种区域具有较好的适应性,将两者结合起来可充分发挥各自的优点．文山提出了一种抛物方程初边值问题的有限元与边界积分的耦合方法,其主要思想是…     

Nonlocal Boundary Conditions for Higher–Order Parabolic Equations

This work deals with the efficient numerical solution of the two–dimensional one–way Helmholtz equation posed on an unbounded domain. In this case one has to introduce artificial boundary conditions to confine the computational domain. Here we construct with the Z –transformation so–called discrete transparent boundary conditions for higher–order parabolic equations schemes. These methods are Padé “Parabolic” approximations of the one–way Helmholtz equation and frequently used in integrated optics and (underwater) acoustics. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analytic and Experimental Studies of the Errors in Numerical Methods for the Valuation of Options

The value of a European option satisfies the Black-Scholes equation with appropriately specified final and boundary conditions.We transform the problem to an initial boundary value problem in dimensionless form.There are two parameters in the coefficients of the resulting linear parabolic partial differential equation.For a range of values of these parameters,the solution of the problem has a boundary or an initial layer.The initial function has a discontinuity in the first-order derivative,which leads to the appearance of an interior layer.We construct analytically the asymptotic solution of the equation in a finite domain.Based on the asymptotic solution we can determine the size of the artificial boundary such that the required solution in a finite domain in x and at the final time is not affected by the boundary.Also,we study computationally the behaviour in the maximum norm of the errors in numerical solutions in cases such that one of the parameters varies from finite (or pretty large) to small values,while the other parameter is fixed and takes either finite (or pretty large) or small values. Crank-Nicolson explicit and implicit schemes using centered or upwind approximations to the derivative are studied.We present numerical computations,which determine experimentally the parameter-uniform rates of convergence.We note that this rate is rather weak,due probably to mixed sources of error such as initial and boundary layers and the discontinuity in the derivative of the solution.     

Composite spectral method for solution of the diffusion equation with specification of energy

Many physical subjects are modeled by nonclassical parabolic boundary value problems with nonlocal boundary conditions replacing the classic boundary conditions. In this article, we introduce a new numerical method for solving the one‐dimensional parabolic equation with nonlocal boundary conditions. The approximate proposed method is based upon the composite spectral functions. The properties of composite spectral functions consisting of terms of orthogonal functions are presented and are utilized to reduce the problem to some algebraic equations. The method is easy to implement and yields very accurate result. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Weak Galerkin finite element method for valuation of American options

We introduce a weak Galerkin finite element method for the valuation of American options governed by the Black-Scholes equation. In order to implement, we need to solve the optimal exercise boundary and then introduce an artificial boundary to make the computational domain bounded. For the optimal exercise boundary, which satisfies a nonlinear Volterra integral equation, it is resolved by a higher-order collocation method based on graded meshes. With the computed optimal exercise boundary, the front-fixing technique is employed to transform the free boundary problem to a one- dimensional parabolic problem in a half infinite area. For the other spatial domain boundary, a perfectly matched layer is used to truncate the unbounded domain and carry out the computation. Finally, the resulting initial-boundary value problems are solved by weak Galerkin finite element method, and numerical examples are provided to illustrate the efficiency of the method.     

A simple and robust boundary treatment for the forced Korteweg–de Vries equation

In this paper, we propose a simple and robust numerical method for the forced Korteweg–de Vries (fKdV) equation which models free surface waves of an incompressible and inviscid fluid flow over a bump. The fKdV equation is defined in an infinite domain. However, to solve the equation numerically we must truncate the infinite domain to a bounded domain by introducing an artificial boundary and imposing boundary conditions there. Due to unsuitable artificial boundary conditions, most wave propagation problems have numerical difficulties (e.g., the truncated computational domain must be large enough or the numerical simulation must be terminated before the wave approaches the artificial boundary for the quality of the numerical solution). To solve this boundary problem, we develop an absorbing non-reflecting boundary treatment which uses outward wave velocity. The basic idea of the proposing algorithm is that we first calculate an outward wave velocity from the solutions at the previous and present time steps and then we obtain a solution at the next time step on the artificial boundary by moving the solution at the present time step with the velocity. And then we update solutions at the next time step inside the domain using the calculated solution on the artificial boundary. Numerical experiments with various initial conditions for the KdV and fKdV equations are presented to illustrate the accuracy and efficiency of our method.     

抛物型初边值问题的自然积分方程及其数值解法

１．引言数值求解无界区域的偏微分方程,自然的处理方式是削去区域的无界部分,即引入一条适当的人工边界ｒ。,将原问题的求解限制在一个适当的有界区域Ｄ内,这样必须在人工边界上引入适当边界的条件．于是很自然地导致这样一个问题：'是否存在一个人工边界条件,使得在这边界条件下,原问题在区域Ｄ内所求得的数值解与原无界区域的解在Ｄ上的限制是完全一致的？"这里我们的着眼点是寻求与原无界区域问题等价的数学形式,以便于数值求解．因为边界元方法可以将区域内的问题转化到区域的边界上去处理,经典的边界元方法常被应用．七十年代…     

A numerical approximation to the solution of an inverse heat conduction problem

The aim of this article is to study the parabolic inverse problem of determination of the leading coefficient in the heat equation with an extra condition at the terminal. After introducing a new variable, we reformulate the problem as a nonclassical parabolic equation along with the initial and boundary conditions. The uniqueness and continuous dependence of the solution upon the data are demonstrated, and then finite difference methods, backward Euler and Crank–Nicolson schemes are studied. The results of some numerical examples are presented to demonstrate the efficiency and the rapid convergence of the methods. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

On an inhomogeneous problem for parabolic-hyperbolic equation

We study a boundary value problem for an inhomogeneous parabolic-hyperbolic equation with a noncharacteristic type change line. Boundary conditions of the first kind are posed on characteristics in the parabolic and hyperbolic parts of the domain where the equation is given, and a condition of the third kind is posed on the noncharacteristic part of the boundary in the parabolic part. First, we study the solvability of an inhomogeneous initial–boundary value problem for a parabolic equation.     

On the combination of some semi‐discretization methods and boundary integral equations for the numerical solution of initial boundary value problems

We consider initial boundary value problems for the homogeneous differential equation of hyperbolic or parabolic type in the unbounded two‐ or three‐dimensional spatial domain with the homogeneous initial conditions and with Dirichlet or Neumann boundary condition. The numerical solution is realized in two steps. At first using the Laguerre transformation or Rothe's method with respect to the time variable the non‐stationary problem is reduced to the sequence of boundary value problems for the non‐homogeneous Helmholtz equation. Further we construct the special integral representation for solutions and obtain the sequence of boundary integral equations (without volume integrals). For the full‐discretization of integral equations we propose some projection methods.     

The operational matrices of Bernstein polynomials for solving the parabolic equation subject to specification of the mass

Some physical problems in science and engineering are modelled by the parabolic partial differential equations with nonlocal boundary specifications. In this paper, a numerical method which employs the Bernstein polynomials basis is implemented to give the approximate solution of a parabolic partial differential equation with boundary integral conditions. The properties of Bernstein polynomials, and the operational matrices for integration, differentiation and the product are introduced and are utilized to reduce the solution of the given parabolic partial differential equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the new technique.     

Boundary-value problem for density-velocity model of collective motion of organisms

The collective motion of organisms is observed at almost all levels of biological systems. In this paper the density-velocity model of the collective motion of organisms is analyzed. This model consists of a system of nonlinear parabolic equations, a forced Burgers equation for velocity and a mass conservation equation for density. These equations are supplemented with the Neumann boundary conditions for the density and the Dirichlet boundary conditions for the velocity. The existence, uniqueness and regularity of solution for the density-velocity problem is proved in a bounded 1D domain. Moreover, a priori estimates for the solutions are established, and existence of an attractor is proved. Finally, some numerical approximations for asymptotical behavior of the density-velocity model are presented.     

A spatial domain decomposition method for parabolic optimal control problems

We present a non-overlapping spatial domain decomposition method for the solution of linear–quadratic parabolic optimal control problems. The spatial domain is decomposed into non-overlapping subdomains. The original parabolic optimal control problem is decomposed into smaller problems posed on space–time cylinder subdomains with auxiliary state and adjoint variables imposed as Dirichlet boundary conditions on the space–time interface boundary. The subdomain problems are coupled through Robin transmission conditions. This leads to a Schur complement equation in which the unknowns are the auxiliary state adjoint variables on the space-time interface boundary. The Schur complement operator is the sum of space–time subdomain Schur complement operators. The application of these subdomain Schur complement operators is equivalent to the solution of an subdomain parabolic optimal control problem. The subdomain Schur complement operators are shown to be invertible and the application of their inverses is equivalent to the solution of a related subdomain parabolic optimal control problem. We introduce a new family of Neumann–Neumann type preconditioners for the Schur complement system including several different coarse grid corrections. We compare the numerical performance of our preconditioners with an alternative approach recently introduced by Benamou.     

半线性抛物型方程奇异摄动问题的数值解

本文讨论具有抛物边界层的半线性抛物型方程奇异摄动问题的数值解法,在非均匀网格上构造了两层非线性差分格式,证明了差分格式是一致收敛的,给出了一些数值例子.     

Transparent boundaries for the parabolic wave equation

The boundary conditions ensuring the full absorption of arbitrary paraxial wave packets, described by the parabolic wave equation, are derived in exact form. This provides a mathematically strict method for reducing an infinite domain to a strip or a cylinder when performing numerical calculations.Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 239, 1997, pp. 211–217.     

三维Poisson方程外问题的高阶局部人工边界条件

１引言假设Ｒ３是一分片光滑的闭曲面．是以为边界的无界区域,=Ｒ3是以为边界的有界区域,并且存在球Ｂ0＝ｘｘＲ0我们考虑下面Ｐｏｉｓｓｏｎ方程的外问题：这里f(x）,ｇ（x）是,上的已知函数,f(x)的支集是紧的,即存在一个球面=x·x=R1,使得x＝xｘＲ１,有ｆx＝０．令＝,则ｆ（ｘ）的支集包含在中,令＝xx＝,表示u在上的外法向微商．用流量为零的条件代替无限远处条件（３）,则我们得到一个新的外问题：我们将分别讨论问题（１）－（３）和（４）－（７）的数值解．由于求解区域的无界性,给数值计算带来了本质性的困难．克服此…     

A non-overlapping domain decomposition dual reciprocity method for solving the forward–backward heat equation in two-dimension

We apply a boundary element dual reciprocity method (DRBEM) to the numerical solution of the forward–backward heat equation in a two-dimensional case. The method is employed for the spatial variable via the fundamental solution of the Laplace equation and the Crank–Nicolson finite difference scheme is utilized to treat the time variable. The physical domain is divided into two non-overlapping subdomains resulting in two standard forward and backward parabolic equations. The subproblems are then treated by the underlying method assuming a virtual boundary in the interface and starting with an initial approximate solution on this boundary followed by updating the solution by an iterative procedure. In addition, we show that the time discrete scheme is unconditionally stable and convergent using the energy method. Furthermore, some computational aspects will be suggested to efficiently deal with the formulation of the proposed method. Finally, two forward–backward problems, for which the exact solution is available, will be numerically solved for two different domains to demonstrate the efficiency of the proposed approach.