Problems on collinear cracks between bonded dissimilar materials under concentrated loads

Following Ref. [6], this paper deals with the problem on collinear cracks between bonded dissimilar materials under a concentrated force and moment at an arbitrary point. Several typical solutions of complex stress functions in closed form are formulated and the stress intensity factors are given. These solutions include a series of results of previous researchers, and redress some errors in the researches of problems containing semi-infinite cracks^[3],[4].     

A dynamic injection operator in a multigrid solution of convection-diffusion equations

A multigrid method is studied for the solution of a linear system resulting from the high-order nine-point discretization of the convection-diffusion equations. The residual injection operator is used as a substitute for the usual full-weighting in the multigrid process. A heuristic analysis is given to obtain a dynamic injection operator that is cost-effective for both diffusion- and convection-dominated problems. Numerical experiments are employed to test the stability and efficiency of the proposed method. © 1998 John Wiley & Sons, Ltd.     

瞬态对流-扩散方程的变分多尺度解法

根据变分多尺度思想,求解了瞬态线性和非线性对流-扩散方程。文中为了简化细尺度方程的求解,忽略了该方程的瞬态性,分别用高阶多项式泡函数(High-order Polynomial Bubble)和自由残量泡(Residual Free Bubble)函数近似细尺度解,进而引入了消除数值伪振荡的稳定化结构。数值算例验证了本文方法的精确性、稳定性和对高Peclet数问题的适应性,证明了上述对细尺度模型的简化是可行的。     

Characteristic Galerkin method for convection-diffusion equations and implicit algorithm using precise integration

This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability. The project sponsored by the State Scientific and Technological Commission of China through “China State Key Project: the Theory and Methodology for Scientific and Engineering Computations with Large Scale”, the National Natural Science Foundation of China and the European Commission Research Project CI1^*CT94-0014.     

Analysis of chebyshev pseudospectral method for multi-dimensional generalized srlw equations

IntroductionAsymmetricregularizedlongwaveequation (SRLWE) 2 x2 -1 u t = x ρ+ 12 u2 ,　 ρ t+ u x=0 (1 )hasbeeninvestigatedinRef.[1 ] .Thesystem (1 )ofequationsisshowntodescribeweaklynonlinearion_acousticwaveandspace_chargewaves.Thehuperbolicsecantsquaredsolitarywaves ,thefourconservationlaws,andsomenumericalresultshavebeenobtainedinRef.[1 ] .Obviously ,eliminatingρin (1 ) ,weobtainaclassofregularizedlongwaveequationutt-uxx+ 12 u2xt-uxxtt =0 . (2 )TheSRLWequationisexplicitlysymmetric…     

H 1 space-time discontinuous finite element method for convection-diffusion equations

An H~1 space-time discontinuous Galerkin (STDG) scheme for convectiondiffusion equations in one spatial dimension is constructed and analyzed. This method is formulated by combining the H~1 Galerkin method and the space-time discontinuous finite element method that is discontinuous in time and continuous in space. The existence and the uniqueness of the approximate solution are proved. The convergence of the scheme is analyzed by using the techniques in the finite difference and finite element methods. An optimal a-priori error estimate in the L~∞ (H~1 ) norm is derived. The numerical exper- iments are presented to verify the theoretical results.     

Local and parallel finite element algorithms for time-dependent convection-diffusion equations

Local and parallel finite element algorithms based on two-grid discretization for the time-dependent convection-diffusion equations are presented. These algorithms are motivated by the observation that, for a solution to the convection-diffusion problem, low frequency components can be approximated well by a relatively coarse grid and high frequency components can be computed on a fine grid by some local and parallel proce- dures. Hence, these local and parallel algorithms only involve one small original problem on the coarse mesh and some correction problems on the local fine grid. One technical tool for the analysis is the local a priori estimates that are also obtained. Some numerical examples are given to support our theoretical analvsis.     

An improved projection scheme applied to pseudospectral methods for the incompressible Navier–Stokes equations

