改进的tanh函数方法与广义变系数KdV和MKdV方程新的精确解

利用改进的tanh函数方法将广义变系数KdV方程和MKdV方程化为一阶变系数非线性常微分方 程组-通过求解这个变系数非线性常微分方程组，获得了广义变系数KdV方程和MKdV方程新的 精确类孤子解、有理形式函数解和三角函数解- 关键词： 改进的tanh函数方法 类孤子解 有理形式函数解 三角函数解     

随机扰动对拟小波方法求解对流扩散方程的影响

引进拟小波方法数值求解对流扩散方程，研究结果表明，计算带宽W有一个极值，当计算带宽W取该极值时，该方程的拟小波解的精度最高，且好于迎风格式。当边界发生随机不等幅扰动时，对于积分时间较长的情况，拟小波格式的效果要稍逊于迎风格式；当边界发生随机等幅扰动时，若计算带宽W取大于等于20的整数时，方程拟小波解的精度与迎风格式相同；当参数受到随机扰动时，W取10时的拟小波解的均方根误差要小于迎风格式；在初值发生随机扰动且计算带宽W取10时，方程的拟小波解的精度最高，好于迎风格式。     

MKdV方程的微扰理论

求得了含微扰的MKdV方程u_t+6u²u_x+u_xxx=εR［u］一级近似解的一般表达式,以及一个特例的具体表达式. 关键词：     

求解Euler方程的隐式无网格算法

研究了求解Eluer方程的稳式无网格算法，用点云离散计算区域，代替通常的网格划分；在当地点云上，引入二次平方极小曲面逼近计算空间导数，用Roe的近似Riemann解确定通量；并用LU-SGS算法求解离散得到的Euler方程稳式时间后差联立方程组，数值模拟了二维翼型跨音速绕流，由于无网格算法区域离散只涉及点云，具有灵活性，适合处理复杂的气动外形。     

高阶Schr(o)dinger方程的高精度辛格式

提出由第三类生成函数法构造高阶Schroedinger方程δu/δt=i(-1)^nδ^2mu/δx^2m的高精度辛格式．首先，给出它的典则Hamilton方程组；然后，成功地克服了本质上是困难的高阶变分导数的计算，并利用第三类生成函数法得到在时间方向具有任意阶精度的半离散方程，进而得到原始方程相关的修正方程的离散形式，最后得到各种精度的辛格式．数值结果表明该格式是有效的，具有高精度及良好的长时间数值行为等特性．     

球形Dirac方程的空间格点求解及假态问题

本文在空间格点上利用虚时间步长方法求解了球形Dirac方程, 着重研究了出现的假态问题. 利用三点数值导数公式离散方程中一阶导数项, 可以证明对于量子数为 κ 和 -κ的单粒子能级能量是完全相同的, 其中一个为物理解, 另一个为假态. 通过在径向Dirac方程中引入Wilson 项, 可以解决假态问题, 得到全部物理解. 文章以 Woods-Saxon 势为例, 考虑 Wilson 项后, 得到与打靶法一致的结果.     

时间分数阶Edwards-Wilkinson方程的标度行为研究

分别采用数值模拟和标度分析方法对1+1维时间分数阶Edwards-Wilkinson方程的标度行为进行研究.利用Caputo分数阶导数数值解求得的生长指数与采用直接标度分析方法得到的结果一致.     

一维双曲守恒律的龙格-库塔控制体积间断有限元方法

给出数值求解一维双曲守恒律方程的新方法——龙格-库塔控制体积间断有限元方法(RKCVDFEM),其中空间离散基于控制体积有限元方法,时间离散基于二阶TVB Runge-Kutta技术,有限元空间选取为分段线性函数空间.理论分析表明,格式具有总变差有界(TVB)的性质,而且空间和时间离散形式上具有二阶精度.数值算例表明,数值解收敛到熵解并且对光滑解的收敛阶是最优的,优于龙格-库塔间断Galerkin方法(RKDGM)的计算结果.     

高阶Schrdinger方程的高精度辛格式

提出由第三类生成函数法构造高阶Schr dinger方程ut=i(-1)m2mux2m的高精度辛格式.首先,给出它的典则Hamilton方程组;然后,成功地克服了本质上是困难的高阶变分导数的计算,并利用第三类生成函数法得到在时间方向具有任意阶精度的半离散方程,进而得到原始方程相关的修正方程的离散形式,最后得到各种精度的辛格式.数值结果表明该格式是有效的,具有高精度及良好的长时间数值行为等特性.     

求解Kuramoto-Sivashinsky方程的平移基无单元Galerkin方法

Kuramoto-Sivashinsky方程是一种可以描述复杂混沌现象的高阶非线性演化方程.方程中高阶导数项的存在,使得传统无单元Galerkin方法采用高次多项式基函数构造形函数时,形函数违背了一致性条件.因此,本文提出了一种采用平移多项式基函数的无单元Galerkin方法.与传统无单元Galerkin方法相比,该方法在方程离散时依然采用Galerkin进行离散,但形函数的构造采用了基于平移多项式基函数的移动最小二乘近似.通过对具有行波解和混沌现象的Kuramoto-Sivashinsky方程的数值模拟,验证了本文方法的有效性.     

