解弹性力学第二类边界积分方程的求积法与分裂外推

利用Sidi奇异求积公式,提出了解曲边多角形域上线性弹性力学第二类边界积分方程的求积法,即离散矩阵的每个元素的生成只需赋值不需计算任何奇异积分.通过估计离散矩阵的特征值和利用Anselone聚紧收敛理论,证明了近似解的收敛性;同时得到了误差的多参数渐近展开式;通过并行地解粗网格上的离散方程,利用分裂外推获得了高精度近似解和后验误差.     

非线性演化方程孤立波的同伦分析法求解

利用同伦分析法求解了Burgers方程，得到了其扭结形孤立波的近似解析解，该解非常接近于相应的精确解.结果表明，同伦分析法可用来求解非线性演化方程的孤立波解.同时，也对所用方法进行了一定扩展，得到了Kadomtsev-Petviashvili(KP)方程的钟形孤立子解.经过扩展后的方法能够更方便地用于求解更多非线性演化方程的高精度近似解析解. 关键词： Burgers方程 同伦分析法 KP方程 孤立波解     

同伦分析法在求解非线性演化方程中的应用

利用同伦分析法求解了(2+1)维改进的 Zakharov-Kuznetsov方程, 得到了它的近似周期解,该解与精确解符合很好. 结果表明，同伦分析法在求解高维非线性演化方程时, 仍然是一种行之有效的方法. 同时,还对该方法进行了一定的扩展， 经过扩展后的方法能够更方便地求解更多非线性演化方程的高精度近似解析解. 关键词： 同伦分析法 改进的 Zakharov-Kuznetsov方程 周期解     

基于高精度解析解的航天器热设计方法

针对微小卫星热设计需求,基于球形双层集总参数模型,解析求解了双节点动态热平衡方程。提出使用多个类谐波近似处理周期性外热流边界的手段,消除外热流近似带来的计算误差。通过单次线性化处理,直接求解动态热平衡方程,得到高精度的动态温度解析解,使用四阶Runge-Kutta方法进行了数值验证。基于高精度解析解,设计了航天器导热–辐射双参数直接调控计算程序,对某小卫星模型进行热设计,对设计结果进行热仿真计算,验证了设计方案的可行性与准确性。     

激光陀螺锁区的谐波测量方法

锁区对速率偏频激光陀螺的性能有重要影响,而高精度的锁区测最依然是一个难题.从激光陀螺闭锁方程出发,通过理论近似的方法得到了锁区与激光陀螺输出信号谐波的近似关系.在此基础卜提出了一种激光陀螺锁区的谐波测量方法,进而通过闭锁方程的严格数值解进行了修正和误差分析.最后利用谐波测量法测量了某一激光陀螺的锁区大小,结果表明精度能够优于5%.     

Fisher方程的有界衰减振荡解

为了研究非线性发展方程的有界衰减振荡解,特选取Fisher方程为例. Fisher方程在描述激发介质的非数值模型(如Belousov-Zhabotinsky (BZ)反应)中, 其解的振幅取负值是有意义的.应用平面动力系统理论,研究了Fisher方程有界行波解存在的条件, 利用LS解法和线性化解法给出了其有界衰减振荡解的近似解析表达式,并进行了误差估计.     

Excel作图和二分法结合解超越方程的一种方法

利用Excel作图和二分法结合解超越方程,既能总体上把握解的概况,又能快速地得到各个解的高精度近似结果,是一种简单、实用的超越方程解法。     

光栅干涉位移传感器信息解算方法研究

面向光栅干涉式位移传感器信号高精度解算的需求,提出了一种基于椭圆拟合补偿及反正切算法的位移解算改进方法。采用傅里叶变换及标准化椭圆参数拟合法对信号误差估计,通过构建线性误差补偿模型对带误差信号进行正交补偿,结合反正切算法对信号进行线性化处理,实现高精度位移解算读出。Matlab仿真结果表明,该方法对信号正交补偿效果显著,线性化处理算法解算位移最大相对误差小于0.2%。为提高光栅干涉式传感器位移解算精度提供了有效途径和方法。     

