New Sixth-Order Compact Schemes for Poisson/Helmholtz Equations

Some new sixth-order compact finite difference schemes for Poisson/Helmholtz equations on rectangular domains in both two- and three-dimensions are developed and analyzed. Different from a few sixth-order compact finite difference schemes in the literature, the finite difference and weight coefficients of the new methods have analytic simple expressions. One of the new ideas is to use a weighted combination of the source term at staggered grid points which is important for grid points near the boundary and avoids partial derivatives of the source term. Furthermore, the new compact schemes are exact for 2D and 3D Poisson equations if the solution is a polynomial less than or equal to 6. The coefficient matrices of the new schemes are $M$-matrices for Helmholtz equations with wave number $K≤0,$ which guarantee the discrete maximum principle and lead to the convergence of the new sixth-order compact schemes. Numerical examples in both 2D and 3D are presented to verify the effectiveness of the proposed schemes.     

二维振荡叶栅非定常粘性流动数值模拟

采用显式四步Runge-Kutta格式，结合Baldwon-Lomax紊流模型求解Navier-Stokes方程，借助运动网格技术，完成了对二维振荡叶栅非定常粘性流动的数值模拟。为了加速求解过程，引入了变系数隐式残差光顺方法，取得了较好效果。数值结果与已公布的数据有很好的一致性。     

面向星敏感器的星模式识别算法

介绍了到目前为止出现的所有面向星敏感器的星模式识别算法,它们是概率统计算法、三角形算法、匹配组算法、网格算法、奇异值分解算法、神经网络算法和遗传算法,并分两组比较了它们的性能,给出了具有指导意义的结论。     

The normal contact of rough surfaces with fine discretization

Summary The numerical treatment of a 3-d frictionless normal contact of rough bodies made of a linear–elastic, isotropic material is considered. First, the governing equations and their discretized forms are presented. Due to the order of the problem, the computation time rapidly increases with the grid refinement. Therefore, common computer programs have to be limited to 50×50 contact points. Several numerical techniques (multigrid, iterative solver, distributive relaxation) are used. Cell-centered grid transformation is applied and grid operators have been established for this purpose. In order to handle upto 1000 × 1000 cells the program CON3D has been developed.     

Natural Vibrations of Ribbed Cylindrical Shells with Low Shear Stiffness

A technique is proposed to determine the natural frequencies of rib-reinforced cylindrical shells with low shear stiffness. The equations of motion based on the Timoshenko model are used. The influence of the transverse shear moduli of the shell and ribs on the minimum natural frequencies and corresponding modes is studied by way of numerical examples. It is shown that the effect of the discrete ribs is different from that predicted by the classical theory of ribbed shells__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 4, pp. 66–74, April 2005.     

一种四面体网格局部自适应方法在有限元中的应用

通过提出一种针对三维四面体网格的区域整形技术，给出了一种简便实用的四面体网格局部自适应加密方法，并将其与流量修正有限元法结合，应用于三维高速流的计算     

建筑风压场数值模拟中压力波动的平抑

在以同位网格为基础的简单流场压力计算中,通常采用动量插值方法来平抑流场中的压力波动现象;但是对于建筑风场等复杂的钝体绕流问题,由该平抑方法得到的收敛风压场仍可能存在小幅波动。为彻底解决同位网格格式下的压力波动,除采用动量插值方法外,本文提出了在压力校正方程的界面流速中添加压力梯度差值项的方法。算例分析表明,该方法计算得到的建筑风压场完全避免了压力波动现象,风压解与试验结果吻合良好。     

平板网架结构分析的超级条元法

本文针对复杂杆系结构提出一种超级有限条的分析方法。这种方法基于半连续半离散思想将各杆件结点自由度归结为整体结构自由度；进一步又将整体自由度简化为少数结线自由度上机计算。这将为复杂结构体系开拓一条简便、实用的工程分析途径。     

双重点最小二乘配点无网格法

配点类无网格法需要计算近似函数的二阶导数,因而在移动最小二乘(MLS)近似中至少要采用二次基函数。本文利用Voronoi图对双重点移动最小二乘近似法进行了改进,建立了基于Voronoi图的双重点移动最小二乘近似(VDG),并利用加权最小二乘法离散微分方程,导出了双重点最小二乘配点无网格法(MD GLS)。该方法将求解域用节点离散,并以节点为生成点建立Voronoi图,取Voronoi多边形的顶点为辅助点。近似函数及其二阶导数的计算过程可分解为两个步骤:首先用场函数节点值拟合辅助点处近似函数的一阶导数,再以辅助点处近似函数的一阶导数值拟合节点处近似函数的二阶导数。由于在每一步中只需计算MLS形函数及其一阶导数,这种近似方法需要较少的影响点和较小的影响域。同时借助于Voronoi结构的优良几何性质,可以快速地搜索影响点。研究表明,与基于MLS的加权最小二乘无网格法(MWLS)相比,这种方法可以显著提高计算效率,并且在精度和收敛性方面也有所改善。     

Transient gas flow simulation using an Adaptive Method of LinesSimulation d'écoulements dans un gazoduc par la méthode adaptative de lignes

Designing natural gas pipelines to safely and efficiently handle unsteady flows, requires knowledge of pressure drop, flowrate and temperature distribution throughout the system. The accurate prediction of these parameters is essential in order to achieve optimum cumulative deliverability, and safe and reliable operation. An Adaptive Method of Lines algorithm is formulated for the solution of Euler system of equations, which fully simulates slow and fast transients. Two test cases present the improvement of the numerical solution from grid adaptation. Good results are obtained both for slow and fast transients simulations proving that the suggested numerical procedure is appropriate for such predictions. To cite this article: E. Tentis et al., C. R. Mecanique 331 (2003).