Navier-Stokes方程的线性化微分求积法

提出一种新的线性化微分求积法(LDQM),将这种目的应用到流函数和涡量形式的Navier-Stokes方程.通过LDQM,非线性方程很容易被解出来,并且容易处理压力的边界条件.为检验本目的,计算了两个数值算例.     

采用微分求积法分析薄板的气动弹性响应

为减轻结构重量，新一代飞行器的结构设计中更多地采用了薄壁结构，在超音速或超高音速飞行条件下，薄板的气动弹性响应相当剧烈。分析薄板的气动弹性响应方法主要有两类：第一类为经典的伽辽金方法；第二类为有限元方法。受薄板形状和边界条件的限制，伽辽金方法能够研究的问题是非常有限的，有限元方法虽然具有普适性，但本质上属于数值方法，其计算精度和收敛性必将受所选单元类型以及数值计算误差的影响。在薄板的气动弹性分析领域，文中采用一种全新的方法，即微分求积方法。     

抛物型方程一类自由边界问题的微分求积区域分裂法

在微分求积法的基础上,结合区域分裂法的优点提出了一种新的数值计算方法——微分求积区域分裂法.数值试验表明,该方法在求解抛物型方程一类初值带有弱奇性的自由边界问题时十分灵活有效.     

基于微分单元的圆柱内辐射导热耦合传热分析

本文充分利用微分单元法不涉及积分、较为稳定、易于实施的特点,进行圆柱坐标系下辐射与导热耦合传热分析.在求解过程中,将热辐射的空间和角度耦合成三维计算域求解,实现空间和角度的高阶离散.针对具有强对流特性的辐射传递方程,采用一种迎风格式可有效抑制数值振荡,实现稳定计算.由于辐射边界条件与导热和对流边界条件不同,它是具有方向性的单向边界条件,且在边界处会存在强烈的间断性和奇异点.为了有效抑制间断并处理奇异点,本文提出双层间断节点方案.通过与文献中圆柱坐标系下辐射与导热耦合传热算例进行对比分析,发现微分单元模型在空间和角度方向上均能实现高分辨率刻画,获得较好的计算精度.     

Galerkin时空耦合谱元法求解声波动方程

提出了时空耦合谱元方法,并将其用于带第一类边界条件的非齐次一维、二维、三维波动方程的求解。分别采用四边形、六面体和超六面体作为计算单元,在每个单元内采用Chebyshev多项式的极值点作为Lagrange插值节点,并且探讨了区域剖分方式对计算精度的影响。时空耦合谱元法能够得到精度很高的数值结果,并且其色散随时间推移是稳定的;当总网格节点数相同时,不同的网格剖分方式所得数值误差不同,当空间方向Chebyshev多项式的阶数较高和时间方向Chebyshev多项式的阶数较低时,得到的数值精度较高;在总节点数相同的情况下,与时间全域方式相比,逐时间子区域方式计算所需要的时间更经济,两种方式可以得到相同的精度。结果表明:时空耦合谱元方法使时空方向精度相匹配,可以提高整体精度;空间方向的Chebyshev多项式对数值精度起主要影响作用;时间子区域方式的采用可以扩大问题的计算区域。      

综合核方法求解辐射输运问题的求积组选取

简要介绍求解辐射输运方程的综合核方法,分析计算误差和收敛性,提出新的求积组和误差修正方法,提高综合核方法的计算精度.通过对基准问题的计算比对表明,采用提出的求积组并通过误差修正,综合核方法在低阶时的结果具有较高的计算精度.     

二阶Y环频率选择表面的设计研究

利用分形结构的自相似性将分形理论应用于频率选择表面(FSS)领域即可使在单屏FSS上具有多频段带通滤波的特性,结合Floquet周期边界条件,采用矩量法研究了二阶Y环分形FSS在不同入射角度下迭代比例因子F及单元排布方式对频率响应特性的影响规律,给出谐振频率的经验估算值.计算及实验结果表明,FSS的谐振频率主要由迭代比例因子及起始单元尺寸决定,而透过率及-3 dB带宽则对排布方式的改变较敏感.实验结果与理论分析一致.     

结合小波变换与微分法改善近红外光谱分析精度

微分法可以有效消除光谱背景和基线漂移,同时会增加光谱噪音;小波变换具有很好的去噪功能,章将微分法和小波变换结合用于重整汽油辛烷值近红外光谱分析。考察了微分噪音对辛烷值分析精度的影响以及小波去噪对微分光谱的噪音扣除以及对辛烷值分析精度改善情况。结果表明,微分光谱可以扣除原始光谱的基线漂移,提高分析精度,同时增加光谱的噪音;噪音对分析精度影响很大。微分光谱经过小波去噪处理后信噪比增加,辛烷值分析精度得到改善。     

微分吸收法测量二极管电压

