叶型优化在非定常流动条件下效果的数值分析

本文对叶型优化在非定常流动条件下的效果进行数值研究。叶型优化基于比较成熟的在定常流动条件下的反问题解法，非定常流动则考虑到动静叶栅之间的相互作用，非定常流动的求解在Ｎ－Ｓ方程的基础上采用分区计算的方法来完成。数值计算的结果表明，定常流动条件下优化得到的叶型在非定常流动条件下同样具有较好的气动性能。     

前缘弯掠斜流转子叶顶间隙流动特性数值分析

本文对前缘弯掠斜流转子叶顶间隙内的流动特性进行了数值分析。结果表明:叶顶间隙气流与主流发生卷吸而生成泄漏涡。泄漏涡作用的区域具有较低的压力分布。在叶片通道内,泄漏涡沿着与转子旋向相反的方向朝相邻叶片的压力面移动。大间隙时的泄漏涡比小间隙时强烈。低流量时泄漏涡的作用区域比高流量时大。在各种流量特性下,叶顶尾缘近吸力面区域都存在着二次间隙流。     

Numerical analysis of turbulent structure in compound meandering open channel by algebraic Reynolds stress model

The turbulent flow in a compound meandering channel with a rectangular cross section is one of the most complicated turbulent flows, because the flow behaviour is influenced by several kinds of forces, including centrifugal forces, pressure‐driven forces and shear stresses generated by momentum transfer between the main channel and the flood plain. Numerical analysis has been performed for the fully developed turbulent flow in a compound meandering open‐channel flow using an algebraic Reynolds stress model. The boundary‐fitted coordinate system is introduced as a method for coordinate transformation in order to set the boundary conditions along the complicated shape of the meandering open channel. The turbulence model consists of transport equations for turbulent energy and dissipation, in conjunction with an algebraic stress model based on the Reynolds stress transport equations. With reference to the pressure–strain term, we have made use of a modified pressure–strain term. The boundary condition of the fluctuating vertical velocity is set to zero not only for the free surface, but also for computational grid points next to the free surface, because experimental results have shown that the fluctuating vertical velocity approaches zero near the free surface. In order to examine the validity of the present numerical method and the turbulent model, the calculated results are compared with experimental data measured by laser Doppler anemometer. In addition, the compound meandering open channel is clarified somewhat based on the calculated results. As a result of the analysis, the present algebraic Reynolds stress model is shown to be able to reasonably predict the turbulent flow in a compound meandering open channel. Copyright © 2005 John Wiley & Sons, Ltd.     

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.

     

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.     

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.     

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

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

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

设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)的关系,得到了它们之间的两个不等式.     

General expressions of peaked traveling wave solutions of CH-γ and CH equations

We use qualitative analysis and numerical simulation to study peaked traveling wave solutions of CH-γ and CH equations. General expressions of peakon and periodic cusp wave solutions are obtained. Some previous results become our special cases.     

A test problem for molecular dynamics integrators

We derive a test problem for evaluating the ability of time-steppingmethods to preserve the statistical properties of systems inmolecular dynamics. We consider a family of deterministic systemsconsisting of a finite number of particles interacting on acompact interval. The particles are given random initial conditionsand interact through instantaneous energy- and momentum-conservingcollisions. As the number of particles, the particle density,and the mean particle speed go to infinity, the trajectory ofa tracer particle is shown to converge to a stationary Gaussianstochastic process. We approximate this system by one describedby a system of ordinary differential equations and provide numericalevidence that it converges to the same stochastic process. Wesimulate the latter system with a variety of numerical integrators,including the symplectic Euler method, a fourth-order Runge-Kuttamethod, and an energyconserving step-and-project method. Weassess the methods' ability to recapture the system's limitingstatistics and observe that symplectic Euler performs significantlybetter than the others for comparable computational expense.