谱元方法求解二维不可压缩Navier-Stokes方程

利用有限元的思想并结合谱方法的精度提出求解偏微分方程的谱元方法,在元素内插值函数使用伪谱Chebyshev逼近,并将此方法应用于求解不可压Navier-Stokes方程,具体求解了二维方腔顶盖驱动流,与公认基准解对比获得了较好的结果。     

基本解方法与边界节点法求解Helmholtz方程的比较研究

基本解方法和边界节点法是基于径向基函数的两种重要无网格边界离散数值技术。针对Helmholtz方程,本文比较研究这两种数值方法在不同计算区域问题上的计算精度、插值矩阵对称性、病态性及计算成本。数值试验结果表明,两种方法都可以有效求解边界数据准确的Helmholtz问题。在数值离散过程中,两种方法都可以通过调整配置点的位...     

关于Orr—Sommerfeld方程的Chebyshev谱方法的讨论

本文讨论了Ｏｒｒ－Ｓｏｍｍｅｒｅｌｄ方程的各种Ｃｈｅｂｙｓｈｅｖ谱离散方法，数值证明了Ｃｈｅｂｙｓｈｅｖ配置法离散Ｏｒｒ－Ｓｏｍｍｅｒｆｅｌｄ方程没有伪谱，并以此构造了适于任意平面平行速度剖面情形，对时间和时空稳定性模式一致有效的无伪谱的离离散方法，其中，对时空稳定性问题本文给出了一种新的迭代法可以快速有效地求出复频率的鞍点，对平面Ｐｏｉｓｅｕｉｌｌｅ流，Ｂｌａｓｉｕｓ边界层流和Ｇａｕｓｓ模型尾迹     

通过Adomian分解法求解二维Helmholtz方程

提出基于Adomian分解法求解二维Helmholtz方程。通过Adomian分解法可以把Helmholtz微分方程和边界条件分别转换成递归代数公式和适用符号计算的简单代数公式。利用边界条件可以很容易得到方程的解析解表达式。Adomian分解法的主要特点在于计算简单快速,并且不需要进行线性化或离散化。最后给出数值实例以验证Adomian分解法求解二维Helmholtz方程的有效性。通过数值计算可以发现,基于Adomian分解法的计算结果非常接近精确解,并且该方法具有良好的收敛性。这表明Adomian分解法能够快速有效求解Helmholtz方程。     

通过Adomian分解法求解二维Helmholtz方程

提出基于Adomian分解法求解二维Helmholtz方程。通过Adomian分解法可以把Helmholtz微分方程和边界条件分别转换成递归代数公式和适用符号计算的简单代数公式。利用边界条件可以很容易得到方程的解析解表达式。Adomian分解法的主要特点在于计算简单快速,并且不需要进行线性化或离散化。最后给出数值实例以验证Adomian分解法求解二维Helmholtz方程的有效性。通过数值计算可以发现,基于Adomian分解法的计算结果非常接近精确解,并且该方法具有良好的收敛性。这表明Adomian分解法能够快速有效求解Helmholtz方程。     

无网格介点法求解Helmholtz方程

作为一种配点型无网格法,无网格介点MIP法具有数值实施简单、计算精度高、运算高效和适用范围广等优点。Helmholtz方程是科学与工程问题中广泛应用的一类特殊方程,因此对MIP法求解此类方程的适用性进行了验证。利用MIP法的d适应性,给出了MIP法求解该方程的两种计算格式。在数值算例中,分别对平面规则域和不规则域上的一般Helmholtz方程,以及轴对称Helmholtz方程进行了数值分析。结果表明,MIP法完全适用于求解Helmholtz方程。而且,MIP法的计算精度和收敛性都优于普通配点法。此外,MIP法的两种计算格式中,L2C0型通常具有更好的计算效果,故建议将该计算格式作为MIP法求解该类方程的标准形式。     

谱元方法求解环形空间内自然对流换热

提出了将谱元方法应用到极坐标系下，利用极坐标系下的谱元方法求解环形空间内自然对流问题。具体求解了原始变量速度和压力的不可压缩Navier-Stokes方程和能量方程，通过在时间方向采用时间分裂方法和空间采用谱元方法对方程进行离散求解，取得了与基准解较一致的计算结果。     

广义概率密度演化方程的Chebyshev拟谱法

