Adaptive finite element analysis of structures under transient dynamic loading using modal superposition

In this paper, the error estimation and adaptive strategy developed for the linear elastodynamic problem under transient dynamic loading based on the Z–Z criterion is utilized for 2D and plate bending problems. An automatic mesh generator based on “growth meshing” is utilized effectively for adaptive mesh refinement. Optimal meshes are obtained iteratively corresponding to the prescribed domain discretization error limit and for a chosen number of basis modes satisfying modal truncation errors. Numerous examples show the effectiveness of the integrated approach in achieving the target accuracy in finite element transient dynamic analysis.     

A dynamic data structure suitable for adaptive mesh refinement in finite element method

This paper describes a dynamic data structure and its implementation, used for an optimum mesh generator. The implementation of this mesh generator was a part of a software package implemented to solve electromagnetic field problems using the finite element method. This mesh generator takes advantage of the Delaunay algorithm, which maximizes the summation of the smallest angles in all triangles and thus creates a mesh that is proved to be an optimum mesh for use in the finite element method. The dynamic data structure is explained and the source code is reviewed. The programs have been written in Pascal programming language.     

A self-organized mesh generator using pattern formation in a reaction–diffusion system

A new type of mesh generator is developed by using a self-organized pattern in a reaction–diffusion system. The system is the Gray–Scott model, which creates a spot pattern in a specific parameter region. The spots correspond to nodes of a mesh. The mesh generator has several advantages: the algorithm is simple and processes to improve the mesh, such as smoothing, (locally) addition, and removal of nodes, are automatically performed by the system.     

Stochastic Domain Decomposition for Time Dependent Adaptive Mesh Generation

The efficient generation of meshes is an important component in the numerical solution of problems in physics and engineering. Of interest are situations where global mesh quality and a tight coupling to the solution of the physical partial differential equation (PDE) is important. We consider parabolic PDE mesh generation and present a method for the construction of adaptive meshes in two spatial dimensions using stochastic domain decomposition that is suitable for an implementation in a multi- or many-core environment. Methods for mesh generation on periodic domains are also provided. The mesh generator is coupled to a time dependent physical PDE and the system is evolved using an alternating solution procedure. The method uses the stochastic representation of the exact solution of a parabolic linear mesh generator to find the location of an adaptive mesh along the (artificial) subdomain interfaces. The deterministic evaluation of the mesh over each subdomain can then be obtained completely independently using the probabilistically computed solutions as boundary conditions. A small scaling study is provided to demonstrate the parallel performance of this stochastic domain decomposition approach to mesh generation. We demonstrate the approach numerically and compare the mesh obtained with the corresponding single domain mesh using a representative mesh quality measure.     

MORTAR型旋转Q_1元的多重网格方法

本文讨论了mortar型旋转Q_1元的多重网格方法.证明了W循环的多重网格法是最优的,即收敛率与网格尺寸及层数无关.同时给出了一种可变的V循环多重网格算法,得到了一个条件数一致有界的预条件子.最后,数值试验验证了我们的理论结果.     

Local residual-type error estimates for adaptive boundary element methods on closed curves

Considering the Galerkin boundary element methods and mesh families with bounded local mesh ratio we derive local error estimates such that the local error is bounded by a local residual together with some global terms which can be expected to be small. If the order of the boundary operator is non-negative and at most two, these estimates show that the local residual is a local error indicator. For the operators of the negative order we obtain the same conclusion if the mesh is β-regular. Our paper improves recent results of Wendland and Yu in several respects.     

Adaptive mesh refinement for the control of cost and quality in finite element analysis

Adaptive strategies are a necessary tool to make finite element analysis applicable to engineering practice. In this paper, attention is restricted to mesh adaptivity. Traditionally, the most common mesh adaptive strategies for linear problems are used to reach a prescribed accuracy. This goal is best met with an h-adaptive scheme in combination with an error estimator. In an industrial context, the aim of the mechanical simulations in engineering design is not only to obtain greatest quality but more often a compromise between the desired quality and the computation cost (CPU time, storage, software, competence, human cost, computer used). In this paper, we propose the use of alternative mesh refinement criteria with an h-adaptive procedure for 3D elastic problems. The alternative mesh refinement criteria (MR) are based on: prescribed number of elements with maximum accuracy, prescribed CPU time with maximum accuracy and prescribed memory size with maximum accuracy. These adaptive strategies are based on a technique of error in constitutive relation (the process could be used with other error estimators) and an efficient adaptive technique which automatically takes into account the steep gradient areas. This work proposes a 3D method of adaptivity with the latest version of the INRIA automatic mesh generator GAMHIC3D.     

