THE APPROXIMATIONS OF THE EXACT BOUNDARY CONDITION AT AN ARTIFICIAL BOUNDARY FOR LINEARIZED INCOMPRESSIBLE VISCOUS FLOWS

1.IntroductionManyproblemsarisinginfluidmechanicsaregiveninanunboundeddomain,suchasfluidflowaroundobstacles.Whencomputingthenumericalsolutionsoftheseproblems,oneoftenintroducesartificialboundariesandsetsupaxtificialboundaryconditionsonthem.Thentheoriginal…     

Positive harmonic functions of finite order in a Denjoy type domain

We introduce a Denjoy type domain and prove that the dimension of the cone of positive harmonic functions of finite order in the domain with vanishing boundary values is one or two, whenever the boundary is included in a certain set.

     

NUMERICAL SOLUTIONS OF PARABOLIC PROBLEMS ON UNBOUNDED 3-D SPATIAL DOMAIN

In this paper,the numerical solutions of heat equation on 3-D unbounded spatial do-main are considered. n artificial boundary Γ is introduced to finite the computationaldomain.On the artificial boundary Γ,the exact boundary condition and a series of approx-imating boundary conditions are derived,which are called artificial boundary conditions.By the exact or approximating boundary condition on the artificial boundary,the originalproblem is reduced to an initial-boundary value problem on the bounded computationaldomain,which is equivalent or approximating to the original problem.The finite differencemethod and finite element method are used to solve the reduced problems on the finitecomputational domain.The numerical results demonstrate that the method given in thispaper is effective and feasible.     

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.     

THE ARTIFICIAL BOUNDARY CONDITION FOR EXTERIOR OSEEN EQUATION IN 2-D SPACE

A finite element method for the solution of Oseen equation in exterior domain is proposed. In this method, a circular artificial boundary is introduced to make the computational domain finite. Then, the exact relation between the normal stress and the prescribed velocity field on the artificial boundary can be obtained analytically. This relation can serve as an boundary condition for the boundary value problem defined on the finite domain bounded by the artificial boundary. Numerical experiment is presented to demonstrate the performance of the method.     

A NON-OVERLAPPING DOMAIN DECOMPOSITION ALGORITHM BASED ON THE NATURAL BOUNDARY REDUCTION FOR WAVE EQUATIONS IN AN UNBOUNDED DOMAIN

In this paper, a new domain decomposition method based on the natural boundary reduction, which solves wave problems over an unbounded domain, is suggestted. An circular artificial boundary is introduced. The original unbounded domain is divided into two subdomains, an internal bounded region and external unbounded region outside the artificial boundary. A Dirichlet-Neumann(D-N) alternating iteration algorithm is constructed. We prove that the algorithm is equavilent to preconditional Richardson iteration method. Numerical studies are performed by finite element method. The numerical results show that the convergence rate of the discrete D-N iteration is independent of the finite element mesh size.     

Optimal Shape Problems for Eigenvalues

ABSTRACT

We consider the first Dirichlet eigenvalue for nonhomogeneous membranes. For given volume we want to find the domain which minimizes this eigenvalue. The problem is formulated as a variational free boundary problem. The optimal domain is characterized as the support of the first eigenfunction. We prove enough regularity for the eigenfunction to conclude that the optimal domain has finite parameter. Finally an overdetermined boundary value problem on the regular part of the free boundary is given.     

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.     

A finite element convergence analysis for 3D Stokes equations in case of variational crimes

We investigate a finite element discretization of the Stokes equations with nonstandard boundary conditions, defined in a bounded three-dimensional domain with a curved, piecewise smooth boundary. For tetrahedral triangulations of this domain we prove, under general assumptions on the discrete problem and without any additional regularity assumptions on the weak solution, that the discrete solutions converge to the weak solution. Examples of appropriate finite element spaces are given.     

A domain decomposition method for linear exterior boundary value problems

