Locking-free DGFEM for elasticity problems in polygons

The h-version of the discontinuous Galerkin finite element method(h-DGFEM) for nearly incompressible linear elasticity problemsin polygons is analysed. It is proved that the scheme is robust(locking-free) with respect to volume locking, even in the absenceof H²-regularity of the solution. Furthermore, it is shown thatan appropriate choice of the finite element meshes leads torobust and optimal algebraic convergence rates of the DGFEMeven if the exact solutions do not belong to H².     

Some Spectral Properties of One Sturm-Liouville Type Problem with Discontinuous Weight

We consider a discontinuous weight Sturm-Liouville equation together with eigenparameter dependent boundary conditions and two supplementary transmission conditions at the point of discontinuity. We extend and generalize some approaches and results of the classic regular Sturm-Liouville problems to the similar problems with discontinuities. In particular, we introduce a special Hilbert space formulation in such a way that the problem under consideration can be interpreted as an eigenvalue problem for a suitable selfadjoint operator, construct the Green’s function and resolvent operator, and derive asymptotic formulas for eigenvalues and normalized eigenfunctions.Original Russian Text Copyright © 2005 Mukhtarov O. Sh. and Kadakal M.__________Translated from Sibirskii Matematicheskii Zhurnal, Vol. 46, No. 4, pp. 860–875, July–August, 2005.     

Comparison of V‐cycle multigrid method for cell‐centered finite difference on triangular meshes

We consider a multigrid algorithm (MG) for the cell centered finite difference scheme (CCFD) on general triangular meshes using a new prolongation operator. This prolongation is designed to solve the diffusion equation with strongly discontinuous coefficient as well as with smooth one. We compare our new prolongation with the natural injection and the weighted operator in Kwak, Kwon, and Lee ( 8 ) and the behaviors of these three prolongation are discussed. Numerical experiments show that (i) for smooth problems, the multigrid with our new prolongation is fastest, the next is the weighted prolongation, and the third is the natural injection; and (ii) for nonsmooth problems, our new prolongation is again fastest, the next is the natural injection, and the third is the weighted prolongation. In conclusion, our new prolongation works better than the natural injection and the weighted operator for both smooth and nonsmooth problems. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

系数密集间断边值问题并行计算分析

研究系数密集间断的边值问题，采用降维方法，使得在数值计算上可实现并行计算，并且收敛性的证明变得很简洁．     

The equicontinuity condition for sequences of solution spaces

The verification of an appropriately stated equicontinuity condition for a sequence of solution spaces is one of the two key points in the theory of the Cauchy problem for equations with singular right-hand sides. We obtain a related sufficient condition. Translated fromMatematicheskie Zametki, Vol. 61, No. 3, pp. 416–423, March, 1997. Translated by V. E. Nazaikinskii     

Simulation of the interaction of light with 3-D metallic nanostructures using a proper orthogonal decomposition-Galerkin reduced-order discontinuous Galerkin time-domain method

In this artice, we report on a reduced-order model (ROM) based on the proper orthogonal decomposition (POD) technique for the system of 3-D time-domain Maxwell's equations coupled to a Drude dispersion model, which is employed to describe the interaction of light with nanometer scale metallic structures. By using the singular value decomposition (SVD) method, the POD basis vectors are extracted offline from the snapshots produced by a high order discontinuous Galerkin time-domain (DGTD) solver. With a Galerkin projection and a second order leap-frog (LF₂) time discretization, a discrete ROM is constructed. The stability condition of the ROM is then analyzed. In particular, when the boundary is a perfect electric conductor condition, the global energy of the ROM is bounded, which is consistent with the characteristics of global energy in the DGTD method. It is shown that the ROM based on Galerkin projection can maintain the stability characteristics of the original high dimensional model. Numerical experiments are presented to verify the accuracy, demonstrate the capabilities of the POD-based ROM and assess its efficiency for 3-D nanophotonic problems.     

Uniform Convergence of Multigrid V-Cycle on Adaptively Refined Finite Element Meshes for Elliptic Problems with Discontinuous Coefficients

The multigrid V-cycle methods for adaptive finite element discretizations of two-dimensional elliptic problems with discontinuous coefficients are considered. Under the conditions that the coefficient is quasi-monotone up to a constant and the meshes are locally refined by using the newest vertex bisection algorithm, some uniform convergence results are proved for the standard multigrid V-cycle algorithm with Gauss-Seidel relaxations performed only on new nodes and their immediate neighbours. The multigrid V-cycle algorithm uses $\mathcal{O}(N)$ operations per iteration and is optimal.     

Flux calculation at the interface between two rock types for two-phase flow in porous media

When simulating two-phase flow in porous media, one has to consider the case where there is a discontinuity in the medium. There relative permeabilities and capillary pressure functions may change and we address the problem of calculating the convective part of the numerical flux at the interface between the two rock types. Several solutions are compared.     

DDA与FEM耦合法在分缝重力坝非线性分析中的应用

混凝土重力坝一般设有数条纵缝。纵缝使坝体的连续性遭到破坏，一般的有限单元法很难处理这样的不连续变形问题。本文首先阐述了DDA(Discontinuous Deformation Analysis)与FEM(Finite Element Method)耦合算法的原理，然后采用作者自行研制的DDA与FEM耦合法程序对一实际重力坝进行了详细的计算分析。重点讨论了纵缝的不同设置形式和缝宽大小对坝体变形和应力分布、坝底面接触应力分布以及坝基面应力分布的影响，为工程设计提供依据。     

一类指数型间断函数的嵌入非连续等参有限元法

在嵌入非连续有限元的基本思想下,提出一类附加位移形函数———指数型间断函数,来模拟由于非连续结构,如裂纹和节理,所导致的位移不连续规律,该附加函数是以到间断处的垂直距离为自变量,且随距离的增大而呈指数衰减的函数.指数型间断函数具有在数学上的便于积分和求导的优点,且比阶梯间断函数更能反映实际破裂后的变形情况.本文用弱解形式推导了嵌入非连续有限元格式,编制了二维4节点和三维8节点的嵌入非连续等参有限元程序,并分别给出了算例.算例表明在模拟裂纹追踪时,指数型间断函数的嵌入非连续等参有限元法可行且有效.