Difference approximations for wave equations via finite elements. I. Construction and performance

Many practical applications of wave equations involve media in which there are interfaces, or discontinuities in material properties. The accurate numerical representation of these interfaces is important in mathematical models. One can develop generalizations of standard finite-difference methods that accommodate sharp interfaces by modifying a straightforward finite-element approach. In two space dimensions, these methods yield explicit, 5-point or 9-point difference schemes that accurately capture reflection, transmission, and refraction at interfaces. The approach also extends readily to the simulation of waves in elastic media. A companion article presents an error analysis for the approach. © 1995 John Wiley & Sons, Inc.     

Prolate spheroidal wave functions on a disc—Integration and approximation of two-dimensional bandlimited functions

We consider the problem of integrating and approximating 2D bandlimited functions restricted to a disc by using 2D prolate spheroidal wave functions (PSWFs). We derive a numerical scheme for the evaluation of the 2D PSWFs on a disc, which is the basis for the numerical implementation of the presented quadrature and approximation schemes. Next, we derive a quadrature formula for bandlimited functions restricted to a disc and give a bound on the integration error. We apply this quadrature to derive an approximation scheme for such functions. We prove a bound on the approximation error and present numerical results that demonstrate the effectiveness of the quadrature and approximation schemes.     

Construction of a Mindlin pseudospectral plate element and evaluating efficiency of the element

In this paper, a Mindlin pseudospectral plate element is constructed to perform static, dynamic, and wave propagation analyses of plate-like structures. Chebyshev polynomials are used as basis functions and Chebyshev–Gauss–Lobatto points are used as grid points. Two integration schemes, i.e., Gauss–Legendre quadrature (GLEQ) and Chebyshev points quadrature (CPQ), are employed independently to form the elemental stiffness matrix of the present element. A lumped elemental mass matrix is generated by only using CPQ due to the discrete orthogonality of Chebyshev polynomials and overlapping of the quadrature points with the grid points. This results in a remarkable reduction of numerical operations in solving the equation of motion for being able to use explicit time integration schemes. Numerical calculations are carried out to investigate the influence of the above two numerical integration schemes in the elemental stiffness formation on the accuracy of static and dynamic response analyses. By comparing with the results of ABAQUS, this study shows that CPQ performs slightly better than GLEQ in various plates with different thicknesses, especially in thick plates. Finally, a one dimensional (1D) and a 2D wave propagation problems are used to demonstrate the efficiency of the present Mindlin pseudospectral plate element.     

Numerical integration schemes for space–time hypersingular integrals in energetic Galerkin BEM

Here we consider exterior Neumann wave propagation problems reformulated in terms of space–time hypersingular boundary integral equations. We deal with quadrature schemes required, in the discretization phase, by the energetic Galerkin boundary element method.     

High-order splitting schemes for semilinear evolution equations

We first derive necessary and sufficient stiff order conditions, up to order four, for exponential splitting schemes applied to semilinear evolution equations. The main idea is to identify the local splitting error as a sum of quadrature errors. The order conditions of the quadrature rules then yield the stiff order conditions in an explicit fashion, similarly to that of Runge–Kutta schemes. Furthermore, the derived stiff conditions coincide with the classical non-stiff conditions. Secondly, we propose an abstract convergence analysis, where the linear part of the vector field is assumed to generate a group or a semigroup and the nonlinear part is assumed to be smooth and to satisfy a set of compatibility requirements. Concrete applications include nonlinear wave equations and diffusion-reaction processes. The convergence analysis also extends to the case where the nonlinear flows in the exponential splitting scheme are approximated by a sufficiently accurate one-step method.     

Rational approximations of trigonometric matrices with application to second-order systems of differential equations

We consider the direct treatment of the second-order system of equations y” (t)+ Ay(t) = tf;(t), such as might arise in finite-element or finite-difference semidiscretizations of the wave equation. We develop the exact solution and some three-term recurrences involving trigonometric matrices. We approximate these trigonometric matrices by rational approximants of Padé type and thus develop a two-parameter family of approximation schemes. We analyze the stability behavior and computational complexity of members of this family and isolate four schemes for numerical experimentation, the results of which we tabulate. We single out as particularly effective the classical Stormer-Numerov method and also a new sixth-order scheme.     

New schemes for the finite-element dynamic analysis of piezoelectric devices