In this paper, we present a domain decomposition method, based on the general theory of Steklov-Poincaré operators, for a class of linear exterior boundary value problems arising in potential theory and heat conductivity. We first use a Dirichlet-to-Neumann mapping, derived from boundary integral equation methods, to transform the exterior problem into an equivalent mixed boundary value problem on a bounded domain. This domain is decomposed into a finite number of annular subregions, and the Dirichlet data on the interfaces is introduced as the unknown of the associated Steklov-Poincaré problem. This problem is solved with the Richardson method by introducing a Dirichlet-Robin-type preconditioner, which yields an iteration-by-subdomains algorithm well suited for parallel computations. The corresponding analysis for the finite element approximations and some numerical experiments are also provided.     

数学物理中的总体和局部场GL方法

为了求解物理化学生物材料和金融中的微分方程,提出了一种总体(Global)和局部(Local)场方法.微分方程的求解区域可以是有限域,无限域,或具曲面边界的部分无限域.其无限域包括有限有界不均匀介质区域.其不均匀介质区域被分划为若干子区域之和.在这含非均匀介质的无限区域,将微分方程的解显式地表示为在若干非均匀介质子区域上和局部子曲面的积分的递归和.把正反算的非线性关系递归地显式化.在无限均匀区域,微分方程的解析解被称为初始总体场.微分方程解的总体场相继地被各个非均匀介质子区域的局部散射场所修正.这种修正过程是一个子域接着另个子域逐步相继地进行的.一旦所有非均匀介质子区域被散射扫描和有限步更新过程全部完成后,微分方程的解就获得了.称其为总体和局部场的方法,简称为GL方法.GL方法完全地不同于有限元及有限差方法,GL方法直接地逐子域地组装逆矩阵而获得解.GL方法无需求解大型矩阵方程,它克服了有限元大型矩阵解的困难.用有限元及有限差方法求解无限域上的微分方程时,人为边界及其上的吸收边界条件是必需的和困难的,人为边界上的吸收边界条件的不精确的反射会降低解的精确度和毁坏反算过程.GL方法又克服了有限元和有限差方法的人为边界的困难.GL方法既不需要任何人为边界又不需要任何吸收边界条件就可以子域接子域逐步精确地求解无限域上的微分方程.有限元和有限差方法都仅仅是数值的方法,GL方法将解析解和数值方法相容地结合起来.提出和证明了三角的格林函数积分方程公式.证明了当子域的直经趋于零时,波动方程的GL方法的数值解收敛于精确解.GL方法解波动方程的误差估计也获得了.求解椭圆型,抛物线型,双曲线型方程的GL模拟计算结果显示出我们的GL方法具有准确,快速,稳定的许多优点.GL方法可以是有网,无网和半网算法.GL方法可广泛应用在三维电磁场,三维弹塑性力学场,地震波场,声波场,流场,量子场等方面.上述三维电磁场等应用领域的GL方法的软件已经由作者研制和发展了。     

An artificial boundary method for the Hull-White model of American interest rate derivatives

The Hull-White (HW) model is a widely used one-factor interest rate model because of its analytical tractability on liquidly traded derivatives, super-calibration ability to the initial term structure and elegant tree-building procedure. As an explicit finite difference scheme, lattice method is subject to some stability criteria, which may deteriorate the computational efficiency for early exercisable derivatives. This paper proposes an artificial boundary method based on the partial differential equations (PDEs) to price interest rate derivatives with early exercise (American) feature under the HW model. We construct conversion factors to extract the market information from the zero-coupon curve and then reduce the infinite computational domain into a finite one by using an artificial boundary on which an exact boundary condition is derived. We then develop an implicit θ-scheme with unconditional stability to solve the PDE in the reduced bounded domain. With a finite computational domain, the optimal exercise strategy can be determined efficiently. Our numerical examples show that the proposed scheme is accurate, robust to the truncation size, and more efficient than the popular lattice method for accurate derivative prices. In addition, the singularity-separating technique is incorporated into the artificial boundary method to enhance accuracy and flexibility of the numerical scheme.     

Partial expansion of a Lipschitz domain and some applications

