带非局部边界条件的热方程的LEGENDRE配置法

对带非局部边界条件的热方程的初边值问题提出了LEGENDRE配置法,并给出其半离散逼近和全离散逼近的稳定性和收敛性分析.数值试验验证了方法的有效性.     

直接配置法求解连续铸钢过程中的最优控制问题

利用描述连续铸钢过程二冷区喷水控制下钢的热传导的半离散化模型 ,我们构造一包含温度梯度约束的最优控制问题 .针对此最优控制问题 ,采用直接配置法进行数值求解 ,得出相应的近似最优控制 .     

KdV方程的Chebyshev-Hermite谱配置法

针对无界区域上Korteweg.-de Vries(KdV)方程构造了时空全离散的ChebyshevHermite谱配置格式,即在空间方向上采用Hermite谱配置方法离散,时间方向上采用Chebyshev谱配置方法离散.提出了一个简单迭代算法,该算法非常适合并行计算.数值结果显示了此算法的有效性.     

关于非线性双曲型方程全离散有限元方法的稳定性和收敛性估计

本文研究非线性双曲型方程混合问题的有限元方法.这类问题的研究,对于非线性振动、渗流力学等实际问题,在理论和实用方面均有价值.关于线性、半线性双曲方程全离散有限元方法及非线性双曲方程半离散有限元方法的收敛性研究,已有[1]—[4].     

二维土壤溶质输运问题的CN有限元法

首先给出二维土壤溶质输运问题时间二阶精度的Crank-Nicolson(CN)时间半离散化格式,然后直接从CN时间半离散化格式出发,建立具有时间二阶精度的全离散化CN有限元格式,并给出CN有限元解的误差分析,最后用数值例子验证全离散化CN有限元格式的优越性.这种方法提高了时间离散的精度,并极大地减少时间方向的迭代步,从而减少实际计算中截断误差的积累,提高计算精度和计算效率.而且方法绕开对空间变量半离散化有限元格式的讨论,使得理论研究更简便.     

非饱和土壤水流问题的CN有限元格式

首先给出二维非饱和土壤水流问题基于Crank-Nicolson(CN)方法的具有时间二阶精度的半离散化格式,然后直接从CN时间半离散化格式出发,建立具有时间二阶精度的全离散化CN有限元格式,并给出误差估计,最后用数值例子说明全离散化CN有限元格式的优越性.这种方法可以绕开关于空间变量的半离散化格式的讨论,提高时间离散的精度,极大地减少时间方向的迭代步,从而减少实际计算中截断误差的积累,提高计算精度和计算效率.     

The Legendre Collocation Spectral Method for the Ground State Solutions of the Bose-Einstein Condensates北大核心CSCD

近年来,有关Bose-Einstein凝聚态基态解的实验研究已经取得了一系列重要的成果.该文在相关研究成果的基础上,首先通过降维和无量纲化方法将Bose-Einstein凝聚态基态解问题转换成能量泛函极值问题,在离散该泛函时,尝试使用Legendre配置谱方法离散该能量泛函的一维和二维情形.其次,对该能量泛函极小值问题进行了数值模拟.最后,通过分析实验数据结果和图像得出,针对非旋转的Bose-Einstein凝聚态的基态解问题可以使用Legendre配置谱方法来求解,且数值结果的误差较小.     

关于人口发展方程半离散算法的研究

利用半离散的方法将人口发展方程的边界条件进行离散,离散后得到两个相应的偏微分方程模型,然后利用算子半群的理论证明了离散后的解都逼近原方程的解,从而证明这种半离散方法是可行的.     

一维Maxwell方程间断解的多区域Legendre-Galerkin谱方法

研究了非一致介质一维Maxwell方程间断问题的多区间Legendre谱方法,建立了半离散和全离散Crank-Nicolson格式,设计了并行算法,并分析了方法的稳定性和收敛性.数值算例显示了多区间Legendre谱方法对于间断问题的有效性.     

凹角区域双曲外问题的精确人工边界条件

研究一类凹角区域双曲型外问题的数值方法.先用Newmark方法对时间进行离散化,在每个时间步求解一个椭圆外问题.然后引入人工边界,并获得精确的人工边界条件.给出半离散化问题的变分问题,证明了变分问题的适定性,并给出了误差估计.最后给出数值例子,以示该方法的可行性与有效性.     

Analysis of a moving collocation method for one-dimensional partial differential equations

A moving collocation method has been shown to be very efficient for the adaptive solution of second- and fourth-order time-dependent partial differential equations and forms the basis for the two robust codes MOVCOL and MOVCOL4.In this paper,the relations between the method and the traditional collocation and finite volume methods are investigated.It is shown that the moving collocation method inherits desirable properties of both methods: the ease of implementation and high-order convergence of the traditional collocation method and the mass conservation of the finite volume method.Convergence of the method in the maximum norm is proven for general linear two-point boundary value problems.Numerical results are given to demonstrate the convergence order of the method.     

