Exponential Time Differencing Method for a Reaction-Diffusion System with Free Boundary

For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.

A Radial Basis Function Method with Improved Accuracy for Fourth Order Boundary Value Problems

Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.

SUFFICIENT AND NECESSARY CONDITIONS ON THE EXISTENCE AND ESTIMATES OF BOUNDARY BLOW-UP SOLUTIONS FOR SINGULAR p-LAPLACIAN EQUATIONS

Let?denote a smooth,bounded domain in R^(N)(N≥2).Suppose that g is a nondecreasing C^(1)positive function and assume that b(x)is continuous and nonnegative inΩ,and that it may be singular on■Ω.In this paper,we provide sufficient and necessary conditions on the existence of boundary blow-up solutions to the p-Laplacian problem△_(p)u=b(x)g(u)for x∈Ω,u(x)→+∞as dist(x,■Ω)→0.The estimates of such solutions are also investigated.Moreover,when b has strong singularity,the nonexistence of boundary blow-up(radial)solutions and infinitely many radial solutions are also considered.

HIGH ACCURACY FINITE VOLUME ELEMENT METHOD FOR TWO-POINT BOUNDARY VALUE PROBLEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.

Implicit finite difference method for fractional percolation equation with Dirichlet and fractional boundary conditions

An implicit finite difference method is developed for a one-dimensional frac- tional percolation equation （FPE） with the Dirichlet and fractional boundary conditions. The stability and convergence are discussed for two special cases, i.e., a continued seep- age flow with a monotone percolation coefficient and a seepage flow with the fractional Neumann boundary condition. The accuracy and efficiency of the method are checked with two numerical examples.

A STUDY ON THE WEIGHT FUNCTION OF THE MOVING LEAST SQUARE APPROXIMATION IN THE LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.

A numerical method based on boundary integral equations and radial basis functions for plane anisotropic thermoelastostatic equations with general variable coefficients

A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable coe?cients. The equations describe the thermoelastic behaviors of nonhomogeneous anisotropic materials with properties that vary smoothly from point to point in space. No restriction is imposed on the spatial variations of the thermoelastic coe?cients as long as all the requirements of the laws of physics are satis?ed. To check the validity and accuracy of the proposed numerical method, some speci?c test problems with known solutions are solved.

Diagonal form fast multipole boundary element method for 2D acoustic problems based on Burton-Miller boundary integral equation formulation and its applications

This paper describes formulation and implementation of the fast multipole boundary element method （FMBEM） for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the FMBEM are described in details. They are moment calculation, moment to moment translation, moment to local translation, and local to local translation. A data structure for the quad-tree construction is proposed which can facilitate implementation. An analytical moment expression is derived, which is more accurate, stable, and efficient than direct numerical computation. Numerical examples are presented to demonstrate the accuracy and efficiency of the FMBEM, and radiation of a 2D vibration rail mode is simulated using the FMBEM.

RESEARCH ON THE COMPANION SOLUTION FOR A THIN PLATE IN THE MESHLESS LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless local boundary integral equation method is a currently developed numerical method, which combines the advantageous features of Galerkin finite element method(GFEM), boundary element method(BEM) and element free Galerkin method(EFGM), and is a truly meshless method possessing wide prospects in engineering applications. The companion solution and all the other formulas required in the meshless local boundary integral equation for a thin plate were presented, in order to make this method apply to solve the thin plate problem.

A Generalized FDTD Method with Absorbing Boundary Condition for Solving a Time-Dependent Linear Schrodinger Equation

The Finite-Difference Time-Domain (FDTD) method is a well-known technique for the analysis of quantum devices. It solves a discretized Schrodinger equation in an iterative process. However, the method provides only a second-order accurate numerical solution and requires that the spatial grid size and time step should satisfy a very restricted condition in order to prevent the numerical solution from diverging. In this article, we present a generalized FDTD method with absorbing boundary condition for solving the one-dimensional (1D) time-dependent Schr?dinger equation and obtain a more relaxed condition for stability. The generalized FDTD scheme is tested by simulating a particle moving in free space and then hitting an energy potential. Numerical results coincide with those obtained based on the theoretical analysis.

THE COUPLING OF BOUNDARY ELEMENT AND FINITE ELEMENT METHODS FOR THE EXTERIOR NONSTATIONARY NAVIER-STOKES EQUATIONS

In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.

FINITE ELEMENT METHOD FOR SOLVING TWO-DIMENSIONAL DIFFUSION-REACTION EQUATIONS OF BOUNDARY LAYER TYPE IN POROUS CATALYST PELLET