An improved projection scheme is proposed and applied to pseudospectral collocation-Chebyshev approximation for the incompressible Navier–Stokes equations. It consists of introducing a correct predictor for the pressure, one which is consistent with a divergence-free velocity field at each time step. The main objective is to allow a time variation of the pressure gradient at boundaries. From different test problems, it is shown that this method, associated with a multistep second-order time scheme, provides a time accuracy of the same order as the temporal scheme used for the pressure, and also improves the prediction of the velocity slip. Moreover, it does not exhibit any numerical boundary layer mentioned as a drawback of fractional steps algorithm, and does not require the use of staggered grids for the velocity and the pressure. Its effectiveness is validated by comparison with a previous time-splitting algorithm proposed by Goda (K. Goda, J. Comput. Phys., 30 , 76–95 (1979)) and implemented by Gresho (P. Gresho, Int. j. numer. methods fluids, 11 , 587–620 (1990)) to finite element approximations. Steady and unsteady solutions for the regularized driven cavity and the rotating cavity submitted to throughflow are also used to assess the efficiency of this algorithm. © 1998 John Wiley & Sons, Ltd.     

一种基于Associated Hermite正交函数求解对流扩散方程的算法

提出了一种基于AH(Associated Hermite)正交基函数求解对流扩散方程的无条件稳定算法。该算法将方程的时间项通过Hermite多项式作为正交基函数进行展开,利用Galerkin方法消除时间变量项,从而导出有限维AH域隐式差分方程,突破了传统显式差分格式稳定性条件的限制,最后通过对AH域展开系数的求解得到该对流扩散方程的数值解。在数值算例中,将该算法与传统显示差分法和交替方向隐式差分法进行对比分析,数值计算结果表明,算法无条件稳定且其计算精度与时间步长无关,对于具有精细结构的对流换热问题,该算法具有明显的效率优势,且保持了较高的精度。     

Modified Laguerre spectral and pseudospectral methods for nonlinear partial differential equations in multiple dimensions

The Laguerre spectral and pseudospectral methods are investigated for multidimensional nonlinear partial differential equations. Some results on the modified Laguerre orthogonal approximation and interpolation are established, which play important roles in the related numerical methods for unbounded domains. As an example, the modified Laguerre spectral and pseudospectral methods are proposed for two-dimensional Logistic equation. The stability and convergence of the suggested schemes are proved. Numerical results demonstrate the high accuracy of these approaches.     

Direct spectral domain decomposition method for 2D incompressible Navier-Stokes equations

An efficient direct spectral domain decomposition method is developed coupled with Chebyshev spectral approximation for the solution of 2D, unsteady and incompressible Navier-Stokes equations in complex geometries. In this numerical approach, the spatial domains of interest are decomposed into several non-overlapping rectangular sub-domains. In each sub-domain, an improved projection scheme with second-order accuracy is used to deal with the coupling of velocity and pressure, and the Chebyshev collocation spectral method (CSM) is adopted to execute the spatial discretization. The influence matrix technique is employed to enforce the continuities of both variables and their normal derivatives between the adjacent sub-domains. The imposing of the Neumann boundary conditions to the Poisson equations of pressure and intermediate variable will result in the indeterminate solution. A new strategy of assuming the Dirichlet boundary conditions on interface and using the first-order normal derivatives as transmission conditions to keep the continuities of variables is proposed to overcome this trouble. Three test cases are used to verify the accuracy and efficiency, and the detailed comparison between the numerical results and the available solutions is done. The results indicate that the present method is efficiency, stability, and accuracy.     

Multidomain collocation techniques for a viscous compressible flow

This paper presents a viscous compressible flow problem to which an equilibrium solution, in terms of density and velocity, can be given implicitly by elementary functions. The corresponding initial boundary value problem is solved by time discretization by the Crank-Nicolson method, Newton linearization and space discretization using multidomain Chebyshev collocation techniques. The physical interval is covered by subintervals of equal length. Each subinterval utilizes the same number of collocation points and each interface consists of one or two points. Six ways of patching are tested. All of them yield solutions with spectral accuracy for a few time steps, but only three are stable in the long run. Details of the density evolution are illustrated.     

Comparison of three spectral methods for the Benjamin-Ono equation: Fourier pseudospectral, rational Christov functions and Gaussian radial basis functions

The Benjamin-Ono equation is especially challenging for numerical methods because (i) it contains the Hilbert transform, a nonlocal integral operator, and (ii) its solitary waves decay only as O(1/|x|²). We compare three different spectral methods for solving this one-space-dimensional equation. The Fourier pseudospectral method is very fast through use of the Fast Fourier Transform (FFT), but requires domain truncation: replacement of the infinite interval by a large but finite domain. Such truncation is unnecessary for a rational basis, but it is simple to evaluate the Hilbert Transform only when the usual rational Chebyshev functions TB_n(x) are replaced by their cousins, the Christov functions; the FFT still applies. Radial basis functions (RBFs) are slow for a given number of grid points N because of the absence of a summation algorithm as fast as the FFT; because RBFs are meshless, however, very flexible grid adaptation is possible.     

