耦合界面计算格式的一维流体力学Lagrange有限点方法

为探索高维多介质流体力学散乱点集上的Lagrange有限点方法,首先对相应一维问题进行研究,提出一种Lagrange有限点方法:在计算区域内(包括物质界面)设置任意离散点集,所有力学量都设在该点集上,在内点和界面点上分别建立离散格式.内点算法为基于Taylor展开的差分方法.界面点算法为显式追踪算法,从定解条件出发,利用Rankine-Hugoniot关系和特征差分方法,计算界面点位置及相应的状态量变化.通过追踪界面点的运动得到物质界面是方法的最大特色.典型算例计算结果与精确解符合很好,验证了算法的合理和有效性.     

Generalized density-functional theory: extended weighted density approaches

A third-order density-functional theory is introduced by an approach that may be used to find density-functional theories to any higher-order accuracy provided only that known homogeneous state correlation functions are utilized as input. It is constructed from the required knowledge of a single weight function at each order. By way of application results are presented for the melting of classical hard spheres using functionals accurate to a third-order functional Taylor series in the homogeneous limit. Within the framework of the modified weighted density approximation, there is a uniform improvement in the solid phase-free energies, pressures and melting parameters, and further improvement also results when these functionals are optimized in a way that utilizes the close packing limit. The sensitivity of the results to existing and proposed models of the third-order direct correlation function is discussed.     

应用遗传算法和序列锥规划进行圆阵的方向图综合

为了解决圆阵方向图旁瓣电平较高的问题,提出了一种圆阵方向图综合方法。该方法将圆阵方向图的峰值旁瓣电平作为目标函数,首先采用改进遗传算法对阵元位置和阵元权值进行联合优化,不仅避免了算法的早熟收敛,而且符合理论意义上的全局最优。其次将遗传算法的优化结果作为初始迭代点,在其附近利用一阶泰勒级数将非凸的圆阵方向图综合问题转化为序列锥规划问题,以便采用凸优化理论进行高效求解。由于该算法增加了寻优操作的后期变异能力,因而有效提高了优化性能。最后仿真表明,阵元数一定的情况下,算法在进一步降低峰值旁瓣电平的同时,可有效减小其动态变化范围,使得圆阵方向图综合性能更优。      

Burton-Miller改进型边界方程的多频计算方法

提出了一种采用Burton-Miller改进型边界积分方程进行多频计算的方法。将Burton-Miller方程中的高奇异积分转化为弱奇异积分形式,获得Burton-Miller改进型边界积分方程;将方程中格林函数进行Taylor级数展开,并把波数从方程中分离出来,从而使随波数变化的计算矩阵表示为波数的矩阵级数形式。数值分析表明,本方法不仅保证了解在全波数范围内的唯一性,并且计算频率点数较多时可以节约大量时间,提高计算效率。      

计算半导体平板微腔自发发射的近似方法

引入新变量,并利用高阶泰勒展开完成半导体平板微腔自发发射的空间积分,由此得到半导体平板微腔TE模式自发发射的近似表达式.在腔长为半个中心波长和高反射率腔面的半导体平板微腔中,结合电子和空穴的费米分布函数,用近似方法计算垂直方向小角度内自发发射谱和总的自发发射谱,分别与数值空间积分基本相同,可以用于计算量子阱平板微腔自发发射谱.     

Directionality of generalized acoustic sensors of arbitrary order

For several decades there has been a great deal of interest in acoustic sensors that can make multiple measurements at a single point in the ocean. The order of such sensors has been defined by linking it to the order of the Taylor series approximation of the pressure field at that point. Following this definition, the pressure, vector, and dyadic sensor is of order zero, one, and two, respectively. For this theoretical study, a multichannel three-dimensional spatial filter is derived for a directional acoustic sensor of arbitrary order. Explicit formulas are found for the filter coefficients that maximize the array gain (directivity index) of the filter as well as an explicit expression for the maximum array gain (directivity index). This process is repeated for the case of a first-order null placed in the direction opposite to the look direction of the multichannel filter. Finally, an example is presented which tracks the array gain and beamwidth of a third-order acoustic sensor as the order of the null is assigned values 0, 1, 2, and 3.     