New finite-element schemes are proposed for investigating harmonic and non-stationary problems for composite elastic and piezoelectric media. These schemes develop the techniques for the finite-element analysis of piezoelectric structures based on symmetric and partitioned matrix algorithms. In order to take account of attenuation in piezoelectric media, a new model is used which extends the Kelvin model for viscoelastic media. It is shown that this model enables the system of finite-element equations to be split into separate scalar equations. The Newmark scheme in a convenient formulation, which does not explicitly use the velocities and accelerations of the nodal degrees of freedom, is employed for the direct integration with respect to time of the finite-element equations of non-stationary problems. The results of numerical experiments are presented which illustrate the effectiveness of the proposed techniques and their implementation in the ACELAN finite-element software package.     

Quadrature formulas for periodic functions and their application to the boundary element method

Two-dimensional and axisymmetric boundary value problems for the Laplace equation in a domain bounded by a closed smooth contour are considered. The problems are reduced to integral equations with a periodic singular kernel, where the period is equal to the length of the contour. Taking into account the periodicity property, high-order accurate quadrature formulas are applied to the integral operator. As a result, the integral equations are reduced to a system of linear algebraic equations. This substantially simplifies the numerical schemes for solving boundary value problems and considerably improves the accuracy of approximation of the integral operator. The boundaries are specified by analytic functions, and the remainder of the quadrature formulas decreases faster than any power of the integration step size. The examples include the two-dimensional potential inviscid circulation flow past a single blade or a grid of blades; the axisymmetric flow past a torus; and free-surface flow problems, such as wave breakdown, standing waves, and the development of Rayleigh-Taylor instability.     

Efficient evaluation of highly oscillatory acoustic scattering surface integrals

We consider the approximation of some highly oscillatory weakly singular surface integrals, arising from boundary integral methods with smooth global basis functions for solving problems of high frequency acoustic scattering by three-dimensional convex obstacles, described globally in spherical coordinates. As the frequency of the incident wave increases, the performance of standard quadrature schemes deteriorates. Naive application of asymptotic schemes also fails due to the weak singularity. We propose here a new scheme based on a combination of an asymptotic approach and exact treatment of singularities in an appropriate coordinate system. For the case of a spherical scatterer we demonstrate via error analysis and numerical results that, provided the observation point is sufficiently far from the shadow boundary, a high level of accuracy can be achieved with a minimal computational cost.     

Numerical elastic-plastic simulation of interlaminar stresses in a notched angle-ply thermoplastic composite laminate

A thermoplastic angle-ply AS4/PEEK laminate with a hole is considered. The interlaminar stresses along the hole edge at different interfaces under uniaxial extension are investigated. According to the symmetries of the structure and loading, a suitable finite-element model is developed. Utilizing a three-dimensional elastic-plastic finite-element procedure elaborated previously, a finite-element modeling of the interlaminar stresses in a thick angle-ply composite laminate is carried out. Based on the interlaminar stresses obtained, the dangerous locations of delamination initiation are predicted. The results obtained indicate that there is some relationship between the dangerous locations and fiber orientations in the adjacent layers, and it maybe inferred that the critical locations are near the regions where the hole edge is tangent to the fiber orientation in the layers adjacent to the interface. The interlaminar stresses at the same interfaces are not sensible to distances from the midplane of the laminate. Very high interlaminar tensile stresses are found to exist on the hole edge at the +25°/+25° or –25°/–25° interfaces, and delaminations can initiate there first. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 45, No. 3, pp. 427-440, May-June, 2009.     

Approximation of differential eigenvalue problems

To solve a one-dimensional second-order differential eigenvalue problem, we use the finite-element method with numerical integration. We analyze the influence of the quadrature formulas used on the error in the approximate eigenvalues and eigenelements. Theoretical results are illustrated by experiments for a model problem.     

Some numerical quadrature schemes of a non-conforming quadrilateral finite element

Numerical quadrature schemes of a non-conforming finite element method for general second order elliptic problems in two dimensional (2-D) and three dimensional (3-D) space are discussed in this paper. We present and analyze some optimal numerical quadrature schemes. One of the schemes contains only three sampling points, which greatly improves the efficiency of numerical computations. The optimal error estimates are derived by using some traditional approaches and techniques. Lastly, some numerical results are provided to verify our theoretical analysis.     

On approximation of discontinuous solutions to the Buckley-Leverett equation