求解非线性偏微分方程的自适应小波精细积分法

以Burgers方程为例,提出了一种求解偏微分方程的自适应多层插值小波配置法,通过引入一种具有插值特性的拟Shannon小波并利用插值小波理论构造了多层自适应插值小波算子,从而在空间实现了偏微分方程的自适应离散.另外,精细时程积分方法和外推法的引入不但有助于提高求解速度和数值结果的精度,而且使时间积分步长的选取具有了自适应性.     

Numerical simulation of the soliton solutions for a complex modified Korteweg—de Vries equation by a finite difference method

In this paper,a Crank-Nicolson-type finite difference method is proposed for computing the soliton solutions of a complex modifed Korteweg de Vries(MKdV)equation(which is equivalent to the Sasa-Satsuma equation)with the vanishing boundary condition.It is proved that such a numerical scheme has the second order accuracy both in space and time,and conserves the mass in the discrete level.Meanwhile,the resuling scheme is shown to be unconditionally stable via the von Nuemann analysis.In addition,an iterative method and the Thomas algorithm are used together to enhance the computational efficiency.In numerical experiments,this method is used to simulate the single-soliton propagation and two-soliton collisions in the complex MKdV equation.The numerical accuracy,mass conservation and linear stability are tested to assess the scheme's performance.     

An improved Hirota bilinear method and new application for a nonlocal integrable complex modified Korteweg-de Vries (MKdV) equation

The Hirota bilinear method has been studied in a lot of local equations, but there are few of works to solve nonlocal equations by Hirota bilinear method. In this letter, we show that the nonlocal integrable complex modified Korteweg-de Vries (MKdV) equation admits multiple complex soliton solutions. A variety of exact solutions including the single bright soliton solutions and two bright soliton solutions are derived via constructing an improved Hirota bilinear method for nonlocal complex MKdV equation. From the gauge equivalence, we can see the difference between the solution of nonlocal integrable complex MKdV equation and the solution of local complex MKdV equation.     

The element-free Galerkin method of numerically solving a regularized long-wave equation

The element-free Galerkin (EFG) method is used in this paper to find the numerical solution to a regularized long-wave (RLW) equation. The Galerkin weak form is adopted to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. The effectiveness of the EFG method of solving the RLW equation is investigated by two numerical examples in this paper.     

Three-dimensional finite-difference lattice Boltzmann model and its application to inviscid compressible flows with shock waves

In this paper, a three-dimensional (3D) finite-difference lattice Boltzmann model for simulating compressible flows with shock waves is developed in the framework of the double-distribution-function approach. In the model, a density distribution function is adopted to model the flow field, while a total energy distribution function is adopted to model the temperature field. The discrete equilibrium density and total energy distribution functions are derived from the Hermite expansions of the continuous equilibrium distribution functions. The discrete velocity set is obtained by choosing the abscissae of a suitable Gauss–Hermite quadrature with sufficient accuracy. In order to capture the shock waves in compressible flows and improve the numerical accuracy and stability, an implicit–explicit finite-difference numerical technique based on the total variation diminishing flux limitation is introduced to solve the discrete kinetic equations. The model is tested by numerical simulations of some typical compressible flows with shock waves ranging from 1D to 3D. The numerical results are found to be in good agreement with the analytical solutions and/or other numerical results reported in the literature.     

非线性偏微分方程边值问题的优化算法研究与应用

针对椭圆类非线性偏微分方程边值问题,以差分法和动态设计变量优化算法为基础,以离散网格点未知函数值为设计变量,以离散网格点的差分方程组构建为复杂程式化形式的目标函数.提出一种求解离散网格点处未知函数值的优化算法.编制了求解未知离散点函数值的通用程序.求解了具体算例.通过与解析解对比,表明了本文提出求解算法的有效性和精确性,将为更复杂工程问题分析提供良好的解决方法. 关键词： 非线性偏微分方程 边值问题 动态设计变量优化算法 程序设计     

An element-free Galerkin(EFG) method for numerical solution of the coupled Schrdinger-KdV equations

The present paper deals with the numerical solution of the coupled Schrdinger-KdV equations using the elementfree Galerkin(EFG) method which is based on the moving least-square approximation.Instead of traditional mesh oriented methods such as the finite difference method(FDM) and the finite element method(FEM),this method needs only scattered nodes in the domain.For this scheme,a variational method is used to obtain discrete equations and the essential boundary conditions are enforced by the penalty method.In numerical experiments,the results are presented and compared with the findings of the finite element method,the radial basis functions method,and an analytical solution to confirm the good accuracy of the presented scheme.     

Meshless analysis of three-dimensional steady-state heat conduction problems

Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson’s equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method (RKPM). A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.     

A Priori Error Estimates of Crank-Nicolson Finite Volume Element Method for Parabolic Optimal Control Problems

In this paper, the Crank-Nicolson linear finite volume element method is applied to solve the distributed optimal control problems governed by a parabolic equation. The optimal convergent order $\mathcal{O}(h^2+k^2)$ is obtained for the numerical solution in a discrete $L^2$-norm. A numerical experiment is presented to test the theoretical result.