点堆中子动力学方程组曲率权重法的解

针对核反应堆安全工程对某些数值计算结果要求较高的精度和正的误差, 以及舰船核反应堆机动性对计算速度的要求, 需要从数学上寻找一种新的数值计算方法, 以满足实际曲线向上凸或向下凹时计算值总是略高于真实值, 且误差不大于欧拉法和改进的欧拉法所得值. 本文研究曲率权重法求解点堆中子动力学方程组, 该方法是在曲率圆法的基础上引入权重的思想来衡量间隔步长上两个曲率对该步长曲率平均值的贡献. 与欧拉法和改进的欧拉法比较, 曲率权重法的计算结果总是能够高于真实值或有正的误差, 且精度和计算速度得到明显提升. 将该方法用于次临界堆阶跃和线性引入反应性时中子密度的求解, 能够快速得到满足计算要求和高精度的数值结果.     

对称旋转不变相干分布式非圆信号二维波达方向估计

在相干分布式非圆信号二维波达方向估计中,利用信号非圆特性可提升估计精度,但现有的低复杂度算法利用泰勒级数近似建立的旋转不变关系会引入额外误差.针对该问题,考虑中心对称的三维立体线阵,提出了一种基于对称旋转不变关系的二维波达方向估计算法.算法首先利用信号非圆特性建立了扩展阵列模型;然后证明了对于任意的中心对称阵列,相干分布源的确定性角信号分布函数矢量具有对称特性,利用此特性在三维立体线阵的三个子阵中分别建立了扩展广义方向矢量的对称旋转不变关系;基于此,通过无须搜索的多项式求根方式分别得到中心方位角和俯仰角估计;最后利用整个阵列广义方向矢量的对称旋转不变关系构造代价函数实现了参数匹配.理论分析和仿真实验表明,相比于现有的低复杂度算法,所提算法避免了泰勒级数近似引入的额外误差,以较小的复杂度代价获得了性能的较大提升.同时,所提算法能够实现三维空间全方位的角度估计.     

任意阶标度分形格分抗与非正则格型标度方程

Carlson分形格电路是分抗的理想逼近情形,但仅具有负半阶运算性能,逼近效益随着电路节次数的增加逐渐降低.虽然可嵌套得到-1/2~n阶(n为大于或等于2的整数)分抗逼近电路,但结构复杂,无法实现任意分数阶运算.通过类比拓展Carlson分形格电路,获得具有高逼近效益的任意实数阶微积算子的分抗逼近电路——标度分形格分抗,并用非正则格型标度方程进行数学描述.分别探讨非正则格型标度方程的近似求解和真实解.通过调节电阻递进比α与电容递进比β的取值,可构造出具有任意运算阶的标度分形格分抗逼近电路.标度拓展极大地提高了标度分形格分抗电路的逼近效益.随着标度因子的增加,负半阶标度分形格分抗的逼近效益逐渐增大并明显高于Carlson分形格分抗.设计了基于五节Carlson分形格分抗与负半阶标度分形格分抗的半阶微分运算电路,并对周期三角波和周期方波信号进行半阶微分运算,实验测试结果与理论分析一致.     

基于FPGA的浮点指数函数算法研究与实现

基于FPGA的核电站仪控设备中涉及大量浮点指数运算,而常用的CORDIC算法和线性逼近法等存在计算范围小、计算精度不高等问题,对FPGA硬件实现指数函数的方法进行研究,并提出一种改进的级数近似法。该方法对输入进行预处理,将输入分解后采用查找表和泰勒级数展开结合的方法,在展开很少项数的情况下快速收敛,发挥查找表法和级数近似法的优势,提高算法的运算精度和效率。在Matlab环境下对改进算法的有效性进行仿真验证,且采用Verilog语言进行编程实现,在Microsemi公司的IGLOO2系列FPGA上进行具体算法性能验证。Matlab仿真和FPGA验证结果均表明,改进的级数近似法能够大幅增大指数函数的自变量输入范围,并提高计算精度。     

