Characteristic-based finite volume scheme for 1D Euler equations

In this paper, a high-order finite-volume scheme is presented for the one- dimensional scalar and inviscid Euler conservation laws. The Simpson＇s quadrature rule is used to achieve high-order accuracy in time. To get the point value of the Simpson＇s quadrature, the characteristic theory is used to obtain the positions of the grid points at each sub-time stage along the characteristic curves, and the third-order and fifth-order central weighted essentially non-oscillatory （CWENO） reconstruction is adopted to estimate the cell point values. Several standard one-dimensional examples are used to verify the high-order accuracy, convergence and capability of capturing shock.     

微分求积单元法在结构工程中的应用

微分求积法（Differential Quadrature Method）是求鳃偏微分方程和积分-微分方程的一种数值方法，该法具有计算简便、精度较高和易于实现等优点。微分求积单元法（Differential Quadrature Element Method）是在微分求积法的基础上结合区域分割和集成规则而形成的一种新的数值计算方法，能通过自适应地选取微分求积网点数目正确模拟构件的刚度和荷载性质，其精度可通过细分单元或增加离散点数目加以提高。微分求积单元法是一种可供选择的、性能优越的数值计算方法。本文将详细论述这一数值方法的基本原理，并通过数值算例说明该方法的应用过程及其优越性，为这一方法在结构工程中的推广应用提供参考。     

轴向均布载荷下压杆稳定问题的DQ解

叙述了微分求积法(differential quadrature method)的一般方法，研究用微分求积法求解在均布轴向载荷下细长杆的稳定问题．通过Newton-Raphson法求解非线性方程组，以及对问题进行线性假设后求解广义特征值方程，得到了精度很高的后屈曲挠度数值和临界载荷数值．与解析解和其他近似解相比，微分求积法具有较高的精度和简便性．     

角度非均匀材料平面V形切口应力奇性分析

提出了一种确定角度非均匀材料平面V形切口尖端应力奇性指数的有效方法。首先,在弹性力学基本方程中引入V形切口尖端位移场的级数渐近展开,建立以位移为特征函数的变系数和非线性微分方程组。然后,采用微分求积法(DQM)求解微分方程组,可得到多阶应力奇性指数及其相对应的特征函数,该法具有公式简单、编程方便、计算量少和精度高等优点,可处理任意开口角度和任意材料组合的V形切口。典型算例验证了微分求积法的有效性和精确性。     

THREE HIGH-ORDER SPLITTING SCHEMES FOR 3D TRANSPORT EQUATION

IntroductionWith the development of modern industry, various pollutants discharge into the air,rivers, lakes and oceans, which makes the environmental qualities worse and has bad effectson the mankind’s health and the sustained development of industry an…     

用DQM求解成层地基一维非线性固结问题

由于多层地基的一维非线性固结问题求解的复杂性，其解析解很难求得。本文基于Davis和Raymond一维非线性固结理论，利用DQM(Differential Quadrature Method)导了初始有效应力沿深度变化、任意边界条件、任意荷载作用下成层地基一维非线性固结的统一表达式，求得了孔压、有效应力和平均固结度的解答。通过解的收敛性分析讨论了DQM解的有效性。由于DQM解对于固结间题各种复杂条件具有统一的矩阵表达式，更便于编程计算和工程应用。最后，用本文解答对三层地基一维非线性固结问题进行了讨论。     

非线性双曲型守恒律的高精度MmB差分格式

构造了一维非线性双曲型守恒律方程的一个高精度、高分辨率的广义G odunov型差分格式。其构造思想是:首先将计算区间划分为若干个互不相交的小区间,再根据精度要求等分小区间,通过各细小区间上的单元平均状态变量,重构各等分小区间交界面上的状态变量,并加以校正;其次,利用近似R iem ann解算子求解细小区间交界面上的数值通量,并结合高阶R unge-K u tta TVD方法进行时间离散,得到了高精度的全离散方法。证明了该格式的Mm B特性。然后,将格式推广到一、二维双曲型守恒方程组情形。最后给出了一、二维Eu ler方程组的几个典型的数值算例,验证了格式的高效性。     

微分求积法的特性及其改进

微分求积法已在科学和工程计算中得到了广泛应用。然而,有关时域微分求积法的数值稳定性、计算精度即阶数等基本特性,仍缺乏系统性的分析结论。依据微分求积法的基本原理,推导证明了微分求积法的权系数矩阵满足V-变换这一重要特性;利用微分求积法和隐式Runge-Kutta法的等值性,证明了时域微分求积法是A-稳定、s级s阶的数值方法。在此基础上,为进一步提高传统微分求积法的计算精度,利用待定系数法和Padé逼近,推导出了一类新的s级2s阶的微分求积法。数值计算对比结果验证了所提出的新微分求积法比传统的微分求积法具有更高的计算精度。     

A CLASS OF COMPACT UPWIND TVD DIFFERENCE SCHEMES

A new method was proposed for constructing total variation diminishing (TVD) upwind schemes in conservation forms. Two limiters were used to prevent non-physical oscillations across discontinuity. Both limiters can ensure the nonlinear compact schemes TVD property. Two compact TVD (CTVD) schemes were tested, one is third-order accuracy, and the other is fifth-order. The performance of the numerical algorithms was assessed by one-dimensional complex waves and Riemann problems, as well as a two-dimensional shock-vortex interaction and a shock-boundary flow interaction. Numerical results show their high-order accuracy and high resolution, and low oscillations across discontinuities.     

Three-dimensional free vibration analysis of functionally graded piezoelectric annular plates on elastic foundations

Three-dimensional free vibration analysis of functionally graded piezoelectric (FGPM) annular plates resting on Pasternak foundations with different boundary conditions is presented. The material properties are assumed to have an exponent-law variation along the thickness. A semi-analytical approach which makes use of state-space method in thickness direction and one-dimensional differential quadrature method in radial direction is utilized to obtain the influences of the Winkler and shearing layer elastic coefficients of the foundations on the non-dimensional natural frequencies of functionally graded piezoelectric annular plates. The analytical solution in the thickness direction can be acquired using the state-space method and approximate solution in the radial direction can be obtained using the one-dimensional differential quadrature method. Numerical results are given to demonstrate the convergency and accuracy of the present method. The influences of the material property graded index, circumferential wave number and thickness of the annular plate on the dynamic behavior are also investigated. Since three-dimensional free vibration analysis of FGPM annular plates on elastic foundations has not been implemented before, the new results can be used as benchmark solutions for future researches.