In this paper, the Lax-Wendroff and “cabaret” schemes for the Buckley-Leverett equation are studied. It is shown that these schemes represent unstable solutions. The choice of an unstable solution depends on the Courant number only. A finite-element version of the “cabaret” scheme is given.     

求解时间分布阶扩散方程的两个高阶有限差分格式

基于复化Simpson公式和复化两点Gauss-Legendre公式，构造了两个求解时间分布阶扩散方程的高阶有限差分格式.不同于以往文献中提出的时间一阶或二阶格式，这两种格式在时间方向都具有三阶精度，而在分布阶和空间方向可达到四阶精度.数值结果表明，两种算法都是稳定且收敛的，从而是有效的.两种格式的收敛速率也通过数值实验进行了验证，并且通过和文献中的算法对比可以得出其更为高效,     

Non-axisymmetric consolidation of poroelastic multilayered materials with anisotropic permeability and compressible constituents

This paper presents an analytical layer-element solution to non-axisymmetric consolidation of multilayered poroelastic materials with anisotropic permeability and compressible constituents. By applying Fourier expansions, Hankel transforms and Laplace transforms to the state variables involved in the governing equations of poroelasticity with respect to the circumferential, radial and time coordinates, respectively, the analytical layer-element (i.e. a symmetric stiffness matrix) is derived, which describes the relationship between the transformed generalized stresses and displacements at the surface (z = 0) and those at an arbitrary depth z, considering the corresponding boundary and continuity conditions at the layer interfaces, the global stiffness matrix of a multilayered system is assembled in the transformed domain. The actual solutions in the physical domain are acquired by applying numerical quadrature schemes for the inversion of the Laplace–Hankel transform. Finally, numerical calculation is presented to investigate the influence of layering and poroelastic material parameters on consolidation process.     

Dissipative/conservative Galerkin method using discrete partial derivatives for nonlinear evolution equations

A new method is proposed for designing Galerkin schemes that retain the energy dissipation or conservation properties of nonlinear evolution equations such as the Cahn–Hilliard equation, the Korteweg–de Vries equation, or the nonlinear Schrödinger equation. In particular, as a special case, dissipative or conservative finite-element schemes can be derived. The key device there is the new concept of discrete partial derivatives. As examples of the application of the present method, dissipative or conservative Galerkin schemes are presented for the three equations with some numerical experiments.     

Mathematical simulation of deformation and wave processes in multilayered structures

The deformation and wave processes induced by collisions of an impactor with deformable layered targets of various configurations are analyzed. The numerical solution of such problems is associated with an adequate treatment of wave processes in a continuous medium, which is an especially difficult task in the case of layered targets. To deal with the former problem, it is proposed to use adaptive Lagrangian triangular meshes. Wave processes are simulated using the grid-characteristic method, which can serve as a basis for algorithms that do not fail near the boundary of the computational domain and at numerous material interfaces. Additionally, hybrid and hybridized grid-characteristic schemes are applied that substantially improve numerical solutions with steep gradients (discontinuous solutions). These methods provide an adequate treatment of wave processes in layered targets (wave reflection and refraction at contact surfaces, secondary-wave interaction, changes in the conditions on these boundaries, etc.).     

Approximate solution of singular integro- differential equations in generalized Hölder spaces

We have elaborated the numerical schemes of collocation methods and mechanical quadrature methods for approximate solution of singular integro- differential equations with kernels of Cauchy type. The equations are defined on the arbitrary smooth closed contours of complex plane. The researched methods are based on Fejér points. Theoretical background of collocation methods and mechanical quadrature methods has been obtained in Generalized Hölder spaces.     

Two very accurate and efficient methods for computing eigenvalues of Sturm–Liouville problems

In this paper, we discuss and compare two useful schemes for Sturm–Liouville eigenvalue problems: differential quadrature method (DQM) and collocation method with sinc functions. These methods are then tested on a few examples and a comparison is made. It is shown that the sinc-collocation method in many instances gives better results.     

Compact high-order finite-element method for elliptic transport problems with variable coefficients

We present a new conforming bilinear Petrov-Galerkin finite-element scheme for elliptic transport problems with variable coefficients. This scheme combines a generalized test function and artificial diffusion to achieve O(h⁴) grid-point accuracy on uniform stencils of 3 × 3 in two dimensions without resorting to the extended stencils of high-order elements. The method is compared with upwind and high-order finite-difference schemes and the standard Galerkin finite-element method for representative test problems. © 1994 John Wiley & Sons, Inc.