给出了微分吸收法测量二极管电压的基本原理和实验结果。利用MCNP程序对轫致辐射-衰减-探测器系统建模,模拟得到了输出剂量与二极管工作电压关系拟合曲线。建立了微分吸收法测量二极管电压测量系统,通过在探测器前端放置不同厚度的吸收片,得到了衰减程度不同的波形。结合理论计算的拟合曲线和实验波形,利用迭代法计算得到了晨光号加速器二极管电压,电压峰值为0.58 MV。和传统方法所测得二极管工作电压进行了比较,结果较为一致。     

电磁场中光子-轴子的转化微分截面

运用费曼微扰方法分别计算了在磁偶极场、电偶极场和均匀静电场及静磁场中光子转化成轴子的非极化微分截面.在电偶极场中,沿光子传播方向及其反方向上的非极化微分截面为零;而在磁偶极场中,在上述方向上通常则具有非零的微分截面,但当光子传播方向平行于磁场偶极距矢量时,该微分截面为零.在均匀的静磁场和均匀静电场中,只有在光子传播方向及其反方向上具有非零的微分截面,但后者小于前者.在轴子质量趋于零的极限条件下,上述过程和光子转化为引力子的过程表现出某些非常类似的性质. 关键词： 轴子 光子 微分截面     

Buckling analysis of nanobeams with exponentially varying stiffness by differential quadrature method

We present the application of differential quadrature(DQ) method for the buckling analysis of nanobeams with exponentially varying stiffness based on four different beam theories of Euler-Bernoulli, Timoshenko, Reddy, and Levison.The formulation is based on the nonlocal elasticity theory of Eringen. New results are presented for the guided and simply supported guided boundary conditions. Numerical results are obtained to investigate the effects of the nonlocal parameter,length-to-height ratio, boundary condition, and nonuniform parameter on the critical buckling load parameter. It is observed that the critical buckling load decreases with increase in the nonlocal parameter while the critical buckling load parameter increases with increase in the length-to-height ratio.     

An efficient differential quadrature methodology for free vibration analysis of arbitrary straight-sided quadrilateral thin plates

In this paper, a new differential quadrature (DQ) methodology is employed to study free vibration of irregular quadrilateral straight-sided thin plates. A four-nodded super element is used to map the irregular physical domain into a square domain in the computational domain. Second order transformation schemes with relative ease and less computation are employed to transform the fourth order governing equation of thin plates between the two domains. The only degree of freedom within the domain is the displacement, whereas along the boundaries, the displacement as well as the second order derivative of the displacement with respect to associated normal co-ordinate variable in computational domain are the two degrees of freedom. Implementing the method, the formulation for the DQ method for the free vibration analysis of plates of straight-sided shapes was presented together with the implementation procedure for the different boundary conditions. To demonstrate the accuracy, convergency and stability of the new methodology, detail studies are made on isotropic plates at acute angles with different geometries, boundary and loading conditions including DQ free-edge boundary condition implementations. Accurate results even with fewer degrees of freedom than for those of comparable numerical algorithms were achieved.     

Simulating non-linear coupled oscillators by an iterative differential quadrature method

An iterative method based on differential quadrature rules is proposed as a new unified frame of resolution for non-linear two-degree-of-freedom systems. Dynamical systems with Duffing-type non-linearity have been considered. Differential quadrature rules have been applied with a careful distribution of sampling points to reduce the governing equation of motion to two second-order non-linear, non-autonomous ordinary differential equations and to solve the time-domain problem. The time domain of the problem is discretized by means of time intervals, with the same distribution of sampling points used to discretize the space domain (which can be seen as a single interval). It will be shown that accurate solutions depend not only on the choice of the distribution of sampling points, but also on the length of the time interval one refers to in the computations. The numerical results, utilized to draw Poincaré maps, are successfully compared with those obtained using the Runge-Kutta method.     

Nonlinear free vibrations of curved double walled carbon nanotubes using differential quadrature method

Nonlinear free vibration analysis of curved double-walled carbon nanotubes (DWNTs) embedded in an elastic medium is studied in this study. Nonlinearities considered are due to large deflection of carbon nanotubes (geometric nonlinearity) and nonlinear interlayer van der Waals forces between inner and outer tubes. The differential quadrature method (DQM) is utilized to discretize the partial differential equations of motion in spatial domain, which resulted in a nonlinear set of algebraic equations of motion. The effect of nonlinearities, different end conditions, initial curvature, and stiffness of the surrounding elastic medium, and vibrational modes on the nonlinear free vibration of DWCNTs is studied. Results show that it is possible to detect different vibration modes occurring at a single vibration frequency when CNTs vibrate in the out-of-phase vibration mode. Moreover, it is observed that boundary conditions have significant effect on the nonlinear natural frequencies of the DWCNT including multiple solutions.     

Element-free Galerkin (EFG) method for analysis of the time-fractional partial differential equations

