Non-stationary subdivision for exponential polynomials reproduction

In this paper we develop a novel approach to construct non-stationary subdivision schemes with a tension control parameter which can reproduce functions in a finite-dimensional subspace of exponential polynomials. The construction process is mainly implemented by solving linear systems for primal and dual subdivision schemes respectively, which are based on different parameterizations. We give the theoretical basis for the existence, uniqueness, and refinement rules of schemes proposed in this paper. The convergence and smoothness of the schemes are analyzed as well. Moreover, conics reproducing schemes are analyzed based on our theory, and a new idea that the tensor parameter ωk of the schemes can be adjusted for conics generation is proposed.     

New numerical methods for the coupled nonlinear Schrödinger equations

In this paper, three numerical schemes with high accuracy for the coupled Schrodinger equations are studied. The conserwtive properties of the schemes are obtained and the plane wave solution is analysised. The split step Runge-Kutta scheme is conditionally stable by linearized analyzed. The split step compact scheme and the split step spectral method are unconditionally stable. The trunction error of the schemes are discussed. The fusion of two solitions colliding with different β is shown in the figures. The numerical experments demonstrate that our algorithms are effective and reliable.     

Concatenating construction of the multisymplectic schemes for 2+1-dimensional sine-Gordon equation

In this paper, taking the 2+1-dimensional sine-Gordon equation as an example, we present the concatenating method to construct the multisymplectic schemes. The-method is to discretizee independently the PDEs in different directions with symplectic schemes, so that the multisymplectic schemes can be constructed by concatenating those symplectic schemes. By this method, we can construct multisymplectic schemes, including some widely used schemes with an accuracy of any order. The numerical simulation on the collisions of solitons are also proposed to illustrate the efficiency of the multisymplectic schemes.     

Parallel difference schemes with interface extrapolation terms for quasi-linear parabolic systems

In this paper some new parallel difference schemes with interface extrapolation terms for a quasi-linear parabolic system of equations are constructed. Two types of time extrapolations are proposed to give the interface values on the interface of sub-domains or the values adjacent to the interface points, so that the unconditional stable parallel schemes with the second accuracy are formed. Without assuming heuristically that the original boundary value problem has the unique smooth vector solution, the existence and uniqueness of the discrete vector solutions of the parallel difference schemes constructed are proved. Moreover the unconditional stability of the parallel difference schemes is justified in the sense of the continuous dependence of the discrete vector solution of the schemes on the discrete known data of the original problems in the discrete W2(2,1) (Q△) norms. Finally the convergence of the discrete vector solutions of the parallel difference schemes with interface extrapolation terms to the unique generalized solution of the original quasi-linear parabolic problem is proved. Numerical results are presented to show the good performance of the parallel schemes, including the unconditional stability, the second accuracy and the high parallelism.     

一维p-Laplacian的加权特征值的最佳估计

