��ά�߽�CIP��ֵ��������Ӧ��

利用函数及其高阶导数值构造五次插值函数近似网格单元内的真实解,改进数值求解双曲类偏微分方程的CIP数值算法。基于之前的一维高阶CIP数值算法思想,不同于利用时间分裂技术,发展了二维高阶CIP数值算法。改进后的算法具有五阶数值精度和显示格式的优点。     

理想弹塑性固体Riemann问题的解结构与初值条件

针对理想弹塑性固体材料的一维Riemann问题，在不考虑真空的情况下，讨论其所有可能存在的解结构，给出每一种解结构下对应的初值条件且证明该系列初值条件的完备性，即任意给定的物理量初值均有且只有一种解结构与之对应.基于该理论，在设计精确或近似理想弹塑性Riemann问题求解器时，可以依据初值条件对任意物理量初值直接判断其对应的解结构，从而提高求解器的精度和效率.数值试验验证了该系列初值条件的正确性和有效性.     

一个新的求解点堆中子动力学方程组的数值方法

使用中子密度一阶泰劳多项式分段近似技术,给出一个新的求解点堆中子动力学方程组的数值方法并采用全隐格式以克服方程组的刚性,同时确保解的必要精度。数值结果表明:在隐式一阶多项式近似下,对合适的反应性输入能够取得足够精确的结果。当反应性给定时,对于求解反应堆动力学问题,能给出一个简法的计算过程。     

高温等离子体相对论性平均原子模型中的半解析波函数

半解析求解平均原子模型方法充分利用了已知精确波函数的解析性质,通过对平均原子模型中势函数的数值拟合,就得到仅含一个数值因子的半解析波函数以及相应的能量本征值.本文列出了等离子体中相对论性平均原子模型的诸方程,特别注意方程求解技术和程序设计中的一些细节.与完全数值解以及其他类似模型得到的数值解进行的比较表明,在较高温度条件下半解析结果的精度是相当高的,求解的效率也很高.此外还对物理模型中某些缺陷进行了分析.     

强阻尼广义sine-Gordon方程特征问题的变分迭代法

研究了在数学、力学中广泛出现的一类非线性强阻尼广义sine-Gordon扰动微分方程问题. 首先, 引入行波变换, 求出退化方程的精确解. 再构造一个泛函, 创建了一个变分迭代算法, 最后, 求出原非线性强阻尼广义sine-Gordon扰动微分方程问题的近似行波解析解. 用变分迭代法可得到的各次近似解, 具有便于求解、精度高等特点. 求得的近似解析解弥补了单纯用数值方法的模拟解的不足.     

第二类边界积分方程Nystrom解的高精度组合方法

第二类边界积分方程常用配置法或Ｇａｌｅｒｋｉｎ法计算，主要困难有：计算积分耗去大量机时；离散方程是满阵且不对称，计算量随剖分精细而急剧增加。本文提出Ｎｙｓｔｒｏｍ近似解的高精度组合法能有效克服上述困难。组合方法是并行地解ｍ个具有ｎ个不同结点的方程组，对得到的ｍ个内点值取算术平均就得到了组合近似，本文证明组合近似精度几乎与解ｍｎ个结点近似方程达到精度同阶，数值结果表明本文方法简单、有效、并且算法高度     

摄动有限体积法重构近似高精度的意义

研讨有限体积(FV)方法重构近似高精度的作用问题.FV方法中积分近似采用中点规则为二阶精度时,重构近似高精度(精度高于二阶)的意义和作用是一个有争议的问题.利用数值摄动技术[1,2]构造了标量输运方程的积分近似为二阶精度、重构近似为任意阶精度的迎风型和中心型摄动有限体积(PFV)格式.迎风PFV格式无条件满足对流有界准则(CBC),中心型PFV格式为正型格式,两者均不会产生数值振荡解.利用PFV格式求解模型方程的数值结果表明:与一阶迎风和二阶中心格式相比,PFV格式精度高、对解的间断分辨率高、稳定性好、雷诺数的适用范围大,数值地"证实"重构近似高精度和PFV格式的实际意义和好处.     

基于拉格朗日方程的含源导热问题分析

