二维US-FDTD方法的数值稳定和色散分析

研究了二维US-FDTD方法的数值稳定性和数值色散特性。通过对增长因子的计算,证明了US-FDTD方法的无条件稳定性。利用增长因子的相位,推导出了US-FDTD方法的数值色散关系式。分析了US-FDTD方法的数值色散误差。数值分析表明,与ADI-FDTD方法一样,数值色散误差仍然是决定US-FDTD时间步长选取的关键因素。同时发现,数值色散受时间步长及网格大小的影响。

Second-order difference scheme for a nonlinear model of wood drying process

A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordinary equation. A difference scheme is derived by the method of reduction of order. First, a new variable is introduced and the original problem is rewritten into a system of the first-order differential equations. Secondly, a difference scheme is constructed for the later problem. The solvability, stability and convergence of the difference scheme are proved by the energy method. The convergence order of the difference scheme is secondorder both in time and in space. A prior error estimate is put forward. The new variable is put aside to reduce the computational cost. A numerical example testifies the theoretical result.

A Class of Two-level High-order Accuracy Explicit Difference Scheme for Solving 3-D Parabolic Partial Differential Equation

A class of two-level high-order accuracy explicit difference scheme for solving 3-D parabolic P.D.E is constructed. Its truncation error is (Δt2+Δx4) and the stability condition is r=Δt/Δx2=Δt/Δy2=Δt/Δz2≤1/6.

Reinforcement strength reduction in FEM for mechanically stabilized earth structures

The factor of safety of mechanically stabilized earth(MSE) structures can be analyzed either using limit equilibrium method(LEM) or strength reduction method(SRM) in finite element/difference method. In LEM, the strengths of the reinforcement members and soils are reduced with the same factor. While using the SRM, only soil strength is reduced during the calculation of the factor of safety. This causes inconsistence in calculating the factor of safety of the MSE structures. To overcome this, an iteration method is proposed to consider the strength reduction of the reinforcements in SRM. The method is demonstrated by using PLAXIS, a finite element software. The results show that the factor of safety converges after a few iterations. The reduction of strength has different effects on the factor of safety depending on the properties of the reinforcements and the soil, and failure modes.

Kreiss矩阵定理的新证明

熟知,Kreiss矩阵定理在差分法稳定性理论中占有十分重要的位置.定理的证明相当复杂,[2]中收入的证明虽经Morton和Schechter作了适当处理,但是被称为证明核心的(R)?(S)的过程仍然繁琐.本文给出一种简单直观的新证明.在证明回路(A)?(R)?(S)?(H)?(A)中,(A)?(R)和(H)?(A)沿用[2]的证明,一并给出.顺便指出,新证明并不影响[2]中对另一重要定理(Buchanan准则)证明的简化。

UNCONDITIONAL STABLE DIFFERENCE METHODS WITH INTRINSIC PARALLELISM FOR SEMILINEAR PARABOLIC SYSTEMS OF DIVERGENCE TYPE

The general finite difference schemes with intrinsic parallelism for the boundary value problem of the semilinear parabolic system of divergence type with bounded measurable coefficients is studied. By the approach of the discrete functional analysis, the existence and uniqueness of the discrete vector solutions of the nonlinear difference system with intrinsic parallelism are proved. Moreover the unconditional stability of the general difference schemes with intrinsic parallelism justified in the sense of the continuous dependence of the discrete vector solution of the difference schemes on the discrete initial data of the original problems in the discrete W_2^(2,1) (Q△) norms. Finally the convergence of the discrete vector solutions of the certain difference schemes with intrinsic parallelism to the unique generalized solution of the original semilinear parabolic problem is proved.

RESIDUAL A POSTERIORI ERROR ESTIMATE OF A NEW TWO-LEVEL METHOD FOR STEADY NAVIER-STOKES EQUATIONS

Residual-based a posteriori error estimate for conforming finite element solutions of incompressible Navier-Stokes equations, which is computed with a new two-level method that is different from Volker John, is derived. A posteriori error estimate contains additional terms in comparison to the estimate for the solution obtained by the standard finite element method. The importance of the additional terms in the error estimates is investigated by studying their asymptotic behavior. For optimal scaled meshes, these bounds are not of higher order than the convergence of discrete solution. The two-level method aims to solve the nonlinear problem on a coarse grid with less computational work, then to solve the linear problem on a fine grid, which is superior to the usual finite element method solving a similar nonlinear problem on the fine grid.

A Novel Approach with Time-Splitting Spectral Technique for the Coupled Schrdinger–Boussinesq Equations Involving Riesz Fractional Derivative

In the present paper the Riesz fractional coupled Schr6dinger-Boussinesq （S-B） equations have been solved by the time-splitting Fourier spectral （TSFS） method. This proposed technique is utilized for discretizing the Schrodinger like equation and further, a pseudospectral discretization has been employed for the Boussinesq-like equation. Apart from that an implicit finite difference approach has also been proposed to compare the results with the solutions obtained from the time-splitting technique. Furthermore, the time-splitting method is proved to be unconditionally stable. The error norms along with the graphical solutions have also been presented here.

An element decomposition method with variance strain stabilization

An element decomposition method with variance strain stabilization(EDM-VSS) is proposed. In the present EDM-VSS, the quadrilateral element is first divided into four sub-triangular cells, and the local strains in sub-triangular cells are obtained using linear interpolation function. For each quadrilateral element, the strain of the whole quadrilateral is the weighted average value of the local strains, which means only one integration point is adopted to construct the stiffness matrix. The stabilization item of the stiffness matrix is constructed by variance of the local strains, which can eliminate the instability of the one-point integration formulation and largely increase the accuracy of the element. Compared with conventional full integration quadrilateral element, the EDM-VSS achieves more accurate results and expends much lower computational cost. More importantly, as no mapping or coordinate transformation is involved in the present EDM-VSS, the restriction on the conventional quadrilateral elements can be removed and problem domain can be discretized in more flexible ways. To verify the accuracy and stability of the present formulation, a number of numerical examples are studied to demonstrate the efficiency of the present EDM-VSS.