Construction algorithms for polynomial lattice rules for multivariate integration

We introduce a new construction algorithm for digital nets for integration in certain weighted tensor product Hilbert spaces. The first weighted Hilbert space we consider is based on Walsh functions. Dick and Pillichshammer calculated the worst-case error for integration using digital nets for this space. Here we extend this result to a special construction method for digital nets based on polynomials over finite fields. This result allows us to find polynomials which yield a small worst-case error by computer search. We prove an upper bound on the worst-case error for digital nets obtained by such a search algorithm which shows that the convergence rate is best possible and that strong tractability holds under some condition on the weights.

We extend the results for the weighted Hilbert space based on Walsh functions to weighted Sobolev spaces. In this case we use randomly digitally shifted digital nets. The construction principle is the same as before, only the worst-case error is slightly different. Again digital nets obtained from our search algorithm yield a worst-case error achieving the optimal rate of convergence and as before strong tractability holds under some condition on the weights. These results show that such a construction of digital nets yields the until now best known results of this kind and that our construction methods are comparable to the construction methods known for lattice rules.

We conclude the article with numerical results comparing the expected worst-case error for randomly digitally shifted digital nets with those for randomly shifted lattice rules.

     

MATLAB/SIMULINK在电力系统仿真中的应用

介绍了Matlab／Simulink仿真工具在电力系统仿真中的应用方法．     

Zero Dissipative, Explicit Numerov-Type Methods for Second Order IVPs with Oscillating Solutions

New explicit, zero dissipative, hybrid Numerov type methods are presented in this paper. We derive these methods using an alternative which avoids the use of costly high accuracy interpolatory nodes. We only need the Taylor expansion at some internal points then. The method is of sixth algebraic order at a cost of seven stages per step while their phase lag order is fourteen. The zero dissipation condition is satisfied, so the methods possess an non empty interval of periodicity. Numerical results over some well known problems in physics and mechanics indicate the superiority of the new method.     

Atomic-scale simulation of nano-grains:structure and diffusion properties

Nanograins are characterized by a typical grain size from lto 100 nm.Mclecular dynamics aimubations have been carried out for the nanograin sphere with the diameters from 1.45 to 10.12 nm.We study the influence of grain size on structure and diffusion properties of the nanograins.The results reveal that as the grain size is reduced,the fraction of grain surface increases significantly,and the surface width is approximately constant；the diffusicn coefficlent is increased sharply,and the relation of the diffusion coefficient and the grain size is close to exponential relation below 10nm.     

Structure of solid monolayers and multilayers of n -hexane on graphite

We present all-atom molecular dynamics simulations ofn-hexane on the basal plane of graphite at monolayer and multilayer coverages. In keeping with experimental data, we find the presence of ordered adsorbed layers both at single monolayer coverage and when the adsorbed layer coexists with excess liquid adsorbate. Using a simulation method that does not impose any particular periodicity on the adsorbed layer, we quantitatively compare our results to the results of neutron diffraction experiments and find a structural transition from a uniaxially incommensurate lattice to a fully commensurate structure on increasing the coverage from a monolayer to a multilayer. The zig-zag backbone planes of all the alkane molecules lie parallel to the graphite surface at the multilayer coverage, while a few molecules are observed to attain the perpendicular orientation at monolayer coverage. Dedicated to Professor C N R Rao on his 70th birthday     

Numerical simulation of thermally and chemically nonequilibrium flows and heat transfer in underexpanded induction plasmatron jets

The results of numerical simulation are presented for thermally and chemically nonequilibrium air plasma flows in a plasmatron discharge channel and underexpanded dissociated and partially ionized air jets flowing past a cylindrical model with a blunt leading edge and cooled copper surface under the experimental conditions realized in a VGU-4 100 kW induction plasmatron (Institute for Problems in Mechanics of the Russian Academy of Sciences) (see, for example, [1, 2]). The nonequilibrium excitation of the vibrational degrees of freedom of the molecules in the modal approximation and the difference between the electron and translational heavy-particle temperatures are taken into account in the calculations. The calculated data on the heat transfer and pressure at the stagnation point are compared with the results obtained within the framework of the thermally equilibrium model. Comparison with the experimental data obtained in the Institute for Problem in Mechanics of the Russian Academy of Sciences (Laboratory for interaction between plasma and radiation and materials) and kindly provided for comparison purposes gives satisfactory agreement.     

可压稠密气固两相流动过程的模拟计算

基于稠密气体分子运动论和颗粒动理学,建立可压稠密气固两相流动模型。采用梯度模拟来考虑气相可压缩性对气相湍流的影响。模拟计算表明气固两相射流速度沿轴向和径向减小,颗粒浓度下降。气固两相射流具有高的颗粒温度,呈现强烈的气固两相湍流流动特性。     

一类复杂循环排队系统的计算机模拟及其在进程调度中的应用

讨论了顾客到达时间和服务时间与等待队列队长有关的排队系统模拟,顾客服务次数不止1次即出现循环的复杂排队系统,引入到达因子、服务台因子和循环因子,在串并联多服务台情形下建立了4类计算机模拟模型,给出该复杂循环排队系统在进程调度中的应用并进行了计算机模拟.     

基于MATLAB的异步电动机线性化控制系统的仿真

从异步电动机在同步旋转坐标系下的状态方程出发，在一定的条件下，对其系数矩阵简化．推出异步电动机的线性化控制模型．并基于该模型得出异步电动机的线性化控制系统．以给定电机为例，对该模型的有效性、响应用MATLAB进行仿真分析．验证的结果说明上述模型具有实用价值．     

张量对称类上诱导线性算子的数值半径

设V是一个n维线性空间,V_x~m(G)为V上的张量对称类.A为V的线性算子T的矩阵,K(A)为V_x~m(G)上的诱导线性算子K(T)的矩阵.本文从K(A)的数值半径Υ(K(A))和可分数值半径Υ_x(K(A))定义出发,研究了Υ(K(A))、Υ_x(K(A))与范数||A||_p(1≤p≤2)、广义矩阵函数d_x~G(A)的关系,得到了它们之间的两个不等式.