We show that a Lipschitz domain can be expanded solely near a part of its boundary, assuming that the part is enclosed by a piecewise C ¹ curve. The expanded domain as well as the extended part are both Lipschitz. We apply this result to prove a regular decomposition of standard vector Sobolev spaces with vanishing traces only on part of the boundary. Another application in the construction of low-regularity projectors into finite element spaces with partial boundary conditions is also indicated.     

Convergence of a difference scheme for the heat equation in a long strip by artificial boundary conditions

The numerical solution of the heat equation on a strip in two dimensions is considered. An artificial boundary is introduced to make the computational domain finite. On the artificial boundary, an exact boundary condition is proposed to reduce the original problem to an initial‐boundary value problem in a finite computational domain. A difference scheme is constructed by the method of reduction of order to solve the problem in the finite computational domain. It is proved that the difference scheme is uniquely solvable, unconditionally stable and convergent with the convergence order 2 in space and order 3/2 in time in an energy norm. A numerical example demonstrates the theoretical results.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

A mapping finite difference model for infinite elastic media

A simple mapping finite difference model is presented for the solution of boundary-value problems in the theory of time-harmonic elastic vibrations. The finite problem domain is condensed by mapping into a smaller finite domain using a suitable coordinate transformation. The field equations and the boundary conditions are also appropriately transformed. The radiation condition at infinity is satisfied through a change of the dependent variable. Finite difference forms of the transformed equations are then solved in the mapped domain, subject to the transformed boundary conditions.     

The boundary integral method for the steady rotating Navier–Stokes equations in exterior domain (I): the existence of solution

In this paper, we apply the boundary integral method to the steady rotating Navier–Stokes equations in exterior domain. Introducing some open ball which decomposes the exterior domain into a finite domain and a infinite domain, we obtain a coupled problem by the steady rotating Navier–Stokes equations in finite domain and a boundary integral equation without using the artificial boundary condition. For the coupled problem, we show the existence of solution in a convex set.     

On the stability and the attraction domain of the stationary solution of a singularly perturbed parabolic equation with a multiple root of the degenerate equation

We construct and justify the asymptotics of a boundary layer solution of a boundary value problem for a singularly perturbed second-order ordinary differential equation for the case in which the degenerate (finite) equation has an identically double root. A specific feature of the asymptotics is the presence of a three-zone boundary layer. The solution of the boundary value problem is a stationary solution of the corresponding parabolic equation. We prove the asymptotic stability of this solution and find its attraction domain.     

Numerical solution for exterior problems

Summary We solve the Helmholtz equation in an exterior domain in the plane. A perfect absorption condition is introduced on a circle which contains the obstacle. This boundary condition is given explicitly by Bessel functions. We use a finite element method in the bounded domain. An explicit formula is used to compute the solution out of the circle. We give an error estimate and we present relevant numerical results.     

On the Maxwell system under impedance boundary conditions with memory

We consider an initial boundary value problem for the system of the Maxwell equations in a bounded domain with smooth boundary on a finite time interval with new boundary conditions with memory. In appropriate function spaces, we define and study the nonselfadjoint operator that is generated by the Maxwell operator under a boundary condition with memory. Using the operator method, we prove an existence and uniqueness theorem for a solution to the initial boundary value problem.     

Analysis of free-edge effects by boundary finite element method

The scaled boundary finite element method is a novel semi-analytical analysis technique that combines the advantages of the finite element method and the boundary element method. Only a part of the boundary of the considered domain has to be discretized but nevertheless the method is solely finite element based. The governing equations are solved in the so-called scaling direction analytically, whereas a finite element approximation of the solution is performed in the circumferential directions, which form the boundary of the considered domain. Thus, the numerical effort can be reduced considerably when handling stress concentration problems such as e.g. the free-edge effect in laminated plates. In order to analyze the free-edge effect in a semi-infinite half plane, some kinematic coupling equations have to be introduced, that not only couple the degrees of freedom on the boundary, but also within the non-discretized domain. The implementation of kinematic coupling equations within the method is presented. Finally, the efficiency of the new approach is shown in some benchmark examples. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)