The present paper deals with the numerical solution of time-fractional partial differential equations using the element-free Galerkin (EFG) method, which is based on the moving least-square approximation. Compared with numerical methods based on meshes, the EFG method for time-fractional partial differential equations needs only scattered nodes instead of meshing the domain of the problem. It neither requires element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. In this method, the first-order time derivative is replaced by the Caputo fractional derivative of order α (0<α ≤1). The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Several numerical examples are presented and the results we obtained are in good agreement with the exact solutions.     

Application of the generalized differential quadrature rule to initial-boundary-value problems

Partial differential equations (PDEs) for the forced vibration of structural beams are solved in this paper using the recently proposed generalized differential quadrature rule (GDQR). The GDQR techniques are first applied to both spatial and time dimensions simultaneously as a whole. No other classical methods are needed in the time dimension. The objective of this paper is to formularize the GDQR expressions and corresponding explicit weighting coefficients, while the derivation of explicit weighting coefficients is one of the most important aspects in the differential quadrature methods. It should be emphasized that the GDQR expressions and weighting coefficients for two-dimensional problems are not a direct application of those for one-dimensional problems, and they are distinctly different for PDEs of different orders. An Euler beam and a Timoshenko beam are employed as examples. Accurate results are obtained. The proposed procedures can be applied to problems in other disciplines of sciences and technology, where the problems may be governed by other PDEs with different orders in the time or spatial dimension.     

Stability Analysis and Order Improvement for Time Domain Differential Quadrature Method

The differential quadrature method has been widely used in scientific and engineering computation. However, for the basic characteristics of time domain differential quadrature method, such as numerical stability and calculation accuracy or order, it is still lack of systematic analysis conclusions. In this paper, according to the principle of differential quadrature method, it has been derived and proved that the weighting coefficients matrix of differential quadrature method meets the importantV-transformation feature. Through the equivalence of the differential quadrature method and the implicit Runge-Kutta method, it has been proved that the differential quadrature method is A-stable and $s$-stage $s$-order method. On this basis, in order to further improve the accuracy of the time domain differential quadrature method, a class of improved differential quadrature method of $s$-stage $2s$-order has been proposed by using undetermined coefficients method and Padéapproximations. The numerical results show that the proposed differential quadrature method is more precise than the traditional differential quadrature method.     

Gauss quadrature based finite temperature Lanczos method

The finite temperature Lanczos method (FTLM), which is an exact diagonalization method intensively used in quantum many-body calculations, is formulated in the framework of orthogonal polynomials and Gauss quadrature. The main idea is to reduce finite temperature static and dynamic quantities into weighted summations related to one- and two-dimensional Gauss quadratures. Then lower order Gauss quadrature, which is generated from Lanczos iteration, can be applied to approximate the initial weighted summation. This framework fills the conceptual gap between FTLM and kernel polynomial method, and makes it easy to apply orthogonal polynomial techniques in the FTLM calculation.     

Conditional quadrature method of moments for kinetic equations

Kinetic equations arise in a wide variety of physical systems and efficient numerical methods are needed for their solution. Moment methods are an important class of approximate models derived from kinetic equations, but require closure to truncate the moment set. In quadrature-based moment methods (QBMM), closure is achieved by inverting a finite set of moments to reconstruct a point distribution from which all unclosed moments (e.g. spatial fluxes) can be related to the finite moment set. In this work, a novel moment-inversion algorithm, based on 1-D adaptive quadrature of conditional velocity moments, is introduced and shown to always yield realizable distribution functions (i.e. non-negative quadrature weights). This conditional quadrature method of moments (CQMOM) can be used to compute exact N-point quadratures for multi-valued solutions (also known as the multi-variate truncated moment problem), and provides optimal approximations of continuous distributions. In order to control numerical errors arising in volume averaging and spatial transport, an adaptive 1-D quadrature algorithm is formulated for use with CQMOM. The use of adaptive CQMOM in the context of QBMM for the solution of kinetic equations is illustrated by applying it to problems involving particle trajectory crossing (i.e. collision-less systems), elastic and inelastic particle–particle collisions, and external forces (i.e. fluid drag).     

Vibration analysis of a circular arch with variable cross-section using differential transformation and generalized differential quadrature

Vibration analysis of circular arches is an important subject in mechanics due to its various applications. In particular, circular arches with variable cross-section have been widely used to satisfy modern architectural and structural requirements. Recently, the generalized differential quadrature method (GDQM) and differential transformation method (DTM) were proposed by Shu and Zhou, respectively. In this study, GDQM and DTM are applied to vibration analysis of circular arches with variable cross-section. The governing equation of motion is derived and the non-dimensional natural frequencies are obtained for various boundary conditions. The concepts of differential transformation and generalized differential quadrature are briefly introduced. The results obtained by these methods are compared with previously published works. GDQM and DTM showed fast convergence, accuracy and validity in solving the vibration problem for circular arches with variable cross-sections.