求解声辐射问题的微分求积法与微分求积单元法

基于无限域中的Helmholtz波动方程,将微分求积法与微分求积单元法应用于二维及三维声辐射问题的求解,对最外层节点施加不同阶数的人工边界条件,区域内使用均匀及非均匀的节点分布方式,分析了节点分布方式及人工边界条件对计算结果的影响,比较了两种数值方法的计算精度。研究结果表明:微分求积法与微分求积单元法,前者精度更高,而后者耗时更少,在频率较低时,具备较高的效费比。人工边界条件对计算结果的影响主要体现在低频段,而节点分布方式的影响主要体现在高频段。非均匀的节点分布方式在不同频段都具备更好的计算精度。      

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.     

Semi-analytical solution for three-dimensional transient response of functionally graded annular plate on a two parameter viscoelastic foundation

The three-dimensional transient analysis of functionally graded annular plates with arbitrary boundary conditions is carried out in this paper. The material properties of the FGM plate are assumed to vary smoothly in an exponential law along the thickness direction. The plate is assumed to rest on a two parameter viscoelastic foundation. A semi-analytical method, which integrates the state space method (SSM), Laplace transform and its inversion, as well as the one-dimensional differential quadrature method (DQM), is proposed to obtain the transient response of the plate. The state space method is used to obtain the analytical solution in the thickness direction. The differential quadrature method is employed to approximate the solution in the radial direction. The Laplace transform and the numerical inversion are used to obtain the solution in time domain. Numerical results show a good agreement between the response histories obtained by the present method and finite element method. The effects of the boundary conditions at the edges, the material graded index, the Winkler and shearing layer elastic coefficients, and the damping coefficient are studied. Numerical examples show that the peak response decreases as the material graded index, the Winkler and shearing layer elastic coefficients, and the damping coefficient increase. The results obtained in this paper can serve as benchmark data in further research.     

二维各向异性磁等离子体的无条件稳定ADE-CNAD-FDTD算法

针对二维各向异性磁等离子体提出一种有效的无条件稳定算法,新算法结合了辅助微分方程(ADE)方法与Crank-Nicolson approximate-decoupling(CNAD)时域有限差分算法仿真各向异性磁等离子体介质。传统的ADE-FDTD方法应用在一维各向异性色散介质具有较高的精度和效率,将提出的新算法ADE-CNAD-FDTD应用到二维各向异性磁等离子体介质中不仅解决了电磁波在具有各向异性和频率色散特性介质中传播的仿真难题,而且去除了CFL稳定性条件。该算法在保留了原有的精度情况下大幅度地提高了计算效率并成为无条件稳定的形式。给出一个算例证明该算法的有效性,通过模拟电磁波在磁等离子体中的传播,仿真结果与传统的ADE-FDTD算法对比,证实了该算法的高效率、无条件稳定性和高精度。     

光子晶体传输特性的时域精细积分法分析

将精细积分法应用于时域有限差分法中,提出了一种求解光子晶体传输特性的时域精细积分法,并对其计算精度及稳定性进行了分析.从一阶麦克斯韦方程出发,在空间上采用Yee元胞进行差分离散,结合吸收边界条件及激励源表达式将方程整理为标准的一阶常微分方程组形式.通过时间步长的精细划分和指数矩阵的加法定理,在时间上利用精细积分法对齐次微分方程进行积分求解,并结合激励向量的特解得到空间离散的场分量,最终通过傅里叶变换求得方程的解.利用时域精细积分法对光子晶体进行了实例计算,并将其结果分别与时域有限差分法和四阶龙格库塔法在精度、稳定性等方面进行了比较,结果表明时域精细积分法具有更高的计算精度,并且克服了时域有限差分法以及四阶龙格库塔法在计算稳定性上对时间步长的限制.提出的方法具有精确、稳定的特点,为光子晶体传输特性的研究提供了一种新的有效的分析方法.     

Transient heat conduction analysis using the MLPG method and modified precise time step integration method

The meshless local Petrov–Galerkin (MLPG) method in conjunction with the modified precise time step integration method in the time domain is proposed for transient heat conduction analysis in this paper. The MLPG method is often referred to as a truly meshless method because it requires no elements or background cells for either field interpolation or background integration. Local weak forms are developed using weighted residual method locally from the partial differential equation of transient heat conduction. In order to simplify the treatment of essential boundary conditions, the natural neighbour interpolation (NNI) is employed for the construction of trial functions. Moreover, the three-node triangular FEM shape functions are taken as test functions to reduce the order of integrands involved in domain integrals. The semi-discrete heat conduction equation is solved numerically with modified precise time step integration method in the time domain. The availability and accuracy of the present method for transient heat conduction analysis are tested through numerical examples.     

