A Collocation Code for Singular Boundary Value Problems in Ordinary Differential Equations

We present a MATLAB package for boundary value problems in ordinary differential equations. Our aim is the efficient numerical solution of systems of ODEs with a singularity of the first kind, but the solver can also be used for regular problems. The basic solution is computed using collocation methods and a new, efficient estimate of the global error is used for adaptive mesh selection. Here, we analyze some of the numerical aspects relevant for the implementation, describe measures to increase the efficiency of the code and compare its performance with the performance of established standard codes for boundary value problems.     

A multistep collocation method for solving advection equations with delay

Advection equations with delay are appeared in the modeling of the dynamics of structured cell populations. In this article, we construct an efficient two-dimensional multistep collocation method for the numerical solution of a class of advection equations with delay. Equations with aftereffect and equations with both aftereffect and retardation of a state variable are considered. Computability of the algorithm and convergence properties of the proposed numerical method are analyzed for solutions in appropriate Sobolev spaces, and it is shown that the proposed scheme enjoys the spectral accuracy. Numerical examples are given and comparison with other existing methods in the literature is made to demonstrate the efficiency, superiority and high accuracy of the presented method.     

A numerical technique based on B-spline for a class of time-fractional diffusion equation

This paper presents an efficient numerical technique for solving a class of time-fractional diffusion equation. The time-fractional derivative is described in the Caputo form. The L1 scheme is used for discretization of Caputo fractional derivative and a collocation approach based on sextic B-spline basis function is employed for discretization of space variable. The unconditional stability of the fully-discrete scheme is analyzed. Two numerical examples are considered to demonstrate the accuracy and applicability of our scheme. The proposed scheme is shown to be sixth order accuracy with respect to space variable and (2 − α)-th order accuracy with respect to time variable, where α is the order of temporal fractional derivative. The numerical results obtained are compared with other existing numerical methods to justify the advantage of present method. The CPU time for the proposed scheme is provided.     

Computational aeroacoustics using the B-spline collocation method

One of the major problems in computational aero-acoustics is the disparity in length scales between the flow field and the acoustic field. As a result, a mapping function is normally used to achieve a non-uniform grid distribution. In this paper, a B-spline collocation method with an arbitrary grid placement capability is proposed. This capability not only allows an optimum grid distribution but also avoids the numerical complexities associated with the mapping function. The B-spline collocation method is applied to the case of spinning co-rotating vortices. The result agrees well with the matched asymptotic solution. To cite this article: R. Widjaja et al., C. R. Mecanique 333 (2005).     

求解Helmholtz方程基于核重构思想的最小二乘配点法

基于核重构思想构造近似函数，将配点法和最小二乘原理相结合对微分方程进行离散， 建立了Helmholtz方程的最小二乘配点格式，并分别研究了Helmholtz方程的波传播问题和 边界层问题. 通过数值算例可以发现，给出的数值计算结果非常接近于精确解，计算精度明显高于SPH 法的数值结果，且随着节点数目的增加，其精确度越来越高，具有良好的收敛性.     

一类弱奇性Volterra积分微分方程的级数展开数值解法

本文考察了一类弱奇性积分微分方程的级数展开数值解法,并给出了相应的收敛性分析.理论分析结果表明,若用已知函数的谱配置多项式逼近已知函数,那么方程的数值解以谱精度逼近方程的真解.数值实验数据也验证了这一理论分析结果.     

建筑结构最优阻尼及最优配置研究

利用阻尼来抗震是建筑物最简单的抗震方法,当在阻尼数目少于楼层数时,为使所布置的阻尼达到最优的抗震效果,本文开发了两种最优阻尼配置方法。第一种方法是根据临界激励法确定计算所有可能布置位置中的临界阻尼,选择最小的临界阻尼位置作为最优配置。第二种方法则是根据给定阻尼,计算最优控制效果以确定最优阻尼配置。然后根据实代码遗传算法,在结构抗震效果基本保证的前提下,寻找最优阻尼以降低结构总的阻尼增量。最后给出了一个计算实例。     

A Spectral Method for the Electrohydrodynamic Flow in a Circular Cylindrical Conduit

This paper presents a combination of the hybrid spectral collocation technique and the spectral homotopy analysis method(SHAM for short) for solving the nonlinear boundary value problem(BVP for short) for the electrohydrodynamic flow of a fluid in an ion drag configuration in a circular cylindrical conduit. The accuracy of the present solution is found to be in excellent agreement with the previously published solution. The authors use an averaged residual error to find the optimal convergence-control parameters. Comparisons are made between SHAM generated results, results from literature and Matlab ode45 generated results, and good agreement is observed.     

Discrete-time orthogonal spline collocation methods for the nonlinear Schrödinger equation with wave operator

In this paper, discrete-time orthogonal spline collocation schemes are proposed for the nonlinear Schrödinger equation with wave operator. These schemes are constructed by using orthogonal spline collocation approaches combined with finite difference methods. The conservative property, the convergence, and the stability of these methods are theoretically analyzed and also verified by extensive numerical experiments. In addition, some interesting phenomena which require further theoretical analysis are discussed numerically.     

Mass and momentum conservation of the least-squares spectral collocation method for the Navier-Stokes equations

From the literature, it is known that the Least-Squares Spectral Element Method (LSSEM) for the stationary Stokes equations performs poorly with respect to mass conservation but compensates this lack by a superior conservation of momentum. Furthermore, it is known that the Least-Squares Spectral Collocation Method (LSSCM) leads to superior conservation of mass and momentum for the stationary Stokes equations. In the present paper, we consider mass and momentum conservation of the LSSCM for time-dependent Stokes and Navier-Stokes equations. We observe that the LSSCM leads to improved conservation of mass (and momentum) for these problems. Furthermore, the LSSCM leads to the well-known time-dependent profiles for the velocity and the pressure profiles. To obtain these results, we use only a few elements, each with high polynomial degree, avoid normal equations for solving the overdetermined linear systems of equations and introduce the Clenshaw-Curtis quadrature rule for imposing the average pressure to be zero. Furthermore, we combined the transformation of Gordon and Hall (transfinite mapping) with the least-squares spectral collocation scheme to discretize the internal flow problems.