用小波—配点法求解一类有奇异性的微分方程

本文用一个含小参数ε的微分方程作为算例，用小波－配点法进行求解，通过不断地对Ｖ空间进行小波的正交分解来提高分辨率，并利用小波的局部化特征，得到了很好的结果。不但看出微分方程的奇异性的存在，而且给出了奇点的位置，也求出了方程的数值解。     

利用第二类Chebyshev小波求二阶常微分方程组边值问题的数值解

给出了一个求二阶常微分方程组边值问题数值解的第二类Chebyshev小波配点法.利用第二类Chebyshev小波积分算子矩阵,将问题转化成代数方程组的运算.数值例子说明了方法的准确性及易操作性.另外,为了表明方法的高精度性和有效性,数值算例结果与解析解,以及运用变分迭代法,B样条配点法,连续遗传算法等得到的结果进行了比较.     

指数时程差分与有理谱配点法求解奇异摄动Burgers-Huxley问题

带小参数ε的Burgers-Huxley方程是一类非线性、非定常奇异摄动初边值问题,本文用指数时程差分与有理谱配点法求其数值解.对空间方向的边界层,用带sinh变换的有理谱配点法便Chebyshev节点在边界层处加密,只需取较少节点即可达到较高精度;时间方向采用指数时程差分与4阶Runge-Kutta法相结合的格式,并用围线积分计算矩阵甬数的方法克服了求解奇异摄动问题时遇到的的数值不稳定堆题.数值实验表明,本文提出的方法在求解左、右边界层和内部层的奇异摄动Burgers-Huxley问题都有较高的精度.     

变系数线性微分方程初值问题数值解的小波方法

通过利用小波尺度函数的正交性并结合配点法 ,本文给出了一种求解变系数线性微分方程初值问题数值解的小波算法 .在一定的假设条件下 ,对算法的收敛性进行了理论分析 .最后 ,我们还给出了一个具体的数值计算例子 .     

求解二维Poisson方程的重心插值配点法

采用重心Lagrange插值配点法计算了二维Poisson方程.采用重心Lagrange插值法构造近似函数,由配点法离散Poisson方程及其边界条件.数值算例表明方法具有理论简单、计算精度高的特点.     

Black-Scholes期权定价方程的自适应小波算法

研究Black-Scholes期权定价方程的自适应算法,对Black-Scholes方程设计插值小波配点离散格式,然后设计自适应算法,该算法能够自动在一个接近最优的网格上找到B-S模型的解,数值试验表明其高效性.     

基于时域配点法的充液储箱系统多模态方程求解

采用时域配点法研究了充液储箱系统多模态方程的稳态周期解.在模型求解过程中,利用牛顿迭代法求解了配点法得到的非线性代数方程组,而牛顿迭代的初值来自谐波平衡法求解得到的低阶谐波近似.数值仿真结果验证了时域配点法的有效性,并验证以二倍激励频率为基频的第二模态的假设形式更为有效.最终通过对比谐波系数数量级提出一种更为简洁有效的模态表达形式.     

瞬态问题稳定性及弥散分析的无网格RKPM配点法研究

介绍了基于强形式的RKPM配点法求解瞬态动力问题的算法,并提出了采用RKPM配点法,配合时间域中心差分求解二阶波动方程的稳定性评价方法,并通过数值算例验证了此方法的正确性．此评价方法可以方便有效地评估出实际计算时的临界时间步长．通过数值算例比较可知,实际算例的计算临界时间步长与本评价方法,所预测的临界时间步长结果非常接近．给出了如何合理地选择RKPM形函数支撑域的建议．最后与径向基函数配点法进行了对比研究．     

SAV/重心插值配点法求解Allen-Cahn方程北大核心CSCD

采用标量辅助变量(scalar auxiliary variable, SAV)方法结合重心插值配点法求解二维Allen-Cahn方程.在时间方向上分别采用Crank-Nicolson格式、二阶向后差分格式离散,空间方向上采用重心Lagrange插值配点法离散,建立了两种无条件能量稳定SAV格式,并给出了重心插值配点格式的逼近性质.数值实验表明:两种SAV配点格式的时间收敛阶为二阶,并满足能量递减规律.与空间采用有限差分法离散对比,重心Lagrange配点格式具有指数收敛的特性.     

用改进的区间样条小波配点法求解Burgers方程

提出了一种用两个一阶导数矩阵的乘积替代二阶导算子矩阵的新算法，用两种不同的区间样条件小波插值算法对含不同参数的Burgers方程进行了验算，结果表明，对于稍大粘性系数的Burgers方程，替代算法的震荡明显小于原算法，因而替代算法有较大的稳定性范围。     

Stabilization of the trial method for the Bernoulli problem in case of prescribed Dirichlet data