A non-uniform grid for triangular differential quadrature

The triangular differential quadrature method based on a non-uniform grid is proposed in the paper. Explicit expressions of the non-uniform grid point coordinates are given and the weighting coefficients of the triangular differential quadrature method are determined with the aid of area coordinates. Two typical examples are presented to testify the effectiveness of the non-uniform grid. It is shown that rapid convergence is achieved under the non-uniform grid in comparison with those from the uniform grid with the same order of approximation.     

结构动力分析的时域配点DQ半解析法

直接从控制微分方程出发,以梁为对象提出了本计算方法。该方法在空间域采用DQ法(DiferentialQuadratureMethod),在时间域取级数,采用时域配点的方式得到求位移场全部待定参数的可解方程组,求解一次线性方程组即可求得全域的响应位移场。     

A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid

A novel accurate numerical model for shallow water equations on sphere have been developed by implementing the high order multi-moment constrained finite volume (MCV) method on the icosahedral geodesic grid. High order reconstructions are conducted cell-wisely by making use of the point values as the unknowns distributed within each triangular cell element. The time evolution equations to update the unknowns are derived from a set of constrained conditions for two types of moments, i.e. the point values on the cell boundary edges and the cell-integrated average. The numerical conservation is rigorously guaranteed. In the present model, all unknowns or computational variables are point values and no numerical quadrature is involved, which particularly benefits the computational accuracy and efficiency in handling the spherical geometry, such as coordinate transformation and curved surface.Numerical formulations of third and fourth order accuracy are presented in detail. The proposed numerical model has been validated by widely used benchmark tests and competitive results are obtained. The present numerical framework provides a promising and practical base for further development of atmospheric and oceanic general circulation models.     

Adaptively deformed mesh based interface method for elliptic equations with discontinuous coefficients

Mesh deformation methods are a versatile strategy for solving partial differential equations (PDEs) with a vast variety of practical applications. However, these methods break down for elliptic PDEs with discontinuous coefficients, namely, elliptic interface problems. For this class of problems, the additional interface jump conditions are required to maintain the well-posedness of the governing equation. Consequently, in order to achieve high accuracy and high order convergence, additional numerical algorithms are required to enforce the interface jump conditions in solving elliptic interface problems. The present work introduces an interface technique based adaptively deformed mesh strategy for resolving elliptic interface problems. We take the advantages of the high accuracy, flexibility and robustness of the matched interface and boundary (MIB) method to construct an adaptively deformed mesh based interface method for elliptic equations with discontinuous coefficients. The proposed method generates deformed meshes in the physical domain and solves the transformed governed equations in the computational domain, which maintains regular Cartesian meshes. The mesh deformation is realized by a mesh transformation PDE, which controls the mesh redistribution by a source term. The source term consists of a monitor function, which builds in mesh contraction rules. Both interface geometry based deformed meshes and solution gradient based deformed meshes are constructed to reduce the L(∞) and L(2) errors in solving elliptic interface problems. The proposed adaptively deformed mesh based interface method is extensively validated by many numerical experiments. Numerical results indicate that the adaptively deformed mesh based interface method outperforms the original MIB method for dealing with elliptic interface problems.     

Some theoretical aspects of elastic wave modeling with a recently developed spectral element method

A spectral element method has been recently developed for solving elastodynamic problems. The numerical solutions are obtained by using the weak formulation of the elastodynamic equation for heterogeneous media, based on the Galerkin approach applied to a partition, in small subdomains, of the original physical domain. In this work, some mathematical aspects of the method and the associated algorithm implementation are systematically investigated. Two kinds of orthogonal basis functions, constructed with Legendre and Chebyshev polynomials, and their related Gauss-Lobatto collocation points are introduced. The related integration formulas are obtained. The standard error estimations and expansion convergence are discussed. An element-by-element pre-conditioned conjugate gradient linear solver in the space domain and a staggered predictor/multi-corrector algorithm in the time integration are used for strong heterogeneous elastic media. As a consequence, neither the global matrices nor the effective force vector is assembled. When analytical formulas are used for the element quadrature, there is even no need for forming element matrix in order to further save memory without losing much in computational efficiency. The element-by-element algorithm uses an optimal tensor product scheme which makes this method much more efficient than finite-element methods from the point of view of both memory storage and computational time requirements. This work is divided into two parts. The first part mainly focuses on theoretical studies with a simple numerical result for the Che-byshev spectral element, and the second part, mainly with the Legendre spectral element, will give the algorithm implementation, numerical accuracy and efficiency analyses, and then the detailed modeling example comparisons of the proposed spectral element method with a pseudo-spectral method, which will be seen in another work by Lin, Wang and Zhang.     