概率密度演化方法(probability density evolution equation, PDEM)为非线性随机结构的动力响应分析提供了新的途径. 通过PDEM获得结构响应概率密度函数(probability density function, PDF)的关键步骤是求解广义概率密度演化方程(generalized probability density evolution equation, GDEE). 对于GDEE的求解通常采用有限差分法, 然而, 由于GDEE是初始条件间断的变系数一阶双曲偏微分方程, 通过有限差分法求解GDEE可能会面临网格敏感性问题、数值色散和数值耗散现象. 文章从全局逼近的角度出发, 基于Chebyshev拟谱法为GDEE构造了全局插值格式, 解决了数值色散、数值耗散以及网格敏感性问题. 考虑GDEE的系数在每个时间步长均为常数, 推导了GDEE在每一个时间步长内时域上的序列矩阵指数解. 由于序列矩阵指数解形式上是解析的, 从而很好地克服了数值稳定性问题. 两个数值算例表明, 通过Chebyshev拟谱法结合时域的序列矩阵指数解求解GDEE得到的结果与精确解以及Monte Carlo模拟的结果非常吻合, 且数值耗散和数值色散现象几乎可以忽略. 此外, 拟谱法具有高效的收敛性且序列矩阵指数解不受CFL (Courant-Friedrichs-Lewy)条件的限制, 因此该方法具有良好的数值稳定性和计算效率.     

Helmholtz方程的微分容积解法

用一种新型的数值技术－－微分容积法（Differential Cubature Method）求解二维Helmholtz方程的边值问题,几个数值算例表明,该方法稳定收敛,并具有较好的数值精度,本文方法适用于求解具有较小波数的Helmholtz方程。     

极坐标系下的Legendre谱元方法求解Poisson-型方程

r=0处的坐标奇异性是求解极坐标下Poisson-型方程的关键。本文提出一种极坐标系下基于Galerkin变分的Legendre谱元方法用于求解圆形区域内的Poisson-型方程,物理区域的径向和周向划分若干单元,计算单元均采用Legendre多项式展开;圆心所在单元的径向使用LGR(Legendre Gauss Radau)积分点,其他单元径向使用LGL(Legendre Gauss Lobatto)积分点,从而避免了极点处1/r坐标奇异性,周向单元均采用LGL积分点。利用区域分解技术,可以避免节点在极点附近聚集;最后求解了多个Dirichlet或Neumann边界条件下的Poisson-型方程算例。数值结果表明,谱元方法具有很高的精度。     

Dynamic analysis of beam-cable coupled systems using Chebyshev spectral element method

The dynamic characteristics of a beam–cable coupled system are investigated using an improved Chebyshev spectral element method in order to observe the effects of adding cables on the beam. The system is modeled as a double Timoshenko beam system interconnected by discrete springs. Utilizing Chebyshev series expansion and meshing the system according to the locations of its connections,numerical results of the natural frequencies and mode shapes are obtained using only a few elements, and the results are validated by comparing them with the results of a finiteelement method. Then the effects of the cable parameters and layout of connections on the natural frequencies and mode shapes of a fixed-pinned beam are studied. The results show that the modes of a beam–cable coupled system can be classified into two types, beam mode and cable mode, according to the dominant deformation. To avoid undesirable vibrations of the cable, its parameters should be controlled in a reasonable range, or the layout of the connections should be optimized.     

An adaptive cell-based domain integration method for treatment of domain integrals in 3D boundary element method for potential and elasticity problems

An adaptive cell-based domain integration method(CDIM) is proposed for the treatment of domain integrals in 3D boundary element method(BEM). The domain integrals are computed in background cells rather than volume elements. The cells are created from the boundary elements based on an adaptive oct-tree structure and no other discretization is needed. Cells containing the boundary elements are subdivided into smaller sub-cells adaptively according to the sizes and levels of the boundary elements; and the sub-cells outside the domain are deleted to obtain the desired accuracy. The method is applied in the 3D potential and elasticity problems in this paper.     

Chebyshev finite spectral method with extended moving grids

A Chebyshev finite spectral method on non-uniform meshes is proposed. An equidistribution scheme for two types of extended moving grids is used to generate grids. One type is designed to provide better resolution for the wave surface, and the other type is for highly variable gradients. The method has high-order accuracy because of the use of the Chebyshev polynomial as the basis function. The polynomial is used to interpolate the values between the two non-uniform meshes from a previous time step to the current time step. To attain high accuracy in the time discretization, the fourth-order Adams-Bashforth-Moulton predictor and corrector scheme is used. To avoid numerical oscillations caused by the dispersion term in the Korteweg-de Vries (KdV) equation, a numerical technique on non-uniform meshes is introduced. The proposed numerical scheme is validated by the applications to the Burgers equation (nonlinear convectiondiffusion problems) and the KdV equation (single solitary and 2-solitary wave problems), where analytical solutions are available for comparisons. Numerical results agree very well with the corresponding analytical solutions in all cases.     