In this paper, we determine the infimum and the supremum of the Dirich-let eigenvalues λn(p) (n = 1,2,…)of the problem t∈ ?[0,T], where 1 < p < ∞, and the weights p are nonnegative and are subject to conditions p(t)dt = M and max(e[0,T] p(t) = H. It is also explained for whatweights p the infimum and the supremum will be attained.     

Global Rank Axioms for Poset Matroids

An excellent introduction to the topic of poset matroids is due to Barnabei, Nicoletti and Pezzoli. In this paper, we investigate the rank axioms for poset matroids; thereby we can characterize poset matroids in a “global” version and a “pseudo-global” version. Some corresponding properties of combinatorial schemes are also obtained.     

EQUIVALENCE OF SEMI-LAGRANGIAN AND LAGRANGE-GALERKIN SCHEMES UNDER CONSTANT ADVECTION SPEED

We compare in this paper two major implementations of large time-step schemes for advection equations, i.e., Semi-Lagrangian and Lagrange-Galerkin techniques. We show that SL schemes are equivalent to exact LG schemes via a suitable definition of the basis functions. In this paper, this equivalence will be proved assuming some simplifying hypoteses, mainly constant advection speed, uniform space grid, symmetry and translation invariance of the cardinal basis functions for interpolation. As a byproduct of this equivalence, we obtain a simpler proof of stability for SL schemes in the constant-coefficient case.     

COMPACT FOURTH-ORDER FINITE DIFFERENCE SCHEMES FOR HELMHOLTZ EQUATION WITH HIGH WAVE NUMBERS

In this paper, two fourth-order accurate compact difference schemes are presented for solving the Helmholtz equation in two space dimensions when the corresponding wave numbers are large. The main idea is to derive and to study a fourth-order accurate compact difference scheme whose leading truncation term, namely, the O（h^4） term, is independent of the wave number and the solution of the Helmholtz equation. The convergence property of the compact schemes are analyzed and the implementation of solving the resulting linear algebraic system based on a FFT approach is considered. Numerical results are presented, which support our theoretical predictions.     

EXPLICITLY SOLVABLE CHARACTERISTIC METHOD FOR CONVECTION DIFFUSION PROBLEMS

In this paper, a new method for numerically solving nonlinear convection-dominated diffusion problems is devised and analysed. The discrete time approximations with time stepping along charactcristics are cstablished and solved in spaces posscssing reproducing kernel functions. At each time step, the exact solution of the approximate problem is given by explicit expression. The computational advantage of this method is that the schemes are absolutely stable, and are explicitly solvable as well. The stability and error estimates are derived. Some numerical results are given.     

HOMOTOPY CONTINUATION METHODS FOR STOCHASTIC TWO-POINT BOUNDARY VALUE PROBLEMS DRIVEN BY ADDITIVE NOISES

In this paper the homotopy continuation method for stochastic two-point boundary value problems driven by additive noises is studied. The existence of the solution of the homotopy equation is proved. Numerical schemes are constructed and error estimates are obtained. Numerical experiments demonstrate the effectiveness of the homotopy continu- ation method over other commonly used methods such as the shooting method.     

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.     

Semi-implicit Runge-Kutta schemes for the Navier-Stokes equations

The stationary Navier-Stokes equations are solved in 2D with semi-implicit Runge-Kutta schemes, where explicit time-integration in the streamwise direction is combined with implicit integration in the body-normal direction. For model problems stability restrictions and convergence properties are studied. Numerical experiments for the flow over a flat plate show that the number of iterations for the semi-implicit schemes is almost independent of the Reynolds number.     

Local and parallel finite element method for solving the biharmonic eigenvalue problem of plate vibration

In this paper, we establish a new local and parallel finite element discrete scheme based on the shifted‐inverse power method for solving the biharmonic eigenvalue problem of plate vibration. We prove the local error estimation of finite element solution for the biharmonic equation/eigenvalue problem and prove the error estimation of approximate solution obtained by the local and parallel scheme. When the diameters of three grids satisfy H⁴ = ?(w²) = ?(h), the approximate solutions obtained by our schemes can achieve the asymptotically optimal accuracy. The numerical experiments show that the computational schemes proposed in this paper are effective to solve the biharmonic eigenvalue problem of plate vibration.     

Conservative numerical methods for the Full von Kármán plate equations

This article is concerned with the numerical solution of the full dynamical von Kármán plate equations for geometrically nonlinear (large‐amplitude) vibration in the simple case of a rectangular plate under periodic boundary conditions. This system is composed of three equations describing the time evolution of the transverse displacement field, as well as the two longitudinal displacements. Particular emphasis is put on developing a family of numerical schemes which, when losses are absent, are exactly energy conserving. The methodology thus extends previous work on the simple von Kármán system, for which longitudinal inertia effects are neglected, resulting in a set of two equations for the transverse displacement and an Airy stress function. Both the semidiscrete (in time) and fully discrete schemes are developed. From the numerical energy conservation property, it is possible to arrive at sufficient conditions for numerical stability, under strongly nonlinear conditions. Simulation results are presented, illustrating various features of plate vibration at high amplitudes, as well as the numerical energy conservation property, using both simple finite difference as well as Fourier spectral discretizations. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1948–1970, 2015     

Multiscale finite element discretizations based on local defect correction for the biharmonic eigenvalue problem of plate buckling

In this paper, we study the multiscale finite element discretizations about the biharmonic eigenvalue problem of plate buckling. On the basis of the work of Dai and Zhou (SIAM J. Numer. Anal. 46[1] [2008] 295‐324), we establish a three‐scale scheme, a multiscale discretization scheme, and the associated parallel version based on local defect correction. We first prove a local priori error estimate of finite element approximations, then give the error estimates of multiscale discretization schemes. Theoretical analysis and numerical experiments indicate that our schemes are suitable and efficient for eigenfunctions with local low smoothness.     

Solution of some plate bending problems using the boundary element method

Some fundamental aspects of the boundary element method of the Kirchhoff theory of thin plate flexure are given. The direct boundary integral equation method with higher conforming properties (using first-order Hermitian interpolation for plate displacement ω, and zero-order Hermitian interpolation for angle of rotation θ, the moment M andthe equivalent shear V) are used for several computational examples. They are: square plate with simply-supported or clamped edges, the same square plate with square central opening and the cantilevered triangular plates. The results of computation as compared with some known experimental and theoritical results showed that the numerical schemes seemed to be satisfactory for the practical applications.     

Analysis of some low-order finite element schemes for Mindlin-Reissner and Kirchhoff plates

Summary We set up a framework for analyzing mixed finite element methods for the plate problem using a mesh dependent energy norm which applies both to the Kirchhoff and to the Mindlin-Reissner formulation of the problem. The analysis techniques are applied to some low order finite element schemes where three degrees of freedom are associated to each vertex of a triangulation of the domain. The schemes proceed from the Mindlin-Reissner formulation with modified shear energy.Dedicated to Professor Ivo Babuka on the occasion of his 60th birthday     

High performance computation for DNS/LES

The paper is focused on high-order compact schemes for direct numerical simulation (DNS) and large eddy simulation (LES) for flow separation, transition, tip vortex, and flow control. A discussion is given for several fundamental issues such as high quality grid generation, high-order schemes for curvilinear coordinates, the CFL condition for complex geometry, and high-order weighted compact schemes for shock capturing and shock–vortex interaction. The computation examples include DNS for K-type and H-type transition, DNS for flow separation and transition around an airfoil with attack angle, control of flow separation by using pulsed jets, and LES simulation for a tip vortex behind the juncture of a wing and flat plate. The computation also shows an almost linear growth in efficiency obtained by using multiple processors.     

弹性约束边界条件下矩形蜂窝夹芯板的自由振动分析

蜂窝夹芯板在飞行器、高速列车等领域有广泛的用途，对其开展振动分析具有明确的科学价值及工程意义.为区别于诸简支等传统约束边界，提出了弹性约束边界下蜂窝夹芯板结构的自由振动特性分析方法.具体来说，首先通过将蜂窝夹芯层等效为各向异性板，将夹芯板问题转变为三层板结构.进一步地，将板结构的位移场函数由改进的二维Fourier级数表示，并基于能量原理的Rayleigh Ritz法得到结构的固有频率和固有振型，理论预测结果与数值模拟分析吻合较好.提出的理论模型可用于系统讨论约束边界对蜂窝夹芯结构自由振动特性的影响，为此类结构的约束方案设计提供理论依据.     

An efficient numerical method for the resolution of the Kirchhoff‐Love dynamic plate equation

We solve numerically the Kirchhoff‐Love dynamic plate equation for an anisotropic heterogeneous material using a spectral method. A mixed velocity‐moment formulation is proposed for the space approximation allowing the use of classical Lagrange finite elements. The benefit of using high order elements is shown through a numerical dispersion analysis. The system resulting from this spatial discretization is solved analytically. Hence this method is particularly efficient for long duration experiments. This time evolution method is compared with explicit and implicit finite differences schemes in terms of accuracy and computation time. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005