On the Effects of Gaps and Uses of Approximation Functions on the Time-Scale Signal Analysis: A Case Study Based on Space Geophysical Events

The presence of gaps is quite common in signals related to space science phenomena. Usually, this presence prevents the direct use of standard time-scale analysis because this analysis needs equally spaced data; it is affected by the time series borders (boundaries), and gaps can cause an increase of internal borders. Numerical approximations can be used to estimate the records whose entries are gaps. However, their use has limitations. In many practical cases, these approximations cannot faithfully reproduce the original signal behaviour. Alternatively, in this work, we compare an adapted wavelet technique (gaped wavelet transform), based on the continuous wavelet transform with Morlet wavelet analysing function, with two other standard approximation methods, namely, spline and Hermite cubic polynomials. This wavelet method does not require an approximation of the data on the gap positions, but it adapts the analysing wavelet function to deal with the gaps. To perform our comparisons, we use 120 magnetic field time series from a well-known space geophysical phenomena and we select and classify their gaps. Then, we analyse the influence of these methods in two time-scale tools. As conclusions, we observe that when the gaps are small (very few points sequentially missing), all the methods work well. However, with large gaps, the adapted wavelet method presents a better performance in the time-scale representation. Nevertheless, the cubic Hermite polynomial approximation is also an option when a reconstruction of the data is also needed, with the price of having a worse time-scale representation than the adapted wavelet method.     

On the use of the finite difference method in a calculation of vibration-rotation energies

The use of the finite difference method to obtain a Taylor series expansion of a potential energy function for a subsequent calculation of the rovibration energies of molecules is considered. A method is proposed that allows the stability of a finite-difference scheme to be increased against the computational inaccuracy upon numerical expansion of a multidimensional potential energy function into a high-order Taylor series. The method is based on the successive elimination of calculated expansion coefficients of a higher order in calculating the lower-order coefficients by the finite difference method. The approach is illustrated for the example of the CO and H₂S molecules.     

基于SIMPLE的紧致方法

通过泰勒展式系数匹配的方法发展了基于非等距网格的有限容积紧致格式,采用延迟修正的方法建立了基于SIMPLE的紧致方法,,该方法能够得到高精度的数值解,增加迭代求解代数方程组的稳定性。对底部加热的方腔内自然对流换热问题进行数值模拟,结果表明,紧致方法比二阶中心差分方法具有更高的精度。     

Cubic dispersion relation for a relativistic backward-waveoscillator

The cubic approximation to the dispersion relation for a relativistic backward-wave oscillator is obtained, and the utility and limits of the approximation are presented. The approximation is obtained by Taylor series expansion of the wave admittance in the dispersion relation for the transverse-magnetic and free-streaming modes of a relativistic, thin, hollow, cylindrical electron beam moving along the axis of a disc-loaded waveguide in a strong axial magnetic field. The resulting cubic dispersion relation yields instability growth rates and frequencies which fall off beyond their maximum more sharply with increasing wavenumber than for the complete dispersion relation. The approximation is found to be quite good near the operating points of contemporary high-power relativistic backward-wave oscillators, namely, for relatively long wavelength and small ratio of Budker's parameter to the relativistic gamma factor of the beam     

二维可压缩流体力学Lagrange有限点方法

提出二维可压缩流体力学问题的拉格朗日有限点方法,将求解区域离散为适当的点集.在每个时间步,每个离散点与其周围适当的五个邻点组成一个基本计算单元.在每个计算单元上,利用有限点方法中的典型微分算子的五点近似公式直接离散流体力学方程中的微分算子,并在每个方程中加上一个人为拉普拉斯粘性项,达到稳定格式的目的.给出时间步长的自动选取算法.数值算例结果验证了算法的有效性,初步展示了其计算大变形流体问题的良好发展潜力.     

A Borel Transform Method for Locating Singularities of Taylor and Fourier Series