求解结构动力响应的卷积型DQ半解析法

通过卷积将原始控制方程构造成包含初始条件的新的具有完整初值问题特征的控制方程.该方程既与Gurtin变分原理一样有合理的数学内涵,又避免了卷积型Gurtin变分原理泛函和计算的繁复.对新的控制方程在时间域取解析函数,在空间域采用离散的DQ法,经对梁的动力响应问题的计算表明,该方法是一种精度好效率高的求解动力响应问题的计算方法.     

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.     

A Stochastic Collocation Method for Delay Differential Equations with Random Input

In this work, we concern with the numerical approach for delay differential equations with random coefficients. We first show that the exact solution of the problem considered admits good regularity in the random space, provided that the given data satisfy some reasonable assumptions. A stochastic collocation method is proposed to approximate the solution in the random space, and we use the Legendre spectral collocation method to solve the resulting deterministic delay differential equations. Convergence property of the proposed method is analyzed. It is shown that the numerical method yields the familiar exponential order of convergence in both the random space and the time space. Numerical examples are given to illustrate the theoretical results.     

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.     

Rapid numerical difference recurrent formula of nonlinear Schrödinger equation and its application

In this paper, a rapid numerical difference recurrent formula, in which it has been taken that the chromatic dispersion and the nonlinearity act together along each fiber segment, is established in the time domain by applying a Maclaurin expansion to the differential form of the nonlinear Schrödinger equation (NLSE) in the frequency domain. The calculated results by using the established formula are contrasted with the known analytical results and the results of the split-step Fourier method (SSFM) and indicated that the rapid numerical difference recurrent formula is very accurate and more reasonable because it abandons an assumption that the dispersive and nonlinear effects can be assumed to act independently as the optical field propagates over each fiber segment. It has been concluded that the established formula in this paper is a scientific, reasonable and effective numerical method for the study of light pulse propagation in a nonlinear optical medium.     

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.     

Two-dimensional thermoelasticity solution for functionally graded thick beams

Two-dimensional thermoelasticity analysis of functionally graded thick beams is presented using the state space method coupled with the technique of differential quadrature. Material properties vary continuously and smoothly through the beam thickness, leading to variable coefficients in the state equation derived from the elasticity equations. Approximate laminate model is employed to translate the state equation into the one with constant coefficients in each layer. To avoid numerical instability, joint coupling matrices are introduced according to the continuity conditions at interfaces in the approximate model. The differential quadrature procedure is applied to discretizing the beam in the axial direction to make easy the treatment of arbitrary end conditions. A simply-supported beam with exponentially varying material properties is considered to validate the present method. Numerical examples are performed to investigate the influences of relative parameters.     

求解对流-扩散-反应问题的改进有限积分法

将求解偏微分方程的有限积分法应用于对流-扩散-反应问题,发现对于非对流占优的对流扩散问题,有限积分法的精度比QUICK法高一个数量级,比传统的有限体积法高两个数量级.处理对流占优的对流-扩散-反应问题时,对流项的离散时引进加权参数,通过调节该参数反映输运的方向性.结果表明这种改进的有限积分法的精度比传统的有限体积法至少高四个数量级,同时明显改进了原来的有限积分法的精度和稳定性.对于对流占优的对流-扩散-反应问题,即使采用粗网格,计算结果也未出现非物理振荡现象,表明改进的有限积分法具有很好的稳定性.     

A study on the multi-freedom broadband impedance model for time-domain simulations

The purpose of this paper is to construct a general broadband impedance model, which is suited for predicting acoustic propagation problems in time domain.A multi-freedom broadband impedance model for sound propagation over impedance surfaces is proposed and the corresponding time domain impedance boundary condition is presented.Basing on the extended Helmholtz resonator,the multi-freedom impedance model is constructed through combing with a sum of rational functions in the form of general complex-conjugate pole-residue pairs and it is proved that the impedance model is well posed.The impedance boundary condition can be implemented into a computational aeroacoustics solver by a recursive convolution technique, which results in a fast and computationally efficient algorithm.The two dimensional and three dimensional benchmark problems are selected to validate the accuracy of the proposed impedance model and time domain simulations.The numerical results are in good agreement with the reference solutions.It is demonstrated that the proposed impedance model can be used to describe the broadband characteristics of acoustic liners,and the corresponding time domain impedance boundary condition is viable and accurate for the prediction of sound propagation over broadband impedance surfaces.