应用基于热质理论的拉格朗日方程对含内热源一维瞬态导热问题进行了近似求解,计算结果与精确解吻合较好;同时指出了现有文献中采用分析力学方法讨论导热问题时存在的某些不足.本文工作表明了基于热质理论用分析力学观点研究导热间题这一方法的有效性,也从一个侧面说明热质和热质能等热学新概念具有一定的合理性.     

单模光纤色散的解析形式

根据波导标量解本征值方程及其递推关系,提出一种利用Gloge关系求解单模光纤中波导色散的理论方法,给出了色散的解析形式,通过分析归一化传输常数的近似解与精确解间的差别论证了这种解析法具有精确求解的计算精度,给出普通单模光纤（G．652）光纤色散的实验数据,并与计算的色散解析解曲线加以比较,二者达到极好的吻合,利用所得到的结果,分析了数值微分法和经验公式的计算精度。     

粘性不可压流的变分多尺度数值模拟

在变分多尺度的理论框架内,将待求解的各个物理量分解到"粗"、"细"两种尺度上.在"细"尺度上采用"泡"函数作为近似函数,通过Petrov-Galerkin方法得到"细"尺度上的近似解;然后引入求解"粗"尺度方程所需的稳定项及与其相适应的稳定化因子;最后运用有限元方法求解"粗"、"细"两种尺度耦合的整体变分多尺度方程,得到有限元近似解.数值算例表明,该处理方法成功地消除了数值求解粘性不可压Navier-Stokes方程过程中,由对流占优和速度-压力失耦引起的数值伪振荡;所引入的稳定化因子适用于结构网格及非结构网格上的数值计算.     

An application of the K-integrals for solving the radiative transfer equation with Fresnel boundary conditions