In this paper,finite element method(FEM)is used to solve two-dimensional diffu-sion-reaction equations of boundary layer type.This kind of equations are usually too complicatedand diffcult to be solved by applying the traditional methods used in chemical engineering becauseof the steep gradients of concentration and temperature.But,these difficulties are easy to be over-comed when the FEM is used.The integraded steps of solving this kind of problems by the FEMare presented in this paper.By applying the FEM to the two actual examples,the conclusion can bereached that the FEM has the advantages of simplicity and good accuracy.

ON FICTITIOUS DOMAIN METHOD FOR THE NUMERICAL SOLUTION TO HEAT CONDUCTION EQUATION WITH DERIVATIVE BOUNDARY CONDITIONS

This paper investigates the stability and convergence of some knowndifference schemes for the numerical solution to heat conduction equation withderivative boundary conditions by the fictitious domain method.The discrete vari-ables at the false mesh points are firstly eliminated from the difference schemes andthe local truncation errors are then analyzed in detail.The stability and convergenceof the schemes are proved by energy method.An improvement is proposed to obtainbetter schemes over the original ones.Several numerical examples and comparisonswith other schemes are presented.

SINGULAR INTEGRAL EQUATIONS AND BOUNDARY ELEMENT METHOD OF CRACKS IN THERMALLY STRESSED PLANAR SOLIDS

Using the method of the boundary integral equation, a set of singular integral equations of the hear transfer problems and the thermo-elastic problems of a crack embedded in a two-dimensional finite body is derived, and then,its numerical method is proposed by the numerical method of the singular integral equations combined with boundary element method. Moreover, the singular nature of temperature gradient field near the crack front is proved by the main-part analysis method of the singular integral equation, and the singular temperature gradients are exactly obtained. Finally, several typical examples calculated.

Improved non-singular local boundary integral equation method

When the source nodes are on the global boundary in the implementation of local boundary integral equation method （LBIEM）, singularities in the local boundary integrals need to be treated specially. In the current paper, local integral equations are adopted for the nodes inside the domain and moving least square approximation （MLSA） for the nodes on the global boundary, thus singularities will not occur in the new al- gorithm. At the same time, approximation errors of boundary integrals are reduced significantly. As applications and numerical tests, Laplace equation and Helmholtz equation problems are considered and excellent numerical results are obtained. Furthermore, when solving the Helmholtz problems, the modified basis functions with wave solutions are adapted to replace the usually-used monomial basis functions. Numerical results show that this treatment is simple and effective and its application is promising in solutions for the wave propagation problem with high wave number.

Simulating regular wave propagation over a submerged bar by boundary element method model and Boussinesq equation model

Numerical models based on the boundary element method and Boussinesq equation are used to simulate the wave transform over a submerged bar for regular waves.In the boundary-element-method model the linear element is used,and the integrals are computed by analytical formulas.The Boussinesq-equation model is the well-known FUNWAVE from the University of Delaware.We compare the numerical free surface displacements with the laboratory data on both gentle slope and steep slope,and find that both the models simulate the wave transform well.We further compute the agreement indexes between the numerical result and laboratory data,and the results support that the boundary-element-method model has a stable good performance,which is due to the fact that its governing equation has no restriction on nonlinearity and dispersion as compared with Boussinesq equation.

A Conceptual Numerical Model of the Wave Equation Using the Complex Variable Boundary Element Method

In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.

Mechanical quadrature methods and extrapolation for solving nonlinear boundary Helmholtz integral equations

This paper presents mechanical quadrature methods （MQMs） for solving nonlinear boundary Helmholtz integral equations. The methods have high accuracy of order O（h3） and low computation complexity. Moreover, the mechanical quadrature methods are simple without computing any singular integration. A nonlinear system is constructed by discretizing the nonlinear boundary integral equations. The stability and convergence of the system are proved based on an asymptotical compact theory and the Stepleman theorem. Using the h3-Richardson extrapolation algorithms （EAs）, the accuracy to the order of O（h5） is improved. To slove the nonlinear system, the Newton iteration is discussed extensively by using the Ostrowski fixed point theorem. The efficiency of the algorithms is illustrated by numerical examples.

A Modified Homotopy Analysis Method for Solving Boundary Layer Equations

A new modification of the Homotopy Analysis Method (HAM) is presented for highly nonlinear ODEs on a semi-infinite domain. The main advantage of the modified HAM is that the number of terms in the series solution can be greatly reduced;meanwhile the accuracy of the solution can be well retained. In this way, much less CPU is needed. Two typical examples are used to illustrate the efficiency of the proposed approach.

The Solution to the Form of Poisson Equation by the Boundary Element Methods

In this paper, the domain integral of the form of Poisson equation is translatedinto complete boundary integral by the fundamental solution of higher-order Laplaceoperator, the dimensions of the problem can be contracted into one. The numericalexamples for Stokes equations show that this method is efficient.