热传导-对流问题的回溯二重水平法

该文给出定常的热传导-对流问题的有限元逼近的一种二重水平方法. 这种二重水平方法包括解一个小的非线性的粗网格系统、一个细网格上的线性Oseen问题和一个粗网格上的线性校正问题. 同时,给出了这种近似解的存在性和收敛性分析.     

A second order isoparametric finite element method for elliptic interface problems

A second order isoparametric finite element method (IPFEM) is proposed for elliptic interface problems. It yields better accuracy than some existing second-order methods, when the coefficients or the flux across the immersed curved interface is discontinuous. Based on an initial Cartesian mesh, a mesh optimization strategy is presented by employing curved boundary elements at the interface, and an incomplete quadratic finite element space is constructed on the optimized mesh. It turns out that the number of curved boundary elements is far less than that of the straight one, and the total degree of freedom is almost the same as the uniform Cartesian mesh. Numerical examples with simple and complicated geometrical interfaces demonstrate the efficiency of the proposed method.     

Error reduction of the adaptive conforming and nonconforming finite element methods with red–green refinement

In this paper, we analyze the convergence of the adaptive conforming and nonconforming $P_1$ finite element methods with red–green refinement based on standard Dörfler marking strategy. Since the mesh after refining is not nested into the one before, the usual Galerkin-orthogonality or quasi-orthogonality for newest vertex bisection does not hold for this case. To overcome such a difficulty, we develop some new quasi-orthogonality instead under certain condition on the initial mesh (Condition A). Consequently, we show convergence of the adaptive methods by establishing the reduction of some total errors. To weaken the condition on the initial mesh, we propose a modified red–green refinement and prove the convergence of the associated adaptive methods under a much weaker condition on the initial mesh (Condition B). Furthermore, we also develop an initial mesh generator which guarantee that all the interior triangles are equilateral triangles (satisfy Condition A) and the other triangles containing at least one vertex on the boundary satisfy Condition B.     

A Lyapunov function for a two-chemical species version of the chemotaxis model

Since the accuracy of finite element solutions of partial differential equations is generally mesh dependent, especially when solutions have singularities and discontinuities, a proper mesh generation is often important and sometimes crucial for an accurate numerical approximation of such problems. In this paper, the mesh transformation method is applied to the boundary value problems of elliptic partial differential equations, and it is proved that the method leads to the optimal finite element solutions. AMS subject classification (2000) 73C50, 65K10, 65N12, 65N30     

A TWO-LEVEL FINITE ELEMENT GALERKIN METHOD FOR THE NONSTATIONARY NAVIER-STOKES EQUATIONS I： SPATIAL DISCRETIZATION

In this article we consider a two-level finite element Galerkin method using mixed finite elements for the two-dimensional nonstationary incompressible Navier-Stokes equations. The method yields a $H^1$-optimal velocity approximation and a $L_2$-optimal pressure approximation. The two-level finite element Galerkin method involves solving one small, nonlinear Navier-Stokes problem on the coarse mesh with mesh size $H$, one linear Stokes problem on the fine mesh with mesh size $h << H$. The algorithm we study produces an approximate solution with the optimal, asymptotic in $h$, accuracy.     

自适应间断Galerkin 有限元的多水平方法

本文讨论在自适应网格上间断Galerkin 有限元离散系统的局部多水平算法. 对于光滑系数和间断系数情形, 利用Schwarz 理论分析了算法的收敛性. 理论和数值试验均说明算法的收敛率与网格层数以及网格尺寸无关. 对强间断系数情形算法是拟最优的, 即收敛率仅与网格层数有关.     

A higher order finite element method with modified graded mesh for singularly perturbed two-parameter problems

This paper presents a modified graded mesh for singularly perturbed two-parameter problems. The mesh is generated recursively using Newton's algorithm and some implicitly defined function. The problem is solved numerically using the finite element method based on higher order polynomials of degree p≥1. We prove parameter uniform convergence of optimal order in ε-weighted energy norm. A test example is taken to compare the proposed graded mesh with others found in the literature.     