Generalized mixed finite element method for 3D elasticity problems

Without applying any stable element techniques in the mixed methods, two simple generalized mixed element(GME) formulations were derived by combining the minimum potential energy principle and Hellinger–Reissner(H–R) variational principle. The main features of the GME formulations are that the common C0-continuous polynomial shape functions for displacement methods are used to express both displacement and stress variables, and the coefficient matrix of these formulations is not only automatically symmetric but also invertible. Hence, the numerical results of the generalized mixed methods based on the GME formulations are stable. Displacement as well as stress results can be obtained directly from the algebraic system for finite element analysis after introducing stress and displacement boundary conditions simultaneously. Numerical examples show that displacement and stress results retain the same accuracy. The results of the noncompatible generalized mixed method proposed herein are more accurate than those of the standard noncompatible displacement method. The noncompatible generalized mixed element is less sensitive to element geometric distortions.     

Spectral method for solving the equal width equation based on Chebyshev polynomials

A spectral solution of the equal width (EW) equation based on the collocation method using Chebyshev polynomials as a basis for the approximate solution has been studied. Test problems, including the migration of a single solitary wave with different amplitudes are used to validate this algorithm which is found to be accurate and efficient. The three invariants of the motion are evaluated to determine the conservation properties of the algorithm. The interaction of two solitary waves is seen to cause the creation of a source for solitary waves. Usually these are of small magnitude, but when the amplitudes of the two interacting waves are opposite, the source produces trains of solitary waves whose amplitudes are of the same order as those of the initial waves. The three invariants of the motion of the interaction of the three positive solitary waves are computed to determine the conservation properties of the system. The temporal evaluation of a Maxwellian initial pulse is then studied. Comparisons are made with the most recent results both for the error norms and the invariant values.     

A Helmholtz pressure equation method for the calculation of unsteady incompressibl eviscous flows

A time-implicit numerical method for solving unsteady incompressible viscous flow problems is introduced. The method is based on introducing intermediate compressibility into a projection scheme to obtain a Helmholtz equation for a pressure-type variable. The intermediate compressibility increases the diagonal dominance of the discretized pressure equation so that the Helmholtz pressure equation is relatively easy to solve numerically. The Helmholtz pressure equation provides an iterative method for satisfying the continuity equation for time-implicit Navier–Stokes algorithms. An iterative scheme is used to simultaneously satisfy, within a given tolerance, the velocity divergence-free condition and momentum equations at each time step. Collocated primitive variables on a non-staggered finite difference mesh are used. The method is applied to an unsteady Taylor problem and unsteady laminar flow past a circular cylinder.     

An efficient numerical method for exterior and interior inverse problems of Helmholtz equation

The generalized pulse-spectrum technique (GPST), an efficient and versatile inversion algorithm, is used with adaptive grids to solve both exterior (scattering) and interior (cavity) boundary-shape inverse problems of two-dimensional Helmholtz equation. Numerical simulations of nontrivial examples are carried out to test the feasibility and to study the general characteristics of GPST without the real measurement data. It is found that GPST does efficiently produce very good results.     

Double-grid Chebyshev spectral elements for acoustic wave modeling

Highly accurate algorithms are needed for modeling wave propagation phenomena in realistic media. The spectral element methods, either based on a Chebyshev or a Legendre polynomial basis, have shown their excellent properties of high accuracy and flexibility in describing complex models outperforming other techniques. In contrast with standard grid methods, which use dense spatial meshes, spectral element methods discretize the computational domain in a very coarse mesh. With constant-property elements, this fact may in some cases reduce seriously the computational efficiency. For instance, if the medium is finely heterogeneous, it may need to be described in a much finer way than the acoustic wave field. The double-grid approach presented in this work is a viable way for overcoming this lack of the method and for handling problems where the medium changes continuously or even sharply on the small scale. The variation in the properties is taken into account by using an independent set of shape functions defined on a temporary local grid in such a way that either the small scale fluctuations are accurately handled, without the need of a global finer grid, and the macroscopic wave field propagation is solved with no loose of computational efficiency.