迎风紧致格式求解Hamilton-Jacobi方程

基于Hamilton-Jacobi(H-J)方程和双曲型守恒律之间的关系,将三阶和五阶迎风紧致格式推广应用于求解H-J方程,建立了高精度的H-J方程求解方法.给出了一维和二维典型数值算例的计算结果,其中包括一个平面激波作用下的Richtmyer Meshkov界面不稳定性问题.数值试验表明,在解的光滑区域该方法具有高精度,而在导数不连续的不光滑区域也获得了比较好的分辨效果.相比于同阶精度的WENO格式,本方法具有更小的数值耗散,从而有利于多尺度复杂流动的模拟中H-J方程的求解.     

单摆系统的振动研究

通过求解变张力弦振动微分方程的边值问题,给出摆线与摆球的质量比为任意值时单摆系统运动的一般解和本征频率满足的方程.利用该方程求得高精度的单摆系统周期数值解,特别是拟合出单摆系统作基频振动时一个范围大、精度高的周期近似公式.同时将理论与实验进行比较,结果二者相符.     

A splitting extrapolation for solving nonlinear elliptic equations with d-quadratic finite elements

Nonlinear elliptic partial differential equations are important to many large scale engineering and science problems. For this kind of equations, this article discusses a splitting extrapolation which possesses a high order of accuracy, a high degree of parallelism, less computational complexity and more flexibility than Richardson extrapolation. According to the problems, some domain decompositions are constructed and some independent mesh parameters are designed. Multi-parameter asymptotic expansions are proved for the errors of approximations. Based on the expansions, splitting extrapolation formulas are developed to compute approximations with high order of accuracy on a globally fine grid. Because these formulas only require us to solve a set of smaller discrete subproblems on different coarser grids in parallel instead of on the globally fine grid, a large scale multidimensional problem is turned into a set of smaller discrete subproblems. Additionally, this method is efficient for solving interface problems.     

Asymptotic solution for a perturbed mechanism of western boundary undercurrents in the Pacific

This paper consider a class of perturbed mechanism for the western boundary undercurrents in the Pacific. The model of generalized governing equations is studied. Using the perturbation method, it constructs the asymptotic solution of the model. And the accuracy of asymptotic solution is proved by the theory of differential inequalities. Thus the uniformly valid asymptotic expansions of the solution are obtained.     

正交各向异性材料切口尖端热弹奇性特征分析

研究正交各向异性平面V形切口,计算其热弹奇性特征.通过引入切口尖端物理场的渐近级数展开式,将应力和热流平衡方程转化为关于奇性指数的特征常微分方程组,采用插值矩阵法求解,获取切口尖端的热流、应力奇性指数和对应的特征角函数.算例表明,该法精度高适应性强.     

双压电材料三维界面端部力电耦合场奇异性

基于切口根部物理场的幂级数渐近展开假设,从三维应力平衡方程和麦克斯韦方程组出发,导出关于双压电材料楔形界面切口端部奇性指数的特征微分方程组,并将切口的力电学边界条件表达为奇性指数和特征角函数的组合,从而将双压电材料楔形界面切口端部奇性指数的计算转化为相应边界条件下常微分方程组特征值的求解,运用插值矩阵法求解界面端部若干阶应力奇性指数和相应特征函数.计算结果与已有结果对比表明本文方法的有效性和具有较高的计算精度.     

Two-body solution of the Newman-Penrose asymptotic equations

The Newman-Penrose nonlinear asymptotic field equations are separated in terms of spin weight spherical harmonics (s.w.s.h.). As an example, the results are used to study the radiation effects on a two-body system. The presence of radiation is manifest through the nonlinear terms in the asymptotic equations. If these terms are assumed to be small, the asymptotic equations can be formally solved by an iteration procedure. For the above example the first step of the iteration procedure is implemented to an accuracy that includes the effects of radiation up to octopole order. The results illustrate the usual internal decay of the orbit as well as an acceleration of the system's center of mass. In favorable cases, the two-body source can reach significant velocities due to the radiation reaction.     