We introduce the reader to an approximate method of solving the transport equation which was developed in the context of neutron thermalisation by Kladnik and Kuscer in 1962 [Kladnik R, Kuscer I. Velocity dependent Milne's problem. Nucl Sci Eng 1962;13:149]. Essentially the method is based upon two special weighted integrals of the one-dimensional transport equation which are valid regardless of the boundary conditions, and any solution must satisfy these integral relationships which are called the K-integrals. To obtain an approximate solution to the transport equation we turn the argument around and insist that any approximate solution must also satisfy the K-integrals. These integrals are particularly useful when the problem under consideration cannot be solved easily by analytic methods. It also has the marked advantage of being applicable to problems where there is energy exchange in a collision and anisotropy of scattering. To establish the feasibility of the method we obtain a number of approximate solutions using the K-integral method for problems to which we have exact analytical solutions. This enables us to validate the method. It is then applied to a new problem that has not yet been solved; namely the calculation of the discontinuity in the scalar intensity at the boundary between two optically dissimilar materials.     

A new approach with piecewise-constant arguments to approximate and numerical solutions of oscillatory problems

This paper is devoted to the development of a novel approximate and numerical method for the solutions of linear and non-linear oscillatory systems, which are common in engineering dynamics. The original physical information included in the governing equations of motion is mostly transferred into the approximate and numerical solutions. Therefore, the approximate and numerical solutions generated by the present method reflect more accurately the characteristics of the motion of the systems. Furthermore, the solutions derived are continuous everywhere with good accuracy and convergence in comparing with Runge-Kutta method. An approximate solution is developed for a linear oscillatory problem and compared with its corresponding exact solution. A non-linear oscillatory problem is also solved numerically and compared with the solutions of Runge-Kutta method. Both the graphical and numerical comparisons are provided in the paper. The accuracy of the approximate and numerical solutions can be controlled as desired by the number of terms in the Taylor series and the value of a single parameter used in the present work. Formulae for numerical computation in solving various linear and non-linear oscillatory problems by the new approach are provided in the paper.     

Ion collection by oblique surfaces of an object in a transversely flowing strongly magnetized plasma

The equations governing a collisionless obliquely flowing plasma around an ion-absorbing object in a strong magnetic field are shown to have an exact analytic solution even for an arbitrary (two-dimensional) object shape, when temperature is uniform, and diffusive transport can be ignored. The solution has an extremely simple geometric embodiment. It shows that the ion collection flux density to a convex body's surface depends only upon the orientation of the surface and provides the theoretical justification and calibration of oblique "Mach probes." The exponential form of this exact solution helps explain the approximate fit of this function to previous numerical solutions.     

Approximate analytic solutions for a generalized Hirota—Satsuma coupled KdV equation and a coupled mKdV equation

<正>This paper applies the variational iteration method to obtain approximate analytic solutions of a generalized Hirota-Satsuma coupled Korteweg-de Vries(KdV) equation and a coupled modified Korteweg-de Vries(mKdV) equation. This method provides a sequence of functions which converges to the exact solution of the problem and is based on the use of the Lagrange multiplier for the identification of optimal values of parameters in a functional.Some examples are given to demonstrate the reliability and convenience of the method and comparisons are made with the exact solutions.     

An Analytical Technique,Based on Natural Transform to Solve Fractional-Order Parabolic Equations

This research article is dedicated to solving fractional-order parabolic equations using an innovative analytical technique. The Adomian decomposition method is well supported by natural transform to establish closed form solutions for targeted problems. The procedure is simple, attractive and is preferred over other methods because it provides a closed form solution for the given problems. The solution graphs are plotted for both integer and fractional-order, which shows that the obtained results are in good contact with the exact solution of the problems. It is also observed that the solution of fractional-order problems are convergent to the solution of integer-order problem. In conclusion, the current technique is an accurate and straightforward approximate method that can be applied to solve other fractional-order partial differential equations.     

Phase reduction models for improving the accuracy of the finite element solution of time-harmonic scattering problems I: General approach and low-order models

This paper introduces a new formulation of high frequency time-harmonic scattering problems in view of a numerical finite element solution. It is well-known that pollution error causes inaccuracies in the finite element solution of short-wave problems. To partially avoid this precision problem, the strategy proposed here consists in firstly numerically computing at a low cost an approximate phase of the exact solution through asymptotic propagative models. Secondly, using this approximate phase, a slowly varying unknown envelope is introduced and is computed using coarser mesh grids. The global procedure is called Phase Reduction. In this first paper, the general theoretical procedure is developed and low-order propagative models are numerically investigated in detail. Improved solutions based on higher order models are discussed showing the potential of the method for further developments.     

Numerical soliton-like solutions of the potential Kadomtsev–Petviashvili equation by the decomposition method

In this Letter we present an Adomian's decomposition method (shortly ADM) for obtaining the numerical soliton-like solutions of the potential Kadomtsev–Petviashvili (shortly PKP) equation. We will prove the convergence of the ADM. We obtain the exact and numerical solitary-wave solutions of the PKP equation for certain initial conditions. Then ADM yields the analytic approximate solution with fast convergence rate and high accuracy through previous works. The numerical solutions are compared with the known analytical solutions.     

Toward Analytic Solution of Nonlinear Differential Difference Equations via Extended Sensitivity Approach

In this paper an efficient computational method based on extending the sensitivity approach(SA) is proposed to find an analytic exact solution of nonlinear differential difference equations.In this manner we avoid solving the nonlinear problem directly.By extension of sensitivity approach for differential difference equations(DDEs),the nonlinear original problem is transformed into infinite linear differential difference equations,which should be solved in a recursive manner.Then the exact solution is determined in the form of infinite terms series and by intercepting series an approximate solution is obtained.Numerical examples are employed to show the effectiveness of the proposed approach.     

Asymmetric and axisymmetric vibrations of finite transversely isotropic circular cylinders

This paper presents an approach for obtaining the exact frequency equations of axisymmetric and asymmetric free vibrations of transversely isotropic circular cylinders. The solution method is based on the three dimensional theory of linear elasticity and uses potential functions. Using this approach, the frequency spectra and vibration mode shapes are plotted for a number of transversely isotropic cylinders. The proposed approach introduces a number of merits compared to earlier approximate and exact solution methods. First, unlike numerically complicated series methods that provide approximate solutions, the proposed approach is exact. Second, combination of scalar functions employed for representing the displacement field is consistent with the physics of the problem. One scalar potential function has been considered for each component of the wave field inside the elastic cylinder. As a result, the solution is systematically divided into coupled and decoupled equations. In addition, by using this approach, there is no need to guess the final of the solution a priori. These merits make the proposed approach suitable for other vibration problems of anisotropic materials.