Given a Taylor series with a finite radius of convergence, its Borel transform defines an entire function. A theorem of Pólya relates the large distance behavior of the Borel transform in different directions to singularities of the original function. With the help of the new asymptotic interpolation method of van der Hoeven, we show that from the knowledge of a large number of Taylor coefficients we can identify precisely the location of such singularities, as well as their type when they are isolated. There is no risk of getting artefacts with this method, which also gives us access to some of the singularities beyond the convergence disk. The method can also be applied to Fourier series of analytic periodic functions and is here tested on various instances constructed from solutions to the Burgers equation. Large precision on scaling exponents (up to twenty accurate digits) can be achieved.     

Explicit one-step time discretizations for discontinuous Galerkin and finite volume schemes based on local predictors

We consider a family of explicit one-step time discretizations for finite volume and discontinuous Galerkin schemes, which is based on a predictor-corrector formulation. The predictor remains local taking into account the time evolution of the data only within the grid cell. Based on a space–time Taylor expansion, this idea is already inherent in the MUSCL finite volume scheme to get second order accuracy in time and was generalized in the context of higher order ENO finite volume schemes. We interpret the space–time Taylor expansion used in this approach as a local predictor and conclude that other space–time approximate solutions of the local Cauchy problem in the grid cell may be applied. Three possibilities are considered in this paper: (1) the classical space–time Taylor expansion, in which time derivatives are obtained from known space-derivatives by the Cauchy–Kovalewsky procedure; (2) a local continuous extension Runge–Kutta scheme and (3) a local space–time Galerkin predictor with a version suitable for stiff source terms. The advantage of the predictor–corrector formulation is that the time evolution is done in one step which establishes optimal locality during the whole time step. This time discretization scheme can be used within all schemes which are based on a piecewise continuous approximation as finite volume schemes, discontinuous Galerkin schemes or the recently proposed reconstructed discontinuous Galerkin or P_NP_M schemes. The implementation of these approaches is described, advantages and disadvantages of different predictors are discussed and numerical results are shown.     

On the time evolution of the surface shape of a charged viscous liquid drop deformed at zero time

The problem of capillary oscillations of the equilibrium spherical shape of a charged viscous incompressible liquid drop is solved in an approximation linear in amplitude of the initial deformation that is represented by a finite sum of axisymmetric modes. In this approximation, the shape of the drop as a function of time, as well as the velocity and pressure fields of the liquid in it, may be represented by infinite series in roots of the dispersion relation and by finite sums in numbers of the initially excited modes. In the cases of low, moderate, and high viscosity, the infinite series in roots of the dispersion relation can be asymptotically correctly replaced by a finite number of terms to find compact analytical expressions that are convenient for further analysis. These expressions can be used for finding higher order approximations in amplitude of the initial deformation.     

An Approximate Approach for Systems of Singular Volterra Integral Equations Based on Taylor Expansion

In this article, an extended Taylor expansion method is proposed to estimate the solution of linear singular Volterra integral equations systems. The method is based on combining the m-th order Taylor polynomial of unknown functions at an arbitrary point and integration method, such that the given system of singular integral equations is converted into a system of linear equations with respect to unknown functions and their derivatives. The required solutions are obtained by solving the resulting linear system. The proposed method gives a very satisfactory solution, which can be performed by any symbolic mathematical packages such as Maple, Mathematica, etc. Our proposed approach provides a significant advantage that the m-th order approximate solutions are equal to exact solutions if the exact solutions are polynomial functions of degree less than or equal to m. We present an error analysis for the proposed method to emphasize its reliability. Six numerical examples are provided to show the accuracy and the efficiency of the suggested scheme for which the exact solutions are known in advance.     

一种高精度的修正Hermite-ENO格式

设计一种基于三单元具有六阶精度的修正Hermite-ENO格式（CHENO）,求解一维双曲守恒律问题.CHENO格式利用有限体积法进行空间离散,在空间层上,使用ENO格式中的Newton差商法自适应选择模板.在重构半节点处的函数值及其一阶导数值时,利用Taylor展开给出修正Hermite插值使其提高到六阶精度,并设计了间断识别法与相应的处理方法以抑制间断处的虚假振荡;在时间层上采用三阶TVD Runge-Kutta法进行函数值及一阶导数值的推进.其主要优点是在达到高阶精度的同时具有紧致性.数值实验表明对一维双曲守恒律问题的求解达到了理论分析结果,是有效可行的.