Convergence analysis for a second-order elliptic equation by a collocation method using scattered points

A collocation method using scattered points applied to a second-order elliptic differential equation is analyzed by establishing a new quadrature formula for the space of the polynomials. We show that a polynomial solution possesses stability and preserves a similar convergence property occurred in the classical high order collocation method.     

Spectral homotopy analysis method and its convergence for solving a class of nonlinear optimal control problems

A combination of the hybrid spectral collocation technique and the homotopy analysis method is used to construct an iteration algorithm for solving a class of nonlinear optimal control problems (NOCPs). In fact, the nonlinear two-point boundary value problem (TPBVP), derived from the Pontryagin’s Maximum Principle (PMP), is solved by spectral homotopy analysis method (SHAM). For the first time, we present here a convergence proof for SHAM. We treat in detail Legendre collocation and Chebyshev collocation. It is indicated that Legendre collocation gives the same numerical results with Chebyshev collocation. Comparisons are made between SHAM, Matlab bvp4c generated results and results from literature such as homotopy perturbation method (HPM), optimal homotopy perturbation method (OHPM) and differential transformations.     

On a spline collocation method for a singularly perturbed problem

We consider a singularly perturbed boundary value problem with two small parameters. The problem is numerically treated by a quadratic spline collocation method. The suitable choice of collocation points provides the inverse monotonicity enabling utilization of barrier function method in the error analysis. Numerical results give justification of the proposed method. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An error analysis and the mesh independence principle for a nonlinear collocation problem

A nonlinear Dirichlet boundary value problem is approximated by an orthogonal spline collocation scheme using piecewise Hermite bicubic functions. Existence, local uniqueness, and error analysis of the collocation solution and convergence of Newton's method are studied. The mesh independence principle for the collocation problem is proved and used to develop an efficient multilevel solution method. Simple techniques are applied for estimating certain discretization and iteration constants that are used in the formulation of a mesh refinement strategy and an efficient multilevel method. Several mesh refinement strategies for solving a test problem are compared numerically. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

A hybrid collocation method for a nonlinear Volterra integral equation with weakly singular kernel

This work is concerned with the numerical solution of a nonlinear weakly singular Volterra integral equation. Owing to the singular behavior of the solution near the origin, the global convergence order of product integration and collocation methods is not optimal. In order to recover the optimal orders a hybrid collocation method is used which combines a non-polynomial approximation on the first subinterval followed by piecewise polynomial collocation on a graded mesh. Some numerical examples are presented which illustrate the theoretical results and the performance of the method. A comparison is made with the standard graded collocation method.     

The use of a stable boundary analogue of the collocation method to approximate the solution of some problems in mechanics

To solve boundary-value problems for elliptic equations, the boundary analogue of the method of least squares is replaced by a boundary analogue of the collocation method. The change is made using a discrete representation of the scalar product in the spaces of functions which, in the case of a smooth boundary, are integrable with their square over the boundary of the region, and of functions which, in the case of a piecewise-linear boundary, are integrable with their square, when weighted, over the boundary of the region. The method used to choose the collocation points which ensure the collocation method to be stable is justified for the case of the Dirichlet problem.     

On a least-squares collocation method for linear differential-algebraic equations

Summary The aim of this note is to extend some results on least-squares collocation methods and to prove the convergence of a least-squares collocation method applied to linear differential-algebraic equations. Some numerical examples are presented.     

Solution of a stochastic Darcy equation by polynomial chaos expansion

This paper deals with solving a boundary value problem for the Darcy equation with a random hydraulic conductivity field.We use an approach based on polynomial chaos expansion in a probability space of input data.We use a probabilistic collocation method to calculate the coefficients of the polynomial chaos expansion. The computational complexity of this algorithm is determined by the order of the polynomial chaos expansion and the number of terms in the Karhunen–Loève expansion. We calculate various Eulerian and Lagrangian statistical characteristics of the flow by the conventional Monte Carlo and probabilistic collocation methods. Our calculations show a significant advantage of the probabilistic collocation method over the directMonte Carlo algorithm.     

A numerical approach for solving initial-boundary value problem describing the process of cooling of a semi-infinite body by radiation

In this article, we have introduced a Taylor collocation method, which is based on collocation method for solving initial-boundary value problem describing the process of cooling of a semi-infinite body by radiation. This method is based on first taking the truncated Taylor expansions of the solution function in the fractional differential equation and then substituting their matrix forms into the equation. Using collocation points, we have the system of nonlinear algebraic equation. Then, we solve the system of nonlinear algebraic equation using Maple 13 and we have the coefficients of Taylor expansion. In addition, numerical results are presented to demonstrate the effectiveness of the proposed method.