旋转Q1非协调元的V循环多重网格法

1．引言近年来,多重网格法已成为行之有效的偏微分方程数值解法,而对非协调元的多重网格法也有众多的研究,例在[1,3］中,作者研究了非协P1元的w循环多重网格法,[10］中,作者研究了＊11s。n非协调元的V循环多重网格法．此外在K豆1,12］中,作者研究了板问题非协调有限元的多重网格法．最近,Rannacher和Turek同构造了所谓的QI非协调元,并用该元离散StokeS问题．而在问中,利用该元来计算晶体,数值效果非常好．同时在同中,作者给出了该元的误差估计和超收敛分析．最近,Chen和oswald同又讨论了该元的多重网格法,并证明了W循环…     

Mesh-dependent stability for finite element approximations of parabolic equations with mass lumping

In this paper, based on some mesh-dependent estimates on the extreme eigenvalues of a general finite element system defined on a simplicial mesh, novel and sharp bounds on the permissible time step size are derived for the mass lumping finite element approximations of parabolic equations. The bounds are dependent not only on the mesh size but also on the mesh shape. These results provide guidance to the stability of numerical solutions of parabolic problems in relation to the unstructured geometric meshing. Numerical experiments on both uniform meshes and adaptive meshes are presented to validate the theoretical analysis.     

Graded Meshes for Higher Order FEM

A singularly perturbed one-dimensional convection-diffusion problem is solved numerically by the finite element method based on higher order polynomials. Numerical solutions are obtained using S-type meshes with special emphasis on meshes which are graded (based on a mesh generating function) in the fine mesh region. Error estimates in the ε-weighted energy norm are proved. We derive an 'optimal' mesh generating function in order to minimize the constant in the error estimate. Two layer-adapted meshes defined by a recursive formulae in the fine mesh region are also considered and a new technique for proving error estimates for these meshes is presented. The aim of the paper is to emphasize the importance of using optimal meshes for higher order finite element methods. Numerical experiments support all theoretical results.     

Application of Nelder-Mead simplex method for unconfined seepage problems

In unconfined seepage problems, the phreatic line resulted from mesh deforming methods is rarely a smooth and continuous curve. The main problem is at the meeting point of the phreatic line with the down stream face of the dam where the phreatic line must be tangent to the seepage face according to the fluid continuity principle. In this paper a mesh deforming finite element method based on Nelder-Mead simplex optimization is presented to solve this problem. The phreatic line is approximated by a 4th degree polynomial and Nelder-Mead simplex method is used to calculate the polynomial’s coefficients minimizing an error function which is introduced based on the conditions on the phreatic line. Tangentiality of the phreatic line to the seepage face is introduced in the solution by a constraint in optimization procedure. The results of the presented method are verified by the results of the nonlinear finite element and other mesh deforming methods.     

Two-level pressure projection finite element methods for Navier-Stokes equations with nonlinear slip boundary conditions

The two-level pressure projection stabilized finite element methods for Navier-Stokes equations with nonlinear slip boundary conditions are investigated in this paper, whose variational formulation is the Navier-Stokes type variational inequality problem of the second kind. Based on the P₁-P₁ triangular element and using the pressure projection stabilized finite element method, we solve a small Navier-Stokes type variational inequality problem on the coarse mesh with mesh size H and solve a large Stokes type variational inequality problem for simple iteration or a large Oseen type variational inequality problem for Oseen iteration on the fine mesh with mesh size h. The error analysis obtained in this paper shows that if h=O(H²), the two-level stabilized methods have the same convergence orders as the usual one-level stabilized finite element methods, which is only solving a large Navier-Stokes type variational inequality problem on the fine mesh. Finally, numerical results are given to verify the theoretical analysis.     

Convergence analysis of the adaptive finite element method with the red-green refinement

In this paper, we analyze the convergence of the adaptive conforming P 1 element method with the red-green refinement. Since the mesh after refining is not nested into the one before, the Galerkin-orthogonality does not hold for this case. To overcome such a difficulty, we prove some quasi-orthogonality instead under some mild condition on the initial mesh (Condition A). Consequently, we show convergence of the adaptive method by establishing the reduction of some total error. To weaken the condition on the...