Artificial compressibility method revisited: Asymptotic numerical method for incompressible Navier–Stokes equations

The artificial compressibility method for the incompressible Navier–Stokes equations is revived as a high order accurate numerical method (fourth order in space and second order in time). Similar to the lattice Boltzmann method, the mesh spacing is linked to the Mach number. An accuracy higher than that of the lattice Boltzmann method is achieved by exploiting the asymptotic behavior of the solution of the artificial compressibility equations for small Mach numbers and the simple lattice structure. An easy method for accelerating the decay of acoustic waves, which deteriorate the quality of the numerical solution, and a simple cure for the checkerboard instability are proposed. The high performance of the scheme is demonstrated not only for the periodic boundary condition but also for the Dirichlet-type boundary condition.     

Self-similar factor approximants for evolution equations and boundary-value problems

The method of self-similar factor approximants is shown to be very convenient for solving different evolution equations and boundary-value problems typical of physical applications. The method is general and simple, being a straightforward two-step procedure. First, the solution to an equation is represented as an asymptotic series in powers of a variable. Second, the series are summed by means of the self-similar factor approximants. The obtained expressions provide highly accurate approximate solutions to the considered equations. In some cases, it is even possible to reconstruct exact solutions for the whole region of variables, starting from asymptotic series for small variables. This can become possible even when the solution is a transcendental function. The method is shown to be more simple and accurate than different variants of perturbation theory with respect to small parameters, being applicable even when these parameters are large. The generality and accuracy of the method are illustrated by a number of evolution equations as well as boundary value problems.     

A low-Mach number fix for Roe’s approximate Riemann solver

We present a low-Mach number fix for Roe’s approximate Riemann solver (LMRoe). As the Mach number Ma tends to zero, solutions to the Euler equations converge to solutions of the incompressible equations. Yet, standard upwind schemes do not reproduce this convergence: the artificial viscosity grows like 1/Ma, leading to a loss of accuracy as Ma → 0. With a discrete asymptotic analysis of the Roe scheme we identify the responsible term: the jump in the normal velocity component ΔU of the Riemann problem. The remedy consists of reducing this term by one order of magnitude in terms of the Mach number. This is achieved by simply multiplying ΔU with the local Mach number. With an asymptotic analysis it is shown that all discrepancies between continuous and discrete asymptotics disappear, while, at the same time, checkerboard modes are suppressed. Low Mach number test cases show, first, that the accuracy of LMRoe is independent of the Mach number, second, that the solution converges to the incompressible limit for Ma → 0 on a fixed mesh and, finally, that the new scheme does not produce pressure checkerboard modes. High speed test cases demonstrate the fall back of the new scheme to the classical Roe scheme at moderate and high Mach numbers.     

Calculation of the linear kernel of the collision integral for the hard-sphere potential

The expansion of a distribution function in spherical harmonics transforms the Boltzmann equation into a system of integro-differential equations with kernels depending only of the magnitudes of velocities. The kernels can be expressed in terms of the sums including the matrix elements (MEs) of the collision integral. The kernels are constructed using new results of ME calculations; analysis of errors is carried out with the help of analytic expressions for kernels, which were derived by Hilbert and Hecke for the hard-sphere model. The concept of generalized matrix elements is introduced and their asymptotic representation is constructed for large values of indices. Analytic expressions for the contribution from MEs with large indices to the kernels are constructed. The high accuracy of the construction of a kernel using MEs is demonstrated.     

On the equivalence of convergent kinetic equations for hot dilute plasmas: II. Generating functions for collision brackets

The generating functions for the collision brackets associated with two alternative convergent kinetic equations are derived for small values of the plasma parameter. It is shown that the first few terms in the asymptotic expansions of these generating functions are identical. Consequently, both kinetic equations give rise to the same transport coefficients in arbitrarily high order of the Chapman-Cowling truncation scheme.