We apply the trial method for the solution of Bernoulli's free boundary problem when the Dirichlet boundary condition is imposed for the solution of the underlying Laplace equation, and the free boundary is updated according to the Neumann boundary condition. The Dirichlet boundary value problem for the Laplacian is solved by an exponentially convergent boundary element method. The update rule for the free boundary is derived from the linearization of the Neumann data around the actual free boundary. With the help of shape sensitivity analysis and Banach's fixed‐point theorem, we shed light on the convergence of the respective trial method. Especially, we derive a stabilized version of this trial method. Numerical examples validate the theoretical findings.Copyright © 2014 John Wiley & Sons, Ltd.     

The Dual Singular Function Method for 2D Boundary Integral Equations

The accuracy of standard boundary element methods for elliptic boundary value problems deteriorates if the boundary of the domain contains corners or if the boundary conditions change along the boundary. Here we first investigate the convergence behaviour of standard spline Galerkin approximation on quasi-uniform meshes for boundary integral equations on polygonal domains. It turns out, that the order of convergence depends on some constant describing the singular behaviour of solutions near corner points of the boundary. In order to recover the full order of convergence for the Galerkin approximation we propose the dual singular function method which is often used for improving the accuracy of finite element methods. The theoretical convergence results are confirmed and illustrated by a numerical example.     

On the coupled NBEM and FEM for a class of nonlinear exterior Dirichlet problem in R~(2*)

In this paper, based on the natural boundary reduction advanced by Feng and Yu, we couple the finite element approach with the natural boundary element method to study the weak solvability and Galerkin approximation of a class of nonlinear exterior boundary value problems. The analysis is mainly based on the variational formulation with constraints. We prove the error estimate of the finite element solution and obtain     

一维抛物型方程一种新的边值问题的解法

在建立目标一维热红外温度模型时,提出了热传导方程的一个新的边值问题;通过比较,选择了GE差分格式,求解方程,并进行了稳定性分析;采用虚拟网格点法处理边界条件,得出了GE格式的完整形式;计算实例表明,分组显示方法更适合此类边值问题的实际计算.     

一种改进边界元法在弹性扭转问题中的应用

本文用一种改进边界元法分析与计算了椭圆截面等直杆的扭转问题.并与边界元法的解进行比较,其结果极为符合.然而,改进边界元法较边界元法所需要的数据量少得多,计算时间也将大大减少了.因此,本文方法对求解Poisson方程问题是一种经济而行之有效的数值计算方法.     

The boundary element method with Lagrangian multipliers

On open surfaces, the energy space of hypersingular operators is a fractional order Sobolev space of order 1/2 with homogeneous Dirichlet boundary condition (along the boundary curve of the surface) in a weak sense. We introduce a boundary element Galerkin method where this boundary condition is incorporated via the use of a Lagrangian multiplier. We prove the quasi‐optimal convergence of this method (it is slightly inferior to the standard conforming method) and underline the theory by a numerical experiment. The approach presented in this article is not meant to be a competitive alternative to the conforming method but rather the basis for nonconforming techniques like the mortar method, to be developed. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

The boundary element formulation for multiparameter structural shape optimization

A problem of optimal shaping of a free and loaded boundary has been formulated in terms of the boundary element method. A maximal stiffness criterion together with the limitation of volume of a body has been used. The problem has been solved with the application of integral optimality conditions being defined on an optimal boundary. Discretization of a boundary and optimality conditions respectively, with the aid of boundary elements, leads to an iterative procedure, from which one can determine the shape parameters generating the position of the optimal boundary.     

An Analogue-difference Method and Application to Induction Motor Models

The paper established a so-called analogue-difference method (ADM) to compute the numerical solutions for boundary value problems of higher-order differential equations, which can be a fundamental method and performs much better than the finite difference method (FDM), even for second-order boundary value problems. Numerical examples and results illustrate the simplicity, efficiency and applicability of the method, which also show that the proposed method has obvious advantages over the methods presented by recent state-of-the-art work for induction motor models.     

Application of reproducing kernel method to third order three-point boundary value problems

This paper investigates the analytical approximate solutions of third order three-point boundary value problems using reproducing kernel method. The solution obtained by using the method takes the form of a convergent series with easily computable components. However, the reproducing kernel method can not be used directly to solve third order three-point boundary value problems, since there is no method of obtaining reproducing kernel satisfying three-point boundary conditions. This paper presents a method for solving reproducing kernel satisfying three-point boundary conditions so that reproducing kernel method can be used to solve third order three-point boundary value problems. Results of numerical examples demonstrate that the method is quite accurate and efficient for singular second order three-point boundary value problems.     

改进的无奇异局部边界积分方程方法

在局部边界积分方程方法中,当源节点位于分析域的整体边界上时,局部边界积分将出现奇异积分问题,这些奇异积分需要做特别的处理．为此,提出了对域内节点采用局部积分方程,而对边界节点直接采用移动最小二乘近似函数引入边界条件来解决奇异积分问题,这同时也解决了对积分边界进行插值引入近似误差的问题．作为应用和数值实验,对Laplace方程和Helmholtz方程问题进行了分析,取得了很好的数值结果．进而,在Helmholtz方程求解中,采用了含波解信息的修正基函数来代替单项式基函数进行近似．数值结果显示,这样处理是简单高效的,在高波数声传播问题的求解中非常具有前景．