An analysis of overlapping Schwarz method for a weakly coupled system of singularly perturbed convection-diffusion equations

In this article, we have developed an overlapping Schwarz method for a weakly coupled system of convection-diffusion equations. The method splits the original domain into two overlapping subdomains. A hybrid difference scheme is proposed in which on the boundary layer region, we use the central finite difference scheme on a uniform mesh, whereas on the nonlayer region, we use the mid-point difference scheme on a uniform mesh. It is shown that the numerical approximations converge in the maximum norm to the exact solution. We have proved that, when appropriate subdomains are used, the method produces almost second-order convergence. Furthermore, it is shown that two iterations are sufficient to achieve the expected accuracy. Numerical examples are presented to support the theoretical results. The main advantage of this method used with the proposed scheme is that it reduces iteration counts very much and easily identifies in which iteration the Schwarz iterate terminates.     

TWO-GRID METHOD FOR CHARACTERISTICS FINITE-ELEMENT SOLUTION OF 2D NONLINEAR CONVECTION-DOMINATED DIFFUSION PROBLEM

Introduction Convection diffusionequationgovernssuchphenomenaastheflowofheatwithina movingfluid,thetransportofdissolvednutrientsorcontaminantswithinthegroundwater,andthetransportofasurfactantortracerwithinanincompressibleoilinapetroleum reservoir.Weconsid…     

Mechanical quadrature methods and extrapolation for solving nonlinear boundary Helmholtz integral equations

This paper presents mechanical quadrature methods (MQMs) for solving nonlinear boundary Helmholtz integral equations. The methods have high accuracy of order O(h ³) and low computation complexity. Moreover, the mechanical quadrature methods are simple without computing any singular integration. A nonlinear system is constructed by discretizing the nonlinear boundary integral equations. The stability and convergence of the system are proved based on an asymptotical compact theory and the Stepleman theorem. Using the h ³-Richardson extrapolation algorithms (EAs), the accuracy to the order of O(h ⁵) is improved. To slove the nonlinear system, the Newton iteration is discussed extensively by using the Ostrowski fixed point theorem. The efficiency of the algorithms is illustrated by numerical examples.     

Taylor expansion method for nonlinear evolution equations

Introduction Thestudyofnonlinearevolutionequationsisafascinatingproblemwhichisattheveryheart oftheunderstandingofmanyimportantproblemsinthenaturalsciences[1,2].Thenonlinear evolutionequationsandtheirnumericalapproximationareveryimportantintheareasof theoreticalmathematicsandcomputationalmathematics.Aninterestingfeatureofthe approximationtheoryofthenonlinearevolutionequationsistheapplicationsofthefunctional analyticmethodstothenumericalapproximationofthenonlinearevolutionequations. Thispaperist…     

An approach for choosing discretization schemes and grid size based on the convection-diffusion equation

A new approach for selecting proper discretization schemes and grid size is presented. This method is based on the convection-diffusion equation and can provide insight for the Navier-Stokes equation. The approach mainly addresses two aspects, i.e., the practical accuracy of diffusion term discretization and the behavior of high wavenumber disturbances. Two criteria are included in this approach. First, numerical diffusion should not affect the theoretical diffusion accuracy near the length scales of interest. This is achieved by requiring numerical diffusion to be smaller than the diffusion discretization error. Second, high wavenumber modes that are much smaller than the length scales of interest should not be amplified. These two criteria provide a range of suitable scheme combinations for convective flux and diffusive flux and an ideal interval for grid spacing. The effects of time discretization on these criteria are briefly discussed.     

The first integral method for solving some important nonlinear partial differential equations

Exact solutions of some important nonlinear partial differential equations are obtained by using the first integral method. The efficiency of the method is demonstrated by applying it for two selected equations.     

Taylor-Galerkin-based spectral element methods for convection-diffusion problems

Several explicit Taylor-Galerkin-based time integration schemes are proposed for the solution of both linear and non-linear convection problems with divergence-free velocity. These schemes are based on second-order Taylor series of the time derivative. The spatial discretization is performed by a high-order Galerkin spectral element method. For convection-diffusion problems an operator-splitting technique is given that decouples the treatment of the convective and diffusive terms. Both problems are then solved using a suitable time scheme. The Taylor-Galerkin methods and the operator-splitting scheme are tested numerically for both convection and convection-diffusion problems.